Differential recognition of highly divergent downy mildew avirulence gene alleles by RPP1 resistance genes from two Arabidopsis lines

The perception of downy mildew avirulence (Arabidopsis thaliana Recognized [ATR]) gene products by matching Arabidopsis thaliana resistance (Recognition of Peronospora parasitica [RPP]) gene products triggers localized cell death (a hypersensitive response) in the host plant, and this inhibits pathogen development. The oomycete pathogen, therefore, is under selection pressure to alter the form of these gene products to prevent detection. That the pathogen maintains these genes indicates that they play a positive role in pathogen survival. Despite significant progress in cloning plant RPP genes and characterizing essential plant components of resistance signaling pathways, little progress has been made in identifying the oomycete molecules that trigger them. Concluding a map-based cloning effort, we have identified an avirulence gene, ATR1NdWsB, that is detected by RPP1 from the Arabidopsis accession Niederzenz in the cytoplasm of host plant cells. We report the cloning of six highly divergent alleles of ATR1NdWsB from eight downy mildew isolates and demonstrate that the ATR1NdWsB alleles are differentially recognized by RPP1 genes from two Arabidopsis accessions (Niederzenz and Wassilewskija). RPP1-Nd recognizes a single allele of ATR1NdWsB; RPP1-WsB also detects this allele plus three additional alleles with divergent sequences. The Emco5 isolate expresses an allele of ATR1NdWsB that is recognized by RPP1-WsB, but the isolate evades detection in planta. Although the Cala2 isolate is recognized by RPP1-WsA, the ATR1NdWsB allele from Cala2 is not, demonstrating that RPP1-WsA detects a novel ATR gene product. Cloning of ATR1NdWsB has highlighted the presence of a highly conserved novel amino acid motif in avirulence proteins from three different oomycetes. The presence of the motif in additional secreted proteins from plant pathogenic oomycetes and its similarity to a host-targeting signal from malaria parasites suggest a conserved role in pathogenicity.     

Arabidopsis RPP4 is a member of the RPP5 multigene family of TIR-NB-LRR genes and confers downy mildew resistance through multiple signalling components

In Arabidopsis, RPP4 confers resistance to Peronospora parasitica (P.p.) races Emoy2 and Emwa1 (downy mildew). We identified RPP4 in Col-0 as a member of the clustered RPP5 multigene family encoding nucleotide-binding leucine-rich repeat proteins with Toll/interleukin-1 receptor domains. RPP4 is the orthologue of RPP5 which, in addition to recognizing P.p. race Noco2, also mediates resistance to Emoy2 and Emwa1. Most differences between RPP4 and RPP5 occur in residues that constitute the TIR domain and in LRR residues that are predicted to confer recognition specificity. RPP4 requires the action of at least 12 defence components, including DTH9, EDS1, PAD4, PAL, PBS2, PBS3, SID1, SID2 and salicylic acid. The ndr1, npr1 and rps5-1 mutations partially compromise RPP4 function in cotyledons but not in true leaves. The identification of RPP4 as a TIR-NB-LRR protein, coupled with its dependence on certain signalling components in true leaves, is consistent with the hypothesis that distinct NB-LRR protein classes differentially signal through EDS1 and NDR1. Our results suggest that RPP4-mediated resistance is developmentally regulated and that in cotyledons there is cross-talk between EDS1 and NDR1 signalling and processes regulating systemic acquired resistance.     

Natural variation reveals key amino acids in a downy mildew effector that alters recognition specificity by an Arabidopsis resistance gene

