Evolution, genomics and epidemiology of Pseudomonas syringae: Challenges in Bacterial Molecular Plant Pathology

Localized recombination drives diversification of killing spectra for phage-derived syringacins

To better understand the potential for antagonistic interactions between members of the same bacterial species, we have surveyed bacteriocin killing activity across a diverse suite of strains of the phytopathogen Pseudomonas syringae. Our data demonstrate that killing activity from phage-derived bacteriocins of P. syringae (R-type syringacins) is widespread. Despite a high overall diversity of bacteriocin activity, strains can broadly be classified into five main killing types and two main sensitivity types. Furthermore, we show that killing activity switches frequently between strains and that switches correlate with localized recombination of two genes that together encode the proteins that specify bacteriocin targeting. Lastly, we demonstrate that phage-derived bacteriocin killing activity can be swapped between strains simply through expression of these two genes in trans. Overall, our study characterizes extensive diversity of killing activity for phage-derived bacteriocins of P. syringae across strains and highlights the power of localized recombination to alter phenotypes that mediate strain interactions during evolution of natural populations and communities.

Recombination of ecologically and evolutionarily significant loci maintains genetic cohesion in the Pseudomonas syringae species complex

Background

Pseudomonas syringae is a highly diverse bacterial species complex capable of causing a wide range of serious diseases on numerous agronomically important crops. We examine the evolutionary relationships of 391 agricultural and environmental strains using whole-genome sequencing and evolutionary genomic analyses.

Results

We describe the phylogenetic distribution of all 77,728 orthologous gene families in the pan-genome, reconstruct the core genome phylogeny using the 2410 core genes, hierarchically cluster the accessory genome, identify the diversity and distribution of type III secretion systems and their effectors, predict ecologically and evolutionary relevant loci, and establish the molecular evolutionary processes operating on gene families. Phylogenetic and recombination analyses reveals that the species complex is subdivided into primary and secondary phylogroups, with the former primarily comprised of agricultural isolates, including all of the well-studied P. syringae strains. In contrast, the secondary phylogroups include numerous environmental isolates. These phylogroups also have levels of genetic diversity typically found among distinct species. An analysis of rates of recombination within and between phylogroups revealed a higher rate of recombination within primary phylogroups than between primary and secondary phylogroups. We also find that “ecologically significant” virulence-associated loci and “evolutionarily significant” loci under positive selection are over-represented among loci that undergo inter-phylogroup genetic exchange.

Conclusions

While inter-phylogroup recombination occurs relatively rarely, it is an important force maintaining the genetic cohesion of the species complex, particularly among primary phylogroup strains. This level of genetic cohesion, and the shared plant-associated niche, argues for considering the primary phylogroups as a single biological species.

Electronic supplementary material

The online version of this article (10.1186/s13059-018-1606-y) contains supplementary material, which is available to authorized users.

Multiple quantitative trait loci contribute to resistance to bacterial canker incited by Pseudomonas syringae pv. actinidiae in kiwifruit (Actinidia chinensis)

Pseudomonas syringae pv. actinidiae (Psa) biovar 3, a virulent, canker-inducing pathogen is an economic threat to the kiwifruit (Actinidia spp.) industry worldwide. The commercially grown diploid (2×) A. chinensis var. chinensis is more susceptible to Psa than tetraploid and hexaploid kiwifruit. However information on the genetic loci modulating Psa resistance in kiwifruit is not available. Here we report mapping of quantitative trait loci (QTLs) regulating resistance to Psa in a diploid kiwifruit population, derived from a cross between an elite Psa-susceptible 'Hort16A' and a resistant male breeding parent P1. Using high-density genetic maps and intensive phenotyping, we identified a single QTL for Psa resistance on Linkage Group (LG) 27 of 'Hort16A' revealing 16-19% phenotypic variance and candidate alleles for susceptibility and resistance at this loci. In addition, six minor QTLs were identified in P1 on distinct LGs, exerting 4-9% variance. Resistance in the F1 population is improved by additive effects from 'Hort16A' and P1 QTLs providing evidence that divergent genetic pathways interact to combat the virulent Psa strain. Two different bioassays further identified new QTLs for tissue-specific responses to Psa. The genetic marker at LG27 QTL was further verified for association with Psa resistance in diploid Actinidia chinensis populations. Transcriptome analysis of Psa-resistant and susceptible genotypes in field revealed hallmarks of basal defense and provided candidate RNA-biomarkers for screening for Psa resistance in greenhouse conditions.

