Protein engineering expands the effector recognition profile of a rice NLR immune receptor
Abstract
Plant nucleotide binding, leucine-rich repeat (NLR) receptors detect pathogen effectors and initiate an immune response. Since their discovery, NLRs have been the focus of protein engineering to improve disease resistance. However, this approach has proven challenging, in part due to their narrow response specificity. Previously, we revealed the structural basis of pathogen recognition by the integrated heavy metal associated (HMA) domain of the rice NLR Pikp (Maqbool et al., 2015). Here, we used structure-guided engineering to expand the response profile of Pikp to variants of the rice blast pathogen effector AVR-Pik. A mutation located within an effector-binding interface of the integrated Pikp–HMA domain increased the binding affinity for AVR-Pik variants in vitro and in vivo. This translates to an expanded cell-death response to AVR-Pik variants previously unrecognized by Pikp in planta. The structures of the engineered Pikp–HMA in complex with AVR-Pik variants revealed the mechanism of expanded recognition. These results provide a proof-of-concept that protein engineering can improve the utility of plant NLR receptors where direct interaction between effectors and NLRs is established, particularly where this interaction occurs via integrated domains.
https://doi.org/10.7554/eLife.47713.001Introduction
Protein engineering offers opportunities to develop new or improved molecular recognition capabilities that have applications in basic research, health and agricultural settings. Protein resurfacing, where the properties of solvent-exposed regions are changed (often by mutation), has been used extensively in diverse areas from antibody engineering for clinical use to the production of more stable, soluble proteins for biotechnology applications (Chapman and McNaughton, 2016).
Intracellular nucleotide binding, leucine-rich repeat (NLR) receptors are key components of plant innate immunity pathways. They recognize the presence or activity of virulence-associated, host-translocated pathogen effector proteins and initiate an immune response (Kourelis and van der Hoorn, 2018; Cesari, 2018). Because they confer resistance to disease, plant NLRs are widely used in crop breeding programs (Dangl et al., 2013). However, the recognition spectra of plant NLRs tend to be very specific, and pathogens may delete detected effectors from their genome or evolve novel effector variants that are not detected by the NLRs to re-establish disease (Yoshida et al., 2016).
The potential of engineering NLRs to overcome these limitations, or to detect new effector activities, is emerging (Rodriguez-Moreno et al., 2017). Gain-of-function random mutagenesis has achieved some success in expanding the activation sensitivity or effector recognition profiles of NLRs (Segretin et al., 2014; Giannakopoulou et al., 2015; Harris et al., 2013). An alternative strategy, in which NLRs perceive protease effectors through their activity on engineered host proteins, can lead to expanded recognition profiles (Carter et al., 2018; Helm et al., 2019; Kim et al., 2016). When detailed knowledge of the direct binding interfaces between an effector and an NLR are known, there is potential for protein resurfacing to modify interactions and thereby impact immune signaling.
Plant NLRs are modular proteins that are defined by their nucleotide-binding (NB-ARC) and leucine-rich repeat (LRR) domains, but that also have either an N-terminal coiled-coil (CC) or a Toll/Interleukin-1/Resistance-protein (TIR) signaling domain (Takken and Goverse, 2012). Many NLRs also contain non-canonical integrated domains (Kroj et al., 2016; Sarris et al., 2016; Bailey et al., 2018). Integrated domains are thought to be derived from ancestral virulence-associated effector targets that directly bind pathogen effectors (or host proteins; Fujisaki et al., 2017) or are modified by them, thereby initiating an immune response (Cesari et al., 2014; Le Roux et al., 2015; Maqbool et al., 2015; Sarris et al., 2015; De la Concepcion et al., 2018). These domains present an exciting target for protein engineering approaches to improve NLR activities. NLRs containing integrated domains (often called the ‘sensor’) typically function in pairs, requiring a second genetically linked NLR (the ‘helper’) for immune signaling (Białas et al., 2018; Eitas and Dangl, 2010).
Two rice NLR pairs, Pik and Pia, contain an integrated Heavy Metal Associated (HMA) domain in their sensor NLR that directly binds effectors from the rice blast pathogen Magnaporthe oryzae (also known as Pycularia oryzae) (Maqbool et al., 2015; Cesari et al., 2013; Ortiz et al., 2017; Guo et al., 2018). The integrated HMA domain in the sensor NLR Pik-1 directly binds the effector AVR-Pik. Co-evolutionary dynamics has driven the emergence of polymorphic Pik-1 HMA domains and AVR-Pik effectors in natural populations, and these display differential disease-resistance phenotypes (Costanzo and Jia, 2010; Kanzaki et al., 2012; Yoshida et al., 2009). The Pikp NLR allele only responds to the effector variant AVR-PikD, but the Pikm allele responds to AVR-PikD, AVR-PikE, and AVR-PikA. These phenotypes can be recapitulated in the model plant Nicotiana benthamiana using a cell-death assay, and are underpinned by differences in effector–receptor binding interfaces that lead to different affinities in vitro (Maqbool et al., 2015; De la Concepcion et al., 2018).
We hypothesized that by combining naturally occurring favorable interactions observed across different interfaces, as defined in different Pik-HMA/AVR-Pik structures (Maqbool et al., 2015; De la Concepcion et al., 2018), we could generate a Pik NLR with improved recognition profiles. Here, we graft an interface from Pikm onto Pikp by mutating two residues in Pikp (Asn261Lys, Lys262Glu), forming PikpNK-KE. This single-site mutation strengthens the cell-death response in N. benthamiana to AVR-PikD, and gains a Pikm-like response to AVR-PikE and AVR-PikA. We show that this gain-of-function phenotype correlates with increased binding affinity of the effectors by the PikpNK-KE–HMA domain in vitro and in vivo, and demonstrate that this mutation results in a Pikm-like structure for PikpNK-KE when in complex with AVR-Pik effectors. Finally, we confirm that the newly engineered interface is responsible for the expanded response of PikpNK-KE by mutation of the effectors.
This study serves as a proof-of-concept for the use of protein resurfacing by targeted mutation to develop plant NLR immune receptors with new capabilities. In the future, such approaches have the potential to improve disease resistance in crops.
Results
Structure-informed engineering expands Pikp-mediated effector recognition in N. benthamiana
By comparing protein interfaces in the structures of Pikp-HMA and Pikm-HMA bound to different AVR-Pik effectors (Maqbool et al., 2015; De la Concepcion et al., 2018), we hypothesized that we could engineer expanded effector recognition capabilities by point mutation of Pikp. We constructed a series of mutations in the previously identified interface 2 and interface 3 regions of Pik-HMA–AVR-Pik structures (De la Concepcion et al., 2018), swapping residues found in Pikm into Pikp (Figure 1A, Figure 1—figure supplement 1). We then screened these mutations for expanded effector recognition by monitoring cell death in a well-established N. benthamiana assay (Maqbool et al., 2015; De la Concepcion et al., 2018). We found that one double mutation in two adjacent amino-acid residues contained within interface 3, Asn261Lys and Lys262Glu (henceforth PikpNK-KE), induced cell death in response to AVR-PikE and AVR-PikA (Figure 1, Figure 1—figure supplement 1A,B). Mutations that cause loss of response to AVR-PikD are most likely to be compromised in effector binding. Western blots confirmed that all of the proteins were expressed in plants (Figure 1—figure supplement 1C).

