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

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

Effect of prenatal arsenic exposure on DNA methylation and leukocyte subpopulations in cord blood

Prenatal arsenic exposure is associated with increased risk of disease in adulthood. This has led to considerable interest in arsenic’s ability to disrupt fetal programming. Many studies report that arsenic exposure alters DNA methylation in whole blood but these studies did not adjust for cell mixture. In this study, we examined the relationship between arsenic in maternal drinking water collected ≤ 16 weeks gestational age and DNA methylation in cord blood (n = 44) adjusting for leukocyte-tagged differentially methylated regions. DNA methylation was quantified using the Infinium HumanMethylation 450 BeadChip array. Recursively partitioned mixture modeling examined the relationship between arsenic and methylation at 473,844 CpG sites. Median arsenic concentration in water was 12 µg/L (range < 1- 510 µg/L). Log₁₀ arsenic was associated with altered DNA methylation across the epigenome (P = 0.002); however, adjusting for leukocyte distributions attenuated this association (P = 0.013). We also observed that arsenic had a strong effect on the distribution of leukocytes in cord blood. In adjusted models, every log₁₀ increase in maternal drinking water arsenic exposure was estimated to increase CD8+ T cells by 7.4% (P = 0.0004) and decrease in CD4+ T cells by 9.2% (P = 0.0002). These results show that prenatal exposure to arsenic had an exposure-dependent effect on specific T cell subpopulations in cord blood and altered DNA methylation in cord blood. Future research is needed to determine if these small changes in DNA methylation alter gene expression or are associated with adverse health effects.     

Regions of variable DNA methylation in human placenta associated with newborn neurobehavior

The placenta regulates the in utero environment and functionally impacts fetal development. Candidate gene studies identified variation in placental DNA methylation is associated with newborn neurologic and behavioral outcomes including movement quality, lethargic behavior, attention, and arousal. We sought to identify novel regions of variable DNA methylation associated with newborn attention, lethargy, quality of movement, and arousal by performing an epigenome-wide association study in 335 infants from a US birth cohort. Methylation status was quantified using the Illumina HumanMethylation450 BeadChip array and associations to newborn outcomes assessed by the NICU Network Neurobehavioral Scales (NNNS) were identified while incorporating established bioinformatics algorithms to control for confounding by cell type composition. Methylation of CpGs within FHIT (cg15970800) and ANKRD11 (cg16710656) demonstrated genome-wide significance (P < 1.8 × 10^?7) in specific associations with infant attention. CpGs whose differential methylation was associated with all 4 neurobehavioral outcomes were common to 50 genes involved in biological processes relating to cellular adhesion and nervous system development. Comprehensive methylation profiling identified relationships between methylation of FHIT and ANKRD11, which have been previously linked to neurodevelopment and behavioral outcomes in genetic association studies. Subtle changes in DNA methylation of these genes within the placenta may impact normal variation of a newborn's ability to alter and track visual and auditory stimuli. Gene ontology analysis suggested that those genes with variable methylation related to these outcomes are over-represented in biological pathways involved in brain development and placental physiology, supportive of our hypothesis for a key role of the placenta in neurobehavioral outcomes.     

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

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

Correlation of DNA methylation patterns to the phenotypic features of Tibetan elite alpinists in extreme hypoxia

High altitude is an extreme environment that imposes hypoxic pressure on physiological processes, and natives living at high altitudes are more adaptive in certain physiological processes. So far, epigenetic modifications under extreme changes in hypoxic pressures are relatively less understood. Here, we recruit 32 Tibetan elite alpinists (TEAs), who have successfully mounted Everest (8848 m) at least five times. Blood samples and physiological phenotypes of TEAs and 32 matched non-alpinist Tibetan volunteers (non-TEAs) are collected for analysis. Genome-wide DNA methylation analysis identifies 23,202 differentially methylated CpGs (P_adj < 0.05, |β| > 0.1) between the two groups. Some differentially methylated CpGs are in hypoxia-related genes such as PPP1R13L, MAP3K7CL, SEPTI-9, and CUL2. In addition, Gene ontology enrichment analysis reveals several inflammation-related pathways. Phenotypic analysis indicates that 12 phenotypes are significantly different between the two groups. In particular, TEAs exhibit higher blood oxygen saturation levels and lower neutrophil count, platelet count, and heart rate. For DNA methylation association analysis, we find that two CpGs (cg16687447, cg06947206) upstream of PTEN were associated with platelet count. In conclusion, extreme hypoxia exposure leads to epigenetic modifications and phenotypic alterations of TEA, providing us clues for exploring the molecular mechanism underlying changes under extreme hypoxia conditions.     

