- Research note
- Open Access
Optimal RNA isolation method and primer design to detect gene knockdown by qPCR when validating Drosophila transgenic RNAi lines
BMC Research Notesvolume 10, Article number: 647 (2017)
RNA interference is employed extensively in Drosophila research to study gene function within a specific cell-type or tissue. Thousands of transgenic Drosophila lines have been generated to express double stranded RNA for gene knockdown; however, no standardized method exists for quantifying their knockdown efficiency. Since antibodies are not available for many proteins, quantitative real-time PCR is often used. Here, we explore how primer design and RNA isolation method can influence detection of gene knockdown using qPCR.
We tested differences in detected gene knockdown efficiency when using purified polyadenylated mRNA or total RNA as templates for cDNA synthesis. We also tested two different primer locations for each gene: one to amplify a region 5′ of the RNAi cut site, and one to amplify a region 3′ of the cut site. Consistently, the strongest gene knockdown was detected when qPCR was performed using 5′ primer sets in combination with mRNA-derived cDNA. Our results indicate that detection of undegraded mRNA cleavage fragments can result in underestimation of true knockdown efficiency for a RNAi construct. Purification of polyadenylated mRNA, combined with primers designed to amplify the non-polyadenylated 5′ mRNA cleavage fragment can avoid this problem.
RNA interference (RNAi) is an important technique for performing loss-of-function experiments in both cell culture and in vivo models. RNAi-induced gene knockdown can give insight into the function of a gene and its encoded protein, within the context of a specific cell-type or tissue. The endogenous RNAi pathway is triggered when a double stranded RNA (dsRNA) molecule is present within a cell . The Dicer-2 protein cleaves dsRNA into short interfering RNA (siRNA) fragments. Argonaute and the RNA-induced silencing complex associate with these siRNA fragments to cleave target mRNA at a sequence-specific site. This cleavage is directed by the complementarity of the siRNA to the mRNA, leaving 5′ and 3′ mRNA fragments that are eventually degraded, and do not contribute to protein production. Current evidence suggests that the 5′ fragment is degraded through a nonstop mRNA decay pathway, while the 3′ fragment is degraded by cellular 5′ to 3′ exonucleases [2, 3].
RNAi is employed extensively in Drosophila research through use of the UAS-Gal4 binary transgenic expression system . Thousands of transgenic lines are available that express dsRNA hairpins under the control of the yeast upstream activation sequence (UAS) enhancer [5, 6]. When these UAS-RNAi lines are crossed to a Gal4 ‘driver’ line, the expression of dsRNA leads to activation of cellular RNAi machinery, which can induce tissue-specific knockdown of a given gene. This system has become an important tool for investigating gene function in specific tissues within the context of a whole organism. Despite the importance and widespread use of Drosophila RNAi reagents, there are still several challenges and limitations to consider. RNAi-mediated gene knockdown is not 100% efficient and knockdown levels can vary considerably between various RNAi lines . This can lead to inconsistencies in phenotypic effects for different RNAi lines that target the same gene. Therefore, gene knockdown should be quantified for different RNAi lines. Western blotting is an ideal technique for verifying gene knockdown; however, this method is not effective for tissue-specific RNAi knockdown. In these cases, one must rely on immunohistochemistry, which is only partially quantitative. Furthermore, specific and effective antibodies are often not available for proteins of interest. They are also expensive and sometimes difficult to generate. Northern blotting is a possibility to detect mRNA knockdown, however, similar to Western blotting it is not useful for tissue specific knockdown, and many labs do not have the facilities to work with radioactive labeled probes. Real-time quantitative PCR (qPCR) can be used as an antibody- and probe-independent method to quantify gene knockdown efficiency, and has the advantage that is more sensitive than both Western and Northern blotting.
In the Drosophila community, considerable effort has been made in the development of a qPCR primer design resource (FlyPrimerBank) that considers the optimal primer sequence, as well as the location of primer amplicons with respect to RNAi reagents . However, there is no standard protocol for quantifying RNAi knockdown, and most RNAi reagents have yet to be tested.
Here, we identify optimal qPCR conditions for testing gene knockdown efficiency in combination with the UAS-Gal4 driven Drosophila RNAi system. Specifically, we identify important factors relating to RNA isolation method and primer amplicon position relative to the RNAi target site on the mRNA. Previous studies performed on cultured cells have reported discrepancies between Western blotting results, which detect protein levels, and qPCR measurement, which detects mRNA levels. In these cases, there is a clear loss of protein that is not reflected by quantification of mRNA [9, 10]. Holmes et al. tested mRNA knockdown detection using several primer sets spanning a single gene and found that 3′ primer sets were not able to detect siRNA mediated knockdown as well as 5′ primer sets . They suggested that this was due to the persistence of the 3′ mRNA cleavage fragment in the single gene that was tested . Detection of the 3′ cleavage fragment would result in overestimation of the amount of functional mRNA present in the cell. We expand upon the study by Holmes et al. by testing several loci using primer sets that amplify the 5′ fragment of cleaved mRNA, in combination with cDNA synthesized from purified polyadenylated mRNA. In principle, only uncleaved functional mRNA transcripts would be detected under these conditions (Fig. 1). By testing a combination of different primer locations and RNA isolation methods at multiple loci, we were able confirm the untested theory first proposed by Holmes et al. , and provide optimal conditions for testing RNAi mediated knockdown using transgenic Drosophila RNAi lines.
Expression of dsRNA using the UAS-Gal4 system activates endogenous cellular RNAi machinery. This leads to cleavage of target mRNA molecules at a specific site, producing 5′ and 3′ mRNA fragments that are eventually degraded (Fig. 1). Only uncleaved mRNA transcripts contribute to protein production; however, mRNA cleavage products that have yet to be degraded are likely represented in cDNA libraries that are used as a template for qPCR. Detection of these fragments would result in overestimation of the amount of functional mRNA present in the cell. The use of primer sets that amplify the 5′ fragment of cleaved mRNA in combination with cDNA synthesized from purified polyadenylated mRNA should circumvent this problem, allowing only for the detection of uncleaved mRNA transcripts (Fig. 1). To test this idea, we performed qPCR using two different primer sets—one that amplified a region on the mRNA transcript 5′ of the siRNA cut site, and the other that amplified from the 3′ mRNA cleavage fragment (Figs. 1, 2). We performed qPCR on cDNA synthesized from total RNA, which contains mRNA transcripts, non-coding RNA, and siRNA-mediated cleavage products, and compared this to qPCR results using cDNA derived from mRNA purified using poly-T beads, which excludes 5′ mRNA cleavage fragments (Fig. 1).
Fly stocks, housing, larva collection
Drosophila UAS-RNAi strains were obtained from the Bloomington Drosophila Stock Center (trr-FBgn0023518, UAS-trr 29563 and UAS-trr 36916; osa-FBgn0261885, UAS-osa 38285; brm-FBgn0000212, UAS-brm 37720 and UAS-brm 31712) and the Vienna Drosophila Resource Center (VDRC) (snr1-FBgn0011715, UAS-snr1 12644). All experiments were performed on standard fly food at 25 °C in 70% humidity on a 12-h light/dark cycle. RNAi stocks were crossed to an in-house UAS-Act-Gal4/CyO-ActGFP driver line, created by balancing y; Act-Gal4/CyO (Bloomington 25374) to a CyO-ActGFP balancer (Bloomington 4533). As a control, this driver line was also crossed to a UAS-mCherry-RNAi line (Bloomington 35785). Third instar larvae containing Act-Gal4 and the UAS-RNAi were selected using a Nightsea fluorescence adaptor (Nightsea, Cat. No SFA-LFS-RB) based on the absence of GFP. Ten larvae of each genotype were pooled into a single tube and flash frozen in liquid nitrogen, followed by storage at − 80 °C. Three biological replicates of each genotype were collected.
Primer design and selection
Two primer sets were designed for each targeted gene: one to amplify a region 5′ of the siRNA cut site on the target mRNA, and the other to amplify a 3′ region (Fig. 2, Additional file 1: Table S1). Primers were selected from FlyPrimerBank , if available, or designed using Primer3 v.4.1.0 [11, 12]. Primers were ordered commercially and validated for efficiency using a cDNA dilution series, according the equation “efficiency = 10[− 1/slope]” .
Total RNA isolation, mRNA isolation, and cDNA synthesis
Larvae were disrupted in QIAzol lysis reagent (Qiagen, Cat. No. 79306) using a pestle and further homogenized using QIAshredder tissue homogenization columns (Qiagen, Cat. No. 79654). Total RNA was isolated using the Qiagen RNeasy Lipid Tissue Mini Kit (Qiagen, Cat. No. 74804) according to the manufacturer’s protocol with on-column DNase digestion (Qiagen, Cat. No. 79254). cDNA from total RNA was synthesized using the SensiFast cDNA Synthesis Kit (Bioline, Cat. No. BIO-65053). Remaining total RNA was used for mRNA isolation using the Oligotex mRNA Mini Kit (Qiagen, Cat. No. 70022) scaled to < 0.25 mg. Samples were eluted twice in 20 μL of Buffer OEB, and cDNA was synthesized using the SensiFast cDNA Synthesis Kit.
Real-time quantitative polymerase chain reaction
Quantitative PCR assays were conducted using the SensiFAST SYBR No-Rox kit (Bioline, Cat. No. BIO-98020). Reaction mixtures were composed according to the manufacturer’s protocol. Reactions were carried out in a Bio-Rad CFX384 Real-Time System using the following cycling conditions: 95 °C for 2 min, followed by 40 cycles at 95 °C for 5 s and 65 °C for 30 s. Three qPCR technical replicates were conducted for each biological replicate. Relative expression was normalized to two reference genes, eIF2Bγ (FBgn0034029) and βCOP (FBgn0008635).
Gene expression was normalized to two reference genes (eIF2Bγ and βCOP) using the ∆∆Ct method, correcting for primer efficiencies according to the method described by Pfaffl . One-tailed t tests were performed using Microsoft Excel (version 15.37) to determine if there was a significant reduction in mRNA level detected for each condition, compared to the UAS-mCherry-RNAi control. To compare differences in mRNA levels between primer set locations and RNA isolation methods, GraphPad Prism 7 was used to perform two-way ANOVA on gene expression values normalized to the UAS-mCherry-RNAi control, with false discovery rate corrections for multiple comparisons.
We expressed Drosophila UAS-RNAi transgenes (UAS-snr1 12644, UAS-brm 31712, UAS-brm 37720, UAS-trr 29563, UAS-trr 36916 and UAS-osa 38285) using the ubiquitous Actin-Gal4 driver line. For all of these lines, Act-Gal4 mediated expression resulted in lethality at the pupal stage, consistent with the function of these genes, which are known to be essential for normal development . This suggests that these RNAi lines induce a strong enough knockdown to produce biological effects, and is congruous with previous studies that suggest the effectiveness of these lines in inducing gene knockdown [15,16,17,18,19].
Despite the lethal phenotype upon expression of UAS-snr1 12644 with Act-Gal4, significant knockdown was only detected when using the 5′ snr1 primer set in combination with purified polyadenylated mRNA (Fig. 3a). For all other conditions—3′ primers, and 5′ primer with total RNA—detected mRNA levels were not significantly different from the control. Furthermore, the 5′ primer set detected a significantly greater knockdown on cDNA derived from purified mRNA (relative expression = 53.4%), as compared to total RNA (70.1%). In keeping with our model (Fig. 1), this result suggests that RNAi cleavage fragments may be masking the true RNAi knockdown efficiency, which could have a serious impact on interpretation of data and planning of experiments.
For both UAS-brm 37720 and UAS-brm 31712, a significant knockdown was detected for nearly all conditions tested when compared to controls (Fig. 3b). However, there were also significant differences in the level of knockdown between the different conditions. 5′ primer sets consistently detected a greater brm knockdown than 3′ primer sets. For both brm RNAi lines tested the strongest knockdown was observed using purified polyadenylated mRNA in combination with a 5′ primer set (56.8% for UAS-brm 37720 and 42.4% for UAS-brm 31712), while the weakest knockdown was observed when using 3′ primer sets in combination with total RNA (90.1% for UAS-brm 37720 and 73.2% for UAS-brm 31712). Again, these results suggest that RNAi cleavage fragments may be masking knockdown under suboptimal testing conditions.
For UAS-trr 29563 and UAS-trr 36916, we again observed the strongest detectable knockdown when using 5′ primers in combination with purified mRNA (Fig. 3c). This was especially striking for UAS-trr 36916, which showed almost threefold greater knockdown when using 5′ primers and mRNA (relative expression = 20.0%) compared to when total RNA was used (72.4 and 64.7% using 5′ and 3′ primer sets, respectively).
Finally, for one RNAi line tested, UAS-osa 38285, there was no significant difference in knockdown detection between the different primer set locations and RNA templates (Fig. 3d), suggesting that these factors are not always critical for optimal knockdown detection.
The data presented in this study suggests that detecting RNAi-mediated knockdown using qPCR is dependent on amplicon location, as well as RNA isolation technique. Our results indicate that the best practice for validating RNA-mediated gene knockdown using qPCR is to use cDNA synthesized from purified polyadenylated mRNA, in combination with a primer set that amplifies from the non-polyadenylated 5′ mRNA cleavage product. This approach ensures that mRNA cleavage products will not be detected by qPCR, providing an accurate estimate of functional mRNA levels.
The degree to which non-functional mRNA cleavage fragments impact knockdown detection is quite variable between different RNAi lines tested. Having only tested six RNAi lines, we cannot predict how universal this phenomenon is, and the mechanisms causing this variation are not clear.
cDNA synthesis, a necessary step before performing qPCR, can also result in potential bias and non-uniform representation of mRNA [22, 23]. In this study we attempted to overcome this limitation by using the Sensifast cDNA synthesis kit, which blends random hexamer primers and anchored oligo dT (deoxythymine) to result in unbiased 3′ and 5′ coverage and reverse transcription of all gene regions.
double stranded RNA
real-time quantitative polymerase chain reaction
short interfering RNA
upstream activation sequence
Vienna Drosophila Resource Centre
Hannon GJ. RNA interference. Nature. 2002;418:244–51.
Hashimoto Y, Takahashi M, Sakota E, Nakamura Y. Nonstop-mRNA decay machinery is involved in the clearance of mRNA 5′-fragments produced by RNAi and NMD in Drosophila melanogaster cells. Biochem Biophys Res Commun. 2017;484:1–7.
Orban TI, Izaurralde E. Decay of mRNAs targeted by RISC requires XRN1, the Ski complex, and the exosome. RNA. 2005;11(4):459–69.
Brand AH, Perrimon N. Targeted gene expression as a means of altering cell fates and generating dominant phenotypes. Development. 1993;118(2):401–15.
Dietzl G, Chen D, Schnorrer F, Su KC, Barinova Y, Fellner M, Gasser B, Kisey K, Oppel S, Scheiblauer S, et al. A genome-wide transgenic RNAi library for conditional gene inactivation in Drosophila. Nature. 2007;448:151–6.
Ni JQ, Zhou R, Czech B, Liu LP, Holderbaum L, Yang-Zhou D, Shim HS, Tao R, Handler D, Karpowicz P, et al. A genome-scale shRNA resource for transgenic RNAi in Drosophila. Nat Methods. 2011;8:405–7.
Boutros M, Ahringer J. The art and design of genetic screens: RNA interference. Nat Rev Genet. 2008;9:554–6.
Hu Y, Sopko R, Foos M, Kelley C, Flockhart I, Ammeux N, Wang X, Perkins L, Perrimon N, Mohr SE. Fly Primer Bank: an online database for Drosophila melanogaster gene expression analysis and knockdown evaluation of RNAi reagents. G3 (Bethesda). 2013;3(9):1607–16.
Holmes K, Williams CM, Chapman EA, Cross MJ. Detection of siRNA induced mRNA silencing by RT-qPCR: considerations for experimental design. BMC Res Notes. 2010;3(3):53.
Maier T, Guell M, Serrano L. Correlation of mRNA and protein in complex biological samples. FEBBS Lett. 2009;583:3966–73.
Koressaar T, Remm M. Enhancements and modifications of primer design program Primer3. Bioinformatics. 2007;23(10):1289–91.
Untergrasser A, Cutcutache I, Koressaar T, Ye J, Faircloth BC, Remm M, Rozen SG. Primer3—new capabilities and interfaces. Nucleic Acids Res. 2012;40(15):e115.
Pfaffl MW. A new mathematical model for relative quantification in real-time RT-qPCR. Nucleic Acids Res. 2001;29:e45.
Eissenberg JC, Shilatifard A. Histone H3 Lysine 4 (H3K4) methylation in development and differentiation. Dev Biol. 2010;339:240–9.
Shi J, Zheug M, Ye Y, Li M, Chen X, Hu X, Sun J, Zhang X, Jland C. Drosophila Brahma complex remodels nucleosome organizations in multiple aspects. Nucleic Acids Res. 2014;42:9730–9.
Waldholm J, Wang Z, Brodin D, Tyagi A, Yu S, Theopold U, Farrants AK, Visa N. SWI/SNF regulates the alternative processing of a specific subset of pre-mRNAs in Drosophila melanogaster. BMC Mol Biol. 2011;12:46.
Herz HM, Mohan M, Garruss AS, Liang K, Takahashia YH, Mickey K, Voets O, Verrijzer CP, Shilatifard A. Enhancer-associated H3K4 monomethylation by Trithorax-related, the Drosophila homolog of mammalian Mll3/Mll4. Genes Dev. 2012;26:2604–20.
Oh H, Slattery M, Ma L, White KP, Mann RS, Irvine KD. Yorkie promotes transcription by recruiting a Histrone methyltransferase complex. Cell Rep. 2014;8(2):449–59.
Xie G, Chen H, Jia D, Shu Z, Palmer WH, Huang YC, Zeng X, Hou SX, Jiao R, Deng WM. The SWI/SNF complex protein Snr1 Is a tumor suppressor in Drosophila imaginal tissues. Cancer Res. 2017;77(4):862–73.
Anderson L, Seilhamer J. A comparison of selected mRNA and protein abundances in human liver. Electrophoresis. 1997;18:533–7.
Greenbaum D, Colangelo C, Williams K, Gerstein M. Comparing protein abundance and mRNA expression levels on a genomic scale. Genome Biol. 2003;4(9):117.
Jeffreys AJ, Wilson V, Neumann R, Keyte J. Amplification of human minisatellites by the polymerase chain reaction: towards DNA fingerprinting of single cells. Nucleic Acids Res. 1988;16(23):10953–71.
Korneev S, Blackshaw S, Davies JA. cDNA libraries from a few neural cells. Progress Neurobiol. 1994;42:339–46.
RLM and TAL assisted with protocol design, maintained the Drosophila lines, made the necessary crosses, collected larvae samples, isolated total and mRNA, conducted qPCR, performed data analysis and statistics, designed figures, and were major contributors to writing the manuscript. MMR designed primers, assisted with protocol design and consulted on the manuscript writing. JMK assisted with protocol design, and was a major contributor to writing of the manuscript. All authors read and approved the final manuscript.
We would like to thank the Transgenic RNAi project, the Bloomington Drosophila Stock Centre, and the Vienna Drosophila Resource Centre for providing Drosophila stocks used in this study.
The authors declare that they have no competing interests.
Availability of data and materials
The datasets generated and/or analyzed during the current study are available from the corresponding author on reasonable request.
Consent for publication
Ethics approval and consent to participate
This work was supported by the Natural Science and Engineering Research Council of Canada (NSERC) and the Canada Research Chairs program. RLM was supported by an NSERC undergraduate student research award.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.