Next Article in Journal
Serological Humoral Immunity Following Natural Infection of Children with High Burden Gastrointestinal Viruses
Next Article in Special Issue
Coinfection of Cotton Plants with Watermelon Mosaic Virus and a Novel Polerovirus in China
Previous Article in Journal
RNA-Seq Analysis of Influenza A Virus-Induced Transcriptional Changes in Mice Lung and Its Possible Implications for the Virus Pathogenicity in Mice
Previous Article in Special Issue
Pepper Mottle Virus and Its Host Interactions: Current State of Knowledge
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Viruses
Volume 13
Issue 10
10.3390/v13102032
viruses-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Expression of an Antiviral Gene GmRUN1 from Soybean Is Regulated via Intron-Mediated Enhancement (IME)

by
Pengfei Diao
1,†,
Hongyu Sun
1,†,
Zhuo Bao
1,
Wenxia Li
1,
Niu Niu
1,
Weimin Li
2 and
Hada Wuriyanghan
1,*
1
Key Laboratory of Forage and Endemic Crop Biology, Ministry of Education, School of Life Sciences, Inner Mongolia University, Hohhot 010070, China
2
Biotechnology Research Institute, Chinese Academy of Agricultural Sciences, Beijing 100081, China
*
Author to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Viruses 2021, 13(10), 2032; https://0-doi-org.brum.beds.ac.uk/10.3390/v13102032
Submission received: 31 August 2021 / Revised: 28 September 2021 / Accepted: 29 September 2021 / Published: 8 October 2021
(This article belongs to the Special Issue State-of-the-Art Plant-Virus Interactions in Asia)
Browse Figures
Versions Notes

Abstract

:
Most of R (resistance) genes encode the protein containing NBS-LRR (nucleotide binding site and leucine-rich repeat) domains. Here, N. benthamiana plants were used for transient expression assays at 3–4 weeks of age. We identified a TNL (TIR-NBS-LRR) encoding gene GmRUN1 that was resistant to both soybean mosaic virus (SMV) and tobacco mosaic virus (TMV). Truncation analysis indicated the importance of all three canonical domains for GmRUN1-mediated antiviral activity. Promoter-GUS analysis showed that GmRUN1 expression is inducible by both salicylic acid (SA) and a transcription factor GmDREB3 via the cis-elements as-1 and ERE (ethylene response element), which are present in its promoter region. Interestingly, GmRUN1 gDNA (genomic DNA) shows higher viral resistance than its cDNA (complementary DNA), indicating the existence of intron-mediated enhancement (IME) for GmRUN1 regulation. We provided evidence that intron2 of GmRUN1 increased the mRNA level of native gene GmRUN1, a soybean antiviral gene SRC7 and also a reporter gene Luciferase, indicating the general transcriptional enhancement of intron2 in different genes. In summary, we identified an antiviral TNL type soybean gene GmRUN1, expression of which was regulated at different layers. The investigation of GmRUN1 gene regulatory network would help to explore the mechanism underlying soybean-SMV interactions.

1. Introduction

Two layers of innate immune systems have evolved to recognize the potential pathogens and initiate an effective defense response. The first type of immune response is initiated by the pattern recognition receptors (PRRs) localized at the plasma membrane [1]. PRR recognizes and responds to evolutionarily conserved pathogen-associated molecular pattern (PAMP), and it is called pattern-triggered immunity (PTI) [2]. Some pathogens secrete virulence effectors to counteract PTI. Plants activate the second type of innate immune system to recognize the virulence effectors, named the effector-triggered immunity (ETI) [3]. ETI is generally mediated by a resistance (R) gene and leads to local necrosis of plants to limit the continuous spread of pathogens, namely hypersensitive response (HR) [4]. Nucleotide-binding site (NBS), leucine-rich repeat (LRR)-containing proteins (NLRs) occupy the largest proportion in plant R proteins [5]. NBS-LRR genes belong to a large gene family, with hundreds of copies in the genome, and are distributed in obvious uneven clusters [6,7]. To date, many NBS-LRR type R genes have been cloned from different plant species [8].
Plant NLR proteins belong to signal transduction ATPases with numeric domains (STAND) superfamily [9]. The central NBS domain performs the function of molecular switch and controls the binding state of ATP/ADP to mediate downstream signal transduction [10,11]. Leucine-rich repeats (LRRs) are ubiquitously present protein domains involved in mediating protein–protein interactions [12,13]. Some studies suggest that LRR motifs can give recognition specificity in plant defense response [14]. According to their different N-terminal structures, these NBS-LRR proteins can be further divided into two categories: TIR-NBS-LRR (TNL) proteins with the homologous domain of toll/interleukin-1 receptor (TIR) and non-TNL (nTNL) proteins [15]. Most nTNL type-R proteins have a coiled coil (CC) structure at the N-terminal, commonly known as CC-NBS-LRR (CNL) type-R protein [16,17]. Generally, a large number of CNL genes are found in all plant genomes, but TNL genes are not identified in monocotyledons [18]. The N-terminal CC or TIR domain can be used as a signal transduction center, which is associated with cellular targets or downstream signaling components of effectors [19].
The expression of plant R gene needs a strict regulation mechanism, and its ability to activate defense signal and trigger immunity depends on its protein level [20]. Overexpression of R gene can lead to autoimmunity and even plant growth retardation [21]. Therefore, the precise regulation of R proteins in plant homeostasis, including transcriptional and translational regulation, is crucial for plant growth and plant disease resistance. The expression of R genes is strictly regulated at multiple steps including transcription, post-transcriptional processing, and transcript turnover [22].
Soybean mosaic virus (SMV) is one of the main members of potyvirus [23]. The infection of SMV causes mosaic, necrosis, and other symptoms in many soybean varieties by means of aphid and seed transmission [24]. Therefore, SMV is a major disease that seriously threatens the yield and quality of soybean [25,26,27,28]. SMV genome is a single-stranded sense RNA which encode eleven functional proteins: P1, HC-Pro, P3, PIPO, 6K1, CI, 6k2, NIa-vpg, NIa-Pro, NIb, and CP [29,30]. Three independent SMV resistance loci, Rsv1, Rsv3, and Rsv4, were identified from soybean. Rsv1 is located on chromosome 13, which may contain one or more members of the NBS-LRR gene family and is highly resistant to most SMV strains except G7 [31,32]. Rsv3 locus was located on chromosome 14 and was resistant to strains G5, G6, and G7 [33]. Rsv4 is resistant to strains v94-5152 and encodes a SMV specific dsRNase [34,35]. In our previous study, we characterized dozens of SMV-responsive NLR genes in the susceptible soybean variety Hefeng25 by transcriptome sequencing [36]. Here, we identified a SMV-resistant soybean gene GmRUN1 encoding a typical TNL protein. GmRUN1 also showed resistance to TMV (tobacco mosaic virus) using transient expression assays in Nicotiana benthamiana, which is the most widely used experimental host in plant virology, due mainly to the large number of diverse plant viruses that can successfully infect it [37]. GmRUN1 genomic DNA (gDNA) is more resistant to these two viruses than its cDNA (complementary DNA), indicating the existence of intron-mediated enhancement (IME) for GmRUN1 regulation. Further analysis showed that intron2 of GmRUN1 was responsible for transcriptional enhancement of GmRUN1. GmRUN1 expression is also inducible by plant hormone SA and a transcription factor GmDREB3.

2. Materials and Methods

2.1. Plant Growth Conditions

Soybean [Glycine max (L.) Merr.] and Nicotiana benthamiana plants were grown in a glasshouse under a 14 h light/10 h dark cycle (24 °C day/22 °C night). Plants were used for transient expression assays at 3–4 weeks of age.

2.2. Construction of Recombinant Vectors

