Skip to main content

Advertisement

Download PDF

Volume 19 Supplement 5

Selected articles from the Biological Ontologies and Knowledge bases workshop 2017

A framework for analyzing DNA methylation data from Illumina Infinium HumanMethylation450 BeadChip

BMC Bioinformatics volume 19, Article number: 115 (2018) Cite this article

Abstract

Background

DNA methylation has been identified to be widely associated to complex diseases. Among biological platforms to profile DNA methylation in human, the Illumina Infinium HumanMethylation450 BeadChip (450K) has been accepted as one of the most efficient technologies. However, challenges exist in analysis of DNA methylation data generated by this technology due to widespread biases.

Results

Here we proposed a generalized framework for evaluating data analysis methods for Illumina 450K array. This framework considers the following steps towards a successful analysis: importing data, quality control, within-array normalization, correcting type bias, detecting differentially methylated probes or regions and biological interpretation.

Conclusions

We evaluated five methods using three real datasets, and proposed outperform methods for the Illumina 450K array data analysis. Minfi and methylumi are optimal choice when analyzing small dataset. BMIQ and RCP are proper to correcting type bias and the normalized result of them can be used to discover DMPs. R package missMethyl is suitable for GO term enrichment analysis and biological interpretation.

Background

DNA methylation is an important epigenetic modification which has shown numerous associations with biological processes and complex diseases such as diabetes, schizophrenia and cancer [1–4]. However, the methylomic landscape in disease pathogenesis has not yet been well characterized, especially in cancer where DNA methylation can be altered dramatically. Interests of exploring the associations between DNA methylation and complex diseases increase in disease studies.

Illumina Infinium HumanMethylation450 (450K) BeadChip array, which covers over 480K CpG sites and targets 96% of CpG islands in human genome [5], has been widely utilized in many large studies, such as The Cancer Genome Atlas (TCGA) and The International Cancer Genome Consortium (ICGC) Project [6]. With the availability of public data resources, a number of methods for analyzing the Illumina 450K array data became rapidly available in the past few years.

Unlike the previous platform Illumina Infinium HumanMethylation27 (27K) BeadChip, in which only one probe type is utilized, the Illumina 450K BeadChip includes two distinct probe types, Infinium I (n=135501) and Infinium II (n=350076) [5]. Each CpG site of Infinium I is targeted by two 50bp probes: one for detecting “methylated” (M) intensity and one for detecting “unmethylated” (U) intensity, whereas each CpG site of Infinium II uses just one probe to distinguish “M” and “U” intensity through different dye colors (green and red), then the β-value, indicating the methylation level of one CpG site, can be computed as β=M/(M+U+α) where α is 100 generally. M-value, M=log2(β/(1−β)), the logit-transformed β-value, is another quantity used in following up analysis.

Owing to the platform design, more loci could be tested simultaneously on a fixed array size. Bibikova et al. [5] report a difference between Infinium I and Infinium II that Infinium II assays demonstrate an average of β-value upward shift for U and downward shift for M (shown in Fig. 1). Dedeurwaerder et al. [7] further evaluated the Illumina 450K BeadChip and reported that the β-values obtained from Infinium II probes had a narrower dynamic range and were less reproducible than those obtained from Infinium I leading to a type design bias. Hence, data preprocessing and normalization is critical for analyzing the Illumina 450K array data. Although many methods and R packages for correcting probe design bias have been proposed [8–14], more attention should be paid on the entire framework for analyzing Illumina 450K array data and selecting different modules on various datasets. Here we present a generalized framework including importing data, quality control, within-array normalization, correction of the type bias, identification of differential methylated probes or regions and biological interpretation, whilst the most popular modules for each step will be introduced, facilitating users to select appropriate module according their needs.

Fig. 1
figure1

Density curves of the β-values of Infinium I and Infinium II of a HPV-HNC sample (GSM937820 in GSE38268). The distribution models of β-values of the two types’ probes are different, where the Infinium II probe’s curve shows a narrower dynamic range

Full size image

Methods

Figure 2 shows the diagram of the framework and the detail of each step will be introduced below.

Fig. 2
figure2

A generalized framework of Illumina 450K array data analysis. There are six steps in the generalized framework leading to a successful methylation data analysis

Full size image

Importing data

