Causes of blood methylomic variation for middle-aged women measured by the HumanMethylation450 array

To address the limitations in current classic twin/family research on the genetic and/or environmental causes of human methylomic variation, we measured blood DNA methylation for 479 women (mean age 56 years) including 66 monozygotic (MZ), 66 dizygotic (DZ) twin pairs and 215 sisters of twins, and 11 random technical duplicates using the HumanMethylation450 array. For each methylation site, we estimated the correlation for pairs of duplicates, MZ twins, DZ twins, and siblings, fitted variance component models by assuming the variation is explained by genetic factors, by shared and individual environmental factors, and by independent measurement error, and assessed the best fitting model. We found that the average (standard deviation) correlations for duplicate, MZ, DZ, and sibling pairs were 0.10 (0.35), 0.07 (0.21), -0.01 (0.14) and -0.04 (0.07). At the genome-wide significance level of 10^?7, 93.3% of sites had no familial correlation, and 5.6%, 0.1%, and 0.2% of sites were correlated for MZ, DZ, and sibling pairs. For 86.4%, 6.9%, and 7.1% of sites, the best fitting model included measurement error only, a genetic component, and at least one environmental component. For the 13.6% of sites influenced by genetic and/or environmental factors, the average proportion of variance explained by environmental factors was greater than that explained by genetic factors (0.41 vs. 0.37, P value <10^?15). Our results are consistent with, for middle-aged woman, blood methylomic variation measured by the HumanMethylation450 array being largely explained by measurement error, and more influenced by environmental factors than by genetic factors.     

Measuring the methylome in clinical samples: Improved processing of the Infinium Human Methylation450 BeadChip Array

The Infinium Human Methylation450 BeadChip Array^TM (Infinium 450K) is an important tool for studying epigenetic patterns associated with disease. This array offers a high-throughput, low cost alternative to more comprehensive sequencing-based methodologies. Here we compare data generated by interrogation of the same seven clinical samples by Infinium 450K and reduced representation bisulfite sequencing (RRBS). This is the largest data set comparing Infinium 450K array to the comprehensive RRBS methodology reported so far. We show good agreement between the two methodologies. A read depth of four or more reads in the RRBS data was sufficient to achieve good agreement with Infinium 450K. However, we observe that intermediate methylation values (20–80%) are more variable between technologies than values at the extremes of the bimodal methylation distribution. We describe careful processing of Infinium 450K data to correct for known limitations and batch effects. Using methodologies proposed by others and newly implemented and combined in this report, agreement of Infinium 450K data with independent techniques can be vastly improved.     

Measuring the methylome in clinical samples: Improved processing of the Infinium Human Methylation450 BeadChip Array

The Infinium Human Methylation450 BeadChip Array^TM (Infinium 450K) is an important tool for studying epigenetic patterns associated with disease. This array offers a high-throughput, low cost alternative to more comprehensive sequencing-based methodologies. Here we compare data generated by interrogation of the same seven clinical samples by Infinium 450K and reduced representation bisulfite sequencing (RRBS). This is the largest data set comparing Infinium 450K array to the comprehensive RRBS methodology reported so far. We show good agreement between the two methodologies. A read depth of four or more reads in the RRBS data was sufficient to achieve good agreement with Infinium 450K. However, we observe that intermediate methylation values (20–80%) are more variable between technologies than values at the extremes of the bimodal methylation distribution. We describe careful processing of Infinium 450K data to correct for known limitations and batch effects. Using methodologies proposed by others and newly implemented and combined in this report, agreement of Infinium 450K data with independent techniques can be vastly improved.     

Validation of DNA methylation profiling in formalin-fixed paraffin-embedded samples using the Infinium HumanMethylation450 Microarray

A formalin-fixed paraffin-embedded (FFPE) sample usually yields highly degraded DNA, which limits the use of techniques requiring high-quality DNA, such as Infinium Methylation microarrays. To overcome this restriction, we have applied an FFPE restoration procedure consisting of DNA repair and ligation processes in a set of paired fresh-frozen (FF) and FFPE samples. We validated the FFPE results in comparison with matched FF samples, enabling us to use FFPE samples on the Infinium HumanMethylation450 Methylation array.     

