A data-driven approach to preprocessing Illumina 450K methylation array data

Springer Science and Business Media LLC - Tập 14 Số 1 - 2013
Ruth Pidsley1, Chloe Wong1, Manuela Volta1, Katie Lunnon1, Jonathan Mill1, Leonard C. Schalkwyk1
1Social, Genetic and Developmental Psychiatry,Institute of Psychiatry, King's College London, De Crespigny Park, London, UK

Tóm tắt

Abstract Background As the most stable and experimentally accessible epigenetic mark, DNA methylation is of great interest to the research community. The landscape of DNA methylation across tissues, through development and in disease pathogenesis is not yet well characterized. Thus there is a need for rapid and cost effective methods for assessing genome-wide levels of DNA methylation. The Illumina Infinium HumanMethylation450 (450K) BeadChip is a very useful addition to the available methods for DNA methylation analysis but its complex design, incorporating two different assay methods, requires careful consideration. Accordingly, several normalization schemes have been published. We have taken advantage of known DNA methylation patterns associated with genomic imprinting and X-chromosome inactivation (XCI), in addition to the performance of SNP genotyping assays present on the array, to derive three independent metrics which we use to test alternative schemes of correction and normalization. These metrics also have potential utility as quality scores for datasets. Results The standard index of DNA methylation at any specific CpG site is β = M/(M + U + 100) where M and U are methylated and unmethylated signal intensities, respectively. Betas (βs) calculated from raw signal intensities (the default GenomeStudio behavior) perform well, but using 11 methylomic datasets we demonstrate that quantile normalization methods produce marked improvement, even in highly consistent data, by all three metrics. The commonly used procedure of normalizing betas is inferior to the separate normalization of M and U, and it is also advantageous to normalize Type I and Type II assays separately. More elaborate manipulation of quantiles proves to be counterproductive. Conclusions Careful selection of preprocessing steps can minimize variance and thus improve statistical power, especially for the detection of the small absolute DNA methylation changes likely associated with complex disease phenotypes. For the convenience of the research community we have created a user-friendly R software package called wateRmelon, downloadable from bioConductor, compatible with the existing methylumi, minfi and IMA packages, that allows others to utilize the same normalization methods and data quality tests on 450K data.

Từ khóa


Tài liệu tham khảo

