No impact of prenatal paracetamol and folic acid exposure on cord blood DNA methylation in children with attention-deficit/hyperactivity disorder

Pharmacoepigenetic studies are important to understand the mechanisms through which medications influence the developing fetus. For instance, we and others have reported associations between prenatal paracetamol exposure and offspring DNA methylation (DNAm). Additionally, folic acid (FA) intake during pregnancy has been associated with DNAm in genes linked to developmental abnormalities. In this study, we aimed to: (i) expand on our previous findings showing differential DNAm associated with long-term prenatal paracetamol exposure in offspring with attention-deficit/hyperactivity disorder (ADHD), and (ii) examine if there is an interaction effect of FA and paracetamol on DNAm in children with ADHD. We used data from the Norwegian Mother, Father and Child Cohort Study (MoBa) and the Medical Birth Registry of Norway (MBRN). We did not identify any impact of paracetamol or any interaction effect of paracetamol and FA on cord blood DNAm in children with ADHD. Our results contribute to the growing literature on prenatal pharmacoepigenetics, but should be replicated in other cohorts. Replication of pharmacoepigenetic studies is essential to ensure robust findings and to increase the clinical relevance of such studies.

1 Introduction

There is an increasing interest in understanding how maternal medication use during pregnancy may affect epigenetic patterns and impact fetal development (Olstad et al., 2021). Epigenetics entails modifications to the DNA which may alter gene expression, without changing the DNA sequence (Rakyan et al., 2011). DNA methylation (DNAm) is the most commonly studied epigenetic modification, whereby a methyl group is attached to cytosine-phosphate-guanine dinucleotides (CpGs) (Rakyan et al., 2011). DNAm is reversible and influenced by both genetics and environmental factors, such as medications (Rakyan et al., 2011; Olstad et al., 2021). The timely establishment of epigenetic patterns is pivotal for normal embryonic and fetal development (Slieker et al., 2015; Luo et al., 2018). During this time, the epigenome is highly plastic, to accommodate the complex process of cellular differentiation, essential to develop functioning organs such as the brain (Slieker et al., 2015; Ciernia and LaSalle, 2016; Luo et al., 2018). Consequently, in fetal life, the epigenome is particularly susceptible to environmental influences, which may contribute to deviations in the normal epigenetic orchestration of cellular differentiation. To this end, pharmacoepigenetic studies are useful to better understand the mechanisms through which medications may impact the developing fetus.

Epidemiological studies have suggested adverse neurodevelopmental outcomes of long-term paracetamol use during pregnancy, including increased risks of attention-deficit/hyperactivity disorder (ADHD) and autism spectrum disorder (ASD) (Bauer et al., 2021). In contrast, folic acid (FA) intake during pregnancy has been associated with reduced ADHD and ASD symptoms in the child (Li et al., 2019). Furthermore, maternal FA intake in pregnancy has been found to reduce the negative impact of prenatal antiepileptic medication exposure on verbal abilities in children at both 1.5, 3, 5, and 8 years old (Husebye et al., 2018; Husebye et al., 2020). These findings suggest FA as a moderator of the association between prenatal medication exposure and child verbal abilities. FA is an essential methyl donor for DNAm. Interestingly, a meta-analysis (Joubert et al., 2016) and a recent FA intervention study (Ondičová et al., 2022), found an association of maternal FA intake during pregnancy with differences in DNAm at genes linked to developmental abnormalities. Several studies have also found associations between prenatal paracetamol exposure and differential DNAm in cord blood (Gervin et al., 2017; Eslamimehr et al., 2022) and placenta (Addo et al., 2019), but none of the differences in DNAm overlapped across the studies (Gervin et al., 2017; Addo et al., 2019; Olstad et al., 2021; Eslamimehr et al., 2022).