In Utero Exposure to Diethylstilbestrol and Blood DNA Methylation in Women Ages 40–59 Years from the Sister Study

In utero exposure to diethylstilbestrol (DES) has been associated with increased risk of adverse health outcomes such as fertility problems and vaginal as well as breast cancer. Animal studies have linked prenatal DES exposure to lasting DNA methylation changes. We investigated genome-wide DNA methylation and in utero DES exposure in a sample of non-Hispanic white women aged 40–59 years from the Sister Study, a large United States cohort study of women with a family history of breast cancer. Using questionnaire information from women and their mothers, we selected 100 women whose mothers reported taking DES while pregnant and 100 control women whose mothers had not taken DES. DNA methylation in blood was measured at 485,577 CpG sites using the Illumina HumanMethylation450 BeadChip. Associations between CpG methylation and DES exposure status were analyzed using robust linear regression with adjustment for blood cell composition and multiple comparisons. Although four CpGs had p<10⁵, after accounting for multiple comparisons using the false discovery rate (FDR), none reached genome-wide significance. In conclusion, adult women exposed to DES in utero had no evidence of large persistent changes in blood DNA methylation.     

DNA methylation profiles in African American prostate cancer patients in relation to disease progression

This study examined whether differential DNA methylation is associated with clinical features of more aggressive disease at diagnosis and prostate cancer recurrence in African American men, who are more likely to die from prostate cancer than other populations. Tumor tissues from 76 African Americans diagnosed with prostate cancer who had radical prostatectomy as their primary treatment were profiled for epigenome-wide DNA methylation levels. Long-term follow-up identified 19 patients with prostate cancer recurrence. Twenty-three CpGs were differentially methylated (FDR q ≤ 0.25, mean methylation difference ≥ 0.10) in patients with vs. without recurrence, including CpGs in GCK, CDKL2, PRDM13, and ZFR2. Methylation differences were also observed between men with metastatic-lethal prostate cancer vs. no recurrence (five CpGs), regional vs. local pathological stage (two CpGs), and higher vs. lower tumor aggressiveness (one CpG). These results indicate that differentially methylated CpG sites identified in tumor tissues of African American men may contribute to prostate cancer aggressiveness.     

Effect of prenatal arsenic exposure on DNA methylation and leukocyte subpopulations in cord blood

Prenatal arsenic exposure is associated with increased risk of disease in adulthood. This has led to considerable interest in arsenic’s ability to disrupt fetal programming. Many studies report that arsenic exposure alters DNA methylation in whole blood but these studies did not adjust for cell mixture. In this study, we examined the relationship between arsenic in maternal drinking water collected ≤ 16 weeks gestational age and DNA methylation in cord blood (n = 44) adjusting for leukocyte-tagged differentially methylated regions. DNA methylation was quantified using the Infinium HumanMethylation 450 BeadChip array. Recursively partitioned mixture modeling examined the relationship between arsenic and methylation at 473,844 CpG sites. Median arsenic concentration in water was 12 µg/L (range < 1- 510 µg/L). Log₁₀ arsenic was associated with altered DNA methylation across the epigenome (P = 0.002); however, adjusting for leukocyte distributions attenuated this association (P = 0.013). We also observed that arsenic had a strong effect on the distribution of leukocytes in cord blood. In adjusted models, every log₁₀ increase in maternal drinking water arsenic exposure was estimated to increase CD8+ T cells by 7.4% (P = 0.0004) and decrease in CD4+ T cells by 9.2% (P = 0.0002). These results show that prenatal exposure to arsenic had an exposure-dependent effect on specific T cell subpopulations in cord blood and altered DNA methylation in cord blood. Future research is needed to determine if these small changes in DNA methylation alter gene expression or are associated with adverse health effects.     

Epigenome-wide differential DNA methylation between HIV-infected and uninfected individuals

Epigenetic control of human immunodeficiency virus-1 (HIV-1) genes is critical for viral integration and latency. However, epigenetic changes in the HIV-1-infected host genome have not been well characterized. Here, we report the first large-scale epigenome-wide association study of DNA methylation for HIV-1 infection. We recruited HIV-infected (n = 261) and uninfected (n = 117) patients from the Veteran Aging Cohort Study (VACS) and all samples were profiled for 485,521 CpG sites in DNA extracted from the blood. After adjusting for cell type and clinical confounders, we identified 20 epigenome-wide significant CpGs for HIV-1 infection. Importantly, 2 CpGs in the promoter of the NLR family, CARD domain containing gene 5 (NLRC5), a key regulator of major histocompatibility complex class I gene expression, showed significantly lower methylation in HIV-infected subjects than in uninfected subjects (cg07839457: t = ?6.03, P_nominal = 4.96 × 10^?9; cg16411857: t = ?7.63, P_nominal = 3.07 × 10^?13). Hypomethylation of these 2 CpGs was replicated in an independent sample (GSE67705: cg07839457: t = ?4.44, P_nominal = 1.61 × 10^?5; cg16411857: t = ?5.90; P = 1.99 × 10^?8). Methylation of these 2 CpGs in NLRC5 was negatively correlated with viral load in the 2 HIV-infected samples (cg07839457: P = 1.8 × 10^?4; cg16411857: P = 0.03 in the VACS; and cg07839457: P = 0.04; cg164111857: P = 0.01 in GSE53840). Our findings demonstrate that differential DNA methylation is associated with HIV infection and suggest the involvement of a novel host gene, NLRC5, in HIV pathogenesis.     

Prenatal antidepressant exposure associated with CYP2E1 DNA methylation change in neonates

Some but not all neonates are affected by prenatal exposure to serotonin reuptake inhibitor antidepressants (SRI) and maternal mood disturbances. Distinguishing the impact of these 2 exposures is challenging and raises critical questions about whether pharmacological, genetic, or epigenetic factors can explain the spectrum of reported outcomes. Using unbiased DNA methylation array measurements followed by a detailed candidate gene approach, we examined whether prenatal SRI exposure was associated with neonatal DNA methylation changes and whether such changes were associated with differences in birth outcomes. Prenatal SRI exposure was first associated with increased DNA methylation status primarily at CYP2E1(β_Non-exposed = 0.06, β_SRI-exposed = 0.30, FDR = 0); however, this finding could not be distinguished from the potential impact of prenatal maternal depressed mood. Then, using pyrosequencing of CYP2E1 regulatory regions in an expanded cohort, higher DNA methylation status—both the mean across 16 CpG sites (P < 0.01) and at each specific CpG site (P < 0.05)—was associated with exposure to lower 3rd trimester maternal depressed mood symptoms only in the SRI-exposed neonates, indicating a maternal mood x SRI exposure interaction. In addition, higher DNA methylation levels at CpG2 (P = 0.04), CpG9 (P = 0.04) and CpG10 (P = 0.02), in the interrogated CYP2E1 region, were associated with increased birth weight independently of prenatal maternal mood, SRI drug exposure, or gestational age at birth. Prenatal SRI antidepressant exposure and maternal depressed mood were associated with altered neonatal CYP2E1 DNA methylation status, which, in turn, appeared to be associated with birth weight.     

Genetic contribution to variation in DNA methylation at maternal smoking-sensitive loci in exposed neonates

