Correcting for intra-experiment variation in Illumina BeadChip data is necessary to generate robust gene-expression profiles

Abstract

BackgroundMicroarray technology is a popular means of producing whole genome transcriptional profiles, however high cost and scarcity of mRNA has led many studies to be conducted based on the analysis of single samples. We exploit the design of the Illumina platform, specifically multiple arrays on each chip, to evaluate intra-experiment technical variation using repeated hybridisations of universal human reference RNA (UHRR) and duplicate hybridisations of primary breast tumour samples from a clinical study.ResultsA clear batch-specific bias was detected in the measured expressions of both the UHRR and clinical samples. This bias was found to persist following standard microarray normalisation techniques. However, when mean-centering or empirical Bayes batch-correction methods (ComBat) were applied to the data, inter-batch variation in the UHRR and clinical samples were greatly reduced. Correlation between replicate UHRR samples improved by two orders of magnitude following batch-correction using ComBat (ranging from 0.9833-0.9991 to 0.9997-0.9999) and increased the consistency of the gene-lists from the duplicate clinical samples, from 11.6% in quantile normalised data to 66.4% in batch-corrected data. The use of UHRR as an inter-batch calibrator provided a small additional benefit when used in conjunction with ComBat, further increasing the agreement between the two gene-lists, up to 74.1%.ConclusionIn the interests of practicalities and cost, these results suggest that single samples can generate reliable data, but only after careful compensation for technical bias in the experiment. We recommend that investigators appreciate the propensity for such variation in the design stages of a microarray experiment and that the use of suitable correction methods become routine during the statistical analysis of the data.

Dexamethasone Stimulated Gene Expression in Peripheral Blood is a Sensitive Marker for Glucocorticoid Receptor Resistance in Depressed Patients

“…Pre-and post-dex samples from the training and the test cohort were processed within a single run each (that is, RNA amplification, hybridization and scanning). To further reduce the described batch effects for Illumina arrays (Kitchen et al, 2010) pre-and post-dex RNAs were hybridized to the same chip for each individual, and cases and controls were randomized across arrays and arrays positions.…”

Individual variation and longitudinal pattern of genome-wide DNA methylation from birth to the first two years of life

“…To eliminate potential batch effect, an empirical Bayes methodology implemented in ComBat software 41 was applied to all data from each of the different chips to effectively remove batch effects in the experiment. 42,49 Analysis of variance (ANOVA) was conducted to evaluate the batch effects for normalized β values and showed that about 96% of CpG sites had significant batch effects before normalization, but this proportion dropped to 0.12% after normalization using ComBat. All of the analyses were based on methylation levels after performing considering gender, two time points, CpG islands and gene structures.…”

A Regulatory T-Cell Gene Signature Is a Specific and Sensitive Biomarker to Identify Children With New-Onset Type 1 Diabetes

“…Before biomarker discovery analysis or hierarchical clustering, batch correction (ComBat [16]), positive control sum normalization, background subtraction (with counts ,1 changed to 1), sum normalization, and log 2 -transformation were performed (Bioconductor [17]). …”