Pseudomonas syringae: enterprising epiphyte and stealthy parasite

Pseudomonas syringae is best known as a plant pathogenic bacterium that causes diseases in a multitude of hosts, and it has been used as a model organism to understand the biology of plant disease. Pathogenic and non-pathogenic isolates of P. syringae are also commonly found living as epiphytes and in the wider environment, including water sources such as rivers and precipitation. Ice-nucleating strains of P. syringae are associated with frost damage to crops. The genomes of numerous strains of P. syringae have been sequenced and molecular genetic studies have elucidated many aspects of this pathogen's interaction with its host plants.

Biological control of Pseudomonas syringae pv. aptata on sugar beet with Bacillus pumilus SS-10.7 and Bacillus amyloliquefaciens (SS-12.6 and SS-38.4) strains

AIM

Assessment of biological control of Pseudomonas syringae pv. aptata using crude lipopeptide extracts (CLEs) of two Bacillus amyloliquefaciens strains (SS-12.6 and SS-38.4) and one Bacillus pumilus strain (SS-10.7).

METHODS AND RESULTS

The minimum inhibitory concentration (MIC) of CLEs and their combinations against the pathogen and potential interaction between the extracts were determined in vitro. The most effective antibacterial activity was achieved with the CLE from B. amyloliquefaciens SS-12.6, with an MIC value of 0·63 mg ml^-1 . Interactions between CLE combinations were mostly indifferent. The biocontrol potential of CLEs, mixtures of CLEs, and cell culture of B. amyloliquefaciens SS-12.6 was tested on sugar beet plants inoculated with P. syringae pv. aptata P53. The best result in inhibiting the appearance of tissue necrosis (up to 92%) was achieved with B. amyloliquefaciens SS-12.6 cell culture.

CONCLUSION

This work demonstrated significant biocontrol potential of the CLE and cell culture of B. amyloliquefaciens SS-12.6 which successfully suppress leaf spot disease severity on sugar beet plants.

SIGNIFICANCE AND IMPACT OF THE STUDY

The findings of biocontrol of sugar beet emerging pathogen will contribute to growers in terms of alternative disease control management. This study represents first assessment of biological control of P. syringae pv. aptata.

Some strains that have converged to infect Prunus spp. trees are members of distinct Pseudomonas syringae genomospecies and ecotypes as revealed by in silico genomic comparison

A complementary taxonomic and population genetic study was performed to delineate genetically and ecologically distinct species within the Pseudomonas syringae complex by assessing 16 strains including pathovar strains that have converged to infect Prunus spp. trees, and two outgroups. Both average nucleotide identity and genome-to-genome distance comparison methods revealed the occurrence of distinct genomospecies, namely 1, 2, 3 and 8 (sensu Gardan et al.), with the latter two being closely related. Strains classified as P. s. pv. morsprunorum clustered into two distinct genomospecies, namely 2 and 8. Both the AdaptML and hierarchical Bayesian analysis of population structure methods highlighted the presence of three ecotypes, and the taxonomically related genomospecies 3 and 8 strains were members of the same ecotype. The distribution of pathogenic and virulence-associated genetic traits among Pseudomonas strains did not reveal any distinct type III secretion system effector or phytotoxin distribution pattern that characterized single genomospecies and strains that infect Prunus spp. The complete WHOP (Woody HOst and Pseudomonas spp.) genomic region and the entire β-ketoadipate gene cluster, including the catBCA operon, were found only in the members of genomospecies 2 and in the two P. s. pv. morsprunorum strains of genomospecies 8. A reduced gene flow between the three ecotypes suggested that point mutations played a larger role during the evolution of the strains than recombination. Our data support the idea that Prunus trees can be infected by different strains of distinct Pseudomonas genomospecies/ecotypes through diverse mechanisms of host colonization and infection. Such strains may represent particular lineages that emerged from environments other than that of the infected plant upon acquiring genetic traits that gave them the ability to cause plant diseases. The complementary assessment of bacterial strains using both taxonomic approaches and methods that reveal ecologically homogeneous populations has proven useful in confirming the cohesion of bacterial clusters.

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.

