Trimester-Specific Associations of Prenatal Lead Exposure With Infant Cord Blood DNA Methylation at Birth.

Gestational exposure to lead (Pb) adversely impacts offspring health through multiple mechanisms, one of which is the alteration of the epigenome including DNA methylation. This study aims to identify differentially methylated CpG sites associated with trimester-specific maternal Pb exposure in umbilical cord blood (UCB) leukocytes. Eighty-nine mother-child dyads from the Early Life Exposure in Mexico to Environmental Toxicants (ELEMENT) longitudinal birth cohorts with available UCB samples were selected for DNA methylation analysis via the Infinium Methylation EPIC BeadChip, which quantifies methylation at >850 000 CpG sites. Maternal blood lead levels (BLLs) during each trimester (T1: 6.56 ± 5.35 µg/dL; T2: 5.93 ± 5.00 µg/dL; T3: 6.09 ± 4.51 µg/dL), bone Pb (patella: 11.8 ± 9.25 µg/g; tibia: 11.8 ± 6.73 µg/g), a measure of cumulative Pb exposure, and UCB Pb (4.86 ± 3.74 µg/dL) were measured. After quality control screening, data from 786 024 CpG sites were used to identify differentially methylated positions (DMPs) and differentially methylated regions (DMRs) by Pb biomarkers using separate linear regression models, controlling for sex and estimated UCB cell-type proportions. We identified 3 DMPs associated with maternal T1 BLL, 2 with T3 BLL, and 2 with tibia bone Pb. We identified one DMR within PDGFRL associated with T1 BLL, one located at chr6:30095136-30095295 with T3 BLL, and one within TRHR with tibia bone Pb (adjusted P-value < .05). Pathway analysis identified 15 overrepresented gene pathways for differential methylation that overlapped among all 3 trimesters with the largest overlap between T1 and T2 (adjusted P-value < .05). Pathways of interest include nodal signaling pathway and neurological system processes. These data provide evidence for differential methylation by prenatal Pb exposure that may be trimester-specific.

Background

Lead (Pb) is a ubiquitous environmental pollutant found in air, soil, water, and food. Human exposures mainly occur through ingestion and inhalation, and Pb can be stored long term in teeth and bones. Pb is a potent neurotoxicant that alters brain development resulting in abnormal cognition and behaviors. Early life Pb exposure also increases the risk for developing a variety of adverse health outcomes later in life, including Alzheimer disease, attention-deficit hyperactivity disorder, cardiovascular disease, and intellectual deficits.[1-5] The Institute for Health Metrics and Evaluation (IHME) estimated in 2017 that Pb exposure accounted for almost 2% of total all-cause mortality and at least 2.5 million years of healthy life lost (disability-adjusted life years [DALYs]) worldwide due to idiopathic developmental intellectual disability,[6] but this may be an underestimate.[7]

The health effects of prenatal Pb exposure can be framed with the Developmental Origins of Health and Disease (DoHAD) theory,[8] which links early-life exposures to the development of disease and adverse health effects later in life. The placenta is not an effective barrier against Pb, and as bone reabsorption increases during pregnancy, Pb is released into the bloodstream from where it can reach the fetus potentially causing abnormalities in development.[9-13] The fetus may be exposed to stored Pb in bone from the mother’s previous exposures nearly decades prior as well as to concurrent Pb exposure circulating through the blood.

It has been hypothesized that toxicant exposures, including Pb, may result in epigenetic modifications during sensitive developmental periods, thereby impacting gene regulation and subsequent development and later life health outcomes. Epigenetic modifications, which include DNA methylation, are mitotically heritable and regulate gene expression without altering the underlying DNA sequence. DNA methylation occurs primarily in cytosine-guanine (CpG) dinucleotides, which are enriched in regions called islands, typically located in promoter regions of genes.[14,15] Establishment of DNA methylation patterns occurs in early gestation and disturbances during this period may result in DNA methylation variations that are propagated through mitosis to new cells and developing organs thus affecting gene expression in a wide range of tissues associated with various developmental physiological processes.[16]