Validation of DNA methylation profiling in formalin-fixed paraffin-embedded samples using the Infinium HumanMethylation450 Microarray

A formalin-fixed paraffin-embedded (FFPE) sample usually yields highly degraded DNA, which limits the use of techniques requiring high-quality DNA, such as Infinium Methylation microarrays. To overcome this restriction, we have applied an FFPE restoration procedure consisting of DNA repair and ligation processes in a set of paired fresh-frozen (FF) and FFPE samples. We validated the FFPE results in comparison with matched FF samples, enabling us to use FFPE samples on the Infinium HumanMethylation450 Methylation array.     

Reducing the risk of false discovery enabling identification of biologically significant genome-wide methylation status using the HumanMethylation450 array

Background

The Illumina HumanMethylation450 BeadChip (HM450K) measures the DNA methylation of 485,512 CpGs in the human genome. The technology relies on hybridization of genomic fragments to probes on the chip. However, certain genomic factors may compromise the ability to measure methylation using the array such as single nucleotide polymorphisms (SNPs), small insertions and deletions (INDELs), repetitive DNA, and regions with reduced genomic complexity. Currently, there is no clear method or pipeline for determining which of the probes on the HM450K bead array should be retained for subsequent analysis in light of these issues.

Results

We comprehensively assessed the effects of SNPs, INDELs, repeats and bisulfite induced reduced genomic complexity by comparing HM450K bead array results with whole genome bisulfite sequencing. We determined which CpG probes provided accurate or noisy signals. From this, we derived a set of high-quality probes that provide unadulterated measurements of DNA methylation.

Conclusions

Our method significantly reduces the risk of false discoveries when using the HM450K bead array, while maximising the power of the array to detect methylation status genome-wide. Additionally, we demonstrate the utility of our method through extraction of biologically relevant epigenetic changes in prostate cancer.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-51) contains supplementary material, which is available to authorized users.     

A systematic study of normalization methods for Infinium 450K methylation data using whole-genome bisulfite sequencing data

DNA methylation plays an important role in disease etiology. The Illumina Infinium HumanMethylation450 (450K) BeadChip is a widely used platform in large-scale epidemiologic studies. This platform can efficiently and simultaneously measure methylation levels at ∼480,000 CpG sites in the human genome in multiple study samples. Due to the intrinsic chip design of 2 types of chemistry probes, data normalization or preprocessing is a critical step to consider before data analysis. To date, numerous methods and pipelines have been developed for this purpose, and some studies have been conducted to evaluate different methods. However, validation studies have often been limited to a small number of CpG sites to reduce the variability in technical replicates. In this study, we measured methylation on a set of samples using both whole-genome bisulfite sequencing (WGBS) and 450K chips. We used WGBS data as a gold standard of true methylation states in cells to compare the performances of 8 normalization methods for 450K data on a genome-wide scale. Analyses on our dataset indicate that the most effective methods are peak-based correction (PBC) and quantile normalization plus β-mixture quantile normalization (QN.BMIQ). To our knowledge, this is the first study to systematically compare existing normalization methods for Illumina 450K data using novel WGBS data. Our results provide a benchmark reference for the analysis of DNA methylation chip data, particularly in white blood cells.     

Genome-scale analysis of DNA methylation in colorectal cancer using Infinium HumanMethylation450 BeadChips

Illumina’s Infinium HumanMethylation450 BeadChip arrays were used to examine genome-wide DNA methylation profiles in 22 sample pairs from colorectal cancer (CRC) and adjacent tissues and 19 colon tissue samples from cancer-free donors. We show that the methylation profiles of tumors and healthy tissue samples can be clearly distinguished from one another and that the main source of methylation variability is associated with disease status. We used different statistical approaches to evaluate the methylation data. In general, at the CpG-site level, we found that common CRC-specific methylation patterns consist of at least 15,667 CpG sites that were significantly different from either adjacent healthy tissue or tissue from cancer-free subjects. Of these sites, 10,342 were hypermethylated in CRC, and 5,325 were hypomethylated. Hypermethylated sites were common in the maximum number of sample pairs and were mostly located in CpG islands, where they were significantly enriched for differentially methylated regions known to be cancer-specific. In contrast, hypomethylated sites were mostly located in CpG shores and were generally sample-specific. Despite the considerable variability in methylation data, we selected a panel of 14 highly robust candidates showing methylation marks in genes SND1, ADHFE1, OPLAH, TLX2, C1orf70, ZFP64, NR5A2, and COL4A. This set was successfully cross-validated using methylation data from 209 CRC samples and 38 healthy tissue samples from The Cancer Genome Atlas consortium (AUC = 0.981 [95% CI: 0.9677–0.9939], sensitivity = 100% and specificity = 82%). In summary, this study reports a large number of loci with novel differential methylation statuses, some of which may serve as candidate markers for diagnostic purposes.     

