Skip to main content

Long-term prenatal exposure to paracetamol is associated with DNA methylation differences in children diagnosed with ADHD

Abstract

Background

Epidemiological studies have shown that long-term exposure to paracetamol during pregnancy is associated with attention-deficit/hyperactivity disorder (ADHD). The mechanism by which paracetamol may modulate the increased risk of developing ADHD is currently unknown. We have conducted an epigenome-wide association study (n = 384 cord blood samples) and investigated whether prenatal exposure to paracetamol is associated with DNA methylation in children diagnosed with ADHD.

Results

Analyses identified significant differences in DNA methylation (n = 6211 CpGs) associated with prenatal exposure to paracetamol for more than 20 days in children diagnosed with ADHD compared to controls. In addition, these samples were differentially methylated compared to samples with ADHD exposed to paracetamol for less than 20 days (n = 2089 CpGs) and not exposed to paracetamol (n = 193 CpGs). Interestingly, several of the top genes ranked according to significance and effect size have been linked to ADHD, neural development, and neurotransmission. Gene ontology analysis revealed enrichment of pathways involved in oxidative stress, neurological processes, and the olfactory sensory system, which have previously been implicated in the etiology of ADHD.

Conclusions

These initial findings suggest that in individuals susceptible to ADHD, prenatal long-term exposure to paracetamol is associated with DNA methylation differences compared to controls.

Background

Attention-deficit/hyperactivity disorder (ADHD) is a common childhood disorder defined by inattention, hyperactivity, and impulsivity, which impair normal functioning and development. It is one of the most prevalent psychiatric disorders of childhood affecting approximately 5% of children worldwide, often continuing into adulthood [1]. ADHD is considered to have a strong genetic basis, with heritability estimates up to 70% [2]. Although candidate gene studies have identified genes associated with ADHD [2], numerous genome-wide association studies (GWAS) [3,4,5,6,7,8,9] and meta-analyses [10] have not found any significant ADHD risk variants [11].

Several epidemiological studies have shown that exposure to a range of putative environmental risk factors are associated with ADHD, particularly during the prenatal and perinatal period. However, these factors can only be regarded as correlates because, at present, convincing evidence that the associations reflect causal relationships have not been established [12]. ADHD is a complex disorder influenced by the interplay of various genetic and environmental factors. The mechanisms by which the environment mediates ADHD susceptibility is not well understood, but epigenetic changes may be involved [13,14,15,16,17].

Recently, five epidemiological studies have revealed a consistent association between long-term exposure to paracetamol (acetaminophen) during pregnancy and ADHD symptoms [18,19,20,21] and adverse neurodevelopment in children [22]. Paracetamol, which is available without prescription, is the most commonly used drug during pregnancy taken by ~50% of pregnant woman to treat fever, headache, and other pain conditions [23]. If the association between prenatal paracetamol exposure and neurodevelopmental disorders is causal, the underlying mechanisms remain to be clarified.

In this study, we aimed to assess whether epigenetic differences are associated with prenatal exposure to paracetamol and development of ADHD. To do this, we have conducted an epigenome-wide association study (EWAS) in samples selected from the Norwegian Mother and Child Cohort (MoBa) and investigated whether long-term prenatal exposure to paracetamol was associated with DNA methylation changes in cord blood and clinical ADHD diagnoses recorded in the Norwegian Patient Registry (NPR) from 2008 to 2014.

Results

Neither paracetamol nor ADHD alone is associated with DNA methylation differences

The samples in this study were selected to enable a systematic investigation of the association between DNA methylation and paracetamol exposure in children with and without ADHD (Table 1). First, we stratified the analyses on paracetamol exposure and ADHD separately to evaluate the associations with DNA methylation. Pairwise analyses of the groups exposed to paracetamol (drug exposure group) or diagnosed with ADHD (ADHD control group) alone compared to control samples revealed no significant differences in DNA methylation (FDR <0.05, data not shown).

Table 1 Sampling and study design

Long-term paracetamol exposure is associated with DNA methylation differences in children with ADHD

Next, we investigate the synergistic effect of prenatal exposure to paracetamol on DNA methylation differences in children with ADHD. To do this, we stratified the analyses on samples exposed to paracetamol (synergistic effect group, Table 1) and controls. Further, we performed comparisons of the synergistic group to the ADHD control group and the drug exposure group to enable interpretation of the associations.

Comparison of the samples exposed to paracetamol with ADHD (synergistic effect group) to each of the other three groups did not identify any differential DNA methylation below the significance threshold (FDR <0.05, data not shown). However, analyses contrasting sporadic and long-term paracetamol use revealed that exposure to paracetamol for ≥20 days was associated with DNA methylation differences in samples with ADHD. These samples showed differential DNA methylation compared to the control group after FDR correction in a large number of CpG sites (n = 6211 CpGs, mean = 0.012, range = −0.28 to 0.20; Fig. 1a, Additional file 1: Table S1). The 6211 CpGs were 90% Winsorized in the long-term exposed samples with ADHD to investigate the potential influence of outliers. The results from this analysis showed that a small proportion of the CpGs (n = 697, 11.2%) were driven by outliers.

Fig. 1
figure 1

