Multi-omics approach reveals dysregulated genes during hESCs neuronal differentiation exposure to paracetamol.

Prenatal paracetamol exposure has been associated with neurodevelopmental outcomes in childhood. Pharmacoepigenetic studies show differences in cord blood DNA methylation between unexposed and paracetamol-exposed neonates, however, causality and impact of long-term prenatal paracetamol exposure on brain development remain unclear. Using a multi-omics approach, we investigated the effects of paracetamol on an in vitro model of early human neurodevelopment. We exposed human embryonic stem cells undergoing neuronal differentiation with paracetamol concentrations corresponding to maternal therapeutic doses. Single-cell RNA-seq and ATAC-seq integration identified paracetamol-induced chromatin opening changes linked to gene expression. Differentially methylated and/or expressed genes were involved in neurotransmission and cell fate determination trajectories. Some genes involved in neuronal injury and development-specific pathways, such as KCNE3, overlapped with differentially methylated genes previously identified in cord blood associated with prenatal paracetamol exposure. Our data suggest that paracetamol may play a causal role in impaired neurodevelopment.

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.

Long-Term Use of Amoxicillin Is Associated with Changes in Gene Expression and DNA Methylation in Patients with Low Back Pain and Modic Changes.

Long-term antibiotics are prescribed for a variety of medical conditions, recently including low back pain with Modic changes. The molecular impact of such treatment is unknown. We conducted longitudinal transcriptome and epigenome analyses in patients (n = 100) receiving amoxicillin treatment or placebo for 100 days in the Antibiotics in Modic Changes (AIM) study. Gene expression and DNA methylation were investigated at a genome-wide level at screening, after 100 days of treatment, and at one-year follow-up. We identified intra-individual longitudinal changes in gene expression and DNA methylation in patients receiving amoxicillin, while few changes were observed in patients receiving placebo. After 100 days of amoxicillin treatment, 28 genes were significantly differentially expressed, including the downregulation of 19 immunoglobulin genes. At one-year follow-up, the expression levels were still not completely restored. The significant changes in DNA methylation (n = 4548 CpGs) were mainly increased methylation levels between 100 days and one-year follow-up. Hence, the effects on gene expression occurred predominantly during treatment, while the effects on DNA methylation occurred after treatment. In conclusion, unrecognized side effects of long-term amoxicillin treatment were revealed, as alterations were observed in both gene expression and DNA methylation that lasted long after the end of treatment.

Does maternal genetic liability to folate deficiency influence the risk of antiseizure medication-associated language impairment and autistic traits in children of women with epilepsy?

BACKGROUND: Prenatal exposure to antiseizure medication (ASM) may lead to low plasma folate concentrations and is associated with impaired neurodevelopment. OBJECTIVES: To examine whether maternal genetic liability to folate deficiency interacts with ASM-associated risk of language impairment and autistic traits in children of women with epilepsy. METHODS: We included children of women with and without epilepsy and with available genetic data enrolled in the Norwegian Mother, Father, and Child Cohort Study. Information on ASM use, folic acid supplement use and dose, dietary folate intake, child autistic traits, and child language impairment was obtained from parent-reported questionnaires. Using logistic regression, we examined the interaction between prenatal ASM exposure and maternal genetic liability to folate deficiency expressed as polygenic risk score of low folate concentrations or maternal rs1801133 genotype (CC or CT/TT) on risk of language impairment or autistic traits. RESULTS: We included 96 children of women with ASM-treated epilepsy, 131 children of women with ASM-untreated epilepsy, and 37,249 children of women without epilepsy. The polygenic risk score of low folate concentrations did not interact with the ASM-associated risk of language impairment or autistic traits in ASM-exposed children of women with epilepsy compared with ASM-unexposed children aged 1.5-8 y. ASM-exposed children had increased risk of adverse neurodevelopment regardless of maternal rs1801133 genotype {adjusted odds ratio [aOR] for language impairment aged 8 y was 2.88 [95% confidence interval (CI): 1.00, 8.26] if CC and aOR 2.88 [95% CI: 1.10, 7.53] if CT/TT genotypes}. In children of women without epilepsy aged 3 y, those with maternal rs1801133 CT/TT compared with CC genotype had increased risk of language impairment (aOR: 1.18; 95% CI: 1.05, 1.34). CONCLUSIONS: In this cohort of pregnant women reporting widespread use of folic acid supplements, maternal genetic liability to folate deficiency did not significantly influence the ASM-associated risk of impaired neurodevelopment.