There are mainly two forms of the Illumina 450K array data: i) raw (*.idat) data which is the direct output from Illumina iScan system and stores intensities for each probe, ii) *.txt data, which is usually got after simple preprocessing and is easier to access. Both file formats can be handled by R packages: *.idat files can be read by illuminaio package [15] and *.txt files can be dealt with minfi [12], wateRmelon [11] et al. We use three datasets for expounding our framework.

Dataset 1: Illumina 450K dataset from Dedeurwaerder et al. (GEO accession number: GSE29290) [7]. Three samples on HCT116WT cell-lines are considered, to evaluate the capability of methods to reduce the replicate variation.

Dataset 2: Illumina 450K dataset of fresh frozen head and neck cancer (HNC) samples form GSE38268 [16], where three of them are HPV+ and other three are HPV- samples.

Dataset 3: Illumina 450K dataset of level 1 methylation data of TCGA KIRC samples, which are *.idat files containing 160 normal samples and 325 tumor samples, to evaluate the efficiency of different modules.

Quality control

After the data imported into R, we would evaluate the quality of data. First, probes displaying a high detected p-value should be filtered out (e.g. > 0.01), while such probes have a β-value of “NA” in *.txt files. It is worth mentioned that different strategies of assigning missing β-value are applied in different modules. In minfi, the β-values are assigned with “NAs” when both M and U intensities are zero, but an additional criteria is that “NAs” will be assigned if either M or U intensities don’t fluoresce above background. Second, there are 850 inbuilt control probes on the 450K array, such as bisulfite conversion I, bisulfite conversion II, extension, hybridization and negative (n=613), which can be used to evaluate the other probes’ intensities. Samples that can’t pass this quality control are excluded in further analysis. Third, probes on Chromosome X or Y should be filtered out to eliminate the impact of sex on differential methylated analysis in many studies. Fourth, Price ME et al. [17] reported that there were 4.3% of Illumina 450K probes containing a known SNP at the targeted site. Such probes should cause problems in inter-sample analysis. Last but not the least, cross-reactive probes on the Illumina 450K array are identified depending on [18], which is particularly problematic because the β-value of such probes is more likely to represent a combination of multiple sites and not the level of initially targeted CpG sites.

Within-array normalization

This step includes background correction and color bias adjustment. Background correction methods have been developed by Triche et al. [19] such as Noob and Normexp, which are based on convolution models and use out-of-band (OOB) probes intensities to measure the background. The lumi R package provides two different methods for eliminating the background. The first one is based on the negative-control probes inbuilt the BeadChip and the second one estimates the background from the density modes of probes intensities. The second one would show bad performance when there are more than two density modes for some sample. The popular methylumi R package [20] proposed a color bias adjustment based on smooth quantile or shift-and-scaling normalization. Globally, these methods seem to improve the data quality in some cases [7].

Correction of the type bias

The type bias is the one that is most crucial to correct as it is the main source decreasing the data quality. There have been several efforts to develop methodologies to correct the probe type bias because of the differences between Infinium I and Infinium II. Because Infinium I probes are more stable and reproducible across different samples, most methods reduce the bias of Infinium II rather than Infinium I probes.

The first method is called peak-based correction (PBC) [7], which rescale the methylation values of Infinium II to the same modes for distribution of methylation values of Infinium I. But this method is sensitive to the shape of β-value density curves and is therefore less robust when the methylation density distribution does not exhibit well-defined peaks.

Touleimat and Tost [8] developed a method called Subset Quantile Normalization (SQN) based on an assumption that the β-values of CpGs form the same biological category should have the same density distribution. They found that the normalization result of using the “relation to CpG” annotation perfectly corrected the bias.

Subset-quantile Within Array Normalization (SWAN) [10] was developed based on the assumption that the β-values distribution should be the same when the probes have the same number of CpGs. But SWAN also alters Infinium I probe data, which increases Infinium I technical variation, and does not seem to improve the data quality when applied to some datasets [21].

Beta Mixture Quantile normalization (BMIQ) [9] method decomposes the density profile of Infinium I and Infinium II probes by fitting a beta-mixture model of three states: unmethylated (default β-value <0.25), hemimethylated (default 0.25 ≤β-value <0.75) and fully methylated (default β-value ≥0.75). Then it uses a quantile normalization to fit β-values distribution of Infinium II to the corresponding β-values distribution of Infinium I. This method does not depend on unceremonious choices of biological characteristics to be used to normalize data. Thus it seems more suitable than other methods. However, some points appear worse after BMIQ correction.