A systematic assessment of normalization approaches for the Infinium 450K methylation platform

The Illumina Infinium HumanMethylation450 BeadChip has emerged as one of the most popular platforms for genome wide profiling of DNA methylation. While the technology is wide-spread, systematic technical biases are believed to be present in the data. For example, this array incorporates two different chemical assays, i.e., Type I and Type II probes, which exhibit different technical characteristics and potentially complicate the computational and statistical analysis. Several normalization methods have been introduced recently to adjust for possible biases. However, there is considerable debate within the field on which normalization procedure should be used and indeed whether normalization is even necessary. Yet despite the importance of the question, there has been little comprehensive comparison of normalization methods. We sought to systematically compare several popular normalization approaches using the Norwegian Mother and Child Cohort Study (MoBa) methylation data set and the technical replicates analyzed with it as a case study. We assessed both the reproducibility between technical replicates following normalization and the effect of normalization on association analysis. Results indicate that the raw data are already highly reproducible, some normalization approaches can slightly improve reproducibility, but other normalization approaches may introduce more variability into the data. Results also suggest that differences in association analysis after applying different normalizations are not large when the signal is strong, but when the signal is more modest, different normalizations can yield very different numbers of findings that meet a weaker statistical significance threshold. Overall, our work provides useful, objective assessment of the effectiveness of key normalization methods.     

A systematic assessment of normalization approaches for the Infinium 450K methylation platform

The Illumina Infinium HumanMethylation450 BeadChip has emerged as one of the most popular platforms for genome wide profiling of DNA methylation. While the technology is wide-spread, systematic technical biases are believed to be present in the data. For example, this array incorporates two different chemical assays, i.e., Type I and Type II probes, which exhibit different technical characteristics and potentially complicate the computational and statistical analysis. Several normalization methods have been introduced recently to adjust for possible biases. However, there is considerable debate within the field on which normalization procedure should be used and indeed whether normalization is even necessary. Yet despite the importance of the question, there has been little comprehensive comparison of normalization methods. We sought to systematically compare several popular normalization approaches using the Norwegian Mother and Child Cohort Study (MoBa) methylation data set and the technical replicates analyzed with it as a case study. We assessed both the reproducibility between technical replicates following normalization and the effect of normalization on association analysis. Results indicate that the raw data are already highly reproducible, some normalization approaches can slightly improve reproducibility, but other normalization approaches may introduce more variability into the data. Results also suggest that differences in association analysis after applying different normalizations are not large when the signal is strong, but when the signal is more modest, different normalizations can yield very different numbers of findings that meet a weaker statistical significance threshold. Overall, our work provides useful, objective assessment of the effectiveness of key normalization methods.     

Novel region discovery method for Infinium 450K DNA methylation data reveals changes associated with aging in muscle and neuronal pathways

We describe a methodology for detecting differentially methylated regions (DMRs) and variably methylated regions (VMRs), in data from Infinium 450K arrays that are very widely used in epigenetic studies. Region detection is more specific than single CpG analysis as it increases the extent of common findings between studies, and is more powerful as it reduces the multiple testing problem inherent in epigenetic whole‐genome association studies (EWAS). In addition, results driven by single erroneous probes are removed. We have used multiple publicly available Infinium 450K data sets to generate a consensus list of DMRs for age, supporting the hypothesis that aging is associated with specific epigenetic modifications. The consensus aging DMRs are significantly enriched for muscle biogenesis pathways. We find a massive increase in VMRs with age and in regions of the genome associated with open chromatin and neurotransmission. Old age VMRs are significantly enriched for neurotransmission pathways. EWAS studies should investigate the role of this interindividual variation in DNA methylation, in the age‐associated diseases of sarcopenia and dementia.     

