- Research article
- Open Access
- Open Peer Review
This article has Open Peer Review reports available.
Selection of suitable reference genes for accurate normalization of gene expression profile studies in non-small cell lung cancer
© Saviozzi et al; licensee BioMed Central Ltd. 2006
Received: 10 February 2006
Accepted: 26 July 2006
Published: 26 July 2006
In real-time RT quantitative PCR (qPCR) the accuracy of normalized data is highly dependent on the reliability of the reference genes (RGs). Failure to use an appropriate control gene for normalization of qPCR data may result in biased gene expression profiles, as well as low precision, so that only gross changes in expression level are declared statistically significant or patterns of expression are erroneously characterized. Therefore, it is essential to determine whether potential RGs are appropriate for specific experimental purposes. Aim of this study was to identify and validate RGs for use in the differentiation of normal and tumor lung expression profiles.
A meta-analysis of lung cancer transcription profiles generated with the GeneChip technology was used to identify five putative RGs. Their consistency and that of seven commonly used RGs was tested by using Taqman probes on 18 paired normal-tumor lung snap-frozen specimens obtained from non-small-cell lung cancer (NSCLC) patients during primary curative resection.
The 12 RGs displayed showed a wide range of Ct values: except for rRNA18S (mean 9.8), the mean values of all the commercial RGs and ESD ranged from 19 to 26, whereas those of the microarray-selected RGs (BTF-3, YAP1, HIST1H2BC, RPL30) exceeded 26. RG expression stability within sample populations and under the experimental conditions (tumour versus normal lung specimens) was evaluated by: (1) descriptive statistic; (2) equivalence test; (3) GeNorm applet. All these approaches indicated that the most stable RGs were POLR2A, rRNA18S, YAP1 and ESD.
These data suggest that POLR2A, rRNA18S, YAP1 and ESD are the most suitable RGs for gene expression profile studies in NSCLC. Furthermore, they highlight the limitations of commercial RGs and indicate that meta-data analysis of genome-wide transcription profiling studies may identify new RGs.
Lung cancer remains the leading cause of cancer-related death in Europe and North America and accounts for nearly 30% of the total [1, 2]. Despite advances in treatment, the prognosis of non-small cell lung cancer (NSCLC) is poor: only 5–15% of patients survive 5 years after diagnosis, mainly in function of the initial stage of the disease . Real-time quantitative PCR (qPCR) is the foremost method of choice for accurate determination of transcripts amounts. It has recently been applied in lung cancer studies to quantify the expression level of predictive and/or prognostic target genes [4, 5].
In gene expression studies, qPCR data need to be normalized to remove nonspecific variability associated with differences in RNA quantity and quality, usually by relative quantification whereby the expression level of the target gene is determined relatively to another gene transcript, the so-called reference gene (RG), and the results are expressed as a target to reference ratio . The reliability of normalized data is highly dependent on RG robustness. An ideal RG should be constitutively expressed and characterized by stable expression in different samples/experimental conditions.
Genes considered to be valid RGs in semi-quantitative techniques (eg. Nothern Blot) may be less reliable when highly sensitive qPCR is used . Identification of more sensitive RGs and their validation within specific biological conditions are thus critical issues in qPCR studies .
Microarray data can be used to identify potential RGs by modeling the expression data to select those with the most stable expression [9–11]. Most lung cancer qPCR studies use commercial RGs, such as glyceraldehyde-3-phosphate dehydrogenase (GAPDH) , β-actin (ACTB) , TATA-binding protein , β-microglobulin , cytochrome oxidase II , 18S ribosomal RNA (rRNA18S) . Their reliability in this contex, however, has not been demonstrated. Furthermore, it is increasingly evident that in many experimental situations the use of GAPDH or ACTB is inappropriate due to their high variability [10, 17–20]. Thus, the aim of the present study was the selection and validation of new RGs for the differentiation of normal and tumor lung specimens. Five putative RGs (ESD, BTF3, HIST1H2BC, RPL30, YAP1) selected from a meta-analysis of geneChip experiments  were validated along with ACTB, GAPDH, RNA18S, PPIA, PGK1, RPLP0 and POLR2A by qPCR assessment of their expression levels on 18 paired normal lung and NSCLC samples. Expression changes between and within these two classes were investigated to define the most stable and comprehensive RGs.
Primary tumour samples and paired normal lung tissues obtained from 18 NSCLC patients (11M, 7F) aged 41–79 (mean 62) during primary curative resection (not preceded by radiotherapy or chemotherapy) at the San Luigi Hospital, Orbassano (Italy) between December 2003 and March 2004 were immediately snap-frozen in liquid nitrogen and stored at -80°C until RNA extraction. The histological subtypes were: 7 adenocarcinomas (ADCA, 38.9%), 10 squamous cell carcinomas (SQCA, 55.6%) and 1 broncoalveolar carcinoma (BAC, 5.5%). Prior written and informed consent was obtained from each patient. The study was approved by the appropriate ethical review board. The guidelines for good clinical practice and the recommendations of the Declaration of Helsinki for biomedical research involving human subjects were followed.
RNA extraction and cDNA synthesis
Total RNA (totRNA) was isolated from lung specimens with the RNeasy 96 Kit and Biorobot 8000 (Qiagen, Hilden, Germany) according to the manufacturer's instructions. RNA was extracted from 15–25 mg and 60–80 mg of tumor and normal lung tissues specimens respectively. To take account of the expression variability due to cellular heterogeneity, totRNA was extracted from a biological duplicate of both normal and tumour specimens. Genomic DNA contaminations were removed by on-column-DNAseI treatment. totRNA was then quantified with an Agilent 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA, USA) and stored at -80°C. Two μgr totRNA were finally retro-transcribed with random hexamer primers and Multiscribe Reverse transcriptase (MMLV) contained in the High Capacity cDNA Archive Kit (Applied Biosystems, Foster City, CA), in accordance with the manufacturer's suggestions.
Gene expression assays.
Assay on Demand
basic transcription factor 3
Esterase D/formylglutathione hydrolase
histone 1, H2bc
ribosomal protein L30
Yes-associated protein 1, 65kDa
polymerase (RNA) II polypeptide A, 220kDa
phosphoglycerate kinase 1
18S ribosomal RNA
ribosomal protein, large, P0
breast cancer 1, early onset
excision repair cross-complementing rodent repair deficiency, complementation group 2
ribonucleotide reductase M2 polypeptide
Meta-analysis of expression data was performed by using affy, MergeMaid and genefilter libraries implemented in Bioconductor . Statistical analyses were performed with R . Gene differential expression was considered significant with an α = 0.05 and p-value <0.05. PCA (Principal Component Analysis, ) and agglomerative hierarchical clustering were performed with TMEV software . Assays on Demand's PCR efficiencies, evaluated by cDNA dilution curves, ranged between 91% (rRNA18S, RPL30) and 100% (POL2RA, ESD) with an average value of 96%. Since PCR efficiencies were high and comparable we used the theoretical value of 2 (100%) for gene expression quantification. So, fold changes in expression levels between the normal and tumor paired samples were evaluated with 2-Δ'Ct . The expression of each target gene was reported as 2-ΔΔCt.
Biological replicates were treated as independent samples in RGs raw Ct distribution analysis and geNorm analysis. Their average Ct value was used in Δ'Ct calculation (i.e for each sample mean Ct_tumor-mean Ct_normal) applied to evaluate the RGs expression stability with equivalence test. We also averaged the ΔCt value (Ct target-Ct reference) of biological replicates to generate the mean ΔΔCt values used in PCA and agglomerative hierarchical clustering according to the Euclidean distance metric and average linkage method.
Five putative RGs were selected from a meta-analysis of a subset of lung tissue transcription profilings . To secure a balanced set of tissues, the analysis focused on 17 normal samples, 18 ADCA (randomly selected from a total of 127), 21 SQCA, 20 pulmonary carcinoids and 6 small-cell lung carcinomas. After RMA intensity calculation  and quantile normalization , an intensity filter was applied to remove all probe sets with intensity ≤100 in all experiments.
Tumour/normal gene expression similarity was measured with the "integrative correlation" (IC) estimation, developed by Parmigiani et al. . This reflects the general consistency of gene expression among the different tissues since it is generated using gene-to-gene correlation obtained by all pairwise correlation for all genes across the tissue groups. A variance/IC filter between each tumor group versus the normal tissues was applied to select probe sets with an inter-sample standard deviation (SD) ≤ 0.5 and IC ≥ 0.2. All the 28 top ranked putative RGs displayed a nearly normal distribution of the residuals representing the expression differences between each tumor sample versus normal samples (data not shown). Within this set we selected five candidates RGs (ESD, BTF3, HIST1H2BC, RPL30, YAP1) involved in different biological processes and for which commercial gene expression assays were available.
RGs expression ranges in lung tissues
RGs expression stability
To allow accurate quantification a RG should own a stable expression among samples while to warrant a trustworthy quantification it should not be modulated by experimental conditions. Ideally to fulfill both conditions RG raw Ct distribution of both normal and tumor populations should thus be characterized by a narrow range of expression with comparable mean values. With reference to stability of expression in both normal and in tumor samples, the lowest RNA transcription ranges (i.e Ct 75°quantile-Ct 25°quantile) were associated withPOLR2A (normal = 0.73; tumor = 0.69) and rRNA18S (normal = 0.69; tumor = 0.93), whereas the highest ranges were detected for GAPDH (normal = 1.8; tumor = 2.2) and RPLP0 (normal = 1.71; tumor = 1.84). Unfortunately, even when equal amount of RNA input are used, retro-transcription reactions are responsible for most of the variation in the experimental determination of mRNA quantities [31, 32]. We are thus unable to estimate how much the technical (i.e. RT efficiency) and the biological (i.e. inter-samples) variability contribute to the measured difference in RNA quantities. Descriptive analysis of raw Ct distributions provides only a rough estimate of RG stability that must be confirmed by other approaches.
To investigate the experimental stability of RG transcription we evaluated the significance of differences in RNA transcription levels between the tumor and normal paired samples, where the latter were used as calibrators in Δ'Ct quantification approach [7, 33]. Since in the two-sided t-test the type II errors (false negatives) are uncontrollable, we performed an equivalence test for dependent samples  to investigate the null hypothesis of inequivalence H0: μT-μN ≤ θL and μT-μN ≥ θU .
Significance of suitable RGs for normalization
This paper describes the first systematic comparison of several conventional and potential RGs and their effectiveness as internal control for relative quantification in lung cancer gene expression profiling. To increase data reliability, we sought to control several variables by (1) using matched pairs of normal and tumor lung samples instead of unpaired samples to minimize inter-individual variations that could arise from genotype profile; (2) analysing biological duplicates of each sample to take into account the expression variability due to cellular heterogeneity in both normal and tumor lung tissues; (3) employing a liquid handling robot for most of the technical procedures (i.e totRNA extraction and qPCR plates assembling), to avoid operator variability; (4) by using high quality certified qPCR reagents and (5) commercial RG probes (Taqman Gene Expression Assay) properly designed to guarantee maximum PCR efficiency and specificity  to allow trustworthy comparison in gene expression stability evaluation. Furthermore, in view of the unsuitability of conventional RGs, such as GAPDH or ACTB, for most of the experimental situations [10, 17–20], we chose a microarray data meta-analysis to identify new potential RGs. Since 2 out of 4 of the best performing RGs (ESD and YAP1) identified in the present study belong to the meta-analysis group, we can confirm the reliability of this strategy in the identification of new candidate RGs within a specific experimental condition. A further advantage is the easier identification of candidate RGs expressed at a comparable order of magnitude to target genes. This might result in an improvement of data normalization since, as highlighted by some researchers [38, 39], at the same Ct level, reference and target experience the same condition and real-time RT-PCR kinetics with respect to polymerase activation, reaction inactivation, stochastic relation between target and primer concentration, and reaction end product inhibition by the generated RT-PCR product. However, the filtering procedures used to search microarray data for stably expressed genes are biased for low expressed genes.
RG expression stability within sample populations and under experimental condition (tumor versus normal lung specimens) has been evaluated by: (1) descriptive statistic; (2) equivalence test; (3) GeNorm applet. All approaches matched and show that for gene expression profiling in NSCLC the most suitable reference genes to be used for normalization are rRNA18S, POLR2A, ESD and YAP1. The effectiveness of POLR2A, the main enzyme in mRNA transcription, is further supported by others . It is characterized by stable expression among different tissues and it is not modulated by treatment with TPA and ionomycin, indicating resistance to cellular activation. The efficacy of rRNA18S, however, has been criticized because: (1) it is highly expressed, representing up to 80% of cellular RNA and (2) it is transcribed by a specific RNA-polymerase . Nevertheless the results of the present study clearly show that there are no significant differences between rRNA18s and the other best performing RGs identified. In addition, rRNA18s has already been recommended for gene expression analysis in non-microdissected kidney biopsies . A recent assessment of RGs for lung cancer studies by Liu et al.  has indicated that GAPDH is the most suitable RG for qPCR studies in NSCLC tissue samples. In our study, however, all the statistical approaches used suggested that GAPDH is far from being a valid RG due to its evident modulation by the neoplastic transformation. Our data are strongly in agreement with previous indication that GAPDH is regulated in response to various stimuli (i.e hypoxia, insulin, mitogen, EGF) and that its transcript is elevated in various cancer tissues [43–45], including lung cancer . The discrepancy with Liu's data could be related to the different ethnical origin of the specimens. Even so, it probably reflects the employment of a different data analysis approach. Liu et al., in fact, defined GAPDH as the optimal RG because its expression was almost equal to the mean expression of the other six endogenous control genes they analyzed. This is not a satisfactory parameter of RG stability because it is closely dependent on the variability of the remaining RGs and could not be used to assess their constancy within investigated samples. Furthermore all their results refer to pooled tumor and normal lung tissues samples and do not consider whether and how high RG expression levels were affected by malignant transformation. The latter assessment is mandatory to warrant the validity of quantification data .
Careful validation of RGs has been repeatedly advocated as mandatory to ensure the accuracy of normalized data in qPCR studies adopting relative quantification. Our results invite the conclusion that for gene expression profiling in NSCLC the RGs most suitable for normalization are rRNA18S, POLR2A, ESD and YAP1. For accurate normalization we recommend their concurrent use since this results in more reliable quantification. Furthermore our results pinpoint the limits of conventional RGs and suggest that relatively simple meta-data analysis of genome wide transcription profiling studies is a useful way of identifying putative RGs.
This work was supported by "Oncongenomic centers" AIRC grant and Fondazione San Paolo. Fellowship for SS was granted by "Programma per la Ricerca Sanitaria 2003", Ministero della Salute.
- Jemal A, Tiwari RC, Murray T, Ghafoor A, Samuels A, Ward E, Feuer EJ, Thun MJ, American Cancer S: Cancer statistics, 2004. CA: a Cancer Journal for Clinicians. 2004, 54 (1): 8-29.Google Scholar
- Mather D, Sullivan SD, Parasuraman TV: Beyond survival: economic analyses of chemotherapy in advanced inoperable NSCLC. Oncology (Huntingt). 1998, 12: 199-209.Google Scholar
- Schiller JH, Harrington D, Belani CP, Langer C, Sandler A, Krook J, Zhu J, Johnson DH, Eastern Cooperative Oncology G: Comparison of four chemotherapy regimens for advanced non-small-cell lung cancer.[see comment]. New England Journal of Medicine. 2002, 346 (2): 92-98. 10.1056/NEJMoa011954.View ArticlePubMedGoogle Scholar
- Bepler G, Sharma S, Cantor A, Gautam A, Haura E, Simon G, Sharma A, Sommers E, Robinson L: RRM1 and PTEN as prognostic parameters for overall and disease-free survival in patients with non-small-cell lung cancer. Journal of Clinical Oncology. 2004, 22 (10): 1878-1885. 10.1200/JCO.2004.12.002.View ArticlePubMedGoogle Scholar
- Rosell R, Felip E, Taron M, Majo J, Mendez P, Sanchez-Ronco M, Queralt C, Sanchez JJ, Maestre J: Gene expression as a predictive marker of outcome in stage IIB-IIIA-IIIB non-small cell lung cancer after induction gemcitabine-based chemotherapy followed by resectional surgery. Clinical Cancer Research. 2004, 10 (12 Pt 2): 4215s-4219s. 10.1158/1078-0432.CCR-040006.View ArticlePubMedGoogle Scholar
- Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods (Duluth). 2001, 25 (4): 402-408.View ArticleGoogle Scholar
- Schmittgen TD, Zakrajsek BA: Effect of experimental treatment on housekeeping gene expression: validation by real-time, quantitative RT-PCR. Journal of Biochemical & Biophysical Methods. 2000, 46 (1–2): 69-81. 10.1016/S0165-022X(00)00129-9.View ArticleGoogle Scholar
- Huggett J, Dheda K, Bustin S, Zumla A: Real-time RT-PCR normalisation; strategies and considerations. Genes & Immunity. 2005, 6 (4): 279-284. 10.1038/sj.gene.6364190.View ArticleGoogle 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 suited for normalization, applied to bladder and colon cancer data sets. Cancer Research. 2004, 64 (15): 5245-5250. 10.1158/0008-5472.CAN-04-0496.View ArticlePubMedGoogle Scholar
- Schmid H, Cohen CD, Henger A, Irrgang S, Schlondorff D, Kretzler M: Validation of endogenous controls for gene expression analysis in microdissected human renal biopsies.[see comment]. Kidney International. 2003, 64 (1): 356-360. 10.1046/j.1523-1755.2003.00074.x.View ArticlePubMedGoogle Scholar
- Warrington JA, Nair A, Mahadevappa M, Tsyganskaya M: Comparison of human adult and fetal expression and identification of 535 housekeeping/maintenance genes. Physiological Genomics. 2000, 2 (3): 143-147.PubMedGoogle Scholar
- Muller-Tidow C, Diederichs S, Bulk E, Pohle T, Steffen B, Schwable J, Plewka S, Thomas M, Metzger R, Schneider PM, et al: Identification of metastasis-associated receptor tyrosine kinases in non-small cell lung cancer. Cancer Research. 2005, 65 (5): 1778-1782. 10.1158/0008-5472.CAN-04-3388.View ArticlePubMedGoogle Scholar
- Brabender J, Metzger R, Salonga D, Danenberg KD, Danenberg PV, Holscher AH, Schneider PM: Comprehensive expression analysis of retinoic acid receptors and retinoid X receptors in non-small cell lung cancer: implications for tumor development and prognosis. Carcinogenesis. 2005, 26 (3): 525-530. 10.1093/carcin/bgi006.View ArticlePubMedGoogle Scholar
- Yuan A, Yu CJ, Shun CT, Luh KT, Kuo SH, Lee YC, Yang PC: Total cyclooxygenase-2 mRNA levels correlate with vascular endothelial growth factor mRNA levels, tumor angiogenesis and prognosis in non-small cell lung cancer patients. International Journal of Cancer. 2005, 115 (4): 545-555. 10.1002/ijc.20898.View ArticlePubMedGoogle Scholar
- Falleni M, Pellegrini C, Marchetti A, Roncalli M, Nosotti M, Palleschi A, Santambrogio L, Coggi G, Bosari S: Quantitative evaluation of the apoptosis regulating genes Survivin, Bcl-2 and Bax in inflammatory and malignant pleural lesions. Lung Cancer. 2005, 48 (2): 211-216. 10.1016/j.lungcan.2004.10.003.View ArticlePubMedGoogle Scholar
- Parekh K, Ramachandran S, Cooper J, Bigner D, Patterson A, Mohanakumar T: Tenascin-C, over expressed in lung cancer down regulates effector functions of tumor infiltrating lymphocytes. Lung Cancer. 2005, 47 (1): 17-29. 10.1016/j.lungcan.2004.05.016.View ArticlePubMedGoogle Scholar
- Bas A, Forsberg G, Hammarstrom S, Hammarstrom ML: Utility of the housekeeping genes 18S rRNA, beta-actin and glyceraldehyde-3-phosphate-dehydrogenase for normalization in real-time quantitative reverse transcriptase-polymerase chain reaction analysis of gene expression in human T lymphocytes. Scandinavian Journal of Immunology. 2004, 59 (6): 566-573. 10.1111/j.0300-9475.2004.01440.x.View ArticlePubMedGoogle Scholar
- Goidin D, Mamessier A, Staquet MJ, Schmitt D, Berthier-Vergnes O: Ribosomal 18S RNA prevails over glyceraldehyde-3-phosphate dehydrogenase and beta-actin genes as internal standard for quantitative comparison of mRNA levels in invasive and noninvasive human melanoma cell subpopulations. Analytical Biochemistry. 2001, 295 (1): 17-21. 10.1006/abio.2001.5171.View ArticlePubMedGoogle Scholar
- Zhong H, Simons JW: Direct comparison of GAPDH, beta-actin, cyclophilin, and 28S rRNA as internal standards for quantifying RNA levels under hypoxia. Biochemical & Biophysical Research Communications. 1999, 259 (3): 523-526. 10.1006/bbrc.1999.0815.View ArticleGoogle Scholar
- Bustin SA: Absolute quantification of mRNA using real-time reverse transcription polymerase chain reaction assays. Journal of Molecular Endocrinology. 2000, 25 (2): 169-193. 10.1677/jme.0.0250169.View ArticlePubMedGoogle Scholar
- Bhattacharjee A, Richards WG, Staunton J, Li C, Monti S, Vasa P, Ladd C, Beheshti J, Bueno R, Gillette M, et al: Classification of human lung carcinomas by mRNA expression profiling reveals distinct adenocarcinoma subclasses. Proceedings of the National Academy of Sciences of the United States of America. 2001, 98 (24): 13790-13795. 10.1073/pnas.191502998.View ArticlePubMedPubMed CentralGoogle Scholar
- Bioconductor. [http://www.bioconductor.org]
- R-project. [http://www.r-project.org]
- Raychaudhuri S, Stuart JM, Altman RB: Principal components analysis to summarize microarray experiments: application to sporulation time series. Pacific Symposium on Biocomputing. 2000, 455-466.Google Scholar
- Saeed AI, Sharov V, White J, Li J, Liang W, Bhagabati N, Braisted J, Klapa M, Currier T, Thiagarajan M, et al: TM4: a free, open-source system for microarray data management and analysis. Biotechniques. 2003, 34 (2): 374-378.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 control genes. Genome Biology. 2002, 3 (7): RESEARCH0034-10.1186/gb-2002-3-7-research0034.View ArticlePubMedPubMed CentralGoogle Scholar
- Gene Quantification Home Page. [http://www.gene-quantification.org]
- Irizarry RA, Hobbs B, Collin F, Beazer-Barclay YD, Antonellis KJ, Scherf U, Speed TP: Exploration, normalization, and summaries of high density oligonucleotide array probe level data. Biostatistics. 2003, 4: 249-264. 10.1093/biostatistics/4.2.249.View ArticlePubMedGoogle Scholar
- Bolstad BM, Irizarry RA, Astrand M, Speed TP: A comparison of normalization methods for high density oligonucleotide array data based on variance and bias. Bioinformatics. 2003, 19 (2): 185-193. 10.1093/bioinformatics/19.2.185.View ArticlePubMedGoogle Scholar
- Parmigiani G, Garrett-Mayer ES, Anbazhagan R, Gabrielson E: A cross-study comparison of gene expression studies for the molecular classification of lung cancer. Clinical Cancer Research. 2004, 10 (9): 2922-2927. 10.1158/1078-0432.CCR-03-0490.View ArticlePubMedGoogle Scholar
- Stahlberg A, Hakansson J, Xian X, Semb H, Kubista M: Properties of the reverse transcription reaction in mRNA quantification. Clinical Chemistry. 2004, 50 (3): 509-515. 10.1373/clinchem.2003.026161.View ArticlePubMedGoogle Scholar
- Curry J, McHale C, Smith MT: Low efficiency of the Moloney murine leukemia virus reverse transcriptase during reverse transcription of rare t(8;21) fusion gene transcripts. Biotechniques. 2002, 32 (4): 768-775.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 (1): 112-114.PubMedGoogle Scholar
- Haller F, Kulle B, Schwager S, Gunawan B, von Heydebreck A, Sultmann H, Fuzesi L: Equivalence test in quantitative reverse transcription polymerase chain reaction: confirmation of reference genes suitable for normalization. Analytical Biochemistry. 2004, 335 (1): 1-9. 10.1016/j.ab.2004.08.024.View ArticlePubMedGoogle Scholar
- McBride GB: Equivalence tests can enhance enviromental, science, and management. Austral & New Zeland J Statist. 1999, 4 (1): 19-29. 10.1111/1467-842X.00058.View ArticleGoogle Scholar
- Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, De Paepe A, Speleman F: GeNorm software manual, update 6 sept 2004. [http://medgenugentbe/~jvdesomp/genorm/]
- Amplification efficiency of TaqMan Gene Expression assays. Application Note Applied Biosystems.Google Scholar
- Bustin SA: Absolute quantification of mRNA using real-time reverse transcription chain reaction assay. J Mol Endocrinol. 2000, 25: 169-193. 10.1677/jme.0.0250169.View ArticlePubMedGoogle Scholar
- Quantification strategies in real-time PCR. Chapter 3 of "A-Z of quantitative PCR".Google Scholar
- Radonic A, Thulke S, Mackay IM, Landt O, Siegert W, Nitsche A: Guideline to reference gene selection for quantitative real-time PCR. Biochemical & Biophysical Research Communications. 2004, 313 (4): 856-862. 10.1016/j.bbrc.2003.11.177.View ArticleGoogle Scholar
- Thellin O, Zorzi W, Lakaye B, De Borman B, Coumans B, Hennen G, Grisar T, Igout A, Heinen E: Housekeeping genes as internal standards: use and limits. Journal of Biotechnology. 1999, 75 (2–3): 291-295. 10.1016/S0168-1656(99)00163-7.View ArticlePubMedGoogle Scholar
- Liu DWCST, Liu HP: Choice of endogenous control for gene expression in nonsmall cell lung cancer. Eur Respir J. 2005, 26 (6): 1002-1008. 10.1183/09031936.05.00050205.View ArticlePubMedGoogle Scholar
- Revillion F, Pawlowski V, Hornez L, Peyrat JP: Glyceraldehyde-3-phosphate dehydrogenase gene expression in human breast cancer. European Journal of Cancer. 2000, 36 (8): 1038-1042. 10.1016/S0959-8049(00)00051-4.View ArticlePubMedGoogle Scholar
- Rondinelli RHEDE, Tricoli JV: Increased glyceraldehyde-3-phosphate dehydrogenase gene expression in late pathological stage human prostate cancer. Prostate Cancer Prostatic Dis. 1997, 1 (2): 66-72. 10.1038/sj.pcan.4500208.View ArticlePubMedGoogle Scholar
- Schek NHBL, Finn OJ: Increased glyceraldehyde-3-phosphate dehydrogenase gene expression in human pancreatic adenocarcinoma. Cancer Research. 1988, 48 (22): 6354-6359.PubMedGoogle Scholar
- Tokunaga K, Nakamura Y, Sakata K, Fujimori K, Ohkubo M, Sawada K, Sakiyama S: Enhanced expression of a glyceraldehyde-3-phosphate dehydrogenase gene in human lung cancers. Cancer Research. 1987, 47 (21): 5616-5619.PubMedGoogle Scholar
- Dheda K, Huggett JF, Chang JS, Kim LU, Bustin SA, Johnson MA, Rook GA, Zumla A: The implications of using an inappropriate reference gene for real-time reverse transcription PCR data normalization. Analytical Biochemistry. 2005, 344 (1): 141-143. 10.1016/j.ab.2005.05.022.View ArticlePubMedGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2407/6/200/prepub
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.