Next Article in Journal
Three-Source Partitioning of Methane Emissions from Paddy Soil: Linkage to Methanogenic Community Structure
Next Article in Special Issue
Analysis of Rice Proteins with DLN Repressor Motif/S
Previous Article in Journal
Dulaglutide Modulates the Development of Tissue-Infiltrating Th1/Th17 Cells and the Pathogenicity of Encephalitogenic Th1 Cells in the Central Nervous System
Previous Article in Special Issue
Transcriptome Analysis Reveals the Accumulation Mechanism of Anthocyanins in Buckwheat (Fagopyrum esculentum Moench) Cotyledons and Flowers
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Induction of Transgene Suppression in Plants via External Application of Synthetic dsRNA

by
Alexandra S. Dubrovina
1,*,
Olga A. Aleynova
1,
Alexander V. Kalachev
2,
Andrey R. Suprun
1,3,
Zlata V. Ogneva
1 and
Konstantin V. Kiselev
1,3
1
Laboratory of Biotechnology, Federal Scientific Center of the East Asia Terrestrial Biodiversity, Far Eastern Branch of the Russian Academy of Sciences, Vladivostok 690022, Russia
2
Laboratory of Embryology, National Scientific Center of Marine Biology, Far Eastern Branch of the Russian Academy of Sciences, Vladivostok 690041, Russia
3
Far Eastern Federal University, The School of Natural Sciences, Vladivostok 690090, Russia
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2019, 20(7), 1585; https://doi.org/10.3390/ijms20071585
Submission received: 1 February 2019 / Revised: 19 March 2019 / Accepted: 26 March 2019 / Published: 29 March 2019
(This article belongs to the Special Issue Transcriptional and Post-transcriptional Gene Regulation in Plants)

Abstract

:
Recent investigations show that exogenously applied small interfering RNAs (siRNA) and long double-stranded RNA (dsRNA) precursors can be taken up and translocated in plants to induce RNA interference (RNAi) in the plant or in its fungal pathogen. The question of whether genes in the plant genome can undergo suppression as a result of exogenous RNA application on plant surface is almost unexplored. This study analyzed whether it is possible to influence transcript levels of transgenes, as more prone sequences to silencing, in Arabidopsis genome by direct exogenous application of target long dsRNAs. The data revealed that in vitro synthesized dsRNAs designed to target the gene coding regions of enhanced green fluorescent protein (EGFP) or neomycin phosphotransferase II (NPTII) suppressed their transcript levels in Arabidopsis. The fact that, simple exogenous application of polynucleotides can affect mRNA levels of plant transgenes, opens new opportunities for the development of new scientific techniques and crop improvement strategies.

1. Introduction

The increasing human population and discussions along the safety of transgenic plants promote the development of new strategies to regulate plant properties without genomic manipulations. Numerous investigations show that it is possible to switch off or decrease expression of particular genes for the regulation of plant stress tolerance, growth, and other processes via induction of RNA interference (RNAi or gene silencing process) [1]. It is known that RNAi serves for the regulation of various processes important for plants, such as growth and development, stress adaptation, or synthesis of biologically active compounds [2,3]. In the course of RNAi, double-stranded RNAs (dsRNAs) are processed by a ribonuclease into small interfering RNA (siRNA) or microRNA (miRNA). These small RNAs are incorporated in the RNA-induced silencing complex that provides cleavage, destabilization, or hindering translation of any homologous mRNAs [4]. As a result, the mRNAs are degraded, or their efficient translation is prevented. The application of RNAi-based approaches, for gene regulation in plants, requires plant genomic modifications or application of the weakened plant viruses.
There are new intriguing messages that synthetic dsRNAs or siRNAs, exogenously applied to the plant surface (direct spray application, mechanical inoculation, loading on clay nanosheets, or using materials promoting RNA adhesion), are capable of entering and spreading into plant vascular system and plant cells [5,6,7,8,9] with the following silencing of the targeted essential genes of infected pathogens and the development of local and system resistance against fungal [5,9,10,11] or viral [6,7,12,13] infections. There is also evidence for processing the external dsRNAs into siRNAs, and induction of RNAi-mediated silencing of the targeted pathogen genes [5,6,7,12,14]. Since RNA molecules can be up-taken and translocated into plant tissues, active studies in this direction are extremely urgent, not only for confirmation and detailed investigation of these facts, but also for investigating the possibility of regulating the expression of plant endogenous genes via this approach. Intriguingly, the patent by Sammons et al. [15] exemplified the possibility of down-regulating transcript levels of plant endogenous genes encoding for herbicide resistance using direct exogenous applications of dsRNA, siRNA, ssRNA, or even DNA molecules. Numata et al. [16] indicated that spraying Arabidopsis with synthetic siRNAs in a complex with a protein carrier induced local loss of anthocyanin pigmentation, probably, due to the silencing of a chalcone synthase gene.
A number of studies provided that plant transgenes are more susceptible to silencing than endogenes [17,18,19,20]. Transgenes usually do not contain introns and 5’/3’-UTRs, which have been shown to contribute to the expression stability of plant endogenes and suppression of RNA silencing [18,19,21]. In addition, transgenes are usually under the control of strong promoters, ensuring a high level of expression and, thus, elevating the possible production of aberrant mRNAs. Aberrant mRNAs have been shown to be perceived by plant RNAi machinery and transcribe into secondary dsRNAs by RNA-directed RNA polymerase 6 (RDR6) leading to amplification of the transgene silencing [22]. The secondary dsRNAs are then converted into secondary siRNAs that map upstream and downstream the primary siRNAs [18]. Therefore, it would be interesting to analyze the influence of external RNAs on transgene activity first. The available studies reported that they failed to regulate the fluorescence of GFP and YFP in transgenic plants by direct siRNA application (e.g. wiping, spraying, injection) without using additional techniques, i.e., a protein carrier or high-pressure spraying [16,23]. Dalakouras et al. [23] reported that spraying GFP-transgenic Nicotiana benthamiana with synthetic GFP-siRNA solutions, under high pressure, using a conventional compressor and an air brush pistol, led to silencing of the GFP fluorescence. The effect was not observed when high pressure was not applied. Numata et al. [16] found that infiltration of in vitro synthesized YFP-siRNA in a complex with carrier peptide into leaves of A. thaliana and poplar resulted in lowering the YFP protein levels and fluorescence observed by confocal laser scanning microscopy. In our study, we aimed to investigate whether it is possible to influence the expression of plant transgenes by simple application of long dsRNA solutions without using additional promoting agents or techniques.

2. Results