Novel imprinted single CpG sites found by global DNA methylation analysis in human parthenogenetic induced pluripotent stem cells

Genomic imprinting is the process of epigenetic modification whereby genes are expressed in a parent-of-origin dependent manner; it plays an important role in normal growth and development. Parthenogenetic embryos contain only the maternal genome. Parthenogenetic embryonic stem cells could be useful for studying imprinted genes. In humans, mature cystic ovarian teratomas originate from parthenogenetic activation of oocytes; they are composed of highly differentiated mature tissues containing all three germ layers. To establish human parthenogenetic induced pluripotent stem cell lines (PgHiPSCs), we generated parthenogenetic fibroblasts from ovarian teratoma tissues. We compared global DNA methylation status of PgHiPSCs with that of biparental human induced pluripotent stem cells by using Illumina Infinium HumanMethylation450 BeadChip array. This analysis identified novel single imprinted CpG sites. We further tested DNA methylation patterns of two of these sites using bisulfite sequencing and described novel candidate imprinted CpG sites. These results confirm that PgHiPSCs are a powerful tool for identifying imprinted genes and investigating their roles in human development and diseases.     

Differential DNA methylation and PM2.5 species in a 450K epigenome-wide association study

Although there is growing evidence that exposure to ambient particulate matter is associated with global DNA methylation and gene-specific methylation, little is known regarding epigenome-wide changes in DNA methylation in relation to particles and, especially, particle components. Using the Illumina Infinium HumanMethylation450 BeadChip, we examined the relationship between one-year moving averages of PM_2.5 species (Al, Ca, Cu, Fe, K, Na, Ni, S, Si, V, and Zn) and DNA methylation at 484,613 CpG probes in a longitudinal cohort that included 646 subjects. Bonferroni correction was applied to adjust for multiple comparisons. Bioinformatics analysis of the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment was also performed. We observed 20 Bonferroni significant (P-value < 9.4× 10^?9) CpGs for Fe, 8 for Ni, and 1 for V. Particularly, methylation at Schlafen Family Member 11 (SLFN11) cg10911913 was positively associated with measured levels of all 3 species. The SLFN11 gene codes for an interferon-induced protein that inhibits retroviruses and sensitizes cancer cells to DNA-damaging agents. Bioinformatics analysis suggests that gene targets may be relevant to pathways including cancers, signal transduction, and cell growth and death. Ours is the first study to examine the epigenome-wide association between ambient particles species and DNA methylation. We found that long-term exposures to specific components of ambient particle pollution, especially particles emitted during oil combustion, were associated with methylation changes in genes relevant to immune responses. Our findings provide insight into potential biologic mechanisms on an epigenetic level.     

High-density array analysis of DNA methylation in Tamoxifen-resistant breast cancer cell lines

Roughly two-thirds of all breast cancers are ERα-positive and can be treated with the antiestrogen, Tamoxifen, however resistance occurs in 33% of women who take the drug for more than 5 y. Aberrant DNA methylation, an epigenetic mechanism that alters gene expression in cancer, is thought to play a role in this resistance. To develop an understanding of Tamoxifen-resistance and identify novel pathways and targets of aberrant methylation, DNA from MCF-7 breast cancer cells and Tamoxifen-resistant derivatives, TMX2–11 and TMX2–28, were analyzed using the Illumina HumanMethylation450 BeadChip. Normalizing against MCF-7 values, ERα-positive TMX2–11 had 4000 hypermethylated sites and ERα-negative TMX2–28 had over 33?000. Analysis of CpG sites altered in both TMX2–11 and TMX2–28 revealed that the Tamoxifen-resistant cell lines share 3000 hypermethylated and 200 hypomethylated CpGs. ZNF350 and MAGED1, two genes hypermethylated in both cell lines, were examined in greater detail. Treatment with 5-aza-2′deoxycitidine caused a significant reduction in promoter methylation of both ZNF350 and MAGED1 and a corresponding increase in expression in TMX2–28. A similar relationship between methylation and expression was not detected in TMX2–11. Our findings are indicative of the variable responses to methylation-targeted breast cancer therapy and highlight the need for biomarkers that accurately predict treatment outcome.     

High-density array analysis of DNA methylation in Tamoxifen-resistant breast cancer cell lines

