Exposure to host resistance mechanisms drives evolution of bacterial virulence in plants

Bacterial pathogenicity to plants and animals has evolved through an arms race of attack and defense. Key players are bacterial effector proteins, which are delivered through the type III secretion system and suppress basal defenses . In plants, varietal resistance to disease is based on recognition of effectors by the products of resistance (R) genes . When recognized, the effector or in this scenario, avirulence (Avr) protein triggers the hypersensitive resistance reaction (HR), which generates antimicrobial conditions . Unfortunately, such gene-for-gene-based resistance commonly fails because of the emergence of virulent strains of the pathogen that no longer trigger the HR . We have followed the emergence of a new virulent pathotype of the halo-blight pathogen Pseudomonas syringae pv. phaseolicola within leaves of a resistant variety of bean. Exposure to the HR led to the selection of strains lacking the avirulence (effector) gene avrPphB (or hopAR1), which triggers defense in varieties with the matching R3 resistance gene. Loss of avrPphB was through deletion of a 106 kb genomic island (PPHGI-1) that shares features with integrative and conjugative elements (ICElands) and also pathogenicity islands (PAIs) in diverse bacteria . We provide a molecular explanation of how exposure to resistance mechanisms in plants drives the evolution of new virulent forms of pathogens.     

Excision from tRNA genes of a large chromosomal region, carrying avrPphB, associated with race change in the bean pathogen, Pseudomonas syringae pv. phaseolicola

Pseudomonas syringae pv. phaseolicola (Pph) race 4 strain 1302A carries avirulence gene avrPphB. Strain RJ3, a sectoral variant from a 1302A culture, exhibited an extended host range in cultivars of bean and soybean resulting from the absence of avrPphB from the RJ3 chromosome. Complementation of RJ3 with avrPphB restored the race 4 phenotype. Both strains showed similar in planta growth in susceptible bean cultivars. Analysis of RJ3 indicated loss of > 40 kb of DNA surrounding avrPphB. Collinearity of the two genomes was determined for the left and right junctions of the deleted avrPphB region; the left junction is approximately 19 kb and the right junction > 20 kb from avrPphB in 1302A. Sequencing revealed that the region containing avrPphB was inserted into a tRNALYS gene, which was re-formed at the right junction in strain 1302A. A putative lysine tRNA pseudogene (PsitRNALYS) was found at the left junction of the insertion. All tRNA genes were in identical orientation in the chromosome. Genes near the left junction exhibited predicted protein homologies with gene products associated with a virulence locus of the periodontal pathogen Actinobacillus actinomycetemcomitans. Specific oligonucleotide primers that differentiate 1302A from RJ3 were designed and used to demonstrate that avrPphB was located in different regions of the chromosome in other strains of Pph. Deletion of a large region of the chromosome containing an avirulence gene represents a new route to race change in Pph.     

Molecular determinants required for the avirulence function of AvrPphB in bean and other plants

The avirulence gene avrPphB from Pseudomonas syringae pv. phaseolicola determines incompatibility, manifested as a hypersensitive reaction (HR), on bean cultivars carrying the R3 resistance gene and also confers avirulence on other plants. The AvrPphB protein carries an embedded consensus myristoylation motif and is cleaved in bacteria and certain plants to yield fragments of about 6 and 28 kDa. We investigated plant recognition and type III translocation determinants in AvrPphB by constructing three N-terminally truncated and two site-directed mutants carrying substitutions in the conserved G63 residue of the myristoylation motif, which lies adjacent to the proteolytic cleavage site. The peptides were either delivered to plant cells by pseudomonads or were expressed transiently in planta via the Agrobacterium tumefaciens or Potato virus X. The 63 amino terminal residues were required for type III-mediated translocation from Pseudomonas strains to the plant, but were partially dispensable for effector recognition following in planta expression. Substitution of the G63 residue resulted in differential HR phenotypes in two different R3 cultivars of bean and abolished effector processing in Pseudomonas strains. Agrobacterium-mediated expression of the mutant proteins elicited HR in resistant bean hosts and in tomato but elicited no reaction in Nicotiana species.     

