# Effector prediction and characterization in the oomycete pathogen Bremia lactucae reveal host-recognized WY domain proteins that lack the canonical RXLR motif

Kelsey Wood, Munir Nur, Juliana D B Gil, Kyle Fletcher, Kim Lakeman, Dasan Gann, Ayumi Gothberg, Tina Khuu, Jennifer Kopetzky, Sanye Naqvi, et al.

<figure><img src="image_1.png" /></figure>

## To cite this version:

Kelsey Wood, Munir Nur, Juliana D B Gil, Kyle Fletcher, Kim Lakeman, et al. Effector prediction and characterization in the oomycete pathogen Bremia lactucae reveal host-recognized WY domain proteins that lack the canonical RXLR motif. PLoS Pathogens, 2020, 16 (10), pp.e1009012. 10.1371/journal.ppat.1009012. hal-03106412

HAL Id: hal-03106412

https://hal.inrae.fr/hal-03106412

Submitted on 11 Jan 2021

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>
RESEARCH ARTICLE

# Effector prediction and characterization in the oomycete pathogen *Bremia lactucae* reveal host-recognized WY domain proteins that lack the canonical RXLR motif

Kelsey J. Wood<sup>1,2</sup>, Munir Nur<sup>1</sup>, Juliana Gil<sup>1,3</sup>, Kyle Fletcher<sup>1</sup>, Kim Lakeman<sup>4</sup>, Dasan Gann<sup>1</sup>, Ayumi Gothberg<sup>1</sup>, Tina Khuu<sup>1</sup>, Jennifer Kopetzky<sup>1</sup>, Sanye Naqvi<sup>1,5</sup>, Archana Pandya<sup>1</sup>, Chi Zhang<sup>1</sup>, Brigitte Maisonneuve<sup>1,5</sup>, Mathieu Pel<sup>4</sup>, Richard Michelmore<sup>1,6*</sup>

1 The Genome Center, University of California, Davis, Davis, California, United States of America,

2 Integrative Genetics & Genomics Graduate Group, University of California, Davis, Davis, California, United States of America, 3 Plant Pathology Graduate Group, University of California, Davis, Davis, California, United States of America, 4 Enza Zaden, Enkhuizen, The Netherlands, 5 INRAE, UR1052, GAFL, Montfavet, France, 6 Departments of Plant Sciences, Molecular & Cellular Biology, Medical Microbiology & Immunology, University of California, Davis, Davis, California, United States of America

* rwmichelmore@ucdavis.edu

<figure><img src="image_3.png" /></figure>

<figure><img src="image_4.png" /></figure>

OPEN ACCESS

**Citation:** Wood KJ, Nur M, Gil J, Fletcher K, Lakeman K, Gann D, et al. (2020) Effector prediction and characterization in the oomycete pathogen *Bremia lactucae* reveal host-recognized WY domain proteins that lack the canonical RXLR motif. PLoS Pathog 16(10): e1009012. https://doi.org/10.1371/journal.ppat.1009012

**Editor:** Mark J. Banfield, John Innes Centre, UNITED KINGDOM

**Received:** June 15, 2020

**Accepted:** September 29, 2020

**Published:** October 26, 2020

**Peer Review History:** PLOS recognizes the benefits of transparency in the peer review process; therefore, we enable the publication of all of the content of peer review and author responses alongside final, published articles. The editorial history of this article is available here: https://doi.org/10.1371/journal.ppat.1009012

**Copyright:** © 2020 Wood et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** The reference genome assembly of *B. lactucae* is available from NCBI, GenBank ID: GCA_004359215.1. RNAseq

## Abstract

Pathogens that infect plants and animals use a diverse arsenal of effector proteins to suppress the host immune system and promote infection. Identification of effectors in pathogen genomes is foundational to understanding mechanisms of pathogenesis, for monitoring field pathogen populations, and for breeding disease resistance. We identified candidate effectors from the lettuce downy mildew pathogen *Bremia lactucae* by searching the predicted proteome for the WY domain, a structural fold found in effectors that has been implicated in immune suppression as well as effector recognition by host resistance proteins. We predicted 55 WY domain containing proteins in the genome of *B. lactucae* and found substantial variation in both sequence and domain architecture. These candidate effectors exhibit several characteristics of pathogen effectors, including an N-terminal signal peptide, lineage specificity, and expression during infection. Unexpectedly, only a minority of *B. lactucae* WY effectors contain the canonical N-terminal RXLR motif, which is a conserved feature in the majority of cytoplasmic effectors reported in *Phytophthora* spp. Functional analysis of 21 effectors containing WY domains revealed 11 that elicited cell death on wild accessions and domesticated lettuce lines containing resistance genes, indicative of recognition of these effectors by the host immune system. Only two of the 11 recognized effectors contained the canonical RXLR motif, suggesting that there has been an evolutionary divergence in sequence motifs between genera; this has major consequences for robust effector prediction in oomycete pathogens.
reads of lettuce cv. Cobham Green infected with *B. lactucae* isolate SF5 is available at NCBI SRA
BioProject PRJNA523226. Sequences of effector proteins are available in File S1. All data for interactions between effectors and individual lines of *Lactuca* spp. are available at http://bremia.ucdavis.edu/BIL/BIL_intersection.php. The sequences of ORFs from the oomycete genomes used for effector prediction, along with scripts used for their analysis, can be found on Github at https://github.com/kelseywood/Bremia-lactucae-WY-effector-characterization.

**Funding:** The work was supported by an NSF Graduate Research Fellowship and a USDA Fellowship #2018-67011-28053 to KW and the NSF/USDA AFRI Microbial Sequencing Program award #2009-65109-05925 to RWM. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

**Competing interests:** The authors have declared that no competing interests exist.

## Author summary

There is a molecular battle that takes place during infection of plants and animals by pathogens. Some of the weapons in the pathogen's arsenal are known as " effectors ;" these are secreted proteins that enter host cells to alter their physiology and suppress the immune system. Effectors can also be a liability for pathogens because potential hosts have evolved ways to recognize these effectors, triggering a defense response leading to localized cell death, preventing the spread of the pathogen. Here we computationally predicted effectors encoded in the genome of *Bremia lactucae*, the pathogen causing the downy mildew disease of lettuce and tested 21 of these proteins in the lab to see what effects they would have on the host plant. Five effectors were demonstrated to suppress the basal immune system of lettuce. Eleven effectors caused programmed cell death when introduced to certain lettuce lines, which indicates recognition of these proteins by the host immune system. In addition to contributing to our understanding of the mechanisms of pathogenesis, this study of effectors facilitates breeding for disease resistant lettuce, which will decrease agricultural reliance on fungicides.

## Introduction

Oomycetes are some of the most devastating pathogens of both plants and animals [1]. Although oomycetes resemble fungi in their filamentous growth and infection structures, they are more closely related to brown algae than to fungi [2]. Notable oomycetes include the plant pathogens causing late blight of potato (*Phytophthora infestans*) [3], sudden oak death (*Phytophthora ramorum*) [4], root rot (*Phytophthora* [5] and *Pythium* spp. [6,7]), white blister rust of Brassica spp. (*Albugo* spp.) [8], and downy mildews (e.g., *Bremia*, *Peronospora*, *Plasmopara* spp.) [9], as well as several important animal pathogens infecting fish (*Saprolegnia* spp.), shell-fish (*Aphanomyces astaci*), and humans and other mammals (*Pythium insidiosum*) [1,10]. Many types of plant and animal pathogens, including oomycetes, secrete proteins known as effectors to promote virulence by manipulating the physiology of the host cells and by suppressing the host immune system [11]. In plants, effectors are also determinants of resistance or susceptibility through their direct or indirect interactions with cognate nucleotide-binding leucine rich repeat proteins (NLRs) encoded by resistance genes [12]. Effectors are secreted from the pathogen and may act extracellularly or they may be translocated into the cytoplasm [11]. One class of translocated effectors from plant pathogenic oomycetes of the class Peronosporales, which includes *Phytophthora* and the downy mildews, are the RXLR effectors, so called for their N-terminal motif usually consisting of arginine (R), followed by any amino acid (X), then followed by leucine (L) and arginine (R) [13]. RXLR effectors also contain an N-terminal signal peptide, which designates them for extracellular transport by way of the endoplasmic reticulum and Golgi apparatus [14]. The RXLR motif is often associated with a downstream EER motif, both of which have been associated with secretion and/or translocation of effectors into the plant cell [14,15]. For some RXLR effectors, such as Avr3a, the RXLR motif has been shown to be cleaved just prior to the EER sequence and therefore plays a role in targeting the protein for translocation rather than directly being involved in uptake into the host cell [16]. The RXLR motif is similar in sequence to the PEXEL motif (RXLX[EDQ]) of the distantly related malaria pathogen (*Plasmodium falciparum*) [17] and the TEXEL motif of *Toxo-plasma gondii* (RRLXX) [18], both of which are required for proteolytic modification in the endoplasmic reticulum and targeting effector proteins for specialized export out of the cell [19].
The downy mildews and the related *Phytophthora* species have different lifestyles, with the downy mildews requiring a living host to complete their lifecycle (obligate biotrophy) and *Phytophthora* spp. using an initial biotrophic phase followed by a necrotrophic phase (hemibiotrophy) [20]; however, their effectors share similar features. Many effectors in the Peronosporales have RXLR and EER motifs, although several alternative sequences to RXLR have been found in downy mildews including RVRN (ATR5 from *Hyaloperonospora arabidopsidis*) [21], QXLR (*Pseudoperonospora cubensis*) [22], GKLR (*Bremia lactucae*) [23,24], and RXLK (*Plasmopara halstedii*) [25]. The C-terminal effector domains of RXLR effectors from *Phytophthora* and downy mildews also share some common sequence motifs and structural features, such as the 24 to 30 amino acid W, Y, and L motifs, which were first identified bioinformatically [26]. Structural analysis on four different RXLR effectors from *Phytophthora infestans* (Avr3a and PexRD2), *P. capsici* (Avr3a11), and the downy mildew pathogen *H. arabidopsidis* (ATR1) revealed that the W and Y motifs form a ~50 amino acid alpha-helical fold that may play a role in protein-protein interactions [27–30]. This effector-associated fold, termed the WY domain after its conserved tryptophan and tyrosine residues, is highly structurally conserved between effectors from multiple Peronosporales species, despite sharing less than 20% sequence similarity across the whole domain [27,31]. The WY domain appears to be specific to the Peronosporales and was predicted to be present in nearly half of the RXLR effectors of *P. infestans* and a quarter of the RXLR effectors in *H. arabidopsidis* [27].

Functional studies of WY domain containing proteins have indicated that certain residues in the WY domain are essential for the immune suppressing functions of *P. sojae* effectors Avr1b [32] and PSR1 [33], and *P. infestans* effectors Avr3a [34] and PexRD2 [35]. Furthermore, mutation of two conserved leucines in the WY domain of PexRD2 disrupted interaction with its target, MAPKKKε, consistent with this domain being important for protein-protein interactions [35]. WY domain containing proteins appear to interact with a variety of host targets, with Avr3a from *P. infestans* targeting the E3 ligase CMPG1 [34], *P. infestans* PexRD54 targeting potato autophagy-related protein ATG8 [36], PsAvh240 from *P. sojae* targeting an aspartic protease [37], and PSR1 and PSR2 from *P. sojae* suppressing host RNA silencing through interactions with the RNA helicase PINP1 and the dsRNA-binding protein DRB4, respectively [38–41]. Mutation analysis of the regions encoding the seven individual WY domains of PSR2 demonstrated differential contributions of each domain to virulence of *P. sojae*, suggesting that the WY domain may act as a module during effector evolution [31]. In addition to its roles in immune suppression, the WY domain has been shown to be important for immune recognition of the effector by nucleotide binding-leucine rich repeat (NB-LRR) resistance proteins [42].

Effector annotation in oomycete genomes has often relied on sequence similarity to known effectors or on prediction of conserved motifs, such as the RXLR motif, or in the case of Crinklers, the LXLFLAK motif [3]. Due to the short length and degeneracy of the RXLR sequence, the motif occurs frequently by chance; therefore, there is a high false positive rate (>50%) using string-based searches [26]. Hidden Markov Model (HMM)-based searches have much lower false positive rates, but the false negative rate may be higher if the genome of interest has diverged significantly from the species used to build the HMM. Downy mildews have a narrow host range and pathogenicity-related proteins are likely to show high lineage-specificity due to co-evolution with their hosts. There is already evidence that downy mildew effectors show divergence from the canonical RXLR motif [21,24,25]; thus, complementary approaches for effector prediction that utilize other conserved features, such as the WY domain, are necessary to fully characterize the repertoire. In addition, a WY domain containing protein from *Plasmopara viticola* lacking the RXLR motif has been found to trigger cell death in grapevine and tobacco species, suggesting that non-RXLR WY proteins may be recognized by resistance
proteins [43]. In addition, due to the importance of the WY domain in effector function [32,34,35], WY domain containing proteins may be better candidates for effectors than those containing only the RXLR motif.

To identify a more complete repertoire of candidate effectors in the reference genome of *B. lactucae* and to test whether the WY domain is informative for predicting effectors in the Pero-nosporales, we searched for this domain using an HMM built from sequences of the WY domain in three *Phytophthora* species [28]. This revealed additional effector candidates that had not been found using an RXLR-based search; similar results were also obtained for other downy mildews and well-studied *Phytophthora* species. These predicted WY proteins from *B. lactucae* had other signatures of oomycete effectors, such as presence of a secretion signal, N-terminal intrinsic disorder, lineage specificity, and expression during infection. A subset of the predicted WY effectors suppressed the host immune system, while others elicited programmed cell death in specific genotypes of lettuce, indicative of recognition by host resistance proteins. Therefore, searches for the WY domains are highly useful for identification of downy mildew effectors lacking the RXLR motif, which was previously considered canonical for effectors of the Peronosporales.

