Genome- and epigenome-wide association studies identify susceptibility of CpG sites and regions for metabolic syndrome in a Korean population

BACKGROUND

While multiple studies have investigated the relationship between metabolic syndrome (MetS) and its related traits (fasting glucose, triglyceride, HDL cholesterol, blood pressure, waist circumference) and DNA methylation, our understanding of the epigenetic mechanisms in MetS remains limited. Therefore, we performed an epigenome-wide meta-analysis of blood DNA methylation to identify differentially methylated probes (DMPs) and differentially methylated regions (DMRs) associated with MetS and its components using two independent cohorts comprising a total of 2,334 participants. We also investigated the specific genetic effects on DNA methylation, identified methylation quantitative trait loci (meQTLs) through genome-wide association studies and further utilized Mendelian randomization (MR) to assess how these meQTLs subsequently influence MetS status.

RESULTS

We identified 40 DMPs and 27 DMRs that are significantly associated with MetS. In addition, we identified many novel DMPs and DMRs underlying inflammatory and steroid hormonal processes. The most significant associations were observed in 3 DMPs (cg19693031, cg26974062, cg02988288) and a DMR (chr1:145440444-145441553) at the TXNIP, which are involved in lipid metabolism. These CpG sites were identified as coregulators of DNA methylation in MetS, TG and FAG levels. We identified a total of 144 cis-meQTLs, out of which only 13 were found to be associated with DMPs for MetS. Among these, we confirmed the identified causal mediators of genetic effects at CpG sites cg01881899 at ABCG1 and cg00021659 at the TANK genes for MetS.

CONCLUSIONS

This study observed whether specific CpGs and methylated regions act independently or are influenced by genetic effects for MetS and its components in the Korean population. These associations between the identified DNA methylation and MetS, along with its individual components, may serve as promising targets for the development of preventive interventions for MetS.

Epigenome-wide association study on ambient PM2.5 exposure in Han Chinese, the NSPT study

Ambient PM2.5 exposure has been recognized as a major health risk and related to aging, cardiovascular, respiratory and neurologic diseases, and cancer. However, underlying mechanism of epigenetic alteration and regulated pathways still remained unclear. The study on methylome effect of PM2.5 exposure was quite limited in Chinese population, and cohort-based study was absent. The study included blood-derived DNA methylation for 3365 Chinese participants from the NSPT cohort. We estimated individual PM2.5 exposure level of short-medium-, medium- and long-term, based on a validated prediction model. We preformed epigenome-wide association studies to estimate the links between PM2.5 exposure and DNA methylation change, as well as stratification and sensitive analysis to examined the robustness of the association models. A systematic review was conducted to obtain the previously published CpGs and examined for replication. We also conducted comparison on the DNA methylation variation corresponding to different time windows. We further conducted gene function analysis and pathway enrichment analysis to reveal related biological response. We identified a total of 177 CpGs and 107 DMRs associated with short-medium-term PM2.5 exposure, at a strict genome-wide significance (P < 5 × 10-8). The effect sizes on most CpGs tended to cease with the exposure of extended time scale. Associated markers and aligned genes were related to aging, immunity, inflammation and carcinogenesis. Enriched pathways were mostly involved in cell cycle and cell division, signal transduction, inflammatory pathway. Our study is the first EWAS on PM2.5 exposure conducted in large-scale Han Chinese cohort and identified associated DNA methylation change on CpGs and regions, as well as related gene functions and pathways.

Alteration of DNA Methylation and Epigenetic Scores Associated With Features of Schizophrenia and Common Variant Genetic Risk

BACKGROUND

Unpacking molecular perturbations associated with features of schizophrenia is a critical step toward understanding phenotypic heterogeneity in this disorder. Recent epigenome-wide association studies have uncovered pervasive dysregulation of DNA methylation in schizophrenia; however, clinical features of the disorder that account for a large proportion of phenotypic variability are relatively underexplored.

METHODS

We comprehensively analyzed patterns of DNA methylation in a cohort of 381 individuals with schizophrenia from the deeply phenotyped Australian Schizophrenia Research Bank. Epigenetic changes were investigated in association with cognitive status, age of onset, treatment resistance, Global Assessment of Functioning scores, and common variant polygenic risk scores for schizophrenia. We subsequently explored alterations within genes previously associated with psychiatric illness, phenome-wide epigenetic covariance, and epigenetic scores.

RESULTS

Epigenome-wide association studies of the 5 primary traits identified 662 suggestively significant (p < 6.72 × 10-5) differentially methylated probes, with a further 432 revealed after controlling for schizophrenia polygenic risk on the remaining 4 traits. Interestingly, we uncovered many probes within genes associated with a variety of psychiatric conditions as well as significant epigenetic covariance with phenotypes and exposures including acute myocardial infarction, C-reactive protein, and lung cancer. Epigenetic scores for treatment-resistant schizophrenia strikingly exhibited association with clozapine administration, while epigenetic proxies of plasma protein expression, such as CCL17, MMP10, and PRG2, were associated with several features of schizophrenia.

CONCLUSIONS

Our findings collectively provide novel evidence suggesting that several features of schizophrenia are associated with alteration of DNA methylation, which may contribute to interindividual phenotypic variation in affected individuals.

Does epigenetic markers of HLA gene show association with coronary artery disease in Indian subjects?

BACKGROUND

DNA methylation, one of the most stable forms of epigenetic modification is associated with the development and progression of coronary artery disease (CAD). Our previously reported study on epigenome-wide microarray analysis showed significantly methylated CpG sites. Top 5 significant CpGs from HLA gene were selected and analysed by Pyrosequencing (PSQ) to determine their association with severity of CAD.

METHODS

Blood samples of 50-age matched angiographically CAD positive male cases with 50 angiographically CAD negative male controls were subjected to lipid profile estimation and PSQ for methylation level analysis. Findings and subgroup analysis were evaluated by Mann-Whitney U; Kruskal-Wallis' rank test and two-way ANOVA by MedCalc (v19.6).

RESULTS

Methylation levels in HLA-DQA1 for cg10217052 was 78.5 (37-85) and 76.5 (24-84); cg09411910 was 81 (72.0 to 93.0) and 81.5 (50.0 to 89.0) in cases and controls respectively. Levels in HLA-DQB1-cg03344051, were 28.88 + 9.41 for cases and 30.36 + 9.37 in controls. For HLA-DRB1-cg07889003, levels in cases and controls were 15.5 (5.00-39.00) and 10.5 (5.00-29.0); while in cg08269402 were 52 (16-65) and 42.5 (17-61) respectively. No association was observed between methylation levels and lipid profile. cg03344051, cg07889003 and cg08269402 were significantly differentiated in double or triple vessel disease (DVD or TVD) as compared to single vessel disease (SVD) suggesting an increase in the extent of methylation with the increase in CAD severity.

CONCLUSION

The present study shows significant increase in the extent of methylation in 3 CpG sites in DVD/TVD cases as compared to SVD cases. Additionally, a novel site, cg07889003 identified in our discovery phase has shown association with the severity of CAD.

Tumor microenvironment deconvolution identifies cell-type-independent aberrant DNA methylation and gene expression in prostate cancer

BACKGROUND

Among men, prostate cancer (PCa) is the second most common cancer and the second leading cause of cancer death. Etiologic factors associated with both prostate carcinogenesis and somatic alterations in tumors are incompletely understood. While genetic variants associated with PCa have been identified, epigenetic alterations in PCa are relatively understudied. To date, DNA methylation (DNAm) and gene expression (GE) in PCa have been investigated; however, these studies did not correct for cell-type proportions of the tumor microenvironment (TME), which could confound results.

METHODS

The data (GSE183040) consisted of DNAm and GE data from both tumor and adjacent non-tumor prostate tissue of 56 patients who underwent radical prostatectomies prior to any treatment. This study builds upon previous studies that examined methylation patterns and GE in PCa patients by using a novel tumor deconvolution approach to identify and correct for cell-type proportions of the TME in its epigenome-wide association study (EWAS) and differential expression analysis (DEA).

RESULTS

The inclusion of cell-type proportions in EWASs and DEAs reduced the scope of significant alterations associated with PCa. We identified 2,093 significantly differentially methylated CpGs (DMC), and 51 genes associated with PCa, including PCA3, SPINK1, and AMACR.

CONCLUSIONS

This work illustrates the importance of correcting for cell types of the TME when performing EWASs and DEAs on PCa samples, and establishes a more confounding-adverse methodology. We identified a more tumor-cell-specific set of altered genes and epigenetic marks that can be further investigated as potential biomarkers of disease or potential therapeutic targets.

Alterations in DNA methylation associate with reduced migraine and headache days after medication withdrawal treatment in chronic migraine patients: a longitudinal study

BACKGROUND

Chronic migraine, a highly disabling migraine subtype, affects nearly 2% of the general population. Understanding migraine chronification is vital for developing better treatment and prevention strategies. An important factor in the chronification of migraine is the overuse of acute headache medication. However, the mechanisms behind the transformation of episodic migraine to chronic migraine and vice versa have not yet been elucidated. We performed a longitudinal epigenome-wide association study to identify DNA methylation (DNAm) changes associated with treatment response in patients with chronic migraine and medication overuse as part of the Chronification and Reversibility of Migraine clinical trial. Blood was taken from patients with chronic migraine (n = 98) at baseline and after a 12-week medication withdrawal period. Treatment responders, patients with ≥ 50% reduction in monthly headache days (MHD), were compared with non-responders to identify DNAm changes associated with treatment response. Similarly, patients with ≥ 50% versus < 50% reduction in monthly migraine days (MMD) were compared.

RESULTS

At the epigenome-wide significant level (p < 9.42 × 10-8), a longitudinal reduction in DNAm at an intronic CpG site (cg14377273) within the HDAC4 gene was associated with MHD response following the withdrawal of acute medication. HDAC4 is highly expressed in the brain, plays a major role in synaptic plasticity, and modulates the expression and release of several neuroinflammation markers which have been implicated in migraine pathophysiology. Investigating whether baseline DNAm associated with treatment response, we identified lower baseline DNAm at a CpG site (cg15205829) within MARK3 that was significantly associated with MMD response at 12 weeks.

CONCLUSIONS

Our findings of a longitudinal reduction in HDAC4 DNAm status associated with treatment response and baseline MARK3 DNAm status as an early biomarker for treatment response, provide support for a role of pathways related to chromatin structure and synaptic plasticity in headache chronification and introduce HDAC4 and MARK3 as novel therapeutic targets.

Impact of tobacco, alcohol, and marijuana on genome-wide DNA methylation and its relationship with hypertension