Another method, called Regression on Correlated Probes (RCP) [14], uses a quantile linear regression model of correlation between pairs of nearby Infinium I and II probes that share the same genomic context to adjust the methylation levels of Infinium II probes. The weakness of RCP is that it may not fit some experimental data leading to a result worse than raw data.

While background is important for measuring absolute methylation levels for single sample/condition experiments, we ignore the background here in this analysis since it can be cancelled out when comparing two conditions. As shown by other studies [22], widely used normalization process, which is based on the assumption that the majority of signals should not change across compared conditions, usually makes mistakes when it was applied to experiments where large portions of signals are differentially expressed. Thus, with a good quality control on the analyzed datasets, we didn’t choose to apply common normalization strategies.

It is also important to remove non-biological variation called batch effects existing between batches and samples. Such batch effects can influence on measurement of global level that could be partially removed through between-sample normalization using principal component analysis.

Identification of DMPs/DMRs

As we mentioned above, the main focus of many methylation studies has been on detecting differentially methylated probes (DMPs) or regions (DMRs) associated with a phenotype. The β-value is the default value for methylation measurement, allowing easy biological interpretation. Another type of value, M-value, is used to express the degree of methylation obtained with Infinium. Due to the heteroscedasticity of β-value, the variance of M-value across the methylation range is approximately constant, so the M-value has better statistical properties. The two types of value are used in different methods.

SQN simply considers a probe as DMP if the absolute value of the difference between β-value medians of paired samples is higher than 0.2:

$$\left|median\left({\beta^{N}_{1},\ldots,\beta^{N}_{n}}\right)-median\left({\beta^{T}_{1},\ldots,\beta^{T}_{n}}\right)\right|\ge0.2 $$

where \(\beta ^{N}_{i}\) and \(\beta ^{T}_{i}\) corresponding β-values of paired normal and tumor samples. The 0.2 threshold represents approximately a difference in methylation level of 20% which can be detected by the Infinium technology with 99% confidence [5]. Then the corresponding differentially methylated gene identities can be obtained from the list of DMPs.

Minfi offers a comprehensive package to analyze Illumina 450K array data, where candidate regions are determined for DMR analysis and locally weighted scatterplot smoothing (LOESS) is adopted to smooth the methylation differences between groups within each determined region. It also can find long-range alterations such as identified hypomethylated blocks [23] based on “open sea” probes. A empirical Bayes moderated t-test is used in limma [24] when sample sizes are less than 10, in which case M-values should be used as they rely much more on Gaussianity assumption [25].

Generally, DMRs are detected by applying various statistical techniques such as Fisher’s exact test [26, 27], t-test [27], Wilcoxon rank sum test [28] or different regression models [29–31].

Biological interpretation

There is a long list of significant CpGs to be interpreted after differential methylation analysis. Using the Infinium annotation file, Illumina 450K probes are classified according to their relations to CGIs and to the closest annotated gene. Regarding their relation to CGIs, the probes are classified into four categories: sites located inside a CGI, sites located in the CGI shores (0-2k bp), sites located in the CGI shelves (2-4k bp) and sites located in the “open sea”. As regards their relation to annotated genes, the sites are categorized as inside the promoter, inside the 5’-UTR region, inside the gene body and inside the 3’-UTR region. Then the significant DMPs can be marked with their related genes.

Performing gene set analysis is a popular way to understand the affected potential gene pathways. Although gene set analysis is well established in gene expression experimental, the research in methylation data is ongoing in different groups. In Illumina 450K array, the numbers of CpGs associated with each gene ranges largely from 1 to 1299 [32]. Genes with larger numbers of probes are more likely to have significant differentially methylated CpGs [33]. With the ontology and knowledgebase developing [34–40], researchers can easily annotate the genes containing DMPs or DMRs to ontology entries, which brings convenience for understanding the function of genes in the pathogenesis of diseases. Obviously, a phenotype is associated with several -omics data, such as mRNA expression and protein expression, which suggests researchers should utilize integrated analysis with multi-dimension data like TCGA project does [41, 42].

Results

Reduce the technical variation and type bias

We evaluated the different modules of correcting methods including SQN, SWAN, Dasen, BMIQ and RCP methods on Dataset 1, described in Table 1. The methylated and unmethylated intensities are imported into R environment and some results are displayed in following figures. For each method, we first plotted the density curves of β-value for three samples in Dataset 1. Then, we computed the standard deviation across the three replicates.