Sequence variations in alleles of the avirulence gene avrPphE.R2 from Pseudomonas syringae pv. phaseolicola lead to loss of recognition of the AvrPphE protein within bean cells and a gain in cultivar-specific virulence

The bean halo blight pathogen, Pseudomonas syringae pv. phaseolicola ( Psph ), is differentiated into nine races based on the presence or absence of five avirulence ( avr ) genes in the bacterium, which interact with corresponding resistance genes, R1–R5 , in Phaseolus vulgaris . The resistance gene R2 is matched by avrPphE , which is located adjacent to the cluster of hrp genes that are required for pathogenicity of Psph . Although only races 2, 4, 5 and 7 are avirulent on cultivars with R2 (inducing the hypersensitive response; HR), homologues of avrPphE are present in all races of Psph . DNA sequencing of avrPphE alleles from races of Psph has demonstrated two routes to virulence: via single basepair changes conferring amino acid substitutions in races 1, 3, 6 and 9 and an insertion of 104 bp in the allele in race 8. We have demonstrated that these base changes are responsible for the difference between virulence and avirulence by generating transconjugants of a virulent race harbouring plasmids expressing the various alleles of avrPphE. Agrobacterium tumefaciens -directed expression of avrPphE from race 4 in bean leaves induced the HR in a resistance gene-specific manner, suggesting that the AvrPphE protein is alone required for HR induction and is recognized within the plant cell. The allele from race 6, which is inactive if expressed in Psph , elicited a weak HR if expressed in planta , whereas the allele from race 1 did not. Our results suggest that the affinity of interaction between AvrPphE homologues and an unknown plant receptor mediates the severity of the plant's response. Mutation of avrPphE alleles did not affect the ability to colonize bean from a low level of inoculum. The avirulence gene avrPphB , which matches the R3 resistance gene, also caused a gene-specific HR following expression in the plant after delivery by A. tumefaciens .     

Identification of three putative signal transduction genes involved in R gene-specified disease resistance in Arabidopsis

The RPS5 disease resistance gene of Arabidopsis mediates recognition of Pseudomonas syringae strains that possess the avirulence gene avrPphB. By screening for loss of RPS5-specified resistance, we identified five pbs (avrPphB susceptible) mutants that represent three different genes. Mutations in PBS1 completely blocked RPS5-mediated resistance, but had little to no effect on resistance specified by other disease resistance genes, suggesting that PBS1 facilitates recognition of the avrPphB protein. The pbs2 mutation dramatically reduced resistance mediated by the RPS5 and RPM1 resistance genes, but had no detectable effect on resistance mediated by RPS4 and had an intermediate effect on RPS2-mediated resistance. The pbs2 mutation also had varying effects on resistance mediated by seven different RPP (recognition of Peronospora parasitica) genes. These data indicate that the PBS2 protein functions in a pathway that is important only to a subset of disease-resistance genes. The pbs3 mutation partially suppressed all four P. syringae-resistance genes (RPS5, RPM1, RPS2, and RPS4), and it had weak-to-intermediate effects on the RPP genes. In addition, the pbs3 mutant allowed higher bacterial growth in response to a virulent strain of P. syringae, indicating that the PBS3 gene product functions in a pathway involved in restricting the spread of both virulent and avirulent pathogens. The pbs mutations are recessive and have been mapped to chromosomes I (pbs2) and V (pbs1 and pbs3).     

Layered basal defenses underlie non-host resistance of Arabidopsis to Pseudomonas syringae pv. phaseolicola