Epigenome-wide DNA methylation association studies have identified highly replicable genomic loci sensitive to maternal smoking during gestation. The role of inter-individual genetic variation in influencing DNA methylation, leading to the possibility of confounding or bias of such associations, has not been assessed. We investigated whether the DNA methylation levels at the top 10 CpG sites previously associated with exposure to maternal smoking during gestation were associated with individual genetic variation at the genome-wide level. Genome-wide association tests between DNA methylation at the top 10 candidate CpG and genome-wide SNPs were performed in 736 case and control participants of the California Childhood Leukemia Study. Three of the strongest maternal-smoking sensitive CpG sites in newborns were significantly associated with SNPs located proximal to each gene: cg18146737 in the GFI1 gene with rs141819830 (P = 8.2×10^?44), cg05575921 in the AHRR gene with rs148405299 (P = 5.3×10^?10), and cg12803068 in the MYO1G gene with rs61087368 (P = 1.3×10^?18). For the GFI1 CpG cg18146737, the underlying genetic variation at rs141819830 confounded the association between maternal smoking and DNA methylation in our data (the regression coefficient changed from ?0.02 [P = 0.139] to ?0.03 [P = 0.015] after including the genotype). Our results suggest that further studies using DNA methylation at cg18146737, cg05575921, or cg12803068 that aim to assess exposure to maternal smoking during gestation should include genotype at the corresponding SNP. New methods are required for adequate and routine inclusion of genotypic influence on DNA methylation in epigenome-wide association studies to control for potential confounding.     

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

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

Genome-wide methylation analysis identifies differentially methylated CpG loci associated with severe obesity in childhood

Childhood obesity is a major public health issue. Here we investigated whether differential DNA methylation was associated with childhood obesity. We studied DNA methylation profiles in whole blood from 78 obese children (mean BMI Z-score: 2.6) and 71 age- and sex-matched controls (mean BMI Z-score: 0.1). DNA samples from obese and control groups were pooled and analyzed using the Infinium HumanMethylation450 BeadChip array. Comparison of the methylation profiles between obese and control subjects revealed 129 differentially methylated CpG (DMCpG) loci associated with 80 unique genes that had a greater than 10% difference in methylation (P-value < 0.05). The top pathways enriched among the DMCpGs included developmental processes, immune system regulation, regulation of cell signaling, and small GTPase-mediated signal transduction. The associations between the methylation of selected DMCpGs with childhood obesity were validated using sodium bisulfite pyrosequencing across loci within the FYN, PIWIL4, and TAOK3 genes in individual subjects. Three CpG loci within FYN were hypermethylated in obese individuals (all P < 0.01), while obesity was associated with lower methylation of CpG loci within PIWIL4 (P = 0.003) and TAOK3 (P = 0.001). After building logistic regression models, we determined that a 1% increase in methylation in TAOK3, multiplicatively decreased the odds of being obese by 0.91 (95% CI: 0.86 – 0.97), and an increase of 1% methylation in FYN CpG3, multiplicatively increased the odds of being obese by 1.03 (95% CI: 0.99 – 1.07). In conclusion, these findings provide evidence that childhood obesity is associated with specific DNA methylation changes in whole blood, which may have utility as biomarkers of obesity risk.     

Whole genome bisulfite sequencing of cell-free DNA and its cellular contributors uncovers placenta hypomethylated domains

BackgroundCirculating cell-free fetal DNA has enabled non-invasive prenatal fetal aneuploidy testing without direct discrimination of the maternal and fetal DNA. Testing may be improved by specifically enriching the sample material for fetal DNA. DNA methylation may allow for such a separation of DNA; however, this depends on knowledge of the methylomes of circulating cell-free DNA and its cellular contributors.ResultsWe perform whole genome bisulfite sequencing on a set of unmatched samples including circulating cell-free DNA from non-pregnant and pregnant female donors and genomic DNA from maternal buffy coat and placenta samples. We find CpG cytosines within longer fragments are more likely to be methylated. Comparison of the methylomes of placenta and non-pregnant circulating cell-free DNA reveal many of the 51,259 identified differentially methylated regions are located in domains exhibiting consistent placenta hypomethylation across millions of consecutive bases. We find these placenta hypomethylated domains are consistently located within regions exhibiting low CpG and gene density. Differentially methylated regions identified when comparing placenta to non-pregnant circulating cell-free DNA are recapitulated in pregnant circulating cell-free DNA, confirming the ability to detect differential methylation in circulating cell-free DNA mixtures.ConclusionsWe generate methylome maps for four sample types at single-base resolution, identify a link between DNA methylation and fragment length in circulating cell-free DNA, identify differentially methylated regions between sample groups, and uncover the presence of megabase-size placenta hypomethylated domains.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-015-0645-x) contains supplementary material, which is available to authorized users.     

