Two NLR immune receptors acquired high-affinity binding to a fungal effector through convergent evolution of their integrated domain

Curation statements for this article:
  • Curated by eLife

    eLife logo

    Evaluation Summary:

    Convergent evolution is often observed in nature, but the molecular mechanisms allowing similar functions to independently emerge are rarely understood. This work determines how the high-affinity recognition of a pathogenic effector produced by the rice blast fungus, Avr-PikD, evolved in the immune receptor Pik-1. The integration of molecular evolution analyses with structure-function biochemical testing is novel to the field and the data quality is exceptional. In addition to advancing knowledge of host-microbe co-evolution, this work is exemplary in its transparency and the breadth of approaches utilized to understand protein evolution, and we expect that this study will provide a conceptual framework for similar studies in the future.

    (This preprint has been reviewed by eLife. We include the public reviews from the reviewers here; the authors also receive private feedback with suggested changes to the manuscript. Reviewer #2 agreed to share their name with the authors.)

This article has been Reviewed by the following groups

Read the full article

Abstract

A subset of plant NLR immune receptors carry unconventional integrated domains in addition to their canonical domain architecture. One example is rice Pik-1 that comprises an integrated heavy metal-associated (HMA) domain. Here, we reconstructed the evolutionary history of Pik-1 and its NLR partner, Pik-2, and tested hypotheses about adaptive evolution of the HMA domain. Phylogenetic analyses revealed that the HMA domain integrated into Pik-1 before Oryzinae speciation over 15 million years ago and has been under diversifying selection. Ancestral sequence reconstruction coupled with functional studies showed that two Pik-1 allelic variants independently evolved from a weakly binding ancestral state to high-affinity binding of the blast fungus effector AVR-PikD. We conclude that for most of its evolutionary history the Pik-1 HMA domain did not sense AVR-PikD, and that different Pik-1 receptors have recently evolved through distinct biochemical paths to produce similar phenotypic outcomes. These findings highlight the dynamic nature of the evolutionary mechanisms underpinning NLR adaptation to plant pathogens.

Article activity feed

  1. Reviewer #3 (Public Review):

    In this manuscript, the authors aim to elucidate the evolutionary history of the paired NLRs Pik-1/Pik-2 in rice. They ask two primary questions:

    (1) When (in evolutionary history) did the paired Pik-1/Pik-2 locus arise and when was the integrated domain integrated into the locus?

    (2) Has the binding affinity of the integrated domain changed over evolutionary time?

    The authors convincingly demonstrate that the integrated domain is undergoing positive selection, that its integration is ancient (~15MYA) and that inferred ancient alleles bind modern AVR-PikD with poor affinity. The subsequent biochemistry experiments and structural analyses identify which residues are important for interactions with AVR-PikD and which allelic combinations induce autoimmunity.

    The biochemical work, while interesting in and of itself for identifying the interacting residues and interactions between domains, was less informative about the evolution of the NLR-effector interaction, and most of the work did not advance our understanding of the questions listed above. The most emphasized biochemistry finding was that of reduced binding affinity of ancestral Pik-1 integrated domain. Specifically, the authors demonstrate that modern AVR-PikD has poor affinity with the ancient Pik-1 integrated domain. From this result the authors infer that ancestral Pik-1 likely bound a different effector. But it was not clear how the authors ruled out binding to an ancient AVR-PikD? I was confused as to why the authors excluded this possibility. Perhaps the authors contend that the absence of the Avr-PikD in other modern blast lineages indicates Avr-PikD is unique to modern rice-infecting M. oryzae. But this modern absence does not preclude Avr-PikD in the ancestral population. Furthermore, changes in binding over time would be the effective null hypothesis in the scenario of coevolving NLR and effector. Their finding seems consistent with expectations of coevolution, a phenomenon that has been widely reported in interactions between NLRs and effectors. The novelty in this manuscript stems from the synthesis of molecular evolution analysis with ancestral state reconstruction and testing.

    Overall this manuscript is exemplary in its integration of biochemical and evolutionary analyses to study plant-pathogen coevolution. While the findings are unsurprising, future emulation of this type of data integration will likely lead to significant insight into the coevolution of plants and their pathogens.

  2. Reviewer #2 (Public Review):

    In this study, Bialas et al. aimed at understanding the evolution of the diversity of Pik-1 immune receptors. First, using phylogenetic and selection analyses they determined that the Pik family of immune receptors is present in multiple grass species, with both Pik-1 and Pik-2 evolving before the radiation of the PACMAD and BOP clades. The author dated the insertion of an HMA domain in a Pik-1 subclade before the radiation of the Oryzinae and detected signs of positive selection on this domain. Using a combination of ancestral sequence reconstructions and biochemistry they determined that two of the extant Pik-1 haplotypes (Pikp-1 and Pikm-1) evolved independently the ability to associate at high affinity with the AVR-PikD effector following two different evolutionary paths. The authors determined that the increased binding correlates at least in one case with the improved ability to induce cell death when co-expressed in tobacco leaves with Pik-2 and AVR-PikD.

    Main strengths:

    The study combines a large diversity of methods to comprehensively address an important question. Despite the large amount of presented data (including a large number of variant names) it was a pleasure to read this very well structured manuscript. The work conducted here by the authors on the ancestral sequence reconstruction, the chimera and the biochemical assays (on two haplotypes!) is impressive and supports a very exciting conclusion. The presentation of all the experimental replicates as supplementary figure is a model of transparency and strengthen the conclusions.

    Weaknesses:

    The conclusions reached by the authors are mostly supported by the presented data, although there are a few points that need to be clarified. The Pik-1 phylogeny (Fig 1A): From the phylogenetic tree presented in Figure 1A it seems that Pik-1 experienced a duplication before the radiation of the BOP and PACMAD clades, with varying patterns of gene retention/loss (for instance loss of both copies in Brachypodium, loss in one clade for maize) and expansion (massive in wheat for instance in the clade where the fusion with the HMA domain did not occur, not in the other). I did not find this point discussed in the manuscript, although this could have an important impact. This would support the hypothesis that the HMA integration occurred before the radiation of the PACMAD clade. A better resolved phylogeny is needed to further test this possibility. In that context, the nomenclature should restrict the Pik-1 name to the actual orthologs, changing the number of Pik-1 per species (in panel 1D for instance).

    In Figure 4C and S13 the Pikp-1 variant I-N11 seems to associate more significantly with AVR-PikD than all the other variants, including I-N2 that was selected for the swap experiments. The reason why I-N2 was selected over other options (including I-N11) should be better explained.

    The correlation between evolution of high-affinity binding to AVR-PikD and the ability to induce immune response should be tested in reconstructed ancestral Pikm-1 variants. The presented data demonstrate nicely the gain of high-affinity binding in Pikm-1, but the impact this may have on the actual immunity function was not tested. It would be important to know whether additional mutations were required or not to turn the ancestral Pik1 into a functional Pikm-1 given that it is the basis for the model proposed in Figure 9. Alternatively, as the result of this experiment would not contradict the model even in absence of immune abilities (it would just add one extra step from high-affinity binding to immune function) the authors could propose this second evolutionary scenario as a supplementary figure.

    The nomenclature used for the Pik variants is not consistent throughout the manuscript, please homogenize as it is not always easy to follow.

    I am not familiar with the besthr R library used for the statistical analyses of the cell death assays, and I am not an expert in biochemistry (SPR, cristal structure) and cannot properly evaluate these aspects of the work.

  3. Reviewer #1 (Public Review):

    This paper was a pleasure to read. It is a tour-de-force study that is well-written, clear, and transparent. The study recounts how the HMA domain became integrated into the Pik NLRs and how it evolved higher affinity binding to a pathogen effector. Strikingly the authors demonstrate adaptability of distinct regions of the HMA:effector interface on two Pik NLRs, driving the convergent evolution of high-affinity binding to the effector. The study furthermore provides a framework for understanding protein evolution in the context of host-microbe interactions. The breadth and depth of the experiments that support the authors conclusions is extraordinary in my view.

  4. Evaluation Summary:

    Convergent evolution is often observed in nature, but the molecular mechanisms allowing similar functions to independently emerge are rarely understood. This work determines how the high-affinity recognition of a pathogenic effector produced by the rice blast fungus, Avr-PikD, evolved in the immune receptor Pik-1. The integration of molecular evolution analyses with structure-function biochemical testing is novel to the field and the data quality is exceptional. In addition to advancing knowledge of host-microbe co-evolution, this work is exemplary in its transparency and the breadth of approaches utilized to understand protein evolution, and we expect that this study will provide a conceptual framework for similar studies in the future.

    (This preprint has been reviewed by eLife. We include the public reviews from the reviewers here; the authors also receive private feedback with suggested changes to the manuscript. Reviewer #2 agreed to share their name with the authors.)

  5. This Zenodo record is a permanently preserved version of a PREreview. You can view the complete PREreview at https://prereview.org/reviews/4537865.

    Abstract

    This is a review of Białas et al. bioRxiv doi: https://doi.org/10.1101/2021.01.26.428286 posted on January 27, 2021. In this paper, the authors study the HMA domain from two NLR immune receptors. The authors identify a common ancestral version for these domains and highlight a convergent evolution that allows the interaction with a common effector.

     

    Summary

    HMA domains are integrated in some NLRs where they confer novel recognition specificities. The rice Pik-1 receptors, are such a NLR family, that recognize Avr-Pik effectors from Magnaporthe oryzae. Using novel Pik-1 orthologs, phylogenetic analyses, algorithms for ancestral sequence reconstruction and functional characterization of chimeric HMAs, this study examines the HMA molecular evolution of the allelic variants Pikp-1 and Pikm-1 to demonstrate that they convergently evolved to recognize Avr-PikD.

     

    Comments

     

    The immune profile in response to Avr-PikD for the majority of the novel Pik orthologs described in this manuscript is not known. However, I would suggest to indicate that Avr-PikD is recognized by rice Pik alleles Pikh, Pikp, Pikm, Piks and Pik* (Kanzaki et al., 2012; De la Conception et al., 2020). Hence, it could be discussed that, while Pikm and Pik* display the typical interface EMVKE, Piks, the closest ortholog of Pikm, is EMAKE and also recognize Avr-PikD.

     

    In this Figure 5A, the legend indicates "Historical mutations […] are shown next to the appropriate nodes." However, the node indicating the L221 mutation (IAQVV to LAQVV) is not at the correct position as both branches deriving from this node harbour the L221 mutation. The LAQVV indication rather describes (i) the branch before this node or (ii) the mutation that appeared at the previous ancient node.

     

    Also, while the historical order of the substitutions places Q228K before V229I, both substitutions are found in the 3 sequences at the bottom of the tree (O. rufipogon W2003, OsPikh-1 and OsPikp-1). Therefore, it appears incorrect to place the two individual substitutions on different nodes that discriminate O. rufipogon W2003 vs the two others.

     

    In Figure 9A, in MKANK-EMVKE (and for the 3 others combinations), it may be more judicious to replace the dash with a slash (as in the text), or maybe even better with an arrow to show the evolution, otherwise it could read like a first domain fused to a second one. The same could be applied to ANK-VKE and AV-VE. Also, as the reader constantly swap between Pikm-1 vs Pikp-1, and between interface 1 vs 2, and between ancHMA vs newHMAs, adding the legend "interface 1" or "interface 2" (as in Fig S26) on top of the 3D structures may also help a reader. Finally, I would also suggest adding, at least, the positions A260/N261/K262, E230, V261/K262/E263 and V231 on the protein models (as it is not so easy to find where are the two K262).

     

    Regarding the amino acid position, your text and description are right, but I had trouble to understand "In both cases, Lys-262 (K-262)" when you described "EMVKE and LKANK" since (i) the lysine residues are in different positions in these two domains and (ii) because the NK part of LKANK is not presented in Fig4A so I presumed that this last lysine would be K-263. I then noticed that the numbers above the ancHMA in Fig4A and Fig7A are different because they each correspond to the presented Pik-HMA (and Pikp and Pikm have a one residue shift (de la Conception et al., 2018)). Therefore, in Fig4A and Fig7A, I think that it would appear judicious to place the numbering below their respective Pik-HMA. In addition, while reporting "Lys-262 is structurally shifted" is correct, I think that it may be pertinent to present an alignment to support this description which explains the 'looping out' of Pikp-HMA.

     

    In Figure 8A, EMANK should be presented in purple for "EM" and green for "ANK" (as in 8B and 8D). As discussed for Figure 5A, EMANK describes the branch before the node. As for EMVKE, the position of the node is correct but few members are not EMVKE (likely because the phylogenetic tree is computed using the complete HMA domains, and not only these five mutations).

     

    In Figure S7, I am not sure if you performed some manual adjustments but I would suggest adjusting TraesCS7D02G007700.1 on the 5' HMA integration site. Some manual adjustments would definitely strengthen the percentage identity.

     

    For this same figure S7, the text refers to LpPik-1 but the figure is LPERR11G24730. It may be a good idea to indicate LpPik-1 (LPERR11G24730) in the text.

     

    I understand that the phylogenetic tree presented in Fig S9A is the annotated version used for Fig 5A and Fig 8A. While the Pik-1 nomenclature is presented in Table S7, I find difficult to find OsPik* in this Fig S9A which is quite important to navigate the Pik-1 phylogeny. Therefore, I would suggest to present OsPik* as Pik*_1_HM048900_1 rather than Pik_1_HM048900_1.

    The authors should pay attention to the OsPik* nomenclature which is presented with too many variations in the different Tables and Figures:

    It reads Pik* in TabS10

    It reads Pik*-1 in Fig2

    It reads Pik-1* in FigS26

    It reads OsPik* in Fig S6

    It reads OsPik-1 in Fig S11

    It reads OsPik*-1 in Fig S3

    It reads Pik_1_HM048900_1 in Fig S9

    It reads OsPik_1* HM048900_1 in Fig S12

     

    In Figure S9, ancHMA I-N2 appears to be the reference sequence. Therefore, in the antepenultimate position, OBART11G23150 should not present an N but a dot.

     

    In Figure S12B, I would suggest to also present the Pikm-1 HMA sequence as a last line in the alignment (even without its own data on probability for marginal reconstruction).

     

    In Figure S26, the top line presents the Pikm-1 EM and the second line the ancHMA MK, but going to the right ancHMA is now on line 1 with IA. Is there a reason to not present ancHMA on line 1?

     

    For homology modelling in Figure S33, which "failed" for MKANK, the authors present no obvious reason for using only Pikm-HMA–AVR-PikD as a model. It may be a good idea to consider using both Pikm-HMA–AVR-PikD (pdb:6FU9) and Pikp-HMA–AVR-PikD (pdb:6G10) as templates. And maybe even other Pik-HMA–AVR-Pik templates.

     

     

    Reviewer

    Freddy Boutrot, Anova-Plus, Evry, France.

     

    References

    De la Conception et al., 2018. Polymorphic residues in rice NLRs expand binding and response to effectors of the blast pathogen. Nature Plants. 4:576-585. https://doi.org/10.1038/s41477-018-0194-x

    De la Conception et al., 2020. The allelic rice immune receptor Pikh confers extended resistance to strains of the blast fungus through a single polymorphism in the effector binding interface. BioRxiv https://www.biorxiv.org/content/10.1101/2020.09.05.284240v1

    Kanzaki et al., 2012. Arms race co-evolution of Magnaporthe oryzae AVR-Pik and rice Pik genes driven by their physical interactions. The Plant Journal. 72:894-907.  https://doi.org/10.1111/j.1365-313X.2012.05110.x