An N-terminal motif in NLR immune receptors is functionally conserved across distantly related plant species

  1. Hiroaki Adachi
  2. Mauricio P Contreras
  3. Adeline Harant
  4. Chih-hang Wu
  5. Lida Derevnina
  6. Toshiyuki Sakai
  7. Cian Duggan
  8. Eleonora Moratto
  9. Tolga O Bozkurt
  10. Abbas Maqbool
  11. Joe Win
  12. Sophien Kamoun  Is a corresponding author
  1. University of East Anglia, Norwich Research Park, United Kingdom
  2. Imperial College London, United Kingdom
12 figures, 1 table and 3 additional files

Figures

Figure 1 with 1 supplement
Transposon-based truncation mutagenesis reveals a short 29 amino-acid region sufficient for NRC4-mediated cell death.

(A) Overview of the strategy for transposon-based C-terminal random truncation of NRC4 proteins. Hairpin Mu-STOP transposon and MuA proteins forming Mu transpososome were used for in vitro …

Figure 1—source data 1

Sequences of NRC4 truncation library.

The Mu-STOP transposon insertion sites were confirmed by PCR amplicon sequencing with Mu-STOP seq Rv primer. The 65 truncate sequences of NRC4 are listed in this file.

https://cdn.elifesciences.org/articles/49956/elife-49956-fig1-data1-v2.txt
Figure 1—figure supplement 1
Images of N. benthamiana leaves expressing truncated NRC4DV::Mu-STOP variants.

The images were taken 7 days after agroinfiltration. ‘DV’, ‘WT’ and ‘TL’ describe autoactive NRC4DV mutant, wild-type NRC4 and the truncation library, respectively. Red characters indicate clones …

Figure 2 with 2 supplements
NRC41-29-YFP induces cell death in Nicotiana benthamiana independently of endogenous NRC4.

(A) Schematic representation of wild-type NRC4-YFP (NRC4WT-YFP) and the variants used for the in planta expression assays. The colour code is: red represents NRC4 1 to 29 amino acid region. (B) NRC41…

Figure 2—figure supplement 1
Knocking out of NRC4a and NRC4b in Nicotiana benthamiana impairs Rpi-blb2-mediated HR cell death.

(A) Schematic representation of sgRNA positions targeting NRC4a and NRC4b. The PAM motifs are marked in blue, and the sequences of sgRNAs are marked in red. (B) Genotyping results of selected T2 nrc4…

Figure 2—figure supplement 2
NRC41-29-YFP cell death is compromised by YFP A206K mutation.

(A) Schematic representation of NRC4DV-YFP, NRC41-29-YFP and the variants used for the in planta expression assays. Arrowheads show A206K mutation site in YFP. The red colour represents NRC4 1 to 29 …

Figure 3 with 2 supplements
NRC4 carries N-terminal sequences that are conserved across distantly related CC-NLRs.

(A) Schematic representation of the different NLR domains used in TRIBE-MCL and phylogenetic analyses. (B) Distribution of plant NLRs across N-terminal domain tribes. The colour codes are: orange …

Figure 3—source data 1

Amino acid sequences of full-length NLRs in the CC-NLR database.

988 NLR sequences used for HMMER analysis are listed.

https://cdn.elifesciences.org/articles/49956/elife-49956-fig3-data1-v2.txt
Figure 3—source data 2

Amino acid sequences of N-terminal domains in the CC-NLR database.

N-terminal domain sequences of 988 proteins used for Tribe-MCL analysis are listed.

https://cdn.elifesciences.org/articles/49956/elife-49956-fig3-data2-v2.txt
Figure 3—source data 3

N-terminal domain tribes of CC-NLRs.

Results of the Tribe-MCL analysis are included in this file.

https://cdn.elifesciences.org/articles/49956/elife-49956-fig3-data3-v2.xlsx
Figure 3—figure supplement 1
Phylogenetic analysis of NLR proteins from dicot and monocot plant species.

NLR proteins were predicted by NLR-parser from N. benthamiana (NbS-), tomato (Solyc-), Arabidopsis (AT-), sugar beet (Bv-), rice (Os-) and barley (HORVU-) proteomes, and were used for the MAFFT …

Figure 3—figure supplement 1—source data 1

Amino acid sequences for CC/TIR-NLR phylogenetic tree.

NB-ARC domain sequences used for phylogenetic analysis are shown with the IDs, N. benthamiana (NbS-), tomato (Solyc-), Arabidopsis (AT-), sugar beet (Bv-), rice (Os-) and barley (HORVU-).

https://cdn.elifesciences.org/articles/49956/elife-49956-fig3-figsupp1-data1-v2.txt
Figure 3—figure supplement 1—source data 2

CC/TIR-NLR phylogenetic tree file.

The phylogenetic tree was saved in newick file format.

https://cdn.elifesciences.org/articles/49956/elife-49956-fig3-figsupp1-data2-v2.txt
Figure 3—figure supplement 2
Phylogenetic analysis of CC-NLR proteins from dicot and monocot plant species.

The phylogenetic tree was constructed with the NB-ARC domain sequences of CC-NLRs as described in Figure 3—figure supplement 2—source data 2. Each leaf is labelled with different colour ranges …

