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.     

Novel exchangeable effector loci associated with the Pseudomonas syringae hrp pathogenicity island: evidence for integron-like assembly from transposed gene cassettes

Pseudomonas syringae strains use a type III secretion system (TTSS) to translocate effector proteins that assist in the parasitism of host plant cells. Some genes for effector proteins are clustered in the exchangeable effector locus (EEL) associated with the hrp pathogenicity island. A polymerase chain reaction-based screen was developed to amplify the EEL from P. syringae strains. Of the 86 strains screened, the EEL was successfully amplified from 30 predominately North American P. syringae pv. syringae strains using hrpK and queA-derived primers and from an additional three strains using hrpL and queA-derived primers. Among the amplified EEL, ten distinct types of EEL were identified that could be classified into six families distinguishable by genetic composition, but other types of EEL may be present in strains isolated in other geographical regions. No linkage with the host range of the source strain was apparent. Gene cassettes carrying conserved flanking, coding, and intergenic sequences, present in different combinations, were identified in the characterized EEL. Six new alleles of known effectors were identified that differed from the homolog in sequence, size, or both of the gene. One of these apparently novel effector proteins, HopPsyB, retained a strongly conserved amino terminus similar to that of HopPsyA, but other regions of the two polypeptides were only weakly similar. hopPsyB was expressed from an apparent operon that included hrpK and a shcA homolog, shcB. Escherichia coli MC4100 expressing the hrp TTSS, ShcB, and HopPsyB elicited the hypersensitive response (HR) in tobacco, consistent with effector production. Indicative of translocation as an effector, P. syringae pv. tomato DC3000 expressing a HopPsyB':'AvrRpt2 fusion elicited the HR in RPS2+ Arabidopsis thaliana. P. syringae pv. tomato DC3000 carrying HopPsyB exhibited slightly enhanced virulence in several Brassica spp. These results are consistent with the hypotheses that the EEL is a source of disparate effectors functioning in pathogenicity of P. syringae strains and that it evolved independently of the hrp pathogenicity island central conserved region, most likely through integron-like assembly of transposed gene cassettes.     

Diversifying selection drives the evolution of the type III secretion system pilus of Pseudomonas syringae

The plant pathogenic bacterium Pseudomonas syringae uses a type III secretion system to inject virulence proteins directly into the cytoplasm of its hosts. The P. syringae type III secretion apparatus is encoded, in part, by the HrpZ operon, which carries the hrpA gene encoding the pilin subunit of the pilus, various components of the structural apparatus, and the HrpZ harpin protein that is believed to produce pores in the host cell membrane. The pilus of the type III system comes into direct contact with the host cell and is, therefore, a likely target of the host's pathogen surveillance systems. We sequenced and analyzed 22 HrpZ operons from P. syringae strains spanning the diversity of the species. Selection analyses, including K(a)/K(s) tests and Tajima's D, revealed strong diversifying selection acting on the hrpA gene. This form of selection enables pathogens to maintain genetic diversity within their populations and is often driven by selection imposed by host defense systems. The HrpZ operon also revealed a single significant recombination event that dramatically changed the evolutionary relationships among P. syringae strains from 2 quite distinct phylogroups. This recombination event appears to have introduced genetic diversity into a clade of strains that may now be undergoing positive selection. The identification of diversifying selection acting on the Hrp pilus across the whole population sample and positive selection within one P. syringae lineage supports a trench warfare coevolutionary model between P. syringae and its plant hosts.     

Changes in race-specific virulence in Pseudomonas syringae pv. phaseolicola are associated with a chimeric transposable element and rare deletion events in a plasmid-borne pathogenicity island

Virulence for bean and soybean is determined by effector genes in a plasmid-borne pathogenicity island (PAI) in race 7 strain 1449B of Pseudomonas syringae pv. phaseolicola. One of the effector genes, avrPphF, confers either pathogenicity, virulence, or avirulence depending on the plant host and is absent from races 2, 3, 4, 6, and 8 of this pathogen. Analysis of cosmid clones and comparison of DNA sequences showed that the absence of avrPphF from strain 1448A is due to deletion of a continuous 9.5-kb fragment. The remainder of the PAI is well conserved in strains 1448A and 1449B. The left junction of the deleted region consists of a chimeric transposable element generated from the fusion of homologs of IS1492 from Pseudomonas putida and IS1090 from Ralstonia eutropha. The borders of the deletion were conserved in 66 P. syringae pv. phaseolicola strains isolated in different countries and representing the five races lacking avrPphF. However, six strains isolated in Spain had a 10.5-kb deletion that extended 1 kb further from the right junction. The perfect conservation of the 28-nucleotide right repeat of the IS1090 homolog in the two deletion types and in the other 47 insertions of the IS1090 homolog in the 1448A genome strongly suggests that the avrPphF deletions were mediated by the activity of the chimeric mobile element. Our data strongly support a clonal origin for the races of P. syringae pv. phaseolicola lacking avrPphF.     