Table 1 Normalization methods for Illumina 450K array data
Full size table

As seen, SQN, Dasen and RCP can significantly reduce the technical variation (Fig. 3) and the same result can be seen base on standard deviation of replicates (Fig. 4). Because BMIQ and RCP do not change the β-values of Infinium I probes, the standard deviation of Infinium I probes of the two methods stays the same as the raw data. SWAN shows the least ability of reducing the variation among replicates.

Fig. 3
figure3

Density curves of β-value of replicates. Density curves of all probes after normalization using different methods are shown in this figure. The separate color represents replicates of the dataset. Raw data shows that there are differences among the replicates, while SQN, Dasen and RCP significantly reduce the technical variation

Full size image
Fig. 4
figure4

Box plot of standard deviation across replicates. The standard deviation between replicates after normalization shows that SQN, Dasen and RCP make the normalized β-values of replicates more similar than other two methods

Full size image

We also plotted the density curves of β-value of Infinium I/II probes for different methods. BMIQ, RCP and SQN show similar performance on the sample GSM815136 such that the Infinium I and Infinium II probes have more similar distribution modes with same local maximum values, while Dasen and SWAN underperformed others regarding removing the type bias (Fig. 5).

Fig. 5
figure5

Density curves of β-value of two types in sample GSM815136 for different methods. Density curves of β-values of Infinium I/II probes are plotted with different colors. BMIQ, RCP and SQN can correct the type bias efficiently on this sample

Full size image

Identify DMPs/DMRs by preprocessed results of different methods

We used the dmpFinder function in minfi R package to evaluate the result of identifying DMPs obtained by the five methods. We selected the probes with p-value of the results obtained after dmpFinder less than 0.05 as DMPs. The number of DMPs of five methods is listed in Table 2. Then the Venn diagram of the detected probes was plotted, shown in Fig. 6. As shown in Fig. 6, BMIQ and RCP methods show larger overlap with other methods, where the ratio of common probes is 43.41% (7370/16977 in BMIQ) and 39.42% (7370/18696 in RCP). Moreover, the most probes identified after BMIQ have intersections with others and only 311 probes (311/16977=1.83%) are identified uniquely. But it is also shown that there are differences among results of different methods, which maybe caused by the models they applying to normalization. Then the GO term enrichment analysis were analyzed with gometh function in missMethyl package. The GO terms of FDR ≤0.05 of five methods were got (see Table 2) and the Venn diagram of GO terms obtained after normalization is shown in Fig. 7. RCP got the highest rate of common (42/76=55.26%) GO terms among the five methods and SQN got the most number (n=199) of GO terms of which 33.67% (67/199) identified by only this method.

Fig. 6
figure6

Venn diagram of DMPs obtained by five methods. After filtering probes with p-value <0.01, we use dmpfinder in minfi to identify DMPs with p-value <0.05. BMIQ and RCP can lead to better results than other three methods

Full size image
Fig. 7
figure7

Venn diagram of GO terms obtained by five methods. There are 42 common GO terms identified by five methods. RCP has the highest rate of common terms (42/76=55.26%). SQN discovers most terms but 33.67% of them are identified by only this method

Full size image
Table 2 The number of obtained DMPs and GO terms by five methods
Full size table

Efficiency of large dataset analysis

We use Dataset3 to evaluate the performance of minfi and meffil [43] in term of importing large data. It took 48 min and less than 3 Gb memory to import data of 485 samples (970 *.idat files) using meffil package on a computer with 8 Gb of RAM and 4 processors while minfi could not run on the same computer. Minfi took two hours and ∼ 23 Gb memory on a server to import the same dataset.

Discussion

In this study we propose a generalized framework for Illumina 450K data analysis. We evaluate five methods for correcting type bias. Analysis of reducing technical replicates showed that different methods optimized different assessment criterion. The summary table of evaluation is shown in Table 3, where the score is set to 1 if there is significant change using the method, otherwise the score is set to 0.5. RCP gets the highest total score based on mentioned criteria.

Table 3 Summary of evaluation of five methods
Full size table

The SQN, Dasen and RCP methods could significantly improve the replicated data quality, while SWAN and BMIQ didn’t show improvement of the replicates. SWAN and Dasen didn’t remove type bias as other methods, which might be due to that models they applying cannot fit the distribution of Infinium I well as other methods.