We and others have provided evidence for the impact of early-life Pb exposure on the epigenome in animal and human studies. In mice, dams were exposed to 1 of 3 doses of Pb acetate in water prior to mating through weaning, and the offspring were followed through adulthood and compared with unexposed mice. Pb was associated with altered DNA methylation in offspring tail tips immediately following cessation of exposure (3 weeks of age) at 2 metastable epialleles.[17] Altered DNA methylation by Pb was also observed in offspring brain (10 months of age) at environmentally responsive retrotransposons.[18] In a human birth cohort, we quantified DNA methylation at 4 genes in 247 umbilical cord blood (UCB) leukocyte DNA samples. Biomarkers of prenatal Pb exposure (maternal tibia: 9.82 [SD =9.45] µg/dL, maternal patella: 14.4 [SD = 14.8] µg/dL, UCB: 6.49 [SD = 3.52] µg/dL) were associated with percent DNA methylation in neonates in long interspersed elements 1 (LINE-1), as well as growth-related genes (IGF2 and HSD11B2) providing evidence for maternal cumulative Pb burden influencing the epigenome of a developing fetus.[19] Finally, the US prospective pregnancy cohort, Project Viva, conducted an epigenome-wide analysis on 268 UCB samples to evaluate the association between prenatal maternal Pb exposure (mid-to-late gestation Pb in erythrocytes averaged 1.22 [SD = 0.63] µg/dL) and DNA methylation.[20] Results identified sex-specific differentially methylated CpG sites by Pb exposure with more found in females (n = 38) than males (n = 2). Mechanistically, metal exposures such as Pb increase oxidative stress and can cause oxidative DNA damage which inhibits the ability of methyltransferases to bind to DNA; this leads to hypomethylation at some loci.[21-23] In addition, Pb exposure has been shown to alter expression and function of epigenetic machinery (ie, DNA methyltransferases [DNMTs]). Several human and animal studies suggest DNA methylation changes from developmental Pb exposure are due to altered DNMT expression and activity.[21,24-27]

Despite mounting evidence for Pb’s effect on the epigenome, no studies have specifically investigated the association between genome-wide DNA methylation and Pb exposure during each of the 3 trimesters of pregnancy, as well as cumulative gestational Pb exposure estimated by maternal patella and tibia measures. We hypothesize that prenatal Pb exposure will be significantly associated with changes in DNA methylation patterns across the genome at birth and that these changes will differ depending on the timing of exposure. A whole genome approach will assist researchers in gaining valuable insight into the effects of prenatal Pb exposure on gene regulation and help to identify potential biomarkers of interest specific to each trimester and cumulative exposure. Here, we performed Pb exposure assessment in multiple biomarkers and DNA methylation profiling in UCB samples from 89 mother-infant pairs enrolled in the Early Life Exposures in Mexico to Environmental Toxicants (ELEMENT) project, a longitudinal birth cohort.

Methods

Study population

The ELEMENT project has used a series of longitudinal birth cohorts to investigate the influence of exposure to Pb and other toxicants—in utero and in childhood—on sensitive periods of development. This project is based on 3 sequentially recruited cohorts comprising 1643 mother-infant pairs, some of whom have been followed for 20+ years.[28] This study uses data and biological samples from the second and third birth cohorts, for which 1530 women were originally enrolled and 1012 mother-infant pairs were followed up after birth. Women were recruited from 1997-2000 and 2001-2005, from the Mexican Social Security Institute hospital in Mexico City. Eligibility and exclusion criteria are as previously described.[28,29] For all recruited mother-child pairs, data collected include sex, gestational age, socioeconomic status, anthropometry data, and other environmental exposures. Families were followed up at multiple time points from infancy through adolescent years. For this study, we selected 97 ELEMENT participants from the second and third cohorts who had archived UCB samples from which we could isolate DNA. The final study comprised 89 mother-infant dyads whose epigenetic data passed downstream quality control (QC). Characteristics of these 89 participants compared with all women of cohorts 2 and 3 that had least one Pb biomarker measurement (n = 1214) and their children are represented in Supplemental Table 1.

At the time of enrollment, all mothers were informed about the study; those who agreed to participate read and signed a letter of informed consent about the original study. The original research protocol and all amendments to the study protocol were approved by the Ethics Committees of the National Institutes of Public Health of Mexico, participating hospitals, and the Internal Review Board at all participating institutions including the University of Michigan.

Pb exposure assessment and genomic DNA isolation

UCB was collected within 12 hours after birth. Cohort 2 blood lead levels (BLL) from each trimester and UCB were measured using inductively coupled plasma mass-spectrometry (ICP-MS, Thermo Finnigan, Bremen, Germany) at the University of California, Santa Cruz, as described previously.[30] Cohort 3 BLLs in maternal venous blood from each trimester and UCB were measured at the Trace Metal Laboratory of the American British Cowdray Hospital using graphite furnace atomic absorption spectrometry (instrument model 3000; PerkinElmer, Norwalk, CT, USA).