Tobacco, alcohol, and marijuana consumption is an important public health problem because of their high use worldwide and their association with the risk of mortality and many health conditions, such as hypertension, which is the commonest risk factor for death throughout the world. A likely pathway of action of substance consumption leading to persistent hypertension is DNA methylation. Here, we evaluated the effects of tobacco, alcohol, and marijuana on DNA methylation in the same cohort (N = 3,424). Three epigenome-wide association studies (EWAS) were assessed in whole blood using the InfiniumHumanMethylationEPIC BeadChip. We also evaluated the mediation of the top CpG sites in the association between substance consumption and hypertension. Our analyses showed 2,569 CpG sites differentially methylated by alcohol drinking and 528 by tobacco smoking. We did not find significant associations with marijuana consumption after correcting for multiple comparisons. We found 61 genes overlapping between alcohol and tobacco that were enriched in biological processes involved in the nervous and cardiovascular systems. In the mediation analysis, we found 66 CpG sites that significantly mediated the effect of alcohol consumption on hypertension. The top alcohol-related CpG site (cg06690548, P-value = 5.9·10^-83) mapped to SLC7A11 strongly mediated 70.5% of the effect of alcohol consumption on hypertension (P-value = 0.006). Our findings suggest that DNA methylation should be considered for new targets in hypertension prevention and management, particularly concerning alcohol consumption. Our data also encourage further research into the use of methylation in blood to study the neurological and cardiovascular effects of substance consumption.

An epigenome-wide association study of child appetitive traits and DNA methylation

The etiology of childhood appetitive traits is poorly understood. Early-life epigenetic processes may be involved in the developmental programming of appetite regulation in childhood. One such process is DNA methylation (DNAm), whereby a methyl group is added to a specific part of DNA, where a cytosine base is next to a guanine base, a CpG site. We meta-analyzed epigenome-wide association studies (EWASs) of cord blood DNAm and early-childhood appetitive traits. Data were from two independent cohorts: the Generation R Study (n = 1,086, Rotterdam, the Netherlands) and the Healthy Start study (n = 236, Colorado, USA). DNAm at autosomal methylation sites in cord blood was measured using the Illumina Infinium HumanMethylation450 BeadChip. Parents reported on their child's food responsiveness, emotional undereating, satiety responsiveness and food fussiness using the Children's Eating Behaviour Questionnaire at age 4-5 years. Multiple regression models were used to examine the association of DNAm (predictor) at the individual site- and regional-level (using DMRff) with each appetitive trait (outcome), adjusting for covariates. Bonferroni-correction was applied to adjust for multiple testing. There were no associations of DNAm and any appetitive trait when examining individual CpG-sites. However, when examining multiple CpGs jointly in so-called differentially methylated regions, we identified 45 associations of DNAm with food responsiveness, 7 associations of DNAm with emotional undereating, 13 associations of DNAm with satiety responsiveness, and 9 associations of DNAm with food fussiness. This study shows that DNAm in the newborn may partially explain variation in appetitive traits expressed in early childhood and provides preliminary support for early programming of child appetitive traits through DNAm. Investigating differential DNAm associated with appetitive traits could be an important first step in identifying biological pathways underlying the development of these behaviors.

Air pollution and epigenetic aging among Black and White women in the US

BACKGROUND

DNA methylation-based measures of biological aging have been associated with air pollution and may link pollutant exposures to aging-related health outcomes. However, evidence is inconsistent and there is little information for Black women.

OBJECTIVE

We examined associations of ambient particulate matter <2.5 μm and <10 μm in diameter (PM2.5 and PM10) and nitrogen dioxide (NO2) with DNA methylation, including epigenetic aging and individual CpG sites, and evaluated whether associations differ between Black and non-Hispanic White (NHW) women.

METHODS

Validated models were used to estimate annual average outdoor residential exposure to PM2.5, PM10, and NO2 in a sample of self-identified Black (n=633) and NHW (n=3493) women residing in the contiguous US. We used sampling-weighted generalized linear regression to examine the effects of pollutants on six epigenetic aging measures (primary: DunedinPACE, GrimAgeAccel, and PhenoAgeAccel; secondary: Horvath intrinsic epigenetic age acceleration [EAA], Hannum extrinsic EAA, and skin & blood EAA) and epigenome-wide associations for individual CpG sites. Wald tests of nested models with and without interaction terms were used to examine effect measure modification by race/ethnicity.

RESULTS

Black participants had higher median air pollution exposure than NHW participants. GrimAgeAccel was associated with both PM10 and NO2 among Black participants, (Q4 versus Q1, PM10: β=1.09, 95% CI: 0.16-2.03; NO2: β=1.01, 95% CI 0.08-1.94) but not NHW participants (p-for-heterogeneity: PM10=0.10, NO2=0.20). In Black participants, we also observed a monotonic exposure-response relationship between NO2 and DunedinPACE (Q4 versus Q1, NO2: β=0.029, 95% CI: 0.004-0.055; p-for-trend=0.03), which was not observed in NHW participants (p-for-heterogeneity=0.09). In the EWAS, pollutants were significantly associated with differential methylation at 19 CpG sites in Black women and one in NHW women.

CONCLUSIONS

In a US-wide cohort study, our findings suggest that air pollution is associated with DNA methylation alterations consistent with higher epigenetic aging among Black, but not NHW, women.

Epigenomic signature of major congenital heart defects in newborns with Down syndrome

BACKGROUND

Congenital heart defects (CHDs) affect approximately half of individuals with Down syndrome (DS), but the molecular reasons for incomplete penetrance are unknown. Previous studies have largely focused on identifying genetic risk factors associated with CHDs in individuals with DS, but comprehensive studies of the contribution of epigenetic marks are lacking. We aimed to identify and characterize DNA methylation differences from newborn dried blood spots (NDBS) of DS individuals with major CHDs compared to DS individuals without CHDs.

METHODS

We used the Illumina EPIC array and whole-genome bisulfite sequencing (WGBS) to quantitate DNA methylation for 86 NDBS samples from the California Biobank Program: (1) 45 DS-CHD (27 female, 18 male) and (2) 41 DS non-CHD (27 female, 14 male). We analyzed global CpG methylation and identified differentially methylated regions (DMRs) in DS-CHD versus DS non-CHD comparisons (both sex-combined and sex-stratified) corrected for sex, age of blood collection, and cell-type proportions. CHD DMRs were analyzed for enrichment in CpG and genic contexts, chromatin states, and histone modifications by genomic coordinates and for gene ontology enrichment by gene mapping. DMRs were also tested in a replication dataset and compared to methylation levels in DS versus typical development (TD) WGBS NDBS samples.

RESULTS

We found global CpG hypomethylation in DS-CHD males compared to DS non-CHD males, which was attributable to elevated levels of nucleated red blood cells and not seen in females. At a regional level, we identified 58, 341, and 3938 CHD-associated DMRs in the Sex Combined, Females Only, and Males Only groups, respectively, and used machine learning algorithms to select 19 Males Only loci that could distinguish CHD from non-CHD. DMRs in all comparisons were enriched for gene exons, CpG islands, and bivalent chromatin and mapped to genes enriched for terms related to cardiac and immune functions. Lastly, a greater percentage of CHD-associated DMRs than background regions were differentially methylated in DS versus TD samples.

CONCLUSIONS

A sex-specific signature of DNA methylation was detected in NDBS of DS-CHD compared to DS non-CHD individuals. This supports the hypothesis that epigenetics can reflect the variability of phenotypes in DS, particularly CHDs.

Fathers' preconception smoking and offspring DNA methylation

BACKGROUND

Experimental studies suggest that exposures may impact respiratory health across generations via epigenetic changes transmitted specifically through male germ cells. Studies in humans are, however, limited. We aim to identify epigenetic marks in offspring associated with father's preconception smoking.

METHODS

We conducted epigenome-wide association studies (EWAS) in the RHINESSA cohort (7-50 years) on father's any preconception smoking (n = 875 offspring) and father's pubertal onset smoking < 15 years (n = 304), using Infinium MethylationEPIC Beadchip arrays, adjusting for offspring age, own smoking and maternal smoking. EWAS of maternal and offspring personal smoking were performed for comparison. Father's smoking-associated dmCpGs were checked in subpopulations of offspring who reported no personal smoking and no maternal smoking exposure.

RESULTS

Father's smoking commencing preconception was associated with methylation of blood DNA in offspring at two cytosine-phosphate-guanine sites (CpGs) (false discovery rate (FDR) < 0.05) in PRR5 and CENPP. Father's pubertal onset smoking was associated with 19 CpGs (FDR < 0.05) mapped to 14 genes (TLR9, DNTT, FAM53B, NCAPG2, PSTPIP2, MBIP, C2orf39, NTRK2, DNAJC14, CDO1, PRAP1, TPCN1, IRS1 and CSF1R). These differentially methylated sites were hypermethylated and associated with promoter regions capable of gene silencing. Some of these sites were associated with offspring outcomes in this cohort including ever-asthma (NTRK2), ever-wheezing (DNAJC14, TPCN1), weight (FAM53B, NTRK2) and BMI (FAM53B, NTRK2) (p < 0.05). Pathway analysis showed enrichment for gene ontology pathways including regulation of gene expression, inflammation and innate immune responses. Father's smoking-associated sites did not overlap with dmCpGs identified in EWAS of personal and maternal smoking (FDR < 0.05), and all sites remained significant (p < 0.05) in analyses of offspring with no personal smoking and no maternal smoking exposure.

CONCLUSION

Father's preconception smoking, particularly in puberty, is associated with offspring DNA methylation, providing evidence that epigenetic mechanisms may underlie epidemiological observations that pubertal paternal smoking increases risk of offspring asthma, low lung function and obesity.

Long-term exposure to ambient fine particulate components and leukocyte epigenome-wide DNA Methylation in older men: the Normative Aging Study

BACKGROUND

Epigenome-wide association studies of ambient fine particulate matter (PM2.5) have been reported. However, few have examined PM2.5 components (PMCs) and sources or included repeated measures. The lack of high-resolution exposure measurements is the key limitation. We hypothesized that significant changes in DNA methylation might vary by PMCs and the sources.

METHODS

We predicted the annual average of 14 PMCs using novel high-resolution exposure models across the contiguous U.S., between 2000-2018. The resolution was 50 m × 50 m in the Greater Boston Area. We also identified PM2.5 sources using positive matrix factorization. We repeatedly collected blood samples and measured leukocyte DNAm with the Illumina HumanMethylation450K BeadChip in the Normative Aging Study. We then used median regression with subject-specific intercepts to estimate the associations between long-term (one-year) exposure to PMCs / PM2.5 sources and DNA methylation at individual cytosine-phosphate-guanine CpG sites. Significant probes were identified by the number of independent degrees of freedom approach, using the number of principal components explaining > 95% of the variation of the DNA methylation data. We also performed regional and pathway analyses to identify significant regions and pathways.

RESULTS

We included 669 men with 1,178 visits between 2000-2013. The subjects had a mean age of 75 years. The identified probes, regions, and pathways varied by PMCs and their sources. For example, iron was associated with 6 probes and 6 regions, whereas nitrate was associated with 15 probes and 3 regions. The identified pathways from biomass burning, coal burning, and heavy fuel oil combustion sources were associated with cancer, inflammation, and cardiovascular diseases, whereas there were no pathways associated with all traffic.

CONCLUSIONS

Our findings showed that the effects of PM2.5 on DNAm varied by its PMCs and sources.

An epigenome-wide association study of child appetitive traits and DNA methylation

