PCR-based expression analysis and identification of microRNAs
|David P Lu1, Rebecca L Read1, David T Humphreys2, Fiona M Battah1, David I K Martin2 and John
E J Rasko1,3*
1Gene and Stem Cell Therapy, Centenary Institute of Cancer Medicine and Cell Biology, University of Sydney, Locked Bag No 6, Newtown 2042, Australia. 2Victor Chang Cardiac Research Institute, Sydney, Australia. 3Sydney Cancer Centre, Royal Prince Alfred Hospital, Sydney, Australia
|*Correspondence to: John EJ Rasko, Email: [email protected], Tel: +612 95656156, Fax: +612 95656101|
|Journal of RNAi and Gene Silencing (2005), 1(1), 44-49|
|© Copyright David P Lu et al|
|(Received 19 May 2005; Revised 16 June 2005; Accepted 16 June 2005; Available online 28 July 2005; Published 12 August 2005)|
microRNAs (miRNAs) are small RNAs that regulate translation and hence control a variety of cellular processes in metazoans. The quantitation and identification of miRNAs has been hampered by their small size and low abundance. Here we describe two robust PCR-based assays of miRNA expression based on the original cloning strategy. The non-quantitative PCR method allows detection and identification of miRNAs and we utilise this method in the discovery of a new miRNA (miR-532) in retinoic acid differentiated P19 cells. The second and quantitative method (QM-RT-PCR) is simple and accurate, and uses commonly available technology. Of particular interest is the specificity of this PCR-based technology compared to hybridisation-based methods including arrays and northern blotting. Here we have shown that a single base pair mismatch in the priming sequence results in a two order of magnitude reduction in the amplification of let-7f. These streamlined methods will complement previously described methods and will facilitate analysis of miRNA expression in rare cell populations where the amount of RNA is limited.
|microRNA, miR-532, RT-PCR, TaqMan, differentiation, hairpin RNA, let-7|
|miRNAs are endogenously encoded and produced small (~22 nucleotides) non-coding RNA molecules that have been implicated in the regulation of cellular differentiation, proliferation, apoptosis and metabolism. They are thought to mediate post-transcriptional gene regulation through interaction with target sequences in the mRNA. The result of this interaction may be translational repression (with partially mismatched targets in the 3 untranslated region) or mRNA cleavage (with fully complementary targets in mRNA) (Kawasaki and Tiara, 2003; Moss et al, 1997; Olsen and Ambros, 1999; Wightman et al, 1993). The very broad phylogenetic distribution of miRNAs suggests that they have important functions, but most of these are yet to be discovered (Grad et al, 2003; Houbaviy et al, 2003; Lagos-Quintana et al, 2001; Lau et al, 2001; Lee and Ambros, 2001; Suh et al, 2004. Analysis of miRNA expression levels has been limited by a lack of sensitive and specific detection methods. To date, detection has typically relied on cloning, northern blotting or solid phase array analysis, which require that large quantities of RNA be available (Kawasaki and Tiara, 2003; Houbaviy et al, 2003; Calin et al, 2002; Chen et al, 2004; Krichevsky et al, 2003; Lee et al, 2003; Thomson et al, 2004). Northern blotting is currently the gold standard for miRNA detection (Ambros et al, 2003) but is also limited by the difficulties of achieving stringent hybridization conditions with these very small molecules. Recently, an elaborate modified Invader assay (Third Wave Technologies, Inc) (Allawi et al, 2004) and an RNA-primed, array-based Klenow enzyme assay (RAKE) (Nelson et al, 2004) were reported to quantify miRNA levels. In addition, a quantitative RTPCR method for the amplification of pre-miRNAs has been developed (Schmittgen et al, 2004). Additional methods for miRNA identification and quantification will provide flexibility and facilitate research in this emerging field. We sought simple and robust methods to detect, quantify, and distinguish miRNAs, particularly when cellular RNA is limited.|
Materials and Methods
|Cell culture, differentiation and RNA extraction|
|HeLa cells were cultured in DMEM supplemented with 10% (v/v) FCS (Sigma, St Louis, Missouri) and P19 mouse embryonic carcinoma cells were cultured in DMEM supplemented with 5% (v/v) FCS, at 37°C in a humidified atmosphere containing 5% (v/v) CO2. The human leukemic cell line HL-60 was obtained from the American Type Culture Collection and maintained in RPMI-1640 medium containing L-glutamine supplemented with 10% (v/v) fetal bovine serum (Gibco, Invitrogen, California, USA) and penicillin/streptomycin (1000U, Gibco) at 37oC in a humidified atmosphere containing 5% (v/v) CO2.|
|HL-60 differentiation experiments were performed using 3.5 x 106 cells in 10ml of medium. Monocyte-macrophage differentiation was induced by the addition of 12-Otetradecanoylphorbol- 13-acetate (TPA; 10-7 M, Sigma, St Louis, MO, USA) to the cultures for up to 5 days. Cells were characterised before and after differentiation induction by flow cytometry, and morphologically with May-Giemsa- Grunwald staining of cells centrifuged onto slides by cytospin (Shandon, Pittsburgh, PA, USA; 350 rpm, 10 min).|
|P19 differentiation experiments were performed in bacterial culture dishes to prevent cell adhesion and encourage the formation of embryoid bodies. Cells were plated in DMEM supplemented with 5% (v/v) FCS and retinoic acid (0.5 μM) for two lots of 48 hr (plates were changed at 48 hr to remove adherent cells). Total RNA was extracted using Trizol according to the manufacturer’s instructions (Gibco Life Technologies), suspended in DEPC-treated water and stored at -70°C. For QM-RT-PCR of HL-60 cells in Figure 3B, RNA was isolated using the RNeasy Mini Kit (a procedure that selectively excludes small RNAs of < 200 nucleotides (RNeasy Mini Handbook, Qiagen)).|
|Adaptors and primers|
|The primers used in this work are listed in Table 1. To amplify miRNAs, a 3 RNA-DNA adaptor was ligated to total RNA (Figure 1). Adaptor-50 or adaptor-64 was used for miRNA-RT-PCR and adaptor-Q was used for QM-RTPCR. To minimize non-specific PCR amplification, the adaptors were designed to be non-homologous to any known genes, as determined by Blast homology searches (www.ncbi.nlm.nih.gov/BLAST). To amplify miRNAs or mRNAs, 5 sequence specific primers were used in conjunction with 3 primers complementary to the adaptor sequence. The expression of L-actin mRNA, detected using gene-specific forward and reverse primers, was used as a positive control. The primers used for sequencing were T7 or pGEM-Trevseq.|
Results and Discussion
|The first method (Figure 1) is a non-quantitative but extremely sensitive reverse transcription polymerase chain reaction (RT-PCR) assay. This miRNA-RT-PCR first requires ligation of a DNA/RNA hybrid oligonucleotide adaptor to total RNA, as originally described for cloning of small RNAs (Elbashir et al, 2001; Lau et al 2001). The subsequent PCR reaction uses a 5 primer that is complementary to ~15 nucleotides at the 5 end of the reverse transcribed target miRNA. The 3 primer is complementary to the adaptor (Figure 1A). This design leaves ~7 ‘unconstrained’ nucleotides (depending on the length of the miRNA) at the 3 end of the specific target miRNA. This portion of the amplified fragment is not determined by either of the primers and so may be used to confirm the identity of the amplified miRNA, by cloning and sequencing the PCR product. This improves on published methods by ligating only one universal adaptor to the 3 end of the miRNA (Grad et al, 2003; Lim et al, 2003; Meister et al, 2004). We tested the ability of miRNA-RT-PCR to detect let-7a in HeLa cells, where its expression has been demonstrated by cloning and northern blotting (Lagos-Quintana et al, 2001). miRNA-RT-PCR produced the predicted 94 base pair fragment which includes the adaptor sequence (Figure 1B). The PCR product was cloned and sequenced; all of eleven sequenced clones were identical to let-7a, although three either lacked the 3 thymidine or contained an additional thymidine at the 3 end (data not shown). The let-7a primer should also amplify let-7b and let-7c, as the three miRNAs differ only in their six 3 bases. let-7b and let-7c were not represented in the sequenced clones, which may indicate that HeLa cells do not express let-7b and let7c, or that they are at least ten fold less abundant than let-7a, consistent with RAKE microarray data (Nelson et al, 2004).|
|We tested the sensitivity of miRNA-RT-PCR by diluting total RNA from HeLa cells and assaying for let-7a. As little as 0.5 ng of total RNA from HeLa cells was sufficient to yield the let-7a product (Figure 1C). This indicates that miRNA-RT-PCR is at least 100-fold more sensitive than a commercial solution hybridization assay [mirVanaTM miRNA Detection Kit, Ambion] and the modified Invader assay (Allawi et al, 2004), both of which are more sensitive than the current ‘gold standard’ northern blotting (Ambros et al, 2003). To further assess sensitivity we used total RNA isolated from 200 HeLa cells and obtained the let-7a product (Figure 1C). This indicates that miRNA-RT-PCR may be useful in detecting miRNA expression in rare cell types, or when total RNA is limiting, such as after laser microdissection.|
|Another application for the miRNA-RT-PCR method is to verify the expression of newly-cloned miRNAs (Figure 2). Newly- cloned sequences, or miRNAs predicted in silico, require evidence that they are actually present in some cell type before they can be classed as a true miRNA (Ambros et al, 2003). We cloned a novel 22 nt miRNA, miR-532 (sequence: CATGCCTTGAGTGTAGGACCGT) from retinoic acid differentiated P19 embryonic carcinoma cells, but could not detect its expression in any tissue by northern blot analysis (data not shown). We assayed total RNA from HeLa cells using miRNA-RT-PCR with a 5 primer specific to this newly identified miRNA, and recovered a PCR product of the correct size (Figure 2A). Sequencing of four clones demonstrated that in each the seven unconstrained nucleotides matched the sequence of the cloned miRNA (Figure 2B), thereby confirming expression of this novel miRNA. The genomic sequence from which this miRNA sequence is derived is predicted to fold into a typical miRNA precursor hairpin structure in both the mouse and human genomes (data not shown). The miRNA-532 sequence is highly conserved in mammals (data not shown) and is located ~400nt upstream of miR- 188. The proximity of miR-532 to miR-188 suggests they are processed from a common pri-miRNA transcript. These data demonstrate that miRNA-RT-PCR can be used to validate novel miRNAs and will be of particular importance as a means of experimentally confirming computationally predicted miRNAs (Xie et al, 2005).|
|Although it is extremely sensitive, miRNA-RT-PCR does not accurately measure or compare miRNA expression levels. To address this limitation, we incorporated a single TaqMan probe and real time PCR technology to create quantitative miRNA (QM)-RT-PCR. The TaqMan probe is derived from a portion of the adaptor and is complementary to the adaptor sequence (Figure 3A). Another TaqMan quantative RT-PCR based method using a self annealing 3’ primer for miRNA detection will be released commercially, however, this method lacks utility in miRNA identification because all bases in the amplified miRNA are constrained by primer sequences (Applied Biosystems, Asia Pacific Advanced Genomics Symposium, 2005). As such sequencing of the product would not provide information on the specific target sequence(s) amplified.|
|We have applied QM-RT-PCR to many different miRNAs in various cell lines; here we demonstrate the key features of the method using let-7f as an informative example. 10-fold serial dilutions of synthetic let-7f were added to a fixed amount of HL-60 cell total RNA depleted of small (<200nt) RNAs. Following ligation and RT, amplification of let-7f was concentration dependent (Figure 3B). A relative measure (2-Ct) of miRNA level was made by incorporating a reference target for all samples (Bustin, 2000), in this case GAPDH mRNA. Expression of miR- 19, a miRNA with little similarity to let-7f, was not detectable in the same experiment (data not shown). To test the specificity of QM-RT-PCR, amplification was performed with an identical reaction mix in which the let- 7f primer was replaced by a primer of the same annealing temperature but with a one base pair mismatch in the 15 5’ priming nucleotides (sequence homologous to let-7a, let-7b and let-7c). All three let-7a/b/c sequences are identical in their 16 5 nucleotides but differ from let-7f by a single base change (Figure 3B box). Use of the mismatched primer reduced amplification by greater than two orders of magnitude, confirming relative specificity under these challenging conditions (Figure 3B). We believe that such a stringent test of specificity has not been previously performed in this way and that the specificity of all arrays including northern blotting and solid phase methods has not been critically tested, as argued elsewhere (Nelson et al, 2004). Nevertheless, QM-RT-PCR would be limited in its ability to distinguish between the very few miRNAs that have 16 nucleotides of 5’ sequence identity, such as let-7a/b/c.|
|To assess the utility of QM-RT-PCR in measuring the temporal regulation of miRNA expression, we examined miRNA miR-21 expression during TPA-induced differentiation to the macrophage lineage in HL-60 cells (Figure 3C) (Kasashima et al, 2004). Kasashima et al. reported a 13-fold increase in miR-21 after 2 days of TPA, measuring abundance by the cumbersome method of counting miR-21 clones. QM-RT-PCR measured an increase in miR-21 expression that exceeded 20-fold after five days (Figure 3D).|
|These two methods provide important additions to the tools available for the study of miRNA. miRNA-RT-PCR is a sensitive miRNA detection method that can be used to confirm the expression and identity of newly cloned miRNAs and uses generally available tools. It is likely this method will find immediate utility in confirming computationally predicted miRNAs. QM-RT-PCR is sensitive, quantitative, and economical: only one universal adaptor and one universal TaqMan probe allows the detection of any miRNA, and the specificity of the reaction for a unique miRNA is determined entirely by the 5 miRNA specific oligonucleotide. This method is therefore well suited to high-throughput strategies for precise quantitation of miRNAs. Clearly, QM-RT-PCR provides a new tool to study temporal patterns of specific miRNA expression in development and cellular differentiation and will complement microarray expression profiling (Nelson et al, 2004; Thomson et al, 2004).|
|We thank members of the Rasko and Martin Laboratories for helpful discussions and Cynthia Ng and Louie Leung for technical assistance. We are grateful to the NHMRC (#183760) and Cure Cancer Australia for supporting this work.|
Statement of Compteting Interests
|The authors declared no competing interests.|
|Allawi HT, Dahlberg JE, Olson, S et al. 2004. Quantification of microRNAs using a modified Invader assay. RNA, 10, 1153-1161.
Ambros V, Bartel B, Bartel DP et al. 2003. A uniform system for miRNA annotation. RNA, 9, 277-279.
Berezikov E, Guryev V, van de Belt J, Wienholds E, Plasterk RH and Cuppen E. 2005.Phylogenetic shadowing and computational identification of human microRNA genes. Cell, 120, 21-24.
Bustin SA. 2000. Absolute quantification of mRNA using realtimereverse transcription polymerase chain reaction assays. J MolEndocrinol, 25, 169-193.
Calin GA, Dumitru CD, Shimizu M, et al. 2002.Frequent deletions and down-regulation of micro-RNA genes miR15 and miR16 at 13q14 in chronic lymphocytic leukemia.ProcNatlAcadSci USA, 99, 15524-15529.
Chen CZ, Li L, Lodish HF and Bartel DP. 2004. MicroRNAs modulate hematopoietic lineage differentiation. Science, 303, 83-86.
Elbashir SM, Lendeckel W and Tuschl T. 2001. RNA interference is mediated by 21- and 22-nucleotide RNAs. Genes Dev, 15, 188-200.
Grad Y, Aach J, Hayes GD et al. 2003.Computational and experimental identification of C. elegans microRNAs.MolCell, 11, 1253-1263.
Houbaviy HB, Murray MF and Sharp PA. 2003. Embryonic stem cell-specific MicroRNAs. Dev Cell, 5, 351-358.
Kasashima K, Nakamura Y and Kozu T. 2004. Altered expression profiles of microRNAs during TPA-induced differentiation of HL-60 cells. BiochemBiophys Res Commun, 322, 403-410.
Kawasaki H and TairaK. 2003.Functional analysis of microRNAs during the retinoic acid-induced neuronal differentiation of human NT2 cells. Nucleic Acids Res Suppl, 243-244. (suppl volume)
Krichevsky AM, King KS, Donahue CP, Khrapko K and KosikKS. 2003. A microRNA array reveals extensive regulation of microRNAs during brain development. RNA, 9,1274-1281.
Lagos-Quintana M, Rauhut R, Lendeckel W and Tuschl T. 2001.
Identification of novel genes coding for small expressed RNAs. Science, 294, 853-858.
Lau NC, Lim LP, Weinstein EG and Bartel DP. 2001. An abundant class of tiny RNAs with probable regulatory roles in Caenorhabditiselegans. Science 294:858-862.
Lee RC and Ambros V. 2001.An extensive class of small RNAs in Caenorhabditiselegans. Science, 294, 862-864.
Lee Y, Ahn C, Han J et al. 2003. The nuclear RNase III Droshainitiates microRNA processing. Nature, 425, 415-419.
Lim LP, Glasner ME, Yekta S, Burge CB and Bartel DP. 2003. Vertebrate microRNAs. Science, 299,1540.
Meister G, Landthaler M, Patkaniowska A, Dorsett Y, Teng G and Tuschl T. 2004. Human Argonaute2 mediates RNA cleavage targeted by miRNAs and siRNAs. Mol Cell, 15, 185-197.
Moss EG, Lee RC and Ambros V. 1997. The cold shock domain protein LIN-28 controls developmental timing in C. elegansand is regulated by the lin-4 RNA. Cell, 88, 637-646.
Nelson PT, Baldwin DA, Scearce LM, Oberholtzer CJ, Tobias JW and Mourelatos Z. 2004.Microarray-based, highthroughputgene expression profiling of microRNAs. Nature Methods, 1, 1-7.
Olsen PH, Ambros V. 1999. The lin-4 regulatory RNA controlsdevelopmental timing in Caenorhabditiselegansby blocking LIN-14 protein synthesis after the initiation of translation. DevBiol, 216, 671-680.
Schmittgen TD, Jiang J, Liu Q and Yang L. 2004.A highthroughputmethod to monitor the expression of microRNA precursors. Nucleic Acids Res, 32, e43.
Suh M-R, Lee Y, Kim JY et al. 2004. Human embryonic stem cells express a unique set of microRNAs. DevBiol, 270, 488-498.
Thomson JM, Parker J, Perou CM and Hammond SM. 2004.A custom microarray platform for analysis of microRNA gene expression. Nature Methods, 1, 1-7.
Wightman B, Ha I and Ruvkun G. 1993. Posttranscriptional regulation of the heterochronic gene lin-14 by lin-4 mediates temporal pattern formation in C. elegans. Cell, 75, 855-862.
Xie X, Lu J, Kulbokas EJ, Golub TR, Mootha V, Lindblad-TohK et al. 2005. Systematic discovery of regulatory motifs in human promoters and 3’ UTRs by comparison of several mammals. Nature, 434, 338-345.