Cytokine Patterns as Predictors of Antibiotic Treatment Effect in Chronic Low Back Pain with Modic Changes: Subgroup Analyses of a Randomized Trial (AIM Study).

Objective: Randomized trials testing the effect of antibiotics for chronic low back pain (LBP) with vertebral bone marrow changes on MRI (Modic changes) report inconsistent results. A proposed explanation is subgroups with low grade discitis where antibiotics are effective, but there is currently no method to identify such subgroups. The objective of the present study was to evaluate whether distinct patterns of serum cytokine levels predict any treatment effect of oral amoxicillin at one-year follow-up in patients with chronic low back pain and Modic changes at the level of a previous lumbar disc herniation. Design: We used data from an overpowered, randomized, placebo-controlled trial (the AIM study) that tested 100 days of oral 750 mg amoxicillin vs placebo three times daily in hospital outpatients with chronic (>6 months) LBP with pain intensity ≥5 on a 0-10 numerical rating scale and Modic changes type 1 (oedema type) or 2 (fatty type). We measured serum levels of 40 inflammatory cytokines at baseline and analysed six predefined potential predictors of treatment effect based on cytokine patterns in 78 randomized patients; three analyses with recursive partitioning, one based on cluster analysis and two based on principal component analyses. The primary outcome was the Roland-Morris Disability Questionnaire score at one-year follow-up in the intention to treat population. The methodology and overall results of the AIM study were published previously. Results: The 78 patients were 25-62 years old and 47 (60%) were women. None of the three recursive partitioning analyses resulted in any suggested subgroups. Of all main analyses, the largest effect estimate (mean difference between antibiotic and placebo groups) was seen in a subgroup not predefined as of main interest (Cluster category 3+4; -2.0, 95% CI: -5.2-1.3, RMDQ points; p-value for interaction 0.54). Conclusion: Patterns of inflammatory serum cytokine levels did not predict treatment effect of amoxicillin in patients with chronic LBP and Modic changes. Clinical Trial Registration Number: ClinicalTrials.gov (identifier: NCT02323412).

Effects of prenatal exposure to (es)citalopram and maternal depression during pregnancy on DNA methylation and child neurodevelopment.

Studies assessing associations between prenatal exposure to antidepressants, maternal depression, and offspring DNA methylation (DNAm) have been inconsistent. Here, we investigated whether prenatal exposure to citalopram or escitalopram ((es)citalopram) and maternal depression is associated with differences in DNAm. Then, we examined if there is an interaction effect of (es)citalopram exposure and DNAm on offspring neurodevelopmental outcomes. Finally, we investigated whether DNAm at birth correlates with neurodevelopmental trajectories in childhood. We analyzed DNAm in cord blood from the Norwegian Mother, Father and Child Cohort Study (MoBa) biobank. MoBa contains questionnaire data on maternal (es)citalopram use and depression during pregnancy and information about child neurodevelopmental outcomes assessed by internationally recognized psychometric tests. In addition, we retrieved ADHD diagnoses from the Norwegian Patient Registry and information on pregnancies from the Medical Birth Registry of Norway. In total, 958 newborn cord blood samples were divided into three groups: (1) prenatal (es)citalopram exposed (n = 306), (2) prenatal maternal depression exposed (n = 308), and (3) propensity score-selected controls (n = 344). Among children exposed to (es)citalopram, there were more ADHD diagnoses and symptoms and delayed communication and psychomotor development. We did not identify differential DNAm associated with (es)citalopram or depression, nor any interaction effects on neurodevelopmental outcomes throughout childhood. Trajectory modeling identified subgroups of children following similar developmental patterns. Some of these subgroups were enriched for children exposed to maternal depression, and some subgroups were associated with differences in DNAm at birth. Interestingly, several of the differentially methylated genes are involved in neuronal processes and development. These results suggest DNAm as a potential predictive molecular marker of later abnormal neurodevelopmental outcomes, but we cannot conclude whether DNAm links prenatal (es)citalopram exposure or maternal depression with child neurodevelopmental outcomes.