To create transgenic Arabidopsis, we introduced pZP-RCS2 [24], bearing the neomycin phosphotransferase II (NPTII), and the enhanced green fluorescent protein (EGFP) genes under the control of the double CaMV 35S promoters (Figure 1a) into Agrobacterium tumefaciens and transformed A. thaliana by floral dip [25]. Three representative independent T3 homozygous lines with high mRNA levels of EGFP and NPTII were chosen for detailed analyses. The quantification of the transgene transcript levels in A. thaliana, performed by quantitative real-time RT-PCR (qRT-PCR), showed that the transgenic lines displayed high levels of NPTII and EGFP expression (Figure 1b,c). Then, we used PCR and in vitro transcription protocol to produce dsRNA molecules of the EGFP and the NPTII genes. The complete EGFP (720 bp) and a large fragment of the NPTII (599 bp out of 798 bp) genes were amplified. The obtained PCR products, containing T7 promoters at both ends, were used as templates for in vitro transcription. For external applications, the synthesized EGFP-dsRNA and NPTII-dsRNA were diluted in water to final concentrations of 0.1, 0.35, and 1 µg/µL. The dsRNAs (100 µL of each dsRNA per individual plant, i.e. 10, 35, and 100 µg) were applied on the surface of four-week-old Arabidopsis by spreading with sterile individual soft brushes. The analysis of concentration effect was performed on the L1 line of A. thaliana and analyzed 7 days post-treatment, indicating that 0.35 µg/µL, resulted in the highest silencing efficiency (Figure 1d,e). Therefore, this concentration was chosen for the following experiments where we studied whether simple exogenous applications of the NPTII-dsRNA and EGFP-dsRNA on plant surface of the three transgenic lines L1, L2, and L3 could lead to changes in transcript levels of the transgenes in comparison with the water-treated control plants 1, 7, and 14 days post-treatment.
qRT-PCR showed that the NPTII transcript levels increased in plants treated with sterile water during cultivation compared with the plants before treatments (Figure 2a). However, the elevation in NPTII transcript levels was lower or was not observed in plants treated with external NPTII-dsRNA (Figure 2a). As for EGFP expression, plants that received EGFP-dsRNA exhibited down-regulation or repression of EGFP expression after the application of the exogenous EGFP-dsRNA, while water-treated plants exhibited either, an increase, or a weak decrease, in EGFP expression (Figure 2b).
We also analyzed NPTII and EGFP expression using primers (Figure 1a; Table S1), designed to align inside the NPTII and EGFP transgene fragments, that have been used for synthesis of the corresponding dsRNAs (Figure 3a,b). Plants that have been treated with dsRNAs, exhibited dramatically high NPTII or EGFP levels, which gradually decreased with time of cultivation (Figure 3a,b). The data indicate that the synthetic NPTII and EGFP dsRNAs have been purified with the total RNA and transcribed into cDNA during reverse transcription. Thus, the NPTII- and EGFP-dsRNAs were stable for at least 7 days of cultivation and were not completely removed from the treated leaves during this time (Figure 3a,b). Figure 3 shows suppression of NPTII in plants that received EGFP-dsRNA, and EGFP in plants that received NPTII-dsRNA, respectively (Figure 3a,b).
Confocal microscopy was used to analyze EGFP fluorescence in L1 of A. thaliana and showed that EGFP fluorescence sharply decreased 1 and 7 days post-treatment with the EGFP dsRNA in the L1 line (Figure 4a). The data are in agreement with the EGFP gene expression analysis, indicating a large decrease in EGFP transcript levels 1 and 7 days post-treatment (Figure 2b). Western blotting showed a sharp decline in EGFP protein level 1 and 7 days post-treatment in L1 (Figure 4b). The results indicate that transgene suppression occurred not only at mRNA but also protein level.
To determine whether EGFP-siRNAs accumulate in the A. thaliana L1 plants, we performed stem-loop RT-PCR as described [26] for small RNA fractions purified from the dsRNA-treated, and water-treated, A. thaliana leaves. This method was originally applied for miRNA detection and then successfully used for siRNA detection in plant tissues [7,8,12]. We used as the target siRNA one of the EGFP-encoding siRNAs (NCBI GB accession: Pr008640320, Pr120179) with a typical structure and with a documented EGFP silencing efficiency [27,28]. Stem-loop RT-PCR followed by cloning and sequencing of the RT-PCR products revealed the presence of the target EGFP-siRNA in the leaves treated with the EGFP-dsRNA, while the siRNA was not detected in the plants before treatment and in the control water-treated plant (Figure 5). The target EGFP-siRNA was also detected in the leaves treated with NPTII-dsRNA, though at a much lower level than after treatment with EGFP-dsRNA (Figure 5).
In case the EGFP-siRNAs were generated by transitive silencing along the transgene sequence after treatment with the NPTII-dsRNA, the read-through transcripts, that include both NPTII and EGFP sequences should be generated in the L1 transgenic plants. Therefore, we analyzed whether the readthrough transcripts, including both the NPTII and EGFP sequences, are accumulated (Figure 6a,b; Figure S1; Table S1). For RT-PCR, we used the cDNAs obtained from the transgenic A. thaliana L1 plants, treated with the NPTII- and EGFP-dsRNAs. The RT-PCR analysis (Figure 6a,b) revealed that there was no amplification in the cDNA probes of the wild-type A. thaliana plants, while a clear amplification product of 2–2.5 kb was detected in the cDNA probes from the transgenic plants and for the positive control (pZP-RCS2 plasmid). DNA sequencing revealed that the amplification product of 2.276 kb contained coding sequences of both transgenes and their regulatory elements (Figure S1). The RT-PCR analysis also indicated that the abundance of the transgene read-through transcripts was lowered after the dsRNA treatments but not after control water treatments (Figure 6a,b).
When analyzing transgene transcription after dsRNA treatment, we noted that the NPTII transcript level was lowered in A. thaliana treated with the EGFP-dsRNA (Figure 2a). Similarly, the EGFP transcript level was down-regulated in the plants treated with NPTII-dsRNA (Figure 2b). In case transitivity and silencing amplification occur, the reduction of NPTII mRNA levels by application of EGFP-dsRNA, and vice versa, can be explained by the secondary siRNA formation. It is also possible that the application of the NPTII- or EGFP-dsRNAs to individual plants could affect expression of the transgenes due to induction and spreading of DNA methylation at the T-DNA region bearing both transgenes. To examine this hypothesis, we analyzed cytosine DNA methylation levels before and after treatment of A. thaliana L1, with EGFP-and NPTII-dsRNAs, using bisulfite sequencing of the dsRNA-targeted transgene regions (Figure 6c,d; Table 1). The data revealed that the extent of DNA methylation considerably increased after the dsRNA treatments.
The analysis of NPTII and EGFP transgene cytosine methylation at different contexts (CG, CHG and CHH) showed that the dsRNA treatments increased methylation in all sequence contexts, while control water treatment did not essentially change methylation status (Table 1). The most pronounced effect for the dsRNA treatments was observed for CHG and CHH methylation elevating it by 1.2-1.3 and 1.3-1.5-fold, respectively. Methylation status at CG positions was considerably affected only for NPTII after NPTII-dsRNA treatment.

3. Discussion