When evaluation focused on detecting DMPs, BMIQ and RCP got more overlapped DMPs and credible GO terms than other methods. It is should pointed out that result may vary largely when using different datasets, which will be validated in further work.

Illumina MethylationEPIC BeadChip [44] microarrays have been used in some project, which contain more probes on a single array. More efficient tools are in urgent need of merging 450K and EPIC array data and the efficiency of analysis should be considered. Minfi has been utilized widely but it cannot handle large dataset on personal computer in our view while the newly package meffil displayed surprising performance.

During the evaluating processing, there were conflicts between R packages, for example, the MethylSet object in wateRmelon and minfi are different because the one in minfi has been updated in the newest version while the wateRmelon still use the previous object constructor. It is should be noticed in case of using different version of them.

Conclusions

It is suggested that the Illumina 450K users should choose proper strategy about importing data, background eliminating, correcting dye bias, correcting the type bias and detecting DMPs or DMRs. When analyzing small dataset, minfi and methylumi are optimal choice to import data and SQN, BMIQ and RCP may be proper to correcting the Infinium I/II bias. R package missMethyl is suitable for GO term enrichment analysis and biological interpretation. In our view, minfi is a proper R package to import data, eliminate background and ENmix package can be used to correct the type bias, then the normalized data should be used in the remaining steps of the framework.