# Results

## Hidden Markov Model-based prediction of the WY domain identifies candidate effectors lacking the RXLR motif in *B. lactucae* and other oomycetes

Our HMM search initially identified 59 candidate WY effectors encoded by the gene models predicted in the *B. lactucae* SF5 assembly [44]; however, four pairs of genes appeared to be allelic based on sequence similarity and read depth, leaving a total of 55 non-redundant genes encoding candidate effectors (Fig 1A; S1 File). N-terminal signal peptides were predicted for 41 of these 55 proteins (Fig 1A), of which two proteins had a predicted transmembrane helix outside of the signal peptide. Several predicted WY containing proteins seemed to be missing their signal peptides due to N-terminal truncation when compared to close relatives. One of these predicted proteins, BLN06, was found to be missing 72 amino acids from the N-terminal sequence, including the signal peptide, in SF5 compared to BL24, which was the source isolate for the cloning of BLN06 reported by Pelgrom *et al.* [45]. BSW04p*, which appears to be allelic to BSW04p based on read depth, also has an N-terminal truncation resulting in loss of the signal peptide. All candidate WY effectors, including those without a signal peptide, are expressed during the lifecycle of *B. lactucae* as evidenced by RNAseq reads (S1 File) with expression confirmed by qRT-PCR for the 21 cloned genes (S1 Fig).

The N-terminal sequences of the 55 predicted WY effectors were examined for the RXLR motif using both HMMs and string searches. One protein was predicted to have an RXLR motif by the HMM and an additional 10 proteins were identified by a string search for [RQGH]XLR, while 33 proteins were predicted to have an EER motif using a string search for [DE][DE][RK] (Fig 1A). To find divergent RXLR-like motifs in the WY effector candidates, we searched for a highly degenerate pattern based on mutational studies of the RXLR motif [35] and natural RXLR variants reported for other downy mildews [5-9]. This revealed an additional 40 proteins with an RXLR-like motif within the first 60 amino acids after the signal peptide matching the pattern [RKHGQ][X]{0,1}[LMYFIV][RNK], many of which also had an EER motif (Fig 1A) (S1 Table). Therefore, the majority of candidate WY effectors in *B. lactucae* have a non-canonical RXLR motif.

In the C-terminal domain, the WY effector candidates had 1-8 WY domains per protein, with a diversity of domain architectures (Fig 1). The WY domain from *B. lactucae* showed high conservation of the characteristic conserved tryptophan (W) residue but appeared to
<figure><img src="image_5.png" /><figcaption>**Fig 1.** WY effector candidates from *B. lactucae* isolate SF5. (A) UPGMA consensus tree of the 59 predicted WY effectors from *B. lactucae* isolate SF5 based on whole protein amino acid sequence alignment using MUSCLE. Four sequences that appear to be allelic are indicated with asterisks next to the sequence name. Branch lengths have been transformed for clarity of presentation; however, bootstrap values (n = 100 iterations) and branch lengths are given for closely related proteins. Signal peptides are shown as circles (red circle for SP; white circle for no SP). RXLR motifs are shown as triangles (black triangle for RXLR ([RQGH]XLR or RXL [QKG]); grey triangle for degenerate RXLR ([RKHGQ][X]{0,1}[LMYFIV][RNK]), white triangle for no RXLR-like sequence) followed by an inverted triangle for EER motifs (yellow triangle for [DE][DE][RK], white triangle for no EER-like sequence). WY domain architecture is shown using blue rectangles, with a black line representing the total length of the protein and rectangle position representing the location of the WY motif as predicted by HMMer 3.1. (B) Sequence logo for the WY domain from *B. lactucae* built from a multiple sequence alignment of WY domains predicted by HMMer 3.1.</figcaption></figure>

https://doi.org/10.1371/journal.ppat.1009012.g001
show equal preference for tyrosine (Y) or phenylalanine (F) for the second characteristic residue of the domain (Fig 1B).

## Many candidate WY effectors across oomycetes lack the RXLR motif

To investigate whether or not there are WY effector candidates that are lacking the canonical RXLR motif in other oomycetes, the predicted open reading frames (ORFs) from several published oomycete genomes were surveyed for the RXLR and EER motifs and the WY domain. Approximately half of the WY proteins predicted in seven downy mildew species lacked the RXLR motif; in six *Phytophthora* spp., the majority of predicted WY proteins had RXLR motifs, but between 9–21% of secreted WY proteins did not contain the RXLR motif (Fig 2; S2 Table). The majority of WY proteins lacking the RXLR motif contained the EER motif (Fig 2; S2 Table). Sequences of RXLR, EER, and WY-containing predicted proteins from all oomycete species analyzed can be found in S2 File. This analysis suggests that the repertoire of candidate WY effectors in other oomycetes may be heavily under-reported.

## Candidate WY effectors have high levels of N-terminal intrinsic disorder

Intrinsic disorder had previously been reported to be a characteristic of the N-terminus of oomycete effectors containing the RXLR motif [46]. Therefore, we investigated whether the predicted degree of structural disorder was a characteristic of candidate WY effectors lacking a

<figure><img src="image_6.png" /><figcaption>**Fig 2.** Distribution of predicted secreted WY effectors with or without RXLR and/or EER motifs in the genomes of downy mildew pathogens and *Phytophthora* species. HMMer was used to search predicted secretomes for each species for the WY domain. For the WY domain containing proteins, the presence of RXLR and EER were determined by searching for [RQGH] XLR within the first 60 amino acids after the signal peptide and [DE][DE][KR] within the first 100 amino acids after the signal peptide. For the non-WY domain containing proteins, the RXLR and EER was determined by searching for a strict "RXLR" with [DE][DE][KR], within the first 60 and 100 amino acids, respectively, plus proteins found by searching for the RXLR-EER domain using the RXLR-EER HMM from [15].</figcaption></figure>

https://doi.org/10.1371/journal.ppat.1009012.g002
<figure><img src="image_7.png" /><figcaption>**Fig 3.** Intrinsic disorder in the first 150 N-terminal amino acid sequences of RXLR and/or WY containing candidate effectors in *B. lactucae*. Proteins were categorized as WY with no RXLR (yellow), WYs with RXLR (black), RXLR+EER (with or without WY, purple), and the total predicted secretome (blue). RXLR and EER motifs were as described in Fig 2. Average positional disorder scores [46] with standard error are shown for each group of proteins.</figcaption></figure>

https://doi.org/10.1371/journal.ppat.1009012.g003

canonical RXLR. The predicted levels of intrinsic disorder were calculated for proteins containing RXLR motifs, for proteins containing WY domains but no RXLR motif, as well as for the entire predicted secretome for comparison. Proteins containing RXLR and/or WY domains had higher levels of intrinsic disorder at their N-termini after the highly ordered signal peptide than the entire set of secreted proteins (Fig 3). Proteins containing a WY domain but lacking a canonical RXLR motif had on average more disordered N-termini than effectors that had RXLR but not a WY domain (Fig 3). The WY domain containing region had higher levels of predicted structure than the RXLR-containing proteins that lacked a WY domain and the secreted proteins as a whole (Fig 3), consistent with the WY domains forming an α-helical bundled structure [28]. This predicted pattern of high intrinsic disorder at the N-terminus and high structure towards the C-terminus was consistently observed in all six downy mildews and six *Phytophthora* species analyzed (S2 Fig). Therefore, a high level of intrinsic disorder is a consistent characteristic of the N-termini of oomycete effectors, regardless of whether they have a canonical RXLR motif; the functional significance of this remains to be investigated.

## Candidate WY effectors show lineage specificity to *B. lactucae*

To evaluate lineage specificity of effectors due to co-evolution of pathogens with their hosts, we used BLAST to identify orthologs in other oomycete species. All of the 41 predicted
<figure><img src="image_8.png" /><figcaption>**Fig 4.** Lineage specificity of predicted secreted WY effectors from *B. lactucae*. Secreted WY proteins predicted in isolate SF5 of *B. lactucae* were used as a query for BLASTp against other oomycete translated ORFs. Isolate C82P24 was queried for *B. lactucae*. The best BLAST hit percent identity for each protein was calculated. The box plot shows the distribution of the best BLAST hits per *B. lactucae* WY protein from each species, with each dot representing an individual data point. No hits were observed to *Albugo laibachii*, *Saprolegnia parasitica*, or *Pythium ultimum*.</figcaption></figure>

https://doi.org/10.1371/journal.ppat.1009012.g004

secreted candidate WY effectors of *B. lactucae* had little sequence similarity to sequences in other genomes with the best BLASTp hit having only 46% identity with a protein from *P. infestans* (Fig 4). Most of the proteins had best-hit identities between 20 to 30%, which is similar to the level of amino acid conservation between WY domains in different effectors [28]. All of the proteins had hits to *B. lactucae* isolate C82P24, with very high identity, indicating conservation of these candidate effectors in *B. lactucae* (Fig 4).

## Candidate WY effectors show distinct subcellular localizations when transiently expressed in lettuce

Eight *B. lactucae* WY effectors were tagged at their N-terminus with either yellow fluorescent protein (YFP) or mCherry and their subcellular localizations were visualized using confocal microscopy (Fig 5). Five of the WY effectors (BLN08, BSW03, BLN06, BSW11, and BSW13) were localized to the cytoplasm and/or plasma membrane and excluded from the nucleus. BSW04p was localized to the nucleus, consistent with its C-terminal nuclear localization signal (NucPred [47] score 0.79). BSW19 showed both nuclear and cytoplasmic localization and also showed localization to round bodies in the cell. BSW14 was excluded from the nucleus and appeared to be localized to round bodies. Full length tagged proteins were observed for all eight tagged WY effectors by Western blot (S3 Fig). Some effectors showed cleavage of
<figure><img src="image_9.png" /><figcaption>**Fig 5.** Confocal microscopy of fluorescently-tagged WY containing candidate effectors from *B. lactucae* transiently expressed in *lettuce*. Panels A-C show N-terminal YFP-tagged effectors (yellow), DAPI nuclear stain (blue), and chlorophyll autofluorescence (red). Panels D-H show N-terminal mCherry-tagged effectors (red) and chlorophyll autofluorescence (blue). Panel I shows free GFP (green) and chlorophyll autofluorescence (red). All squares are 200 x 200 µM.</figcaption></figure>

https://doi.org/10.1371/journal.ppat.1009012.g005

mCherry or YFP on the Western blot (S3 Fig); however, since their localization patterns are distinct from free GFP (Fig 5), this degradation likely occurred during protein extraction.

**Several candidate WY effectors suppress PTI-induced ROS activity**

We tested candidate WY effectors for their ability to suppress PAMP-triggered immunity (PTI). The cloned candidate WY effectors were transiently expressed in *Nicotiana benthami-ana* and the level of reactive oxygen species (ROS) production induced by flg22 was measured. Several effectors significantly suppressed ROS production to a similar extent as two known bacterial suppressors of PTI (Fig 6). Interestingly, the addition of N-terminal fluorescent tags on some of the effectors influenced their ability to suppress PTI. For BSW04p and BSW13, the
<figure><img src="image_10.png" /><figcaption>**Fig 6.** Pattern-triggered immunity suppression ability of candidate WY containing effectors from *B. lactucae* as measured by *flg22-triggered reactive oxygen species burst*. Bar plot showing averages and standard error in PTI response, measured by the area under the curve (AUC) of luminescence during the 40-minute ROS assay, for leaf discs expressing each effector relative to leaf discs expressing the GFP control. DsRed was used as an additional negative control. *P. syringae* Hop effectors were used as positive controls for PTI suppression (purple bars). Tagged effectors are indicated by mCh (mCherry, red bars) or YFP (yellow fluorescent protein, yellow bars). Asterisks indicate PTI responses that are statistically different from the GFP control (t-test, *p* < 0.05).</figcaption></figure>

https://doi.org/10.1371/journal.ppat.1009012.g006

addition of YFP or mCherry tags increased the suppression ability, possibly by providing higher expression or increased protein stability; for BSW14 and BSW19, the addition of YFP and mCherry tags reduced the ability of these effectors to suppress PTI (Fig 6). All fluorescently-tagged proteins were confirmed to be full length by Western blotting, although all effectors, except for BSW04p, also showed a band representing cleaved mCherry or YFP (S4 Fig). The level of induction of ROS by flg22 was significantly higher with some effectors; however, there was no induction of ROS in the absence of flg22. Therefore, the results from this assay suggest that at least five effectors can suppress PTI and some may actually increase PTI responsiveness.

## Multiple candidate WY effectors trigger cell death in certain lettuce genotypes