Arabidopsis is a non-host for Pseudomonas syringae pv. phaseolicola NPS3121 (Pph), a bacterial pathogen of bean. Pph does not induce a hypersensitive response in Arabidopsis. Here we show that Arabidopsis instead resists Pph with multi-layered basal defense. Our approach was: (i) to identify defense readouts induced by Pph; (ii) to determine whether mutations in known Arabidopsis defense genes disrupt Pph-induced defense signaling; (iii) to determine whether heterologous type III effectors from pathogens of Arabidopsis suppress Pph-induced defense signaling, and (iv) to ascertain how basal defenses contribute to resistance against Pph by individually or multiply disrupting defense signaling pathways with mutations and heterologous type III effectors. We demonstrate that Pph elicits a minimum of three basal defense-signaling pathways in Arabidopsis. These pathways have unique readouts, including PR-1 protein accumulation and morphologically distinct types of callose deposition. Further, they require distinct defense genes, including PMR4, RAR1, SID2, NPR1, and PAD4 . Finally, they are suppressed differentially by heterologous type III effectors, including AvrRpm1 and HopM1. Pph growth is enhanced only when multiple defense pathways are disrupted. For example, mutation of NPR1 or SID2 combined with the action of AvrRpm1 and HopM1 renders Arabidopsis highly susceptible to Pph. Thus, non-host resistance of Arabidopsis to Pph is based on multiple, individually effective layers of basal defense.     

Transgenic tobacco plants expressing the bacterial mc gene resist virus infection

A bacterial rnc gene coding for a double-stranded RNA-dependent RNase III endoribonuclease and a mutant, rnc70, were expressed in tobacco plants. The RNase III protein produced in the transgenic plants was the same size as the bacterial protein. Expression of the wild-type gene could cause stunting in some plant lines, but not in others. Expression of the mutant protein did not affect normal growth and development of the transgenic plants. Transgenic plants of the R1 and R2 generations, expressing the wild type, as well as a mutant protein, were resistant to infection by three disparate RNA plant viruses with a divided genome but not against two viruses with a single-stranded RNA genome. Introduction of the rnc gene in crop plants may provide resistance to economically important virus diseases.     

Crystal structures of the type III effector protein AvrPphF and its chaperone reveal residues required for plant pathogenesis

The avrPphF locus from Pseudomonas syringae pv. phaseolicola, the causative agent of bean halo-blight disease, encodes proteins which either enhance virulence on susceptible hosts or elicit defense responses on hosts carrying the R1 resistance gene. Here we present the crystal structures of the two proteins from the avrPphF operon. The structure of AvrPphF ORF1 is strikingly reminiscent of type III chaperones from bacterial pathogens of animals, indicating structural conservation of these specialized chaperones, despite high sequence divergence. The AvrPphF ORF2 effector adopts a novel "mushroom"-like structure containing "head" and "stalk" subdomains. The head subdomain possesses limited structural homology to the catalytic domain of bacterial ADP-ribosyltransferases (ADP-RTs), though no ADP-RT activity was detected for AvrPphF ORF2 in standard assays. Nonetheless, this structural similarity identified two clusters of conserved surface-exposed residues important for both virulence mediated by AvrPphF ORF2 and recognition of this effector by bean plants expressing the R1 resistance gene.     

Arabidopsis GH3-LIKE DEFENSE GENE 1 is required for accumulation of salicylic acid, activation of defense responses and resistance to Pseudomonas syringae

In Arabidopsis, the GH3-like gene family consists of 19 members, several of which have been shown to adenylate the plant hormones jasmonic acid, indole acetic acid and salicylic acid (SA). In some cases, this adenylation has been shown to catalyze hormone conjugation to amino acids. Here we report molecular characterization of the GH3-LIKE DEFENSE GENE 1 (GDG1), a member of the GH3-like gene family, and show that GDG1 is an important component of SA-mediated defense against the bacterial pathogen Pseudomonas syringae. Expression of GDG1 is induced earlier and to a higher level in response to avirulent pathogens compared to virulent pathogens. gdg1 null mutants are compromised in several pathogen defense responses, including activation of defense genes and resistance against virulent and avirulent bacterial pathogens. Accumulation of free and glucoside-conjugated SA (SAG) in response to pathogen infection is compromised in gdg1 mutants. All defense-related phenotypes of gdg1 can be rescued by external application of SA, suggesting that gdg1 mutants are defective in the SA-mediated defense pathway(s) and that GDG1 functions upstream of SA. Our results suggest that GDG1 contributes to both basal and resistance gene-mediated inducible defenses against P. syringae (and possibly other pathogens) by playing a critical role in regulating the levels of pathogen-inducible SA. GDG1 is allelic to the PBS3 (avrPphB susceptible) gene.     

