Herbivorous insects independently evolved salivary effectors to regulate plant immunity by destabilizing the malectin-LRR RLP NtRLP4

  1. Xin Wang
  2. Jia-Bao Lu
  3. Yi-Zhe Wang
  4. Xu-Hong Zhou
  5. Jian-Ping Chen
  6. Chuan-Xi Zhang
  7. Jun-Min Li  Is a corresponding author
  8. Hai-Jian Huang  Is a corresponding author
  1. State Key Laboratory for Quality and Safety of Agro-Products, Key Laboratory of Biotechnology in Plant Protection of MARA, Zhejiang Key Laboratory of Green Plant Protection, Institute of Plant Virology, Ningbo University, China
  2. Institute of Applied Ecology, Fujian Agriculture and Forestry University, China
  3. Office of Science and Technology, Yunnan University of Chinese Medicine, China
6 figures and 2 additional files

Figures

Figure 1 with 4 supplements
BtRDP is a salivary protein and secreted into plants.

(a) Detection of BtRDP in Nicotiana tabacum plants. The untreated and Bemisia tabaci-infested tobacco plants, as well as the salivary gland samples, are collected for western blotting assays. Asterisk indicates non-specific binding. Rubisco staining (RbcL) is conducted to visualize the amount of sample loading. (b, c) Detection of BtRDP in different tissues by western blotting and quantitative real-time PCR (qRT-PCR) assays. Ov, ovary; SG, salivary gland; Ca, carcass; FB, fat body. Silver staining and anti-actin serum are used to visualize the sample loading. (d) Expression patterns of BtRDP in different development stages. Ny, nymph; Ps, pseudopupa; Fe, female; Ma, male. B. tabaci tubulin is used as an internal control. The relative quantitative method (2−ΔΔCt) is used to evaluate the quantitative variation. Data are presented as mean values ± SEM (n = 3 independent biological replicates). Different lowercase letters indicate statistically significant differences at the p < 0.05 level according to one-way ANOVA test followed by Tukey’s multiple comparisons test. Immunohistochemical staining of BtRDP in salivary glands (e) and salivary sheath secreted from B. tabaci (f). The salivary gland and its nearby tissues are dissected and incubated with anti-BtRDP serum or pre-immune serum conjugated with Alexa Fluor 488 NHS Ester (green) and actin dye phalloidin–rhodamine (red). The nucleus is stained with 4′,6-diamidino-2-phenylindole (DAPI, blue). PSG, principal salivary gland; ASG, accessory salivary gland. The lower images represent the enlarged images of the boxed area in the upper image. Experiments are repeated thrice for (a, b), while twice for (e, f). Similar results are observed, and representative images are displayed.

Figure 1—source data 1

PowerPoint file containing original western blots for Figure 1a, indicating the relevant bands and treatments.

https://cdn.elifesciences.org/articles/108737/elife-108737-fig1-data1-v1.zip
Figure 1—source data 2

Original files for western blot analysis displayed in Figure 1a.

https://cdn.elifesciences.org/articles/108737/elife-108737-fig1-data2-v1.zip
Figure 1—source data 3

PowerPoint file containing original western blots for Figure 1b, indicating the relevant bands and treatments.

https://cdn.elifesciences.org/articles/108737/elife-108737-fig1-data3-v1.zip
Figure 1—source data 4

Original files for western blot analysis displayed in Figure 1b.

https://cdn.elifesciences.org/articles/108737/elife-108737-fig1-data4-v1.zip
Figure 1—figure supplement 1
Characteristic of BtRDP.

(a) Deduced amino acid sequence of BtRDP. Arrow indicates the signal peptide cleavage site. (b–g) Mass spectrums of the identified unique peptides in watery saliva of Bemisia tabaci.

Figure 1—figure supplement 2
Analysis of insect RDP.

(a) Phylogenetic tree of RDP sequences from Aleyrodidae species. The trees are constructed by RAxML v0.9.0 using the maximum likelihood method with 1000 bootstrap replicates. Nodes with bootstrap values greater than 50 are displayed. (b) Sequence alignments of RDPs. The amino acid sequence of RDPs used for phylogenetic tree construction is aligned using ClusterX software. Black shades indicate the conserved regions.

Figure 1—figure supplement 3
Expression patterns of Bemisia tabaci BtRDP.

Transcripts per million (TPM) expression values of BtRDP in different tissues (a) and at different developmental stages (b) are determined based on the transcriptomic data. Data are presented as mean values ± SEM (n = 3 independent biological replicates).

Figure 1—figure supplement 4
Expression patterns of BtActin and Bt18s rRNA.

The transcript level of each gene in different tissues (a) and at different developmental stages (b) is quantified by quantitative real-time PCR (qRT-PCR). Bemisia tabaci tubulin is used as an internal control. The relative quantitative method (2−ΔΔCt) is used to evaluate the quantitative variation. Ov, ovary; SG, salivary gland; Ca, carcass; FB, fat body; Ny, nymph; Ps, pseudopupa; Fe, female; Ma, male. Data are presented as mean values ± SEM (n = 3 independent biological replicates). Different lowercase letters indicate statistically significant differences at p < 0.05 level according to one-way ANOVA test followed by Tukey’s multiple comparisons test.

Figure 2 with 5 supplements
Effects of BtRDP on Bemisia tabaci.