Bone Pb was measured in maternal left patella (trabecular bone) and mid-shaft of the left tibia (cortical bone) 1 month post-partum as an indicator of cumulative Pb exposure during pregnancy using a spot-source 109Cd K-shell X-ray fluorescence (K-XRF) instrument. The technical specifications and validation of this instrument are described in detail elsewhere.[31] Analysis of means and standard deviations once weekly for QC and calibration measures did not disclose any significant shifts in precision or accuracy. Tibia and patella bone Pb values below the limit of detection (ie, negative values) were dropped from subsequent analyses resulting in 73 patella bone Pb and 46 tibia bone Pb measures.

DNA was isolated from UCB leukocytes using Qiagen kits and standard protocols for blood DNA isolation. Nucleic acid yield and purity were assessed first using a NanoDrop spectrophotometer (ThermoFisher Scientific), and double-stranded DNA was also quantified via a Qubit fluorometer. All DNA samples were stored at −80°C until later use.

Epigenetic analysis

All genomic DNA samples were bisulfite converted (500 ng) using the Zymo EZ DNA Methylation kit (Zymo Research) and the recommended incubation conditions and methods for downstream Infinium analysis. In short, this treatment converts unmethylated cytosines to uracils while leaving methylated cytosines intact. DNA methylation was quantified at 863 904 CpG sites following hybridization to the Infinium MethylationEPIC BeadChip (Illumina).[32] Beadchips were processed and scanned on the Illumina iScan at the University of Michigan (UM) Advanced Genomics Core using the Infinium Methylation EPIC kit.

We used a preprocessing pipeline for the resulting raw EPIC data based on functional normalization[33] that uses the “minfi” package in R.[34] The pipeline included several steps: loading raw image files, linking to phenotypic information, dropping probes with single-nucleotide polymorphisms (SNPs) within query sites, removing poor-quality samples (>5% failed probes) and poorly detected probes (P-value > 1e−6 when comparing to background), normalizing the data to correct for background signal and dye bias, checking and correcting for batch effects (by slide and array), and removing CpG sites on the X and Y chromosomes. After QC and normalization, 786 024 CpG sites and 89 samples passed all QC measures and were subsequently used in statistical analysis.

Estimates of cell-type composition (T lymphocytes [CD4T, CD8T], B cells, NK cells, monocytes, granulocytes, and nucleated red blood cells [nRBCs]) for each sample was performed using an established method based on UCB cell-type-specific differentially methylated regions. Cell-type proportions were adjusted for in downstream statistical analysis to avoid confounding by cell-type composition.[35] Beta-values and M-values (log2-transformed beta-values) were then calculated from the processed data and used in statistical analysis.[36]

Statistical analysis

All statistical analyses were performed in the R Project for Statistical Computing (version 3.4.3). Analyses were performed to identify differentially methylated loci by Pb exposure. Pb biomarkers tested were maternal BLLs at each trimester, UCB Pb levels, and Pb levels in maternal patella and tibia, which were treated as non-transformed continuous variables. Bone and UCB Pb levels were normally distributed. Trimester-specific Pb measures were right-skewed yet all biomarker concentrations were biologically plausible and relevant to this study population. As such the relationships between maternal BLLs and DNA methylation were tested with and without BLL outliers. We performed bivariate analyses relating Pb biomarkers and covariates (eg, cell-type proportions, gestational age, maternal age, maternal education, household income category) to identify potential confounders. We also performed singular value decomposition to identify covariates associated with principal components of the DNA methylation data. We then chose potential confounding variables to include in the final statistical models of site-specific DNA methylation data that were associated with both Pb and DNA methylation. We also retained several variables that had a large effect on DNA methylation but were not associated with Pb (ie, infant sex, cell type proportions). In models with BLLs, offspring sex and estimated cell-type compositions of granulocytes and nRBC[35] were included. Models with bone Pb measures as the predictor of interest additionally included maternal age and cohort. Other covariates available in this study population, including gestational age, were not considered as confounders because they were neither associated with prenatal Pb exposure nor DNA methylation in our study sample.