Decreased serum concentrations of antiseizure medications in children with drug resistant epilepsy following treatment with ketogenic diet.

OBJECTIVE: To examine the potential influence of a ketogenic diet on serum concentrations of antiseizure medications (ASMs) in children with drug resistant epilepsy. METHODS: We investigated the serum concentrations of ASMs in 25 children with drug resistant epilepsy, 2-13 years of age, treated with a classical ketogenic diet for 12 weeks. The patients were recruited from the National Centre for Epilepsy from August 15th, 2017, to January 24th, 2022. Changes in ASM serum concentrations were analyzed using a mixed effect model analysis. Significance level was set at P < 0.05 for all comparisons. RESULTS: The participants used 12 different ASMs during the study. The mean number of ASMs was 2.4 (±SD 0.7). None of the participants changed the type or dose of the ASMs during the intervention period. The serum concentrations of clobazam (n = 9, P = 0.002), desmethylclobazam (n = 9, P = 0.010), and lamotrigine (n = 6, P = 0.016) decreased significantly during the dietary treatment. The analytes with the largest reduction in serum concentration after 12 weeks of dietary treatment were clobazam (mean change -38%) and desmethylclobazam (mean change -37%). We found no significant change in the serum concentrations of levetiracetam, topiramate, and valproic acid. SIGNIFICANCE: We identified a significant decrease in the serum concentrations of clobazam, desmethylclobazam, and lamotrigine following a 12-week ketogenic diet intervention in children with drug resistant epilepsy. An unintended decrease in the serum concentrations of ASMs may render the patient prone to seizures. Measurements of ASM serum concentrations might be useful in patients on a ketogenic diet, especially in patients with lack of efficacy of the dietary treatment.

A multi-omics approach to visualize early neuronal differentiation from hESCs in 4D.

Neuronal differentiation of pluripotent stem cells is an established method to study physiology, disease, and medication safety. However, the sequence of events in human neuronal differentiation and the ability of in vitro models to recapitulate early brain development are poorly understood. We developed a protocol optimized for the study of early human brain development and neuropharmacological applications. We comprehensively characterized gene expression and epigenetic profiles at four timepoints, because the cells differentiate from embryonic stem cells towards a heterogeneous population of progenitors, immature and mature neurons bearing telencephalic signatures. A multi-omics roadmap of neuronal differentiation, combined with searchable interactive gene analysis tools, allows for extensive exploration of early neuronal development and the effect of medications.

Low reliability of DNA methylation across Illumina Infinium platforms in cord blood: implications for replication studies and meta-analyses of prenatal exposures.

BACKGROUND: There is an increasing interest in the role of epigenetics in epidemiology, but the emerging research field faces several critical biological and technical challenges. In particular, recent studies have shown poor correlation of measured DNA methylation (DNAm) levels within and across Illumina Infinium platforms in various tissues. In this study, we have investigated concordance between 450 k and EPIC Infinium platforms in cord blood. We could not replicate our previous findings on the association of prenatal paracetamol exposure with cord blood DNAm, which prompted an investigation of cross-platform DNAm differences. RESULTS: This study is based on two DNAm data sets from cord blood samples selected from the Norwegian Mother, Father and Child Cohort Study (MoBa). DNAm of one data set was measured using the 450 k platform and the other data set was measured using the EPIC platform. Initial analyses of the EPIC data could not replicate any of our previous significant findings in the 450 k data on associations between prenatal paracetamol exposure and cord blood DNAm. A subset of the samples (n = 17) was included in both data sets, which enabled analyses of technical sources potentially contributing to the negative replication. Analyses of these 17 samples with repeated measurements revealed high per-sample correlations ([Formula: see text] 0.99), but low per-CpG correlations ([Formula: see text] ≈ 0.24) between the platforms. 1.7% of the CpGs exhibited a mean DNAm difference across platforms > 0.1. Furthermore, only 26.7% of the CpGs exhibited a moderate or better cross-platform reliability (intra-class correlation coefficient ≥ 0.5). CONCLUSION: The observations of low cross-platform probe correlation and reliability corroborate previous reports in other tissues. Our study cannot determine the origin of the differences between platforms. Nevertheless, it emulates the setting in studies using data from multiple Infinium platforms, often analysed several years apart. Therefore, the findings may have important implications for future epigenome-wide association studies (EWASs), in replication, meta-analyses and longitudinal studies. Cognisance and transparency of the challenges related to cross-platform studies may enhance the interpretation, replicability and validity of EWAS results both in cord blood and other tissues, ultimately improving the clinical relevance of epigenetic epidemiology.