Childhood appetitive traits are consistently associated with obesity risk, and yet their etiology is poorly understood. Appetitive traits are complex phenotypes which are hypothesized to be influenced by both genetic and environmental factors, as well as their interactions. Early-life epigenetic processes, such as DNA methylation (DNAm), may be involved in the developmental programming of appetite regulation in childhood. In the current study, we meta-analyzed epigenome-wide association studies (EWASs) of cord blood DNAm and early-childhood appetitive traits. Data were from two independent cohorts: the Generation R Study (n=1,086, Rotterdam, the Netherlands) and the Healthy Start study (n=236, Colorado, USA). DNAm at autosomal methylation sites in cord blood was measured using the Illumina Infinium HumanMethylation450 BeadChip. Parents reported on their child's food responsiveness, emotional undereating, satiety responsiveness and food fussiness using the Children's Eating Behaviour Questionnaire at age 4-5 years. Multiple regression models were used to examine the association of DNAm (predictor) at the individual site- and regional-level (using DMRff) with each appetitive trait (outcome), adjusting for covariates. Bonferroni-correction was applied to adjust for multiple testing. There were no associations of DNAm and any appetitive trait at the individual site-level. However, at the regional level, we identified 45 associations of DNAm with food responsiveness, 7 associations of DNAm with emotional undereating, 13 associations of DNAm with satiety responsiveness, and 9 associations of DNAm with food fussiness. This study shows that DNAm in the newborn may partially explain variation in appetitive traits expressed in early childhood and provides preliminary support for early programming of child appetitive traits through DNAm. Investigating differential DNAm associated with appetitive traits could be an important first step in identifying biological pathways underlying the development of these behaviors.

Epigenetic marks associated with gestational diabetes mellitus across two time points during pregnancy

An adverse intrauterine or periconceptional environment, such as hyperglycemia during pregnancy, can affect the DNA methylation pattern both in mothers and their offspring. In this study, we explored the epigenetic profile in maternal peripheral blood samples through pregnancy to find potential epigenetic biomarkers for gestational diabetes mellitus (GDM), as well as candidate genes involved in GDM development. We performed an epigenome-wide association study in maternal peripheral blood samples in 32 pregnant women (16 with GDM and 16 non-GDM) at pregnancy week 24-28 and 36-38. Biochemical, anthropometric, and obstetrical variables were collected from all the participants. The main results were validated in an independent cohort with different ethnic origin (European = 307; South Asians = 165). Two hundred and seventy-two CpGs sites remained significantly different between GDM and non-GDM pregnant women across two time points during pregnancy. The significant CpG sites were related to pathways associated with type I diabetes mellitus, insulin resistance and secretion. Cg01459453 (SELP gene) was the most differentiated in the GDM group versus non-GDM (73.6 vs. 60.9, p = 1.06E-11; FDR = 7.87E-06). Three CpG sites (cg01459453, cg15329406, and cg04095097) were able to discriminate between GDM cases and controls (AUC = 1; p = 1.26E-09). Three differentially methylated positions (DMPs) were replicated in an independent cohort. To conclude, epigenetic marks during pregnancy differed between GDM cases and controls suggesting a role for these genes in GDM development. Three CpGs were able to discriminate GDM and non-GDM groups with high specificity and sensitivity, which may be biomarker candidates for diagnosis or prediction of GDM.

DNA methylation is differentially associated with glycemic outcomes by different types of weight-loss interventions: an epigenome-wide association study

BACKGROUND

Alterations in DNA methylation (DNAm) have been reported to be a mechanism by which bariatric surgeries resulted in considerable metabolic improvements. Previous studies have mostly focused on change in DNAm following weight-loss interventions, yet whether DNAm prior to intervention can explain the variability in glycemic outcomes has not been investigated. Here, we aim to examine whether baseline DNAm is differentially associated with glycemic outcomes induced by different types of weight-loss interventions.

METHODS

Participants were 75 adults with severe obesity who underwent non-surgical intensive medical intervention (IMI), adjustable gastric band (BAND) or Roux-en-Y gastric bypass (RYGB) (n = 25 each). Changes in fasting plasma glucose (FPG) and glycated hemoglobin (HbA1c) were measured at 1-year after intervention. DNAm was quantified by Illumina 450 K arrays in baseline peripheral blood DNA. Epigenome-wide association studies were performed to identify CpG probes that modify the effects of different weight-loss interventions on glycemic outcomes, i.e., changes in FPG and HbA1c, by including an interaction term between types of intervention and DNAm. Models were adjusted for weight loss and baseline clinical factors.

RESULTS

Baseline DNAm levels at 3216 and 117 CpGs were differentially associated with changes in FPG and HbA1c, respectively, when comparing RYGB versus IMI. Of these, 79 CpGs were significant for both FPG and HbA1c. The identified genes are enriched in adaptive thermogenesis, temperature homeostasis and regulation of cell population proliferation. Additionally, DNAm at 6 CpGs was differentially associated with changes in HbA1c when comparing RYGB versus BAND.

CONCLUSIONS

Baseline DNAm is differentially associated with glycemic outcomes in response to different types of weight-loss interventions, independent of weight loss and other clinical factors. Such findings provided initial evidence that baseline DNAm levels may serve as potential biomarkers predictive of differential glycemic outcomes in response to different types of weight-loss interventions.

Epigenome-wide association study identifies Vogt-Koyanagi-Harada disease-associated DNA methylation loci in Chinese

DNA methylation is one of the important epigenetic mechanisms for modulating gene expression. By performing a genome-wide methylation association analysis of whole peripheral blood from 60 Vogt-Koyanagi-Harada disease (VKH) patients and 60 healthy controls, we depicted the global DNA methylation status of VKH disease. Further pyrosequencing validation in 160 patients and 159 controls identified 3 aberrant CpG sites in HLA gene regions including cg04026937 and cg18052547 (located in HLA-DRB1 region), and cg13778567 (HLA-DQA1). We also identified 9 aberrant CpG sites in non-HLA gene regions including cg13979407, cg21075643, cg24290586, cg10135747 and cg22707857 (BTNL2), cg22155039 (NOTCH4), cg02605387 (TNXB), cg06255004 (AGPAT2) and cg18855195 (RIBC2). Increased mRNA levels of BTNL2, NOTCH4 and TNXB were identified in VKH patients when compared with healthy controls, consistent with the hypomethylated CpG status in these gene regions. Moreover, seven aberrantly methylated CpG sites may serve as a diagnostic marker for VKH disease (AUC = 84.95%, 95%CI: 79.49%-90.41%).

Epigenomic profiling of isolated blood cell types reveals highly specific B cell smoking signatures and links to disease risk

BACKGROUND

Tobacco smoking alters the DNA methylation profiles of immune cells which may underpin some of the pathogenesis of smoking-associated diseases. To link smoking-driven epigenetic effects in specific immune cell types with disease risk, we isolated six leukocyte subtypes, CD14+ monocytes, CD15+ granulocytes, CD19+ B cells, CD4+ T cells, CD8+ T cells, and CD56+ natural killer cells, from whole blood of 67 healthy adult smokers and 74 nonsmokers for epigenome-wide association study (EWAS) using Illumina 450k and EPIC methylation arrays.

RESULTS

Numbers of smoking-associated differentially methylated sites (smCpGs) at genome-wide significance (p < 1.2 × 10-7) varied widely across cell types, from 5 smCpGs in CD8+ T cells to 111 smCpGs in CD19+ B cells. We found unique smoking effects in each cell type, some of which were not apparent in whole blood. Methylation-based deconvolution to estimate B cell subtypes revealed that smokers had 7.2% (p = 0.033) less naïve B cells. Adjusting for naïve and memory B cell proportions in EWAS and RNA-seq allowed the identification of genes enriched for B cell activation-related cytokine signaling pathways, Th1/Th2 responses, and hematopoietic cancers. Integrating with large-scale public datasets, 62 smCpGs were among CpGs associated with health-relevant EWASs. Furthermore, 74 smCpGs had reproducible methylation quantitative trait loci single nucleotide polymorphisms (SNPs) that were in complete linkage disequilibrium with genome-wide association study SNPs, associating with lung function, disease risks, and other traits.

CONCLUSIONS

We observed blood cell-type-specific smCpGs, a naïve-to-memory shift among B cells, and by integrating genome-wide datasets, we identified their potential links to disease risks and health traits.

Entorhinal cortex epigenome-wide association study highlights four novel loci showing differential methylation in Alzheimer's disease

BACKGROUND

Studies on DNA methylation (DNAm) in Alzheimer's disease (AD) have recently highlighted several genomic loci showing association with disease onset and progression.

METHODS

Here, we conducted an epigenome-wide association study (EWAS) using DNAm profiles in entorhinal cortex (EC) from 149 AD patients and control brains and combined these with two previously published EC datasets by meta-analysis (total n = 337).

RESULTS

We identified 12 cytosine-phosphate-guanine (CpG) sites showing epigenome-wide significant association with either case-control status or Braak's tau-staging. Four of these CpGs, located in proximity to CNFN/LIPE, TENT5A, PALD1/PRF1, and DIRAS1, represent novel findings. Integrating DNAm levels with RNA sequencing-based mRNA expression data generated in the same individuals showed significant DNAm-mRNA correlations for 6 of the 12 significant CpGs. Lastly, by calculating rates of epigenetic age acceleration using two recently proposed "epigenetic clock" estimators we found a significant association with accelerated epigenetic aging in the brains of AD patients vs. controls.

CONCLUSION

In summary, our study represents the hitherto most comprehensive EWAS in AD using EC and highlights several novel differentially methylated loci with potential effects on gene expression.

The X-factor in ART: does the use of assisted reproductive technologies influence DNA methylation on the X chromosome?

BACKGROUND

Assisted reproductive technologies (ART) may perturb DNA methylation (DNAm) in early embryonic development. Although a handful of epigenome-wide association studies of ART have been published, none have investigated CpGs on the X chromosome. To bridge this knowledge gap, we leveraged one of the largest collections of mother-father-newborn trios of ART and non-ART (natural) conceptions to date to investigate sex-specific DNAm differences on the X chromosome. The discovery cohort consisted of 982 ART and 963 non-ART trios from the Norwegian Mother, Father, and Child Cohort Study (MoBa). To verify our results from the MoBa cohort, we used an external cohort of 149 ART and 58 non-ART neonates from the Australian 'Clinical review of the Health of adults conceived following Assisted Reproductive Technologies' (CHART) study. The Illumina EPIC array was used to measure DNAm in both datasets. In the MoBa cohort, we performed a set of X-chromosome-wide association studies ('XWASs' hereafter) to search for sex-specific DNAm differences between ART and non-ART newborns. We tested several models to investigate the influence of various confounders, including parental DNAm. We also searched for differentially methylated regions (DMRs) and regions of co-methylation flanking the most significant CpGs. Additionally, we ran an analogous model to our main model on the external CHART dataset.

RESULTS

In the MoBa cohort, we found more differentially methylated CpGs and DMRs in girls than boys. Most of the associations persisted after controlling for parental DNAm and other confounders. Many of the significant CpGs and DMRs were in gene-promoter regions, and several of the genes linked to these CpGs are expressed in tissues relevant for both ART and sex (testis, placenta, and fallopian tube). We found no support for parental DNAm-dependent features as an explanation for the observed associations in the newborns. The most significant CpG in the boys-only analysis was in UBE2DNL, which is expressed in testes but with unknown function. The most significant CpGs in the girls-only analysis were in EIF2S3 and AMOT. These three loci also displayed differential DNAm in the CHART cohort.

CONCLUSIONS