Comparative genomics of Spiraeoideae-infecting Erwinia amylovora strains provides novel insight to genetic diversity and identifies the genetic basis of a low-virulence strain

Erwinia amylovora is the causal agent of fire blight, one of the most devastating diseases of apple and pear. Erwinia amylovora is thought to have originated in North America and has now spread to at least 50 countries worldwide. An understanding of the diversity of the pathogen population and the transmission to different geographical regions is important for the future mitigation of this disease. In this research, we performed an expanded comparative genomic study of the Spiraeoideae-infecting (SI) E. amylovora population in North America and Europe. We discovered that, although still highly homogeneous, the genetic diversity of 30 E. amylovora genomes examined was about 30 times higher than previously determined. These isolates belong to four distinct clades, three of which display geographical clustering and one of which contains strains from various geographical locations ('Widely Prevalent' clade). Furthermore, we revealed that strains from the Widely Prevalent clade displayed a higher level of recombination with strains from a clade strictly from the eastern USA, which suggests that the Widely Prevalent clade probably originated from the eastern USA before it spread to other locations. Finally, we detected variations in virulence in the SI E. amylovora strains on immature pear, and identified the genetic basis of one of the low-virulence strains as being caused by a single nucleotide polymorphism in hfq, a gene encoding an important virulence regulator. Our results provide insights into the population structure, distribution and evolution of SI E. amylovora in North America and Europe.

Defining essential processes in plant pathogenesis with Pseudomonas syringae pv. tomato DC3000 disarmed polymutants and a subset of key type III effectors

Pseudomonas syringae pv. tomato DC3000 and its derivatives cause disease in tomato, Arabidopsis and Nicotiana benthamiana. The primary virulence factors include a repertoire of 29 effector proteins injected into plant cells by the type III secretion system and the phytotoxin coronatine. The complete repertoire of effector genes and key coronatine biosynthesis genes have been progressively deleted and minimally reassembled to reconstitute basic pathogenic ability in N. benthamiana, and in Arabidopsis plants that have mutations in target genes that mimic effector actions. This approach and molecular studies of effector activities and plant immune system targets have highlighted a small subset of effectors that contribute to essential processes in pathogenesis. Most notably, HopM1 and AvrE1 redundantly promote an aqueous apoplastic environment, and AvrPtoB and AvrPto redundantly block early immune responses, two conditions that are sufficient for substantial bacterial growth in planta. In addition, disarmed DC3000 polymutants have been used to identify the individual effectors responsible for specific activities of the complete repertoire and to more effectively study effector domains, effector interplay and effector actions on host targets. Such work has revealed that AvrPtoB suppresses cell death elicitation in N. benthamiana that is triggered by another effector in the DC3000 repertoire, highlighting an important aspect of effector interplay in native repertoires. Disarmed DC3000 polymutants support the natural delivery of test effectors and infection readouts that more accurately reveal effector functions in key pathogenesis processes, and enable the identification of effectors with similar activities from a broad range of other pathogens that also defeat plants with cytoplasmic effectors.

Pseudomonas syringae: what it takes to be a pathogen