In this study, the transgenic adult plants of A. thaliana, overexpressing the NPTII and EGFP, have been treated with long dsRNAs targeting these transgenes. Taken together, the results show that the exogenous application of dsRNA spread with soft brushes suppressed NPTII and EGFP transcript levels in A. thaliana. The current data also indicated that applying 35 µg of the transgene-encoding dsRNA per four-week-old plant is the most efficient amount, due to the combination of effectiveness and lower cost of synthesis. Thus, simple foliar applications of dsRNAs complimentary to commonly used transgenes lead to their repression. The detected dsRNA-induced down-regulation of transgene transcript levels suggests that the exogenously applied dsRNAs were uptaken by plants, and initiated RNAi-mediated silencing of the target mRNA transcripts. In our study, the detection of EGFP-derived siRNA in the dsRNA-treated plants, indicated that the EGFP-encoding dsRNAs were processed into siRNAs, and the observed transgene silencing effect was based on RNAi-mediated silencing. The detection of the EGFP-derived siRNA after NPTII-dsRNA treatment suggests that NPTII silencing was accompanied by the phenomenon of transitivity [18,19]. Transgenes have been shown to be more prone for transitivity, which is the spreading of silencing beyond the initial target site. Transgenes are usually under the control of strong promoters, ensuring a high level of expression and, thus, elevating possible production of aberrant mRNAs, e.g., truncated and/or read-through transcripts, leading to the amplification of the transgene silencing [18,19,22]. In our study, we detected transgene read-through transcripts, including both the NPTII and EGFP coding sequences and the regulatory elements, in all the dsRNA-treated and non-treated L1 transgenic A. thaliana plants. According to the data obtained, the silencing transitivity occurred before and after treatments in the transgenic A. thaliana plants and was enhanced after the dsRNA treatments. Surprisingly, we also noted that the NPTII mRNA level was suppressed in A. thaliana treated with the EGFP-dsRNA (Figure 2a). Similarly, the EGFP mRNA level was down-regulated in the plants treated with NPTII-dsRNA (Figure 2b). If transitivity occurs, the reduction of NPTII mRNA levels by the application of EGFP-dsRNA, and vice versa, can be explained by secondary siRNA formation. Further experiments, including small RNA sequencing, are needed to examine the whole diversity of the transgene-siRNAs generated after external dsRNA application.
In our study, we also found that the level of cytosine DNA methylation considerably increased after the dsRNA treatments at the dsRNA-targeted protein-coding regions of the transgenes. RNA-directed DNA methylation (RdDM), a process where siRNAs induce de novo DNA methylation of the target homologous DNA sequences, and the DNA sequences upstream/downstream of the target sites [19,29,30], is supposed to be responsible for methylation induced at the protein-coding regions of post-transcriptionally silenced genes [31,32]. Methylation at the gene coding regions was also presumed to contribute to PTGS reinforcement. In our study, DNA methylation at the coding sequences of the NPTII and EGFP transgenes homologous to the applied dsRNAs could reflect the influence of 24nt siRNAs, a RdDM “footprint”, arising as a byproduct of unprecise dicing. Therefore, it is possible that the application of the NPTII- or EGFP-encoding dsRNAs could promote transgene mRNA degradation and also affect transcription of the transgenes and/or heterochromatin formation, due to the induction and spreading of DNA methylation at the T-DNA region bearing both transgenes. Although only sporadic investigations report on the role of methylation in gene body regions for gene activity, there were at least some studies that documented the suppression of gene expression and heterochromatinization, due to coding region methylation [33,34].
Understanding the mechanisms for recognition and uptake of extracellular nucleic acids by plants has only begun in recent years [35]. Several studies documented active plant perception of extra-celluar RNAs originating from plant viruses and pathogenic bacteria and acting as genuine pathogen-associated molecular patterns (PAMPs) and microbe-associated molecular patterns (MAMPs) to induce the pattern-triggered plant immune signaling [36,37]. MAMPs and PAMPs are recognized by pattern recognition receptors (PRRs). Plant PRRs are either, surface-localized receptor kinases, or receptor-like proteins, containing various ligand-binding ectodomains [38]. The molecular mechanisms, and the specific PRRs responsible for extracellular nucleic acid recognition and uptake, are generally uncharacterized. Niehl et al. [37] showed that exogenously applied dsRNAs act in plants via a somatic embryogenesis receptor-like kinase 1 (SERK1), and not through antiviral dicer-like (DCL) proteins, and suggested SERK1 as a potential dsRNA receptor. It is likely that transgene-encoding dsRNAs applied in this study and dsRNAs used in other related studies [5,6,9,10,11,12,13,39] could also act via PRRs. Further research is needed to shed light on the mechanisms of extracellular dsRNA uptake, local and systemic dsRNA movements, and processing in plant tissues and cells. It is possible that exogenously applied dsRNA molecules could also affect transcript levels of plant endogenous genes related to various signaling pathways. Active studies in this direction could promote development of new alternatives to genetically modified crops and new appealing methodologies that can be used in plant functional studies.

4. Materials and Methods

4.1. Generation of Transgenic Arabidopsis

Plants (Arabidopsis thaliana ecotype Columbia L.) were grown in pots filled with commercially available rich soil in a controlled environmental chamber at 22 °C (Sanyo MLR-352, Panasonic, Osaka, Japan), kept on a 16/8 h day/night cycle at a light intensity of ~120 μmol m−2 s−1.
To generate the EGFP- and NPTII-overexpressing plants, we used the binary plasmid construction pZP-RCS2-EGFP-NPTII (Figure 1a) that was kindly provided by Professor Alexander Krichevsky (State University of New York, Stony Brook, NY, USA) [24]. The construction carried the NPTII and EGFP genes under the control of the double CaMV 35S promoters. A. thaliana was transformed by floral-dip as described [25]. The transgenic lines used in this study were T3 homozygous plants with single copy insertion.

4.2. dsRNA Synthesis and Application

dsRNA of EGFP and NPTII was synthesized using MEGAscript RNAi Kit (ThermoFisher Scientific, Waltham, MA, USA). The T7 promoter sequence was introduced into both the 5’ and 3’ ends of EGFP and NPTII by PCR. The primers are shown in Table S1 and Figure 1a. The obtained PCR products were used as templates for in vitro transcription and dsRNA synthesis following the manufacturer’s protocol. The resultant dsRNAs were analyzed by gel electrophoresis and were spectrophotometrically analyzed to estimate the amount, purity, and integrity of the dsRNAs. The dsRNAs were applied on the surface of four-week-old Arabidopsis by spreading with individual soft brushes (natural pony hair) sterilized by autoclaving. We did not observe wounding visually and using confocal microscopy. All the treatments were performed at 9:00–9:30 pm.

4.3. Total RNA Extraction, Reverse Transcription, and qRT-PCR

The leaves of A. thaliana for nucleic acid isolation were collected from the same four-week-old individual plant at all time points (before, 1, 7 and 14 days post-treatments) for each type of treatment. One typical adult leaf was used for the nucleic acid isolation each time. RNA isolation was carried out from the A. thaliana leaves before treatment, 1, 7, and 14 days post-treatments. The isolation of total RNA was performed using the cetyltrimethylammonium bromide (CTAB)-based protocol and complementary DNAs were synthesized as described [40]. The reverse transcription products were amplified by PCR and verified on the absence for DNA contamination using primers listed in Table S1. The qRT-PCRs were performed with EvaGreen Real-time PCR (Biotium, Hayward, CA, USA) as described in [25], using cDNAs of NPTII, EGFP, and two internal controls (GAPDH and UBQ) selected in previous studies as relevant reference genes for qRT-PCRs on Arabidopsis [41]. The expression was calculated by the 2−ΔΔCT method [42]. All GenBank accession numbers and primers are listed in Table S1.

4.4. Laser Scanning Microscopy

