- Methodology article
- Open Access
RNase If -treated quantitative PCR for dsRNA quantitation of RNAi trait in genetically modified crops
BMC Biotechnologyvolume 18, Article number: 3 (2018)
RNA interference (RNAi) technology has been widely used to knockdown target genes via post-transcriptional silencing. In plants, RNAi is used as an effective tool with diverse applications being developed such as resistance against insects, fungi, viruses, and metabolism manipulation. To develop genetically modified (GM) RNAi traits for insect control, a transgene is created and composed of an inversely-repeated sequence of the target gene with a spacer region inserted between the repeats. The transgene design is subject to form a self-complementary hairpin RNA (hpRNA) and the active molecules are > 60 bp doubled-stranded RNA (dsRNA) derived from the hpRNA. However, in some cases, an undesirable intermediate such as single-stranded RNA (ssRNA) may be formed, which is not an active molecule. The aforementioned characteristics of RNAi traits lead to increase the challenges for RNAi-derived dsRNA quantitation.
To quantify the dsRNA and distinguish it from the ssRNA in transgenic maize, an analytical tool is required to be able to effectively quantify dsRNA which contains a strong secondary structure. Herein, we develop a modified qRT-PCR method (abbreviated as RNase If -qPCR) coupled with a ssRNA preferred endonuclease (i.e., RNase If). This method enables the precise measurement of the active molecules (i.e., dsRNA) derived from RNAi traits of GM crops and separately quantifies the dsRNA from ssRNA. Notably, we also demonstrate that the RNase If -qPCR is comparable to a hybridization-based method (Quantigene Plex 2.0).
To our best knowledge, this is the first report of a method combining RNase If with modified qRT-PCR protocol. The method represents a reliable analytical tool to quantify dsRNA for GM RNAi crops. It provides a cost-effective and feasible analytical tool for general molecular laboratory without using additional equipment for other methods. The RNase If -qPCR method demonstrates high sensitivity (to 0.001 pg/ μL of dsRNA), precision and accuracy. In this report, we demonstrated the deployment of this method to characterize the RNAi events carrying v-ATPase C in maize during trait development process. The method can be utilized in any application which requires the dsRNA quantification such as double-stranded RNA virus or sprayable dsRNA as herbicide.
RNA interference (RNAi) mediates a posttranscriptional gene-silencing (PTGS) mechanism via the process of exogenous double-stranded RNA (dsRNA) silencing complementary to endogenous mRNA in the cell . Due to its high specificity in knockdown gene expression, RNAi has been developed into a widely used tool in plants for a variety of applications such as nematode resistance , insect resistance , fungi resistance , viral resistance , or parasitic weed resistance . The use of RNAi as a pest control tool was first discovered and tested using the non-parasitic nematode Caenorhabditis elegans  and has been extensively investigated in the past decades for the development of a new mode of action for genetic modified (GM) crops to control pests in addition to conventional Bt (Bacillus thuringiensis) traits (detailed reviews in [8, 9]).
The RNAi approach for insect control becomes a practical application in crop protection through the expression of hairpin RNAs (hpRNA) in GM crops to downregulate essential genes in target insects. Transgenic plants producing hpRNAs to control insect pests were first demonstrated in western corn rootworm (WCR) (Diabrotica virgifera)  and cotton bollworm (Helicoverpa armigera) . Both studies used a transgenic RNAi approach to produce plant material for oral ingestion. This allowed the dsRNA to silence the essential genes of the targeted insects for effective insect control. The key attributes of a successful RNAi approach for insect control are: (1) Identification of essential genes and (2) sufficient amount of dsRNA for delivery to the insect.
In the study by Mao et al. (2007) plants engineered to express dsRNAs against the bollworm CYP6AE14 gene showed effective larvae growth inhibition . More interestingly, when the transgene was introduced into Arabidopsis dicer mutants (e.g. dcl2, dcl3 and dcl4), it resulted in greater abundance of longer dsRNA, CYP6AE14 gene in the insect is repressed more effectively by transgenic plant materials. This result suggests that a minimal length of dsRNA is required to achieve the optimal efficiency of repression of targeted genes in pests. Moreover, two recent studies demonstrated the length of dsRNA is critical to achieve the desired efficacy. In the study described by Bolognesi et al., different sequence lengths of dsRNAs against Dv Snf7 were tested in bioassays, a size cut-off of approximately 60 bp dsRNA is required to display efficacy against corn rootworm . Another study reported by Li et al. also showed dsRNA against Dv v-ATPase C with at least 60 bp in length resulted in high levels of larval mortality . In contrast, feeding the pooled 21-bp siRNAs that have homology to v-ATPase C do not cause any mortality in both WCR larvae and adults.
Interestingly, the gene encoding RNA-dependent RNA polymerase (RdRP), required for the secondary siRNA amplification and transitive RNAi effects, is absent in insects  but is present in nematodes and plants [14, 15]. The lack of the RdRP activity may implicate that any effect of RNAi is a direct result from oral ingestion of dsRNA. A recent study displays evidence that the RNAi effect in WCR is not attributed to the transitive RNAi. In contrast, the siRNA distribution is solely restricted to sequences within the initial target sequence regions . Taken together, these results indicate that the sequence length and initial concentration of transgene-derived dsRNAs for oral ingestion are key attributes for insect control.
Since the length and concentration of dsRNA are critical attributes to achieve the efficacy for controlling WCR, a reliable analytical tool to quantify dsRNA is required for RNAi trait development. The approach of detecting dsRNAs can be challenging due to its strong secondary structure derived from self-complementary inverted repeats (IRs). The secondary structure of the RNA molecules is often unpredicted and is determined by specific base pairing interactions within the encoded nucleotide sequence . In addition, due to the design of transgene cassette, transgene may form dsRNA variants; some RNA variants may lead to intermediate forms of ssRNA, leading to non-efficacious molecules.
Studies in viral biology have developed some methods to characterize viral dsRNA such as immunofluorescence analysis or conventional RT-PCR. Using the monoclonal dsRNA-specific mouse antibody, the dsRNA may be detected via immunoblot or immunohistochemistry [18, 19]. Moreover, studies in viral biology use high (e.g., 96 °C) and low temperature (e.g., 70 °C) to denature the dsRNA and ssRNA, respectively . Recently, a hybridization-based method using the QuantiGene technology (Affymetrix) has been reported for effective dsRNA quantitation . In this report, we developed a modified reverse transcription PCR (RT-PCR) by coupling the RNase If -treatment RNA prior to cDNA conversion (abbreviation: RNaseIf -qPCR) which has proven to be both cost effective and feasible. RNase If is an endonuclease which has preferential activity to digest ssRNA over dsRNA. RNaseIf -qPCR allows users to precisely distinguish dsRNA from ssRNA and effectively assay the dsRNA in a high-throughput manner. The reported method here may be easily adopted for any applications in the emerging field of RNAi-based insect pest control with the need for dsRNA quantitation.
Synthetic double-stranded RNA and single-stranded RNA of Dv v-ATPase C
In this study, we used RNAi carrying sequence targeting vacuolar H + -ATPase (V-ATPase) subunit C in western corn rootworm (WCR) (Diabrotica virgifera) for our testing. V-ATPase is a highly conserved enzyme complex that acidifies intracellular organelles through pumping protons across the plasma membranes and consists of catalytic V1 and membrane-bound V0 two subcomplexes . V-ATPase C in WCR is one of the effective dsRNA target. Previous study has demonstrated the efficacy of dsRNA targeting Dv v-ATPase C . To validate the method, a 149-bp sequence (nucleotide positions 9-157 of Dv v-ATPase C open reading frame (ORF)) was selected for dsRNA synthesis . The dsRNA was synthesized (agroRNA, Seoul, Korea). For ssRNA, a 125-bp sequence overlapping with the 149-bp dsRNA region was selected for ssRNA synthesis. The ssRNA was synthesized through IDT DNA Technologies (Iowa). One hundred nanogram synthetic Dv v-ATPase C dsRNA was treated with DNase (Ambion, Life Technologies, Austin TX, USA) to eliminate DNA contamination. After inactivation of the DNase, the dsRNA was then diluted to 100 pg/μL. The synthetic ssRNA, without further DNase treatment, was diluted to 100 pg/μL with RNase free H2O. The 10-point standard dilution and mixtures of dsRNA, ssRNA and wild-type maize B104 leaf RNA were created at the different ratios as specified in Tables 1, 2 and 3.
Transgenic Dv v-ATPase C RNAi maize materials
The transgenic maize events producing Dv v-ATPase C hairpin RNA (hpRNA) were created based on the vector previously described . Briefly, the binary vector carrying a 149-bp sequence (nucleotide positions 9-157 of Dv v-ATPase C open reading frame (ORF)) of Dv v-ATPase C sequence was constructed. The binary vector was designed to express a hpRNA with 149-bp invert-repeated sequence with a loop sequence in between the inverted repeats. Transcription of the v-ATPase C mRNA was driven by the maize uniqutin 1 promoter and terminated by a 3′ untranslated region (UTR) from a maize peroxidase 5 gene (ZmPer5 3’UTR). The binary vector was then transformed into Agrobacterium tumefaciens. The binary vector was transformed into maize via Agrobacterium-mediated transformation of immature embryos isolated from the inbred line, Zea mays c.v. B104 following the standard method .
RNase If -qPCR
Plant RNA extractions were carried out using the Norgen Plant/Fungi RNA purification kit (Norgen Biotek, Thorold, ON, Canada). Four leaf punches of maize plants at v4 stages were collected on ice into a 96-well collection plate. Samples were then ground in lysis buffer for 5 min using a 3.2 mm stainless steel bead. The supernatant was separated from the plant debris and precipitated with 100% ETOH. The subsequent binding and wash steps followed the manufacturer’s procedure with an on-column DNase I (Norgen) treatment for 15 min. Once eluted, the isolated RNA was then quantified on the NanoDrop-8000 (Thermo Fisher Scientific, Waltham, MA, USA) and normalized to 50 ng/μL.
RNase If treatment
The isolated RNA was digested with RNase If (New England Biolabs, Ipswich, MA, USA) per the manufacturer’s protocol at 37 °C for 10 min followed by heat inactivation at 70 °C for 20 min. The digested samples were purified using the RNA Clean-up and concentration kit (Norgen) following the manufacturer’s protocol. The RNase If treated and purified dsRNA were then used for cDNA synthesis.
First-strand cDNA synthesis
High capacity cDNA RT kit (Life Technologies, Carlsbad, CA, USA) was used for qRT-PCR analysis with the following modifications. (1) The dsRNA is pre-incubated at 95 °C and 70 °C for 5 min in the presence of the random hexamers and flash frozen with ice. (2) Multiscribe™ reverse transcriptase was then added to the pre-incubated RNA:hexamers mixture, then follow the manufacture protocol for cDNA conversion.
For the synthetic dsRNA/ssRNA mixture or leaf total RNA from transgenic plants, an aliquot of RNA was treated by RNase If to digest the ssRNA while another aliquot without RNase If treatment containing both the ssRNA and dsRNA was set aside. The same proportion of RNase If -treated and non-treated RNA was used for downstream cDNA conversion, respectively with the pre-incubation of RNA and random hexamers described above. For example, 1 μg of total RNA was digested via RNase If and eluted in 50 μL after clean-up. Five microliter (~ 100 ng) of RNaseIf-treated RNA was then used for cDNA conversion (both 95 °C and 70 °C); in a parallel set, 100 ng of untreated RNA was directly used for cDNA conversion (both 95 °C and 70 °C).
After the cDNA conversion using the RNA from different treatments, the ensuing cDNA was quantified via qRT-PCR using the LightCycler®480 system (Roche Diagnostics, Indianapolis, IN, USA) with a hydrolysis probe. The assays were run for the target gene (i.e., Dv v-ATPase C in this study) and the Zea mays TIP41 as the reference gene. For cDNA converted from the RNA without being treated by RNase If, the assays were performed in biplex reaction (i.e., Dv v-ATPase C and TIP41). For cDNA from RNA treated with RNase If, only Dv v-ATPase C was assayed in the singleplex reaction. To normalize the dsRNA, a separate reaction using the same proportion of RNA without RNase If treatment was used, as the TIP41 signal from a separate non RNase If - treated RNA was used to normalize dsRNA. To prepare the assays, LightCycler®480 Probes Master (Roche) was prepared at 1× final concentration in a 10 μL volume containing 0.4 μM of each primer and 0.2 μM of each probe. TIP41: (Forward primer: TGAGGGTAATGCCAACTGGTT; Reverse primer: GCAATGTAACCGAGTGT CTCTCAA and Probe: HEX-TTTTTGGCTTAGAGTTGATGGTGTACTGATGA-BHQ1). Dv v-ATPase C (Forward primer: GGGACACGATGAATAATATGACAAGTA; Reverse primer: CCAACTTTCAAGTCCGGGATT and Probe: FAM-CAACCAACTACAAGTTTC-MGB).
A two-step qRT-PCR program was performed with 40 cycles of the denaturation step at 95 °C for 10 s and extension step at 58 °C for 40 s with fluorescence acquisition. Cp values, the point at which the florescence signal crosses the background threshold using the fit points algorithm (LightCycler® software release 1.5) and the relative transcript level (RTL) was calculated using 2^ (−∆Ct) where ∆Ct = Ct (target gene) - Ct (reference gene) to obtain the relative quantitation over the endogenous gene.
Sample preparation, digestion & denaturation
Synthetic ssRNA and dsRNA was diluted to 100 pg with 5 μg of wild-type maize B104 leaf total RNA. An aliquot of both was treated with RNase If to remove the ssRNA and another aliquot was left un-treated. The samples were then precipitated and resuspended in 20 μL of RNA sample loading buffer (Sigma-Aldrich, St. Louis, MO, USA). Five microliter (5 ng) of 0.3-1.5Kb DIG-labeled RNA molecular weight marker III (Roche) was aliquoted and added to 15 μL of RNA sample loading buffer. Samples and marker were denatured at 70 ͦ C for 15 min and immediately placed on ice.
Gel preparation and running
A 1.5% denaturing gel was created using Roche Mp agarose (Roche) with the addition of 6.5% formaldehyde in 1× MESA buffer (MOPS-EDTA-sodium acetate). Gel was cast and run in electrophoresis equipment capable of running samples at 65 V/30 mA for 2 h and 30 min.
Following electrophoresis, the gel was rinsed in 2X SSC for 15 min. The gel was passively transferred to a nylon membrane (Millipore, Billerica, MA, USA) overnight at RT, using 10X SSC as the transfer buffer. Following the transfer, the membrane was rinsed in 2X SSC for 5 min, subsequently crosslinked (Stratagene, San Diego, CA), and allowed to dry at RT for up to 2 days.
Membrane hybridization, processing, and size determination
The membrane was prehybridized in UltraHyb buffer (Ambion) for 1-2 h at 65 °C. Blots described were probed with a 149-nt antisense single-stranded RNA complementary to the sense-strand RNA of the v-ATPase C hpRNA labelled with DIG digoxigenin-11-UTP using the Roche DIG RNA labeling Kit (Roche) according to manufacturer’s protocols . Probes at a concentration of 80 ng/ml of hybridization solution were denatured at 95 °C for 5 min and placed on ice for at least 2 min before adding to fresh hybridization solution. Hybridization occurred overnight at a temperature of 65 °C in hybridization tubes. Following hybridization, the blots were subjected to the DIG immunological washes and detection (DIG wash and block buffer set, Roche) following manufacturer’s instructions. The wrapped blot was put up film for 1-20 min, developed and scanned.
QuantiGene Plex 2.0 assay
The QuantiGene Plex 2.0 assay (QGP) assay (Affymetrix, Thermo Fisher Scientific, Waltham, MA, USA) used in this study was the custom triplex panel QGP Corn 41,669 3 Plx carrying three plex assay of Dv v-ATPase C, Elf1α and TIP like-41 (TIP41). One hundred twenty-five nanogram of total RNA from the same RNA source that was used for RNase If -qPCR was analyzed as assay input for QGP assay. The QGP assay was performed following the standard protocol of QGP except with the additional step to incubate RNA with probe sets. Briefly, the purified RNA was incubated with the probe sets at 95 °C for 5 min prior to mixing with the capture beads for overnight hybridization . A typical probe set for a single target mRNA consists of a family of four or more different capture extenders (CEs), label extenders (LEs) and blocking probes (BPs). After the overnight hybridization, the analytes were washed of excess probes and remaining purified RNA, and the signal amplification reagents, consisting of the preamplifier (PreAmp), amplifiers (Amp), and label probes (LPs), were sequentially hybridized to the captured target RNAs via specific Luminex beads. The LP also includes a biotin molecule, which in turn is a binding site for the final signal amplification reagent, streptavidin-conjugated R-phycoerythrin (SAPE). Last, the Luminex 200 Flow Cytometer instrument was used to read the resulting fluorescence signal associated with individual capture beads.
The QGP signal is reported as median fluorescence intensity (MFI) and is proportional to the number of target RNA molecules extracted from the purified RNA of the samples. MFI was calculated by measuring signals on 35-100 beads/target. MFI values obtained from background blank wells with no target RNA were subtracted from MFI of each target reading. The relative quantification of Dv v-ATPase C was calculated by dividing the Dv v-ATPase C (MFI ATPase C - MFI background) to TIP41 (MFI TIP41 - MFI background).
The RNAi transgene design is composed of inverted-repeat sequences of the target gene inserted into the genome with a spacer loop region connecting the repeats. The mRNA transcripts from such transgenes are subject to form a self-complementary hairpin RNA (hpRNA) structure with a based-paired stem encoded by the inverted repeats and single stranded loop sequence. However, in addition to the hpRNA, the mRNA transcripts may sometimes form an unpredicted intermediate such as partial single-stranded mRNA (ssRNA). Because of the strong secondary structure of the hpRNA and the potential for unexpected ssRNA formation, challenges for reliable RNAi analytics during the development of RNAi trait could be encountered. Since the dsRNA strands are the active molecules for efficacy [11, 12], an analytical method which distinguishes dsRNA from ssRNA and effectively quantifies the dsRNA molecules is required for a RNAi trait development. To provide a reliable method to quantify the active RNA molecule (i.e., dsRNA longer than 60 bp) of an RNAi trait in transgenic plants, we developed the RNase If -qPCR assay. Figure 1 shows the assay overview scheme for the RNase If -qPCR assay (see protocol in methods). Briefly, RNA is first extracted from the plant samples. Then, the isolated RNA is aliquoted and treated with RNase If endonuclease to remove ssRNA, while an aliquot of RNA without RNase If was set aside. After the RNase If treatment, the remaining dsRNA (Fig. 1b) is further purified and converted to cDNA. For RNA aliquots treated with RNase If, the RNA is pre-incubated with random hexamers at both 95 °C and 70 °C, and respective samples are converted to cDNA using standard protocols (Fig. 1c and d). The cDNA is then quantified using quantitative real-time PCR (qRT-PCR). The qRT-PCR result from the 95 °C pre-incubated cDNA conversion represents quantification of only dsRNA (abbreviation R95 RTL(relative transcript level)). The qRT-PCR result from the 70 °C pre-incubated cDNA conversion (abbreviation R70 RTL) serves as the negative control for RNase If treatment. If the RNase If treatment was effective, these samples (R70 RTL) should yield no detection as all the ssRNA should have been digested and the 70 °C is insufficient to open the dsRNA duplex for annealing of the random hexamers inside the dsRNA region being assayed. In addition to the RNase If -treated RNA, an equal amount of the non- RNase If -treated RNA is used directly for cDNA conversion with random hexamers at the same 95 °C and 70 °C pre-incubation. The result from the 70 °C pre-incubated cDNA conversion detects and quantifies the ssRNA present in the samples (abbreviation 70 RTL), which resembles regular qRT-PCR protocol. The 70 RTL and R95 RTL result may serve as a direct comparison between RNase If -qPCR and regular qPCR. The result from the 95 °C pre-incubated cDNA conversion represents the total of both the dsRNA and ssRNA present in the sample (abbreviation 95 RTL). Both ssRNA and dsRNA RTL were calculated using the endogenous gene quantified by 70 RTL treatment.
RNase If -qPCR assay demonstrates dilutional linearity, accuracy and precision
Since the whole procedure involves multiple steps, we first performed linearity testing to evaluate if the dsRNA or ssRNA analyzed by this method can display appropriate linearity. For this test, synthetic dsRNA and ssRNA were used for assay validation. A 149-bp sequence which targets the Dv v-ATPase C gene was selected for dsRNA synthesis. An ssRNA with overlapping region to the dsRNA was also synthesized (Fig. 2a). To confirm the secondary structure of dsRNA, a northern blot was used (Fig. 2b). The synthetic dsRNA and ssRNA were mixed with wild-type (WT) maize leaf RNA, separately or together, and then treated with RNase If to remove the ssRNA. RNA blot analysis for the ssRNA revealed no bands after the RNase If treatment, thus confirming complete ssRNA digestion via RNase If treatment. ssRNA was not detected after RNase If treatment, suggesting the capability of ssRNA digestion via RNase If (lane1 and 2). True to expectation, the dsRNA remained detected after the RNase If treatment with the RNA blot (lane3 and 4). This confirms the dsRNA duplex structure and validates the specificity of the RNase If treatment. The mixed dsRNA and ssRNA were still detected after the RNase If treatment, suggesting the dsRNA digestion is not affected in the presence of ssRNA (lane5 and 6).
To test for dilutional linearity, a 10-point dilution series of dsRNA and ssRNA were spiked into WT maize RNA, respectively. The samples were treated through the RNase If -qPCR assay (the corresponding dsRNA/ssRNA input in the final qRT-PCR is listed in Table 1). In the 10-point serial dilution of dsRNA, the R95 RTL to detect the dsRNA displayed dilutional linearity (Fig. 3a). The R70 RTL displayed no signal, indicating an effective RNase If digestion (Table 1). Much like for the dsRNA, the 10-point serial dilution of ssRNA, the 70 RTL to detect the ssRNA also displayed exceptional linearity (R2 0.99) (Fig. 3b).
Using the 10-point dilution of synthetic dsRNA and ssRNA as outlined above, the post-extraction spike recovery was assessed, which may be inferred for the accuracy of the assay. The amount of 10-point dilution (~ 0.001 – 0.667 pg/ μL in the final qRT-PCR input), the corresponding RTL (R95 RTL, R70RTL and 70RTL) and % recovery are presented in Table 1. The accuracy of the assay is indicated as percent recovery that was calculated as a percentage of the back-calculated measurement against the standard curve over the theoretical amount. As shown in Table 1, the % recovery for each standard point is within the acceptance range (80 – 120%) except a few points slightly above (e.g. 121%) and below (75%). Altogether, the linearity testing suggests that the RNase If -qPCR assay demonstrates dilutional linearity with acceptable accuracy over a range of concentration. Also, the method demonstrates exceptional sensitivity. The limit of quantitation (LOQ) for Dv v-ATPase C dsRNA was determined from the lowest amount from the 10-point dilution as ~ 0.001 pg/ μL.
To further assess the repeatability of the assay (i.e., inter-assay precision), we created 10-point dilution of dsRNA and repeated the assays three times; each by different analysts. The Cp values from qRT-PCR of each standard point was captured and the coefficient of variation (% CV) was calculated as inter-assay precision using the Cp values of each standard point observed from three assay runs (Table 2 and Additional file 1: Figure S1). As shown in Table 2, average % CV among three independent runs are ~ 5% (acceptance criteria 20%), suggesting the acceptable repeatability for this assay. Based on the aforementioned tests in linearity, accuracy and precision, the results are in compliance with ICH guideline for validation analytical procedures .
RNase If -qPCR assay determines the concordant ratios of dsRNA and ssRNA mixtures
The RNase If -qPCR assay is developed not only to effectively quantify dsRNA but also as a means to distinguish dsRNA from ssRNA. To test the assay for this capability, we mixed distinct ratios of dsRNA and ssRNA together with WT maize leaf RNA to resemble the varied amounts of dsRNA and ssRNA possible in transgenic plants and performed the RNase If -qPCR assay to test the accuracy of this assay. The dsRNA and ssRNA were measured from each distinct ratios of mixtures and the result are shown in Table 3. The nine different ratios of dsRNA and ssRNA mixtures, arranging from 1:5 to 10:1, were analyzed via RNase If -qPCR assay with at least two independent assays. As shown in Table 3, the RNase If -qPCR assay was able to determine the correct ratios of dsRNA and ssRNA for each mixture (Chi-square analysis showed no significant difference between the expected and measured ratios), demonstrating that the assay is capable of accurately measuring and distinguishing the dsRNA from ssRNA.
The practical applicability of RNase If -qPCR assay to characterize transgenic maize materials
To demonstrate the practical applicability of this assay for deployment in trait development process, we tested transgenic maize containing an RNAi transgene. Transgenic events carrying Dv v-ATPase C RNAi were produced using the method described (details in materials and methods). RNA was isolated from leaf tissues of distinct T0 transgenic events, and were analyzed via the RNase If -qPCR assay to quantify the dsRNA and ssRNA RTL, respectively (Fig. 4). Within the distinct T0 transgenic events, different ratios of dsRNA and ssRNA were identified using this assay. As shown in Fig. 4, using this method, we could identify the varied expression levels of Dv v-ATPase C including both dsRNA and ssRNA: with some events having high abundance of dsRNA, while some exhibiting low expression, which is often observed in regular protein-coding trait a s well due to inserted genomic location. Importantly, using this method, we could distinguish the abundance of active molecules (i.e., dsRNA) from ssRNA, .
To further investigate the dosage effects using this method, we deployed the assay in a segregating S1 population including homozygous and hemizygous transgenic events. Distinct transgenic events from different constructs were assayed. In each transgenic event, the homozygous plants showed ~ 1.6-2.0 fold increase in dsRNA expression as compared to the corresponding hemizygous plants, suggesting that the assay may quantify the dsRNA with an accurate dose response (Fig. 5). Altogether, the assay as deployed on transgenic plants indicates that the RNase If -qPCR assay delivers reliable quantification to the active molecules of an RNAi transgene.
Comparison of RNase If -qPCR assay quantification and QuantiGene Plex 2.0 RNA assay
The utility of QuantiGene Plex 2.0 (QGP) platform for dsRNA detection has been previously shown . The QuantiGene Plex 2.0 RNA assay is a platform that utilizes hybridization-based technology in which the target mRNAs are captured via custom designed probe sets coupled with magnetic microbeads. The detection of the signal results from the amplification of the captured oligonucleotides-RNA complex via branched DNA technology. To evaluate the concordance of QGP and the RNase If -qPCR, a custom QGP assay was designed based on the Dv v-ATPase C RNAi transgene in maize and the endogenous gene (e.g. TIP41). Purified leaf RNA from different transgenic maize events expressing Dv v-ATPase C RNAi transgene was used for cross comparison. RNase If -qPCR displays high concordance with QGP via the RTL measurements of the Dv v-ATPase C dsRNA expression among distinct transgenic events (Fig. 6). Regression analysis illustrates that RNase If -qPCR data had a significantly high correlation with QGP data: with the correlation coefficient R and slope for RNase If –qPCR versus QGP of 0.89 (p < 0.0001) and 1.07, respectively. Together, it suggests that the RNase If -qPCR assay can serve as an alternative analytical tool for dsRNA detection for RNAi trait in GM crops.
The application of RNAi in plants has been extensively explored in the past decades. Considering the application of RNAi for insect control was successfully demonstrated , this may lead to revolutionary benefits to agriculture for effective insect pest control with a distinct mode of action. Consequently, a reliable and precise analytical tool for RNAi is essential for RNAi trait development. In this study, we described a new modified qRT-PCR assay coupled with RNase If treatment. This new assay allows for the exclusive detection of the active RNAi molecules present in transgenes (i.e., dsRNA; > 60 bp). Through the assay validation process, we demonstrated that the assay displays dilutional linearity, precision and accuracy. Importantly, it can accurately distinguish between dsRNA from ssRNA generating concordant dsRNA-to-ssRNA ratios. This allows for the accurate quantitation of the active molecules (i.e., dsRNA) vs. non-efficacious intermediate (ssRNA). Applied to transgenic maize plants, the assay identifies events containing a significant level of Dv v-ATPase C RNAs which are present in single-stranded form (ssRNA). Without the ability to differentiate dsRNA from ssRNA, results may be misleading with inaccurate levels of active RNAi molecules (i.e., dsRNA). Thus, the described method here may serve as an effective tool to quantify dsRNA for RNAi traits in GM crops. In addition to the utilization of this method to characterize RNAi traits in crops, the reported method may be used for potential GM detection of the active ingredient as RNA-based approaches.
The success of this new dsRNA assay involves two critical steps: (1) Effective RNase If treatment and (2) 95 °C pre-incubation with random hexamers prior to cDNA conversion. Because it is important to ensure effective RNase If digestion, a QC step using 70 °C pre-incubation (R70 RTL) is important to evaluate if there are any undigested ssRNA remaining. Secondly, due to the strong intramolecular binding affinity of dsRNA, a 95 °C pre-incubation helps to ensure the dsRNA is open and accessible to the random hexamers prior to cDNA conversion. In addition to those two steps, because the RNase If -treated RNA lacks ssRNA, the RTL is normalized to the reference gene from a separate preparation produced without RNase If -treated RNA. Such cross normalization is not unprecedented used for qRT-PCR . For traditional qRT-PCR using SYBR® Green I, cross normalization is required for relative quantification due to the singleplex nature of SYBR® Green I, and the results are still considered valid. In this assay, the equal proportions of RNA input for both the RNase If and non-RNase If cDNA conversions is particularly important in order to obtain reliable results due to the fact the RNA samples are treated separately. Overall, the method represents a cost-effective and feasible solution to quantify and differentiate dsRNA and may be leveraged into any application which requires for dsRNA quantitation.
In conclusion, the RNase If -qPCR assay developed here provides an alternative method to effectively quantify the active molecules (i.e., dsRNA) of RNAi traits in GM crops. The modifications to conventional qPCR methods with the aid of a specific endonuclease and modified cDNA conversion protocol without additional equipment other than those used for conventional qPCR makes this method both feasible and cost-effective for all laboratories. This assay demonstrates equivalence to another hybridization-based method (i.e., QGP). Overall, the simplicity of this method makes it an alternative method for any application with the need of dsRNA quantitation such as the study of dsRNA virus in microbes.
Polymerase chain reaction
QuantiGene Plex 2.0
- RNase If -qPCR:
RNase If -treated quantitative real-time PCR
Western corn rootworm
Baulcombe D. RNA silencing in plants. Nature. 2004;431(7006):356–63.
Sindhu AS, Maier TR, Mitchum MG, Hussey RS, Davis EL, Baum TJ. Effective and specific in planta RNAi in cyst nematodes: expression interference of four parasitism genes reduces parasitic success. J Exp Bot. 2009;60(1):315–24.
Baum JA, Bogaert T, Clinton W, Heck GR, Feldmann P, Ilagan O, Johnson S, Plaetinck G, Munyikwa T, Pleau M, Vaughn T, Roberts J. Control of coleopteran insect pests through RNA interference. Nat Biotech. 2007;25(11):1322–6.
Eschen-Lippold L, Landgraf R, Smolka U, Schulze S, Heilmann M, Heilmann I, Hause G, Rosahl S. Activation of defense against Phytophthora Infestans in potato by down-regulation of syntaxin gene expression. New Phytol. 2012;193(4):985–96.
Krubphachaya P, Juricek M, Kertbundit S. Induction of RNA-mediated resistance to papaya ringspot virus type W. BMB Rep. 2007;40(3):404–11.
Aly R, Cholakh H, Joel DM, Leibman D, Steinitz B, Zelcer A, Naglis A, Yarden O, Gal-On A. Gene silencing of mannose 6-phosphate reductase in the parasitic weed Orobanche Aegyptiaca through the production of homologous dsRNA sequences in the host plant. Plant Biotechnol J. 2009;7(6):487–98.
Fire A, Xu S, Montgomery MK, Kostas SA, Driver SE, Mello CC. Potent and specific genetic interference by double-stranded RNA in Caenorhabditis Elegans. Nature. 1998;391(6669):806–11.
Katoch R, Sethi A, Thakur N, Murdock L. RNAi for insect control: current perspective and future challenges. Appl Biochem Biotechnol. 2013;171(4):847–73.
Yu N, Christiaens O, Liu J, Niu J, Cappelle K, Caccia S, Huvenne H, Smagghe G. Delivery of dsRNA for RNAi in insects: an overview and future directions. Insect Sci. 2013;20(1):4–14.
Mao Y-B, Cai W-J, Wang J-W, Hong G-J, Tao X-Y, Wang L-J, Huang Y-P, Chen X-Y. Silencing a cotton bollworm P450 monooxygenase gene by plant-mediated RNAi impairs larval tolerance of gossypol. Nat Biotech. 2007;25(11):1307–13.
Bolognesi R, Ramaseshadri P, Anderson J, Bachman P, Clinton W, Flannagan R, Ilagan O, Lawrence C, Levine S, Moar W, Mueller G, Tan J, Uffman J, Wiggins E, Heck G, Segers G. Characterizing the mechanism of action of double-stranded RNA activity against western corn rootworm (Diabrotica Virgifera Virgifera LeConte). PLoS One. 2012;7(10):e47534.
Li H, Khajuria C, Rangasamy M, Gandra P, Fitter M, Geng C, Woosely A, Hasler J, Schulenberg G, Worden S. Long dsRNA but not siRNA initiates RNAi in western corn rootworm larvae and adults. J Appl Entomol. 2015;139(6):432–5.
Tomoyasu Y, Miller SC, Tomita S, Schoppmeier M, Grossmann D, Bucher G. Exploring systemic RNA interference in insects: a genome-wide survey for RNAi genes in Tribolium. Genome Biol. 2008;9(1):R10.
Xie Z, Fan B, Chen C, Chen Z. An important role of an inducible RNA-dependent RNA polymerase in plant antiviral defense. Proc Natl Acad Sci U S A. 2001;98(11):6516–21.
Pak J, Fire A. Distinct populations of primary and secondary effectors during RNAi in C. Elegans. Science. 2007;315(5809):241–4.
Li H, Bowling AJ, Gandra P, Rangasamy M, Pence HE, McEwan RE, Khajuria C, Siegfried BD, Narva KE. Systemic RNAi in western corn rootworm, Diabrotica Virgifera Virgifera, does not involve transitive pathways. Insect Sci 2016. https://doi.org/10.1111/1744-7917.12382.
Tinoco I, Bustamante C. How RNA folds. J Mol Biol. 1999;293(2):271–81.
Lukács N. Detection of virus infection in plants and differentiation between coexisting viruses by monoclonal antibodies to double-stranded RNA. J Virol Methods. 1994;47(3):255–72.
Weber F, Wagner V, Rasmussen SB, Hartmann R, Paludan SR. Double-stranded RNA is produced by positive-strand RNA viruses and DNA viruses but not in detectable amounts by negative-strand RNA viruses. J Virol. 2006;80(10):5059–64.
Antunes TFS, Amaral RJV, Ventura JA, Godinho MT, Amaral JG, Souza FO, Zerbini PA, Zerbini FM, Fernandes PMB. The dsRNA virus papaya meleira virus and an ssRNA virus are associated with papaya sticky disease. PLoS One. 2016;11(5):e0155240.
Armstrong TA, Chen H, Ziegler TE, Iyadurai KR, Gao A-G, Wang Y, Song Z, Tian Q, Zhang Q, Ward JM, Segers GC, Heck GR, Staub JM. Quantification of Transgene-derived double-stranded RNA in plants using the QuantiGene nucleic acid detection platform. J Agric Food Chem. 2013;61(51):12557–64.
Forgac M. Vacuolar ATPases: rotary proton pumps in physiology and pathophysiology. Nat Rev Mol Cell Biol. 2007;8(11):917–29.
Beringer J, Chen W, Garton R, Sardesai N, Wang P-H, Zhou N, Gupta M, Wu H. Comparison of the impact of viral and plant-derived promoters regulating selectable marker gene on maize transformation and transgene expression. Plant Cell Rep. 2017;36(4):519–28.
Q2(R1) I. Validation of analytical procedures: text and methodology. 2005. http://www.ich.org/products/guidelines.html.
Wang Y, Zhu W, Levy DE. Nuclear and cytoplasmic mRNA quantification by SYBR green based real-time RT-PCR. Methods. 2006;39(4):356–62.
We would like to thank Dr. Marcelo German (Dow AgroSciences LLC) for critical review and comments of the manuscript. Ryan Lee (Dow AgroSciences LLC) for coordination of field sampling. Nitigna Patel (Dow AgroSciences LLC) for high throughput PCR assays. Zhifang Gao (Dow AgroSciences LLC) for initial discussion and advice.
This work was supported by Dow AgroSciences LLC.
Availability of data and materials
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
RNase If -qPCR assay inter-assay precision. A 10-point dilution of dsRNA was analyzed by different analysts for at least three times. Cp and log2 of dsRNA amount for each dilution were shown. (DOCX 42 kb)