Genes that co-localized with the significant CpGs and DMRs associated with ART are implicated in several key biological processes (e.g., neurodevelopment) and disorders (e.g., intellectual disability and autism). These connections are particularly compelling in light of previous findings indicating that neurodevelopmental outcomes differ in ART-conceived children compared to those naturally conceived.

Epigenome-wide association study of plasma lipids in West Africans: the RODAM study

BACKGROUND

DNA-methylation has been associated with plasma lipid concentration in populations of diverse ethnic backgrounds, but epigenome-wide association studies (EWAS) in West-Africans are lacking. The aim of this study was to identify DNA-methylation loci associated with plasma lipids in Ghanaians.

METHODS

We conducted an EWAS using Illumina 450k DNA-methylation array profiles of extracted DNA from 663 Ghanaian participants. Differentially methylated positions (DMPs) were examined for association with plasma total cholesterol (TC), LDL-cholesterol, HDL-cholesterol, and triglycerides concentrations using linear regression models adjusted for age, sex, body mass index, diabetes mellitus, and technical covariates. Findings were replicated in independent cohorts of different ethnicities.

FINDINGS

We identified one significantly associated DMP with triglycerides (cg19693031 annotated to TXNIP, regression coefficient beta -0.26, false discovery rate adjusted p-value 0.001), which replicated in-silico in South African Batswana, African American, and European populations. From the top five DMPs with the lowest nominal p-values, two additional DMPs for triglycerides (CPT1A, ABCG1), two DMPs for LDL-cholesterol (EPSTI1, cg13781819), and one for TC (TXNIP) replicated. With the exception of EPSTI1, these loci are involved in lipid transport/metabolism or are known GWAS-associated loci. The top 5 DMPs per lipid trait explained 9.5% in the variance of TC, 8.3% in LDL-cholesterol, 6.1% in HDL-cholesterol, and 11.0% in triglycerides.

INTERPRETATION

The top DMPs identified in this study are in loci that play a role in lipid metabolism across populations, including West-Africans. Future studies including larger sample size, longitudinal study design and translational research is needed to increase our understanding on the epigenetic regulation of lipid metabolism among West-African populations.

FUNDING

European Commission under the Framework Programme (grant number: 278901).

Associations of DNA Methylation With Behavioral Problems, Gray Matter Volumes, and Negative Life Events Across Adolescence: Evidence From the Longitudinal IMAGEN Study

BACKGROUND

Negative life events (NLEs) increase the risk for externalizing behaviors (EBs) and internalizing behaviors (IBs) in adolescence and adult psychopathology. DNA methylation associated with behavioral problems may reflect this risk and long-lasting effects of NLEs.

METHODS

To identify consistent associations between blood DNA methylation and EBs or IBs across adolescence, we conducted longitudinal epigenome-wide association studies (EWASs) using data from the IMAGEN cohort, collected at ages 14 and 19 years (n = 506). Significant findings were validated in a separate subsample (n = 823). Methylation risk scores were generated by 10-fold cross-validation and further tested for their associations with gray matter volumes and NLEs.

RESULTS

No significant findings were obtained for the IB-EWAS. The EB-EWAS identified a genome-wide significant locus in a gene linked to attention-deficit/hyperactivity disorder (ADHD) (IQSEC1, cg01460382; p = 1.26 × 10-8). Other most significant CpG sites were near ADHD-related genes and enriched for genes regulating tumor necrosis factor and interferon-γ signaling, highlighting the relevance of EB-EWAS findings for ADHD. Analyses with the EB methylation risk scores suggested that it partly reflected comorbidity with IBs in late adolescence. Specific to EBs, EB methylation risk scores correlated with smaller gray matter volumes in medial orbitofrontal and anterior/middle cingulate cortices, brain regions known to associate with ADHD and conduct problems. Longitudinal mediation analyses indicated that EB-related DNA methylation were more likely the outcomes of problematic behaviors accentuated by NLEs, and less likely the epigenetic bases of such behaviors.

CONCLUSIONS

Our findings suggest that novel epigenetic mechanisms through which NLEs exert short and longer-term effects on behavior may contribute to ADHD.

Epigenome-wide DNA methylation in externalizing behaviours: A review and combined analysis

DNA methylation (DNAm) is one of the most frequently studied epigenetic mechanisms facilitating the interplay of genomic and environmental factors, which can contribute to externalizing behaviours and related psychiatric disorders. Previous epigenome-wide association studies (EWAS) for externalizing behaviours have been limited in sample size, and, therefore, candidate genes and biomarkers with robust evidence are still lacking. We 1) performed a systematic literature review of EWAS of attention-deficit/hyperactivity disorder (ADHD)- and aggression-related behaviours conducted in peripheral tissue and cord blood and 2) combined the most strongly associated DNAm sites observed in individual studies (p < 10^-3) to identify candidate genes and biological systems for ADHD and aggressive behaviours. We observed enrichment for neuronal processes and neuronal cell marker genes for ADHD. Astrocyte and granulocytes cell markers among genes annotated to DNAm sites were relevant for both ADHD and aggression-related behaviours. Only 1 % of the most significant epigenetic findings for ADHD/ADHD symptoms were likely to be directly explained by genetic factors involved in ADHD. Finally, we discuss how the field would greatly benefit from larger sample sizes and harmonization of assessment instruments.

Epigenome-wide DNA methylation in leukocytes and toenail metals: The normative aging study

BACKGROUND

Environmental metal exposures have been associated with multiple deleterious health endpoints. DNA methylation (DNAm) may provide insight into the mechanisms underlying these relationships. Toenail metals are non-invasive biomarkers, reflecting a medium-term time exposure window.

OBJECTIVES

This study examined variation in leukocyte DNAm and toenail arsenic (As), cadmium (Cd), lead (Pb), manganese (Mn), and mercury (Hg) among elderly men in the Normative Aging Study, a longitudinal cohort.

METHODS

We repeatedly collected samples of blood and toenail clippings. We measured DNAm in leukocytes with the Illumina HumanMethylation450 K BeadChip. We first performed median regression to evaluate the effects of each individual toenail metal on DNAm at three levels: individual cytosine-phosphate-guanine (CpG) sites, regions, and pathways. Then, we applied a Bayesian kernel machine regression (BKMR) to assess the joint and individual effects of metal mixtures on DNAm. Significant CpGs were identified using a multiple testing correction based on the independent degrees of freedom approach for correlated outcomes. The approach considers the effective degrees of freedom in the DNAm data using the principal components that explain >95% variation of the data.

RESULTS

We included 564 subjects (754 visits) between 1999 and 2013. The numbers of significantly differentially methylated CpG sites, regions, and pathways varied by metals. For example, we found six significant pathways for As, three for Cd, and one for Mn. The As-associated pathways were associated with cancer (e.g., skin cancer) and cardiovascular disease, whereas the Cd-associated pathways were related to lung cancer. Metal mixtures were also associated with 47 significant CpG sites, as well as pathways, mainly related to cancer and cardiovascular disease.

CONCLUSIONS

This study provides an approach to understanding the potential epigenetic mechanisms underlying observed relations between toenail metals and adverse health endpoints.

Systemic interindividual epigenetic variation in humans is associated with transposable elements and under strong genetic control

BACKGROUND

Genetic variants can modulate phenotypic outcomes via epigenetic intermediates, for example at methylation quantitative trait loci (mQTL). We present the first large-scale assessment of mQTL at human genomic regions selected for interindividual variation in CpG methylation, which we call correlated regions of systemic interindividual variation (CoRSIVs). These can be assayed in blood DNA and do not reflect interindividual variation in cellular composition.

RESULTS

We use target-capture bisulfite sequencing to assess DNA methylation at 4086 CoRSIVs in multiple tissues from each of 188 donors in the NIH Gene-Tissue Expression (GTEx) program. At CoRSIVs, DNA methylation in peripheral blood correlates with methylation and gene expression in internal organs. We also discover unprecedented mQTL at these regions. Genetic influences on CoRSIV methylation are extremely strong (median R²=0.76), cumulatively comprising over 70-fold more human mQTL than detected in the most powerful previous study. Moreover, mQTL beta coefficients at CoRSIVs are highly skewed (i.e., the major allele predicts higher methylation). Both surprising findings are independently validated in a cohort of 47 non-GTEx individuals. Genomic regions flanking CoRSIVs show long-range enrichments for LINE-1 and LTR transposable elements; the skewed beta coefficients may therefore reflect evolutionary selection of genetic variants that promote their methylation and silencing. Analyses of GWAS summary statistics show that mQTL polymorphisms at CoRSIVs are associated with metabolic and other classes of disease.

CONCLUSIONS

A focus on systemic interindividual epigenetic variants, clearly enhanced in mQTL content, should likewise benefit studies attempting to link human epigenetic variation to the risk of disease.

PM2.5-related DNA methylation and the association with lung function in non-smokers

PM_2.5 exposure leads to lung function alteration. The potential pathway underlying above association, especially the role of DNA methylation is unclear. The objectives of this study are to evaluate the associations of personal PM_2.5 concentrations with DNA methylation at the epigenome-wide level, and investigate how PM_2.5-related DNA methylation affects lung function. A total of 402 observations of non-smokers were selected from the Wuhan-Zhuhai cohort. PM_2.5 exposure was estimated through a model established in the same population. Blood DNA methylation levels were determined through Illumina Infinium MethylationEPIC BeadChips. Lung function was tested through spirometry on the day of blood sampling. The associations of PM_2.5 exposure with DNA methylation and DNA methylation with lung function were determined through linear mixed models. Ten PM_2.5-related CpG sites (mapped to 7 different genes) were observed with false discovery rate <0.05. Methylation levels of cg24821877, cg24862131, cg23530876, cg11149743 and cg10781276 were positively associated with PM_2.5 concentrations. While methylation levels of cg10314909, cg08968107, cg18362281, cg24663971 and cg17834632 were negatively associated with PM_2.5 concentrations. The top CpG was cg24663971 (P = 1.51✕10^-9). Among the above 10 sites, significantly positive associations of methylation levels of cg24663971 with FVC%pred and FEV₁%pred, and cg10314909 with FVC, FVC%pred, and FEV₁%pred were observed. Age had modification effect on the associations between cg24663971 methylation and FVC%pred, and the associations were more obvious among participants with age ≥58 years. In conclusion, PM_2.5 exposure was associated with DNA methylation, and PM_2.5-related DNA methylation was associated with lung function among Wuhan urban non-smokers.

Wildfire-related PM2.5 and DNA methylation: An Australian twin and family study

BACKGROUND

Wildfire-related fine particulate matter (PM_2.5) has many adverse health impacts, but its impacts on human epigenome are unknown. We aimed to evaluate the associations between long-term exposure to wildfire-related PM_2.5 and blood DNA methylation, and whether the associations differ from those with non-wildfire-related PM_2.5.

METHODS

We studied 479 Australian women comprising 132 twin pairs and 215 of their sisters. Blood-derived DNA methylation was measured using the HumanMethylation450 BeadChip array. Data on 3-year (year of blood collection and previous two years) average wildfire-related and non-wildfire-related PM_2.5 at 0.01°×0.01° spatial resolution were created by combining information from satellite observations, chemical transport models, and ground-based observations. Exposure data were linked to each participant's home address, assuming the address did not change during the exposure window. For DNA methylation of each cytosine-guanine dinucleotide (CpG), and for global DNA methylation represented by the average of all measured CpGs or CpGs in repetitive elements, we evaluated their associations with wildfire- or non-wildfire-related PM_2.5 using a within-sibship analysis controlling for factors shared between siblings and other important covariates. Differentially methylated regions (DMRs) were defined by comb-p and DMRcate.