Structure-informed engineering expands Pikp-mediated effector recognition in N. benthamiana.
(A) Sequence alignment of the Pikp-1 and Pikm-1 HMA domains. Secondary structure features of the HMA fold are shown above, and the residues that are located at binding interfaces are as colored. Key residues from interface 2 and interface 3 involved in this study are highlighted in red. (B) Representative leaf images showing Pikp- (left) or Pikp-1NK-KE (right)-mediated cell death in response to AVR-Pik variants as autofluorescence under UV light. (C) Autofluorescence intensity is scored as previously described (Maqbool et al., 2015; De la Concepcion et al., 2018). Cell death assay scores are represented as dot plots for Pikp and PikpNK-KE (blue and purple, respectively). For each sample, all of the data points are represented as dots with a distinct color for each of the three biological replicates; these dots are plotted around the cell death score for visualization purposes. The size of the centre dot at each cell death value is directly proportional to the number of replicates in the sample with that score. The total number of repeats was 80. Data for Pikp have been previously shown (De la Concepcion et al., 2018), but was acquired at the same time as those for PikpNK-KE. The estimation methods used to visualize differences in the data sets are shown in Figure 1—figure supplement 3. (D) Conductivity measurements showing ion leakage as a quantitative measure of cell death. The centre line represents the median, the box limits are the upper and lower quartiles, the whiskers extend to the largest value within Q1 – 1.5x the interquartile range (IQR) and the smallest value within Q3 + 1.5x IQR. All the data points are shown as dots with distinct colors for each biological replicate. For each experiment, six biological replicates with 5 or 10 internal repeats were performed (total data points = 40). ‘p’ is the p-value obtained from statistical analysis and Tukey’s HSD (honestly significant difference) test.
-
Figure 1—source data 1
Cell death scoring data used in the preparation of Figure 1C.
- https://doi.org/10.7554/eLife.47713.009
-
Figure 1—source data 2
Conductivity measurements used in the preparation of Figure 1D and Figure 1—figure supplement 2B.
- https://doi.org/10.7554/eLife.47713.010
We subsequently focused on this double mutant, and independently repeated the cell-death assay to ensure its robustness (Figure 1B,C, Figure 1—figure supplement 3). When looking at the intensity of cell death mediated by PikpNK-KE, as for the Pikm allele (De la Concepcion et al., 2018), we observed an effector hierarchy of AVR-PikD > AVR-PikE > AVR-PikA . PikpNK-KE shows a comparable, but elevated (not statistically significant in the case of AVR-PikE), response to effector variants when compared to that of Pikm (Figure 1—figure supplements 2 and 4). The cell-death response to AVR-PikD mediated by PikpNK-KE is similar but elevated when compared to that mediated by Pikp (Figure 1B,C). PikpNK-KE does not show a response to the stealthy AVR-PikC variant. To obtain a quantitative measure of cell death, we performed ion-leakage assays (Figure 1D, Figure 1—figure supplement 2, Supplementary file 1). The results of these assays correlate well with the cell-death index scores based on autofluorescence described above.
We conclude that the single Asn261Lys/Lys262Glu (PikpNK-KE) mutation at interface 3 in the Pikp NLR expands this protein’s recognition profile to include the effector variants AVR-PikE and AVR-PikA, similar to that observed for Pikm.
The engineered PikpNK-KE-HMA mutant shows increased binding to effector variants in vivo and in vitro
We used yeast-2-hybrid (Y2H) and surface plasmon resonance (SPR) to determine whether the expanded PikpNK-KE cell-death response in N. benthamiana correlates with increased binding affinity of the PikpNK-KE-HMA domain for AVR-Pik effectors.
As AVR-PikE and AVR-PikA showed some interaction with Pikp-HMA using these approaches (Maqbool et al., 2015; De la Concepcion et al., 2018), we tested interactions with PikpNK-KE-HMA side-by-side with interactions with wild-type Pikp-HMA. Using Y2H, we observed increases in growth and in blue coloration (both indicative of protein–protein interaction) for Pikp-HMANK-KE with effectors AVR-PikE and, particularly, AVR-PikA when compared with Pikp-HMA (Figure 2A). This was accentuated with more stringent conditions (imposed by increasing concentration of Aureobasidin A). Unexpectedly, we also observed limited yeast growth for Pikp-HMANK-KE with AVR-PikC at the lower stringency (Figure 2A). The unrelated M. oryzae effector AVR-Pii was used as a negative control. Expression of all proteins was confirmed in yeast (Figure 2—figure supplement 1).

PikpNK-KE shows increased binding to effector variants in vivo and in vitro when compared to wild type Pikp.
(A) Yeast-two-hybrid assay of Pikp-HMA and Pikp-HMANK-KE with AVR-Pik alleles. Control plates for yeast growth are on the left, with quadruple dropout media supplemented with X-α-gal and increasing concentrations of Aureobasidin A on the right for each combination of HMA/AVR-Pik. The unrelated M. oryzae effector AVR-Pii was used as a negative control. Growth and the development of blue coloration in the selection plate are both indicative of protein–protein interaction. HMA domains were fused to the GAL4 DNA binding domain, and AVR-Pik alleles to the GAL4 activator domain. Each experiment was repeated a minimum of three times, with similar results. (B) Box plots showing %Rmax, as measured by surface plasmon resonance, for Pikp-HMA and Pikp-HMANK-KE with the AVR-Pik effectors alleles at an HMA concentration of 40 nM. Pikp-HMA and Pikp-HMANK-KE are represented by blue and purple boxes, respectively. The centre line within each box represents the median, the box limits are the upper and lower quartiles, the whiskers extend to the largest value within Q1 – 1.5 × the interquartile range (IQR) and the smallest value within Q3 + 1.5 × IQR. All of the data points are represented as dots with distinct colors for each biological replicate. For each experiment, three biological replicates with three internal repeats were performed. ‘p’ is the p-value obtained from statistical analysis and Tukey’s HSD. For results of experiments with 4 nM and 100 nM HMA protein concentration, see Figure 2—figure supplement 2. (C) Co-immunoprecipitation of full length Pikp-1 and Pikp-1NK-KE with AVR-Pik variants. N-terminally 4xMyc tagged AVR-Pik effectors were transiently co-expressed with Pikp-1:6xHis3xFLAG (left) or with Pikp-1NK-KE:6xHis3xFLAG (right) in N. benthamiana. Immunoprecipitates (IPs) obtained with anti-FLAG antiserum, and total protein extracts, were probed with appropriate antisera. The dashed line indicates a crop site on the same blot used to compose the figure. Each experiment was repeated at least three times, with similar results. The asterisks mark the Pik-1 band. PS = Ponceau Stain.
-
Figure 2—source data 1
Surface plasmon resonance measurements used in the preparation of Figure 2B.
- https://doi.org/10.7554/eLife.47713.018
Next, we produced the Pikp-HMANK-KE domain protein via overexpression in E. coli and purified it to homogeneity using well-established procedures for these domains (see Materials and methods; Maqbool et al., 2015; De la Concepcion et al., 2018). Using SPR, we measured the binding affinity of the Pikp-HMANK-KE domain (alongside both wild-type Pikp-HMA and Pikm-HMA) to AVR-Pik effectors (Figure 2B, Figure 2—figure supplements 2 and 3). Response units (RU) were measured following the injection of Pik-HMAs at three different concentrations, after capturing AVR-Pik effectors on a Biacore NTA chip. RUs were then normalized to the theoretical maximum response (Rmax), assuming a 2:1 interaction model for Pikp-HMA and Pikp-HMANK-KE, and a 1:1 interaction model for Pikm-HMA, as previously described (De la Concepcion et al., 2018). These data showed an increased binding of Pikp-HMANK-KE to AVR-PikE and AVR-PikA compared to that of wild-type Pikp (Figure 2B, Figure 2—figure supplement 2, Supplementary file 2). The binding of Pikp-HMANK-KE to the AVR-Pik effectors was also greater than that of Pikm-HMA (Figure 2—figure supplement 3, Supplementary file 2), correlating with the results of cell-death assays (Figure 1—figure supplement 2). In both cases, statistical analysis has been carried out for the 40 nM analyte (Pik-HMA) data as a representative concentration. Although neither Pikp-HMA nor Pikm-HMA domains show binding to AVR-PikC by SPR, we observe a gain-of-binding of this effector variant with Pikp-HMANK-KE (Figure 2B, Figure 2—figure supplements 2 and 3, Supplementary file 2).
These results show that the Pikp-HMANK-KE mutant has a higher binding affinity for effectors AVR-PikE and AVR-PikA than does the wild-type protein. This suggests that the increased binding affinity to the HMA domain correlates with the expanded cell-death response in planta (Figure 1B).
The engineered PikpNK-KE mutant expands the association of full-length Pik-1 with effector variants in planta
In addition to interaction with the isolated HMA domain, we tested whether the Asn261Lys/Lys262Glu mutant could expand effector variant binding in the context of the full-length NLR. After generating the mutant in the full-length protein, we co-expressed either Pikp-1 or Pikp-1NK-KE with the AVR-Pik effector variants in N. benthamiana, followed by immunoprecipitation and western blotting to determine effector association.
AVR-PikD shows a robust association with Pikp-1. However, although we observe limited binding to the isolated Pikp-HMA domain in Y2H and SPR screens, we did not detect association of AVR-PikE and AVR-PikA with the full-length Pikp-1 in planta (Figure 2C). By contrast, we observe clear association of AVR-PikE and AVR-PikA with the Pikp-1NK-KE mutant, albeit with reduced intensity compared to AVR-PikD, correlating with the hierarchical cell-death response observed in planta (Figure 1C).
We also detect a very low level of association between full-length Pikp-1NK-KE and AVR-PikC (Figure 2C). However, co-expression of Pikp-1NK-KE and AVR-PikC does not result in macroscopic cell death in N. benthamiana (Figure 1C).
These results show that effector variant association with full-length Pikp-1 and Pikp-1NK-KE correlates with the in planta cell-death response (Figure 1C).
The effector-binding interface in the PikpNK-KE mutant adopts a Pikm-like conformation
Having established that the PikpNK-KE mutant displays an expanded effector recognition profile compared to wild-type Pikp, we sort to determine the structural basis of this activity. To this end, we determined the crystal structures of Pikp-HMANK-KE bound to AVR-PikD, and to AVR-PikE. We obtained samples of Pikp-HMANK-KE/AVR-PikD and Pikp-HMANK-KE/AVR-PikE complexes by co-expression in E. coli (described in the Materials and methods and by De la Concepcion et al., 2018). Each complex was crystallized (see Materials and methods) and X-ray diffraction data were collected at the Diamond Light Source (Oxford, UK) to 1.6 Å and 1.85 Å resolution, respectively. The details of X-ray data collection, structure solution, and completion are given in the Materials and methods and in Table 1.
Data collection and refinement statistics
https://doi.org/10.7554/eLife.47713.019PikpNK-KE–AVR-PikD | PikpNK-KE–AVR-PikE | |
---|---|---|
Data collection statistics | ||
Wavelength (Å) | 0.9763 | 0.9763 |
Space group | P 21 21 21 | P 21 21 21 |
Cell dimensions | ||
a, b, c (Å) | 29.79, 65.33, 75.86 | 66.46, 80.70, 105.58 |
Resolution (Å)* | 32.80–1.60 (1.63–1.60) | 29.50–1.85 (1.89–1.85) |
Rmerge (%)# | 8.1 (97.1) | 5.2 (75.1) |
I/σI# | 16.1 (2.6) | 31.0 (4.1) |
Completeness (%)# | 100 (100) | 99.8 (97.8) |
Unique reflections# | 20,294 (978) | 49337 (2963) |
Redundancy# | 12.8 (13.3) | 18.3 (17.8) |
CC(1/2) (%)# | 99.9 (86.6) | 100 (95.2) |
Refinement and model statistics | ||
Resolution (Å) | 32.82–1.60 (1.64–1.60) | 29.52–1.85 (1.90–1.85) |
Rwork/Rfree (%)^ | 19.7/23.2 (25.5/27.3) | 18.6/23.0 (29.1/35.0) |
No. atoms (Protein) | 1277 | 3604 |
B-factors (Protein) | 25.6 | 39.7 |
R.m.s. deviations^ | ||
Bond lengths (Å) | 0.009 | 0.012 |
Bond angles (°) | 1.5 | 1.4 |
Ramachandran plot (%)** | ||
Favored | 98.1 | 97.3 |
Outliers | 0 | 0.2 |
MolProbity Score | 1.41 (93th percentile) | 1.59 (91st percentile) |
-
*The highest resolution shell is shown in parenthesis.
#As calculated by Aimless, ^As calculated by Refmac5, **As calculated by MolProbity.
The overall architecture of these complexes is the same as that observed for all Pik-HMA/AVR-Pik effector structures. A key interaction at interface 3, one of the previously defined Pik-HMA–AVR-Pik interfaces (De la Concepcion et al., 2018), involves a lysine residue (Lys262 in Pikp and Pikm) that forms intimate contacts within a pocket on the effector surface (Figure 3). In order to position this lysine in the effector pocket, Pikp has to loop-out regions adjacent to this residue, compromising the packing at the interface (De la Concepcion et al., 2018, Figure 3A [left panel], B [left panel], C and D). By contrast, in Pikm, where the position of the lysine is shifted one residue to the N-terminus, no looping-out is required to locate the lysine within the pocket (Figure 3A [right panel], B [right panel], C and D). In the PikpNK-KE mutant, the position of this key lysine is shifted one residue to the N-terminus compared to that in the wild-type, so that it occupies the same position in the PikpNK-KE sequence as in the Pikm sequence. In the crystal structures of Pikp-HMANK-KE in complex with either AVR-PikD or AVR-PikE, we see that this region of the HMA adopts a Pikm-like conformation (Figure 3A [middle panel], B [middle panel], C and D), with no looping-out of the preceding structure. This confirms that in making the PikpNK-KE mutant, we have resurfaced Pikp to have a more robust, Pikm-like interface in this region.