Differentially methylated positions (DMPs) by Pb were identified using the “limma” package.[37,38] Briefly, linear regression models were run for each CpG site using M-values, which are logit-transformed beta values representing the proportion of methylation at a given CpG site, with one Pb exposure biomarker as the predictor of interest and the aforementioned covariates included. Limma calculates fold-change in methylation, t-statistics, and P-values by applying empirical Bayes smoothing to the standard errors. To account for multiple testing, we used the false discovery rate (FDR) method[39] and consider DMPs with corrected P-value (q-value) < .05 to be statistically significant. For comparison purposes between studies, we also report estimates (β) from the same regression analysis with betas (non-transformed methylation values) as the outcome.

We used the statistical package DMRcate to test associations between each Pb exposure variable and differential methylation at the regional level across the entire genome.[40] The Gaussian kernel bandwidth lambda values used was the default of 1000, and 2 CpG sites were considered the minimum consecutive loci to be included in a differentially methylated region (DMR). Differentially methylated regions with FDR < 0.05 were statistically significant.

We used LRpath, a multiple-comparison enrichment testing tool that uses logistic regression to test for gene sets that have higher significance values by raw P-value than expected at random,[41] to identify gene pathways enriched for hyper- or hypo-methylation by Pb exposure. It considers CpG sites in known genes, their respective log fold-changes by Pb from linear regression models, and associated P-values.

Results

Descriptive statistics of ELEMENT study data

Among the 89 mothers included in this analysis, 41 (46%) had male infants (Table 1). The mean age of women at offspring birth was 26.4 (SD = 4.8) years. The average gestational age at birth was 39.0 (SD = 1.1) weeks. The mean (SD) maternal BLLs averaged over all 3 trimesters was 6.18 (4.51) µg/dL, with the first trimester (T1) mean of 6.56 (5.35) µg/dL, second trimester (T2) mean of 5.93 (5.00) µg/dL, and third trimester (T3) mean of 6.09 (4.51) µg/dL. Maternal BLLs between the trimesters were highly correlated (r > 0.63) according to Spearman rank-order correlation test. Average Pb concentration in patella was 9.42 (10.28) µg/g and in tibia was 7.51 (9.46) µg/g. Umbilical cord blood Pb levels were measured to represent offspring Pb exposure at birth and averaged 4.86 (3.74) µg/dL. Geometric means of BLL among women in this study were more than 4 times higher than women of childbearing ages within a similar time period in the United States according to NHANES (1999-2002 and 2003-2006).[42] Sources of Pb exposure in Mexico City during this time included the use of traditional lead-glazed ceramics for cooking, and Pb in gasoline which was not phased out in full until 1997.[43] Pregnant women in countries such as Malaysia, India, Iran, Malta, and Bangladesh continue to have exposures similar or up to 2½ times greater on average than that of the women in this study.[44-48]

Table 1.

Characteristics of ELEMENT mother-infant pairs with UCB DNA methylation data.

Characteristics	No.	Mean ± SD or N (%)	Range
Mothers
Age (y)	89	26.4 ± 4.81	18.0-37.0
Blood lead (µg/dL)
First trimester (T1)	69	6.56 ± 5.35	0.90-35.8
Second trimester (T2)	74	5.93 ± 5.00	0.80-38.2
Third trimester (T3)	76	6.09 ± 4.51	0.90-34.0
Average all trimesters	77	6.18 ± 4.51	1.17-33.1
Bone lead (µg/g)
Patella	73	11.8 ± 9.25	0.20-42.0
Tibia	46	11.8 ± 6.73	0.40-28.8
Household income	77
Lowest		12 (15.6)
Low-middle		31 (40.3)
Middle		20 (26.0)
Middle-high		9 (11.6)
Highest		5 (6.5)
Maternal education (total years)	89	11.0 ± 2.39	3.00-17.0
Children
UCB lead (µg/dL)	86	4.86 ± 3.74	0.00-19.5
Gestational age (wk)	87	39.0 ± 1.09	36.0-41.0
Male (%)	89	41 (46.0)

Abbreviations: UCB, umbilical cord blood.

We compared participants included in this analysis with all ELEMENT mother-infant pairs from the same cohorts (2 and 3). Demographics and Pb biomarker concentrations were not statistically different between the subset and the entire study (Supplemental Table 1) with the exception of gestational age which was slightly higher in the subset with DNA methylation data.

CpG site-specific analysis and pathways analysis

