Complete nucleotide sequence and analysis of pPSR1 (72,601 bp), a pPT23A-family plasmid from Pseudomonas syringae pv. syringae A2

On the use of antibiotics to control plant pathogenic bacteria: a genetic and genomic perspective

Despite growing attention, antibiotics (such as streptomycin, oxytetracycline or kasugamycin) are still used worldwide for the control of major bacterial plant diseases. This raises concerns on their potential, yet unknown impact on antibiotic and multidrug resistances and the spread of their genetic determinants among bacterial pathogens. Antibiotic resistance genes (ARGs) have been identified in plant pathogenic bacteria (PPB), with streptomycin resistance genes being the most commonly reported. Therefore, the contribution of mobile genetic elements (MGEs) to their spread among PPB, as well as their ability to transfer to other bacteria, need to be further explored. The only well-documented example of ARGs vector in PPB, Tn5393 and its highly similar variants (carrying streptomycin resistance genes), is concerning because of its presence outside PPB, in Salmonella enterica and Klebsiella pneumoniae, two major human pathogens. Although its structure among PPB is still relatively simple, in human- and animal-associated bacteria, Tn5393 has evolved into complex associations with other MGEs and ARGs. This review sheds light on ARGs and MGEs associated with PPB, but also investigates the potential role of antibiotic use in resistance selection in plant-associated bacteria.

Cooperative virulence via the collective action of secreted pathogen effectors

Although virulence is typically attributed to single pathogenic strains, here we investigated whether effectors secreted by a population of non-virulent strains could function as public goods to enable the emergence of collective virulence. We disaggregated the 36 type III effectors of the phytopathogenic bacterium Pseudomonas syringae strain PtoDC3000 into a 'metaclone' of 36 coisogenic strains, each carrying a single effector in an effectorless background. Each coisogenic strain was individually unfit, but the metaclone was collectively as virulent as the wild-type strain on Arabidopsis thaliana, suggesting that effectors can drive the emergence of cooperation-based virulence through their public action. We show that independently evolved effector suits can equally drive this cooperative behaviour by transferring the effector alleles native to the strain PmaES4326 into the conspecific but divergent strain PtoDC3000. Finally, we transferred the disaggregated PtoDC3000 effector arsenal into Pseudomonas fluorescens and show that their cooperative action was sufficient to convert this rhizosphere-inhabiting beneficial bacterium into a phyllosphere pathogen. These results emphasize the importance of microbial community interactions and expand the ecological scale at which disease may be attributed.

The toxic guardians - multiple toxin-antitoxin systems provide stability, avoid deletions and maintain virulence genes of Pseudomonas syringae virulence plasmids

Background

Pseudomonas syringae is a γ-proteobacterium causing economically relevant diseases in practically all cultivated plants. Most isolates of this pathogen contain native plasmids collectively carrying many pathogenicity and virulence genes. However, P. syringae is generally an opportunistic pathogen primarily inhabiting environmental reservoirs, which could exert a low selective pressure for virulence plasmids. Additionally, these plasmids usually contain a large proportion of repeated sequences, which could compromise plasmid integrity. Therefore, the identification of plasmid stability determinants and mechanisms to preserve virulence genes is essential to understand the evolution of this pathogen and its adaptability to agroecosystems.

Results

The three virulence plasmids of P. syringae pv. savastanoi NCPPB 3335 contain from one to seven functional stability determinants, including three highly active toxin-antitoxin systems (TA) in both pPsv48A and pPsv48C. The TA systems reduced loss frequency of pPsv48A by two orders of magnitude, whereas one of the two replicons of pPsv48C likely confers stable inheritance by itself. Notably, inactivation of the TA systems from pPsv48C exposed the plasmid to high-frequency deletions promoted by mobile genetic elements. Thus, recombination between two copies of MITEPsy2 caused the deletion of an 8.3 kb fragment, with a frequency of 3.8 ± 0.3 × 10^− 3. Likewise, one-ended transposition of IS801 generated plasmids containing deletions of variable size, with a frequency of 5.5 ± 2.1 × 10^− 4, of which 80% had lost virulence gene idi. These deletion derivatives were stably maintained in the population by replication mediated by repJ, which is adjacent to IS801. IS801 also promoted deletions in plasmid pPsv48A, either by recombination or one-ended transposition. In all cases, functional TA systems contributed significantly to reduce the occurrence of plasmid deletions in vivo.

Conclusions

Virulence plasmids from P. syringae harbour a diverse array of stability determinants with a variable contribution to plasmid persistence. Importantly, we showed that multiple plasmid-borne TA systems have a prominent role in preserving plasmid integrity and ensuring the maintenance of virulence genes in free-living conditions. This strategy is likely widespread amongst native plasmids of P. syringae and other bacteria.

Electronic supplementary material

The online version of this article (10.1186/s13100-019-0149-4) contains supplementary material, which is available to authorized users.

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.

Plasmid Replicons from Pseudomonas Are Natural Chimeras of Functional, Exchangeable Modules

