- Research article
- Open Access
Evaluation of microarray-based DNA methylation measurement using technical replicates: the Atherosclerosis Risk In Communities (ARIC) Study
© Bose et al.; licensee BioMed Central Ltd. 2014
Received: 27 January 2014
Accepted: 8 September 2014
Published: 19 September 2014
DNA methylation is a widely studied epigenetic phenomenon; alterations in methylation patterns influence human phenotypes and risk of disease. As part of the Atherosclerosis Risk in Communities (ARIC) study, the Illumina Infinium HumanMethylation450 (HM450) BeadChip was used to measure DNA methylation in peripheral blood obtained from ~3000 African American study participants. Over 480,000 cytosine-guanine (CpG) dinucleotide sites were surveyed on the HM450 BeadChip. To evaluate the impact of technical variation, 265 technical replicates from 130 participants were included in the study.
For each CpG site, we calculated the intraclass correlation coefficient (ICC) to compare variation of methylation levels within- and between-replicate pairs, ranging between 0 and 1. We modeled the distribution of ICC as a mixture of censored or truncated normal and normal distributions using an EM algorithm. The CpG sites were clustered into low- and high-reliability groups, according to the calculated posterior probabilities. We also demonstrated the performance of this clustering when applied to a study of association between methylation levels and smoking status of individuals. For the CpG sites showing genome-wide significant association with smoking status, most (~96%) were seen from sites in the high reliability cluster.
We suggest that CpG sites with low ICC may be excluded from subsequent association analyses, or extra caution needs to be taken for associations at such sites.
DNA methylation is one of the most commonly occurring epigenetic phenomena in the human genome. It is one of the major regulators of gene transcription and plays a vital role in many cellular processes. In the last decade, numerous studies have shown that abnormal methylation patterns are linked to phenotypic differences and development of disease [1–6].
Recent technological advances have provided multiple platforms for systematically interrogating DNA methylation variation across the genome . Among them, the Illumina HumanMethylation450 BeadChip (HM450) (Illumina, Inc.) is a new-generation array constituting a major extension of the previous Infinium HumanMethylation27 BeadChip (HM27) (Illumina, Inc.) and can be used to assess the methylation condition of more than 480,000 cytosines distributed over the entire genome. Sandoval et al.  validated the HM450 chip by showing that methylation patterns measured in colorectal cancer cell lines and normal mucosa were consistent with those found in bisulfite genomic sequencing. Other recent studies have extensively evaluated data generated from this chip and developed data processing and analysis pipelines [9–11].
Similar to other microarray experiments (e.g. RNA expression), it is important to evaluate the impact of technical variation in the measurement. A well-known source of bias for epigenome-wide association studies (EWASs) of DNA methylation is the so-called “batch effect”, which is largely caused by technical differences from one chip to the next. Although many statistical methods have been proposed to correct for batch effects [12–17], no approach has been generally accepted. Alternatively, it will be useful to consider statistical measures that can quantify the extent to which the measured methylation level at a specific CpG site is affected by technical errors. For CpG sites with large inter-individual variation, it may be reasonable to assume that technical differences will have relatively low impact at these sites. However, the HM450 chip contains a large number of CpG sites with little inter-individual variation in methylation levels, for which it is critical to consider additional statistics to evaluate the impact of technical errors.
In experiments, technical replicates are often included which can be used to evaluate the consistency of measurement. Meng et al.  suggested using technical replicates to identify and exclude “non-variable” CpG sites, at which the technical “noise” outperforms true biological variation. Their data came from the Illumina GoldenGate methylation array and consisted of 311 samples assayed at 1,505 sites. However, no comprehensive work has been done for the HM450 chip.
In the Atherosclerosis Risk in Communities (ARIC) study, the HM450 chip was used for genome wide profiling of DNA methylation in 2,873 African-American individuals assayed at 485,577 CpG sites, with 130 replicate pairs or triplets included. We evaluated the consistency of methylation measurement at each CpG site based on the intraclass correlation coefficient (ICC), which compares the within- and between-replicate variations. We next modeled the ICC values using a mixture model approach to classify CpG sites on the HM450 chip into high and low reliability clusters. We demonstrated the relationship of this classification to the results of an association study between methylation level and smoking status of ARIC study participants.
The aim of this work was to apply a statistical procedure to identify CpG sites which give consistent values for methylation levels in DNA extracted from peripheral blood using the HM450 chip (or a similar methylation profiling array), thereby providing insight into the performance of this chip across CpG sites. We hope that our results will facilitate subsequent EWAS analyses.
Based on the ICC, we quantified the reliability of methylation measurements on the HM450 chip using technical replicates included in the analysis of DNA methylation in ARIC study participants. The CpG sites were then classified into low and high reliability groups by modeling the distribution of the ICC. We further demonstrated the relationship of measurement reliability to the associations between methylation and smoking status in the ARIC study.
Distribution of ICC
Classification of CpG sites using mixture models
Fitted models for ICC distribution
Censored normal mixture
Truncated normal mixture
Maximized log likelihood
Data application: association with smoking status
Next, we investigated how association patterns between phenotypes and methylation levels vary across different reliability clusters using smoking status from the ARIC Study as an illustrative example. Complete information on methylation levels, smoking status, and covariates was available for 2,835 African American participants. The average age was 56.7 years (standard deviation, SD, 5.9 years), 63.4% were female, and 25.6% were current smokers. 874 CpG methylation sites displayed genome wide significant association (p < 10-7), with current smoking status, or 723 CpG sites after excluding those with at least one common SNP within 50 base pairs. Among the 723 associated sites, only 31 were from the low reliability cluster, and none from the zero-ICC cluster.
Recent technological advances have provided high-throughput array for systematically interrogating DNA methylation across the genome. This allows investigators to evaluate regions of the genome in which variation in DNA methylation may influence gene expression and ultimately disease risk. One of the potential problems that could threaten the validity of methylation-phenotype association results is the variability of technical measurements (e.g., batch effects). Here we used technical replicates to assess reliability of methylation measurement on CpG sites interrogated by the HM450 chip. Using a mixture model approach, we classified CpG sites into low and high reliability groups. The ICC value can serve as a measure to quantify the impact of technical errors, especially batch effects. In addition, CpG sites with low reliability could be excluded in subsequent analyses to improve the power of EWASs by reducing the burden of multiple hypothesis tests.
In our analysis of 265 technical replicates available as part of the ARIC methylation data, we classified CpG sites on the HM450 chip by a mixture of truncated normal and normal distributions, with ~44% sites in the cluster of high reliability. For CpG sites with high ICC values, the variation between independent biological samples is relatively large compared to the variation due to technical error. We expect statistical tests performed for the sites with high ICC will have more power to detect association between methylation level and phenotype of interest than for the sites with low ICC. In the follow-up EWAS of smoking status, we showed that a majority (96%) of genome-wide significant associations are from the high reliability cluster, even though only 44% of all sites were in this cluster.
The ICC provides a measure of reliability for the methylation measurement. In EWASs, it is well known that batch effects can threaten the validity of association results . Commonly used methods for batch effects correction include empirical Bayes methods , surrogate variables [15, 16], and linear mixed-effects (LMM) models. However, it is difficult to completely control or remove batch effects in association tests. Alternatively, we can use ICC to quantify to what extent the batch effects can affect the methylation measures. The probes for CpG sites with low ICC are expected to be more vulnerable to technical variability including batch effects. Similar to the ICC for replicates, we calculated ICC for chip effects (ICCchip), which is the proportion of variation in measured methylation levels due to chip difference in the experiment. Additional file 1: Figure S2 shows that when ICCchip is close to 1, most of the variation in methylation measures is due to technical “errors”. Reliability of methylation measures is therefore low at these CpG sites, with ICC for replicates close to 0. On the other hand, when ICC for replicates is close to 1, technical variability has minimal contribution to the methylation measures, and little batch effects are observed (ICCchip ≈ 0).
In addition to the EWAS of smoking status using LMM, we performed the EWAS using a linear regression model without adjusting for chip and chip position (results not shown). There were 1,783 CpG sites significantly associated with smoking status, after excluding sites with low pass rate or with at least one common SNP within 50 base pairs. Among these sites, there was no site in the zero-ICC cluster, 253 sites from the low reliability cluster, and 1,530 sites from the high reliability cluster. All except 9 sites identified by the LMM were also genome-wide significant in this model. The number of significant sites without adjustment for chip and chip position was approximately 2.5-fold higher than we observed in the EWAS with those adjustments. We hypothesize that some of the additional associations are false positives due to technical errors. We also compared our results to published EWASs of smoking studies [23–27], which were carried out in individuals of European and African ancestry. We were able to validate ~80% of the reported signals with genome-wide significance. Of these, 32 sites were found within the body of the aryl-hydrocarbon receptor repressor (AHRR) gene, with a 15% decrease in the mean methylation level of cp05575921 in current smokers compared to non-smokers (p = 6.2 × 10-196) as previously shown in the other studies.
The ICC value is determined by two variance components, for variance between independent samples and for variance within replicates of the same sample. When sample characteristics (e.g., cell type composition, sample ethnicity, or genetic variants) are controlled for, becomes smaller which leads to smaller ICC than the current estimates (in Additional file 1: Figure S3). When batch effects are controlled for, both and can be reduced and the ICC can increase or decrease at different CpG sites (in Additional file 1: Figure S4).
The empirical distribution of C we observed is similar to that of Meng et al. , in which an earlier version of the Illumina GoldenGate array having 1,505 CpG sites was used. In their analysis, Pearson product–moment correlation was primarily used between replicate pairs, although ICC was also introduced. Our analysis updates the information for the HM450 chip, and focuses on ICC which allows more than two replicates per biological sample.
Meng et al. observed that, compared to CpG sites not in CpG islands, a higher proportion of sites in CpG islands showed low inter-individual variation in methylation levels relative to within-individual variation (“non-variable”), consistent with a previous finding by Bock et al. . Our results further confirmed this observation. In contrast, Dedeurwaerder et al.  studied technical replicates of a human colon cancer cell line and showed that the probe-wise variation of methylation levels were greater for Infinium II probes than those for Infinium I probes. However, their findings were limited to within-individual variation, with no comparison to inter-individual variation. In addition, given the small number of replicates, the impact of batch effects was minimal in their study. Therefore the results were not directly comparable.
From the distribution of ICC values we noted that the distribution of high ICC component is slightly left skewed, with few CpG sites having ICC close to 1. Because of the normal distribution assumption for this component in the mixture model, this skewness may lead to an underestimate for the corresponding variance parameter and under-fit the distribution of ICC in the middle (0.2-0.5), in both of the truncated and censored normal mixture models. The direct consequence is that CpG sites with moderate ICC values will have a small posterior probability of being in the high reliability cluster. As a possible remedy, we considered three-component mixture models, with an additional normal distribution component to model the moderate ICC values. However, the model fit was poorer than the 2-component models (results not shown).
To classify the CpG sites into different reliability groups, Meng et al.  chose a cutoff of 0.5 for the posterior probability and suggested to exclude the “non-variable” sites in the low reliability group. In contrast, we set a rather small cutoff of 0.01 to include more sites with moderate ICC into the high reliability cluster, given the discussion above. In the ARIC study, the low reliability cluster contributes only ~4% of significant associations in the EWAS of smoking status, with only two associated sites having ICC < 0.1. We believe that validation using other approaches is especially important for sites with low ICC, for example, bisulphite pyrosequencing and replication in an external validation cohort. In this specific example, restricting analysis to CpG sites in the high reliability cluster would result in the loss of only a few significant association signals and could improve the power of EWAS by reducing the number of hypothesis tests by more than a half. To demonstrate, we calculated the statistical power at different effect sizes (%variance explained), based on a sample size of 2,500 and genome-wide significance level determined using Bonferroni correction. We assume 50% CpG sites will be excluded due to low ICC. We observe power increase of ~5% under this setting (in Additional file 1: Figure S5). For other studies, investigators may decide their own cutoff for the ICC. Even if all CpG sites on the chip are to be included in the association analysis, the ICC can still serve as a measure of how likely the association results are to be affected by technical errors.
Here, we aimed to describe the distribution of ICC using a statistical approach. Our results clearly approximate the observed distribution. However, it is unknown what biological or chemical factors explain the difference in performance of the corresponding probes on the chip between the low (including zero-ICC) and high reliability clusters. It is also unknown whether our results can be generalized to other populations or other types of tissues. The methylation profiling in ARIC samples are from white blood cells, which is a mixture of multiple cell types, and the participants are African Americans. It has been shown that both cell type composition  and sample ethnicity  can heavily influence DNA methylation pattern. It is unknown whether our results will still hold when a purified cell type, e.g., monocytes or neutrophils, or different tissue is used, or when samples from different ethnic groups are included. As an example, Etcheverry et al.  carried out an EWAS to identify CpG sites differentially methylated between glioblastoma and control samples using brain tissues. Among the 616 CpG sites identified, 584 are interrogated on the HM450 BeadChip. We examined the ICC values for these 584 sites in our data, and found 486 sites to belong to the high reliability group (ICC > 0.37), 88 to belong to the low reliability group, and 10 to belong to the zero-ICC group. Enrichment of significant associations in the high-ICC group in that study suggests that our results may have reasonable generality.
In addition, we used the “raw” methylation measures produced by the GenomeStudio software with minimal normalization. Other normalization procedures or use of M-values  instead of β-values may result in different ICC values. However, the same analytical approach used here can be applied in other studies, if technical replicates are included.
Another limitation of the study is that we did not take into account the estimation uncertainty when estimating the ICC, which can lead to different confidence intervals for this reliability measure. We will consider a weighted estimation approach in future studies, and re-classify the CpG sites on the HM450 chip.
We examined the reliability of methylation measurements using the latest HM450 chip, and demonstrated that the CpG sites assayed on this array can be potentially classified according to different levels of reliability. We also evaluated the impact of measurement reliability on results of EWAS. The biological differences between the clusters of CpG sites need to be further investigated. The estimated ICC values for all CpG sites on the HM450 chip are available in Additional file 2. We hope that our results can provide additional guidance on inclusion/exclusion of CpG sites for future EWASs using the HM450 chip, and our analysis approach can be generalized to other types of methylation arrays.
To evaluate the reliability of methylation measurement assayed on the HM450 chip, an intraclass correlation coefficient (ICC) can be calculated for each CpG site using the technical replicates. We aim to classify the CpG sites into multiple reliability groups, by modeling the distribution of ICC values using a mixture distribution.
The ARIC methylation data
The Atherosclerosis Risk in Communities (ARIC) Study is a prospective cohort study of cardiovascular disease risk in four U.S. communities . Between 1987 and 1989, 7,082 men and 8,710 women aged 45–64 years were recruited from Forsyth County, North Carolina; Jackson, Mississippi (African Americans only); suburban Minneapolis, Minnesota; and Washington County, Maryland. The ARIC Study protocol was approved by the institutional review board of each participating university. Participants underwent a baseline clinical examination (Visit 1) and four subsequent follow-up clinical exams (Visits 2–5). Written informed consent was obtained for each clinic exam, including that for genetic studies.
Genomic DNA was extracted from peripheral blood leukocyte samples using the Gentra Puregene Blood Kit (Qiagen) according to the manufacturer’s instructions (http://www.qiagen.com). Bisulphite conversion of 1 ug genomic DNA was performed using the EZ-96 DNA Methylation Kit (Deep Well Format) (Zymo Research) according to the manufacturer’s instructions (http://www.zymoresearch.com). Bisulphite conversion efficiency was determined by PCR amplification of the converted DNA before proceeding with methylation analyses on the Illumina platform using Zymo Research’s Universal Methylated Human DNA Standard and Control Primers.
Bisulphite-converted DNA from 2,905 African American participants at Visit 2 (1990–92; n = 2,504) or Visit 3 (1993–95; n = 441) were measured for methylation status using the HM450 chip. The degree of methylation is determined using Illumina GenomeStudio 2011.1, Methylation module 1.9.0 software. The methylation score for each CpG was represented as a beta (β) value calculated by dividing the fluorescence intensity of the methylated allele by the sum of the intensities of the methylated allele and unmethylated allele. β-values may take any value between 0 (non-methylated) and 1 (completely methylated). Background subtraction was conducted with the GenomeStudio software using built-in negative control bead types on the array. An average normalization was also applied in GenomeStudio to minimize scanner-to-scanner variation using ~90 normalization probe pairs included on the array which target housekeeping regions with no underlying CpG sites. These probes are used to independently calculate normalization values in the green and red channels so that all samples have the same average intensity.
Individuals (n = 32) were excluded from analysis if a pass rate for the DNA sample for the study participant was less than 99% (probes with a detection p-value >0.01/all probes on the array). Probes on the HM450 chip for which the pass rate was less than 99% (sample with a detection p-value >0.01 at probe/all samples) were not analyzed (n = 11,789).
Technical replicates were included for 130 samples (total n = 265 with 5 samples replicated 3 times), from which the ICC values are calculated for all probes. All except two of the replicate pairs were distributed on different chips.
ICC estimation and modeling
which takes values in [0,1] and measures the extent to which the measurements in a replicate-set resemble each other. An ICC of 1 indicates perfect measurement accuracy at the CpG site, and 0 implies little variation between independent samples but large measurement error between technical replicates. Thus, we can classify sites into clusters of varying degrees of reliability based on ICC values.
where l(ij) indexes the batch (chip) where replicate i of sample j is located so bl(ij) is the corresponding batch effect, and X jp is the pth covariate (characteristic of biological sample), such as age or sex, for sample j. Comparing this model to the random-effects ANOVA model above, the within replicate set variation is decomposed into variances of batch-specific measurement error (bl(ij)) and other technical error (e ij ), and corresponds to the between replicate set (biological sample) variation. When the measure of methylation is vulnerable to batch effects at a given CpG site, the variance of bl(ij) is large, leading to small ICC.
where δ i = 1 for the censored observations (i.e., ICC = 0) and 0 otherwise, p is the mixing proportion for the first cluster, the set of parameters, and Φ(⋅) and ϕ(⋅) denote the cumulative and probability distribution function of standard normal distribution, respectively.
where p0 and n0 denote the sample proportion and number of sites with ICC = 0, respectively, and p1 = p(1 - p0) and p2(1 - p)(1 - p0).
The two models were then compared using the maximized log likelihood values. The CpG sites interrogated by the HM450 chip can be classified based on the corresponding posterior probabilities, , computed from either of the two models.
Data application: association with smoking status
To demonstrate how association patterns between methylation levels at CpG sites and phenotypic characteristics of individuals vary across the clustered sites, we investigated the cross-sectional association of methylation levels and smoking status in the ARIC Study. We fit a linear mixed-effects (LMM) model of methylation on self-reported smoking status (never/former smoker vs. current smoker), adjusting for age, sex, body mass index, self-reported alcohol consumption (never/former vs. current drinker), the top principal components of ancestry derived from genotype data, visit number, field center, and chip row position as fixed covariates. The values of the covariates were from the specific visit at which the participant’s sample was obtained. The batch (chip) effects were modeled as random effects. Probes with at least one single nucleotide polymorphism (SNP) within 50 base pairs of the CpG site and minor allele frequency (MAF) > .05 based on all samples in the 1000 Genomes project (http://www.1000genomes.org/) were excluded . We used p < 10-7 as genome-wide significance based on a Bonferroni correction, and summarize the number of significant associations by the low and high reliability group, classified as above.
The Atherosclerosis Risk in Communities Study is carried out as a collaborative study supported by National Heart, Lung, and Blood Institute contracts (HHSN268201100005C, HHSN268201100006C, HHSN268201100007C, HHSN268201100008C, HHSN268201100009C, HHSN268201100010C, HHSN268201100011C, and HHSN268201100012C). The authors thank the staff and participants of the ARIC study for their important contributions. Methylation profiling of the ARIC samples was also supported by the National Institutes of Health (NIH) American Recovery and Reinvestment Act of 2009 (ARRA) Building on GWAS for NHLBI-diseases: the U.S. CHARGE consortium (5RC2HL102419) (PI: E. Boerwinkle). This work was carried out in part using computing resources at the University of Minnesota Supercomputing Institute. The authors have no conflict of interest.
- Jones PA, Takai D: The role of DNA methylation in mammalian epigenetics. Science. 2001, 293 (5532): 1068-1070. 10.1126/science.1063852.View ArticlePubMedGoogle Scholar
- Robertson KD: DNA methylation and human disease. Nat Rev Genet. 2005, 6 (8): 597-610.View ArticlePubMedGoogle Scholar
- Scarano MI, Strazzullo M, Matarazzo MR, D’Esposito M: DNA methylation 40 years later: Its role in human health and disease. J Cell Physiol. 2005, 204 (1): 21-35. 10.1002/jcp.20280.View ArticlePubMedGoogle Scholar
- Heyn H, Esteller M: DNA methylation profiling in the clinic: applications and challenges. Nat Rev Genet. 2012, 13 (10): 679-692. 10.1038/nrg3270.View ArticlePubMedGoogle Scholar
- Qiu P, Zhang L: Identification of markers associated with global changes in DNA methylation regulation in cancers. BMC Bioinformatics. 2012, 13 (Suppl 13): S7-10.1186/1471-2105-13-S13-S7.View ArticlePubMed CentralPubMedGoogle Scholar
- Liu J, Chen J, Ehrlich S, Walton E, White T, Perrone-Bizzozero N, Bustillo J, Turner JA, Calhoun VD: Methylation patterns in whole blood correlate with symptoms in schizophrenia patients. Schizophr Bull. 2013, sbt080-Google Scholar
- Laird PW: Principles and challenges of genomewide DNA methylation analysis. Nat Rev Genet. 2010, 11 (3): 191-203. 10.1038/nrg2732.View ArticlePubMedGoogle Scholar
- Sandoval J, Heyn H, Moran S, Serra-Musach J, Pujana MA, Bibikova M, Esteller M: Validation of a DNA methylation microarray for 450,000 CpG sites in the human genome. Epigenetics. 2011, 6 (6): 692-702. 10.4161/epi.6.6.16196.View ArticlePubMedGoogle Scholar
- Marabita F, Almgren M, Lindholm ME, Ruhrmann S, Fagerstrom-Billai F, Jagodic M, Sundberg CJ, Ekstrom TJ, Teschendorff AE, Tegner J, Gomez-Cabrero D: An evaluation of analysis pipelines for DNA methylation profiling using the illumina human Methylation450 BeadChip platform. Epigenetics. 2013, 8 (3): 333-346. 10.4161/epi.24008.View ArticlePubMed CentralPubMedGoogle Scholar
- Dedeurwaerder S, Defrance M, Calonne E, Denis H, Sotiriou C, Fuks F: Evaluation of the Infinium Methylation 450 K technology. Epigenomics. 2011, 3 (6): 771-784. 10.2217/epi.11.105.View ArticlePubMedGoogle Scholar
- Touleimat N, Tost J: Complete pipeline for Infinium((R)) Human Methylation 450 K BeadChip data processing using subset quantile normalization for accurate DNA methylation estimation. Epigenomics. 2012, 4 (3): 325-341. 10.2217/epi.12.21.View ArticlePubMedGoogle Scholar
- Sun Z, Chai HS, Wu Y, White WM, Donkena KV, Klein CJ, Garovic VD, Therneau TM, Kocher JP: Batch effect correction for genome-wide methylation data with Illumina Infinium platform. BMC Med Genet. 2011, 4: 84-Google Scholar
- Chen C, Grennan K, Badner J, Zhang D, Gershon E, Jin L, Liu C: Removing batch effects in analysis of expression microarray data: an evaluation of six batch adjustment methods. PLoS One. 2011, 6 (2): e17238-10.1371/journal.pone.0017238.View ArticlePubMed CentralPubMedGoogle Scholar
- Leek JT, Scharpf RB, Bravo HC, Simcha D, Langmead B, Johnson WE, Geman D, Baggerly K, Irizarry RA: Tackling the widespread and critical impact of batch effects in high-throughput data. Nat Rev Genet. 2010, 11 (10): 733-739. 10.1038/nrg2825.View ArticlePubMedGoogle Scholar
- Leek JT, Johnson WE, Parker HS, Jaffe AE, Storey JD: The sva package for removing batch effects and other unwanted variation in high-throughput experiments. Bioinformatics. 2012, 28 (6): 882-883. 10.1093/bioinformatics/bts034.View ArticlePubMed CentralPubMedGoogle Scholar
- Teschendorff AE, Zhuang J, Widschwendter M: Independent surrogate variable analysis to deconvolve confounding factors in large-scale microarray profiling studies. Bioinformatics. 2011, 27 (11): 1496-1505. 10.1093/bioinformatics/btr171.View ArticlePubMedGoogle Scholar
- Johnson WE, Li C, Rabinovic A: Adjusting batch effects in microarray expression data using empirical Bayes methods. Biostatistics. 2007, 8 (1): 118-127. 10.1093/biostatistics/kxj037.View ArticlePubMedGoogle Scholar
- Meng H, Joyce AR, Adkins DE, Basu P, Jia Y, Li G, Sengupta TK, Zedler BK, Murrelle EL, van den Oord EJ: A statistical method for excluding non-variable CpG sites in high-throughput DNA methylation profiling. BMC Bioinformatics. 2010, 11: 227-10.1186/1471-2105-11-227.View ArticlePubMed CentralPubMedGoogle Scholar
- Etcheverry A, Aubry M, de Tayrac M, Vauleon E, Boniface R, Guenot F, Saikali S, Hamlat A, Riffaud L, Menei P, Quillien V, Mosser J: DNA methylation in glioblastoma: impact on gene expression and clinical outcome. BMC Genomics. 2010, 11: 701-10.1186/1471-2164-11-701.View ArticlePubMed CentralPubMedGoogle Scholar
- Bibikova M, Le J, Barnes B, Saedinia-Melnyk S, Zhou L, Shen R, Gunderson KL: Genome-wide DNA methylation profiling using Infinium(R) assay. Epigenomics. 2009, 1 (1): 177-200. 10.2217/epi.09.14.View ArticlePubMedGoogle Scholar
- Chen YA, 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-209. 10.4161/epi.23470.View ArticlePubMed CentralPubMedGoogle Scholar
- Barfield RT, Almli LM, Kilaru V, Smith AK, Mercer KB, Duncan R, Klengel T, Mehta D, Binder EB, Epstein MP, Ressler KJ, Conneely KN: Accounting for population stratification in DNA methylation studies. Genet Epidemiol. 2014, 38 (3): 231-241. 10.1002/gepi.21789.View ArticlePubMed CentralPubMedGoogle Scholar
- Breitling LP, Yang R, Korn B, Burwinkel B, Brenner H: Tobacco-smoking-related differential DNA methylation: 27 K discovery and replication. Am J Hum Genet. 2011, 88 (4): 450-457. 10.1016/j.ajhg.2011.03.003.View ArticlePubMed CentralPubMedGoogle Scholar
- Wan ES, Qiu W, Baccarelli A, Carey VJ, Bacherman H, Rennard SI, Agusti A, Anderson W, Lomas DA, Demeo DL: Cigarette smoking behaviors and time since quitting are associated with differential DNA methylation across the human genome. Hum Mol Genet. 2012, 21 (13): 3073-3082. 10.1093/hmg/dds135.View ArticlePubMed CentralPubMedGoogle Scholar
- Joubert BR, Haberg SE, Nilsen RM, Wang X, Vollset SE, Murphy SK, Huang Z, Hoyo C, Midttun O, Cupul-Uicab LA, Ueland PM, Wu MC, Nystad W, Bell DA, Peddada SD, London SJ: 450 K epigenome-wide scan identifies differential DNA methylation in newborns related to maternal smoking during pregnancy. Environ Health Perspect. 2012, 120 (10): 1425-1431. 10.1289/ehp.1205412.View ArticlePubMed CentralPubMedGoogle Scholar
- Shenker NS, Polidoro S, van Veldhoven K, Sacerdote C, Ricceri F, Birrell MA, Belvisi MG, Brown R, Vineis P, Flanagan JM: Epigenome-wide association study in the European Prospective Investigation into Cancer and nutrition (EPIC-Turin) identifies novel genetic loci associated with smoking. Hum Mol Genet. 2013, 22 (5): 843-851. 10.1093/hmg/dds488.View ArticlePubMedGoogle Scholar
- Monick MM, Beach SR, Plume J, Sears R, Gerrard M, Brody GH, Philibert RA: Coordinated changes in AHRR methylation in lymphoblasts and pulmonary macrophages from smokers. American journal of medical genetics Part B, Neuropsychiatric genetics: the official publication of the International Society of Psychiatric Genetics. 2012, 159B (2): 141-151. 10.1002/ajmg.b.32021.View ArticleGoogle Scholar
- Bock C, Walter J, Paulsen M, Lengauer T: Inter-individual variation of DNA methylation and its implications for large-scale epigenome mapping. Nucleic Acids Res. 2008, 36 (10): e55-10.1093/nar/gkn122.View ArticlePubMed CentralPubMedGoogle Scholar
- Jaffe AE, Irizarry RA: Accounting for cellular heterogeneity is critical in epigenome-wide association studies. Genome Biol. 2014, 15 (2): R31-10.1186/gb-2014-15-2-r31.View ArticlePubMed CentralPubMedGoogle Scholar
- Irizarry RA, Ladd-Acosta C, Carvalho B, Wu H, Brandenburg SA, Jeddeloh JA, Wen B, Feinberg AP: Comprehensive high-throughput arrays for relative methylation (CHARM). Genome Res. 2008, 18 (5): 780-790. 10.1101/gr.7301508.View ArticlePubMed CentralPubMedGoogle Scholar
- The ARIC investigators: The Atherosclerosis Risk in Communities (ARIC) Study: design and objectives. Am J Epidemiol. 1989, 129 (4): 687-702.Google Scholar
- Dempster A, Laird N, Rubin D: Maximum likelihood from incomplete data via the E-M algorithm. Journal of the Royal Statistical Society (B). 1977, 39: 1-38.Google Scholar
- Lee G, Scott C: EM algorithms for multivariate Gaussian mixture models with truncated and censored data. Computational Statistics & Data Analysis. 2012, 56 (9): 2816-2829. 10.1016/j.csda.2012.03.003.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. 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.