Roadmap to new virulence determinants in Pseudomonas syringae: insights from comparative genomics and genome organization

Systematic comparison of the current repertoire of virulence-associated genes for three Pseudomonas syringae strains with complete genome sequences, P. syringae pv. tomato DC3,000, P. syringae pv. phaseolicola 1448A, and P. syringae pv. syringae B728a, is prompted by recent advances in virulence factor identification in P. syringae and other bacteria. Among these are genes linked to epiphytic fitness, plant- and insect-active toxins, secretion pathways, and virulence regulators, all reflected in the recently updated DC3,000 genome annotation. Distribution of virulence genes in relation to P. syringae genome organization was analyzed to distinguish patterns of conservation among genomes and association between genes and mobile genetic elements. Variable regions were identified on the basis of deviation in sequence composition and gaps in syntenic alignment among the three genomes. Mapping gene location relative to the genome structure revealed strong segregation of the HrpL regulon with variable genome regions (VR), divergent distribution patterns for toxin genes depending on association with plant or insect pathogenesis, and patterns of distribution for other virulence genes that highlight potential sources of strain-to-strain differences in host interaction. Distribution of VR among other sequenced bacterial genomes was analyzed and future plans for characterization of this potential reservoir of virulence genes are discussed.     

Identification of a novel Pseudomonas syringae Psy61 effector with virulence and avirulence functions by a HrpL-dependent promoter-trap assay

The hrp pathogenicity island of Pseudomonas syringae encodes a type III secretion system (TTSS) that translocates effectors into plant cells. Most genes encoding effectors are dispersed in the P. syringae genome. Regardless of location, all are regulated coordinately by the alternative sigma factor HrpL. An HrpL-dependent promoter-trap assay was developed to screen genomic libraries of P. syringae strains for promoters whose activity in Escherichia coli is dependent on an inducible hrpL construct. Twenty-two HrpL-dependent promoter fragments were isolated from P. syringae Psy61 that included promoters for known HrpL-dependent genes. One fragment also was isolated that shared no similarity with known genes but retained a near consensus HrpL-dependent promoter. The sequence of the region revealed a 375-amino acid open reading frame encoding a 40.5-kDa product that was designated HopPsyL. HopPsyL was structurally similar to other secreted effectors and carried a putative chloroplast-targeting signal and two predicted transmembrane domains. HopPsyL':'AvrRpt2 fusions were translocated into host cells via the P. syringae pv. tomato DC3000 hrp TTSS. A hopPsyL::kan mutant of Psy61 exhibited strongly reduced virulence in Phaseolus vulgaris cv. Kentucky Wonder, but did not appear to act as a defense response suppressor. The ectopically expressed gene reduced the virulence of Pseudomonas syringae DC3000 transformants in Arabidopsis thaliana Col-0. The gene was shown to be conserved in 6 of 10 P. syringae pv. syringae strains but was not detected in 35 strains of other pathovars. HopPsyL appears to be a novel TTSS-dependent effector that functions as a host-species-specific virulence factor in Psy61.     

Pseudomonas syringae type III effector repertoires: last words in endless arguments

Many plant pathogens subvert host immunity by injecting compositionally diverse but functionally similar repertoires of cytoplasmic effector proteins. The bacterial pathogen Pseudomonas syringae is a model for exploring the functional structure of such repertoires. The pangenome of P. syringae encodes 57 families of effectors injected by the type III secretion system. Distribution of effector genes among phylogenetically diverse strains reveals a small set of core effectors targeting antimicrobial vesicle trafficking and a much larger set of variable effectors targeting kinase-based recognition processes. Complete disassembly of the 28-effector repertoire of a model strain and reassembly of a minimal functional repertoire reveals the importance of simultaneously attacking both processes. These observations, coupled with growing knowledge of effector targets in plants, support a model for coevolving molecular dialogs between effector repertoires and plant immune systems that emphasizes mutually-driven expansion of the components governing recognition.     

Pto- and Prf-mediated recognition of AvrPto and AvrPtoB restricts the ability of diverse pseudomonas syringae pathovars to infect tomato