The leaves were mounted in distilled water in a Petri dish and were observed under a Zeiss LSM 780 confocal microscope (Far East Centre of Electron Microscopy, National Scientific Center of Marine Biology, Far Eastern Branch of the Russian Academy of Sciences, Vladivostok, Russia) equipped with a Plan-Apochromat 20×/0.8 objective. The excitation wavelength of an argon laser was set at 488 nm and EGFP emission was detected using a 520–525 nm band pass filter. To quantify EGFP fluorescence 10 Z-stacks were captured for each leaf (5 Z-stacks from the top and 5 Z-stacks from bottom side of the leaf). For fluorescence intensity quantification, Z-stacks were processed using ImageJ 1.51w image processing software [43,44]. Each Z-stack was imported into ImageJ and z-projected using sum projection. The projected images reached the threshold using the method from Otsu [45], and the area and mean fluorescence of the foreground (signal of EGFP) and the background along the entire image were measured. The corrected EGFP fluorescence was calculated according to the formula: Corrected EGFP fluorescence = mean fluorescence of the foreground − (area of the foreground × mean fluorescence of the background).

4.5. Protein Extraction and Western Blotting

Proteins were isolated from A. thaliana by grinding 30 mg of frozen leaf tissue in 300 µL of 0.1 M sodium phosphate buffer pH 8.0 [46] with 4% sodium dodecyl sulfate, followed by ultrasound treatment for 30 s at maximum amplitude, using Q55 Qsonica sonicator (Qsonica LLC, Newtown, CT, USA). After centrifugation (16,000× g for 20 min at +4 °C), the protein concentration was quantified using Qubit®Protein Assay Kit and Qubit 2.0 fluorometer (Invitrogen, Carlsbad, CA, USA). A quantity of 5 µg of the protein extracts were separated by SDS–PAGE in 12% polyacrylamide gels using a Helicon VE-10 (Helicon, Moscow, Russia) at 200 V and were electro-transferred to polyvinylidene difluoride (PVDF) membrane. After blocking for 1 h in PBS (1.7 mM KH2PO4, 5.2 mM Na2HPO4, 150 mM NaCl) with 2.5% nonfat dry milk at room temperature, membranes were incubated overnight at +5 °C with anti-GFP antibody (ab6556 Abcam, Cambridge, United Kingdom, USA) at 1:1000 dilution. Following three washings with PBS containing 0.1% Tween 20 and 2.5% nonfat dry milk at room temperature, the membrane was incubated with secondary antibody (ab6112 Abcam, USA) at 1:2000 dilution for 1 h. After three washings, the membranes were visualized by treatment with 40 mL water solution containing 10 mg 3,3′-diaminobenzidine, 4 mL EtOH, 2 mL of 1 M Tris-HCl (pH 7.6), and 25 µL H2O2. As a positive control, we used total proteins purified from E. coli Rosetta(DE3)pLysS transformed with pET28a-EGFP (Merck, Darmstadt, Germany). The EGFP from pZP-RCS2-EGFP-NPTII plasmid was PCR-amplified using primers (Table S1) containing BamH I and Hind III restriction sites, The PCR products were digested with the restriction enzymes (SibEnzyme, Novosibirsk, Russia), and the resulting DNA fragments were ligated to pET28a.

4.6. DNA Extraction and Bisulfite Sequencing

Sampling of the leaf material for DNA extraction was performed as described above. Extraction of genomic DNA and bisulfite sequencing were performed as described [47]. After DNA conversion, a 352-bp NPTII and a 349-bp EGFP fragments were amplified using primers listed in Table S1. The level of C to T transitions in the converted PCR products was greater than 95%. We sequenced 18 clones for each DNA region from the 2 biological replicates (9 clones per each individual plant).

4.7. Detection of EGFP-siRNA by Stem-Loop RT-PCR and DNA Sequencing

