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.     

Genetic variation in DNMT3B and increased global DNA methylation is associated with suicide attempts in psychiatric patients

Recently, a significant epigenetic component in the pathology of suicide has been realized. Here we investigate candidate functional SNPs in epigenetic‐regulatory genes, DNMT1 and DNMT3B, for association with suicide attempt (SA) among patients with co‐existing psychiatric illness. In addition, global DNA methylation levels [5‐methyl cytosine (5‐mC%)] between SA and psychiatric controls were quantified using the Methylflash Methylated DNA Quantification Kit. DNA was obtained from blood of 79 suicide attempters and 80 non‐attempters, assessed for DSM‐IV Axis I disorders. Functional SNPs were selected for each gene (DNMT1; n = 7, DNMT3B; n = 10), and genotyped. A SNP (rs2424932) residing in the 3′ UTR of the DNMT3B gene was associated with SA compared with a non‐attempter control group (P = 0.001; Chi‐squared test, Bonferroni adjusted P value = 0.02). Moreover, haplotype analysis identified a DNMT3B haplotype which differed between cases and controls, however this association did not hold after Bonferroni correction (P = 0.01, Bonferroni adjusted P value = 0.56). Global methylation analysis showed that psychiatric patients with a history of SA had significantly higher levels of global DNA methylation compared with controls (P = 0.018, Student's t‐test). In conclusion, this is the first report investigating polymorphisms in DNMT genes and global DNA methylation quantification in SA risk. Preliminary findings suggest that allelic variability in DNMT3B may be relevant to the underlying diathesis for suicidal acts and our findings support the hypothesis that aberrant DNA methylation profiles may contribute to the biology of suicidal acts. Thus, analysis of global DNA hypermethylation in blood may represent a biomarker for increased SA risk in psychiatric patients.     

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.     

DNA methylation in embryonic stem cells

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

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.     

Genome‐wide variation in DNA methylation is associated with stress resilience and plumage brightness in a wild bird

Individuals often differ in their ability to cope with challenging environmental and social conditions. Evidence from model systems suggests that patterns of DNA methylation are associated with variation in coping ability. These associations could arise directly if methylation plays a role in controlling the physiological response to stressors by, among other things, regulating the release of glucocorticoids in response to challenges. Alternatively, the association could arise indirectly if methylation and resilience have a common cause, such as early‐life conditions. In either case, methylation might act as a biomarker for coping ability. At present, however, relatively little is known about whether variation in methylation is associated with organismal performance and resilience under natural conditions. We studied genome‐wide patterns of DNA methylation in free‐living female tree swallows (Tachycineta bicolor) using methylated DNA immunoprecipitation (MeDIP) and a tree swallow genome that was assembled for this study. We identified areas of the genome that were differentially methylated with respect to social signal expression (breast brightness) and physiological traits (ability to terminate the glucocorticoid stress response through negative feedback). We also asked whether methylation predicted resilience to a subsequent experimentally imposed challenge. Individuals with brighter breast plumage and higher stress resilience had lower methylation at differentially methylated regions across the genome. Thus, widespread differences in methylation predicted both social signal expression and the response to future challenges under natural conditions. These results have implications for predicting individual differences in resilience, and for understanding the mechanistic basis of resilience and its environmental and social mediators.     

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.     

Expression of various genes is controlled by DNA methylation during mammalian development

Despite thousands of articles about 5-methylcytosine (m(5)C) residues in vertebrate DNA, there is still controversy concerning the role of genomic m(5)C in normal vertebrate development. Inverse correlations between expression and methylation are seen for many gene regulatory regions [Heard et al., 1997; Attwood et al., 2002; Plass and Soloway, 2002] although much vertebrate DNA methylation is in repeated sequences [Ehrlich et al., 1982]. At the heart of this debate is whether vertebrate DNA methylation has mainly a protective role in limiting expression of foreign DNA elements and endogenous transposons [Walsh and Bestor, 1999] or also is important in the regulation of the expression of diverse vertebrate genes involved in differentiation [Attwood et al., 2002]. Enough thorough studies have now been reported to show that many tissue- or development-specific changes in methylation at vertebrate promoters, enhancers, or insulators regulate expression and are not simply inconsequential byproducts of expression differences. One line of evidence comes from mutants with inherited alterations in genes encoding DNA methyltransferases and from rodents or humans with somatically acquired changes in DNA methylation that illustrate the disease-producing effects of abnormal methylation. Another type of evidence derives from studies of in vivo correlations between tissue-specific changes in DNA methylation and gene expression coupled with experiments demonstrating cause-and-effect associations between DNA hyper- or hypomethylation and gene expression. In this review, I summarize some of the strong evidence from both types of studies. Taken together, these studies demonstrate that DNA methylation in mammals modulates expression of many genes during development, causing major changes in or important fine-tuning of expression. Also, I discuss previously established and newly hypothesized mechanisms for this epigenetic control.     

DNA methylation in imprinted genes IGF2 and GNASXL is associated with prenatal maternal stress

Epigenetic regulation of imprinted genes during embryonic development is influenced by the prenatal environment. Our aim was to examine the effect of maternal emotional stress and cortisol levels during pregnancy on methylation of imprinted genes, insulin‐like growth factor 2 (IGF2) and guanine nucleotide‐binding protein, alpha stimulating extra‐large (GNASXL), using umbilical cord blood DNA. Maternal depressed mood (Edinburgh Depression Scale; EDS), pregnancy‐related anxiety questionnaire (PRAQ) and cortisol day profiles were assessed throughout pregnancy. At birth, a cord blood sample (n = 80) was taken to study DNA methylation of IGF2 DMR0 (differentially methylated region), IGF2 anti‐sense (IGF2AS) and GNASXL using Sequenom Epi TYPER. Linear mixed models were used to examine the relationship between DNA methylation and maternal stress, while correcting for confounders. We also studied the association of DNA methylation with the child ponderal index at birth. We found a cytosine–guanine dinucleotide (CpG)‐specific association of PRAQ subscales with IGF2 DMR0 (CpG5, P < 0.0001) and GNASXL (CpG11, P = 0.0003), while IGF2AS was associated with maternal EDS scores (CpG33, P = 0.0003) and cortisol levels (CpG33, P = 0.0006; CpG37‐38, P = 0.0005). However, there was no association of methylation with ponderal index at birth. In conclusion, maternal stress during pregnancy, as defined by cortisol measurements, EDS and PRAQ scores, is associated with DNA methylation of imprinted genes IGF2 and GNASXL. Our results provide further evidence that prenatal adversity can influence imprinted gene methylation, although future studies are needed to unravel the exact mechanisms.     

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.