Methylation

Data Generation Protocols Data Analysis Protocols
HpaII tiny fragment Enrichment by Ligation-mediated PCR (HELP) Assay with NimbleGen Arrays (Roche) ALL P2ALL symbol
Infinium HumanMethylation27 BeadChip (Illumina) AMLAML symbol
Infinium HumanMethylation450 BeadChip (Illumina) AMLAML symbol , CCSKKidney tumor symbol , NBLNeuroblastoma symbol , OSOsteosarcoma symbol , WTKidney tumor symbol
Infinium MethylationEPIC BeadChip Kit (Illumina) ALALALL symbol

HpaII tiny fragment Enrichment by Ligation-mediated PCR (HELP) Assay with NimbleGen Arrays (Roche) for Acute Lymphoblastic Leukemia Phase 2 (ALL P2)

*Protocol performed at Weill Cornell Medical College.

DNA was extracted using QIAGEN QIAamp DNA Mini Kit according to manufacturer’s protocol at St. Jude’s Children’s Research Hospital.         

Nucleic acid labeling, hybridization and array scanning protocols were used according to NimbleGen (Roche) manufacturer’s protocols (see NimbleGen Dual-color DNA Labeling Kit, NimbleGen Hybridization System 4 and NimbleGen MS 200 Microarray Scanner respectively).                 

Normalization data transformation protocols were carried out at as follows: Median normalized log2 ratio of signal intensity of HpaII and MspI, as detailed in Thompson et al, Bioinformatics 2008;24:1161-7. Software: NimbleScan version 2.5.26, Value: median normalized log2-transformed HpaII/MspI ratios, "NA": if MspI signal intensity < 1 mean absolute deviation (MAD) above the median of random probe signals.                  

Transformation of Level 2 data matrix into per sample Level 3 files with probe set gene annotations added was performed at NCI Center for Bioinformatics and Information Technology.

Infinium HumanMethylation27 Bead Chip (Illumina) for Acute Myeloid Leukemia (AML)

*Protocols performed at the Johns Hopkins University.