The molecular basis underlying the ability of pathogens to infect certain plant species and not others is largely unknown. Pseudomonas syringae is a useful model species for investigating this phenomenon because it comprises more than 50 pathovars which have narrow host range specificities. Tomato (Solanum lycopersicum) is a host for P. syringae pv. tomato, the causative agent of bacterial speck disease, but is considered a nonhost for other P. syringae pathovars. Host resistance in tomato to bacterial speck disease is conferred by the Pto protein kinase which acts in concert with the Prf nucleotide-binding lucine-rich repeat protein to recognize P. syringae pv. tomato strains expressing the type III effectors AvrPto or AvrPtoB (HopAB2). The Pto and Prf genes were isolated from the wild tomato species S. pimpinellifolium and functional alleles of both of these genes now are known to exist in many species of tomato and in other Solanaceous species. Here, we extend earlier reports that avrPto and avrPtoB genes are widely distributed among pathovars of P. syringae which are considered nonhost pathogens of tomato. This observation prompted us to examine the possibility that recognition of these type III effectors by Pto or Prf might contribute to the inability of many P. syringae pathovars to infect tomato species. We show that 10 strains from presumed nonhost P. syringae pathovars are able to grow and cause pathovar-unique disease symptoms in tomato leaves lacking Pto or Prf, although they did not reach the population levels or cause symptoms as severe as a control P. syringae pv. tomato strain. Seven of these strains were found to express avrPto or avrPtoB. The AvrPto- and AvrPtoB-expressing strains elicited disease resistance on tomato leaves expressing Pto and Prf. Thus, a gene-for-gene recognition event may contribute to host range restriction of many P. syringae pathovars on tomato species. Furthermore, we conclude that the diverse disease symptoms caused by different Pseudomonas pathogens on their normal plant hosts are due largely to the array of virulence factors expressed by each pathovar and not to specific molecular or morphological attributes of the plant host.     

Cheating, trade-offs and the evolution of aggressiveness in a natural pathogen population

The evolutionary dynamics of pathogens are critically important for disease outcomes, prevalence and emergence. In this study we investigate ecological conditions that may promote the long-term maintenance of virulence polymorphisms in pathogen populations. Recent theory predicts that evolution towards increased virulence can be reversed if less-aggressive social 'cheats' exploit more aggressive 'cooperator' pathogens. However, there is no evidence that social exploitation operates within natural pathogen populations. We show that for the bacterium Pseudomonas syringae, major polymorphisms for pathogenicity are maintained at unexpectedly high frequencies in populations infecting the host Arabidopsis thaliana. Experiments reveal that less-aggressive strains substantially increase their growth potential in mixed infections and have a fitness advantage in non-host environments. These results suggest that niche differentiation can contribute to the maintenance of virulence polymorphisms, and that both within-host and between-host growth rates modulate cheating and cooperation in P. syringae populations.     

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.     

hrp genes of Pseudomonas solanacearum are homologous to pathogenicity determinants of animal pathogenic bacteria and are conserved among plant pathogenic bacteria.

The majority of bacterial plant diseases are caused by members of three bacterial genera, Pseudomonas, Xanthomonas, and Erwinia. The identification and characterization of mutants that have lost the abilities to provoke disease symptoms on a compatible host and to induce a defensive hypersensitive reaction (HR) on an incompatible host have led to the discovery of clusters of hrp genes (hypersensitive reaction and pathogenicity) in phytopathogenic bacteria from each of these genera. Here, we report that predicted protein sequences of three hrp genes from Pseudomonas solanacearum show remarkable sequence similarity to key virulence determinants of animal pathogenic bacteria of the genus Yersinia. We also demonstrate DNA homologies between P. solanacearum hrp genes and hrp gene clusters of P. syringae pv. phaseolicola, Xanthomonas campestris pv. campestris, and Erwinia amylovora. By comparing the role of the Yersinia determinants in the control of the extracellular production of proteins required for pathogenicity, we propose that hrp genes code for an export system that might be conserved among many diverse bacterial pathogens of plants and animals but that is distinct from the general export pathway.     

Contributions of the effector gene hopQ1-1 to differences in host range between Pseudomonas syringae pv. phaseolicola and P. syringae pv. tabaci

