Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
Comparative Study
. 2007 Jun 1:8:140.
doi: 10.1186/1471-2164-8-140.

Selection of DDX5 as a novel internal control for Q-RT-PCR from microarray data using a block bootstrap re-sampling scheme

Affiliations
Comparative Study

Selection of DDX5 as a novel internal control for Q-RT-PCR from microarray data using a block bootstrap re-sampling scheme

Li-Jen Su et al. BMC Genomics. .

Abstract

Background: The development of microarrays permits us to monitor transcriptomes on a genome-wide scale. To validate microarray measurements, quantitative-real time-reverse transcription PCR (Q-RT-PCR) is one of the most robust and commonly used approaches. The new challenge in gene quantification analysis is how to explicitly incorporate statistical estimation in such studies. In the realm of statistical analysis, the various available methods of the probe level normalization for microarray analysis may result in distinctly different target selections and variation in the scores for the correlation between microarray and Q-RT-PCR. Moreover, it remains a major challenge to identify a proper internal control for Q-RT-PCR when confirming microarray measurements.

Results: Sixty-six Affymetrix microarray slides using lung adenocarcinoma tissue RNAs were analyzed by a statistical re-sampling method in order to detect genes with minimal variation in gene expression. By this approach, we identified DDX5 as a novel internal control for Q-RT-PCR. Twenty-three genes, which were differentially expressed between adjacent normal and tumor samples, were selected and analyzed using 24 paired lung adenocarcinoma samples by Q-RT-PCR using two internal controls, DDX5 and GAPDH. The percentage correlation between Q-RT-PCR and microarray were 70% and 48% by using DDX5 and GAPDH as internal controls, respectively.

Conclusion: Together, these quantification strategies for Q-RT-PCR data processing procedure, which focused on minimal variation, ought to significantly facilitate internal control evaluation and selection for Q-RT-PCR when corroborating microarray data.

PubMed Disclaimer

Figures

Figure 1
Figure 1
Flow chart for prioritization of potential internal controls.
Figure 2
Figure 2
Bootstrap box plots of the gene expression intensity of various internal controls. (A) The box plot results show the best 14 internal control candidates, all of which exhibited consistent expression intensity in the NHRI lung adenocarcinoma microarray dataset for each re-sampling process. Moreover, also included are 10 well-known Q-RT-PCR internal controls contained in 23 probe sets on the HG-U133A chip. These are shown as #1–23 in x-axis. The detailed probe set characteristics were shown in Table 1. Except ABCF1, BHLHB2 and LAPTM4A, the gene expression intensities of top 12 internal control candidates, GAPDH, and ACTB from the Boston (B) and the Ann Arbor lung cancer datasets (C) were also compared. DDX5: (DEAD (Asp-Glu-Ala-Asp) box polypeptide 5), PKM2: (pyruvate kinase, muscle), BHLHB2: (basic helix-loop-helix domain containing, class B, 2), GLO1: (glyoxalase I), LAPTM4A: (lysosomal-associated protein transmembrane 4 alpha), SET: (SET translocation (myeloid leukemia-associated)), CLTC: (clathrin, heavy chain (Hc)), MSN: (MSN/ALK fusion; moesin/anaplastic lymphoma kinase fusion protein), ABCF1: (ATP-binding cassette, sub-family F (GCN20), member 1), EPHB3: (EPH receptor B3), CCL5: (chemokine (C-C motif) ligand 5), PTPN21: (protein tyrosine phosphatase, non-receptor type 21), DDR1: (discoidin domain receptor family, member 1), 1–4: ACTB (actin, beta), 5–6: B2M (beta-2-microglobulin), 7–12: GAPDH (glyceraldehyde-3-phosphate dehydrogenase), 13: HMBS (hydroxymethylbilane synthase), 14: HPRT1 (hypoxanthine phosphoribosyltransferase 1), 15–19: RPL13A (ribosomal protein L13a), 20: RPL32 (ribosomal protein L32), 21: SDHA (succinate dehydrogenase complex, subunit A, flavoprotein (Fp)), 22: UBC (ubiquitin C), 23: YWHAZ (tyrosine 3-monooxygenase/tryptophan 5-monooxygenase activation protein, zeta polypeptide).
Figure 3
Figure 3
Bootstrap replicates of the descriptive statistics related to the variation in gene expression. In total, 27 pairs of adjacent normal and tumor samples with 12 un-paired samples were used in this analysis. A total of 39 blocks of microarray samples were used for the block bootstrap. The re-sampling process was repeated 1,000 times to obtain 1,000 bootstrap replicates of the minimum (green color), first quartile (red color), median (black color), third quartile (blue color) and maximum (cyan color) expression levels for each gene. Each result was ranked by the order of medians. The bootstrap replicates of all five statistics for DDX5 expression remain roughly constants, but those for GAPDH expression vary greatly.
Figure 4
Figure 4
Gene expression patterns of DDX5 and GAPDH from 66 Affymetrix microarray chips using three different probe level quantile normalizations. The gene expression patterns of two internal controls, DDX5 (A) and GAPDH (B: 212581_x_at and C: M33197_3_at) from Affymetrix chips by MAS5 (marked by diamonds), RMA (marked by squares) and GC-RMA (marked by triangles) probe level quantile normalizations. Normalization was performed per chip and per gene using GeneSpring® 7.3 software. The expression levels of these two internal controls were related to the median of the intensities on the 66 chips. The DDX5 expression patterns in each chip did not significantly alter compared to greater variation in GAPDH.
Figure 5
Figure 5
Box plots of DDX5 relative expression patterns exhibit small variation across various microarray datasets. The gene expression patterns of DDX5 were obtained from other microarray databases. These datasets were from 84 cancer cell lines (NCI60), which were classified into 8 different cancer cell types and other cell lines (A). GAPDH #5: M33197_3_at; GAPDH #10: 212581_x_at. The box plot results indicated that DDX5 exhibited only small variation across the various NCI60 cell types. For the lung cancer cell lines, DDX5, CLTC and MSN all gave lower variances than ACTB and GAPDH (B). Both ACTB (#1–4) and GAPDH (#5–10) are contained in 10 probe sets on the HG-U133A chip and are shown as #1–10 on x-axis and in Table 1. The variation of DDX5 was also smaller for the HeLa cell cycle dataset, lung cancer dataset and HCC dataset from the Stanford Microarray Database (C). BR: breast cancer, CO: colon cancer, GL: glioblastoma, KI: Kidney, LE: leukemia, LU: lung cancer, ME: melanoma and OV: ovarian cancer.

References

    1. Jain KK. Applications of biochips: from diagnostics to personalized medicine. Curr Opin Drug Discov Devel. 2004;7:285–289. - PubMed
    1. Lipshutz RJ, Fodor SP, Gingeras TR, Lockhart DJ. High density synthetic oligonucleotide arrays. Nat Genet. 1999;21:20–24. doi: 10.1038/4447. - DOI - PubMed
    1. Bustin SA. Absolute quantification of mRNA using real-time reverse transcription polymerase chain reaction assays. J Mol Endocrinol. 2000;25:169–193. doi: 10.1677/jme.0.0250169. - DOI - PubMed
    1. Yun JJ, Heisler LE, Hwang. 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 Res. 2006;34:e85. doi: 10.1093/nar/gkl400. - DOI - PMC - PubMed
    1. Lee LG, Connell CR, Bloch W. Allelic discrimination by nick-translation PCR with fluorogenic probes. Nucleic Acids Res. 1993;21:3761–3766. doi: 10.1093/nar/21.16.3761. - DOI - PMC - PubMed

Publication types

MeSH terms