(a) Treating B. tabaci with dsBtRDP significantly reduces the transcript and protein level of target gene. The dsGFP-treated B. tabaci is used as a control. Effects of BtRDP knockdown on insect reproduction (b) and feeding behavior (c). Electrical penetration graph (EPG) is used to monitor the insect feeding behavior, which can be classified into nonpenetration (np), pathway duration (C), phloem salivation (E1), and phloem ingestion (E2) phases. All EPG recordings are performed for 8 hr. Typical EPG waveforms are displayed in Figure 2—figure supplement 3. (d–f) Insect performance on transgenic plants overexpressing BtRDP (oeBtRDP). (d) Detection of BtRDP in transgenic Nicotiana tabacum overexpressing a complete coding region of BtRDP. The empty vector (EV) plants are used as the control. The flag tag is fused to the C-terminal ends of recombinant proteins. (e) Comparison of insect reproduction on oeBtRDP#1 and EV plants. Five B. tabaci individuals are confined to indicated plants for 3 days, and the oviposited eggs are counted. (f) Attraction of oeBtRDP#1 and EV leaves to B. tabaci in a two-choice equipment. A group of 40 female B. tabaci are released into a device containing oeBtRDP#1 and EV leaves. The number of insects settling on each leaf is counted. After 48 hr, the number of eggs on each leaf is counted. (g–i) Insect performance on transgenic plants overexpressing BtRDP-sp (oeBtRDP-sp). (g) Detection of BtRDP-sp in transgenic N. tabacum overexpressing BtRDP without a signal peptide. The insect reproduction (h) and settlement (i) on oeBtRDP-sp transgenic plants are recorded. The EPG data are first checked for normality and homogeneity of variance, and data not fitting a normal distribution are subjected to log10 transformation. Data are presented as mean ± SEM. Statistical significance in panels (ac), (e), (f), (h), and (i) was assessed using two-tailed unpaired Student’s t-tests. ***p < 0.001; **p < 0.01; *p < 0.05; ns, not significant. Western blotting assays are repeated thrice with similar results.

Figure 2—source data 1

PowerPoint file containing original western blots for Figure 2a, indicating the relevant bands and treatments.

https://cdn.elifesciences.org/articles/108737/elife-108737-fig2-data1-v1.zip
Figure 2—source data 2

Original files for western blot analysis displayed in Figure 2a.

https://cdn.elifesciences.org/articles/108737/elife-108737-fig2-data2-v1.zip
Figure 2—source data 3

PowerPoint file containing original western blots for Figure 2d, indicating the relevant bands and treatments.

https://cdn.elifesciences.org/articles/108737/elife-108737-fig2-data3-v1.zip
Figure 2—source data 4

Original files for western blot analysis displayed in Figure 2d.

https://cdn.elifesciences.org/articles/108737/elife-108737-fig2-data4-v1.zip
Figure 2—source data 5

PowerPoint file containing original western blots for Figure 2g, indicating the relevant bands and treatments.

https://cdn.elifesciences.org/articles/108737/elife-108737-fig2-data5-v1.zip
Figure 2—source data 6

Original files for western blot analysis displayed in Figure 2g.

https://cdn.elifesciences.org/articles/108737/elife-108737-fig2-data6-v1.zip
Figure 2—figure supplement 1
BtRDP is efficiently and specifically suppressed by dsBtRDP.

(a) Relative transcript level of BtRDP, BtFTSP, BtE3, and Bt56 after dsRNA treatments. Bemisia tabaci is treated with dsGFP and dsBtRDP. The expression pattern of each gene is determined 4 days post-treatment using quantitative real-time PCR (qRT-PCR) methods. B. tabaci tubulin is used as an internal control. The relative quantitative method (2−ΔΔCt) is used to evaluate the quantitative variation. Data are presented as mean values ± SEM (n = 3 independent biological replicates). (b) Protein level of BtRDP after dsRNA treatments. Actin is used as an internal control. Two independent biological replicates are displayed. (c) Immunohistochemical staining of BtRDP after dsRNA treatment. The salivary gland (SG) and its nearby tissues are collected from dsGFP- and dsBtRDP-treated B. tabaci. The samples are incubated with anti-BtRDP serum conjugated with Alexa Fluor 488 NHS Ester (green) and actin dye phalloidin–rhodamine (red), and examined by Leica SP8. The nucleus is stained with 4′,6-diamidino-2-phenylindole (DAPI, blue). The lower image in each treatment represents the enlarged images of the boxed area in the upper image. Experiments are repeated twice, with each inspecting more than 10 salivary glands.

Figure 2—figure supplement 1—source data 1

PowerPoint file containing original western blots for Figure 2—figure supplement 1b, indicating the relevant bands and treatments.

https://cdn.elifesciences.org/articles/108737/elife-108737-fig2-figsupp1-data1-v1.zip
Figure 2—figure supplement 1—source data 2

Original files for western blot analysis displayed in Figure 2—figure supplement 1b.

https://cdn.elifesciences.org/articles/108737/elife-108737-fig2-figsupp1-data2-v1.zip
Figure 2—figure supplement 2
Effects of dsRNA treatment on insect survivorship and salivary sheath formation.

(a, b) Effects of dsRNA treatment on salivary sheath formation. Newly emerged Bemisia tabaci adults are injected with dsGFP and dsBtRDP. The dsRNA-treated B. tabaci are fed on artificial diets. The salivary sheaths left on parafilm are inspected by scanning electron microscopy (a), and the length of salivary sheath is measured from the top to base of salivary sheath (b). Data are presented as mean values ± SEM. Twenty salivary sheaths from each treatment are measured. p-value is determined by two-tailed unpaired Student’s t-test. *p < 0.05. (c) Effects of dsRNA treatment on insect survivorship. A group of 20–30 B. tabaci are placed in a leaf cage, and their mortality is recorded for 10 consecutive days. Three independent biological replications are performed. Differences in survivorship between the two treatments are tested by log-rank test. ns, not significant.

Figure 2—figure supplement 3
Typical electrical penetration graph (EPG) waveforms for Bemisia tabaci feeding on Nicotiana tabacum.

The B. tabaci feeding behavior can be classified into nonpenetration (np), pathway duration (C), phloem salivation (E1), and phloem ingestion (E2).

Figure 2—figure supplement 4
Effects of empty vector (EV) transgenic plants on Bemisia tabaci.

A two-choice experiment is performed to investigate the attraction of wild type (WT) and EV transgenic plants to B. tabaci. A group of 40 female B. tabaci are released into a device containing EV and WT leaves. The number of insects settling on each leaf is counted at 3, 6, 12, 24, 36, and 48 hr (a). After 48 hr, the number of eggs on each leaf is counted (b). Data are presented as mean values ± SEM. Ten independent biological replicates are performed. p-values are determined by two-tailed unpaired Student’s t-test. ns, not significant.

