Batch effect correction for genome-wide methylation data with Illumina Infinium platform
Open Access
- 16 December 2011
- journal article
- Published by Springer Science and Business Media LLC in BMC Medical Genomics
- Vol. 4 (1), 84
- https://doi.org/10.1186/1755-8794-4-84
Abstract
Background Genome-wide methylation profiling has led to more comprehensive insights into gene regulation mechanisms and potential therapeutic targets. Illumina Human Methylation BeadChip is one of the most commonly used genome-wide methylation platforms. Similar to other microarray experiments, methylation data is susceptible to various technical artifacts, particularly batch effects. To date, little attention has been given to issues related to normalization and batch effect correction for this kind of data. Methods We evaluated three common normalization approaches and investigated their performance in batch effect removal using three datasets with different degrees of batch effects generated from HumanMethylation27 platform: quantile normalization at average β value (QNβ); two step quantile normalization at probe signals implemented in "lumi" package of R (lumi); and quantile normalization of A and B signal separately (ABnorm). Subsequent Empirical Bayes (EB) batch adjustment was also evaluated. Results Each normalization could remove a portion of batch effects and their effectiveness differed depending on the severity of batch effects in a dataset. For the dataset with minor batch effects (Dataset 1), normalization alone appeared adequate and "lumi" showed the best performance. However, all methods left substantial batch effects intact in the datasets with obvious batch effects and further correction was necessary. Without any correction, 50 and 66 percent of CpGs were associated with batch effects in Dataset 2 and 3, respectively. After QNβ, lumi or ABnorm, the number of CpGs associated with batch effects were reduced to 24, 32, and 26 percent for Dataset 2; and 37, 46, and 35 percent for Dataset 3, respectively. Additional EB correction effectively removed such remaining non-biological effects. More importantly, the two-step procedure almost tripled the numbers of CpGs associated with the outcome of interest for the two datasets. Conclusion Genome-wide methylation data from Infinium Methylation BeadChip can be susceptible to batch effects with profound impacts on downstream analyses and conclusions. Normalization can reduce part but not all batch effects. EB correction along with normalization is recommended for effective batch effect removal.Keywords
This publication has 33 references indexed in Scilit:
- DNA methylation profiling reveals novel biomarkers and important roles for DNA methyltransferases in prostate cancerGenome Research, 2011
- Impaired hydroxylation of 5-methylcytosine in myeloid cancers with mutant TET2Nature, 2010
- A comparison of batch effect removal methods for enhancement of prediction performance using MAQC-II microarray gene expression dataThe Pharmacogenomics Journal, 2010
- Genome-wide DNA methylation analysis in cohesin mutant human cell linesNucleic Acids Research, 2010
- Supervised normalization of microarraysBioinformatics, 2010
- Age-dependent DNA methylation of genes that are suppressed in stem cells is a hallmark of cancerGenome Research, 2010
- Human aging-associated DNA hypermethylation occurs preferentially at bivalent chromatin domainsGenome Research, 2010
- Epigenetic profiling of somatic tissues from human autopsy specimens identifies tissue- and individual-specific DNA methylation patternsHuman Molecular Genetics, 2009
- lumi: a pipeline for processing Illumina microarrayBioinformatics, 2008
- The power and the promise of DNA methylation markersNature Reviews Cancer, 2003