Figure 3—figure supplement 2—source data 1

Amino acid sequences for CC-NLR phylogenetic tree.

NB-ARC domain sequences used for phylogenetic analysis are shown with the IDs, N. benthamiana (NbS-), tomato (Solyc-), Arabidopsis (AT-), sugar beet (Bv-), rice (Os-) and barley (HORVU-).

https://cdn.elifesciences.org/articles/49956/elife-49956-fig3-figsupp2-data1-v2.txt
Figure 3—figure supplement 2—source data 2

CC-NLR phylogenetic tree file.

The phylogenetic tree was saved in newick file format.

https://cdn.elifesciences.org/articles/49956/elife-49956-fig3-figsupp2-data2-v2.txt
Figure 4 with 2 supplements
The MADA motif is a conserved unit at the very N-terminus of NRC4 and ZAR1.

(A) Schematic representation of a classical CC-NLR protein highlighting the position of the MADA motif. Consensus sequence pattern of the MADA motif identified by MEME along with an alignment of …

Figure 4—source data 1

Output of the HMMER search using the MADA motif HMM against tomato and Arabidopsis proteomes.

HMM scores are listed with the IDs, tomato (Solyc-) and Arabidopsis (AT-), and annotation information.

https://cdn.elifesciences.org/articles/49956/elife-49956-fig4-data1-v2.xlsx
Figure 4—source data 2

Amino acid sequences of the MADA motif.

The sequences were extracted from MEME output against N-terminal domain Tribe 2 and were used to build the MADA motif HMM.

https://cdn.elifesciences.org/articles/49956/elife-49956-fig4-data2-v2.xls
Figure 4—figure supplement 1
CC-NLRs have conserved protein sequence patterns in the beginning of the N-terminal domains.

Consensus sequence patterns in N-terminal domains were identified by MEME from 226 Tribe 1, 102 Tribe 2, 83 Tribe 3 and 59 Tribe 4 members. Motif logos describe the N-terminal consensus patterns …

Figure 4—figure supplement 2
N terminus of NRC4 possesses a consensus pattern coined MADA motif.

(A) Consensus sequence of the MADA motif. The MADA motif logo was generated by MEME from 87 N-terminal domains of Tribe 2 members. (B) Graphical representation of the MADA HMM used to screen …

Figure 5 with 1 supplement
The MADA motif is conserved in ~20% of CC-NLRs.

(A) Schematic representation of a classical CC-NLR protein highlighting the regions used for HMMER searches (MADA-HMM) and for TRIBE-MCL. (B) Occurrence of MADA/MADAL-CC-NLRs in representative …

Figure 5—source data 1

Output of the HMMER search using the MADA motif HMM against the CC-NLR database.

HMM scores of the predicted MADA motifs are listed by IDs, N. benthamiana (NbS-), tomato (Solyc-), Arabidopsis (AT-), sugar beet (Bv-), rice (Os-) and barley (HORVU-), with Tribe-MCL result, the start (‘MADA_strat’) and end (‘MADA_end’) positions of the MADA motifs in the CC-NLRs.

https://cdn.elifesciences.org/articles/49956/elife-49956-fig5-data1-v2.xlsx
Figure 5—source data 2

List of the predicted Arabidopsis MADA-CC-NLRs.

The IDs are listed with the HMM score.

https://cdn.elifesciences.org/articles/49956/elife-49956-fig5-data2-v2.xlsx
Figure 5—figure supplement 1
Bar graph of MADA/MADAL-CC-NLRs according to HMM score.

(A) HMM score bar graph for CC-NLR database (988 proteins). MADA/MADAL-CC-NLRs from HMMER analysis are shown in orange and blue, respectively. (B) HMM score bar graph with plant species information. …

NRC-dependent sensors (NRC-S) do not have the MADA motif.

(A) Distribution of NRCs (NRC-H) and NRC-dependent sensors (NRC-S) across N-terminal domain tribes of CC-NLRs. Individual NLR members of the NRC superclade were classified based on phylogenetic …

Figure 6—source data 1

HMM scores of NRC-superclade proteins.

HMM scores are listed by IDs, N. benthamiana (NbS-), tomato (Solyc-) and sugar beet (Bv-) with Tribe-MCL result, the start (‘MADA_strat’) position of the MADA motifs and NRC clade information (‘NRC-H’ and ‘NRC-S’).

https://cdn.elifesciences.org/articles/49956/elife-49956-fig6-data1-v2.xlsx
Figure 7 with 1 supplement
L9, V10 and L14 triple mutation impairs cell death activity of autoimmune NRC4DV.

(A) Schematic representation of NRC4 and the mutated sites in the MADA motif. Mutated sites and substituted residues are shown as red characters in the NRC4 sequence alignment. (B) Cell death …

Figure 7—figure supplement 1
NRC41-29-YFP cell death is compromised by L9, V10 and L14 triple mutation.

(A) Schematic representation of NRC41-29-YFP, ZAR11-144-YFP and the variants used for the in planta expression assays. Arrowheads show triple alanine mutation sites in NRC41-29 and ZAR11-144. The …

