Recommendations for the design and analysis of epigenome-wide association studies

Advances in epigenome-wide association studies for common diseases

Epigenome-wide association studies (EWASs) provide a systematic approach to uncovering epigenetic variants underlying common diseases. Discoveries have shed light on novel molecular mechanisms of disease and enabled the application of epigenetic variants as biomarkers. Here, we highlight the recent advances in this emerging line of research and discuss key challenges for current and future studies.

Dietary fat quality impacts genome-wide DNA methylation patterns in a cross-sectional study of Greek preadolescents

The type and the amount of dietary fat have a significant influence on the metabolic pathways involved in the development of obesity, metabolic syndrome, diabetes type 2 and cardiovascular diseases. However, it is unknown to what extent this modulation is achieved through DNA methylation. We assessed the effects of cholesterol intake, the proportion of energy intake derived from fat, the ratio of polyunsaturated fatty acids (PUFA) to saturated fatty acids (SFA), the ratio of monounsaturated fatty acids (MUFA) to SFA, and the ratio of MUFA+PUFA to SFA on genome-wide DNA methylation patterns in normal-weight and obese children. We determined the genome-wide methylation profile in the blood of 69 Greek preadolescents (∼10 years old) as well as their dietary intake for two consecutive weekdays and one weekend day. The methylation levels of one CpG island shore and four sites were significantly correlated with total fat intake. The methylation levels of 2 islands, 11 island shores and 16 sites were significantly correlated with PUFA/SFA; of 9 islands, 26 island shores and 158 sites with MUFA/SFA; and of 10 islands, 40 island shores and 130 sites with (MUFA+PUFA)/SFA. We found significant gene enrichment in 34 pathways for PUFA/SFA, including the leptin pathway, and a significant enrichment in 5 pathways for (MUFA+PUFA)/SFA. Our results suggest that specific changes in DNA methylation may have an important role in the mechanisms involved in the physiological responses to different types of dietary fat.

Variation of DNA methylation in candidate age-related targets on the mitochondrial-telomere axis in cord blood and placenta

BACKGROUND

Epigenetics is tissue-specific and potentially even cell-specific, but little information is available from human reproductive studies about the concordance of DNA methylation patterns in cord blood and placenta, as well as within-placenta variations. We evaluated methylation levels at promoter regions of candidate genes in biological ageing pathways (SIRT1, TP53, PPARG, PPARGC1A, and TFAM), a subtelomeric region (D4Z4) and the mitochondrial genome (MT-RNR1, D-loop).

METHODS

Ninety individuals were randomly chosen from the ENVIRONAGE birth cohort to evaluate methylation concordance between cord blood and placenta using highly quantitative bisulfite-PCR pyrosequencing. In a subset of nineteen individuals, a more extensive sampling scheme was performed to examine within-placenta variation.

RESULTS

The DNA methylation levels of the subtelomeric region and mitochondrial genome showed concordance between cord blood and placenta with correlation coefficients ranging from r = 0.31 to 0.43, p ≤ 0.005, and also between the maternal and foetal sides of placental tissue (r = 0.53 to 0.72, p ≤ 0.05). For the majority of targets, an agreement in methylation levels between four foetal biopsies was found (with intra-class correlation coefficients ranging from 0.16 to 0.72), indicating small within-placenta variation.

CONCLUSIONS

The methylation levels of the subtelomeric region (D4Z4) and mitochondrial genome (MT-RNR1, D-loop) showed concordance between cord blood and placenta, suggesting a common epigenetic signature of these targets between tissues. Concordance was lacking between the other genes that were studied. In placental tissue, methylation patterns of most targets on the mitochondrial-telomere axis were not strongly influenced by sample location.

Grasping nettles: cellular heterogeneity and other confounders in epigenome-wide association studies

Platform technologies for measurement of CpG methylation at multiple loci across the genome have made ambitious epigenome-wide association studies affordable and practicable. In contrast to genetic studies, which estimate the effects of structural changes in DNA, and transcriptomic studies, which measure genomic outputs, epigenetic studies can access states of regulation of genome function in particular cells and in response to specific stimuli. Although many factors complicate the interpretation of epigenetic variation in human disease, cell-specific methylation patterns and the cellular heterogeneity present in peripheral blood and tissue biopsies are anticipated to cause the most problems. In this review, we suggest that the difficulties may be exaggerated and we explore how cellular heterogeneity may be embraced with appropriate study designs and analytical tools. We further suggest that systematic mapping of the loci influenced by age, sex and genetic polymorphisms will bring important biological insights as well as improved control of epigenome-wide association studies.

An integrated platform for bovine DNA methylome analysis suitable for small samples

Background

Oocytes and early embryos contain minute amounts of DNA, RNA and proteins, making the study of early mammalian development highly challenging. The study of the embryo epigenome, in particular the DNA methylome, has been made accessible thanks to the possibility of amplifying specific sequences according to their initial methylation status. This paper describes a novel platform dedicated to the genome-wide study of bovine DNA methylation, including a complete pipeline for data analysis and visualization. The platform allows processing and integrating of DNA methylome and transcriptome data from the same sample. Procedures were optimized for genome-wide analysis of 10 ng of DNA (10 bovine blastocysts). Bovine sperm and blastocysts were compared as a test of platform capability.

Results

The hypermethylation of bovine sperm DNA compared to the embryo genome was confirmed. Differentially methylated regions were distributed across various classes of bovine sperm genomic feature including primarily promoter, intronic and exonic regions, non-CpG-island regions (shore, shelf and open-sea) and CpG islands with low-to-intermediate CpG density. The blastocyst genome bore more methylation marks than sperm DNA only in CpG islands with high CpG density. Long-terminal-repeat retrotransposons (LTR), LINE and SINE were more methylated in sperm DNA, as were low-complexity repetitive elements in blastocysts.

Conclusions

This is the first early embryo compatible genome-wide epigenetics platform for bovine. Such platforms should improve the study of the potential epigenetic risks of assisted reproductive technologies (ART), the establishment sequence of embryonic cell lines and potential deviations in both gene expression and DNA methylation capable of having long-term impact.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-451) contains supplementary material, which is available to authorized users.

Infinium monkeys: Infinium 450K array for the Cynomolgus macaque (Macaca fascicularis)

The Infinium Human Methylation450 BeadChip Array (Infinium 450K) is a robust and cost-efficient survey of genome-wide DNA methylation patterns. Macaca fascicularis (Cynomolgus macaque) is an important disease model; however, its genome sequence is only recently published, and few tools exist to interrogate the molecular state of Cynomolgus macaque tissues. Although the Infinium 450K is a hybridization array designed to the human genome, the relative conservation between the macaque and human genomes makes its use in macaques feasible. Here, we used the Infinium 450K array to assay DNA methylation in 11 macaque muscle biopsies. We showed that probe hybridization efficiency was related to the degree of sequence identity between the human probes and the macaque genome sequence. Approximately 61% of the Human Infinium 450K probes could be reliably mapped to the Cynomolgus macaque genome and contain a CpG site of interest. We also compared the Infinium 450K data to reduced representation bisulfite sequencing data generated on the same samples and found a high level of concordance between the two independent methodologies, which can be further improved by filtering for probe sequence identity and mismatch location. We conclude that the Infinium 450K array can be used to measure the DNA methylome of Cynomolgus macaque tissues using the provided filters. We also provide a pipeline for validation of the array in other species using a simple BLAST-based sequence identify filter.

Epigenome-wide association study reveals longitudinally stable DNA methylation differences in CD4+ T cells from children with IgE-mediated food allergy

Food allergy is mediated by a combination of genetic and environmental risk factors, potentially mediated by epigenetic mechanisms. CD4+ T-cells are key drivers of the allergic response, and may therefore harbor epigenetic variation in association with the disease phenotype. Here we retrospectively examined genome-wide DNA methylation profiles (~450,000 CpGs) from CD4+ T-cells on a birth cohort of 12 children with IgE-mediated food allergy diagnosed at 12-months, and 12 non-allergic controls. DNA samples were available at two time points, birth and 12-months.

CASE

control comparisons of CD4+ methylation profiles identified 179 differentially methylated probes (DMP) at 12-months and 136 DMP at birth (FDR-adjusted P value<0.05, delta β>0.1). Approximately 30% of DMPs were coincident with previously annotated SNPs. A total of 92 [corrected] allergy-associated non-SNP DMPs were present at birth when individuals were initially disease-free, potentially implicating these loci in the causal pathway. Pathway analysis of differentially methylated genes identified several MAP kinase signaling molecules. Mass spectrometry was used to validate 15 CpG sites at 3 candidate genes. Combined analysis of differential methylation with gene expression profiles revealed gene expression differences at some but not all allergy associated differentially methylated genes. Thus, dysregulation of DNA methylation at MAPK signaling-associated genes during early CD4+ T-cell development may contribute to suboptimal T-lymphocyte responses in early childhood associated with the development of food allergy.

The effect of genotype and in utero environment on interindividual variation in neonate DNA methylomes

Integrating the genotype with epigenetic marks holds the promise of better understanding the biology that underlies the complex interactions of inherited and environmental components that define the developmental origins of a range of disorders. The quality of the in utero environment significantly influences health over the lifecourse. Epigenetics, and in particular DNA methylation marks, have been postulated as a mechanism for the enduring effects of the prenatal environment. Accordingly, neonate methylomes contain molecular memory of the individual in utero experience. However, interindividual variation in methylation can also be a consequence of DNA sequence polymorphisms that result in methylation quantitative trait loci (methQTLs) and, potentially, the interaction between fixed genetic variation and environmental influences. We surveyed the genotypes and DNA methylomes of 237 neonates and found 1423 punctuate regions of the methylome that were highly variable across individuals, termed variably methylated regions (VMRs), against a backdrop of homogeneity. MethQTLs were readily detected in neonatal methylomes, and genotype alone best explained ∼25% of the VMRs. We found that the best explanation for 75% of VMRs was the interaction of genotype with different in utero environments, including maternal smoking, maternal depression, maternal BMI, infant birth weight, gestational age, and birth order. Our study sheds new light on the complex relationship between biological inheritance as represented by genotype and individual prenatal experience and suggests the importance of considering both fixed genetic variation and environmental factors in interpreting epigenetic variation.

Possible molecular mechanisms linking air pollution and asthma in children

BACKGROUND

Air pollution has many effects on the health of both adults and children, but children's vulnerability is unique. The aim of this review is to discuss the possible molecular mechanisms linking air pollution and asthma in children, also taking into account their genetic and epigenetic characteristics.

RESULTS

Air pollutants appear able to induce airway inflammation and increase asthma morbidity in children. A better definition of mechanisms related to pollution-induced airway inflammation in asthmatic children is needed in order to find new clinical and therapeutic strategies for preventing the exacerbation of asthma. Moreover, reducing pollution-induced oxidative stress and consequent lung injury could decrease children's susceptibility to air pollution. This would be extremely useful not only for the asthmatic children who seem to have a genetic susceptibility to oxidative stress, but also for the healthy population. In addition, epigenetics seems to have a role in the lung damage induced by air pollution. Finally, a number of epidemiological studies have demonstrated that exposure to common air pollutants plays a role in the susceptibility to, and severity of respiratory infections.

CONCLUSIONS

Air pollution has many negative effects on pediatric health and it is recognised as a serious health hazard. There seems to be an association of air pollution with an increased risk of asthma exacerbations and acute respiratory infections. However, further studies are needed in order to clarify the specific mechanism of action of different air pollutants, identify genetic polymorphisms that modify airway responses to pollution, and investigate the effectiveness of new preventive and/or therapeutic approaches for subjects with low antioxidant enzyme levels. Moreover, as that epigenetic changes are inheritable during cell division and may be transmitted to subsequent generations, it is very important to clarify the role of epigenetics in the relationship between air pollution and lung disease in asthmatic and healthy children.

Pyrosequencing Evaluation of Widely Available Bisulfite Conversion Methods: Considerations for Application

Introduction

Bisulfite treatment of DNA introduces methylation-dependent sequence changes through selective chemical conversion of nonmethylated cytosine to uracil and serves as pretreatment step for the majority of DNA methylation analysis methods.

Methods

We have evaluated the conversion performance of five of the most commonly used bisulfite treatment kits [MethylDetector (Active Motif), Epitect+ (Qiagen), Zymo Methylation, Zymo Gold and Zymo Lightning (all from Zymo Research)] by pyrosequencing four different regions with variable methylation levels, including: a repetitive element (ALUSX), a gene with low levels of methylation (IL6ST), an imprinted gene expected to be approximately 50% methylated (IGF2), and a fully methylated gene (ST3GAL2). In addition, we have studied the influence of duration (3 vs. 16 h) and type (fixed temperature vs. cycling program) of incubation protocol on the conversion efficiency of each evaluated kit.

Results

All kits produced similar conversion rates of ALUSX, IGF2 and ST3GAL2, while the conversion of the low methylated IL6ST gene was variable between kits. The Zymo kits were highly consistent in their performance even when different protocols of incubation were applied, generating full conversion at the low methylated gene IL6; this was not true for the MethylDetector and Epitect+ kits. However, long-cycling incubation could produce similar conversion rates for the same locus in combination with Active Motif and Qiagen kits.