RESULTS

The 3-year average wildfire-related PM_2.5 (range: 0.3 to 7.6 µg/m³_, mean: 1.6 µg/m³) was negatively, but not significantly (p-values greater than 0.05) associated with all seven global DNA methylation measures. There were 26 CpGs and 33 DMRs associated with wildfire-related PM_2.5 (Bonferroni adjusted p-value < 0.05) mapped to 47 genes enriched for pathways related to inflammatory regulation and platelet activation. These genes have been related to many human diseases or phenotypes e.g., cancer, mental disorders, diabetes, obesity, asthma, blood pressure. These CpGs, DMRs and enriched pathways did not overlap with the 1 CpG and 7 DMRs associated with non-wildfire-related PM_2.5.

CONCLUSIONS

Long-term exposure to wildfire-related PM_2.5 was associated with various blood DNA methylation signatures in Australian women, and these were distinct from those associated with non-wildfire-related PM_2.5.

Whole-blood methylation signatures are associated with and accurately classify multiple sclerosis disease severity

BACKGROUND

The variation in multiple sclerosis (MS) disease severity is incompletely explained by genetics, suggesting genetic and environmental interactions are involved. Moreover, the lack of prognostic biomarkers makes it difficult for clinicians to optimise care. DNA methylation is one epigenetic mechanism by which gene-environment interactions can be assessed. Here, we aimed to identify DNA methylation patterns associated with mild and severe relapse-onset MS (RMS) and to test the utility of methylation as a predictive biomarker.

METHODS

We conducted an epigenome-wide association study between 235 females with mild (n = 119) or severe (n = 116) with RMS. Methylation was measured with the Illumina methylationEPIC array and analysed using logistic regression. To generate hypotheses about the functional consequence of differential methylation, we conducted gene set enrichment analysis using ToppGene. We compared the accuracy of three machine learning models in classifying disease severity: (1) clinical data available at baseline (age at onset and first symptoms) built using elastic net (EN) regression, (2) methylation data using EN regression and (3) a weighted methylation risk score of differentially methylated positions (DMPs) from the main analysis using logistic regression. We used a conservative 70:30 test:train split for classification modelling. A false discovery rate threshold of 0.05 was used to assess statistical significance.

RESULTS

Females with mild or severe RMS had 1472 DMPs in whole blood (839 hypermethylated, 633 hypomethylated in the severe group). Differential methylation was enriched in genes related to neuronal cellular compartments and processes, and B-cell receptor signalling. Whole-blood methylation levels at 1708 correlated CpG sites classified disease severity more accurately (machine learning model 2, AUC = 0.91) than clinical data (model 1, AUC = 0.74) or the wMRS (model 3, AUC = 0.77). Of the 1708 selected CpGs, 100 overlapped with DMPs from the main analysis at the gene level. These overlapping genes were enriched in neuron projection and dendrite extension, lending support to our finding that neuronal processes, rather than immune processes, are implicated in disease severity.

CONCLUSION

RMS disease severity is associated with whole-blood methylation at genes related to neuronal structure and function. Moreover, correlated whole-blood methylation patterns can assign disease severity in females with RMS more accurately than clinical data available at diagnosis.

An explorative epigenome-wide association study of plasma renin and aldosterone concentration in a Ghanaian population: the RODAM study

BACKGROUND

The epigenetic regulation of the renin-angiotensin-aldosterone system (RAAS) potentially plays a role in the pathophysiology underlying the high burden of hypertension in sub-Saharan Africans (SSA). Here we report the first epigenome-wide association study (EWAS) of plasma renin and aldosterone concentrations and the aldosterone-to-renin ratio (ARR).

METHODS

Epigenome-wide DNA methylation was measured using the Illumina 450K array on whole blood samples of 68 Ghanaians. Differentially methylated positions (DMPs) were assessed for plasma renin concentration, aldosterone, and ARR using linear regression models adjusted for age, sex, body mass index, diabetes mellitus, hypertension, and technical covariates. Additionally, we extracted methylation loci previously associated with hypertension, kidney function, or that were annotated to RAAS-related genes and associated these with renin and aldosterone concentration.

RESULTS

We identified one DMP for renin, ten DMPs for aldosterone, and one DMP associated with ARR. Top DMPs were annotated to the PTPRN2, SKIL, and KCNT1 genes, which have been reported in relation to cardiometabolic risk factors, atherosclerosis, and sodium-potassium handling. Moreover, EWAS loci previously associated with hypertension, kidney function, or RAAS-related genes were also associated with renin, aldosterone, and ARR.

CONCLUSION

In this first EWAS on RAAS hormones, we identified DMPs associated with renin, aldosterone, and ARR in a SSA population. These findings are a first step in understanding the role of DNA methylation in regulation of the RAAS in general and in a SSA population specifically. Replication and translational studies are needed to establish the role of these DMPs in the hypertension burden in SSA populations.

Maternal blood pressure associates with placental DNA methylation both directly and through alterations in cell-type composition

BACKGROUND

Maternal blood pressure levels reflect cardiovascular adaptation to pregnancy and proper maternal-fetal exchanges through the placenta and are very sensitive to numerous environmental stressors. Maternal hypertension during pregnancy has been associated with impaired placental functions and with an increased risk for children to suffer from cardiovascular and respiratory diseases later on. Investigating changes in placental DNA methylation levels and cell-type composition in association with maternal blood pressure could help elucidate its relationships with placental and fetal development.

METHODS

Taking advantage of a large cohort of 666 participants, we investigated the association between epigenome-wide DNA methylation patterns in the placenta, measured using the Infinium HumanMethylation450 BeadChip, placental cell-type composition, estimated in silico, and repeated measurements of maternal steady and pulsatile blood pressure indicators during pregnancy.

RESULTS

At the site-specific level, no significant association was found between maternal blood pressure and DNA methylation levels after correction for multiple testing (false discovery rate < 0.05), but 5 out of 24 previously found CpG associations were replicated (p-value < 0.05). At the regional level, our analyses highlighted 64 differentially methylated regions significantly associated with at least one blood pressure component, including 35 regions associated with mean arterial pressure levels during late pregnancy. These regions were found enriched for genes implicated in lung development and diseases. Further mediation analyses show that a significant part of the association between steady blood pressure-but not pulsatile pressure-and placental methylation can be explained by alterations in placental cell-type composition. In particular, elevated blood pressure levels are associated with a decrease in the ratio between mesenchymal stromal cells and syncytiotrophoblasts, even in the absence of preeclampsia.

CONCLUSIONS

This study provides the first evidence that the association between maternal steady blood pressure during pregnancy and placental DNA methylation is both direct and partly explained by changes in cell-type composition. These results could hint at molecular mechanisms linking maternal hypertension to lung development and early origins of childhood respiratory problems and at the importance of controlling maternal blood pressure during pregnancy.

Gaseous air pollutants and DNA methylation in a methylome-wide association study of an ethnically and environmentally diverse population of U.S. adults

Epigenetic mechanisms may underlie air pollution-health outcome associations. We estimated gaseous air pollutant-DNA methylation (DNAm) associations using twelve subpopulations within Women's Health Initiative (WHI) and Atherosclerosis Risk in Communities (ARIC) cohorts (n = 8397; mean age 61.3 years; 83% female; 46% African-American, 46% European-American, 8% Hispanic/Latino). We used geocoded participant address-specific mean ambient carbon monoxide (CO), nitrogen oxides (NO2; NOx), ozone (O3), and sulfur dioxide (SO2) concentrations estimated over the 2-, 7-, 28-, and 365-day periods before collection of blood samples used to generate Illumina 450 k array leukocyte DNAm measurements. We estimated methylome-wide, subpopulation- and race/ethnicity-stratified pollutant-DNAm associations in multi-level, linear mixed-effects models adjusted for sociodemographic, behavioral, meteorological, and technical covariates. We combined stratum-specific estimates in inverse variance-weighted meta-analyses and characterized significant associations (false discovery rate; FDR<0.05) at Cytosine-phosphate-Guanine (CpG) sites without among-strata heterogeneity (PCochran's Q > 0.05). We attempted replication in the Cooperative Health Research in Region of Augsburg (KORA) study and Normative Aging Study (NAS). We observed a -0.3 (95% CI: -0.4, -0.2) unit decrease in percent DNAm per interquartile range (IQR, 7.3 ppb) increase in 28-day mean NO2 concentration at cg01885635 (chromosome 3; regulatory region 290 bp upstream from ZNF621; FDR = 0.03). At intragenic sites cg21849932 (chromosome 20; LIME1; intron 3) and cg05353869 (chromosome 11; KLHL35; exon 2), we observed a -0.3 (95% CI: -0.4, -0.2) unit decrease (FDR = 0.04) and a 1.2 (95% CI: 0.7, 1.7) unit increase (FDR = 0.04), respectively, in percent DNAm per IQR (17.6 ppb) increase in 7-day mean ozone concentration. Results were not fully replicated in KORA and NAS. We identified three CpG sites potentially susceptible to gaseous air pollution-induced DNAm changes near genes relevant for cardiovascular and lung disease. Further harmonized investigations with a range of gaseous pollutants and averaging durations are needed to determine the effect of gaseous air pollutants on DNA methylation and ultimately gene expression.

Cell Type-Specific Signal Analysis in Epigenome-Wide Association Studies

Hundreds of epigenome-wide association studies (EWAS) have been performed, successfully identifying replicated epigenomic signals in processes such as aging and smoking. Despite this progress, it remains a major challenge in EWAS to detect both cell type-specific and cell type confounding effects impacting study results. One way to identify these effects is through eFORGE (experimentally derived Functional element Overlap analysis of ReGions from EWAS), a published tool that uses 815 datasets from large-scale mapping studies to detect enriched tissues, cell types, and genomic regions. Here, I show that eFORGE analysis can be extended to EWAS differentially variable positions (DVPs), identifying target cell types and tissues. In addition, I also show that eFORGE tissue-specific enrichment can be detected for sites below EWAS significance threshold. I develop on these and other analysis examples, extending our knowledge of eFORGE cell type- and tissue-specific enrichment results for different EWAS.

Epigenome-wide association study of bronchopulmonary dysplasia in preterm infants: results from the discovery-BPD program

BACKGROUND

Bronchopulmonary dysplasia (BPD) is a lung disease in premature infants caused by therapeutic oxygen supplemental and characterized by impaired pulmonary development which persists into later life. While advances in neonatal care have improved survival rates of premature infants, cases of BPD have been increasing with limited therapeutic options for prevention and treatment. This study was designed to explore the relationship between gestational age (GA), birth weight, and estimated blood cell-type composition in premature infants and to elucidate early epigenetic biomarkers associated with BPD.

METHODS

Cord blood DNA from preterm neonates that went on to develop BPD (n = 14) or not (non-BPD, n = 93) was applied to Illumina 450 K methylation arrays. Blood cell-type compositions were estimated using DNA methylation profiles. Multivariable robust regression analysis elucidated CpGs associated with BPD risk. cDNA microarray analysis of cord blood RNA identified differentially expressed genes in neonates who later developed BPD.

