Downy Mildew effector HaRxL21 interacts with the transcriptional repressor TOPLESS to promote pathogen susceptibility

Hyaloperonospora arabidopsidis (Hpa) is an oomycete pathogen causing Arabidopsis downy mildew. Effector proteins secreted from the pathogen into the plant play key roles in promoting infection by suppressing plant immunity and manipulating the host to the pathogen’s advantage. One class of oomycete effectors share a conserved ‘RxLR’ motif critical for their translocation into the host cell. Here we characterize the interaction between an RxLR effector, HaRxL21 (RxL21), and the Arabidopsis transcriptional co-repressor Topless (TPL). We establish that RxL21 and TPL interact via an EAR motif at the C-terminus of the effector, mimicking the host plant mechanism for recruiting TPL to sites of transcriptional repression. We show that this motif, and hence interaction with TPL, is necessary for the virulence function of the effector. Furthermore, we provide evidence that RxL21 uses the interaction with TPL, and its close relative TPL-related 1, to repress plant immunity and enhance host susceptibility to both biotrophic and necrotrophic pathogens.

Introduction of TPL for its repressive activity in low auxin conditions [33,34]. Repression of jasmonate signalling relies on interaction between TPL and the EAR motif-containing Novel Interactor of JAZ (NINJA) [25,37]. Interaction with TPL and transcriptional repressor activity of NINJA were abolished when Leu residues in the EAR motif were mutated to Ala [25].
EAR motifs have now been identified in many plant proteins involved in development, stress and defence [38][39][40]. So far, a few pathogen effectors mimicking this EAR motif have been described although little is known about subsequent corepressor recruitment. For example, the Xanthomonas campestris type III effector XopD has been found to contain two tandemly repeated EAR motifs downstream of a DNA binding domain which are required for XopD-dependent virulence in tomato [41,42]. Effectors from the XopD superfamily that contain conserved EAR motifs have been found in Xanthomonas, Acidovorax and Pseudomonas species [41]. In addition, a conserved EAR motif is required for virulence in the Ralstonia solanacearum effector PopP2 although again no interaction with known Arabidopsis corepressors has been identified so far [43]. There are however, examples of pathogen effectors interacting with TPL family members including interaction of the Melampsora larici-populina effector MLP124017 with TPR4 [44].
The aim of our work was to determine how RxL21 is manipulating Arabidopsis to promote infection by Hpa. We demonstrate that RxL21 interacts in planta with the Arabidopsis corepressors TPL and TPR1 via an EAR motif and that this interaction is essential for RxL21 virulence activity against both Hpa and a necrotrophic plant pathogen. We find there is cooccurrence of RxL21 and TPR1 binding sites on promoter regions of a set of TPR1-repressed defence-related genes, suggesting that RxL21 virulence function involves perturbation of TPL/ TPR1 transcriptional repression during mobilization of host immunity.

RxL21 is conserved across multiple Hpa isolates and contains a C-terminal EAR motif
RxL21 is a 45 kDa effector protein identified from the genome of Hpa [16]. It contains a 'RLLR-DEER' motif at the N-terminus and an EAR (LxLxL) motif (amino acid sequence LMLTL) at the C-terminus. Alleles of RxL21 have been found in Hpa isolates Cala2, Emco5, Emoy2, Emwa1, Hind2, Maks9, Noks1 and Waco9 ([45] and BioProject PRJNA298674); alignment of their amino acid sequences is shown in S1 Fig. The signal peptide cleavage site was predicted to be between position 16 and 17 (SignalP-5.0). The RxLR-DEER motif is conserved across all alleles and with the exception of Noks1 (truncated due to Serine 197 changed to a Stop codon) the EAR motif at the C-terminus is also conserved between all aligned alleles. RxL21 is expressed in conidiospores of virulent (Waco9) and avirulent (Emoy2) Hpa isolates during Col-0 infection, as well as in Waco9 1, 3 and 5 days post inoculation [46] consistent with it playing a role in virulence.

RxL21 expression in planta causes enhanced susceptibility to both biotrophic and necrotrophic pathogens
Constitutive expression of RxL21 (from Hpa isolate Emoy2) in Arabidopsis has been shown to enhance growth of Hpa isolate Noco2 and a Pst isolate impaired in the suppression of early immune responses (Pst DC300 ΔavrPto/ΔavrPtoB) [16]. In addition, RxL21 increased bacterial growth when delivered into Arabidopsis via the type-three secretion system of Pst DC3000-LUX [16]. Here, Arabidopsis plants overexpressing RxL21 [16] were screened for susceptibility to Hpa isolates Noks1 and Maks9, and the necrotrophic pathogen Botrytis cinerea.
RxL21 lines were compared to both Col-0 wild type (WT) and lines expressing 35S::GUS (Col-0 GUS) which had been transformed and selected alongside the effector lines. Presence of the effector conferred enhanced susceptibility to both obligate biotroph Hpa isolates in two independent transgenic lines (RxL21a, b) compared to both Col-0 and Col-0 GUS, measured by total sporangiophores per seedling at 4 days post inoculation (dpi) (Fig 1A). RxL21 expression also resulted in increased lesion size caused by B. cinerea infection compared to controls ( Fig  1B). Hence, it appears that the RxL21 effector is targeting a mechanism essential for a full defence response against both biotrophic and necrotrophic pathogens.

