The Renilla luciferase gene as a reference gene for normalization of gene expression in transiently transfected cells
© Jiwaji et al; licensee BioMed Central Ltd. 2010
Received: 3 August 2010
Accepted: 31 December 2010
Published: 31 December 2010
The importance of appropriate normalization controls in quantitative real-time polymerase chain reaction (qPCR) experiments has become more apparent as the number of biological studies using this methodology has increased. In developing a system to study gene expression from transiently transfected plasmids, it became clear that normalization using chromosomally encoded genes is not ideal, at it does not take into account the transfection efficiency and the significantly lower expression levels of the plasmids. We have developed and validated a normalization method for qPCR using a co-transfected plasmid.
The best chromosomal gene for normalization in the presence of the transcriptional activators used in this study, cadmium, dexamethasone, forskolin and phorbol-12-myristate 13-acetate was first identified. qPCR data was analyzed using geNorm, Normfinder and BestKeeper. Each software application was found to rank the normalization controls differently with no clear correlation. Including a co-transfected plasmid encoding the Renilla luciferase gene (Rluc) in this analysis showed that its calculated stability was not as good as the optimised chromosomal genes, most likely as a result of the lower expression levels and transfection variability. Finally, we validated these analyses by testing two chromosomal genes (B2M and ActB) and a co-transfected gene (Rluc) under biological conditions. When analyzing co-transfected plasmids, Rluc normalization gave the smallest errors compared to the chromosomal reference genes.
Our data demonstrates that transfected Rluc is the most appropriate normalization reference gene for transient transfection qPCR analysis; it significantly reduces the standard deviation within biological experiments as it takes into account the transfection efficiencies and has easily controllable expression levels. This improves reproducibility, data validity and most importantly, enables accurate interpretation of qPCR data.
The analysis of gene expression using quantitative real-time polymerase chain reaction (qPCR) has become increasingly important as biological research has focused on developing insights into the complex regulatory networks that exist within cells . qPCR is often the assay of choice as it is sensitive and reproducible; it allows the simultaneous analysis of gene expression in a number of different samples and as a result of the high dynamic range, this technique is suitable even when only a few cells are available. The speed of analysis and the potential for automation and multiplexing makes qPCR an attractive technique for the analysis of gene expression [2–4].
Unfortunately, problems attributed to the biological and technical variability which can occur between the different steps of the experimental procedures, are associated with the qPCR assay. The technical variables include the amount of starting materials in the reactions, the quality of the RNA samples and the efficiency of the enzymatic steps (i.e. reverse transcription and PCR) [5, 6]. The biological variables include the differences in the levels of transcriptional expression of genes between tissues and cell types . To take into account these variations, internal reference genes are often used to normalize the qPCR data [8, 9]. Ideally, the internal reference gene should be expressed at levels comparable to the gene of interest and the levels of expression of the gene selected as the internal reference should not vary between the samples and treatments selected for analysis [8, 10]. The selection of the most appropriate internal reference gene serves to decrease the error both within the experiment and between biological experiments [8, 9, 11]. In addition, it allows valid analyses of qPCR data to be conducted . This makes the selection of the internal reference an important factor in the design of a qPCR experiment. This becomes particularly important when the sequence targeted for analysis is being transiently transfected into cells; however validated methods for this type of experiment are not currently available.
In a living cell, it is unlikely that the transcription of any gene is resistant to changes in the cell cycle or in the levels of nutrients. It is therefore important that in the selection of the reference gene, the candidate genes should be regulated at a minimal level. A number of studies have shown that the classical internal reference genes, such as the glyceraldehyde-3-phosphate dehydrogenase gene (GAPDH), are not always the most suitable reference genes , and that the levels of GAPDH mRNA fluctuate in the cell. This is understandable considering the many pathways in which this protein is involved, including endocytosis, translational control, export of nuclear tRNA, DNA replication and repair, apoptosis and glycolysis . GAPDH was originally selected as a normalization reference as the gene encoded a protein with a 'housekeeping' function. It was not until later that the role of GAPDH in the cell was more fully understood, and thus its potential unsuitability as an internal reference. It is, therefore, important that the design of a biological study includes the evaluation of potential internal reference genes, and that the most appropriate reference genes are selected. The availability of software applications such as geNorm, Normfinder and BestKeeper, that use statistical methods to select the most appropriate internal reference genes, make this task easier [13–15].
Most of the traditional internal reference genes are chromosomal genes. The use of a chromosomal internal reference takes into account all the technical and biological variables that are present within the experiment bar one. None of the internal reference genes, as long as the gene was present on the chromosomal DNA, would take into account the variation in the transfection efficiency between samples that had been transiently transfected with plasmid DNA encoding a gene of interest. Also, choosing a genomically encoded gene with an appropriate expression level is challenging. We have adapted a technique often applied to the normalization of reporter enzyme activity to qPCR, and have demonstrated that using a co-transfected plasmid-encoded Renilla reniformis luciferase gene (Rluc) as an internal reference is more appropriate for the normalization of qPCR data in transiently transfected systems.
Results and Discussion
We are developing a plasmid vector to study transcriptional activation derived from cis-acting regulatory elements of interest in transiently transfected mammalian cell lines. We planned to validate the putative cis-acting regulatory elements using a reporter enzyme system that could be adapted to allow a high-throughput analysis of the samples of interest. Traditional reporter enzyme assays usually use the activity of one reporter enzyme to normalize the activity of another enzyme of interest when two plasmids, each separately encoding one of the enzymes, have been transiently transfected into cells. One of the most common systems is the firefly luciferase-Renilla luciferase enzyme system, where the enzyme activity of Renilla luciferase (encoded by the control plasmid) is used to normalize the firefly luciferase enzyme activity (present on the experimental plasmid).
A complete workflow for the firefly luciferase-Renilla luciferase system from the plasmid constructs used to transfect the cells to the reporter enzyme assay system used to quantify the reporter enzyme activity has been commercially developed (Promega). This uses one of three control plasmid constructs with different promoters upstream of the Rluc gene . pRL-SV40, which encodes the early simian virus 40 (SV40) enhancer/promoter region, and pRL-CMV, which encodes the cytomegalovirus (CMV) immediate early enhancer/promoter region both provide high levels of Renilla luciferase enzyme expression. In contrast, pRL-TK encodes the weaker herpes simplex virus (HSV) thymidine kinase (TK) promoter, resulting in lower levels of Renilla luciferase enzyme expression.
For the studies described here the most appropriate system would have a promoter upstream of Renilla luciferase with relatively high expression levels such that a low concentration of the control plasmid could be used in the transient transfections. It would also be advantageous to use a promoter that did not encode transcription factor binding sequences. All three control plasmids provided as part of the commercial kit have their limitations. Breuning et al.  and Svensson et al.  reported that the CMV promoter was susceptible to up-regulation in mammalian cells under conditions of cell stress. This is not ideal since in cis-regulatory element studies it is likely that some of the treatments used to activate the cis-regulatory elements will result in cell stress. A number of reports have also highlighted that HSV TK promoter is susceptible to altered regulation under a range of experimental conditions [19–21]. Ho et al.  performed a comparison between pRL-TK and pRL-SV40 and stated that of the two plasmids, pRL-SV40 was less susceptible to non-specific activation by the GATA transcription factor than pRL-TK. As this project aims to study cis-acting regulatory elements, we wanted to select the control plasmid that was least affected by undesired transcriptional regulation so we selected pRL-SV40 as the control plasmid.
HEK293 cells were transiently tranfected with pRL-SV40 (50 ng) and an equivalent concentration of our reporter system plasmids pMN85 or pMN86. The concentration of pMN85/pMN86 was determined in an experiment similar to the one that was conducted for pRL-SV40 (data not shown). The control plasmid, pRL-SV40, is 3.7 kb in size compared to pMN85/pMN86, both of which are 5.1 kb in size. Yin et al.  reported that there was an inverse correlation between the plasmid size and the transfection efficiency of that plasmid. As both the reporter system plasmids are of comparable size (pMN86 is 82 bp larger than pMN85), we determined that the transfection efficiencies between these plasmids would not be significantly different. In addition, despite the differences in the transfection efficiencies of the reporter system plasmids and the control plasmid, the transfection efficiencies of the respective plasmids should be comparable between samples allowing us to normalize the data derived from the reporter system plasmids to that of pRL-SV40, the control plasmid.
Cd-dependent induction of MRE-directed expression of the firefly luciferase reporter enzyme activity
Normalized Fold Induction
1.00 ± 0.00
212.77 ± 8.92
In order to quantify the transcriptional activity of the induction of the Cd responsive MRE sequence using qPCR, we needed to find an appropriate normalization reference for these experiments. We planned to test the induction of cis-acting regulatory elements in the presence of a variety of compounds including metals (i.e. Cd), steroids (i.e. dexamethasone), cAMP activators (i.e. forskolin) and activators of protein kinase C (i.e. phorbol-12-myristate 13-acetate (TPA)). Hence, the reference gene selected as a normalization control had to be stably expressed and unaffected by these inducers.
Traditionally, a chromosomal gene is selected as a reference gene for normalization. To this end, we tested a panel of 11 chromosomal genes as internal controls for qPCR normalization using polymerase chain reaction (PCR) primers and fluorescent probes provided in the geNorm Housekeeping Selection Kit (PrimerDesign Ltd). The panel included the genes for 18S rRNA (18S), beta-actin (ActB), ATP synthase (Atp5B), beta-2-microglobulin (B2M), cytochrome c-1 (Cyc1), isoform 2 of the eukaryotic translation initiation factor 4A (Eif4A2), GAPDH, ribosomal protein L13a (RPL13A), subunit A of the succinate dehydrogenase complex (SDHA), DNA topoisomerase I (TOP1) and the zeta polypeptide of the tyrosine-3-monooxygenase/tryptophan-5-monooxygenase activation protein (YWHAZ). The expression of these 11 chromosomal genes was measured in HEK293 cells that had been co-transfected with pMN85 or pMN86 and pRL-SV40 in the absence and presence of Cd, dexamethasone, forskolin and TPA. The data was compared using three software applications; geNorm, Normfinder and BestKeeper (websites are shown in the Methods section).
Ranking of chromosomal reference genes based on geNorm, Normfinder and BestKeeper analyses
Normfinder is also a Visual Basic Application and it too assigns a stability value to the candidate reference genes. In contrast to geNorm, Normfinder uses a model-based approach to provide a value for the two most stable reference genes with the least intra- and inter-group variation . The robustness of the Normfinder approach for the selection of reference genes for normalization was demonstrated by Andersen et al. . As for the geNorm analysis, the experimental Cq values were converted to a linear scale of expression using the delta Ct method. This data loaded into Normfinder in two ways- (1) with the data separated into groups corresponding to treatment with the different activators, and (2) with all the treatment groups combined. In Normfinder, stability is expressed as a stability value in arbitrary units. Analysis of our data (Table 2) indicated that the RPL13A chromosomal gene was the most stable (stability value = 0.195 ± 3.090) and that again SDHA was the least stable gene (stability value = 31.894 ± 10.087).
BestKeeper is an Excel-based tool that calculates the variability in the expression of the pool of reference genes by analyzing the Cq values directly and defining variability by the standard deviation (SD) and coefficient of variance (CV) . All genes with SD values of greater than 1 are considered to be unsuitable reference genes. The gene with the lowest SD is considered to be the most suitable reference gene. BestKeeper then performs a comparative analysis based on the pair-wise correlations of the all the candidate reference genes to each other and generates a BestKeeper index. The comparison of this BestKeeper index to each candidate reference gene, using pair-wise correlation analyses, results in each candidate reference gene being assigned a value for the Pearson correlation coefficient (r) and the probability (p). The most suitable reference gene is that with an r value closest to 1.0, with the p value providing an indication of the significance of the r value. The most unstable gene in the two previous analyses, SDHA, was not included in the BestKeeper analysis. The BestKeeper analysis (Table 2) of the potential reference genes highlighted B2M, Eif4A2 and RPL13A as being the most suitable normalization reference genes (r = 0.974, p = 0.001) and TOP1 as the least suitable (r = 0.704, p = 0.119).
Comparison of the ranking of B2M and ActB, the selected chromosomal reference genes, and Rluc, present on a transiently co-transfected plasmid, based on geNorm, Normfinder and BestKeeper analyses
The stability values for Rluc were generally lower than for the chromosomal genes (0.63 (geNorm), 1.236 (Normfinder) and 0.694 (BestKeeper), as might be expected, since this reflects both the problem of variable transfection efficiency, and the lower levels of expression, meaning that more rounds of PCR are needed, accentuating any variable transfection efficiencies. However, the scores were still better than some, and close to the majority, of the generally used chromosomal genes, suggesting that this could still be considered as a candidate system for normalization [24–26]. According to the geNorm manual [http://medgen.ugent.be/~jvdesomp/genorm/geNorm_manual.pdf], genes with M≥1.5 are deemed unsuitable as reference genes. Based on this cut-off point, Rluc is not excluded as a suitable reference gene. Unfortunately, neither Andersen et al.  nor Pfaffl et al.  have reported guidelines on the range of acceptable stability values for the output of the Normfinder or BestKeeper applications, respectively [14, 15]. This makes it difficult to determine the suitability of Rluc as a normalization control based on the values of the output of the Normfinder and BestKeeper applications. The output of all three bioinformatic analyses indicated that Rluc was not the most suitable gene for use as a reference control.
Cq values obtained for Fluc, Rluc, B2M and ActB in qPCR reactions
Cq value ± SD
27.12 ± 2.76
27.98 ± 0.44
19.92 ± 0.42
19.08 ± 0.52
Copies of Fluc normalized to the copies of the reference genes Rluc, B2M and ActB
Normalized copies (× 10-3)
Normalized copies (× 10-3)
Induction of Fluc expression in the presence of Cd
Fluc/Rluc ± SD (% SD)
Fluc/B2M ± SD (% SD)
Fluc/ActB ± SD (% SD)
1.34 ± 0.42
1.07 ± 0.27
1.35 ± 0.28
61.42 ± 10.12 (16%)
71.38 ± 16.24 (23%)
74.92 ± 25.80 (35%)
Credible Intervals for Fluc expression in the presence of Cd
95% Credible intervals for fold induction
Mean fold induction
A good control for normalization in qPCR experiments should be stably expressed at levels similar to those of the genes of interest, and should take into account all factors that could affect expression levels in the experiment, including transfection efficiency. While we observed that the chromosomal internal control genes B2M and ActB were more stably expressed than transfected Rluc under the experimental conditions, the levels of expression of B2M and ActB were much higher than those of genes transiently transfected into the cell. The Rluc gene was expressed at levels similar to the gene of interest and was demonstrated by a number of biological analyses to be the most appropriate control, ensuring that differences in the transfection efficiencies between the samples were taken into account, thus significantly improving the reproducibility and the validity of the qPCR experiments.
General Culture Conditions
Escherichia coli DH5α cells  were used for the construction, screening and propagation of plasmid constructs. E. coli DH5α cells were cultured in lysogeny broth (LB) or on LB agar containing 100 μg/mL ampicillin as described by Sezonov et al. . HEK293 cells (LGC Standards, Middlesex, UK; ATCC Number CRL-1573)  were maintained in Dulbecco's Modified Eagle Medium (DMEM, Invitrogen Ltd, Paisley, UK; 41966-029) supplemented with 4 mM L-glutamine (Sigma-Aldrich, Dorset, UK; G7513) and 10% foetal bovine serum (FBS; Biosera, East Sussex, UK; S1810/500) at 37°C in an atmosphere that contained 5% CO2 as per the supplier's propagation instructions.
Chromosomal reference genes
The geNorm Housekeeping Gene Selection Kit (PrimerDesign Ltd, Southampton, UK; ge-PP-12-hu) was used as the panel of chromosomal reference genes. It included forward primers, Fluorescein (FAM)-labelled fluorescent probes and reverse primers for the chromosomal genes for 18S, ActB, Atp5B, B2M, Cyc1, Eif4A2, GAPDH, RPL13A, SDHA, TOP1 and YWHAZ.
Description of plasmids pMN85 and pMN86
Plasmid pMN85 was derived from pGL3 Basic (Promega Corporation, Madison, WI, USA; E1751). DNA encoding (1) a multiple cloning cassette into which transcription factor binding sites could be inserted for analysis and (2) a PTK was inserted between the Kpn I and Hind III restriction enzyme sites upstream of Fluc in pGL3 generating pMN85. In pMN86 (Figure 2), five copies of a nucleotide sequence sensitive to induction in the presence of Cd (MRE; CGG GTG CGC CCG GCC CGA) were inserted between the Eco RI and Eco RV restriction enzyme sites upstream of PTK.
Transfection of HEK293 cells
Plasmid DNA was prepared using the PureYield Plasmid Midiprep System (Promega Corporation, Madison, WI, USA; A2495). 1 × 106 HEK293 cells and, unless otherwise stated, a total of 100 ng DNA (50 ng pRL-SV40 (Promega Corporation, Madison, WI, USA; E2231) and 50 ng of either pMN85 or pMN86) was used in each transfection. HEK293 cells were transfected using Genejuice (Novagen supplied by Merck Chemicals Ltd, Nottingham, UK; 70967-6) as recommended by the manufacturers instructions. Transfected HEK293 cells were incubated for 24 hours and subsequently induced with 12.5 μM CdCl2 (Sigma-Aldrich, Dorset, UK; 439800), 2.5 μM dexamethasone (Sigma-Aldrich, Dorset, UK; D8893), 25 μM forskolin (Sigma-Aldrich, Dorset, UK; F6886) or 0.25 μM TPA (Sigma-Aldrich, Dorset, UK; P1585) for 4 hours. Each experiment contained n = 6 biological replicates and experiments were repeated at least three times. Results between the experiments were comparable.
Reporter enzyme assays
The firefly and Renilla luciferase activities in 20 μL HEK293 cell lysate were quantified with the dual luciferase reporter assay system as per the manufacturer's instructions (Promega Corporation, Madison, WI, USA; E1910) in a Luminoskan Ascent (Thermo Fisher Scientific, Waltham, Ma, USA) luminometer. The firefly and Renilla luciferase enzyme activities were measured for each biological sample. The firefly luciferase enzyme activity was normalized to the Renilla luciferase enzyme activity. The Cd-treated samples were then normalized to untreated samples yielding values for the fold induction of the firefly luciferase enzyme activity directed by the MRE in the presence of Cd.
RNA purification and cDNA synthesis
Total RNA was prepared from 1 × 106 HEK293 cells using the miRNeasy mini kit (Qiagen Ltd, West Sussex, UK; 217004), mRNA was isolated from 4 μg of total RNA with 25 μL of the Dynabeads mRNA Purification Kit (Invitrogen Ltd, Paisley, UK; 610-06) and cDNA produced from the purified mRNA with 100 ng random hexamers (Invitrogen Ltd, Paisley, UK; 48190011) or 500 μg/mL oligo dT(12-18) primer (Invitrogen Ltd, Paisley, UK; 18418-012), 2.5 mM of each dNTP (Promega Corporation, Madison, WI, USA; U1240), 40 U RNaseOUT (Invitrogen Ltd, Paisley, UK; 10777-019), and 100 U Superscript II enzyme (Invitrogen Ltd, Paisley, UK; 18064014) in a final volume of 20 μL. The cDNA synthesis reaction was incubated at 25°C for 10 minutes, 42°C for 60 minutes and finally, at 70°C for 15 minutes.
Nucleotide sequences of the primers used in the qPCR reactions and the length of the amplicons generated
Transcript (amplicon length)
Sequence (5' to 3')
Fluc (137 nts)
Rluc (140 nts)
B2M (125 nts)
ActB (142 nts)
Analysis of the qPCR data for the reference genes
Cq values obtained for HEK293 cells transiently transfected with pMN85/pMN86 and pRL-SV40 were grouped into 5 sets: (1) Untreated samples, (2) + Cd, (3) + dexamethasone, (4) + forskolin and (5) + TPA. Each subset contained n = 6 biological replicates. These sets of data were analyzed with three applications- geNorm http://medgen.ugent.be/~jvdesomp/genorm/, Normfinder http://www.mdl.dk/publicationsnormfinder.htm and BestKeeper http://www.gene-quantification.com/bestkeeper.html#download as described by the documentation for the software applications.
Generation of standard curves for Fluc, Rluc, B2M and ActB
Linear dsDNA was used to generate standard curves for qPCR, after Hou et al.  highlighted the benefits of using linear DNA in preference to circular DNA. The Fluc, Rluc, B2M and ActB PCR products generated in the qPCR reactions above were cloned into the pGEM-T-Easy plasmid using the pGEM-T-Easy Vector System I (Promega Corporation, Madison, WI, USA; A1360). Ligations were transformed into Subcloning Efficiency™ DH5α™ Competent Cells (Invitrogen Ltd, Paisley, UK; 18265-017). Plasmid constructs were screened for the desired insert using the GoTaq® Flexi DNA Polymerase (Promega Corporation, Madison, WI, USA; M8291) in PCR reactions that contained 2 mM MgCl2, 25 ng DNA, 0.2 mM of each dNTP (Promega Corporation, Madison, WI, USA; U1240), 1 μM gene specific forward primer (Table 8), 1 μM gene specific reverse primer (Table 8) and 0.5 U GoTaq DNA polymerase in a final volume of 25 μL. The PCR conditions were: (1) 92°C for 2 minutes, 1 cycle; (2) 92°C for 30 seconds, 52°C for 30 seconds, 72°C for 1 minute, 30 cycles and (3) 72°C for 2 minutes, 1 cycle. The DNA in the PCR reactions was separated on a 1% Seakem LE Agarose gel (Lonza Biologics plc, Slough, Berkshire, UK; 50004) containing 0.01% Sybr Safe DNA gel stain (Invitrogen Ltd, Paisley, UK; S33102). Constructs that showed a positive PCR reaction were selected for validation by DNA sequencing.
The efficiencies and errors of the qPCR standard curves as calculated by the Lightcycler 480 software (Roche Diagnostics Ltd, Burgess Hill, UK)
qPCR data analysis
Cq values obtained for HEK293 cells transiently transfected with pMN85/pMN86 and pRL-SV40 were grouped into 4 sets: (1) pMN85 - Untreated, (2) pMN85 - + Cd, (3) pMN86 - Untreated, (4) pMN86 - + Cd. The Cq values were converted to copy numbers using the Fluc, Rluc, B2M and ActB standard curves. The copy number of Fluc was normalized to the copy number of Rluc, B2M or ActB. The copy number of the normalized Cd-treated samples were then compared to the copy number of the normalized untreated samples for pMN85 and pMN86 yielding values for the fold induction of the Fluc directed by the MRE in the presence of Cd.
Statistical analysis of qPCR data
The data were analysed using JAGS , a Gibbs sampler for hierarchical models and Coda , a tool for examining Markov Chain Monte Carlo runs. The data were modelled using a log-normal distribution, with uninformative normal priors for the mean parameter and uninformative gamma priors for the precision. For each set of data, 10 independent runs were used, each with a burn in of 20000 iterations and 100000 samples taken. From these samples, 95% credible intervals were calculated for the posterior mean and standard deviation of the fold induction (Table 7).
List of Abbreviations
- 18S :
18S rRNA gene
- ActB :
beta actin gene
- Atp5B :
ATP synthase gene
- amp :
β-lactamase gene conferring resistance to ampicillin
- B2M :
Black Hole Quencher 1
coefficient of variance
Dulbecco's Modified Eagle Medium
efficiency of the qPCR reaction
- Eif4A2 :
eukaryotic translation initiation factor 4A, isoform 2 gene
foetal bovine serum
- Fluc :
firefly luciferase gene
- GAPDH :
glyceraldehyde-3-phosphate dehydrogenase gene
Human Embryonic Kidney 293
- HMBS :
hydroxymethylbilane synthase gene
herpes simplex virus
- M value:
geNorm measure of expression stability
metal responsive element
- ori :
origin of replication
polymerase chain reaction
thymidine kinase promoter
quantitative real-time polymerase chain reaction
BestKeeper coefficient of correlation
- Rluc Renilla :
- RPL13A :
ribosomal protein L13a gene
- SDHA :
succinate dehydrogenase complex, subunit A gene
simian virus 40
- TOP1 :
DNA topoisomerase I gene
- YWHAZ :
tyrosine-3-monooxygenase/tryptophan-5-monooxygenase activation protein (zeta polypeptide) gene.
We would like to thank Sarah Cumming for discussions during the preparation of this manuscript and the EPSRC for funding this study through a Basic Technology Grant (EP/E0327545/1).
- Ginzinger DG: Gene quantification using real-time quantitative PCR: an emerging technology hits the mainstream. Experimental Hematology. 2002, 30: 503-512. 10.1016/S0301-472X(02)00806-8View ArticlePubMedGoogle Scholar
- Higuchi R, Fockler C, Dollinger G, Watson R: Kinetic PCR analysis: real-time monitoring of DNA amplification reactions. Biotechnology. 1993, 11: 1026-1030. 10.1038/nbt0993-1026View ArticlePubMedGoogle Scholar
- Heid CA, Stevens J, Lival KJ, Williams PM: Real time quantitative PCR. Genome Research. 1996, 6: 986-994. 10.1101/gr.6.10.986View ArticlePubMedGoogle Scholar
- Nolan T, Hands RE, Bustin SA: Quantification of mRNA using real-time RT-PCR. Nature Protocols. 2006, 1: 1559-1582. 10.1038/nprot.2006.236View ArticlePubMedGoogle Scholar
- Bustin SA: Quantification of mRNA using real-time reverse transcription PCR (RT-PCR): trends and problems. Journal of Molecular Endocrinology. 2002, 29: 23-39. 10.1677/jme.0.0290023View ArticlePubMedGoogle Scholar
- Bustin SA: Absolute quantification of mRNA using real-time reverse transcription polymerase chain reaction. Journal of Molecular Endocrinology. 2000, 25: 169-193. 10.1677/jme.0.0250169View ArticlePubMedGoogle Scholar
- Aerts JL, Gonzalesm MI, Topalian SL: Selection of appropriate reference genes to assess expression of tumor antigens using real-time PCR. Biotechniques. 2004, 36: 84-91.PubMedGoogle Scholar
- Dheda K, Huggett JF, Bustin SA, Johnson MA, Rook G, Zumla A: Validation of housekeeping genes for normalizing RNA expression in real-time PCR. Biotechniques. 2004, 37: 112-119.PubMedGoogle Scholar
- Radonic A, Thulke S, Mackay IM, Landt O, Siegert W, Nitsche A: Guideline to reference gene selection for quantitative real-time PCR. Biochemical and Biophysical Research Communications. 2004, 313: 856-862. 10.1016/j.bbrc.2003.11.177View ArticlePubMedGoogle Scholar
- Dheda K, Huggett JF, Chang JS, Kim LU, Bustin SA, Johnson MA, Rook GAW, Zumla A: The implications of using an inappropriate reference gene for real-time reverse transcription PCR data normalization. Analytical Biochemistry. 2005, 344: 141-143. 10.1016/j.ab.2005.05.022View ArticlePubMedGoogle Scholar
- Huggett J, Dheda K, Bustin S, Zumla A: Real-time RT-PCR normalization; strategies and considerations. Genes and Immunity. 2005, 6: 279-284. 10.1038/sj.gene.6364190View ArticlePubMedGoogle Scholar
- Suzuki T, Higgins PJ, Crawford DR: Control selection for RNA quantitation. Biotechniques. 2000, 29: 332-337.PubMedGoogle Scholar
- Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, De Paepe A, Speleman F: Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal reference genes. Genome Biology. 2002, 3 (7): research0034.1-0034.11. 10.1186/gb-2002-3-7-research0034. 10.1186/gb-2002-3-7-research0034View ArticleGoogle Scholar
- Pfaffl MW, Tichopad A, Prgomet C, Neuvians TP: Determination of stable housekeeping genes, differentially regulated target genes and sample integrity: BestKeeper - Excel-based tool using pair-wise correlations. Biotechnology Letters. 2004, 26: 509-515. 10.1023/B:BILE.0000019559.84305.47View ArticlePubMedGoogle Scholar
- Andersen CL, Jensen JL, Orntoft TF: Normalization of real-time quantitative reverse transcription PCR data: a model based variance estimation approach to identify genes sioted for normalization, applied to bladder and colon cancer data sets. Cancer Research. 2004, 64: 5245-5250. 10.1158/0008-5472.CAN-04-0496View ArticlePubMedGoogle Scholar
- pRL Renilla Luciferase Reporter Vectors. Technical Bulletin No. TB550 Promega Corporation, Madison WI USA. 2008
- Bruening W, Giasson B, Mushynski W, Durham HD: Activation of stress-activated MAP protein kinases up-regulates expression of transgenes driven by the cytomegalovirus immediate/early promoter. Nucleic Acids Research. 1998, 26: 486-489. 10.1093/nar/26.2.486View ArticlePubMedPubMed CentralGoogle Scholar
- Svensson RU, Barnes JM, Rokhlin OW, Cohen MB, Henry MD: Chemotherapeutic Agents Up-regulate the Cytomegalovirus Promoter: Implications for Bioluminescence Imaging of Tumor Response to Therapy. Cancer Res. 2007, 67: 10445-10454. 10.1158/0008-5472.CAN-07-1955View ArticlePubMedGoogle Scholar
- Thavathiru E, Das GM: Activation of pRL-TK by 12S E1A oncoprotein: drawbacks of using an internal reference reporter in transcription assays. Biotechniques. 2001, 31: 528-532.PubMedGoogle Scholar
- Matuszyk J, Ziolo E, Cebrat M, Kockel I, Strzadala L: Nurr1 affects pRL-TK but not phRG-B internal control plasmid in genetic reporter system. Biochemical and Biophysical Research Communications. 2002, 294: 1036-1039. 10.1016/S0006-291X(02)00601-0View ArticlePubMedGoogle Scholar
- Ho CKM, Strauss JF: Activation of the control reporter plasmids pRL-TK and pRL-SV40 by multiple GATA transcription factors can lead to aberrant normalization of transfection efficiency. BMC Biotechnology. 2004, 4: 10-14. 10.1186/1472-6750-4-10View ArticlePubMedPubMed CentralGoogle Scholar
- Yin W, Xiang P, Li Q: Investigations of the effectof DNA size in transient transfection assay using dual luciferase system. Analytical Biochemistry. 2005, 346: 289-294. 10.1016/j.ab.2005.08.029View ArticlePubMedGoogle Scholar
- Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2-ΔΔCt method. Methods. 2001, 25: 402-408. 10.1006/meth.2001.1262View ArticlePubMedGoogle Scholar
- Perez R, Tupac-Yupanqui I, Dunner S: Evaluation of suitable reference genes for gene expression studies in bovine muscular tissue. BMC Molecular Biology. 2008, 9: 79- 10.1186/1471-2199-9-79View ArticlePubMedPubMed CentralGoogle Scholar
- Valente V, Teixeira SA, Neder L, Okamoto OK, Oba-Shinjo SM, Marie SKN, Scrideli CA, Paco-Larson ML, Carlotti CG: Selection of suitable housekeeping genes for expression analysis in glioblastoma using quantitative RT-PCR. BMC Molecular Biology. 2009, 10: 17- 10.1186/1471-2199-10-17View ArticlePubMedPubMed CentralGoogle Scholar
- Lallemant B, Evrard A, Combescure C, Chapuis H, Chambon G, Raynal C, Reynaud C, Sabra O, Joubert D, Hollande F, Lallemant JG, Lumbroso S, Brouillet JP: Reference gene selection for head and neck squamous cell carcinoma gene expression studies. BMC Molecular Biology. 2009, 10: 78- 10.1186/1471-2199-10-78View ArticlePubMedPubMed CentralGoogle Scholar
- Frost P, Nilsen F: Validation of reference genes for transcription profiling in the salmon louse Lepeophtheirus salmonis by quantitative real-time PCR. Vet Parasitology. 2003, 118: 169-174. 10.1016/j.vetpar.2003.09.020. 10.1016/j.vetpar.2003.09.020View ArticleGoogle Scholar
- Hanahan D: Studies on transformation of Escherichia coli with plasmids. Journal of Molecular Biology. 1983, 166: 557-580. 10.1016/S0022-2836(83)80284-8View ArticlePubMedGoogle Scholar
- Sezonov G, Joseleau-Petit D, D'Ari R: Escherichia coli physiology in Luria-Bertani broth. Journal of Bacteriology. 2007, 189: 8746-8749. 10.1128/JB.01368-07View ArticlePubMedPubMed CentralGoogle Scholar
- Graham FL, Smiley J, Russell WC, Nairn R: Characteristics of a human cell line transformed by DNA from human adenovirus Type 5. Journal of General Virology. 1977, 36: 59-72. 10.1099/0022-1317-36-1-59View ArticlePubMedGoogle Scholar
- Hou Y, Zhang H, Miranda L, Lin S: Serious overestimation in quantitative PCR by circular (supercoiled) plasmid standard: microalgal pcna as the model gene. PLoS ONE. 2010, 5: e9545- 10.1371/journal.pone.0009545View ArticlePubMedPubMed CentralGoogle Scholar
- Whelan JA, Russell NB, Whelan MA: A method for using absolute quantification of cDNA using real-time PCR. Journal of Immunological Methods. 2003, 278: 261-269. 10.1016/S0022-1759(03)00223-0View ArticlePubMedGoogle Scholar
- Yun JJ, Heisler LE, Hwang IIL, Wilkins O, Lau SK, Hyrcza M, Jayabalasingham B, Jin J, McLaurin J, Tsao MS, Der SD: Genomic DNA functions as a universal external standard in quantitative real-time PCR. Nucleic Acids Research. 2006, 34: e85- 10.1093/nar/gkl400View ArticlePubMedPubMed CentralGoogle Scholar
- Plummer M: JAGS: A program for analysis of Bayesian graphical models using Gibbs sampling. 2003, Proceedings of the 3rd International Workshop on Distributed Statistical ComputingGoogle Scholar
- Plummer M, Best N, Cowles K, Vines K: CODA: Convergence Diagnosis and Output Analysis for MCMC. R News. 2006, 6: 11-17.Google Scholar