Reducing the risk of false discovery enabling identification of biologically significant genome-wide methylation status using the HumanMethylation450 array

BACKGROUND

The Illumina HumanMethylation450 BeadChip (HM450K) measures the DNA methylation of 485,512 CpGs in the human genome. The technology relies on hybridization of genomic fragments to probes on the chip. However, certain genomic factors may compromise the ability to measure methylation using the array such as single nucleotide polymorphisms (SNPs), small insertions and deletions (INDELs), repetitive DNA, and regions with reduced genomic complexity. Currently, there is no clear method or pipeline for determining which of the probes on the HM450K bead array should be retained for subsequent analysis in light of these issues.

RESULTS

We comprehensively assessed the effects of SNPs, INDELs, repeats and bisulfite induced reduced genomic complexity by comparing HM450K bead array results with whole genome bisulfite sequencing. We determined which CpG probes provided accurate or noisy signals. From this, we derived a set of high-quality probes that provide unadulterated measurements of DNA methylation.

CONCLUSIONS

Our method significantly reduces the risk of false discoveries when using the HM450K bead array, while maximising the power of the array to detect methylation status genome-wide. Additionally, we demonstrate the utility of our method through extraction of biologically relevant epigenetic changes in prostate cancer.

Epigenome-wide association of liver methylation patterns and complex metabolic traits in mice

Heritable epigenetic factors can contribute to complex disease etiology. Here we examine the contribution of DNA methylation to complex traits that are precursors to heart disease, diabetes, and osteoporosis. We profiled DNA methylation in the liver using bisulfite sequencing in 90 mouse inbred strains, genome-wide expression levels, proteomics, metabolomics, and 68 clinical traits and performed epigenome-wide association studies (EWAS). We found associations with numerous clinical traits including bone density, insulin resistance, expression, and protein and metabolite levels. A large proportion of associations were unique to EWAS and were not identified using GWAS. Methylation levels were regulated by genetics largely in cis, but we also found evidence of trans regulation, and we demonstrate that genetic variation in the methionine synthase reductase gene Mtrr affects methylation of hundreds of CpGs throughout the genome. Our results indicate that natural variation in methylation levels contributes to the etiology of complex clinical traits.

An epigenome-wide association study meta-analysis of educational attainment

The epigenome is associated with biological factors, such as disease status, and environmental factors, such as smoking, alcohol consumption and body mass index. Although there is a widespread perception that environmental influences on the epigenome are pervasive and profound, there has been little evidence to date in humans with respect to environmental factors that are biologically distal. Here we provide evidence on the associations between epigenetic modifications-in our case, CpG methylation-and educational attainment (EA), a biologically distal environmental factor that is arguably among the most important life-shaping experiences for individuals. Specifically, we report the results of an epigenome-wide association study meta-analysis of EA based on data from 27 cohort studies with a total of 10 767 individuals. We find nine CpG probes significantly associated with EA. However, robustness analyses show that all nine probes have previously been found to be associated with smoking. Only two associations remain when we perform a sensitivity analysis in the subset of never-smokers, and these two probes are known to be strongly associated with maternal smoking during pregnancy, and thus their association with EA could be due to correlation between EA and maternal smoking. Moreover, the effect sizes of the associations with EA are far smaller than the known associations with the biologically proximal environmental factors alcohol consumption, body mass index, smoking and maternal smoking during pregnancy. Follow-up analyses that combine the effects of many probes also point to small methylation associations with EA that are highly correlated with the combined effects of smoking. If our findings regarding EA can be generalized to other biologically distal environmental factors, then they cast doubt on the hypothesis that such factors have large effects on the epigenome.

Multi-OMICS analyses of frailty and chronic widespread musculoskeletal pain suggest involvement of shared neurological pathways

Chronic widespread musculoskeletal pain (CWP) and frailty are prevalent conditions in older people. We have shown previously that interindividual variation in frailty and CWP is genetically determined. We also reported an association of frailty and CWP caused by shared genetic and common environmental factors. The aim of this study was to use omic approaches to identify molecular genetic factors underlying the heritability of frailty and its genetic correlation with CWP. Frailty was quantified through the Rockwood Frailty Index (FI) as a proportion of deficits from 33 binary health deficit questions in 3626 female twins. Common widespread pain was assessed using a screening questionnaire. OMICS analysis included 305 metabolites and whole-genome (>2.5 × 10 SNPs) and epigenome (∼1 × 10 MeDIP-seq regions) assessments performed on fasting blood samples. Using family-based statistical analyses, including path analysis, we examined how FI scores were related to molecular genetic factors and to CWP, taking into account known risk factors such as fat mass and smoking. Frailty Index was significantly correlated with 51 metabolites after correction for multiple testing, with 20 metabolites having P-values between 2.1 × 10 and 4.0 × 10. Three metabolites (uridine, C-glycosyl tryptophan, and N-acetyl glycine) were statistically independent and thought to exert a direct effect on FI. Epiandrosterone sulphate, previously shown to be highly inversely associated with CWP, was found to exert an indirect influence on FI. Bioinformatics analysis of genome-wide association study and EWAS showed that FI and its covariation with CWP was through genomic regions involved in neurological pathways. Neurological pathway involvement accounts for the associated conditions of aging CWP and FI.