Conclusions

The selection of a long-cycling protocol during conversion permits standardization of protocols, improving the reproducibility of methylation estimates across laboratories for gene-specific, genome-wide and bisulfite-based sequencing analyses.

Assessment of RainDrop BS-seq as a method for large-scale, targeted bisulfite sequencing

We present a systematic assessment of RainDrop BS-seq, a novel method for large-scale, targeted bisulfite sequencing using microdroplet-based PCR amplification coupled with next-generation sequencing. We compared DNA methylation levels at 498 target loci (1001 PCR amplicons) in human whole blood, osteosarcoma cells and an archived tumor tissue sample. We assessed the ability of RainDrop BS-seq to accurately measure DNA methylation over a range of DNA quantities (from 10 to 1500 ng), both with and without whole-genome amplification (WGA) following bisulfite conversion. DNA methylation profiles generated using at least 100 ng correlated well (median R = 0.92) with those generated on Illumina Infinium HumanMethylation450 BeadChips, currently the platform of choice for epigenome-wide association studies (EWAS). WGA allowed for testing of samples with a starting DNA amount of 10 and 50 ng, although a reduced correlation was observed (median R = 0.79). We conclude that RainDrop BS-seq is suitable for measuring DNA methylation levels using nanogram quantities of DNA, and can be used to study candidate epigenetic biomarker loci in an accurate and high-throughput manner, paving the way for its application to routine clinical diagnostics.

Array-based DNA methylation profiling in male infertility reveals allele-specific DNA methylation in PIWIL1 and PIWIL2

OBJECTIVE

To identify CpG sites differentially methylated in peripheral blood of men with idiopathic infertility due to impaired spermatogenesis as compared with fertile controls.

DESIGN

DNA methylation profiling on peripheral blood samples using the HumanMethylation450 BeadChip (Illumina) in patients and controls, single-nucleotide polymorphism (SNP) typing by Sanger sequencing.

SETTING

University institute in cooperation with genetic and infertility clinics.

PATIENT(S)

30 infertile men with normal CFTR and AZF tests and karyotype, and 10 fertile male controls.

INTERVENTION(S)

None.

MAIN OUTCOME MEASURE(S)

DNA methylation levels at CpG sites.

RESULT(S)

We identified 471 CpGs (287 genes) as differentially methylated between patients and controls. These were significantly enriched for the gene ontology functions MHC class II receptor activity and piwi-interacting (piRNA) binding. The latter was associated with two methylation-sensitive SNPs in the genes PIWIL1 and PIWIL2, respectively, which showed significant allele distribution skewing in the infertile cohort. We found that 445 (94.5%) of 471 differentially methylated CpGs were associated with SNPs, but 26 (15 genes) were not genomically templated, including the ENO1, MTA2, BRSK2, and LBX2 genes previously associated with fertility and spermatogenesis.

CONCLUSION(S)

Our study identifies surrogate DNA methylation markers for idiopathic infertility in peripheral blood and suggests that allele-specific DNA methylation differences at regulatory sites of genes involved in piRNA regulation are associated with disturbed spermatogenesis.

Improved reporting of DNA methylation data derived from studies of the human placenta

Epigenetic variation is increasingly hypothesized as a mechanism underlying the effect of the in utero environment on long-term postnatal health; however, there is currently little clear data to support this in humans. A number of biological and technical factors provide challenges for the design of clinical epigenetic studies: from the type of cells or tissues that are available to the large range of predicted confounders that may impact findings. The human placenta, in addition to other neonatal tissues and whole blood, is commonly sampled for the study of epigenetic modifications. However there is little conformity for the most appropriate methods for study design, data analysis, and importantly, data interpretation. Here we present general recommendations for the reporting of DNA methylation in biological samples, with specific focus on the placenta. We outline key guidelines for: (1) placental sampling, (2) data analysis and presentation, and (3) interpretation of DNA methylation data. We emphasize the need to consider methodological noise, increase statistical power and to ensure appropriate adjustment for biological covariates. Finally, we highlight that epigenetic changes may be non-pathological and not necessarily translate into disease-associated changes. Improved reporting of DNA methylation data will be critical to identify epigenetic-based effects and to better understand the full phenotypic impact of these widely-reported epigenomic changes.