RxL21 interacts with the CTLH domain of TPL via its EAR motif
It has been previously reported that RxL21 interacts with TPL using the matrix yeast two hybrid technique (Y2H) [17]. To confirm that this interaction could occur in planta, we first determined that RxL21 and TPL co-localise to the nucleus when transiently expressed in Nicotiana benthamiana leaves (Fig 2A). We next investigated the specific interaction motifs of both partners, which have not been previously reported. Two EAR-mutation constructs (RxL21ΔEAR1; Δ402-409 and RxL21ΔEAR2; Δ360-409) were used to determine whether the amino acids flanking the EAR motif are necessary for interaction with TPL ( Fig 2B). As expected, full-length RxL21 showed a positive interaction with TPL via activation of histidine (GAL1::HIS3) and adenine (GAL2::ADE2) reporter genes. The screen was performed in both directions with RxL21 and TPL fused to both the activation domain (AD) or the DNA binding domain (DB). Selective plates containing 3AT (a competitive inhibitor of the HIS3 gene product) were used as an additional control for increased stringency. Binding affinity appeared to be influenced by the direction of cloning with activation of the adenine reporter gene (and high stringency activation of the histidine reporter) only detected when TPL was fused to the DB domain. Deletion of either the initial N terminal sequence (RxL21ΔN) or the RxLR-DEER motif (RxL21ΔRxLR) did not prevent interaction with TPL. However, deletion of the RxL21 EAR motif abolished the interaction with TPL (Fig 2B), as did deletion of the CTLH domain   Interacts with the EAR domain of RxL21 in yeast. RxL21 was cloned without the signal peptide and lacking the EAR motif (RxL21ΔEAR1; Δ402-409 and HaRxL21ΔEAR2; Δ360-409), without the N-terminus (RxL21ΔN; Δ1-36) and RxLR-DEER motif (RxL21ΔRxLR; Δ1-51) and combinations of the above. Amino acid locations of deletion constructs are indicated. Deleting the EAR motif at the C-terminus of RxL21 abolishes interaction with TPL by Y2H. Successful mating is indicated by growth on media lacking Leucine and Tryptophan (-LW). Growth on media additionally lacking Histidine (-LWH) indicates GAL1::HIS3 reporter gene activation and protein-protein interaction. 3AT is used to increase stringency in a concentration dependent manner. Growth on media additionally lacking Adenine (-LWHA) indicates activation of the GAL2::ADE2 reporter gene. AD; activation domain. DB; DNA binding domain. (C) Deleting the CTLH domain of TPL (TPLΔCTLH; Δ25-91) abolishes interaction with RxL21 in yeast. Interaction (indicated by growth on-LWH +1mM 3AT) was observed between RxL21 and TPL but not between RxL21 and TPLΔCTLH. All Y2H was repeated at least twice with similar results. (D) TPL interacts with RxL21 and the interaction requires the EAR motif in planta. TPL:HA, TPLΔCTLH:HA together with RxL21 and RxL21 EAR motif mutants RxL21ΔEAR, RxL21-FMFTF and RxL21-AMATA (under control of an estradiol-inducible promoter and N-terminal Myc tag) were transiently expressed in N. benthamiana leaves. RxL21 expression was induced by 30μM β-estradiol 24 hr prior to harvesting. c-myc beads were used for immunoprecipitation (IP). HA antibody was used for TPL and TPLΔCTLH immunoblots (IB) and c-myc antibody was used to detect RxL21, RxL21ΔEAR and the respective EAR motif mutants. https://doi.org/10.1371/journal.ppat.1008835.g002

PLOS PATHOGENS
Pathogen effector HaRxL21 targets a host transcriptional repressor of TPL ( Fig 2C). Hence, both the EAR motif of RxL21 and the CTLH domain of TPL are necessary for direct protein-protein interaction. The mechanism of interaction between the effector and its host target (TPL) is mimicking the mechanism by which plant proteins (such as NINJA [25]) interact with TPL and recruit it to a complex for transcriptional repression.

Leucine residues in the EAR motif are critical for interaction with TPL
The three Leucine (Leu/L) residues (Leu389, Leu391, Leu393) in the EAR motif of RxL21 were mutated individually, in pairs or all together to Phenylalanine (Phe/F), Alanine (Ala/A) and Isoleucine (Ile/I). Substitution to Phe was chosen because previously Leu to Phe mutation of Leu354, Leu356 and Leu358 in the Arabidopsis repressor SUPERMAN prevented corepressor activity [47]. In addition, the three Leu residues in RxL21 were mutated to Ala (as it lacks functional side chains) [48] and Ile because it is the most similar amino acid to Leu [49].
Mutagenesis of all three Leu residues to Phe, Ala or Ile was found to abolish interaction of RxL21 with TPL by Y2H ( S2 Fig). Mutation of the 'x' residues within the EAR motif between the Leu residues (Methionine390 and Threonine392) had no effect on the interaction. Mutation of the first Leu residue (Leu389) or any pair of Leu residues in the EAR motif was sufficient to abolish RxL21-TPL interaction (S2A Fig). This analysis was repeated using TPL from N. benthamiana (NbTPL). RxL21 interacted with NbTPL and the same residues in the RxL21 EAR motif were necessary for interaction with NbTPL as for interaction with AtTPL (S2B Fig). These data are consistent with previous observations that mutation of any two Leucine residues within the DLELRL hexapeptide of SUPERMAN is sufficient to abolish repression of transcription [47].

RxL21 and TPL interact in planta
Bimolecular fluorescence complementation (Split YFP) was used to verify the interaction in planta. RxL21 was cloned into a vector with an N-terminal fused YFP fragment to minimise steric hindrance to the C-terminal EAR motif. Agrobacterium tumefaciens harbouring the N and C fragments of E-YFP fused to RxL21 and TPL were infiltrated into N. benthamiana and imaged using confocal microscopy. YFP fluorescence was detected in the nucleus, but no fluorescence was observed when RxL21ΔEAR was used instead of the full-length effector (S3 Fig).
The interaction between RxL21 and TPL was confirmed by co-immunoprecipitation (Co-IP) using transient expression in N. benthamiana. A Myc tag was fused to the N terminus of RxL21 (Myc-RxL21) and co-expressed with C-terminal HA-tagged TPL (TPL-HA). Immunoprecipitation of Myc-RxL21 using an anti-Myc antibody resulted in Co-IP of TPL-HA ( Fig  2D) confirming direct protein-protein interaction occurs in planta. Furthermore, Myc-tagged RxL21ΔEAR was unable to pull down TPL, and TPLΔCTLH was not immunoprecipitated by Myc-RxL21, demonstrating that, as in Y2H, the EAR motif and CTLH domain are required for RxL21-TPL interaction. Variants of RxL21 in which the Leu residues in the EAR motif were mutated to Phe or Ala (RxL21-FMFTF or RxL21-AMATA) were also unable to pull down TPL, confirming our observations in yeast that the Leu residues in the EAR motif of RxL21 are necessary for interaction with TPL.

An EAR motif is not sufficient to mediate interaction with TPL
In a previous large-scale study, no other Hpa effectors were found to interact with TPL [17]. The amino acid sequence of 134 RxLR (RxL), 491 RxL-like (RxLL) and 20 Crinkler (CRN) candidate effector proteins in the Hpa genome [15,45] were searched for the presence of the LxLxL EAR motif. The predicted signal peptide cleavage site in the effectors was then