To generate cDNA clones of soybean GmRUN1 and GmDREB3 genes, total RNA was isolated from soybean cv. Hefeng25 leaves using the TRIzol reagent (Invitrogen, Waltham, MA, USA, Cat#15596026), and cDNA was generated using the GoScript reverse-transcription system (Promega, Madison, WI, USA, Cat#A5001) following the manufacturers’ instructions, and the sequences were then amplified from this cDNA using the primers listed in Supplemental Table S1 with PrimeSTAR® GXL DNA Polymerase (Takara, Kusatsu, Japan, Cat#R050A). For full-length GmRUN1 gene, PCR product was cloned into the binary vector pBI121 digested by SmaI and SacΙ using the In-Fusion HD cloning kit (Takara, Cat#639650). For GmDREB3 gene, PCR product was cloned into the TA cloning vector using the pMD19-T vector cloning kit (Takara, Cat#3271) and then recombined with the pMD1-T7 vector digested by BamHI and XhoΙ using the ClonExpress Ultra one-step cloning kit (Vazyme, Nanjing, China, Cat#C115-01). For truncated domain fragments of GmRUN1 gene, three different domains TIR (1–165 aa), NBS (166–474 aa), LRR (475–1088 aa) were identified by SMART website (http://smart.embl-heidelberg.de/ (Accessed: 1 October 2021)), and the PCR product was cloned into the pCB301-2μ-HDV vector linearized by PCR.
To generate truncated promotor clones of GmRUN1 gene, genomic DNA was isolated from soybean cv. Hefeng25 leaves using the CTAB method according to the manufacturer’s instructions, and the GmRUN1 truncated promoter regions were amplified and cloned into the pBI121 vector by replacing the CaMV 35S promoter before the β-glucuronidase (GUS) gene.
To generate exogenous intron insert clones of GmRUN1 gene, the Luciferase gene was cloned from pGWB435-LUC (GenBank No. AB294455.1) and then recombined into the entry vector pHSG299. Four introns of GmRUN1 gene were amplified from the genomic DNA and cloned into recombined vector pHSG299-LUC digested by ScaI site. The recombined entry clones were then recombined with the binary vector pCambia1300 for luciferase assay.

2.3. Transient Expression and Virus Inoculation

Agrobacterium tumefaciens GV3101 strains carrying recombinant binary vectors wereused to infiltrate N. benthamiana leaves. Liquid cultures of all Agrobacterium strains were initially grown at 28 °C with agitation in Luria–Bertani (LB) media supplemented with the appropriate antibiotics. The bacterial cells were pelleted by centrifugation at 8000 rpm for 1 min, resuspended in infiltration buffer (10 mM MgCl2, 10 mM MES, and 200 μM AS) and adjusted to an appropriate OD600 for infiltration. PJL24 which carried GFP (green fluorescent protein) in TMV genome was used as infectious clones for verification of GmRUN1 resistance, and GFP fluorescence was detected by a handheld long-wave (365 nm) UV lamp. SMV-N1 strain was used to infect N. benthamiana leaves with mechanical inoculation, the SMV-infected leaves under quartz sand grinding in 1×phosphate buffer was daubed to infiltrated site using a writing brush.

2.4. Promoter Analysis and In Situ GUS Activity Assay

Promoter elements were predicted for the 3000 bp genomic sequence upstream of GmRUN1 gene by PlantCare website (http://bioinformatics.psb.ugent.be/webtools/plantcare/html/ (Accessed: 1 October 2021)). The promoter-GUS recombinant vectors, alone or together with the recombinant vectors expressing GmDREB3 proteins, were infiltrated into N. benthamiana. SMV and SA were also applied to the leaves for induction analysis by rub and spray, respectively, and the leaf discs (1 cm in diameter) were cut at 2 dpi (day post inoculation) for GUS staining with X-Gluc as the substrate according to the literature [38].

2.5. Luciferase Reporter Assay and Fluorescence Quantitative Analysis

Agrobacterium strain GV3101 carrying LUC gene inserted different GmRUN1 introns that recombined in pCambia1300 expression vector were infiltrated into N. benthamiana leaves. After 24–48 h, 20 μL 0.5 mM D-luciferin was applied to infiltrated site in the dark. After dark treatment for 3–5 min, the LUC expression was detected by CCD imaging system (Berthold Technologies, Bad Wildbad, Germany, LB 985) with IndiGO software at 560 nm for exposure time of 1–3 min, and the fluorescence intensity was visualized to assay the effect of different GmRUN1 introns. The infiltrated sites of leaves without imaging in the same batch were cut and quick freezing in liquid nitrogen for qRT-PCR analysis. To quantitatively analyze the level of GmRUN1 resistance, Gel-Pro analyzer software was used for the detection of TMV-GFP intensity in the green channel, and the data were imported to GraphPad Prism7 software to analyze statistically significant and draw graph.

2.6. Quantitative Real-Time PCR (qRT-PCR) Analysis

Total RNA extraction and reverse transcription were performed as described above, and the quantitative RT-PCR using gene-specific primers was carried out subsequently in an Analytikjena qPCR instrument using TransStrat® Tip Green qPCR SuperMix (TransGen Biotech, Beijing, China, Cat#AQ141). The data were normalized to ACTIN expression by the cycle threshold (CT) 2-ΔΔCT method according to the literature [39] and analyzed by Graphpad Prism7 software. All experiments were repeated at least three times. Primers used in this study are listed in Table S1.

3. Results

3.1. GmRUN1 cDNA Shows Partial Resistance to TMV and SMV in N. Benthamiana Transient Expression Assay

Previously, we identified NBS-LRR family genes that were involved in SMV–soybean interactions [36]. In the present study, we characterized one such gene with sequence ID of XM_006592417.3 (gene locus: G12g132200) and designated it as GmRUN1 as it was homologous with Vitis rotundifolia RUN1 (resistance to Uncinula Necator 1) gene (Figure 1A,B) [40,41]. ORF (open-reading frame) sequence or genomic DNA sequence of GmRUN1 was amplified from soybean cDNA or genomic DNA and were ligated into binary vectors pCB301 and pBI121 to obtain recombinant overexpression vectors. We used N. benthamiana transient expression system to investigate the role of GmRUN1 for SMV resistance. In the parallel experiment, TMV-GFP infectious clone (pJL24) was also used to investigate the antiviral role of GmRUN1 toward TMV [42]. The above recombinant vectors were transformed into Agrobacterium strain GV3101 and transiently expressed in tobacco leaves to detect their resistance to TMV and SMV. In this study, a TMV-resistant tobacco gene N [43] and a SMV/TMV resistant soybean gene SRC7 (unpublished data from the same lab) were used as positive control, while Agrobacterium carrying empty vector pBI121 was used as negative control. We defined almost no GFP fluorescence (for TMV) or severe hypersensitive response (HR) (for SMV) as full resistance, weak GFP fluorescence or mild HR as partial resistance, and strong GFP fluorescence or no HR as no resistance. As a result, GmRUN1 genomic DNA showed full resistance to both TMV and SMV, while its cDNA displayed only partial resistance when it was expressed from either of the binary vectors (Figure 1C–E; Table 1 and Table 2).

3.2. Three Canonical Domains of GmRUN1 Are Indispensable for Its Antiviral Activity

GmRUN1 is a typical TNL (TIR-NBS-LRR) protein, containing three domains, including a typical N-terminal TIR domain, central NBS domain, and C-terminal LRR domain (Figure 2A). To investigate the importance of these domains, we made truncation analysis. We expressed five truncations for GmRUN1, including GmRUN1TIR, GmRUN1NBS, GmRUN1LRR, GmRUN1TN (abbreviation for TIR-NBS of GmRUN1), and GmRUN1NL (abbreviation for NBS-LRR of GmRUN1) (Figure 2A). Truncation analysis showed that all of the TIR, NBS, and LRR domains were indispensable for GmRUN1 antiviral activity, as the deletion of any of them abolished its resistance to both TMV and SMV (Figure 2B–D; Table 3 and Table 4).

3.3. GmRUN1 Has IME Phenomena

As GmRUN1 gDNA is resistant to TMV and SMV and its cDNA without the introns has been proved to be partially resistant, we speculate that the introns of GmRUN1 might be involved in its antiviral activity. GmRUN1 gDNA contains four introns with the sizes of 6738, 243, 134, and 102 nt, and all of them follow the “GT-AG” rule (Figure 1A). Some introns were reported to enhance gene expression, and this phenomenon was named intron-mediated enhancement (IME) [44,45]. To determine the possible IME of four GmRUN1 introns, we inserted them at the +166 site of luciferase (LUC) reporter gene (Figure 3A). Interestingly, the insertion of any of four introns from GmRUN1 abolished LUC signal (Figure 3B). We speculate that the insertion of these introns may disturb the LUC ORF and subsequent protein expression. The PCR amplicons still contained the introns when total RNA was extracted from the inoculated leaves, and LUC cDNA was amplified, indicating that the introns were not correctly spliced from LUC mRNA (Figure 3C). Therefore, these introns led to insertion mutation in LUC gene, showing that the introns of one gene may not be spliced normally when was inserted into other genes. However, we observed that LUC mRNA level was significantly upregulated by 3.3 times upon insertion of intron2, showing that intron2 has IME effects at transcriptional level (Figure 3D). We then replaced the intron of soybean antiviral gene SRC7 with GmRUN1 intron2 (Figure 4A), and the mRNA level of SRC7 was increased by 2.6 folds, further demonstrating the IME effect of GmRUN1 intron2 at the transcriptional level (Figure 4B). Furthermore, the insertion of GmRUN1 intron2 did not abolish SRC7 antiviral activity, indicating that GmRUN1 intron2 was correctly spliced from SRC7 mRNA (Figure 4C; Table 5). It also indicates that the splicing of intron depends on inserted genes. SRC7 was homologous to GmRUN1; therefore, it is reasonable that GmRUN1 intron2 can be spliced in SRC7 but not in nonhomologous LUC gene. To further demonstrate the importance of these introns, we made GmRUN1 expression constructs with truncations in different introns (Figure 5A). The depletion of intron1 increased while further deletion of intron2 decreased antiviral activity, indicating the enhancement of intron2 for antiviral activity (Figure 5B, Table 6). Taken together, GmRUN1 intron2 has an IME effect.

3.4. GmRUN1 Expression Is Transcriptionally Regulated by SA

Having confirmed the IME effect of intron2 on GmRUN1 expression, we next intended to clone the GmRUN1 promoter and further examine the transcriptional regulation of GmRUN1. We used GUS reporter to assess promoter activity. The promoter region of GmRUN1 was amplified from soybean genomic DNA and was cloned into binary vector pBI121 to obtain Pro:GUS reporter vectors. We made four different Pro:GUS constructs, harboring different lengths of the promoter region, namely Pro2415:GUS, Pro2592:GUS, Pro2237:GUS, and Pro2414: GUS (Figure 6A). When these Pro:GUS constructs are transiently expressed in N. benthamiana, they did not show any GUS signal, demonstrating that they possessed very low basal transcriptional activity (Figure 6B). The infection of SMV did not induce the GUS expression in any of the Pro:GUS constructs, indicating that GmRUN1 expression is not inducible by SMV infection (Figure 6B). Salicylic acid (SA) is a well-known defense hormone which is generally implicated in plant immunity against plant viruses including SMV [36]. Treatment with MeSA, an analog of SA, elevated GUS expression for the Pro2592:GUS and Pro2414:GUS reporters but not in Pro2415:GUS and Pro2237:GUS reporters (Figure 6B).

3.5. GmRUN1 Expression Is Transcriptionally Induced by Transcription Factor GmDREB3

Based on the above results, we further analyzed cis-regulatory element present in GmRUN1 promoter using PlantCare website (http://bioinformatics.psb.ugent.be/webtools/plantcare/html/ (Accessed: 1 October 2021)). Interestingly, we found two cis-elements immediately upstream of GmRUN1 ORF which were present in SA-inducible constructs Pro2592:GUS and Pro2414:GUS, while lacking in SA noninducible Pro2415:GUS and Pro2237:GUS constructs (Figure 6A). One cis element is the SA-responsive as-1, which might contribute to SA induction of GmRUN1 promoter. Another cis element is the ERF binding element ERE, which is close to as-1 element on GmRUN1 promoter (Table 7). Previously, we showed that GmRUN1 expression is repressed upon SMV infection in soybean [36]. From our RNA-seq data, we also found that GmDREB3 is significantly repressed by SMV infection, and GmDREB3 expression pattern is positively correlated with that of GmRUN1 (Figure 7A). The data lead us to the assumption that GmDREB3 might be a positive regulator of GmRUN1. Coexpression of GmDREB3 elevated GUS expression for the Pro2592:GUS and Pro2414:GUS reporters but not in Pro2415:GUS and Pro2237:GUS reporters, further demonstrating positive regulation of GmRUN1 expression by GmDREB3 via the ERE cis element at GmRUN1 promoter (Figure 7B).

4. Discussion

In this study, we identified an antiviral gene GmRUN1 from soybean. GmRUN1 encodes a typical TIR-NBS-LRR protein and gives resistance to both a Potyvirus SMV and a Tobamovirus TMV. Truncation analysis showed that the all of the canonical domains TIR, NBS, and LRR were necessary for GmRUN1 antiviral activity. In some reports, TIR-NBS domains but not LRR domain were sufficient to trigger immune responses, especially in transient expression assay. For example, Arabidopsis powdery mildew resistance gene TN2 and autoimmune-related gene CHS1 encode functional TIR-NBS proteins [46,47]. Furthermore, TIR-NBS genes were broadly reported in variety of plant species, including the leguminous plants such as soybean and common bean [48,49]. Therefore, the antiviral mechanism of GmRUN1 might be different from those of TIR-NBS genes and deserves further investigation.
Overexpression of GmRUN1 showed HR upon infection with SMV. Cell-death-triggering activity of R proteins should be under strict control so as to trigger timely immune response only upon pathogen infection and to also avoid fitness costs at pathogen-free conditions [50]. Therefore, multiple layers of regulation at transcriptional, post-transcriptional, and protein activity levels exist for expression control of R genes [22]. The expression of GmRUN1 is inducible by the treatment with major defense hormone SA. Promoter truncation analysis narrowed down the SA-responsive element to a −326~−336 nt (10 nt) region, where we identified a SA-responsive cis element as-1, which is most likely contributed to SA induction on GmRUN1 expression [51,52]. SA is an important defense hormone which contributes to the immunity against various pathogens, especially plant viruses, and leads to systemic acquired resistance (SAR) [53]. For example, SA was reported to be involved in defense response mediated by R genes such as tobacco N and Arabidopsis AtTN10 [48,54]. Besides this as-1 element, we also found a ERE element in the SA-responsive promoter region. A putative transcription factor GmDREB3 upregulates GmRUN1 promoter activity via ERE-dependent manner, and the expression of GmRUN1 is positively correlated with GmDREB3 expression, demonstrating that GmRUN1 expression is transcriptionally induced by GmDREB3. The effects of DREB on R gene expression and involvement in plant immunity were also observed in other plant species [55,56,57].
It has been shown that some introns had positive regulatory roles on gene expression, but some others possessed inhibitory effects [58]. In 1987, Callis et al. first discovered that introns can mediate the enhancement of gene expression in maize cells [44]. Subsequently, the phenomenon of IME was observed in mammals, nematodes, and yeasts [59,60,61]. Some endogenous introns can compensate for the low-level expression driven by the weak promoter [62]. IME is also associated with specific sequence motifs, such as TTNGATYTG and CGATT [63]. Several introns were also shown to have both promoter activity and enhancer function [64]. In most studies, IME was attributable to increased mRNA accumulation at the transcriptional level; however, there are also some data that provide evidence of enhanced translation [65,66]. Here, we provided evidence that intron2 of GmRUN1 increased mRNA level of native gene GmRUN1, a homologous gene SRC7 and also a reporter gene Luciferase, indicating the general transcriptional enhancement of this intron2 in different genes. However, we also showed that correct splicing of intron2 might depend on the inserted genes. The in-depth study of IME phenomenon provides a basis for the wide application of functional introns.

5. Conclusions

Plant viruses pose threats to agriculturally important crops. SMV is a major pathogen of soybean and causes heavy yield losses worldwide. Although several R genes have been cloned from multiple host varieties, soybean–SMV interactions are still elusive. In this study, we identified a TNL (TIR-NBS-LRR)-type antiviral gene GmRUN1 from soybean. GmRUN1 genomic DNA showed full resistance to both TMV and SMV, while its cDNA displayed only partial resistance, and three canonical domains of GmRUN1 are indispensable for its antiviral activity. GmRUN1 is spectacular as its expression is regulated at multiple layers, such as SA induction, GmDREB3 transcriptional activation, and IME from intron2. GmRUN1 represents a novel SMV-resistance gene and deserves further functional study.

Supplementary Materials

The following are available online at www.mdpi.com/article/10.3390/v13102032/s1, Table S1: Primers used in this study.

Author Contributions

Conceptualization, H.W.; data curation, P.D. and H.S.; formal analysis, P.D., H.S. and Z.B.; funding acquisition, H.W.; investigation, P.D., H.S., Z.B., W.L. (Wenxia Li) and N.N.; methodology, H.S.; project administration, H.W.; supervision, W.L. (Weimin Li) and H.W.; visualization, P.D.; writing—original draft, P.D. and H.S.; writing—review and editing, H.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the science and technology program of Inner Mongolia Autonomous Region (2020GG0045, 201802061), innovation guidance project of Inner Mongolia Autonomous Region (KCBJ2018001), Natural Science Foundation of Inner Mongolia (2017MS0330), and National Natural Science Foundation of China (grant No. 31860055, 31460036).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

We thank Bryce Falk for providing pJL24 (TMV-GFP) vector and Zhenghe Li from Zhejiang University for providing pCB301-2µ-HDV vector.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zipfel, C. Pattern-recognition receptors in plant innate immunity. Curr. Opin. Immunol. 2008, 20, 10–16. [Google Scholar] [CrossRef]
  2. Boller, T.; He, S.Y. Innate immunity in plants: An arms race between pattern recognition receptors in plants and effectors in microbial pathogens. Science 2009, 324, 742–744. [Google Scholar] [CrossRef] [Green Version]
  3. Jones, J.D.G.; Dangl, J.L. The plant immune system. Nature 2006, 444, 323–329. [Google Scholar] [CrossRef] [Green Version]
  4. Staskawicz, B.J.; Ausubel, F.M.; Baker, B.J.; Ellis, J.G.; Jones, J. Molecular genetics of plant disease resistance. Science 1995, 268, 661–667. [Google Scholar] [CrossRef] [Green Version]
  5. Kourelis, J.; van der Hoorn, R.A.L. Defended to the Nines: 25 years of resistance gene cloning identifies nine mechanisms for r protein function. Plant Cell 2018, 30, 285–299. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Luo, S.; Zhang, Y.; Hu, Q.; Chen, J.; Li, K.; Lu, C.; Liu, H.; Wang, W.; Kuang, H. Dynamic nucleotide-binding site and leu-cine-rich repeat-encoding genes in the grass family. Plant Physiol. 2012, 159, 197–210. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. He, L.; Du, C.; Covaleda, L.; Xu, Z.; Robinson, A.F.; Yu, J.Z.; Kohel, R.J.; Zhang, H.-B. Cloning, characterization, and evolution of the nbs-lrr-encoding resistance gene analogue family in polyploid cotton (Gossypium hirsutum L.). Mol. Plant-Microbe. Interact. 2004, 17, 1234–1241. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Hulbert, S.H.; Webb, C.A.; Smith, S.M.; Sun, Q. Resistance gene complexes: Evolution and utilization. Annu. Rev. Phytopathol. 2001, 39, 285–312. [Google Scholar] [CrossRef] [PubMed]
  9. Lelpe, D.D.; Koonin, E.V.; Aravind, L. STAND, a class of P-loop NTPases including animal and plant regulators of pro-grammed cell death: Multiple, complex domain architectures, unusual phyletic patterns, and evolution by horizontal gene transfer. J. Mol. Biol. 2004, 343, 1–28. [Google Scholar]
  10. Walker, J.E.; Saraste, M.; Runswick, M.J.; Gay, N.J. Distantly related sequences in the alpha subunits and beta subunits of atp synthase myosin kinases and other atp requiring enzymes and a common nucleotide binding fold. EMBO J. 1982, 1, 945–952. [Google Scholar] [CrossRef]
  11. Tameling, W.I.L.; Elzinga, S.D.J.; Darmin, P.S.; Vossen, J.H.; Takken, F.; Haring, M.A.; Cornelissen, B.J.C. The tomato R gene products I-2 and Mi-1 are functional atp binding proteins with atpase activity. Plant Cell 2002, 14, 2929–2939. [Google Scholar] [CrossRef] [PubMed]
  12. Shanmugam, V. Role of extracytoplasmic leucine rich repeat proteins in plant defence mechanisms. Microbiol. Res. 2005, 160, 83–94. [Google Scholar] [CrossRef]
  13. Velásquez, A.C.; Nomura, K.; Cooper, M.D.; Herrin, B.R.; He, S.Y. Leucine-rich-repeat-containing variable lymphocyte receptors as modules to target plant-expressed proteins. Plant Methods 2017, 13, 29. [Google Scholar] [CrossRef] [Green Version]
  14. Ellis, J.; Dodds, P.; Pryor, T. Structure, function and evolution of plant disease resistance genes. Curr. Opin. Plant Biol. 2000, 3, 278–284. [Google Scholar] [CrossRef]
  15. Eitas, T.K.; Dangl, J.L. NB-LRR proteins: Pairs, pieces, perception, partners, and pathways. Curr. Opin. Plant Biol. 2010, 13, 472–477. [Google Scholar] [CrossRef] [Green Version]
  16. Meyers, B.C.; Dickerman, A.W.; Michelmore, R.W.; Sivaramakrishnan, S.; Sobral, B.W.; Young, N.D. Plant disease resistance genes encode members of an ancient and diverse protein family within the nucleotide-binding superfamily. Plant J. 1999, 20, 317–332. [Google Scholar] [CrossRef] [PubMed]
  17. Meyers, B.; Kozik, A.; Griego, A.; Kuang, H.; Michelmore, R.W. Genome-wide analysis of NBS-LRR–Encoding genes in arabidopsis. Plant Cell 2003, 15, 809–834. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Yue, J.-X.; Meyers, B.C.; Chen, J.-Q.; Tian, D.; Yang, S. Tracing the origin and evolutionary history of plant nucleotide-binding site-leucine-rich repeat (NBS-LRR) genes. New Phytol. 2012, 193, 1049–1063. [Google Scholar] [CrossRef]
  19. Mukhtar, M.S.; Carvunis, A.-R.; Dreze, M.; Epple, P.; Steinbrenner, J.; Moore, J.; Tasan, M.; Galli, M.; Hao, T.; Nishimura, M.T.; et al. Independently Evolved Virulence Effectors Converge onto Hubs in a Plant Immune System Network. Science 2011, 333, 596–601. [Google Scholar] [CrossRef] [Green Version]
  20. Bieri, S.; Mauch, S.; Shen, Q.H.; Peart, J.; Devoto, A.; Casais, C.; Ceron, F.; Schulze, S.; Steinbiss, H.H.; Shirasu, K.; et al. RAR1 positively controls steady state levels of barley MLA resistance proteins and enables sufficient MLA6 ac-cumulation for effective resistance. Plant Cell 2004, 16, 3480–3495. [Google Scholar] [CrossRef] [Green Version]
  21. Mackey, D.; Belkhadir, Y.; Alonso, J.; Ecker, J.; Dangl, J.L. Arabidopsis RIN4 is a target of the type III virulence effector AvrRpt2 and modulates RPS2-mediated resistance. Cell 2003, 112, 379–389. [Google Scholar] [CrossRef] [Green Version]
  22. Lai, Y.; Eulgem, T. Transcript-level expression control of plant NLR genes. Mol. Plant Pathol. 2018, 19, 1267–1281. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Mayo, M.A.; Pringle, C.R. Virus taxonomy-1997. J. Gen. Virol. 1998, 79, 649–657. [Google Scholar] [CrossRef]
  24. Yang, Y.; Lin, J.; Zheng, G.; Zhang, M.; Zhi, H. Recombinant soybean mosaic virus is prevalent in Chinese soybean fields. Arch. Virol. 2014, 159, 1793–1796. [Google Scholar] [CrossRef]
  25. Cho, E.K.; Goodman, R.M. Strains of soybean mosaic virus: Classification based on virulence in resistant soybean cultivars. Phytopathology 1979, 69, 467–470. [Google Scholar] [CrossRef]
  26. Cho, E.K.; Chung, B.J.; Lee, S.H. Studies on identification and classification of soybean virus diseases in Korea part 2 etiology of a necrotic disease of Glycine-max. Plant Dis. Rep. 1977, 61, 313–317. [Google Scholar]
  27. Saruta, M.; Kikuchi, A.; Okabe, A.; Sasaya, T. Molecular characterization of A2 and D strains of Soybean mosaic virus, which caused a recent virus outbreak in soybean cultivar Sachiyutaka in Chugoku and Shikoku regions of Japan. J. Gen. Plant Pathol. 2005, 71, 431–435. [Google Scholar] [CrossRef]
  28. Li, K.; Yang, Q.H.; Zhi, H.J.; Gai, J.Y. Identification and distribution of soybean mosaic virus strains in Southern China. Plant Dis. 2010, 94, 351–357. [Google Scholar] [CrossRef] [Green Version]
  29. Eggenberger, A.L.; Stark, D.M.; Beachy, R.N. The nucleotide sequence of a soybean mosaic virus coat protein-coding region and its expression in escherichia coli, agrobacterium tumefaciens and tobacco callus. J. Gen. Virol. 1989, 70, 1853–1860. [Google Scholar] [CrossRef]
  30. Liu, J.-Z.; Fang, Y.; Pang, H. The current status of the Soybean-Soybean Mosaic Virus (SMV) Pathosystem. Front. Microbiol. 2016, 7, 1906. [Google Scholar] [CrossRef]
  31. Zhang, C.; Grosic, S.; Whitham, S.A.; Hill, J.H. The requirement of multiple defense genes in soybean Rsv1-mediated extreme resistance to soybean mosaic virus. Mol. Plant-Microbe Interact. 2012, 25, 1307–1313. [Google Scholar] [CrossRef]
  32. Hajimorad, M.R.; Hill, J.H. Rsv1-mediated resistance against soybean mosaic virus-N is hypersensitive re-sponse-independent at inoculation site, but has the potential to initiate a hypersensitive response-like mechanism. Mol. Plant-Microbe Interact. 2001, 14, 587–598. [Google Scholar] [CrossRef] [Green Version]
  33. Chen, P.; Buss, G.R.; Tolin, S.A. Resistance to soybean mosaic virus conferred by two independent dominant genes in PI 486355. J. Hered. 1993, 84, 25–28. [Google Scholar] [CrossRef]
  34. Ishibashi, K.; Saruta, M.; Shimizu, T.; Shu, M.; Anai, T.; Komatsu, K.; Yamada, N.; Katayose, Y.; Ishikawa, M.; Ishimoto, M.; et al. Soybean antiviral immunity conferred by dsRNase targets the viral replication complex. Nat. Commun. 2019, 10, 1–10. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Ilut, D.C.; Lipka, A.E.; Jeong, N.; Bae, D.N.; Kim, D.H.; Kim, J.H.; Redekar, N.; Yang, K.; Park, W.; Kang, S.-T.; et al. Identification of haplotypes at the Rsv4 genomic region in soybean associated with durable resistance to soybean mosaic virus. Theor. Appl. Genet. 2016, 129, 453–468. [Google Scholar] [CrossRef] [PubMed]
  36. Zhao, Q.; Li, H.; Sun, H.; Li, A.; Liu, S.; Yu, R.; Cui, X.; Zhang, D.; Wuriyanghan, H.; Su, H. Salicylic acid and broad spectrum of NBS-LRR family genes are involved in SMV-soybean interactions. Plant Physiol. Biochem. 2018, 123, 132–140. [Google Scholar] [CrossRef]
  37. Goodin, M.M.; Zaitlin, D.; Naidu, R.A.; Lommel, S.A. Nicotiana benthamiana: Its history and future as a model for plant–pathogen interactions. Mol. Plant-Microbe Interact. 2014, 21, 1015–1026. [Google Scholar] [CrossRef] [Green Version]
  38. Jefferson, R.A.; Kavanagh, T.A.; Bevan, M.W. GUS fusions: Beta-glucuronidase as a sensitive and versatile gene fusion marker in higher plants. EMBO J. 1987, 6, 3901–3908. [Google Scholar] [CrossRef]
  39. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2(T) (-Delta Delta C) method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef]
  40. Feechan, A.; Kocsis, M.; Riaz, S.; Zhang, W.; Gadoury, D.M.; Walker, M.A.; Dry, I.B.; Reisch, B.; Cadle-Davidson, L. Strategies for RUN1 deployment using RUN2 and REN2 to manage grapevine powdery mildew informed by studies of race specificity. Phytopathology 2015, 105, 1104–1113. [Google Scholar] [CrossRef] [Green Version]
  41. Pauquet, J.; Bouquet, A.; This, P.; Adam-Blondon, A.-F. Establishment of a local map of AFLP markers around the powdery mildew resistance gene Run1 in grapevine and assessment of their usefulness for marker assisted selection. Theor. Appl. Genet. 2001, 103, 1201–1210. [Google Scholar] [CrossRef]
  42. Lindbo, J.A. TRBO: A high-efficiency tobacco mosaic virus RNA-based overexpression vector. Plant Physiol. 2007, 145, 1232–1240. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Whitham, S.; Dinesh-Kumar, S.P.; Choi, D.; Hehl, R.; Corr, C.; Baker, B. The product of the tobacco mosaic virus resistance gene N: Similarity to toll and the interleukin-1 receptor. Cell 1994, 78, 1101–1115. [Google Scholar] [CrossRef]
  44. Callis, J.; Fromm, M.; Walbot, V. Introns increase gene expression in cultured maize cells. Genes Dev. 1987, 1, 1183–1200. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Rose, A.B. Requirements for intron-mediated enhancement of gene expression in Arabidopsis. RNA 2002, 8, 1444–1453. [Google Scholar] [CrossRef] [Green Version]
  46. Zhao, T.; Rui, L.; Li, J.; Nishimura, M.; Vogel, J.; Liu, N.; Liu, S.; Zhao, Y.; Dangl, J.L.; Tang, D. A Truncated NLR protein, TIR-NBS2, is required for activated defense responses in the exo70B1 Mutant. PLoS Genet. 2015, 11, e1004945. [Google Scholar] [CrossRef]
  47. Zhang, Y.; Wang, Y.; Liu, J.; Ding, Y.; Wang, S.; Zhang, X.; Liu, Y.; Yang, S. Temperature-dependent autoimmunity mediated by chs1 requires its neighboring TNL gene SOC 3. New Phytol. 2017, 213, 1330–1345. [Google Scholar] [CrossRef]
  48. Nandety, R.S.; Caplan, J.L.; Cavanaugh, K.; Perroud, B.; Wroblewski, T.; Michelmore, R.W.; Meyers, B.C. The role of TIR-NBS and TIR-X proteins in plant basal defense responses. Plant Physiol. 2013, 162, 1459–1472. [Google Scholar] [CrossRef] [Green Version]
  49. Neupane, S.; Ma, Q.; Mathew, F.M.; Varenhorst, A.J.; Andersen, E.J.; Nepal, M.P. Evolutionary divergence of TNL disease-resistant proteins in soybean (Glycine max) and common bean (Phaseolus vulgaris). Biochem. Genet. 2018, 56, 397–422. [Google Scholar] [CrossRef]
  50. Karasov, T.L.; Chae, E.; Herman, J.J.; Bergelson, J. Mechanisms to mitigate the trade-off between growth and defense. Plant Cell 2017, 29, 666–680. [Google Scholar] [CrossRef] [Green Version]
  51. Garreton, V.; Carpinelli, J.; Jordana, X.; Holuigue, L. The as-1 promoter element is an oxidative stress-responsive element and salicylic acid activates it via oxidative species. Plant Physiol. 2002, 130, 1516–1526. [Google Scholar] [CrossRef] [Green Version]
  52. Grüner, R.; Strompen, G.; Pfitzner, A.J.P.; Pfitzner, U.M. Salicylic acid and the hypersensitive response initiate distinct signal transduction pathways in tobacco that converge on the as-1-like element of the PR-1a promoter. JBIC J. Biol. Inorg. Chem. 2003, 270, 4876–4886. [Google Scholar]
  53. Zhou, J.-M.; Zhang, Y. Plant Immunity: Danger Perception and Signaling. Cell 2020, 181, 978–989. [Google Scholar] [CrossRef] [PubMed]
  54. Yang, T.; Xu, Z.P.; Lv, R.; Zhu, L.S.; Peng, Q.D.; Qiu, L.; Xi, D.H. N gene enhances resistance to Chili veinal mottle virus and hypersensitivity to salt stress in tobacco. J. Plant Physiol. 2018, 230, 92–100. [Google Scholar] [CrossRef]
  55. Gutterson, N.; Reuber, T.L. Regulation of disease resistance pathways by AP2/ERF transcription factors. Curr. Opin. Plant Biol. 2004, 7, 465–471. [Google Scholar] [CrossRef]
  56. Huang, P.-Y.; Catinot, J.; Zimmerli, L. Ethylene response factors in Arabidopsis immunity. J. Exp. Bot. 2015, 67, 1231–1241. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Wu, J.; Deng, Y.; Hu, J.; Jin, C.; Zhu, X.; Li, D. Genome-wide analyses of direct target genes of an ERF11 transcription factor involved in plant defense against bacterial pathogens. Biochem. Biophys. Res. Commun. 2020, 532, 76–81. [Google Scholar] [CrossRef] [PubMed]
  58. Hormuzdi, S.G.; Penttinen, R.; Jaenisch, R.; Bornstein, P. A Gene-targeting approach identifies a function for the first intron in expression of the alpha 1 (I) collagen gene. Mol. Cell. Biol. 1998, 18, 3368–3375. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  59. Furger, A.; O’Sullivan, J.M.; Binnie, A.; Lee, B.A.; Proudfoot, N.J. Promoter proximal splice sites enhance transcription. Genes Dev. 2002, 16, 2792–2799. [Google Scholar] [CrossRef] [Green Version]
  60. Okkema, P.G.; Harrison, S.W.; Plunger, V.; Aryana, A.; Fire, A. Sequence requirements for myosin gene expression and regulation in Caenorhabditis elegans. Genetics 1993, 135, 385–404. [Google Scholar] [CrossRef]
  61. Moabbi, A.M.; Agarwal, N.; El Kaderi, B.; Ansari, A. Role for gene looping in intron-mediated enhancement of transcription. Proc. Natl. Acad. Sci. USA 2012, 109, 8505–8510. [Google Scholar] [CrossRef] [Green Version]
  62. Morello, L.; Giani, S.; Troina, F.; Breviario, D. Testing the IMEter on rice introns and other aspects of intron-mediated enhancement of gene expression. J. Exp. Bot. 2010, 62, 533–544. [Google Scholar] [CrossRef] [PubMed]
  63. Parra, G.; Bradnam, K.; Rose, A.B.; Korf, I. Comparative and functional analysis of intron-mediated enhancement signals reveals conserved features among plants. Nucleic Acids Res. 2011, 39, 5328–5337. [Google Scholar] [CrossRef] [Green Version]
  64. Liao, L.; Ning, G.; Liu, C.; Zhang, W.; Bao, M. The intron from the 5′-UTR of the FBP11 gene in petunia displays promoter- and enhancer-like functions. Sci. Hortic. 2013, 154, 96–101. [Google Scholar] [CrossRef]
  65. Nott, A.; Le Hir, H.; Moore, M.J. Splicing enhances translation in mammalian cells: An additional function of the exon junction complex. Genes Dev. 2004, 18, 210–222. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Samadder, P.; Sivamani, E.; Lu, J.; Li, X.; Qu, R. Transcriptional and post-transcriptional enhancement of gene ex-pression by the 5′ UTR intron of rice rubi3 gene in transgenic rice cells. Mol. Genet. Genom. 2008, 279, 429–439. [Google Scholar] [CrossRef] [PubMed]
Viruses 13 02032 g001 550
Figure 1. Antiviral activity of GmRUN1. (A) Gene architecture of GmRUN1. Block and line indicate exon and intron, respectively. (B) Evolutionary analyses of GmRUN1. The evolutionary tree was built using the neighbor-joining method conducted in MEGA7. All positions with less than 50% site coverage were eliminated. (C) Transient expression assay for GmRUN1 antiviral activity for TMV. N. benthamiana leaves were infiltrated with Agrobacterium tumefaciens GV3101 inocula (OD600 = 1.0) carrying different recombinant vectors and co-infected with TMV-GFP. GFP was visualized under hand-held UV lamp (Wavelength = 365 nm) at 5 dpi (days post infiltration). N: tobacco N protein. SRC7: SMV resistance cluster 7. EV: empty vector. GmRUN1 cDNA-1: GmRUN1 cDNA expressed from pCB301 vector. GmRUN1 cDNA-2: GmRUN1 cDNA expressed from pBI211 vector. (D) Fluorescence quantification of GmRUN1 transient expression assay. TMV-GFP intensity was analyzed by Gel-Pro analyzer software and normalized against positive control (N). (E) Transient expression assay for GmRUN1 antiviral activity for SMV.
Figure 1. Antiviral activity of GmRUN1. (A) Gene architecture of GmRUN1. Block and line indicate exon and intron, respectively. (B) Evolutionary analyses of GmRUN1. The evolutionary tree was built using the neighbor-joining method conducted in MEGA7. All positions with less than 50% site coverage were eliminated. (C) Transient expression assay for GmRUN1 antiviral activity for TMV. N. benthamiana leaves were infiltrated with Agrobacterium tumefaciens GV3101 inocula (OD600 = 1.0) carrying different recombinant vectors and co-infected with TMV-GFP. GFP was visualized under hand-held UV lamp (Wavelength = 365 nm) at 5 dpi (days post infiltration). N: tobacco N protein. SRC7: SMV resistance cluster 7. EV: empty vector. GmRUN1 cDNA-1: GmRUN1 cDNA expressed from pCB301 vector. GmRUN1 cDNA-2: GmRUN1 cDNA expressed from pBI211 vector. (D) Fluorescence quantification of GmRUN1 transient expression assay. TMV-GFP intensity was analyzed by Gel-Pro analyzer software and normalized against positive control (N). (E) Transient expression assay for GmRUN1 antiviral activity for SMV.
Viruses 13 02032 g001
Viruses 13 02032 g002 550
Figure 2. Truncation analysis of GmRUN1. (A) Domain architecture of GmRUN1. Yellow, green, and blue boxes indicate TIR, NBS, and LRR domains, respectively. TN: TIR-NBS, NL: NBS-LRR, TNL: TIR-NBS-LRR. (B,D) Transient expression assay of different domains for antiviral activity. N. benthamiana leaves were infiltrated with Agrobacterium GV3101 inocula (OD600 = 1.0) carrying different recombinant vectors and co-infected with TMV-GFP (B) or SMV (D). GFP was visualized under hand-held UV lamp (Wavelength = 365 nm) at 5 dpi (days post infiltration). (C) Fluorescence quantification of truncated GmRUN1 transient expression assay. TMV-GFP intensity was analyzed by Gel-Pro analyzer software and normalized against positive control (GmRUN1 gDNA).
Figure 2. Truncation analysis of GmRUN1. (A) Domain architecture of GmRUN1. Yellow, green, and blue boxes indicate TIR, NBS, and LRR domains, respectively. TN: TIR-NBS, NL: NBS-LRR, TNL: TIR-NBS-LRR. (B,D) Transient expression assay of different domains for antiviral activity. N. benthamiana leaves were infiltrated with Agrobacterium GV3101 inocula (OD600 = 1.0) carrying different recombinant vectors and co-infected with TMV-GFP (B) or SMV (D). GFP was visualized under hand-held UV lamp (Wavelength = 365 nm) at 5 dpi (days post infiltration). (C) Fluorescence quantification of truncated GmRUN1 transient expression assay. TMV-GFP intensity was analyzed by Gel-Pro analyzer software and normalized against positive control (GmRUN1 gDNA).
Viruses 13 02032 g002
Viruses 13 02032 g003 550
Figure 3. Influence of GmRUN1 introns insertions on LUC reporter gene expression. (A) Schematic diagram of LUC recombinant vector construction (top panel) and transient expression (bottom panel). (B) Luciferase reporter assay of GmRUN1 introns. LUC recombinant vectors with or without different GmRUN1 introns were transiently expressed in N. benthamiana leaves. Images were taken using a Berthold camera 24 and 30 h after infiltration. (C) Semiquantitative PCR for transcription assay of LUC inserted different GmRUN1 introns. Lines: (1), LUC cDNA, 1653 bp; (2), LUC-intron2 cDNA, 1896 bp; (3), LUC-intron3 cDNA, 1787 bp; (4), LUC-intron4 cDNA, 1755 bp; (5–8), N. benthamiana ACTIN of line 1~4; and (9) Negative control. (D) qRT-PCR assays of LUC mRNA level after inserted with GmRUN1 introns. Error bars show the SD between biological replicates performed (n = 3), and Tukey’s multiple comparisons test was performed between samples in different groups.
Figure 3. Influence of GmRUN1 introns insertions on LUC reporter gene expression. (A) Schematic diagram of LUC recombinant vector construction (top panel) and transient expression (bottom panel). (B) Luciferase reporter assay of GmRUN1 introns. LUC recombinant vectors with or without different GmRUN1 introns were transiently expressed in N. benthamiana leaves. Images were taken using a Berthold camera 24 and 30 h after infiltration. (C) Semiquantitative PCR for transcription assay of LUC inserted different GmRUN1 introns. Lines: (1), LUC cDNA, 1653 bp; (2), LUC-intron2 cDNA, 1896 bp; (3), LUC-intron3 cDNA, 1787 bp; (4), LUC-intron4 cDNA, 1755 bp; (5–8), N. benthamiana ACTIN of line 1~4; and (9) Negative control. (D) qRT-PCR assays of LUC mRNA level after inserted with GmRUN1 introns. Error bars show the SD between biological replicates performed (n = 3), and Tukey’s multiple comparisons test was performed between samples in different groups.
Viruses 13 02032 g003
Viruses 13 02032 g004 550
Figure 4. GmRUN1 intron2 enhances SRC7 expression. (A) Gene architecture of SRC7. Block and line indicate exon and intron, respectively. (B) qRT-PCR assays of SRC7 mRNA level after replaced by GmRUN1 intron2. Error bars show the SD between biological replicates performed (n = 3), and Tukey’s multiple comparisons test was performed between samples in different groups. (C) Transient expression assay for recombinant SRC7-intron2 antiviral activity. N. benthamiana leaves were infiltrated with Agrobacterium tumefaciens GV3101 inocula carrying different recombinant vectors (OD600 = 0.01) and co-infected with TMV-GFP. N: tobacco N protein. SRC7: SMV resistance cluster 7. SRC7-intron2: SRC7 intron replaced by GmRUN1 intron2. EV: empty vector.
Figure 4. GmRUN1 intron2 enhances SRC7 expression. (A) Gene architecture of SRC7. Block and line indicate exon and intron, respectively. (B) qRT-PCR assays of SRC7 mRNA level after replaced by GmRUN1 intron2. Error bars show the SD between biological replicates performed (n = 3), and Tukey’s multiple comparisons test was performed between samples in different groups. (C) Transient expression assay for recombinant SRC7-intron2 antiviral activity. N. benthamiana leaves were infiltrated with Agrobacterium tumefaciens GV3101 inocula carrying different recombinant vectors (OD600 = 0.01) and co-infected with TMV-GFP. N: tobacco N protein. SRC7: SMV resistance cluster 7. SRC7-intron2: SRC7 intron replaced by GmRUN1 intron2. EV: empty vector.
Viruses 13 02032 g004
Viruses 13 02032 g005 550
Figure 5. Antiviral activity of GmRUN1 truncations in different introns. (A) Schematic diagram of truncated vector construction with different GmRUN1 introns. Block and line indicate exon and intron, respectively. (B) Transient expression assay of truncated GmRUN1 introns for antiviral activity. N. benthamiana leaves were infiltrated with Agrobacterium tumefaciens GV3101 inocula (OD600 = 0.5) carrying different recombinant vectors and co-infected with TMV-GFP. GFP was visualized under hand-held UV lamp (Wavelength = 365 nm) at 5 dpi (days post infiltration).
Figure 5. Antiviral activity of GmRUN1 truncations in different introns. (A) Schematic diagram of truncated vector construction with different GmRUN1 introns. Block and line indicate exon and intron, respectively. (B) Transient expression assay of truncated GmRUN1 introns for antiviral activity. N. benthamiana leaves were infiltrated with Agrobacterium tumefaciens GV3101 inocula (OD600 = 0.5) carrying different recombinant vectors and co-infected with TMV-GFP. GFP was visualized under hand-held UV lamp (Wavelength = 365 nm) at 5 dpi (days post infiltration).
Viruses 13 02032 g005
Viruses 13 02032 g006 550
Figure 6. Transcriptional regulation of GmRUN1 promoter. (A) Schematic diagram of truncated vector construction with different region of GmRUN1 promoter. Colored dashed lines and boxes indicate cis-acting regulatory element predicted by PlantCare database. (B) GUS activity assay of different region of GmRUN1 promoter. GUS activity was detected at 3 days post SA or SMV induction. EV: pBI121 empty vector. Pro1~4: Pro2415:GUS, Pro2592:GUS, Pro2237:GUS, and Pro2414:GUS. EVΔGUS: pBI121 empty vector removed GUS gene.
Figure 6. Transcriptional regulation of GmRUN1 promoter. (A) Schematic diagram of truncated vector construction with different region of GmRUN1 promoter. Colored dashed lines and boxes indicate cis-acting regulatory element predicted by PlantCare database. (B) GUS activity assay of different region of GmRUN1 promoter. GUS activity was detected at 3 days post SA or SMV induction. EV: pBI121 empty vector. Pro1~4: Pro2415:GUS, Pro2592:GUS, Pro2237:GUS, and Pro2414:GUS. EVΔGUS: pBI121 empty vector removed GUS gene.
Viruses 13 02032 g006
Viruses 13 02032 g007 550
Figure 7. GmDREB3 positively regulates GmRUN1. (A) Venn diagrams of predicted TFs transcription level consistent with GmRUN1 at 1/5/10 dpi (left panel). Red number indicates GmDREB3 which shows consistent expression pattern with GmRUN1 in graph at 1/5/10 dpi (right panel). (B) GUS activity assay of different region of GmRUN1 promoter coexpressed with GmDREB3. GUS activity was detected at 3 dpi. EV: pBI121 empty vector. Pro1~4: Pro2415:GUS, Pro2592:GUS, Pro2237:GUS, and Pro2414:GUS. EVΔGUS: pBI121 empty vector removed GUS gene.
Figure 7. GmDREB3 positively regulates GmRUN1. (A) Venn diagrams of predicted TFs transcription level consistent with GmRUN1 at 1/5/10 dpi (left panel). Red number indicates GmDREB3 which shows consistent expression pattern with GmRUN1 in graph at 1/5/10 dpi (right panel). (B) GUS activity assay of different region of GmRUN1 promoter coexpressed with GmDREB3. GUS activity was detected at 3 dpi. EV: pBI121 empty vector. Pro1~4: Pro2415:GUS, Pro2592:GUS, Pro2237:GUS, and Pro2414:GUS. EVΔGUS: pBI121 empty vector removed GUS gene.
Viruses 13 02032 g007
Table
Table 1. Phenotypic statistics of TMV appearance upon transient expression of different genes.
Table 1. Phenotypic statistics of TMV appearance upon transient expression of different genes.
Gene NameFull ResistanceNo ResistancePartial Resistance
pBI121-N30/30 (100%)0/30 (0%)0/30 (0%)
pBI121-SRC730/30 (100%)0/30 (0%)0/30 (0%)
pBI121-GmRUN1 gDNA24/30 (80%)0/30 (0%)6/30 (20%)
pCB301-GmRUN1 cDNA0/30 (0%)13/30 (43%)17/30 (57%)
pBI121-GmRUN1 cDNA0/30 (0%)24/30 (80%)6/30 (20%)
pBI1210/30 (0%)30/30 (100%)0/30 (0%)
Table
Table 2. Phenotypic statistics of SMV upon transient expression of different genes.
Table 2. Phenotypic statistics of SMV upon transient expression of different genes.
Gene NameHR/Total Leaves
pBI121-GmRUN1 gDNA14/30 (47%)
pCB301-GmRUN1 cDNA9/30 (30%)
pBI121-GmRUN1 cDNA7/30 (23%)
pCB3010/30 (0%)
Table
Table 3. Phenotypic statistics of TMV appearance upon transient expression of different genes.
Table 3. Phenotypic statistics of TMV appearance upon transient expression of different genes.
Gene NameFull ResistanceNo ResistancePartial Resistance
pBI121-GmRUN1 gDNA17/25 (68%)0/25 (0%)8/25 (32%)
pCB301-GmRUN1TIR0/25 (0%)25/25 (100%)0/25 (0%)
pCB301-GmRUN1NBS0/25 (0%)25/25 (100%)0/25 (0%)
pCB301- GmRUN1LRR0/25 (0%)25/25 (100%)0/25 (0%)
pCB301-GmRUN1TN0/25 (0%)25/25 (100%)0/25 (0%)
pCB301-GmRUN1NL0/25 (0%)25/25 (100%)0/25 (0%)
pCB301-GmRUN1TNL0/25 (0%)6/25 (24%)19/25 (76%)
pCB3010/25 (0%)25/25 (100%)0/25 (0%)
Table
Table 4. Phenotypic statistics of SMV upon transient expression of different genes.
Table 4. Phenotypic statistics of SMV upon transient expression of different genes.
Gene NameHR/Total Leaves
pBI121-GmRUN1 gDNA10/25 (40%)
pCB301-GmRUN1TIR1/25 (4%)
pCB301-GmRUN1NBS0/25 (0%)
pCB301-GmRUN1LRR0/25 (0%)
pCB301-GmRUN1TN1/25 (4%)
pCB301-GmRUN1NL0/25 (0%)
pCB301-GmRUN1TNL6/25 (24%)
pCB3010/25 (0%)
Table
Table 5. Phenotypic statistics of TMV appearance upon transient expression of different genes.
Table 5. Phenotypic statistics of TMV appearance upon transient expression of different genes.
Gene NameFull ResistanceNo ResistancePartial Resistance
pCambia1300-SRC7-intron215/20 (75%)0/20 (0%)5/20 (25%)
pCambia1300-SRC77/20 (35%)0/20 (0%)13/20 (65%)
pCambia1300-N8/20 (40%)0/20 (0%)12/20 (60%)
pCambia13000/20 (0%)20/20 (100%)0/20 (0%)
Table
Table 6. Phenotypic statistics of TMV appearance upon transient expression of different genes.
Table 6. Phenotypic statistics of TMV appearance upon transient expression of different genes.
Gene NameFull ResistanceNo ResistancePartial Resistance
pCB301-gDNA9/35 (26%)0/35 (0%)26/35 (74%)
pCB301-GmRUN1-ExIn-115/35 (43%)0/35 (0%)20/35 (57%)
pCB301-GmRUN1-ExIn-52/35 (6%)0/35 (0%)33/35 (94%)
pCB301-GmRUN1-ExIn-112/35 (6%)5/35 (14%)28/35 (80%)
pCB301-cDNA0/35 (0%)17/35 (49%)18/35 (51%)
pCB3010/35 (0%)35/35 (100%)0/35 (0%)
Table
Table 7. Promoter elements and their functions of GmRUN1.
Table 7. Promoter elements and their functions of GmRUN1.
Promoter ElementSequenceFunction
TATA-boxTATAACore promoter element around –30 of transcription start
CAAT-boxCCAATCommon cis-acting element in promoter and enhancer regions
CGTCA-motifCGTCACis-acting regulatory element involved in the MeJA-responsiveness
TGACG-motifTGACGCis-acting regulatory element involved in the MeJA-responsiveness
as-1TGACGCis-acting element related to salicylic acid induction
TCA-elementTCAGAAGAGGCis-acting element involved in salicylic acid responsiveness
TCATCATCTTCATUnknown functional element
ABREACGTGCis-acting element involved in the abscisic acid responsiveness
EREATTTTAAACis-acting element involved in ethylene response
AREAAACCACis-acting regulatory element essential for the anaerobic induction
LTRCCGAAACis-acting element involved in low-temperature responsiveness
MBSCAACTGMYB binding site involved in drought inducibility
W-boxTTGACCCis-acting element involved in disease resistance
CAT-boxGCCACTCis-acting regulatory element related to meristem expression
AT1-motifAATTATTTTTTATTPart of a light-responsive module
GT1-motifGGTTAATLight-responsive element
TCT-motifTCTTAC Part of a light-responsive element
G-BoxCACGTGCis-acting regulatory element involved in light responsiveness
MREAACCTAAMYB binding site involved in light responsiveness
Box 4ATTAATPart of a conserved DNA module involved in light responsiveness
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Diao, P.; Sun, H.; Bao, Z.; Li, W.; Niu, N.; Li, W.; Wuriyanghan, H. Expression of an Antiviral Gene GmRUN1 from Soybean Is Regulated via Intron-Mediated Enhancement (IME). Viruses 2021, 13, 2032. https://0-doi-org.brum.beds.ac.uk/10.3390/v13102032

AMA Style

Diao P, Sun H, Bao Z, Li W, Niu N, Li W, Wuriyanghan H. Expression of an Antiviral Gene GmRUN1 from Soybean Is Regulated via Intron-Mediated Enhancement (IME). Viruses. 2021; 13(10):2032. https://0-doi-org.brum.beds.ac.uk/10.3390/v13102032

Chicago/Turabian Style

Diao, Pengfei, Hongyu Sun, Zhuo Bao, Wenxia Li, Niu Niu, Weimin Li, and Hada Wuriyanghan. 2021. "Expression of an Antiviral Gene GmRUN1 from Soybean Is Regulated via Intron-Mediated Enhancement (IME)" Viruses 13, no. 10: 2032. https://0-doi-org.brum.beds.ac.uk/10.3390/v13102032

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop