Report on the Infinium 450k Methylation Array Analysis Workshop: April 20, 2012 UCL,London, UK

A new platform for DNA methylome analysis is Illumina''s Infinium HumanMethylation450. This technology is an extension of the previous HumanMethylation27 BeadChip and allows the methylation status of 12 samples per chip and 4 to 8 chips (total of 48 to 96 samples) to be assessed simultaneously for more than 480,000 cytosines across the genome. The platform incorporates two different probe types using different assay designs (InfiniumI and InfiniumII). Although this has allowed the assessment of more CpG sites, it has also introduced technical variation between the two probe types, which has complicated the analysis process. Many groups are working on normalization methods and analysis pipelines while many others are struggling to make sense of their new data sets. This motivated the organization of a meeting held at University College London that focused solely on the analysis methods and problems related to this new platform. The meeting was attended by 125 computational and bench scientists from 11 countries. There were 10 speakers, a small poster session and a discussion session.     

Report on the 2nd Annual Infinium Humanmethylation450 Array Workshop

The Illumina Infinium HumanMethylation450 BeadChip – the successor to their hugely popular HumanMethylation27 BeadChip – is arguably the most prevalent platform for large-scale studies of DNA methylome analysis. After the success of last year’s meeting¹ that discussed initial analysis strategies for this then-new platform, this year’s meeting (held at Queen Mary, University of London) included the presentation of now established pipelines and normalization methods for data analysis, as well as some exciting tools for down-stream analysis. The importance of defining cell composition was a new topic mentioned by most speakers. The epigenome varies between cell types and insuring that methylation differences are related to sample treatment and not a differing cell population is essential. The meeting was attended by 215 computational and bench scientists from 18 countries. There were 11 speakers, a small poster session, and a discussion session. Talks were recorded and are now freely available at http://www.illumina.com/applications/epigenetics/array-based_methylation_analysis/methylation-array-analysis-education.ilmn     

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.     

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.     

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.     

Agreement in DNA methylation levels from the Illumina 450K array across batches,tissues, and time

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.     

Integration of silicon-via electrodes with different recording characteristics on a glass microprobe using a glass reflowing process

Electrodes on planar type microelectromechanical system (MEMS) microprobes mainly record neurons on the top-side of probe shaft (called a top-side electrode). However, it is often necessary to record neurons other than those on the top-side of the probe shaft. This study uses the glass reflowing technique to embed silicon-vias in a glass probe to implement a microprobe capable of recording neurons around the shaft. The proposed technology makes it possible to fabricate, distribute, and integrate four types of electrodes on the shaft: top-side, back-side, double-side, and sidewall electrodes. These electrodes have different recording characteristics. The in vitro and in vivo (using crayfish and rat brain) experiments in this study shows that the top-side and back-side electrodes are respectively more sensitive to neurons on the top-side and back-side of the probe shaft. In contrast, signals recorded by double-side electrode and sidewall electrode are equally sensitive to neurons around the probe shaft. This study enables the implementation and integration of these four types of electrodes, meeting the requirements of various neural applications.     

Histochemical classification based upon reaction types and its application to carbohydrate histochemistry

Summary A method of classification is presented, which divides histochemical visualization reactions into categories based on general reaction types. This scheme is dependent upon the reaction between two elements, the substrate and the probe. The substrate represents a tissue component(s) with one or more reactive groups that can combine directly with the probe. Alternatively, the substrate reactive groups are chemically modified (activation) with a suitable reagent before reaction with the probe. Probes are of three types: those that yield a coloured product, those that result in a colourless product, and those that produce a coloured product only after a further reaction.Methods used in carbohydrate histochemistry are divided into one, two and three probe reactions. Two probe reactions are further subdivided into sequences involving one or two coloured products (one and two dye sequences); three probe reactions into sequences involving one, two or three coloured products (one, two and three dye sequences). This classification permits the rationalization and organization of methods, and provides a framework for the examination of existing methods and the development of new ones.     

A novel in situ probe for oxygen uptake rate measurement in mammalian cell cultures

The newly developed in situ oxygen uptake rate (in situ OUR) probe presented in this article is based on the in situ microscope technology platform. It is designed to measure the oxygen uptake rate (OUR) of mammalian cells, an important parameter for metabolic flux analysis, inside a reactor (in situ) and in real-time. The system isolates a known volume of cell culture from the bulk inside the bioreactor, monitors the oxygen consumption over time, and releases the sample again. The sample is mixed during the measurement with a new agitation system to keep the cells in suspension and prevent oxygen concentration gradients. The OUR measurement system also doubles as a standard dissolved oxygen (DO) probe for process monitoring when it is not performing OUR measurements. It can be equipped with two different types of optical sensors (i.e., DO, pH) simultaneously or a conventional polarographic DO-probe (Clark type). This new probe was successfully tested in baby hamster kidney perfusion cell cultures.     

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.     

Sequence-matched probes produce increased cross-platform consistency and more reproducible biological results in microarray-based gene expression measurements

Cancer derived microarray data sets are routinely produced by various platforms that are either commercially available or manufactured by academic groups. The fundamental difference in their probe selection strategies holds the promise that identical observations produced by more than one platform prove to be more robust when validated by biology. However, cross-platform comparison requires matching corresponding probe sets. We are introducing here sequence-based matching of probes instead of gene identifier-based matching. We analyzed breast cancer cell line derived RNA aliquots using Agilent cDNA and Affymetrix oligonucleotide microarray platforms to assess the advantage of this method. We show, that at different levels of the analysis, including gene expression ratios and difference calls, cross-platform consistency is significantly improved by sequence- based matching. We also present evidence that sequence-based probe matching produces more consistent results when comparing similar biological data sets obtained by different microarray platforms. This strategy allowed a more efficient transfer of classification of breast cancer samples between data sets produced by cDNA microarray and Affymetrix gene-chip platforms.     

Genome-wide DNA methylation profiling of recurrent and non-recurrent chordomas

Chordomas are an aggressive rare type of malignant bone tumors arising from the remnant of the notochord. Chordomas occur mainly in vertebral bones and account for 1–4% of malignant bone tumors. Management and treatment of chordomas are difficult as they are resistant to conventional chemotherapy; therefore, they are mainly treated with surgery and radiation therapy. In this study, we performed DNA methylation profiling of 26 chordomas and normal nucleus pulposus samples plus UCH-1 chordoma cell line using the Illumina Infinium HumanMethylation450 BeadChips. Combined bisulfite restriction analysis and bisulfite sequencing was used to confirm the methylation data. Gene expression was analyzed using RT-PCR before and after 5-aza-2’-deoxycytidine (5-azaDC) treatment of chordoma cell lines. Analysis of the HumanMethylation450 BeadChip data led to the identification of 8,819 loci (2.9%) that were significantly differentially methylated (>0.2 average β-value difference) between chordomas and nucleus pulposus samples (adjusted P < 0.05). Among these, 5,868 probes (66.5%) were hypomethylated, compared to 2,951 (33.5%) loci that were hypermethylated in chordomas compared to controls. From the 2,951 differentially hypermethylated probes, 33.3% were localized in the promoter region (982 probes) and, among these, 104 probes showed cancer-specific hypermethylation. Ingenuity Pathway Analysis indicates that the cancer-specific differentially methylated loci are involved in various networks including cancer disease, nervous system development and function, cell death and survival, cellular growth, cellular development, and proliferation. Furthermore, we identified a subset of probes that were differentially methylated between recurrent and non-recurrent chordomas. BeadChip methylation data was confirmed for these genes and gene expression was shown to be upregulated in methylated chordoma cell lines after treatment with 5-azaDC. Understanding epigenetic changes in chordomas may provide insights into chordoma tumorigenesis and development of epigenetic biomarkers.     

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.     

Epigenetic landscape correlates with genetic subtype but does not predict outcome in childhood acute lymphoblastic leukemia

Although children with acute lymphoblastic leukemia (ALL) generally have a good outcome, some patients do relapse and survival following relapse is poor. Altered DNA methylation is highly prevalent in ALL and raises the possibility that DNA methylation-based biomarkers could predict patient outcome. In this study, genome-wide methylation analysis, using the Illumina Infinium HumanMethylation450 BeadChip platform, was carried out on 52 diagnostic patient samples from 4 genetic subtypes [ETV6-RUNX1, high hyperdiploidy (HeH), TCF3-PBX1 and dic(9;20)(p11–13;q11)] in a 1:1 case-control design with patients who went on to relapse (as cases) and patients achieving long-term remission (as controls). Pyrosequencing assays for selected loci were used to confirm the array-generated data. Non-negative matrix factorization consensus clustering readily clustered samples according to genetic subgroups and gene enrichment pathway analysis suggested that this is in part driven by epigenetic disruption of subtype specific signaling pathways. Multiple bioinformatics approaches (including bump hunting and individual locus analysis) were used to identify CpG sites or regions associated with outcome. However, no associations with relapse were identified. Our data revealed that ETV6-RUNX1 and dic(9;20) subtypes were mostly associated with hypermethylation; conversely, TCF3-PBX1 and HeH were associated with hypomethylation. We observed significant enrichment of the neuroactive ligand-receptor interaction pathway in TCF3-PBX1 as well as an enrichment of genes involved in immunity and infection pathways in ETV6-RUNX1 subtype. Taken together, our results suggest that altered DNA methylation may have differential impacts in distinct ALL genetic subtypes.