Group-wise differences in DNA methylation. Scatter plots of DNA methylation differences between a long-term exposed samples with ADHD (n = 19) and controls (n = 96), b long-term exposed samples with ADHD (n = 19) and sporadically exposed samples with ADHD (n = 77), and c long-term exposed samples with ADHD (n = 19) and unexposed samples with ADHD (n = 96). Each point represents a CpG site (n = 390,787), and significant CpGs are marked in red (BH FDR <0.05). d Venn diagram summarizing the overlap of CpGs that were differentially methylated between the three groups. e Example box plots showing the DNA methylation levels at two representative CpGs that were differentially methylated between long-term exposed samples with ADHD and unexposed samples with ADHD

The long-term exposed samples with ADHD were also differentially methylated compared to exposed samples with ADHD in both the sporadic (n = 2089 CpGs, mean = 0.025, range = −0.13 to 0.17; Fig. 1b, Additional file 2: Table S2) and unexposed samples with ADHD (n = 193 CpGs, mean = 0.029, range = −0.11 to 0.17; Fig. 1c, Additional file 3: Table S3). There was substantial overlap between the differentially methylated CpGs identified by these three comparisons (Fig. 1d). Quantile-quantile (Q-Q) plots of the observed versus expected p values and corresponding histograms of p values for the three comparisons above show strong inflation (Additional file 4: Figure S1). Here, the observed quantiles are consistently higher than their expected values under the null hypothesis of no association between long-term exposure to paracetamol and DNA methylation differences in children with ADHD to each of the three comparisons (lambda ranging from 1.7 to 2.2). We performed surrogate variable analysis (SVA) to detect possible hidden covariates with a global effect on DNA methylation between the groups. No significantly associated surrogate variables were detected after correction for the covariates already included in the statistical model (i.e., cell-type composition using a suitable reference panel based on cord blood, gestational age, maternal age, smoking, and bisulfite conversion). Hence, these results suggest that there are no unmeasured covariates with a strong confounding effect on the significantly associated DNA methylation differences.

We could not detect any effect due to use in different trimesters, but common for the samples in the high paracetamol exposure group was that they were exposed to paracetamol in at least two of the three trimesters, usually across all trimesters.

The majority of the significant CpGs between groups involved very small differences in DNA methylation. To increase the likelihood of biological relevance of the differences in DNA methylation being related to ADHD and paracetamol exposure, we chose to focus the downstream analyses on the CpGs displaying both statistical significance and a large effect size (combined score, see the “Methods” section). Interestingly, top ranked CpGs overlapped genes previously linked to ADHD (SGTB [24] and CADM1 [25]). In addition, genes involved in neural development (REST [26], ZNRF1 [27], RabGEF1 [28], and CUX1 [29]), neurotransmission (SYN2 [30] and DDAH1 [31]; Fig. 1e (left)), noradrenergic neurotransmitter system (PHOX2A), voltage channels (KCNJ8 and CACNG8 [32]), and peroximal lipid metabolism (PHYH; Fig. 1e (right)) clustered among the top ranked genes, all of which have been implicated in neurodevelopmental disorders. Regional analyses did not identify any significant differential DNA methylation between groups.

Enrichment of pathways involved in oxidative stress and the olfactory system

Gene ontology (GO) analyses of the top 100 significantly differentially methylated genes ranked on the combined score identified enrichment of pathways involved in oxidative stress among the genes showing increased DNA methylation (Additional files 1, 2, and 3: Tables S1–3). This enrichment was more prominent among the differentially methylated genes between samples with ADHD, but discordant for paracetamol exposure (Table 2). In addition, receptor-mediated endocytosis, vesicle-mediated transport, and bicarbonate transport pathways, which are involved in the neurotransmission, were enriched. Genes showing decreased DNA methylation were enriched for pathways related to the olfactory system and neurological system processes. The olfactory system and perception of smell were highly and consistently enriched among the top 100 genes identified in all the three comparisons (Table 2).

Table 2 Gene ontology results based on top 100 ranked according to significance and effect size

Discussion

We have studied DNA methylation differences associated with prenatal exposure to paracetamol and ADHD in cord blood samples selected from a large prospective birth cohort. This is the first study to investigate the recent pharmacoepidemiological findings of increased risk of ADHD after long-term prenatal exposure to paracetamol. In the pairwise comparisons of the groups, only the group with long-term exposure to paracetamol and ADHD showed differential DNA methylation. These results indicate a dosage effect of paracetamol on the observed differences in DNA methylation. This is in agreement with the epidemiological link identified in five large recent studies [18,19,20,21,22]. The same effect has been observed in adult mice where long-term exposure to paracetamol, but not ibuprofen [33], during neonatal brain development was associated with cognitive functions [34].

Several studies, based on selected candidate genes [13, 14, 16] and two recent EWAS [15, 17], have revealed an association between DNA methylation differences and ADHD symptoms. No differential DNA methylation associated with ADHD alone was identified in our study. Several factors involving variations in the sample size, criteria underlying the ADHD phenotype, study design, type of tissue, time point, and confounding factors are likely to explain the lack of replication in our study (or between any of the other studies). Most of the studies are based on either peripheral whole-blood or saliva samples from older children. Hence, the associated differences in DNA methylation are likely to also involve different environmental influences during childhood and having tissue-specific roles. The recent EWAS by Walton et al. [17] is perhaps the most comparable to our study in terms of the time point and tissue investigated. Cord blood, which is sampled at the time of birth, reflects prenatal environmental influences on the development of ADHD. However, we did not replicate any of the 13 CpGs reported in this study.