RPP13, a member of the cytoplasmic class of disease resistance genes, encodes one of the most variable Arabidopsis proteins so far identified. This variability is matched in ATR13, the protein from the oomycete downy mildew pathogen Hyaloperonospora parasitica recognized by RPP13, suggesting that these proteins are involved in tight reciprocal coevolution. ATR13 exhibits five domains: an N-terminal signal peptide, an RXLR motif, a heptad leucine/isoleucine repeat, an 11-amino-acid repeated sequence and a C-terminal domain. We show that the conserved RXLR-containing domain is dispensable for ATR13-mediated recognition, consistent with its role in transport into the plant cytoplasm. Sequencing ATR13 from 16 isolates of H. parasitica revealed high levels of amino acid diversity across the entire protein. The leucines/isoleucines of the heptad leucine repeat were conserved, and mutation of particular leucine or isoleucine residues altered recognition by RPP13. Natural variation has not exploited this route to detection avoidance, suggesting a key role of this domain in pathogenicity. The extensive variation in the 11-amino-acid repeat units did not affect RPP13 recognition. Domain swap analysis showed that recognition specificity lay in the C-terminal domain of ATR13. Variation analyses combined with functional assays allowed the identification of four amino acid positions that may play a role in recognition specificity. Site-directed mutagenesis confirmed that a threonine residue is absolutely required for RPP13 recognition and that recognition can be modulated by the presence of either an arginine or glutamic acid at other sites. Mutations in these three amino acids had no effect on the interaction of ATR13 with a resistance gene unlinked to RPP13, consistent with their critical role in determining RPP13-Nd recognition specificity.     

Structural elucidation and functional characterization of the Hyaloperonospora arabidopsidis effector protein ATR13

The oomycete Hyaloperonospora arabidopsidis (Hpa) is the causal agent of downy mildew on the model plant Arabidopsis thaliana and has been adapted as a model system to investigate pathogen virulence strategies and plant disease resistance mechanisms. Recognition of Hpa infection occurs when plant resistance proteins (R-genes) detect the presence or activity of pathogen-derived protein effectors delivered to the plant host. This study examines the Hpa effector ATR13 Emco5 and its recognition by RPP13-Nd, the cognate R-gene that triggers programmed cell death (HR) in the presence of recognized ATR13 variants. Herein, we use NMR to solve the backbone structure of ATR13 Emco5, revealing both a helical domain and a disordered internal loop. Additionally, we use site-directed and random mutagenesis to identify several amino acid residues involved in the recognition response conferred by RPP13-Nd. Using our structure as a scaffold, we map these residues to one of two surface-exposed patches of residues under diversifying selection. Exploring possible roles of the disordered region within the ATR13 structure, we perform domain swapping experiments and identify a peptide sequence involved in nucleolar localization. We conclude that ATR13 is a highly dynamic protein with no clear structural homologues that contains two surface-exposed patches of polymorphism, only one of which is involved in RPP13-Nd recognition specificity.     

Maintenance of genetic variation in plants and pathogens involves complex networks of gene-for-gene interactions

The RPP13 [recognition of Hyaloperonospora arabidopsidis (previously known as Peronospora parasitica )] resistance ( R ) gene in Arabidopsis thaliana exhibits the highest reported level of sequence diversity among known R genes. Consistent with a co-evolutionary model, the matching effector protein ATR13 ( A. thaliana -recognized) from H. arabidopsidis reveals extreme levels of allelic diversity. We isolated 23 new RPP13 sequences from a UK metapopulation, giving a total of 47 when combined with previous studies. We used these in functional studies of the A. thaliana accessions for their resistance response to 16 isolates of H. arabidopsidis . We characterized the molecular basis of recognition by the expression of the corresponding ATR13 genes from these 16 isolates in these host accessions. This allowed the determination of which alleles of RPP13 were responsible for pathogen recognition and whether recognition was dependent on the RPP13/ATR13 combination. Linking our functional studies with phylogenetic analysis, we determined that: (i) the recognition of ATR13 is mediated by alleles in just a single RPP13 clade; (ii) RPP13 alleles in other clades have evolved the ability to detect other pathogen ATR protein(s); and (iii) at least one gene, unlinked to RPP13 in A. thaliana, detects a different subgroup of ATR13 alleles.     

The downy mildew effector proteins ATR1 and ATR13 promote disease susceptibility in Arabidopsis thaliana