We found that DNA methylation levels at 3 CpG sites demonstrated significant associations (FDR < 0.05) with T1 BLLs, 1 CpG with T3 BLLs, and 2 with tibia bone Pb levels (Table 2). There were no statistically significant DMPs in analyses with maternal BLLs in T2, UCB Pb, or maternal patella bone Pb. Of the 3 statistically significant DMPs with the maternal BLLs during T1, the sites mapped to RAB5A, EXT1, and a non-genic region (chr8:12615485). The DMPs in RAB5A and EXT1 had less DNA methylation (−0.022% and −0.024%, respectively) per 1 µg/dL increase in maternal T1 BLL. In contrast, the non-genic region (chr8:12615485) was hypermethylated by Pb (0.35% per µg/dL of T1 BLL). The DMP by T3 maternal BLLs was hypomethylated (−0.52% per µg/dL BLL) and in a non-genic region (chr1:160839299). Maternal tibia bone Pb-associated DMPs were hypermethylated in LNFN1 (0.37% per µg/g bone Pb) and hypomethylated in a non-genic region located at chr6:792305 (−0.21% per µg/g bone Pb).

Table 2.

Statistically significant DMPs (q < 0.05) by first-trimester maternal BLL, third-trimester maternal BLL, and maternal tibia Pb.

Pb biomarker	Probe ID	Gene name	Chr	Pos	β estimate	P-value	Q-value
T1 BLL	cg17138393	RAB5A	chr3	19988887	−0.000221	5.46E-08	0.0227
	cg03390844		chr8	12615485	0.00348	5.78E-08	0.0227
	cg00984923	EXT1	chr8	119124069	−0.000235	1.49E-07	0.0390
T3 BLL	cg01328348		chr1	160839299	−0.00519	1.81E-08	0.0142
Tibia Pb	cg00002033	LRFN1	chr19	39798481	0.00367	3.57E-08	0.0281
	cg03463208		chr6	792305	−0.00205	7.23E-08	0.0284

Abbreviations: BLL, blood lead level; Chr, chromosome; DMP, differentially methylated position; pos, base pair position from hg19.

Results from the analysis are shown for CpG sites associated with Pb biomarker levels below q-value of 0.05 from models adjusted for sex and estimated cell-type proportions (granulocytes and nucleated red blood cells). The P-values and q-values are obtained from the analyses using methylation expressed as M-values as the outcome variable, and reported association estimates (β) are from analysis using methylation beta values as the outcome variable.

LRpath analysis was conducted to test for gene sets enriched among continuous significance values (raw P-values) by each Pb biomarker. After multiple-comparisons correction, there were 109 gene sets enriched with T1 BLLs, 76 with T2 BLLs, 28 with T3 BLLs, and 12 with UCB Pb levels, as well as 15 with patella and 68 with tibia bone measures that met a q-value cut-off of 5% (Supplemental Table 2). Differently methylated genes associated with Pb exposure across all 3 trimesters and at birth were more likely to be enriched in gene pathways involved in detection of chemical stimulus, sensory perception and smell, immunoglobin binding, transmembrane receptor activity, and olfactory receptor activity. Four pathways were statistically significant (q-value < 0.05) among all blood Pb variables, with 35 uniquely overlapping between T1 and T2, 4 uniquely overlapping between T2 and T3, none matching between T1 and T3 only, and 2 matching between T1 and at birth UCB (Figure 1). There were only 4 matching pathways when comparing results from models of tibia and patella bone Pb measures (Supplemental Figure 1).

Figure 1.

Venn diagram of functional annotations (q-value < 0.05) associated with maternal Pb exposure compared across all 3 trimesters (first trimester, T1; second trimester, T2; third trimester, T3) and umbilical cord blood (UCB).

Significant pathways from Gene ontology (GO) and Kyoto encyclopedia of genes and genomes (KEGG) pathways related to neurological development included neurological system processes associated with T1 (q-value = 2.87E10−6) and T2 (q-value = 1.77E10−3) BLLs and tibia bone Pb (q-value = 5.65E10−9). Nodal system pathway was also associated with T1 (q-value = 1.40E10−3) and T2 (q-value = 1.36E10−3) BLLs. Neurotransmitter receptor activity genes were associated with tibia bone Pb (q-value = 7.84E10−3).

Regional analysis