Roughly two-thirds of all breast cancers are ERα-positive and can be treated with the antiestrogen, Tamoxifen, however resistance occurs in 33% of women who take the drug for more than 5 y. Aberrant DNA methylation, an epigenetic mechanism that alters gene expression in cancer, is thought to play a role in this resistance. To develop an understanding of Tamoxifen-resistance and identify novel pathways and targets of aberrant methylation, DNA from MCF-7 breast cancer cells and Tamoxifen-resistant derivatives, TMX2–11 and TMX2–28, were analyzed using the Illumina HumanMethylation450 BeadChip. Normalizing against MCF-7 values, ERα-positive TMX2–11 had 4000 hypermethylated sites and ERα-negative TMX2–28 had over 33 000. Analysis of CpG sites altered in both TMX2–11 and TMX2–28 revealed that the Tamoxifen-resistant cell lines share 3000 hypermethylated and 200 hypomethylated CpGs. ZNF350 and MAGED1, two genes hypermethylated in both cell lines, were examined in greater detail. Treatment with 5-aza-2′deoxycitidine caused a significant reduction in promoter methylation of both ZNF350 and MAGED1 and a corresponding increase in expression in TMX2–28. A similar relationship between methylation and expression was not detected in TMX2–11. Our findings are indicative of the variable responses to methylation-targeted breast cancer therapy and highlight the need for biomarkers that accurately predict treatment outcome.     

DNA methylation in embryonic stem cells

Embryonic stem cells (ESCs) are pluripotent, self‐renewing cells. These cells can be used in applications such as cell therapy, drug development, disease modeling, and the study of cellular differentiation. Investigating the interplay of epigenetics, genetics, and gene expression in control of pluripotence and differentiation could give important insights on how these cells function. One of the best known epigenetic factors is DNA methylation, which is a major mechanism for regulation of gene expression. This phenomenon is mostly seen in imprinted genes and X‐chromosome inactivation where DNA methylation of promoter regions leads to repression of gene expression. Differential DNA methylation of pluripotence‐associated genes such as Nanog and Oct4/Pou5f1 has been observed between pluripotent and differentiated cells. It is clear that tight regulation of DNA methylation is necessary for normal development. As more associations between aberrant DNA methylation and disease are reported, the demand for high‐throughput approaches for DNA methylation analysis has increased. In this article, we highlight these methods and discuss recent DNA methylation studies on ESCs. J. Cell. Biochem. 109: 1–6, 2010. © 2009 Wiley‐Liss, Inc.     

Methyl CpG Binding Domain Ultra‐Sequencing: a novel method for identifying inter‐individual and cell‐type‐specific variation in DNA methylation

Experience‐dependent changes in DNA methylation can exert profound effects on neuronal function and behaviour. A single learning event can induce a variety of DNA modifications within the neuronal genome, some of which may be common to all individuals experiencing the event, whereas others may occur in a subset of individuals. Variations in experience‐induced DNA methylation may subsequently confer increased vulnerability or resilience to the development of neuropsychiatric disorders. However, the detection of experience‐dependent changes in DNA methylation in the brain has been hindered by the interrogation of heterogeneous cell populations, regional differences in epigenetic states and the use of pooled tissue obtained from multiple individuals. Methyl CpG Binding Domain Ultra‐Sequencing (MBD Ultra‐Seq) overcomes current limitations on genome‐wide epigenetic profiling by incorporating fluorescence‐activated cell sorting and sample‐specific barcoding to examine cell‐type‐specific CpG methylation in discrete brain regions of individuals. We demonstrate the value of this method by characterizing differences in 5‐methylcytosine (5mC) in neurons and non‐neurons of the ventromedial prefrontal cortex of individual adult C57BL/6 mice, using as little as 50 ng of genomic DNA per sample. We find that the neuronal methylome is characterized by greater CpG methylation as well as the enrichment of 5mC within intergenic loci. In conclusion, MBD Ultra‐Seq is a robust method for detecting DNA methylation in neurons derived from discrete brain regions of individual animals. This protocol will facilitate the detection of experience‐dependent changes in DNA methylation in a variety of behavioural paradigms and help identify aberrant experience‐induced DNA methylation that may underlie risk and resiliency to neuropsychiatric disease.     