Figure 2—figure supplement 5
Effects of BtRDP and BtRDP-sp overexpression on Bemisia tabaci.

Insect performance on oeBtRDP#2 (a, b) and oeBtRDP-sp#2 (c, d) transgenic plants is tested. (a, c) Comparison of insect reproduction on empty vector (EV) and oeBtRDP#2/oeBtRDP-sp#2 transgenic plants. Five B. tabaci individuals are confined to indicated plants for 3 days, and the oviposited eggs are counted. Twelve independent biological replicates are performed. (b, d) Attraction of EV and oeBtRDP#2/oeBtRDP-sp#2 leaves to B. tabaci in a two-choice equipment. A group of 40 female B. tabaci are released into a device containing indicated leaves. The number of insects settling on each leaf is counted at 3, 6, 12, 24, 36, and 48 hr. After 48 hr, the number of eggs on each leaf is counted. Data are presented as mean values ± SEM. Ten independent biological replicates are performed. p-values are determined by two-tailed unpaired Student’s t-test. ***p < 0.001.

Figure 3 with 4 supplements
NtRLP4 interacts with BtRDP and confers plant resistance to Bemisia tabaci.

(a) Domain organization of NtRLP4. NtRLP4 contains a predicted N-terminal signal peptide (SP), a malectin-like domain, an LRR domain, and a transmembrane domain (TM). (b) The protein level of NtRLP4 in response to B. tabaci infestation. The untreated Nicotiana tabacum is used as a negative control. Four independent biological replicates are performed, and the representative images are displayed. The band density is measured using ImageJ. The density values from four biological replicates are calculated and the mean value in the controls is set at 1.0. (c) Yeast two-hybrid assays showing the interaction between BtRDP and NtRLP4. The different combinations of constructs are transformed into yeast cells and are grown on the selective medium SD/-Trp/-Leu (DDO), and the interactions are tested with SD/-Trp/-Leu/-His/-Ade (QDO). (d) Co-immunoprecipitation assay showing the interaction between BtRDP and NtRLP4. Total proteins are extracted from N. benthamiana leaves transiently co-expressing NtRLP4-myc/NtCf9-myc with BtFTSP-flag/BtRDP-flag. All genes are expressed with a complete coding region, and the myc/flag tags are fused at the C-terminus. Precipitation is performed using flag beads. The samples are probed with anti-flag and anti-myc antibodies for immunoblotting analysis. (e–h) Analysis of transgenic N. tabacum overexpressing NtRLP4 (oeRLP). (e) Detection of NtRLP4 level in oeRLP plants. The empty vector (EV) plant is used as a control. Two independent oeRLP lines are selected. The samples are probed with an anti-NtRLP4 antibody. Rubisco staining (RbcL) is used to visualize the amount of sample loading. (f) Attraction of oeRLP#1 and EV leaves to B. tabaci in a two-choice equipment. A group of 40 female B. tabaci are released into a device containing oeRLP#1 and EV leaves. The number of insects settling on each leaf is counted at each time point. After 48 hr, the number of eggs on each leaf is counted. (g) Comparison of insect reproduction on oeRLP#1 and EV plants. Five B. tabaci individuals are confined to indicated plants for 3 days, and the oviposited eggs are counted. (h) Relative transcript levels of salicylic acid (SA)- and jasmonic acid (JA)-associated genes in oeRLP#1 and EV plants. PAL, phenylalanine ammonia lyase; NPR1, nonexpressor of pathogenesis-related protein 1; FAD7, fatty acid desaturase 7; PDF1.2, plant defensin 1.2. Two independent biological replicates are performed in (d) and (e). Data in (fh) are presented as mean values ± SEM. For insect bioassays in (f) and (g), n = 10 independent biological replicates. For quantitative real-time PCR (qRT-PCR) in (h), n = 3 independent biological replicates. p-values are determined by two-tailed unpaired Student’s t-test. ***p < 0.001; **p < 0.01; *p < 0.05; ns, not significant.

Figure 3—source data 1

PowerPoint file containing original western blots for Figure 3b, indicating the relevant bands and treatments.

https://cdn.elifesciences.org/articles/108737/elife-108737-fig3-data1-v1.zip
Figure 3—source data 2

Original files for western blot analysis displayed in Figure 3b.

https://cdn.elifesciences.org/articles/108737/elife-108737-fig3-data2-v1.zip
Figure 3—source data 3

PowerPoint file containing original western blots for Figure 3d, indicating the relevant bands and treatments.

https://cdn.elifesciences.org/articles/108737/elife-108737-fig3-data3-v1.zip
Figure 3—source data 4

Original files for western blot analysis displayed in Figure 3d.

https://cdn.elifesciences.org/articles/108737/elife-108737-fig3-data4-v1.zip
Figure 3—source data 5

PowerPoint file containing original western blots for Figure 3e, indicating the relevant bands and treatments.

https://cdn.elifesciences.org/articles/108737/elife-108737-fig3-data5-v1.zip
Figure 3—source data 6

Original files for western blot analysis displayed in Figure 3e.

https://cdn.elifesciences.org/articles/108737/elife-108737-fig3-data6-v1.zip
Figure 3—figure supplement 1
Expression patterns of NtRLP4 in response to Bemisia tabaci infestation.

Relative transcript level of NtRLP4 in response to B. tabaci infestation is determined by quantitative real-time PCR (qRT-PCR). Data are presented as mean values ± SEM. Different lowercase letters indicate statistically significant differences at p < 0.05 level according to one-way ANOVA test followed by Tukey’s multiple comparisons test. Three independent biological replicates are performed.

Figure 3—figure supplement 2
Sequence alignments of SlCf9 and NtCf9.

The amino acid sequence of SlCf9 and NtCf9 is aligned using ClusterX software. Black shades indicate the conserved regions.

Figure 3—figure supplement 3
BtRDP interacts with NtRLP4.