DNA methylation in human lipid metabolism and related diseases

PURPOSE OF REVIEW

It is becoming increasingly evident that epigenetic mechanisms, particularly DNA methylation, play a role in the regulation of blood lipid levels and lipid metabolism-linked phenotypes and diseases.

RECENT FINDINGS

Recent genome-wide methylation and candidate gene studies of blood lipids have highlighted several robustly replicated methylation markers across different ethnicities. Furthermore, many of these lipid-related CpG sites associated with blood lipids are also linked to lipid-related phenotypes and diseases. Integrating epigenome-wide association studies (EWAS) data with other layers of molecular data such as genetics or the transcriptome, accompanied by relevant statistical methods (e.g. Mendelian randomization), provides evidence for causal relationships. Recent data suggest that epigenetic changes can be consequences rather than causes of dyslipidemia. There is sparse information on many lipid classes and disorders of lipid metabolism, and also on the interplay of DNA methylation with other epigenetic layers such as histone modifications and regulatory RNAs.

SUMMARY

The current review provides a literature overview of epigenetic modifications in lipid metabolism and other lipid-related phenotypes and diseases focusing on EWAS of DNA methylation from January 2016 to September 2017. Recent studies strongly support the importance of epigenetic modifications, such as DNA methylation, in lipid metabolism and related diseases for relevant biological insights, reliable biomarkers, and even future therapeutics.

Placental DNA methylation changes associated with maternal prepregnancy BMI and gestational weight gain

BACKGROUND

Maternal obesity prior to or during pregnancy influences fetal growth, predisposing the offspring to increased risk for obesity across the life course. Placental epigenetic mechanisms may underlie these associations. We conducted an epigenome-wide association study to identify placental DNA methylation changes associated with maternal prepregnancy body mass index (BMI) and rate of gestational weight gain at first (GWG1), second (GWG2), and third trimester (GWG3).

METHOD

Participants of the NICHD Fetal Growth Studies with genome-wide placental DNA methylation (n = 301) and gene expression (n = 75) data were included. Multivariable-adjusted regression models were used to test the associations of 1 kg/m2 increase in prepregnancy BMI or 1 kg/week increase in GWG with DNA methylation levels. Genes harboring top differentially methylated CpGs (FDR P < 0.05) were evaluated for placental gene expression. We assessed whether DNA methylation sites known to be associated with BMI in child or adult tissues, were also associated with maternal prepregnancy BMI in placenta.

RESULTS

Prepregnancy BMI was associated with DNA methylation at cg14568196[EGFL7], cg15339142[VETZ], and cg02301019[AC092377.1] (FDR P < 0.05, P ranging from 1.4 × 10-10 to 1.7 × 10-9). GWG1 or GWG2 was associated with DNA methylation at cg17918270[MYT1L], cg20735365[DLX5], and cg17451688[SLC35F3] (FDR P < 0.05, P ranging from 6.4 × 10-10 to 1.2 × 10-8). Both prepregnancy BMI and DNA methylation at cg1456819 [EGFL7] were negatively correlated with EGFL7 expression in placenta (P < 0.05). Several CpGs previously implicated in obesity traits in children and adults were associated with prepregnancy BMI in placenta. Functional annotations revealed that EGFL7 is highly expressed in placenta and the differentially methylated CpG sites near EGFL7 and VEZT were cis-meQTL targets in blood.

CONCLUSIONS

We identified placental DNA methylation changes at novel loci associated with prepregnancy BMI and GWG. The overlap between CpGs associated with obesity traits in placenta and other tissues in children and adults suggests that epigenetic mechanisms in placenta may give insights to early origins of obesity.

Epigenome-wide association study for glyphosate induced transgenerational sperm DNA methylation and histone retention epigenetic biomarkers for disease

The herbicide glyphosate has been shown to promote the epigenetic transgenerational inheritance of pathology and disease in subsequent great-grand offspring (F3 generation). This generational toxicology suggests the impacts of environmental exposures need to assess subsequent generations. The current study was designed to identify epigenetic biomarkers for glyphosate-induced transgenerational diseases using an epigenome-wide association study (EWAS). Following transient glyphosate exposure of gestating female rats (F0 generation), during the developmental period of gonadal sex determination, the subsequent transgenerational F3 generation, with no direct exposure, were aged to 1 year and animals with specific pathologies identified. The pathologies investigated included prostate disease, kidney disease, obesity, and presence of multiple disease. The sperm were collected from the glyphosate lineage males with only an individual disease and used to identify specific differential DNA methylation regions (DMRs) and the differential histone retention sites (DHRs) associated with that pathology. Unique signatures of DMRs and DHRs for each pathology were identified for the specific diseases. Interestingly, at a lower statistical threshold overlapping sets of DMRs and DHRs were identified that were common for all the pathologies. This is one of the first observations that sperm histone retention can potentially act as a biomarker for specific diseases. The DMR and DHR associated genes were identified and correlated with known pathology specific-associated genes. Observations indicate transgenerational epigenetic biomarkers of disease pathology can be identified in the sperm that appear to assess disease susceptibility. These biomarkers suggest epigenetic diagnostics could potentially be used to facilitate preventative medicine.