EDS1 in tomato is required for resistance mediated by TIR-class R genes and the receptor-like R gene Ve

In tobacco and other Solanaceae species, the tobacco N gene confers resistance to tobacco mosaic virus (TMV), and leads to induction of standard defense and resistance responses. Here, we report the use of N-transgenic tomato to identify a fast-neutron mutant, sun1-1 (suppressor of N), that is defective in N-mediated resistance. Induction of salicylic acid (SA) and expression of pathogenesis-related (PR) genes, each signatures of systemic acquired resistance, are both dramatically suppressed in sun1-1 plants after TMV treatment compared to wild-type plants. Application of exogenous SA restores PR gene expression, indicating that SUN1 acts upstream of SA. Upon challenge with additional pathogens, we found that the sun1-1 mutation impairs resistance mediated by certain resistance (R) genes, (Bs4, I, and Ve), but not others (Mi-1). In addition, sun1-1 plants exhibit enhanced susceptibility to TMV, as well as to virulent pathogens. sun1-1 has been identified as an EDS1 homolog present on chromosome 6 of tomato. The discovery of enhanced susceptibility in the sun1-1 (Le_eds1-1) mutant plant, which contrasts to reports in Nicotiana benthamiana using virus-induced gene silencing, provides evidence that the intersection of R gene-mediated pathways with general resistance pathways is conserved in a Solanaceous species. In tomato, EDS1 is important for mediating resistance to a broad range of pathogens (viral, bacterial, and fungal pathogens), yet shows specificity in the class of R genes that it affects (TIR-NBS-LRR as opposed to CC-NBS-LRR). In addition, a requirement for EDS1 for Ve-mediated resistance in tomato exposes that the receptor-like R gene class may also require EDS1.     

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.     

Molecular evolution of virulence in natural field strains of Xanthomonas campestris pv. vesicatoria

The avrBs2 avirulence gene of the bacterial plant pathogen Xanthomonas campestris pv. vesicatoria triggers disease resistance in pepper plants containing the Bs2 resistance gene and contributes to bacterial virulence on susceptible host plants. We studied the effects of the pepper Bs2 gene on the evolution of avrBs2 by characterizing the molecular basis for virulence of 20 X. campestris pv. vesicatoria field strains that were isolated from disease spots on previously resistant Bs2 pepper plants. All field strains tested were complemented by a wild-type copy of avrBs2 in their ability to trigger disease resistance on Bs2 plants. DNA sequencing revealed four mutant alleles of avrBs2, two of which consisted of insertions or deletions of 5 nucleotides in a repetitive region of avrBs2. The other two avrBs2 alleles were characterized by point mutations with resulting single amino acid changes (R403P or A410D). We generated isogenic X. campestris pv. vesicatoria strains by chromosomal avrBs2 gene exchange to study the effects of these mutations on the dual functions of avrBs2 in enhancing bacterial virulence and inducing plant resistance by in planta bacterial growth experiments. The deletion of 5 nucleotides led to loss of avrBs2-induced resistance on Bs2 pepper plants and abolition of avrBs2-mediated enhancement of fitness on susceptible plants. Significantly, the point mutations led to minimal reduction in virulence function of avrBs2 on susceptible pepper plants, with either minimal (R403P allele) or an intermediate level of (A410D allele) triggering of resistance on Bs2 plants. Consistent with the divergent selection pressures on avrBs2 exerted by the Bs2 resistance gene, our results show that avrBs2 is evolving to decrease detection by the Bs2 gene while at the same time maintaining its virulence function.     

Distribution and sequence analysis of a family of type ill-dependent effectors correlate with the phylogeny of Ralstonia solanacearum strains