(a) Yeast two-hybrid assays showing the interaction between NtRLP4 and different whitefly salivary proteins (BtFTSP-sp, BtSP16.3-sp, and BtSP37.4-sp, without signal peptides), as well as BtRDP-sp and truncated NtRLP4. The different combinations of constructs are transformed into yeast cells and are grown on the selective medium SD/-Trp/-Leu (DDO), and the interactions are tested with SD/-Trp/-Leu/-His/-Ade (QDO). (b) Expression of three control salivary proteins at the protein level. The target bands are indicated by asterisks. (c) Bimolecular fluorescence complementation assay showing the specific interaction between BtRDP and NtRLP4. YFP fluorescence is observed when co-expressing the N-terminal nYFP tag fused NtRLP4 and N-terminal cYFP tag fused BtRDP-sp. Bar = 40 μm. (d) Co-immunoprecipitation (Co-IP) assays on oeBtRDP transgenic plants. oeBtRDP transgenic plants that overexpress BtRDP-flag are incubated with anti-flag beads. Endogenous NtRLP4 can be immunoprecipitated by BtRDP-flag in oeBtRDP transgenic plants, while the flag tag alone in EV plants fails to immunoprecipitate endogenous NtRLP4. Rubisco staining (RbcL) is used to visualize the amount of sample loading.

Figure 3—figure supplement 3—source data 1

PowerPoint file containing original western blots for Figure 3—figure supplement 3b, indicating the relevant bands and treatments.

https://cdn.elifesciences.org/articles/108737/elife-108737-fig3-figsupp3-data1-v1.zip
Figure 3—figure supplement 3—source data 2

Original files for western blot analysis displayed in Figure 3—figure supplement 3b.

https://cdn.elifesciences.org/articles/108737/elife-108737-fig3-figsupp3-data2-v1.zip
Figure 3—figure supplement 3—source data 3

PowerPoint file containing original western blots for Figure 3—figure supplement 3d, indicating the relevant bands and treatments.

https://cdn.elifesciences.org/articles/108737/elife-108737-fig3-figsupp3-data3-v1.zip
Figure 3—figure supplement 3—source data 4

Original files for western blot analysis displayed in Figure 3—figure supplement 3d.

https://cdn.elifesciences.org/articles/108737/elife-108737-fig3-figsupp3-data4-v1.zip
Figure 3—figure supplement 4
NtSOBIR1 interacts with NtRLP4 but not BtRDP.

(a) Domain organization of NtSOBIR1. NtSOBIR1 contains a predicted N-terminal signal peptide (SP), a leucine-rich repeat (LRR) domain, a transmembrane (TM) domain, and a kinase domain. (b) Co-immunoprecipitation assay showing the interaction between NtSOBIR1 and NtRLP4. (c) Co-immunoprecipitation assay showing the interaction between NtRLP4 and BtRDP, but not NtSOBIR1 and BtRDP. The complete coding region of NtSOBIR1 and NtRLP4 is fused with flag and myc tags at the C-terminal ends. Total proteins are extracted from N. benthamiana leaves co-expressing NtRLP4-myc with GFP-flag or NtSOBIR1-flag. Precipitation is performed using flag beads. The samples are probed with anti-flag and anti-myc antibodies for immunoblot analysis. Experiments are repeated twice with similar results.

Figure 3—figure supplement 4—source data 1

PowerPoint file containing original western blots for Figure 3—figure supplement 4b, indicating the relevant bands and treatments.

https://cdn.elifesciences.org/articles/108737/elife-108737-fig3-figsupp4-data1-v1.zip
Figure 3—figure supplement 4—source data 2

Original files for western blot analysis displayed in Figure 3—figure supplement 4b.

https://cdn.elifesciences.org/articles/108737/elife-108737-fig3-figsupp4-data2-v1.zip
Figure 3—figure supplement 4—source data 3

PowerPoint file containing original western blots for Figure 3—figure supplement 4d, indicating the relevant bands and treatments.

https://cdn.elifesciences.org/articles/108737/elife-108737-fig3-figsupp4-data3-v1.zip
Figure 3—figure supplement 4—source data 4

Original files for western blot analysis displayed in Figure 3—figure supplement 4d.

https://cdn.elifesciences.org/articles/108737/elife-108737-fig3-figsupp4-data4-v1.zip
Figure 4 with 10 supplements
BtRDP suppresses plant defenses by promoting NtRLP4 degradation.

(a) Effects of dsBtRDP suppression on insect reproduction (n = 10 independent biological replicates) when feeding on empty vector (EV) and NtRLP4-silenced (RNAi-RLP) transgenic Nicotiana tabacum. (b) Transient overexpressing BtRDP-mCherry attenuates H2O2 accumulation caused by NtRLP4-GFP overexpression. N. tabacum co-expressing RFP-mCherry and NtRLP4-GFP is used as a control. All genes are expressed with a complete coding region. The experiment is repeated five times with similar results. (c) Relative transcript level of salicylic acid (SA)- and jasmonic acid (JA)-associated genes in BtRDP-mCherry/NtRLP4-GFP and RFP-mCherry/NtRLP4-GFP plants (n = 3 independent biological replicates). PAL, phenylalanine ammonia lyase; NPR1, nonexpressor of pathogenesis-related protein 1; FAD7, fatty acid desaturase 7; PDF1.2, plant defensin 1.2. (d) Comparison of insect reproduction (n = 16 independent biological replicates) on transient-expressed BtRDP-mCherry/NtRLP4-GFP and RFP-mCherry/NtRLP4-GFP plants. (e) Attraction of BtRDP-mCherry/NtRLP4-GFP and RFP-mCherry/NtRLP4-GFP leaves to B. tabaci in a two-choice equipment (n = 16 independent biological replicates). The number of eggs on each leaf is counted at 48 hr post insect release. (f) Effects of BtRDP and BtRDP-sp overexpression on NtRLP4 accumulation. Transgenic oeBtRDP (left) and oeBtRDP-sp (right) plants are probed with anti-NtRLP4 and anti-flag antibodies for immunoblotting analysis. The experiments are repeated twice with similar results. (g) Degradation of NtRLP4 by BtRDP in N. benthamiana leaves. NtRLP4-myc and BtRDP-flag are transiently co-expressed in N. benthamiana plants through Agrobacterium infiltration. Agrobacterium carrying NtCf9-myc and GUS-flag is used as negative controls. (h) Effects of 26S proteasome inhibitor (MG132) on NtRLP4 accumulation. Co-infiltrated leaves are treated with MG132 at 24 hr post-injection. The samples are probed with anti-flag and anti-myc antibodies for immunoblot analysis. Rubisco staining (RbcL) is conducted to visualize the amount of sample loading. The small triangle indicates the different concentrations of Agrobacterium (OD600 = 0.05, 0.3, and 1.0). (i) NtRLP4 is ubiquitinated in planta. NtRLP4-myc is co-expressed transiently with HA-UBQ in N. benthamiana leaves. Extracted total proteins are immunoprecipitated by anti-myc beads and immunoblotted with anti-myc or anti-HA antibody. Experiments are repeated thrice for (g) and (h), while twice for (i). Band density is measured using ImageJ. The density values from three biological replicates are calculated and the mean value in the first lane is set at 1.0. Data are presented as mean values ± SEM. p-values are determined by two-tailed unpaired Student’s t-test. ***p < 0.001; **p < 0.01; *p < 0.05; ns, not significant.