The downy mildew (Hyaloperonospora parasitica) effector proteins ATR1 and ATR13 trigger RPP1-Nd/WsB- and RPP13-Nd-dependent resistance, respectively, in Arabidopsis thaliana. To better understand the functions of these effectors during compatible and incompatible interactions of H. parasitica isolates on Arabidopsis accessions, we developed a novel delivery system using Pseudomonas syringae type III secretion via fusions of ATRs to the N terminus of the P. syringae effector protein, AvrRPS4. ATR1 and ATR13 both triggered the hypersensitive response (HR) and resistance to bacterial pathogens in Arabidopsis carrying RPP1-Nd/WsB or RPP13-Nd, respectively, when delivered from P. syringae pv tomato (Pst) DC3000. In addition, multiple alleles of ATR1 and ATR13 confer enhanced virulence to Pst DC3000 on susceptible Arabidopsis accessions. We conclude that ATR1 and ATR13 positively contribute to pathogen virulence inside host cells. Two ATR13 alleles suppressed bacterial PAMP (for Pathogen-Associated Molecular Patterns)-triggered callose deposition in susceptible Arabidopsis when delivered by DC3000 DeltaCEL mutants. Furthermore, expression of another allele of ATR13 in plant cells suppressed PAMP-triggered reactive oxygen species production in addition to callose deposition. Intriguingly, although Wassilewskija (Ws-0) is highly susceptible to H. parasitica isolate Emco5, ATR13Emco5 when delivered by Pst DC3000 triggered localized immunity, including HR, on Ws-0. We suggest that an additional H. parasitica Emco5 effector might suppress ATR13-triggered immunity.     

Multiple candidate effectors from the oomycete pathogen Hyaloperonospora arabidopsidis suppress host plant immunity

Oomycete pathogens cause diverse plant diseases. To successfully colonize their hosts, they deliver a suite of effector proteins that can attenuate plant defenses. In the oomycete downy mildews, effectors carry a signal peptide and an RxLR motif. Hyaloperonospora arabidopsidis (Hpa) causes downy mildew on the model plant Arabidopsis thaliana (Arabidopsis). We investigated if candidate effectors predicted in the genome sequence of Hpa isolate Emoy2 (HaRxLs) were able to manipulate host defenses in different Arabidopsis accessions. We developed a rapid and sensitive screening method to test HaRxLs by delivering them via the bacterial type-three secretion system (TTSS) of Pseudomonas syringae pv tomato DC3000-LUX (Pst-LUX) and assessing changes in Pst-LUX growth in planta on 12 Arabidopsis accessions. The majority (~70%) of the 64 candidates tested positively contributed to Pst-LUX growth on more than one accession indicating that Hpa virulence likely involves multiple effectors with weak accession-specific effects. Further screening with a Pst mutant (ΔCEL) showed that HaRxLs that allow enhanced Pst-LUX growth usually suppress callose deposition, a hallmark of pathogen-associated molecular pattern (PAMP)-triggered immunity (PTI). We found that HaRxLs are rarely strong avirulence determinants. Although some decreased Pst-LUX growth in particular accessions, none activated macroscopic cell death. Fewer HaRxLs conferred enhanced Pst growth on turnip, a non-host for Hpa, while several reduced it, consistent with the idea that turnip's non-host resistance against Hpa could involve a combination of recognized HaRxLs and ineffective HaRxLs. We verified our results by constitutively expressing in Arabidopsis a sub-set of HaRxLs. Several transgenic lines showed increased susceptibility to Hpa and attenuation of Arabidopsis PTI responses, confirming the HaRxLs' role in Hpa virulence. This study shows TTSS screening system provides a useful tool to test whether candidate effectors from eukaryotic pathogens can suppress/trigger plant defense mechanisms and to rank their effectiveness prior to subsequent mechanistic investigation.     

Computational prediction and molecular characterization of an oomycete effector and the cognate Arabidopsis resistance gene