Figure 8 with 2 supplements
L9E, L13E and L17E single mutations impair cell death activity of autoimmune NRC4DV.

(A) Schematic representation of NRC4 and the glutamic acid (E) mutant scan of the MADA motif. Mutated sites are shown as red characters in the NRC4 sequence. (B) Cell death observed in N. benthamiana

Figure 8—figure supplement 1
Alanine mutants do not compromise HR cell death triggered by autoactive NRC4.

(A) Schematic representation of NRC4 and the alanine (A) mutant scan of the MADA motif. Mutated sites are shown as red characters in the NRC4 sequence. (B) Cell death observed in N. benthamiana

Figure 8—figure supplement 2
Mapping loss of function mutations on N-terminal α helices of NRC4.

(A) Cartoon representation of N-terminal α helices of NRC4 resistosome (zoom in grey box of Figure 4B). (B, C) N-terminal α helices are rotated 90 degrees and mutated amino acids are shown as stick …

NRC41-29-YFP forms MADA- and YFP-dependent puncta.

(A) Subcellular localization of NRC41-29-YFP and the mutant proteins in N. benthamiana epidermal cells. N. benthamiana leaves expressing YFP, YFPA206K, NRC41-29-YFP, NRC41-29-YFPA206K and NRC41-29L9E

Figure 10 with 1 supplement
First 17 amino acids of NRC4 can be functionally replaced by the N-terminus of other MADA/MADAL-CC-NLRs.

(A) Alignment of the N-terminal region of the MADA/MADAL-CC-NLRs. Key residues for cell death activity identified in Figure 8 are marked as red characters with asterisks in the sequence alignment. …

Figure 10—figure supplement 1
Quantification of cell death triggered by NRC4 MADA motif chimeras.

(A) Box plots showing cell death intensity scored as an HR index based on three independent experiments. Statistical differences among the samples were analysed with Tukey’s HSD test (p<0.01). (B) …

The chimeric protein ZAR11-17-NRC4 complements NRC4 function in Rpi-blb2-mediated resistance.

(A) Schematic representation of NRC4 complementation assay for Rpi-blb2-mediated resistance. Wild-type and the variants of NRC4 were co-expressed with RFP-Rpiblb2 in wild-type or nrc4a/b_9.1.3 N. …

Evolution of NLRs from singletons to networks.

We propose that the N-terminal MADA motif/α1 helix has emerged early in the evolution of CC-NLRs and has remained constrained throughout time as singletons evolved from multifunctional proteins into …

Tables

Key resources table
Reagent type (species) or resourceDesignationSource or referenceIdentifiersAdditional information
Genetic reagent (Nicotiana benthamiana)NRC4-KO N. benthamiana
(nrc4a/b_9.1.3 and nrc4a/b_1.2.1)
This paperMaterials and methods: Generation of N. benthamiana
nrc4a/b CRISPR/Cas9
mutants
Recombinant
DNA reagent
pGEM::Mu-STOPThis paperMaterials and methods:
Mu-STOP in vitro
transposition
Commercial assay, kitMutation Generation
System Kit
Thermo FisherCat #: F-701Materials and methods:
Mu-STOP in vitro
transposition
Gene
(Solanum lycopersicum)
Tomato
genome
sequence
(Tomato ITAG release 2.40)
Sol Genomics
Network
(https://solgenomics.net/)
Materials and methods: Bioinformatic and phylogenetic analyses
Gene
(N. benthamiana)
N. benthamiana genome sequence (N. benthamiana Genome v0.4.4)Sol Genomics
Network (https://solgenomics.net/)
Materials and methods: Bioinformatic and phylogenetic analyses
Gene (Arabidopsis thaliana)Arabidopsis
genome
sequence (Araport11)
https://www.araport.org/Materials and methods: Bioinformatic and
phylogenetic analyses
Gene
(Beta vulgaris)
Sugar beet
genome sequence
(RefBeet-1.2)
http://bvseq.molgen.mpg.de/index.shtmlMaterials and methods: Bioinformatic and phylogenetic analyses
Gene
(Oryza sativa)
Rice genome
sequence
(Rice Gene Models
in Release 7)
http://rice.plantbiology.msu.edu/Materials and methods: Bioinformatic and
phylogenetic analyses
Gene
(Hordeum vulgare)
Barley genome
sequence (IBSC_v2)
https://www.barleygenome.org.uk/Materials and methods: Bioinformatic and
phylogenetic analyses
Other3D structure
of ZAR1
Protein Data Bank6J5TMaterials and methods:
Structure homology modelling

Additional files

Supplementary file 1

Primers used for generating NRC4 variants by Golden Gate cloning.

https://cdn.elifesciences.org/articles/49956/elife-49956-supp1-v2.pptx
Supplementary file 2

The MADA-HMM for HMMER analysis.

This MADA-HMM was used for searching MADA-CC-NLRs from CC-NLR database (Figure 3—source data 1).

https://cdn.elifesciences.org/articles/49956/elife-49956-supp2-v2.hmm
Transparent reporting form
https://cdn.elifesciences.org/articles/49956/elife-49956-transrepform-v2.pdf

Download links