Gestational diabetes mellitus epigenetically affects genes predominantly involved in metabolic diseases

Offspring exposed to gestational diabetes mellitus (GDM) have an increased risk for chronic diseases, and one promising mechanism for fetal metabolic programming is epigenetics. Therefore, we postulated that GDM exposure impacts the offspring’s methylome and used an epigenomic approach to explore this hypothesis. Placenta and cord blood samples were obtained from 44 newborns, including 30 exposed to GDM. Women were recruited at first trimester of pregnancy and followed until delivery. GDM was assessed after a 75-g oral glucose tolerance test at 24–28 weeks of pregnancy. DNA methylation was measured at > 485,000 CpG sites (Infinium HumanMethylation450 BeadChips). Ingenuity Pathway Analysis was conducted to identify metabolic pathways epigenetically affected by GDM. Our results showed that 3,271 and 3,758 genes in placenta and cord blood, respectively, were potentially differentially methylated between samples exposed or not to GDM (p-values down to 1 × 10⁻⁰⁶; none reached the genome-wide significance levels), with more than 25% (n = 1,029) being common to both tissues. Mean DNA methylation differences between groups were 5.7 ± 3.2% and 3.4 ± 1.9% for placenta and cord blood, respectively. These genes were likely involved in the metabolic diseases pathway (up to 115 genes (11%), p-values for pathways = 1.9 × 10⁻¹³ < p < 4.0 × 10⁻⁰³; including diabetes mellitus p = 4.3 × 10⁻¹¹). Among the differentially methylated genes, 326 in placenta and 117 in cord blood were also associated with newborn weight. Our results therefore suggest that GDM has epigenetic effects on genes preferentially involved in the metabolic diseases pathway, with consequences on fetal growth and development, and provide supportive evidence that DNA methylation is involved in fetal metabolic programming.     

Circadian gene methylation profiles are associated with obesity,metabolic disturbances and carbohydrate intake

The circadian clock regulates the daily rhythms of several physiological and behavioral processes. Disruptions in clock genes have been associated with obesity and related comorbidities. This study aimed to analyze the association of DNA methylation signatures at circadian rhythm pathway genes with body mass index (BMI), metabolic profiles and dietary intakes. DNA methylation profiling was determined by microarray in white blood cells from 474 adults from the Methyl Epigenome Network Association (MENA) project. Kyoto Encyclopedia of Genes and Genomes database was used to identify the genes integrating the circadian rhythm pathway. Network enrichment analyses were performed with the PathDIP platform. Associations between circadian methylation patterns with anthropometric measurements, the metabolic profile, clinical data and dietary intakes were analyzed. DNA methylation patterns of nine CpG sites at six circadian rhythm pathway genes were strongly correlated with BMI (false discovery rates <0.0001). These CpGs encompassed cg09578018 (RORA), cg20406576 (PRKAG2), cg10059324 (PER3), cg01180628 (BHLHE40), cg23871860 (FBXL3), cg16964728 (RORA), cg14129040 (CREB1), cg07012178 (PRKAG2) and cg24061580 (PRKAG2). Interestingly, network enrichment analyses revealed that the six BMI-associated genes statistically contributed to the regulation of the circadian rhythm pathway (p = 1.9E-10). In addition, methylation signatures at cg09578018 (RORA), cg24061580 (PRKAG2), cg01180628 (BHLHE40) and cg10059324 (PER3) also correlated with insulin resistance (p < 0.0001) and mean arterial blood pressure (p < 0.0001). Furthermore, relevant correlations (p < 0.05) between methylation at cg09578018 (RORA) and cg01180628 (BHLHE40) with total energy and carbohydrate intakes were found. This investigation revealed potential associations of DNA methylation profiles at circadian genes with obesity, metabolic disturbances and carbohydrate intake, with potential impact on weight homeostasis.