In Ralstonia solanacearum, we previously have reported on the characterization of popP1 and popP2 genes. These genes encode type III-dependent pathogenicity effectors related to the large family of AvrRxv/YopJ cysteine proteases that are shared among pathogens of plants and animals. In this study, we identify a third gene, named popP3, that is inactivated in the genome sequence of strain GMI1000 by insertion of a copy of the insertion sequence ISRso13. The three popP genes are localized on two large chromosomal pathogenicity islands, with popP1 and popP2 being present on the same island. Phylogenic analysis demonstrated that the PopP2 and PopP3 proteins are clearly distinct from other effectors of this family previously characterized in plant and animal pathogens. Analysis of the distribution and allelic variations of the three genes in 30 strains representative of the biodiversity of R. solanacearum established that popP genes are distributed widely among strains from two of the three phyla previously defined on the basis of the structure of the core genome. Sequencing of the popP genes from the different strains revealed limited allelic variations at the three loci but did not show evidence of recombination between the popP genes. Limited allelic variation together with occurrence of insertion sequences within or in the close vicinity of popP genes and the presence of gene duplications in these pathogenicity islands suggest that genomic rearrangements might be a major evolutionary driving force controlling evolution of the genes encoded in these regions. The implications of these observations in terms of bacterial evolution, gene acquisition, and horizontal gene transfers are discussed.     

Genomic and Gene-Expression Comparisons among Phage-Resistant Type-IV Pilus Mutants of Pseudomonas syringae pathovar phaseolicola

Pseudomonas syringae pv. phaseolicola (Pph) is a significant bacterial pathogen of agricultural crops, and phage Φ6 and other members of the dsRNA virus family Cystoviridae undergo lytic (virulent) infection of Pph, using the type IV pilus as the initial site of cellular attachment. Despite the popularity of Pph/phage Φ6 as a model system in evolutionary biology, Pph resistance to phage Φ6 remains poorly characterized. To investigate differences between phage Φ6 resistant Pph strains, we examined genomic and gene expression variation among three bacterial genotypes that differ in the number of type IV pili expressed per cell: ordinary (wild-type), non-piliated, and super-piliated. Genome sequencing of non-piliated and super-piliated Pph identified few mutations that separate these genotypes from wild type Pph–and none present in genes known to be directly involved in type IV pilus expression. Expression analysis revealed that 81.1% of gene ontology (GO) terms up-regulated in the non-piliated strain were down-regulated in the super-piliated strain. This differential expression is particularly prevalent in genes associated with respiration—specifically genes in the tricarboxylic acid cycle (TCA) cycle, aerobic respiration, and acetyl-CoA metabolism. The expression patterns of the TCA pathway appear to be generally up and down-regulated, in non-piliated and super-piliated Pph respectively. As pilus retraction is mediated by an ATP motor, loss of retraction ability might lead to a lower energy draw on the bacterial cell, leading to a different energy balance than wild type. The lower metabolic rate of the super-piliated strain is potentially a result of its loss of ability to retract.     

A mutation within the leucine-rich repeat domain of the Arabidopsis disease resistance gene RPS5 partially suppresses multiple bacterial and downy mildew resistance genes.

Recognition of pathogens by plants is mediated by several distinct families of functionally variable but structurally related disease resistance (R) genes. The largest family is defined by the presence of a putative nucleotide binding domain and 12 to 21 leucine-rich repeats (LRRs). The function of these LRRs has not been defined, but they are speculated to bind pathogen-derived ligands. We have isolated a mutation in the Arabidopsis RPS5 gene that indicates that the LRR region may interact with other plant proteins. The rps5-1 mutation causes a glutamate-to-lysine substitution in the third LRR and partially compromises the function of several R genes that confer bacterial and downy mildew resistance. The third LRR is relatively well conserved, and we speculate that it may interact with a signal transduction component shared by multiple R gene pathways.     

A Small GTPase of the Rab Family Is Required for Root Hair Formation and Preinfection Stages of the Common Bean–Rhizobium Symbiotic Association