Based on these pharmacoepidemiological and -epigenetic findings, we postulated that DNAm could be a plausible biological mechanism of action for the paracetamol-ADHD associations (Gervin et al., 2017). Therefore, we have previously conducted an epigenome-wide association study (EWAS) on the association of prenatal paracetamol exposure with DNAm in children with ADHD (Gervin et al., 2017). In our previous study, we found that DNAm at more than 6,000 CpGs was associated with prenatal paracetamol exposure only in children later developing ADHD (Gervin et al., 2017). To strengthen these findings, we wanted to replicate and expand on our previous study, both by increasing the number of samples and CpGs, and by exploring whether FA may impact a potential effect of paracetamol on DNAm. To this end, we aimed to: (i) expand on our previous findings showing epigenome-wide differences in DNAm associated with long-term prenatal paracetamol exposure (≥20 days) in children with ADHD (Gervin et al., 2017), and (ii) examine if there is an interaction effect of FA and paracetamol on DNAm in children with ADHD. To do this, we selected umbilical cord blood samples from the Norwegian Mother, Father and Child Cohort Study (MoBa), which contains information on maternal use of paracetamol and FA during pregnancy, and conducted an EWAS.

2 Materials and methods

This study is based on umbilical cord blood samples from MoBa, conducted by the Norwegian Institute of Public Health (NIPH) (Magnus et al., 2006; Magnus et al., 2016). MoBa is a prospective, population-based birth cohort (n = 114,500 children, n = 95,200 mothers and n = 75,200 fathers), including births in Norway between 1998 and 2008 (Magnus et al., 2006; Magnus et al., 2016). In 40.6% of pregnancies, parents consented to participate. Throughout pregnancy and childhood, the participants complete multiple questionnaires. MoBa also includes a biobank with approximately 90,000 blood samples collected from both parents during pregnancy, and from mother and child (umbilical cord) at birth (Paltiel et al., 2014). This study is based on Data Version 8 released by MoBa in 2015. All MoBa participants have provided a written, informed consent to participate in the cohort and can retract their consent at any time. The establishment of MoBa and initial data collection was based on a license from the Norwegian Data Protection Agency and approval from The Regional Committees for Medical and Health Research Ethics. The MoBa cohort is currently regulated by the Norwegian Health Registry Act. All the analyzed data were de-identified, and the linking of MoBa to the relevant health registries was handled by NIPH and the respective registries. Our study was approved by the Regional Committee for Medical Research Ethics South East Norway.

The selection of individuals for the study was based on observational data from MoBa questionnaires Q1 (gestational weeks 0–13), Q3 (gestational weeks 13–29) and Q4 (gestational week 30 to delivery), the Medical Birth Registry of Norway (MBRN) and the Norwegian Patient Registry (NPR). MBRN is a national health registry containing information about all births in Norway and NPR contains diagnoses asserted by specialists in governmental hospitals and outpatient clinics. MoBa was linked to NPR and MBRN using the personal 11-digit identification number unique to every permanent resident of Norway. We included 261 pregnancies from MoBa, divided into three groups: (i) children with ADHD prenatally exposed to paracetamol ≥20 days (exposed group; n = 61), (ii) children with ADHD and unexposed to paracetamol (ADHD-controls; n = 100), and (iii) children without ADHD and unexposed to paracetamol (population controls; n = 100). We did not include a paracetamol-control group (i.e., children prenatally exposed to paracetamol ≥20 days without ADHD), as in our prior study we did not identify any significant differences in DNAm between this group and the exposed group (Gervin et al., 2017). Pre-term births (<37 weeks) and twins were excluded. The definitions of the different measures and the selection of covariates are presented in the Supplementary Methods (Supplementary Methods).

