Yaoyao Jiang1,2
Weiran Zheng3,4Jing Li3
Peng Liu2Kaili Zhong2
Peng Jin4Miaoze Xu2
Jian Yang2Jianping Chen1,2,3*- 1College of Plant Protection, Fujian Agriculture and Forestry University, Fuzhou, China
- 2State Key Laboratory for Managing Biotic and Chemical Treats to the Quality and Safety of Agro-Products, Institute of Plant Virology, Ningbo University, Ningbo, China
- 3State Key Laboratory Breeding Base for Zhejiang Sustainable Pest and Disease Control, Zhejiang Provincial Key Laboratory of Plant Virology, Institute of Virology and Biotechnology, Zhejiang Academy of Agricultural Sciences, Hangzhou, China
- 4College of Plant Protection, Hunan Agricultural University, Changsha, China
WRKY transcription factors play important roles in plants, including responses to stress; however, our understanding of the function of WRKY genes in plant responses to viral infection remains limited. In this study, we investigate the role of NbWRKY40 in Nicotiana benthamiana resistance to tomato mosaic virus (ToMV). NbWRKY40 is significantly downregulated by ToMV infection, and subcellular localization analysis indicates that NbWRKY40 is targeted to the nucleus. In addition, NbWRKY40 activates W-box-dependent transcription in plants and shows transcriptional activation in yeast cells. Overexpressing NbWRKY40 (OEWRKY40) inhibits ToMV infection, whereas NbWRKY40 silencing confers susceptibility. The level of salicylic acid (SA) is significantly higher in OEWRKY40 plants compared with that of wild-type plants. In addition, transcript levels of the SA-biosynthesis gene (ICS1) and SA-signaling genes (PR1b and PR2) are dramatically higher in OEWRKY40 plants than in the control but lower in NbWRKY40-silenced plants than in the control. Furthermore, electrophoretic mobility shift assays show that NbWRKY40 can bind the W-box element of ICS1. Callose staining reveals that the plasmodesmata is decreased in OEWRKY40 plants but increased in NbWRKY40-silenced plants. Exogenous application of SA also reduces viral accumulation in NbWRKY40-silenced plants infected with ToMV. RT-qPCR indicates that NbWRKY40 does not affect the replication of ToMV in protoplasts. Collectively, our findings suggest that NbWRKY40 likely regulates anti-ToMV resistance by regulating the expression of SA, resulting in the deposition of callose at the neck of plasmodesmata, which inhibits viral movement.
Introduction
Plants, as sessile organisms, are often seriously affected by biotic stressors and, thus, have evolved a wide array of sophisticated mechanisms to mitigate the detrimental effects of these stressors, many of which involve transcription factors (TFs) (Miller et al., 2017). Indeed, studies show that TFs are essential components of a plant’s stress response, and they function by binding specific cis-acting elements to regulate the expression of genes containing such elements in their promoters (Jakoby et al., 2002; Pandey and Somssich, 2009; Tiwari et al., 2012; Li J. et al., 2019).
The WRKY TFs are plant-specific transcriptional regulators involved in a variety of signaling pathways (Rushton et al., 2010; Song et al., 2016). WRKY TFs that contain WRKY domains and zinc-finger motifs can be categorized into five groups: I, IIa–IIb, IIc, IId–IIe, and III (Liang et al., 2017). All WRKY factors contain a WRKYGQK sequence motif and exhibit a high binding affinity to W-box sequences [TTGAC(C/T)], which are found upstream of many defense-related genes (Rushton et al., 2010) as well in the promoters of WRKY TFs themselves. Therefore, WRKY TFs possess regulatory, auto-regulatory, and cross-regulatory properties (Eulgem et al., 1999). Indeed, since the first WRKY gene SPF1 was identified in sweet potato (Ishiguro and Nakamura, 1994), many WRKY TFs have been reported to be transcriptional regulators that are involved in complex and interconnected transcriptional networks in plants (Rushton et al., 2010).
Tremendous progress has been made in elucidating the function of WRKY proteins associated with abiotic and biotic stress responses. AtWRKY25, AtWRKY26, OsWRKY11, and OsWRKY30 are associated with heat-, salt-, and drought-tolerance (Wu et al., 2009; Li et al., 2011; Shen et al., 2012; Scarpeci et al., 2013). OsWRKY42 is associated with high temperature- and salinity-tolerance (Pillai et al., 2018), and OsWRKY76 functions as a negative regulator in responses to cold stress and blast disease (Yokotani et al., 2013). AtWRKY70 is reported to function at the intersection of the salicylic acid (SA)- and jasmonic acid (JA)-mediated defense-signaling pathways during resistance against specific bacterial and fungal pathogens (Li et al., 2004). In addition, AtWRKY46 is reported to improve basal resistance against Pseudomonas syringae through its overlapping and synergetic functions with AtWRKY70 and AtWRKY53 (Hu et al., 2012). Many WRKY proteins are also involved in plant responses to viral infection. For example, CaWRKYd, a Capsicum annuum gene, is reported to contribute to tobacco mosaic virus (TMV)-mediated hypersensitive-response (HR) apoptosis by regulating downstream gene expression (Huh et al., 2012). However, βC1 proteins encoded by begomoviruses have been shown to bind to WRKY20, interfering with the biosynthesis and accumulation of glucosinolates in plants as well as several defensive factors. This not only benefits the begomovirus, but also whiteflies (i.e., the vector of begomoviruses), and deters non-vector insects (Zhao et al., 2019).
Several WRKY genes are reported to participate in various plant hormone-mediated signaling pathways (Ishihama and Yoshioka, 2012; Nan and Gao, 2019). GhWRKY15-overexpressing plants exhibit enhanced resistance to TMV and cucumber mosaic virus (CMV) through the regulation of ROS signaling pathways (Yu et al., 2012). In Arabidopsis, AtWRKY8 is reported to mediate crosstalk between ABA and ethylene (ET) signaling pathways by directly regulating the expression of abscisic acid insensitive 4 (ABI4), 1-aminocyclopropane-1-carboxylic acid synthase 6 (ACS6), and ethylene response factor 104 (ERF104), thereby enhancing the defense response against TMV-cg (Chen et al., 2013).
Salicylic acid (SA) is a key plant defense hormone that has a remarkable impact on plant defense against various pathogens (Zhang and Li, 2019). Pathogen-induced SA is synthesized either from isochorismate synthase or via the phenylpropanoid biosynthesis pathway (Wildermuth et al., 2001; Ullah et al., 2019). SA was first reported as an inducer of plant disease resistance against TMV in tobacco (White, 1979). Since then, more and more studies show that SA plays an important role as a signaling molecule in plant defense responses against viral infection (Zhou et al., 2014; Yang L. et al., 2016; Li T. et al., 2019; Shaw et al., 2019). SA-mediated defense responses are usually associated with the induction of a number of pathogenesis-related (PR) proteins, which are usually considered to be markers for SA-mediated resistance to viral infection (Whitham et al., 2003; Love et al., 2005). An increasing body of evidence suggests that Non-expressor of PR1 (NPR1), NPR3, and NPR4 serve as SA receptors participating in two opposite signaling pathways downstream of SA (Wu et al., 2012; Ding et al., 2018).
Previous studies also show that WRKY proteins participate in SA-mediated plant immunity. For example, in Arabidopsis, 49 of the 74 WRKY genes were modulated when plants were treated with SA (Dong et al., 2003). In Arabidopsis, AtWRKY18 is a positive regulator of defense-related gene expression and is reported to increase resistance to the bacterial pathogen P. syringae (Chen and Chen, 2002). AtWRKY18, AtWRKY40, and AtWRKY60 are structurally similar and somewhat functionally redundant, and double or triple mutants of these proteins are reported to exhibit reduced resistance to P. syringae but greater resistance to Botrytis cinerea (Xu et al., 2006). Furthermore, the resistance of NbWRKY40-silenced plants to Phytophthora parasitica and B. cinerea can be induced in Nicotiana benthamiana through the SA-mediated signaling pathway (Ma et al., 2016). Interestingly, W-boxes are also found in the promoters of SA-biosynthesis genes and SA-signaling genes, such as NPR1, PR2, PR10, and Isochorismate synthase 1 (ICS1) (Yu et al., 2001; Li et al., 2004; van Verk et al., 2011). However, numerous relationships between WRKY proteins, SA, and viral infections remain largely unclear in plants.
The aim of this study is to analyze the regulation pattern of NbWRKY40 cloned from Nicotiana benthamiana in response to tomato mosaic virus (ToMV) infection and to identify related genes in order to elucidate potential antiviral mechanisms and to determine the involvement of SA-mediated pathways. Overall, the results of our study provide molecular evidence regarding the potential contributions of WRKY TFs to antiviral defense signaling in plants.
Materials and Methods
Plant Materials, Growth Conditions, and Treatments
Nicotiana benthamiana plants were grown in a greenhouse (26 ± 1°C, 60% RH) under long-day conditions (16 h light/8 h dark). Agrobacterium cultures that contained GFP-tagged recombinant ToMV (ToMV-GFP) derivatives were used to infect N. benthamiana plants (Liu et al., 2014). Cultures were shaken overnight at 28°C, pelleted, resuspended in an induction buffer (1 M MgCl2 and 10 mM MES, pH = 5.6, and 100 mM acetosyringone), diluted to an OD600 of 0.6, incubated at room temperature for 3 h, and then used to infiltrate the fourth leaves of 21-day-old plants. At 3 days postinoculation (dpi), plants were photographed under a 100-W handheld long-wave ultraviolet lamp using a digital camera (Canon EOS R, Tokyo, Japan). Some of the 21-day-old seedlings were also sprayed with an SA buffer (500 μM; Sigma-Aldrich, St. Louis, MO, United States). At appropriate times, total RNA was extracted from SA-treated leaves and stored at −80°C for future analysis. The process was repeated at least three times for each treatment.
NbWRKY40 Identification and Analysis
The open-reading-frame sequence of NbWRKY40 (Niben101Scf04944g05002.1) was downloaded from NCBI1. Homologs from other species were identified using an NCBI database sequence matching tool (blastn), and sequences with high similarity were downloaded. DNAMAN version 6.0 (Lynnon Biosoft, Quebec, Canada) was used to generate a multiple sequence alignment. MEGA 7.0 was used to perform a phylogenetic analysis. The neighbor-joining method was used to construct a phylogenetic tree with bootstrap values of 1000 (Kumar et al., 2016).
For NbWRKY40 gene cloning, total RNA was extracted from N. benthamiana seedlings using TRIzol reagent (Invitrogen, Carlsbad, CA, United States) and then treated with DNase (Invitrogen, Carlsbad, CA, United States). The full-length NbWRKY40 sequence was amplified using reverse-transcription PCR (RT-PCR), specific primers (Supplementary Table 1), and super-fidelity DNA polymerase (Phanta Max, Vazyme Biotech Co., Ltd., Nanjing, China). The purified PCR product was then inserted into a plasmid vector pGEM-T easy and transformed into E. coli (pGEM-T Easy vector; Promega, Madison, WI, United States), and the full NbWRKY40 was sequenced from positive clones.
Quantitative RT-PCR Analysis
RNA was extracted from leaves or protoplasts using TRIzol reagent (Invitrogen, Carlsbad, CA, United States) and then treated with DNase. cDNA was generated by reverse transcription using a First Strand cDNA Synthesis Kit (TOYOBO, Osaka, Japan), diluted 1:20, and used as a template for quantitative real-time PCR (RT-qPCR). RT-qPCR was performed using an AceQ RT-qPCR SYBR Green Master Mix (Vazyme, Nanjing, China) with the ubiquitin-conjugating enzyme (UBC) as an internal reference gene (Supplementary Table 1). At least three biological replicates, each with three technical replicates, were used for each treatment. Relative expression levels were calculated using the comparative 2–ΔΔCt method (Willems et al., 2008). The primers used for RT-qPCR are listed in Supplementary Table 1.
Generation of N. benthamiana Transgenic Lines
Specific primers (Supplementary Table 1) were used to obtain the full-length cDNA of NbWRKY40. pCV-GFP, which was constructed in our laboratory, was used in this study (Lu et al., 2011). Next, PCR products were double digested with BamHI and SacI and inserted into the BamHI and SacI sites of pCV-His in which the His tag was used to replace the GFP tag of pCV-GFP to create the fusion construct pCV-NbWRKY40-His. The recombinant clone was then transferred into Agrobacterium tumefaciens GV3101 for transformation of N. benthamiana. NbWRKY40 transgenic plants were confirmed by PCR amplification with NbWRKY40-specific primers until T2 generation.
Transcriptional Activation Activity Analysis
Full-length NbWRKY40 cDNA was generated using primers pGBKT7-NbWRKY40-F and pGBKT7-NbWRKY40-R (Supplementary Table 1), double digested using EcoRI and BamHI, and inserted into the EcoRI and BamHI sites of pGBKT7. The resulting fusion construct pGBKT7-NbWRKY40 as well as the pGBKT7-Solanum lycopersicum Abscisic acid insensitive 3-F (SlABI3-F) vector (positive control) (Gao et al., 2013, 2020) and the empty pGBKT7 vector (negative control) were then transformed into yeast receptor state Y2HGold using the PEG-LiAC method (Yang et al., 2010). The transformants were screened by streaking on SD medium lacking tryptophan (SD/-Trp) and on SD/-Trp medium with X-alpha-gal (SD/-Trp-X-α-Gal) with two replicates of each transformant on each plate. The pCAMBIA1300-35Smini-GUS recombinant plasmid, which was used as a reporter plasmid, was constructed by synthesizing three tandem W-box sequences (5′-CGTTGACCGTTGACCGAGTTGACCTTTTTA- 3′), ligating the fragment to the N-terminus of the CaMV 35S minimal promoter, and then substituting the resulting W-box-35S mini promoter for the CaMV 35S promoter in the pCAMBIA1300 vector. A mutant W-box (mW-box-35S mini, 5′-CGTAGACGGTAGACGGAGTAGACGTTTTTA-3′) was used as a negative control. Meanwhile, the pCB1300-NbWRKY40 recombinant plasmid, which was used as an effector plasmid, was constructed by amplifying the full NbWRKY40 cDNA, double digesting the product using EcoRI and BamHI, and then ligating the fragment to the EcoRI and BamHI sites of the pCAMBIA1300 vector. The effector and reporter plasmids were cotransformed into the A. tumefaciens strain GV3101 as described by (Yang et al. (2000), and the resulting strain was used to infiltrate N. benthamiana (Yang et al., 2000). At the same time, the effector SlWRKY8 and the reporter plasmids (W–box–35Smini–GUS) were cotransformed into the A. tumefaciens strain GV3101 as a positive control. GUS histochemical analysis was performed as described previously (Chen et al., 2013).
Subcellular Localization of NbWRKY40
To observe the subcellular localization pattern of NbWRKY40, the full–length NbWRKY40 protein was amplified using specific primers (Supplementary Table 1) and ligated to the N–terminus of GFP in the pCV–GFP to construct the plasmid pCV–NbWRKY40. Then, both pCV–NbWRKY40 and pCV–GFP were transformed individually into A. tumefaciens strain GV3101 by electroporation, and the fourth leaves of 21–day–old N. benthamiana plants were infiltrated with dual–transformed Agrobacterium. At 2 dpi, infiltrated leaves were collected and evaluated for GFP fluorescence under a Leica TCS SP5 confocal laser scanning microscope (Leica Microsystems, Heidelberg, Germany). Digital images were captured and postacquisition image processing was performed using Adobe Photoshop version 7.0 (Adobe Systems, Inc., San Jose, CA, United States).
Hormone Treatment
Twenty-one-day-old N. benthamiana seedlings were sprayed with either 500 μM SA (Sigma–Aldrich, PCode Number 101998016) diluted in sterile distilled water containing 0.1% Triton X–100 or 0.1% Triton X–100 alone (control). After 12 h, samples were inoculated using ToMV–GFP. At 4 dpi, N. benthamiana plants were photographed under long–wavelength UV light using a Canon digital camera, and leaves were harvested for further analysis.
Hormone Extraction and Analysis
Levels of ABA, SA, and JA extracted from plants were measured and analyzed by Zoonbio Biotechnology Co., Ltd. (Nanjing, China) as described previously (Forcat et al., 2008; Fu et al., 2012). N. benthamiana samples of approximately 1 g were ground in a precooled mortar containing 10 mL of an extraction buffer consisting of isopropanol/hydrochloric acid. Extracts were then shaken at 4°C for 30 min. Next, dichloromethane (20 mL) was added, and the samples were shaken again at 4°C for 30 min and centrifuged at 14,000 g at 4°C for 3 min. The organic phase was extracted and dried under liquid nitrogen. The pellets were dissolved in 150 mL methanol (0.1% methane acid) and filtered through a 0.22–mm filter membrane. The purified products were analyzed by performing high–performance liquid chromatography–tandem mass spectrometry (HPLC–MS/MS) using the following parameters: injection volume, 2 mL; spray voltage, 4500 V; air curtain pressure, 15 psi; atomizer pressure, 65 psi; auxiliary gas pressure, 70 psi; and atomization temperature, 400°C.
Virus–Induced Gene Silencing (VIGS)
Approximately 300 bp partial fragments of the NbWRKY40 coding sequence were RT–PCR–amplified using total RNA from N. benthamiana as a template using specific primers (Supplementary Table 1), double–digested using BamHI and SmaI, and then inserted into the TRV–RNA2 expression vector. Both the resulting TRV:NbWRKY40 and TRV–RNA1 vectors were electroporated into A. tumefaciens GV3101 to knock down NbWRKY40 expression as described previously (Yang et al., 2018). The empty TRV–RNA2 and TRV–RNA1 (TRV:00) were used to generate a negative control.
Protein Extraction and Analysis
Western blot assays were performed as described by Yang J. et al. (2016). Briefly, total protein was extracted using a lysis buffer (100 mM Tris–HCl, pH 8.8, 60% SDS, 2% β–mercaptoethanol) and protease inhibitor cocktail tablets (1 tablet per 50 mL buffer; Roche), and extracted proteins were separated using 12.5% SDS–PAGE and transferred to nitrocellulose membranes (Life Technologies), which were blocked for 2 h using TBST buffer (150 mM NaCl, 10 mM Tris–HCl, pH 8.0, 0.05% Tween–20) that contained 10% powdered milk. To detect GFP expression, the membranes were incubated for 4 h at room temperature with anti–GFP polyclonal antibodies (monoclonal antibody; 1:2,000 dilution; Quanshijin, Beijing, China), washed three times with TBST buffer, and then incubated with 1:5000 secondary antibody (Sigma) in TBST. Images were captured using Molecular Image ChemiDoc Touch (Bio–Rad, Hercules, United States), and the GFP signal was quantified using ImageJ.
Analysis of the cis–Regulatory Element of SA–Related Genes
In this study, 2000–bp sequences upstream of the translational start sites of the SA–related genes were considered to be promoter sequences. PlantCARE software2 was used to predict the cis-regulatory elements based on these promoter sequences.
Plasmodesmal Callose Staining
Nicotiana benthamiana epidermal cells injected with aniline blue dissolved in a sodium phosphate buffer at pH 7.5 were incubated in the dark for 5 min. Then, the injected leaf tissue was dissected, washed with sterile water, and observed under a Leica TCS SP5 confocal laser scanning microscope (Leica Microsystems, Heidelberg, Germany).
N. benthamiana Protoplast Isolation and Transfection
Nicotiana benthamiana protoplasts were isolated from 21-day-old seedlings. In short, the young leaves were chopped up and immersed in enzyme solution (0.5 M mannitol, 1.5% cellulose RS (Yakult Honsha, Tokyo, Japan), 0.75% macerozyme R10 (Yakult Honsha), 1 mM CaCl2, and 0.1% BSA). The mixture was incubated on a shaking incubator (60 rpm) in the dark at room temperature for 4–5 h and then filtered through Miracloth. The protoplasts were pelleted by centrifugation at 200 g for 5 min and then resuspended in an equal volume of W5 solution (154 mM NaCl, 125 mM CaCl2, 5 mM KCl and 1.5 mM MES, adjusted to pH 5.7), followed by centrifugation and resuspension in MMG solution (0.4 M mannitol, 15 mM MgCl2 and 4.7 mM MES, adjusted to pH 5.7). Plasmid DNA (10 or 20 μg) was added to the protoplast solution and transfected with 40% polyethylene glycol (PEG) solution [40% PEG4000 (Sigma-Aldrich, St. Louis, MO, United States, PCode Number, 102078930), 0.4 M mannitol, and 100 mM Ca(NO3)2] at room temperature for 20 min. W5 solution was gradually added to dilute the PEG solution and then discarded. The transfected protoplasts were incubated overnight at room temperature and then observed under a confocal microscope.
Electrophoretic Mobility Shift Assay (EMSA)
For the EMSA, specific promoter fragments of SA-related genes containing the W-box were synthesized as biotin end labels. The unlabeled W-box oligonucleotide served as a competitor. The stabilized streptavidin-horseradish peroxidase conjugate was used for super shift identification. The assay was performed using a LightShift® Chemiluminescent EMSA Kit (Thermo Scientific, Waltham, MA, United States) according to the manufacturer’s instructions.
Results
NbWRKY40 Cloning and Expression Characterization During ToMV Infection
A previous study shows that NbWKRY40 is involved in the N. benthamiana response to inoculation with Phytophthora parasitica or Botrytis cinerea (Ma et al., 2016). In this study, RNA-seq analysis revealed significantly lower levels of NbWRKY40 expression in N. benthamiana after infection with ToMV than in non-infected control plants. Therefore, in order to investigate the effect of ToMV infection on the relative expression of NbWRKY40 in N. benthamiana, we performed RT-qPCR. The results show that the lowest expression levels were detected in inoculated leaves at 12 dpi (0.23-fold) and in systemic leaves at 10 dpi (0.20-fold) compared to controls (Figures 1A,B). These results indicate that NbWRKY40 expression is suppressed during ToMV infection, potentially to promote ToMV infection.
Figure 1. NbWRKY40 cloning and expression characterization during ToMV infection. (A) Time course of NbWRKY40 levels in inoculated leaves following infection with ToMV. (B) Time course of NbWRKY40 levels in systemic leaves following infection with ToMV. The relative expression of NbWRKY40 in ToMV-inoculated and systemic leaves was measured using quantitative real-time polymerase chain reaction and the 2–ΔΔCt method with Nbubiquitin as an internal control. Values and error bars represent the mean ± the SD (n = 3 with three technical replicates for each biological replicate). **P < 0.01, n.s., not significant, based on Dunnett’s test. (C) Sequence alignment of Oryza sativa (Os), Arabidopsis thaliana (At), Nicotiana sylvestris (Ns), and Nicotiana benthamiana (Nb) WRKY proteins. Identical amino acids are highlighted in black. (D) Phylogenetic tree of Arabidopsis thaliana (At), Nicotiana sylvestris (Ns), Nicotiana benthamiana (Nb), Oryza sativa (Os), and Zanthoxylum bungeanum (Zb) WRKY proteins. The red triangle indicates NbWRKY40. The phylogenetic tree was constructed using the neighbor-joining method in MEGA 7.0. Numbers at nodes indicate bootstrap values based on 1000 resamplings.
To further investigate the function of NbWRKY40 in plant resistance, we analyzed the cDNA of NbWRKY40 (Niben101Scf04944g05002.1) isolated from N. benthamiana. RT-PCR confirmed that the cDNA of NbWRKY40 was 768 bp in length and represented the complete open reading frame. The predicted NbWRKY40 protein was 254 amino acids in length with a molecular weight of 28.95 kDa. Amino acid sequence alignment revealed similarity (37.82–95.29%) to homologs from other species, including Nicotiana sylvestris WRKY40 (NsWRKY40, 95.29%), Arabidopsis thaliana WRKY60 (AtWRKY60, 46.85%), AtWRKY40 (44.89%), AtWRKY18 (44.74%), and Oryza sativa WRKY76 (OsWRKY76, 37.82%) (Figure 1A). Similar to other members of the WRKY group II family, NbWRKY40 also contained a C2H2 motif (C-X5-C-X23-H-X1-H) and a conserved WRKYGQK core sequence (Figure 1C). To demonstrate the evolutionary relationship between NbWRKY40 and WRKY proteins of O. sativa, A. thaliana, and N. sylvestris, a phylogenetic tree was constructed using the neighbor-joining method in MEGA 7.0. The resulting phylogenetic tree indicated that NbWRKY40 was most similar to IIa subgroup proteins, such as AtWRKY60, AtWRKY18, and AtWRKY40 (Figure 1C). Furthermore, we also performed amino acid sequence alignment with homologs NbWRKY40a, NbWRKY40b, NbWRKY40c, NbWRKY40d, and NbWRKY40e, which were analyzed in a previous study (Ma et al., 2016). The sequence alignment showed that the NbWRKY40 cloned in this study shared 100% sequence identity with NbWRKY40e and exhibited 40.87, 29.91, 39.22, and 44.07% sequence identities with NbWRKY40a, NbWRKY40b, NbWRKY40c, and NbWRKY40d, respectively (Supplementary Figure 1).
NbWRKY40 Activity in Yeast and Plants
In order to determine the transcriptional activity of the NbWRKY40 protein, a yeast one-hybrid assay was performed. Yeast cells transformed with pGBKT7-NbWRKY40, the empty pGBKT7 vector (negative control), or the abscisic acid insensitive pGBKT7-SlABI3-F vector (positive control) were all able to form white colonies on the SD/-Trp medium (Figure 2A). However, when transformants were streaked on the SD/-Trp/X-α-Gal medium, only the yeast transformant carrying pGBKT7-SlABI3-F was able to form colonies that turned blue, whereas transformants containing an empty pGBKT7 vector or pGBKT7-NbWRKY40 formed white colonies (Figure 2A). These results suggest that the NbWRKY40 protein did not possess transcriptional activation activity in yeast.
Figure 2. Transcriptional activation ability of NbWRKY40. (A) Transcriptional activation ability of NbWRKY40 in yeast cells. (B) Triple tandem repeats of the W-box and mW-box. (C) Structures of reporter and effector constructs. (D) In vivo histochemical analysis of GUS activity in cotransfected Nicotiana benthamiana leaves. GUS staining was analyzed in the fourth leaves of 28-day-old plants that were Agrobacterium-infiltrated with the reporter and effector at an OD600 of 0.6. (E) Electrophoretic mobility shift assay (EMSA) analysis showing that NbWRKY40 binds to the W-box promoter in vitro. The “shift band” black arrow indicates the binding of NbWRKY40 to the biotin-labeled W-box promoter. The “+” indicates the presence of the corresponding component, whereas “–” indicates the absence of the corresponding component.
As previously reported, WRKY TFs can bind the cis-element W-box to regulate gene transcription activity (Rushton et al., 2010). Based on this hypothesis, a transient co-expression assay was performed to determine whether the transcription of genes containing the cis-element W-box in their promoter sequences could be activated by NbWRKY40 in plant cells (Figures 2B,C). Tobacco leaves that were cotransformed with W-box-35Smini-GUS and 35S:NbWRKY40 exhibited strong GUS staining like that observed for leaves that were cotransformed with W-box-35Smini-GUS and 35S:SlWRKY8, which was used as a positive control (Gao et al., 2020), whereas those leaves transformed with only the reporter vector (W-box-35Smini-GUS) exhibited much less staining (Figure 2D). No staining was observed in leaves inoculated with mW-box-35Smini-GUS and/or 35S:NbWRKY40 (Figure 2D). These results indicate that GUS transcription can be activated by the binding of NbWRKY40 to W-box motifs.
To further confirm the binding of NbWRKY40 to the W-box element, an EMSA was performed. The analysis indicates that NbWRKY40 was bound to the W-box, but this did not occur in the presence of excess unlabeled W-box (Figure 2E). GST alone was used as a negative control to confirm protein-DNA specificity. Taken together, these results indicate that NbWRKY40 functions as a TF that activates the expression of genes that contain the W-box element in their respective promoter sequences.
Effect of ToMV on NbWRKY40 Localization
To examine the subcellular localization of NbWRKY40, a fusion plasmid (NbWRKY40-GFP) containing a C-terminal GFP tag was constructed and introduced into N. benthamiana epidermal cells using Agrobacterium infiltration. At 2 dpi, NbWRKY40-GFP was mainly localized in the nuclei as demonstrated by co-localization with DAPI (Figure 3A). Next, to determine whether the localization of NbWRKY40 could be affected by ToMV infection, NbWRKY40-GFP was co-infiltrated with ToMV into N. benthamiana epidermal cells by Agrobacterium infiltration. At 2 dpi, strong GFP fluorescence was observed in the nuclei of plants infiltrated with either NbWRKY40-GFP alone or co-infiltration with ToMV (Figure 3A). However, western blot assays revealed that NbWRKY40 expression levels were lower in plants under ToMV stress than those in the mock-inoculated controls (Figures 3B,C).
Figure 3. Subcellular localization of NbWRKY40 in Nicotiana benthamiana following infection with ToMV. (A) Subcellular localization of NbWRKY40 in response to ToMV infection. Agrobacterium infiltration was used to introduce recombinant plasmid NbWRKY40-GFP + ToMV into healthy N. benthamiana epidermal cells. Healthy leaves infiltrated with an Agrobacterium culture containing NbWRKY40-GFP or infiltrated with an Agrobacterium culture that only contained GFP acted as controls. The confocal microscopy images were captured 2 dpi under bright-field fluorescence to show cell morphology; under dark field to show green fluorescence, indicating localization of the NbWRKY40 protein, and blue fluorescence, indicating nuclei stained blue by 4,6-diamidino-2-phenyl-indole dihydrochloride (DAPI); and under combination fluorescence to show the three images merged. Scale bar, 100 μm. (B) Immunoblot of total protein extracted from the leaves of ToMV-inoculated (+ToMV) and mock-inoculated (–ToMV) plants. Anti-GFP was used to detect GFP fractions. Equal amounts of protein were used for immunoblotting and for staining with Coomassie blue (CBB). (C) Immunoblot of nuclear protein extracted from the leaves of ToMV-inoculated and mock-inoculated plants. Anti-GFP was used to detect GFP fractions. Equal amounts of protein were used for immunoblotting and for staining with CBB.
Effect of NbWRKY40 Expression on ToMV Accumulation
To further investigate the relationship between NbWRKY40 expression and ToMV infection, transgenic N. benthamiana plants overexpressing NbWRKY40 lines 23, 24, and 50 (L-23, L-24, and L-50) were constructed and evaluated using RT-qPCR and western blot analysis (Supplementary Figures 2A,B). As expected, higher levels of NbWRKY40 accumulated in transgenic lines than in WT plants. Next, to determine the role of NbWRKY40 in the plant defense response to ToMV infection, the responses of transgenic and WT plants to ToMV stresses were investigated. L-23, L-24, L-50, and WT plants were inoculated with ToMV-GFP. At 7 dpi, prominent areas of GFP fluorescence were observed in systemic leaves of infiltrated plants. Interestingly, less GFP fluorescence was observed in systemic leaves of L-23, L-24, and L-50 plants infiltrated with ToMV-GFP than in WT plants infiltrated with ToMV-GFP (Figure 4A). Western blot analysis also indicates that lower levels of GFP had accumulated in systemic leaves of L-23, L-24, and L-50 plants inoculated with ToMV-GFP than in the WT (Figure 4B). Transient silencing of NbWRKY40 was achieved using TRV vector-based VIGS. At 7 dpi, plants infiltrated with TRV1 + TRV2-NbWRKY40 were confirmed by RT-qPCR (Supplementary Figure 2C). As such, three VIGS (#2, #5, and #8) and vector control (TRV:00) plants were selected for further study (i.e., challenged with ToMV-GFP). At 7 dpi, more GFP fluorescence was observed in systemic leaves of #2, #5, and #8 plants than in those of TRV:00 plants (Figure 4D). Western blot analysis also confirmed that higher levels of GFP had accumulated in the #2, #5, and #8 plants than in plants co-infiltrated with ToMV-GFP and TRV:00 (Figure 4E). Moreover, RT-qPCR analysis also confirmed that the expression level of the ToMV CP was lower in L-23, L-24, and L-50 plants than in WT plants after ToMV infection (Figure 4C), but much higher in #2, #5, and #8 plants (Figure 4F). These findings clearly demonstrate that NbWRKY40 expression levels were negatively related to ToMV accumulation levels in N. benthamiana.
Figure 4. Expression level of NbWRKY40 affects ToMV accumulation in Nicotiana benthamiana. (A) Systemic ToMV-GFP fluorescence in wild-type (WT) and T2-homozygous transgenic (L-23, L-24, and L-50) plants at 7 days postinfiltration (dpi). (B) Western blot analysis of GFP accumulation in ToMV-GFP-inoculated WT, L-23, L-24, and L-50 plants at 7 dpi. Coomassie blue (CBB)-stained rubisco gel and ImageJ (United States National Institutes of Health, http://rsb.info.nih.gov/nih-image/) were used to determine protein loading. (C) Relative expression level of ToMV-CP in NbWRKY40 overexpressing plants. (D) Systemic ToMV-GFP fluorescence of TRV:00- and TRV:WRKY40 (#2, #5, #8)-treated plants at 7 dpi. (E) Western blot analysis of GFP expression at 7 dpi. CBB-stained rubisco gel and ImageJ were used to determine protein loading. (F) Relative expression level of ToMV-CP in NbWRKY40-silenced plants. The relative expression level of ToMV-CP was calculated using the 2–ΔΔCt method. The expression level of the Nbubiquitin gene in N. benthamiana was used as an internal control. Values and error bars represent the mean ± SD of three independent biological replicates with three technical replicates per sample. **P < 0.01 based on Student’s t-test.
Effect of NbWRKY40 on SA-Related Genes and the Deposition of PD Callose in Plants
WRKY genes are reported to participate in plant hormone-mediated signaling pathways (Ishihama and Yoshioka, 2012; Nan and Gao, 2019). To assess the effect of NbWRKY40 on plant hormone levels, the ABA, SA, and JA content levels of L-23 and WT plant tissues were compared. L-23 plants contained significantly higher levels of SA than WT plants but similar levels of ABA and JA (Figure 5A). To confirm this finding, RT-qPCR was used to measure the expression levels of SA-related genes. The expression levels of ICS1, PR1b, and PR2 were significantly higher in L-23 plants than in WT plants (Figure 5B) but were significantly lower in NbWRKY40-silenced plants than in WT plants (Figure 5C).
Figure 5. Effect of NbWRKY40 on salicylic acid signaling. (A) Hormone levels in the leaves of wild-type (WT) and T2-homozygous transgenic (L-23) Nicotiana benthamiana. ABA, abscisic acid; SA, salicylic acid; JA, jasmonic acid. (B) Expression of the SA synthesis genes, ICS1 (isochorismate synthase 1), SAMT (SA methyl transferase), and SABP2 (salicylic acid-binding protein 2), and the SA-dependent signaling-related genes, PR2 and PR1b (pathogenesis-related protein) in WT and L-23 plants. (C) Expression of SA-related genes in VIGS and control N. benthamiana plants. (D) Aniline blue staining of WT-, L-23, TRV:00- and TRV:NbWRKY40-treated leaves revealing the callose in plasmodesmata (PDs) and guard cells. (E) Transcript levels of the callose biosynthesis gene callose synthase 1 (Cals1) in L-23 plants and WT plants (E) and in TRV-NbWRKY40 plants and TRV:00 plants (F). (G) Relative transcript levels of viral RNA in protoplasts of WT and L-23 plants inoculated with ToMV. (H) Relative transcript levels of viral RNA in protoplasts of TRV:00 and TRV:NbWRKY40 plants inoculated with ToMV. (I) Binding of WRKY40 to ICS1 promoter fragments. The expression of Nbubiquitin was used as an internal control. Values and error bars indicate means ± SD (n = 3 with three technical replicates for each biological replicate; **P < 0.01, n.s., not significant, based on Student’s t-test).
Callose deposition at plasmodesmata (PD) is regulated by SA (Wang et al., 2013; Cui and Lee, 2016). We used aniline blue staining to assess the amount of callose deposited at PDs. A greater amount of callose was deposited on PDs in L-23 plant leaves than on PDs in WT plant leaves; however, only a very small amount of callose was deposited on PDs in leaves of NbWRKY40-silenced plants compared with that in TRV:00-treated, control plants (Figure 5D and Supplementary Figure 3). To understand this phenomenon, we used RT-qPCR to analyze the relative expression of the callose biosynthesis gene callose synthase 1 (Cals1). The transcript level of Cals1 in L-23 plants was significantly higher than that in WT plants, whereas the expression of CalS1 was significantly lower in TRV-NbWRKY40 plants than in TRV:00 plants (Figures 5E,F).
To confirm the possible role of NbWRKY40 in ToMV replication, we conducted protoplast transfection assays. Protoplasts were isolated from N. benthamiana and then transfected with ToMV-GFP. At 48 h post transfection, GFP fluorescence was observed using confocal microscopy (Supplementary Figure 4). RT-qPCR was performed to monitor viral (−) RNA accumulation in protoplasts infiltrated with ToMV-GFP. The relative expression levels of viral RNA and GFP fluorescence in WT, L-23, TRV:00, and TRV-NbWRKY40 plants indicates that NbWRKY40 did not affect the replication of ToMV in protoplasts (Figures 5G,H).
WRKY protein is generally thought to bind to the consensus W-box sequence TTGAC (C/T). In order to identify whether SA-related genes contain the cis-regulatory element W-box, the PlantCARE tool was used to analyze the 2000-bp sequence upstream of the putative translation start site of the SA-related genes. As shown in Supplementary Figure 5, the SA-related gene NbICS1 contains the WRKY identification area of the W-box, which is located at 1673 to 1678 bp upstream of the ATG start codon. As a first step toward the characterization of WRKY40 binding sites in the NbICS1 gene promoter, we prepared 84-bp promoter fragments that contained a TTGAC (C/T) core sequence in the center. After biotin labeling, the fragments were assayed for their ability to bind NbWRKY40 through EMSA. The results show that NbWRKY40 protein can recognize and bind to the biotin-labeled probes of NbICS1 to cause mobility changes (Figure 5I). These findings indicate that the NbWRKY40 protein can bind to the W-box in the promoter region of NbICS1.
Effect of SA on ToMV Infection
Salicylic acid treatment affects the expression of a large number of WRKY genes in Arabidopsis, suggesting an important role for WRKY-mediated transcriptional control in gene expression (Dong et al., 2003). We used RT-qPCR to investigate the effect of SA treatment on NbWRKY40 expression. At 2 h after SA treatment, the expression level of NbWRKY40 was 12.3-fold higher than that in the control plants (Figure 6A). In addition, to investigate whether an exogenous application of SA can reduce viral accumulation in NbWRKY40-silenced plants infected with ToMV, NbWRKY40-silenced seedlings were treated with 500 μM SA or 0.1% triton X-100 (which acted as a negative control) before inoculation with ToMV-GFP. At 3 dpi, stronger GFP fluorescence was observed in the inoculated leaves of plants sprayed with 0.1% Triton X-100 than in those sprayed with SA (Figure 6B), and western blot analysis indicated that lower levels of GFP had accumulated in SA-treated leaves than in controls (Figure 6C). At 7 dpi, SA-treated plants exhibited weaker GFP fluorescence and lower levels of GFP accumulation than control plants (Figures 6D,E), which indicates that SA pretreatment can reduce the susceptibility of NbWRKY40-silenced plants to ToMV. Moreover, to determine whether the localization of NbWRKY40 could be affected by SA treatment, N. benthamiana was pretreated with SA and then infiltrated with NbWRKY40-GFP. At 2 dpi, the NbWRKY40-GFP protein mainly accumulated in the nucleus with stronger GFP fluorescence observed in SA-treated leaves than in control leaves (Figure 6F). Western blot assays also indicated that the protein expression level of NbWRKY40-GFP was higher in plants pretreated with SA than in those that did not receive the SA treatment (Figures 6G,H).
Figure 6. Effect of exogenous salicylic acid (SA) application on viral infection and subcellular localization of NbWRKY40. (A) Relative expression of NbWRKY40 in response to SA treatment (B) GFP fluorescence of NbWRKY40-silenced leaves inoculated with ToMV-GFP at 3 days postinfiltration (dpi) following SA treatment is less intense than that following 0.1% Triton X-100 treatment alone. (C) Western blot analysis of GFP accumulation in ToMV-GFP-inoculated leaves at 3 dpi. Coomassie brilliant blue (CBB)-stained rubisco gel and ImageJ (United States National Institutes of Health, http://rsb.info.nih.gov/nih-image/) were used to determine protein loadings. (D) GFP fluorescence of systemic leaves of NbWRKY40-silenced plants inoculated with ToMV-GFP at 7 dpi following SA treatment is less intense than that following 0.1% Triton X-100 treatment alone. (E) Western blot analysis of GFP accumulation in non-infected leaves of ToMV-GFP-inoculated plants at 7 dpi. CBB-stained rubisco gel and ImageJ (United States National Institutes of Health, http://rsb.info.nih.gov/nih-image/) were used to determine protein loadings. (F) Subcellular localization of NbWRKY40 protein in response to SA treatment. The recombinant plasmid NbWRKY40-GFP was introduced into Nicotiana benthamiana epidermal cells using Agrobacterium infiltration. Leaves infiltrated with an Agrobacterium culture that only contained GFP acted as controls. The confocal microscopy images were captured under bright-field fluorescence to show cell morphology; under dark field to show green fluorescence, indicating localization of the NbWRKY40 protein, and blue fluorescence, indicating nuclei stained blue by 4,6-diamidino-2-phenyl-indole dihydrochloride (DAPI); and under combination fluorescence to show the three images merged. Scale bar, 100 μm. (G) Immunoblot of total protein extracted from the leaves of SA- and mock-treated plants. Anti-GFP was used to detect GFP fractions. Equal amounts of protein were used for immunoblotting and for staining with CBB. (H) Immunoblot of nuclear protein extracted from the leaves of SA- and mock-treated plants. Anti-GFP was used to detect GFP fractions. Equal amounts of protein were used for immunoblotting and for staining with CBB. **P < 0.01.
Discussion
WRKY TFs comprise one of the largest families of plant transcription factors, and many WRKY genes have been reported to play important roles in host defense mechanisms (Huh et al., 2012; Ishihama and Yoshioka, 2012). However, the functional roles of WRKYs and their involvement in defense mechanisms against pathogen infection remain unclear. In this study, we cloned a WRKY TF gene, NbWRKY40, from N. benthamiana. Multiple sequence alignments indicated that NbWRKY40 was similar to five other WRKY proteins (OsWRKY76, AtWRKY60, AtWRKY40, AtWRKY18, and NsWRKY40; Figure 1C), and phylogenetic analysis indicated that NbWRKY40 was most similar to subgroup IIa WRKY proteins (e.g., AtWRKY40, AtWRKY18, AtWRKY60, and OsWRKY76; Figure 1D). Among the identified homologs, OsWRKY76 is a negative regulator of defense-related metabolite biosynthesis and affects SA production through its participation in the phenylpropanoid pathway (Liang et al., 2017), which suggests that NbWRKY40 may also be involved in plant defense mechanisms.
Consistent with the domain of WRKY proteins, NbWRKY40 contained both the conserved WRKY motif and a zinc-finger motif in the central region of the genome (Figure 1C). Our confocal microscopy observations indicate that NbWRKY40 targets the nucleus exclusively, at least in the leaf epidermal cells of N. benthamiana (Figure 3A). To determine whether NbWRKY40 is a transcriptional activator or repressor, we also performed Agrobacterium-mediated transient co-expression analyses, which reveal that NbWRKY40 can activate GUS transcription (Figures 2B–D). Furthermore, the binding of NbWRKY40 to the W-box motif is confirmed by EMSA (Figure 2E). Taken together, these results suggest that NbWRKY40 functions as a transcription activator by specifically binding the W-box motif, thereby regulating gene expression.
Previous studies show that the transcripts of WRKY genes can be strongly induced by several different pathogens, including Xanthomonas oryzae, Ralstonia solanacearum, Coniothyrium diplodiella, and CMV (Shi et al., 2014; Hwang et al., 2016; Zhang et al., 2016; Zou et al., 2019). However, NbWRKY40 expression was significantly reduced in ToMV-infected plants (Figures 1A,B). Generally, viral infection also affects the subcellular distribution of the host protein, leading to the loss of normal functions. For example, rice stripe virus (RSV) infection affects the distribution pattern of OsHSP20 in rice cells (Li et al., 2015). Here, we demonstrate that ToMV infection did not affect the localization of NbWYKY40 in the nucleus but did inhibit expression levels of NbWYKY40 in the plant cells (Figure 3A). Furthermore, the constitutive expression of NbWRKY40 in transgenic plants enhanced the resistance of plants to ToMV infection, whereas NbWRKY40-silenced plants were more susceptible to infection (Figure 4). Similarly, many WRKY TFs have been confirmed as positive regulators in responses to viral infection, such as SlWRKY8, AtWRKY30, and CaWRKYd (Dang et al., 2014; Zou et al., 2019; Gao et al., 2020). Taken together, our findings suggest that NbWRKY40 plays a positive role in regulating the defense response of N. benthamiana against ToMV infection.
WRKY TFs are known to be associated with the SA-, JA-, and ABA-mediated signaling pathways (Dong et al., 2003). In our study, the SA level was dramatically higher in NbWRKY40 transgenic plants than in WT (Figure 5A). Pathogen-induced SA is predominantly biosynthesized from chorismate by ICS1, followed by the induction of PR gene expression (Wildermuth et al., 2001; Whitham et al., 2003; Love et al., 2005; Hao et al., 2018). Our RT-qPCR analyses show that expression levels of ICS1, PR2, and PR1b were significantly more upregulated in the L-23 plants than in WT plants, whereas these genes were downregulated when NbWRKY40 was silenced (Figures 5B,C). Previous studies report that SA can enhance plant antiviral defense through inhibiting the movement and replication of several plant viruses (Murphy and Carr, 2002; Wong et al., 2002). In this study, we observed that both locally and systemically infected leaves of NbWRKY40-silenced plants pretreated with SA were less susceptible to ToMV infection (Figure 6), which suggests that SA is negatively related to ToMV accumulation in N. benthamiana. Furthermore, RT-qPCR analyses show that NbWRKY40 did not affect the replication of ToMV in protoplasts, indicating that NbWRKY40 affects viral infection, possibly through inhibiting the movement of ToMV. The movement of plant viruses between cells is controlled by the deposition of callose at the neck of PD during virus infection (Li et al., 2012; Cui et al., 2018). In Arabidopsis, SA has been shown to participate in the regulation of callose deposition at PDs by affecting the expression of Cals1 and Cals8 (Wang et al., 2013; Cui and Lee, 2016). Here, we find that overexpression of NbWRKY40 did increase SA levels (Figure 5A). Furthermore, PD callose deposition was significantly higher in L-23 plants than in WT plants; however, levels were significantly lower in NbWRKY40-silenced plants than in WT plants (Figure 5F). These results indicate that PD callose deposition levels in N. benthamiana plants were related to the accumulation of SA controlled by NbWRKY40. WRKY TFs can regulate gene expression by binding the W-box element at the promoter of genes. WRKY28 and WRKY46 bind to the W-box at the promoters of ICS1, leading to increased ICS1 expression in Arabidopsis protoplasts (van Verk et al., 2011). Consistent with these previous findings, our results also indicate that NbWRKY40 functions as an activator through binding the W-box element of ICS1 promoters (Figures 2, 5) to promote SA synthesis. In addition, the relative expression of NbWRKY40 transcripts was induced by SA treatment (Figure 6A). The inoculation of Capsicum annum with R. solanacearum has been reported to induce CaCDPK15 and indirectly activate downstream CaWRKY40, which, in turn, potentiates CaCDPK15 expression, creating a positive feedback loop. It is thought that this positive feedback loop amplifies defense signaling against R. solanacearum infections and efficiently activates strong plant immunity (Shen et al., 2016). Thus, we speculate that NbWRKY40 plays a role in a positive feedback loop for SA synthesis to amplify defense signaling. Overall our findings suggest that the NbWRKY40 restricts ToMV infection, possibly through regulating the expression of SA, resulting in the deposition of callose at the neck of PD to inhibit viral movement.
Conclusion
In this study, we identify a WRKY transcription factor in N. benthamiana, namely NbWRKY40, that is induced by ToMV infection. NbWRKY40 overexpression improved resistance to ToMV stress, whereas knockdown of the gene enhanced the susceptibility of NbWRKY-silenced plants to ToMV infection. Analysis of molecular mechanisms involved in enhanced resistance to ToMV revealed the extensive roles of NbWRKY40 in upregulating the expression of SA-synthesis and/or -responsive genes. Callose staining reveals that PD of overexpression WRKY40 plants were less permeable than those of WT plants, whereas those of NbWRKY40-silenced plants were more permeable than WT. Furthermore, our analyses show that NbWRKY40 did not affect the replication of ToMV in protoplasts. Taken together, our findings suggest that NbWRKY40 may mediate resistance to ToMV, possibly through regulating the expression of SA, resulting in the deposition of callose at the neck of PD to inhibit viral movement.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.
Author Contributions
JY and JC conceived the project and designed the experiments. YJ, WZ, and JY conducted the experiments with assistance from JL, PL, KZ, PJ, and MX. All authors analyzed and discussed the results. YJ, JY, and JC wrote the manuscript.
Funding
This work was funded by the National Key R&D Plan in China (2018YFD0200408, 2018YFD0200507, and 2017YFD-0201701), the National Natural Science Foundation of China (31901954), the Natural Science Foundation of Ningbo City (2019A610415 and 2019A610410), the National Key Project for Research on Transgenic Biology (2016ZX08002-001), the China Agriculture Research System from the Ministry of Agriculture of the P.R. China (CARS-03), and K.C. Wong Magna Funding in Ningbo University.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
We are very grateful to Prof. Jie Zhou (Zhejiang Academy of Agricultural Sciences, Hangzhou, China) for kindly providing us with the pCAMBIA1300-35Smini-GUS vector. We thank Professor Qian-Sheng Liao (College of Life Science, Zhejiang SCI-Tech University, Hangzhou, China) for providing the ToMV-GFP vector.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2020.603518/full#supplementary-material
Supplementary Figure 1 | Alignment of the amino acid sequences of the NbWRKY40 protein with the NbWRKY40a (NbS00048676g0011.1), NbWRKY40b (NbS00010134g0012.1), NbWRKY40c (NbS00051158g0004.1), NbWRKY40d (NbS00001817g0013.1), and NbWRKY40e (NbS00036148g0004.1).
Supplementary Figure 2 | Assessment of transgene and silencing efficiency in Nicotiana benthamiana. (A) Western blot analysis of transgenic lines (L-23, L-24, and L-50). (B) Relative expression levels of NbWRKY40 in transgenic lines (L-23, L-24, and L-50) and wild-type (WT) plants. (C) Relative expression levels of NbWRKY40 in VIGS and control N. benthamiana plants determined by RT-qPCR.
Supplementary Figure 3 | Relative intensity of aniline blue staining of WT-, L- 23-, TRV: 00-, and TRV:NbWRKY40-treated leaves to assess the amount of callose deposited at PDs and guard cells. Error bars represent the SD of the means of three biological repeats. A two-sample unequal variance directional t-test was used to test the significance of the difference (∗∗P-value < 0.01).
Supplementary Figure 4 | Confocal microscopy images of GFP fluorescence in Nicotiana benthamiana protoplasts 48 h post transfection with ToMV-GFP. Protoplasts were isolated from the transgenic NbWRKY40 overexpression line 23 (L-23), wild type (WT), transiently NbWRKY40 silenced plants (TRV:NbWRKY40), and N. benthamiana inoculated with TRV:00. Scale bar, 100 μm.
Supplementary Figure 5 | Schematic of the W-box element of SA-related gene promoters. Promoter sequences (2000 bp) upstream of genes were chosen for cis-regulatory element analysis using the PlantCARE online tool (http://www.dna.affrc.go.jp/PLACE/).
Supplementary Table 1 | Primers used in this study.
Footnotes
References
Chen, C., and Chen, Z. (2002). Potentiation of developmentally regulated plant defense response by AtWRKY18, a pathogen induced Arabidopsis transcription factor. Plant Physiol. 129, 706–716. doi: 10.11104/pp.001057
Chen, L., Zhang, L., Li, D., Wang, F., and Yu, D. (2013). WRKY8 transcription factor functions in the TMV-cg defense response by mediating both abscisic acid and ethylene signaling in Arabidopsis. Proc. Natl. Acad. Sci. U.S.A. 110, E1963–E1971. doi: 10.1073/pnas.1221347110
Cui, W., and Lee, J. Y. (2016). Arabidopsis callose synthases CalS1/8 regulate plasmodesmal permeability during stress. Nat. Plants 2:16034. doi: 10.1038/nplants.2016.34
Cui, X., Lu, L., Wang, Y., Yuan, X., and Chen, X. (2018). The interaction of soybean reticulon homology domain protein (GmRHP) with Soybean mosaic virus encoded P3 contributes to the viral infection. Biochem. Biophys. Res. Commun. 495, 2105–2110. doi: 10.1016/j.bbrc.2017.12.043
Dang, F., Wang, Y., She, J., Lei, Y., Liu, Z., Eulgem, T., et al. (2014). Overexpression of CaWRKY27, a subgroup IIe WRKY transcription factor of Capsicum annuum, positively regulates tobacco resistance to Ralstonia solanacearum infection. Physiol. Plant 150, 397–411. doi: 10.1111/ppl.12093
Ding, Y., Sun, T., Ao, K., Peng, Y., Zhang, Y., Li, X., et al. (2018). Opposite roles of salicylic acid receptors NPR1 and NPR3/NPR4 in transcriptional regulation of plant immunity. Cell 173, 1454–1467–15. doi: 10.1016/j.cell.2018.03.044
Dong, J., Chen, C., and Chen, Z. (2003). Expression profiles of the Arabidopsis WRKY gene superfamily during plant defense response. Plant Mol. Biol. 51, 21–37. doi: 10.1023/a:1020780022549
Eulgem, T., Rushton, P. J., Schmelzer, E., Hahlbrock, K., and Somssich, I. E. (1999). Early nuclear events in plant defence signalling: rapid gene activation by WRKY transcription factors. EMBO J. 18, 4689–4699. doi: 10.1093/emboj/18.17.4689
Forcat, S., Bennett, M. H., Mansfield, J. W., and Grant, M. R. (2008). A rapid and robust method for simultaneously measuring changes in the phytohormones ABA, JA and SA in plants following biotic and abiotic stress. Plant Methods 4:16. doi: 10.1186/1746-4811-4-16
Fu, J., Chu, J., Sun, X., Wang, J., and Yan, C. (2012). Simple, rapid, and simultaneous assay of multiple carboxyl containing phytohormones in wounded tomatoes by UPLC-MS/MS using single SPE purification and isotope dilution. Anal. Sci. 28, 1081–1087. doi: 10.2116/analsci.28.1081
Gao, Y., Liu, J., Zhang, Z., Sun, X., Zhang, N., Fan, J., et al. (2013). Functional characterization of two alternatively spliced transcripts of tomato ABSCISIC ACID INSENSITIVE3 (ABI3) gene. Plant Mol. Biol. 82, 131–145. doi: 10.1007/s11103-013-0044-1
Gao, Y. F., Liu, J. K., Yang, F. M., Zhang, G. Y., Wang, D., Zhang, L., et al. (2020). The WRKY transcription factor WRKY8 promotes resistance to pathogen infection and mediates drought and salt stress tolerance in Solanum lycopersicum. Physiol. Plant 168, 98–117. doi: 10.1111/ppl.12978
Hao, Q., Wang, W., Han, X., Wu, J., Lyu, B., Chen, F., et al. (2018). Isochorismate-based salicylic acid biosynthesis confers basal resistance to Fusarium graminearum in barley. Mol. Plant Pathol. 19, 1995–2010. doi: 10.1111/mpp.12675
Hu, Y., Dong, Q., and Yu, D. (2012). Arabidopsis WRKY46 coordinates with WRKY70 and WRKY53 in basal resistance against pathogen Pseudomonas syringae. Plant Sci. 185-186, 288–297. doi: 10.1016/j.plantsci.2011.12.003
Huh, S. U., Choi, L. M., Lee, G. J., Kim, Y. J., and Paek, K. H. (2012). Capsicum annuum WRKY transcription factor d (CaWRKYd) regulates hypersensitive response and defense response upon Tobacco mosaic virus infection. Plant Sci. 197, 50–58. doi: 10.1016/j.plantsci.2012.08.013
Hwang, S. H., Kwon, S. I., Jang, J. Y., Fang, I. L., Lee, H., Choi, C., et al. (2016). OsWRKY51, a rice transcription factor, functions as a positive regulator in defense response against Xanthomonas oryzae pv. oryzae. Plant Cell Rep. 35, 1975–1985. doi: 10.1007/s00299-016-2012-0
Ishiguro, S., and Nakamura, K. (1994). Characterization of a cDNA encoding a novel DNA-binding protein, SPF1, that recognizes SP8 sequences in the 5′ upstream regions of genes coding for sporamin and beta-amylase from sweet potato. Mol. Gen. Genet. 244, 563–571. doi: 10.1007/BF00282746
Ishihama, N., and Yoshioka, H. (2012). Post-translational regulation of WRKY transcription factors in plant immunity. Curr. Opin. Plant Biol. 15, 431–437. doi: 10.1016/j.pbi.2012.02.003
Jakoby, M., Weisshaar, B., Droge-Laser, W., Vicente-Carbajosa, J., Tiedemann, J., Kroj, T., et al. (2002). bZIP transcription factors in Arabidopsis. Trends Plant Sci. 7, 106–111. doi: 10.1016/s1360-1385(01)02223-3
Kumar, S., Stecher, G., and Tamura, K. (2016). MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol. Biol. Evol. 33, 1870–1874. doi: 10.1093/molbev/msw054
Li, J., Brader, G., and Palva, E. T. (2004). The WRKY70 transcription factor: a node of convergence for jasmonate-mediated and salicylate-mediated signals in plant defense. Plant Cell 16, 319–331. doi: 10.1105/tpc.016980
Li, J., Han, G., Sun, C., and Sui, N. (2019). Research advances of MYB transcription factors in plant stress resistance and breeding. Plant Signal Behav. 14:1613131. doi: 10.1080/15592324.2019.1613131
Li, J., Xiang, C. Y., Yang, J., Chen, J. P., and Zhang, H. M. (2015). Interaction of HSP20 with a viral RdRp changes its subcellular localization and distribution pattern in plants. Sci. Rep. 5:14016. doi: 10.1038/srep14016
Li, S., Fu, Q., Chen, L., Huang, W., and Yu, D. (2011). Arabidopsis thaliana WRKY25, WRKY26, and WRKY33 coordinate induction of plant thermotolerance. Planta 233, 1237–1252. doi: 10.1007/s00425-011-1375-2
Li, T., Huang, Y., Xu, Z. S., Wang, F., and Xiong, A. S. (2019). Salicylic acid-induced differential resistance to the Tomato yellow leaf curl virus among resistant and susceptible tomato cultivars. BMC Plant Biol. 19:173. doi: 10.1186/s12870-019-1784-0
Li, W., Zhao, Y., Liu, C., Yao, G., Wu, S., Hou, C., et al. (2012). Callose deposition at plasmodesmata is a critical factor in restricting the cell-to-cell movement of Soybean mosaic virus. Plant Cell Rep. 31, 905–916. doi: 10.1007/s00299-011-1211-y
Liang, X., Chen, X., Li, C., Fan, J., and Guo, Z. (2017). Metabolic and transcriptional alternations for defense by interfering OsWRKY62 and OsWRKY76 transcriptions in rice. Sci. Rep. 7:2474. doi: 10.1038/s41598-017-02643-x
Liu, N., Cheng, X., Lu, R., and Liao, Q. (2014). Construction of Agrobacterium-mediated infectious clone and its expression of Tomato mosaic virus (ToMV) Strain. J. Agric. Biotechnol. 22, 1027–1034.
Love, A. J., Yun, B. W., Laval, V., Loake, G. J., and Milner, J. J. (2005). Cauliflower mosaic virus, a compatible pathogen of Arabidopsis, engages three distinct defense-signaling pathways and activates rapid systemic generation of reactive oxygen species. Plant Physiol. 139, 935–948. doi: 10.1104/pp.105.066803
Lu, Y., Yan, F., Guo, W., Zheng, H., Lin, L., Peng, J., et al. (2011). Garlic virus X 11-kDa protein granules move within the cytoplasm and traffic a host protein normally found in the nucleolus. Mol. Plant Pathol. 12, 666–676. doi: 10.1111/j.1364-3703.2010.00699
Ma, L., Zhang, X., Dou, D., and Cai, C. (2016). Functional analysis of NbWRKY40 transcription factors of Nicotiana benthamiana. Acta Phytopathol. Sin. 46, 791–802.
Miller, R. N., Costa Alves, G. S., and Van Sluys, M. A. (2017). Plant immunity: unravelling the complexity of plant responses to biotic stresses. Ann. Bot. 119, 681–687. doi: 10.1093/aob/mcw284
Murphy, A. M., and Carr, J. P. (2002). Salicylic acid has cell-specific effects on tobacco mosaic virus replication and cell-to-cell movement. Plant Physiol. 128, 552–563. doi: 10.1104/pp.010688
Nan, H., and Gao, L. Z. (2019). Genome-wide analysis of WRKY Genes and their response to hormone and mechanic stresses in carrot. Front. Genet. 10:363. doi: 10.3389/fgene.2019.00363
Pandey, S. P., and Somssich, I. E. (2009). The role of WRKY transcription factors in plant immunity. Plant Physiol. 150, 1648–1655. doi: 10.1104/pp.109.138990
Pillai, S. E., Kumar, C., Patel, H. K., and Sonti, R. V. (2018). Overexpression of a cell wall damage induced transcription factor, OsWRKY42, leads to enhanced callose deposition and tolerance to salt stress but does not enhance tolerance to bacterial infection. BMC Plant Biol. 18:177. doi: 10.1186/s12870-018-1391-5
Rushton, P. J., Somssich, I. E., Ringler, P., and Shen, Q. J. (2010). WRKY transcription factors. Trends Plant Sci. 15, 247–258. doi: 10.1016/j.tplants.2010.02.006
Scarpeci, T. E., Zanor, M. I., Mueller-Roeber, B., and Valle, E. M. (2013). Overexpression of AtWRKY30 enhances abiotic stress tolerance during early growth stages in Arabidopsis thaliana. Plant Mol. Biol. 83, 265–277. doi: 10.1007/s11103-013-0090-8
Shaw, J., Yu, C., Makhotenko, A. V., Makarova, S. S., Love, A. J., Kalinina, N. O., et al. (2019). Interaction of a plant virus protein with the signature Cajal body protein coilin facilitates salicylic acid-mediated plant defence responses. New Phytol. 224, 439–453. doi: 10.1111/nph.15994
Shen, H., Liu, C., Zhang, Y., Meng, X., Zhou, X., Chu, C., et al. (2012). OsWRKY30 is activated by MAP kinases to confer drought tolerance in rice. Plant Mol. Biol. 80, 241–253. doi: 10.1007/s11103-012-9941-y
Shen, L., Yang, S., Yang, T., Liang, J., Cheng, W., Wen, J., et al. (2016). CaCDPK15 positively regulates pepper responses to Ralstonia solanacearum inoculation and forms a positive-feedback loop with CaWRKY40 to amplify defense signaling. Sci. Rep. 6:22439. doi: 10.1038/srep22439
Shi, W., Hao, L., Li, J., Liu, D., Guo, X., and Li, H. (2014). The Gossypium hirsutum WRKY gene GhWRKY39-1 promotes pathogen infection defense responses and mediates salt stress tolerance in transgenic Nicotiana benthamiana. Plant Cell Rep. 33, 483–498. doi: 10.1007/s00299-013-1548-5
Song, L., Huang, S. C., Wise, A., Castanon, R., Nery, J. R., Chen, H., et al. (2016). A transcription factor hierarchy defines an environmental stress response network. Science 354:aag1550. doi: 10.1126/science.aag1550
Tiwari, S. B., Belachew, A., Ma, S. F., Young, M., Ade, J., Shen, Y., et al. (2012). The EDLL motif: a potent plant transcriptional activation domain from AP2/ERF transcription factors. Plant J. 70, 855–865. doi: 10.1111/j.1365-313X.2012.04935.x
Ullah, C., Tsai, C. J., Unsicker, S. B., Xue, L., Reichelt, M., Gershenzon, J., et al. (2019). Salicylic acid activates poplar defense against the biotrophic rust fungus Melampsora larici-populina via increased biosynthesis of catechin and proanthocyanidins. New Phytol. 221, 960–975. doi: 10.1111/nph.15396
van Verk, M. C., Bol, J. F., and Linthorst, H. J. (2011). WRKY transcription factors involved in activation of SA biosynthesis genes. BMC Plant Biol. 11:89. doi: 10.1186/1471-2229-11-89
Wang, X., Sager, R., Cui, W., Zhang, C., Lu, H., and Lee, J. Y. (2013). Salicylic acid regulates Plasmodesmata closure during innate immune responses in Arabidopsis. Plant Cell 25, 2315–2329. doi: 10.1105/tpc.113.110676
White, R. F. (1979). Acetylsalicylic acid (aspirin) induces resistance to tobacco mosaic virus in tobacco. Virology 99, 410–412. doi: 10.1016/0042-6822(79)90019-9
Whitham, S. A., Quan, S., Chang, H. S., Cooper, B., Estes, B., Zhu, T., et al. (2003). Diverse RNA viruses elicit the expression of common sets of genes in susceptible Arabidopsis thaliana plants. Plant J. 33, 271–283. doi: 10.1046/j.1365-313x.2003.01625.x
Wildermuth, M. C., Dewdney, J., Wu, G., and Ausubel, F. M. (2001). Isochorismate synthase is required to synthesize salicylic acid for plant defence. Nature 414, 562–565. doi: 10.1038/35107108
Willems, E., Leyns, L., and Vandesompele, J. (2008). Standardization of real-time PCR gene expression data from independent biological replicates. Anal. Biochem. 379, 127–129. doi: 10.1016/j.ab.2008.04.036
Wong, C. E., Carson, R. A., and Carr, J. P. (2002). Chemically induced virus resistance in Arabidopsis thaliana is independent of pathogenesis-related protein expression and the NPR1 gene. Mol. Plant Microbe Interact. 15, 75–81. doi: 10.1094/MPMI.2002.15.1.75
Wu, X., Shiroto, Y., Kishitani, S., Ito, Y., and Toriyama, K. (2009). Enhanced heat and drought tolerance in transgenic rice seedlings overexpressing OsWRKY11 under the control of HSP101 promoter. Plant Cell Rep. 28, 21–30. doi: 10.1007/s00299-008-0614-x
Wu, Y., Zhang, D., Chu, J. Y., Boyle, P., Wang, Y., Brindle, I. D., et al. (2012). The Arabidopsis NPR1 protein is a receptor for the plant defense hormone salicylic acid. Cell Rep. 1, 639–647. doi: 10.1016/j.celrep.2012.05.008
Xu, X., Chen, C., Fan, B., and Chen, Z. (2006). Physical and functional interactions between pathogen-induced Arabidopsis WRKY18, WRKY40, and WRKY60 transcription factors. Plant Cell 18, 1310–1326. doi: 10.1105/tpc.105.037523
Yang, J., Zhang, F., Xie, L., Song, X. J., Li, J., Chen, J. P., et al. (2016). Functional identification of two minor capsid proteins from Chinese wheat mosaic virus using its infectious full-length cDNA clones. J. Gen. Virol. 97, 2441–2450. doi: 10.1099/jgv.0.000532
Yang, J., Zhang, T. Y., Liao, Q. S., He, L., Li, J., Zhang, H. M., et al. (2018). Chinese wheat mosaic virus-induced gene silencing in monocots and dicots at low temperature. Front. Plant Sci. 9:1627. doi: 10.3389/fpls.2018.01627
Yang, L., Ji, W., Zhu, Y., Gao, P., Li, Y., Cai, H., et al. (2010). GsCBRLK, a calcium/calmodulin-binding receptor-like kinase, is a positive regulator of plant tolerance to salt and ABA stress. J. Exp. Bot. 61, 2519–2533. doi: 10.1093/jxb/erq084
Yang, L., Xu, Y., Liu, Y., Meng, D., Jin, T., and Zhou, X. (2016). HC-Pro viral suppressor from tobacco vein banding mosaic virus interferes with DNA methylation and activates the salicylic acid pathway. Virology 497, 244–250. doi: 10.1016/j.virol.2016.07.024
Yang, Y., Li, R., and Qi, M. (2000). In vivo analysis of plant promoters and transcription factors by agroinfiltration of tobacco leaves. Plant J. 22, 543–551. doi: 10.1046/j.1365-313x.2000.00760.x
Yokotani, N., Sato, Y., Tanabe, S., Chujo, T., Shimizu, T., Okada, K., et al. (2013). WRKY76 is a rice transcriptional repressor playing opposite roles in blast disease resistance and cold stress tolerance. J. Exp. Bot. 64, 5085–5097. doi: 10.1093/jxb/ert298
Yu, D., Chen, C., and Chen, Z. (2001). Evidence for an important role of WRKY DNA binding proteins in the regulation of NPR1 gene expression. Plant Cell 13, 1527–1540. doi: 10.1105/tpc.010115
Yu, F., Huaxia, Y., Lu, W., Wu, C., Cao, X., and Guo, X. (2012). GhWRKY15, a member of the WRKY transcription factor family identified from cotton (Gossypium hirsutum L.), is involved in disease resistance and plant development. BMC Plant Biol. 12:144. doi: 10.1186/1471-2229-12-144
Zhang, Y., and Li, X. (2019). Salicylic acid: biosynthesis, perception, and contributions to plant immunity. Curr. Opin. Plant Biol. 50, 29–36. doi: 10.1016/j.pbi.2019.02.004
Zhang, Y., Yu, H., Yang, X., Li, Q., Ling, J., Wang, H., et al. (2016). CsWRKY46, a WRKY transcription factor from cucumber, confers cold resistance in transgenic-plant by regulating a set of cold-stress responsive genes in an ABA-dependent manner. Plant Physiol. Biochem. 108, 478–487. doi: 10.1016/j.plaphy.2016.08.013
Zhao, P., Yao, X., Cai, C., Li, R., Du, J., Sun, Y., et al. (2019). Viruses mobilize plant immunity to deter nonvector insect herbivores. Sci. Adv. 5:eaav9801. doi: 10.1126/sciadv.aav9801
Zhou, T., Murphy, A. M., Lewsey, M. G., Westwood, J. H., Zhang, H. M., Gonzalez, I., et al. (2014). Domains of the cucumber mosaic virus 2b silencing suppressor protein affecting inhibition of salicylic acid-induced resistance and priming of salicylic acid accumulation during infection. J. Gen. Virol. 95(Pt 6), 1408–1413. doi: 10.1099/vir.0.063461-0
Keywords: NbWRKY40, salicylic acid, VIGS, tomato mosaic virus, transcription factor
Citation: Jiang Y, Zheng W, Li J, Liu P, Zhong K, Jin P, Xu M, Yang J and Chen J (2021) NbWRKY40 Positively Regulates the Response of Nicotiana benthamiana to Tomato Mosaic Virus via Salicylic Acid Signaling. Front. Plant Sci. 11:603518. doi: 10.3389/fpls.2020.603518
Received: 07 September 2020; Accepted: 01 December 2020;
Published: 15 January 2021.
Edited by:
Zuhua He, Center for Excellence in Molecular Plant Sciences, Chinese Academy of Sciences, ChinaReviewed by:
Honghui Lin, Sichuan University, ChinaKen Komatsu, Tokyo University of Agriculture and Technology, Japan
Copyright © 2021 Jiang, Zheng, Li, Liu, Zhong, Jin, Xu, Yang and Chen. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jianping Chen, jpchen2001@126.com