PLOS PATHOGENS
Pathogen effector HaRxL21 targets a host transcriptional repressor identified using SignalP [50] and any effectors containing an LxLxL motif within the predicted signal peptide were subsequently excluded from further analysis.
The LxLxL motif was found to be present in 16 RxLs and 35 RxLLs (S1 Table). Of these, 2 RxLs and 4 RxLLs were found to contain multiple LxLxL motifs. The position of the LxLxL motif was noted, since previous reports of the EAR motif in functionally characterized transcriptional repressors is often at the N or C terminus [25,51,52]. The LxLxL motif of 5 RxLs (including RxL21) and 7 RxLLs was found to be within 35 amino acids of the C-terminus of the protein. No RxLs or RxLLs were found to have a LxLxL motif at the N terminus. We compared LxLxL-containing effectors against available subcellular localization data [53] and found, in addition to RxL21, HaRxL48, HaRxL55, HaRxLL100 and HaRxLL470 were nuclear localised in N. benthamiana [53], suggesting a possible functional role for the EAR motif in these effectors in manipulating host transcription. Ten of the LxLxL-containing Hpa effectors (including three of the four nuclear-localised ones) were tested using Y2H. None of the effectors showed an interaction with TPL (S4 Fig) indicating that the presence of an EAR motif is not sufficient to mediate an interaction with TPL and that RxL21 may have a unique function.