For the analysis, we retrieved EGFP-siRNAs from the NCBI siRNA database (Available online: https://www.ncbi.nlm.nih.gov/projects/genome/rnai/) and based our stem-loop primer construction on a EGFP-siRNA (Accession: Pr008640320, Pr120179) with typical structure deposited for the cloning vector pEGFP-N1 (accession U55762), which bears identical EGFP coding sequence to pZP-RCS2-EGFP-NPTII (used in our work). For isolation of low molecular weight (LMW) RNA fraction, we began isolation of total RNA using the CTAB-based protocol as described [40]. After incubation of 1 mL of the water phase, containing total RNA mixed, with 250 µL of 10 M lithium chloride at 4 °C overnight, we separated low-molecular weight (LMW) RNA by centrifugation at 16,000× g for 20 min, followed by the collection of the aqueous phase, containing LMW RNAs. Then, the LMW RNA was precipitated as described [48]. The RT reaction (20 µL) was performed with 0.2 µL of EGFP-siRNA stem-loop RT primer (100 µM) and ThermoScript reverse transcriptase (Invitrogen, Life Technologies) using the conditions of 16 °C for 30 min, 30 °C for 30 s, 42 °C for 30 s, 50 °C for 1 s for 60 cycles, and 85 °C for 5 min according to published pulsed RT protocol and primer design instructions [26]. The targeted EGFP-siRNA was PCR amplified as described [49] with the 25 µL of the PCR mixture containing 0.1 µL of EGFP-specific forward and universal reverse primers (100 µM) designed as described [26]. The amplification conditions consisted of one cycle at 95 °C for 2 min followed by 40 cycles of 95 °C for 10 s, 58 for 10 s, and 72 °C for 20 s. The primer sequences are listed in Table S1. The PCR products were isolated from agarose gels using a Cleanup Mini Kit (Eurogene, Moscow, Russia), subcloned, and sequenced as described [49].

Supplementary Materials

Supplementary materials can be found at https://www.mdpi.com/1422-0067/20/7/1585/s1. Table S1 Primers used for amplification of the NPTII and EGFP transgenes in RT-PCR, qRT-PCRs, and bisulfite sequencing. Fig. S1 DNA sequencing of the amplification product including coding sequences of both the EGFP and NPTII transgenes and their regulatory elements.

Author Contributions

A.S.D. and K.V.K. performed research design, experiments, data analysis, and paper preparation. O.A.A., A.R.S., and Z.V.O. were involved in plant management, RNA isolation, and performed qRT-PCRs. A.V.K. performed confocal microscopy.

Funding

This work was supported by the Russian Science Foundation (17-74-10083), except for laser scanning microscopy that was supported by the Far Eastern Center for Electron Microscopy (National Scientific Center of Marine Biology, Russian Academy of Sciences).

Acknowledgments

The authors are grateful to Kirill Shefer for the technical assistance in microscopic studies.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kamthan, A.; Chaudhuri, A.; Kamthan, M.; Datta, A. Small RNAs in plants: Recent development and application for crop improvement. Front. Plant Sci. 2015, 6, 208. [Google Scholar] [CrossRef]
  2. Borges, F.; Martienssen, R.A. The expanding world of small RNAs in plants. Nat. Rev. Mol. Cell Biol. 2015, 16, 727–741. [Google Scholar] [CrossRef] [Green Version]
  3. Singh, A.; Gautam, V.; Singh, S.; Sarkar Das, S.; Verma, S.; Mishra, V.; Mukherjee, S.; Sarkar, A.K. Plant small RNAs: advancement in the understanding of biogenesis and role in plant development. Planta 2018, 248, 545–558. [Google Scholar] [CrossRef] [PubMed]
  4. Wilson, R.C.; Doudna, J.A. Molecular mechanisms of RNA interference. Annu. Rev. Biophys. 2013, 42, 217–239. [Google Scholar] [CrossRef]
  5. Koch, A.; Biedenkopf, D.; Furch, A.; Weber, L.; Rossbach, O.; Abdellatef, E.; Linicus, L.; Johannsmeier, J.; Jelonek, L.; Goesmann, A.; et al. An RNAi-based control of Fusarium graminearum infections through spraying of long dsRNAs involves a plant passage and is controlled by the fungal silencing machinery. PLoS Pathog. 2016, 12, e1005901. [Google Scholar] [CrossRef] [PubMed]
  6. Mitter, N.; Worrall, E.A.; Robinson, K.E.; Li, P.; Jain, R.G.; Taochy, C.; Fletcher, S.J.; Carroll, B.J.; Lu, G.Q.; Xu, Z.P. Clay nanosheets for topical delivery of RNAi for sustained protection against plant viruses. Nat. Plants 2017, 3, 16207. [Google Scholar] [CrossRef] [PubMed]
  7. Konakalla, N.C.; Kaldis, A.; Berbati, M.; Masarapu, H.; Voloudakis, A.E. Exogenous application of double-stranded RNA molecules from TMV p126 and CP genes confers resistance against TMV in tobacco. Planta 2016, 244, 961–969. [Google Scholar] [CrossRef] [PubMed]
  8. Gogoi, A.; Sarmah, N.; Kaldis, A.; Perdikis, D.; Voloudakis, A. Plant insects and mites uptake double-stranded RNA upon its exogenous application on tomato leaves. Planta 2017, 246, 1233–1241. [Google Scholar] [CrossRef]
  9. Song, X.S.; Gu, K.X.; Duan, X.X.; Xiao, X.M.; Hou, Y.P.; Duan, Y.B.; Wang, J.X.; Zhou, M.G. A myosin5 dsRNA that reduces the fungicide resistance and pathogenicity of Fusarium asiaticum. Pest. Biochem. Physiol. 2018, 150, 1–9. [Google Scholar] [CrossRef]
  10. McLoughlin, A.G.; Wytinck, N.; Walker, P.L.; Girard, I.J.; Rashid, K.Y.; de Kievit, T.; Fernando, W.G.D.; Whyard, S.; Belmonte, M.F. Identification and application of exogenous dsRNA confers plant protection against Sclerotinia sclerotiorum and Botrytis cinerea. Sci. Rep. 2018, 9, 7320. [Google Scholar] [CrossRef] [PubMed]
  11. Gu, K.X.; Song, X.S.; Xiao, X.M.; Duan, X.X.; Wang, J.X.; Duan, Y.B.; Hou, Y.P.; Zhou, M.G. A β2-tubulin dsRNA derived from Fusarium asiaticum confers plant resistance to multiple phytopathogens and reduces fungicide resistance. Pestic. Biochem. Physiol. 2019, 153, 36–46. [Google Scholar] [CrossRef]
  12. Kaldis, A.; Berbati, M.; Melita, O.; Reppa, C.; Holeva, M.; Otten, P.; Voloudakis, A. Exogenously applied dsRNA molecules deriving from the Zucchini yellow mosaic virus (ZYMV) genome move systemically and protect cucurbits against ZYMV. Mol. Plant. Pathol. 2018, 19, 883–895. [Google Scholar] [CrossRef]
  13. Niehl, A.; Soininen, M.; Poranen, M.M.; Heinlein, M. Synthetic biology approach for plant protection using dsRNA. Plant Biotechnol. J. 2018, 16, 1679–1687. [Google Scholar] [CrossRef] [PubMed]
  14. Song, X.S.; Gu, K.X.; Duan, X.X.; Xiao, X.M.; Hou, Y.P.; Duan, Y.B.; Wang, J.X.; Yu, N.; Zhou, M.G. Secondary amplification of siRNA machinery limits the application of spray-induced gene silencing. Mol. Plant Pathol. 2018, 19, 2543–2560. [Google Scholar] [CrossRef]
  15. Sammons, R.; Ivashuta, S.; Liu, H.; Wang, D.; Feng, P.; Kouranov, A.; Andersen, S. Polynucleotide Molecules for Gene Regulation in Plants. U.S. Patent 2011/0296556 A1, 1 December 2011. [Google Scholar]
  16. Numata, K.; Ohtani, M.; Yoshizumi, T.; Demura, T.; Kodama, Y. Local gene silencing in plants via synthetic dsRNA and carrier peptide. Plant Biotechnol. J. 2014, 12, 1027–1034. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Vermeersch, L.; De Winne, N.; Nolf, J.; Bleys, A.; Kovařík, A.; Depicker, A. Transitive RNA silencing signals induce cytosine methylation of a transgenic but not an endogenous target. Plant J. 2013, 74, 867–879. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Vermeersch, L.; De Winne, N.; Depicker, A. Introns reduce transitivity proportionally to their length, suggesting that silencing spreads along the pre-mRNA. Plant J. 2010, 64, 392–401. [Google Scholar] [CrossRef] [Green Version]
  19. Dadami, E.; Dalakouras, A.; Zwiebel, M.; Krczal, G.; Wassenegger, M. An endogene-resembling transgene is resistant to DNA methylation and systemic silencing. RNA Biol. 2014, 11, 934–941. [Google Scholar] [CrossRef] [Green Version]
  20. Koscianska, E.; Kalantidis, K.; Wypijewski, K.; Sadowski, J.; Tabler, M. Analysis of RNA silencing in agroinfiltrated leaves of Nicotiana benthamiana and Nicotiana tabacum. Plant Mol. Biol. 2005, 59, 647–661. [Google Scholar] [CrossRef]
  21. Christie, M.; Croft, L.J.; Carroll, B.J. Intron splicing suppresses RNA silencing in Arabidopsis. Plant J. 2011, 68, 159–167. [Google Scholar] [CrossRef] [Green Version]
  22. Luo, Z.; Chen, Z. Improperly terminated, unpolyadenylated mRNA of sense transgenes is targeted by RDR6-mediated RNA silencing in Arabidopsis. Plant Cell 2007, 19, 943–958. [Google Scholar] [CrossRef] [PubMed]
  23. Dalakouras, A.; Wassenegger, M.; McMillan, J.N.; Cardoza, V.; Maegele, I.; Dadami, E.; Runne, M.; Krczal, G.; Wassenegger, M. Induction of silencing in plants by high-pressure spraying of in vitro-synthesized small RNAs. Front. Plant Sci. 2016, 7, 1327. [Google Scholar] [CrossRef]
  24. Tzfira, T.; Tian, G.W.; Lacroix, B.; Vyas, S.; Li, J.; Leitner-Dagan, Y.; Krichevsky, A.; Taylor, T.; Vainstein, A.; Citovsky, V. pSAT vectors: a modular series of plasmids for autofluorescent protein tagging and expression of multiple genes in plants. Plant Mol. Biol. 2005, 57, 503–516. [Google Scholar] [CrossRef]
  25. Dubrovina, A.S.; Kiselev, K.V.; Khristenko, V.S.; Aleynova, O.A. VaCPK20, a calcium-dependent protein kinase gene of wild grapevine Vitis amurensis Rupr., mediates cold and drought stress tolerance. J. Plant Physiol. 2015, 185, 1–12. [Google Scholar] [CrossRef]
  26. Varkonyi-Gasic, E.; Wu, R.; Wood, M.; Walton, E.F.; Hellens, R.P. Protocol: A highly sensitive RT-PCR method for detection and quantification of microRNAs. Plant Methods 2007, 3, 12. [Google Scholar] [CrossRef]
  27. Werner, N.S.; Windoffer, R.; Strnad, P.; Grund, C.; Leube, R.E.; Magin, T.M. Epidermolysis bullosa simplex-type mutations alter the dynamics of the keratin cytoskeleton and reveal a contribution of actin to the transport of keratin subunits. Mol. Biol. Cell 2004, 15, 990–1002. [Google Scholar] [CrossRef]
  28. Gruber, J.; Manninga, H.; Tuschl, T.; Osborn, M.; Weber, K. Specific RNAi mediated gene knockdown in zebrafish cell lines. RNA Biol. 2005, 2, 101–105. [Google Scholar] [CrossRef] [PubMed]
  29. Movahedi, A.; Sun, W.; Zhang, J.; Wu, X.; Mousavi, M.; Mohammadi, K.; Yin, T.; Zhuge, Q. RNA-directed DNA methylation in plants. Plant Cell Rep. 2015, 34, 1857–1862. [Google Scholar] [CrossRef]
  30. Cuerda-Gil, D.; Slotkin, R.K. Non-canonical RNA-directed DNA methylation. Nat. Plants 2016, 2, 16163. [Google Scholar] [CrossRef] [PubMed]
  31. Jones, L.; Hamilton, A.J.; Voinnet, O.; Thomas, C.L.; Maule, A.J.; Baulcombe, D.C. RNA-DNA interactions and DNA methylation in post-transcriptional gene silencing. Plant Cell 1999, 11, 2291–2301. [Google Scholar] [CrossRef]
  32. Morel, J.B.; Mourrain, P.; Béclin, C.; Vaucheret, H. DNA methylation and chromatin structure affect transcriptional and post-transcriptional transgene silencing in Arabidopsis. Curr. Biol. 2000, 10, 1591–1594. [Google Scholar] [CrossRef] [Green Version]
  33. Hohn, T.; Corsten, S.; Rieke, S.; Muller, M.; Rothney, H. Methylation of coding region alone inhibits gene expression in plant protoplasts. Proc. Natl. Acad. Sci. USA 1996, 93, 8334–8339. [Google Scholar] [CrossRef] [PubMed]
  34. Chawla, R.; Nicholson, S.J.; Folta, K.M.; Srivastava, V. Transgene-induced silencing of Arabidopsis phytochrome A gene via exonic methylation. Plant J. 2007, 52, 1105–1118. [Google Scholar] [CrossRef] [Green Version]
  35. Bhat, A.; Ryu, C.M. Plant perceptions of extracellular DNA and RNA. Mol. Plant 2016, 9, 956–958. [Google Scholar] [CrossRef] [PubMed]
  36. Lee, B.; Park, Y.S.; Lee, S.; Song, G.C.; Ryu, C.M. Bacterial RNAs activate innate immunity in Arabidopsis. New Phytol. 2016, 209, 785–797. [Google Scholar] [CrossRef]
  37. Niehl, A.; Wyrsch, I.; Boller, T.; Heinlein, M. Double-stranded RNAs induce a pattern-triggered immune signaling pathway in plants. New Phytol. 2016, 211, 1008–1019. [Google Scholar] [CrossRef] [Green Version]
  38. Zipfel, C. Plant pattern-recognition receptors. Trends Immunol. 2014, 35, 345–351. [Google Scholar] [CrossRef] [PubMed]
  39. Wang, M.; Weiberg, A.; Lin, F.M.; Thomma, B.P.; Huang, H.D.; Jin, H. Bidirectional cross-kingdom RNAi and fungal uptake of external RNAs confer plant protection. Nat. Plants 2016, 2, 16151. [Google Scholar] [CrossRef] [PubMed]
  40. Kiselev, K.V.; Dubrovina, A.S.; Shumakova, O.A.; Karetin, Y.A.; Manyakhin, A.Y. Structure and expression profiling of a novel calcium-dependent protein kinase gene, CDPK3a, in leaves, stems, grapes, and cell cultures of wild-growing grapevine Vitis amurensis Rupr. Plant Cell Rep. 2013, 32, 431–442. [Google Scholar] [CrossRef]
  41. Czechowski, T.; Stitt, M.; Altmann, T.; Udvardi, M.K.; Scheible, W.R. Genome-wide identification and testing of superior reference genes for transcript normalization in Arabidopsis. Plant Physiol. 2005, 139, 5–17. [Google Scholar] [CrossRef]
  42. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef]
  43. Schindelin, J.; Arganda-Carreras, I.; Frise, E.; Kaynig, V.; Longair, M.; Pietzsch, T.; Preibisch, S.; Rueden, C.; Saalfeld, S.; Schmid, B.; et al. Fiji: An open-source platform for biological-image analysis. Nat. Methods 2012, 9, 676–682. [Google Scholar] [CrossRef] [PubMed]
  44. Schneider, C.; Rasband, W.; Eliceiri, K. NIH image to ImageJ: 25 years of image analysis. Nat. Methods 2012, 9, 671–675. [Google Scholar] [CrossRef] [PubMed]
  45. Otsu, N. A threshold selection method from gray-level histograms. IEEE Trans Sys. Man. Cybern. 1979, 9, 62–66. [Google Scholar] [CrossRef]
  46. Sambrook, J.F.; Russell, D.W. Molecular Cloning: A Laboratory Manual, 3rd ed.; Cold Spring Harbor Laboratory Press: Cold Spring Harbor, NY, USA, 2001; 2100p. [Google Scholar]
  47. Ogneva, Z.V.; Dubrovina, A.S.; Kiselev, K.V. Age-associated alterations in DNA methylation and expression of methyltransferase and demethylase genes in Arabidopsis thaliana. Biol. Plant. 2016, 60, 628–634. [Google Scholar] [CrossRef]
  48. Peng, J.; Xia, Z.; Chen, L.; Shi, M.; Pu, J.; Guo, J.; Fan, Z. Rapid and efficient isolation of high-quality small RNAs from recalcitrant plant species rich in polyphenols and polysaccharides. PLoS ONE 2014, 9, e95687. [Google Scholar] [CrossRef] [PubMed]
  49. Kiselev, K.V.; Aleynova, O.A.; Grigorchuk, V.P.; Dubrovina, A.S. Stilbene accumulation and expression of stilbene biosynthesis pathway genes in wild grapevine Vitis amurensis Rupr. Planta 2017, 245, 151–159. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Characterization of transformed A. thaliana lines, and investigation of the dsRNA concentration effect. (a) Schematic representation of the T-DNA region from the pZP-RCS2-EGFP-NPTII vector [24]. 2x35S—the double 35S promoter of the cauliflower mosaic virus (CaMV); EGFP—the enhanced green fluorescent protein (EGFP) gene; NPTII—the neomycin phosphotransferase II (NPTII) gene; Tnos—nopaline synthase terminator. D1s, D1a, D2s, D2a—primers for dsRNA synthesis; R1s, R1a, R3s, R3a—primers for qRT-PCR estimation of transgene expression after dsRNA treatments; R2s, R2a, R4s, R4a—primers designed to align inside the transgene fragments that have been used for synthesis of the corresponding dsRNAs; M1s, M1a, M2s, M2a—primers for bisulfite DNA sequencing. (b,c) Quantification the NPTII and EGFP mRNAs in the untreated transgenic A. thaliana. RNA was extracted from the wild type (WT) and transformed A. thaliana lines (L1, L2, L3). (d,e) Quantification the NPTII and EGFP mRNAs (log2 fold change) in L1 of A. thaliana in response to external application of NPTII-dsRNA, and EGFP-dsRNA 7 days post-treatment, respectively. 10, 35, 100—the synthesized EGFP-dsRNA and NPTII-dsRNA were diluted in water to concentrations of 0.1, 0.35, and 1 µg/µL (100 µL per plant). Control—sterile filtered water. qRT-PCR data are presented as mean ± SE. *, **—significantly different from the untreated plants at P ≤ 0.05 and 0.01, respectively, according to the Student’s t-test.