The 21 cloned, untagged WY candidate effectors were screened for their ability to elicit cell death in 215 different accessions of wild and cultivated lettuce (S3 File) using Agrobacterium-mediated transient expression. These accessions collectively express the majority of the known *Dm* genes as well as new resistance factors [48]. Eleven of the 21 WY effectors (BLR13, BLR36, BLN06, BLN08, BSW03, BSW04m, BSW04p, BSW11, BSW13, BSW14, and BSW19) triggered chlorosis or necrosis in one or more accessions of lettuce (Fig 7). Alleles of the same effector showed similar patterns of necrosis and/or chlorosis except for BSW11 by UCDM14 and PI491178 (Fig 7). For BLN06, we found that the allele cloned from SF5, which lacks a signal peptide, caused necrosis in NunDm17 and chlorosis in RYZ2146 (S5 Fig), as was previously
<table><thead><tr><th>Genotype</th><th>species</th><th>BLR13</th><th>BLR36</th><th>BLN06</th><th>BLN08</th><th>BSW03</th><th>BSW04m</th><th>BSW04m/p</th><th>BSW04p</th><th>BSW11a</th><th>BSW11b</th><th>BSW13</th><th>BSW14</th><th>BSW19</th><th>GFP</th><th>+control</th></tr></thead><tbody><tr><td>CGN14312</td><td>indica</td><td>1.3</td><td>3.0</td><td>3.3</td><td></td><td></td><td></td><td></td><td></td><td>1.5</td><td>1.5</td><td>1.5</td><td>2.8</td><td>1.0</td><td>0.8</td><td>3.9</td></tr><tr><td>CGN14316</td><td>indica</td><td>0.8</td><td>1.5</td><td>0.0</td><td>0.5</td><td>0.4</td><td>1.5</td><td>2.3</td><td>1.0</td><td>1.0</td><td>1.8</td><td>1.6</td><td>0.6</td><td>0.0</td><td>0.4</td><td>3.8</td></tr><tr><td>25 acc.</td><td>saligna</td><td>0.0</td><td>0.4</td><td>0.0</td><td>3.2</td><td>0.2</td><td>0.2</td><td>0.3</td><td>0.2</td><td>0.0</td><td>0.0</td><td>0.0</td><td>3.6</td><td>0.1</td><td>0.1</td><td>3.9</td></tr><tr><td>11 acc.</td><td>saligna</td><td>0.0</td><td>2.9</td><td>0.0</td><td>3.6</td><td>0.4</td><td>0.3</td><td>0.5</td><td>0.3</td><td>0.0</td><td>0.0</td><td>0.1</td><td>4.0</td><td>0.0</td><td>0.2</td><td>4.0</td></tr><tr><td>CGN10888</td><td>saligna</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.8</td><td>0.0</td><td>0.0</td><td>0.0</td><td>1.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>4.0</td><td>0.0</td><td>0.0</td><td>4.0</td></tr><tr><td>CGN11341</td><td>saligna</td><td>0.0</td><td>0.0</td><td>0.0</td><td>1.6</td><td>0.5</td><td>0.3</td><td></td><td>0.4</td><td>0.0</td><td>0.0</td><td>0.0</td><td>3.8</td><td>0.0</td><td>0.0</td><td>3.8</td></tr><tr><td>CGN13326</td><td>saligna</td><td>0.0</td><td>1.7</td><td>0.0</td><td>3.7</td><td>0.0</td><td>0.3</td><td>0.5</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>4.0</td><td>0.0</td><td>0.2</td><td>4.0</td></tr><tr><td>CGN13327</td><td>saligna</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>4.0</td><td>0.0</td><td>0.0</td><td>4.0</td></tr><tr><td>CGN13330</td><td>saligna</td><td>0.0</td><td>0.3</td><td>0.3</td><td>1.3</td><td>0.5</td><td>0.6</td><td>1.0</td><td>0.3</td><td>0.2</td><td>0.0</td><td>0.4</td><td>4.0</td><td>0.0</td><td>0.3</td><td>4.0</td></tr><tr><td>CGN15716</td><td>saligna</td><td>0.0</td><td>0.0</td><td>0.0</td><td></td><td></td><td></td><td></td><td></td><td>0.0</td><td>0.0</td><td>0.0</td><td>4.0</td><td>0.0</td><td>0.0</td><td>4.0</td></tr><tr><td>CGN5282</td><td>saligna</td><td>0.0</td><td>0.0</td><td>0.0</td><td>2.5</td><td>0.0</td><td>0.0</td><td>1.5</td><td>0.0</td><td>0.0</td><td>0.0</td><td>4.0</td><td>4.0</td><td>0.0</td><td>0.0</td><td>4.0</td></tr><tr><td>CGN5309</td><td>saligna</td><td>0.0</td><td>0.8</td><td>2.5</td><td>3.7</td><td>0.3</td><td>0.6</td><td>0.3</td><td>1.0</td><td>0.0</td><td>0.0</td><td>4.0</td><td>4.0</td><td>0.0</td><td>0.3</td><td>4.0</td></tr><tr><td>CGN5314</td><td>saligna</td><td>0.0</td><td>0.8</td><td>3.3</td><td>3.8</td><td>0.8</td><td>0.8</td><td>0.8</td><td>1.5</td><td>0.0</td><td>0.0</td><td>0.0</td><td>4.0</td><td>0.0</td><td>0.4</td><td>4.0</td></tr><tr><td>CGN5315</td><td>saligna</td><td>0.0</td><td>0.8</td><td>1.8</td><td>3.9</td><td>0.3</td><td>1.0</td><td>0.8</td><td>1.9</td><td>0.0</td><td>0.0</td><td>0.0</td><td>4.0</td><td>0.0</td><td>0.0</td><td>4.0</td></tr><tr><td>CGN5318</td><td>saligna</td><td>0.0</td><td>1.3</td><td>0.0</td><td>3.5</td><td>0.8</td><td>0.3</td><td>2.0</td><td>0.3</td><td>0.0</td><td>0.0</td><td>0.0</td><td>4.0</td><td>0.0</td><td>0.1</td><td>3.7</td></tr><tr><td>CGN5322</td><td>saligna</td><td>0.0</td><td>0.0</td><td>0.0</td><td>3.5</td><td>0.2</td><td>0.1</td><td>1.0</td><td>2.4</td><td>0.0</td><td>0.0</td><td>0.0</td><td>4.0</td><td>0.0</td><td>0.1</td><td>4.0</td></tr><tr><td>CGN5329</td><td>saligna</td><td>0.0</td><td>0.0</td><td>0.0</td><td>2.4</td><td>0.1</td><td>0.9</td><td>1.0</td><td>1.8</td><td>0.0</td><td>0.0</td><td>0.0</td><td>2.8</td><td>2.5</td><td>0.3</td><td>3.9</td></tr><tr><td>CGN5330</td><td>saligna</td><td>0.0</td><td>0.0</td><td>0.0</td><td>3.3</td><td>0.0</td><td>0.6</td><td>1.0</td><td>3.1</td><td>0.0</td><td>0.0</td><td>0.0</td><td>3.5</td><td>0.0</td><td>0.1</td><td>4.0</td></tr><tr><td>LB108</td><td>saligna</td><td>0.0</td><td>0.0</td><td>0.0</td><td>1.8</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td></td><td>0.0</td><td>2.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>4.0</td></tr><tr><td>PI491000</td><td>saligna</td><td>0.0</td><td>1.0</td><td>0.0</td><td>4.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.3</td><td>0.0</td><td>1.5</td><td>1.6</td><td>0.0</td><td>0.0</td><td>0.0</td><td>3.9</td></tr><tr><td>PI491207</td><td>saligna</td><td>0.0</td><td>0.0</td><td>0.0</td><td>4.0</td><td>1.1</td><td>2.0</td><td>2.0</td><td>0.5</td><td>0.0</td><td>0.0</td><td>0.0</td><td>4.0</td><td>0.0</td><td>0.2</td><td>4.0</td></tr><tr><td>PI491208</td><td>saligna</td><td>0.0</td><td>0.0</td><td>0.0</td><td>3.6</td><td>0.6</td><td>2.0</td><td>3.6</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>4.0</td><td>0.0</td><td>0.7</td><td>4.0</td></tr><tr><td>Blonde Lente</td><td>sativa</td><td>0.5</td><td>0.5</td><td>0.0</td><td>0.0</td><td>0.5</td><td>0.5</td><td>1.0</td><td>0.0</td><td>2.0</td><td>2.7</td><td>0.5</td><td>0.0</td><td>0.0</td><td>0.2</td><td>4.0</td></tr><tr><td>Capitan</td><td>sativa</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.6</td><td>0.2</td><td>0.2</td><td>3.3</td><td>3.5</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.2</td><td>3.8</td></tr><tr><td>Capsule</td><td>sativa</td><td>0.0</td><td>0.0</td><td>0.0</td><td></td><td></td><td></td><td></td><td></td><td>0.0</td><td>0.0</td><td>0.0</td><td>4.0</td><td>0.0</td><td>0.0</td><td>4.0</td></tr><tr><td>Femke</td><td>sativa</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>2.4</td><td>3.8</td><td>3.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>4.0</td></tr><tr><td>Fenston</td><td>sativa</td><td>0.0</td><td>0.0</td><td>1.0</td><td>0.0</td><td>0.6</td><td>1.3</td><td>1.5</td><td>2.8</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.4</td><td>3.6</td></tr><tr><td>FrRsal-1</td><td>sativa</td><td>0.0</td><td>1.5</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.1</td><td>3.4</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>3.8</td></tr><tr><td>Ninja</td><td>sativa</td><td>0.0</td><td>1.2</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.2</td><td>2.0</td><td>3.4</td><td>0.0</td><td>0.6</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.2</td><td>4.0</td></tr><tr><td>Pennlake</td><td>sativa</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>2.1</td><td>0.0</td><td>0.4</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>3.5</td></tr><tr><td>PI491226</td><td>sativa</td><td>0.0</td><td>0.0</td><td>0.0</td><td>3.7</td><td>0.6</td><td>0.5</td><td>1.0</td><td>1.5</td><td>0.0</td><td>0.0</td><td>0.0</td><td>4.0</td><td>0.0</td><td>0.1</td><td>3.9</td></tr><tr><td>RYZ2164</td><td>sativa</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>3.3</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>4.0</td></tr><tr><td>Salvius</td><td>sativa</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>2.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>3.6</td></tr><tr><td>UCDM14</td><td>sativa</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.1</td><td>0.1</td><td>0.0</td><td>0.4</td><td>0.0</td><td>3.1</td><td>0.0</td><td>0.0</td><td>0.2</td><td>0.0</td><td>0.4</td><td>3.8</td></tr><tr><td>Versaï</td><td>sativa</td><td>0.0</td><td>0.0</td><td>0.4</td><td>0.0</td><td>0.0</td><td>0.7</td><td></td><td>0.0</td><td>3.2</td><td>2.3</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>4.0</td></tr><tr><td>ViAE</td><td>sativa</td><td>0.0</td><td>0.0</td><td>0.0</td><td>3.4</td><td>3.5</td><td>3.9</td><td>2.0</td><td>3.8</td><td>0.0</td><td>0.0</td><td>0.0</td><td>2.8</td><td>0.0</td><td>0.1</td><td>3.7</td></tr><tr><td>ViCQ</td><td>sativa</td><td>1.0</td><td>0.0</td><td>0.0</td><td>3.9</td><td>2.8</td><td>3.8</td><td>3.3</td><td>3.4</td><td>0.0</td><td>0.0</td><td>0.0</td><td>3.2</td><td>0.0</td><td>0.1</td><td>3.8</td></tr><tr><td>Arm999</td><td>serriola</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>3.6</td><td>4.0</td><td>1.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>4.0</td></tr><tr><td>CGN14280</td><td>serriola</td><td>2.2</td><td>0.0</td><td>0.5</td><td>0.3</td><td>0.0</td><td>0.0</td><td>0.3</td><td>0.4</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.3</td><td>3.3</td></tr><tr><td>ISR-380</td><td>serriola</td><td>0.5</td><td>0.5</td><td>1.0</td><td>0.0</td><td>1.0</td><td>0.5</td><td>0.5</td><td>0.5</td><td></td><td>0.3</td><td>1.0</td><td>2.5</td><td>0.4</td><td>4.0</td><td>4.0</td></tr><tr><td>LS102</td><td>serriola</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.3</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>1.0</td><td>2.3</td><td>1.3</td><td>0.2</td><td>4.0</td></tr><tr><td>PI491108</td><td>serriola</td><td>0.0</td><td>0.0</td><td>0.0</td><td>2.1</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>3.3</td><td>0.0</td><td>0.0</td><td>4.0</td></tr><tr><td>PI491178</td><td>serriola</td><td>0.0</td><td>0.0</td><td>0.0</td><td></td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>4.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>4.0</td></tr><tr><td>CGN14305</td><td>virosa</td><td>0.0</td><td>0.0</td><td>0.0</td><td>1.5</td><td>0.7</td><td>0.9</td><td>1.5</td><td>1.4</td><td>0.0</td><td>0.0</td><td>0.0</td><td>2.3</td><td>0.0</td><td>0.4</td><td>3.6</td></tr><tr><td>CGN4683</td><td>virosa</td><td></td><td></td><td></td><td>3.5</td><td>0.0</td><td>0.0</td><td></td><td>0.0</td><td></td><td></td><td></td><td></td><td>0.0</td><td>0.0</td><td>4.0</td></tr><tr><td>CGN5333</td><td>virosa</td><td>0.0</td><td>0.0</td><td>0.5</td><td>0.0</td><td>0.1</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>2.3</td><td>1.5</td><td>0.0</td><td>3.9</td></tr><tr><td>CGN9365</td><td>virosa</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.2</td><td>0.0</td><td>0.8</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>3.8</td><td>0.0</td><td>0.0</td><td>4.0</td></tr><tr><td>LS238</td><td>virosa</td><td>0.0</td><td>0.0</td><td></td><td>1.8</td><td>2.0</td><td>3.9</td><td>0.0</td><td>3.6</td><td>0.0</td><td>0.0</td><td>0.0</td><td>2.0</td><td>0.0</td><td>0.0</td><td>3.4</td></tr><tr><td>LS241</td><td>virosa</td><td>0.0</td><td>0.0</td><td>0.0</td><td>2.5</td><td>2.3</td><td>3.9</td><td>0.0</td><td>3.2</td><td>0.0</td><td>0.0</td><td>0.0</td><td>3.5</td><td>0.0</td><td>0.0</td><td>3.1</td></tr></tbody></table>

