- Research
- Open access
- Published:
Gestational exposure to environmental chemicals and epigenetic alterations in the placenta and cord blood mononuclear cells
Epigenetics Communications volume 4, Article number: 4 (2024)
Abstract
Background
Exposure to environmental chemicals such as phthalates, phenols, and polycyclic aromatic hydrocarbons (PAHs) during pregnancy can increase the risk of adverse newborn outcomes. We explored the associations between maternal exposure to select environmental chemicals and DNA methylation in cord blood mononuclear cells (CBMC) and placental tissue (maternal and fetal sides) to identify potential mechanisms underlying these associations.
Method
This study included 75 pregnant individuals who planned to give birth at the University of Cincinnati Hospital between 2014 and 2017. Maternal urine samples during the delivery visit were collected and analyzed for 37 biomarkers of phenols (12), phthalates (13), phthalate replacements (4), and PAHs (8). Cord blood and placenta tissue (maternal and fetal sides) were also collected to measure the DNA methylation intensities using the Infinium HumanMethylation450K BeadChip. We used linear regression, adjusting for potential confounders, to assess CpG-specific methylation changes in CBMC (n = 54) and placenta [fetal (n = 67) and maternal (n = 68) sides] associated with gestational chemical exposures (29 of 37 biomarkers measured in this study). To account for multiple testing, we used a false discovery rate q-values < 0.05 and presented results by limiting results with a genomic inflation factor of 1±0.5. Additionally, gene set enrichment analysis was conducted using the Kyoto Encyclopedia of Genes and Genomics pathways.
Results
Among the 29 chemical biomarkers assessed for differential methylation, maternal concentrations of PAH metabolites (1-hydroxynaphthalene, 2-hydroxyfluorene, 4-hydroxyphenanthrene, 1-hydroxypyrene), monocarboxyisononyl phthalate, mono-3-carboxypropyl phthalate, and bisphenol A were associated with altered methylation in placenta (maternal or fetal side). Among exposure biomarkers associated with epigenetic changes, 1-hydroxynaphthalene, and mono-3-carboxypropyl phthalate were consistently associated with differential CpG methylation in the placenta. Gene enrichment analysis indicated that maternal 1-hydroxynaphthalene was associated with lipid metabolism and cellular processes of the placenta. Additionally, mono-3-carboxypropyl phthalate was associated with organismal systems and genetic information processing of the placenta.
Conclusion
Among the 29 chemical biomarkers assessed during delivery, 1-hydroxynaphthalene and mono-3-carboxypropyl phthalate were associated with DNA methylation in the placenta.
Background
Pregnant individuals are exposed to a spectrum of environmental chemicals that may increase their risk of adverse maternal or fetal health outcomes [1]. Among these chemicals, metabolites of polycyclic aromatic hydrocarbons (PAHs), phenols, and phthalates are frequently detected in the U.S. pregnant population [1]. Despite the non-persistent nature of these chemicals, exposures during pregnancy have been associated with adverse birth outcomes [2,3,4,5,6,7,8]. Exposure to exogenous endocrine-disrupting chemicals (EDC) such as phenols, phthalates, and PAHs (present in the environment, food, and consumer products) may result in homeostatic, reproductive, and developmental dysregulation [9,10,11].
Environmental chemical exposures during pregnancy have the potential to influence up to three generations that encompass the mother, developing child, and developing gametes within the developing embryo/fetus [12,13,14]. Exposure to EDCs during fetal development may result in adverse health outcomes during adulthood, and these environmental insults can be observed from the alterations in DNA methylation [12]. PAHs, phthalates, and phenols, exhibit endocrine-disrupting properties by mimicking or altering the synthesis/metabolism/transport/action/excretion of hormones (thyroid, insulin, estrogen, androgens), which in turn could disrupt hormonal physiology related to epigenetic pathways [9, 15, 16]. Several studies evaluated the associations between gestational exposure to phthalates, phenols, and cord blood global methylation, suggesting inconsistent patterns in CpG-specific associations [17,18,19,20,21,22].
However, a limited number of studies have examined placenta methylation changes in relation to gestational exposure to environmental chemicals [23]. A few of these studies highlighted associations between gestational phenol/phthalate biomarkers and alterations in the global DNA methylations(DNAm), which are limited to the fetal side of the placenta [22]. Given the importance of the placenta, which plays a critical role in nutrient regulation, gas/waste exchange, and hormone secretion, and the relevance of placental epigenetic alterations to newborn phenotyping, this study focuses on cord blood mononuclear cells (CBMC) and the placenta methylation [23,24,25,26]. In the current study, we quantified the global DNAm of CBMC and maternal/fetal side of placenta changes associated with non-persistent EDC biomarkers to highlight potential developmental origins of diseases.
Materials and methods
Study participants
We enrolled 75 pregnant individuals between August 2014 and September 2017 from a labor and delivery unit at the University of Cincinnati Medical Center, Cincinnati, Ohio. The Institutional Review Board at the University of Cincinnati reviewed and approved the research protocol. The study participants were women with singleton pregnancies, aged between 18 and 45 years, who were admitted for childbirth and provided written consent for their involvement in the study at enrollment. Individuals with a medical history involving diabetes, thyroid disorders, cardiovascular, renal, hepatic conditions, cancers affecting pregnancy, or severe maternal or fetal complications at the time of recruitment were excluded from participation.
We collected 30 mL of maternal urine during the hospital visit for childbirth. Additionally, 10 mL of cord blood and placenta tissue from the maternal and fetal sides were collected. Detailed information on extracting the maternal and fetal side of the placenta tissue is available in the supplement information (SM1). The biospecimens were stored at -80 °C and shipped to labs on dry ice overnight; biospecimens were kept frozen until analysis. Demographic, medical, tobacco use during pregnancy, pre-pregnancy body mass index (BMI), and pregnancy outcome data for both the mother and infant were obtained through questionnaires and a review of medical records.
Chemical exposure assessment
The urinary concentrations of eight PAH metabolites (1-hydroxynaphthalene, 2-hydroxynaphthalene, 2-hydroxyfluorene, 1-hydroxyphenanthrene, 2,3-hydroxyphenanthrene, 4-hydroxyphenanthrene, 9-hydroxyphenanthrene, and 1-hydroxypyrene) were quantified at NSF International – Applied Research Center in Ann Arbor, MI using ultra-high-performance liquid chromatography tandem mass spectrometry [27]. The assays were implemented using a Thermo Scientific Transcend TXII Turbulent Flow system coupled with a Thermo Scientific Vantage triple quadrupole mass spectrometer using multiple reaction monitoring in negative mode.
We also quantified urinary concentrations of 29 biomarkers of triclocarban, phenols, parabens, phthalates, di(2-ethylhexyl) terephthalate, and di(isononyl) cyclohexane-1,2-dicarboxylate using on-line solid-phase extraction coupled with high-performance liquid chromatography isotope dilution tandem mass spectrometry [28, 29]. Concentrations of these biomarkers were quantified at the Centers for Disease Control and Prevention (CDC), Atlanta, GA. The involvement of the CDC laboratory did not constitute engagement in human-participant research.
The limit of detection (LOD) for each chemical biomarker is listed in Table 1. Biomarker concentrations below LOD were imputed with LOD divided by the square root of two [30]. Among the 37 chemical biomarkers quantified from the maternal urine, we excluded eight biomarkers. Specifically, we excluded four phthalate/phthalate replacement biomarkers (mono-2-ethyl-5-carboxypentyl phthalate, mono-2-ethyl-5-hydroxyhexyl phthalate, mono-2-ethyl-5-oxohexyl phthalate, and mono-2-ethylhexyl phthalate) and three parabens (ethyl paraben, methyl paraben, propyl paraben) detected in all study participants due to their potential exposure from medical products during labor and delivery. Additionally, we excluded mono-isononyl phthalate, which was detected in 6% of participants, as outlined in Table 1.
DNA methylation, preprocessing, and cell deconvolution
We extracted DNA using a DNeasy Blood and Tissue kit with RNase A, following the manufacturer’s protocols. The DNA was then prepared for the array following the Illumina Infinium HD Assay Methylation Protocol Guide. We randomly assigned DNA samples to well plates within a chip, as visualized in Figure S1, to reduce potential batch effects. Bisulfite conversion of the DNA was performed using the EZ DNA Methylation Kit (Zymo Research, Irvine, CA), and DNA methylation levels were measured using Infinium HumanMethylation450K BeadChip (Illumina, San Diego, CA) by the Genomics, Epigenomics and Sequencing Core at the University of Cincinnati.
Raw methylation data were processed following the steps recommended by Wilhelm-Benartzi et al. using the minfi Bioconductor R-package [31, 32]. Quality control steps include removing samples with > 5% high detection p-values (p-value > 1e-7; n = 0) and probes that failed in one or more samples using the same p-value cutoff [33]. Background correction with dye-bias normalization was performed using the normal-exponential out-of-band (noob) procedure [34]. Additionally, bias resulting from type-2 probe values was normalized using the beta-mixture quantile normalization method [35] and followed by potential batch effect correction using an empirical Bayes framework available through the ComBat function [36]. Additionally, the CpG probes located on sex chromosomes, enriched at single-nucleotide polymorphisms (SNPs), and those prone to cross-reactivity were excluded [37]. Following these preprocessing steps, we included 418,997 CpGs from CBMC, 415,604 from the fetal side of the placenta, and 412,460 from the maternal side of the placenta for further downstream analysis. We estimated the proportions of five cell types [T-lymphocytes (CD8T, CD4T), NK cells (CD56+), B-cells (CD19+), and monocytes (CD14+)] using CBMC methylation intensities using Houseman’s reference-based (cord blood – HM450) method [38, 39]. Additionally, we estimated the proportions of six cell types [trophoblasts, syncytiotrophoblast, stromal, Hofbauer, endothelial, and nucleated red blood cells (nRBC)] based on placenta methylation intensities using a reference-based (third-trimester placenta) approach [40, 41].
These CpG sites were annotated to genes using the University of California, Santa Cruz Genome Browser (Human GRCh37/hg19) database [42, 43]. In subsequent analysis, we used M-values (log2 ratio of the intensities of methylated versus unmethylated probe) to enhance statistical validity [44].
Statistical analysis
Differential methylated CpG positions
We identified differentially methylated CpGs using the 22 chemical biomarkers (urine creatinine adjusted and log10-transformed) measured in maternal urine as independent variables and the M-values (methylation intensity) of CpGs as dependent variables. Additionally, the associations between 7 maternal chemical biomarkers and the methylation intensities were assessed, considering chemical biomarkers as binary variables (detected or not-detected) due to a relatively lower detection frequency (detected among 20–40% of participants) among the study participants. These associations were adjusted for covariates based on prior literature, including maternal age, race, education, tobacco use during pregnancy, BMI, child sex, and cell type proportions [45, 46]. Associations between 22 biomarker concentrations measured on a continuous scale and M-value per CpG site were assessed using multiple linear regression [47]. Additionally, the association between seven biomarkers considered as binary variables and M-values per CpG site were tested using the ANOVA test [47]. Among the 75 subjects in this study, we limited our sample size to 54 while analyzing CBMC samples, 67 fetal side placenta, and 68 maternal side placenta, due to the availability of DNA methylation (DNAm) and chemical biomarker data. To handle the false positives as a result of multiple tests, we considered CpG sites with Benjamini-Hochberg’s false discovery rate (FDR) at the 5% level (q-values < 0.05) as statistically significant [48]. To expand on biological phenotype interpretations, we performed downstream analysis using CpGs with FDR q-value < 0.05 and for chemical biomarkers with genomic inflation factor (λ) 1 ± 0.5, to minimize false positives [49]. The inflation factor is the ratio between the median of the observed and expected chi-squared statistic [50]. Although the λ value close to one reflects no evidence of inflation (false positives), we considered a liberal threshold of 1 ± 0.5, due to the exploratory study objective [51].
Gene enrichment, pathways, and phenotype inference
Gene enrichment analysis was performed assuming Wallenius noncentral hypergeometric distribution, using the missMethyl R-package considering the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways [50]. Additionally, we performed a sensitivity analysis using the 500 CpGs with relatively lowest q-values for gene enrichment and highlighted the pathways with q-values < 0.1.
The analyses and preprocessing throughout this article were performed using R-packages: minfi 1.44, sva 3.46, wateRmelon 2.4, planet 1.8, CpGassoc 2.6, missMethyl 1.32, R 4.2 [52].
Associations between maternal biomarkers and estimated cell composition
We additionally explored potential associations between maternal chemical biomarkers and cell compositions estimated using DNAm intensities from CBMC/placenta. We tested for normality distribution assumption for estimated cell proportions using the Shapiro-Wilk test [53]. These associations were assessed using a multiple linear regression approach adjusting for maternal age, race, education, tobacco use during pregnancy, BMI, and child sex. The beta coefficients with a corresponding 95% confidence interval without a null value were considered statistically significant. Additionally, we performed a sensitivity analysis by transforming both the maternal urinary chemical biomarker concentrations and cell proportions to log10 scale to address potential non-linear trends of the outcome variables.
Results
This study involved pregnant individuals with a median age of 29 years (interquartile range (IQR): 25–32) at the time they gave birth, and the majority (83.3%) had previous pregnancies (Table S1). About half of the subjects were non-Hispanic Black females (47.2%). The median BMI of the study participants was 25.9 kg/m2 (IQR: 22.5–30.5), and 17% self-reported smoking tobacco during their pregnancy. Most study participants (75%) attained less than a 4-year college education, and 64% had a household income below the median for Cincinnati, OH. All the participants included in this study had full-term pregnancies, with gestational ages ranging from 37 to 41 weeks. We observed correlations between maternal urinary biomarkers that varied between − 0.29 and 0.78 (Figure S2). Overall, the biomarkers within a chemical class showed relatively strong positive correlations.
Differential methylation of CpG sites
Among the epigenetic-wide associations (EWAS) assessed with 22 chemical biomarkers, we observed statistically significant associations with differential methylation of CpG sites for 7 biomarkers (Table S2). These biomarkers encompass PAHs (1-hydroxynaphthalene, 2-hydroxynaphthalene, 2-hydroxyfluorene, 2,3-hydroxyphenanthrene, 4-hydroxyphenanthrene, 1-hydroxypyrene); phenols (bisphenol A); and phthalates (mono-3-carboxypropyl phthalate [MCPP]), mono-2-ethyl-5-hydroxyhexyl terephthalate, mono-2-ethyl-5-carboxypentyl terephthalate, mono carboxy-isononyl phthalate). Notably, 2 phthalate metabolites (mono-2-ethyl-5-carboxypentyl terephthalate, and mono-2-ethyl-5-hydroxyhexyl terephthalate) associated with maternal side CpG methylation; 2,3-hydroxyphenanthrene association with fetal side placenta methylation; and 2-hydroxynaphthalene, MCPP associations with CBMC methylation were not highlighted in this study due to higher lambda values (genomic inflation factor).
Among the 22 biomarkers assessed using the linear regression approach, we observed a relatively higher number (n = 46) of CpGs from the fetal side of the placenta associated with 2-hydroxyfluorene concentrations (Fig. 1). Among these associations, most (11 of 46) were related to hypomethylation of CpG sites that are within 200 base pairs of the promotor [TSS (transcription start site)] region of genes—followed by 1-hydroxypyrene associated with 31 CpGs from the fetal side of the placenta. Nine of the 31 CpG sites associated with 1-hydroxypyrene were hypomethylated and within 200 base pairs of the TSS region of a gene. There were 23 CpG sites from the fetal side of the placenta that were commonly associated with 1-hydroxypyrene and 2-hydroxyfluorene (Figs. 2 and 3). Similarly, a higher number of CpG sites (n = 10) from the maternal side of the placenta were commonly associated with 1-hydroxypyrene and 2-hydroxyfluorene.
Summary of CpG-specific associations (q-value < 0.05) with gestational chemical biomarkers. Horizontal bars represent the total number of CpGs associated with chemical biomarkers by methylation sample source. MP-maternal side placenta, FP-fetal side placenta. Vertical bars represent the number of CpG associations that overlap across and are exclusive to a chemical biomarker
Association between PAH biomarkers concentrations and CpG methylation intensities. This figure covers 1-hydroxynaphthalene and 2-hydroxyfluorene. Methylation intensities from fetal-side placenta (FP) and maternal-side placenta (MP). The x-axis represents the beta coefficient values, and the y-axis represents raw p-values on the -log10 scale. The CpG ID, chromosome position, and UCSC Reference Gene names were labeled. CpG-specific associations with q-values < 0.05 were color-coded in red (hypomethylation) and blue (hypermethylation)
Gestational urinary concentrations of 2-hydroxyfluorene were associated with hypomethylation of seven CpG sites from the fetal side of the placenta, which are located at the promotor region of genes zinc finger protein 345 C (ZNF354C) and coiled-coil domain containing 63 (CCDC63). The CpG sites closest to the ZNF354C gene were in the TSS1500 region and were labeled as cg19761115, cg10864596, and cg12623648. For the CCDC63 gene, the closest CpG sites were in the TSS200 region and were labeled as cg05924445, cg10995082, cg01805574, and cg02855409.
Furthermore, we observed associations between 2-hydroxyfluorene and hypomethylation of 4 of 16 CpG sites (cg03434149, cg16796590, cg13181052, cg10390736) from the maternal side of the placenta that are within 200 base pairs of the TSS region of major facilitator superfamily domain containing 2 A (MFSD2A) gene. Additionally, there were associations between 1-hydroxypyrene and hypomethylation of 5 of 16 CpG sites (cg19694978, cg03434149, cg16796590, cg13181052, cg10390736) located within 200 base pairs of the TSS region of MFSD2A gene.
Highlighting seven biomarkers assessed using the ANOVA test, we observed that MCPP was associated with differential methylation of the placenta. About 130 CpGs from the maternal side of the placenta and 69 CpGs from the fetal site of the placenta were associated with maternal MCPP concentrations. Among these CpGs associated with MCPP, 34 from the maternal side of the placenta and 15 from the fetal side of the placenta were located in the TSS regions. Here, we highlighted the results with more than a CpG located in the TSS region and associated with MCPP concentrations. Among the CpGs associated with MCPP on the maternal side of the placenta, four specific CpGs (cg05180121, cg17114573, cg08216792, cg22492943) were identified within 200 base pairs of the TSS region of the DNA Ligase 4 (LIG4) gene. On the fetal side of the placenta, two CpGs (cg08231348, cg11792516) associated with MCPP were located within 200 base pairs of the TSS region of the hook microtubule tethering protein (HOOK1) gene, and two CpGs (cg02788400, cg05505803) located within 200 base pairs of the TSS region of the DAZ (deleted in azoospermia) interacting zinc finger protein (DZIP1) gene.
Enrichment of Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways
None of the KEGG pathways enriched with FDR q-values below 0.1. Therefore, we highlighted the pathways based on raw p-values less than 0.01. Our findings from the enrichment analysis suggest that 1-hydroxynaphthalene urinary concentrations were associated with 6 KEGG pathways using CpG methylation intensities from the maternal and fetal side of the placenta (Table 2). Of these, there were four pathways enriched for CpGs from the fetal side of the placenta that represent lipid metabolism (fatty acid biosynthesis, degradation, and metabolism) and cellular process (ferroptosis). At the same time, two pathways were enriched while using CpGs from the maternal side of the placenta that represent viral infectious disease (Herpes simplex virus-1) and genetic information replication or repair (non-homologous end-joining). In contrast, results from our sensitivity analysis using the top 500 affected CpG sites suggest associations between CpGs from the maternal side placenta and 1-hydroxynaphthalene were through fatty acid metabolism pathway (Table S3).
There were seven CpG sites (cg25448281, cg11768167, cg25307371, cg01969938, cg25044186, cg05902997, cg14056357) from the maternal side of the placenta that were associated with 1-hydroxynaphthalene. Most of these CpGs (6 of 7) associated with 1-hydroxynaphthalene were hypomethylated. These hypomethylated CpGs are located closer to the promotor region of LEPROPT/LEPR (cg25307371), LIG4 (cg01969938), SFRS12 (cg05902997), and ZNF229 (cg14056357) genes, except for cg25448281, which is located at the end of a protein-coding sequence (3’UTR) of FARP1 gene.
There were relatively lower (n = 2; cg20707765, cg25307371) CpG site-specific associations from the fetal side of the placenta that correspond to leptin receptor overlapping transcript (LEPROT/LEPR) and acyl-CoA synthetase long-chain family member 5 (ACSL5), genes. Urinary concentrations of 1-hydroxynaphthalene were associated with hypermethylation of cg20707765, located within 1500 base pairs of the TSS region of gene ACSL5. At the same time, 1-hydroxynaphthalene was associated with hypomethylation of cg25307371, located within 1500 base pairs of gene LEPROT/LEPR.
Furthermore, CpGs from the maternal side of the placenta, associated with MCPP concentrations, were associated with genetic information processing (non-homologous end-joining [NHEJ]) and organismal systems (neutrophil extracellular trap formation). The NHEJ could be attributable to the four CpGs associated with MCPP from the maternal side of the placenta, positioned proximally to the TSS of the LIG4 gene. Additionally, enrichment of the neutrophil extracellular trap formation could be related to the CpG (cg23207054) from the maternal side of the placenta that is located within 200 base pairs of the colony-stimulating factor 3 (CSF3) gene and associated with the MCPP concentrations.
Associations between maternal chemical biomarkers and estimated cell composition
Most of these cell proportions estimated from the DNA methylation intensities showed normal distributions, except for NK-cells from CBMC, Hofbauer, and nRBC cells from the fetal side of the placenta that are slightly skewed (Figure S3). We observed a weak association pattern between chemical biomarkers and estimated cell compositions in CBMC/placenta (Figure S4-A). Among the cells estimated using DNAm intensities in CBMC, we observed statistically significant negative associations between biomarkers of PAHs and B-cells. Similarly, we observed statistically significant negative associations between biomarkers of PAHs and nRBC cells estimated from the DNAm intensities measured from the maternal side of the placenta. However, we also observed statistically significant positive associations between phthalate biomarkers and CD4T cells from CBMC as well as Trophoblasts from the maternal side of the placenta. Moreover, transforming the cell proportions on a log10 scale, we identified additional negative associations between 2-hydroxyfluorene, 2,3-hydroxyphenanthrene, monobenzyl phthalate, and B-cell proportions estimated from CBMC; 1-hydroxyphenanthrene and nRBCs from maternal placenta; monobenzyl phthalate and stromal cell proportions estimated from maternal placenta (Figure S4-B). The results presented here for NK cells from cord blood, nRBC from the fetal placenta, and Hofbauer cells from the maternal/fetal placenta require caution in interpretation since the cell proportions were not normally distributed (Table S4).
Discussion and conclusions
Among the 29 environmental chemical biomarkers assessed for EWAS, 11 biomarkers (6 PAHs, 4 phthalates, and one phenol) were associated with differential CpG methylation levels. Across the sample media (CBMC, fetal side placenta, and maternal side placenta), the differential methylation associations were prominent with the fetal side of the placenta. Our gene enrichment analysis results indicate associations between 1-hydroxynaphthalene, MCPP, and placental metabolism pathways. Our results highlight plausible mechanisms linking exposure to environmental chemicals during pregnancy and placenta development.
Focusing on the 1-hydroxynaphthalene and CpG methylation while using the fetal side of the placenta, two CpG sites are closest to LEPROT/LEPR (cg25307371), and ACSL5 (cg20707765) genes that were differentially methylated. The LEPROT/LEPR gene is associated with obesity and ACSL5 – obesity and gastrointestinal conditions [54, 55]. There are studies focused on the associations between PAH metabolites and DNA methylation, suggesting altered global methylation patterns, epigenetic aging, chronic diseases (pulmonary, cardiovascular), endocrine disruption, nervous system disorders, and cancer [56,57,58,59,60]. Additionally, gestational measurements of the PAH metabolites were associated with methylation of ACSL3 from cord white blood cells, which was found to be a potential biomarker for environmentally-related asthma [61].
Summarizing our findings from the gene enrichment analysis, we observed associations between differentially methylated CpGs from the fetal side of the placenta related to 1-hydroxynaphthalene and lipid metabolism and cellular process. Alterations in lipid metabolism pathways may result in adverse metabolic conditions. Highlighting on the MCPP and CpG methylation from the maternal side of the placenta, there were 34 differentially methylated CpGs in the TSS region representing 31 distinct genes. Furthermore, our findings from the enrichment analysis suggested associations between CpGs from the maternal side of the placenta associated with MCPP and genetic information processing (NHEJ) and organismal systems (neutrophil extracellular trap formation) pathways. The enrichment of these pathways could be relevant to the differential methylation of four CpGs that are located near the TSS region of the LIG4 gene and a CpG that is located near the TSS region of the CSF3 gene. Previous studies have implicated that the LIG4 gene is associated with phenotypic conditions such as metatarsus adductus, abnormality of chromosome stability, abnormality in bone marrow cell morphology, acute lymphoblastic leukemia, and thrombocytopenia [54, 55]. Additionally, the CSF3 gene is associated with the production and differentiation of granulocytes and macrophages [54, 55].
In contrast to our primary findings, which revealed associations between 1-hydroxynaphthalene and changes in lipid metabolism with the fetal-side placenta, our sensitivity analysis, using 500 CpGs with the lowest p-values, indicated associations between 1-hydroxynaphthalene and altered fatty acid metabolism with the maternal-side placenta. Overall, we observed associations between 1-hydroxynaphthalene, mono-3-carboxypropyl phthalate and differential placenta methylation, which could lead to alterations in fatty acid metabolism, genetic information processing and organismal systems. Similar results have been presented in the literature, suggesting associations between PAH metabolites and metabolic syndrome, cardiovascular conditions, and systemic inflammation [62, 63]. Additionally, a mice-based study suggested prenatal exposure to PAHs and their associated multi-generational (offspring and grand-offspring) epigenetic effects [64]. Furthermore, Jednak et al. identified associations between maternal MCPP exposures and placenta methylation, identifying differential methylation regions of the fibroblast growth factor 12 (FGF12) (chr3:192445514–192445539), and the ADP-ribosyltransferase 5 (ART5) (chr11:3663491–3663843) genes, suggesting potential alterations in protein function [65]. However, our study’s findings did not align with those reported by Jednak et al., as they measured phthalate biomarkers at ~ 22–29 weeks of pregnancy, contrasting with our study’s measurements at delivery.
Most of the existing literature is supported by the use of whole blood or cord blood to identify global methylation changes associated with organic pollutants, bisphenol A, and polybrominated diphenyl ethers [23]. Our study contributes by exploring potential early-life changes in DNA methylation associated with a spectrum of environmental exposures, using placental tissue. However, we did not find consistent patterns while comparing our results with the existing literature. These inconsistencies could relate to heterogeneity in the sample media to estimate methylation intensities. Although we attempted to minimize false discoveries using the q-value thresholds while identifying differentially methylated CpGs, false positives may still influence our results. Our results are limited to the discovery set analysis and would benefit from a replication study to validate our findings. Additionally, we focused on non-persistent chemical biomarkers measured at a single point in time during the delivery visit, which may not adequately represent exposure during pregnancy and may introduce exposure misclassification bias. The relatively strong correlation pattern between maternal urinary chemical biomarkers may introduce bias due to multi-collinearity.
Furthermore, our findings could be influenced by bias due to unmeasured confounding. Another limitation of our study is the use of the HM450 array that contains substantially fewer type-II probes than compared to the Infinium MethylationEPICv1/v2 chip, which may lead to an underrepresentation of our underlying findings [66, 67]. However, this limitation could be to our advantage to minimize false positives from multiple testing given with our study sample size.
In summary, we observed associations between gestational urinary concentrations of select PAH, phthalate, and phenol biomarkers and differential methylation patterns in the placenta (maternal or fetal). Our findings from the gene enrichment analysis suggest that gestational exposure to some PAH and phthalate metabolites may influence newborn health. However, to confirm these findings, further replication studies are necessary. Additionally, conducting comprehensive multi-omics studies, integrating exposome, epigenetics, transcriptomics, and phenotype data in a larger sample, could enhance our understanding of the biological pathways that underlie the associations between exposure and phenotype.
Data availability
The HumanMethylation450K BeadChip data discussed in this publication have been deposited in NCBI’s Gene Expression Omnibus and are accessible through GEO Series accession number GSE269983.
References
Woodruff TJ, Zota AR, Schwartz JM. Environmental Chemicals in pregnant women in the United States: NHANES 2003–2004. Environ Health Perspect. 2011;119(6):878–85.
Wolff MS, Engel SM, Berkowitz GS, Ye X, Silva MJ, Zhu C, et al. Prenatal phenol and phthalate exposures and birth outcomes. Environ Health Perspect. 2008;116(8):1092–7.
Engel SM, Zhu C, Berkowitz GS, Calafat AM, Silva MJ, Miodovnik A, Wolff MS. Prenatal phthalate exposure and performance on the neonatal behavioral Assessment Scale in a multiethnic birth cohort. Neurotoxicology. 2009;30(4):522–8.
Nidens N, Krönke A, Jurkutat A, Schlingmann M, Poulain T, Nüchter M, et al. Associations of prenatal exposure to phthalates and one phthalate substitute with anthropometric measures in early life: results from the German LIFE child cohort study. Best Pract Res Clin Endocrinol Metab. 2021;35(5):101532.
Perera FP, Whyatt RM, Jedrychowski W, Rauh V, Manchester D, Santella RM, Ottman R. Recent developments in Molecular Epidemiology: a study of the effects of Environmental Polycyclic Aromatic Hydrocarbons on Birth outcomes in Poland. Am J Epidemiol. 1998;147(3):309–14.
Jedrychowski WA, Majewska R, Spengler JD, Camann D, Roen EL, Perera FP. Prenatal exposure to fine particles and polycyclic aromatic hydrocarbons and birth outcomes: a two-pollutant approach. Int Arch Occup Environ Health. 2017;90(3):255–64.
Welch BM, Keil AP, Buckley JP, Calafat AM, Christenbury KE, Engel SM, et al. Associations between prenatal urinary biomarkers of Phthalate exposure and Preterm Birth: a pooled study of 16 US cohorts. JAMA Pediatr. 2022;176(9):895–905.
van den Dries MA, Keil AP, Tiemeier H, Pronk A, Spaan S, Santos S, et al. Prenatal exposure to Nonpersistent Chemical Mixtures and fetal growth: a Population-based study. Environ Health Perspect. 2021;129(11):117008.
Zhang Y, Dong S, Wang H, Tao S, Kiyama R. Biological impact of environmental polycyclic aromatic hydrocarbons (ePAHs) as endocrine disruptors. Environ Pollut. 2016;213:809–24.
Diamanti-Kandarakis E, Bourguignon JP, Giudice LC, Hauser R, Prins GS, Soto AM, et al. Endocrine-disrupting chemicals: an endocrine Society scientific statement. Endocr Rev. 2009;30(4):293–342.
Duh-Leong C, Maffini MV, Kassotis CD, Vandenberg LN, Trasande L. The regulation of endocrine-disrupting chemicals to minimize their impact on health. Nat Reviews Endocrinol. 2023;19(10):600–14.
Ho S-M, Cheong A, Adgent MA, Veevers J, Suen AA, Tam NNC, et al. Environmental factors, epigenetics, and developmental origin of reproductive disorders. Reprod Toxicol. 2017;68:85–104.
Breton CV, Marsit CJ, Faustman E, Nadeau K, Goodrich JM, Dolinoy DC, et al. Small-magnitude effect sizes in epigenetic end points are important in children’s Environmental Health studies: the Children’s Environmental Health and Disease Prevention Research Center’s epigenetics Working Group. Environ Health Perspect. 2017;125(4):511–26.
Perera F, Herbstman J. Prenatal environmental exposures, epigenetics, and disease. Reprod Toxicol. 2011;31(3):363–73.
American Academy of Pediatrics Council on Environmental Health. Endocrine disrupting chemicals. In: Etzel RA, Balk SJ, editors. Pediatric Environmental Health. Itasca, IL: American Academy of Pediatrics; 2019. pp. 511–20.
Baccarelli A, Pesatori AC, Bertazzi PA. Occupational and environmental agents as endocrine disruptors: experimental and human evidence. J Endocrinol Investig. 2000;23(11):771–81.
Sol CM, Gaylord A, Santos S, Jaddoe VWV, Felix JF, Trasande L. Fetal exposure to phthalates and bisphenols and DNA methylation at birth: the Generation R Study. Clin Epigenetics. 2022;14(1):125.
Lee J, Kim J, Zinia SS, Park J, Won S, Kim WJ. Prenatal phthalate exposure and cord blood DNA methylation. Sci Rep. 2023;13(1):7046.
Jedynak P, Tost J, Calafat AM, Bourova-Flin E, Broséus L, Busato F, et al. Pregnancy exposure to phthalates and DNA methylation in male placenta — an epigenome-wide association study. Environ Int. 2022;160:107054.
Miura R, Araki A, Minatoya M, Miyake K, Chen ML, Kobayashi S, et al. An epigenome-wide analysis of cord blood DNA methylation reveals sex-specific effect of exposure to bisphenol A. Sci Rep. 2019;9(1):12369.
Solomon O, Yousefi P, Huen K, Gunier RB, Escudero-Fung M, Barcellos LF, et al. Prenatal phthalate exposure and altered patterns of DNA methylation in cord blood. Environ Mol Mutagen. 2017;58(6):398–410.
Jedynak P, Tost J, Calafat AM, Bourova-Flin E, Busato F, Forhan A, et al. Pregnancy exposure to synthetic phenols and placental DNA methylation — an epigenome-wide association study in male infants from the EDEN cohort. Environ Pollut. 2021;290:118024.
Mortillo M, Marsit CJ. Select early-Life Environmental exposures and DNA methylation in the Placenta. Curr Environ Health Rep. 2023;10(1):22–34.
Vlahos A, Mansell T, Saffery R, Novakovic B. Human placental methylome in the interplay of adverse placental health, environmental exposure, and pregnancy outcome. PLoS Genet. 2019;15(8):e1008236.
Marsit CJ. Placental epigenetics in children’s Environmental Health. Semin Reprod Med. 2016;34(1):36–41.
Hogg K, Price EM, Hanna CW, Robinson WP. Prenatal and perinatal environmental influences on the human fetal and placental epigenome. Clin Pharmacol Ther. 2012;92(6):716–26.
Onyemauwa F, Rappaport SM, Sobus JR, Gajdosová D, Wu R, Waidyanatha S. Using liquid chromatography-tandem mass spectrometry to quantify monohydroxylated metabolites of polycyclic aromatic hydrocarbons in urine. J Chromatogr B Analyt Technol Biomed Life Sci. 2009;877(11–12):1117–25.
Ye X, Kuklenyik Z, Needham LL, Calafat AM. Automated On-Line column-switching HPLC-MS/MS method with peak focusing for the determination of Nine Environmental Phenols in urine. Anal Chem. 2005;77(16):5407–13.
Silva MJ, Samandar E, Preau JL, Reidy JA, Needham LL, Calafat AM. Quantification of 22 phthalate metabolites in human urine. J Chromatogr B. 2007;860(1):106–12.
Lubin JH, Colt JS, Camann D, Davis S, Cerhan JR, Severson RK, et al. Epidemiologic evaluation of measurement data in the presence of detection limits. Environ Health Perspect. 2004;112(17):1691–6.
Wilhelm-Benartzi CS, Koestler DC, Karagas MR, Flanagan JM, Christensen BC, Kelsey KT, et al. Review of processing and analysis methods for DNA methylation array data. Br J Cancer. 2013;109(6):1394–402.
Aryee MJ, Jaffe AE, Corrada-Bravo H, Ladd-Acosta C, Feinberg AP, Hansen KD, Irizarry RA. Minfi: a flexible and comprehensive Bioconductor package for the analysis of Infinium DNA methylation microarrays. Bioinformatics. 2014;30(10):1363–9.
Liu Y, Eliot MN, Papandonatos GD, Kelsey KT, Fore R, Langevin S, et al. Gestational perfluoroalkyl substance exposure and DNA methylation at Birth and 12 years of age: a longitudinal Epigenome-wide Association study. Environ Health Perspect. 2022;130(3):037005.
Triche TJ Jr, Weisenberger DJ, Van Den Berg D, Laird PW, Siegmund KD. Low-level processing of Illumina Infinium DNA methylation BeadArrays. Nucleic Acids Res. 2013;41(7):e90–e.
Teschendorff AE, Marabita F, Lechner M, Bartlett T, Tegner J, Gomez-Cabrero D, Beck S. A beta-mixture quantile normalization method for correcting probe design bias in Illumina Infinium 450 k DNA methylation data. Bioinformatics. 2013;29(2):189–96.
Johnson WE, Li C, Rabinovic A. Adjusting batch effects in microarray expression data using empirical Bayes methods. Biostatistics. 2007;8(1):118–27.
Zhou W, Laird PW, Shen H. Comprehensive characterization, annotation and innovative use of Infinium DNA methylation BeadChip probes. Nucleic Acids Res. 2017;45(4):e22.
Houseman EA, Accomando WP, Koestler DC, Christensen BC, Marsit CJ, Nelson HH, et al. DNA methylation arrays as surrogate measures of cell mixture distribution. BMC Bioinformatics. 2012;13(1):86.
Salas LA, Koestler DC, Butler RA, Hansen HM, Wiencke JK, Kelsey KT, Christensen BC. An optimized library for reference-based deconvolution of whole-blood biospecimens assayed using the Illumina HumanMethylationEPIC BeadArray. Genome Biol. 2018;19(1):64.
Yuan V, Hui D, Yin Y, Peñaherrera MS, Beristain AG, Robinson WP. Cell-specific characterization of the placental methylome. BMC Genomics. 2021;22(1):6.
Yuan V. Placental DNA methylation analysis tools. 1.8 ed: Bioconductor; 2023.
Haeussler M, Zweig AS, Tyner C, Speir ML, Rosenbloom KR, Raney BJ, et al. The UCSC Genome Browser database: 2019 update. Nucleic Acids Res. 2018;47(D1):D853–8.
Hansen KD. IlluminaHumanMethylation450kanno. ilmn12. hg19: annotation for Illumina’s 450k methylation arrays. R Package Version. 2016;06:0.
Du P, Zhang X, Huang C-C, Jafari N, Kibbe WA, Hou L, Lin SM. Comparison of Beta-value and M-value methods for quantifying methylation levels by microarray analysis. BMC Bioinformatics. 2010;11(1):587.
Si J, Meir AY, Hong X, Wang G, Huang W, Pearson C, et al. Maternal pre-pregnancy BMI, offspring epigenome-wide DNA methylation, and childhood obesity: findings from the Boston Birth Cohort. BMC Med. 2023;21(1):317.
Fernandez-Jimenez N, Fore R, Cilleros-Portet A, Lepeule J, Perron P, Kvist T, et al. A meta-analysis of pre-pregnancy maternal body mass index and placental DNA methylation identifies 27 CpG sites with implications for mother-child health. Commun Biology. 2022;5(1):1313.
Barfield RT, Kilaru V, Smith AK, Conneely KN. CpGassoc: an R function for analysis of DNA methylation microarray data. Bioinformatics. 2012;28(9):1280–1.
Benjamini Y, Drai D, Elmer G, Kafkafi N, Golani I. Controlling the false discovery rate in behavior genetics research. Behav Brain Res. 2001;125(1):279–84.
Devlin B, Roeder K. Genomic Control for Association Studies. Biometrics. 1999;55(4):997–1004.
Phipson B, Maksimovic J, Oshlack A. missMethyl: an R package for analyzing data from Illumina’s HumanMethylation450 platform. Bioinformatics. 2016;32(2):286–8.
Williams CJ, Li Z, Harvey N, Lea RA, Gurd BJ, Bonafiglia JT, et al. Genome wide association study of response to interval and continuous exercise training: the Predict-HIIT study. J Biomed Sci. 2021;28(1):37.
Core Team R. R. R: a language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing; 2022.
Shapiro SS, Wilk MB. An analysis of variance for normality (complete samples). Biometrika. 1995;52(3):591.
Stelzer G, Rosen N, Plaschkes I, Zimmerman S, Twik M, Fishilevich S et al. The GeneCards Suite: From Gene Data Mining to Disease Genome Sequence Analyses. Curr Protoc Bioinformatics. 2016;54:1.30.1-1.3.
Safran M, Rosen N, Twik M, BarShir R, Stein TI, Dahary D, et al. The GeneCards suite. In: Abugessaisa I, Kasukawa T, editors. Practical guide to Life Science Databases. Singapore: Springer Nature Singapore; 2021. pp. 27–56.
Li J, Zhu X, Yu K, Jiang H, Zhang Y, Wang B, et al. Exposure to Polycyclic Aromatic Hydrocarbons and accelerated DNA methylation aging. Environ Health Perspect. 2018;126(6):067005.
Herbstman JB, Tang D, Zhu D, Qu L, Sjödin A, Li Z, et al. Prenatal Exposure to Polycyclic Aromatic Hydrocarbons, Benzo[< i > a]pyrene–DNA Adducts, and Genomic DNA Methylation in Cord Blood. Environmental Health Perspectives. 2012;120(5):733-8.
Martin EM, Fry RC. Environmental influences on the Epigenome: exposure- Associated DNA methylation in human populations. Annu Rev Public Health. 2018;39:309–33.
Yang P, Gong YJ, Cao WC, Wang RX, Wang YX, Liu C, et al. Prenatal urinary polycyclic aromatic hydrocarbon metabolites, global DNA methylation in cord blood, and birth outcomes: a cohort study in China. Environ Pollut. 2018;234:396–405.
Das DN, Ravi N. Influences of polycyclic aromatic hydrocarbon on the epigenome toxicity and its applicability in human health risk assessment. Environ Res. 2022;213:113677.
Perera F, Tang W-y, Herbstman J, Tang D, Levin L, Miller R, Ho S-m. Relation of DNA methylation of 5′-CpG island of ACSL3 to Transplacental exposure to Airborne Polycyclic Aromatic Hydrocarbons and Childhood Asthma. PLoS ONE. 2009;4(2):e4488.
Li K, Yin R, Wang Y, Zhao D. Associations between exposure to polycyclic aromatic hydrocarbons and metabolic syndrome in U.S. adolescents: cross-sectional results from the National Health and Nutrition Examination Survey (2003–2016) data. Environ Res. 2021;202:111747.
Li W, Chen D, Peng Y, Lu Z, Wang D. Association of polycyclic aromatic hydrocarbons with systemic inflammation and metabolic syndrome and its components. Obes (Silver Spring). 2023;31(5):1392–401.
Lichtiger L, Jezioro J, Rivera J, McDonald JD, Terry MB, Sahay D, Miller RL. Prenatal airborne polycyclic aromatic hydrocarbon exposure, altered regulation of peroxisome proliferator-activated receptor gamma (ppar)γ, and links with mammary cancer. Environ Res. 2023;231(Pt 2):116213.
Jedynak P, Tost J, Calafat AM, Bourova-Flin E, Broséus L, Busato F, Forhan A, Heude B, Jakobi M, Schwartz J, Slama R, Vaiman D, Lepeule J, Philippat C. Pregnancy exposure to phthalates and DNA methylation in male placenta - an epigenome-wide association study. Environ Int. 2022;160:107054.
Pidsley R, Zotenko E, Peters TJ, Lawrence MG, Risbridger GP, Molloy P, et al. Critical evaluation of the Illumina MethylationEPIC BeadChip microarray for whole-genome DNA methylation profiling. Genome Biol. 2016;17(1):208.
Kaur D, Lee SM, Goldberg D, Spix NJ, Hinoue T, Li H-T, et al. Comprehensive evaluation of the Infinium human MethylationEPIC v2 BeadChip. Epigenetics Commun. 2023;3(1):6.
Acknowledgements
We acknowledge the staff assistance from the Department of Obstetrics and Gynecology at the University of Cincinnati in participant recruitment. We are thankful to our study participants for their time and efforts.
Funding
This research was funded by pilot projects from the National Institute of Environmental Health Sciences (NIEHS), supported by the Center of Environmental Genetics (P30ES006096) at the University of Cincinnati and the Center of Excellence in Environmental Toxicology (P30ES013508) at the University of Pennsylvania. This project was also supported by grants from the Veteran Affairs [VA-I01BX005395 (SMH, YKL), VA‐IK6BX006182 (SMH)], NIEHS [R01ES032675 (SMH, YKL), R01ES028277 (AC), R01ES033054 (AC), R01ES032836 (JMB)], Department of Defense [DoD‐W81XWH‐22‐1‐0152 (SMH, YKL, AC)], and National Science Foundation [RII Track‐2 FEC, Award #2217824 (SMH, YKL)].
Author information
Authors and Affiliations
Contributions
JP, JB, EDF, SMH, YLK, SH, and AC – conceptualization and methodology. EDF – collected biospecimens. AMC and JCB – quantified maternal urinary biomarkers except PAHs. XZ – performed the methylation. AMV, SK, ZP and AC – acquired the data. JB, SMH, YLK, SH, AMV, and AC – acquired funding to conduct this study. JP – performed the statistical analysis under the guidance of AC, JB, SMH, YLK, and SH. JP – writing the original draft. All authors reviewed, edited, and approved the final version of the manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
All procedures performed in this study involving human subjects were in compliance with the ethical standards of the Institutional Review Board at the University of Cincinnati. Informed consent was obtained from all the subjects included in this study.
Consent for publication
All authors approve the manuscript for publication.
Competing interests
The authors declare no competing interests.
Disclaimer
The findings and conclusions in this manuscript are those of the authors and do not necessarily represent the official position of the Centers for Disease Control and Prevention (CDC). Use of trade names is for identification only and does not imply endorsement by the CDC, the Public Health Service, or the US Department of Health and Human Services.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Electronic supplementary material
Below is the link to the electronic supplementary material.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Puvvula, J., Braun, J.M., DeFranco, E.A. et al. Gestational exposure to environmental chemicals and epigenetic alterations in the placenta and cord blood mononuclear cells. Epigenetics Commun. 4, 4 (2024). https://doi.org/10.1186/s43682-024-00027-7
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s43682-024-00027-7