Figure 1. Characterization of transformed A. thaliana lines, and investigation of the dsRNA concentration effect. (a) Schematic representation of the T-DNA region from the pZP-RCS2-EGFP-NPTII vector [24]. 2x35S—the double 35S promoter of the cauliflower mosaic virus (CaMV); EGFP—the enhanced green fluorescent protein (EGFP) gene; NPTII—the neomycin phosphotransferase II (NPTII) gene; Tnos—nopaline synthase terminator. D1s, D1a, D2s, D2a—primers for dsRNA synthesis; R1s, R1a, R3s, R3a—primers for qRT-PCR estimation of transgene expression after dsRNA treatments; R2s, R2a, R4s, R4a—primers designed to align inside the transgene fragments that have been used for synthesis of the corresponding dsRNAs; M1s, M1a, M2s, M2a—primers for bisulfite DNA sequencing. (b,c) Quantification the NPTII and EGFP mRNAs in the untreated transgenic A. thaliana. RNA was extracted from the wild type (WT) and transformed A. thaliana lines (L1, L2, L3). (d,e) Quantification the NPTII and EGFP mRNAs (log2 fold change) in L1 of A. thaliana in response to external application of NPTII-dsRNA, and EGFP-dsRNA 7 days post-treatment, respectively. 10, 35, 100—the synthesized EGFP-dsRNA and NPTII-dsRNA were diluted in water to concentrations of 0.1, 0.35, and 1 µg/µL (100 µL per plant). Control—sterile filtered water. qRT-PCR data are presented as mean ± SE. *, **—significantly different from the untreated plants at P ≤ 0.05 and 0.01, respectively, according to the Student’s t-test.
Ijms 20 01585 g001
Figure 2. dsRNA-induced EGFP and NPTII transgene suppression in four-week-old Arabidopsis thaliana. The expression levels (log2 fold change) of the NPTII (a) and EGFP (b) in response to external application of sterile filtered water (Control), NPTII-dsRNA, and EGFP-dsRNAs relative to NPTII and EGFP levels in the plants before treatments. The dsRNA were diluted in water to 0.35 µg/µL (100 µL per plant). The RNA was isolated before, 1, 7, and 14 days post-treatment from an individual plant of each line per treatment (four independent experiments). qRT-PCR data are presented as mean ± SE. *, **—significantly different from the untreated plants at P ≤ 0.05 and 0.01, respectively, according to the Student’s t-test.
Figure 2. dsRNA-induced EGFP and NPTII transgene suppression in four-week-old Arabidopsis thaliana. The expression levels (log2 fold change) of the NPTII (a) and EGFP (b) in response to external application of sterile filtered water (Control), NPTII-dsRNA, and EGFP-dsRNAs relative to NPTII and EGFP levels in the plants before treatments. The dsRNA were diluted in water to 0.35 µg/µL (100 µL per plant). The RNA was isolated before, 1, 7, and 14 days post-treatment from an individual plant of each line per treatment (four independent experiments). qRT-PCR data are presented as mean ± SE. *, **—significantly different from the untreated plants at P ≤ 0.05 and 0.01, respectively, according to the Student’s t-test.
Ijms 20 01585 g002
Figure 3. The analysis of NPTII (a) and EGFP (b) transgene transcript levels (log2 fold change) using primers designed to align inside the NPTII and EGFP transgene fragments, which have been used for synthesis of the corresponding dsRNAs, in the dsRNA-treated four-week-old Arabidopsis thaliana relative to that in the plants before treatments. One plant of the L1, L2, and L3 lines for each treatment was treated with NPTII-dsRNA, EGFP-dsRNAs or filtered sterile water (Control). The dsRNA were diluted in water to 0.35 µg/µL (100 µL per plant). In an experiment, the RNA was isolated before, 1, 7, and 14 days post-treatment from each plant (three independent experiments). qRT-PCR data are presented as mean ± SE. *, **—significantly different from the plants before treatment at P ≤ 0.05 and 0.01, respectively, according to the Student’s t-test.
Figure 3. The analysis of NPTII (a) and EGFP (b) transgene transcript levels (log2 fold change) using primers designed to align inside the NPTII and EGFP transgene fragments, which have been used for synthesis of the corresponding dsRNAs, in the dsRNA-treated four-week-old Arabidopsis thaliana relative to that in the plants before treatments. One plant of the L1, L2, and L3 lines for each treatment was treated with NPTII-dsRNA, EGFP-dsRNAs or filtered sterile water (Control). The dsRNA were diluted in water to 0.35 µg/µL (100 µL per plant). In an experiment, the RNA was isolated before, 1, 7, and 14 days post-treatment from each plant (three independent experiments). qRT-PCR data are presented as mean ± SE. *, **—significantly different from the plants before treatment at P ≤ 0.05 and 0.01, respectively, according to the Student’s t-test.
Ijms 20 01585 g003
Figure 4. The effect of external dsRNA application on the EGFP fluorescence and EGFP protein levels in four-week-old Arabidopsis plants of the L1 line. EGFP silencing efficiency was visualized 1, 7, and 14 days post-treatment of L1 rosettes using confocal microscopy (a) and western blotting (b). Pc—total protein purified from E. coli Rosetta(DE3)pLysS transformed with pET28a-EGFP. Values are the mean ± SE. ** ‒ significantly different from the negative control at P ≤ 0.01, respectively, according to the Student’s t-test.
Figure 4. The effect of external dsRNA application on the EGFP fluorescence and EGFP protein levels in four-week-old Arabidopsis plants of the L1 line. EGFP silencing efficiency was visualized 1, 7, and 14 days post-treatment of L1 rosettes using confocal microscopy (a) and western blotting (b). Pc—total protein purified from E. coli Rosetta(DE3)pLysS transformed with pET28a-EGFP. Values are the mean ± SE. ** ‒ significantly different from the negative control at P ≤ 0.01, respectively, according to the Student’s t-test.
Ijms 20 01585 g004
Figure 5. Detection of EGFP-siRNA (NCBI siRNA database acession: Pr008640320 or Pr120179) in the dsRNA-treated and water-treated A. thaliana L1 leaves using stem-loop RT-PCR before treatment and 7 days post application. (a) Agarose gel electrophoresis of the stem-loop RT-PCR products. M—low molecular weight DNA ladder; Nc—negative control (RT-PCR mixture without RNA). (b) Sequences of the cloned RT-PCR product and predicted stem-loop structure of the target EGFP-siRNA (sense strand) with primer positions. Dashed lines depict position of the stem-loop RT primer; highlighted in green and blue are the forward and reverse primer positions for stem-loop RT-PCR; highlighted in yellow is the EGFP-siRNA sense strand.
Figure 5. Detection of EGFP-siRNA (NCBI siRNA database acession: Pr008640320 or Pr120179) in the dsRNA-treated and water-treated A. thaliana L1 leaves using stem-loop RT-PCR before treatment and 7 days post application. (a) Agarose gel electrophoresis of the stem-loop RT-PCR products. M—low molecular weight DNA ladder; Nc—negative control (RT-PCR mixture without RNA). (b) Sequences of the cloned RT-PCR product and predicted stem-loop structure of the target EGFP-siRNA (sense strand) with primer positions. Dashed lines depict position of the stem-loop RT primer; highlighted in green and blue are the forward and reverse primer positions for stem-loop RT-PCR; highlighted in yellow is the EGFP-siRNA sense strand.
Ijms 20 01585 g005
Figure 6. Detection of the read-through transcripts, including both the EGFP and NPTII sequences, and analysis of the total cytosine DNA methylation, within the NPTII and EGFP transgenes in the transgenic L1 line of A. thaliana, before and 7 days after foliar application of dsRNAs or filtered sterile water. (a,b) Agarose gel electrophoresis of the RT-PCR products. C—control treatment (sterile filtered water). 10, 35, 100—treatment with the dsRNAs diluted in water to concentrations of 0.1, 0.35, and 1 µg/µL (100 µL per plant). M—low molecular weight DNA ladder; Nc—negative control (RT-PCR mixture without RNA); Pc—positive control (pZP-RCS2 plasmid). (c,d) Bisulfite sequencing of the DNA purified from the dsRNA- and water-treated plants. The dsRNAs were diluted in water to 0.35 µg/µL (100 µL per plant). The data are presented as mean ± SE. Significant differences at *—P < 0.05 and **—P < 0.01 versus values of cytosine methylation before treatment according to the Student’s t-test. The total number of cytosines was regarded as 100% in each DNA region.
Figure 6. Detection of the read-through transcripts, including both the EGFP and NPTII sequences, and analysis of the total cytosine DNA methylation, within the NPTII and EGFP transgenes in the transgenic L1 line of A. thaliana, before and 7 days after foliar application of dsRNAs or filtered sterile water. (a,b) Agarose gel electrophoresis of the RT-PCR products. C—control treatment (sterile filtered water). 10, 35, 100—treatment with the dsRNAs diluted in water to concentrations of 0.1, 0.35, and 1 µg/µL (100 µL per plant). M—low molecular weight DNA ladder; Nc—negative control (RT-PCR mixture without RNA); Pc—positive control (pZP-RCS2 plasmid). (c,d) Bisulfite sequencing of the DNA purified from the dsRNA- and water-treated plants. The dsRNAs were diluted in water to 0.35 µg/µL (100 µL per plant). The data are presented as mean ± SE. Significant differences at *—P < 0.05 and **—P < 0.01 versus values of cytosine methylation before treatment according to the Student’s t-test. The total number of cytosines was regarded as 100% in each DNA region.
Ijms 20 01585 g006
Table 1. Methylation status in the NPTII and EGFP transgene protein-coding regions in the transgenic L1 line of A. thaliana before and 7 days after foliar application of NPTII-dsRNA, EGFP-dsRNAs, or filtered sterile water.
Table 1. Methylation status in the NPTII and EGFP transgene protein-coding regions in the transgenic L1 line of A. thaliana before and 7 days after foliar application of NPTII-dsRNA, EGFP-dsRNAs, or filtered sterile water.
GeneTreatment CGCHGCHH
NPTIIWater (before)
Water (7 dpa)
65.6 ± 2.160.1 ± 3.355.9 ± 3.8
65.3 ± 3.263.3 ± 2.459.3 ± 4.7
NPTII-dsRNA (before)
NPTII-dsRNA (7 dpa)
68.5 ± 4.460.7 ± 3.051.6 ± 4.5
82.5 ± 4.9 **79.2 ± 3.4 **77.5 ± 3.6 **
EGFP-dsRNA (before)
EGFP-dsRNA (7 dpa)
65.3 ± 3.465.0 ± 3.755.6 ± 2.4
80.7 ± 4.378.2 ± 2.978.1 ± 2.2 **
EGFPWater (before)
Water (7 dpa)
62.1 ± 2.459.2 ± 3.154.7 ± 3.9
62.4 ± 4.259.5 ± 2.857.9 ± 5.2
NPTII-dsRNA (before)
NPTII-dsRNA (7 dpa)
65.2 ± 5.962.2 ± 3.359.3 ± 4.3
75.8 ± 5.775.9 ± 4.6 *73.9 ± 3.4 **
EGFP-dsRNA (before)
EGFP-dsRNA (7 dpa)
59.9 ± 4.256.1 ± 3.854.1 ± 3.2
74.9 ± 4.875.1 ± 4.7 **73.1 ± 3.0 **
The data are presented as the mean ± SE, * P < 0.05; ** P < 0.01 versus the values measured before treatment.