To study the role of type III-secreted effectors in the host adaptation of the tobacco ( Nicotiana sp.) pathogen Pseudomonas syringae pv. tabaci , a selection of seven strains was first characterized by multilocus sequence typing (MLST) to determine their phylogenetic affinity. MLST revealed that all strains represented a tight phylogenetic group and that the most closely related strain with a completely sequenced genome was the bean ( Phaseolus vulgaris ) pathogen P. syringae pv. phaseolicola 1448A. Using primers designed to 21 P. syringae pv. phaseolicola 1448A effector genes, it was determined that P. syringae pv. phaseolicola 1448A shared at least 10 effectors with all tested P. syringae pv. tabaci strains. Six of the 11 effectors that failed to amplify from P. syringae pv. tabaci strains were individually expressed in one P. syringae pv. tabaci strain. Although five effectors had no effect on phenotype, growth in planta and disease severity of the transgenic P. syringae pv. tabaci expressing hopQ1-1 _Pph1448A were significantly increased in bean, but reduced in tobacco. We conclude that hopQ1-1 has been retained in P. syringae pv. phaseolicola 1448A, as this effector suppresses immunity in bean, whereas hopQ1-1 is missing from P. syringae pv. tabaci strains because it triggers defences in Nicotiana spp. This provides evidence that fine-tuning effector repertoires during host adaptation lead to a concomitant reduction in virulence in non-host species.     

Pseudomonas syringae exchangeable effector loci: sequence diversity in representative pathovars and virulence function in P. syringae pv. syringae B728a

Pseudomonas syringae is a plant pathogen whose pathogenicity and host specificity are thought to be determined by Hop/Avr effector proteins injected into plant cells by a type III secretion system. P. syringae pv. syringae B728a, which causes brown spot of bean, is a particularly well-studied strain. The type III secretion system in P. syringae is encoded by hrp (hypersensitive response and pathogenicity) and hrc (hrp conserved) genes, which are clustered in a pathogenicity island with a tripartite structure such that the hrp/hrc genes are flanked by a conserved effector locus and an exchangeable effector locus (EEL). The EELs of P. syringae pv. syringae B728a, P. syringae strain 61, and P. syringae pv. tomato DC3000 differ in size and effector gene composition; the EEL of P. syringae pv. syringae B728a is the largest and most complex. The three putative effector proteins encoded by the P. syringae pv. syringae B728a EEL--HopPsyC, HopPsyE, and HopPsyV--were demonstrated to be secreted in an Hrp-dependent manner in culture. Heterologous expression of hopPsyC, hopPsyE, and hopPsyV in P. syringae pv. tabaci induced the hypersensitive response in tobacco leaves, demonstrating avirulence activity in a nonhost plant. Deletion of the P. syringae pv. syringae B728a EEL strongly reduced virulence in host bean leaves. EELs from nine additional strains representing nine P. syringae pathovars were isolated and sequenced. Homologs of avrPphE (e.g., hopPsyE) and hopPsyA were particularly common. Comparative analyses of these effector genes and hrpK (which flanks the EEL) suggest that the EEL effector genes were acquired by horizontal transfer after the acquisition of the hrp/hrc gene cluster but before the divergence of modern pathovars and that some EELs underwent transpositions yielding effector exchanges or point mutations producing effector pseudogenes after their acquisition.     

Comparative genomics of host-specific virulence in Pseudomonas syringae

While much study has gone into characterizing virulence factors that play a general role in disease, less work has been directed at identifying pathogen factors that act in a host-specific manner. Understanding these factors will help reveal the variety of mechanisms used by pathogens to suppress or avoid host defenses. We identified candidate Pseudomonas syringae host-specific virulence genes by searching for genes whose distribution among natural P. syringae isolates was statistically associated with hosts of isolation. We analyzed 91 strains isolated from 39 plant hosts by DNA microarray-based comparative genomic hybridization against an array containing 353 virulence-associated (VA) genes, including 53 type III secretion system effectors (T3SEs). We identified individual genes and gene profiles that were significantly associated with strains isolated from cauliflower, Chinese cabbage, soybean, rice, and tomato. We also identified specific horizontal gene acquisition events associated with host shifts by mapping the array data onto the core genome phylogeny of the species. This study provides the largest suite of candidate host-specificity factors from any pathogen, suggests that there are multiple ways in which P. syringae isolates can adapt to the same host, and provides insight into the evolutionary mechanisms underlying host adaptation.     

Evolution of microbial virulence: the benefits of stress

Although genome sequencing of microbial pathogens has shed light on the evolution of virulence, the drivers of the gain and loss of genes and of pathogenicity islands (gene clusters), which contribute to the emergence of new disease outbreaks, are unclear. Recent experiments with the bean pathogen Pseudomonas syringae pv. phaseolicola illustrate how exposure to resistance mechanisms acts as the driving force for genome reorganization. Here we argue that the antimicrobial conditions generated by host defences can accelerate the generation of genome rearrangements that provide selective advantages to the invading microbe. Similar exposure to environmental stress outside the host could also drive the horizontal gene transfer that has led to the evolution of pathogenicity towards both animals and plants.     