Improved reproducibility in genome-wide DNA methylation analysis for PAXgene-fixed samples compared with restored formalin-fixed and paraffin-embedded DNA

Formalin fixation has been the standard method for conservation of clinical specimens for decades. However, a major drawback is the high degradation of nucleic acids, which complicates its use in genome-wide analyses. Unbiased identification of biomarkers, however, requires genome-wide studies, precluding the use of the valuable archives of specimens with long-term follow-up data. Therefore, restoration protocols for DNA from formalin-fixed and paraffin-embedded (FFPE) samples have been developed, although they are cost-intensive and time-consuming. An alternative to FFPE and snap-freezing is the PAXgene Tissue System, developed for simultaneous preservation of morphology, proteins, and nucleic acids. In the current study, we compared the performance of DNA from either PAXgene or formalin-fixed tissues to snap-frozen material for genome-wide DNA methylation analysis using the Illumina 450K BeadChip. Quantitative DNA methylation analysis demonstrated that the methylation profile in PAXgene-fixed tissues showed, in comparison with restored FFPE samples, a higher concordance with the profile detected in frozen samples. We demonstrate, for the first time, that DNA from PAXgene conserved tissue performs better compared with restored FFPE DNA in genome-wide DNA methylation analysis. In addition, DNA from PAXgene tissue can be directly used on the array without prior restoration, rendering the analytical process significantly more time- and cost-effective.     

Comparison of DNA methylation measured by Illumina 450K and EPIC BeadChips in blood of newborns and 14-year-old children

Analysis of DNA methylation helps to understand the effects of environmental exposures as well as the role of epigenetics in human health. Illumina, Inc. recently replaced the HumanMethylation450 BeadChip (450K) with the EPIC BeadChip, which nearly doubles the measured CpG sites to >850,000. Although the new chip uses the same underlying technology, it is important to establish if data between the two platforms are comparable within cohorts and for meta-analyses. DNA methylation was assessed by 450K and EPIC using whole blood from newborn (n = 109) and 14-year-old (n = 86) participants of the Center for the Health Assessment of Mothers and Children of Salinas. The overall per-sample correlations were very high (r >0.99), although many individual CpG sites, especially those with low variance of methylation, had lower correlations (median r = 0.24). There was also a small subset of CpGs with large mean methylation β-value differences between platforms, in both the newborn and 14-year datasets. However, estimates of cell type proportion prediction by 450K and EPIC were highly correlated at both ages. Finally, differentially methylated positions between boys and girls replicated very well by both platforms in newborns and older children. These findings are encouraging for application of combined data from EPIC and 450K platforms for birth cohorts and other population studies. These data in children corroborate recent comparisons of the two BeadChips in adults and in cancer cell lines. However, researchers should be cautious when characterizing individual CpG sites and consider independent methods for validation of significant hits.     

An evaluation of analysis pipelines for DNA methylation profiling using the Illumina HumanMethylation450 BeadChip platform

The proper identification of differentially methylated CpGs is central in most epigenetic studies. The Illumina HumanMethylation450 BeadChip is widely used to quantify DNA methylation; nevertheless, the design of an appropriate analysis pipeline faces severe challenges due to the convolution of biological and technical variability and the presence of a signal bias between Infinium I and II probe design types. Despite recent attempts to investigate how to analyze DNA methylation data with such an array design, it has not been possible to perform a comprehensive comparison between different bioinformatics pipelines due to the lack of appropriate data sets having both large sample size and sufficient number of technical replicates. Here we perform such a comparative analysis, targeting the problems of reducing the technical variability, eliminating the probe design bias and reducing the batch effect by exploiting two unpublished data sets, which included technical replicates and were profiled for DNA methylation either on peripheral blood, monocytes or muscle biopsies. We evaluated the performance of different analysis pipelines and demonstrated that: (1) it is critical to correct for the probe design type, since the amplitude of the measured methylation change depends on the underlying chemistry; (2) the effect of different normalization schemes is mixed, and the most effective method in our hands were quantile normalization and Beta Mixture Quantile dilation (BMIQ); (3) it is beneficial to correct for batch effects. In conclusion, our comparative analysis using a comprehensive data set suggests an efficient pipeline for proper identification of differentially methylated CpGs using the Illumina 450K arrays.