Optimization of cDNA microarrays procedures using criteria that do not rely on external standards
- Torunn Bruland†1Email author,
- Endre Anderssen†1, 2,
- Berit Doseth1,
- Hallgeir Bergum2,
- Vidar Beisvag1, 2 and
- Astrid Lægreid1, 2
© Bruland et al; licensee BioMed Central Ltd. 2007
Received: 05 March 2007
Accepted: 18 October 2007
Published: 18 October 2007
The measurement of gene expression using microarray technology is a complicated process in which a large number of factors can be varied. Due to the lack of standard calibration samples such as are used in traditional chemical analysis it may be a problem to evaluate whether changes done to the microarray procedure actually improve the identification of truly differentially expressed genes. The purpose of the present work is to report the optimization of several steps in the microarray process both in laboratory practices and in data processing using criteria that do not rely on external standards.
We performed a cDNA microarry experiment including RNA from samples with high expected differential gene expression termed "high contrasts" (rat cell lines AR42J and NRK52E) compared to self-self hybridization, and optimized a pipeline to maximize the number of genes found to be differentially expressed in the "high contrasts" RNA samples by estimating the false discovery rate (FDR) using a null distribution obtained from the self-self experiment. The proposed high-contrast versus self-self method (HCSSM) requires only four microarrays per evaluation. The effects of blocking reagent dose, filtering, and background corrections methodologies were investigated. In our experiments a dose of 250 ng LNA (locked nucleic acid) dT blocker, no background correction and weight based filtering gave the largest number of differentially expressed genes. The choice of background correction method had a stronger impact on the estimated number of differentially expressed genes than the choice of filtering method. Cross platform microarray (Illumina) analysis was used to validate that the increase in the number of differentially expressed genes found by HCSSM was real.
The results show that HCSSM can be a useful and simple approach to optimize microarray procedures without including external standards. Our optimizing method is highly applicable to both long oligo-probe microarrays which have become commonly used for well characterized organisms such as man, mouse and rat, as well as to cDNA microarrays which are still of importance for organisms with incomplete genome sequence information such as many bacteria, plants and fish.
Gene-expression microarrays are widely used for large scale studies in cells, and are emerging as promising tools in clinical diagnosis, with potential impact on the assessment of prognosis and choice of treatment . Microarray data are also increasingly becoming essential for data driven model building in system biology approaches [2, 3]. However, the measurement of gene expression using e.g. cDNA microarray technology is a complex process that involves several steps, which have not yet been fully optimized. In a typical cDNA microarray experiment, total RNA from two biological sources in which gene expression levels are to be compared, is labelled with two different fluorescent dyes and hybridized to an array of spotted cDNA probes. After hybridization and scanning the data are further pre-processed. This often includes image analysis, filtering of data and background subtraction before normalization and transformation of data. Relative gene expression levels are given as ratios of intensities of fluorescent emission from labelled the RNAs'. Many different approaches and methods can be used at each stage of a microarray experiment, and there is still no definitive consensus about which methods to choose (reviewed and discussed in ).
Sample quality, labeling protocol, hybridization conditions, scanning protocols and image acquisition as well as the different stages in the microarray data analysis pipeline can all contribute to the overall uncertainty of the conclusions drawn. There has therefore recently been an emerging focus on the need for universally applicable standards, reference materials and analytic guidelines to assist in the standardization of microarray experiments [4–10]. In order to choose or optimize microarray methods several criteria can be used. The use of external RNA spike-in controls is the most recommended approach for technology assessment and optimization. However, the recommended use of spike-in RNA might not always be feasible. For optimal application several controls and probes must be used and the spike-in controls must be representative of the endogenous RNA with respect to e.g. length, sequence characteristics, melting temperature and cross-hybridization risk . On some arrays suitable spike-in controls are not printed, and printing a sufficient number of spikes can leave insufficient room for the genes of biological interest. Other valuable tools are calibrated reference RNA samples such as those developed in the MAQC project  or mixed tissue RNA samples developed by Thompson et al. [11, 12]. However, it is not always possible to have RNA samples with known differences in the expression level of many genes from the species of interest.
There is a large number of steps in the microarray production, hybridization, scanning and data analysis that have adjustable parameters that may be optimized. We will show how HCSSM can be used to optimize data analysis and hybridization protocols. As mentioned above, the standard microarray data analysis pipeline consists of image analysis, filtering of data, background correction, normalization and identification of differentially expressed genes. Recent studies have focused on the normalization step and evaluated a large number of normalization methods for cDNA microarray data [15, 16]. We have therefore chosen to test the HCSSM on filtering and background correction, as these are fields where there is still considerable debate in the literature [4, 17].
Effects of level of filtration on the number of differentially expressed genes
It is common in microarray analysis to reduce the impact of spots that are malformed or have intensities outside the linear range of the scanner. Such spots are either removed, or weighted down using some measure of spot quality. Filtering has been shown to introduce bias in microarray studies , but has also been shown to significantly reduce variation in self-self hybridizations . However, as far as we know, no work has shown how the bias versus variance trade-off in practice influences the ability to detect differentially expressed genes. Therefore, we tested a number of different filtering methods of increasing complexity to evaluate how they influence the ability to find differentially expressed genes in a study.
Often filtering consists of removing spots with intensities below a (arbitrary) threshold . We have however found no investigations into the effect of such filtering. Several filtering methods based on spot quality statistics have been reported [13, 18, 20]. These quality measures can be used as filters by setting a spot quality threshold, but can also be used as weights to reduce the impact of low quality spots.
Where G and R are the green and red signal intensities, σ R , σ G are the pixel standard deviations in the red and green channel respectively.
The delta term is a small number added to stabilize spots with very small uncertainties.
Effects of background correction on the number of differentially expressed genes
The rationale of background correction is that the e.g. observed green signal intensity (G) is a linear combination of a true signal coming from the labelled RNA (G t ) and a background signal (G b ) so that:
G = G t + G b
More accurate log ratios can then be obtained by correcting the observed signal by subtracting the background. The problem then is to correctly estimate the background signal for each spot in the image analysis phase, and correcting this estimate to satisfy prior criteria e.g. background smoothness and non-negativity of the true signal. Background correction in itself may also destabilize results by increasing variance , and may decrease the ability to detect differentially expressed genes. Therefore, we use the HCSSM to compare the ability to find differentially expressed genes when no background correction is used versus a non negative background correction method .
Effect of choice and dose of dT blockers in cDNA microarray experiments
The HCSS-method can be applied to the laboratory procedures as well as to the data processing steps in microarray research. One major source of bias in cDNA microarrays is cross-hybridization of the labelled RNA to non-target homologous probe sequences on the array [22–24]. A substantial proportion of the non-specific cross-hybridization signal is due to poly(dA)-poly(dT) cross hybridization since poly(dT)-containing molecules produced during labelling of poly(dA) tails of sample RNAs by reverse transcription can bind promiscuously to the poly(dA) stretch of cDNA probe spots. These poly(dA) hybridization signals reduce the ability to detect differentially expressed genes. Thus, minimization of poly(dA) hybridization is a major challenge in cDNA microarray experiments. To reduce the poly(dA) signals it is common practice to add a blocker like synthetic poly(dA) in the probe or poly(dT) in the hybridization mixture [25–27]. Recently, an LNA dT Blocker containing Locked Nucleic Acid (LNA) nucleotides  has been introduced. There is currently no standardized protocol available that can be used to evaluate and optimize blocking procedures with respect to optimizing the observed number of differentially expressed genes. In the present study we therefore compared the effect of poly(dA)40–60 and different doses of LNA dT blocker on the number of differentially expressed genes, and used the HCSSM to optimize identification of truly differentially expressed genes.
Results and discussion
The choice and dose of blocking reagent influences the number of differentially expressed genes estimated
In order to evaluate how different conditions designed to block poly(dA)-poly(dT) cross hybridization affect the ability to detect differentially expressed genes, we applied the HCSS-method to analysis of microarray experiments performed in the presence of either poly(dA)40–60 or varying amounts of LNA. Poly(dA)40–60 binds to the poly(dT) segment of the labelled target in the hybridization solution and thus competes with the poly(dA) segment in the microarray cDNA probes. LNA dT blocker contains Locked Nucleic Acid (LNA) nucleotides  at key positions within the (dT) synthetic strand, and is designed to block poly (dA) sequences present in the microarray cDNA probes and prevent them from hybridizing to poly(dT) segment of the labelled targets.
Background correction method has a stronger impact on the estimated number of differentially expressed genes than filtering
Validation of cDNA microarray data generated by HCSSM
Selected GO categories associated with genes differentially expressed in the cell lines AR42J versus NRK52E.
response to stimulus
response to stress
Cross platform microarray analysis
Since RT-PCR is still considered the "golden standard" for gene expression measurements, we further used SYBR green-based quantitative real-time PCR to validate the relative gene expression for a few of the 1478 additional genes identified with optimal conditions. Commonly, a 2 fold change is reported as the cut-off below which microarray and qRT-PCR data begin to loose correlation, although recent reports observed significant correlation where genes exhibited 1.4 fold change or higher [34, 35]. We therefore selected seven genes in the fold ratio range 4.8 to 1.8 (log2 ratio range 2.2 to 0.9 on the cDNA arrays). Results are shown in Additional files 3 and 4. Description of PCR primers and qRT-PCR protocol are shown in Additional file 5. Five of the selected genes (Ica1, c-fos, Btg2, Uhrf and Ube2b) were detected as differentially expressed on both microarray platforms. Two genes (Hoxa2b and Irfrd) were only identified on the the cDNA platform. Five genes with fold ratio > 2 were verified by qRT-PCR (see Additional files 3 and 4). Two genes (Hoxa2 and Ube2b) with fold ratio < 2 were not verified. Hox2b was excluded from the results due to technical problems (see Additional file 4). The gene Ube2b was expressed in the opposite direction relative to the cDNA microarray, but in the same direction as on the Illumina platform (higher in AR42J than in NRK52E on Illumina, lower on cDNA). Consistent with other reports [35, 36], the microarray ratios were compressed compared to qRT-PCR results (see Additional file 3). Overall, the qRT-PCR measurement corresponded 100 % (5/5) with the Illumina data and 71 % (5/7) with the cDNA microarray data. The latter confirmation rate is in accordance with cDNA microarray results of others (, and references therein).
Some of the discrepancy between cDNA microarry versus Illumina and qRT-PCR may be due to errors (reported by e.g. ) in the IMAGE collection used for probes in our cDNA microarrays. We thus consider the results from qRT-PCR and the Illumina platform more likely to be "true data" in the present study. In spite of the known errors on cDNA platforms, cross platform analysis confirmed 72 % of the results. Validation of a few selected genes by qRT-PCR supported this finding. Taken together, functional annotation analysis and the alternative methods for validation of the cDNA micoarry data indicated that the increase in the number of differentially expressed genes found by HCSSM was mainly real.
In the present study we have used the HCSSM in cDNA microarray control experiments to evaluate blocking conditions, filtration and background correction. We found that the choice and dose of poly(dT)-poly(dA) blocking agents during hybridization influenced the estimated number of differentially expressed genes. When HCSSM was used to evaluate steps in the computational analysis of microarray data we found that no background correction and weight based filtering gave the largest estimated numbers of differentially expressed genes. The background correction method, however, has a stronger impact on the numbers of differentially expressed genes found than filtering. The results show that HCSSM may be a useful and simple approach to optimize cDNA microarrays procedures without including external standards.
AR42J (rat pancreatic acinar cell derived, ATCC, Rockville, MD) were grown in Dulbecco's modified Eagle's medium (DMEM) with 4.5 g/l glucose (Invitrogen, Carlsbad, CA), 15 % fetal calf serum (Euroclone Ltd, Devon, UK), 1 mM sodium pyruvate, 0.1 mg/ml L-glutamine (Invitrogen, Carlsbad, CA), 10 U/ml penicillin/streptomycin (Invitrogen, Carlsbad, CA), and 1 μ g/ml fungizone (Invitrogen Carlsbad, CA). NRK52E cells (rat kidney epithelial cells, ATCC, Rockville, MD) were grown in DMEM with 4.5 g/l glucose (Invitrogen, Carlsbad, CA) supplemented with 5 % FCS (Euroclone Ltd), 0.1 mg/ml L-glutamine (Invitrogen, Carlsbad, CA) and 10 U/ml penicillin/streptomycin (Invitrogen, Carlsbad, CA). All cells were grown at 37°C in a humidified atmosphere at 5 % CO2.
Treatment of cells and isolation of RNA
AR42J and NRK52E cells were cultured in 75 cm2 culture flasks for 72 h until confluence was reached. Total RNA was isolated using RNeasy® Midi kit (Qiagen, Germantown, MD) according to the manufacturer's instruction. RNA was quantified and assessed for purity by measurement of OD260 and OD280 using a UV fiberoptic spectrophotometer (Nanodrop Technologies, Rockland, DE) and was qualitatively assessed by measurement of relative 28S and 18S ribosomal band intensities using a Bioanalyzer and RNA LabChip capillary gel electrophoresis assay (Agilent Technologies, Palo Alto, CA). Aliquots were kept frozen at -80°C until further processing.
cDNA Microarray procedure
cDNA Microarrays were manufactured by the Norwegian Microarray Consortium  using 15000 cDNA rat probes from Research Genetics (IMAGE collection) printed in duplicates on Corning CMT Gaps II slides (Corning Inc. NY). The probes were dissolved in 50 % DMSO to ensure high printing quality. The microarray slides were UV cross linked at 300 mJ in order to fix the DNA to the glass slides.
Labelling was performed by reverse transcription of total RNA from the cell lines (3 μ g each) in the presence of primers containing the capture sequence for subsequent hybridization of Cy3- and Cy5-labelled dendrimers, using the Genisphere 3DNA Array 350 Expression Array Detection kit (Genisphere, Hatfield, PA) as described in the manufacturer's protocol. Different combinations of blocking reagents were added to the hybridization mixture: Mouse COT-1 DNA (Life Technologies) was added to all mixtures (0.1 μ g/μ g RNA) in order to block hybridization of repetitive elements. Different concentrations (0.1–4.0 μ l = 25–1000 ng) of LNA dT Blocker (Genisphere, Hatfield, PA) and 1 μ l = 1000 ng Poly(dA)40–60 (Stratagen Spotreporter; La Jolla, CA) were used as illustrated in Figure 2.
Hybridization was done in a humidified hybridization chamber (Corning Inc., NY) at 60°C for 14–15 h in a total hybridization volume of 60 μ l. Post hybridization washes were done for 15 min at 55°C with 2× SSC and 0.2 % SDS, for 10 min at room temperature with 2× SSC, and finally for 10 min at room temperature with 0.2× SSC. After the washing, the Cy3- and Cy5-labelled dendrimers were hybridized to the capture sequence at the reverse transcribed sample RNA at 60°C for 3 h as described in the manufacturer's protocol. After dendrimer hybridization, the slides were washed as described above and dried by centrifugation at 1500 rpm for 5 min.
Scanning and image analysis of cDNA microarrays
The slides were scanned at a resolution of 10 μ m by use of Packard Bioscience Scanarray Express HT scanner (Packard BioScience, Billerica, MA). A laser power of 100 % was used, and excitation of Cy3 and Cy5 was performed at a wavelength of 532 nm and 635 nm, respectively. Signals were detected by use of photomultiplier tubes (PMTs) with two channels. In order to maximizing the certainty of the weakest spots, the PMT voltage was adjusted to keep the background intensity between 200 and 400 (mean spot intensity) in each channel as suggested in . The GenePix 5.0 image analysis software (Axon Instruments, Inc., Union City, CA) was used for spot segmentation and intensity calculations. Spots and regions with high unspecific binding of dye or dust particles were manually flagged and excluded from the analysis.
Data analysis of cDNA microarrays
Modified T-scores and false discovery rate (FDR) values were calculated as described under Applications in the Background section for all combinations of filter, background corrections, and blocking treatments.
A progression of increasingly complex filter functions was evaluated.
Coarse filter: Spots flagged as bad or missing during the image analysis are removed.
Medium filter: Also removes spots with intensities below 200 and spots with more than 70 % saturation in both channels.
Fine filter: Also removes spots that are small, i.e. spots with diameters less than 20.
Uncertain filter: Removes spots flagged by GenePix 5.0 as "bad" or "missing" and spots with deltaM > 3.
The delta term is a small number added to stabilize spots with very small uncertainties.
Background correction methods
Three different background corrections strategies were applied:None were only foreground signal is used to calculate ratios, Edwards which is a subtraction method that ensures positive values  and dampened Edwards were a small number (e.g. 20) is added to the corrected signal to avoid extreme ratios in those cases where the background corrected intensity would be very low.
After filtering and background correction, all arrays were normalized using loess to remove intensity dependent variations in the ratios . Differentially expressed genes were found at the 5 % FDR level using the methodology described in the Background section.
Cross-platform microarray analysis
The Illumina Gene Expression system was used for cross-platform analysis of differentially expressed genes in AR42J cells versus NRK52E cells. RNA amplifications and hybridization were performed at the Finnish DNA Microarray Centre at Turku Centre of Biotechnology . Briefly, aliquots from the same RNA as in the cDNA microarray experiments were amplified with Ambion's Illumina® TotalPrep RNA Amplification kit (cat no AMIL1791) using 400 ng of total RNA as input material. The in vitro transcription (IVT) amplification that incorporated biotin-labelled nucleotides was performed overnight (14 h) at 37°C. After the amplifications the cRNA concentrations where checked with NanoDrop ND-1000 and cRNA quality was controlled by BioRad's Experion electrophoresis station.
A total of 750 ng of each biotin-labelled cRNA sample was hybridized to Illumina's Sentrix® RaRef-12-v1 Expression BeadChips at 58°C overnight (17 h) according to the manufacturer's instructions. The hybridized biotinylated cRNA was detected with 1 μ g/ml Cyanine3-streptavidine, (GE Healthcare Biosciences; cat no PA43001) and the Beadchips were scanned with Illumina BeadArray Reader (Factor = 1, PMT = 521, Filter = 100 %). Numerical results were extracted with Bead Studio v188.8.131.52 without any normalization or background subtraction.
The Illumina data were analysed using a modified T-test  with p-values adjusted for false discovery rate , and genes with adjusted p-values < 0.05 were taken as significant. Entrez Gene IDs were used as common identifiers for the microarray platforms. The genes were identified according to UniGen build #163 in GeneTools [29, 30].
Database submission of microarray data
The microarray data were prepared according to minimum information about a microarray experiment (MIAME) recommendations  and deposited in the Array Express . Detailed information about the microarray designs (platforms) and raw data files from the experiments are accessible in ArrayExpress by use of these accession numbers; A-MEXP-358 (cDNA platform), E-TABM-202 (cDNA experiment) and E-TABM-307 (Illumina experiment).
This work was supported by The Norwegian Research Council, The Norwegian Cancer Association, The Cancer Research Fund at St. Olav's Hospital, and the scientific collaboration between NTNU and St. Olav's University Hospital (Samarbeidsorganet). The cDNA microarray work was done with support from the national microarray technology platform (NMC), which is supported by the functional genomic program (FUGE). The Illumina work was performed at the Finnish DNA Microarray Centre at Turku Centre of Biotechnology .
We thank Bob Getts from Genisphere for helpful discussion on optimizing and evaluation of blocking conditions, Kristine Misund and Tonje S. Steigedal from NTNU for designing qRT-PCR primers and Anne Kristensen from NTNU for skilful technical assistance.
- Sandvik AK, Alsberg BK, Norsett KG, Yadetie F, Waldum HL, Laegreid A: Gene expression analysis and clinical diagnosis. Clin Chim Acta. 2006, 363 (1-2): 157-164. 10.1016/j.cccn.2005.05.046.PubMedView ArticleGoogle Scholar
- Ge H, Walhout AJM, Vidal M: Integrating `omic' information: a bridge between genomics and systems biology. Trends in Genetics. 2003, 19 (10): 551-560. 10.1016/j.tig.2003.08.009.PubMedView ArticleGoogle Scholar
- Kitano H: Systems biology: a brief overview. Science. 2002, 295 (5560): 1662-1664. 10.1126/science.1069492.PubMedView ArticleGoogle Scholar
- Allison DB, Cui X, Page GP, Sabripour M: Microarray data analysis: from disarray to consolidation and consensus. Nat Rev Genet. 2006, 7 (1): 55-65. 10.1038/nrg1749.PubMedView ArticleGoogle Scholar
- Andersen MT, Foy CA: The development of microarray standards. Anal Bioanal Chem. 2005, 381 (1): 87-89. 10.1007/s00216-004-2825-5.PubMedView ArticleGoogle Scholar
- Novoradovskaya N, Whitfield ML, Basehore LS, Novoradovsky A, Pesich R, Usary J, Karaca M, Wong WK, Aprelikova O, Fero M, Perou CM, Botstein D, Braman J: Universal Reference RNA as a standard for microarray experiments. BMC Genomics. 2004, 5 (1): 20-10.1186/1471-2164-5-20.PubMed CentralPubMedView ArticleGoogle Scholar
- van Bakel H, Holstege FC: In control: systematic assessment of microarray performance. EMBO Rep. 2004, 5 (10): 964-969. 10.1038/sj.embor.7400253.PubMed CentralPubMedView ArticleGoogle Scholar
- Kawasaki ES: The End of the Microarray Tower of Babel: Will Universal Standards Lead the Way?. J Biomol Tech. 2006, 17 (3): 200-206.PubMed CentralPubMedGoogle Scholar
- Brazma A, Krestyaninova M, Sarkans U: Standards for systems biology. Nat Rev Genet. 2006, 7 (8): 593-605. 10.1038/nrg1922.PubMedView ArticleGoogle Scholar
- Shi L, Reid LH, Jones WD, Shippy R, Warrington JA, Baker SC, Collins PJ, de Longueville F, Kawasaki ES, Lee KY, Luo Y, Sun YA, Willey JC, Setterquist RA, Fischer GM, Tong W, Dragan YP, Dix DJ, Frueh FW, Goodsaid FM, Herman D, Jensen RV, Johnson CD, Lobenhofer EK, Puri RK, Schrf U, Thierry-Mieg J, Wang C, Wilson M, Wolber PK, Zhang L, Amur S, Bao W, Barbacioru CC, Lucas AB, Bertholet V, Boysen C, Bromley B, Brown D, Brunner A, Canales R, Cao XM, Cebula TA, Chen JJ, Cheng J, Chu TM, Chudin E, Corson J, Corton JC, Croner LJ, Davies C, Davison TS, Delenstarr G, Deng X, Dorris D, Eklund AC, Fan XH, Fang H, Fulmer-Smentek S, Fuscoe JC, Gallagher K, Ge W, Guo L, Guo X, Hager J, Haje PK, Han J, Han T, Harbottle HC, Harris SC, Hatchwell E, Hauser CA, Hester S, Hong H, Hurban P, Jackson SA, Ji H, Knight CR, Kuo WP, LeClerc JE, Levy S, Li QZ, Liu C, Liu Y, Lombardi MJ, Ma Y, Magnuson SR, Maqsodi B, McDaniel T, Mei N, Myklebost O, Ning B, Novoradovskaya N, Orr MS, Osborn TW, Papallo A, Patterson TA, Perkins RG, Peters EH, Peterson R, Philips KL, Pine PS, Pusztai L, Qian F, Ren H, Rosen M, Rosenzweig BA, Samaha RR, Schena M, Schroth GP, Shchegrova S, Smith DD, Staedtler F, Su Z, Sun H, Szallasi Z, Tezak Z, Thierry-Mieg D, Thompson KL, Tikhonova I, Turpaz Y, Vallanat B, Van C, Walker SJ, Wang SJ, Wang Y, Wolfinger R, Wong A, Wu J, Xiao C, Xie Q, Xu J, Yang W, Zhang L, Zhong S, Zong Y, Slikker W: The MicroArray Quality Control (MAQC) project shows inter- and intraplatform reproducibility of gene expression measurements. Nat Biotechnol. 2006, 24 (9): 1151-1161. 10.1038/nbt1239.PubMedView ArticleGoogle Scholar
- Han T, Melvin CD, Shi L, Branham WS, Moland CL, Pine PS, Thompson KL, Fuscoe JC: Improvement in the Reproducibility and Accuracy of DNA Microarray Quantification by Optimizing Hybridization Conditions. BMC Bioinformatics. 2006, 7 Suppl 2: S17-10.1186/1471-2105-7-S2-S17.PubMedView ArticleGoogle Scholar
- Thompson KL, Rosenzweig BA, Pine PS, Retief J, Turpaz Y, Afshari CA, Hamadeh HK, Damore MA, Boedigheimer M, Blomme E, Ciurlionis R, Waring JF, Fuscoe JC, Paules R, Tucker CJ, Fare T, Coffey EM, He Y, Collins PJ, Jarnagin K, Fujimoto S, Ganter B, Kiser G, Kaysser-Kranich T, Sina J, Sistare FD: Use of a mixed tissue RNA design for performance assessments on multiple microarray formats. Nucleic Acids Res. 2005, 33 (22): e187-10.1093/nar/gni186.PubMed CentralPubMedView ArticleGoogle Scholar
- Hautaniemi S, Edgren H, Vesanen P, Wolf M, Jarvinen AK, Yli-Harja O, Astola J, Kallioniemi O, Monni O: A novel strategy for microarray quality control using Bayesian networks. Bioinformatics. 2003, 19 (16): 2031-2038. 10.1093/bioinformatics/btg275.PubMedView ArticleGoogle Scholar
- Smyth GK: Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol. 2004, 3 (Article3 ):
- Wu W, Xing EP, Myers C, Mian IS, Bissell MJ: Evaluation of normalization methods for cDNA microarray data by k-NN classification. BMC Bioinformatics. 2005, 6: 191-10.1186/1471-2105-6-191.PubMed CentralPubMedView ArticleGoogle Scholar
- Fujita A, Sato J, Rodrigues L, Ferreira C, Sogayar M: Evaluating different methods of microarray data normalization. BMC Bioinformatics. 2006, 7 (1): 469-10.1186/1471-2105-7-469.PubMed CentralPubMedView ArticleGoogle Scholar
- Qin LX, Kerr KF: Empirical evaluation of data transformations and ranking statistics for microarray analysis. Nucleic Acids Res. 2004, 32 (18): 5471-5479. 10.1093/nar/gkh866.PubMed CentralPubMedView ArticleGoogle Scholar
- Tom BD, Gilks WR, Brooke-Powell ET, Ajioka JW: Quality determination and the repair of poor quality spots in array experiments. BMC Bioinformatics. 2005, 6: 234-10.1186/1471-2105-6-234.PubMed CentralPubMedView ArticleGoogle Scholar
- Quackenbush J: Microarray data normalization and transformation. Nat Genet. 2002, 32 Suppl: 496-501. 10.1038/ng1032.PubMedView ArticleGoogle Scholar
- Bylesjo M, Eriksson D, Sjodin A, Sjostrom M, Jansson S, Antti H, Trygg J: MASQOT: a method for cDNA microarray spot quality control. BMC Bioinformatics. 2005, 6: 250-10.1186/1471-2105-6-250.PubMed CentralPubMedView ArticleGoogle Scholar
- Edwards D: Non-linear normalization and background correction in one-channel cDNA microarray studies. Bioinformatics. 2003, 19 (7): 825-833. 10.1093/bioinformatics/btg083.PubMedView ArticleGoogle Scholar
- Flikka K, Yadetie F, Laegreid A, Jonassen I: XHM: a system for detection of potential cross hybridizations in DNA microarrays. BMC Bioinformatics. 2004, 5: 117-10.1186/1471-2105-5-117.PubMed CentralPubMedView ArticleGoogle Scholar
- Chuaqui RF, Bonner RF, Best CJ, Gillespie JW, Flaig MJ, Hewitt SM, Phillips JL, Krizman DB, Tangrea MA, Ahram M, Linehan WM, Knezevic V, Emmert-Buck MR: Post-analysis follow-up and validation of microarray experiments. Nat Genet. 2002, 32 Suppl: 509-514. 10.1038/ng1034.PubMedView ArticleGoogle Scholar
- Chen Y, Chou CC, Lu X, Slate E, Peck K, Xu W, Voit E, Almeida J: A multivariate prediction model for microarray cross-hybridization. BMC Bioinformatics. 2006, 7 (1): 101-10.1186/1471-2105-7-101.PubMed CentralPubMedView ArticleGoogle Scholar
- Handley D, Serban N, Peters D, O'Doherty R, Field M, Wasserman L, Spirtes P, Scheines R, Glymour C: Evidence of systematic expressed sequence tag IMAGE clone cross-hybridization on cDNA microarrays. Genomics. 2004, 83 (6): 1169-1175. 10.1016/j.ygeno.2003.12.010.PubMedView ArticleGoogle Scholar
- Nguyen C, Rocha D, Granjeaud S, Baldit M, Bernard K, Naquet P, Jordan BR: Differential gene expression in the murine thymus assayed by quantitative hybridization of arrayed cDNA clones. Genomics. 1995, 29 (1): 207-216. 10.1006/geno.1995.1233.PubMedView ArticleGoogle Scholar
- Pan SJ, Rigney DR, Ivy JL: Outliers involving the Poly(A) effect among highly-expressed genes in microarrays. BMC Genomics. 2002, 3 (1): 35-10.1186/1471-2164-3-35.PubMed CentralPubMedView ArticleGoogle Scholar
- Vester B, Wengel J: LNA (locked nucleic acid): high-affinity targeting of complementary RNA and DNA. Biochemistry. 2004, 43 (42): 13233-13241. 10.1021/bi0485732.PubMedView ArticleGoogle Scholar
- GeneTools - a collection of resources for functional annotations . [http://www.genetools.microarray.ntnu.no/common/intro.php]
- Beisvag V, Junge FK, Bergum H, Jolsum L, Lydersen S, Gunther CC, Ramampiaro H, Langaas M, Sandvik AK, Laegreid A: GeneTools--application for functional annotation and statistical hypothesis testing. BMC Bioinformatics. 2006, 7: 470-10.1186/1471-2105-7-470.PubMed CentralPubMedView ArticleGoogle Scholar
- Cheadle C, Becker KG, Cho-Chung YS, Nesterova M, Watkins T, Wood W, Prabhu V, Barnes KC: A rapid method for microarray cross platform comparisons using gene expression signatures. Mol Cell Probes. 2007, 21 (1): 35-46. 10.1016/j.mcp.2006.07.004.PubMedView ArticleGoogle Scholar
- Larkin JE, Frank BC, Gavras H, Sultana R, Quackenbush J: Independence and reproducibility across microarray platforms. Nat Methods. 2005, 2 (5): 337-344. 10.1038/nmeth757.PubMedView ArticleGoogle Scholar
- Bosotti R, Locatelli G, Healy S, Scacheri E, Sartori L, Mercurio C, Calogero R, Isacchi A: Cross platform microarray analysis for robust identification of differentially expressed genes. BMC Bioinformatics. 2007, 8 Suppl 1: S5-10.1186/1471-2105-8-S1-S5.PubMedView ArticleGoogle Scholar
- Morey JS, Ryan JC, Van Dolah FM: Microarray validation: factors influencing correlation between oligonucleotide microarrays and real-time PCR. Biol Proced Online. 2006, 8: 175-193. 10.1251/bpo126.PubMed CentralPubMedView ArticleGoogle Scholar
- Wang Y, Barbacioru C, Hyland F, Xiao W, Hunkapiller KL, Blake J, Chan F, Gonzalez C, Zhang L, Samaha RR: Large scale real-time PCR validation on gene expression measurements from two commercial long-oligonucleotide microarrays. BMC Genomics. 2006, 7: 59-10.1186/1471-2164-7-59.PubMed CentralPubMedView ArticleGoogle Scholar
- Dallas PB, Gottardo NG, Firth MJ, Beesley AH, Hoffmann K, Terry PA, Freitas JR, Boag JM, Cummings AJ, Kees UR: Gene expression levels assessed by oligonucleotide microarray analysis and quantitative real-time RT-PCR -- how well do they correlate?. BMC Genomics. 2005, 6 (1): 59-10.1186/1471-2164-6-59.PubMed CentralPubMedView ArticleGoogle Scholar
- Hofsli E, Thommesen L, Yadetie F, Langaas M, Kusnierczyk W, Falkmer U, Sandvik AK, Laegreid A: Identification of novel growth factor-responsive genes in neuroendocrine gastrointestinal tumour cells. Br J Cancer. 2005, 92 (8): 1506-1516. 10.1038/sj.bjc.6602535.PubMed CentralPubMedView ArticleGoogle Scholar
- Halgren RG, Fielden MR, Fong CJ, Zacharewski TR: Assessment of clone identity and sequence fidelity for 1189 IMAGE cDNA clones. Nucleic Acids Res. 2001, 29 (2): 582-588. 10.1093/nar/29.2.582.PubMed CentralPubMedView ArticleGoogle Scholar
- The Norwegian Microarray Consortium. [http://www.microarray.no/]
- Lyng H, Badiee A, Svendsrud DH, Hovig E, Myklebost O, Stokke T: Profound influence of microarray scanner characteristics on gene expression ratios: analysis and procedure for correction. BMC Genomics. 2004, 5 (1): 10-10.1186/1471-2164-5-10.PubMed CentralPubMedView ArticleGoogle Scholar
- Yang YH, Dudoit S, Luu P, Lin DM, Peng V, Ngai J, Speed TP: Normalization for cDNA microarray data: a robust composite method addressing single and multiple slide systematic variation. Nucl Acids Res. 2002, 30 (4): e15-10.1093/nar/30.4.e15.PubMed CentralPubMedView ArticleGoogle Scholar
- Finnish DNA Microarray Centre. [http://microarrays.btk.fi/]
- Benjamini Y: Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing. Journal of the Royal Statistical Society. 1995, 57 (1): 289-Google Scholar
- Brazma A, Hingamp P, Quackenbush J, Sherlock G, Spellman P, Stoeckert C, Aach J, Ansorge W, Ball CA, Causton HC, Gaasterland T, Glenisson P, Holstege FC, Kim IF, Markowitz V, Matese JC, Parkinson H, Robinson A, Sarkans U, Schulze-Kremer S, Stewart J, Taylor R, Vilo J, Vingron M: Minimum information about a microarray experiment (MIAME)-toward standards for microarray data. Nat Genet. 2001, 29 (4): 365-371. 10.1038/ng1201-365.PubMedView ArticleGoogle Scholar
- Parkinson H, Sarkans U, Shojatalab M, Abeygunawardena N, Contrino S, Coulson R, Farne A, Lara GG, Holloway E, Kapushesky M, Lilja P, Mukherjee G, Oezcimen A, Rayner T, Rocca-Serra P, Sharma A, Sansone S, Brazma A: ArrayExpress--a public repository for microarray gene expression data at the EBI. Nucleic Acids Res. 2005, 33 (Database issue): D553-5. 10.1093/nar/gki056.PubMed CentralPubMedView ArticleGoogle Scholar
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.