The Pseudomonas syringae type III effector HopAM1 enhances virulence on water-stressed plants

Pseudomonas syringae strains deliver diverse type III effector proteins into host cells, where they can act as virulence factors. Although the functions of the majority of type III effectors are unknown, several have been shown to interfere with plant basal defense mechanisms. Type III effectors also could contribute to bacterial virulence by enhancing nutrient uptake and pathogen adaptation to the environment of the host plant. We demonstrate that the type III effector HopAM1 (formerly known as AvrPpiB) enhances the virulence of a weak pathogen in plants that are grown under drought stress. This is the first report of a type III effector that aids pathogen adaptation to water availability in the host plant. Expression of HopAM1 makes transgenic Ws-0 Arabidopsis hypersensitive to abscisic acid (ABA) for stomatal closure and germination arrest. Conditional expression of HopAM1 in Arabidopsis also suppresses basal defenses. ABA responses overlap with defense responses and ABA has been shown to suppress defense against P. syringae pathogens. We propose that HopAM1 aids P. syringae virulence by manipulation of ABA responses that suppress defense responses. In addition, host ABA responses enhanced by type III delivery of HopAM1 protect developing bacterial colonies inside leaves from osmotic stress.     

Using Weakly Conserved Motifs Hidden in Secretion Signals to Identify Type-III Effectors from Bacterial Pathogen Genomes

Background

As one of the most important virulence factor types in gram-negative pathogenic bacteria, type-III effectors (TTEs) play a crucial role in pathogen-host interactions by directly influencing immune signaling pathways within host cells. Based on the hypothesis that type-III secretion signals may be comprised of some weakly conserved sequence motifs, here we used profile-based amino acid pair information to develop an accurate TTE predictor.

Results

For a TTE or non-TTE, we first used a hidden Markov model-based sequence searching method (i.e., HHblits) to detect its weakly homologous sequences and extracted the profile-based k-spaced amino acid pair composition (HH-CKSAAP) from the N-terminal sequences. In the next step, the feature vector HH-CKSAAP was used to train a linear support vector machine model, which we designate as BEAN (Bacterial Effector ANalyzer). We compared our method with four existing TTE predictors through an independent test set, and our method revealed improved performance. Furthermore, we listed the most predictive amino acid pairs according to their weights in the established classification model. Evolutionary analysis shows that predictive amino acid pairs tend to be more conserved. Some predictive amino acid pairs also show significantly different position distributions between TTEs and non-TTEs. These analyses confirmed that some weakly conserved sequence motifs may play important roles in type-III secretion signals. Finally, we also used BEAN to scan one plant pathogen genome and showed that BEAN can be used for genome-wide TTE identification. The webserver and stand-alone version of BEAN are available at http://protein.cau.edu.cn:8080/bean/.     

Closing the circle on the discovery of genes encoding Hrp regulon members and type III secretion system effectors in the genomes of three model Pseudomonas syringae strains

Pseudomonas syringae strains translocate large and distinct collections of effector proteins into plant cells via the type III secretion system (T3SS). Mutations in T3SS-encoding hrp genes are unable to elicit the hypersensitive response or pathogenesis in nonhost and host plants, respectively. Mutations in individual effectors lack strong phenotypes, which has impeded their discovery. P. syringae effectors are designated Hop (Hrp outer protein) or Avr (avirulence) proteins. Some Hop proteins are considered to be extracellular T3SS helpers acting at the plant-bacterium interface. Identification of complete sets of effectors and related proteins has been enabled by the application of bioinformatic and high-throughput experimental techniques to the complete genome sequences of three model strains: P. syringae pv. tomato DC3000, P. syringae pv. phaseolicola 1448A, and P. syringae pv. syringae B728a. Several recent papers, including three in this issue of Molecular Plant-Microbe Interactions, address the effector inventories of these strains. These studies establish that active effector genes in P. syringae are expressed by the HrpL alternative sigma factor and can be predicted on the basis of cis Hrp promoter sequences and N-terminal amino-acid patterns. Among the three strains analyzed, P. syringae pv. tomato DC3000 has the largest effector inventory and P. syringae pv. syringae B728a has the smallest. Each strain has several effector genes that appear inactive. Only five of the 46 effector families that are represented in these three strains have an active member in all of the strains. Web-based community resources for managing and sharing growing information on these complex effector arsenals should help future efforts to understand how effectors promote P. syringae virulence.