Genome-wide decrease in DNA methylation in adults with epilepsy treated with modified ketogenic diet: A prospective study.

OBJECTIVE: The aim of this study was to investigate the impact of the modified ketogenic diet on DNA methylation in adults with epilepsy. METHODS: In this prospective study, we investigated the genome-wide DNA methylation in whole blood in 58 adults with epilepsy treated with the modified ketogenic for 12 weeks. Patients were recruited from the National Center for Epilepsy, Norway, from March 1, 2011 to February 28, 2017. DNA methylation was analyzed using the Illumina Infinium MethylationEPIC BeadChip array. Analysis of variance and paired t-test were used to identify differentially methylated loci after 4 and 12 weeks of dietary treatment. A false discovery rate approach with a significance threshold of <5% was used to adjust for multiple comparisons. RESULTS: We observed a genome-wide decrease in DNA methylation, both globally and at specific sites, after 4 and 12 weeks of dietary treatment. A substantial share of the differentially methylated positions (CpGs) were annotated to genes associated with epilepsy (n = 7), lipid metabolism (n = 8), and transcriptional regulation (n = 10). Furthermore, five of the identified genes were related to inositol phosphate metabolism, which may represent a possible mechanism by which the ketogenic diet attenuates seizures. SIGNIFICANCE: A better understanding of the modified ketogenic diet's influence at the molecular level may be the key to unraveling the mechanisms by which the diet can ameliorate seizures and possibly to identifying novel therapeutic targets for epilepsy.

Prenatal medication exposure and epigenetic outcomes: a systematic literature review and recommendations for prenatal pharmacoepigenetic studies.

When used during pregnancy, analgesics and psychotropics pass the placenta to enter the foetal circulation and may induce epigenetic modifications. Where such modifications occur and whether they disrupt normal foetal developme nt, are currently unanswered questions. This field of prenatal pharmacoepigenetics has received increasing attention, with several studies reporting associations between in utero medication exposure and offspring epigenetic outcomes. Nevertheless, no recent systematic review of the literature is available. Therefore, the objectives of this review were to (i) provide an overview of the literature on the association of prenatal exposure to psychotropics a nd analgesics with epigenetic outcomes, and (ii) suggest recommendations for future studies within prenatal pharmacoepigenetics. We performed systematic literature searches in five databases. The eligible studies assessed human prenatal exposure to psychotropics or analgesics, with epigenetic analyses of offspring tissue as an outcome. We identified 18 eligible studies including 4,419 neonates exposed to either antidepressants, antiepileptic drugs, paracetamol, acetylsalicylic acid, or methadone. The epigenetic outcome in all studies was DNA methylation in cord blood, placental tissue or buccal cells. Although most studies found significant differences in DNA methylation upon medication exposure, almost no differences were persistent across studies for similar medications and sequencing methods. The reviewed studies were challenging to compare due to poor transparency in reporting, and heterogeneous methodology, design, genome coverage, and statistical modelling. We propose 10 recommendations for future prenatal pharmacoepigenetic studies considering both epidemiological and epigenetic perspectives. These recommendations may improve the quality, comparability, and clinical relevance of such studies. PROSPERO registration ID: CRD42020166675.