<table><thead><tr><th>Color Key</th><th>Score</th></tr></thead><tbody><tr><td>No reaction</td><td>0</td></tr><tr><td>Mild chlorosis</td><td>1</td></tr><tr><td>Chlorosis</td><td>2</td></tr><tr><td>Mild Necrosis</td><td>3</td></tr><tr><td>Necrosis</td><td>4</td></tr><tr><td>No data</td><td>-</td></tr></tbody></table>
**Fig 7.** Agroinfiltration results from the screen of diverse genotypes of lettuce for recognition of candidate WY containing effectors from *B. lactucae* indicated by necrosis and/or chlorosis. Candidate effectors were expressed in lettuce leaves using Agroinfiltration and leaves were scored for their reaction four to seven days post-infiltration. Qualitative scores were converted into numeric scores for data analysis: 0 = no reaction (green), 1 = mild chlorosis (light green), 2 = chlorosis (yellow), 3 = mild necrosis or mixture of chlorosis and necrosis (pink), and 4 = full necrosis (magenta). The figure shows average scores for each effector on each genotype; only accessions and effectors that had a necrotic or chlorotic response are shown. GFP was used for a negative control and either HopM1 or PsojNIP was used as a positive control; only leaves that showed an appropriate response to controls were scored. The two effectors containing RXLR motifs are highlighted in dark grey, while effectors with no RXLR motif are highlighted in blue. Scores for all accessions screened and sample sizes can be found at http://bremia.ucdavis.edu/BIL/BIL_intersection.php. Two groups of 28 and 11 accessions of *L. maligna* that had identical reactions are shown together; the individual accessions making up each group are shown in S3 File.

https://doi.org/10.1371/journal.ppat.1009012.g007

found for full length BLN06 from the isolate BL24 [45]. BLN06 from SF5 did not cause necrosis or chlorosis on LS102 despite being expressed (S5 Fig). BLN06 from BL24 causes chlorosis in LS102 [45]; however, without simultaneous comparisons with a clone of BLN06 from BL24, the reasons for this difference are unclear.

Two effectors, BLN08 and BSW14, triggered cell death on many accessions of *L. maligna*. Cell death was observed for additional effectors, up to four per genotype, on several accessions of *L. maligna* (Fig 7; S2 File). Two lines of *L. sativa* tested, ViAE and ViCQ, which both have introgressions from the wild lettuce species *L. virosa* [49], were observed to have cell death in response to five effectors (BLN08, BSW03, BSW04m, BSW04p, and BSW14) as well as the chimeric sequence BSW04m/p (Fig 8). BSW04m/p was cloned using the same primers as BSW04p but did not align to a single gene in the reference genome and was therefore determined to be a chimera between the paralogous sequences BSW04p and BSW04m. The *L. virosa* accessions that were the resistance donors for ViAE and ViCQ also were observed to have cell death in response to these five effectors, but not the chimeric sequence BSW04m/p (Fig 8). To assess whether effector triggered cell death is genetically linked to resistance to *B. lactucae* in ViAE, we tested 67 F<sub>8,9</sub> recombinant inbred lines (RILs) of Cobham Green x ViAE [49] using agroinfiltration of the five effectors and by screening these lines for resistance to *B. lactucae*. There were 43 resistant RILs, 19 susceptible RILs, and five segregating RILs. All five effectors triggered cell death in all the resistant RILs and none triggered cell death in the susceptible RILs (Fig 9), demonstrating genetic linkage between effector triggered cell death and resistance to *B. lactucae*. Effector response was also found to be segregating in the five RILs that segregated for resistance to *B. lactucae*.

Many *L. sativa* lines and a few genotypes of *L. maligna* and *L. serriola* were observed to have necrosis or yellowing in response to the chimeric effector BSW04m/p. Some of the genotypes that had effector triggered cell death caused by BSW04m/p also had cell death in response to BSW04p, but not the paralog BSW04m, which shares 91% amino acid identity with BSW04p. Three lines in the original screen, Capitan, Ninja, and Femke, share the resistance gene *Dm11*; therefore, we tested two additional lines, Fila and Mondian, which also contain *Dm11*, for recognition of BSW04p. These lines similarly had effector triggered cell death caused by BSW04m/p and BSW04p, but not BSW04m. Since BSW04p causes cell death on multiple varieties that contain *Dm11*, BSW04p is a candidate for the protein encoded by the *Avr11* gene.

## A region in the C-terminus of BSW04p is responsible for cell death on *L. sativa* cv. Ninja

To determine the region of the BSW04p sequence responsible for cell death on Ninja, we generated four domain-swapped chimeric proteins between BSW04m (which does not cause cell death on Ninja, despite being expressed (S6 Fig)) and BSW04p (which causes cell death on Ninja) at different locations along the protein (Fig 10C). Only the chimeras with the C-
<figure><img src="image_11.png" /><figcaption>Fig 8. Example agroinfiltration results of reactions of *L. sativa* ViAE and ViCQ and their respective progenitor R gene donors *L. virosa* LS241 and LS238 to candidate WY containing effectors from *B. lactucae*. Photos are representative of a typical leaf. GFP and HopM1 were used as negative and positive controls for necrosis, respectively. Leaves were collected five days post-infiltration and the first column for each accession shows the uncleared leaf tissue; the second column shows leaf tissue cleared in ethanol to make necrosis (brown areas) more visible.</figcaption></figure>

https://doi.org/10.1371/journal.ppat.1009012.g008

terminal domain of BSW04p elicited cell death on Ninja (Fig 10). Sequence comparisons of all the isoforms tested (both the domain swaps and alleles obtained through the original cloning experiment) narrowed down the causative region for cell death on Ninja between residues 405 and 458, which contains five amino acid polymorphisms between cell death triggering and non-cell death triggering sequences (Fig 10D). All of the other C-terminal polymorphisms outside of this region were found to be present in the unrecognized BSW04m alleles BSW04m-10C or BSW04m-4C. The causative region is located between two predicted WY domains (Fig 10C).

Discussion

Effectors play a critical role in interactions between pathogens and their hosts. Accurate prediction and annotation of effector repertoires is foundational for functional genomics studies of pathogens. Due to their economic importance in agriculture, an increasing number of
<figure><img src="image_12.png" /><figcaption>**Fig 9.** Agroinfiltration scores of the ViAE x Cobham Green RILs in response to the five effectors recognized by ViAE. Scores as in Fig 7. Candidate effectors were expressed in lettuce leaves using Agroinfiltration and leaves were scored for their reaction three to five days post-infiltration. GFP was used for a negative control and HopM1 used as a positive control for all leaves; only leaves showing the appropriate response to both were scored. Sample sizes for each group are 35 plants of Cobham Green, 45 plants of ViAE, 18 lines of susceptible RILs (88 plants), and 44 lines of resistant RILs (205 plants). Bar plot shows average scores with standard error bars.</figcaption></figure>

https://doi.org/10.1371/journal.ppat.1009012.g009

*Phytophthora* and downy mildew pathogens are being sequenced. Prior to this study, the majority of cloned avirulence genes have encoded proteins with a RXLR motif [50]. This was consistent with most bioinformatic studies using the RXLR motif as their primary, or only, criterion to identify candidate effectors. Our results reveal that there are many candidate effectors containing the effector-related WY domain that lack the canonical RXLR motif, especially in the downy mildews, but also in *Phytophthora* spp. In *B. lactucae*, we have shown that these non-RXLR, WY domain containing effectors show characteristics of RXLR effectors, such as N-terminal intrinsic disorder, lineage specificity, and expression during infection. Furthermore, some of these proteins can act as suppressors of the host immune system and some trigger cell death in certain lettuce lines. Therefore, despite lacking the canonical RXLR-motif, WY domain containing proteins can have both virulence and avirulence activities. Consequently, numerous candidate effector genes are likely to have been missed in the genomic
<figure><img src="image_13.png" /><figcaption>**Fig 10.** (A) Agroinfiltration scores for the BSW04 paralogs and domain swaps. One leaf per plant was infiltrated and scores were averaged across N replicates. Scoring and color scheme as in Fig 7. (B) Example agroinfiltration results from BSW04m-BSW04p domain swaps. Infiltrations were done in *L. sativa* cv. Ninja and scored 5 days post infiltration. (C) Sequence alignment of BSW04m, BSW04p, and the domain swaps. Black lines on sequence alignment indicate amino acid differences from BSW04m. Orange lines indicate recombination breakpoints used in generating the domain swaps. Blue squares show the predicted WY domains. Proteins with sequences that cause necrosis on Ninja are colored magenta; those that do not are colored aqua. (D) Sequence alignment of the putative recognized region of BSW04p for recognized and unrecognized clones of BSW04.</figcaption></figure>

https://doi.org/10.1371/journal.ppat.1009012.g010

analysis of other species within the Peronosporales. Bioinformatics pipelines for predicting effectors in both *Phytophthora* spp. and downy mildew pathogens should include an HMM for the WY domain because it has a low false positive rate (Table 1) and may discover more high confidence effectors than searches for the RXLR motif alone.

Some WY effectors were initially predicted to not have a signal peptide and therefore not to be secreted; this could have been due to misannotation or may reflect biological reality. Several

<table><caption>**Table 1.** False positive rates for HMM and string searches across 14 oomycete genomes.</caption><thead><tr><th>Motif</th><th>String/HMM</th><th>False Positive Rate (Average +SD)</th><th>Range</th></tr></thead><tbody><tr><td>RXLR</td><td>RXLR</td><td>45.2% ± 3.4%</td><td>32–54%</td></tr><tr><td>RXLR</td><td>[RQGH]XLR</td><td>55.2% ± 3.1%</td><td>44–61%</td></tr><tr><td>RXLR</td><td>Whisson HMM [15]</td><td>1.5% ± 0.8%</td><td>0.1–3.3%</td></tr><tr><td>EER</td><td>[DE][DE][KR]</td><td>49.4% ± 2.6%</td><td>41–62%</td></tr><tr><td>RXLR-EER</td><td>[RQGH]XLR + [DE][DE][KR]</td><td>15.4% ± 4.7%</td><td>7.0–23%</td></tr><tr><td>Degenerate RXLR</td><td>[RKHGQ][X]{0,1}[LMYFIV][RNK]</td><td>74% ± 0.8%</td><td>71–77%</td></tr><tr><td>WY</td><td>Boutemy HMM [28]</td><td>0.05% ± 0.09%</td><td>0–0.3%</td></tr></tbody></table>

https://doi.org/10.1371/journal.ppat.1009012.t001
appeared to be misannotated since they had a predicted signal peptide starting at the second start codon of the gene model. Therefore, studies that rely on the presence of a signal peptide when using predicted ORFs may be missing true effectors with misannotated start codons. Some WY effectors had clearly lost the signal peptide due to N-terminal deletion, such as BLN06 (relative to the isolate BL24 [45]) and BSW04m (relative to its paralog BSW04p). Hypothetically, signal peptide loss may be an evolutionary strategy for the pathogen to evade recognition of an effector while retaining the effector function encoded in the genome for future deployment.

Although only 10% of WY effectors from *B. lactucae* contained the canonical RXLR motif, nearly all N-terminal sequences contained a degenerate RXLR motif; whether or not this can function similarly to the canonical RXLR motif is yet to be determined. Functional studies such as those done for effectors from *Phytophthora* [51] would be informative to ascertain which of these degenerate motifs can function similarly to the canonical RXLR motif in protein secretion [16]. However, at present, the degenerate RXLR regular expression used in this paper should not be used on its own to search for genes encoding RXLR proteins due to the high false positive rate (>70%; Table 1). WY effectors lacking the non-canonical RXLR motif had features similar to proteins containing the canonical RXLR motif such as high N-terminal intrinsic disorder and the presence of an EER motif. N-terminal intrinsic disorder has been predicted for RXLR effectors of *Phytophthora* [46] and is also a common feature of bacterial effectors [52]. The biological significance of these features remains to be determined; however, it may be that these intrinsically disordered domains are important in post-translational modification or protein-protein interactions, as is the case for intrinsically disordered regions in other organisms [53,54].

The number of WY domains per protein varied considerably, from one to eight in *B. lactucae*. Many pathogen effectors exhibit rapid evolution and divergence due to the selective pressure of evolving host targets and host resistance proteins [26]. Duplication of domains may allow for evolution of novel effector functions or for evasion of host recognition while retaining function [31]. Variation in the number of repeated domains is reminiscent of NB-LRR proteins, which recognize effectors (or their activity) and also have leucine-rich repeat domains [55,56]. Duplication of regions encoding WY domains may have happened within a gene through replication errors or between different genes through recombination. The genomic sequences of the repeats can be analyzed in multiple isolates to reveal the origins of these duplications and domain expansion or contraction as has been done for NB-LRR proteins [57].

Analysis of the amino acid sequences of the WY domain in *B. lactucae* showed that it may be better considered a W[Y/F] domain, due to the equal preference for phenylalanine and tyrosine for the second characteristic residue. This substitution has a BLOSUM score of 3, indicating that it is fairly common. Both tyrosine and phenylalanine are aromatic amino acids; however, the hydroxyl group on tyrosine makes it slightly bulkier, more polar, and introduces a potential phosphorylation site. Variation in this region is not uncommon in other oomycete species: the fifth of the seven domains in Psr2 of *P. sojae* has an F instead of a Y and is similar in structure to the single domain of ATR1 in *H. arabidopsidis* that has a cysteine at the Y position [31]. Further structural characterization is needed to reveal whether these substitutions alter tertiary structure and biological function.