Plasmids are a main factor for the evolution of bacteria through horizontal gene exchange, including the dissemination of pathogenicity genes, resistance to antibiotics and degradation of pollutants. Their capacity to duplicate is dependent on their replication determinants (replicon), which also define their bacterial host range and the inability to coexist with related replicons. We characterize a second replicon from the virulence plasmid pPsv48C, from Pseudomonas syringae pv. savastanoi, which appears to be a natural chimera between the gene encoding a newly described replication protein and a putative replication control region present in the widespread family of PFP virulence plasmids. We present extensive evidence of this type of chimerism in structurally similar replicons from species of Pseudomonas, including environmental bacteria as well as plant, animal and human pathogens. We establish that these replicons consist of two functional modules corresponding to putative control (REx-C module) and replication (REx-R module) regions. These modules are functionally separable, do not show specificity for each other, and are dynamically exchanged among replicons of four distinct plasmid families. Only the REx-C module displays strong incompatibility, which is overcome by a few nucleotide changes clustered in a stem-and-loop structure of a putative antisense RNA. Additionally, a REx-C module from pPsv48C conferred replication ability to a non-replicative chromosomal DNA region containing features associated to replicons. Thus, the organization of plasmid replicons as independent and exchangeable functional modules is likely facilitating rapid replicon evolution, fostering their diversification and survival, besides allowing the potential co-option of appropriate genes into novel replicons and the artificial construction of new replicon specificities.

Comparative genomics of Pseudomonas syringae pv. syringae strains B301D and HS191 and insights into intrapathovar traits associated with plant pathogenesis

Pseudomonas syringae pv. syringae is a common plant-associated bacterium that causes diseases of both monocot and dicot plants worldwide. To help delineate traits critical to adaptation and survival in the plant environment, we generated complete genome sequences of P. syringae pv. syringae strains B301D and HS191, which represent dicot and monocot strains with distinct host specificities. Intrapathovar comparisons of the B301D (6.09 Mb) and HS191 (5.95 Mb plus a 52 kb pCG131 plasmid) genomes to the previously sequenced B728a genome demonstrated that the shared genes encompass about 83% of each genome, and include genes for siderophore biosynthesis, osmotolerance, and extracellular polysaccharide production. Between 7% and 12% of the genes are unique among the genomes, and most of the unique gene regions carry transposons, phage elements, or IS elements associated with horizontal gene transfer. Differences are observed in the type III effector composition for the three strains that likely influences host range. The HS191 genome had the largest number at 25 of effector genes, and seven effector genes are specific to this monocot strain. Toxin production is another major trait associated with virulence of P. syringae pv. syringae, and HS191 is distinguished by genes for production of syringopeptin SP25 and mangotoxin.