Modeling dependency structures in 450k DNA methylation data.

MOTIVATION: DNA methylation has been shown to be spatially dependent across chromosomes. Previous studies have focused on the influence of genomic context on the dependency structure, while not considering differences in dependency structure between individuals. RESULTS: We modeled spatial dependency with a flexible framework to quantify the dependency structure, focusing on inter-individual differences by exploring the association between dependency parameters and technical and biological variables. The model was applied to a subset of the Finnish Twin Cohort study (N = 1611 individuals). The estimates of the dependency parameters varied considerably across individuals, but were generally consistent across chromosomes within individuals. The variation in dependency parameters was associated with bisulfite conversion plate, zygosity, sex and age. The age differences presumably reflect accumulated environmental exposures and/or accumulated small methylation differences caused by stochastic mitotic events, establishing recognizable, individual patterns more strongly seen in older individuals. AVAILABILITY AND IMPLEMENTATION: The twin dataset used in the current study are located in the Biobank of the National Institute for Health and Welfare, Finland. All the biobanked data are publicly available for use by qualified researchers following a standardized application procedure (https://thl.fi/en/web/thl-biobank/for-researchers). A R-script for fitting the dependency structure to publicly available DNA methylation data with the software used in this article is provided in supplementary data. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Methotrexate Treatment of Newly Diagnosed RA Patients Is Associated With DNA Methylation Differences at Genes Relevant for Disease Pathogenesis and Pharmacological Action.

Background: Methotrexate (MTX) is the first line treatment of rheumatoid arthritis (RA), and methylation changes in bulk T cells have been reported after treatment with MTX. We have investigated cell-type specific DNA methylation changes across the genome in naïve and memory CD4+ T cells before and after MTX treatment of RA patients. DNA methylation profiles of newly diagnosed RA patients (N=9) were assessed by reduced representation bisulfite sequencing. Results: We found that MTX treatment significantly influenced DNA methylation levels at multiple CpG sites in both cell populations. Interestingly, we identified differentially methylated sites annotated to two genes; TRIM15 and SORC2, previously reported to predict treatment outcome in RA patients when measured in bulk T cells. Furthermore, several of the genes, including STAT3, annotated to the significant CpG sites are relevant for RA susceptibility or the action of MTX. Conclusion: We detected CpG sites that were associated with MTX treatment in CD4+ naïve and memory T cells isolated from RA patients. Several of these sites overlap genetic regions previously associated with RA risk and MTX treatment outcome.

Macrophage migration inhibitory factor: a potential biomarker for chronic low back pain in patients with Modic changes.

BACKGROUND: Low back pain (LBP) is a leading cause of disability worldwide, but the aetiology remains poorly understood. Finding relevant biomarkers may lead to better understanding of disease mechanisms. Patients with vertebral endplate bone marrow lesions visualised on MRI as Modic changes (MCs) have been proposed as a distinct LBP phenotype, and inflammatory mediators may be involved in the development of MCs. OBJECTIVES: To identify possible serum biomarkers for LBP in patients with MCs. METHODS: In this case control study serum levels of 40 cytokines were compared between patients with LBP and MC type 1 (n=46) or type 2 (n=37) and healthy controls (n=50). RESULTS: Analyses identified significantly higher levels of six out of 40 cytokines in the MC type 1 group (MC1), and five in the MC type 2 group (MC2) compared with healthy controls. Six cytokines were moderately correlated with pain. Principal component analyses revealed clustering and separation of patients with LBP and controls, capturing 40.8% of the total variance, with 10 cytokines contributing to the separation. Macrophage migration inhibitory factor (MIF) alone accounted for 92% of the total contribution. Further, receiver operating characteristics analysis revealed that MIF showed an acceptable ability to distinguish between patients and controls (area under the curve=0.79). CONCLUSIONS: These results suggest that cytokines may play a role in LBP with MCs. The clinical significance of the findings is unknown. MIF strongly contributed to clustering of patients with LBP with MCs and controls, and might be a biomarker for MCs. Ultimately, these results may guide future research on novel treatments for this patient group.

Genetic Susceptibility to Drug Teratogenicity: A Systematic Literature Review.

Since the 1960s, drugs have been known to cause teratogenic effects in humans. Such teratogenicity has been postulated to be influenced by genetics. The aim of this review was to provide an overview of the current knowledge on genetic susceptibility to drug teratogenicity in humans and reflect on future directions within the field of genetic teratology. We focused on 12 drugs and drug classes with evidence of teratogenic action, as well as 29 drugs and drug classes with conflicting evidence of fetal safety in humans. An extensive literature search was performed in the PubMed and EMBASE databases using terms related to the drugs of interest, congenital anomalies and fetal development abnormalities, and genetic variation and susceptibility. A total of 29 studies were included in the final data extraction. The eligible studies were published between 1999 and 2020 in 10 different countries, and comprised 28 candidate gene and 1 whole-exome sequencing studies. The sample sizes ranged from 20 to 9,774 individuals. Several drugs were investigated, including antidepressants (nine studies), thalidomide (seven studies), antiepileptic drugs (five studies), glucocorticoids (four studies), acetaminophen (two studies), and sex hormones (estrogens, one study; 17-alpha hydroxyprogesterone caproate, one study). The main neonatal phenotypic outcomes included perinatal complications, cardiovascular congenital anomalies, and neurodevelopmental outcomes. The review demonstrated that studies on genetic teratology are generally small, heterogeneous, and exhibit inconsistent results. The most convincing findings were genetic variants in SLC6A4, MTHFR, and NR3C1, which were associated with drug teratogenicity by antidepressants, antiepileptics, and glucocorticoids, respectively. Notably, this review demonstrated the large knowledge gap regarding genetic susceptibility to drug teratogenicity, emphasizing the need for further efforts in the field. Future studies may be improved by increasing the sample size and applying genome-wide approaches to promote the interpretation of results. Such studies could support the clinical implementation of genetic screening to provide safer drug use in pregnant women in need of drugs.

Rheumatoid Arthritis Patients, Both Newly Diagnosed and Methotrexate Treated, Show More DNA Methylation Differences in CD4+ Memory Than in CD4+ Naïve T Cells.

Background: Differences in DNA methylation have been reported in B and T lymphocyte populations, including CD4+ T cells, isolated from rheumatoid arthritis (RA) patients when compared to healthy controls. CD4+ T cells are a heterogeneous cell type with subpopulations displaying distinct DNA methylation patterns. In this study, we investigated DNA methylation using reduced representation bisulfite sequencing in two CD4+ T cell populations (CD4+ memory and naïve cells) in three groups: newly diagnosed, disease modifying antirheumatic drugs (DMARD) naïve RA patients (N = 11), methotrexate (MTX) treated RA patients (N = 18), and healthy controls (N = 9) matched for age, gender and smoking status. Results: Analyses of these data revealed significantly more differentially methylated positions (DMPs) in CD4+ memory than in CD4+ naïve T cells (904 vs. 19 DMPs) in RA patients compared to controls. The majority of DMPs (72%) identified in newly diagnosed and DMARD naïve RA patients with active disease showed increased DNA methylation (39 DMPs), whereas most DMPs (80%) identified in the MTX treated RA patients in remission displayed decreased DNA methylation (694 DMPs). Interestingly, we also found that about one third of the 101 known RA risk loci overlapped (±500 kb) with the DMPs. Notably, introns of the UBASH3A gene harbor both the lead RA risk SNP and two DMPs in CD4+ memory T cells. Conclusion: Our results suggest that RA associated DNA methylation differences vary between the two T cell subsets, but are also influenced by RA characteristics such as disease activity, disease duration and/or MTX treatment.

Systematic evaluation and validation of reference and library selection methods for deconvolution of cord blood DNA methylation data.

BACKGROUND: Umbilical cord blood (UCB) is commonly used in epigenome-wide association studies of prenatal exposures. Accounting for cell type composition is critical in such studies as it reduces confounding due to the cell specificity of DNA methylation (DNAm). In the absence of cell sorting information, statistical methods can be applied to deconvolve heterogeneous cell mixtures. Among these methods, reference-based approaches leverage age-appropriate cell-specific DNAm profiles to estimate cellular composition. In UCB, four reference datasets comprising DNAm signatures profiled in purified cell populations have been published using the Illumina 450 K and EPIC arrays. These datasets are biologically and technically different, and currently, there is no consensus on how to best apply them. Here, we systematically evaluate and compare these datasets and provide recommendations for reference-based UCB deconvolution. RESULTS: We first evaluated the four reference datasets to ascertain both the purity of the samples and the potential cell cross-contamination. We filtered samples and combined datasets to obtain a joint UCB reference. We selected deconvolution libraries using two different approaches: automatic selection using the top differentially methylated probes from the function pickCompProbes in minfi and a standardized library selected using the IDOL (Identifying Optimal Libraries) iterative algorithm. We compared the performance of each reference separately and in combination, using the two approaches for reference library selection, and validated the results in an independent cohort (Generation R Study, n = 191) with matched Fluorescence-Activated Cell Sorting measured cell counts. Strict filtering and combination of the references significantly improved the accuracy and efficiency of cell type estimates. Ultimately, the IDOL library outperformed the library from the automatic selection method implemented in pickCompProbes. CONCLUSION: These results have important implications for epigenetic studies in UCB as implementing this method will optimally reduce confounding due to cellular heterogeneity. This work provides guidelines for future reference-based UCB deconvolution and establishes a framework for combining reference datasets in other tissues.

Exploring the influence from whole blood DNA extraction methods on Infinium 450K DNA methylation.

Genome-wide DNA methylation studies are becoming increasingly important in unraveling the epigenetic basis of cell biology, aging and human conditions. The aim of the present study was to explore whether different methods for extracting DNA from whole blood can affect DNA methylation outcome, potentially confounding DNA methylation studies. DNA was isolated from healthy blood donors (n = 10) using three different extraction methods (i.e. two automatic extractions methods based on magnetic beads or isopropanol precipitation, and manual organic extraction). DNA methylation was analyzed using the Infinium HumanMethylation450 Bead Chip (Infinium 450K) (n = 30 samples in total), which is a frequently used method in genome-wide DNA methylation analyses. Overall, the different extraction methods did not have a significant impact on the global DNA methylation patterns. However, DNA methylation differences between organic extraction and each of the automated methods were in general larger than differences between the two automated extraction methods. No CpG sites or regions reached genome-wide significance when testing for differential methylation between extraction methods. Although this study is based on a small sample, these results suggest that extraction method is unlikely to confound Infinium 450K methylation analysis in whole blood.

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

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.

Differentially Methylated DNA Regions in Monozygotic Twin Pairs Discordant for Rheumatoid Arthritis: An Epigenome-Wide Study.

OBJECTIVES: In an explorative epigenome-wide association study (EWAS) to search for gene independent, differentially methylated DNA positions and regions (DMRs) associated with rheumatoid arthritis (RA) by studying monozygotic (MZ) twin pairs discordant for RA. METHODS: Genomic DNA was isolated from whole blood samples from 28 MZ twin pairs discordant for RA. DNA methylation was measured using the HumanMethylation450 BeadChips. Smoking, anti-cyclic citrullinated peptide antibodies, and immunosuppressive treatment were included as covariates. Pathway analysis was performed using GREAT. RESULTS: Smoking was significantly associated with hypomethylation of a DMR overlapping the promoter region of the RNF5 and the AGPAT1, which are implicated in inflammation and autoimmunity, whereas DMARD treatment induced hypermethylation of the same region. Additionally, the promotor region of both S100A6 and EFCAB4B were hypomethylated, and both genes have previously been associated with RA. We replicated several candidate genes identified in a previous EWAS in treatment-naïve RA singletons. Gene-set analysis indicated the involvement of immunologic signatures and cancer-related pathways in RA. CONCLUSION: We identified several differentially methylated regions associated with RA, which may represent environmental effects or consequences of the disease and plausible biological pathways pertinent to the pathogenesis of RA.