Epigenome-wide meta-analysis of PTSD symptom severity in three military cohorts implicates DNA methylation changes in genes involved in immune system and oxidative stress

Epigenetic factors modify the effects of environmental factors on biological outcomes. Identification of epigenetic changes that associate with PTSD is therefore a crucial step in deciphering mechanisms of risk and resilience. In this study, our goal is to identify epigenetic signatures associated with PTSD symptom severity (PTSS) and changes in PTSS over time, using whole blood DNA methylation (DNAm) data (MethylationEPIC BeadChip) of military personnel prior to and following combat deployment. A total of 429 subjects (858 samples across 2 time points) from three male military cohorts were included in the analyses. We conducted two different meta-analyses to answer two different scientific questions: one to identify a DNAm profile of PTSS using a random effects model including both time points for each subject, and the other to identify a DNAm profile of change in PTSS conditioned on pre-deployment DNAm. Four CpGs near four genes (F2R, CNPY2, BAIAP2L1, and TBXAS1) and 88 differentially methylated regions (DMRs) were associated with PTSS. Change in PTSS after deployment was associated with 15 DMRs, of those 2 DMRs near OTUD5 and ELF4 were also associated with PTSS. Notably, three PTSS-associated CpGs near F2R, BAIAP2L1 and TBXAS1 also showed nominal evidence of association with change in PTSS. This study, which identifies PTSD-associated changes in genes involved in oxidative stress and immune system, provides novel evidence that epigenetic differences are associated with PTSS.

Novel epigenetic biomarkers for hematopoietic cancer found in twins

BACKGROUND AND PURPOSE

This article aims to identify epigenetic markers and detect early development of hematopoietic malignancies through an epigenome wide association study of DNA methylation data.

MATERIALS AND METHODS

This register-based study includes 1,085 Danish twins with 31 hematopoietic malignancies and methylation levels from 450,154 5'-C-phospate-G-3' (CpG) sites. Associations between methylation levels and incidence of hematopoietic malignancy is studied through time-to-event regression. The matched case-cotwin design, where one twin has a malignancy and the cotwin does not, is applied to enhance control for unmeasured shared confounding and false discoveries. Predictive performance is validated in the independent Older Finnish Twin Cohort.

RESULTS AND INTERPRETATION

We identified 67 epigenetic markers for hematopoietic malignancies of which 12 are linked to genes associated with hematologic malignancies. For some markers, we discovered a 2-3-fold relative risk difference for high versus low methylation. The identification of these 67 sites enabled the formation of a predictor demonstrating a cross-validated time-varying area under the curve (AUC) of 92% 3 years after individual blood sampling and persistent performance above 70% up to 6 years after blood sampling. This predictive performance was to a large extent recovered in the validation sample showing an overall Harrell's C of 73%. In conclusion, from a large population representative twin study on hematopoietic cancers, novel epigenetic markers were identified that may prove useful for early diagnosis.

Newborn adiposity is associated with cord blood DNA methylation at IGF1R and KLF7

OBJECTIVE

This study aimed to identify whether cord blood DNA methylation at specific loci is associated with neonatal adiposity, a key risk factor for childhood obesity.

METHODS

An epigenome-wide association study was conducted using the Hyperglycemia and Adverse Pregnancy Outcome (HAPO) Study as a discovery sample. Linear regression models adjusted for maternal and offspring covariates and cell counts were used to analyze associations between neonatal adiposity as measured by sum of three skinfold thicknesses and cord blood DNA methylation. Assays were performed with Illumina EPIC arrays (791,359 CpG sites after quality control). Replication was performed in an independent cohort, Genetics of Glucose regulation in Gestation and Growth (Gen3G).

RESULTS

In 2740 HAPO samples, significant associations were identified at 89 CpG sites after accounting for multiple testing (Bonferroni-adjusted p < 0.05). Replication analyses conducted in 139 Gen3G participants confirmed associations for seven CpG sites. These included IGF1R, which encodes a transmembrane receptor involved in cell growth and survival that binds insulin-like growth factor I and insulin, and KLF7, which encodes a regulator of cell proliferation and inhibitor of adipogenesis; both are key regulators of growth during fetal life.

CONCLUSIONS

These findings support epigenetic mechanisms in the developmental origins of neonatal adiposity and as potential biomarkers of metabolic disease risk.