We tested for DMRs by Pb exposure for each trimester using DMRcate, which combines data from at least 2 consecutive CpG sites. DMRcate identified 1 region of differential methylation with T1 BLLs, 1 region with T3 BLLs, and 1 region with tibia bone Pb levels. No regions were identified with T2 BLLs, UCB Pb, or patella bone Pb. The region identified in T1 that was differentially methylated included 4 CpG sites in a region that is located between chr8:17433625-17433761 and is within the promoter of PDGFRL (Figure 2). The region identified in T3 included 12 CpG sites located at chr6:30095136-30095295 (Figure 3). Tibia bone Pb was associated with a region with 4 CpG sites within the promoter of TRHR (Figure 4).

Figure 2.

Differentially methylated region (DMR) by maternal blood Pb concentrations in T1 located within 200 base pairs of the transcription start site of PDGFRL (chr8:17433625-17433761). Linear regression modeling, treating Pb exposure as a continuous variable, adjusted for sex and estimated cell-type proportions. DNA methylation (beta value) is plotted for each sample at 6 CpG sites in the region included 4 sites of the statistically significant DMR. Exposure quartile ranges are represented from light red to dark red with quartile 4 (Q4), representing the top 25% most Pb exposed during T1, being the most methylated at each of the 4 CpG sites within the DMR.

Figure 3.

DMR by maternal blood Pb concentrations in T3 at chr6:30095136-30095295 includes 12 CpG sites within an intergenic region with no known regulatory function. This DMR was selected from a linear regression modeling, treating Pb exposure as a continuous variable, adjusting for sex and estimated cell-type proportions. DNA methylation (beta value) is plotted for each sample at the 12 CpG sites in the DMR. Exposure quartile ranges are represented from light red to dark red with Q4, representing the top 25% most Pb exposed during T3, being the most methylated at each of the 12 CpG sites within the DMR.

Figure 4.

DMR by maternal tibia bone Pb located within 800 base pairs of the transcription start site of TRHR (chr8:110098835-110098870). The DMR was identified by a linear regression model, treating Pb exposure as a continuous variable, and adjusting for sex and estimated cell-type proportions, maternal age, and cohort. DNA methylation (beta value) is plotted for each sample for the 4 CpG sites in the DMR. Exposure quartile ranges are represented from light red to dark red with Q4, representing the top 25% most Pb exposed when cumulative Pb is measured in tibia, being the most methylated at each of the 4 CpG sites within the DMR.

Discussion

We performed an epigenome-wide DNA methylation analysis using DNA extracted from UCB leukocytes and the Illumina MethylationEPIC beadchip in 89 mother-infant pairs to identify differentially methylated genes by trimester-specific and cumulative Pb exposure biomarkers. To our knowledge, this is the first genome-wide DNA methylation study to investigate the association of DNA methylation at birth with prenatal Pb exposure at each of the 3 trimesters of pregnancy, as well as cumulative prenatal Pb exposure in maternal patella and tibia bones. We found that maternal BLLs during T1 and T3, as well as maternal tibia bone Pb levels, were associated with differential methylation at several CpG sites in UCB.

Epidemiological studies of Pb-exposed pregnant women have reported that BLLs below 5 µg/dL can have adverse impacts on neurodevelopment, manifesting for example as decreased scores in the neonatal behavioral neurological assessment test.[49] Maternal BLLs under 10 µg/dL have been associated with a 3.5- to 6-point decrease in the Mental Development Index, a measure of IQ using the Bayley Scales of Infant Development, in 24-month-old infants (n = 146).[29] Although the inverse association between prenatal Pb exposure and neurodevelopment is known, whether epigenetic perturbation is one of the mechanisms underlying this association is not known.

In the epigenome-wide analysis performed here, the statistically significant DMPs and DMRs associated with biomarkers of prenatal Pb exposure annotated to several genes involved in neurodevelopment or neurologic function. In this study, DNA methylation levels at 3 CpG sites were significantly associated with maternal T1 BLLs, for which one is in a DNase I hypersensitivity site (DHS) within the first exon of RAB5A, a Ras-related protein. DNase I hypersensitivity site are regulatory regions of the genome often located at transcription start sites where the DNA is accessible for transcription factor binding. RAB5A is an important mediator of endocytic pathways due to its role in intracellular membrane trafficking and fusion, and localizing early endosomes.[50-53] Evidence from an in vitro study indicates that loss of RAB5A activity may perturb endosomal dynamics and may underlie neuronal dysfunction and degeneration in various motor neuron and neurodegenerative diseases.[54] The second CpG site associated with T1 BLLs is located at chr8:12615485, a DHS that is also classified to be within an enhancer region (according to FANTOM5)[55] but unclassified to date in terms of the specific gene(s) it may regulate. The third CpG site is in a DHS, a suspected enhancer, and is within 200 base pairs of the transcription start site of EXT1, exostosin glycosyltransferase 1. EXT1 has a primary role in the biosynthesis of heparan sulfate, which is known to regulate various biological processes such as cell proliferation, growth factor signaling, and embryonic development, specifically mammalian neuronal development.[56-59]