Bisulfite conversion of genomic DNA was performed with EZ DNA methylation Kit (Zymo Research, Irvine, CA, #D5002) following the manufacturer’s protocol with modifications for the Infinium Methylation Assay. Briefly, one microgram of genomic DNA was mixed with 5 µl of Dilution Buffer and incubated at 37°C for 15 minutes and then mixed with 100 µl of conversion reagent prepared as instructed in the protocol. Mixtures were incubated in a thermocycler for 16 cycles at 95°C for 30 seconds and 50°C for 60 minutes. Bisulfite-converted DNA samples were loaded onto the provided 96-column plates for desulphonation, washing and elution. The concentration of bisulfite-converted, eluted DNA was measured by UV-absorbance using a NanoDrop-1000 (Thermo Fisher Scientific, Waltham, MA). Bisulfite-converted genomic DNA was analyzed using the Infinium Human Methylation27 Beadchip Kit (Illumina, San Diego, CA, #WG-311-1202). DNA amplification, fragmentation, array hybridization, extension and staining were performed with reagents provided in the kit according to the manufacturer’s protocol (Illumina Infinium II Methylation Assay, #WG-901-2701). Briefly, 4 µl of bisulfite-converted genomic DNA at a minimum concentration of 20 ng/µL) was added to 0.8 ml 96-well storage plate (Thermo Fisher Scientific), denatured in 0.014N sodium hydroxide, neutralized and then amplified for 20-24 hours at 37°C. Samples were fragmented at 37°C for 60 minutes and precipitated in isopropanol. Re-suspended samples were denatured in a 96-well plate heat block at 95°C for 20 minutes. 15 µl of each sample was loaded onto a 12-sample BeadChip, assembled in the hybridization chamber as instructed by the manufacturer and incubated at 48°C for 16-20 hours.  Following hybridization, the BeadChips were washed and assembled in a fluid flow-through station for primer-extension reaction and staining with reagents and buffers provided.

Polymer-coated BeadChips were scanned in an iScan scanner (Illumina) using Inf Methylation mode.

Signal intensity and Beta value data were extracted using the Methylation module of GenomeStudio (Illumina, v2011.1) software following the Methylation analysis pipeline without normalization or background subtraction using BeadChip content descriptors provided by the manufacturer (HumanMethylation27_270596_v.1.2.bpm).

Summary beta values for each locus with annotations for Illumina probe name, gene symbol, chromosome and CpG position (UCSC hg18).

Data were normalized using functional normalization (funnorm) as implemented in the minfi package and summarized as beta values [M /(M+U)] with annotation at each locus for Illumina probe name, gene symbol, chromosome and CpG position (UCSC hg19).  Probes having an annotated SNV within the CpG or SBE site are masked as NA across all samples.  Probes where the non-detection probability was >  0.01 are masked as NA for individual samples.

Infinium HumanMethylation450 Bead Chip (Illumina) for Acute Myeloid Leukemia (AML)

*Protocols performed at the Phoenix Children’s Hospital.

Bisulfite conversion of genomic DNA was performed with EZ DNA methylation Kit (Zymo Research, Irvine, CA, #D5002) following the manufacturer’s protocol with modifications for the Infinium Methylation Assay. Briefly, one microgram of genomic DNA was mixed with 5 µl of Dilution Buffer and incubated at 37°C for 15 minutes and then mixed with 100 µl of conversion reagent prepared as instructed in the protocol. Mixtures were incubated in a thermocycler for 16 cycles at 95°C for 30 seconds and 50°C for 60 minutes. Bisulfite-converted DNA samples were loaded onto the provided 96-column plates for desulphonation, washing and elution. Bisulfite-converted genomic DNA was analyzed using the Infinium Human Methylation450K Beadchip Kit (Illumina, San Diego, CA, #WG-314-1001). DNA amplification, fragmentation, array hybridization, extension and staining were performed with reagents provided in the kit according to the manufacturer’s protocol (Illumina Infinium II Methylation Assay, #WG-901-2701). Briefly, 4 µl of bisulfite-converted genomic DNA was added to 0.8 ml 96-well storage plate (Thermo Fisher Scientific), denatured in 0.014N sodium hydroxide, neutralized and then amplified for 20-24 hours at 37°C. Samples were fragmented at 37°C for 60 minutes and precipitated in isopropanol. Re-suspended samples were denatured in a 96-well plate heat block at 95°C for 20 minutes. 15 µl of each sample was loaded onto a 12-sample BeadChip, assembled in the hybridization chamber as instructed by the manufacturer and incubated at 48°C for 16-20 hours.  Following hybridization, the BeadChips were washed and assembled in a fluid flow-through station for primer-extension reaction and staining with reagents and buffers provided.

Polymer-coated BeadChips were scanned in an iScan scanner (Illumina) using Inf Methylation mode.

Methylated and unmethylated signal intensity and detection p-values were extracted after background correction and dye-bias equalization by normal-exponential convolution ('noob') as implemented in the minfi package.

Summary beta values for each locus with annotations for Illumina probe name, gene symbol, chromosome and CpG position (UCSC hg18).

Data were normalized using functional normalization (funnorm) as implemented in the minfi package and summarized as beta values [M /(M+U)] with annotation at each locus for Illumina probe name, gene symbol, chromosome and CpG position (UCSC hg19).  Probes having an annotated SNV within the CpG or SBE site are masked as NA across all samples.  Probes where the non-detection probability was > 0.01 are masked as NA for individual samples.

Infinium HumanMethylation450 Bead Chip (Illumina) for Clear Cell Sarcoma of the Kidney (CCSK)

*Protocol performed at Ann & Robert H. Lurie Children's Hospital.

DNA was extracted from tumor samples at Nationwide Children's BioPathology Center (BPC) by using the standard BPC protocol. The DNA samples were analyzed by Pico green to verify gDNA concentration, spectrophotometry to verify DNA purity, and gel electrophoresis to verify DNA quality.  DNA samples (1.5 ug) diluted in nuclease-free water were provided to the Northwestern University Genomics Core in 96-well plate format for Illumina 450K DNA methylation analysis. Randomly selected samples were tested by Northwestern University to verify that the correct concentration had been provided.                

Nucleic acid hybridization and labeling were performed according to the manufacturer's protocol for the Illumina 450K array at the Northwestern University Genomics Core Facility. Nucleic acid labeling is completed after the hybridization step with Illumina Infinium 450K arrays.   

The array scanning protocol was performed according to the manufacturer's protocol for the Illumina 450K array at the Northwestern University Genomics Core Facility.                 

Raw data files (1 red and 1 green .idat file per sample and 1 .sdf file from each array, which included 12 samples per array) were processed at the Northwestern University Genomics Core Facility by BeadStudio software. The following subtables were generated by BeadStudio and were downloaded in .txt format (Level 2 data): (1) the Sample Methylation Profile .txt, (2) the Control Profile .txt, and (3) the Control Probe Profile .txt. Several quality control steps were used for these data. Samples were subjected to the internal quality controls in the Bioconductor lumi package. Samples were subjected to a color balance check using the Bioconductor lumi package. Gender analysis was performed to increase our confidence that the data correctly corresponded to the expected sample. Because one of the X chromosomes is heavily methylated in females, the density of X-chromosome methylation is a good indicator of gender. Unsupervised hierarchical clustering based on all methylated regions on the X-chromosome was performed using average-linkage clustering with CLUSTER and the results were displayed with TREEVIEW. The tumors were assigned a gender based on their clustering in the tumor dendogram, which was checked against the known gender of the sample.               

The Sample Methylation Profile text, which included information for all of the samples, was broken down at the DCC into a single .txt file (level 2 data) per sample containing the following columns: (1) the sample ID, (2) probe name, (3) AVG_Beta value, (4) gene Symbol, (5) chromosome, and (6) position.

Infinium HumanMethylation450 Bead Chip (Illumina) for Neuroblastoma (NBL)

*Protocols performed at the USC Epigenome Center.

Labeling, hybridization and scanning protocols were performed following the manufacturer’s protocol using the Infinium Human Methylation450K Beadchip Kit (Illumina, San Diego, CA, #WG-314-1001).

Level 2 data contain background-corrected methylated (M) and unmethylated (U) summary intensities as extracted by the methylumi package.  Non-detection probabilities (P-values) were computed as the minimum of the two values (one per allele) for the empirical cumulative density function of the negative control probes in the appropriate color channel. Background correction is performed via normal-exponential deconvolution using out-of-band probes (Triche, Jr. et al, Nucl. Acids Res 2013). Multiple-batch archives have the intensities in each of the two channels multiplicatively scaled to match a reference sample (sample with R/G ratio of normalization control probes closest to 1.0.).

Level 3 data contain derived summary measures (beta values: M/(M+U) for each interrogated locus) with annotations (based on Illumina's manifest on GEO, GPL13534) for gene symbol, chromosome (UCSC hg19, Feb 2009), and CpG/CpH coordinate (UCSC hg19, Feb 2009). Probes having a common SNP (common SNP is a SNP with Minor Allele Frequency > 1% as defined by the UCSC snp135common track) within 10bp of the interrogated CpG site or having 15bp from the interrogated CpG site overlap with a REPEAT element (as defined by RepeatMasker and Tandem Repeat Finder Masks based on UCSC hg19, Feb 2009) are masked as NA across all samples, and probes with a non-detection probability (P-value) greater than 0.05 in a given sample are masked as NA on that chip. Probes that are mapped to multiple sites on hg19 are annotated as NA for chromosome and 0 for CpG/CpH coordinate.

Infinium HumanMethylation450 Bead Chip (Illumina) for Osteosarcoma (OS)

*Protocols performed at the Phoenix Children’s Hospital and Baylor College of Medicine.

Bisulfite conversion of genomic DNA was performed with EZ-96 DNA methylation Kit (Zymo Research, Irvine, CA, #D5002) following the manufacturer’s protocol with modifications for the Infinium Methylation Assay. Briefly, one microgram of genomic DNA was mixed with 5 µl of Dilution Buffer and incubated at 37°C for 15 minutes and then mixed with 100 µl of conversion reagent prepared as instructed in the protocol. Mixtures were incubated in a thermocycler for 16 cycles at 95°C for 30 seconds and 50°C for 60 minutes. Bisulfite-converted DNA samples were loaded onto the provided 96-column plates for desulphonation, washing and elution. Bisulfite-converted genomic DNA was analyzed using the Infinium Human Methylation450K Beadchip Kit (Illumina, San Diego, CA, #WG-314-1001). DNA amplification, fragmentation, array hybridization, extension and staining were performed with reagents provided in the kit according to the manufacturer’s protocol (Illumina Infinium II Methylation Assay, #WG-901-2701). Briefly, 4 µl of bisulfite-converted genomic DNA was added to 0.8 ml 96-well storage plate (Thermo Fisher Scientific), denatured in 0.014N sodium hydroxide, neutralized and then amplified for 20-24 hours at 37°C. Samples were fragmented at 37°C for 60 minutes and precipitated in isopropanol. Re-suspended samples were denatured in a 96-well plate heat block at 95°C for 20 minutes. 15 µl of each sample was loaded onto a 12-sample BeadChip, assembled in the hybridization chamber as instructed by the manufacturer and incubated at 48°C for 16-20 hours.  Following hybridization, the BeadChips were washed and assembled in a fluid flow-through station for primer-extension reaction and staining with reagents and buffers provided.

Polymer-coated BeadChips were scanned using Illumina iScan technology which outputs data in the format of IDAT files.  These are then used retrieve the probe intensities and calculate the beta-values. Raw unmethylated and methylated intensities were background corrected using out-of-band correction.

Probe intensities were then color corrected using Lumi's dye bias correction algorithm.  Beta-values were calculated from probe intensities and corrected for probe bias using the beta mixture quantile dilation (BMIQ) normalization method.

Infinium HumanMethylation450 Bead Chip (Illumina) for Wilms Tumor (WT)

*Protocol performed at Ann & Robert H. Lurie Children's Hospital.

DNA was extracted from tumor samples at Nationwide Children's BioPathology Center (BPC) by using the standard BPC protocol. The DNA samples were analyzed by Pico green to verify gDNA concentration, spectrophotometry to verify DNA purity, and gel electrophoresis to verify DNA quality.  DNA samples (1.5 ug) diluted in nuclease-free water were provided to the Northwestern University Genomics Core in 96-well plate format for Illumina 450K DNA methylation analysis. Randomly selected samples were tested by Northwestern University Genomics Core to verify that the correct concentration had been provided.                

Nucleic acid hybridization and labeling were performed according to the manufacturer's protocol for the Illumina 450K array at the Northwestern University Genomics Core Facility. Nucleic acid labeling is completed after the hybridization step with Illumina Infinium 450K arrays.   

The array scanning protocol was performed according to the manufacturer's protocol for the Illumina 450K array at the Northwestern University Genomics Core Facility.                 

Raw data files (1 red and 1 green .idat file per sample and 1 .sdf file from each array, which included 12 samples per array) were processed at the Northwestern University Genomics Core Facility by BeadStudio software. The following subtables were generated by BeadStudio and were downloaded in .txt format (Level 2 data): (1) the Sample Methylation Profile .txt, (2) the Control Profile .txt, and (3) the Control Probe Profile .txt. Several quality control steps were used for these data. Samples were subjected to the internal quality controls in the Bioconductor lumi package. Samples were subjected to a color balance check using the Bioconductor lumi package. Gender analysis was performed to increase our confidence that the data correctly corresponded to the expected sample. Because one of the X chromosomes is heavily methylated in females, the density of X-chromosome methylation is a good indicator of gender. Unsupervised hierarchical clustering based on all methylated regions on the X-chromosome was performed using average-linkage clustering with CLUSTER and the results were displayed with TREEVIEW. The tumors were assigned a gender based on their clustering in the tumor dendogram, which was checked against the known gender of the sample. The X-chromosome methylation profile corresponded to gender in the majority (~95%) of the samples; the other samples did not show gender-specific patterns.                 

The Sample Methylation Profile text, which included information for all of the samples, was broken down at the DCC into a single .txt file (level 2 data) per sample containing the following columns: (1) the sample ID, (2) probe name, (3) AVG_Beta value, (4) gene Symbol, (5) chromosome, and (6) position.

Infinium MethylationEPIC BeadChip Kit (Illumina) for Acute Leukemia of Ambiguous Lineage (ALAL)

*Protocol performed at St. Jude Children's Research Hospital.

Raw data from the Infinium MethylationEPIC BeadChip Kit (Illumina Inc.) were analyzed using the ChAMP1 R package.

In general, the raw *.idat files were imported through “minfi” method2 and then the following filters were applied to exclude the probes: 1) with detection P-value above 0.01 in one or more samples; 2) with beadcount <3 in at least 5% of samples; 3) as non-CpG probes, 4) identified as SNPs3; 5) aligned to multiple locations4 and 6) on the X or Y chromosome. After filtering, “BMIQ” normalization from ChAMP package was used as the author suggested to calculate methylation beta values. Batch effect was observed by the singular value decomposition method5 and adjusted by ComBat normalization method6.

References

  1. Morris TJ, et al. (2014). ChAMP: 450k Chip Analysis Methylation Pipeline. Bioinformatics. 30(3):428-30. (PMID: 24336642)
  2. Aryee MJ, et al. (2014). Minfi: a flexible and comprehensive Bioconductor package for the analysis of Infinium DNA methylation microarrays. Bioinformatics. 30(10):1363-9. (PMID: 24478339)
  3. Zhou W, et al. (2017). Comprehensive characterization, annotation and innovative use of Infinium DNA methylation BeadChip probes. Nucleic Acids Res. 45(4):e22. (PMID: 27924034)
  4. Nordlund J, et al. (2013). Genome-wide signatures of differential DNA methylation in pediatric acute lymphoblastic leukemia. Genome Biol. 14(9):r105. (PMID: 24063430)
  5. Teschendorff AE, et al. (2009). An epigenetic signature in peripheral blood predicts active ovarian cancer. PloS One. 4(12):e8274. (PMID: 20019873)
  6. Johnson WE, et al. (2007). Adjusting batch effects in microarray expression data using empirical Bayes methods. Biostatistics. 8(1):118-27. (PMID: 16632515)
Last updated: May 30, 2019