Hyaloperonospora arabidopsidis (Hpa) is an obligate biotroph oomycete pathogen of the model plant Arabidopsis thaliana and contains a large set of effector proteins that are translocated to the host to exert virulence functions or trigger immune responses. These effectors are characterized by conserved amino-terminal translocation sequences and highly divergent carboxyl-terminal functional domains. The availability of the Hpa genome sequence allowed the computational prediction of effectors and the development of effector delivery systems enabled validation of the predicted effectors in Arabidopsis. In this study, we identified a novel effector ATR39-1 by computational methods, which was found to trigger a resistance response in the Arabidopsis ecotype Weiningen (Wei-0). The allelic variant of this effector, ATR39-2, is not recognized, and two amino acid residues were identified and shown to be critical for this loss of recognition. The resistance protein responsible for recognition of the ATR39-1 effector in Arabidopsis is RPP39 and was identified by map-based cloning. RPP39 is a member of the CC-NBS-LRR family of resistance proteins and requires the signaling gene NDR1 for full activity. Recognition of ATR39-1 in Wei-0 does not inhibit growth of Hpa strains expressing the effector, suggesting complex mechanisms of pathogen evasion of recognition, and is similar to what has been shown in several other cases of plant-oomycete interactions. Identification of this resistance gene/effector pair adds to our knowledge of plant resistance mechanisms and provides the basis for further functional analyses.     

Arabidopsis downy mildew resistance gene RPP27 encodes a receptor-like protein similar to CLAVATA2 and tomato Cf-9

The Arabidopsis Ler-RPP27 gene confers AtSgt1b-independent resistance to downy mildew (Peronospora parasitica) isolate Hiks1. The RPP27 locus was mapped to a four-bacterial artificial chromosome interval on chromosome 1 from genetic analysis of a cross between the enhanced susceptibility mutant Col-edm1 (Col-sgt1) and Landsberg erecta (Ler-0). A Cf-like candidate gene in this interval was PCR amplified from Ler-0 and transformed into mutant Col-rpp7.1 plants. Homozygous transgenic lines conferred resistance to Hiks1 and at least four Ler-0 avirulent/Columbia-0 (Col-0) virulent isolates of downy mildew pathogen. A full-length RPP27 cDNA was isolated, and analysis of the deduced amino acid sequences showed that the gene encodes a receptor-like protein (RLP) with a distinct domain structure, composed of a signal peptide followed by extracellular Leu-rich repeats, a membrane spanning region, and a short cytoplasmic carboxyl domain. RPP27 is the first RLP-encoding gene to be implicated in disease resistance in Arabidopsis, enabling the deployment of Arabidopsis techniques to investigate the mechanisms of RLP function. Homology searches of the Arabidopsis genome, using the RPP27, Cf-9, and Cf-2 protein sequences as a starting point, identify 59 RLPs, including the already known CLAVATA2 and TOO MANY MOUTHS genes. A combination of sequence and phylogenetic analysis of these predicted RLPs reveals conserved structural features of the family.     

White rust (Albugo candida) resistance loci on three Arabidopsis chromosomes are closely linked to downy mildew (Peronospora parasitica) resistance loci

Two accessions of Arabidopsis thaliana (Ksk-1 and Ksk-2) were used to identify and map three loci ( RAC1 , RAC2 and RAC3 ) of genes that confer R esistance to A lbugo c andida (white rust). The phenotypes associated with these genes were classified as either FN (necrotic f lecks on upper surface of cotyledons and n o blisters) for RAC2 and RAC3, or FYN ( f lecks surrounded by y ellowing and n o blisters) for RAC1 . Both phenotypes exhibited rapid death of host cells penetrated by the parasite (hypersensitive response), with callose deposition commonly encasing the haustorium. F₆ recombinant inbred lines were produced specifically for the purpose of mapping each RAC locus relative to molecular markers. Dominant resistance at the locus RAC1 in Ksk-1 was previously mapped to chromosome 1 between RFLP markers m253 and m254, and co-segregating with a downy mildew resistance specificity RPP9 in the accession Wei-0. We report here a fine-scale map interval and co-segregating markers for this locus, which in turn enabled mapping of a previously unnoticed source of resistance in Ksk-1 designated RAC3 that exhibits an FN phenotype hyperstatic to the FYN phenotype of RAC1. RAC3 is closely linked to the RPP8/HRT on chromosome 5, a locus which contains specificities for resistance to downy mildew and turnip crinkle virus. Recombinant inbreds also enabled mapping of recessive resistance at RAC2 in Ksk-2 to the bottom arm of chromosome 3, in the 6 cM interval between two downy mildew resistance loci ( RPP1 and RPP13 ) .     