DNAm was assessed using the Infinium HumanMethylation EPIC BeadChip (Illumina). Samples were randomly allocated to sample plates and beadchips, and processed as previously described (Gervin et al., 2017). All analyses were performed in the R programming language (http://www.r-project.org/). The DNAm data was quality controlled using the quality control module implemented in RnBeads (v. 2.8.1) (Assenov et al., 2014; Müller et al., 2019). First, samples with >5% low-quality CpGs or low bisulfite intensity were removed (0 samples). CpGs with >5% low-quality values were also removed (n = 8,947). Low-quality probes either exhibited a detection p-value > 10⁻⁶ or a bead count <3. We performed background correction with the ENmix exponential-truncated-normal out-of-band (oob) method (Xu et al., 2016), dye bias correction with RELIC (REgression on Logarithm of Internal Control probes) (Xu et al., 2017) and probe-type correction with RCP (Regression of Correlated Probes) (Niu et al., 2016). We then removed probes with SNPs overlapping with the CpG interrogation site or the nucleotide extension site (n = 29,176), cross-reactive probes (n = 14,921) (McCartney et al., 2016; Pidsley et al., 2016) and probes on the sex chromosomes (n = 17,532). The final data set included 795,515 probes and 261 samples.

We used β values (ratio of methylated signal to total signal) for visualizations and M values for statistical tests (Du et al., 2010). Principal component analysis (PCA) on the DNAm data was used to test the strength of association between covariates and DNAm variation (tests applied on categorical and continuous variables included one-way analysis of variance [ANOVA] and Spearman’s correlation test, respectively). To identify differentially methylated CpGs associated with paracetamol, we fit linear regression models onto the mean DNAm differences in limma (Smyth et al., 2005). We included the CD8⁺ T cell proportion as a covariate in all models. Interaction was assessed by including an interaction term of paracetamol and FA in the model. We pairwise compared the exposed group to the ADHD-control and the population control groups. All comparisons were adjusted for multiple testing (false discovery rate [FDR] <0.05) (Benjamini, 2010). We performed a surrogate variable analysis to examine unmeasured sources of variation in DNAm (Leek et al., 2012). Further details on the statistical analyses are presented in the Supplementary methods (Supplementary Methods).

3 Results

To enable a systematic replication and expansion of our previous findings (Gervin et al., 2017), we selected three groups: one exposed group, consisting of children with ADHD prenatally exposed to paracetamol ≥20 days (n = 61), and two control groups, the ADHD-control group, including children with ADHD and unexposed to paracetamol (n = 100), and the population control group, including children without ADHD and unexposed to paracetamol (n = 100). The samples were selected from MoBa, the same cohort as our previous study was based on, and 17 of the 19 exposed samples in the original study were also included in the present study. We selected model covariates by assessing their contribution to DNAm variation using PCA (Supplementary Methods). These analyses identified that the estimated CD8⁺ T cell proportion was significantly associated with principal components (PCs) 2 and 3, and differed between the comparison groups (Supplementary Figure S1; Supplementary Tables S1, S2), and therefore, this covariate was included in the regression models.

We ran three different models to assess the impact of paracetamol and FA on DNAm in children with ADHD (Table 1): (i) a crude model to examine whether prenatal paracetamol exposure was associated with DNAm in the cord blood of children with ADHD, (ii) an adjusted model where we adjusted the crude model for FA intake during pregnancy to investigate the influence of FA status on the variance in DNAm, and (iii) an interaction model where we examined any interaction effect of FA and paracetamol on DNAm. We ran the three models comparing the exposed group to each of the two control groups separately, and also ran sensitivity analyses excluding the CD8⁺ T cell proportion from the models, which did not change the results (Supplementary Figure S2).

TABLE 1

TABLE 1. Overview of the three models used to assess the association of prenatal paracetamol and folic acid (FA) exposure with cord blood DNA methylation (DNAm) in children with attention-deficit/hyperactivity disorder (ADHD). We ran each model separately for the exposed group compared to each of the two control groups.

To test the association of DNAm with prenatal paracetamol exposure and ADHD, we ran linear regression models comparing the exposed group to the ADHD-control group and the population control group (the crude model). This analysis did not identify any significant differences between the exposed group and either of the two control groups (Figure 1A; Supplementary Figure S3A; Supplementary Table S3), and thereby did not replicate our previous findings. We then assessed whether FA could influence the association between paracetamol exposure and DNAm by running an adjusted model, adjusting for FA intake during pregnancy. This analysis did not reveal any significant CpGs and the p values were largely similar to those observed in the crude model (Figure 1B; Supplementary Figure S3B). Finally, to understand whether there is an interaction effect of paracetamol and FA intake on DNAm in children with ADHD, we ran an interaction model where the group and FA exposure variables interacted. This analysis neither revealed any significant CpGs (Figure 1C; Supplementary Figure S3C).

FIGURE 1

FIGURE 1. Volcano plots comparing mean differences in DNAm (Δβ) for each CpG between the exposed group and either the ADHD-control group or the population control group. (A) Crude model, only adjusted for CD8⁺ T cell proportion. (B) Adjusted model, adjusting for FA intake during pregnancy and CD8⁺ T cell proportion. (C) Interaction model, assessing any interaction effect of prenatal paracetamol and FA exposure on DNAm. Each CpG is plotted by the -log₁₀ of the p values against the mean per-CpG difference in DNAm (Δβ) of the compared groups. In (C) mean Δβ reflects the interaction term (i.e., mean Δβ_{no FA}–Δβ_FA). The red dotted line indicates the FDR-adjusted p value significance threshold (FDR < 0.05).

For the crude model, we also explored whether the difference in results between our present and previous study could be explained by differences in normalization methods, selection of covariates or number of probes on the microarrays (EPIC vs. 450 k, respectively). Conducting our analysis using the same normalization methods (i.e., noob and BMIQ) and/or the same covariates as in our previous study (i.e., sex, smoking, maternal age, gestational age, estimated cell type composition, BeadChip and bisulfite conversion plate), we still find no significant CpGs (Supplementary Table S4). Yet, some CpGs exhibit FDR-adjusted p values closer to significance (p ≥ 0.13), but none of these overlap with the significant findings in our previous study. Finally, to assess whether the heightened burden of multiple testing associated with the increased number of probes on the EPIC microarray could explain our null findings, we also ran the crude model only on the probes overlapping between the 450 k and EPIC microarrays (Supplementary Table S5). This analysis did not reveal any CpGs significantly associated with prenatal paracetamol exposure in children with ADHD.

4 Discussion

Overall, we identified no impact of paracetamol nor interaction of paracetamol and FA on cord blood DNAm in children with ADHD in our data. Results from the crude model did not identify any significant differences in DNAm between the exposed group and either of the two control groups. Consequently, we did not replicate our previous results (Gervin et al., 2017). These findings were surprising, as the present study is based on a larger number of samples from the same cohort (MoBa), using the same study design and similar inclusion criteria as in the original study. However, our studies are performed 5 years apart and methods have evolved, including the introduction of the Illumina Infinium EPIC platform and novel analysis methods, such as normalization and cell type deconvolution procedures. While we could not explain our replication failure solely by differing DNAm measurements between platforms (Olstad et al., 2022), we cannot exclude that other aspects of the analysis may have contributed to the lack of replication.

To explore if differences in the analysis between our present and previous study contributed to our non-replication, we examined whether varying normalization methods, covariates and number of probes on the microarrays would change the results. We did not identify any CpGs significantly associated with prenatal paracetamol exposure in children with ADHD in these analyses.

Our findings also differ from two recent studies on prenatal paracetamol exposure and DNAm, which found differential DNAm at some CpGs associated with paracetamol (Addo et al., 2019; Eslamimehr et al., 2022). However, comparison to these studies is difficult for several reasons, including differing study designs, exposure definitions, phenotypes, tissue types and analysis pipelines. There are also other challenges when comparing current prenatal pharmacoepigenetic studies (Olstad et al., 2021), which may contribute to the lack of replication and overlap of findings.

When assessing whether FA contributed to the variance in DNAm or if there was an interaction effect of paracetamol and FA on DNAm, we did not find any significant CpGs, suggesting that there is no effect of paracetamol and FA intake on DNAm in children with ADHD. Previous studies have reported differential DNAm associated with maternal FA intake during pregnancy in genes linked to developmental abnormalities (Joubert et al., 2016; Ondičová et al., 2022). However, these studies examined FA intake as the main exposure rather than as an interacting factor with medication use (Joubert et al., 2016; Ondičová et al., 2022), and therefore, are not directly comparable to the current study.

While we have improved multiple aspects of our previous study (Gervin et al., 2017), there are some limitations to the current study. The sample sizes of the groups are relatively small, albeit larger than in our previous study. Although the MoBa cohort is one of the world’s largest prospective birth cohorts, only 61 children with a clinical ADHD diagnosis were exposed to paracetamol ≥20 days. Additionally, we did not assess the dose and timing of FA intake during pregnancy. However, in Norway, the recommended FA supplement dosage is 400 μg/day, starting 1 month prior to pregnancy and lasting throughout the first trimester of pregnancy.

In conclusion, this study did not replicate previous findings in MoBa or other studies investigating the influence of paracetamol on DNAm, and did not identify any interaction effect of paracetamol and FA on DNAm in children with ADHD. Our results contribute to the growing literature on prenatal pharmacoepigenetics, but should be replicated in other cohorts. Replication of pharmacoepigenetic studies is essential to ensure robust findings and to increase the clinical relevance of such studies.

Data availability statement

The datasets presented in this article are available from the Norwegian Mother, Father and Child Cohort Study, but restrictions apply to the availability of these data and so are not publicly available. However, data are available from the authors upon reasonable request and with permission from the Norwegian Mother, Father and Child Cohort Study. Requests to access the datasets should be directed to EWO, e.w.olstad@farmasi.uio.no.

Ethics statement

The studies involving human participants were reviewed and approved by the Regional Committees for Medical and Health Research Ethics (REC). Written informed consent to participate in this study was provided by the participants’ legal guardian/next of kin.

Author contributions

EWO, HN, RL, and KG conceived the idea of and designed the study. EWO conducted the analyses. EWO generated plots and tables and drafted the first version of the paper. EWO, HN, RL, and KG all revised the paper. All authors contributed to the article and approved the submitted version.

Funding

This work was supported by the European Research Council (ERC) Starting Grant “DrugsInPregnancy” (grant number: 678033) and the Research Council of Norway (NFR) “Medications during pregnancy: effects on neurodevelopmental disorders in childhood and epigenetic outcomes” (grant number: 241117). MoBa is supported by the Norwegian Ministry of Health and the Ministry of Education and Research, NIH/NIEHS (contract number: NO1-ES-75558), NIH/NINDS (grant numbers: UO1 NS 047537-01 and UO1 NS047537-06A1), and the Norwegian Research Council/FUGE (grant numbers; 151918/S10). The Norwegian Mother, Father and Child Cohort Study is supported by the Norwegian Ministry of Health and Care Services and the Ministry of Education and Research.

Acknowledgments

We are grateful to all the participating families in Norway taking part in this on-going cohort study.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2023.1204879/full#supplementary-material

References

Addo, K. A., Bulka, C., Dhingra, R., Santos, H. P., Smeester, L., O’Shea, T. M., et al. (2019). Acetaminophen use during pregnancy and DNA methylation in the placenta of the extremely low gestational age newborn (ELGAN) cohort. Environ. Epigenetics 5, dvz010. doi:10.1093/eep/dvz010

PubMed Abstract | CrossRef Full Text | Google Scholar

Assenov, Y., Müller, F., Lutsik, P., Walter, J., Lengauer, T., and Bock, C. (2014). Comprehensive analysis of DNA methylation data with RnBeads. Nat. Methods 11, 1138–1140. doi:10.1038/nmeth.3115

PubMed Abstract | CrossRef Full Text | Google Scholar

Bauer, A. Z., Swan, S. H., Kriebel, D., Liew, Z., Taylor, H. S., Bornehag, C. G., et al. (2021). Paracetamol use during pregnancy — A call for precautionary action. Nat. Rev. Endocrinol. 17, 757–766. doi:10.1038/s41574-021-00553-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Benjamini, Y. (2010). Discovering the false discovery rate. J. R. Stat. Soc. Ser. B 72, 405–416. doi:10.1111/j.1467-9868.2010.00746.x

CrossRef Full Text | Google Scholar

Ciernia, A. V., and LaSalle, J. (2016). The landscape of DNA methylation amid a perfect storm of autism aetiologies. Nat. Rev. Neurosci. 17, 411–423. doi:10.1038/nrn.2016.41

PubMed Abstract | CrossRef Full Text | Google Scholar

Du, P., Zhang, X., Huang, C-C., Jafari, N., Kibbe, W. A., Hou, L., et al. (2010). Comparison of Beta-value and M-value methods for quantifying methylation levels by microarray analysis. BMC Bioinforma. 11, 587. doi:10.1186/1471-2105-11-587

PubMed Abstract | CrossRef Full Text | Google Scholar

Eslamimehr, S., Jones, A. D., Anthony, T. M., Arshad, S. H., Holloway, J. W., Ewart, S., et al. (2022). Association of prenatal acetaminophen use and acetaminophen metabolites with DNA methylation of newborns: Analysis of two consecutive generations of the isle of wight birth cohort. Environ. Epigenetics 8, dvac002–10. doi:10.1093/eep/dvac002

CrossRef Full Text | Google Scholar

Gervin, K., Nordeng, H., Ystrom, E., Reichborn-kjennerud, T., and Lyle, R. (2017). Long-term prenatal exposure to paracetamol is associated with DNA methylation differences in children diagnosed with ADHD. Clin. Epigenetics 9, 77. doi:10.1186/s13148-017-0376-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Husebye, E. S. N., Gilhus, N. E., Riedel, B., Spigset, O., Daltveit, A. K., and Bjørk, M. H. (2018). Verbal abilities in children of mothers with epilepsy: Association to maternal folate status. Neurology 91, e811–e821. doi:10.1212/WNL.0000000000006073

PubMed Abstract | CrossRef Full Text | Google Scholar

Husebye, E. S. N., Gilhus, N. E., Spigset, O., Daltveit, A. K., and Bjørk, M. H. (2020). Language impairment in children aged 5 and 8 years after antiepileptic drug exposure in utero – The Norwegian mother and child cohort study. Eur. J. Neurol. 27, 667–675. doi:10.1111/ene.14140

PubMed Abstract | CrossRef Full Text | Google Scholar

Joubert, B. R., Den Dekker, H. T., Felix, J. F., Bohlin, J., Ligthart, S., Beckett, E., et al. (2016). Maternal plasma folate impacts differential DNA methylation in an epigenome-wide meta-analysis of newborns. Nat. Commun. 7, 10577. doi:10.1038/ncomms10577

PubMed Abstract | CrossRef Full Text | Google Scholar

Leek, J. T., Johnson, W. E., Parker, H. S., Jaffe, A. E., and Storey, J. D. (2012). The SVA package for removing batch effects and other unwanted variation in high-throughput experiments. Bioinformatics 28, 882–883. doi:10.1093/bioinformatics/bts034

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, M., Francis, E., Hinkle, S. N., Ajjarapu, A. S., and Zhang, C. (2019). Preconception and prenatal nutrition and neurodevelopmental disorders: A systematic review and meta-analysis. Nutr 11, 1628. doi:10.3390/nu11071628

CrossRef Full Text | Google Scholar

Luo, C., Hajkova, P., and Ecker, J. R. (2018). Dynamic DNA methylation: In the right place at the right time. Science 361, 1336–1340. doi:10.1126/science.aat6806

PubMed Abstract | CrossRef Full Text | Google Scholar

Magnus, P., Birke, C., Vejrup, K., Haugan, A., Alsaker, E., Daltveit, A. K., et al. (2016). Cohort profile update: The Norwegian mother and child cohort study (MoBa). Int. J. Epidemiol. 45, 382–388. doi:10.1093/ije/dyw029

PubMed Abstract | CrossRef Full Text | Google Scholar

Magnus, P., Irgens, L. M., Haug, K., Nystad, W., Skjærven, R., Stoltenberg, C., et al. (2006). Cohort profile: The Norwegian mother and child cohort study (MoBa). Int. J. Epidemiol. 35, 1146–1150. doi:10.1093/ije/dyl170

PubMed Abstract | CrossRef Full Text | Google Scholar

McCartney, D. L., Walker, R. M., Morris, S. W., McIntosh, A. M., Porteous, D. J., and Evans, K. L. (2016). Identification of polymorphic and off-target probe binding sites on the Illumina Infinium MethylationEPIC BeadChip. Genomics Data 9, 22–24. doi:10.1016/j.gdata.2016.05.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Müller, F., Scherer, M., Assenov, Y., Lutsik, P., Walter, J., Lengauer, T., et al. (2019). RnBeads 2.0: Comprehensive analysis of DNA methylation data. Genome Biol. 20, 55. doi:10.1186/s13059-019-1664-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Niu, L., Xu, Z., and Taylor, J. A. (2016). Rcp: A novel probe design bias correction method for Illumina methylation BeadChip. Bioinformatics 32, 2659–2663. doi:10.1093/bioinformatics/btw285

PubMed Abstract | CrossRef Full Text | Google Scholar

Olstad, E. W., Nordeng, H. M. E., and Gervin, K. (2021). Prenatal medication exposure and epigenetic outcomes: A systematic literature review and recommendations for prenatal pharmacoepigenetic studies. Epigenetics 17, 357–380. doi:10.1080/15592294.2021.1903376

PubMed Abstract | CrossRef Full Text | Google Scholar

Olstad, E. W., Nordeng, H. M. E., Sandve, G. K., Lyle, R., and Gervin, K. (2022). Low reliability of DNA methylation across Illumina Infinium platforms in cord blood: Implications for replication studies and meta-analyses of prenatal exposures. Clin. Epigenetics 14, 80. doi:10.1186/s13148-022-01299-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Ondičová, M., Irwin, R. E., Thursby, S. J., Hilman, L., Caffrey, A., Cassidy, T., et al. (2022). Folic acid intervention during pregnancy alters DNA methylation, affecting neural target genes through two distinct mechanisms. Clin. Epigenetics 14, 63–20. doi:10.1186/s13148-022-01282-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Paltiel, L., Haugan, A., Skjerden, T., Harbak, K., Bækken, S., Stensrud, N. K., et al. (2014). The biobank of the Norwegian mother and child cohort study – present status. Nor. Epidemiol. 24, 29–35.

CrossRef Full Text | Google Scholar

Pidsley, R., Zotenko, E., Peters, T. J., Lawrence, M. G., Risbridger, G. P., Molloy, P., et al. (2016). Critical evaluation of the Illumina MethylationEPIC BeadChip microarray for whole-genome DNA methylation profiling. Genome Biol. 17, 208. doi:10.1186/s13059-016-1066-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Rakyan, V. K., Down, T. A., Balding, D. J., and Beck, S. (2011). Epigenome-wide association studies for common human diseases. Nat. Rev. Genet. 12, 529–541. doi:10.1038/nrg3000

PubMed Abstract | CrossRef Full Text | Google Scholar

Slieker, R. C., Roost, M. S., van Iperen, L., Suchiman, H. E. D., Tobi, E. W., Carlotti, F., et al. (2015). DNA methylation landscapes of human fetal development. PLOS Genet. 11, e1005583. doi:10.1371/journal.pgen.1005583

PubMed Abstract | CrossRef Full Text | Google Scholar

Smyth, G. K. (2005). “limma: Linear models for microarray data,” in Bioinforma comput biol solut using R bioconductor. Editors R. Gentleman, V. J. Carey, W. Huber, R. A. Irizarry, and S. Dudoit (New York, NY: Springer), 397–420.

Google Scholar

Xu, Z., Langie, S. A. S., De Boever, P., Taylor, J. A., and Niu, L. (2017). Relic: A novel dye-bias correction method for Illumina methylation BeadChip. BMC Genomics 18, 4. doi:10.1186/s12864-016-3426-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, Z., Niu, L., Li, L., and Taylor, J. A. (2016). ENmix: A novel background correction method for Illumina HumanMethylation450 BeadChip. Nucleic Acids Res. 44, e20. doi:10.1093/nar/gkv907

PubMed Abstract | CrossRef Full Text | Google Scholar