Pseudomonas syringae is one of the best-studied plant pathogens and serves as a model for understanding host-microorganism interactions, bacterial virulence mechanisms and host adaptation of pathogens as well as microbial evolution, ecology and epidemiology. Comparative genomic studies have identified key genomic features that contribute to P. syringae virulence. P. syringae has evolved two main virulence strategies: suppression of host immunity and creation of an aqueous apoplast to form its niche in the phyllosphere. In addition, external environmental conditions such as humidity profoundly influence infection. P. syringae may serve as an excellent model to understand virulence and also of how pathogenic microorganisms integrate environmental conditions and plant microbiota to become ecologically robust and diverse pathogens of the plant kingdom.

Entomopathogenicity to Two Hemipteran Insects Is Common but Variable across Epiphytic Pseudomonas syringae Strains

Strains of the well-studied plant pathogen Pseudomonas syringae show large differences in their ability to colonize plants epiphytically and to inflict damage to hosts. Additionally, P. syringae can infect some sap-sucking insects and at least one P. syringae strain is highly virulent to insects, causing death to most individuals within as few as 4 days and growing to high population densities within insect hosts. The likelihood of agricultural pest insects coming into contact with transient populations of P. syringae while feeding on plants is high, yet the ecological implications of these interactions are currently not well understood as virulence has not been tested across a wide range of strains. To investigate virulence differences across strains we exposed the sweet potato whitefly, Bemisia tabaci, and the pea aphid, Acyrthosiphon pisum, both of which are cosmopolitan agricultural pests, to 12 P. syringae strains. We used oral inoculations with bacteria suspended in artificial diet in order to assay virulence while controlling for other variables such as differences in epiphytic growth ability. Generally, patterns of pathogenicity remain consistent across the two species of hemipteran insects, with bacterial strains from phylogroup II, or genomospecies 1, causing the highest rate of mortality with up to 86% of individuals dead after 72 h post infection. The rate of mortality is highly variable across strains, some significantly different from negative control treatments and others showing no discernable difference. Interestingly, one of the most pathogenic strains to both aphids and whiteflies (Cit7) is thought to be non-pathogenic on plants. We also found Cit7 to establish the highest epiphytic population after 48 h on fava beans. Between the nine P. syringae strains tested for epiphytic ability there is also much variation, but epiphytic ability was positively correlated with pathogenicity to insects, suggesting that the two traits may be linked and that strains likely to be found on plants may often be entomopathogenic. Our study highlights that there may be a use for epiphytic bacteria in the biological control of insect crop pests. It also suggests that interactions with epiphytic bacteria could be evolutionary and ecological drivers for hemipteran insects.

Clarification of Taxonomic Status within the Pseudomonas syringae Species Group Based on a Phylogenomic Analysis

The Pseudomonas syringae phylogenetic group comprises 15 recognized bacterial species and more than 60 pathovars. The classification and identification of strains is relevant for practical reasons but also for understanding the epidemiology and ecology of this group of plant pathogenic bacteria. Genome-based taxonomic analyses have been introduced recently to clarify the taxonomy of the whole genus. A set of 139 draft and complete genome sequences of strains belonging to all species of the P. syringae group available in public databases were analyzed, together with the genomes of closely related species used as outgroups. Comparative genomics based on the genome sequences of the species type strains in the group allowed the delineation of phylogenomic species and demonstrated that a high proportion of strains included in the study are misclassified. Furthermore, representatives of at least 7 putative novel species were detected. It was also confirmed that P. ficuserectae, P. meliae, and P. savastanoi are later synonyms of P. amygdali and that “P. coronafaciens” should be revived as a nomenspecies.

Recombination of Virulence Genes in Divergent Acidovorax avenae Strains That Infect a Common Host

