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.     

MB-COMT promoter DNA methylation is associated with working-memory processing in schizophrenia patients and healthy controls

Many genetic studies report mixed results both for the associations between COMT polymorphisms and schizophrenia and for the effects of COMT variants on common intermediate phenotypes of the disorder. Reasons for this may include small genetic effect sizes and the modulation of environmental influences. To improve our understanding of the role of COMT in the disease etiology, we investigated the effect of DNA methylation in the MB-COMT promoter on neural activity in the dorsolateral prefrontal cortex during working memory processing as measured by fMRI - an intermediate phenotype for schizophrenia. Imaging and epigenetic data were measured in 102 healthy controls and 82 schizophrenia patients of the Mind Clinical Imaging Consortium (MCIC) study of schizophrenia. Neural activity during the Sternberg Item Recognition Paradigm was acquired with either a 3T Siemens Trio or 1.5T Siemens Sonata and analyzed using the FMRIB Software Library (FSL). DNA methylation measurements were derived from cryo-conserved blood samples. We found a positive association between MB-COMT promoter methylation and neural activity in the left dorsolateral prefrontal cortex in a model using a region-of-interest approach and could confirm this finding in a whole-brain model. This effect was independent of disease status. Analyzing the effect of MB-COMT promoter DNA methylation on a neuroimaging phenotype can provide further evidence for the importance of COMT and epigenetic risk mechanisms in schizophrenia. The latter may represent trans-regulatory or environmental risk factors that can be measured using brain-based intermediate phenotypes.     

MB-COMT promoter DNA methylation is associated with working-memory processing in schizophrenia patients and healthy controls

Many genetic studies report mixed results both for the associations between COMT polymorphisms and schizophrenia and for the effects of COMT variants on common intermediate phenotypes of the disorder. Reasons for this may include small genetic effect sizes and the modulation of environmental influences. To improve our understanding of the role of COMT in the disease etiology, we investigated the effect of DNA methylation in the MB-COMT promoter on neural activity in the dorsolateral prefrontal cortex during working memory processing as measured by fMRI - an intermediate phenotype for schizophrenia. Imaging and epigenetic data were measured in 102 healthy controls and 82 schizophrenia patients of the Mind Clinical Imaging Consortium (MCIC) study of schizophrenia. Neural activity during the Sternberg Item Recognition Paradigm was acquired with either a 3T Siemens Trio or 1.5T Siemens Sonata and analyzed using the FMRIB Software Library (FSL). DNA methylation measurements were derived from cryo-conserved blood samples. We found a positive association between MB-COMT promoter methylation and neural activity in the left dorsolateral prefrontal cortex in a model using a region-of-interest approach and could confirm this finding in a whole-brain model. This effect was independent of disease status. Analyzing the effect of MB-COMT promoter DNA methylation on a neuroimaging phenotype can provide further evidence for the importance of COMT and epigenetic risk mechanisms in schizophrenia. The latter may represent trans-regulatory or environmental risk factors that can be measured using brain-based intermediate phenotypes.     

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.     

Ultra-sensitive immunodetection of 5'methyl cytosine for DNA methylation analysis on oligonucleotide microarrays.

For the determination of methylation levels in genomic regulatory DNA sequences a high-sensitive assay for detecting 5'methyl-cytosines (5'mC) in non-bisulfite-treated DNA has been established. The system is designed for the application of immunofluorescence using a monoclonal antibody that specifically recognizes 5'mC in single-stranded DNA hybridized to oligonucleotide microarrays. For assay readout an ultra-sensitive fluorescence scanner with submicrometer resolution was used. To minimize autofluorescence 150-microm thin glass slides with an aldehyde-functionalized surface were developed. These methodological improvements allowed the detection of 5'mC in synthetic oligonucleotides hybridized to microarrays with atto molar analytical sensitivity. Using enzymatic fragmented genomic DNA from myeloid leukemia tumor cell lines differences in the methylation status of gene regulatory sequences for E-cadherin, p15/CDKN2b and p16/CDKN2a were demonstrated. Thus, this novel technique can potentially be used for DNA methylation analysis in various scientific fields.     

Agglomerates of aberrant DNA methylation are associated with toxicant-induced malignant transformation

Epigenetic dysfunction is a known contributor in carcinogenesis, and is emerging as a mechanism involved in toxicant-induced malignant transformation for environmental carcinogens such as arsenicals or cadmium. In addition to aberrant DNA methylation of single genes, another manifestation of epigenetic dysfunction in cancer is agglomerative DNA methylation, which can participate in long-range epigenetic silencing that targets many neighboring genes and has been shown to occur in several types of clinical cancers. Using in vitro model systems of toxicant-induced malignant transformation, we found hundreds of aberrant DNA methylation events that emerge during malignant transformation, some of which occur in an agglomerative fashion. In an arsenite-transformed prostate epithelial cell line, the protocadherin (PCDH), HOXC and HOXD gene family clusters are targeted for agglomerative DNA methylation. The agglomerative DNA methylation changes induced by arsenicals appear to be common and clinically relevant events, since they occur in other human cancer cell lines and models of malignant transformation, as well as clinical cancer specimens. Aberrant DNA methylation in general occurred more often within histone H3 lysine-27 trimethylation stem cell domains. We found a striking association between enrichment of histone H3 lysine-9 trimethylation stem cell domains and toxicant-induced agglomerative DNA methylation, suggesting these epigenetic modifications may become aberrantly linked during malignant transformation. In summary, we found an association between toxicant-induced malignant transformation and agglomerative DNA methylation, which lends further support to the hypothesis that epigenetic dysfunction plays an important role in toxicant-induced malignant transformation.     

Agglomerates of aberrant DNA methylation are associated with toxicant-induced malignant transformation

Epigenetic dysfunction is a known contributor in carcinogenesis, and is emerging as a mechanism involved in toxicant-induced malignant transformation for environmental carcinogens such as arsenicals or cadmium. In addition to aberrant DNA methylation of single genes, another manifestation of epigenetic dysfunction in cancer is agglomerative DNA methylation, which can participate in long-range epigenetic silencing that targets many neighboring genes and has been shown to occur in several types of clinical cancers. Using in vitro model systems of toxicant-induced malignant transformation, we found hundreds of aberrant DNA methylation events that emerge during malignant transformation, some of which occur in an agglomerative fashion. In an arsenite-transformed prostate epithelial cell line, the protocadherin (PCDH), HOXC and HOXD gene family clusters are targeted for agglomerative DNA methylation. The agglomerative DNA methylation changes induced by arsenicals appear to be common and clinically relevant events, since they occur in other human cancer cell lines and models of malignant transformation, as well as clinical cancer specimens. Aberrant DNA methylation in general occurred more often within histone H3 lysine-27 trimethylation stem cell domains. We found a striking association between enrichment of histone H3 lysine-9 trimethylation stem cell domains and toxicant-induced agglomerative DNA methylation, suggesting these epigenetic modifications may become aberrantly linked during malignant transformation. In summary, we found an association between toxicant-induced malignant transformation and agglomerative DNA methylation, which lends further support to the hypothesis that epigenetic dysfunction plays an important role in toxicant-induced malignant transformation.     

Distinctive patterns of DNA methylation associated with Parkinson disease

Parkinson disease (PD) is a multifactorial neurodegenerative disorder with high incidence in the elderly, where environmental and genetic factors are involved in etiology. In addition, epigenetic mechanisms, including deregulation of DNA methylation have been recently associated to PD. As accurate diagnosis cannot be achieved pre-mortem, identification of early pathological changes is crucial to enable therapeutic interventions before major neuropathological damage occurs. Here we investigated genome-wide DNA methylation in brain and blood samples from PD patients and observed a distinctive pattern of methylation involving many genes previously associated to PD, therefore supporting the role of epigenetic alterations as a molecular mechanism in neurodegeneration. Importantly, we identified concordant methylation alterations in brain and blood, suggesting that blood might hold promise as a surrogate for brain tissue to detect DNA methylation in PD and as a source for biomarker discovery.     

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.     

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.     

Identification of a new locus and validation of previously reported loci showing differential methylation associated with smoking. The REGICOR study

Smoking increases the risk of many diseases and could act through changes in DNA methylation patterns. The aims of this study were to determine the association between smoking and DNA methylation throughout the genome at cytosine-phosphate-guanine (CpG) site level and genomic regions. A discovery cross-sectional epigenome-wide association study nested in the follow-up of the REGICOR cohort was designed and included 645 individuals. Blood DNA methylation was assessed using the Illumina HumanMethylation450 BeadChip. Smoking status was self-reported using a standardized questionnaire. We identified 66 differentially methylated CpG sites associated with smoking, located in 38 genes. In most of these CpG sites, we observed a trend among those quitting smoking to recover methylation levels typical of never smokers. A CpG site located in a novel smoking-associated gene (cg06394460 in LNX2) was hypomethylated in current smokers. Moreover, we validated two previously reported CpG sites (cg05886626 in THBS1, and cg24838345 in MTSS1) for their potential relation to atherosclerosis and cancer diseases, using several different approaches: CpG site methylation, gene expression, and plasma protein level determinations. Smoking was also associated with higher THBS1 gene expression but with lower levels of thrombospondin-1 in plasma. Finally, we identified differential methylation regions in 13 genes and in four non-coding RNAs. In summary, this study replicated previous findings and identified and validated a new CpG site located in LNX2 associated with smoking.     

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.     

Locus-specific DNA methylation in the placenta is associated with levels of pro-inflammatory proteins in cord blood and they are both independently affected by maternal smoking during pregnancy

We investigated the impact of maternal smoking during pregnancy on placental DNA methylation and how this may mediate the association between maternal smoking and pro-inflammatory proteins in cord blood. The study population consisted of 27 individuals exposed to maternal smoking throughout pregnancy, 32 individuals exposed during a proportion of the pregnancy, and 61 unexposed individuals. Methylation of 11 regions within 6 genes in placenta tissue was assessed by pyrosequencing. Levels of 7 pro-inflammatory proteins in cord blood were assessed by electrochemiluminescence. Differential methylation was observed in the CYP1A1 promoter and AHRR gene body regions between women who smoked throughout pregnancy and non-smokers on the fetal-side of the placenta and in the GFI1 promoter between women who quit smoking while pregnant and non-smokers on the maternal-side of the placenta. Maternal smoking resulted in elevated levels of IL-8 protein in cord blood, which was not mediated by DNA methylation of our candidate regions at either the maternal or the fetal side of the placenta. Placental DNA methylation was associated with levels of inflammatory proteins in cord blood. Our observations suggest that maternal smoking during pregnancy affects both placental DNA methylation and the neonate's immune response.     

Imputation of missing covariate values in epigenome-wide analysis of DNA methylation data

DNA methylation is a widely studied epigenetic mechanism and alterations in methylation patterns may be involved in the development of common diseases. Unlike inherited changes in genetic sequence, variation in site-specific methylation varies by tissue, developmental stage, and disease status, and may be impacted by aging and exposure to environmental factors, such as diet or smoking. These non-genetic factors are typically included in epigenome-wide association studies (EWAS) because they may be confounding factors to the association between methylation and disease. However, missing values in these variables can lead to reduced sample size and decrease the statistical power of EWAS. We propose a site selection and multiple imputation (MI) method to impute missing covariate values and to perform association tests in EWAS. Then, we compare this method to an alternative projection-based method. Through simulations, we show that the MI-based method is slightly conservative, but provides consistent estimates for effect size. We also illustrate these methods with data from the Atherosclerosis Risk in Communities (ARIC) study to carry out an EWAS between methylation levels and smoking status, in which missing cell type compositions and white blood cell counts are imputed.     

Differential DNA methylation in umbilical cord blood of infants exposed to mercury and arsenic in utero

Mercury and arsenic are known developmental toxicants. Prenatal exposures are associated with adverse childhood health outcomes that could be in part mediated by epigenetic alterations that may also contribute to altered immune profiles. In this study, we examined the association between prenatal mercury exposure on both DNA methylation and white blood cell composition of cord blood, and evaluated the interaction with prenatal arsenic exposure. A total of 138 mother-infant pairs with postpartum maternal toenail mercury, prenatal urinary arsenic concentrations, and newborn cord blood were assessed using the Illumina Infinium Methylation450 array. White blood cell composition was inferred from DNA methylation measurements. A doubling in toenail mercury concentration was associated with a 2.5% decrease (95% CI: 5.0%, 1.0%) in the estimated monocyte proportion. An increase of 3.5% (95% CI: 1.0, 7.0) in B-cell proportion was observed for females only. Among the top 100 CpGs associated with toenail mercury levels (ranked on P-value), there was a significant enrichment of loci located in North shore regions of CpG islands (P = 0.049), and the majority of these loci were hypermethylated (85%). Among the top 100 CpGs for the interaction between arsenic and mercury, there was a greater than expected proportion of loci located in CpG islands (P = 0.045) and in South shore regions (P = 0.009) and all of these loci were hypermethylated. This work supports the hypothesis that mercury may be contributing to epigenetic variability and immune cell proportion changes, and suggests that in utero exposure to mercury and arsenic, even at low levels, may interact to impact the epigenome.     

Differentiation of trophoblast lineage is associated with DNA methylation and demethylation.

Our previous study has shown that the placenta and kidney had different genomic methylation patterns regarding CpG island loci detected by restriction landmark genomic scanning (RLGS). To investigate whether differentiation involves changes in DNA methylation, we analyzed the rat Rcho-1 cell line, which retains trophoblast cell features and differentiates from stem cells into trophoblast giant cells in vitro. By RLGS, a total of 1,232 spots were identified in the Rcho-1 stem and differentiated giant cells. Four spots (0.3%) were detected only in giant cells, implying that the loci were originally methylated, but became demethylated during differentiation. Another four spots (0.3%) were detected only in stem cells, implying that these loci, originally unmethylated, became methylated during differentiation. DNAs from three loci that became methylated during differentiation were cloned and sequenced. All showed high homologies with expressed sequence tags (ESTs) or with genomic DNA of other species, suggesting that these loci are biologically important. Thus, the eight differentially methylated loci should be good tools to study epigenetic modification specific to differentiation of trophoblast giant cells.     

DNA methylation dynamics during the mammalian life cycle

DNA methylation is dynamically remodelled during the mammalian life cycle through distinct phases of reprogramming and de novo methylation. These events enable the acquisition of cellular potential followed by the maintenance of lineage-restricted cell identity, respectively, a process that defines the life cycle through successive generations. DNA methylation contributes to the epigenetic regulation of many key developmental processes including genomic imprinting, X-inactivation, genome stability and gene regulation. Emerging sequencing technologies have led to recent insights into the dynamic distribution of DNA methylation during development and the role of this epigenetic mark within distinct genomic contexts, such as at promoters, exons or imprinted control regions. Additionally, there is a better understanding of the mechanistic basis of DNA demethylation during epigenetic reprogramming in primordial germ cells and during pre-implantation development. Here, we discuss our current understanding of the developmental roles and dynamics of this key epigenetic system.