RESULTS

The development of BPD and the need for oxygen supplementation were strongly associated with GA (BPD, p < 1.0E-04; O2 supplementation, p < 1.0E-09) and birth weight (BPD, p < 1.0E-02; O2 supplementation, p < 1.0E-07). The estimated nucleated red blood cell (NRBC) percent was negatively associated with birth weight and GA, positively associated with hypomethylation of the tobacco smoke exposure biomarker cg05575921, and high-NRBC blood samples displayed a hypomethylation profile. Epigenome-wide association study (EWAS) identified 38 (Bonferroni) and 275 (false discovery rate 1%) differentially methylated CpGs associated with BPD. BPD-associated CpGs in cord blood were enriched for lung maturation and hematopoiesis pathways. Stochastic epigenetic mutation burden at birth was significantly elevated among those who developed BPD (adjusted p = 0.02). Transcriptome changes in cord blood cells reflected cell cycle, development, and pulmonary disorder events in BPD.

CONCLUSIONS

While results must be interpreted with caution because of the small size of this study, NRBC content strongly impacted DNA methylation profiles in preterm cord blood and EWAS analysis revealed potential insights into biological pathways involved in BPD pathogenesis.

Clinical Implication of Smoking-Related Aryl-Hydrocarbon Receptor Repressor (AHRR) Hypomethylation in Japanese Adults

BACKGROUND

Tobacco smoking is a leading preventable cause of morbidity and mortality worldwide; still, the success rate of smoking cessation is low in general. From the viewpoint of public health and clinical care, an objective biomarker of long-term smoking behavior is sought.Methods and Results:This study assessed DNA methylation as a biomarker of smoking in a hospital setting through a combination of molecular approaches including genetic, DNA methylation and mRNA expression analyses. First, in an epigenome-wide association study involving Japanese individuals with chronic cardiovascular disease (n=94), genome-wide significant smoking association was identified at 2 CpG sites on chromosome 5, with the strongest signal at cg05575921 located in intron 3 of the aryl-hydrocarbon receptor repressor (AHRR) gene. Highly significant (P<1×10^-27) smoking-cg05575921 association was validated in 2 additional panels (n=339 and n=300). For the relationship of cg05575921 methylation extent with time after smoking cessation and cumulative cigarette consumption among former smokers, smoking-related hypomethylation was found to remain for ≥20 years after smoking cessation and to be affected by multiple factors, such ascis-interaction of genetic variation. There was a significant inverse correlation (P=0.0005) between cg05575921 methylation extent andAHRRmRNA expression.

CONCLUSIONS

The present study results support that reversion of AHRR hypomethylation can be a quantifiable biomarker for progress in and observance of smoking cessation, although some methodological points need to be considered.

DNA Methylation Imputation Across Platforms

In this chapter, we will provide a review on imputation in the context of DNA methylation, specifically focusing on a penalized functional regression (PFR) method we have previously developed. We will start with a brief review of DNA methylation, genomic and epigenomic contexts where imputation has proven beneficial in practice, and statistical or computational methods proposed for DNA methylation in the recent literature (Subheading 1). The rest of the chapter (Subheadings 2-4) will provide a detailed review of our PFR method proposed for across-platform imputation, which incorporates nonlocal information using a penalized functional regression framework. Subheading 2 introduces commonly employed technologies for DNA methylation measurement and describes the real dataset we have used in the development of our method: the acute myeloid leukemia (AML) dataset from The Cancer Genome Atlas (TCGA) project. Subheading 3 comprehensively reviews our method, encompassing data harmonization prior to model building, the actual building of penalized functional regression model, post-imputation quality filter, and imputation quality assessment. Subheading 4 shows the performance of our method in both simulation and the TCGA AML dataset, demonstrating that our penalized functional regression model is a valuable across-platform imputation tool for DNA methylation data, particularly because of its ability to boost statistical power for subsequent epigenome-wide association study. Finally, Subheading 5 provides future perspectives on imputation for DNA methylation data.

Short- and intermediate-term exposure to ambient fine particulate elements and leukocyte epigenome-wide DNA methylation in older men: the Normative Aging Study

BACKGROUND

Several epigenome-wide association studies (EWAS) of ambient particulate matter with aerodynamic diameter ≤ 2.5 µm (PM_2.5) have been reported. However, EWAS of PM_2.5 elements (PEs), reflecting different emission sources, are very limited.

OBJECTIVES

We performed EWAS of short- and intermediate-term exposure to PM_2.5 and 13 PEs. We hypothesized that significant changes in DNAm may vary by PM_2.5 mass and its elements.

METHODS

We repeatedly collected blood samples in the Normative Aging Study and measured leukocyte DNA methylation (DNAm) with the Illumina HumanMethylation450K BeadChip. We collected daily PM_2.5 and 13 PEs at a fixed central site. To estimate the associations between each PE and DNAm at individual cytosine-phosphate-guanine (CpG) sites, we incorporated a distributed-lag (0-27 d) term in the setting of median regression with subject-specific intercept and examined cumulative lag associations. We also accounted for selection bias due to loss to follow-up and mortality prior to enrollment. Significantly differentially methylated probes (DMPs) were identified using Bonferroni correction for multiple testing. We further conducted regional and pathway analyses to identify significantly differentially methylated regions (DMRs) and pathways.

RESULTS

We included 695 men with 1,266 visits between 1999 and 2013. The subjects had a mean age of 75 years. The significant DMPs, DMRs, and pathways varied by to PM_2.5 total mass and PEs. For example, PM_2.5 total mass was associated with 2,717 DMPs and 10,470 DMRs whereas Pb was associated with 3,173 DMPs and 637 DMRs. The identified pathways by PM_2.5 mass were mostly involved in mood disorders, neuroplasticity, immunity, and inflammation, whereas the pathways associated with motor vehicles (BC, Cu, Pb, and Zn) were related with cardiovascular disease and cancer (e.g., "PPARs signaling").

CONCLUSIONS

PM_2.5 and PE were associated with methylation changes at multiple probes and along multiple pathways, in ways that varied by particle components.

Ionizing radiations induce shared epigenomic signatures unraveling adaptive mechanisms of cancerous cell lines with or without methionine dependency

BACKGROUND

Although radiation therapy represents a core cancer treatment modality, its efficacy is hampered by radioresistance. The effect of ionizing radiations (IRs) is well known regarding their ability to induce genetic alterations; however, their impact on the epigenome landscape in cancer, notably at the CpG dinucleotide resolution, remains to be further deciphered. In addition, no evidence is available regarding the effect of IRs on the DNA methylome profile according to the methionine dependency phenotype, which represents a hallmark of metabolic adaptation in cancer.

METHODS

We used a case-control study design with a fractionated irradiation regimen on four cancerous cell lines representative of HCC (HepG2), melanoma (MeWo and MeWo-LC1, which exhibit opposed methionine dependency phenotypes), and glioblastoma (U251). We performed high-resolution genome-wide DNA methylome profiling using the MethylationEPIC BeadChip on baseline conditions, irradiated cell lines (cumulative dose of 10 Gy), and non-irradiated counterparts. We performed epigenome-wide association studies to assess the effect of IRs and methionine-dependency-oriented analysis by carrying out epigenome-wide conditional logistic regression. We looked for epigenome signatures at the locus and single-probe (CpG dinucleotide) levels and through enrichment analyses of gene ontologies (GO). The EpiMet project was registered under the ID#AAP-BMS_003_211.

RESULTS

EWASs revealed shared GO annotation pathways associated with increased methylation signatures for several biological processes in response to IRs, including blood circulation, plasma membrane-bounded cell projection organization, cell projection organization, multicellular organismal process, developmental process, and animal organ morphogenesis. Epigenome-wide conditional logistic regression analysis on the methionine dependency phenotype highlighted several epigenome signatures related to cell cycle and division and responses to IR and ultraviolet light.

CONCLUSIONS

IRs generated a variation in the methylation level of a high number of CpG probes with shared biological pathways, including those associated with cell cycle and division, responses to IRs, sustained angiogenesis, tissue invasion, and metastasis. These results provide insight on shared adaptive mechanisms of the epigenome in cancerous cell lines in response to IR. Future experiments should focus on the tryptic association between IRs, the initiation of a radioresistance phenotype, and their interaction with methionine dependency as a hallmark of metabolic adaptation in cancer.

Impact of BMI and waist circumference on epigenome-wide DNA methylation and identification of epigenetic biomarkers in blood: an EWAS in multi-ethnic Asian individuals

BACKGROUND

The prevalence of obesity and its related chronic diseases have been increasing especially in Asian countries. Obesity-related genetic variants have been identified, but these explain little of the variation in BMI. Recent studies reported associations between DNA methylation and obesity, mostly in non-Asian populations.

METHODS

We performed an epigenome-wide association study (EWAS) on general adiposity (body mass index, BMI) and abdominal adiposity (waist circumference, WC) in 409 multi-ethnic Asian individuals and replicated BMI and waist-associated DNA methylation CpGs identified in other populations. The cross-lagged panel model and Mendelian randomization were used to assess the temporal relationship between methylation and BMI. The temporal relationship between the identified CpGs and inflammation and metabolic markers was also examined.

RESULTS

EWAS identified 116 DNA methylation CpGs independently associated with BMI and eight independently associated with WC at false discovery rate PFDR < 0.05 in 409 Asian samples. We replicated 110 BMI-associated CpGs previously reported in Europeans and identified six novel BMI-associated CpGs and two novel WC-associated CpGs. We observed high consistency in association direction of effect compared to studies in other populations. Causal relationship analyses indicated that BMI was more likely to be the cause of DNA methylation alteration, rather than the consequence. The causal analyses using BMI-associated methylation risk score also suggested that higher levels of the inflammation marker IL-6 were likely the consequence of methylation change.

CONCLUSION

Our study provides evidence of an association between obesity and DNA methylation in multi-ethnic Asians and suggests that obesity can drive methylation change. The results also suggested possible causal influence that obesity-related methylation changes might have on inflammation and lipoprotein levels.

DNA methylation signatures of incident coronary heart disease: findings from epigenome-wide association studies

Coronary heart disease (CHD) is a type of cardiovascular disease (CVD) that affects the coronary arteries, which provide oxygenated blood to the heart. It is a major cause of mortality worldwide. Various prediction methods have been developed to assess the likelihood of developing CHD, including those based on clinical features and genetic variation. Recent epigenome-wide studies have identified DNA methylation signatures associated with the development of CHD, indicating that DNA methylation may play a role in predicting future CHD. This narrative review summarises recent findings from DNA methylation studies of incident CHD (iCHD) events from epigenome-wide association studies (EWASs). The results suggest that DNA methylation signatures may identify new mechanisms involved in CHD progression and could prove a useful adjunct for the prediction of future CHD.

Ambient temperature and genome-wide DNA methylation: A twin and family study in Australia

Little is known about the association between ambient temperature and DNA methylation, which is a potential biological process through which ambient temperature affects health. This study aimed to evaluate the association between ambient temperature and DNA methylation across human genome. We included 479 Australian women, including 132 twin pairs and 215 sisters of these twins. Blood-derived DNA methylation was measured using the HumanMethylation450 BeadChip array. Data on average ambient temperature during eight different exposure windows [lag0d (the blood draw day), lag0-7d (the current day and previous seven days prior to blood draw), lag0-14d, lag0-21d, lag0-28d, lag0-90d, lag0-180d, and lag0-365d)] was linked to each participant's home address. For each cytosine-guanine dinucleotide (CpG), we evaluated the association between its methylation level and temperature using generalized estimating equations (GEE), adjusting for important covariates. We used comb-p and DMRcate to identify differentially methylated regions (DMRs). We identified 31 CpGs at which blood DNA methylation were significantly associated with ambient temperature with false discovery rate [FDR] < 0.05. There were 82 significant DMRs identified by both comb-p (Sidak p-value < 0.01) and DMRcate (FDR < 0.01). Most of these CpGs and DMRs only showed association with temperature during one specific exposure window. These CpGs and DMRs were mapped to 85 genes. These related genes have been related to many human chronic diseases or phenotypes (e.g., diabetes, arthritis, breast cancer, depression, asthma, body height) in previous studies. The signals of short-term windows (lag0d and lag0-21d) showed enrichment in biological processes related to cell adhesion. In conclusion, short-, medium-, and long-term exposures to ambient temperature were all associated with blood DNA methylation, but the target genomic loci varied by exposure window. These differential methylation signals may serve as potential biomarkers to understand the health impacts of temperature.

Residential surrounding greenness and DNA methylation: An epigenome-wide association study

BACKGROUND

DNA methylation is a potential biological mechanism through which residential greenness affects health, but little is known about its association with greenness and whether the association could be modified by genetic background. We aimed to evaluate the association between surrounding greenness and genome-wide DNA methylation and potential gene-greenness interaction effects on DNA methylation.

METHODS

We measured blood-derived DNA methylation using the HumanMethylation450 BeadChip array (Illumina) for 479 Australian women, including 66 monozygotic, 66 dizygotic twin pairs, and 215 sisters of these twins. Surrounding greenness was represented by Normalized Difference Vegetation Index (NDVI) and Enhanced Vegetation Index (EVI) within 300, 500, 1000 or 2000 m surrounding participants' home addresses. For each cytosine-guanine dinucleotide (CpG), the associations between its methylation level and NDVI or EVI were evaluated by generalized estimating equations, after adjusting for age, education, marital status, area-level socioeconomic status, smoking behavior, cell-type proportions, and familial clustering. We used comb-p and DMRcate to identify significant differentially methylated regions (DMRs). For each significant CpG, we evaluated the interaction effects of greenness and single-nucleotide polymorphisms (SNPs) within ±1 Mb window on its methylation level.

RESULTS

We found associations between surrounding greenness and blood DNA methylation for one CpG (cg04720477, mapped to the promoter region of CNP gene) with false discovery rate [FDR] < 0.05, and for another 9 CpGs with 0.05 ≤ FDR < 0.10. For two of these CpGs, we found 33 SNPs significantly (FDR < 0.05) modified the greenness-methylation association. There were 35 significant DMRs related to surrounding greenness that were identified by both comb-p (Sidak p-value < 0.01) and DMRcate (FDR < 0.01). Those CpGs and DMRs were mapped to genes related to many human diseases, such as mental health disorders and neoplasms as well as nutritional and metabolic diseases.

CONCLUSIONS

Surrounding greenness was associated with blood DNA methylation of many loci across human genome, and this association could be modified by genetic variations.

Detecting cord blood cell type-specific epigenetic associations with gestational diabetes mellitus and early childhood growth

BACKGROUND

Epigenome-wide association studies (EWAS) have provided opportunities to understand the role of epigenetic mechanisms in development and pathophysiology of many chronic diseases. However, an important limitation of conventional EWAS is that profiles of epigenetic variability are often obtained in samples of mixed cell types. Here, we aim to assess whether changes in cord blood DNA methylation (DNAm) associated with gestational diabetes mellitus (GDM) exposure and early childhood growth markers occur in a cell type-specific manner.

RESULTS

We analyzed 275 cord blood samples collected at delivery from a prospective pre-birth cohort with genome-wide DNAm profiled by the Illumina MethylationEPIC array. We estimated proportions of seven common cell types in each sample using a cord blood-specific DNAm reference panel. Leveraging a recently developed approach named CellDMC, we performed cell type-specific EWAS to identify CpG loci significantly associated with GDM, or 3-year-old body mass index (BMI) z-score. A total of 1410 CpG loci displayed significant cell type-specific differences in methylation level between 23 GDM cases and 252 controls with a false discovery rate < 0.05. Gene Ontology enrichment analysis indicated that LDL transportation emerged from CpG specifically identified from B-cells DNAm analyses and the mitogen-activated protein kinase pathway emerged from CpG specifically identified from natural killer cells DNAm analyses. In addition, we identified four and six loci associated with 3-year-old BMI z-score that were specific to CD8+ T-cells and monocytes, respectively. By performing genome-wide permutation tests, we validated that most of our detected signals had low false positive rates.

CONCLUSION

Compared to conventional EWAS adjusting for the effects of cell type heterogeneity, the proposed approach based on cell type-specific EWAS could provide additional biologically meaningful associations between CpG methylation, prenatal maternal GDM or 3-year-old BMI. With careful validation, these findings may provide new insights into the pathogenesis, programming, and consequences of related childhood metabolic dysregulation. Therefore, we propose that cell type-specific analyses are worth cautious explorations.

Pre-diagnostic DNA methylation patterns differ according to mammographic breast density amongst women who subsequently develop breast cancer: a case-only study in the EPIC-Florence cohort

PURPOSE

Mammographic breast density (MBD) is a marker of increased breast cancer (BC) risk, yet much remains to be clarified about the underlying mechanisms. We investigated whether DNA methylation patterns differ between high- vs. low-MBD women who developed BC during an 8.9-year median follow-up in the Florence section of the European Prospective Investigation into Cancer and Nutrition.

METHODS

We analysed 96 pairs of women with BC arising on high- vs. low-MBD breasts (BI-RADS category III-IV vs. I). DNA methylation was determined on pre-diagnostic blood samples using the Illumina Infinium MethylationEPIC BeadChip assay. The statistical analysis was conducted by performing an epigenome-wide association study (EWAS), by searching differentially methylated regions (DMRs) in gene promoters (followed by functional enrichment and gene annotation analysis); and through a "candidate pathways" approach focusing on pre-defined inflammation-related pathways.

RESULTS

In EWAS, no single CpG site was differentially methylated between high- and low-MBD women after correction for multiple testing. A total of 140 DMRs were identified, of which 131 were hyper- and 9 hypo-methylated amongst high-MBD women. These DMRs encompassed an annotation cluster of 35 genes coding for proteins implicated in transcription regulation and DNA binding. The "apoptosis signalling" was the only inflammation-related candidate pathway differentially methylated between high- and low-MBD women.

CONCLUSION

Pre-diagnostic methylation patterns differ between high- vs. low-MBD women who subsequently develop BC, particularly, in genes involved in the regulation of DNA transcription and cell apoptosis. Our study provides novel clues about the mechanisms linking MBD and BC.

Sex-specific DNA methylation differences in Alzheimer's disease pathology

Sex is an important factor that contributes to the clinical and biological heterogeneities in Alzheimer's disease (AD), but the regulatory mechanisms underlying sex disparity in AD are still not well understood. DNA methylation is an important epigenetic modification that regulates gene transcription and is known to be involved in AD. We performed the first large-scale sex-specific meta-analysis of DNA methylation differences in AD neuropathology, by re-analyzing four recent epigenome-wide association studies totaling more than 1000 postmortem prefrontal cortex brain samples using a uniform analytical pipeline. For each cohort, we employed two complementary analytical strategies, a sex-stratified analysis that examined methylation-Braak stage associations in male and female samples separately, and a sex-by-Braak stage interaction analysis that compared the magnitude of these associations between different sexes. Our analysis uncovered 14 novel CpGs, mapped to genes such as TMEM39A and TNXB that are associated with the AD Braak stage in a sex-specific manner. TMEM39A is known to be involved in inflammation, dysregulated type I interferon responses, and other immune processes. TNXB encodes tenascin proteins, which are extracellular matrix glycoproteins demonstrated to modulate synaptic plasticity in the brain. Moreover, for many previously implicated genes in AD neuropathology, such as MBP and AZU1, our analysis provided the new insights that they were predominately driven by effects in only one sex. These sex-specific DNA methylation differences were enriched in divergent biological processes such as integrin activation in females and complement activation in males. Our study implicated multiple new loci and biological processes that affected AD neuropathology in a sex-specific manner.

DNA methylation biomarkers of myocardial infarction and cardiovascular disease

BACKGROUND

The epigenetic landscape underlying cardiovascular disease (CVD) is not completely understood and the clinical value of the identified biomarkers is still limited. We aimed to identify differentially methylated loci associated with acute myocardial infarction (AMI) and assess their validity as predictive and causal biomarkers.

RESULTS

We designed a case-control, two-stage, epigenome-wide association study on AMI (n_discovery = 391, n_validation = 204). DNA methylation was assessed using the Infinium MethylationEPIC BeadChip. We performed a fixed-effects meta-analysis of the two samples. 34 CpGs were associated with AMI. Only 12 of them were available in two independent cohort studies (n ~ 1800 and n ~ 2500) with incident coronary and cardiovascular disease (CHD and CVD, respectively). The Infinium HumanMethylation450 BeadChip was used in those two studies. Four of the 12 CpGs were validated in association with incident CHD: AHRR-mapping cg05575921, PTCD2-mapping cg25769469, intergenic cg21566642 and MPO-mapping cg04988978. We then assessed whether methylation risk scores based on those CpGs improved the predictive capacity of the Framingham risk function, but they did not. Finally, we aimed to study the causality of those associations using a Mendelian randomization approach but only one of the CpGs had a genetic influence and therefore the results were not conclusive.

CONCLUSIONS

We have identified 34 CpGs related to AMI. These loci highlight the relevance of smoking, lipid metabolism, and inflammation in the biological mechanisms related to AMI. Four were additionally associated with incident CHD and CVD but did not provide additional predictive information.

Persistent variations of blood DNA methylation associated with treatment exposures and risk for cardiometabolic outcomes in long-term survivors of childhood cancer in the St. Jude Lifetime Cohort

BACKGROUND

It is well-established that cancer treatment substantially increases the risk of long-term adverse health outcomes among childhood cancer survivors. However, there is limited research on the underlying mechanisms. To elucidate the pathophysiology and a possible causal pathway from treatment exposures to cardiometabolic conditions, we conducted epigenome-wide association studies (EWAS) to identify the DNA methylation (DNAm) sites associated with cancer treatment exposures and examined whether treatment-associated DNAm sites mediate associations between specific treatments and cardiometabolic conditions.

METHODS

We included 2052 survivors (median age 33.7 years) of European ancestry from the St. Jude Lifetime Cohort Study, a retrospective hospital-based study with prospective clinical follow-up. Cumulative doses of chemotherapy and region-specific radiation were abstracted from medical records. Seven cardiometabolic conditions were clinically assessed. DNAm profile was measured using MethylationEPIC BeadChip with blood-derived DNA.

RESULTS