Maternal T1 Pb exposure was also associated with differential methylation in UCB at the regional level, identified by analyzing consecutive probes. We observed a region of 4 CpG sites located at chr8:17433625-17433761 within 200 base pairs of the transcription start site of PDGFRL, platelet-derived growth factor receptor-like protein. This gene encodes a known tumor suppressor that inhibits the growth of colorectal cancer cells in vitro; the function of this gene in development, if any, is not yet known.[60]

We observed one association between DNA methylation and T3 BLL. The CpG site associated with T3 BLLs is located at chr1:160839299 in an intergenic region within a DHS. In the regional analysis, there was one DMR that was statistically associated with maternal T3 Pb exposure. It annotated to chr6: 30095136-30095295 with 12 methylation sites within an intergenic region with no known regulatory function.

DNA methylation levels at 2 CpG sites were significantly associated with tibia bone Pb, a measure of cumulative Pb exposure throughout and before pregnancy. The first is within a DHS located in the gene LRFN1, Leucine Rich Repeat and Fibronectin Type III Domain Containing 1. LRFN proteins promote neurite outgrowth and synapse formation in hippocampal neurons and have a role in the developing nervous system with involvement in the regulation and conservation of excitatory synapses.[61-63] The second is located within a CpG island at chr6:792305 of no known regulatory function to date. We identified one region of differential methylation associated with tibia bone Pb located in the TRHR gene, thyrotropin-releasing hormone receptor, on chr8:110098835-110098870 within 800 base pairs of the transcription start site. This receptor signals synthesis, secretion, and bioactivity of the thyroid stimulating hormone within the pituitary gland in the brain by activating the phosphatidylinositol-calcium-protein kinase C transduction pathway.[64,65] We did not identify associations between patella Pb and DNA methylation. Pb in tibia has a residence time of 25 to 30 years[66] while patella Pb has a half-life ranging from months to years.[67] While both serve as proxies for gestational Pb exposure, differences in half-life and bioavailability (ie, following bone turnover) may influence the relationship between maternal bone Pb and fetal exposure dose throughout pregnancy.

We hypothesized that the greatest number of DMPs by Pb would be identified for T1, yet only a few significant results were obtained for T1, along with T3 and bone Pb. During early gestation, there are 2 major events influencing epigenetic programming: (1) post-fertilization there is a massive wave of demethylation to trigger embryonic developmental; (2) following, there is a subsequent re-methylation that shapes the epigenome in first cell lineages.[68,69] These events occur during early T1, and as such environmental perturbation during T1 could lead to epigenetic reprogramming that is propagated across all germ layers and tissues. Differentiation, which is critical for nervous system development, also occurs during T1.[70] Throughout the remainder of fetal development, dividing somatic cells are methylated via maintenance epigenetic machinery.[16] Thus, while gestational exposures are expected to have the most dramatic effects on the DNA methylome across tissues, subtle changes can still occur in mid- to late-gestation that may be tissue-specific.

Our study is one of a growing body of literature that shows associations between gestational Pb exposure and offspring DNA methylation in rodent models and humans.[17-20,71-73] While the DMPs and DMRs identified in this study by T1 and T3 maternal BLL and maternal bone Pb were not the same as those previously reported in other cohorts, collectively these studies show that Pb exposure and DNA methylation vary by timing of exposure, sex, and population. A rural Bangladesh birth cohort study with high environmental Pb exposure performed epigenome-wide DNA methylation profiling using the Infinium HumanMethylation450K BeadChip (450K) on UCB and assessed its association with prenatal Pb exposure measured in maternal urine at 8 weeks’ gestation and in maternal erythrocytes at 14 weeks’ gestation (n = 127).[74] This study identified 9 loci associated with urine Pb and 2 loci with erythrocyte Pb. Another study, a US prospective pregnancy cohort, Project Viva, with relatively low levels of Pb exposure (mean maternal erythrocyte Pb, 1.22 ± 0.63 μg/dL) conducted an epigenome-wide analysis using 450K on 268 UCB samples to evaluate the association between prenatal Pb exposure measured on average at 27.9 weeks’ gestation and DNA methylation.[20] Their results identified 4 CpG sites in all newborns, as well as sex-specific CpG sites with more found in females (n = 38) than males (n = 2), associated with prenatal Pb. Differences in findings from these cohorts and ELEMENT could be due to differences in exposure levels between the cohorts, participant racial/ethnic background, sex differences, timing of Pb exposure assessment, and statistical power to detect true associations. Even so, our study adds to the previous evidence that prenatal Pb exposure is associated with altered DNA methylation with likely small effect sizes, and this potential mechanism of toxicity or biomarker of response to exposure should be studied in larger cohorts and/or in consortiums.