Figure 4—source data 1

PowerPoint file containing original western blots for Figure 4f, indicating the relevant bands and treatments.

https://cdn.elifesciences.org/articles/108737/elife-108737-fig4-data1-v1.zip
Figure 4—source data 2

Original files for western blot analysis displayed in Figure 4f.

https://cdn.elifesciences.org/articles/108737/elife-108737-fig4-data2-v1.zip
Figure 4—source data 3

PowerPoint file containing original western blots for Figure 4g, indicating the relevant bands and treatments.

https://cdn.elifesciences.org/articles/108737/elife-108737-fig4-data3-v1.zip
Figure 4—source data 4

Original files for western blot analysis displayed in Figure 4g.

https://cdn.elifesciences.org/articles/108737/elife-108737-fig4-data4-v1.zip
Figure 4—source data 5

PowerPoint file containing original western blots for Figure 4h, indicating the relevant bands and treatments.

https://cdn.elifesciences.org/articles/108737/elife-108737-fig4-data5-v1.zip
Figure 4—source data 6

Original files for western blot analysis displayed in Figure 4h.

https://cdn.elifesciences.org/articles/108737/elife-108737-fig4-data6-v1.zip
Figure 4—source data 7

PowerPoint file containing original western blots for Figure 4i, indicating the relevant bands and treatments.

https://cdn.elifesciences.org/articles/108737/elife-108737-fig4-data7-v1.zip
Figure 4—source data 8

Original files for western blot analysis displayed in Figure 4i.

https://cdn.elifesciences.org/articles/108737/elife-108737-fig4-data8-v1.zip
Figure 4—figure supplement 1
Effect of NtRLP4 overexpression on Bemisia tabaci performance.

(a) Attraction of empty vector (EV) and oeNtRLP4#2 (oeRLP#2) transgenic Nicotiana tabacum to B. tabaci in a two-choice equipment. A group of 40 female B. tabaci are released into a device containing oeRLP#2 and EV leaves. The number of insects settling on each leaf is counted at 3, 6, 12, 24, 36, and 48 hr. After 48 hr, the number of eggs on each leaf is counted. (b) Comparison of insect reproduction on EV and oeRLP#2 transgenic plants. Five B. tabaci individuals are confined to indicated plants for 3 days, and the oviposited eggs are counted. (c) Attraction of leaves transiently overexpressing GFP and NtRLP4-GFP to B. tabaci in a two-choice equipment. (d) Comparison of insect reproduction on N. tabacum leaves transiently overexpressing GFP and NtRLP4-GFP. (e) Relative transcript level of salicylic acid (SA)- and jasmonic acid (JA)-associated genes in N. tabacum leaves transiently overexpressing GFP and NtRLP4-GFP. PAL, phenylalanine ammonia lyase; NPR1, nonexpressor of pathogenesis-related protein 1; FAD7, fatty acid desaturase 7; PDF1.2, plant defensin 1.2. (f) Transiently overexpressing NtRLP4-GFP induces H2O2 accumulation. For bioassays in (ad), 10 independent biological replicates are performed. For qPCR analysis in (e), three independent biological replicates are performed. Data are presented as mean values ± SEM. p-values are determined by two-tailed unpaired Student’s t-test. ***p < 0.001. The experiment in (f) is repeated five times with similar results.

Figure 4—figure supplement 2
Transcriptomic comparison of empty vector (EV) and oeNtRLP4#1 (oeRLP#1) transgenic plants.

(a) Principal component analysis (PCA) of gene expression patterns in EV and oeRLP#1 transgenic Nicotiana tabacum plants. The first two principal components (PC1 and PC2) based on transcriptomic results are shown. Four independent biological replicates are performed. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis of differentially expressed genes (DEGs) that significantly upregulated (b) and downregulated (c). Enriched p-values are calculated according to one-sided hypergeometric test using TBtools software. Expression patterns of DEGs that are annotated as RLK/RLP (d) or WRKY transcription factor (e). Gene expression patterns are illustrated by a heatmap. The max transcripts per million value of each gene is set as 100.

Figure 4—figure supplement 3
Effect of NtSOBIR1 overexpression and silencing on Bemisia tabaci performance.

(a) Validation of NtSOBIR1 transient overexpression using quantitative real-time PCR (qRT-PCR) (n = 4 biological replicates). (b) NtSOBIR1 transient overexpression induces cell death phenotype 5 days post agro-injection. (c) Comparison of insect reproduction on Nicotiana tabacum transiently overexpressing GFP and NtSOBIR1-GFP. Five B. tabaci individuals are confined to indicated plants for 3 days, and the oviposited eggs are counted. (d) Attraction of GFP and NtSOBIR1-GFP leaves to B. tabaci in a two-choice equipment. A group of 40 female B. tabaci are released into a device containing GFP and NtSOBIR1-GFP leaves. The number of insects settling on each leaf is counted at 3, 6, 12, 24, 36, and 48 hr (left). After 48 hr, the number of eggs on each leaf is counted (right). (e) Silencing efficiency of hairpin NtSOBIR1 using qRT-PCR (n = 3 independent biological replicates). (f) Transiently silencing NtSOBIR1 affects growth of infiltrated leaves. (g) Comparison of insect reproduction on Nicotiana tabacum leaves infiltrated with hairpin empty vector (EV) and NtSOBIR1. (h) Attraction of hairpin EV and NtSOBIR1-treated leaves to B. tabaci in a two-choice equipment. The experiments in (b) and (f) are repeated 10 times with similar results. In (c), (d), (g), and (h), 10 independent biological replicates are performed. p-values are determined by two-tailed unpaired Student’s t-test. ***p < 0.001; **p < 0.01.