In oomycetes, WY domain containing effectors have been shown to have several functions including PTI suppression [58]. Here we show five WY effectors from *B. lactucae*, four of which lack the canonical RXLR motif, are able to suppress the host immune system by interfering with pathogen-triggered production of ROS. Suppression of host defenses is critical to the survival of *B. lactucae* and therefore it is not surprising that multiple effectors target the basal
immune system. Identification of the host targets of these effectors will determine which steps
in the signal transduction cascade are modulated by each effector or may reveal candidates for
susceptibility genes in the host that are required for successful proliferation of B. lactucae.

Effectors are powerful tools for the discovery and characterization of host resistance genes [59]. Eleven of the *B. lactucae* effectors tested caused necrosis or severe chlorosis on one or more lettuce lines. Necrosis and/or chlorosis is characteristic of recognition of effectors by R genes. This reaction could also be due to the non-specific toxic effects of effectors; however, the fact that necrosis only occurs in specific lettuce genotypes is consistent with a specific inter-action with R genes rather than a general toxicity. R genes that recognize effectors can be iden-tified by map-based cloning using segregating F<sub>2:3</sub> and RIL populations. Loci for recognition of five *B. lactucae* effectors (BLG01, BLG03, BLN08, BLR31, and BLR38) have been successfully mapped in other studies [24,45,60]. Two of these effectors, BLG01 and BLN08, have been shown to be recognized broadly by the non-host *L. saligna* [24,60]. Our study confirms the results for BLN08 and revealed BSW14 as an additional effector recognized by *L. saligna*. BLN08 and BSW14 share little sequence similarity though both contain WY domains. Non-host resistance in *L. saligna* is clearly complex [61], but these results indicate that it is mediated in part by recognition of multiple effectors.