A large number of differentially methylated CpGs were identified in this study, particularly between samples with ADHD and long-term exposure to paracetamol and controls (unexposed without ADHD). The majority involved small differences in DNA methylation. Sensitivity analyses showed that outliers in this group influenced a small proportion (11.2%) of the differentially methylated CpGs. It is not possible to differentiate the influence of technical outliers from extreme DNA methylation values reflecting paracetamol dosage-effect (high number of days of exposure) associations. To increase biological relevance, the analyses and interpretation of the results were based on the genes displaying significant and the largest differences in DNA methylation. Several of the top ranked genes have previously been linked to not only ADHD but also neural development and neurotransmission. These are categories involved in the brain development and function and are also enriched among genes genetically associated with ADHD [35]. ADHD is a complex disorder, which involves many genes interacting within multiple pathways, also epigenetically dysregulated genes. Pathways are interconnected, and a defect in one may lead to dysfunction in others, potentially leading to a high number of epigenetically dysregulated genes. The top 100 genes with the largest DNA methylation differences in our study were enriched in pathways involved in oxidative stress, neurological system processes, and the olfactory system. Several studies have detected altered levels of oxidative stress in different tissues in ADHD patients, including the blood [36], saliva [37], and serum [38]. This has also been observed in autism spectrum disorders [39]. This could suggest that oxidative stress is a peripheral biomarker of ADHD or that ADHD is a systemic abnormality with consequences in the brain. The brain is particularly susceptible to oxidative stress due to a high level of oxygen utilization and a modest antioxidant defense [40]. Interestingly, oxidative damage and redox imbalance have also been suggested as one of the paracetamol toxicity mechanisms [41]. Further, pathways related to the olfactory system were highly enriched among the genes displaying decreased DNA methylation. The olfactory system is responsible for the detection of smell. Activated olfactory receptors induce a nerve impulse, which is ultimately transmitted to the brain. Deficits in the function of this system are a common feature in several neurodevelopmental disorders, including ADHD [42]. Specifically, increased odor sensitivity has been found in both children [43] and adults [44] with ADHD.

Interestingly, we identified more differentially methylated positions when comparing the long-term exposed samples to sporadic exposed samples than to unexposed samples with ADHD. Whereas this observation is difficult to explain, it could reflect the fact that dose-response relationships are usually not linear. Unfortunately, it is difficult to model and monitor a dose effect and timing during pregnancy due to the way the MoBa questionnaires were designed and completed. Future studies including a larger sample should investigate whether this trend remains.

The results presented in this study should be interpreted in light of several limitations. Although this study is based on one of the largest prospective birth cohorts in the world, the number of samples possessing both clinical ADHD diagnoses and long-term prenatal exposure to paracetamol was relatively small. Hence, the association between paracetamol and ADHD should be replicated in a larger set of samples. We report DNA methylation differences in cord blood samples, and whether these differences also reflect changes in the brain needs to be established. Although the analyses and results were focused on differentially methylated genes with large effect sizes, the effect on gene expression needs to be explored. Heritable impulsive traits are associated with paracetamol use during pregnancy [45]. Future studies should investigate the link found here while controlling for intergenerational transmission of impulsive traits.

There are also several strengths that should be taken into account. We used a clinical diagnosis of ADHD given by specialists obtained through linkage to the NPR. Although registry-based diagnoses of ADHD has been used as outcome in a number of previous Scandinavian studies on ADHD, including the abovementioned Danish birth cohort study of paracetamol use during pregnancy and ADHD [18], they have not been validated. The study is based on cord blood samples collected at birth, thereby eliminating the confounding of ADHD stimulant drugs, which are known to influence DNA methylation [46]. Moreover, as samples with maternal co-medication were excluded, there was no influence of other psychotropics or analgesics on DNA methylation. In addition, other known confounders were included as covariates. Importantly, correction for cord blood cell-type composition between groups was done based on a suitable cord blood reference data set, which we recently established [47].

Conclusion

This study identified altered DNA methylation differences at genes involved in oxidative stress, neural transmission, and olfactory sensory pathways associated with long-term exposure to paracetamol during pregnancy in children diagnosed with ADHD. These results suggest that in individuals susceptible to ADHD, prenatal long-term exposure to paracetamol is associated with DNA methylation changes. That is, individuals susceptible to ADHD respond differently compared to controls to long-term paracetamol exposure during development. This is compatible with the epidemiological evidence of the increased risk of developing ADHD with long-term paracetamol exposure. Further, these results lend novel insights into the etiology of ADHD and may serve as disease biomarkers in blood. Future developments should also investigate the functional impact of these DNA methylation differences on neural development in model systems.

Methods

Study population