Analysis of expressed sequence tags derived from pea leaves infected by Peronospora viciae f. sp. pisi

A cDNA library was constructed from field pea leaves infected by the downy mildew pathogen, Peronospora viciae f. sp. pisi, using a suppression subtractive hybridisation approach. The library consists of 399 expressed sequence tags, from which 207 unisequences were obtained after sequence assembly. Of the unisequences, six were shown to be of Peronospora viciae f. sp. pisi origin. The remaining unisequences were subjected to gene ontology analysis and their functions were predicted in silico. Eleven of these unisequences (representing 24 clones) shared significant sequence similarities with Arabidopsis genes known to be involved in downy mildew resistance, including the well‐characterised genes RPP5, RPP6 and RPP27. Expression analysis of five selected unisequences by real‐time PCR indicated that all five were up‐regulated during downy mildew pathogenesis, suggesting a significant role for these genes in the host response to downy mildew infection.     

454 Genome sequencing of Pseudoperonospora cubensis reveals effector proteins with a QXLR translocation motif

Pseudoperonospora cubensis is a biotrophic oomycete pathogen that causes downy mildew of cucurbits, a devastating foliar disease threatening cucurbit production worldwide. We sequenced P. cubensis genomic DNA using 454 pyrosequencing and obtained random genomic sequences covering approximately 14% of the genome, thus providing the first set of useful genomic sequence information for P. cubensis. Using bioinformatics approaches, we identified 32 putative RXLR effector proteins. Interestingly, we also identified 29 secreted peptides with high similarity to RXLR effectors at the N-terminal translocation domain, yet containing an R-to-Q substitution in the first residue of the translocation motif. Among these, a family of QXLR-containing proteins, designated as PcQNE, was confirmed to have a functional signal peptide and was further characterized as being localized in the plant nucleus. Internalization of secreted PcQNE into plant cells requires the QXLR-EER motif. This family has a large number of near-identical copies within the P. cubensis genome, is under diversifying selection at the C-terminal domain, and is upregulated during infection of plants, all of which are common characteristics of characterized oomycete effectors. Taken together, the data suggest that PcQNE are bona fide effector proteins with a QXLR translocation motif, and QXLR effectors are prevalent in P. cubensis. Furthermore, the massive duplication of PcQNE suggests that they might play pivotal roles in pathogen fitness and pathogenicity.     

A genetic interval and physical contig spanning the Peronospora parasitica (At) avirulence gene locus ATR1Nd

In Peronospora parasitica (At) (downy mildew), the genetic determinants of cultivar-specific recognition by Arabidopsis thaliana are the Arabidopsis thaliana-recognised (ATR) avirulence genes. We describe the identification of 10 amplified fragment length polymorphism (AFLP) markers that define a genetic mapping interval for the ATR1Nd avirulence allele, the presence of which is perceived by the RPP1Nd resistance gene. Furthermore, we have constructed a P. parasitica (At) bacterial artificial chromosome (BAC) library comprising over 630Mb of cloned DNA. We have isolated 16 overlapping clones from the BAC library that form a contig spanning the genetic interval. BAC sequence-derived markers and a total mapping population of 311 F(2) individuals were used to refine the ATR1Nd locus to a 1cM interval that is represented by four BAC clones and spans less than 250kb of DNA. This work demonstrates that map-based cloning techniques are feasible in this organism and provides the critical foundations for cloning ATR1Nd using such a strategy.     