Bacterial etiolation and decline (BED), caused by Acidovorax avenae, is an emerging disease of creeping bentgrass on golf courses in the United States. We performed the first comprehensive analysis of A. avenae on a nationwide collection of turfgrass- and maize-pathogenic A. avenae. Surprisingly, our results reveal that the turfgrass-pathogenic A. avenae in North America are not only highly divergent but also belong to two distinct phylogroups. Both phylogroups specifically infect turfgrass but are more closely related to maize pathogens than to each other. This suggests that, although the disease is only recently reported, it has likely been infecting turfgrass for a long time. To identify a genetic basis for the host specificity, we searched for genes closely related among turfgrass strains but distantly related to their homologs from maize strains. We found a cluster of 11 such genes generated by three ancient recombination events within the type III secretion system (T3SS) pathogenicity island. Ever since the recombination, the cluster has been conserved by strong purifying selection, hinting at its selective importance. Together our analyses suggest that BED is an ancient disease that may owe its host specificity to a highly conserved cluster of 11 T3SS genes.

Similar levels of gene content variation observed for Pseudomonas syringae populations extracted from single and multiple host species

Bacterial strains of the same species collected from different hosts frequently exhibit differences in gene content. In the ubiquitous plant pathogen Pseudomonas syringae, more than 30% of genes encoded by each strain are not conserved among strains colonizing other host species. Although they are often implicated in host specificity, the role of this large fraction of the genome in host-specific adaptation is largely unexplored. Here, we sought to relate variation in gene content between strains infecting different species to variation that persists between strains on the same host. We fully sequenced a collection of P. syringae strains collected from wild Arabidopsis thaliana populations in the Midwestern United States. We then compared patterns of variation observed in gene content within these A. thaliana-isolated strains to previously published P. syringae sequence from strains collected on a diversity of crop species. We find that strains collected from the same host, A. thaliana, differ in gene content by 21%, 2/3 the level of gene content variation observed across strains collected from different hosts. Furthermore, the frequency with which specific genes are present among strains collected within the same host and among strains collected from different hosts is highly correlated. This implies that most gene content variation is maintained irrespective of host association. At the same time, we identify specific genes whose presence is important for P. syringae's ability to flourish within A. thaliana. Specifically, the A. thaliana strains uniquely share a genomic island encoding toxins active against plants and surrounding microbes, suggesting a role for microbe-microbe interactions in dictating the abundance within this host. Overall, our results demonstrate that while variation in the presence of specific genes can affect the success of a pathogen within its host, the majority of gene content variation is not strongly associated with patterns of host use.

The Scientific, Economic, and Social Impacts of the New Zealand Outbreak of Bacterial Canker of Kiwifruit (Pseudomonas syringae pv. actinidiae)

The introduction of Pseudomonas syringae pv. actinidiae (Psa) severely damaged the New Zealand kiwifruit industry, which in 2010 was based on only two cultivars. Despite an extraordinarily quick and strong response by industry, government, and scientists to minimize the economic and social impacts, the economic consequences of this outbreak were severe. Although our understanding of Psa epidemiology and control methods increased substantively over the past six years, the kiwifruit industry largely recovered because of the introduction of a less-susceptible yellow-fleshed cultivar. The New Zealand population of Psa is clonal but has evolved rapidly since its introduction by exchanging mobile genetic elements, including integrative conjugative elements (ICEs), with the local bacterial populations. In some cases, this has led to copper resistance. It is currently believed that the center of origin of the pathogen is Japan or Korea, but biovar 3, which is responsible for the global outbreak, originated in China.

Adaptation, specialization, and coevolution within phytobiomes

Growth patterns of individual plants and evolutionary trajectories of plant communities are intimately linked with and are critically affected by host-associated microbiomes. Research across systems has begun to shed light on how these phytobiomes are established under laboratory and natural conditions, and have cultivated hope that a better understanding of the governing principles for host-microbe interactions can guide attempts to engineer microbiomes to boost agricultural yields. One important, yet relatively understudied, parameter in regards to phytobiome membership is the degree to which specialization and coevolution between plant species and microbes provides structure to these communities. In this article, I provide an overview of mechanisms enabling adaptation and specialization of phytobiome communities to host plants as well as the potential for plants themselves to recruit and cultivate beneficial interactions. I further explore the possibility of host-beneficial microbe coevolution and suggest particular situations that could promote the evolution of such close-knit partnerships. It is my hope that this overview will encourage future experiments that can begin to fill in this black box of ecological and evolutionary interactions across phytobiomes.