Two types of genetic carrier, the IncP genomic island and the novel IncP-1β plasmid, for the aac(2')-IIa gene that confers kasugamycin resistance in Acidovorax avenae ssp. avenae

A unique aminoglycoside antibiotic, kasugamycin (KSM), has been used to control many plant bacterial and fungal diseases in several countries. The emergence of KSM-resistant Acidovorax avenae ssp. avenae and Burkholderia glumae, which cause rice bacterial brown stripe and rice bacterial grain and seedling rot, respectively, is a serious threat for the effective control of these diseases. Previously, we have identified the aac(2')-IIa gene, encoding a KSM 2'-N-acetyltransferase, from both KSM-resistant pathogens. Although all KSM-resistant isolates from both species possess the aac(2')-IIa gene, only A. avenae strain 83 showed higher resistance than other strains. In this research, kinetic analysis indicates that an amino acid substitution from serine to threonine at position 146 of AAC(2')-IIa in strain 83 is not involved in this increased resistance. Whole draft genome analysis of A. avenae 83 shows that the aac(2')-IIa gene is carried by the novel IncP-1β plasmid pAAA83, whereas the genetic carrier of other strains, the IncP genomic island, is inserted into their chromosomes. The difference in the nucleotides of the promoter region of aac(2')-IIa between strain 83 and other strains indicates an additional transcription start site and results in the increased transcription of aac(2')-IIa in strain 83. Moreover, biological characterization of pAAA83 demonstrates that it can be transferred by conjugation and maintained in the host cells. These results demonstrate that acquisition of the aac(2')-IIa gene takes place in at least two ways and that the gene module, which includes aac(2')-IIa and the downstream gene, may be an important unit for the dissemination of antibiotic resistance.

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 rulB gene of plasmid pWW0 is a hotspot for the site-specific insertion of integron-like elements found in the chromosomes of environmental Pseudomonas fluorescens group bacteria

The rulAB operon of Pseudomonas spp. confers fitness traits on the host and has been suggested to be a hotspot for insertion of mobile elements that carry avirulence genes. Here, for the first time, we show that rulB on plasmid pWW0 is a hotspot for the active site-specific integration of related integron-like elements (ILEs) found in six environmental pseudomonads (strains FH1–FH6). Integration into rulB on pWW0 occurred at position 6488 generating a 3 bp direct repeat. ILEs from FH1 and FH5 were 9403 bp in length and contained eight open reading frames (ORFs), while the ILE from FH4 was 16 233 bp in length and contained 16 ORFs. In all three ILEs, the first 5.1 kb (containing ORFs 1–4) were structurally conserved and contained three predicted site-specific recombinases/integrases and a tetR homologue. Downstream of these resided ORFs of the ‘variable side’ with structural and sequence similarity to those encoding survival traits on the fitness enhancing plasmid pGRT1 (ILE_FH1 and ILE_FH5) and the NR-II virulence region of genomic island PAGI-5 (ILE_FH4). Collectively, these ILEs share features with the previously described type III protein secretion system effector ILEs and are considered important to host survival and transfer of fitness enhancing and (a)virulence genes between bacteria.

pA506, a conjugative plasmid of the plant epiphyte Pseudomonas fluorescens A506

Conjugative plasmids are known to facilitate the acquisition and dispersal of genes contributing to the fitness of Pseudomonas spp. Here, we report the characterization of pA506, the 57-kb conjugative plasmid of Pseudomonas fluorescens A506, a plant epiphyte used in the United States for the biological control of fire blight disease of pear and apple. Twenty-nine of the 67 open reading frames (ORFs) of pA506 have putative functions in conjugation, including a type IV secretion system related to that of MOBP6 family plasmids and a gene cluster for type IV pili. We demonstrate that pA506 is self-transmissible via conjugation between A506 and strains of Pseudomonas spp. or the Enterobacteriaceae. The origin of vegetative replication (oriV) of pA506 is typical of those in pPT23A family plasmids, which are present in many pathovars of Pseudomonas syringae, but pA506 lacks repA, a defining locus for pPT23A plasmids, and has a novel partitioning region. We selected a plasmid-cured derivative of A506 and compared it to the wild type to identify plasmid-encoded phenotypes. pA506 conferred UV resistance, presumably due to the plasmid-borne rulAB genes, but did not influence epiphytic fitness of A506 on pear or apple blossoms in the field. pA506 does not appear to confer resistance to antibiotics or other toxic elements. Based on the conjugative nature of pA506 and the large number of its genes that are shared with plasmids from diverse groups of environmental bacteria, the plasmid is likely to serve as a vehicle for genetic exchange between A506 and its coinhabitants on plant surfaces.

Sequence determination and analysis of three plasmids of Pseudomonas sp. GLE121, a psychrophile isolated from surface ice of Ecology Glacier (Antarctica)

Pseudomonas sp. GLE121 (a psychrophilic Antarctic strain) carries three plasmids: pGLE121P1 (6899 bp), pGLE121P2 (8330 bp) and pGLE121P3 (39,583 bp). Plasmids pGLE121P1 and pGLE121P2 show significant sequence similarity to members of the IncP-9 and IncP-7 incompatibility groups, respectively, while the largest replicon, pGLE121P3, is highly related to plasmid pNCPPB880-40 of Pseudomonas syringae pathovar tomato NCPPB880. All three plasmids have a narrow host range, limited to members of the genus Pseudomonas. Plasmid pGLE121P3 encodes a conjugal transfer system, while pGLE121P1 carries only a putative MOB module, conserved in many mobilizable plasmids. Plasmid pGLE121P3 contains an additional load of genetic information, including a pair of genes with homology to the rulAB operon, responsible for ultraviolet radiation (UVR) tolerance. Given the increasing UV exposure in Antarctic regions, the expression of these genes is likely to be an important adaptive response.

Genetic analysis of streptomycin-resistant (Sm(R)) strains of Erwinia amylovora suggests that dissemination of two genotypes is responsible for the current distribution of Sm(R) E. amylovora in Michigan

Streptomycin-resistant (Sm(R)) strains of the fire blight pathogen Erwinia amylovora were first isolated in southwest Michigan in 1991. Since that time, resistant strains have progressed northward to other apple-producing regions in the state. A total of 98.7% of Sm(R) strains isolated between 2003 and 2009 in Michigan harbored the strA-strB genes on transposon Tn5393. strA and strB encode phosphotransferase enzymes that modify streptomycin to a nonbactericidal form. Mutational resistance to streptomycin, caused by a point mutation-mediated target-site alteration of the ribosomal S12 protein, occurred in 1.3% of E. amylovora strains from Michigan. Tn5393 was originally introduced to E. amylovora on the plasmid pEa34; thus, the first Sm(R) strains isolated contained both pEa34 and the ubiquitous nonconjugative plasmid pEA29. More recently, we have observed Sm(R) strains in which Tn5393 is present on pEA29, suggesting that the transposon has moved via transposition from pEa34 to pEA29. Almost all of the strains containing Tn5393 on pEA29 had lost pEa34. Of 210 pEA29::Tn5393 plasmids examined, the transposon was inserted at either nucleotide position 1,515 or 17,527. Both of these positions were in noncoding regions of pEA29. Comparative sequencing of the housekeeping genes groEL and potentially variable sequences on pEA29 was done in an attempt to genetically distinguish Sm(R) strains from streptomycin-sensitive (Sm(S)) strains isolated in Michigan. Only 1 nucleotide difference within the total 2,660 bp sequenced from each strain was observed in 2 of 29 strains; multiple sequence differences were observed between the Michigan strains and E. amylovora control strains isolated in the western United States or from Rubus spp. Alterations in virulence observable using an immature pear fruit assay were detected in three of eight Sm(R) strains examined. Our current genetic data indicate that only two Sm(R) strain genotypes (strains containing pEA29::Tn5393 with Tn5393 inserted at either nucleotide position 1,515 or 17,527 on the plasmid) are responsible for the dissemination of Tn5393-encoded streptomycin resistance in Michigan, and that the Sm(R) and Sm(S) strains in Michigan compose a homogenous group.

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

Bacterial blight, caused by Pseudomonas savastanoi pv. glycinea (Psg), is a common disease of soybean. In an effort to compare a current field isolate with one isolated in the early 1960s, the genomes of two Psg strains, race 4 and B076, were sequenced using 454 pyrosequencing. The genomes of both Psg strains share more than 4,900 highly conserved genes, indicating very low genetic diversity between Psg genomes. Though conserved, genome rearrangements and recombination events occur commonly within the two Psg genomes. When compared to each other, 437 and 163 specific genes were identified in B076 and race 4, respectively. Most specific genes are plasmid-borne, indicating that acquisition and maintenance of plasmids may represent a major mechanism to change the genetic composition of the genome and even acquire new virulence factors. Type three secretion gene clusters of Psg strains are near identical with that of P. savastanoi pv. phaseolicola (Pph) strain 1448A and they shared 20 common effector genes. Furthermore, the coronatine biosynthetic cluster is present on a large plasmid in strain B076, but not in race 4. In silico subtractive hybridization-based comparative genomic analyses with nine sequenced phytopathogenic pseudomonads identified dozens of specific islands (SIs), and revealed that the genomes of Psg strains are more similar to those belonging to the same genomospecies such as Pph 1448A than to other phytopathogenic pseudomonads. The number of highly conserved genes (core genome) among them decreased dramatically when more genomes were included in the subtraction, suggesting the diversification of pseudomonads, and further indicating the genome heterogeneity among pseudomonads. However, the number of specific genes did not change significantly, suggesting these genes are indeed specific in Psg genomes. These results reinforce the idea of a species complex of P. syringae and support the reclassification of P. syringae into different species.

The Pseudomonas syringae effector protein HopZ1a suppresses effector-triggered immunity

*The Pseudomonas syringae pv syringae type III effector HopZ1a is a member of the HopZ effector family of cysteine-proteases that triggers immunity in Arabidopsis. This immunity is dependent on HopZ1a cysteine-protease activity, and independent of known resistance genes. We have previously shown that HopZ1a-triggered immunity is partially additive to that triggered by AvrRpt2. These partially additive effects could be caused by at least two mechanisms: their signalling pathways share a common element(s), or one effector interferes with the response triggered by the other. *Here, we investigate the molecular basis for the partially additive effect displayed by AvrRpt2- and HopZ1a-triggered immunities, by analysing competitive indices, hypersensitive response and symptom induction, PR-1 accumulation, expression of PR genes, and systemic acquired resistance (SAR) induction. *Partially additive effects between these defence responses require HopZ1a cysteine-protease activity, and also take place between HopZ1a and AvrRps4 or AvrRpm1-triggered responses. We establish that HopZ1a-triggered immunity is independent of salicylic acid (SA), EDS1, jasmonic acid (JA) and ethylene (ET)-dependent pathways, and show that HopZ1a suppresses the induction of PR-1 and PR-5 associated with P. syringae pv tomato (Pto)-triggered effector-triggered immunity (ETI)-like defences, AvrRpt2-triggered immunity, and Pto or Pto (avrRpt2) activation of SAR, and that suppression requires HopZ1a cysteine-protease activity. *Our results indicate that HopZ1a triggers an unusual resistance independent of known pathways and suppresses SA and EDS1-dependent resistance.

Catalytic domain of the diversified Pseudomonas syringae type III effector HopZ1 determines the allelic specificity in plant hosts

The type III secretion systems (T3SS) and secreted effectors (T3SEs) are essential virulence factors in Gram-negative bacteria. During the arms race, plants have evolved resistance (R) genes to detect specific T3SEs and activate defence responses. However, this immunity can be efficiently defeated by the pathogens through effector evolution. HopZ1 of the plant pathogen Pseudomonas syringae is a member of the widely distributed YopJ T3SE family. Three alleles are known to be present in P. syringae, with HopZ1a most resembling the ancestral allelic form. In this study, molecular mechanisms underlying the sequence diversification-enabled HopZ1 allelic specificity is investigated. Using domain shuffling experiments, we present evidence showing that a central domain upstream of the conserved catalytic cysteine residue determines HopZ1 recognition specificity. Random and targeted mutagenesis identified three amino acids involved in HopZ1 allelic specificity. Particularly, the exchange of cysteine141 in HopZ1a with lysine137 at the corresponding position in HopZ1b abolished HopZ1a recognition in soybean. This position is under strong positive selection, suggesting that the cysteine/lysine mutation might be a key step driving the evolution of HopZ1. Our data support a model in which sequence diversification imposed by the plant R gene-associated immunity has driven HopZ1 evolution by allowing allele-specific substrate-binding.

Allelic variants of the Pseudomonas syringae type III effector HopZ1 are differentially recognized by plant resistance systems

The bacterial plant pathogen Pseudomonas syringae depends on the type III secretion system and type III-secreted effectors to cause disease in plants. HopZ is a diverse family of type III effectors widely distributed in P. syringae isolates. Among the HopZ homologs, HopZ1 is ancient to P. syringae and has been shown to be under strong positive selection driven by plant resistance-imposed selective pressure. Here, we characterized the virulence and avirulence functions of the three HopZ1 alleles in soybean and Nicotiana benthamiana. In soybean, HopZ1 alleles have distinct functions: HopZ1a triggers defense response, HopZ1b promotes bacterial growth, and HopZ1c has no observable effect. In N. benthamiana, HopZ1a and HopZ1b both induce plant defense responses. However, they appear to trigger different resistance pathways, evidenced by two major differences between HopZ1a- and HopZ1b-triggered hypersensitive response (HR): i) the putative N-acylation sites had no effect on HopZ1a-triggered cell death, whereas it greatly enhanced HopZ1b-triggered cell death; and ii) the HopZ1b-triggered HR, but not the HopZ1a-triggered HR, was suppressed by another HopZ homolog, HopZ3. We previously demonstrated that HopZ1a most resembled the ancestral allelic form of HopZ1; therefore, this new evidence suggested that differentiated resistance systems have evolved in plant hosts to adapt to HopZ1 diversification in P. syringae.

Identification of new type III effectors and analysis of the plant response by competitive index

In recent years, many efforts have been directed towards the identification of new type III-secreted effectors, and the completion of the secretomes of several Pseudomonas syringae pathovars. Several functional and bioinformatic screenings have been used to search for candidates, which have been tested for translocation into the plant cell, an essential criterion for the identification of new type III effector proteins. The most common translocation assay is based on the use of DeltaAvrRpt2 as a reporter. When fused to a type III effector protein, DeltaAvrRpt2 is translocated and elicits a hypersensitive response in leaves of Arabidopsis thaliana expressing the RPS2 resistance protein. This approach has been used widely and has allowed the identification of a considerable number of new effectors in a fast and convenient manner. However, as the hypersensitive response is a semi-quantitative assay, and the conditions do not resemble those occurring in nature, effectors with low expression or translocation efficiency could fail to translocate sufficient DeltaAvrRpt2 to trigger a clear hypersensitive response. In keeping with these limitations, this test has failed to detect some true effectors that have been confirmed as such by other means. In order to increase the sensitivity of this method, we have developed a modification of the DeltaAvrRpt2-based translocation assay using a competitive index in mixed infection to monitor the limitation of growth associated with the induction of the hypersensitive response. We have tested several effector candidates from P. syringae pv. phaseolicola and other P. syringae pathovars, and have compared the results obtained by our competitive index translocation assay with those obtained by standard hypersensitive response assays. We have identified six type III secretion system-translocated proteins using this approach, five of which failed to be identified by hypersensitive response assays. In addition, we have analysed the defence response triggered by one of these effectors using competitive index assays.

De novo assembly using low-coverage short read sequence data from the rice pathogen Pseudomonas syringae pv. oryzae

We developed a novel approach for de novo genome assembly using only sequence data from high-throughput short read sequencing technologies. By combining data generated from 454 Life Sciences (Roche) and Illumina (formerly known as Solexa sequencing) sequencing platforms, we reliably assembled genomes into large scaffolds at a fraction of the traditional cost and without use of a reference sequence. We applied this method to two isolates of the phytopathogenic bacteria Pseudomonas syringae. Sequencing and reassembly of the well-studied tomato and Arabidopsis pathogen, Pto(DC3000), facilitated development and testing of our method. Sequencing of a distantly related rice pathogen, Por(1_)(6), demonstrated our method's efficacy for de novo assembly of novel genomes. Our assembly of Por(1_6) yielded an N50 scaffold size of 531,821 bp with >75% of the predicted genome covered by scaffolds over 100,000 bp. One of the critical phenotypic differences between strains of P. syringae is the range of plant hosts they infect. This is largely determined by their complement of type III effector proteins. The genome of Por(1_6) is the first sequenced for a P. syringae isolate that is a pathogen of monocots, and, as might be predicted, its complement of type III effectors differs substantially from the previously sequenced isolates of this species. The genome of Por(1_6) helps to define an expansion of the P. syringae pan-genome, a corresponding contraction of the core genome, and a further diversification of the type III effector complement for this important plant pathogen species.

Global genomic analysis of Pseudomonas savastanoi pv. savastanoi plasmids

Pseudomonas savastanoi pv. savastanoi strains harbor native plasmids belonging to the pPT23A plasmid family (PFPs) which are detected in all pathovars of the related species Pseudomonas syringae examined and contribute to the ecological and pathogenic fitness of their host. However, there is a general lack of information about the gene content of P. savastanoi pv. savastanoi plasmids and their role in the interaction of this pathogen with olive plants. We designed a DNA macroarray containing 135 plasmid-borne P. syringae genes to conduct a global genetic analysis of 32 plasmids obtained from 10 P. savastanoi pv. savastanoi strains. Hybridization results revealed that the number of PFPs per strain varied from one to four. Additionally, most strains contained at least one plasmid (designated non-PFP) that did not hybridize to the repA gene of pPT23A. Only three PFPs contained genes involved in the biosynthesis of the virulence factor indole-3-acetic acid (iaaM, iaaH, and iaaL). In contrast, ptz, a gene involved in the biosynthesis of cytokinins, was found in five PFPs and one non-PFP. Genes encoding a type IV secretion system (T4SS), type IVA, were found in both PFPs and non-PFPs; however, type IVB genes were found only on PFPs. Nine plasmids encoded both T4SSs, whereas seven other plasmids carried none of these genes. Most PFPs and non-PFPs hybridized to at least one putative type III secretion system effector gene and to a variety of additional genes encoding known P. syringae virulence factors and one or more insertion sequence transposase genes. These results indicate that non-PFPs may contribute to the virulence and fitness of the P. savastanoi pv. savastanoi host. The overall gene content of P. savastanoi pv. savastanoi plasmids, with their repeated information, mosaic arrangement, and insertion sequences, suggests a possible role in adaptation to a changing environment.

Genomics of IncP-1 antibiotic resistance plasmids isolated from wastewater treatment plants provides evidence for a widely accessible drug resistance gene pool

The dramatic spread of antibiotic resistance is a crisis in the treatment of infectious diseases that affect humans. Several studies suggest that wastewater treatment plants (WWTP) are reservoirs for diverse mobile antibiotic resistance elements. This review summarizes findings derived from genomic analysis of IncP-1 resistance plasmids isolated from WWTP bacteria. Plasmids that belong to the IncP-1 group are self-transmissible, and transfer to and replicate in a wide range of hosts. Their backbone functions are described with respect to their impact on vegetative replication, stable maintenance and inheritance, mobility and plasmid control. Accessory genetic modules, mainly representing mobile genetic elements, are integrated in-between functional plasmid backbone modules. These elements carry determinants conferring resistance to nearly all clinically relevant antimicrobial drug classes, to heavy metals, and quaternary ammonium compounds used as disinfectants. All plasmids analysed here contain integrons that potentially facilitate integration, exchange and dissemination of resistance gene cassettes. Comparative genomics of accessory modules located on plasmids from WWTP and corresponding modules previously identified in other bacterial genomes revealed that animal, human and plant pathogens and other bacteria isolated from different habitats share a common pool of resistance determinants.

Plasmids and rickettsial evolution: insight from Rickettsia felis

BACKGROUND

The genome sequence of Rickettsia felis revealed a number of rickettsial genetic anomalies that likely contribute not only to a large genome size relative to other rickettsiae, but also to phenotypic oddities that have confounded the categorization of R. felis as either typhus group (TG) or spotted fever group (SFG) rickettsiae. Most intriguing was the first report from rickettsiae of a conjugative plasmid (pRF) that contains 68 putative open reading frames, several of which are predicted to encode proteins with high similarity to conjugative machinery in other plasmid-containing bacteria.

METHODOLOGY/PRINCIPAL FINDINGS

Using phylogeny estimation, we determined the mode of inheritance of pRF genes relative to conserved rickettsial chromosomal genes. Phylogenies of chromosomal genes were in agreement with other published rickettsial trees. However, phylogenies including pRF genes yielded different topologies and suggest a close relationship between pRF and ancestral group (AG) rickettsiae, including the recently completed genome of R. bellii str. RML369-C. This relatedness is further supported by the distribution of pRF genes across other rickettsiae, as 10 pRF genes (or inactive derivatives) also occur in AG (but not SFG) rickettsiae, with five of these genes characteristic of typical plasmids. Detailed characterization of pRF genes resulted in two novel findings: the identification of oriV and replication termination regions, and the likelihood that a second proposed plasmid, pRFdelta, is an artifact of the original genome assembly.

CONCLUSION/SIGNIFICANCE

Altogether, we propose a new rickettsial classification scheme with the addition of a fourth lineage, transitional group (TRG) rickettsiae, that is unique from TG and SFG rickettsiae and harbors genes from possible exchanges with AG rickettsiae via conjugation. We offer insight into the evolution of a plastic plasmid system in rickettsiae, including the role plasmids may have played in the acquirement of virulence traits in pathogenic strains, and the likely origin of plasmids within the rickettsial tree.

Phylogenetic analysis of the pPT23A plasmid family of Pseudomonas syringae

The pPT23A plasmid family of Pseudomonas syringae contains members that contribute to the ecological and pathogenic fitness of their P. syringae hosts. In an effort to understand the evolution of these plasmids and their hosts, we undertook a comparative analysis of the phylogeny of plasmid genes and that of conserved chromosomal genes from P. syringae. In total, comparative sequence and phylogenetic analyses were done utilizing 47 pPT23A family plasmids (PFPs) from 16 pathovars belonging to six genomospecies. Our results showed that the plasmid replication gene (repA), the only gene currently known to be distributed among all the PFPs, had a phylogeny that was distinct from that of the P. syringae hosts of these plasmids and from those of other individual genes on PFPs. The phylogenies of two housekeeping chromosomal genes, those for DNA gyrase B subunit (gyrB) and primary sigma factor (rpoD), however, were strongly associated with genomospecies of P. syringae. Based on the results from this study, we conclude that the pPT23A plasmid family represents a dynamic genome that is mobile among P. syringae pathovars.

Genomic insights into the contribution of phytopathogenic bacterial plasmids to the evolutionary history of their hosts

Plasmids are common residents of phytopathogenic bacteria and contribute significantly to host evolution in a multi-faceted manner. Plasmids tend to encode determinants of virulence and ecological fitness that can enhance adaptation to a specific niche or can influence niche expansion. Many of these determinants appear to have been acquired from other bacteria via horizontal transfer, illustrating an important function of plasmids in the acquisition of sequences that enable rapid evolution. These genes can ultimately be delivered to the host chromosome through plasmid integration events, thus stabilizing important acquired determinants within the genome. Most plasmids characterized in phytopathogenic bacteria are self-transmissible and possess suites of genes encoding type IV secretion systems. In addition, the phytopathogenic bacterial plasmid "mobilome" includes insertion sequence and other transposable elements that contribute to the movement of sequences within and between genomes. Possession of mosaic and ever-changing plasmids allows phytopathogenic bacteria to maintain a dynamic, flexible genome and possible advantage in host-pathogen and other environmental interactions that belies the concept of plasmids as apparently selfish genetic elements.

Type III effector diversification via both pathoadaptation and horizontal transfer in response to a coevolutionary arms race

The concept of the coevolutionary arms race holds a central position in our understanding of pathogen–host interactions. Here we identify the molecular mechanisms and follow the stepwise progression of an arms race in a natural system. We show how the evolution and function of the HopZ family of type III secreted effector proteins carried by the plant pathogen Pseudomonas syringae are influenced by a coevolutionary arms race between pathogen and host. We surveyed 96 isolates of P. syringae and identified three homologs (HopZ1, HopZ2, and HopZ3) distributed among ∼45% of the strains. All alleles were sequenced and their expression was confirmed. Evolutionary analyses determined that the diverse HopZ1 homologs are ancestral to P. syringae, and have diverged via pathoadaptive mutational changes into three functional and two degenerate forms, while HopZ2 and HopZ3 have been brought into P. syringae via horizontal transfer from other ecologically similar bacteria. A PAML selection analysis revealed that the C terminus of HopZ1 is under strong positive selection. Despite the extensive genetic variation observed in this family, all three homologs have cysteine–protease activity, although their substrate specificity may vary. The introduction of the ancestral hopZ1 allele into strains harboring alternate alleles results in a resistance protein-mediated defense response in their respective hosts, which is not observed with the endogenous allele. These data indicate that the P. syringae HopZ family has undergone allelic diversification via both pathoadaptive mutational changes and horizontal transfer in response to selection imposed by the host defense system. This genetic diversity permits the pathogen to avoid host defenses while still maintaining a virulence-associated protease, thereby allowing it to thrive on its current host, while simultaneously impacting its host range.

Pathogens and their hosts impose reciprocal selective pressures on each other, such that the improvement of one selects for the improvement of the other. Pathogens that are able to evolve increasingly effective methods of attacking their hosts select for hosts that are able to mount increasingly effective defenses against pathogen attack. This coevolutionary interaction is commonly referred to as an arms race, or the Red Queen Principle, taken from Lewis Carroll's Through the Looking Glass, and What Alice Found There, in which Alice and the Red Queen had to run as fast as they could simply to stay in the same place. Many pathogenic bacteria rely on specialized virulence proteins, called type III secreted effectors (T3SEs), to cause disease. These proteins are injected into the cells of the host, and often act to disrupt the host defense response. This study shows how the HopZ family of T3SEs in the pathogen Pseudomonas syringae evolves in response to coevolutionary selective pressures imposed by its plant hosts. The authors identify the version of the hopZ gene that is most similar to the one carried by the ancestral strain, and then show how this version has been modified by mutation and selection in response to the host defense systems. They also identify genes related to hopZ from other species that were brought into P. syringae presumably in response to this same host-imposed selective pressure. Finally, the authors show how the genetic diversity in this gene family permits the pathogen to avoid host defenses while still maintaining an important virulence-associated function. This study provides a clearer picture of the molecular interactions that drive coevolutionary interactions, and insight into how ecological processes play out at the molecular and evolutionary scale.

Anatomy of the replication origin of plasmid ColE2-P9

The plasmid ColE2-P9 origin is a 32-bp region which is specifically recognized by the plasmid-specified Rep protein to initiate DNA replication. We analyzed the structural and functional organization of the ColE2 origin by using various derivatives carrying deletions and single-base-pair substitutions. The origin may be divided into three subregions: subregion I, which is important for stable binding of the Rep protein; subregion II, which is important for binding of the Rep protein and for initiation of DNA replication; and subregion III, which is important for DNA replication but apparently not for binding of the Rep protein. The Rep protein might recognize three specific DNA elements in subregions I and II. The relative transformation frequency of the autonomously replicating plasmids carrying deletions in subregion I is lower, and nevertheless the copy numbers of these plasmids in host bacteria are higher than those of the wild-type plasmid. Efficient and stable binding of the Rep protein to the origin might be important for the replication efficiency to be at the normal (low) level. Subregion II might be essential for interaction with the catalytic domain of the Rep protein for primer RNA synthesis. The 8-bp sequence across the border of subregions II and III, including the primer sequence, is conserved in the (putative) origins of many plasmids, the putative Rep proteins of which are related to the ColE2-P9 Rep protein. Subregion III might be required for a step that is necessary after Rep protein binding has taken place.

A DNA polymerase V homologue encoded by TOL plasmid pWW0 confers evolutionary fitness on Pseudomonas putida under conditions of environmental stress

Plasmids in conjunction with other mobile elements such as transposons are major players in the genetic adaptation of bacteria in response to changes in environment. Here we show that a large catabolic TOL plasmid, pWW0, from Pseudomonas putida carries genes (rulAB genes) encoding an error-prone DNA polymerase Pol V homologue which increase the survival of bacteria under conditions of accumulation of DNA damage. A study of population dynamics in stationary phase revealed that the presence of pWW0-derived rulAB genes in the bacterial genome allows the expression of a strong growth advantage in stationary phase (GASP) phenotype of P. putida. When rulAB-carrying cells from an 8-day-old culture were mixed with Pol V-negative cells from a 1-day-old culture, cells derived from the aged culture out-competed cells from the nonaged culture and overtook the whole culture. At the same time, bacteria from an aged culture lacking the rulAB genes were only partially able to out-compete cells from a fresh overnight culture of the parental P. putida strain. Thus, in addition to conferring resistance to DNA damage, the plasmid-encoded Pol V genes significantly increase the evolutionary fitness of bacteria during prolonged nutritional starvation of a P. putida population. The results of our study indicate that RecA is involved in the control of expression of the pWW0-encoded Pol V.

Comparative genomic analysis of the pPT23A plasmid family of Pseudomonas syringae

Members of the pPT23A plasmid family of Pseudomonas syringae play an important role in the interaction of this bacterial pathogen with host plants. Complete sequence analysis of several pPT23A family plasmids (PFPs) has provided a glimpse of the gene content and virulence function of these plasmids. We constructed a macroarray containing 161 genes to estimate and compare the gene contents of 23 newly analyzed and eight known PFPs from 12 pathovars of P. syringae, which belong to four genomospecies. Hybridization results revealed that PFPs could be distinguished by the type IV secretion system (T4SS) encoded and separated into four groups. Twelve PFPs along with pPSR1 from P. syringae pv. syringae, pPh1448B from P. syringae pv. phaseolicola, and pPMA4326A from P. syringae pv. maculicola encoded a type IVA T4SS (VirB-VirD4 conjugative system), whereas 10 PFPs along with pDC3000A and pDC3000B from P. syringae pv. tomato encoded a type IVB T4SS (tra system). Two plasmids encoded both T4SSs, whereas six other plasmids carried none or only a few genes of either the type IVA or type IVB secretion system. Most PFPs hybridized to more than one putative type III secretion system effector gene and to a variety of additional genes encoding known P. syringae virulence factors. The overall gene contents of individual PFPs were more similar among plasmids within each of the four groups based on T4SS genes; however, a number of genes, encoding plasmid-specific functions or hypothetical proteins, were shared among plasmids from different T4SS groups. The only gene shared by all PFPs in this study was the repA gene, which encoded sequences with 87 to 99% amino acid identityamong 25 sequences examined. We proposed a model to illustrate the evolution and gene acquisition of the pPT23A plasmid family. To our knowledge, this is the first such attempt to conduct a global genetic analysis of this important plasmid family.

Nucleotide sequences, genetic organization, and distribution of pEU30 and pEL60 from Erwinia amylovora

The nucleotide sequences, genetic organization, and distribution of plasmids pEU30 (30,314 bp) and pEL60 (60,145 bp) from the plant pathogen Erwinia amylovora are described. The newly characterized pEU30 and pEL60 plasmids inhabited strains isolated in the western United States and Lebanon, respectively. The gene content of pEU30 resembled plasmids found in plant-associated bacteria, while that of pEL60 was most similar to IncL/M plasmids inhabiting enteric bacteria.

Lineage-specific regions in Pseudomonas syringae pv tomato DC3000

SUMMARY Comparative analyses of the chromosome of Pseudomonas syringae pv tomato DC3000 with the finished, complete genomes of Pseudomonas aeruginosa PAO1, an animal pathogen, and the non-pathogenic soil inhabitant Pseudomonas putida KT2440 revealed a high degree of sequence conservation in genes involved in 'housekeeping functions'. However, divergence is present among these three fluorescent pseudomonads, yielding 'suites' of species-specific genes that may provide the genetic basis for adaptation to an ecological niche and lifestyle. For DC3000, 1053 genes located on the chromosome were specific to DC3000 and not present in PAO1 or KT2440. The majority of these DC3000-specific genes either lack a known function or are mobile genetic elements. However, these genes do share features among themselves such as association with regions of atypical trinucleotides, unusual G+C content and localization within large tracts of DC3000-specific sequence, suggestive of lateral gene transfer events. Indeed, a comparison of syntenic blocks among these three complete Pseudomonas genomes revealed that a substantial portion (533) of the DC3000-specific chromosomal genes (1053) were located in lineage-specific regions (defined as being larger than 2 kb and enriched in mobile genetic elements and/or genes specific to DC3000 in this three-way comparison). A large proportion of mobile genetic elements (199 of 318 genes; 63%), which are highly enriched in DC3000, were present within such regions. Similarly, most of the genes encoding type III secretion system virulence effectors were located in lineage-specific regions. Consistent with the plasticity of the DC3000 genome, a putative chromosomal inversion mediated by identical copies of ISPsy6 involving 2838 kb (44%) of the DC3000 genome was detected. These data suggest that a substantial portion of the differentiation of DC3000, a plant pathogen, from an animal pathogen and a soil inhabitant has involved transfer of a large number of novel genes coupled with amplification of mobile genetic elements.

Nucleotide sequence and evolution of the five-plasmid complement of the phytopathogen Pseudomonas syringae pv. maculicola ES4326

Plasmids are transmissible, extrachromosomal genetic elements that are often responsible for environmental or host-specific adaptations. In order to identify the forces driving the evolution of these important molecules, we determined the complete nucleotide sequence of the five-plasmid complement of the radish and Arabidopsis pathogen Pseudomonas syringae pv. maculicola ES4326 and conducted an intraspecific comparative genomic analysis. To date, this is the most complex fully sequenced plasmid complement of any gram-negative bacterium. The plasmid complement comprises two pPT23A-like replicons, pPMA4326A (46,697 bp) and pPMA4326B (40,110 bp); a pPS10-like replicon, pPMA4326C (8,244 bp); and two atypical, replicase-deficient replicons, pPMA4326D (4,833 bp) and pPMA4326E (4,217 bp). A complete type IV secretion system is found on pPMA4326A, while the type III secreted effector hopPmaA is present on pPMA4326B. The region around hopPmaA includes a shorter hopPmaA homolog, insertion sequence (IS) elements, and a three-element cassette composed of a resolvase, an integrase, and an exeA gene that is also present in several human pathogens. We have also identified a novel genetic element (E622) that is present on all but the smallest plasmid (pPMA4326E) that has features of an IS element but lacks an identifiable transposase. This element is associated with virulence-related genes found in a wide range of P. syringae strains. Comparative genomic analyses of these and other P. syringae plasmids suggest a role for recombination and integrative elements in driving plasmid evolution.