Bibikova M, Le J, Barnes B: Genome-wide DNA methylation profiling using Infinium assay. Epigenomics. 2009, 1: 177-200. 10.2217/epi.09.14. [http://www.ingentaconnect.com/content/fm/epi/2009/00000001/00000001/art00019]

Bibikova M, Barnes B, Tsan C, Ho V, Klotzle B, Le JM, Delano D, Zhang L, Schroth GP, Gunderson KL, Fan JB, Shen R: High density DNA methylation array with single CpG site resolution. Genomics. 2011, 98 (4): 288-295. 10.1016/j.ygeno.2011.07.007. [http://www.ncbi.nlm.nih.gov/pubmed/21839163]

Sandoval J, Heyn Ha, Moran S, Serra-Musach J, Pujana MA, Bibikova M, Esteller M: Validation of a DNA methylation microarray for 450,000 CpG sites in the human genome. Epigenetics. 2011, 6 (6): 692-702. 10.4161/epi.6.6.16196. [http://www.landesbioscience.com/journals/epigenetics/article/16196/]

Du P, Zhang X, Huang CC, Jafari N, Kibbe WA, Hou L, Lin SM: Comparison of Beta-value and M-value methods for quantifying methylation levels by microarray analysis. BMC Bioinformatics. 2010, 11: 587-10.1186/1471-2105-11-587.

Dedeurwaerder S, Defrance M, Calonne E, Sotiriou C, Fuks F: Evaluation of the Infinium Methylation 450K technology. Epigenetics. 2011, 3 (6): 771-784.

Maksimovic J, Gordon L, Oshlack A: SWAN: Subset quantile within-array normalization for Illumina Infinium HumanMethylation450 BeadChips. Genome Biol. 2012, 13 (6): R44-10.1186/gb-2012-13-6-r44. [http://www.ncbi.nlm.nih.gov/pubmed/22703947]

Touleimat N, Tost J: Complete pipeline for Infinium ®, Human Methylation 450K BeadChip data processing using subset quantile normalization for accurate DNA methylation estimation. Epigenomics. 2012, 4: 325-341. 10.2217/epi.12.21. [http://www.futuremedicine.com/doi/abs/10.2217/epi.12.21]

Carén H, Djos A, Nethander M, Sjöberg RM, Kogner P, Enström C, Nilsson S, Martinsson T: Identification of epigenetically regulated genes that predict patient outcome in neuroblastoma. BMC Cancer. 2011, 11: 66-10.1186/1471-2407-11-66. [http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=3045360&tool=pmcentrez&rendertype=abstract]

Dempster EL, Pidsley R, Schalkwyk LC, Owens S, Georgiades A, Kane F, Kalidindi S, Picchioni M, Kravariti E, Toulopoulou T, Murray RM, Mill J: Disease-associated epigenetic changes in monozygotic twins discordant for schizophrenia and bipolar disorder. Human Mol Genet. 2011, 20 (24): 4786-4796. 10.1093/hmg/ddr416.

Kinoshita M, Numata S, Tajima A, Shimodera S, Ono S, Imamura A, Iga JI, Watanabe S, Kikuchi K, Kubo H, Nakataki M, Sumitani S, Imoto I, Okazaki Y, Ohmori T: DNA methylation signatures of peripheral leukocytes in Schizophrenia. Neuromolecular Med. 2012, [http://www.ncbi.nlm.nih.gov/pubmed/22961555]

Rakyan VK, Beyan H, Down Ta, Hawa MI, Maslau S, Aden D, Daunay A, Busato F, Mein Ca, Manfras B, Dias KRM, Bell CG, Tost J, Boehm BO, Beck S, Leslie RD: Identification of type 1 diabetes-associated DNA methylation variable positions that precede disease diagnosis. PLoS Genet. 2011, 7 (9): e1002300-10.1371/journal.pgen.1002300. [http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=3183089&tool=pmcentrez&rendertype=abstract]

Sun Z, Chai HS, Wu Y, White WM, Donkena KV, Klein CJ, Garovic VD, Therneau TM, Kocher JPa: Batch effect correction for genome-wide methylation data with Illumina Infinium platform. BMC Med Genomics. 2011, 4: 84-10.1186/1755-8794-4-84. [http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=3265417&tool=pmcentrez&rendertype=abstract]

Irizarry Ra, Hobbs B, Collin F, Beazer-Barclay YD, Antonellis KJ, Scherf U, Speed TP: Exploration, normalization, and summaries of high density oligonucleotide array probe level data. Biostatistics (Oxford, England). 2003, 4 (2): 249-264. 10.1093/biostatistics/4.2.249. [http://www.ncbi.nlm.nih.gov/pubmed/12925520]

Roessler J, Ammerpohl O, Gutwein J, Hasemeier B, Anwar SL, Kreipe H, Lehmann U: Quantitative cross-validation and content analysis of the 450k DNA methylation array from Illumina, Inc. BMC Res Notes. 2012, 5: 210-10.1186/1756-0500-5-210. [http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=3420245&tool=pmcentrez&rendertype=abstract]

Schulz R, Woodfine K, Menheniott TR, Bourc’his D, Bestor T, Oakey RJ: WAMIDEX: a web atlas of murine genomic imprinting and differential expression. Epigenetics. 2008, 3 (2): 89-96. 10.4161/epi.3.2.5900. [http://www.landesbioscience.com/journals/epi/article/5900/]

Schulz R: Wamidex, accessed 26 Jan2012. [https://atlas.genetics.kcl.ac.uk/]

Chen Ya, Choufani S, Ferreira JC, Grafodatskaya D, Butcher DT, Weksberg R: Sequence overlap between autosomal and sex-linked probes on the Illumina HumanMethylation27 microarray. Genomics. 2011, 97 (4): 214-222. 10.1016/j.ygeno.2010.12.004. [http://www.ncbi.nlm.nih.gov/pubmed/21211562]

Davis S, Du P, Bilke S, Triche JrT, Bootwalla M: Methylumi: Handle Illumina Methylation Data. 2012, R package version 2.2.0, .

Hansen KD, Aryee M: Minfi: Analyze Illumina’s 450k Methylation Arrays. R package version 1.2.0

Wang D, Yan L, Hu Q, Sucheston LE, Higgins MJ, Ambrosone CB, Johnson CS, Smiraglia DJ, Liu S: IMA: an R package for high-throughput analysis of Illumina’s 450K Infinium methylation data. Bioinformatics (Oxford, England). 2012, 28 (5): 729-730. 10.1093/bioinformatics/bts013. [http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=3289916&tool=pmcentrez&rendertype=abstract]

Teschendorff AE, Marabita F, Lechner M, Bartlett T, Tegner J, Gomez-Cabrero D, Beck S: A Beta-Mixture quantile normalisation method for correcting probe design bias in Illumina Infinium 450k DNA methylation data. Bioinformatics. 2012, 29 (2): 189-96. [http://www.ncbi.nlm.nih.gov/pmc/articles/PMC3546795/]

Heijmans BT, Mill J: Commentary: The seven plagues of epigenetic epidemiology. Int J Epidemiol. 2012, 41: 74-78. 10.1093/ije/dyr225. [http://www.ncbi.nlm.nih.gov/pubmed/22269254]

R Development Core Team: Vienna: R Foundation for Statistical Computing; 2012.: , ISBN 3-900051-07-0 [http://www.R-project.org/]

Gentleman RC, Carey VJ, Bates DM: Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004, 5 (R80): [http://genomebiology.com/2004/5/10/R80]

Smyth GK: Limma: linear models for microarray data. Bioinformatics and Computational Biology Solutions using R and Bioconductor. Edited by: Gentleman R CareyVDudoit S. 2005, New York: Springer, 397-420.

Sing T, Sander O, Beerenwinkel N, Lengauer T: 2009, R package version 1.0-4 [http://CRAN.R-project.org/package=ROCR.]

Thông tin
Thông tin xuất bản

Springer Science and Business Media LLC

Tập 14 Số 1

Thông tin tác giả