References

  1. 1

    Koga Y, Pelizzola M, Cheng E, Krauthammer M, Sznol M, Ariyan S, Narayan D, Molinaro AM, Halaban R, Weissman SM. Genome-wide screen of promoter methylation identifies novel markers in melanoma. Genome Res. 2009; 19(8):1462–70.

  2. 2

    Teng M, Balch C, Liu Y, Li M, Huang TH, Wang Y, Nephew KP, Li L. The influence of cis-regulatory elements on dna methylation fidelity. PloS ONE. 2012; 7(3):32928.

  3. 3

    Esteller M. Cancer epigenomics: Dna methylomes and histone-modification maps. Nat Rev Genet. 2007; 8(4):286.

  4. 4

    Irizarry RA, Ladd-Acosta C, Wen B, Wu Z, Montano C, Onyango P, Cui H, Gabo K, Rongione M, Webster M, et al. The human colon cancer methylome shows similar hypo-and hypermethylation at conserved tissue-specific cpg island shores. Nat Genet. 2009; 41(2):178–86.

  5. 5

    Bibikova M, Barnes B, Tsan C, Ho V, Klotzle B, Le JM, Delano D, Zhang L, Schroth GP, Gunderson KL, et al. High density dna methylation array with single cpg site resolution. Genomics. 2011; 98(4):288–95.

  6. 6

    Consortium ICG, et al. International network of cancer genome projects. Nature. 2010; 464(7291):993.

  7. 7

    Dedeurwaerder S, Defrance M, Calonne E, Denis H, Sotiriou C, Fuks F. Evaluation of the infinium methylation 450k technology. Epigenomics. 2011; 3(6):771–84.

  8. 8

    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(3):325–41.

  9. 9

    Teschendorff AE, Marabita F, Lechner M, Bartlett T, Tegner J, Gomez-Cabrero D, Beck S. A beta-mixture quantile normalization method for correcting probe design bias in illumina infinium 450 k dna methylation data. Bioinformatics. 2012; 29(2):189–96.

  10. 10

    Maksimovic J, Gordon L, Oshlack A. Swan: Subset-quantile within array normalization for illumina infinium humanmethylation450 beadchips. Genome Biol. 2012; 13(6):44.

  11. 11

    Pidsley R, Wong CC, Volta M, Lunnon K, Mill J, Schalkwyk LC. A data-driven approach to preprocessing illumina 450k methylation array data. BMC Genomics. 2013; 14(1):293.

  12. 12

    Aryee MJ, Jaffe AE, Corrada-Bravo H, Ladd-Acosta C, Feinberg AP, Hansen KD, Irizarry RA. Minfi: a flexible and comprehensive bioconductor package for the analysis of infinium dna methylation microarrays. Bioinformatics. 2014; 30(10):1363–9.

  13. 13

    Xu Z, Niu L, Li L, Taylor JA. Enmix: a novel background correction method for illumina humanmethylation450 beadchip. Nucleic Acids Res. 2015; 44(3):20.

  14. 14

    Niu L, Xu Z, Taylor JA. Rcp: a novel probe design bias correction method for illumina methylation beadchip. Bioinformatics. 2016; 32(17):2659–63.

  15. 15

    Smith ML, Baggerly KA, Bengtsson H, Ritchie ME, Hansen KD. illuminaio: An open source idat parsing tool for illumina microarrays. F1000Research. 2013; 2:264.

  16. 16

    Lechner M, Fenton T, West J, Wilson G, Feber A, Henderson S, Thirlwell C, Dibra HK, Jay A, Butcher L, et al. Identification and functional validation of hpv-mediated hypermethylation in head and neck squamous cell carcinoma. Genome Med. 2013; 5(2):15.

  17. 17

    Price EM, Cotton AM, Lam LL, Farré P., Emberly E, Brown CJ, Robinson WP, Kobor MS. Additional annotation enhances potential for biologically-relevant analysis of the illumina infinium humanmethylation450 beadchip array. Epigenetics Chromatin. 2013; 6(1):4.

  18. 18

    Chen Y-a, Lemire M, Choufani S, Butcher DT, Grafodatskaya D, Zanke BW, Gallinger S, Hudson TJ, Weksberg R. Discovery of cross-reactive probes and polymorphic cpgs in the illumina infinium humanmethylation450 microarray. Epigenetics. 2013; 8(2):203–9.

  19. 19

    Triche Jr TJ, Weisenberger DJ, Van Den Berg D, Laird PW, Siegmund KD. Low-level processing of illumina infinium dna methylation beadarrays. Nucleic Acids Res. 2013; 41(7):90.

  20. 20

    Davis S, Du P, Bilke S, Triche T, Bootwalla M. methylumi: Handle illumina methylation data. R Package version 2.0. 2014.

  21. 21

    Dedeurwaerder S, Defrance M, Bizet M, Calonne E, Bontempi G, Fuks F. A comprehensive overview of infinium humanmethylation450 data processing. Brief Bioinforma. 2013; 15(6):929–41.

  22. 22

    Robinson MD, Oshlack A. A scaling normalization method for differential expression analysis of rna-seq data. Genome Biol. 2010; 11(3):25.

  23. 23

    Berman BP, Weisenberger DJ, Aman JF, Hinoue T, Ramjan Z, Liu Y, Noushmehr H, Lange CP, van Dijk CM, Tollenaar RA, et al. Regions of focal dna hypermethylation and long-range hypomethylation in colorectal cancer coincide with nuclear lamina-associated domains. Nat Genet. 2012; 44(1):40–6.

  24. 24

    Smyth GK. Limma: linear models for microarray data. In: Bioinformatics and Computational Biology Solutions Using R and Bioconductor. New York: Springer: 2005. p. 397–420.

  25. 25

    Zhuang J, Widschwendter M, Teschendorff AE. A comparison of feature selection and classification methods in dna methylation studies using the illumina infinium platform. BMC Bioinformatics. 2012; 13(1):59.

  26. 26

    Stockwell PA, Chatterjee A, Rodger EJ, Morison IM. Dmap: differential methylation analysis package for rrbs and wgbs data. Bioinformatics. 2014; 30(13):1814–22.

  27. 27

    Warden CD, Lee H, Tompkins JD, Li X, Wang C, Riggs AD, Yu H, Jove R, Yuan Y-C. Cohcap: an integrative genomic pipeline for single-nucleotide resolution dna methylation analysis. Nucleic Acids Res. 2013; 41(11):117.

  28. 28

    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. 2012; 28(5):729–30.

  29. 29

    Hansen KD, Langmead B, Irizarry RA. Bsmooth: from whole genome bisulfite sequencing reads to differentially methylated regions. Genome Biol. 2012; 13(10):83.

  30. 30

    Park Y, Figueroa ME, Rozek LS, Sartor MA. Methylsig: a whole genome dna methylation analysis pipeline. Bioinformatics. 2014; 30(17):2414–22.

  31. 31

    Hebestreit K, Dugas M, Klein H-U. Detection of significantly differentially methylated regions in targeted bisulfite sequencing data. Bioinformatics. 2013; 29(13):1647–53.

  32. 32

    Phipson B, Maksimovic J, Oshlack A. missmethyl: an r package for analyzing data from illumina’s humanmethylation450 platform. Bioinformatics. 2015; 32(2):286–8.

  33. 33

    Geeleher P, Hartnett L, Egan LJ, Golden A, Raja Ali RA, Seoighe C. Gene-set analysis is severely biased when applied to genome-wide methylation data. Bioinformatics. 2013; 29(15):1851–7.

  34. 34

    Cheng L, Jiang Y, Wang Z, Shi H, Sun J, Yang H, Zhang S, Hu Y, Zhou M. Dissim: an online system for exploring significant similar diseases and exhibiting potential therapeutic drugs. Sci Rep. 2016; 6:30024.

  35. 35

    Cheng L, Sun J, Xu W, Dong L, Hu Y, Zhou M. Oahg: an integrated resource for annotating human genes with multi-level ontologies. Sci Rep. 2016; 6:34820.

  36. 36

    Peng J, Lu J, Shang X, Chen J. Identifying consistent disease subnetworks using dnet. Methods. 2017; 131:104–10.

  37. 37

    Cheng L, Yang H, Zhao H, Pei X, Shi H, Sun J, Zhang Y, Wang Z, Zhou M. Metsigdis: a manually curated resource for the metabolic signatures of diseases. Brief Bioinforma. 2017;:bbx103.

  38. 38

    Peng J, Wang H, Lu J, Hui W, Wang Y, Shang X. Identifying term relations cross different gene ontology categories. BMC Bioinforma. 2017; 18(16):573.

  39. 39

    Peng J, Xue H, Shao Y, Shang X, Wang Y, Chen J. A novel method to measure the semantic similarity of hpo terms. Int J Data Min Bioinforma. 2017; 17(2):173–88.

  40. 40

    Peng J, Zhang X, Hui W, Lu J, Li Q, Shang X. Improving the measurement of semantic similarity by combining gene ontology and co-functional network: a random walk based approach. BMC Syst Biol. 2018;:12(Suppl2). In press.

  41. 41

    Network CGAR, et al. Integrated genomic analyses of ovarian carcinoma. Nature. 2011; 474(7353):609.

  42. 42

    Network CGAR, et al. Integrated genomic characterization of endometrial carcinoma. Nature. 2013; 497(7447):67.

  43. 43

    Min J, Hemani G, Smith GD, Relton CL, Suderman M. Meffil: efficient normalisation and analysis of very large dna methylation samples. bioRxiv. 2017:125963.

  44. 44

    Pidsley R, Zotenko E, Peters TJ, Lawrence MG, Risbridger GP, Molloy P, Van Djik S, Muhlhausler B, Stirzaker C, Clark SJ. Critical evaluation of the illumina methylationepic beadchip microarray for whole-genome dna methylation profiling. Genome Biol. 2016; 17(1):208.