The PikpNK-KE-HMA mutant adopts a Pikm-like conformation at the effector-binding interface.
Schematic view of the different conformations adopted by Pikp-HMA, Pikp-HMANK-KE and Pikm-HMA at interface 3 when in complex with AVR-PikD or AVR-PikE. In each panel, the effector is shown as sticks with the molecular surface also shown and colored as labeled. Pik-HMA residues are colored as labeled and shown in the Cα-worm with side-chain representation. (A) Schematic of Pikp-HMA (left), Pikp-HMANK-KE (middle) and Pikm-HMA (right) bound to AVR-PikD. Important residues in the HMA–effector interaction are labeled as shown. (B) Schematic of HMA residues as for panel (A), but bound to AVR-PikE. (C) Superposition showing Pikp-HMA, Pikp-HMANK-KE and Pikm-HMA chains (colored in blue, purple and yellow, respectively) bound to AVR-PikD. For clarity, only the Lys-261/262 side chain is shown. (D) Superposition as described before, but bound to AVR-PikE.
We found only limited structural perturbations at either of the other previously defined interfaces (interface 1 or 2 [De la Concepcion et al., 2018]) resulting from the binding of the AVR-PikD or AVR-PikE effectors to Pikp-HMA or Pikp-HMANK-KE (Figure 3—figure supplement 1). We therefore conclude that the effects of the PikpNK-KE mutant on protein function are mediated by altered interactions at interface 3.
Mutation in AVR-Pik effectors at the engineered binding interface impacts the in planta response and in vivo binding
To further confirm that the engineered binding interface is responsible for the expanded recognition of AVR-PikE and AVR-PikA by PikpNK-KE, we used mutants in the effectors at interface 2 (AVR-PikDH46E) and interface 3 (AVR-PikD,E,AE53R), which have previously been shown to impact the interactions and in planta responses of wild-type NLR alleles (De la Concepcion et al., 2018). We tested whether these mutants affected the cell-death response in N. benthamiana, the interactions between effectors and Pikp-HMANK-KE (using Y2H), and the interactions between effectors and full-length PikpNK-KE (using in planta co-immunoprecipitation [co-IP]).
First, we investigated the impact of mutation at interface 2 using the AVR-PikDH46E mutant. We found that cell death in N. benthamiana is essentially blocked when co-expressing either Pikp or PikpNK-KE with this mutant (AVR-PikDH46E), suggesting that the engineered NLR is still reliant on this interface for response (Figure 4A, Figure 4—figure supplement 1). Intriguingly, Y2H shows that the AVR-PikDH46E mutant displays some interaction with Pikp-HMANK-KE (Figure 4B), similar to this mutant’s interaction with Pikm-HMA (De la Concepcion et al., 2018), although this interaction is barely observed by co-IP with the full-length NLR (Figure 4C).