Functional Analysis of Hyaloperonospora arabidopsidis RXLR Effectors

The biotrophic plant pathogen Hyaloperonospora arabidopsidis produces a set of putative effector proteins that contain the conserved RXLR motif. For most of these RXLR proteins the role during infection is unknown. Thirteen RXLR proteins from H. arabidopsidis strain Waco9 were analyzed for sequence similarities and tested for a role in virulence. The thirteen RXLR proteins displayed conserved N-termini and this N-terminal conservation was also found in the 134 predicted RXLR genes from the genome of H. arabidopsidis strain Emoy2. To investigate the effects of single RXLR effector proteins on plant defense responses, thirteen H. arabidopsidis Waco9 RXLR genes were expressed in Arabidopsis thaliana. Subsequently, these plants were screened for altered susceptibility to the oomycetes H. arabidopsidis and Phytophthora capsici, and the bacterial pathogen Pseudomonas syringae. Additionally, the effect of the RXLR proteins on flg22-triggered basal immune responses was assessed. Multifactorial analysis of results collated from all experiments revealed that, except for RXLR20, all RXLR effector proteins tested affected plant immunity. For RXLR9 this was confirmed using a P. syringae ΔCEL-mediated effector delivery system. Together, the results show that many H. arabidopsidis RXLR effectors have small effects on the plant immune response, suggesting that suppression of host immunity by this biotrophic pathogen is likely to be caused by the combined actions of effectors.     

RPP13 is a simple locus in Arabidopsis thaliana for alleles that specify downy mildew resistance to different avirulence determinants in Peronospora parasitica

Disease resistance (R) genes are found in plants as either simple (single allelic series) loci, or more frequently as complex loci of tandemly repeated genes. These different loci are likely to be under similar evolutionary forces from pathogens, but the contrast between them suggests important differences in mechanisms associated with DNA structure and recombination that generate and maintain R gene diversity. The RPP13 locus in Arabidopsis represents an important paradigm for studying the evolution of an R gene at a simple locus. The RPP13 allele from the accession Nd-1, designated RPP13-Nd, confers resistance to five different isolates of the biotrophic oomycete, Peronospora parasitica (causal agent of downy mildew), and encodes an NBS-LRR type R protein with a putative amino-terminal leucine zipper. The RPP13-Rld allele, cloned from the accession Rld-2, encodes a different specificity. Comparison of three RPP13 alleles revealed a high rate of amino acid divergence within the LRR domain, less than 80% identity overall, compared to the remainder of the protein (> 95% identity). We also found evidence for positive selection in the LRR domain for amino acid diversification outside the core conserved beta-strand/beta-turn motif, suggesting that more of the LRR structure is available for interaction with target molecules than has previously been reported for other R gene products. Furthermore, an amino acid sequence (LLRVLDL) identical in an LRR among RPP13 alleles is conserved in other LZ NBS-LRR type R proteins, suggesting functional significance.     

Identification of non-TIR-NBS-LRR markers linked to the <Emphasis Type="Italic">Pl5</Emphasis>/<Emphasis Type="Italic">Pl8</Emphasis> locus for resistance to downy mildew in sunflower

The resistance of sunflower, Helianthus annuus L., to downy mildew, caused by Plasmopara halstedii, is conferred by major genes denoted by Pl. Using degenerate and specific primers, 16 different resistance gene analogs (RGAs) have been cloned and sequenced. Sequence comparison and Southern-blot analysis distinguished six classes of RGA. Two of these classes correspond to TIR-NBS-LRR sequences while the remaining four classes correspond to the non-TIR-NBS-LRR type of resistance genes. The genetic mapping of these RGAs on two segregating F2 populations showed that the non-TIR-NBS-LRR RGAs are clustered and linked to the Pl5/ Pl8 locus for resistance to downy mildew in sunflower. These and other results indicate that different Pl loci conferring resistance to the same pathogen races may contain different sequences.