The samples were obtained from MoBa, which is a large prospective population-based pregnancy cohort study conducted by the Norwegian Institute of Public Health [48]. Participants were recruited from all over Norway between 1999 and 2008. The women consented to participation in 40.6% of the pregnancies. The cohort now includes 114,500 children, 95,200 mothers, and 75,200 fathers. Blood samples were obtained from both parents during pregnancy and from mothers and children (umbilical cord) at birth from approximately 90,000 participants [49, 50]. This study is based on Data Version 8 released by MoBa in 2015. The MoBa cohort was linked to the Medical Birth Registry of Norway (MBRN) and the Norwegian Patient Registry (NPR) using the women’s personal 11-digit identification number.

Measures

Prenatal use of paracetamol

Information about paracetamol use (ATC code: N02BE01) was based on three questionnaires (Q1: from conception to gestational week 18, Q3: to gestational week 30, and Q4: to delivery). In order to minimize potential confounding by gestational age, only term children were included in the study (≥37 weeks). Samples with no reported co-medication together with paracetamol were selected to be certain that the number of days of use corresponded to paracetamol use. We defined long-term use as use of paracetamol for 20 days or more during pregnancy based on previous studies [18,19,20,21,22]. Moreover, trimester of paracetamol use was recorded and assessed in sub-analyses.

Measure of ADHD

Information about the children’s ADHD diagnoses from 2008 to 2014 was obtained from the NPR. From 2008, all hospitals owned and funded by the government and outpatient clinics report individual-level diagnoses given by specialists according to the 10th revision of the International Classification of Disease (ICD-10). Children registered with an ICD-10 diagnosis of hyperkinetic disorder (HKD) (F90.0, F90.1, F90.8, or F90.9) between 2008 and 2013 were identified and defined as having ADHD. HKD requires the combination of inattentive and hyperactive symptoms and thus corresponds to the combined subtype in the Diagnostic and Statistical Manual (DSM) system [51, 52].

Sampling and study design

The participants included in this study (n = 384) were selected to allow a cohort study design with 96 samples in each of the four categories (Table 1): (1) unexposed to paracetamol without ADHD (randomly selected control group), (2) exposed to paracetamol without ADHD (drug exposure group), (3) unexposed to paracetamol with ADHD (ADHD control group), and (4) exposed to paracetamol with ADHD (synergistic effect group). Selection criteria for paracetamol exposure and the clinical ADHD diagnose are described above. Potential covariates from MoBa and MBRN included infant sex, gestational age at delivery, maternal age, smoking, and alcohol consumption during pregnancy.

DNA methylation analysis

Microarray preprocessing and quality control