While the mechanism by which Pb perturbs DNA methylation is not entirely known, several studies provide evidence for direct or indirect changes in DNMT expression and activity by Pb. DNA methyltransferases add methyl groups to CpG sites; if DNMTs fail to add methyl marks during cellular replication, DNA can become passively demethylated.[75,76] Pb exposure has been shown to increase reactive oxygen species which inhibits the DNMT binding to DNA causing hypomethylation at some but not all loci.[21-23,77,78] In addition, several studies have provided evidence for Pb exposure effects on DNMT expression and function which could result in either increases or decreases in methylation.[21,24-27]

Our study had several strengths including Pb exposure assessment at multiple time periods during gestation, epigenome-wide assessment of DNA methylation, a prospective study design, and rich data on key covariates. At the same time, our study has limitations. First, due to the sample size, we were unable to perform a sex-specific analysis to determine how sex may influence the relationship between DNA methylation and developmental Pb exposure. We were also underpowered to detect all “true” differentially methylated genes by Pb exposure, especially those with small effect sizes. Second, although epigenomic changes in early gestation would be expected to propagate across all tissues, it is still important to consider tissue specificity when conducting differential methylation studies. While brain is the primary tissue of interest in terms of Pb toxicity, blood-based epigenetic measures are typically necessary in longitudinal epidemiological studies. Cross-tissue studies should be conducted to better understand the biological connection between our findings in UCB with DNA methylation in the brain.[79] The Toxicant Exposures and Responses by Genomic and Epigenomic Regulators of Transcription (TaRGET II) consortium is currently promoting efforts toward understanding the role of environment in disease susceptibility as a function of epigenomic perturbations across target and surrogate tissue using mouse models (ie, blood and brain), which will provide insights on cross-tissue, inter-individual variation of epigenomic marks.[80] Finally, although DNA methylation changes would be expected to alter gene expression, we are unable to test this hypothesis; storage conditions of archived UCB samples in this cohort did not allow for RNA isolation.

In conclusion, we found that T1, T3, and the cumulative measure in tibia of gestational Pb exposure were associated with several statistically significant changes in DNA methylation. In total, 5 of the 6 CpG sites identified within the study were located within DHS which are important functional regions for gene regulation. There was no statistically significant differential methylation by T2 BLL, cord blood Pb, and patella bone Pb. Our findings are of importance because prenatal Pb exposure has been previously associated with adverse neurodevelopmental outcomes, but whether epigenetic mechanisms contribute to these long-term effects is not well characterized. Pb still remains a widespread environmental health problem. Dietary changes and micronutrient supplementation during pregnancy could play an essential role in neutralizing epigenomic perturbations due to environmental Pb exposure,[81-83] and the concentration of Pb in circulation can also be reduced by dietary supplementation (ie, to calcium).[84] Many public health promotion efforts are centered on primary prevention of Pb exposure in early childhood even though gestation is also a critical developmental period when Pb can affect long-term health. Our results suggest that prenatal Pb exposure may modify DNA methylation profiles at birth, and this should be explored as one potential mechanism underlying Pb’s neurodevelopmental effects.

Click here for additional data file.

Supplemental material, Supplemental_Figure_and_table for Trimester-Specific Associations of Prenatal Lead Exposure With Infant Cord Blood DNA Methylation at Birth by Christine A Rygiel, Dana C Dolinoy, Wei Perng, Tamara R Jones, Maritsa Solano, Howard Hu, Martha M Téllez-Rojo, Karen E Peterson and Jaclyn M Goodrich in Epigenetics Insights