By performing multiple treatment-specific EWAS, we identified 935 5'-cytosine-phosphate-guanine-3' (CpG) sites mapped to 538 genes/regions associated with one or more cancer treatments at the epigenome-wide significance level (p < 9 × 10-8). Among the treatment-associated CpGs, 8 were associated with obesity, 63 with hypercholesterolemia, and 17 with hypertriglyceridemia (false discovery rate-adjusted p < 0.05). We observed substantial mediation by methylation at four independent CpGs (cg06963130, cg21922478, cg22976567, cg07403981) for the association between abdominal field radiotherapy (abdominal-RT) and risk of hypercholesterolemia (70.3%) and by methylation at three CpGs (cg19634849, cg13552692, cg09853238) for the association between abdominal-RT and hypertriglyceridemia (54.6%). In addition, three CpGs (cg26572901, cg12715065, cg21163477) partially mediated the association between brain-RT and obesity with a 32.9% mediation effect, and two CpGs mediated the association between corticosteroids and obesity (cg22351187, 14.2%) and between brain-RT and hypertriglyceridemia (cg13360224, 10.5%). Notably, several mediator CpGs reside in the proximity of well-established dyslipidemia genes: cg21922478 (ITGA1) and cg22976567 (LMNA).

CONCLUSIONS

In childhood cancer survivors, cancer treatment exposures are associated with DNAm patterns present decades following the exposure. Treatment-associated DNAm sites may mediate the causal pathway from specific treatment exposures to certain cardiometabolic conditions, suggesting the utility of DNAm sites as risk predictors and potential mechanistic targets for future intervention studies.

Nonlinear ridge regression improves cell-type-specific differential expression analysis

BACKGROUND

Epigenome-wide association studies (EWAS) and differential gene expression analyses are generally performed on tissue samples, which consist of multiple cell types. Cell-type-specific effects of a trait, such as disease, on the omics expression are of interest but difficult or costly to measure experimentally. By measuring omics data for the bulk tissue, cell type composition of a sample can be inferred statistically. Subsequently, cell-type-specific effects are estimated by linear regression that includes terms representing the interaction between the cell type proportions and the trait. This approach involves two issues, scaling and multicollinearity.

RESULTS

First, although cell composition is analyzed in linear scale, differential methylation/expression is analyzed suitably in the logit/log scale. To simultaneously analyze two scales, we applied nonlinear regression. Second, we show that the interaction terms are highly collinear, which is obstructive to ordinary regression. To cope with the multicollinearity, we applied ridge regularization. In simulated data, nonlinear ridge regression attained well-balanced sensitivity, specificity and precision. Marginal model attained the lowest precision and highest sensitivity and was the only algorithm to detect weak signal in real data.

CONCLUSION

Nonlinear ridge regression performed cell-type-specific association test on bulk omics data with well-balanced performance. The omicwas package for R implements nonlinear ridge regression for cell-type-specific EWAS, differential gene expression and QTL analyses. The software is freely available from https://github.com/fumi-github/omicwas.

DNA methylation mediates the effect of maternal smoking on offspring birthweight: a birth cohort study of multi-ethnic US mother-newborn pairs

Background

Maternal smoking affects more than half a million pregnancies each year in the US and is known to result in fetal growth restriction as measured by lower birthweight and its associated long-term consequences. Maternal smoking also has been linked to altered fetal DNA methylation (DNAm). However, what remains largely unexplored is whether these DNAm alterations are merely markers of smoking exposure or if they also have implications for health outcomes. This study tested the hypothesis that fetal DNAm mediates the effect of maternal smoking on newborn birthweight.

Methods

This study included mother–newborn pairs from a US predominantly urban, low-income multi-ethnic birth cohort. DNAm in cord blood were determined using the Illumina Infinium MethylationEPIC BeadChip. After standard quality control and normalization procedures, an epigenome-wide association study (EWAS) of maternal smoking was performed using linear regression models, controlling for maternal age, education, race, parity, pre-pregnancy body mass index, alcohol consumption, gestational age, maternal pregestational/gestational diabetes, child sex, cord blood cell compositions and batch effects. To quantify the degree to which cord DNAm mediates the smoking-birthweight association, the VanderWeele-Vansteelandt approach for single mediator and structural equational model for multiple mediators were used, adjusting for pertinent covariates.

Results

The study included 954 mother–newborn pairs. Among mothers, 165 (17.3%) ever smoked before or during pregnancy. Newborns with smoking exposure had on average 258 g lower birthweight than newborns without exposure (P < 0.001). Using a false discovery rate (FDR) < 0.05 as the significance cutoff, the EWAS identified 38 differentially methylated CpG sites associated with maternal smoking. Of those, 17 CpG sites were mapped to previously reported genes: GFI1, AHRR, CYP1A1, and CNTNAP2; 8 of those, located in the first three genes, were Bonferroni significantly associated with newborn birthweight and mediated the smoking-birthweight association. The combined mediation effect of the three genes explained 67.8% of the smoking-birthweight association.

Conclusions

Our study not only lends further support that maternal smoking alters fetal DNAm in a multiethnic population, but also suggests that fetal DNAm substantially mediates the maternal smoking–birthweight association. Our findings, if further validated, indicate that DNAm modification is likely an important pathway by which maternal smoking impairs fetal growth and, perhaps, even long-term health outcomes.

An epigenome-wide association study identifies multiple DNA methylation markers of exposure to endocrine disruptors

BACKGROUND

Exposure to environmental endocrine disrupting chemicals (EDCs) may play an important role in the epidemic of metabolic diseases. Epigenetic alterations may functionally link EDCs with gene expression and metabolic traits.

OBJECTIVES

We aimed to evaluate metabolic-related effects of the exposure to endocrine disruptors including five parabens, three bisphenols, and 13 metabolites of nine phthalates as measured in 24-hour urine on epigenome-wide DNA methylation.

METHODS

A blood-based epigenome-wide association study was performed in 622 participants from the Lifelines DEEP cohort using Illumina Infinium HumanMethylation450 methylation data and EDC excretions in 24-hour urine. Out of the 21 EDCs, 13 compounds were detected in >75% of the samples and, together with bisphenol F, were included in these analyses. Furthermore, we explored the putative function of identified methylation markers and their correlations with metabolic traits.

RESULTS

We found 20 differentially methylated cytosine-phosphate-guanines (CpGs) associated with 10 EDCs at suggestive p-value < 1 × 10^-6, of which four, associated with MEHP and MEHHP, were genome-wide significant (Bonferroni-corrected p-value < 1.19 × 10^-7). Nine out of 20 CpGs were significantly associated with at least one of the tested metabolic traits, such as fasting glucose, glycated hemoglobin, blood lipids, and/or blood pressure. 18 out of 20 EDC-associated CpGs were annotated to genes functionally related to metabolic syndrome, hypertension, obesity, type 2 diabetes, insulin resistance and glycemic traits.

CONCLUSIONS

The identified DNA methylation markers for exposure to the most common EDCs provide suggestive mechanism underlying the contributions of EDCs to metabolic health. Follow-up studies are needed to unravel the causality of EDC-induced methylation changes in metabolic alterations.

Cord blood DNA methylome in newborns later diagnosed with autism spectrum disorder reflects early dysregulation of neurodevelopmental and X-linked genes

BACKGROUND

Autism spectrum disorder (ASD) is a neurodevelopmental disorder with complex heritability and higher prevalence in males. The neonatal epigenome has the potential to reflect past interactions between genetic and environmental factors during early development and influence future health outcomes.

METHODS

We performed whole-genome bisulfite sequencing of 152 umbilical cord blood samples from the MARBLES and EARLI high-familial risk prospective cohorts to identify an epigenomic signature of ASD at birth. Samples were split into discovery and replication sets and stratified by sex, and their DNA methylation profiles were tested for differentially methylated regions (DMRs) between ASD and typically developing control cord blood samples. DMRs were mapped to genes and assessed for enrichment in gene function, tissue expression, chromosome location, and overlap with prior ASD studies. DMR coordinates were tested for enrichment in chromatin states and transcription factor binding motifs. Results were compared between discovery and replication sets and between males and females.

RESULTS

We identified DMRs stratified by sex that discriminated ASD from control cord blood samples in discovery and replication sets. At a region level, 7 DMRs in males and 31 DMRs in females replicated across two independent groups of subjects, while 537 DMR genes in males and 1762 DMR genes in females replicated by gene association. These DMR genes were significantly enriched for brain and embryonic expression, X chromosome location, and identification in prior epigenetic studies of ASD in post-mortem brain. In males and females, autosomal ASD DMRs were significantly enriched for promoter and bivalent chromatin states across most cell types, while sex differences were observed for X-linked ASD DMRs. Lastly, these DMRs identified in cord blood were significantly enriched for binding sites of methyl-sensitive transcription factors relevant to fetal brain development.

CONCLUSIONS

At birth, prior to the diagnosis of ASD, a distinct DNA methylation signature was detected in cord blood over regulatory regions and genes relevant to early fetal neurodevelopment. Differential cord methylation in ASD supports the developmental and sex-biased etiology of ASD and provides novel insights for early diagnosis and therapy.

Immediate and durable effects of maternal tobacco consumption alter placental DNA methylation in enhancer and imprinted gene-containing regions

BACKGROUND

Although exposure to cigarette smoking during pregnancy has been associated with alterations of DNA methylation in the cord blood or placental cells, whether such exposure before pregnancy could induce epigenetic alterations in the placenta of former smokers has never been investigated.

METHODS

Our approach combined the analysis of placenta epigenomic (ENCODE) data with newly generated DNA methylation data obtained from 568 pregnant women, the largest cohort to date, either actively smoking during their pregnancy or formerly exposed to tobacco smoking.

RESULTS

This strategy resulted in several major findings. First, among the 203 differentially methylated regions (DMRs) identified by the epigenome-wide association study, 152 showed "reversible" alterations of DNA methylation, only present in the placenta of current smokers, whereas 26 were also found altered in former smokers, whose placenta had not been exposed directly to cigarette smoking. Although the absolute methylation changes were smaller than those observed in other contexts, such as in some congenital diseases, the observed alterations were consistent within each DMR. This observation was further supported by a demethylation of LINE-1 sequences in the placentas of both current (beta-coefficient (β) (95% confidence interval (CI)), - 0.004 (- 0.008; 0.001)) and former smokers (β (95% CI), - 0.006 (- 0.011; - 0.001)) compared to nonsmokers. Second, the 203 DMRs were enriched in epigenetic marks corresponding to enhancer regions, including monomethylation of lysine 4 and acetylation of lysine 27 of histone H3 (respectively H3K4me1 and H3K27ac). Third, smoking-associated DMRs were also found near and/or overlapping 10 imprinted genes containing regions (corresponding to 16 genes), notably including the NNAT, SGCE/PEG10, and H19/MIR675 loci.

CONCLUSIONS

Our results pointing towards genomic regions containing the imprinted genes as well as enhancers as preferential targets suggest mechanisms by which tobacco could directly impact the fetus and future child. The persistence of significant DNA methylation changes in the placenta of former smokers supports the hypothesis of an "epigenetic memory" of exposure to cigarette smoking before pregnancy. This observation not only is conceptually revolutionary, but these results also bring crucial information in terms of public health concerning potential long-term detrimental effects of smoking in women.