Mutation of AVR-Pik effectors at the engineered binding interface compromises binding and response.
(A) (Left) A representative leaf image showing Pikp-1NK-KE-mediated cell death induced by AVR-Pik variants and mutants as autofluorescence under UV light (the AVR-PikDH46E mutant is located at interface 2, whereas the AVR-PikDE53R, AVR-PikEE53R, and AVR-PikAE53R mutants are located at interface 3). Autofluorescence intensity is scored as in Figure 1. (Right) PikpNK-KE cell-death assay quantification in the form of dot plots. For each sample, the data points are represented as dots with a distinct color for each of the three biological replicates; these dots are plotted around the cell death score for visualization purposes. The size of the central dot at each cell death value is proportional to the number of replicates of the sample with that score. The number of repeats was 90. (B) Yeast-two-hybrid assay of Pikp-HMANK-KE with AVR-Pik variants and mutants. Control plate for yeast growth is on the left with quadruple dropout media supplemented with X-α-gal and increasing concentrations of Aureobasidin A on the right for each combination of HMA/AVR-Pik. The unrelated M. oryzae effector AVR-Pii was used as a negative control. Growth and the development of blue coloration in the selection plate are both indicative of protein–protein interaction. HMA domains were fused to the GAL4 DNA binding domain, and AVR-Pik alleles to the GAL4 activator domain. Each experiment was repeated a minimum of three times, with similar results. (C) Co-immunoprecipitation of full-length Pikp-1NK-KE with AVR-Pik variants and mutants. N-terminally 4xMyc tagged AVR-Pik effectors were transiently co-expressed with Pikp-1NK-KE:6xHis3xFLAG in N. benthamiana leaves. Immunoprecipitates (IPs) obtained with anti-FLAG antiserum, and total protein extracts, were probed with appropriate antisera. Each experiment was repeated at least three times, with similar results. The asterisks mark the Pik-1 band. PS = Ponceau Stain.
-
Figure 4—source data 1
Cell-death scoring data used in the preparation of Figure 4A.
- https://doi.org/10.7554/eLife.47713.025
Second, we investigated the impact of mutations at interface 3 using the Glu53Arg (E53R) mutant in AVR-PikD, AVR-PikE and AVR-PikA. We found that the AVR-PikDE53R mutant has essentially no effect on recognition of the effector by PikpNK-KE in N. benthamiana, and little or no effect on the interaction with Pikp-HMANK-KE or full-length PikpNK-KE (Figure 4A,B,C, Figure 4—figure supplement 1). By contrast, the equivalent mutation in AVR-PikE and AVR-PikA restricted the cell-death response in N. benthamiana, reduced the binding to Pikp-HMANK-KE in Y2H (as shown by the reduced blue coloration) and produced a less intense band for the effector following PikpNK-KE co-IP (Figure 4A,B,C). The expression of all proteins in yeast was confirmed by western blot (Figure 4—figure supplement 2).
These results provide evidence to show that while interactions across interface 2 remain important for the PikpNK-KE interaction with AVR-Pik effectors, it is the altered interaction at interface 3, as observed in the structures, that is responsible for the expanded recognition profile of this engineered mutant.
Discussion
Plants, including food crops, are under continuous threat from pathogens and pests, and new solutions to control disease are required. Although largely elusive to date, engineered plant NLR-type intracellular immune receptors have the potential to improve disease-resistance breeding (Dangl et al., 2013; Rodriguez-Moreno et al., 2017). NLR-integrated domains are a particularly attractive target for protein engineering because they interact directly with pathogen effectors (or host effector targets). Further, where tested, their binding affinities in vitro correlate with in planta immunity phenotypes (Maqbool et al., 2015; De la Concepcion et al., 2018; Guo et al., 2018), allowing biochemical and structural techniques to inform NLR design directly.
Here, we show that the recognition profile of the rice NLR Pikp can be expanded to include different AVR-Pik variants by engineering the binding interface between these proteins. This strengthens the hypothesis that tighter binding affinity between effectors and integrated HMA domains correlates with increased immune signaling in plants. This was previously shown for both natural alleles of Pik (Maqbool et al., 2015; De la Concepcion et al., 2018) and for Pia (Guo et al., 2018), but is now also shown for an engineered NLR. We propose this may be a general model for integrated domains that directly bind effectors.
Natural variation in Pik NLRs has given rise to different effector recognition profiles, and contributions from different binding interfaces were suggested to underpin this phenotype (De la Concepcion et al., 2018). In particular, a more favorable interaction at one interface (interface 3) in Pikm, compared to that in Pikp, was concluded to have evolved to compensate for changes in binding at a different site (interface 2). Here, through mutation of residues in Pikp (forming PikpNK-KE), we have combined favorable interfaces from Pikp and Pikm into a single protein. This has resulted in an expanded recognition phenotype for effector binding and response in planta. In vitro, the binding of PikpNK-KE-HMA to effector variants is consistently higher than that of Pikm-HMA (Figure 2—figure supplement 3). This is likely to be underpinned by differences in a cluster of key residues within interface 2 (spatially equivalent residues Val222, Lys228, and Glu230 in Pikp, and Ala223, Gln229, and Val231 in Pikm) that form more energetically favorable contacts, including water-mediated hydrogen bonds. However, despite a positive trend (which is significant in some cases), the levels of cell death mediated by Pikm and PikpNK-KE in planta are comparable (Figure 1—figure supplements 2 and 4).
Although the PikpNK-KE mutant did not deliver a cell-death response in N. benthamiana to the stealthy AVR-PikC effector variant, it did show a gain-of-binding in vitro, as well as weak binding in Y2H and in in planta co-IP. It is surprising that the gain-of-binding for AVR-PikC in vitro, unlike that for AVR-PikA, does not correlate directly with gain of in vivo interactions or cell-death response in N. benthamiana. We hypothesize that this gain-of-binding for AVR-PikC observed in vitro is not sufficiently robust, especially in the context of the full-length NLR, for binding or for triggering immune signaling. Nevertheless, this work sets the scene for future interface engineering experiments that may further improve the response profiles of Pik NLRs to currently unrecognized effector variants. It also suggests that future work should test the disease resistance profile of M. oryzae strains carrying the different effector variants in rice expressing the engineered receptor.
The integrated HMA domain in the NLR RGA5 (the sensor of the Pia NLR pair in rice) binds to M. oryzae effectors AVR1-CO39 and AVR-Pia via a different interface, and it has been suggested that these binding sites are mutually exclusive (Guo et al., 2018). This raises the possibility that an HMA domain could be engineered to bind and respond to multiple effectors (Guo et al., 2018). Recently, the Pikp-HMA domain was shown to interact with AVR-Pia at the same interface as that used by the RGA5-HMA domain, and this probably underpins the partial resistance of rice plants to M. oryzae expressing AVR-Pia (Varden et al., 2019). This presents a starting point for the use of Pikp as a chassis for such studies. Although it remains to be seen whether any such resurfaced HMA domain can bind to multiple effectors, these studies suggest that this strategy has potential as a novel approach.
Plant breeding is required to provide new genetic solutions to disease resistance in crops. This is necessary to limit the environmental and social damage caused by pesticides, and to deal with changes in climate and in the globalization of agriculture that result in the spread of pathogens and pests into new environments (Islam et al., 2016; Michelmore et al., 2017; Savary et al., 2019). Classical breeding for disease resistance has been limited by issues such as linkage drag and hybrid incompatibility, as also seen in model plant species (Bomblies et al., 2007). Novel molecular approaches such as engineering ‘decoys’ (Kim et al., 2016) and protein resurfacing, as described here, combined with modern transformation (Altpeter et al., 2016) and breeding pipelines (Watson et al., 2018), offer the opportunity for more targeted approaches to breeding for disease resistance. These strategies will complement other emerging technologies in NLR identification (Arora et al., 2019) and NLR stacking (Dangl et al., 2013) as methods to develop improved crops for the future.
Accession codes
Protein structures, and the data used to derive them, have been deposited at the Protein DataBank (PDB) with accession codes 6R8K (Pikp-HMANK-KE/AVR-PikD) and 6R8M (Pikp-HMANK-KE/AVR-PikE).
Materials and methods
Gene cloning
Request a detailed protocolFor in vitro studies, Pikp-HMANK-KE (encompassing residues 186 to 263) was amplified from wild-type Pikp-HMA by introducing the mutations in the reverse primer, followed by cloning into pOPINM (Berrow et al., 2007). The wild-type Pikp-HMA, Pikm-HMA, and AVR-Pik expression constructs used in this study are as described in De la Concepcion et al. (2018).
For Y2H, we cloned Pikp-HMANK-KE (as above) into pGBKT7 using an In-Fusion cloning kit (Takara Bio USA), following the manufacturer’s protocol. The wild-type Pikp-HMA domain in pGBKT7 and AVR-Pik effector variants in pGADT7 were generated as described in De la Concepcion et al. (2018).
For protein expression in planta, the Pikp-HMANK-KE domain was generated using site-directed mutagenesis by introducing the mutations in the reverse primer. This domain was then assembled into a full-length construct using Golden Gate cloning (Engler et al., 2008) and into the plasmid pICH47742 with a C-terminal 6xHis/3xFLAG tag. Expression was driven by the Agrobacterium tumefaciens Mas promoter and terminator. Full-length Pikp-1, Pikp-2, and AVR-Pik variants were generated as described in De la Concepcion et al. (2018). All DNA constructs were verified by sequencing.
Expression and purification of proteins for in vitro binding studies
Request a detailed protocolpOPINM constructs encoding Pikp-HMA, Pikm-HMA and Pikp-HMANK-KE were produced in E. coli SHuffle cells (Lobstein et al., 2012), using the protocol described in De la Concepcion et al. (2018). Cell cultures were grown in autoinduction media (Studier, 2005) at 30°C for 5–7 hr and then at 16°C overnight. Cells were harvested by centrifugation and re-suspended in 50 mM Tris-HCl (pH 7.5), 500 mM NaCl, 50 mM glycine, 5% (vol/vol) glycerol, and 20 mM imidazole supplemented with EDTA-free protease inhibitor tablets (Roche). Cells were sonicated and, following centrifugation at 40,000xg for 30 min, the clarified lysate was applied to a Ni2+-NTA column connected to an AKTA Xpress purification system (GE Healthcare). Proteins were step-eluted with elution buffer (50 mM Tris-HCl [pH7.5], 500 mM NaCl, 50 mM glycine, 5% (vol/vol) glycerol, and 500 mM imidazole) and directly injected onto a Superdex 75 26/60 gel filtration column pre-equilibrated with 20 mM HEPES (pH 7.5) and 150 mM NaCl. Purification tags were removed by incubation with 3C protease (10 μg/mg fusion protein) followed by passing through tandem Ni2+-NTA and MBP Trap HP columns (GE Healthcare). The flow-through was concentrated as appropriate and loaded onto a Superdex 75 26/60 gel filtration column for final purification and buffer exchange into 20 mM HEPES (pH 7.5) and 150 mM NaCl.
AVR-Pik effectors, with either a 3C protease-cleavable N-terminal SUMO or a MBP tag and with a non-cleavable C-terminal 6xHis tag, were produced in and purified from E. coli SHuffle cells as previously described (Maqbool et al., 2015; De la Concepcion et al., 2018). All protein concentrations were determined using a Direct Detect Infrared Spectrometer (Merck).
Co-expression and purification of Pik-HMA and AVR-Pik effectors for crystallization
Request a detailed protocolPikp-HMANK-KE was co-expressed with AVR-PikD or AVR-PikE effectors in E. coli SHuffle cells following co-transformation of pOPINM:Pikp-HMANK-KE and pOPINA:AVR-PikD/E (which were prepared as described in De la Concepcion et al., 2018). Cells were grown in autoinduction media (supplemented with both carbenicillin and kanamycin), harvested, and processed as described in De la Concepcion et al. (2018). Protein concentrations were measured using a Direct Detect Infrared Spectrometer (Merck).
Protein–protein interaction: yeast-two-hybrid analyses
Request a detailed protocolTo detect protein–protein interactions between Pik-HMAs and AVR-Pik effectors in a yeast two-hybrid system, we used the Matchmaker Gold System (Takara Bio USA). Plasmid DNA encoding Pikp-HMANK-KE in pGBKT7, generated in this study, was co-transformed into chemically competent Y2HGold cells (Takara Bio, USA) with the individual AVR-Pik variants or mutants in pGADT7 , as described previously (De la Concepcion et al., 2018). Single colonies grown on selection plates were inoculated in 5 ml of SD-Leu-Trp overnight at 30 °C. Saturated culture was then used to make serial dilutions of OD600 1, 1−1, 1−2, and 1−3, respectively. 5 μl of each dilution was then spotted on a SD-Leu-Trp plate as a growth control, and on a SD-Leu-Trp-Ade-His plate containing X-α-gal and supplemented with Aureobasidin A (Takara Bio, USA). Plates were imaged after incubation for 60–72 hr at 30 °C. Each experiment was repeated a minimum of three times, with similar results.
To confirm protein expression in yeast, total protein extracts from transformed colonies were produced by boiling the cells for 10 min in LDS Runblue sample buffer. Samples were centrifugated and the supernatant was subjected to SDS-PAGE gels prior to western blotting. The membranes were probed with anti-GAL4 DNA-BD (Sigma) for the HMA domains in pGBKT7 and with the anti-GAL4 activation domain (Sigma) antibodies for the AVR-Pik effectors in pGADT7.
Protein–protein interaction: surface plasmon resonance
View detailed protocolSurface plasmon resonance (SPR) experiments were performed on a Biacore T200 system (GE Healthcare) using an NTA sensor chip (GE Healthcare). The system was maintained at 25°C, and a flow rate of 30 μl/min was used. All proteins were prepared in SPR running buffer (20 mM HEPES [pH 7.5], 860 mM NaCl, 0.1% Tween 20). C-terminally 6xHis-tag AVR-Pik variants were immobilized on the chip, giving a response of 200 ± 100. The sensor chip was regenerated between each cycle with an injection of 30 μl of 350 mM EDTA.
For all the assays, the level of binding was expressed as a percentage of the theoretical maximum response (Rmax) normalized for the amount of ligand immobilized on the chip. The cycling conditions were the same as those used in De la Concepcion et al. (2018). For each measurement, in addition to subtracting the response in the reference cell, a further buffer-only subtraction was made to correct for bulk refractive index changes or machine effects (Myszka, 1999). SPR data were exported and plotted using R v3.4.3 (https://www.r-project.org/) and the function ggplot2 (Wickham, 2009). Each experiment was repeated a minimum of three times, including internal repeats, with similar results. The proteins used came from three independent preparations for the HMA domains and two independent preparations of the AVR-Pik effectors.
Protein–protein interaction: in planta co-immunoprecipitation (co-IP)
Request a detailed protocolTransient gene-expression in planta for Co-IP was performed by delivering T-DNA constructs within the A. tumefaciens GV3101 strain into 4-week-old N. benthamiana plants grown at 22–25°C with high-light intensity. A. tumefaciens strains carrying Pikp-1 or Pikp −1NK-KE were mixed with strains carrying the AVR-Pik effector, at OD600 0.2 each, in agroinfiltration medium (10 mM MgCl2 and 10 mM 2-(N-morpholine)-ethanesulfonic acid (MES) [pH5.6]), supplemented with 150 µM acetosyringone. For detection of complexes in planta, leaf tissue was collected 3 days post infiltration (dpi), frozen, and ground to fine powder in liquid nitrogen using a pestle and mortar. Leaf powder was mixed with two times weight/volume ice-cold extraction buffer (10% glycerol, 25 mM Tris [pH 7.5], 1 mM EDTA, 150 mM NaCl, 2% w/v PVPP, 10 mM DTT, 1x protease inhibitor cocktail [Sigma], 0.1% Tween 20 [Sigma]), centrifuged at 4,200 g/4°C for 20–30 min, and the supernatant was passed through a 0.45 μm Minisart syringe filter. The presence of each protein in the input was determined by SDS-PAGE and western blot. Pik-1 and AVR-Pik effectors were detected by probing the membrane with anti-FLAG M2 antibody (SIGMA) and anti c-Myc monoclonal antibody (Santa Cruz), respectively. For immunoprecipitation, 1.5 ml of filtered plant extract was incubated with 30 μl of M2 anti-FLAG resin (Sigma) in a rotatory mixer at 4°C. After three hours, the resin was pelleted (800 g, 1 min) and the supernatant removed. The pellet was washed and resuspended in 1 ml of IP buffer (10% glycerol, 25 mM Tris [pH 7.5], 1 mM EDTA, 150 mM NaCl, 0.1% Tween 20 [Sigma]) and pelleted again by centrifugation as before. Washing steps were repeated five times. Finally, 30 μl of LDS Runblue sample buffer was added to the agarose and incubated for 10 min at 70°C. The resin was pelleted again, and the supernatant loaded onto SDS-PAGE gels prior to western blotting. Membranes were probed with anti-FLAG M2 (Sigma) and anti c-Myc (Santa Cruz) monoclonal antibodies.
N. benthamiana cell-death assays
Request a detailed protocolA. tumefaciens GV3101 vectors carrying Pikp-1, Pikm-1, or Pikp-1NK-KE were resuspended in induction media (10 mM MES [pH 5.6], 10 mM MgCl2 and 150 µM acetosyringone) and mixed with Pikp-2 (or Pikm-2 in the Pikm cell-death assay), AVR-Pik effectors, and P19 at OD600 0.4, 0.4, 0.6 and 0.1, respectively. Four-week-old N. benthamiana leaves were infiltrated using a needleless syringe. Leaves were collected at 5 dpi to measure UV autofluorescence (a proxy for cell death) or ion leakage.
Cell-death scoring: UV autofluorescence
Request a detailed protocolDetached leaves were imaged at 5 dpi from the abaxial side of the leaves for UV fluorescence images. Photos were taken using a Nikon D4 camera with a 60 mm macro lens, ISO set 1600 and exposure ~10 secs at F14. The filter was a Kodak Wratten No.8 and white balance was set to 6250 degrees Kelvin. Blak-Ray longwave (365 nm) B-100AP spot light lamps were moved around the subject during the exposure to give an even illumination. Images shown are representative of three independent experiments, with internal repeats. The cell death index used for scoring was as presented previously (Maqbool et al., 2015). Dotplots were generated using R v3.4.3 (https://www.r-project.org/) and the graphic package ggplot2 (Wickham, H., 2009). The size of the centre dot at each cell death value is directly proportional to the number of replicates in the sample with that score. All individual data points are represented as dots.
Cell-death scoring: ion leakage
Request a detailed protocolFor ion-leakage quantification, plants were infiltrated with the relevant constructs on two different leaves. At 5 dpi, leaves were detached and two leaf discs with a diameter of 8 mm (one disc per leaf spot) were collected and floated in 1.5 mL of Milli-Q water. Conductivity (µS/cm) was measured immediately after transferring the leaf disc to water (time zero) using a LAQUAtwin EC-33 conductivity meter (Horiba UK Ltd). Leaf discs were then incubated for 6 h at room temperature with gentle shaking before measuring the final conductivity. The assay was carried out in six biological replicates with a total of 40 technical replicates (2 discs x 40 plants). Conductivity data for each sample were exported and plotted using R v3.4.3 (https://www.r-project.org/) and the function ggplot2 (Wickham, H., 2009).
Crystallization, data collection and structure solution
Request a detailed protocolFor crystallization, Pikp-HMANK-KE in complex with AVR-PikD or AVR-PikE were concentrated following gel filtration. Sitting drop vapor diffusion crystallization trials were set up in 96 well plates, using an Oryx nano robot (Douglas Instruments, United Kingdom). Plates were incubated at 20°C, and crystals typically appeared after 24–48 hr. For data collection, all crystals were harvested from the Morpheus HT-96 screen (Molecular Dimensions), and snap-frozen in liquid nitrogen. Crystals used for data collection appeared from the following conditions: (i) Pikp-HMANK-KE/AVR-PikD (10 mg/ml), Morpheus HT-96 condition D4 (0.12 M alcohols [0.2 M 1,6-hexanediol, 0.2 M 1-butanol, 0.2 M 1,2-propanediol, 0.2 M 2-propanol, 0.2 M 1,4-butanediol, and 0.2 M 1,3-propanediol], 0.1 M Buffer system 1 [1 M imidazole, MES monohydrate (acid) (pH 6.5)], and 50% v/v precipitant mix 4 [25%v/v MPD, 25% v/v PEG 1000, 25% v/v PEG3350]); (ii) Pikp-HMANK-KE/AVR-PikE (15 mg/ml), Morpheus HT-96 condition A8 (0.06M divalents [0.3 M magnesium chloride hexahydrate, 0.3 M calcium chloride dihydrate]), 0.1M Buffer system 2 (sodium HEPES; MOPS [acid] [pH 7.5]), 37.5% v/v Precipitant mix 4 (25%v/v MPD, 25% v/v PEG 1000, 25% v/v PEG3350).
X-ray data sets were collected at the Diamond Light Source using beamline i03 (Oxford, UK). The data were processed using the xia2 pipeline (Winter, 2010) and CCP4 (Winn et al., 2011). The structures were solved by molecular replacement using PHASER (McCoy et al., 2007) and the coordinates of AVR-PikD and a monomer of Pikp-HMA from PDB entry 6G10. The final structures were obtained through iterative cycles of manual rebuilding and refinement using COOT (Emsley et al., 2010) and REFMAC5 (Murshudov et al., 2011), as implemented in CCP4 (Winn et al., 2011). Structures were validated using the tools provided in COOT and MOLPROBITY (Chen et al., 2010).
Statistical analyses
Request a detailed protocolQualitative cell-death scoring from autofluorescence was analyzed using estimation methods (Ho et al., 2019) and visualized with estimation graphics using the besthr R library (MacLean, 2019). Briefly, in this process all autofluorescence (cell-death) scores in samples under comparison were ranked, irrespective of sample. The mean ranks of the control and test sample were taken and a bootstrap process was begun on ranked test data, in which samples of equal size to the experiment were replaced and the mean rank calculated. After 1000 bootstrap samples, rank means were calculated, a distribution of the mean ranks was drawn and its 2.5 and 97.5 quantiles calculated. If the mean of the control data is outside of these boundaries, the control and test means were considered to be different. Quantitative data from ion-leakage (cell-death) and SPR assays were analyzed by preparing a linear mixed effects model of sample on ion leakage/SPR and post-hoc comparisons performed for sample contrasts using Tukey’s HSD method in the R package nlme (Pinheiro and Bates, 2019) and in lsmeans (Lenth, 2016).
Data availability
Protein structures, and the data used to derive these, have been deposited at the Protein DataBank (PDB) with accession codes 6R8K (Pikp-HMANK-KE/AVR-PikD) and 6R8M (Pikp-HMANK-KE/AVR-PikE).
-
Protein Data BankID 6R8K. Complex of rice blast (Magnaporthe oryzae) effector protein AVR-PikD with an engineered HMA domain of Pikp-1 from rice (Oryza sativa).
-
Protein Data BankID 6R8M. Complex of rice blast (Magnaporthe oryzae) effector protein AVR-PikE with an engineered HMA domain of Pikp-1 from rice (Oryza sativa).
References
-
Advancing crop transformation in the era of genome editingThe Plant Cell 28:1510–1520.https://doi.org/10.1105/tpc.16.00196
-
Lessons in effector and NLR biology of Plant-Microbe systemsMolecular Plant-Microbe Interactions 31:34–45.https://doi.org/10.1094/MPMI-08-17-0196-FI
-
Convergent evolution of effector protease recognition by Arabidopsis and barleyMolecular Plant-Microbe Interactions : MPMI 32:550–565.https://doi.org/10.1094/MPMI-07-18-0202-FI
-
Scratching the surface: resurfacing proteins to endow new properties and functionCell Chemical Biology 23:543–553.https://doi.org/10.1016/j.chembiol.2016.04.010
-
MolProbity: all-atom structure validation for macromolecular crystallographyActa Crystallographica Section D Biological Crystallography 66:12–21.https://doi.org/10.1107/S0907444909042073
-
NB-LRR proteins: pairs, pieces, perception, partners, and pathwaysCurrent Opinion in Plant Biology 13:472–477.https://doi.org/10.1016/j.pbi.2010.04.007
-
Features and development of cootActa Crystallographica. Section D, Biological Crystallography 66:486–501.https://doi.org/10.1107/S0907444910007493
-
Engineering a decoy substrate in soybean to enable recognition of the soybean mosaic virus NIa proteaseMolecular Plant-Microbe Interactions 32:760–769.https://doi.org/10.1094/MPMI-12-18-0324-R
-
Arms race co-evolution of Magnaporthe oryzae AVR-Pik and rice pik genes driven by their physical interactionsThe Plant Journal : For Cell and Molecular Biology 72:894–907.https://doi.org/10.1111/j.1365-313X.2012.05110.x
-
Least-Squares means: the R package lsmeansJournal of Statistical Software 69:1–33.https://doi.org/10.18637/jss.v069.i01
-
Phaser crystallographic softwareJournal of Applied Crystallography 40:658–674.https://doi.org/10.1107/S0021889807021206
-
Foundational and translational research opportunities to improve plant healthMolecular Plant-Microbe Interactions 30:515–516.https://doi.org/10.1094/MPMI-01-17-0010-CR
-
REFMAC5 for the refinement of macromolecular crystal structuresActa Crystallographica Section D Biological Crystallography 67:355–367.https://doi.org/10.1107/S0907444911001314
-
Improving biosensor analysisJournal of Molecular Recognition 12:279–284.https://doi.org/10.1002/(SICI)1099-1352(199909/10)12:5<279::AID-JMR473>3.0.CO;2-3
-
Transfer and engineering of immune receptors to improve recognition capacities in cropsCurrent Opinion in Plant Biology 38:42–49.https://doi.org/10.1016/j.pbi.2017.04.010
-
The global burden of pathogens and pests on major food cropsNature Ecology & Evolution 3:430–439.https://doi.org/10.1038/s41559-018-0793-y
-
Single amino acid mutations in the potato immune receptor R3a expand response to Phytophthora effectorsMolecular Plant-Microbe Interactions 27:624–637.https://doi.org/10.1094/MPMI-02-14-0040-R
-
Protein production by auto-induction in high density shaking culturesProtein Expression and Purification 41:207–234.https://doi.org/10.1016/j.pep.2005.01.016
-
How to build a pathogen detector: structural basis of NB-LRR functionCurrent Opinion in Plant Biology 15:375–384.https://doi.org/10.1016/j.pbi.2012.05.001
-
Bookggplot2: Elegant Graphics for Data AnalysisNew York, United States: Springer-Verlag.https://doi.org/10.1007/978-0-387-98141-3
-
Overview of the CCP4 suite and current developmentsActa Crystallographica. Section D, Biological Crystallography 67:235–242.https://doi.org/10.1107/S0907444910045749
-
xia2 : an expert system for macromolecular crystallography data reductionJournal of Applied Crystallography 43:186–190.https://doi.org/10.1107/S0021889809045701
Decision letter
-
Thorsten NürnbergerReviewing Editor; University of Tübingen, Germany
-
Detlef WeigelSenior Editor; Max Planck Institute for Developmental Biology, Germany
-
Thorsten NürnbergerReviewer; University of Tübingen, Germany
In the interests of transparency, eLife includes the editorial decision letter and accompanying author responses. A lightly edited version of the letter sent to the authors after peer review is shown, indicating the most substantive concerns; minor comments are not usually included.
Thank you for submitting your article "Protein engineering expands the effector recognition profile of a rice NLR immune receptor" for consideration by eLife. Your article has been reviewed by three peer reviewers, including Thorsten Nürnberger as the Reviewing Editor and Reviewer #1, and the evaluation has been overseen by Detlef Weigel as the Senior Editor.
The reviewers have discussed the reviews with one another and the Reviewing Editor has drafted this decision to help you prepare a revised submission.
Summary:
Engineering of plant NLR-type immune receptors for expanded effector recognition is an exciting and promising tool for future molecular breeding approaches. All reviewers agreed that this study provides a substantial step forward in our understanding on how to achieve this goal. While highly supportive of your work, the reviewers request some revisions, as listed below.
Essential revisions:
The quality of the Y2H data is not sufficient. Here, growth on triple or quadruple dropout medium supplemented with different 3AT concentrations that accentuate growth differences is required. Likewise, cell death phenotypes observed upon overexpression in N. benthamiana require proper quantification and statistical evaluation.
Reviewer #1:
This is a really nice and technically challenging complement to a number of papers of this laboratory. There is no doubt that it would make a nice contribution to eLife. My only concern is about quantification of cell death in this study. Here, the authors may resort to ion leakage measurements that are certainly better quantifiable than the cell death index chosen here.
Reviewer #2:
The manuscript at hand introduces an engineered variant of the rice NLR protein Pikp, the design of which was inspired by another rice NLR named Pikm. The authors argue that this Pikp variant, named PikpNK-KE binds the Magnaporthe oryzae effector variants AVR-PikA and AVR-PikE in addition to the cognate ligand AVR-PikD of wild-type Pikp. Binding of PikpNK-KE to AVR-PikA and AVR-PikE was shown in various ways including yeast-two-hybrid assay (Y2H), cell death assay in N. benthamiana, co-immunoprecipitation (Co-IP) and surface plasmon resonance (SPR) (the latter two using recombinant protein from E. coli, HMA domain only). In addition, the authors determined two crystal structures of the HMA domains of PikpNK-KE in complex with AVR-PikD and of PikpNK-KE in complex with AVR-PikE at 1.6 Å and 1.85 Å, respectively. The authors compare these crystal structures with previously published crystal structures of wild-type Pikp and Pikm bound to AVR-PikD and AVR-PikE. They conclude with using this structural information to develop effector variants that are not bound and unrecognized by the PikpNK-KE receptor, to which end they used cell death assay, Y2H and Co-IP.
As few studies exist that succeeded at engineering NLRs with novel specificity, the article is highly relevant to scientist within the field of plant immunity and beyond. I believe the study can be helpful and serve as a template to scientists who have similar endeavors of designing novel immune receptors. The novelty of the resolved structures lies in the mechanistic insight into effector binding of an engineered NLR.
For those reasons, I recommend publication of the manuscript in eLife.
However, major points need to be addressed before publication which are listed as follows:
1) Figure 2A and B: Why binding of AVR-PikC in case of PikpNK-KE? Requires more substantial discussion with special focus on differences in yeast, tobacco and in vitro.
2) Figure 2—figure supplement 3: Discuss (by using crystal structure) why binding of pHMANK-KE is consistently higher than mHMA.
3) Figure 4B: Control missing. Why not show mutation of N46 in AVR-PikE and corresponding residue in AVR-PikA?
Reviewer #3:
This article describes structure-guided engineering of the Pikp NLR immune receptor from rice that recognizes the effector AVR-PikD from the blast fungus Magnaporthe oryzae. Engineering of the integrated HMA decoy domain of Pikm1 domain was performed to expand its recognition spectrum to the AVR-Pik allelic variants AVR-PikA and AVR-PikE. For this, the work builds on previous work on structural and functional characterization of the HMA domains of Pikm1 and PikP1.
The manuscript is very interesting, extremely well written and documented and provides a first example of structure guided engineering of an NLR immune receptor for extended recognition specificity. However, several weaknesses require revision:
- An important drawback of the study is that there is no analysis of the engineered Pikp-1NK-KE mutant in rice plants to test whether it confers indeed an extended recognition specificity and resistance spectrum in transgenic rice.
- Another important weakness are the Yeast two hybrid data that are not convincing (cf specific comments).
- The N. benthamiana "cell death" assay should be improved by measuring fluorescence intensity on infiltrated spots and performing statistical analysis on these quantitative results.
- co-IP experiments with the isolated HMA domain (wt and mutant) would strengthen the study (c.f. specific comments).
Specific Comments:
Results and data, images and graphs "monitoring cell death in a well-established N. benthamiana assay": all images intended to show N. benthamiana cell death show autofluorescence, not cell death. Graphs are also based on autofluorescence induction and not cell death. Previous publications suggest that the cell death response is rather weak. Description of the experiments should be more accurate here and throughout the text. In addition, I do not understand why auto fluorescence is not quantified using image analysis software but interpreted visually (with a 6 level scoring system). Quantitative measurements combined with robust statistical analysis should be performed rather than scoring and qualitative interpretation of results from scoring. Differences between AVRPik-C and AVRPik-A on Pikp-1NK-KE on the photo in Figure 1B seems quite subtle in particular if one assumes that this is among the most differential images that were obtained.
Figure 1—figure supplement 2 -Statistical test to conclude on difference is missing. Unclear what is meant by "qualitatively stronger response".
Figure 2A and 4B – Y2H experiments provide only weak support for the author's claims (in certain cases no support at all): growth differences between the different yeast clones are, in most cases, rather subtle and often don’t correlate with the activity and in vitro binding data. E.g. in Figure 2A, Pikp-HMA/ARVPikD and Pikp-HMA/AVR-PikA grow identically and Pikp-HMA/AVR-PikE even grow better. Pikp-HMA/AVR-PikE growth is stronger than Pikp-HMA NK/KE/AVR-PikE. In Figure 4B, all clone seem to grow the same. Differences in color are not always easy to appreciate and rely on subjective interpretation. I would like to see growth on triple or quadruple dropout medium supplemented with different 3AT concentrations that accentuate growth differences and allow to unambiguously differentiate between growth strength of different yeast clones. If authors prefer to rely on β-Galactosidase activity, they should provide quantitative data from enzyme assays. Also apparent inconsistencies between Y2H and SPR or N. benthamiana cell death assays should be addressed in the text. For Pikp-HMANK/KE, a negative control (unrelated effector) should be provided since it seems to interact with all AVR-Pik constructs that were tested including some that are expected to not interact or to interact only very weakly (e.g. AVRPikD-H46E). It should be indicated on the figure that HMAs are fused to BD and Avr's to AD.
Figure 2C and Figure 4C: It would be interesting to see co-IP experiments with the isolated HMA domain (wt and mutant) to compare, on the one hand, to the other experiments with the isolated HMA domain (SPR and Y2H) and on the other to full-length Pikp1 variants.
Figure 2B and Figure 2—figure supplements 2 and 3 – SPR experiments: Binding of C is similar to binding of A and E. Statistical analysis should be performed to determine if there are differences. Why do authors not calculate affinities from SPR and discuss differences in affinities? One question to address would be, can the weak differences in affinity between C and A (or E) explain recognition specificity?
Figure 3—figure supplement 1: It looks to me as if the QtPISA interaction radar for Pikp-HMANK/KE/AVR-PikE does not support that binding is biologically relevant. Area is highly similar to the one for Pikp-HMA/AVR-PikE. Please comment on this point.
Introduction, fourth paragraph: recognition does not necessarily rely on binding. Can be posttranslational modification (e.g. RRS1) or other mechanism.
Subsection “The engineered PikpNK-KE mutant expands association of full-length Pik-1 to effector variants in planta”: Replace "interact" by "associate with".
Subsection “The engineered PikpNK-KE mutant expands association of full-length Pik-1 to effector variants in planta”, last paragraph: Replace "binding" by "association".
Varden et al., 2019: This study is not peer-reviewed or properly published (Bioarchive paper) and therefore not be cited.
https://doi.org/10.7554/eLife.47713.035Author response
Essential revisions:
The quality of the Y2H data is not sufficient. Here, growth on triple or quadruple dropout medium supplemented with different 3AT concentrations that accentuate growth differences is required. Likewise, cell death phenotypes observed upon overexpression in N. benthamiana require proper quantification and statistical evaluation.
Further details are provided below in the responses to each of the reviewers’ comments.
In brief, we have repeated our Y2H assays with increasing concentrations of Aureobasidin A, the antibiotic used in the Matchmaker Gold Y2H system (more stringent than 3AT in other systems). We screened different concentrations from 0.2µg/ml to 1µg/ml (0.2 µg/ml being the standard concentration in the Matchmaker protocol), our new figures show results for Aureobasidin A concentrations of 0.2, 0.4 and 0.8µg/ml. We hope you agree that this has improved the quality of the Y2H data (making it more robust), and that the results are clearer. Further, we have added a negative control, AVR-Pii, an unrelated effector from M. oryzae.
For the cell death assays, we do not think it appropriate to quantify the autofluorescence from the images as during acquisition the UV lamp is moved on a gantry and uniform light distribution is not guaranteed to allow this level of accuracy (we have added details to the Materials and methods about image acquisition). Therefore, we have chosen to develop the ion leakage assay to quantify cell death as an additional experiment. This has the added advantage of being another independent measure of immunity-related readouts in this study. As can be seen, the ion leakage assay compliments the scoring from autofluorescence very well. This includes a comparison between Pikp and Pikm that was the subject of our previous manuscript published in Nature Plants (direct comparison not included here, although is available).
Regarding statistics, as the data from the cell death assays is categorical, we use estimation methods to highlight differences. For the quantitative data (ion leakage assays and SPR binding data), we have used statistical analysis and Tukey’s HSD, and adjusted the text of the manuscript to take the results into account where necessary. In recognition of his contribution to this part of the manuscript, we have added an author, Dan MacLean, to the manuscript.
Reviewer #1:
This is a really nice and technically challenging complement to a number of papers of this laboratory. There is no doubt that it would make a nice contribution to eLife. My only concern is about quantification of cell death in this study. here, the authors may resort to ion leakage measurements that are certainly better quantifiable than the cell death index chosen here.
In response to this reviewer, and others, we have now performed ion leakage experiments from N. benthamiana tissue. The results are presented in the new Figure 1D, Figure 1—figure supplement 2, and Supplementary file 1. The main text and Materials and methods have been updated appropriately. As the cell death scoring (by autofluorescence) and ion leakage assays correlate very well for the experiments with the AVR-Pik variants, we did not extend the ion leakage assays to include the effector mutants described in Figure 4.
Reviewer #2:
[…] Major points need to be addressed before publication which are listed as follows:
1) Figure 2A and B: Why binding of AVR-PikC in case of PikpNK-KE? Requires more substantial discussion with special focus on differences in yeast, tobacco and in vitro.
We have edited the discussion concerning the binding of AVR-PikC to PikpNK-KE as measured in vitro (SPR), when there is only weak apparent binding by Y2H (accentuated in new results) and by co-IP in planta. We were surprised by the SPR result (especially in comparison to AVR-PikA). We think it is best interpreted as a difference between working in vitro with the purified Pik-HMA domain only vs. in vivo with the HMA domains (tagged with the GAL4-DNA binding domain, yeast) and with the full-length NLRs (in plants).
2) Figure 2—figure supplement 3: Discuss (by using crystal structure) why binding of pHMANK-KE is consistently higher than mHMA.
We have added a section to the Discussion on this point, and also put the in vitro binding data in the context of the in planta results.
3) Figure 4B: Control missing. Why not show mutation of N46 in AVR-PikE and corresponding residue in AVR-PikA?
An N46E mutation in AVR-PikE produces the same protein as AVR-PikDH46E, so this is already included in the experiment. While it is indeed the case that we have not included the N46E mutation in the AVR-PikA background, we feel this is an unnecessary control. Given our other data, it is highly likely that this mutation would have the same effect as in the AVR-PikD background.
Reviewer #3:
[…] The manuscript is very interesting, extremely well written and documented and provides a first example of structure guided engineering of an NLR immune receptor for extended recognition specificity. However, several weaknesses require revision:
- An important drawback of the study is that there is no analysis of the engineered Pikp-1NK-KE mutant in rice plants to test whether it confers indeed an extended recognition specificity and resistance spectrum in transgenic rice.
We appreciate the comment, but we consider this work is beyond the scope of the manuscript.
- Another important weakness are the Yeast two hybrid data that are not convincing (cf specific comments).
Please see comments below, and the responses elsewhere to address the quality of the Y2H.
- The N. benthamiana "cell death" assay should be improved by measuring fluorescence intensity on infiltrated spots and performing statistical analysis on these quantitative results.
As mentioned above, we don’t think that the way we measure the autofluorescence lends itself to quantification per say. We have performed ion leakage assays to strengthen the conclusions and present those new results here, along with statistical analyses (including using estimation methods to analyse the scored autofluorescence data).
- co-IP experiments with the isolated HMA domain (wt and mutant) would strengthen the study (c.f. specific comments).
Please see comments below.
Specific Comments:
Results and data, images and graphs "monitoring cell death in a well-established N. benthamiana assay": all images intended to show N. benthamiana cell death show autofluorescence, not cell death. Graphs are also based on autofluorescence induction and not cell death. Previous publications suggest that the cell death response is rather weak. Description of the experiments should be more accurate here and throughout the text. In addition, I do not understand why auto fluorescence is not quantified using image analysis software but interpreted visually (with a 6 level scoring system). Quantitative measurements combined with robust statistical analysis should be performed rather than scoring and qualitative interpretation of results from scoring. Differences between AVRPik-C and AVRPik-A on Pikp-1NK-KE on the photo in Figure 1B seems quite subtle in particular if one assumes that this is among the most differential images that were obtained.
The experiments are now described in more detail in the main text. As mentioned previously, our method for imaging fluorescence is not standardised in terms of light distribution, meaning it is not strictly useful for quantitative scoring. Having said this, we now include using estimation methods for analysis of this data in the revised manuscript. More importantly, as mentioned above, we have also performed ion leakage assays to add a quantitative measure of the cell death to our results, including statistical analysis.
Regarding the image in Figure 1B. We have not selectively picked one of the most differential images to show, it is just one example of the many repeats we have conducted to ensure we are reporting a robust phenotype. We feel it does show a clear difference in response of AVR-PikA vs. AVR-PikC.
Figure 1—figure supplement 2 -Statistical test to conclude on difference is missing. Unclear what is meant by "qualitatively stronger response".
Having performed the estimation methods analysis, we have revised our wording describing this data.
Figure 2A and 4B – Y2H experiments provide only weak support for the author's claims (in certain cases no support at all): growth differences between the different yeast clones are, in most cases, rather subtle and do often not correlate with the activity and in vitro binding data. E.g. in Figure 2A, Pikp-HMA/ARVPikD and Pikp-HMA/AVR-PikA grow identically and Pikp-HMA/AVR-PikE even grow better. Pikp-HMA/AVR-PikE growth is stronger than Pikp-HMA NK/KE/AVR-PikE. In Figure 4B, all clone seem to grow the same. Differences in color are not always easy to appreciate and rely on subjective interpretation. I would like to see growth on triple or quadruple dropout medium supplemented with different 3AT concentrations that accentuate growth differences and allow to unambiguously differentiate between growth strength of different yeast clones. If authors prefer to rely on β-Galactosidase activity, they should provide quantitative data from enzyme assays. Also apparent inconsistencies between Y2H and SPR or N. benthamiana cell death assays should be addressed in the text. For Pikp-HMA NK/KE, a negative control (unrelated effector) should be provided since it seems to interact with all AVR-Pik constructs that were tested including some that are expected to not interact or to interact only very weakly (e.g. AVRPikD-H46E). It should be indicated on the figure that HMAs are fused to BD and Avr's to AD.
Thanks for the suggestions. We have repeated our Y2H assays (with some modifications), and made several changes in the new figures:
1) We repeated the assays adding increasing concentrations of Aureobasidin A (we chose this antibiotic over 3AT because is more stringent and used with the Matchmaker Gold system). We hope you agree that this has improved the quality of the Y2H data (making it more robust), and that the results are clearer.
2) As suggested, we have added an unrelated rice blast effector, AVR-Pii, as a negative control in the Y2H.
3) We had indicated the fusion (AD or BD) in supplementary Western blots regarding yeast protein expression in the original submission. We have now also added this information to the relevant main manuscript figure legends.
4) As also requested by reviewer 2, we have added extra discussion to the text concerning the binding of AVR-PikC to PikpNK-KE as measured in vitro (SPR), when there is only weak binding as measured by Y2H and very limited interaction by co-IP in planta.
Figure 2C and Figure 4C: It would be interesting to see co-IP experiments with the isolated HMA domain (wt and mutant) to compare, on the one hand, to the other experiments with the isolated HMA domain (SPR and Y2H) and on the other to full-length Pikp1 variants.
We thought about such experiments in earlier stages of the project. However, we decided to focus on the full-length NLR when expressing in planta, as this is more similar to the native state of the interaction.
Figure 2B and Figure 2—figure supplements 2 and 3 – SPR experiments: Binding of C is similar to binding of A and E. Statistical analysis should be performed to determine if there are differences. Why do authors not calculate affinities from SPR and discuss differences in affinities? One question to address would be, can the weak differences in affinity between C and A (or E) explain recognition specificity?
As mentioned above, we have added extra discussion to the text concerning the binding of AVR-PikC to PikpNK-KE as measured in vitro (SPR), when there is only weak binding as measured by Y2H and very limited interaction by co-IP in planta. Statistical analysis for the SPR data has now been included and changes to the text made where relevant.
Determining absolute, accurate, binding affinities in SPR requires additional concentrations of the analyte, and saturated binding to be achieved. As for our previously published work (Nature Plants), we were most interested in comparative binding analysis between effector variants and mutants for which the Rmax analysis is sufficient. Further, we struggle to work with high (saturating) concentrations of some HMA domains in SPR, including Pikp-HMA and Pikp-HMANK-KE,as they can “stick” non-specifically to the chip surface and affect data quality at elevated concentrations. For these reasons we favour the Rmax analysis.
Figure 3—figure supplement 1: It looks to me as if the QtPISA interaction radar for Pikp-HMANK/KE/AVR-PikE does not support that binding is biologically relevant. Area is highly similar to the one for Pikp-HMA/AVR-PikE. Please comment on this point.
Having reviewed the interaction radars presented, we agree with the reviewer that these are not useful in these cases and we have removed them. The interaction radars are a simplistic view of the interfaces, and here they prove to be confusing.
Introduction, fourth paragraph: recognition does not necessarily rely on binding. Can be posttranslational modification (e.g. RRS1) or other mechanism.
We have added “or be modified by them” which covers this valid point.
Subsection “The engineered PikpNK-KE mutant expands association of full-length Pik-1 to effector variants in planta”: Replace "interact" by "associate with".
Replaced “interaction” by “association”, also in the subsection heading “The engineered PikpNK-KE mutant expands association of full-length Pik-1 to effector variants in planta”.
Subsection “The engineered PikpNK-KE mutant expands association of full-length Pik-1 to effector variants in planta”, last paragraph: Replace "binding" by "association".
Replaced.
Varden et al., 2019: This study is not peer-reviewed or properly published (Bioarchive paper) and therefore not be cited.
Citation updated as the paper is now accepted at JBC.
https://doi.org/10.7554/eLife.47713.036Article and author information
Author details
Funding
Biotechnology and Biological Sciences Research Council (BB/J004553)
- Kamoun S
- Banfield MJ
Biotechnology and Biological Sciences Research Council (BB/P012574)
- Kamoun S
- Banfield MJ
Biotechnology and Biological Sciences Research Council (BB/M02198X)
- Franceschetti M
- Kamoun S
- Banfield MJ
H2020 European Research Council (743165)
- Kamoun S
- Banfield MJ
John Innes Foundation
- De la Concepcion JC
- Franceschetti M
- Banfield MJ
Gatsby Charitable Foundation
- Sophien Kamoun
Japan Society for the Promotion of Science (15H05779)
- Ryohei Terauchi
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
This work was supported by the BBSRC (grants BB/J004553, BB/P012574, BB/M02198X), the ERC (proposal 743165), the John Innes Foundation, the Gatsby Charitable Foundation, and JSPS KAKENHI 15H05779. We thank the Diamond Light Source, UK (beamline i03 under proposal MX13467) for access to X-ray data collection facilities. We also thank David Lawson and Clare Stevenson (JIC X-ray Crystallography/Biophysical Analysis Platform) for help with protein-structure determination and SPR, and Andrew Davis and Phil Robinson (JIC Bioimaging facilities) for photography.
Senior Editor
- Detlef Weigel, Max Planck Institute for Developmental Biology, Germany
Reviewing Editor
- Thorsten Nürnberger, University of Tübingen, Germany
Reviewer
- Thorsten Nürnberger, University of Tübingen, Germany
Version history
- Received: April 17, 2019
- Accepted: September 17, 2019
- Accepted Manuscript published: September 19, 2019 (version 1)
- Version of Record published: September 30, 2019 (version 2)
Copyright
© 2019, De la Concepcion et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 4,627
- Page views
-
- 793
- Downloads
-
- 78
- Citations
Article citation count generated by polling the highest count across the following sources: Crossref, Scopus, PubMed Central.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Ecology
- Plant Biology
Global agro-biodiversity has resulted from processes of plant migration and agricultural adoption. Although critically affecting current diversity, crop diffusion from Classical antiquity to the Middle Ages is poorly researched, overshadowed by studies on that of prehistoric periods. A new archaeobotanical dataset from three Negev Highland desert sites demonstrates the first millennium CE&'s significance for long-term agricultural change in southwest Asia. This enables evaluation of the 'Islamic Green Revolution' (IGR) thesis compared to 'Roman Agricultural Diffusion' (RAD), and both versus crop diffusion during and since the Neolithic. Among the finds, some of the earliest aubergine (Solanum melongena) seeds in the Levant represent the proposed IGR. Several other identified economic plants, including two unprecedented in Levantine archaeobotany-jujube (Ziziphus jujuba/mauritiana) and white lupine (Lupinus albus)-implicate RAD as the greater force for crop migrations. Altogether the evidence supports a gradualist model for Holocene-wide crop diffusion, within which the first millennium CE contributed more to global agricultural diversity than any earlier period.
-
- Microbiology and Infectious Disease
- Plant Biology
Purinergic signaling activated by extracellular nucleotides and their derivative nucleosides trigger sophisticated signaling networks. The outcome of these pathways determine the capacity of the organism to survive under challenging conditions. Both extracellular ATP (eATP) and Adenosine (eAdo) act as primary messengers in mammals, essential for immunosuppressive responses. Despite the clear role of eATP as a plant damage-associated molecular pattern, the function of its nucleoside, eAdo, and of the eAdo/eATP balance in plant stress response remain to be fully elucidated. This is particularly relevant in the context of plant-microbe interaction, where the intruder manipulates the extracellular matrix. Here, we identify Ado as a main molecule secreted by the vascular fungus Fusarium oxysporum. We show that eAdo modulates the plant's susceptibility to fungal colonization by altering the eATP-mediated apoplastic pH homeostasis, an essential physiological player during the infection of this pathogen. Our work indicates that plant pathogens actively imbalance the apoplastic eAdo/eATP levels as a virulence mechanism.