Statistical methods for classification of 5hmC levels based on the Illumina Inifinium HumanMethylation450 (450k) array data, under the paired bisulfite (BS) and oxidative bisulfite (oxBS) treatment.
Ontology highlight
ABSTRACT:
SUBMITTER: Slynko A
PROVIDER: S-EPMC6563990 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:DNA methylation analysis has become an integral part of biomedical research. For high-throughput applications such as epigenome-wide association studies, the Infinium HumanMethylation450 (450K) BeadChip is currently the platform of choice. However, BeadChip processing relies on traditional bisulfite (BS) based protocols which cannot discriminate between 5-methylcytosine (5mC) and 5-hydroxymethylcytosine (5hmC). Here, we report the adaptation of the recently developed oxidative bisulfite (oxBS) chemistry to specifically detect both 5mC and 5hmC in a single workflow using 450K BeadChips, termed oxBS-450K. Supported by validation using mass spectrometry and pyrosequencing, we demonstrate reproducible (R(2)>0.99) detection of 5hmC in human brain tissue using the optimised oxBS-450K protocol described here.
| S-EPMC4304834 | biostudies-literature
Project description:DNA methylation is the most widely studied epigenetic mark and is known to be essential to normal development and frequently disrupted in disease. The Illumina HumanMethylation450 BeadChip assays the methylation status of CpGs at 485,577 sites across the genome. Here we present Subset-quantile Within Array Normalization (SWAN), a new method that substantially improves the results from this platform by reducing technical variation within and between arrays. SWAN is available in the minfi Bioconductor package.
| S-EPMC3446316 | biostudies-literature
Project description:BackgroundAs 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.ResultsThe 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.ConclusionsCareful 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.
| S-EPMC3769145 | biostudies-literature
Project description:Motivation5-Methylcytosine (5mC) and 5-hydroxymethylcytosine (5hmC) are important epigenetic regulators of gene expression. 5mC and 5hmC levels can be computationally inferred at single base resolution using sequencing or array data from paired DNA samples that have undergone bisulfite and oxidative bisulfite conversion. Current estimation methods have been shown to produce irregular estimates of 5hmC level or are extremely computation intensive.ResultsWe developed an efficient method oxBS-MLE based on binomial modeling of paired bisulfite and oxidative bisulfite data from sequencing or array analysis. Evaluation in several datasets showed that it outperformed alternative methods in estimate accuracy and computation speed.Availability and implementationoxBS-MLE is implemented in Bioconductor package ENmix.Contactniulg@ucmail.uc.eduSupplementary information: Supplementary data are available at Bioinformatics online.
| S-EPMC5181539 | biostudies-literature
Project description:Due to their relatively low-cost per sample and broad, gene-centric coverage of CpGs across the human genome, Illumina's 450k arrays are widely used in large scale differential methylation studies. However, by their very nature, large studies are particularly susceptible to the effects of unwanted variation. The effects of unwanted variation have been extensively documented in gene expression array studies and numerous methods have been developed to mitigate these effects. However, there has been much less research focused on the appropriate methodology to use for accounting for unwanted variation in methylation array studies. Here we present a novel 2-stage approach using RUV-inverse in a differential methylation analysis of 450k data and show that it outperforms existing methods.
| S-EPMC4652745 | biostudies-literature
Project description:The Infinium 450K Methylation array is an established tool for measuring methylation. However, the bisulfite (BS) reaction commonly used with the 450K array cannot distinguish between 5-methylcytosine (5mC) and 5-hydroxymethylcytosine (5hmC). The oxidative-bisulfite assay disambiguates 5mC and 5hmC. We describe the use of oxBS in conjunction with the 450K array (oxBS-array) to analyse 5hmC/5mC in cerebellum DNA. The "methylation" level derived by the BS reaction is the combined level of 5mC and 5hmC at a given base, while the oxBS reaction gives the level of 5mC alone. The level of 5hmC is derived by subtracting the oxBS level from the BS level. Here we present an analysis method that distinguishes genuine positive levels of 5hmC at levels as low as 3%. We performed four replicates of the same sample of cerebellum and found a high level of reproducibility (average r for BS = 98.3, and average r for oxBS = 96.8). In total, 114,734 probes showed a significant positive measurement for 5hmC. The range at which we were able to distinguish 5hmC occupancy was between 3% and 42%. In order to investigate the effects of multiple replicates on 5hmC detection we also simulated fewer replicates and found that decreasing the number of replicates to two reduced the number of positive probes identified by > 50%. We validated our results using qPCR in conjunction with glucosylation of 5hmC sites followed by MspI digestion and we found good concordance with the array estimates (r = 0.94). This experiment provides a map of 5hmC in the cerebellum and a robust dataset for use as a standard in future 5hmC analyses. We also provide a novel method for validating the presence of 5hmC at low levels, and highlight some of the pitfalls associated with measuring 5hmC and 5mC.
| S-EPMC4338296 | biostudies-literature
Project description:BACKGROUND: Measurement of genome-wide DNA methylation (DNAm) has become an important avenue for investigating potential physiologically-relevant epigenetic changes. Illumina Infinium (Illumina, San Diego, CA, USA) is a commercially available microarray suite used to measure DNAm at many sites throughout the genome. However, it has been suggested that a subset of array probes may give misleading results due to issues related to probe design. To facilitate biologically significant data interpretation, we set out to enhance probe annotation of the newest Infinium array, the HumanMethylation450 BeadChip (450 k), with >485,000 probes covering 99% of Reference Sequence (RefSeq) genes (National Center for Biotechnology Information (NCBI), Bethesda, MD, USA). Annotation that was added or expanded on includes: 1) documented SNPs in the probe target, 2) probe binding specificity, 3) CpG classification of target sites and 4) gene feature classification of target sites. RESULTS: Probes with documented SNPs at the target CpG (4.3% of probes) were associated with increased within-tissue variation in DNAm. An example of a probe with a SNP at the target CpG demonstrated how sample genotype can confound the measurement of DNAm. Additionally, 8.6% of probes mapped to multiple locations in silico. Measurements from these non-specific probes likely represent a combination of DNAm from multiple genomic sites. The expanded biological annotation demonstrated that based on DNAm, grouping probes by an alternative high-density and intermediate-density CpG island classification provided a distinctive pattern of DNAm. Finally, variable enrichment for differentially methylated probes was noted across CpG classes and gene feature groups, dependant on the tissues that were compared. CONCLUSION: DNAm arrays offer a high-throughput approach for which careful consideration of probe content should be utilized to better understand the biological processes affected. Probes containing SNPs and non-specific probes may affect the assessment of DNAm using the 450 k array. Additionally, probe classification by CpG enrichment classes and to a lesser extent gene feature groups resulted in distinct patterns of DNAm. Thus, we recommend that compromised probes be removed from analyses and that the genomic context of DNAm is considered in studies deciphering the biological meaning of Illumina 450 k array data.
| S-EPMC3740789 | biostudies-literature
Project description:BACKGROUND:DNA methylation has been recognized as a key mechanism in cell differentiation. Various studies have compared tissues to characterize epigenetically regulated genomic regions, but due to differences in study design and focus there still is no consensus as to the annotation of genomic regions predominantly involved in tissue-specific methylation. We used a new algorithm to identify and annotate tissue-specific differentially methylated regions (tDMRs) from Illumina 450k chip data for four peripheral tissues (blood, saliva, buccal swabs and hair follicles) and six internal tissues (liver, muscle, pancreas, subcutaneous fat, omentum and spleen with matched blood samples). RESULTS:The majority of tDMRs, in both relative and absolute terms, occurred in CpG-poor regions. Further analysis revealed that these regions were associated with alternative transcription events (alternative first exons, mutually exclusive exons and cassette exons). Only a minority of tDMRs mapped to gene-body CpG islands (13%) or CpG islands shores (25%) suggesting a less prominent role for these regions than indicated previously. Implementation of ENCODE annotations showed enrichment of tDMRs in DNase hypersensitive sites and transcription factor binding sites. Despite the predominance of tissue differences, inter-individual differences in DNA methylation in internal tissues were correlated with those for blood for a subset of CpG sites in a locus- and tissue-specific manner. CONCLUSIONS:We conclude that tDMRs preferentially occur in CpG-poor regions and are associated with alternative transcription. Furthermore, our data suggest the utility of creating an atlas cataloguing variably methylated regions in internal tissues that correlate to DNA methylation measured in easy accessible peripheral tissues.
| S-EPMC3750594 | biostudies-literature
Project description:BACKGROUND: The newly released 450k DNA methylation array from Illumina, Inc. offers the possibility to analyze more than 480,000 individual CpG sites in a user friendly standardized format. In this study the relationship between the ?-values provided by the Illumina, Inc. array for each individual CpG dinucleotide and the quantitative methylation levels obtained by pyrosequencing were analyzed. In addition, the representation of microRNA genes and imprinted loci on the Illumina, Inc. array was assessed in detail. Genomic DNA from 4 human breast cancer cell lines (IPH-926, HCC1937, MDA-MB-134, PMC42) and 18 human breast cancer specimens as well as 4 normal mammary epithelial fractions was analyzed on 450k DNA methylation arrays. The ?-values for 692 individual CpG sites from 62 different genes were cross-validated using conventional quantitative pyrosequencing. FINDINGS: The newly released 450k methylation array from Illumina, Inc. shows a high concordance with quantitative pyrosequencing if identical CpG sites are analyzed in cell lines (Spearman r?=?0.88, p???0.0001), which is somewhat reduced in primary tumor specimens (Spearman r?=?0.86, p???0.0001). 80.7% of the CpG sites show an absolute difference in methylation level of less than 15 percentage points. If different CpG sites in the same CpG islands are targeted the concordance is lower (r?=?0.83 in cell lines and r?=?0.7 in primary tumors). The number of CpG sites representing microRNA genes and imprinted loci is very heterogeneous (range: 1 - 70 CpG sites for microRNAs and 1 - 288 for imprinted loci). CONCLUSIONS: The newly released 450k methylation array from Illumina, Inc. provides a genome-wide quantitative representation of DNA methylation aberrations in a convenient format. Overall, the congruence with pyrosequencing data is very good. However, for individual loci one should be careful to translate the ?-values directly into percent methylation levels.
| S-EPMC3420245 | biostudies-literature
Project description:Epigenome-wide association studies (EWAS) have focused primarily on DNA methylation as a chemically stable and functional epigenetic modification. However, the stability and accuracy of the measurement of methylation in different tissues and extraction types is still being actively studied, and the longitudinal stability of DNA methylation in commonly studied peripheral tissues is of great interest. Here, we used data from two studies, three tissue types, and multiple time points to assess the stability of DNA methylation measured with the Illumina Infinium HumanMethylation450 BeadChip array. Redundancy analysis enabled visual assessment of agreement of replicate samples overall and showed good agreement after removing effects of tissue type, age, and sex. At the probe level, analysis of variance contrasts separating technical and biological replicates clearly showed better agreement between technical replicates versus longitudinal samples, and suggested increased stability for buccal cells versus blood or blood spots. Intraclass correlations (ICCs) demonstrated that inter-individual variability is of similar magnitude to within-sample variability at many probes; however, as inter-individual variability increased, so did ICC. Furthermore, we were able to demonstrate decreasing agreement in methylation levels with time, despite a maximal sampling interval of only 576 days. Finally, at 6 popular candidate genes, there was a large range of stability across probes. Our findings highlight important sources of technical and biological variation in DNA methylation across different tissues over time. These data will help to inform longitudinal sampling strategies of future EWAS.
| S-EPMC5837078 | biostudies-literature