Figure 4—figure supplement 4
Effects of RLP silencing on Bemisia tabaci performance.

(a) Detection of NtRLP4 protein in empty vector (EV) and NtRLP4-silenced (RNAi-RLP) transgenic Nicotiana tabacum. Rubisco staining (RbcL) is conducted to visualize the amount of sample loading. (b, c) Attraction of EV and RNAi-RLP plants to B. tabaci in a two-choice equipment. Two NtRLP4-silenced plants (RNAi-RLP#1 in b; RNAi-RLP#2 in c) are assayed. A group of 40 female B. tabaci are released into a device containing two leaves. The number of insects settling on each leaf is counted at 3, 6, 12, 24, 36, and 48 hr. After 48 hr, the number of eggs on each leaf is counted. Ten independent biological replicates are performed. Data are presented as mean values ± SEM. p-values are determined by two-tailed unpaired Student’s t-test. ***p < 0.001.

Figure 4—figure supplement 4—source data 1

PowerPoint file containing original western blots for Figure 4—figure supplement 4a, indicating the relevant bands and treatments.

https://cdn.elifesciences.org/articles/108737/elife-108737-fig4-figsupp4-data1-v1.zip
Figure 4—figure supplement 4—source data 2

Original files for western blot analysis displayed in Figure 4—figure supplement 4a.

https://cdn.elifesciences.org/articles/108737/elife-108737-fig4-figsupp4-data2-v1.zip
Figure 4—figure supplement 5
Influence of BtRDP on NtRLP4 accumulation by fluorescent analysis.

The NtRLP4-GFP is transiently co-expressed with red fluorescent protein (RFP)-mCherry, BtFTSP-mCherry, or BtRDP-mCherry via agroinfiltration. The samples are imaged by confocal microscopy at 48 hr post-injection. RFP-mCherry and BtFTSP-mCherry are used as negative controls. Three independent biological replicates are performed, and three representative images are taken in each biological replicate. Fluorescence intensity is measured using ImageJ. The intensity values from three biological replicates are calculated and the mean value in the upper left image is set at 1.0. Representative fluorescence images are displayed. The small triangle indicates the different concentrations (OD600 = 0.1 and 1.0) of Agrobacterium. Bar = 40 μm.

Figure 4—figure supplement 6
Effects of BtRDP on suppressing NtRLP4-associated plant defenses.

(a) Comparison of insect reproduction on Nicotiana tabacum plants transiently overexpressing GFP along and BtRDP-mCherry/NtRLP4-GFP. Five Bemisia tabaci individuals are confined to indicated plants for 3 days, and the oviposited eggs are counted. (b) Attraction of GFP- and BtRDP-mCherry/NtRLP4-GFP-expressed leaves to B. tabaci in a two-choice equipment. A group of 40 female B. tabaci are released into a device containing two leaves. The number of insects settling on each leaf is counted at 3, 6, 12, 24, 36, and 48 hr. After 48 hr, the number of eggs on each leaf is counted. Data are presented as mean values ± SEM. Ten independent biological replicates are performed. p-values are determined by two-tailed unpaired Student’s t-test. *p < 0.05; ns, not significant.

Figure 4—figure supplement 7
Degradation of NtRLP4 by purified BtRDP-sp in Nicotiana benthamiana leaves.

Recombinant BtRDP-sp-his and GFP-his proteins are expressed in Escherichia coli and purified by Ni-NTA. N. benthamiana plants overexpressing NtRLP4-myc or NtCf9-myc are then infiltrated with different concentrations of purified BtRDP-sp-his and GFP-his. The samples are probed with anti-myc antibodies for immunoblot analysis. Rubisco staining (RbcL) is conducted to visualize the amount of sample loading. Coomassie brilliant blue (CBB) staining is conducted to visualize the amount of recombinant BtRDP-sp-his and GFP-his proteins. Experiments are repeated three times with similar results. Band density is measured using ImageJ. The density values from three biological replicates are calculated and the mean value in the first lane is set at 1.0. Data are presented as mean values ± SEM.

Figure 4—figure supplement 7—source data 1

PowerPoint file containing original western blots for Figure 4—figure supplement 7, indicating the relevant bands and treatments.

https://cdn.elifesciences.org/articles/108737/elife-108737-fig4-figsupp7-data1-v1.zip
Figure 4—figure supplement 7—source data 2

Original files for western blot analysis displayed in Figure 4—figure supplement 7.

https://cdn.elifesciences.org/articles/108737/elife-108737-fig4-figsupp7-data2-v1.zip
Figure 4—figure supplement 8
Effect of BtRDP on NtRLP4 and NtSOBIR1.

NtRLP4-myc and NtSOBIR1-flag were transiently co-expressed with different concentrations of BtRDP-flag in Nicotiana benthamiana plants through Agrobacterium infiltration. The samples are probed with anti-flag and anti-myc antibodies for immunoblot analysis. Rubisco staining (RbcL) is conducted to visualize the amount of sample loading. The amount of BtRDP was controlled by infiltrating different concentrations of Agrobacterium (OD600 = 0.1, 0.3, and 1.0).

Figure 4—figure supplement 8—source data 1

PowerPoint file containing original western blots for Figure 4—figure supplement 8, indicating the relevant bands and treatments.

https://cdn.elifesciences.org/articles/108737/elife-108737-fig4-figsupp8-data1-v1.zip
Figure 4—figure supplement 8—source data 2

Original files for western blot analysis displayed in Figure 4—figure supplement 8.

https://cdn.elifesciences.org/articles/108737/elife-108737-fig4-figsupp8-data2-v1.zip
Figure 4—figure supplement 9
Effect of BtRDP on the transcript level of NtRLP4.

NtRLP4-myc is agroinfiltrated together with different concentrations of BtRDP-flag or GFP-flag into Nicotiana benthamiana leaves. The transcript level of NtRLP4 is determined by quantitative real-time PCR (qRT-PCR). Data are presented as mean values ± SEM. Three independent biological replicates are performed. The same lowercase letters ‘a’ indicate no statistically significant differences at p < 0.05 according to the one-way ANOVA test followed by Tukey’s multiple comparisons test.

Figure 4—figure supplement 10
Effects of autophagy inhibitor on NtRLP4 accumulation.

NtRLP4-myc and BtRDP-flag were transiently co-expressed in Nicotiana benthamiana plants through Agrobacterium infiltration. Co-infiltrated leaves are treated with autophagy inhibitor bafilomycin A1 (BAF) and 3-Methyladenine (3-MA) at 24 hr post-injection. The samples are probed with anti-flag and anti-myc antibodies for immunoblot analysis. Rubisco staining (RbcL) is conducted to visualize the amount of sample loading. The small triangle indicates the different concentrations of Agrobacterium (OD600 = 0.1 and 1.0). Experiments are repeated three times with similar results. Band density is measured using ImageJ. The density values from three biological replicates are calculated and the mean value in the first lane is set at 1.0. Data are presented as mean values ± SEM. p-values are determined by two-tailed unpaired Student’s t-test. ***p < 0.001; **p < 0.01; *p < 0.05; ns, not significant.

Figure 4—figure supplement 10—source data 1

PowerPoint file containing original western blots for Figure 4—figure supplement 10, indicating the relevant bands and treatments.

https://cdn.elifesciences.org/articles/108737/elife-108737-fig4-figsupp10-data1-v1.zip
Figure 4—figure supplement 10—source data 2

Original files for western blot analysis displayed in Figure 4—figure supplement 10.

https://cdn.elifesciences.org/articles/108737/elife-108737-fig4-figsupp10-data2-v1.zip
Figure 5 with 5 supplements
Rice RLP4 is targeted by salivary protein NlSP104 from Nilaparvata lugens.

(a) Domain organization of Oryza sativa RLP4 (OsRLP4). (b, c) Yeast two-hybrid and Co-IP assays showing the interaction between OsRLP4 and NlSP104. All genes are expressed with a complete coding region in Co-IP assays, while NlSP104 without a signal peptide is used in yeast two-hybrid assays. Experiments in (b) are repeated twice with the similar result. (d) The expression patterns of NlSP104 in different N. lugens tissues. Te, testis; Ov, ovary; SG, salivary gland; Ca, carcass; FB, fat body. Data are presented as mean values ± SEM (n = 3 independent biological replicates). Different lowercase letters indicate statistically significant differences at p < 0.05 level according to one-way ANOVA test followed by Tukey’s multiple comparisons test. Effects of dsRNA treatment on insect survivorship (e), reproduction (f), and honeydew excretion (g). For survivorship analysis, a group of 30 N. lugens are reared in a cage. Three independent biological replications are performed. Differences in survivorship between the two treatments are tested by log-rank test. ns, not significant. For reproduction analysis, n = 10 and 13 individuals are tested in dsGFP- and dsNlSP104-treatment, respectively. The sterile females are excluded from data analysis. For honeydew analysis, n = 10 independent biological replicates. p-values in (f) and (g) are determined by two-tailed unpaired Student’s t-test. ***p < 0.001; **p < 0.01. (h) Effect of NlSP104 on the accumulation of transient-expressed OsRLP4. OsRLP4-myc is agro-injected into Nicotiana benthamiana together with different concentrations of NlSP104-flag and NlSP7-flag. NtCf9-myc is used as a negative control. The small triangle indicates the different concentrations (OD600 = 0.05, 0.3, and 1.0) of Agrobacterium. Rubisco staining (RbcL) is conducted to visualize the amount of sample loading. Experiments are repeated three times with similar results. Band density is measured using ImageJ. The density values from three biological replicates are calculated with the mean value in the first lane being set at 1.0.

Figure 5—source data 1

PowerPoint file containing original western blots for Figure 5b, indicating the relevant bands and treatments.

https://cdn.elifesciences.org/articles/108737/elife-108737-fig5-data1-v1.zip
Figure 5—source data 2

Original files for western blot analysis displayed in Figure 5b.

https://cdn.elifesciences.org/articles/108737/elife-108737-fig5-data2-v1.zip
Figure 5—source data 3

PowerPoint file containing original western blots for Figure 5h, indicating the relevant bands and treatments.

https://cdn.elifesciences.org/articles/108737/elife-108737-fig5-data3-v1.zip
Figure 5—source data 4

Original files for western blot analysis displayed in Figure 5h.

https://cdn.elifesciences.org/articles/108737/elife-108737-fig5-data4-v1.zip
Figure 5—figure supplement 1
Sequence alignment and phylogenetic tree of RLP4 homologs.

(a) The amino acid sequence of RLP4s from Nicotiana tabacum, N. benthamiana, Solanum lycopersicum, and Oryza sativa are aligned using ClusterX software. Black shades indicate the conserved regions. (b) Phylogenetic analysis of RLP4 and Cf9-associated genes used in this study. The unrooted phylogenetic trees are constructed by RAxML v0.9.0 using the maximum likelihood method with 1000 bootstrap replicates.

Figure 5—figure supplement 2
Influence of BtRDP on RLP4 homology in Solanum lycopersicum.

(a) Domain organization of S. lycopersicum RLP4 (SlRLP4). SlRLP4 contains a predicted N-terminal signal peptide (SP), a malectin-like domain, an LRR domain, and a transmembrane (TM) domain. (b, c) Yeast two-hybrid and co-immunoprecipitation (Co-IP) assays showing the interaction between BtRDP and SlRLP4. In (b), BtRDP is expressed without a signal peptide (BtRDP-sp), while SlRLP4 and SlCf9 are expressed without a signal peptide and transmembrane domain (SlRLP4(23–546) and SlCf9(24–524)). (d) Effect of BtRDP on the accumulation of SlRLP4. SlRLP4-myc and SlCf9-myc are agro-injected together with different concentrations of BtRDP-flag or GFP-flag. In (c, d), the complete coding region of BtRDP, SlRLP4, and SfCf9 is fused with flag or myc tags at C-terminal ends, respectively. The small triangle indicates the different concentrations (OD600 = 0.05, 0.3, and 1.0) of Agrobacterium. Rubisco staining (RbcL) is conducted to visualize the amount of sample loading. Experiments are repeated three times with similar results. Band density is measured using ImageJ. The density values from three biological replicates are calculated. For SlRLP4-myc bands (lanes 1–6), lane 1 is set at 1.0. For SlCf9-myc bands (lanes 7–9), lane 7 is set at 1.0.

Figure 5—figure supplement 2—source data 1

PowerPoint file containing original western blots for Figure 4—figure supplement 7c, indicating the relevant bands and treatments.

https://cdn.elifesciences.org/articles/108737/elife-108737-fig5-figsupp2-data1-v1.zip
Figure 5—figure supplement 2—source data 2

Original files for western blot analysis displayed in Figure 4—figure supplement 7c.

https://cdn.elifesciences.org/articles/108737/elife-108737-fig5-figsupp2-data2-v1.zip
Figure 5—figure supplement 2—source data 3

PowerPoint file containing original western blots for Figure 4—figure supplement 7d, indicating the relevant bands and treatments.

https://cdn.elifesciences.org/articles/108737/elife-108737-fig5-figsupp2-data3-v1.zip
Figure 5—figure supplement 2—source data 4

Original files for western blot analysis displayed in Figure 4—figure supplement 7d.

https://cdn.elifesciences.org/articles/108737/elife-108737-fig5-figsupp2-data4-v1.zip
Figure 5—figure supplement 3
Yeast two-hybrid assays showing the interaction between OsRLP4 and salivary proteins from planthopper species.

OsRLP4(29–551) (OsRLP4 without signal peptides and transmembrane domains) is fused to pGADT7 vector, while salivary proteins without signal peptides are fused to pGBKT7 vector, respectively. The different combinations of constructs are transformed into yeast cells and are grown on the selective medium SD/-Trp/-Leu (DDO), and the interactions are tested with SD/-Trp/-Leu/-His/-Ade (QDO). Twenty salivary proteins are selected. The accession number of each protein is listed as follows: (1) MF278694.1, (2) XP_039291719.1, (3) KU365967.1, (4) MF278706.1, (5) MF278711.1, (6) MF278714.1, (7) MF278715.1, (8) MF278720.1, (9) XP_022195702.1, (10) XP_022196818.2, (11) XP_039284081.1, (12) XP_022192221.2, (13) XP_022207944.2, (14) KT764973, (15) RZF44823.1, (16) RZF42644.1, (17) RZF33006.1, (18) RZF48570.1, (19) RZF42817.1, and (20) RZF33751.1.

Figure 5—figure supplement 4
Three-dimensional structure of BtRDP and NlSP104.

AlphaFold2 is used to predict the protein structure. The structures are colored as a spectrum from N-terminus (blue) to C-terminus (red).

Figure 5—figure supplement 5
Expression patterns of Nilaparvata lugens NlSP104.

Transcripts per million (TPM) expression values of NlSP104 in different tissues are determined based on the transcriptomic data. The head that contains salivary glands exhibits the highest NlSP104 expression in male and female N. lugens. Data are presented as mean values ± SEM (n = 3 independent biological replicates).

The proposed model for the suppression of receptor-like protein (RLP)-mediated plant defenses by salivary effectors.

Host plants employ pattern-recognition receptors (PRRs) to detect various damage-associated molecular patterns (DAMPs) and herbivore-associated molecular patterns (HAMPs) triggered by insect feeding. The RLP4/SOBIR1 complex plays a vital role in initiating pattern-triggered immunity (PTI), including H2O2 burst, upregulation of jasmonic acid (JA), and downregulation of salicylic acid (SA), which hinders insect feeding. The whitefly Bemisia tabaci and planthopper Nilaparvata lugens independently evolved salivary proteins that targeted plant RLP4. B. tabaci salivary sheath protein BtRDP interacts with the leucine-rich repeat (LRR) domain of RLP4 from Nicotiana tabaci and Solanum lycopersicum, while N. lugens NlSP104 targets both the LRR domain and the malectin-like domain of Oryza sativa RLP4. These interactions promote the ubiquitin-dependent degradation of RLP4, thereby disrupting the stability of the RLP4/SOBIR1 complex. The presence of salivary effectors causes a hormonal shift and suppresses the H2O2 burst, finally favoring insect feeding.

Additional files

Supplementary file 1

Sequences, accessions, and expression patterns of associated genes in this study.

(A) Identification of RDP and SP101 homologs in insect species. (B) Proteins from a Nicotiana benthamiana cDNA library screened by yeast two-hybrid using BtRDP as a bait. (C) Differentially expressed genes between empty vector (EV) and oeRLP#1 transgenic plant. (D) Primers used in this study.

https://cdn.elifesciences.org/articles/108737/elife-108737-supp1-v1.docx
MDAR checklist
https://cdn.elifesciences.org/articles/108737/elife-108737-mdarchecklist1-v1.pdf

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

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)

  1. Xin Wang
  2. Jia-Bao Lu
  3. Yi-Zhe Wang
  4. Xu-Hong Zhou
  5. Jian-Ping Chen
  6. Chuan-Xi Zhang
  7. Jun-Min Li
  8. Hai-Jian Huang
(2026)
Herbivorous insects independently evolved salivary effectors to regulate plant immunity by destabilizing the malectin-LRR RLP NtRLP4
eLife 14:RP108737.
https://doi.org/10.7554/eLife.108737.4