Knots Untie: Molecular Determinants Involved in Knot Formation Induced by Pseudomonas savastanoi in Woody Hosts

The study of the molecular basis of tree diseases is lately receiving a renewed attention, especially with the emerging perception that pathogens require specific pathogenicity and virulence factors to successfully colonize woody hosts. Pathosystems involving woody plants are notoriously difficult to study, although the use of model bacterial strains together with genetically homogeneous micropropagated plant material is providing a significant impetus to our understanding of the molecular determinants leading to disease. The gammaproteobacterium Pseudomonas savastanoi belongs to the intensively studied Pseudomonas syringae complex, and includes three pathogenic lineages causing tumorous overgrowths (knots) in diverse economically relevant trees and shrubs. As it occurs with many other bacteria, pathogenicity of P. savastanoi is dependent on a type III secretion system, which is accompanied by a core set of at least 20 effector genes shared among strains isolated from olive, oleander, and ash. The induction of knots of wild-type size requires that the pathogen maintains adequate levels of diverse metabolites, including the phytohormones indole-3-acetic acid and cytokinins, as well as cyclic-di-GMP, some of which can also regulate the expression of other pathogenicity and virulence genes and participate in bacterial competitiveness. In a remarkable example of social networking, quorum sensing molecules allow for the communication among P. savastanoi and other members of the knot microbiome, while at the same time are essential for tumor formation. Additionally, a distinguishing feature of bacteria from the P. syringae complex isolated from woody organs is the possession of a 15 kb genomic island (WHOP) carrying four operons and three other genes involved in degradation of phenolic compounds. Two of these operons mediate the catabolism of anthranilate and catechol and, together with another operon, are required for the induction of full-size tumors in woody hosts, but not in non-woody micropropagated plants. The use of transposon mutagenesis also uncovered a treasure trove of additional P. savastanoi genes affecting virulence and participating in diverse bacterial processes. Although there is still much to be learned on what makes a bacterium a successful pathogen of trees, we are already untying the knots.

Suppression of Plant Immune Responses by the Pseudomonas savastanoi pv. savastanoi NCPPB 3335 Type III Effector Tyrosine Phosphatases HopAO1 and HopAO2

The effector repertoire of the olive pathogen P. savastanoi pv. savastanoi NCPPB 3335 includes two members of the HopAO effector family, one of the most diverse T3E families of the P. syringae complex. The study described here explores the phylogeny of these dissimilar members, HopAO1 and HopAO2, among the complex and reveals their activities as immune defense suppressors. Although HopAO1 is predominantly encoded by phylogroup 3 strains isolated from woody organs of woody hosts, both HopAO1 and HopAO2 are phylogenetically clustered according to the woody/herbaceous nature of their host of isolation, suggesting host specialization of the HopAO family across the P. syringae complex. HopAO1 and HopAO2 translocate into plant cells and show hrpL-dependent expression, which allows their classification as actively deployed type III effectors. Our data also show that HopAO1 and HopAO2 possess phosphatase activity, a hallmark of the members of this family. Both of them exert an inhibitory effect on early plant defense responses, such as ROS production and callose deposition, and are able to suppress ETI responses induced by the effectorless polymutant of P. syringae pv. tomato DC3000 (DC3000D28E) in Nicotiana. Moreover, we demonstrate that a ΔhopAO1 mutant of P. savastanoi NCPBB 3335 exhibits a reduced fitness and virulence in olive plants, which supports the relevance of this effector during the interaction of this strain with its host plants. This work contributes to the field with the first report regarding functional analysis of HopAO homologs encoded by P. syringae or P. savastanoi strains isolated from woody hosts.