The EAR domain is necessary for RxL21 virulence function
To determine whether interaction with TPL is required for RxL21 virulence function, homozygous transgenic plants were generated expressing N-terminal HA tagged RxL21 or RxL21ΔEAR in a Col-4 background under the control of the 35S promoter (HA:21#1/2 or HA:21ΔEAR#1/2). Expression of the transgenes was verified by qPCR and western blot (S5 Fig). We noted that protein expression levels in all transgenic lines were very low perhaps due to a detrimental effect of effector accumulation. However, given the effector mRNA accumulation, we challenged these lines, compared to Col-4 WT and HA:GFP (35S promoter, Col-4 background) control plants, to Hpa infection. HA:21 lines exhibited enhanced susceptibility whereas HA:21ΔEAR lines did not differ in susceptibility from WT and GFP controls (Fig 3A). Similarly, compared to HA:GFP controls, both HA:21 lines showed significantly enhanced lesion area at 72 hpi after inoculation with B. cinerea whereas HA:21ΔEAR lines did not ( Fig  3B). Enhanced sporulation was also clearly seen on leaves from the HA:21#1 line compared to Col-4 WT, HA:GFP controls and HA:21ΔEAR#1 line (S6 Fig). These data demonstrate that the EAR motif of RxL21 (and hence most likely interaction with TPL/TPR1) is necessary for enhanced susceptibility to biotrophic and necrotrophic pathogens.
We also generated lines expressing RxL21 and RxL21ΔEAR under an estradiol-inducible promoter and with an N-terminal Myc tag in Col-0 background (Est:21#2/6/9 and Est:21ΔEAR#1/2/3). Expression of RxL21 after estradiol induction was verified by qPCR and protein accumulation (before and after induction) was assessed by western blot (S7 Fig). After estradiol induction, Est:21 lines showed enhanced susceptibility to Hpa ( Fig 3C) and B. cinerea infection (Fig 3D). In contrast, Est:21ΔEAR lines did not show enhanced susceptibility to either pathogen, confirming that the EAR motif is necessary for the virulence function of RxL21, and that RxL21 expression shortly before infection is sufficient to enhance host susceptibility.

TPL family members function in defence against multiple pathogens
It has previously been shown that TPL, TPR1 and TPR4 function in defence against Pst, with redundancy observed between TPL and TPR1 [20]. We screened tpr1-tpl-tpr4 triple knockout plants using Hpa isolate Noks1. The tpr1-tpl-tpr4 plants were found to be significantly more susceptible (Fig 4A) than both Col-0 WT and Col-0 GUS control plants. In addition, tpr1-tpl-tpr4 plants showed significantly increased lesion area after infection with B. cinerea (Fig 4B)

PLOS PATHOGENS
Pathogen effector HaRxL21 targets a host transcriptional repressor and visual symptoms of sporulation ( Fig 4C). These data suggest that transcriptional regulation by TPL and family members is a crucial part of immune signalling against multiple pathogens with different lifestyles.

PLOS PATHOGENS
Pathogen effector HaRxL21 targets a host transcriptional repressor

RxL21 also interacts with TPR1
TPL/TPR interactions can exhibit both redundancy (eg. all family members interact with multiple IAAs) and specificity (eg. several ERF TFs interact with 1 member or a subset of the TPL/ TPR family [23]). We tested interaction between all four TPRs and RxL21 using Y2H. In addition to TPL, RxL21 was found to interact with TPR1 (S8A Fig), the most closely related family member to TPL with 95% similarity at the amino acid level [20]. As with TPL, this interaction was abolished when the EAR motif of RxL21 was deleted. The RxL21-TPR1 interaction was also confirmed by Co-IP (S8B Fig

RxL21 interaction with TPL affects host transcription
Both the nuclear localisation and protein-protein interaction with TPL within the nucleus suggest that RxL21 plays a role in manipulating host transcription. We observed that when GFP:: RxL21 is transiently expressed in N. benthamiana, treating with DAPI before imaging results in co-localisation of GFP and DAPI (S9 Fig). DAPI has been shown to accumulate with high levels of heterochromatin [54] and therefore tightly bound chromatin in a repressed state. We therefore speculated that RxL21 co-localisation with DAPI is indicative of an effector mechanism involving transcriptional repression via interaction with TPL.

PLOS PATHOGENS
Pathogen effector HaRxL21 targets a host transcriptional repressor Consequently, we performed RNA-sequencing to compare gene expression in HA:21 and HA:21ΔEAR lines, before and after induction of PTI with flg22; a conserved stretch of 22 amino acids from bacterial flagellin that is recognised by the receptor FLS2 [55]. We compared two independent HA:21 lines with two HA:21ΔEAR lines, 2 h after both mock treatment and flg22 induction (full details of lines and treatments used for each sample is included in S2 Table). The total number of 75 bp paired-end reads generated was 676,059,072 across 24 samples. Transcript abundance was calculated by pseudoalignment of reads to AtRTD2 (Arabidopsis reference transcript dataset) [56] using Kallisto [57]. Genes with low expression across all samples were removed and TMM normalization was performed [58]. RNA count data after filtering and normalisation is provided in Sheet B in S2 Table. A PCA plot was used to visualise the data quality and variance between conditions, using the average read counts of biological replicates (S10 Fig); samples from independent transgenic lines (HA:21#1/2 or HA:21ΔEAR#1/2) were treated as biological replicates. 49% of variance between the samples can be explained by PC1 (flg22 treatment) however separation between RxL21 and RxL21ΔEAR lines is also observed. Differentially expressed genes (DEGs) were defined as having a log 2 fold change of 1 or greater, with an adjusted p-value cut-off (after false discovery rate correction) of 0.05. DEGs are listed in S3 Table; the total number of DEGs for each comparison is shown in Fig 5A.

RxL21 does not perturb the overall PTI response
The scale of transcriptional reprogramming in response to flg22 treatment is similar in both the HA:21 and HA:21ΔEAR lines, and indeed the vast majority of the flg22-induced DEGs are conserved across the two genotypes (2681 genes, with direction of differential expression also conserved: 591 downregulated and 2090 upregulated).  Table; Pearson correlation coefficients of 0.841 and 0.844 respectively). Hence the majority of the flg22 response is intact in HA:21 lines and RxL21-induced pathogen susceptibility is not due to large-scale interference with the PTI response.

RxL21 causes differential expression of defence-related genes
There were only 417 genes differentially expressed (DE) between RxL21 and RxL21ΔEAR transgenic lines under mock conditions, and 240 after flg22 treatment (Fig 5A; S3 Table). 113 genes were DE in both of these conditions, with the direction of differential expression conserved ( Fig 5B). We analysed the DEGs in two groups: those apparent under mock conditions (417 genes, 316 up and 101 down) and the subset of DEGs specific to flg22 treatment (127; 68 up and 59 down) (Fig 5B). In both the genes downregulated and upregulated by RxL21 compared to RxL21ΔEAR (under mock conditions) we found significant over-representation of Gene Ontology (GO) terms involved in response to biotic stimulus/defence and secondary metabolism. Enzymatic activity including oxidoreductase and hydrolase activity, and components of the endomembrane system were also enriched across the DEGs in both mock and flg22 treatments (S5 Table).
A significant proportion of the genes whose expression is perturbed by the RxL21 effector are associated with plant immunity (S6 Table). 82 genes differentially expressed in response to flg22 in WT plants [59] are mis-regulated by the effector, either before or after flg22 treatment. 213 genes differentially expressed in response to RxL21 are differentially expressed during Pst infection (DC3000 and/or DC3000HrpA compared to mock) [13], and 142 during B. cinerea infection [60].

PLOS PATHOGENS
Pathogen effector HaRxL21 targets a host transcriptional repressor We identified specific cases where RxL21 appeared to be suppressing the host defence response. Fifteen genes that are induced in response to flg22 treatment in WT Arabidopsis [59] were downregulated by the full-length effector compared to control after flg22 induction (including AVRPPHB SUSCEPTIBLE 3, CRK24, CRK38 and MYB85). The Pst data set profiles gene expression in both virulent Pst DC3000 (capable of secreting effector proteins directly For Est:21 and Est:21ΔEAR, data were obtained from 2 independent lines each with 3 biological replicates and expression was normalised to Arabidopsis tubulin gene. Black circles are individual data points and bars denote the mean ± SE of target gene expression. Letters indicate significant differences (P < 0.05) (One-way ANOVA with Tukey's honest significance difference). https://doi.org/10.1371/journal.ppat.1008835.g005

PLOS PATHOGENS
Pathogen effector HaRxL21 targets a host transcriptional repressor into plant cells via the Type III secretion system) as well as a strain DC3000HrpA-which lacks the Type III secretion system. Hence, we can distinguish between host response and pathogen manipulation of gene expression via effector proteins. 20 defence genes that are suppressed by Pst effectors were regulated in the same manner by RxL21 (i.e. both Pst effectors and RxL21 downregulated or upregulated expression compared to respective controls). Furthermore, 18 genes that are specifically expressed in response to Pst effectors (and not part of the host PTI response) were also differentially expressed in response to RxL21. These data indicate not only that RxL21 DEGs are involved in plant immunity, but also that shared mechanisms of host manipulation exist between RxL21 from Hpa and Pst effectors.
We quantified expression of 8 DEGs using qPCR in an independent set of lines; estradiolinducible lines Est:21, Est:21ΔEAR and Col-0 control. We selected eight of the genes downregulated by RxL21 compared to RxL21ΔEAR; these genes code for two TFs (WRKY63 and NAC019), three receptor-like proteins, two calmodulins and a protein of unknown function (Sheet B in S6 Table). Seven out of the eight genes were significantly downregulated in Est:21 compared to Est:21ΔEAR lines (Fig 5C). Six of these genes also showed reduced expression in Est:21 compared to wildtype Col-0. NAC019 did not differ between Col-0 and Est:21 but was reduced in Est:21 compared to Est:21ΔEAR, and WRKY63 showed an intermediate expression level in Est:21ΔEAR between Est:21 and Col-0. In general, the qPCR shows a similar pattern of log 2 fold change for the eight selected genes between RxL21 and RxL21ΔEAR lines compared to RNA-seq results (S11 Fig).

RxL21-repressed genes are enriched for TPR1 binding targets
We hypothesised that RxL21 modulates host gene expression via interaction with TPL/TPR1. RxL21 could recruit these corepressors to new sites on the genome (by either binding to the DNA itself or to TFs), or RxL21 might sustain TPL/TPR1 repression when it would normally be relieved during infection. Using iDNA-Prot [61] we found no evidence to suggest that RxL21 contains a DNA binding motif. Promoters of DEGs under mock conditions were significantly enriched for WRKY TF binding motifs in both up and down-regulated genes (S7 Table). In the DEGs only evident after flg22 treatment, we found significant enrichment for MYB TF binding motifs (as well as a WRKY binding motif) in the upregulated genes and CAMTA (Calmodulin-binding transcription activator) binding motifs in the genes downregulated by RxL21 (S7 Table). This suggests that RxL21 may exert at least some of its effects via modulating TPL/TPR interaction with WRKY, MYB and CAMTA TFs.
We also compared DEGs from our RNAseq analysis to chromatin immunoprecipitation (ChIP)-seq data revealing TPR1 target/bound genes from pTPR1:TPR1:GFP Col plants (Griebel et al. [62], Fig 6A). Only genes repressed by RxL21 compared to RxL21ΔEAR (both under mock and flg22 treatment) were enriched for TPR1 targets (Sheet A, S8 Table) with binding observed immediately upstream of the transcription start site (TSS) (Fig 6A, purple and orange  lines). Hence, for at least some of the RxL21 DEGs, RxL21 appears to impair expression by interaction with TPR1 upstream of the TSS. As expected, the level of TPR1 ChIP signal in RxL21-repressed DEGs was lower than in the group of all genes defined as TPR1 targets from the ChIP-seq data. The 20 genes repressed by RxL21 and identified as TPR1 targets encode several known positive regulators of immunity against biotrophic pathogens (PBS3, ICS1 and RLP20) and regulators of abiotic stress tolerance (WRKY46, WRKY63, NAC019) (S8 Table). The identification of WRKY46 and 63 as likely direct targets of RxL21 suppression via TPR1, and the enrichment of WRKY binding motifs in the promoters of RxL21 DEGs, suggests that these two TFs could be one key mechanism underlying RxL21-induced pathogen susceptibility.
Another key mechanism appears to be mis-regulation of salicylic acid (SA) signalling with

PLOS PATHOGENS
Pathogen effector HaRxL21 targets a host transcriptional repressor RxL21/TPR1 targets including genes for two enzymes required for SA accumulation (ICS1 and PBS3).
Genes that do not respond to the effector (RxL21 control in Fig 6A; Sheet B in S8 Table) do not show enrichment for TPR1 binding, while interestingly, genes induced by RxL21 (compared to RxL21ΔEAR in mock or after flg22 treatment; Fig 6A, yellow and brown lines) demonstrate depletion of TPR1 binding around the TSS. This suggests that promoters of genes upregulated by RxL21 are not TPR1 targets and mis-regulation of these by the effector is likely to be an indirect effect rather than through mis-placing bound TPR1.
To confirm the observed overlap of RxL21 repressed genes with TPR1 binding sites, we performed independent ChIP-qPCR in Est:21 and Est:21ΔEAR lines (Fig 6B) on seven genes (marked in bold in S8 Table). Promoter regions encompassing 500 bp upstream and 500 bp down-stream of the TSS for five of these genes (NAC019, AED1, STMP6, AT5G44574 and HR4) were enriched approximately 5-fold in RxL21 immunoprecipitated samples compared to RxL21ΔEAR samples. The promoter regions of CRK38 and WRKY63 didn't show strong ChIP-qPCR of RxL21-repressed genes. Two-week old seedlings, overexpressing RxL21 or RxL21ΔEAR with an Nterminal myc tag and under the control of an estradiol inducible promoter, were treated with 30 μM β-estradiol to induce RxL21 expression 18 hrs before harvesting and cross-linking with 1% formaldehyde. ChIP assays were performed with an anti-Myc antibody. In the ChIP-qPCR, the enrichment of immunoprecipitated DNA was normalized by the percent input method (signals obtained from ChIP samples were divided by signals obtained from an input sample). Error bar represents ± SE of four technical repeats. Arabidopsis Actin 2 was included as a control but no amplification was detected after 40 cycles. The experiment was repeated with similar results. https://doi.org/10.1371/journal.ppat.1008835.g006

PLOS PATHOGENS
Pathogen effector HaRxL21 targets a host transcriptional repressor enrichment, however we noted that CRK38 was only repressed by RxL21 after flg22 treatment in our data set. We also amplified AtActin2 as a negative control which was not detectable in ChIP samples after 40 qPCR cycles. Promoter regions of these five genes can, therefore, be directly bound by both RxL21 (and not RxL21ΔEAR) and TPR1. Furthermore, a lack of enrichment in the ChIP-qPCR assays in RxL21ΔEAR samples suggests that RxL21 binds to these promoter regions via TPR1 (maintaining repression when it would normally be relieved) rather than recruiting TPR1 to new locations.

RxL21 mimics a host gene regulation mechanism
Manipulation of host transcription by effector proteins is a key mechanism by which plant pathogens are able to render the plant a more hospitable environment for colonization. In this study we provide evidence that an effector from the oomycete pathogen Hpa mimics a mechanism by which plants regulate transcription throughout their development and in response to stress; the EAR motif is a repression motif that occurs in plant transcriptional repressors and mediates interaction with members of the TPL/TPR corepressor family. This mechanism of TPL recruitment to transcriptional repression complexes is highly conserved across the plant kingdom [26,34,36]. The Hpa effector RxL21 contains a C-terminal EAR motif which we have shown to be responsible for interaction with the CTLH domain of the transcriptional corepressor TPL in both in vitro yeast assays and in planta (Co-IP and Split YFP; Figs 2 and S3). This interaction is specific, and is abolished even when Leu residues in the EAR motif are substituted with the similar amino acid Ile (Fig 2 and S2 Fig). A pivotal role for TPL family members in regulation of immune responses has been previously implicated [20,63] and Arabidopsis plants lacking expression of TPL and its closest homologues, TPR1 and TPR4, show enhanced susceptibility to Pst, Hpa and B. cinerea (Zhu et al., 2010 and Fig 4). Interestingly, the effector RxL21 also interacts with TPR1 ( S8 Fig) and, when expressed in planta, enhances susceptibility of Arabidopsis to a biotrophic (Hpa), hemibiotrophic (Pst) and necrotrophic pathogen (B. cinerea) (Figs 1, 3 and S6). There are a few examples of EAR motifs in other pathogen effectors -the bacterial XopD effector family [41] and an effector from fungi that contains an EAR motif and interacts with TPR4 [44]. Crucially we show that the C-terminal EAR motif (and thus interaction with TPL/TPR1) is required for enhanced pathogen susceptibility provided by the RxL21 effector (Fig 3).

RxL21 alters expression of a subset of the host defence response
By performing RNA-seq comparing the transcriptional effects of expressing RxL21 compared to the effector lacking the EAR motif (RxL21ΔEAR) we were able to identify specific DEGs that result from RxL21-TPL/TPR1 interaction (Fig 5). Importantly, there is no large-scale change in gene expression in response to presence of the effector (either with or without flg22 treatment to induce PTI), yet effector expression still delivers a striking susceptibility enhancement to the host plant. This is similar to the situation with Pst where blocking of effector delivery through the Type III Secretion System only perturbed the expression of 872 PTI-regulated genes with the vast majority (5350) showing no change in expression or amplification of usual PTI-regulation [13]. These findings highlight the ability of pathogen effectors to target (and hence identify) specific components of the plant defence response with great effect.
The genes mis-regulated by RxL21 include several important components of plant defence responses. These include down regulation of several receptor-like protein (RLP)genes; RLP6, RLP22, RLP23 and RLP35, with RLP20 DE specifically after flg22 treatment. RLPs are the second largest family of Arabidopsis leucine rich repeat-containing receptors; several members

PLOS PATHOGENS
Pathogen effector HaRxL21 targets a host transcriptional repressor have been shown to function as PAMP receptors and it is becoming increasingly clear that RLPs play a critical role in pathogen recognition [64]. Consistent with this, many RLP genes (including RLP6, 22, 23 and 35) show increased expression during pathogen infection and/or hormone treatment [64,65]. Several RLPs are involved in fungal and oomycete resistance [66][67][68][69][70]. In addition, several other genes downregulated by RxL21 are known components of the defence response. These include genes that impact resistance against biotrophic pathogens such as accelerated cell death 6 (ACD6) and Cysteine-rich receptor-like kinase CRK13 [71] that positively regulate SA signalling [72], as well as Calmodulin-like CML41 that regulates flg22-induced stomatal closure [73]. Furthermore, miR825 target AT3G04220 [74] and the lipid transfer protein AZI3 [75] contribute to defence against necrotrophs and the lectin receptor kinase LecRK4.1 positively regulates Arabidopsis PTI and resistance against biotrophic and necrotrophic pathogens [76]. Reduction in expression of these genes is likely to contribute to the observed enhanced susceptibility of RxL21 over-expressing plants to both biotrophic and necrotrophic pathogens.
Twenty genes that are direct TPR1 targets are repressed by RxL21, seven of which were specific to flg22 treatment. This includes tetraspanin9 (TET9) which was repressed 2.6-fold by RxL21 compared to RxL21ΔEAR. TET8/9 promote the formation of exosome-like extracellular vesicles to deliver host sRNAs into B. cinerea to decrease fungal virulence [77]. TET9 accumulates around fungal infection sites after B. cinerea infection and tet9 loss-of-function mutants display weak but consistent enhanced susceptibility towards B. cinerea [77], leading us to speculate that TET9 repression could contribute to the B. cinerea susceptibility we see in RxL21-expressing lines.
Two WRKY TFs (WRKY46 and WRKY63) are direct TPR1 targets and repressed in RxL21 expressing plants in comparison to the RxL21ΔEAR variant. WRKY TFs have been shown to regulate defence responses against biotrophic and necrotrophic pathogens [78,79] and WRKY46 overexpressing plants are more resistant to Pst [80]. Moreover, non-host resistance against Erwinia amylovora is regulated by WRKY46 and WRKY54 via EDS1 [81] and WRKY46 is a transcriptional activator of PBS3 expression [82]. Consistent with this, PBS3 expression is also reduced in RxL21 compared to RxL21ΔEAR over-expressing plants.
Recently, it was shown that PBS3 binds to and protects EDS1 from proteasomal digestion [83], hence reductions in PBS3 expression could lead to lower levels of EDS1 protein and enhanced susceptibility to biotrophic pathogens [84]. Both WRKY46 and 63 bind to the promoter of the SA biosynthesis gene SID2 with increased levels of SA, correlated with increased levels of expression of these and several other WRKY TFs [85]. Consistent with this, SID2 (also called ICS1) is downregulated by RxL21 after flg22 treatment, however this may also be due to direct repression via association with TPR1 as SID2 is also a target of TPR1.
WRKY TF binding motifs were overrepresented in the promoters of RxL21 DEGs implying that many of the DEGs are downstream targets of WRKY TFs. It is possible that WRKY46 and WRKY63 are directly repressed by RxL21 which in turn leads to changes in expression of their target genes. We observed overrepresentation of CAMTA motifs in genes that were specifically downregulated by RxL21 in an EAR-dependent manner upon flg22 treatment, indicating that target genes of CAMTA TFs are showing differential expression. There is no evidence to date that CAMTA TFs are TPL/TPR targets (Causier, Ashworth, Guo, & Davies, 2012a) but CAMTA TFs play a role in immune regulation through suppressing pathogen-responsive genes and could therefore be important targets for RxL21 manipulation via TPL/TPR1 [86][87][88].

Recruitment of RxL21 to TPL/TPR1 transcriptional complexes
We have shown that the suppression of immunity by RxL21 depends on its EAR domain, and hence most likely through modifying actions of TPL and TPR1. Several scenarios are possible.

PLOS PATHOGENS
Pathogen effector HaRxL21 targets a host transcriptional repressor On the one hand, RxL21 could interact with a TF and mediate subsequent TPL/TPR1 binding (to novel sites) in a manner similar to NINJA or JAZ proteins [25]. Evidence for an interaction of RxL21 with the TF TCP14 was seen in Y2H experiments [17] and later it was shown that TCP14 shifts RxL21 into subnuclear foci [18]. However, if RxL21 was binding TFs independently of TPL/TPR1 then our ChIP-PCR experiments would be expected to show a similar enrichment in the RxL21 and RxL21ΔEAR at the tested loci. An alternative model is that RxL21 interferes with the repression (or lifting of repression) of existing TPL/TPR1 targets. We provide evidence to show that, at least for some genes, RxL21 appears to bind to TPL/ TPR1 within transcriptional complexes at plant gene promoters and prevent transcriptional de-repression (Fig 6). First, there is a significant over-representation in TPR1 binding sites upstream of genes that are repressed by full-length RxL21 compared to RxL21ΔEAR (where interaction with TPL/TPR1 is lost by deletion of the EAR motif). Second, promoters of several RxL21-repressed genes were shown to be binding targets of TPR1 and RxL21 (and not of RxL21ΔEAR). This indicates that in at least these cases, RxL21 likely binds to TPR1 protein already present at the gene promoter (as there is no binding of RxL21 to these promoters without the EAR motif) and subsequently exerts its activity on the TPL/TPR1 complex.
The TPL family of proteins are closely related to Groucho / Tup1 proteins in animals and fungi [89] which, like TPL, have been shown to form tetrameric structures. Tetramerization was suggested as a mechanism to enable recruitment of multiple TFs to a single complex [30] and binding of EAR motif peptides does not prevent tetramerization [30,90,91]. It is therefore plausible that RxL21 is able to bind TPL/TPR1 via its EAR motif while other epitopes of the TPL oligomer are binding other proteins such as TFs within the transcriptional complex. How TPL/TPR1-bound RxL21 behaves is not known. However, TPR1 activity was shown to be regulated by the (SUMO) E3 ligase SIZ1 [92], perhaps via SUMOylation inhibiting the corepressor activity of TPR1 by preventing its interaction with HISTONE DEACETYLASE 19 (HDA19). It is possible that RxL21 is shielding the SUMO attachment sites K282 and K721 in TPR1 and/or preventing (SUMO) E3 ligase activity, thereby enhancing the TPL/TPR corepressor activity.
There is no evidence that the only mechanism of RxL21 action is through maintained repression of direct TPL/TPR1 targets and it is worth noting that many DEGs were upregulated in the presence of RxL21. In addition to immune suppression, pathogen effectors are known to directly manipulate plant gene expression, and hence physiology, to aid infection [93,94]. During Pst infection, more than 1500 genes were specifically upregulated in response to Pst effectors, and we observed an overlap (15 genes) between these genes and those upregulated by RxL21. However, it is possible that upregulated genes are downstream targets of TFs or other regulators targeted by RxL21. This appears to be the most likely explanation given that the RxL21 upregulated genes are not enriched for TPR1 binding in wildtype plants and the enrichment for WRKY TF binding motifs in the promoter regions of these genes.
Remarkably, in planta expression of RxL21 alone can increase susceptibility to pathogens with varying lifestyle and virulence strategies. To our knowledge there are very few, if any, effectors that exhibit this activity. RxL21 is one of several effectors that mimic the EAR motif for transcriptional repression, and appears to actively initiate and/or maintain repression of gene expression mediated by TPL/TPR1. We show here that the RxL21 EAR motif is essential for its virulence function, and for modifying the expression of key host defence genes. Future interrogations will be to examine the RxL21 mode of action on TPL/TPR1 transcriptional repression complexes and determine which effector-manipulated defences ultimately result in enhanced susceptibility of the host plant.

PLOS PATHOGENS
Pathogen effector HaRxL21 targets a host transcriptional repressor

Yeast-2-Hybrid
Yeast-2-Hybrid (Y2H) screening was performed as described in Dreze et al. [96]. Briefly, yeast strains Y8800 and Y8930 harbouring AD-X and DB-X constructs respectively were mated on yeast-extract, peptone, dextrose, Adenine (YPDA) medium. Yeast was replica plated onto Synthetic complete (SC) media lacking combinations of Leu, Trp, His and Ade. 3AT (3-amino-1,2,4-triazole) was added to selective plates lacking His to increase stringency as a competitive inhibitor of the HIS3 gene product. Plates were cleaned using sterile velvets after 1 day and imaged after 4 days. To determine that constructs were not autoactivating, AD-X and BD-X constructs were mated with Y9830 and Y8800 harbouring empty pDB and pAD respectively and any autoactivating constructs were removed from further screening. To generate EAR motif mutants, site-directed mutagenesis (SDM) was performed using a QuikChange II Site-Directed Mutagenesis Kit (Agilent, Santa Clara, US). Primer sequences for SDM are included in S9 Table. Plant material 35S::HaRxL21 lines (RxL21a/b) are described in Fabro et al. [16]. To generate 35S::HA:: HaRxL21 lines in Arabidopsis, RxL21 and RxL21ΔEAR were cloned into pEarleygate201 [97] and transformed into Arabidopsis ecotype Col-4 using floral dipping {Clough:1998vw}. Independent transformants were selected on BASTA (Bayer CropScience, Wolfenbüttel, Germany) until homozygous. Western Blotting was performed on 14 day old seedlings to determine protein expression using anti-HA high affinity antibody (Roche, Penzberg, Germany). Primer sequences for RxL21 and RxL21 deletion variants are in S9 Table. To generate estradiol inducible lines, RxL21 and RxL21ΔEAR were cloned into the pER8 plasmid {Zuo:2000} and transformed into Arabidopsis (Col-0 background) via floral dip. Independent homozygous lines were selected on hygromycin (Invitrogen, Carlsbad, US).

Pathogen Assays
Hpa infection screens were performed as described in Tomé et al. [98]. Briefly, plants were sown in a randomized block design within the inner modules of p40 seed trays at a density of approximately 30 seedlings per module. Modules at the periphery of the tray were sown with WT. Plants were grown under a 12 h photoperiod, 20˚C, 60% humidity. Hpa isolates Noks1 and Maks9 were maintained on Col-0 plants and inoculated onto 2 week old seedlings at 3x10 4 spores/mL. After infection, plants were covered with a lid and sealed. Sporangiophore counts were performed using a dissecting microscope at 4 dpi. Botrytis infection screens were performed as described in Windram et al. [60]. Briefly, B. cinerea (strain pepper) was maintained on sterile tinned apricot halves in petri dishes, kept in the dark at 25˚C. Detached Arabidopsis leaves were placed on 0.8% agar and inoculated with a 10 μL droplet of spore suspension at 1x10 5 spores/mL in 50% sterile grape juice. Lesions were imaged at 24, 48, 72 and 96 h post infection (hpi). For Botrytis infections on estradiol-inducible lines, plants were grown at 10 h photoperiod, 16˚C at 60% humidity. Gene expression was induced by 30 μM β-

PLOS PATHOGENS
Pathogen effector HaRxL21 targets a host transcriptional repressor estradiol, 18 h before infection. Spores were washed from the surface of the B. cinerea culture plate using potato dextrose broth (PDB) medium and leaves were infected using a 5 μL droplet at a concentration of 1x10 4 spores/mL.

Localisation and Split YFP
For localisation, pK7WGF2 (N-terminal GFP fusion; [99]) and pB7WGR2 (N-terminal RFP fusion; [100]) vectors were used. BiFP (Bimolecular Fluorescence complementation in Planta) vectors were used for generation of C-or N-YFP fusion constructs [101]. Expression vectors were transformed into A. tumefaciens strain GV3101 and cultured overnight at 28˚C. A. tumefaciens harbouring expression constructs, along with P19 [102] were co-infiltrated into N. benthamiana. Leaves were imaged 48 h after infiltration by laser scanning confocal microscopy.

Co-immunoprecipitation
For Co-immunoprecipitation (Co-IP), TPL and RxL21 were expressed in pCsVMV-HA3-N-1300 (C-terminal HA tag) and Per8 (N-terminal myc tag) vectors. Proteins were expressed transiently in N. benthamiana. Leaves of about four-week-old plants were infiltrated with A. tumefaciens (OD600 = 0.2) in infiltration buffer (10 mM MgCl 2 , 10 mM MES, and 100 μM acetosyringone). After 24 hpi, leaves were sprayed with 30 μM estradiol to induce the expression of RxL21 and subsequent mutants. At 48 hpi, approximately 2 g of tissue from the infiltrated area was collected and frozen with liquid N 2 , then ground into powder using a mortar and pestle. All of the steps were carried out either on ice or in a 4˚C cold room. About three volumes of nuclei extraction buffer (NEB) (20 mM Hepes pH 8, 250mM Sucrose, 1mM Magnesium Chloride, 5mM Potassium Chloride, 40% (v/v) glycerol, 0.25% Triton X-100, 0.1 mM PMSF, Protease inhibitor cocktail (Roche cOmplete), were added to each sample. The resuspended samples were filtered using miracloth and centrifuged at 3,320 g for 15 min at 4˚C. The pellet was subsequently resuspended and washed with NEB in 2 mL microcentrifuge tubes. Washed pellet was resuspended in lysis buffer (10 mM Tris pH7.5, 0.5 mM EDTA, 1 mM PMSF, 1% Protease Inhibitor, 150 mM Sodium Chloride, 0.5% Igepal) and sonicated to break the nuclei. After sonication, the supernatant was centrifuged again to remove additional debris, and the supernatant was used as input for IP. Afterwards, the μMACS kit of magnetic microbeads, conjugated to an anti-c-myc monoclonal antibody (130-091-123; Miltenyi Biotech), was used for IP according to the manufacturer's instructions. Eluted proteins were analysed by western blot using HRP-conjugated anti-myc antibody (130-092-113; Miltenyi Biotec), or an anti-HA antibody (3F10; Roche).

RNA-Seq
Sample preparation. Arabidopsis seedlings were grown on ½ MS agar with sucrose for 14 days at 22˚C with photoperiod of 8 h light and 16 h dark, and 150 mmol photons � m -2 × s -1 . Seedlings were then transferred to 6 well plates containing 5 mL ½ MS liquid media (approximately 8 seedlings per well) and left to rest overnight. Media was replaced with 5 mL fresh ½ MS liquid media containing 100 nM flg22 or 5 mL ½ MS control. Each sample consisted of total pooled seedlings from a well, these were harvested 2 h post induction, briefly dried on tissue and flash frozen in liquid nitrogen. For each of 35S::HA::RxL21 and 35S::HA::RxL21ΔEAR, two independently transformed Arabidopsis lines were used (HA:21#1/2 and HA:21ΔEAR#1/ 2), each with 3 biological replicates. RNA was extracted from Arabidopsis tissue using Trizol (Invitrogen) and cleaned up using a RNeasy Mini Kit (Qiagen, Hilden, Germany) including on-column DNase treatment using a RNase-free DNase kit (Qiagen). Libraries were made using the NEBNext Ultra Directional RNA library prep kit. Each library was sequenced on two lanes of Illumina HiSeq 4000 generating 75 bp paired end reads. The sequencing files were examined using fastQC version 0.11.5.
Quantification of transcripts and DE gene analysis. Transcript abundances were generated using Kallisto version 0.44.0 [57] using an index generated using Arabidopsis Thaliana Reference Transcript Dataset 2 [56]. The number of bootstraps was set to 100. DEGs were generated using the 3D-RNA-seq app [103]. Sequencing replicates from the two HiSeq 4000 lanes were merged. Data was pre-processed by filtering to remove genes which did not meet the following criteria; 1) An expressed transcript must have at least 3 out of 23 samples with CPM (count per million read) �1 and 2) An expressed gene must have at least one expressed transcript. One sample was removed from further analysis (RxL21ΔEAR.flg22 line 2, biorep 3) due to outlying sequencing depth. RNA-seq read counts were normalised with Trimmed Mean of M-values method [58]. Models of expression contrasts (RxL21.flg22 vs RxL21.mock, RxL21ΔEAR.flg22 vs RxL21ΔEAR.mock, RxL21.mock vs RxL21ΔEAR.mock and RxL21.flg22 vs RxL21ΔEAR.flg22) were fitted using Limma Voom [104,105]. Genes were significantly DE in the contrast groups if they had BH adjusted p-value < 0.05 and log 2 -fold-change �1. We performed Go-Term analysis for biological process, cellular component and molecular function using AgriGO ( [106], http://bioinfo.cau.edu.cn/agriGO/). Significance of GO enrichment was determined with FDR < 0.05. Known Arabidopsis TF DNA-binding motifs were retrieved from CIS-DB version 1.02 [107], and those described in Franco-Zorilla et al., [108]. Promoter sequences defined as the 500 bp upstream of the transcription start site were retrieved from https://www.arabidopsis.org/ (Araport 11 annotation). Motif occurrences were determined using FIMO [109], and promoters defined as positive for a motif if it had at least one match with a P value < 10 −4 . Motif enrichment was assessed using the hypergeometric distribution against the background of all genes. P values < 0.001 were considered significant to allow for multiple testing. RNA seq data is available online at NCBI under Bioproject ID: PRJNA622757.

qPCR
Total RNA extraction was performed with TRIzol (Invitrogen) following the manufacturer's instructions. cDNA was synthesized by using the RevertAid cDNA synthesis kit (Thermofisher Scientific, Waltham, US). Quantitative PCR was performed in the Applied Biosystems 7500 FAST real-time PCR system using SYBR Green JumpStart Taq ReadyMix (Sigma-Aldrich, St Louis, US). Transcript levels of target genes were determined via the 2-ΔΔCt method [110], normalized to the amount of Arabidopsis tubulin 4 (AT5G44340) transcript. For expression analysis of HA:RxL21 and HA: 21ΔEAR lines, expression was normalised to Actin2 (AT3G18780) and UBQ5 (AT3G62250) transcript levels. Primer sequences are detailed in S9 Table.

Chromatin immunoprecipitation (ChIP) PCR
ChIP PCR was performed as described in Gendrel et al. [111]. Briefly, 2 week old seedlings were harvested and cross-linked in 1% formaldehyde (Sigma-Aldrich) solution under vacuum for 15 min. The isolated chromatin complex was resuspended in Nuclei lysis buffer 50 mM Tris-HCl, 10m M EDTA, % SDS and one tablet PI (Roche cOmplete) and sheared by sonication (Bioruptor, Diagenode, Ougrée, Belgium) to reduce the average DNA fragment size to around 500 bp. The sonicated chromatin complex was diluted in ChIP dilution buffer (16.7 mM Tris-HCl, 1.2 mM EDTA, 167 mM NaCl, 1.1% Triton X-100 and one tablet PI) and immunoprecipitated by anti-myc antibody (ab9132; ChIP-grade, Abcam, Cambridge, UK).

PLOS PATHOGENS
Pathogen effector HaRxL21 targets a host transcriptional repressor After reverse cross-linking, the immunoprecipitated DNA was extracted by using equal amounts of phenol/chloroform/isoamyl alcohol and precipitated with 100% EtOH, 1/10 volume of 3 M Sodium acetate and, 10 mg/mL glycogen. DNA was resuspended in Milli-Q water and analyzed by qPCR with gene specific primers. Primer sequences are shown in S9 Table. Input % in IP samples was calculated by (100 � 2^(Adjusted input-IP)). AtActin2 (AT3G18780) was included as a control.
Supporting information S1 Data. Excel spreadsheet containing, in separate sheets, the underlying numerical data and statistical analysis for Figures/Figure panels 1, 3, 4, 5c, 6b Table. Over-represented GO categories in differentially expressed genes between RxL21 and RxL21ΔEAR in both mock and flg22 conditions. Singular enrichment analysis of GO terms was performed using AgriGO. The cutoff P-value after false discovery rate correction is  Table. A) Overlap between genes associated with TPR1 binding sites and genes differentially regulated in Arabidopsis plants expressing the RxL21 effector. "RxL21 induced" or "RxL21 repressed" indicates genes which show significantly higher or lower expression respectively in HA:21-expressing plants compared to HA:21ΔEAR-expressing plants. Number of genes for each comparison is shown with overlap to TPR1 targets in brackets. P values indicate significance of overlap between each group of genes with TPR1 targets. Genes selected for verification by ChIP-PCR are indicated in bold (Column D). Gene descriptions taken from Araport 11. B) The control set of 150 genes which show no differential expression between HA:21 and HA:21ΔEAR. (XLSX) S9

PLOS PATHOGENS
Pathogen effector HaRxL21 targets a host transcriptional repressor