Bisulfite conversion of 500-ng cord blood DNA was done using the EZ-96 DNA Methylation-Gold Kit (Zymo Research) according to the manufacturer’s instructions. The samples were randomly located on 96-well plates to minimize potential batch effects related to bisulfite conversion. DNA methylation status was assessed using the Infinium HumanMethylation 450 BeadChip (Illumina). All analyses were carried out using the R programming language (http://www.r-project.org/), and the raw data were preprocessed using the approach implemented in RnBeads v.1.2.1 [53] in two steps (before and after normalization). First, cross-reactive probes [54] (n = 29,233), probes with overlapping SNPs in any of the bases in the target sequence (n = 41,930), and probes and samples with unreliable measurements (detection p values >0.01) (n = 11,680 probes and 7 samples) were removed. Background subtraction was performed using noob.methylumi [55], and β values (ratio of methylated signal divided by the sum of the methylated and unmethylated signal) were normalized using BMIQ [56]. Finally, probes located on the sex chromosomes (n = 9489) and non-CpG probes (n = 2458) were removed, resulting in a final data set consisting of 390,787 probes and 377 samples. Cell-type proportions (CD8+ and CD4+ T-lymphocytes, natural killer cells, B cells, monocytes, and granulocytes) for each sample were estimated using the Houseman reference-based approach [57] and a suitable validated cord blood reference data set (R package FlowSorted.CordBloodNorway.450K), which we have recently established [47].

Differential DNA methylation analyses

All statistical tests were performed on the M values (log2 of the β values), which are more statistically valid for differential DNA methylation analyses than β values [58]. To identify differentially methylated positions associated with paracetamol and/or ADHD, we fitted a linear regression model using the limma package [59] and the mean DNA methylation differences.

We performed three sets of analyses to explore the relationship between prenatal exposure to paracetamol and DNA methylation differences in children with ADHD (Table 1). First, we stratified the analyses on paracetamol exposure and ADHD separately to evaluate the associations with DNA methylation. Second, to investigate whether prenatal exposure to paracetamol is associated with differences in DNA methylation in children with ADHD, we compared exposed samples with ADHD (synergistic effect group) to controls. Theoretically, significant associations from this comparison reflect both associations with ADHD and paracetamol. Hence, we also performed comparisons of this group to unexposed samples with ADHD (ADHD control group) as well as exposed samples without ADHD (drug exposure group) for interpretation purposes. Third, we divided the exposed samples with ADHD (synergistic effect group) into sporadic (6–19 days, n = 77 samples) or long-term (≥20 days, n = 19 samples) groups and used similar comparison groups as above. This enabled exploration of the epidemiological evidence that long-term exposure during pregnancy is associated with development of ADHD [18,19,20,21,22].

All analyses were conducted on the CpG site and region levels. Regional analysis was performed on predefined genomic regions (5 kb tiles, genes, promoters and CpG islands). Briefly, aggregated p values for each region were obtained by combining the uncorrected p values for each CpG in the region [53]. To adjust for multiple testing, a false discovery rate (FDR) cutoff of less than 5% was used for genome-wide significance by using the method of Benjamini and Hochberg (BH) [60]. Sex, smoking, maternal age, gestational age, estimated cell-type composition, BeadChip, and bisulfite conversion plate were included as covariates in the statistical model.

Analysis was also done based on a combined score of statistical significance and effect size to increase the biological relevance of the observed differences in DNA methylation as implemented in RnBeads [53]. Briefly, (1) the absolute difference in group means, (2) the quotient in group means, and (3) the p values from the linear model were ranked for all sites. Then, the rank for each of the three measures was aggregated into a combined score. The differentially methylated CpGs and regions were ranked according to this combined score.

Surrogate variable analyses

Surrogate variable analyses were performed to search for unmeasured sources of variation in DNA methylation confounding DNA methylation differences between groups. Specifically, surrogate variables were estimated directly from the DNA methylation data using the sva R package [61] with default parameters and group comparisons as targets.

Gene ontology analysis

Gene ontology (GO) analysis was performed using the GOstats package based on the top 100 differentially methylated genes ranked according to the combined score (effect size). Both over- and under-representation of GO terms were tested with a standard hypergeometric test among genes showing differential DNA methylation.

Abbreviations

ADHD:

Attention-deficit/hyperactivity disorder

BH:

Benjamini and Hochberg

DSM:

Diagnostic and statistical manual

EWAS:

Epigenome-wide association study

FDR:

False discovery rate

GO:

Gene ontology

GWAS:

Genome-wide association study

HKD:

Hyperkinetic disorder

ICD-10:

10th revision of the international classification of disease

MBRN:

Medical birth registry of Norway

NPR:

Norwegian patient registry

SVA:

Surrogate variable analysis

References

  1. Polanczyk GV, Salum GA, Sugaya LS, Caye A, Rohde LA. Annual research review: a meta-analysis of the worldwide prevalence of mental disorders in children and adolescents. J Child Psychol Psychiatry. 2015;56:345–65.

    Article  PubMed  Google Scholar 

  2. Faraone SV, Perlis RH, Doyle AE, Smoller JW, Goralnick JJ, Holmgren MA, et al. Molecular genetics of attention-deficit/hyperactivity disorder. Biol Psychiatry. 2005;57:1313–23.

    Article  CAS  PubMed  Google Scholar 

  3. Lesch K-P, Timmesfeld N, Renner TJ, Halperin R, Röser C, Nguyen TT, et al. Molecular genetics of adult ADHD: converging evidence from genome-wide association and extended pedigree linkage studies. J Neural Transm Springer Vienna. 2008;115:1573–85.

    Article  CAS  Google Scholar 

  4. Neale BM, Su J, Anney R, Franke B, Zhou K, Maller JB, et al. Genome-wide association scan of attention deficit hyperactivity disorder. American journal of medical genetics. Part B, neuropsychiatric genetics: the official publication of the International Society of Psychiatric Genetics. NIH Public Access. 2008;147B:1337–44.

    CAS  Google Scholar 

  5. Mick E, Todorov A, Smalley S, Hu X, Loo S, Todd RD, et al. Family-based genome-wide association scan of attention-deficit/hyperactivity disorder. Journal of the American Academy of Child & Adolescent Psychiatry. 2010;49:898–905.e3.

    Article  Google Scholar 

  6. Neale BM, Medland S, Ripke S, Anney RJL, Asherson P, Buitelaar J, et al. Case-control genome-wide association study of attention-deficit/hyperactivity disorder. Journal of the American Academy of Child & Adolescent Psychiatry. 2010;49:906–20.

    Article  Google Scholar 

  7. Hinney A, Scherag A, Jarick I, Albayrak Ö, Pütter C, Pechlivanis S, et al. Genome-wide association study in German patients with attention deficit/hyperactivity disorder. American journal of medical genetics part B: neuropsychiatric genetics. Wiley Subscription Services, Inc A Wiley Company. 2011;156:888–97.

    Google Scholar 

  8. Stergiakouli E, Hamshere M, Holmans P, Langley K, Zaharieva I, deCODE Genetics, et al. Investigating the contribution of common genetic variants to the risk and pathogenesis of ADHD. AJP. 2012;169:186–94.

    Article  Google Scholar 

  9. Zayats T, Athanasiu L, Sønderby IE, Djurovic S, Westlye LT, Tamnes CK, et al. Genome-wide analysis of attention deficit hyperactivity disorder in Norway. Yao Y-G, editor. PLoS One. Public Library of Science (PLoS); 2015; 10: e0122501.

  10. Neale BM, Medland SE, Ripke S, Asherson P, Franke B, Lesch K-P, et al. Meta-analysis of genome-wide association studies of attention-deficit/hyperactivity disorder. Journal of the American Academy of Child & Adolescent Psychiatry. 2010;49:884–97.

    Article  Google Scholar 

  11. Akutagava-Martins GC, Rohde LA, Hutz MH. Genetics of attention-deficit/hyperactivity disorder: an update. Expert Review of Neurotherapeutics Taylor & Francis. 2016;16:145–56.

    Article  CAS  Google Scholar 

  12. Thapar A, Cooper M. Attention deficit hyperactivity disorder. Lancet; . Elsevier 2016; 387:1240–1250.

  13. van Mil NH, Steegers-Theunissen RPM, Bouwland-Both MI, Verbiest MMPJ, Rijlaarsdam J, Hofman A, et al. DNA methylation profiles at birth and child ADHD symptoms. J Psychiatr Res. 2014;49:51–9.

    Article  PubMed  Google Scholar 

  14. Xu Y, Chen X-T, Luo M, Tang Y, Zhang G, Wu D, et al. Multiple epigenetic factors predict the attention deficit/hyperactivity disorder among the Chinese Han children. J Psychiatr Res. 2015;64:40–50.

    Article  PubMed  Google Scholar 

  15. Wilmot B, Fry R, Smeester L, Musser ED, Mill J, Nigg JT. Methylomic analysis of salivary DNA in childhood ADHD identifies altered DNA methylation in VIPR2. J Child Psychol Psychiatry. 2015;57:152–60.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Park S, Lee JM, Kim JW, Cho DY, Yun HJ, Han DH, et al. Associations between serotonin transporter gene (SLC6A4) methylation and clinical characteristics and cortical thickness in children with ADHD. Psychol Med. 2015;45:3009–17.

    Article  CAS  PubMed  Google Scholar 

  17. Walton E, Pingault JB, Cecil CAM, Gaunt TR, Relton CL, Mill J, et al. Epigenetic profiling of ADHD symptoms trajectories: a prospective, methylome-wide study. Mol Psychiatry. 2017;2:250–6.

    Article  Google Scholar 

  18. Liew Z, Ritz B, Rebordosa C, Lee P-C, Olsen J. Acetaminophen use during pregnancy, behavioral problems, and hyperkinetic disorders. JAMA Pediatr. 2014;168:313–20.

    Article  PubMed  Google Scholar 

  19. Thompson JMD, Waldie KE, Wall CR, Murphy R, Mitchell EA, the ABC study group. Associations between acetaminophen use during pregnancy and ADHD symptoms measured at ages 7 and 11 years. Hashimoto K, editor. PLoS One. 2014;9:e108210.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Avella-Garcia CB, Julvez J, Fortuny J, Rebordosa C, García-Esteban R, Galán IR, et al. Acetaminophen use in pregnancy and neurodevelopment: attention function and autism spectrum symptoms. Int J Epidemiol; Oxford University Press 2016; 6:1987-1996.

  21. Stergiakouli E, Thapar A, Smith GD. Association of acetaminophen use during pregnancy with behavioral problems in childhood: evidence against confounding. JAMA Pediatr. 2016;10:964–70.

    Article  Google Scholar 

  22. Brandlistuen RE, Ystrom E, Nulman I, Koren G, Nordeng H. Prenatal paracetamol exposure and child neurodevelopment: a sibling-controlled cohort study. Int J Epidemiol. 2014;42:1702–13.

    Article  Google Scholar 

  23. Lupattelli A, Spigset O, Twigg MJ, Zagorodnikova K, Mårdby AC, Moretti ME, et al. Medication use in pregnancy: a cross-sectional, multinational web-based study. BMJ Open. 2014;4:e004365.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Elia J, Glessner JT, Wang K, Takahashi N, Shtir CJ, Hadley D, et al. Genome-wide copy number variation study associates metabotropic glutamate receptor gene networks with attention deficit hyperactivity disorder. - PubMed-NCBI. Nat Genet. 2011;44:78–84.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Kitagishi Y, Minami A, Nakanishi A, Ogura Y, Matsuda S. Neuron membrane trafficking and protein kinases involved in autism and ADHD. International journal of molecular sciences. Multidisciplinary Digital Publishing Institute (MDPI). 2015;16:3095–115.

    CAS  Google Scholar 

  26. Ballas N, Grunseich C, Lu DD, Speh JC, Mandel G. REST and its corepressors mediate plasticity of neuronal gene chromatin throughout neurogenesis. Cell; Elsevier 2005; 121:645–657.

  27. Wakatsuki S, Furuno A, Ohshima M, Araki T. Oxidative stress-dependent phosphorylation activates ZNRF1 to induce neuronal/axonal degeneration. J Cell Biol; Rockefeller Univ Press 2015; 211:881–896.

  28. Tam S-Y, Lilla JN, Chen C-C, Kalesnikoff J, Tsai M. RabGEF1/Rabex-5 regulates TrkA-mediated neurite outgrowth and NMDA-induced Signaling activation in NGF-differentiated PC12 cells. Obukhov AG, editor. PLoS One; Public Library of Science 2015; 10: e0142935.

  29. Cubelos B, Sebastián-Serrano A, Beccari L, Calcagnotto ME, Cisneros E, Kim S, et al. Cux1 and Cux2 regulate dendritic branching, spine morphology, and synapses of the upper layer neurons of the cortex. Neuron; Elsevier 2010; 66:523–535.

  30. Medrihan L, Cesca F, Raimondi A, Lignani G, Baldelli P, Benfenati F. Synapsin II desynchronizes neurotransmitter release at inhibitory synapses by interacting with presynaptic calcium channels. Nat Comms. 2013;4:1512.

    Article  Google Scholar 

  31. Calabrese V, Mancuso C, Calvani M, Rizzarelli E, Butterfield DA, Stella AMG. Nitric oxide in the central nervous system: neuroprotection versus neurotoxicity. Nat Rev Neurosci. 2007;8:766–75.

    Article  CAS  PubMed  Google Scholar 

  32. Imbrici P, Camerino DC, Tricarico D. Major channels involved in neuropsychiatric disorders and therapeutic perspectives. Front Genet Frontiers Media SA. 2013;4:76.

    Google Scholar 

  33. Philippot G, Nyberg F, Gordh T, Fredriksson A, Viberg H. Short-term exposure and long-term consequences of neonatal exposure to Δ9-tetrahydrocannabinol (THC) and ibuprofen in mice. Behav Brain Res. 2016;307:137–44.

    Article  CAS  PubMed  Google Scholar 

  34. Viberg H, Eriksson P, Gordh T, Fredriksson A. Paracetamol (acetaminophen) administration during neonatal brain development affects cognitive function and alters its analgesic and anxiolytic response in adult male mice. Toxicol Sci; Oxford University Press 2014; 138:139–147.

  35. Mooney MA, McWeeney SK, Faraone SV, Hinney A, Hebebrand J, IMAGE2 Consortium, et al. Pathway analysis in attention deficit hyperactivity disorder: an ensemble approach. Am J Med Genet B Neuropsychiatr Genet. 2016;6:815–26.

    Article  Google Scholar 

  36. Guney E, Cetin FH, Alisik M, Tunca H, Tas Torun Y, Iseri E, et al. Attention deficit hyperactivity disorder and oxidative stress: a short term follow up study. Psychiatry Res. 2015;229:310–7.

    Article  CAS  PubMed  Google Scholar 

  37. Archana E, Pai P, Prabhu BK, Shenoy RP, Prabhu K, Rao A. Altered biochemical parameters in saliva of pediatric attention deficit hyperactivity disorder. Neurochem Res Springer US. 2011;37:330–4.

    Article  Google Scholar 

  38. Ceylan MF, Sener S, Bayraktar AC, Kavutcu M. Changes in oxidative stress and cellular immunity serum markers in attention-deficit/hyperactivity disorder. Psychiatry Clin Neurosci. 2012;66:220–6.

    Article  CAS  PubMed  Google Scholar 

  39. Goldani AAS, Downs SR, Widjaja F, Lawton B, Hendren RL. Biomarkers in Autism. Frontiers in Psychiatry. Frontiers Media SA; 2014; 5:950.

  40. Popa-Wagner A, Mitran S, Sivanesan S, Chang E, Buga A-M. ROS and brain diseases: the good, the bad, and the ugly. Oxidative Med Cell Longev. 2013;2013:1–14.

    Article  Google Scholar 

  41. Masutani H. Oxidative stress and redox imbalance in acetaminophen toxicity. The Pharmacogenomics Journal. 2001;1:165–6.

    Article  CAS  PubMed  Google Scholar 

  42. Ghanizadeh A, Bahrani M, Miri R, Sahraian A. Smell identification function in children with attention deficit hyperactivity disorder. Psychiatry Investigation Korean Neuropsychiatric Association. 2012;9:150–3.

    Google Scholar 

  43. Romanos M, Renner TJ, Schecklmann M, Hummel B, Roos M, Mering von C, et al. Improved odor sensitivity in attention-deficit/hyperactivity disorder. Biol Psychiatry. 2008;64:938–40.

    Article  CAS  PubMed  Google Scholar 

  44. Weiland R, Macht M, Ellgring H, Groß-Lesch S, Lesch K-P, Pauli P. Olfactory and gustatory sensitivity in adults with attention-deficit/hyperactivity disorder. Atten Defic Hyperact Disord. 2011;3:53–60.

    Article  PubMed  Google Scholar 

  45. Ystrom E, Vollrath ME, Nordeng H. Effects of personality on use of medications, alcohol, and cigarettes during pregnancy. Eur J Clin Pharmacol; Springer-Verlag 2011; 68:845–851.

  46. Kalda A, Zharkovsky A. Epigenetic mechanisms of Psychostimulant-induced addiction. Int Rev Neurobiol. 2015:85–105.

  47. Gervin K, Page CM, Aass HCD, Jansen MA, Fjeldstad HE, Andreassen BK, et al. Cell type specific DNA methylation in cord blood: a 450K-reference data set and cell count-based validation of estimated cell type composition. Epigenetics. 2016;9:690–8.

    Article  Google Scholar 

  48. Magnus P, Birke C, Vejrup K, Haugan A, Alsaker E, Daltveit AK, et al. Cohort profile update: the Norwegian mother and child cohort study (MoBa). Int J Epidemiol; . Oxford University Press 2016:1-7.

  49. Rønningen KS, Paltiel L, Meltzer HM, Nordhagen R, Lie KK, Hovengen R, et al. The biobank of the Norwegian mother and child cohort study: a resource for the next 100 years. Eur J Epidemiol. 2006;21:619–25.

    Article  PubMed  PubMed Central  Google Scholar 

  50. Paltiel L, Haugan A, Skjerden T, Harbak K, Bækken S, Stensrud NK, et al. The biobank of the Norwegian mother and child cohort study—present status. Norsk Epidemiologi. 2014;24:29–35.

  51. Taylor E, Schachar R, Thorley G, Wieselberg HM, Everitt B, Rutter M. Which boys respond to stimulant medication? A controlled trial of methylphenidate in boys with disruptive behaviour. Psychol Med; Cambridge University Press 1987; 17:121–143.

  52. Thapar A, Pine DS, Leckman JF, Scott S, Snowling MJ, Taylor E. Rutter’s child and adolescent psychiatry. Chichester, UK: Wiley; 2015.

  53. Assenov Y, Müller F, Lutsik P, Walter J, Lengauer T, Bock C. Comprehensive analysis of DNA methylation data with RnBeads. Nat Meth Nature Publishing Group. 2014;11:1138–40.

    Article  CAS  Google Scholar 

  54. Chen Y-A, 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. 2014;8:203–9.

    Article  Google Scholar 

  55. 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:e90.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  56. Teschendorff AE, Marabita F, Lechner M, Bartlett T, Tegner J, Gomez-Cabrero D, et al. A beta-mixture quantile normalization method for correcting probe design bias in Illumina Infinium 450 k DNA methylation data. Bioinformatics. 2013;29:189–96.

    Article  CAS  PubMed  Google Scholar 

  57. Houseman E, Accomando WP, Koestler DC, Christensen BC, Marsit CJ, Nelson HH, et al. DNA methylation arrays as surrogate measures of cell mixture distribution. BMC Bioinformatics. 2012;13:86.

    Article  PubMed  PubMed Central  Google Scholar 

  58. Du P, Zhang X, Huang C-C, Jafari N, Kibbe WA, Hou L, et al. Comparison of Beta-value and M-value methods for quantifying methylation levels by microarray analysis. BMC Bioinformatics BioMed Central. 2010;11:587.

    Article  CAS  Google Scholar 

  59. 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:e47.

    Article  PubMed  PubMed Central  Google Scholar 

  60. Benjamini Y. Discovering the false discovery rate. Journal of the Royal Statistical Society: series B (statistical methodology); Blackwell Publishing Ltd; 2010; 72:405–416.

  61. Leek JT, Johnson WE, Parker HS, Jaffe AE, Storey JD. The sva package for removing batch effects and other unwanted variation in high-throughput experiments. Bioinformatics; . Oxford University Press 2012; 28:882–883.

Download references

Acknowledgements

The Infinium 450K HumanMethylation BeadChips were processed at the Norwegian Microarray Consortium (Oslo University Hospital, Radium). We are grateful to all the participating families in Norway who take part in this on-going MoBa cohort study.

Funding

The study was funded by the Norwegian Research Council (grant number 32251), the Southern and Eastern Norway Regional Health Authority (grant number 39671), and the European Research Council Starting Grant (grant number 678033). The Norwegian Mother and Child Cohort Study is supported by the Norwegian Ministry of Health and the Ministry of Education and Research, NIH/NIEHS (contract number N01-ES-75558), and NIH/NINDS (grant number 1: UO1 NS 047537–01 and grant number 2: UO1 NS 047537-06A1).

Availability of data and materials

The data that support the findings in this study are available from the Norwegian Mother and Child Cohort, but restrictions apply to the availability of these data, which were used under license for the current study and so are not publicly available. Data are however available from the authors upon reasonable request and with permission from the Norwegian Mother and Child Cohort.

Authors’ contributions

KG, HN, and RL conceived and designed the experiments. KG performed the experiments. KG and RL performed statistical analysis. KG, HN, EY, TR-K, and RL analyzed the data. KG, HN, EY, TR-K, and RL wrote the paper. All authors read and approved the final manuscript.

Author information

Authors and Affiliations

Authors

Corresponding authors

Correspondence to Kristina Gervin or Robert Lyle.

Ethics declarations

Ethics approval and consent to participate

This work was approved by the National Institute of Public Health and Southern and Eastern Norway Regional Health Authority. Written informed consent was obtained from parents of all subjects in MoBa.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1: Table S1.

Significantly methylated CpGs between long-term exposed samples with ADHD and controls. (XLSX 1646 kb)

Additional file 2: Table S2.

Significantly methylated CpGs between long-term exposed samples with ADHD and sporadically exposed samples with ADHD. (XLSX 585 kb)

Additional file 3: Table S3.

Significantly methylated CpGs between long-term exposed samples with ADHD and unexposed samples with ADHD. (XLSX 101 kb)

Additional file 4: Figure S1.

Enrichment of small p values associated with differences in DNA methylation in long-term exposed children with ADHD. Q-Q plots of the observed versus expected p values from the comparisons DNA methylation in cord blood from long-term exposed children with ADHD (synergy group, ≥20 days) to controls (blue, lambda = 2.2), sporadically exposed children with ADHD (red, lambda = 1.8) and unexposed children with ADHD (green, lambda = 1.7). The corresponding histograms of the nominal p values for each of the three comparisons display an enrichment of small p values. (PDF 887 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Gervin, K., Nordeng, H., Ystrom, E. et al. Long-term prenatal exposure to paracetamol is associated with DNA methylation differences in children diagnosed with ADHD. Clin Epigenet 9, 77 (2017). https://doi.org/10.1186/s13148-017-0376-9

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13148-017-0376-9

Keywords