Share and Cite

MDPI and ACS Style

Dubrovina, A.S.; Aleynova, O.A.; Kalachev, A.V.; Suprun, A.R.; Ogneva, Z.V.; Kiselev, K.V. Induction of Transgene Suppression in Plants via External Application of Synthetic dsRNA. Int. J. Mol. Sci. 2019, 20, 1585. https://doi.org/10.3390/ijms20071585

AMA Style

Dubrovina AS, Aleynova OA, Kalachev AV, Suprun AR, Ogneva ZV, Kiselev KV. Induction of Transgene Suppression in Plants via External Application of Synthetic dsRNA. International Journal of Molecular Sciences. 2019; 20(7):1585. https://doi.org/10.3390/ijms20071585

Chicago/Turabian Style

Dubrovina, Alexandra S., Olga A. Aleynova, Alexander V. Kalachev, Andrey R. Suprun, Zlata V. Ogneva, and Konstantin V. Kiselev. 2019. "Induction of Transgene Suppression in Plants via External Application of Synthetic dsRNA" International Journal of Molecular Sciences 20, no. 7: 1585. https://doi.org/10.3390/ijms20071585

APA Style

Dubrovina, A. S., Aleynova, O. A., Kalachev, A. V., Suprun, A. R., Ogneva, Z. V., & Kiselev, K. V. (2019). Induction of Transgene Suppression in Plants via External Application of Synthetic dsRNA. International Journal of Molecular Sciences, 20(7), 1585. https://doi.org/10.3390/ijms20071585

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