Download references

Acknowledgements

None declared.

Funding

Publication costs were funded by the Major State Research Development Program of China [No:2016YFC1202302], the fundamental research funds for the central universities (grant no. HIT.NSRIF.201652), the National Natural Science Foundation of China (No:61571152) and the National High-tech R&D Program of China (863 Program) [Nos:2015AA020101, 2015AA020108].

Availability of data and materials

Dataset 1 and Dataset 2 analyzed during the current study are available in the GEO repository, GSE29290, GSE38268. Dataset 3 is available in the GDC repository TCGA-KIRC

About this supplement

This article has been published as part of BMC Bioinformatics Volume 19 Supplement 5, 2018: Selected articles from the Biological Ontologies and Knowledge bases workshop 2017. The full contents of the supplement are available online at https://doi.org/https://bmcbioinformatics.biomedcentral.com/articles/supplements/volume-19-supplement-5.

Author information

Affiliations

  1. School of Computer Science and Technology, Harbin Institute of Technology, Harbin, 150001, China
    • Zhenxing Wang
    • , XiaoLiang Wu
    •  & Yadong Wang
Authors
  1. Zhenxing Wang
    View author publications
    You can also search for this author in
  2. XiaoLiang Wu
    View author publications
    You can also search for this author in
  3. Yadong Wang
    View author publications
    You can also search for this author in

Contributions

ZXW performed analyses and wrote the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Yadong Wang.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver(http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Wang, Z., Wu, X. & Wang, Y. A framework for analyzing DNA methylation data from Illumina Infinium HumanMethylation450 BeadChip. BMC Bioinformatics 19, 115 (2018). https://doi.org/10.1186/s12859-018-2096-3

Download citation

Keywords

  • DNA methylation
  • Illumina 450K
  • Normalization
  • Ontology interpretation

BMC Bioinformatics

ISSN: 1471-2105

Contact us