Prenatal vitamin intake in first month of pregnancy and DNA methylation in cord blood and placenta in two prospective cohorts
Epigenetics & Chromatin volume 15, Article number: 28 (2022)
Prenatal vitamin use is recommended before and during pregnancies for normal fetal development. Prenatal vitamins do not have a standard formulation, but many contain calcium, folic acid, iodine, iron, omega-3 fatty acids, zinc, and vitamins A, B6, B12, and D, and usually they contain higher concentrations of folic acid and iron than regular multivitamins in the US Nutrient levels can impact epigenetic factors such as DNA methylation, but relationships between maternal prenatal vitamin use and DNA methylation have been relatively understudied. We examined use of prenatal vitamins in the first month of pregnancy in relation to cord blood and placenta DNA methylation in two prospective pregnancy cohorts: the Early Autism Risk Longitudinal Investigation (EARLI) and Markers of Autism Risk Learning Early Signs (MARBLES) studies.
In placenta, prenatal vitamin intake was marginally associated with −0.52% (95% CI −1.04, 0.01) lower mean array-wide DNA methylation in EARLI, and associated with −0.60% (−1.08, −0.13) lower mean array-wide DNA methylation in MARBLES. There was little consistency in the associations between prenatal vitamin intake and single DNA methylation site effect estimates across cohorts and tissues, with only a few overlapping sites with correlated effect estimates. However, the single DNA methylation sites with p-value < 0.01 (EARLI cord nCpGs = 4068, EARLI placenta nCpGs = 3647, MARBLES cord nCpGs = 4068, MARBLES placenta nCpGs = 9563) were consistently enriched in neuronal developmental pathways.
Together, our findings suggest that prenatal vitamin intake in the first month of pregnancy may be related to lower placental global DNA methylation and related to DNA methylation in brain-related pathways in both placenta and cord blood.
Vitamins and minerals are critical for normal fetal development. The World Health Organization recommends supplementation during pregnancy with iron, folic acid, vitamin A, calcium, and iodine . The American College of Obstetricians and Gynecologists additionally recommends supplementation with choline and vitamins B6, B12, C, and D . Prenatal vitamins do not have a standard formulation, but most contain calcium, iodine, omega-3 fatty acids, zinc, and vitamins A and D as well as more iron and B vitamins, and about twice as much folic acid compared to multivitamins [3,4,5,6]. In the US, prenatal vitamin use among pregnant people is estimated to be between 78 and 92%, with 55–60% reporting use in the first trimester [7,8,9,10]. In another study of the EARLI cohort, 59.7% reported prenatal vitamin use in the first month . Despite the use of these supplements, a recent study found that a significant number of pregnant people in the US still do not meet the recommended nutrient intake levels . Deficiencies in these nutrients are associated with multiple disorders including anemia and preeclampsia in the parent and impaired neurodevelopment, neural tube defects, and recurrent wheezing in the child [4, 12,13,14]. Understanding the molecular implications of these nutrients on fetal development is critical.
Nutrient levels can impact epigenetic factors such as DNA methylation. Relationships between prenatal vitamin use and DNA methylation have been relatively understudied, as most prior research focused on individual nutrients. For example, folic acid supplementation is associated with DNA methylation differences in cord blood at both differentially methylated regions  and differentially methylated positions in an epigenome-wide meta-analysis (n = 1988) . Another epigenome-wide association study examining vitamin D levels in two cohorts (n = 1416) found no association with cord blood DNA methylation . Maternal dietary intake of three dietary patterns was not associated with DNA methylation in placental tissue, implying that minor nutritional deficiencies were not associated with differences in DNA methylation (n = 573) . Conversely, a different study found that prenatal vitamin supplementation was associated with cord blood DNA methylation (n = 130) . A systematic review of randomized control trials of maternal micronutrient supplementation and DNA methylation in cord blood, blood spots, placental tissue, and buccal swabs found inconsistent results, but these studies varied in exposure timing and dose, types of micronutrients, sample size, and analytical methods . Additional studies of the associations between prenatal vitamin use and DNA methylation are needed.
Cord blood and placental tissues are useful sources of information to understand fetal development. Epigenetic factors such as DNA methylation play a role in normal and abnormal development . DNA methylation patterns from cord blood and placental tissues may be used as biomarkers of in utero exposures or to predict future health . DNA methylation has tissue-specific patterns, and tissues can have varying sensitivity to nutrient changes. Both cord blood and placental tissue develop early in gestation which represents an important window for early prenatal vitamin supplementation. The umbilical cord begins to form during the fourth week of gestation and blood is flowing in it by the fifth week, but the structure is not fully developed until the 12th week . Umbilical cord blood contains a higher level of hematopoietic stem cells compared to adult blood and includes differentiated cell types such as B cells, natural killer cells, T cells, monocytes, granulocytes, and nucleated red blood cells [24, 25]. The placenta is the site of gas, nutrient, and waste exchange between the parent and fetus and has metabolic and endocrine functions. The trophoblastic cells of the placenta are of fetal origin, but these are in close contact with the maternal decidua and blood vessels through the chorionic villi . The placenta starts to develop from the trophectoderm layer of the blastocyst following implantation in the maternal endometrium, around 6–7 days . Following zygote formation epigenetic reprogramming occurs, making early pregnancy a potentially important exposure window for DNA methylation impacts . The associations between prenatal vitamin use in the first month of pregnancy and DNA methylation in cord blood and placental tissue are still unclear. Here, we examined use of prenatal vitamin use during the first month of pregnancy and DNA methylation in cord blood and placenta in two prospective pregnancy cohorts. We analyzed array-wide mean DNA methylation, as well as associations at single DNA methylation sites and tested for enriched gene ontology processes.
Study sample descriptive statistics
For cord blood, variables of interest (prenatal vitamin use, maternal education, maternal age, genetic principal components, sex, gestational age, estimated cell proportions) were available for 113 samples (66.5% of total) in EARLI and 201 samples (82.7% of total) in MARBLES. For placenta, variables of interest were available for 88 samples (69.3% of total) in EARLI and 70 samples (77.8% of total) in MARBLES. Analyses were done separately in the four cohort/tissue groups: EARLI cord blood, EARLI placenta, MARBLES cord blood, and MARBLES placenta. In all cohort/tissue groups, DNA methylation principal components were associated with sex and estimated cell proportions (p < 0.001 in related to at least one of principal components one through three). Laboratory batch was associated with DNA methylation principal components in EARLI cord (p = 0.002 with principal component 2) and MARBLES cord (p = 0.002 with principal component 2, p < 0.001 with principal component 5) (Additional file 1: Figure S3). In cord blood, cell type composition of monocytes (p < 0.001) and nRBC (p < 0.001) differed between cohorts. In placenta, composition estimates of Hofbauer cells (p < 0.001) and nRBC (p = 0.003) differed between cohorts.
Among subjects with cord blood samples, use of prenatal vitamins in the first month of pregnancy was associated with maternal education in both cohorts. In EARLI, 69.0% of mothers who took prenatal vitamins in the first month had a college degree, compared to 42.9% of mothers who did not take prenatal vitamins. In MARBLES, 63.7% of mothers who took prenatal vitamins in the first month of pregnancy had at least a college degree, while 34.3% of mothers who did not take prenatal vitamins had a college degree (Table 1). Among the subset with placenta samples, this same trend was less pronounced in EARLI (64.2% mothers who took prenatal vitamins had college degree versus 51.4% among mothers who did not take prenatal vitamins), and present in MARBLES (69.7% taking prenatal vitamins in first month had college degree compared to 43.2% among those without prenatal vitamin use) (Table 2). By the latter half of pregnancy, the vast majority of mothers (> 75%) were taking prenatal vitamins (Additional file 1: Table S1). Estimated cell type did not differ by prenatal vitamin status, and maternal pre-pregnancy BMI did not differ but prenatal vitamin use.
Mean DNA methylation differences with prenatal vitamin use
Mean array-wide DNA methylation differed by prenatal vitamin use in the first month of pregnancy in placenta (Additional file 1: Figure S4). In EARLI, prenatal vitamin use was marginally associated with −0.52% (95% confidence interval: −1.04, 0.01) lower mean array-wide placenta DNA methylation. In MARBLES, prenatal vitamin use in the first month was associated with −0.60% (95% confidence interval: −1.08, −0.13) lower mean array-wide placenta DNA methylation. Differences by prenatal vitamin use were smaller in magnitude for island regions. In cord blood there were no differences in mean array-wide DNA methylation by prenatal vitamin use in either cohort.
Prenatal vitamin epigenome-wide association results
Across all analyses, no single CpG site met the more stringent significance threshold (p-value < 10–7). In EARLI placenta tissue 3442 sites had a nominal association (p-value < 0.01). Of these sites, 94.4% had lower DNA methylation with prenatal vitamin use in the first pregnancy month and these sites had an average of −3.9 percent lower DNA methylation (Fig. 1). The top site associated with prenatal vitamin use in placenta in EARLI was cg24700222 associated with COL24A1 gene (effect estimate = −7.3, p-value = 2.0 × 10–6) (Additional file 2: Table S2). In MARBLES placental tissue, 9216 sites had p-value < 0.01, and 96.4% had lower DNA methylation with prenatal vitamin use and these sites had an average of −4.2 percent lower DNA methylation. The top site associated with prenatal vitamin use in placenta in MARBLES was cg00711959 (effect estimate = −11.0, p-value = 3.3 × 10–6), which is located on chromosome 2 and is not annotated to a gene (Additional file 3: Table S3).
In cord blood, among nominally 3273 associated sites in EARLI, 80.5% had lower DNA methylation with prenatal vitamin use and these sites had an average of 1.8 percent lower DNA methylation. The top site associated with prenatal vitamin use in cord blood in EARLI was cg18452703 (effect estimate = −3.0, p-value = 2.7 × 10–6), located on chromosome 10 and not annotated to a gene. (Additional file 4: Table S4). In MARBLES cord blood, 58.3% of the 2,348 nominally associated sites had lower DNA methylation for prenatal vitamin use and these sites had an average 1.0 percent lower DNA methylation. The top site associated with prenatal vitamin use in cord blood in MARBLES was cg04551619 (effect estimate = 0.56, p-value = 1.8 × 10–5), which is located on chromosome 1 and is not annotated to a gene (Additional file 5: Table S5).
For the DNA methylation sites commonly measured across tissues and cohorts (n = 422,383CpGs), overall correlations among effect estimates for the associations of prenatal vitamin use and DNA methylation were weak to modest (Fig. 2). The strongest correlations in effect estimates were between MARBLES placenta (measured on EPIC) and EARLI placenta (measured on 450k) (r = 0.17, p-value < 2.2 × 10–16), and between EARLI cord (measured on 450k) and EARLI placenta (r = 0.17, p-value < 2.2 × 10–16). In contrast, correlation of cord blood effect estimates between the cohorts was non-existent (r = 0.02, p-value < 2.2 × 10–16), and correlation between cord blood and placenta in MARBLES was similarly weak (r = 0.06, p-value < 2.2 × 10–16).
We restricted to the DNA methylation sites associated (p-value < 0.01) with prenatal vitamin use and calculated overlap and correlation across tissues and cohorts. In EARLI cord blood, 4068 CpGs had p-value < 0.01, EARLI placenta had 3647 CpG sites, MARBLES cord had 4,025 CpGs, and MARBLES placenta had 9563 CpGs. Sites reaching this threshold were largely tissue and cohort specific (Fig. 3A). The largest levels of overlap mirrored the correlation patterns. Between EARLI and MARBLES placenta, 101 CpG sites had p-value < 0.01, and these 101 sites had correlation r = 0.80 (Fig. 3B). The 101 CpG sites had effect estimates in the same direction in MARBLES and EARLI placenta, except for one site, with 99 of them having negative effect estimates (Fig. 4). In contrast, comparing cord blood between MARBLES and EARLI there were 20 CpGs with p-value < 0.01 in both (Fig. 3B), with more modest correlation (r = 0.50) and mixed directions of effect (Fig. 4). Between cord blood and placenta in EARLI, 63 sites had p-value < 0.01 in both, with those CpG effect estimates having correlation r = 0.71. In MARBLES, 31 CpGs had p-value < 0.01 in both cord blood and placenta, with effect estimates having correlation r = 0.59 (Fig. 3B). EARLI had more consistency in direction of effect between cord blood and placenta than MARBLES (Additional file 1: Figure S5).
We compared single site results with a previous study on maternal plasma folate during pregnancy and cord blood methylation (data were accessed through dbGaP, accession number: phs001059.v1.p1.c1) . In that study, 7219 CpG sites had p-value < 0.01. Looking only at sites overlapping between the 450k and EPIC chips, the number of sites with p-value < 0.01 was 3365 in EARLI placenta, 3771 in EARLI cord, 4625 in MARBLES placenta, and 2342 in MARBLES cord. The level of overlap was minimal but more than expected at random between the maternal plasma folate study and our prenatal vitamin usage study, with 72 CpGs overlapping with EARLI placenta results (Fisher test p = 0.007), 131 CpGs with EARLI cord (Fisher test p < 0.001), 96 CpGs with MARBLES placenta (Fisher test p = 0.009), and 45 CpGs with MARBLES cord (Fisher test p = 0.10).
Gene ontologies enriched in prenatal vitamin use-related DNA methylation
CpG sites associated with prenatal vitamin use (p-value < 0.01) in the epigenome-wide regression models were used for gene ontology analysis (EARLI cord nCpGs = 4068, EARLI placenta nCpGs = 3647, MARBLES cord nCpGs = 4025, MARBLES placenta nCpGs = 9563). Across all four tissue and cohort combinations, the top pathways by rank sum were neuron and development-related pathways, such as nervous system development, neuron differentiation, and cell projection morphogenesis (Table 3, Additional file 6: Table S6). In cord blood, the top pathways by rank sum across the two cohorts included largely development-related pathways, such as sensory organ development, neuron differentiation, and embryo development. In MARBLES cord blood, no pathways had FDR < 1.0, and the top pathway was protein catabolic process in the vacuole (Additional file 7: Table S7). In EARLI cord blood, no pathways had FDR < 1.0, and the top pathway was nervous system development (Additional file 8: Table S8). In placenta, the top pathways by rank sum across the two cohorts were neuron and synaptic signaling pathways, including neurogenesis, and chemical synaptic transmission. In MARBLES placenta, the top pathway was neuron development (FDR = 0.14) (Additional file 9: Table S9). In EARLI placenta, the top pathways were regulation of bone remodeling (FDR = 0.48) and system development (FDR = 0.48) (Additional file 10: Table S10). Comparing the pathways with unadjusted p-value < 0.01, mirroring single site comparisons most pathways were unique to tissue/cohort (Additional file 1: Figure S6). A total of eight pathways, including neuron differentiation and chemical synaptic transmission, had p < 0.01 in all but MARBLES cord blood.
The top 1000 CpG sites associated with prenatal vitamin use in the first month of pregnancy were examined for enrichment in chromatin state marker signatures. Sites associated with prenatal vitamin use in EARLI cord blood were enriched in the repressed polycomb chromatin state, followed by weak repressed polycomb, bivalent enhancer, and bivalent/poised transcript start sequence markers (Additional file 1: Figure S7).
Replication testing with whole genome bisulfite sequencing (WGBS)
Data for WGBS were available on 63 samples [33 with prenatal vitamin use] in EARLI cord blood, 91 samples [39 with prenatal vitamin use] in MARBLES placenta, 45 samples [17 with prenatal vitamin use] in MARBLES cord on HiSeq 4000, and 42 samples [17 with prenatal vitamin use] in MARBLES cord on HiSeq X Ten. Differentially methylated regions (DMRs) with p-value < 0.05 in each are summarized in Additional file 11: Table S11. Overall, a majority of DMRs showed lower DNA methylation with prenatal vitamin use, and agreement in the direction of effect was strongest in placenta, where there were 803 array sites within 5 kb of DMRs identified in WGBS, and 66.6% had the same direction of effect across the two measures (Additional file 1: Table S12).
In two prospective pregnancy cohorts, we found that prenatal vitamin intake in the first month of pregnancy was related to lower average DNA methylation in placenta and cord blood. The magnitude of this association was strongest in placenta. We observed little consistency in the associations between prenatal vitamin intake and single DNA methylation site effect estimates across cohorts and tissues, with only a few overlapping sites with correlated effect estimates. However, the single DNA methylation sites associated with prenatal vitamin use were consistently enriched in neuron developmental pathways. Together these findings suggest that prenatal vitamin intake may be related to placental global DNA methylation and related to DNA methylation in brain-related pathways in both placenta and cord blood.
Previous research has examined DNA methylation and prenatal vitamin supplementation, or components of prenatal multivitamins, at different supplementation time points. Vitamin B12 is a component of prenatal vitamins, and prenatal deficiency leads to adverse infant hematologic and neurologic outcomes . In terms of DNA methylation, a cohort study of cord blood samples found that higher maternal serum vitamin B12 concentrations (collected at mean 10.6 weeks of gestation) were associated with lower global cord blood DNA methylation, using bisulfite pyrosequencing . Another study of 516 cord blood samples observed periconceptional vitamin B12 intake was not associated with long interspersed nuclear element-1 (LINE-1) DNA methylation, often used as a proxy for global methylation . Placentas from rats who were treated with excess folate and were vitamin B12 deficient had lower global DNA methylation levels compared to controls and rats with normal folate + B12 deficiency, but when they were given omega-3, the DNA methylation returned to the control group level . In humans, mothers with high folate and low vitamin B12 during early pregnancy had significantly lower cord blood DNA methylation as measured by LC–MS/MS . These previous studies were heterogeneous in species, sample type, and DNA methylation measurement methods. The conflicting results highlight the need for more comprehensive analyses of the associations between early pregnancy vitamin supplementation and DNA methylation, considering the intake of multiple vitamins.
Some prenatal vitamins also contain omega-3 fatty acids. The omega-3 fatty acid, docosahexaenoic acid (DHA) and polyunsaturated fatty acids are essential during pregnancy for fetal development . A study of DHA, given at 20+ weeks of gestation, found no difference in global DNA methylation of blood spots and leukocytes between the children from the treatment and control groups . They did identify 21 differentially methylated regions with most regions showing lower methylation levels in the treatment group compared to the control. The direction of this finding matches ours where a majority of our differentially methylated sites (p-value < 0.01) across both cohorts and tissue types associated lower methylation with prenatal vitamin use. A randomized control trial of polyunsaturated fatty acid supplementation from weeks 18–22 to birth observed no difference in the cord blood global DNA methylation between the non-smoking treatment and control groups . The review by Andraos et al. (2018) summarizes the DNA methylation results from randomized control trials of nutritional supplementations during pregnancy, finding overall that micronutrient supplementation does not substantially affect offspring DNA methylation in cord blood, blood spots, placental tissue, and buccal swabs, though lack of standardized methods complicate comparing results of studies .
Folate and folic acid are components of prenatal vitamins. In a randomized control trial of maternal folic acid supplementation in the second and third trimesters, supplementation was associated with lower DNA methylation in cord blood at LINE-1 compared to the control group [37, 38]. A meta-analysis of two cohorts of 1,988 infants identified 443 DNA methylation sites in cord blood associated (FDR < 0.05) with maternal plasma folate levels measured during pregnancy (median 18 weeks gestation, and median 12.9 weeks gestation in the two cohorts) . Of these sites, 416 (94%) had lower DNA methylation with higher plasma folate levels, suggesting a similar direction of association to what we observed. Using gene ontologies, prenatal vitamin-associated DNA methylation sites were enriched for developmental and neuronal pathways in both our cohorts and tissues. The meta-analysis also observed an enrichment of developmental and neurodevelopmental pathways . Although we did not observe correlation of associations at individual DNA methylation sites to this previous study, our findings were consistent with respect to overall trends of lower DNA methylation and pathway enrichment.
Prenatal vitamins overlap with multivitamins in their nutrient components, although prenatal vitamins typically contain more folic acid than multivitamins (about twice as much) . Previously in the EARLI cohort, multivitamin usage in the 3 months before pregnancy in mothers without the MTHFR allele variant was associated with higher global cord blood DNA methylation . However, similar to our current findings in cord blood, the previous EARLI study observed no association between prenatal vitamin use in the 3 months prior to pregnancy and global cord blood DNA methylation (p < 0.05). These analyses differed in their exposure window (prior paper: 3 months prior to pregnancy versus current paper: first month of pregnancy).
Differences in DNA methylation have been associated with other early-life exposures to environmental factors, in addition to prenatal vitamin use [39,40,41]. In this study, we found larger differences in DNA methylation and more consistency between cohorts in placenta compared to cord blood. This finding is likely due to known DNA methylation characteristics of placenta, including global hypomethylation relative to somatic tissues and the presence of partially methylated domains . In placenta, methylation of gene bodies is predictive of active expression and exposure to home/garden pesticides was previously shown to be associated with higher global DNA methylation in MARBLES placenta samples [43, 44]. Combined, these findings suggest that the lower global DNA methylation levels associated with prenatal vitamin use may be reflecting a more quiescent genome with reduced activation of environmentally responsive genes that may be related to the known effects of nutrients in reducing oxidative stress . We observed surprisingly low correlation in cord blood effect estimates between the two cohorts. This could suggest a lack of robust single site effects, but batch effects may have contributed to reduced correlation. Within the EARLI study, we observed higher correlation between the placenta and cord blood DNA methylation associations with prenatal vitamin intake, than we did in the MARBLES study. The EARLI placenta and cord blood DNA methylation measures were conducted at the same time and on the same laboratory plates. In MARBLES, the cord blood and placenta samples were run at different times and on different laboratory plates. Cell type proportions are a large driver of variability in DNA methylation measures in most tissues. Environmental exposures can even be highly associated with cell type proportions , and cell type proportions can mediate the effects of environmental exposures on DNA methylation . In our study, cell composition did not differ by prenatal vitamin status, but we did observe differences in the cell composition of tissues by cohort.
The prior studies reviewed above and the current study have great heterogeneity with respect to the type, method, and timing of prenatal vitamin-related exposure measures. As described, the associations vary, but more studies have reported lower DNA methylation with pregnancy vitamin use. Most environmental epigenetic studies have used the Illumina BeadArrays for DNA methylation measures, though the tissue measured at birth varies. Replication testing or meta-analyses are necessary to determine the reproducibility of these findings and few studies to date have included more than one cohort. These steps will require harmonization of exposure and DNA methylation measures across studies. Once reproducibility and specificity are determined, these DNA methylation patterns may be used as biomarkers of exposure or may predict future health .
There were a number of strengths to this study. We analyzed two prospective pregnancy cohorts and two tissue types with 201 total samples from EARLI and 271 total samples from MARBLES. The prospective design minimized recall bias of prenatal vitamin use. Furthermore, we had consistent exposure measures including using data from questionnaires, consistent cleaning of the vitamin questionnaires, and consistent exposure timing during pregnancy. The results were based on a diverse subset of participants from five study sites in the US. We employed rigorous methods to preprocess and analyze the DNA methylation data. Our analysis across the two cohorts and two tissue types showed similar gene ontology pathways. Finally, we compared the array-based DNA methylation single site results to the WGBS results.
There were several limitations to our study. First, the samples used for DNA methylation were from one time point at birth, so long-term differences in DNA methylation were not assessed. We examined prenatal vitamin use as a yes/no response for any use in the first pregnancy month; effects may differ by frequency of intake, nutrient composition, and nutrient dose. We did not account for differences in underlying nutritional deficiencies which may have affected our results. Future studies could assess nutritional intake in addition to prenatal vitamins. Exposure to other environmental chemicals may lower the DNA methylation levels of adults, children, and infants [30, 48,49,50,51,52]. It is possible that differences in chemical exposures may account for or mask some of the differences across DNA methylation related to prenatal vitamin use. While there was little positional overlap in effects between array-based DNA methylation and WGBS, we observed similar directions of effect across measurement methods in placenta. This finding may have been due to a greater overlap in these samples between measured on array compared to WGBS. Data on method of delivery, known to affect DNA methylation , were incomplete in EARLI, so we did not adjust for that information in our regression analyses for both cohorts. Additional caution should be applied when interpreting our findings, since prenatal vitamin use was associated with higher educational attainment, and there may be unmeasured confounding with such factors as improved prenatal care, improved diet, and other factors. Studies with available measures in these areas will be highly valuable. Future studies could examine a larger cohort, compare these results to a general population cohort, or conduct a meta-analysis of multiple prenatal vitamin studies. Future studies may also consider examining the placenta, to explore whether it is more sensitive to effects of supplementation, as our study suggests.
We found that prenatal vitamin use in the first month of pregnancy was associated with lower DNA methylation, particularly in the placenta. Prenatal vitamin use is recommended before and during pregnancies for normal fetal development. Given its importance, additional research is needed to understand the underlying biological mechanisms of development. By demonstrating an association between prenatal vitamin intake and DNA methylation at birth, we lay the foundation for DNA methylation as a biomarker of prenatal vitamin exposure. While promising, these findings also highlight the need for larger studies in this area with standardized, well measured, and longitudinal prenatal vitamin exposure measures for testing with epigenetic marks and fetal and postnatal growth and heath.
The Early Autism Risk Longitudinal Investigation (EARLI) and Markers of Autism Risk Learning Early Signs (MARBLES) studies are enriched risk prospective pregnancy cohorts studying autism etiology [54, 55]. The EARLI study was reviewed and approved by Human Subjects Institutional Review Boards (IRBs) from each of the four study sites (Johns Hopkins University, Drexel University, University of California Davis, and Kaiser Permanente Northern California). The MARBLES protocol was reviewed and approved by the Human Subjects IRB from University of California Davis. Secondary data analysis for this manuscript were approved by the Human Subjects IRB for the University of Michigan. Both studies recruited mothers of children with clinically confirmed ASD who were early in a subsequent pregnancy or were trying to become pregnant. In EARLI there were 232 mothers with a subsequent sibling born through this study between November 2009 and March 2012. In MARBLES there were 389 enrolled mothers that gave birth to 425 subsequent siblings between December 1, 2006 and July 1, 2016.
Covariate and exposure assessment
Demographics, behaviors, and medical history were all collected longitudinally via maternal self-report questionnaire. In these questionnaires, mothers were asked if they used prenatal vitamins for each month of pregnancy (yes/no). Data for the first month of pregnancy and for 3 months prior to pregnancy were collected at study enrollment.
Sample collection and processing
In EARLI, biospecimens including cord blood and placenta, were collected and archived for 213 births. Full thickness placental tissue from a central cotyledon was collected. Sterile punch biopsy forceps were used to extract placental samples from the maternal and fetal sides. Whole cord blood was also collected at delivery. Samples were transported to the Johns Hopkins Biological Repository (JHBR) for aliquoting and archiving (-80ºC). Placental DNA was extracted with the DNeasy Tissue Kit (Qiagen), and cord blood DNA was extracted using the DNA Midi kit (Qiagen, Valencia, CA). DNA was quantified using the Nanodrop (ThermoFisher Waltham, MA) and normalized DNA aliquots were sent to the Center for Inherited Disease Research (Johns Hopkins University). DNA samples were bisulfite treated and cleaned using the EZ DNA methylation gold kit (Zymo Research, Irvine, CA) according to manufacturer’s instructions. DNA was plated randomly and was assayed on the Infinium HumanMethylation450 BeadChip (Illumina, San Diego, CA) . Methylation control gradients and between-plate repeated tissue controls were used.
In MARBLES, placental tissues and cord blood were collected at delivery and immediately processed and frozen. The MARBLES study used orientation to the umbilical cord to ensure that all placenta samples were isolated from the chorionic villus from the fetal side of the placenta. Placental and cord blood samples were stored at −80ºC in the UC Davis repository. Cord blood and placenta samples were processed for methylation measures. Placenta DNA was extracted with Gentra Puregene kit (Qiagen) and cord blood DNA was extracted using the DNA Midi kit (Qiagen, Valencia, CA). Samples were bisulfite treated and cleaned using the EZ DNA methylation gold kit (Zymo Research, Irvine, CA). DNA was plated randomly and assayed on the Infinium HumanMethylationEPIC BeadChip (Illumina, San Diego, CA) at the Johns Hopkins SNP Center, a shared lab and informatics operation with the Center for Inherited Disease Research (Johns Hopkins University). DNA methylation control gradients and between-plate repeated tissue controls were used.
DNA methylation processing
For all methylation samples, we used the minfi library (version 1.30.0) in R (version 3.6) to process raw Illumina image files into noob background corrected methylation values [57, 58]. In EARLI, cord blood and placenta samples were run on the 450k array together in two batches, and thus preprocessed together. Samples from multiple births (cord blood n = 2 samples, placenta n = 6 samples), as well as samples with discordant DNA methylation predicted sex and observed infant sex were removed (cord blood n = 3, placenta n = 1). Probes with failed detection P-value (> 0.01) in > 5% of samples were removed (n = 661), as were probes documented as cross-reactive (n = 29,153) . Y-chromosome probes (n = 48) were dropped from analysis. There were 170 EARLI cord blood samples, and 127 EARLI placenta samples with 455,650 probes that passed DNA methylation quality control.
In MARBLES, placenta and cord blood samples were run on the EPIC array at different times and preprocessed separately. First, we dropped cord blood samples from multiple births (cord blood n = 8 samples). Samples that had mismatched predicted sex were dropped (cord n = 3). For siblings not from multiple births, all but one sibling was dropped (cord n = 13). Probes were dropped if they had detection-p (p > 0.01) failure in greater than 5% of samples (n = 4630). Cross-reactive probes (n = 42,967) and Y chromosome probes were dropped from analysis (n = 379) . There were 243 MARBLES cord blood samples with 817,883 probes that passed DNA methylation quality control. Second, no placenta samples had mismatched predicted sex. There were no samples from multiple births, and all but one sample from siblings were dropped (placenta n = 2). Probes that failed detection-p in > 5% of samples (n = 1,699), cross-reactive probes (n = 43,068), and remaining Y-chromosome probes were dropped from analysis (n = 84). There were 90 MARBLES placenta samples with 821,008 probes that passed quality control. Sample exclusion is summarized in Additional file 1: Figure S1 and CpG probe exclusion is summarized in Additional file 1: Figure S2.
In cord blood samples, cell type (CD8+ T-cell, CD4+ T-cell, natural killer cell, B-cell, monocyte, granulocyte, and nucleated red blood cell) proportions were estimated using a combined reference panel with the IDOL method . In placenta samples, EpiDISH  was used to predict proportions of placenta cell types using a reference panel from the planet package: trophoblasts, stromal cells, Hofbauer cells, endothelial cells, nucleated red blood cells, and syncytiotrophoblasts . Mean DNA methylation per person was calculated as the mean across all probes . Mean DNA methylation restricted to probes in genomic regions (CpG island, shore, shelf, or open sea) were also computed. We used Illumina’s annotation of CpG sites to assign genomic regions (CpG island, CpG shore, CpG shelf, open sea) [64, 65].
Genetics data processing
In EARLI, genetic data were measured using the Omni5 + exome array (Illumina) at the John Hopkins University Center of Inherited Disease Research (CIDR). Data on 4.6 million single nucleotide polymorphisms (SNPs) were generated for 841 EARLI family biosamples (including maternal, paternal, proband, and infant samples) from 254 families and 18 HapMap control samples. Samples were processed together, but only data from infants with cord blood or placenta methylation were used. No samples had missing genotypes at > 3% of probes, or excess heterozygosity or homozygosity [4 standard deviations]. Probes were removed if they had technical problems flagged by CIDR or missing genomic location information. Single nucleotide polymorphisms (SNPs) with minor allele frequencies > 5% were removed if they had a missingness rate > 5%, and SNPs with minor allele frequency < 5% were removed if they had a missingness rate > 1%. There were 2.5 million clean SNPs for 827 samples, which were merged with the 1000 genomes project (1000GP, version 5) data  and principal components for genetic ancestry were computed.
In MARBLES, SNPs on 643 infant and mother samples from 234 families were genotyped using the Illumina Mega array at the John Hopkins University Center of Inherited Disease Research (CIDR). Maternal and infant samples were processed together, but only data from infants with cord blood or placenta methylation measures were used. We again applied stringent quality control criteria  to the raw 1.75 million genotypes to remove low quality SNPs and samples. Our criteria include removal of samples with call rate < 98%, sex discrepancy, and relatedness (pi-hat < 0.18) to non-familial samples. We also filtered SNPs with call rates < 95%, excess hetero- or homozygosity, and minor allele frequency (MAF) < 5%. After quality control, 620 samples and 758 thousand SNPs remained. Principal components were calculated on genotype data, and these principal components were used to adjust for genetic ancestry in models.
Study sample descriptive statistics were calculated for each of the four cohort/tissue groups. For continuous covariates (maternal age at delivery, gestational age, estimated cell proportions), we calculated mean and standard deviation. For categorical covariates (maternal education, infant sex, infant race/ethnicity), we provided number and frequency. We tested for differences in covariates and prenatal vitamin status using t-tests for continuous covariates, and Chi-square tests for categorical covariates. Because of systematic differences in measures, in terms of platform and timing/batches, we chose not to pool the EARLI and MARBLES samples and instead to conduct analyses separately within each sample. We conducted principal component analysis on the DNA methylation data and evaluated principal component associations with covariates. These associations were visualized with a heatmap of p-values.
In multivariable linear regression analyses, first, we examined array-wide mean DNA methylation differences by prenatal vitamin intake in the first month of pregnancy. Regression models were adjusted for infant sex, maternal age, gestational age, maternal education, and genetic ancestry principal components. To provide flexibility in conceptualization of cell type proportions we used estimated cell proportions as terms in regression models. Since cell composition estimates sum up to 100%, to avoid collinearity issues in models, we did not use all predicted cell types in models. For placenta, syncytiotrophoblast and Hofbauer proportions were used, while in cord blood granulocyte and nucleated red blood cell proportions were used. Batch was also adjusted for as a covariate. In the EARLI cohort, samples were measured in two runs, and an indicator variable for run was used for adjustment. MARBLES cord blood was run in a single batch, and sample plate (of which there were three) was used as a covariate for adjustment. MARBLES placenta was run in a single batch on a single plate, and no batch covariate was used. Since smoking has known impacts on DNA methylation, and smoking rates in both cohorts are very low, we excluded mothers with smoking during pregnancy in this analysis (number dropped: MARBLES cord n = 10, MARBLES placenta n = 4, EARLI cord n = 6, EARLI cord n = 6). We visualized regression coefficients and 95% confidence intervals using forest plots.
Next, we performed epigenome-wide association analyses by examining single CpG site differential DNA methylation. We fit parallel linear models for each probe. Models were again adjusted for infant sex, maternal age, gestational age, maternal education, genetic ancestry PCs, and estimated cell proportions. Regression and empirical Bayes standard error moderation were performed using the limma package . We visualized findings using volcano plots of effect estimates and −log10 (p-values). For sites reaching a nominal p-value threshold (p < 0.05), we calculated the proportion of sites that had higher DNA methylation with prenatal vitamin intake and the proportion of sites with lower DNA methylation. To compare pairwise results across cohort and tissues, we examined Pearson correlations of effect estimates from these regression models, across all sites in common between the 450k and EPIC methylation arrays. We also focused on CpG sites that had p-value < 0.01 in multiple cohort/tissues, examining the overlap of such sites with an upset plot and the Pearson correlation of overlapping sites. For sites prioritized in multiple cohorts/tissues, we also used scatter plots to visualize the effect estimates.
We tested enrichment for gene ontology biological processes using the missMethyl package . As input to missMethyl, CpG sites with p-value < 0.01 in the epigenome-wide regression models were used. We ranked the gene ontologies by significance, then computed a rank sum by adding the ranks across the four cohort/tissue groups. In addition, we tested for enrichment of chromatin state types using eFORGE 2.0 . The top 1000 CpG sites for each cohort/tissue analysis was input into the eFORGE site, with appropriate array platform chosen (450k for EARLI, EPIC for MARBLES), and Consolidated Roadmap Epigenomics—All 15-state marks and 1 kb window proximity options, with other options set at defaults.
In EARLI, whole genome bisulfite sequencing (WGBS) data were available on 63 cord blood samples sequenced on the HiSeq X [51 overlapping with array methylation data]. Sample processing and WGBS quality control and alignment for cord blood samples  and placenta samples [72, 73] have been previously discussed. In MARBLES, WGBS data were available for 91 placenta samples sequenced on HiSeq X Ten (89 overlapping with the array methylation data), 45 cord blood samples sequenced on HiSeq 4000 [30 overlapping with array methylation data], and 42 cord blood samples sequenced on the HiSeq X [35 overlapping with array methylation data].
Raw sequencing reads were preprocessed, mapped to human genome, and converted to CpG methylation count matrices with CpG_Me using default parameters [74,75,76]. Reads were trimmed for adapters and methylation bias, aligned to the reference genome, and filtered for PCR duplicates. Methylation counts at all sites were extracted to Bismark cytosine methylation reports. The CpG_Me workflow incorporates Trim Galore, Bismark, Bowtie2, SAMtools, and MultiQC [75, 77,78,79,80].
Differentially methylated regions (DMRs) were identified between prenatal vitamin intake during the first month of pregnancy with adjustment for sex and 10 permutation tests using DMRichR . The DMR analysis utilizes a smoothing and weighting algorithm to weight regions with high coverage and low variation. Permutation testing was performed on pooled null distribution to generate empirical p-values as significant DMRs. The DMRichR pipeline utilized dmrseq and bsseq algorithms [81, 82]. To evaluate consistency with array results, we identified array probes within 5 kb of the DMR and examined concordance in estimated directions of effect of those CpG probes and the DMR.
Availability of data and materials
Data are available through the National Database for Autism Research (NDAR) accession numbers for EARLI (1600) and MARBLES (2462). Code files are available through GitHub (https://github.com/bakulskilab).
Early Autism Risk Longitudinal Investigation
Differentially methylated regions
Long Interspersed Nuclear Element-1
Markers of Autism Risk Learning Early Signs
World Health Organization. Essential Nutrition Actions: improving maternal, newborn, infant and young child health and nutrition. Geneva: World Health Organization; 2013.
The American College of Obstetricians and Gynecologists. Nutrition During Pregnancy [Internet]. https://www.acog.org/womens-health/faqs/nutrition-during-pregnancy?utm_source=redirect&utm_medium=web&utm_campaign=int#:~:text=When you are pregnant you,first 12 weeks of pregnancy. Accessed 9 Mar 2021
Ottney A, Lebeau L. A “secret shopper” survey of community pharmacist prenatal vitamin recommendations. Jam Pharm Assoc. 2021. https://doi.org/10.1016/j.japh.2021.07.015.
Duerbeck NB, Dowling DD, Duerbeck JM. Prenatal vitamins: what is in the bottle? Obstet Gynecol Surv. 2014;69(12):777.
Schmidt RJ, Hansen RL, Hartiala J, Allayee H, Schmidt LC, Tancredi DJ, et al. Prenatal vitamins, one-carbon metabolism gene variants, and risk for autism. Epidemiology. 2011;22(4):476–85.
Brieger KK, Bakulski KM, Pearce CL, Baylin A, Dou JF, Feinberg JI, et al. The association of prenatal vitamins and folic acid supplement intake with odds of autism spectrum disorder in a high-risk sibling cohort, the early autism risk longitudinal investigation (EARLI). J Autism Dev Disord. 2021. https://doi.org/10.1007/s10803-021-05110-9.
Branum AM, Bailey R, Singer BJ. Dietary supplement use and folate status during pregnancy in the United States. J Nutr. 2013;143(4):486–92.
Aronsson CA, Vehik K, Yang J, Uusitalo U, Hay K, Joslowski G, et al. Use of dietary supplements in pregnant women in relation to sociodemographic factors—a report from The Environmental Determinants of Diabetes in the Young (TEDDY) study. Public Health Nutr. 2013;16(8):1390–402.
Sullivan KM, Ford ES, Azrak MF, Mokdad AH. Multivitamin use in pregnant and nonpregnant women: results from the behavioral risk factor surveillance system. Public Health Rep. 2009;124(3):384–90.
Saldanha LG, Dwyer JT, Andrews KW, Brown LL, Costello RB, Ershow AG, et al. Is nutrient content and other label information for prescription prenatal supplements different from nonprescription products? J Acad Nutr Diet. 2017;117(9):1429–36.
Bailey RL, Pac SG, Fulgoni VL III, Reidy KC, Catalano PM. Estimation of total usual dietary intakes of pregnant women in the United States. JAMA Netw Open. 2019;2(6):e195967–e195967.
Sfakianaki AK. Prenatal vitamins: a review of the literature on benefits and risks of various nutrient supplements. Formulary. 2013;48(2):77–82.
Patel A, Lee SY, Stagnaro-Green A, MacKay D, Wong AW, Pearce EN. Iodine content of the best-selling United States adult and prenatal multivitamin preparations. Thyroid. 2018;29(1):124–7.
Vahdaninia M, Mackenzie H, Helps S, Dean T. Prenatal intake of vitamins and allergic outcomes in the offspring: a systematic review and meta-analysis. J Allergy Clin Immunol Pract. 2017;5(3):771-778.e5.
Hoyo C, Murtha AP, Schildkraut JM, Jirtle RL, Demark-Wahnefried W, Forman MR, et al. Methylation variation at IGF2 differentially methylated regions and maternal folic acid use before and during pregnancy. Epigenetics. 2011;6(7):928–36.
Joubert BR, den Dekker HT, Felix JF, Bohlin J, Ligthart S, Beckett E, et al. Maternal plasma folate impacts differential DNA methylation in an epigenome-wide meta-analysis of newborns. Nat Commun. 2016;7(1):10577.
Suderman M, Stene LC, Bohlin J, Page CM, Holvik K, Parr CL, et al. 25-Hydroxyvitamin D in pregnancy and genome wide cord blood DNA methylation in two pregnancy cohorts (MoBa and ALSPAC). J Steroid Biochem Mol Biol. 2016;159:102–9.
Lecorguillé M, Charles MA, Lepeule J, Lioret S, de Lauzon-Guillain B, Forhan A, et al. Association between dietary patterns reflecting one-carbon metabolism nutrients intake before pregnancy and placental DNA methylation. Epigenetics. 2021;31:1–16.
Bakulski KM, Dou JF, Feinberg JI, Brieger KK, Croen LA, Hertz-Picciotto I, et al. Prenatal multivitamin use and MTHFR genotype are associated with newborn cord blood DNA methylation. Int J Environ Res Public Health. 2020;17:9190.
Andraos S, de Seymour JV, O’Sullivan JM, Kussmann M. The impact of nutritional interventions in pregnant women on dna methylation patterns of the offspring: a systematic review. Mol Nutr Food Res. 2018;62(24):1800034.
Robertson KD. DNA methylation and human disease. Nat Rev Genet. 2005;6:597–610.
James P, Sajjadi S, Tomar AS, Saffari A, Fall CHD, Prentice AM, et al. Candidate genes linking maternal nutrient exposure to offspring health via DNA methylation: a review of existing evidence in humans with specific focus on one-carbon metabolism. Int J Epidemiol. 2018;47(6):1910–37.
Spurway J, Logan P, Pak S. The development, structure and blood flow within the umbilical cord with particular reference to the venous system. Australas J Ultrasound Med. 2012;15(3):97–102.
Broxmeyer HE, Douglas GW, Hangoc G, Cooper S, Bard J, English D, et al. Human umbilical cord blood as a potential source of transplantable hematopoietic stem/progenitor cells. Proc Natl Acad Sci. 1989;86(10):3828–32.
Bakulski KM, Feinberg JI, Andrews SV, Yang J, McKenney Brown SL., S, et al. DNA methylation of cord blood cell types: Applications for mixed cell birth studies. Epigenetics. 2016;11(5):354–62.
Gude NM, Roberts CT, Kalionis B, King RG. Growth and function of the normal human placenta. Thromb Res. 2004;114(5):397–407.
Kaufmann P, Frank HG. Placental Development. In: Polin RA, Fox WW, Abman SHBTF, editors. Fetal and Neonatal Physiology. Philadelphia: WB Saunders; 2004.
Faulk C, Dolinoy DC. Timing is everything: the when and how of environmentally induced changes in the epigenome of animals. Epigenetics. 2011;6(7):791–7.
Finkelstein JL, Layden AJ, Stover PJ. Vitamin B-12 and perinatal health. Adv Nutr. 2015;6(5):552–63.
McKay JA, Groom A, Potter C, Coneyworth LJ, Ford D, Mathers JC, et al. Genetic and non-genetic influences during pregnancy on infant global and site specific DNA methylation: role for folate gene variants and vitamin B12. PloS one. 2012;7(3):e33290–e33290.
Boeke CE, Baccarelli A, Kleinman KP, Burris HH, Litonjua AA, Rifas-Shiman SL, et al. Gestational intake of methyl donors and global LINE-1 DNA methylation in maternal and cord blood: Prospective results from a folate-replete population. Epigenetics. 2012;7(3):253–60.
Kulkarni A, Dangat K, Kale A, Sable P, Chavan-Gautam P, Joshi S. Effects of altered maternal folic acid, vitamin B12 and docosahexaenoic acid on placental global DNA methylation patterns in Wistar rats. PLoS ONE. 2011;6(3): e17706.
Plumptre L, Tammen SA, Sohn KJ, Masih SP, Visentin CE, Aufreiter S, et al. Maternal and cord blood folate concentrations are inversely associated with fetal DNA Hydroxymethylation, but not DNA Methylation, in a cohort of pregnant Canadian women. J Nutr. 2020;150(2):202–11.
Coletta JM, Bell SJ, Roman AS. Omega-3 fatty acids and pregnancy. Rev Obstet Gynecol. 2010;3(4):163–71.
van Dijk SJ, Zhou J, Peters TJ, Buckley M, Sutcliffe B, Oytam Y, et al. Effect of prenatal DHA supplementation on the infant epigenome: results from a randomized controlled trial. Clin Epigenetics. 2016;8(1):114.
Lee HS, Barraza-Villarreal A, Hernandez-Vargas H, Sly PD, Biessy C, Ramakrishnan U, et al. Modulation of DNA methylation states and infant immune system by dietary supplementation with ω-3 PUFA during pregnancy in an intervention study. Am J Clin Nutr. 2013;98(2):480–7.
Caffrey A, Irwin RE, McNulty H, Strain JJ, Lees-Murdock DJ, McNulty BA, et al. Gene-specific DNA methylation in newborns in response to folic acid supplementation during the second and third trimesters of pregnancy: epigenetic analysis from a randomized controlled trial. Am J Clin Nutr. 2018;107(4):566–75.
Yang AS, Estécio MRH, Doshi K, Kondo Y, Tajara EH, Issa JJ. A simple method for estimating global DNA methylation using bisulfite PCR of repetitive DNA elements. Nucleic Acids Res. 2004;32(3):e38–e38.
Sen A, Heredia N, Senut MC, Hess M, Land S, Qu W, et al. Early life lead exposure causes gender-specific changes in the DNA methylation profile of DNA extracted from dried blood spots. Epigenomics. 2015;7(3):379–93.
Demetriou CA, van Veldhoven K, Relton C, Stringhini S, Kyriacou K, Vineis P. Biological embedding of early-life exposures and disease risk in humans: a role for DNA methylation. Eur J Clin Invest. 2015;45(3):303–32.
Mitchell C, Schneper LM, Notterman DA. DNA methylation, early life environment, and health outcomes. Pediatr Res. 2016;79(1):212–9.
Schroeder DI, Blair JD, Lott P, Yu HOK, Hong D, Crary F, et al. The human placenta methylome. Proc Natl Acad Sci USA. 2013;110(15):6037–42.
Schroeder DI, Jayashankar K, Douglas KC, Thirkill TL, York D, Dickinson PJ, et al. Early developmental and evolutionary origins of gene body DNA methylation patterns in Mammalian placentas. PLoS Genet. 2015;11(8):e1005442–e1005442.
Schmidt RJ, Schroeder DI, Crary-Dooley FK, Barkoski JM, Tancredi DJ, Walker CK, et al. Self-reported pregnancy exposures and placental DNA methylation in the MARBLES prospective autism sibling study. Environ Epigenetics. 2016;2(4):dvw024.
Rodríguez-Cano AM, Calzada-Mendoza CC, Estrada-Gutierrez G, Mendoza-Ortega JA, Perichart-Perera O. Nutrients, mitochondrial function, and perinatal health. Nutrients. 2020;12(7):2166.
Middleton LYM, Nguyen VK, Dou J, Park SK, Colacino JA, Bakulski KM. Environmental chemical-wide associations with immune biomarkers in the US: A cross-sectional analysis. medRxiv [Internet]. 2022; https://www.medrxiv.org/content/early/2022/04/25/2022.03.22.22272789
Campbell KA, Colacino JA, Park SK, Bakulski KM. Cell types in environmental epigenetic studies: Biological and epidemiological frameworks. Curr Environ Health Rep. 2020;7(3):185–97.
Rusiecki JA, Andrea B, Valentina B, Letizia T, Moore LE, Bonefeld-Jorgensen EC. Global DNA hypomethylation is associated with high serum-persistent organic pollutants in Greenlandic Inuit. Environ Health Perspect. 2008;116(11):1547–52.
Keon-Yeop K, Dong-Sun K, Sung-Kook L, In-Kyu L, Jung-Ho K, Yoon-Seok C, et al. Association of low-dose exposure to persistent organic pollutants with global DNA hypomethylation in healthy Koreans. Environ Health Perspect. 2010;118(3):370–4.
Bollati V, Baccarelli A, Hou L, Bonzini M, Fustinoni S, Cavallo D, et al. Changes in DNA methylation patterns in subjects exposed to low-dose benzene. Cancer Res. 2007;67(3):876–80.
Baccarelli A, Wright RO, Bollati V, Tarantini L, Litonjua AA, Suh HH, et al. Rapid DNA methylation changes after exposure to traffic particles. Am J Respir Crit Care Med. 2009;179(7):572–8.
Breton CV, Byun HM, Wenten M, Pan F, Yang A, Gilliland FD. Prenatal tobacco smoke exposure affects global and gene-specific DNA methylation. Am J Respir Crit Care Med. 2009;180(5):462–7.
Schlinzig T, Johansson S, Gunnar A, Ekström TJ, Norman M. Epigenetic modulation at birth—altered DNA-methylation in white blood cells after Caesarean section. Acta Paediatr. 2009;98(7):1096–9.
Newschaffer CJ, Croen LA, Fallin MD, Hertz-Picciotto I, Nguyen DV, Lee NL, et al. Infant siblings and the investigation of autism risk factors. J Neurodev Disord. 2012;4(1):7.
Hertz-Picciotto I, Schmidt RJ, Walker CK, Bennett DH, Oliver M, Shedd-Wise KM, et al. A prospective study of environmental exposures and early biomarkers in autism spectrum disorder: Design, protocols, and preliminary data from the MARBLES study. Environ Health Perspect. 2018;126(11):117004.
Bibikova M, Barnes B, Tsan C, Ho V, Klotzle B, Le JM, et al. High density DNA methylation array with single CpG site resolution. Genomics. 2011;98(4):288–95.
Aryee MJ, Jaffe AE, Corrada-Bravo H, Ladd-Acosta C, Feinberg AP, Hansen KD, et al. Minfi: A flexible and comprehensive bioconductor package for the analysis of infinium DNA methylation microarrays. Bioinformatics. 2014;30(10):1363–9.
Triche TJ, 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.
Chen Y, Lemire M, Choufani S, Butcher DT, Grafodatskaya D, Zanke BW, et al. Discovery of cross-reactive probes and polymorphic CpGs in the Illumina Infinium HumanMethylation450 microarray. Epigenetics. 2013;8(2):203–9.
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.
Gervin K, Salas LA, Bakulski KM, van Zelm MC, Koestler DC, Wiencke JK, et al. Systematic evaluation and validation of reference and library selection methods for deconvolution of cord blood DNA methylation data. Clin Epigenetics. 2019;11(1):125.
Teschendorff AE, Breeze CE, Zheng SC, Beck S. A comparison of reference-based algorithms for correcting cell-type heterogeneity in Epigenome-Wide Association Studies. BMC Bioinformatics. 2017;18(1):105.
Yuan V, Hui D, Yin Y, Penaherrera MS, Beristain AG, Robinson WP. Cell-specific characterization of the placental methylome. BMC Genomics. 2021;22(1):6.
Hansen KD. IlluminaHumanMethylation450kanno.ilmn12.hg19: Annotation for Illumina’s 450k methylation arrays. 2016.
Hansen KD. IlluminaHumanMethylationEPICanno.ilm10b2.hg19: Annotation for Illumina’s EPIC methylation arrays. 2016.
Abecasis GR, Auton A, Brooks LD, DePristo MA, Durbin RM, Handsaker RE, et al. An integrated map of genetic variation from 1,092 human genomes. Nature. 2012;491(7422):56–65.
Anderson CA, Pettersson FH, Clarke GM, Cardon LR, Morris AP, Zondervan KT. Data quality control in genetic case-control association studies. Nat Protoc. 2010;5(9):1564–73.
Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, et al. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;43(7): e47.
Phipson B, Maksimovic J, Oshlack A. missMethyl: an R package for analyzing data from Illumina’s HumanMethylation450 platform. Bioinformatics. 2016;32(2):286–8.
Breeze CE, Reynolds AP, van Dongen J, Dunham I, Lazar J, Neph S, et al. eFORGE v.20: updated analysis of cell type-specific signal in epigenomic data. Bioinformatics. 2019;35(22):4767–9.
Mordaunt CE, Jianu JM, Laufer BI, Zhu Y, Hwang H, Dunaway KW, et al. Cord blood DNA methylome in newborns later diagnosed with autism spectrum disorder reflects early dysregulation of neurodevelopmental and X-linked genes. Genome Med. 2020;12(1):88.
Zhu Y, Mordaunt CE, Yasui DH, Marathe R, Coulson RL, Dunaway KW, et al. Placental DNA methylation levels at CYP2E1 and IRS2 are associated with child outcome in a prospective autism study. Hum Mol Genet. 2019;28(16):2659–74.
Zhu Y, Gomez JA, Laufer BI, Mordaunt CE, Mouat JS, Soto DC, et al. Placental methylome reveals a 22q13.33 brain regulatory gene locus associated with autism. Genome Biol. 2022;23(1):46.
Laufer BI, Hwang H, Jianu JM, Mordaunt CE, Korf IF, Hertz-Picciotto I, et al. Low-pass whole genome bisulfite sequencing of neonatal dried blood spots identifies a role for RUNX1 in Down syndrome DNA methylation profiles. Hum Mol Genet. 2021;29(21):3465–76.
Krueger F, Andrews SR. Bismark: a flexible aligner and methylation caller for Bisulfite-Seq applications. Bioinformatics. 2011;27(11):1571–2.
Coulson RL, Yasui DH, Dunaway KW, Laufer BI, Vogel Ciernia A, Zhu Y, et al. Snord116-dependent diurnal rhythm of DNA methylation in mouse cortex. Nat Commun. 2018;9(1):1616.
Krueger F. 2015. Trim Galore!: A wrapper tool around Cutadapt and FastQC to consistently apply quality and adapter trimming to FastQ files. Babraham Inst.
Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012;9(4):357–9.
Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al. The Sequence Alignment/Map format and SAMtools. Bioinformatics. 2009;25(16):2078–9.
Ewels P, Magnusson M, Lundin S, Käller M. MultiQC: summarize analysis results for multiple tools and samples in a single report. Bioinformatics. 2016;32(19):3047–8.
Korthauer K, Chakraborty S, Benjamini Y, Irizarry RA. Detection and accurate false discovery rate control of differentially methylated regions from whole genome bisulfite sequencing. Biostatistics. 2019;20(3):367–83.
Hansen KD, Langmead B, Irizarry RA. BSmooth: from whole genome bisulfite sequencing reads to differentially methylated regions. Genome Biol. 2012;13(10):R83.
Nutrient and DNA methylation measures and analyses were supported by the National Institutes of Health (R01 ES025574, PI: Schmidt). Cord blood DNA methylation was also funded by R01 ES025531, PI: Fallin. Funding for the EARLI study was provided by the National Institutes of Health (R01ES016443, PI: Newschaffer) and Autism Speaks (003953 PI: Newschaffer). The MARBLES study was funded by National Institutes of Health grants (R24ES028533, R01ES020392, R01ES028089, R01ES020392, P01ES011269) and an EPA STAR grant (#RD-83329201). Support for this research was also provided by the National Institutes of Health grants (P30 ES017885 and R24 ES030893, PI: Fallin). Mr. Dou and Dr. Bakulski were also supported by grants (R01 ES025531, PI: Fallin; R01 AG067592, MPI: Bakulski; R01 MD013299). Dr. Zhu and Dr. LaSalle were supported by R01 ES029213. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.
Ethics approval and consent to participate
The EARLI study was reviewed and approved by Human Subjects Institutional Review Boards (IRBs) from each of the four study sites (Johns Hopkins University, Drexel University, University of California Davis, and Kaiser Permanente Northern California). The MARBLES protocol was reviewed and approved by the Human Subjects IRB from University of California Davis. The University of Michigan Institutional Review Board (HUM00116291) approved these secondary data analyses.
Consent for publication
All co-authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Prenatal vitamin intake in first month of pregnancy and DNA methylation in cord blood and placenta. This supplementary material contains 7 supplementary figures and 8 supplementary table. An additional 10 large supplementary tables are provided as comma separated value tables.
Single site DNA methylation array adjusted regression results: Use of prenatal vitamins in the first month of pregnancy in placenta of EARLI cohort.
Single site DNA methylation array adjusted regression results: Use of prenatal vitamins in the first month of pregnancy in placenta of MARBLES cohort.
Single site DNA methylation array adjusted regression results: Use of prenatal vitamins in the first month of pregnancy in cord blood of EARLI cohort.
Single site DNA methylation array adjusted regression results: Use of prenatal vitamins in the first month of pregnancy in cord blood of MARBES cohort.
Gene ontology ranking by p-value for EARLI placenta, EARLI cord blood, MARBLES placenta, and MARBLES cord blood, ordered by rank sum across the four cohort/tissue groups.
Gene ontology enrichment results for MARBLES cord blood.
Gene ontology enrichment results for EARLI cord blood.
Gene ontology enrichment results for MARBLES placenta.
Gene ontology enrichment results for EARLI placenta.
Differentially methylated region analysis for whole genome bisulfite sequencing data.
About this article
Cite this article
Dou, J.F., Middleton, L.Y.M., Zhu, Y. et al. Prenatal vitamin intake in first month of pregnancy and DNA methylation in cord blood and placenta in two prospective cohorts. Epigenetics & Chromatin 15, 28 (2022). https://doi.org/10.1186/s13072-022-00460-9
- DNA methylation
- Prenatal vitamins
- Pregnancy cohort
- Cord blood