Legume plants are able to establish a symbiotic relationship with soil bacteria from the genus Rhizobium, leading to the formation of nitrogen-fixing root nodules. Successful nodulation requires both the formation of infection threads (ITs) in the root epidermis and the activation of cell division in the cortex to form the nodule primordium. This study describes the characterization of RabA2, a common bean (Phaseolus vulgaris) cDNA previously isolated as differentially expressed in root hairs infected with Rhizobium etli, which encodes a protein highly similar to small GTPases of the RabA2 subfamily. This gene is expressed in roots, particularly in root hairs, where the protein was found to be associated with vesicles that move along the cell. The role of this gene during nodulation has been studied in common bean transgenic roots using a reverse genetic approach. Examination of root morphology in RabA2 RNA interference (RNAi) plants revealed that the number and length of the root hairs were severely reduced in these plants. Upon inoculation with R. etli, nodulation was completely impaired and no induction of early nodulation genes (ENODs), such as ERN1, ENOD40, and Hap5, was detected in silenced hairy roots. Moreover, RabA2 RNAi plants failed to induce root hair deformation and to initiate ITs, indicating that morphological changes that precede bacterial infection are compromised in these plants. We propose that RabA2 acts in polar growth of root hairs and is required for reorientation of the root hair growth axis during bacterial infection.     

Constitutive disease resistance requires EDS1 in the Arabidopsis mutants cpr1 and cpr6 and is partially EDS1-dependent in cpr5

The systemic acquired resistance (SAR) response in Arabidopsis is characterized by the accumulation of salicylic acid (SA), expression of the pathogenesis-related (PR) genes, and enhanced resistance to virulent bacterial and oomycete pathogens. The cpr (constitutive expressor of PR genes) mutants express all three SAR phenotypes. In addition, cpr5 and cpr6 induce expression of PDF1.2, a defense-related gene associated with activation of the jasmonate/ethylene-mediated resistance pathways. cpr5 also forms spontaneous lesions. In contrast, the eds1 (enhanced disease susceptibility) mutation abolishes race-specific resistance conferred by a major subclass of resistance (R) gene products in response to avirulent pathogens. eds1 plants also exhibit increased susceptibility to virulent pathogens. Epistasis experiments were designed to explore the relationship between the cpr- and EDS1-mediated resistance pathways. We found that a null eds1 mutation suppresses the disease resistance phenotypes of both cpr1 and cpr6. In contrast, eds1 only partially suppresses resistance in cpr5, leading us to conclude that cpr5 expresses both EDS1-dependent and EDS1-independent components of plant disease resistance. Although eds1 does not prevent lesion formation on cpr5 leaves, it alters their appearance and reduces their spread. This phenotypic difference is associated with increased pathogen colonization of cpr5 eds1 plants compared to cpr5. The data allow us to place EDS1 as a necessary downstream component of cpr1- and cpr6-mediated responses, but suggest a more complex relationship between EDS1 and cpr5 in plant defense.     

Tomato MAPKKKε is a positive regulator of cell-death signaling networks associated with plant immunity

Mitogen-activated protein (MAP) kinase cascades are fundamental components of the signaling pathways associated with plant immunity. Despite the large number of MAP kinase kinase kinases (MAPKKK) encoded in the plant genome, only very few of them have an assigned function. Here, we identified MAPKKK gene of tomato (Solanum lycopersicum), SIMAPKKKε, which is required for hypersensitive response cell death and disease resistance against Gram-negative bacterial pathogens. Silencing of SIMAPKKKε compromised tomato resistance to Xanthomonas campestris and Pseudomonas syringae strains, resulting in the appearance of disease symptoms and enhanced bacterial growth. In addition, silencing of NbMAPKKKε in Nicotiana benthamiana plants significantly inhibited the cell death triggered by expression of different R gene/effector gene pairs. Conversely, overexpression of either the full-length SIMAPKKKε gene or its kinase domain in N. benthamiana leaves caused pathogen-independent activation of cell death that required an intact kinase catalytic domain. Moreover, by suppressing the expression of various MAPKK and MAPK genes and overexpressing the SIMAPKKKε kinase domain, we identified a signaling cascade acting downstream of SIMAPKKKε that includes MEK2, WIPK and SIPK. Additional epistasis experiments revealed that SIPKK functions as a negative regulator of SIMAPKKKε-mediated cell death. Our results provide evidence that SIMAPKKKε is a signaling molecule that positively regulates cell death networks associated with plant immunity.