Resistance to *B. lactucae* is present in some accessions of *L. virosa* [62,63]; however, its genetic basis and efficacy this has been little studied due to sexual incompatibility between *L. virosa* and other *Lactuca* spp. The lettuce lines ViAE and ViCQ contain introgressed disease resistance from *L. virosa*, which was achieved using embryo rescue followed by extensive back-crossing to *L. sativa* [49]. We have shown that ViAE and ViCQ recognize five WY domain containing effectors from *B. lactucae* and are therefore expected to be highly resistant to *B. lactucae* in the field; this is consistent with resistance of ViAE to all 162 isolates tested from California (https://bremia.ucdavis.edu/bremia_database.php). For ViAE, we show that these effector recognitions are genetically linked to each other, and importantly, completely genetically cosegregate with resistance to *B. lactucae*. Experiments are underway to determine the molecular basis of this resistance.

The non-RXLR, WY containing effector BSW04p is a candidate for Avr11 due to recognition by multiple lettuce lines containing *Dm11*. Confirmation that Avr11 encodes BSW04p requires assessing the nucleotide and expression level polymorphism of the effector in virulent and avirulent isolates of *B. lactucae*, as well as cloning the cognate R gene and testing whether it confers resistance to isolates expressing Avr11. Interestingly, the *Dm11* containing lines that recognize BSW04p do not recognize the paralog BSW04m, which elicits cell death in other lines of lettuce, such as ViAE and one accession of *L. serriola*. Using domain swap experiments, we narrowed down the region responsible for cell death in Ninja to the C-terminal half of BSW04p, in a region located in between two WY domains; because WY proteins often contain multiple equally spaced WY repeats, it is possible that this region also contains a WY domain that was not predicted by the HMM [31]. Structural characterization is needed to determine whether this region forms a WY-like fold and what structural consequences the amino acid differences have between BSW04p and BSW04m.

Not all RXLR candidate effectors have WY domains and the presence of a WY domain is not required for the avirulence phenotype of all effectors; for example, WY domains were not predicted for *B. lactucae* effectors BLG01 and BLG03, yet they elicit an immune response in lettuce [24]. Structural elucidation of an RXLR effector lacking the WY motif *H. arabidopsidis* ATR13 revealed that it contained a helical fold that was distinct from the WY fold [64]. It would be informative to identify and compare the protein structures of additional RXLR effectors that lack WY domains to determine whether there are other conserved C-terminal domains that may be involved in effector functions.
RXLR and WY effectors provide tools for monitoring pathogen populations and effector-
driven resistance breeding. Analysis of phenotypically diverse isolates from multiple locations
will allow the characterization of individual effector repertoires as well as the development of
the pan-repertoire for the pathogen species. Cloned effectors also will be highly instrumental
for identification of their cognate resistance genes as well as effector targets in the host. In
addition, screens for resistance using transient expression of individual effectors will allow the
pyramiding of resistance genes with different specificities that will maximize the evolutionary
hurdle for the pathogen to become virulent. Ultimately, knowledge of effector repertoires will
allow data-driven deployment of resistance genes leading to more durable disease resistance
[65].

**Materials and methods**

Effector prediction

**Signal peptide and WY domain prediction.** To search for the WY domain, an HMM was built using HMMer v3.1 [66] based on 721 amino acid sequences of WY domains in RXLR motif containing proteins from *P. infestans*, *P. sojae*, and *P. ramorum* obtained from Boutemy *et al.* [28]. The sequence described by the HMM spans the WY domain in the crystal structures of Avr3a11 and PexRD2 [27]. Candidate WY effectors were predicted using this HMM to search translated predicted ORF sequences (>80 amino acids) as well as gene models in the draft genome of *B. lactucae* isolate SF5 [44]. Sequences with a positive HMM bit score were considered to be putative WY domain effectors as in [28]. We also used an HMM built from predicted *B. lactucae* WY domains, but found that it did not identify any additional candidates compared to the *Phytophthora* based HMM. Signal peptide prediction was performed on candidate WY effectors using SignalP v 4.1 [67] and PhobiusSP [68]. Default settings were altered for SignalP v 4.1 to have sensitivity similar to SignalP v 3.0. Output was compared between SignalP v 4.1 sensitive and the combination of SignalP4.0 + SignalP 3.0, and identical results were obtained from the two methods. Gene models were more accurate for predicting signal peptides than translated ORFs—in part due to misannotated start codons upstream of the probable true start codon and signal peptides in the ORFs. However, on several occasions the gene model was missing a signal peptide found in the ORF; these gene models were manually updated to reflect this likely true start site. SignalP v 4.1 in sensitive mode was better able to predict signal peptides in proteins with misannotated start codons compared to SignalP v 4.0 or v 5.0.

**RXLR and EER prediction.** A combination of string searches and HMMs were used to search for the RXLR motif in oomycete predicted WY proteins. The following strings were used based on variants observed in downy mildews: [RQGH]XLR for RXLR and [DE][DE] [KR] for EER. An HMM from Whisson *et al.* [15] was also tested, although this did not reveal any additional RXLR motifs. To search for additional non-canonical RXLR motifs in WY candidates, a highly degenerate string of [RKHGQ][X]{0,1}[LMYFIV][RNK] was used. Searches were limited to the first 100 amino acids of the N-terminus of the predicted protein sequence.

**Estimation of false positive rates for effector prediction.** We determined false positive rates for each motif (Table 1) by analyzing multiple permutations of the non-redundant secretome using the same pipeline as described above for effector prediction. At least 10 random permutations of the sequence space were created using the MEME fasta-shuffle-letters program (with a kmer size of 1) with peptides starting after the cleavage site identified by SignalP v4.1. The false positive rate for each motif was estimated as the average frequency of detection in the permutated sequences divided by the observed frequency in the original sequences.
## Sequence comparison

To generate a neighbor-joining tree of the WY effector candidates, whole protein amino acid alignments were performed using MUSCLE 3.8.425 [69] implemented in Geneious 11.0.5 (http://www.geneious.com) with default settings. A UPGMA tree was built in Geneious with bootstrap resampling (100 replicates). The resulting tree was annotated using the Interactive Tree of Life [70] and stylistically refined in the GNU Image Manipulation Program (http://www.gimp.org). To build a sequence logo for the WY domain from *B. lactucae*, a multiple sequence alignment of WY domains was built using MUSCLE 3.8.425 [69], manually corrected for alignment errors, and the sequence logo generated using Weblogo (http://weblogo.berkeley.edu/logo.cgi).

## Prediction of intrinsic disorder

PONDR VSL2 [71] was used to calculate levels of intrinsic disorder for groups of candidate effectors, grouped by the presence of an RXLR-like motif, EER-like motif, and/or WY domain. PONDR VSL2 is a meta-predictor of protein disorder that utilizes neural networks and amino acid context to give a weighted disorder score at each amino acid position. The average sequence disorder was calculated at each amino acid position in each group of effectors and aligned on the first amino acid. Plots were generated from the average positional disorder scores calculated from all peptides in a group. The entire non-redundant, predicted secretome was used as a reference for comparison.

## Lineage specificity

To determine if the candidate WY proteins are unique to *B. lactucae*, BLASTp [72] was used to search for orthologs in other oomycete species. BLASTp-based sequence comparisons with an e-value threshold of 0.01 were performed against the following oomycete genomes: *Albugo laibachii* [8], *H. arabidopsidis* [73], *P. capsici* [74], *Phytophthora cinnamomi* [75], *P. infestans* [3], *Phytophthora parasitica* [76], *P. ramorum* [5], *P. sojae* [5], *Pythium ultimum* [7], *Saprolegnia parasitica* [77], *Pseudoperonospora cubensis* [22], *Plasmopara halstedii* [25], *Plasmopara viticola* JL-7-2 [78] and INRA-PV221 [79], *Peronospora tabacina* [80], and *Peronospora effusa* [81].

## Gateway cloning of effector candidates

Twenty-one WY effectors that differed in their number of WY domains were selected at random for cloning and functional analysis. Genes were cloned from amplified genomic DNA in an attempt to capture both alleles of each effector from the homokaryotic isolate SF5 [44]. In order to capture additional allelic diversity, some effectors were also cloned from the heterokaryotic isolate C82P24 [44]. Effector candidates without their signal peptide were cloned into the pEarleyGate100 (pEG100) vector for plant expression [82] using Gateway cloning (Thermo Fisher Scientific). Candidate effector genes were amplified from genomic DNA isolated from spores of *B. lactucae* isolate SF5 or C82P24 using Phusion High Fidelity Polymerase (Thermo Fisher Scientific) and primers (S4 File) that amplified each ORF after the predicted signal peptide cleavage site. The Kozak sequence (ACCATG) was added to the forward primer for correct translational initiation. PCR products were purified using precipitation with either polyethylene glycol (PEG) or homemade solid phase reversible immobilization (SPRI) beads to remove primers and primer-dimers, recombined into pDONR207 using BP Clonase II (Thermo Fisher Scientific), and transformed into chemically competent *Escherichia coli* DH10B cells. The resulting entry clones were sequenced using primers designed to
pDONR207 to confirm gene identity and identify alleles. Two alleles were obtained for many effectors; in addition, more than two alleles were obtained for several genes with the additional alleles representing chimeras between the two haplotypes of isolate SF5, which is diploid and homokaryotic. The chimeric sequences show distinct recombination breakpoints and may have arisen either during PCR or Gateway-based recombination. Clones of all unique sequences for each effector were retained because they could be informative for dissecting the sequence basis of host recognition. This resulted in multiple distinct clones (wildtype alleles plus chimeric sequences) for some effectors.

The entry clones were then recombined into pEG100 using LR Clonase II and transformed into *E. coli* to generate expression clones. The expression clones were then transformed into *Agrobacterium tumefaciens* strain C58rif+ using electroporation. Kanamycin-resistant colonies were confirmed for the transgene using gene-specific primers. For the domain-swap experiments, chimeras were generated using PCR-fusion/Gateway cloning using overlapping primers as previously described [83]. For localization experiments, effectors were cloned without their signal peptide either as described above into pEG104 [31], resulting in an N-terminal YFP fusion, or they were cloned into the Gateway compatible binary vector pl2B-KAN-LjUB1-GW-tHSP with mCherry added using overlapping primers [83]. The vector pl2B-KAN-LjUB1-GW-tHSP was constructed using Golden Gate assembly [84] with EC50507 T-DNA backbone vector (https://www.ensa.ac.uk/resources/) with a Gateway acceptor site added (including ccdb) and a pNOS-Kan-tNOS cassette added.

## Gene expression analysis

Expression of the candidate WY effectors in *B. lactucae* was evaluated using RNA-seq reads from a seven-day timecourse experiment of *B. lactucae* infecting lettuce seedlings (https://www.ncbi.nlm.nih.gov/bioproject/PRJNA523226). Reads were mapped to all gene models using SAMtools [85] and then the candidate WY gene models were extracted and counted using BEDtools multicov [86]. Raw read counts are given in S1 File. To confirm expression of all cloned candidate WY effectors, specific quantitative real-time PCR (qRT-PCR) primers were designed for each transcript for SYBR Green based quantification [87]. Five seedlings per replicate of *L. sativa* cv. Cobham Green infected with *B. lactucae* isolate SF5 were collected and frozen at 0, 2, 4, and 6 days post infection with uninfected seedlings also taken as a negative control. Total RNA was extracted from seedlings using the RNeasy Plant Kit (Qiagen). Samples were treated with DNase to remove genomic DNA and RNA was reverse transcribed into cDNA using RevertAid RT (Thermo Fisher Scientific). qRT-PCR was performed using PerfeCTa SYBR Green SuperMix (Quantabio) on a 7900HT Fast Real-Time PCR machine (Applied Biosystems).

## Agrobacterium-mediated transient assays

Agroinfiltration and transient expression experiments were performed on *N. benthamiana* and lettuce using conditions optimized for lettuce [88]. Cultures of *A. tumefaciens* strain C58Rif+ were started from glycerol stocks in Luria-Bertani (LB) media supplemented with 50 µg/µL each rifamycin and kanamycin and grown for 16–18 hours at 30°C. For infiltration, cultures were pelleted, washed, and resuspended in 10 mM MgCl₂ to OD = 0.3–0.5.

## Subcellular localization prediction and characterization

Subcellular localization was assessed for nine cloned candidate WY effectors using Agroinfiltration as described above. Three to five days post-Agroinfiltration, lettuce leaves expressing N-terminal YFP or mCherry fusions of effectors were examined for subcellular localization by
confocal microscopy. Nuclear staining was performed by incubating cut leaf tissue in 18 nM DAPI in water for at least five minutes. Imaging was performed on a Zeiss LSM 710 laser scanning confocal microscope with a 40x objective lens. Nuclear localization was predicted using NucPred [47] and nuclear localization motifs were predicted using LOCALIZER [89].

## Western blotting

Leaf tissue was collected from lettuce or *N. benthamiana* at 2–3 days post infiltration. Proteins were extracted from leaf tissue using GTEN buffer (10% glycerol, 25 mM Tris pH 7.5, 1 mM EDTA, 150 mM NaCl) supplemented with 2% PVPP, 10 mM DTT, 1x protease inhibitor cocktail, and 0.1% Tween-20 [90]. Total protein (20 µg) was run on a 4–12% Bis-Tris SDS-PAGE gel at 200V for 30 minutes using MOPS SDS running buffer. Proteins were transferred to a PVDF membrane using wet transfer and were probed using anti-GFP (Roche 1181446001) or anti-mCherry (BioVision 5993–100) antibodies. Secondary HRP-conjugated antibodies were used for visualization on x-ray film using ECL reagent (Pierce).

## PTI suppression assay

Effectors were expressed in five-week-old plants of *N. benthamiana* using Agroinfiltration as described above. Two days post-infiltration, two leaf discs (3.8 mm) were taken using a cork borer from each infiltration site away from leaf veins. Leaf discs were floated abaxial side up in 200 µL of distilled water in a 96-well white assay plate and incubated at room temperature for 24 hrs. To measure suppression of ROS production [91], the water was removed and 100 µL of assay solution was added to the leaf discs, which contained 17 µg/mL luminol (Sigma-Aldrich, St Louis, MO) and 10 µg/mL horseradish peroxidase type 6A (Sigma-Aldrich). One of the two paired leaf discs was exposed to 100 nM of flg22 peptide (flg+) and the other was not exposed to any elicitor as a control for endogenous ROS production. Three published bacterial suppressors of PTI from *P. syringae* pv tomato DC3000 [92] were used as positive controls for suppression of ROS production: HopS2 [93], HopX1, and HopT1 [94]; GFP and DsRed were used as negative controls. Luminescence was measured on a FilterMax F5 Plate Reader (Molecular Devices, San Jose, CA) promptly after adding in the assay solution and was measured every two minutes for a total of 40 minutes. For the untagged proteins, the assay was repeated on two separate occasions with a total of 15 replicates for the flg22+ treatment and eight replicates of the flg22- treatment. For the fluorescently-tagged proteins, the assay was repeated three times with a total of 16 replicates for both the flg22+ and the flg22- treatments. For quantification of total ROS burst, the area under the curve for luminescence was calculated from 2–40 minutes in R version 3.5.1 [95] using the ‘auc’ function with liner interpolation in the MESS package [96].

## Lettuce germplasm screen for effector recognition

Agroinfiltration was performed on lettuce as described above. The youngest fully expanded leaves of 3 to 4-week-old greenhouse grown lettuce plants were infiltrated with the *A. tumefaciens* cultures using a needleless syringe. Leaves were examined four to five days post-infiltration for signs of macroscopic cell death, indicative of immune recognition by host resistance proteins. *A. tumefaciens* containing pEG100:35S-GFP was used as a negative control; *A. tumefaciens* containing pBAV139:35S-HopM1 from *Pseudomonas syringae* [92] or pB7WG2:35S-PsojNIP from *P. sojae* [97] were used as positive controls. To control for false negatives, leaves that did not show cell death in response to the positive control were excluded from the analysis. To control for false positives, leaves that showed cell death in response to GFP control were also excluded from the analysis. In addition, any cell death observed in the
initial screening was confirmed by repeating infiltrations on at least four plants for effectors found to cause cell death.

## B. lactucae disease assays

Lettuce seedlings were screened for resistance or susceptibility to *B. lactucae* as described previously [49]. Briefly, cotyledons of week-old lettuce seedlings were sprayed with spores from *B. lactucae* isolate Bl: 23 at a concentration of 1–2 x 10⁵ spores/mL and scored for presence/absence of sporulation at 7 and 14 days post inoculation.

## Supporting information

S1 Table. RXLR-like sequences observed within 100 amino acids of the N-terminus of *B. lactucae* WY proteins, sorted by first amino acid. Sequences that occur in more than one protein are highlighted in red.
(DOCX)

S2 Table. Number of predicted secreted WY effectors with or without RXLR and/or EER motifs in the genomes of downy mildew pathogens and *Phytophthora* species.
(DOCX)

S1 Fig. Quantitative real-time PCR results for candidate effectors from *B. lactucae*. (A) Normalized expression levels for each effector transcript at 0, 2, 4, and 6 days post infection (dpi) and uninfected lettuce as a control. *B. lactucae* B-tubulin (BlBtub) was used as a house-keeping gene to quantify total amount of *B. lactucae* RNA on infected seedlings and *L. sativa* actin (LsAct) was used to quantify the amount of lettuce RNA. Transcript levels normalized to B-tubulin using the 2^ΔCt^ method (2(^Ct effector-Ct b-tubulin))^[98]. (B) Raw Ct values for each effector transcript at 0, 2, 4, and 6 days post infection (dpi) and uninfected lettuce as a control. *B. lactucae* B-tubulin (BlBtub) was used as a housekeeping gene to quantify total amount of *B. lactucae* RNA on infected seedlings and *L. sativa* actin (LsAct) was used to quantify the amount of lettuce RNA.
(TIFF)

S2 Fig. Intrinsic disorder of first 150 amino acids of RXLR and WY containing candidate effectors predicted from other oomycete genomes. The x-axis shows the amino acid position; the y-axis shows the intrinsic disorder (VSL2 PONDR) score. Proteins were categorized as WY with no RXLR (yellow), WYs with RXLR (black), RXLR+EER (with or without WY, purple), and the total predicted secretome (blue) using the same color scheme as Fig 2. RXLR and EER motifs were as described in Figs 2 and 3. Average positional disorder scores are shown for each class of proteins.
(TIFF)

S3 Fig. Stability of fluorescently tagged WY effectors transiently expressed in lettuce. (A) Anti-mCherry (mCh) Western blot. (B) Anti-YFP Western blot. (C) Expected sizes for each of the fusion proteins.
(TIFF)

S4 Fig. Stability of fluorescently tagged WY effectors transiently expressed in *Nicotiana benthamiana*. (A) Anti-mCherry (mCh) Western blot. (B) Anti-YFP Western blot. (C) Expected sizes for each of the fusion proteins.
(TIFF)
**S5 Fig. Agroinfiltration results from BLN06-SF5 and sequence comparison between BLN06, SF5, and BL24.** (A) Agroinfiltration scores from various isoforms of BLN06 from *B. lactucae* isolate SF5. Scores as in Fig 7. One leaf per plant was infiltrated and scores were averaged across 10 replicates. All leaves were confirmed to be expressing mCherry (mCh)-BLN06 by confocal microscopy. (B) Example confocal microscopy of N-terminal mCherry-tagged BLN06 in LS102, RYZ2164, and NunDm17 lettuce genotypes. (C) Sequence comparison between BLN06 isoforms from isolates SF5 and BL24 [45]. BLN06 was initially cloned from SF5 with the first 28 amino acids removed in case there was a cryptic signal peptide. For the mCherry-BLN06 fusions, these 28 amino acids were included in the protein sequence. Both alleles of BLN06 were cloned from SF5 (BLN06-SF5-1 and BLN06-SF5-2).
(TIFF)

**S6 Fig. BSW04m is expressed in Ninja but does not cause cell death.** (A) Agroinfiltration scores from *L. sativa* cv Ninja three days post infiltration. GFP is used as a negative control, HopM1 is used as a positive control for necrosis. (B) Representative photo of infiltrated Ninja leaves before clearing (left) and after clearing (right) with ethanol to visualize necrosis (brown areas). (C) Example confocal microscopy result of mCherry-BSW04m (red) transiently expressed in Ninja. All four leaves showed similar expression and nuclear localization. Chlorophyll autofluorescence is shown in blue. The fusion protein is intact as measured by Western blot (S3 Fig).
(TIFF)

**S1 File. Sequences and NCBI reference numbers of the 59 WY proteins predicted from the *B. lactucae* SF5 genome assembly.**
(XLSX)

**S2 File. Sequences of predicted translated ORFs from 14 oomycete genomes with results from RXLR and EER motif searches and WY HMM searches on each sequence.**
(XLSX)

**S3 File. Lettuce genotypes tested for recognition of *B. lactucae* effectors.**
(XLSX)

**S4 File. Primers used in Gateway Cloning of the *B. lactucae* effectors and for qRT-PCR.**
(XLSX)

**Acknowledgments**

We thank Keri Cavanaugh, Alyssa Schweickert, Catherine Lopez, Amber Robbins, Natalie Hamada, Jesus Cardenas, Pauline Sanders, and Sebastian Reyes Chin Wo (all UC Davis) for wet lab and computational assistance. We thank Anne Giesbers for the plasmid pL2B-KAN-LjUB1-GW-tHSP. We thank Gitta Coaker (UC Davis) for assistance with the ROS assay. We thank Ales Lebeda (Palacký University, Olomouc, Czech Republic) and Alex Beharav (The Hebrew University, Tel Aviv, Israel) for lines of *L. saligna*. We thank Lien Bertier, Anne Giesbers, and Elizabeth Georgian for helpful comments on earlier versions of this manuscript.

**Author Contributions**

**Conceptualization:** Kelsey J. Wood, Richard Michelmore.

**Formal analysis:** Kelsey J. Wood, Munir Nur, Juliana Gil, Kyle Fletcher, Kim Lakeman, Ayumi Gothberg, Tina Khuu, Jennifer Kopetzky, Archana Pandya.
**Funding acquisition:** Kelsey J. Wood, Richard Michelmore.

**Investigation:** Kelsey J. Wood, Munir Nur, Juliana Gil, Kim Lakeman, Dasan Gann, Ayumi Gothberg, Tina Khuu, Sanye Naqvi, Archana Pandya, Chi Zhang, Brigitte Maisonneuve.

**Supervision:** Kelsey J. Wood, Mathieu Pel, Richard Michelmore.

**Visualization:** Kelsey J. Wood.

**Writing – original draft:** Kelsey J. Wood, Richard Michelmore.

**Writing – review & editing:** Kelsey J. Wood, Richard Michelmore.

## References

1. Derevnina L, Petre B, Kellner R, Dagdas YF, Sarowar MN, Giannakopoulou A, et al. Emerging oomycete threats to plants and animals. Philos Trans R Soc B Biol Sci. 2016; 371. https://doi.org/10.1098/rstb.2015.0459 PMID: 28080985

2. Baldauf SL. The Deep Roots of Eukaryotes. 2003; 300:1703–1707. https://doi.org/10.1126/science.1085544 PMID: 12805537

3. Haas BJ, Kamoun S, Zody MC, Jiang RHY, Handsaker RE, Cano LM, et al. Genome sequence and analysis of the Irish potato famine pathogen *Phytophthora infestans*. Nature. 2009; 461:393–398. https://doi.org/10.1038/nature08358 PMID: 19741609

4. Rizzo DM, Garbelotto M, Hansen EM. *Phytophthora ramorum*: Integrative Research and Management of an Emerging Pathogen in California and Oregon Forests. Annu Rev Phytopathol. 2005; 43:309–335. https://doi.org/10.1146/annurev.phyto.42.040803.140418 PMID: 16078887

5. Tyler BM, Tripathy S, Zhang X, Dehal P, Jiang RHY, Aerts A, et al. *Phytophthora* genome sequences uncover evolutionary origins and mechanisms of pathogenesis. Science. 2006; 313:1261–1266. https://doi.org/10.1126/science.1128796 PMID: 16946064

6. Adhikari BN, Hamilton JP, Zerillo MM, Tisserat N, Lévesque CA, Buell CR. Comparative genomics reveals insight into virulence strategies of plant pathogenic oomycetes. PLoS One. 2013; 8:e75072. https://doi.org/10.1371/journal.pone.0075072 PMID: 24124466

7. Lévesque CA, Brouwer H, Cano L, Hamilton JP, Holt C, Huitema E, et al. Genome sequence of the necrotrophic plant pathogen *Pythium ultimum* reveals original pathogenicity mechanisms and effector repertoire. Genome Biol. 2010; 11:1–22. https://doi.org/10.1186/gb-2010-11-7-r73 PMID: 20626842

8. Kemen E, Gardiner A, Schultz-Larsen T, Kemen AC, Balmuth AL, Robert-Seilaniantz A, et al. Gene gain and loss during evolution of obligate parasitism in the white rust pathogen of *Arabidopsis thaliana*. PLoS Biol. 2011; 9:e1001094. https://doi.org/10.1371/journal.pbio.1001094 PMID: 21750662

9. Thines M, Choi Y. Evolution, Diversity, and Taxonomy of the Peronosporaceae, with Focus on the Genus *Peronospora*. Phytopathology. 2015; 106:6–18. https://doi.org/10.1094/PHYTO-05-15-0127-RVW PMID: 26649784

10. Phillips AJ, Anderson VL, Robertson EJ, Secombes CJ, van West P. New insights into animal pathogenic oomycetes. Trends Microbiol. 2007; 16:13–19. https://doi.org/10.1016/j.tim.2007.10.013 PMID: 18096392

11. Bozkurt TO, Schornack S, Banfield MJ, Kamoun S. Oomycetes, effectors, and all that jazz. Curr Opin Plant Biol. Elsevier Ltd; 2012; 15:483–92. https://doi.org/10.1016/j.pbi.2012.03.008 PMID: 22483402

12. Jones JDG, Dangl JL. The plant immune system. Nature. 2006; 444:323–329. https://doi.org/10.1038/nature05286 PMID: 17108957

13. Rehmany AP, Gordon A, Rose LE, Allen RL, Armstrong MR, Whisson SC, et al. Differential recognition of highly divergent downy mildew avirulence gene alleles by RPP1 resistance genes from two *Arabidopsis* lines. Plant Cell. 2005; 17:1839–1850. https://doi.org/10.1105/tpc.105.031807 PMID: 15894715

14. Birch PRJ, Armstrong M, Bos J, Boevink P, Gilroy EM, Taylor RM, et al. Towards understanding the virulence functions of RXLR effectors of the oomycete plant pathogen Phytophthora infestans. J Exp Bot. 2009; 60:1133–1140. https://doi.org/10.1093/jxb/ern353 PMID: 19204033

15. Whisson SC, Boevink PC, Moleleki L, Avrova AO, Morales JG, Gilroy EM, et al. A translocation signal for delivery of oomycete effector proteins into host plant cells. Nature. 2007; 450:115–118. https://doi.org/10.1038/nature06203 PMID: 17914356

16. Wawra S, Trusch F, Matena A, Apostolakis K, Linne U, Zhukov I, et al. The RxLR Motif of the Host Targeting Effector AVR3a of *Phytophthora infestans* Is Cleaved Before Secretion. Plant Cell. 2017; 29: tpc.00552.2016. https://doi.org/10.1105/tpc.16.00552 PMID: 28522546
17. Hiller NL, Hiller NL, Bhattacharjee S, Ooij C Van, Liolios K, Harrison T, et al. A Host-Targeting Signal in Virulence Proteins Reveals a Secretome in Malarial Infection. Science (80-). 2012; 1934:1934–1938.
https://doi.org/10.1126/science.1102737 PMID: 15591203

18. Coffey MJ, Sleebs BE, Uboldi AD, Garnham A, Franco M, Marino ND, et al. An aspartyl protease defines a novel pathway for export of *Toxoplasma* proteins into the host cell. Elife. 2015; 4:1–34. https://doi.org/10.7554/eLife.10809 PMID: 26576949

19. Boddey JA, O'Neill MT, Lopaticki S, Carvalho TG, Hodder AN, Nebl T, et al. Export of malaria proteins requires co-translational processing of the PEXEL motif independent of phosphatidylinositol-3-phosphate binding. Nat Commun. 2016; 7. https://doi.org/10.1038/ncomms10470 PMID: 26832821

20. Thines M. Bridging the gulf: *Phytophthora* and downy mildews are connected by rare grass parasites. Ausubel FM, editor. PLoS One. Public Library of Science; 2009; 4:e4790. https://doi.org/10.1371/journal.pone.0004790 PMID: 19274081

21. Bailey K, Cevik V, Holton N, Byrne-Richardson J, Sohn KH, Coates M, et al. Molecular cloning of ATR5 (Emoy2) from *Hyaloperonospora arabidopsidis*, an avirulence determinant that triggers RPP5-mediated defense in *Arabidopsis*. Mol Plant Microbe Interact. 2011; 24:827–838. https://doi.org/10.1094/MPMI-12-10-0278 PMID: 21361788

22. Tian M, Win J, Savory E, Burkhardt A, Held M, Brandizzi F, et al. 454 Genome sequencing of *Pseudo-peronospora cubensis* reveals effector proteins with a QXLR translocation motif. Mol Plant Microbe Interact. The American Phytopathological Society; 2011; 24:543–53. https://doi.org/10.1094/MPMI-08-10-0185 PMID: 21261462

23. Stassen J, Seidl M. Effector identification in the lettuce downy mildew *Bremia lactucae* by massively parallel transcriptome sequencing. Mol Plant Pathol. 2012; 13:719–31. https://doi.org/10.1111/j.1364-3703.2011.00780.x PMID: 22293108

24. Stassen JHM, Boer E den, Vergeer PWJ, Andel A, Ellendorff U, Pelgrom K, et al. Specific In Planta Recognition of Two GKLR Proteins of the Downy Mildew *Bremia lactucae* Revealed in a Large Effector Screen in Lettuce. Mol Plant Microbe Interact. The American Phytopathological Society; 2013; 26:1259–1270. https://doi.org/10.1094/MPMI-05-13-0142-R PMID: 23883357

25. Sharma R, Xia X, Cano LM, Evangelisti E, Kemen E, Judelson H, et al. Genome analyses of the sun-flower pathogen *Plasmopara halstedii* provide insights into effector evolution in downy mildews and Phytophthora. BMC Genomics. 2015; 16:741. https://doi.org/10.1186/s12864-015-1904-7 PMID: 26438312

26. Jiang RHY, Tripathy S, Govers F, Tyler BM. RXLR effector reservoir in two *Phytophthora* species is dominated by a single rapidly evolving superfamily with more than 700 members. Proc Natl Acad Sci U S A. 2008; 105:4874–4879. https://doi.org/10.1073/pnas.0709303105 PMID: 18344324

27. Win J, Krasileva K V, Kamoun S, Shirasu K, Staskawicz BJ, Banfield MJ. Sequence divergent RXLR effectors share a structural fold conserved across plant pathogenic oomycete species. PLoS Pathog. 2012; 8:e1002400. https://doi.org/10.1371/journal.ppat.1002400 PMID: 22253591

28. Boutemy L, King S, Win J, Hughes R, Clarke T, Blumenschein T, et al. Structures of *Phytophthora* RXLR effector proteins a conserved but adaptable fold underpins functional diversity. J Biol Chem. 2011; 286:35834–42. https://doi.org/10.1074/jbc.M111.262303 PMID: 21813644

29. Chou S, Krasileva K V, Holton JM, Steinbrenner AD, Alber T, Staskawicz BJ. *Hyaloperonospora arabidopsidis* ATR1 effector is a repeat protein with distributed recognition surfaces. Proc Natl Acad Sci U S A. 2011; 108:13323–13328. https://doi.org/10.1073/pnas.1109791108 PMID: 21788488

30. Yaeno T, Li H, Chaparro-Garcia A, Schornack S, Koshiba S, Watanabe S, et al. Phosphatidylinositol monophosphate-binding interface in the oomycete RXLR effector AVR3a is required for its stability in host cells to modulate plant immunity. Proc Natl Acad Sci. 2011; 108:14682–14687. https://doi.org/10.1073/pnas.1106002108 PMID: 21821794

31. He J, Ye W, Choi DS, Wu B, Zhai Y, Guo B, et al. Structural analysis of *Phytophthora* suppressor of RNA silencing 2 (PSR2) reveals a conserved modular fold contributing to virulence. Proc Natl Acad Sci. 2019; 2:201819481. https://doi.org/10.1073/pnas.1819481116 PMID: 30926664

32. Dou D, Kale SD, Wang X, Chen Y, Wang Q, Wang X, et al. Conserved C-terminal motifs required for avirulence and suppression of cell death by *Phytophthora sojae* effector Avr1b. Plant Cell. 2008; 20:1118–1133. https://doi.org/10.1105/tpc.107.057067 PMID: 18390593

33. Zhang P, Jia Y, Shi J, Chen C, Ye W, Wang Y, et al. The WY domain in the *Phytophthora* effector PSR1 is required for infection and RNA silencing suppression activity. New Phytol. 2019; 223:839–852. https://doi.org/10.1111/nph.15836 PMID: 30963588

34. Bos JIB, Armstrong MR, Gilroy EM, Boevink PC, Hein I, Taylor RM, et al. *Phytophthora infestans* effector AVR3a is essential for virulence and manipulates plant immunity by stabilizing host E3 ligase CMPG1. Proc Natl Acad Sci U S A. 2010; 107:9909–14. https://doi.org/10.1073/pnas.0914408107 PMID: 20457921
35. King SRF, McLellan H, Boevink PC, Armstrong MR, Bukharova T, Sukarta O, et al. *Phytophthora infestans* RXLR Effector PexRD2 Interacts with Host MAPKKK(varepsilon) to Suppress Plant Immune Signaling. Plant Cell. 2014; 26:1345-59. https://doi.org/10.1105/tpc.113.120055 PMID: 24632534

36. Maqbool A, Hughes RK, Dagdas YF, Tregidgo N, Zess E, Belhaj K, et al. Structural basis of host Autophagy-related protein 8 (ATG8) binding by the Irish potato famine pathogen effector protein PexRD54. J Biol Chem. 2016; 8:jbc.M116.744995. https://doi.org/10.1074/jbc.M116.744995 PMID: 27458016

37. Guo B, Wang H, Yang B, Jiang W, Jing M, Li H, et al. *Phytophthora sojae* effector PsAvh240 inhibits a host aspartic protease secretion to promote infection. Mol Plant. Chinese Society for Plant Biology; 2019; https://doi.org/10.1016/j.molp.2019.01.017 PMID: 30703565

38. Qiao Y, Liu L, Xiong Q, Flores C, Wong J, Shi J, et al. Oomycete pathogens encode RNA silencing suppressors. Nat Genet. Nature Publishing Group; 2013; 45:330-3. https://doi.org/10.1038/ng.2525 PMID: 23377181

39. Xiong Q, Ye W, Choi D, Wong J, Qiao Y, Tao K, et al. Phytophthora Suppressor of RNA Silencing 2 Is a Conserved RxLR Effector that Promotes Infection in Soybean and *Arabidopsis thaliana*. Mol plant-microbe Interact. 2014; 27:1379-1389. Available: http://apsjournals.apsnet.org/doi/abs/10.1094/MPMI-06-14-0190-R PMID: 25387135

40. Hou Y, Zhai Y, Feng L, Karimi HZ, Rutter BD, Zeng L, et al. A *Phytophthora* Effector Suppresses Trans-Kingdom RNAi to Promote Disease Susceptibility. Cell Host Microbe. 2019; 25:153-165.e5. https://doi.org/10.1016/j.chom.2018.11.007 PMID: 30595554

41. Qiao Y, Shi J, Zhai Y, Hou Y, Ma W. *Phytophthora* effector targets a novel component of small RNA pathway in plants to promote infection. Proc Natl Acad Sci. 2015; 112:5850-5855. https://doi.org/10.1073/pnas.1421475112 PMID: 25902521

42. Goritschnig S, Steinbrenner AD, Grunwald DJ, Staskawicz BJ. Structurally distinct *Arabidopsis thaliana* NLR immune receptors recognize tandem WY domains of an oomycete effector. New Phytol. 2016; 210:984-96. https://doi.org/10.1111/nph.13823 PMID: 26725254

43. Combier M, Evangelisti E, Piron M-C, Rengel D, Legrand L, Shenhav L, et al. A secreted WY-domain-containing protein present in European isolates of the oomycete *Plasmopara viticola* induces cell death in grapevine and tobacco species. PLoS One. 2019; 14:e0220184. https://doi.org/10.1371/journal.pone.0220184 PMID: 31356604

44. Fletcher K, Gil J, Bertier LD, Kenefick A, Wood KJ, Zhang L, et al. Genomic signatures of heterokaryosis in the oomycete pathogen *Bremia lactucae*. Nat Commun. 2019; 10:2645. https://doi.org/10.1038/s41467-019-10550-0 PMID: 31201315

45. Pelgrom AJE, Eikelhof J, Elberse J, Meisrimler CN, Raedts R, Klein J, et al. Recognition of lettuce downy mildew effector BLR38 in *Lactuca serriola* LS102 requires two unlinked loci. Molecular Plant Pathology. 2018;240-253. https://doi.org/10.1111/mpp.12751 PMID: 30251420

46. Shen D, Li Q, Sun P, Zhang M, Dou D. Intrinsic disorder is a common structural characteristic of RxLR effectors in oomycete pathogens. Fungal Biol. 2017; 121:911-919. https://doi.org/10.1016/j.funbio.2017.07.005 PMID: 29029698

47. Brameier M, Krings A, MacCallum RM. NucPred—Predicting nuclear localization of proteins. Bioinformatics. 2007; 23:1159-1160. https://doi.org/10.1093/bioinformatics/btm066 PMID: 17332022

48. Parra L, Maisonneuve B, Lebeda A, Schut J, Christopoulou M, Jeuken M, et al. Rationalization of genes for resistance to *Bremia lactucae* in lettuce. Euphytica. 2016; 210:309-326. https://doi.org/10.1007/s10681-016-1687-1

49. Maisonneuve B. *Lactuca virosa*, a source of disease resistance genes for lettuce breeding: results and difficulties for gene introgression. Eucarpia leafy Veg '03. 2003; 2003:61-67. Available: http://www.leafyvegetables.nl/download/12_061-067_maisonneuve.pdf

50. Anderson RG, Deb D, Fedkenheuer K, Mcdowell JM. Recent Progress in RXLR Effector Research. Mol Plant Microbe Interact. 2015; 28:1063-1072. https://doi.org/10.1094/MPMI-01-15-0022-CR PMID: 26125490

51. Wang S, Boevink PC, Welsh L, Zhang R, Whisson SC, Birch PRJ. Delivery of cytoplasmic and apoplast-tic effectors from *Phytophthora infestans* haustoria by distinct secretion pathways. New Phytol. 2017; 205-215. https://doi.org/10.1111/nph.14696 PMID: 28758684

52. Marin M, Uversky VN, Ott T. Intrinsic Disorder in Pathogen Effectors: Protein Flexibility as an Evolutionary Hallmark in a Molecular Arms Race. Plant Cell. 2013; 25:3153-3157. https://doi.org/10.1105/tpc.113.116319 PMID: 24038649

53. Uversky VN. Intrinsic disorder-based protein interactions and their modulators. Curr Pharm Des. 2013; 19:4191-4193. Available: http://www.ncbi.nlm.nih.gov/pubmed/23170892 https://doi.org/10.2174/1381612811319230005 PMID: 23170892
54. Iakoucheva LM, Brown CJ, Lawson JD, Obradović Z, Dunker AK. Intrinsic Disorder in Cell-signaling and Cancer-associated Proteins. J Mol Biol. 2002; 323:573-584. https://doi.org/10.1016/S0022-2836(02)00969-5 PMID: 12381310

55. Meyers BC, Kozik A, Griego A, Kuang H, Michelmore RW. Genome-Wide Analysis of NBS-LRR- Encoding Genes in *Arabidopsis*. 2003; 15:809-834. https://doi.org/10.1105/tpc.009308 PMID: 12671079

56. Hofberger JA, Jones JDG. A Novel Approach for Multi-Domain and Multi-Gene Family Identification Provides Insights into Evolutionary Dynamics of Disease Resistance Genes in Core Eudicot Plants. BMC Genomics. 2014; 31:1-35. https://doi.org/10.1186/1471-2164-15-966 PMID: 25380807

57. Kuang H, Woo S-S, Meyers BC, Nevo E, Michelmore RW. Multiple genetic processes result in hetero- geneous rates of evolution within the major cluster disease resistance genes in lettuce. Plant Cell. 2004; 16:2870-94. https://doi.org/10.1105/tpc.104.025502 PMID: 15494555

58. He Q, McLellan H, Hughes RK, Boevink PC, Armstrong M, Lu Y, et al. *Phytophthora infestans* effector SFI 3 targets potato UBK to suppress early immune transcriptional responses. New Phytol. John Wiley & Sons, Ltd (10.1111); 2018; nph.15635. https://doi.org/10.1111/nph.15635 PMID: 30536576

59. Vleeshouwers VGAA, Oliver RP. Effectors as tools in disease resistance breeding against biotrophic, hemibiotrophic, and necrotrophic plant pathogens. Mol Plant Microbe Interact. 2014; 27:196-206. https://doi.org/10.1094/MPMI-10-13-0313-IA PMID: 24405032

60. Giesbers AKJ, Pelgrom AJE, Visser RGF, Niks RE, Van den Ackerveken G, Jeuken MJW. Effector- mediated discovery of a novel resistance gene against *Bremia lactucae* in a nonhost lettuce species. New Phytol. 2017; https://doi.org/10.1111/nph.14741 PMID: 28833168

61. Giesbers AKJ, Boer E Den, Braspenning DNJ, Bouten TPH, Specken JW, van Kaauwen MPW, et al. Bidirectional backcrosses between wild and cultivated lettuce identify loci involved in nonhost resistance to downy mildew. Theor Appl Genet. 2018; 1-16. https://doi.org/10.1007/s00122-018-3112-8 PMID: 29802449

62. M. Bonnier FJ, Reinink K, Groenwold R. Genetic Analysis of Lactuca Accessions with New Major Gene Resistance to Lettuce Downy Mildew. Phytopathology. 1994. p. 462. https://doi.org/10.1094/Phyto-84- 462

63. Bonnier FJM, Reinink K, Groenwold R. New sources of major gene resistance in *Lactuca* to *Bremia lactucae*. Euphytica. 1992; 61:203-211. https://doi.org/10.1007/BF00039659

64. Leonelli L, Pelton J, Schoeffler A, Dahlbeck D, Berger J, Wemmer DE, et al. Structural elucidation and functional characterization of the *Hyaloperonospora arabidopsidis* effector protein ATR13. PLoS Pathog. 2011; 7:e1002428. https://doi.org/10.1371/journal.ppat.1002428 PMID: 22194684

65. Michelmore RW, Christopoulou M, Caldwell KS. Impacts of resistance gene genetics, function, and evolution on a durable future. Annu Rev Phytopathol. 2013; 51:291-319. https://doi.org/10.1146/annurev- phyto-082712-102334 PMID: 23682913

66. Eddy SR. Accelerated Profile HMM Searches. PLoS Computational Biology. 2011. p. e1002195. https://doi.org/10.1371/journal.pcbi.1002195 PMID: 22039361

67. Petersen TN, Brunak S, von Heijne G, Nielsen H. SignalP 4.0: discriminating signal peptides from trans-membrane regions. Nat Methods. Nature Publishing Group; 2011; 8:785-786. https://doi.org/10.1038/ nmeth.1701 PMID: 21959131

68. Käll L, Krogh A, Sonnhammer ELL. Advantages of combined transmembrane topology and signal peptide prediction-the Phobius web server. Nucleic Acids Res. 2007; 35:429-432. https://doi.org/10.1093/ nar/gkm256 PMID: 17483518

69. Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004; 32:1792-7. https://doi.org/10.1093/nar/gkh340 PMID: 15034147

70. Letunic I, Bork P. Interactive Tree Of Life (iTOL): an online tool for phylogenetic tree display and annotation. Bioinformatics. 2007; 23:127-8. https://doi.org/10.1093/bioinformatics/btl529 PMID: 17050570

71. Peng K, Radivojac P, Vucetic S, Dunker AK, Obradovic Z. Length-dependent prediction of protein in intrinsic disorder. BMC Bioinformatics. 2006; 7:1-17. https://doi.org/10.1186/1471-2105-7-1 PMID: 16393334

72. Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990; 215:403-410. https://doi.org/10.1016/S0022-2836(05)80360-2 PMID: 2231712

73. Baxter L, Tripathy S, Ishaque N, Boot N, Cabral A, Kemen E, et al. Signatures of adaptation to obligate biotrophy in the *Hyaloperonospora arabidopsidis* genome. Science (80-). 2010; 330:1549-1551. https://doi.org/10.1126/science.1195203 PMID: 21148394

74. Lamour KH, Mudge J, Gobena D, Hurtado-Gonzales OP, Schmutz J, Kuo A, et al. Genome sequencing and mapping reveal loss of heterozygosity as a mechanism for rapid adaptation in the vegetable
pathogen *Phytophthora capsici*. Mol Plant Microbe Interact. 2012; 25:1350–60. https://doi.org/10.1094/MPMI-02-12-0028-R PMID: 22712506

75. Institute JG. *Phytophthora cinnamomi* var. *cinnamomi* CBS 144.22 [Internet]. 2013. Available: https://gold.jgi.doe.gov/project?id=17783

76. Broad Institute. *Phytophthora parasitica* INRA-310 Genome sequencing and assembly. [Internet]. 2018 [cited 6 Jun 2018]. Available: https://data.nal.usda.gov/dataset/phytophthora-parasitica-inra-310-genome-sequencing-and-assembly

77. Jiang RHY, de Bruijn I, Haas BJ, Belmonte R, Löbach L, Christie J, et al. Distinctive Expansion of Potentia Virulence Genes in the Genome of the Oomycete Fish Pathogen *Saprolegnia parasitica*. PLoS Genet. 2013; 9:e1003272. https://doi.org/10.1371/journal.pgen.1003272 PMID: 23785293

78. Yin L, An Y, Qu J, Li X, Zhang Y, Dry I, et al. Genome sequence of *Plasmopara viticola* and insight into the pathogenic mechanism. Sci Rep. 2017; 7:1–12. https://doi.org/10.1038/s41598-016-0028-x PMID: 28127051

79. Dussert Y, Couture C, Mazet ID, Gouzy J, Piron M-C, Kuchly C, et al. Genome assembly and annotation of *Plasmopara viticola*, the grapevine downy mildew pathogen [Internet]. Portail Data Inra; 2019. https://doi.org/10.15454/4NYHD6

80. Derevnina L, Reyes-Chin Wo S, Martin F, Wood K, Froenicke L, Spring O, et al. Genome Sequence and Architecture of the Tobacco Downy Mildew Pathogen, *Peronospora tabacina*. Mol plant-microbe Interact. 2015; 28:1198–1215. https://doi.org/10.1094/MPMI-05-15-0112-R PMID: 26196322

81. Fletcher K, Klosterman SJ, Derevnina L, Martin F, Bertier LD, Koike S, et al. Comparative genomics of downy mildews reveals potential adaptations to biotrophy. BMC Genomics. BMC Genomics; 2018; 19:8–10. https://doi.org/10.1186/s12864-017-4386-y PMID: 29298683

82. Earley KW, Haag JR, Pontes O, Opper K, Juehne T, Song K, et al. Gateway-compatible vectors for plant functional genomics and proteomics. Plant J. 2006; 45:616–629. https://doi.org/10.1111/j.1365-313X.2005.02617.x PMID: 16441352

83. Atanassov II, Atanassov II, Etchells JP, Turner SR. A simple, flexible and efficient PCR-fusion/Gateway cloning procedure for gene fusion, site-directed mutagenesis, short sequence insertion and domain deletions and swaps. Plant Methods. 2009; 5. https://doi.org/10.1186/1746-4811-5-14 PMID: 19863796

84. Engler C, Marillonnet S. Golden Gate Cloning. In: Valla S, Lale R, editors. DNA Cloning and Assembly Methods. Totowa, NJ: Humana Press; 2014. pp. 119–131. https://doi.org/10.1007/978-1-62703-764-8_9 PMID: 24395361

85. Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al. The Sequence Alignment/Map format and SAMtools. Bioinformatics. 2009; 25:2078–2079. https://doi.org/10.1093/bioinformatics/btp352 PMID: 19505943

86. Quinlan AR, Hall IM. BEDTools: A flexible suite of utilities for comparing genomic features. Bioinformatics. 2010; 26:841–842. https://doi.org/10.1093/bioinformatics/btq033 PMID: 20110278

87. Thornton B, Basu C. Real-time PCR (qPCR) primer design using free online software. Biochem Mol Biol Educ. 2011; 39:145–154. https://doi.org/10.1002/bmb.20461 PMID: 21445907

88. Wroblewski T, Tomczak A, Michelmore R. Optimization of *Agrobacterium*-mediated transient assays of gene expression in lettuce, tomato and *Arabidopsis*. Plant Biotechnol J. 2005; 3:259–273. Available: http://www.ncbi.nlm.nih.gov/pubmed/17173625 https://doi.org/10.1111/j.1467-7652.2005.00123.x PMID: 17173625

89. Sperschneider J, Catanzariti AM, Deboer K, Petre B, Gardiner DM, Singh KB, et al. LOCALIZER: Subcellular localization prediction of both plant and effector proteins in the plant cell. Sci Rep. Nature Publishing Group; 2017; 7:1–14. https://doi.org/10.1038/s41598-016-0028-x PMID: 28127051

90. Win J, Kamoun S, Jones AME. Purification of effector-target protein complexes via transient expression in *Nicotiana benthamiana*. Methods Mol Biol. Humana Press; 2011; 712:181–94. https://doi.org/10.1007/978-1-61737-998-7_15 PMID: 21359809

91. Felix G, Duran JD, Volko S, Boller T. Plants have a sensitive perception system for the most conserved domain of bacterial flagellin. Plant J. 1999; 18:265–276. https://doi.org/10.1046/j.1365-313x.1999.00265.x PMID: 10377992

92. Wroblewski T, Caldwell KS, Piskurewicz U, Cavanaugh KA, Xu H, Kozik A, et al. Comparative large-scale analysis of interactions between several crop species and the effector repertoires from multiple pathovars of *Pseudomonas* and *Ralstonia*. Plant Physiol. 2009; 150:1733–49. https://doi.org/10.1104/pp.109.140251 PMID: 19571308

93. Jamir Y, Guo M, Oh HS, Petnicki-Ocwieja T, Chen S, Tang X, et al. Identification of *Pseudomonas syringae* type III effectors that can suppress programmed cell death in plants and yeast. Plant J. 2004; 37:554–565. https://doi.org/10.1046/j.1365-313x.2003.01982.x PMID: 14756767
94. Gimenez-Ibanez S, Hann DR, Chang JH, Segonzac C, Boller T, Rathjen JP. Differential suppression of Nicotiana benthamiana innate immune responses by transiently expressed *Pseudomonas syringae* type III effectors. Front Plant Sci. 2018; 9:1–19. https://doi.org/10.3389/fpls.2018.00001 PMID: 29410674

95. Team RC. R: A language and environment for statistical computing. [Internet]. Vienna, Austria.: R Foundation for Statistical Computing; 2018. Available: http://www.r-project.org/

96. Ekström CT. MESS: Miscellaneous Esoteric Statistical Scripts [Internet]. 2019. Available: https://cran.r-project.org/package=MESS

97. Cabral A, Oome S, Sander N, Küfner I, Nürnberger T, Van den Ackerveken G. Nontoxic Nep1-like proteins of the downy mildew pathogen *Hyaloperonospora arabidopsidis*: repression of necrosis-inducing activity by a surface-exposed region. Mol Plant-Microbe Interact. 2012; 25:697–708. https://doi.org/10.1094/MPMI-10-11-0269 PMID: 22235872

98. Pfaffl MW. Relative quantification. Real-time PCR. 2007; 64–82.