Accuracy and calibration of commercial oligonucleotide and custom cDNA microarrays.

We compared the accuracy of microarray measurements obtained with oligonucleotide arrays (GeneChip, Affymetrix) with a laboratory-developed cDNA array by assaying test RNA samples from an experiment using a paradigm known to regulate many genes measured on both arrays. We selected 47 genes represented on both arrays, including both known regulated and unregulated transcripts, and established reference relative expression measurements for these genes in the test RNA samples using quantitative reverse transcriptase real-time PCR (QRTPCR) assays. The validity of the reproducible (average coefficient of variation = 11.8%) QRTPCR measurements were established through application of a new mathematical model. The performance of both array platforms in identifying regulated and non-regulated genes was identical. With either platform, 16 of 17 definitely regulated genes were correctly identified, and no definitely unregulated transcript was falsely identified as regulated. Accuracy of the fold-change measurements obtained with each platform was assessed by determining measurement bias. Both platforms consistently underestimate the relative changes in mRNA expression between experimental and control samples. The bias observed with cDNA arrays was predictable for fold-changes <250-fold by QRTPCR and could be corrected by the calibration function F(c) = F(a(cDNA))(q), where F(a(cDNA)) is the microarray-determined fold-change comparing experimental with control samples, q is the correction factor and F(c) is the calibrated value. The bias observed with the commercial oligonucleotide arrays was less predictable and calibration was unfeasible. Following calibration, fold-change measurements generated by custom cDNA arrays were more accurate than those obtained by commercial oligonucleotide arrays. Our study demonstrates systematic bias of microarray measurements and identifies a calibration function that improves the accuracy of cDNA array data.

[1]  P. Chomczyński,et al.  Single-step method of RNA isolation by acid guanidinium thiocyanate-phenol-chloroform extraction. , 1987, Analytical biochemistry.

[2]  Calvin L. Williams,et al.  Modern Applied Statistics with S-Plus , 1997 .

[3]  P. Brown,et al.  Yeast microarrays for genome wide parallel genetic and gene expression analysis. , 1997, Proceedings of the National Academy of Sciences of the United States of America.

[4]  J. Mesirov,et al.  Molecular classification of cancer: class discovery and class prediction by gene expression monitoring. , 1999, Science.

[5]  L. Overbergh,et al.  Quantification of murine cytokine mRNAs using real time quantitative reverse transcriptase PCR. , 1999, Cytokine.

[6]  J. R. Koehler,et al.  Modern Applied Statistics with S-Plus. , 1996 .

[7]  D. Botstein,et al.  The transcriptional program in the response of human fibroblasts to serum. , 1999, Science.

[8]  D. Botstein,et al.  Exploring the new world of the genome with DNA microarrays , 1999, Nature Genetics.

[9]  S. P. Fodor,et al.  High density synthetic oligonucleotide arrays , 1999, Nature Genetics.

[10]  S. Walker,et al.  Quantitative RT-PCR : Pitfalls and Potential , 1999 .

[11]  M. Ko,et al.  Genome-wide expression profiling of mid-gestation placenta and embryo using a 15,000 mouse developmental cDNA microarray. , 2000, Proceedings of the National Academy of Sciences of the United States of America.

[12]  Ash A. Alizadeh,et al.  Distinct types of diffuse large B-cell lymphoma identified by gene expression profiling , 2000, Nature.

[13]  S. Bustin Absolute quantification of mRNA using real-time reverse transcription polymerase chain reaction assays. , 2000, Journal of molecular endocrinology.

[14]  N. Walker Real‐time and quantitative PCR: applications to mechanism‐based toxicology , 2001, Journal of biochemical and molecular toxicology.

[15]  K. Miura,et al.  Quantitative assessment of DNA microarrays--comparison with Northern blot analyses. , 2001, Genomics.

[16]  S. Sealfon,et al.  Gonadotropin-releasing Hormone Receptor-coupled Gene Network Organization * , 2001 .

[17]  M. Pfaffl,et al.  A new mathematical model for relative quantification in real-time RT-PCR. , 2001, Nucleic acids research.

[18]  R. Nuttall,et al.  An evaluation of the performance of cDNA microarrays for detecting changes in global mRNA expression. , 2001, Nucleic acids research.

[19]  S. Sealfon,et al.  Monitoring G-protein-coupled receptor signaling with DNA microarrays and real-time polymerase chain reaction. , 2002, Methods in enzymology.