Differential DNA Methylation Research Articles

Evaluation of epigenetic methylation biomarkers for the detection of colorectal cancer using droplet digital PCR

Colorectal cancer (CRC) is the third most common cancer worldwide. Screening programs allow early diagnosis and have improved the clinical management of this disease. Aberrant DNA methylation is increasingly being explored as potential biomarkers for many types of cancers. In this study we investigate the methylation of ten target genes in 105 CRC and paired normal adjacent colonic tissue samples using a MethylLight droplet digital PCR (ML-ddPCR) assay. Receiver operator characteristic (ROC) curves were used to determine the diagnostic performance of all target genes individually and in combination. All 515 different combinations of genes showed significantly higher levels of methylation in CRC tissue. The combination of multiple target genes into a single test generally resulted in greater diagnostic accuracy when compared to single target genes. Our data confirms that ML-ddPCR is able to reliably detect significant differences in DNA methylation between CRC tissue and normal adjacent colonic tissue in a specific selection of target genes.

Sex‐specific DNA methylation differences associated with Alzheimer’s disease

AbstractBackgroundSex has increasingly been recognized as an important factor contributing to the heterogeneity in Alzheimer’s disease (AD).MethodWe performed a sex‐specific meta‐analysis of two large epigenome‐wide association studies (ADNI, AIBL) in the blood, with a total of 633 and 651 samples in females and males, respecitively. For each dataset, we adjusted covariate variables age, sex, batch, and immune cell‐type proportions. The Inverse‐variance fixed‐effects meta‐analysis model was then used to prioritize the most consistent DNAm differences across datasets.ResultIn female samples, 2 CpGs, mapped to the PRRC2A and RPS8 genes, reached the 5% false discovery rate (FDR). No CpGs reached 5% FDR in male samples analysis. At a more relaxed significance threshold of P < 1×10−5, an additional 22 CpGs and 5 CpGs were identified in female samples and male samples, respectively. For these 29 AD‐associated CpGs, the odds ratios ranged from 0.843 to 1.335 in females and 0.993 to 1.058 in males. The majority of these CpGs were hypermethylated in AD subjects (23 CpGs), located outside CpG islands or shores (27 CpGs), or in distal regions located greater than 2k bp from the TSS (27 CpGs). Only 2 of these 29 CpGs were located in gene promoters at SLC5A8 and C16orf89.At 5% Sidak corrected P‐value, our DMR analysis using comb‐p software identified 49 distinct differentially methylated regions (DMRs) in female and male samples each. The median numbers of CpGs in these DMRs are 5 CpGs in females and 4 CpGs in males. Among the top 10 most significant DMRs, about half of them (6 in females, 5 in males) were hypermethylated in AD. In females, 5 out of the 10 DMRs were mapped to promoter regions of the NNAT, CCDC169‐SOHLH2, ARHGEF15, LPAR5, and ZNF595 genes. In males, 4 out of the 10 DMRs were mapped to promoter regions of the MCCC1, PM20D1, PRR19, TTC23, and LRRC28 genes.ConclusionAs sex is a strong factor underlying phenotypic variability in AD, the results of our study are particularly relevant for a better understanding AD pathophysiology. They also provide a valuable resource for sex‐specific biomarkers and drug targets in AD.

DNA methylation differences in noncoding regions in ER negative breast tumors between Black and White women.

Incidence of estrogen receptor (ER)-negative breast cancer, an aggressive tumor subtype associated with worse prognosis, is higher among African American/Black women than other US racial and ethnic groups. The reasons for this disparity remain poorly understood but may be partially explained by differences in the epigenetic landscape. We previously conducted genome-wide DNA methylation profiling of ER- breast tumors from Black and White women and identified a large number of differentially methylated loci (DML) by race. Our initial analysis focused on DML mapping to protein-coding genes. In this study, motivated by increasing appreciation for the biological importance of the non-protein coding genome, we focused on 96 DMLs mapping to intergenic and noncoding RNA regions, using paired Illumina Infinium Human Methylation 450K array and RNA-seq data to assess the relationship between CpG methylation and RNA expression of genes located up to 1Mb away from the CpG site. Twenty-three (23) DMLs were significantly correlated with the expression of 36 genes (FDR<0.05), with some DMLs associated with the expression of single gene and others associated with more than one gene. One DML (cg20401567), hypermethylated in ER- tumors from Black versus White women, mapped to a putative enhancer/super-enhancer element located 1.3 Kb downstream of HOXB2. Increased methylation at this CpG correlated with decreased expression of HOXB2 (Rho=-0.74, FDR<0.001) and other HOXB/HOXB-AS genes. Analysis of an independent set of 207 ER- breast cancers from TCGA similarly confirmed hypermethylation at cg20401567 and reduced HOXB2 expression in tumors from Black versus White women (Rho=-0.75, FDR<0.001). Our findings indicate that epigenetic differences in ER- tumors between Black and White women are linked to altered gene expression and may hold functional significance in breast cancer pathogenesis.

DNA methylation and gene expression changes in mouse mammary tissue during successive lactations: part II – the impact of lactation rank

ABSTRACT Mastitis is among the main reasons women cease breastfeeding. In farm animals, mastitis results in significant economic losses and the premature culling of some animals. Nevertheless, the effect of inflammation on the mammary gland is not completely understood. This article discusses the changes to DNA methylation in mouse mammary tissue caused by lipopolysaccharide-induced inflammation after in vivo intramammary challenges and the differences in DNA methylation between 1st and 2nd lactations. Lactation rank induces 981 differential methylations of cytosines (DMCs) in mammary tissue. Inflammation in 1st lactation compared to inflammation in 2nd lactation results in the identification of 964 DMCs. When comparing inflammation in 1st vs. 2nd lactations with previous inflammation history, 2590 DMCs were identified. Moreover, Fluidigm PCR data show changes in the expression of several genes related to mammary function, epigenetic regulation, and the immune response. We show that the epigenetic regulation of two successive physiological lactations is not the same in terms of DNA methylation and that the effect of lactation rank on DNA methylation is stronger than that of the onset of inflammation. The conditions presented here show that few DMCs are shared between comparisons, suggesting a specific epigenetic response depending on lactation rank, the presence of inflammation, and even whether the cells had previously suffered inflammation. In the long term, this information could lead to a better understanding of the epigenetic regulation of lactation in both physiological and pathological conditions. Abbreviations: RRBS, reduced representation bisulphite sequencing; RT-qPCR, real-time quantitative polymerase chain reaction; MEC, mammary epithelial cells; MaSC, mammary stem cell; TSS, transcription start site; TTS, transcription termination site; UTR, untranslated region; SINE, short interspersed nuclear element; LINE, long interspersed nuclear element; CGI, CpG island; DEG, differentially expressed gene; DMC, differentially methylated cytosine; DMR, differentially methylated region; GO term, gene ontology term; MF, molecular function; BP, biological process

Metabolome plasticity in 241 Arabidopsis thaliana accessions reveals evolutionary cold adaptation processes.

Acclimation and adaptation of metabolism to a changing environment are key processes for plant survival and reproductive success. In the present study, 241 natural accessions of Arabidopsis (Arabidopsis thaliana) were grown under two different temperature regimes, 16 °C and 6 °C, and growth parameters were recorded, together with metabolite profiles, to investigate the natural genome × environment effects on metabolome variation. The plasticity of metabolism, which was captured by metabolic distance measures, varied considerably between accessions. Both relative growth rates and metabolic distances were predictable by the underlying natural genetic variation of accessions. Applying machine learning methods, climatic variables of the original growth habitats were tested for their predictive power of natural metabolic variation among accessions. We found specifically habitat temperature during the first quarter of the year to be the best predictor of the plasticity of primary metabolism, indicating habitat temperature as the causal driver of evolutionary cold adaptation processes. Analyses of epigenome- and genome-wide associations revealed accession-specific differential DNA-methylation levels as potentially linked to the metabolome and identified FUMARASE2 as strongly associated with cold adaptation in Arabidopsis accessions. These findings were supported by calculations of the biochemical Jacobian matrix based on variance and covariance of metabolomics data, which revealed that growth under low temperatures most substantially affects the accession-specific plasticity of fumarate and sugar metabolism. Our findings indicate that the plasticity of metabolic regulation is predictable from the genome and epigenome and driven evolutionarily by Arabidopsis growth habitats.

Traffic-related air pollution and supplemental folic acid intake in relation to DNA methylation in granulosa cells

BackgroundHigher exposure to traffic-related air pollution (TRAP) is related to lower fertility, with specific adverse effects on the ovary. Folic acid may attenuate these effects. Our goal was to explore the relation of TRAP exposure and supplemental folic acid intake with epigenetic aging and CpG-specific DNA methylation (DNAm) in granulosa cells (GC). Our study included 61 women undergoing ovarian stimulation at a fertility center (2005–2015). DNAm levels were profiled in GC using the Infinium MethylationEPIC BeadChip. TRAP was defined using a spatiotemporal model to estimate residence-based nitrogen dioxide (NO2) exposure. Supplemental folic acid intake was measured with a validated food frequency questionnaire. We used linear regression to evaluate whether NO2 or supplemental folic acid was associated with epigenetic age acceleration according to the Pan-tissue, mural GC, and GrimAge clocks or DNAm across the genome adjusting for potential confounders and accounting for multiple testing with a false discovery rate < 0.1.ResultsThere were no associations between NO2 or supplemental folic acid intake and epigenetic age acceleration of GC. NO2 and supplemental folic acid were associated with 9 and 11 differentially methylated CpG sites. Among these CpGs, only cg07287107 exhibited a significant interaction (p-value = 0.037). In women with low supplemental folic acid, high NO2 exposure was associated with 1.7% higher DNAm. There was no association between NO2 and DNAm in women with high supplemental folic acid. The genes annotated to the top 250 NO2-associated CpGs were enriched for carbohydrate and protein metabolism, postsynaptic potential and dendrite development, and membrane components and exocytosis. The genes annotated to the top 250 supplemental folic acid-associated CpGs were enriched for estrous cycle, learning, cognition, synaptic organization and transmission, and size and composition of neuronal cell bodies.ConclusionsWe found no associations between NO2, supplemental folic acid, and DNAm age acceleration of GC. However, there were 20 differentially methylated CpGs and multiple enriched GO terms associated with both exposures suggesting that differences in GC DNAm could be a plausible mechanism underlying the effects of TRAP and supplemental folic acid on ovarian function.

Saliva biopsy: Detecting the difference of EBV DNA methylation in the diagnosis of nasopharyngeal carcinoma.

Saliva sampling is a non-invasive method, and could be performed by donors themselves. However, there are few studies reporting biomarkers in saliva in the diagnosis of NPC. A total of 987 salivary samples were used in this study. First, EBV DNA methylation was profiled by capture sequencing in the discovery cohort (n = 36). Second, a q-PCR based method was developed and five representative EBV DNA CpG sites (11 029 bp, 45 849 bp, 57 945 bp, 66 226 bp and 128 102 bp) were selected and quantified to obtain the methylated density in the validation cohort1 (n = 801). Third, a validation cohort2 (n = 108) was used to further verify the differences of EBV methylation in saliva. A significant increase of EBV methylation was found in NPC patients compared with controls. The methylated score of EBV genome obtained by capture sequencing could distinguish NPC from controls (sensitivity 90%, specificity 100%). Further, the methylated density of EBV DNA CpG sites revealed by q-PCR showed a good diagnostic performance. The sensitivity and specificity of detecting a single CpG site (11 029 bp) could reach 75.4% and 99.7% in the validation cohort1, and 78.2% and 100% in the validation cohort2. Besides, the methylated density of the CpG site was found to decrease below the COV in NPC patients after therapy, and increase above the COV after recurrence. Our study provides an appealing alternative for the non-invasive detection of NPC without clinical setting. It paves the way for conducting a home-based large-scale screening in the future.

Multi-omic analysis in normal colon organoids highlights MSH4 as a novel marker of defective mismatch repair in Lynch syndrome and microsatellite instability.

Lynch syndrome (LS) is a hereditary condition that increases the risk of colorectal (CRC) and extracolonic cancers that exhibit microsatellite instability (MSI-H). MSI-H is driven by defective mismatch repair (dMMR), and approximately 15% of nonhereditary CRCs also exhibit MSI-H. Here, we aimed to better define mechanisms underlying tumor initiation in LS and MSI-H cancers through multi-omic analyses of LS normal colon organoids and MSI-H tumors. Right (n = 35) and left (n = 23) colon organoids generated from normal colon biopsies at routine colonoscopy of LS and healthy individuals were subjected to Illumina EPIC array. Differentially methylated region (DMR) analysis was performed by DMRcate. RNA-sequencing (n = 16) and bisulfite-sequencing (n = 15) were performed on a subset of right colon organoids. CRISPR-cas9-mediated editing of MMR genes in colon organoids of healthy individuals was followed by quantitative PCR of MSH4. The relationship between MSH4 expression and tumor mutational burden was further explored in three independent tumor data sets. We identified a hypermethylated region of MSH4 in both the right and left colon organoids of LS versus healthy controls, which we validated using bisulfite-sequencing. DMR analysis in three gastrointestinal and one endometrial data set revealed that this region was also hypermethylated in MSI-H versus microsatellite stable (MSS) tumors. MSH4 expression was increased in colon organoids of LS versus healthy subjects and in publicly available MSI-H versus MSS tumors across four RNA-seq and four microarray data sets. CRISPR-cas9 editing of MLH1 and MSH2, but not MSH6, in normal colon organoids significantly increased MSH4 expression. MSH4 expression was significantly associated with tumor mutational burden in three publicly available data sets. Our findings implicate DNA methylation and gene expression differences of MSH4 as a marker of dMMR and as a potential novel biomarker of LS. Our study of LS colon organoids supports the hypothesis that dMMR exists in the colons of LS subjects prior to CRC.

Identification of Cell-Specific Differential DNA Methylation Associated With Methotrexate Treatment Response in Rheumatoid Arthritis.

We undertook this study to estimate changes in cell-specific DNA methylation (DNAm) associated with methotrexate (MTX) response using whole blood samples collected from rheumatoid arthritis (RA) patients before and after initiation of MTX treatment. Patients included in this study were from the Rheumatoid Arthritis Medication Study (n=66) and the University of California San Francisco Rheumatoid Arthritis study (n=11). All patients met the American College of Rheumatology RA classification criteria. Blood samples were collected at baseline and following treatment. Disease Activity Scores in 28 joints using the C-reactive protein level were collected at baseline and after 3-6 months of treatment with MTX. Methylation profiles were generated using the Illumina Infinium HumanMethylation450 and MethylationEPIC v1.0 BeadChip arrays using DNA from whole blood. MTX response was defined using the EULAR response criteria (responders showed good/moderate response; nonresponders showed no response). Differentially methylated positions were identified using the Limma software package and Tensor Composition Analysis, which is a method for identifying cell-specific differential DNAm at the CpG level from tissue-level ("bulk") data. Differentially methylated regions were identified using Comb-p software. We found evidence of differential global methylation between treatment response groups. Further, we found patterns of cell-specific differential global methylation associated with MTX response. After correction for multiple testing, 1 differentially methylated position was associated with differential DNAm between responders and nonresponders at baseline in CD4+ T cells, CD8+ T cells, and natural killer cells. Thirty-nine cell-specific differentially methylated regions associated with MTX treatment response were identified. There were no significant findings in analyses of whole blood samples. We identified cell-specific changes in DNAm that were associated with MTX treatment response in RA patients. Future studies of DNAm and MTX treatment response should include measurements of DNAm from sorted cells.

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.

Differential effect of panic on the DNA methylation of the glucocorticoid receptor gene exon 1F in chronic subjective tinnitus with distress.

Tinnitus can be regarded as a chronic stressor, leading to dysregulation of the hypothalamic-pituitary-adrenal (HPA) axis. There is important comorbidity with anxiety, particularly panic, potentially associated with differences in HPA axis functioning and methylation patterns of HPA axis-related genes. This study examines DNA methylation of the glucocorticoid receptor gene ( NR3C1 ) exon 1F in adults with chronic subjective tinnitus and the possible differential effect of panic. In a well characterized tinnitus sample ( n = 22, half of which had co-occurring panic attacks), and unaffected controls ( n = 31) methylation patterns of the CpG sites were determined using pyrosequencing and compared between groups through linear mixed models. Gene expression was determined using quantitative PCR on mRNA. Comparing the combined tinnitus groups to the control group, no DNA methylation differences were observed; however, the tinnitus group with panic attacks showed consistently higher mean methylation values across all CpGs compared to the tinnitus-only and the control group ( P = 0.03 following Tukey correction), which became even more pronounced when accounting for childhood trauma ( P = 0.012). Moreover, a significant positive correlation was found between methylation of the CpG7 site and the Beck Anxiety Inventory total score ( P = 0.001) in the total population. NR3C1 -1F expression was not significantly different between the three groups. Panic is associated with higher DNA methylation of the NR3C1 exon 1F in adults with chronic subjective tinnitus, consistent with the reduced negative glucocorticoid feedback and HPA axis hyperfunction observed in individuals with panic disorder.

In utero exposure to diethylstilbestrol and blood DNA methylation in adult women: Results from a meta-analysis of two cohort studies

BackgroundPrenatal exposure to diethylstilbestrol (DES) is associated with several adverse health outcomes. Animal studies have shown associations between prenatal DES exposure and DNA methylation. ObjectiveThe aim of this study was to explore blood DNA methylation in women exposed and unexposed to DES in utero. MethodsSixty women (40 exposed and 20 unexposed) in the National Cancer Institute's Combined DES Cohort Study and 199 women (99 exposed and 100 unexposed women) in the Sister Study Cohort were included in this analysis. Within each study, robust linear regression models were used to assess associations between DES exposure and blood DNA methylation. Study-specific associations were combined using fixed-effect meta-analysis with inverse variance weights. Our analysis focused on CpG sites located within nine candidate genes identified in animal models. We further explored whether in utero DES exposure was associated with age acceleration. ResultsBlood DNA methylation levels at 10 CpG sites in six of the nine candidate genes were statistically significantly associated with prenatal DES exposure (P < 0.05) in this meta-analysis. Genes included EGF, EMB, EGFR, WNT11, FOS, and TGFB1, which are related to cell proliferation and differentiation. The most statistically significant CpG site was cg19830739 in gene EGF, and it was associated with lower methylation levels in women prenatally exposed to DES compared with those not exposed (P < 0.0001; false discovery rate<0.05). The association between prenatal DES exposure in utero and age acceleration was not statistically significant (P = 0.07 for meta-analyzed results). ConclusionsThere are few opportunities to investigate the effects of prenatal DES exposure. These findings suggest that in utero DES exposure may be associated with differential blood DNA methylation levels, which could mediate the increased risk of several adverse health outcomes observed in exposed women. Our findings need further evaluation using larger data sets.

Identification of a novel PARP4 gene promoter CpG locus associated with cisplatin chemoresistance

The protein family of poly (ADP-ribose) polymerases (PARPs) is comprised of multifunctional nuclear enzymes. Several PARP inhibitors have been developed as new anticancer drugs to combat resistance to chemotherapy. Herein, we characterized PARP4 mRNA expression profiles in cisplatin-sensitive and cisplatin-resistant ovarian cancer cell lines. PARP4 mRNA expression was significantly upregulated in cisplatin-resistant ovarian cancer cell lines, and this upregulation was associated with hypomethylation of specific cytosine-phosphate-guanine (CpG) sites (cg18582260 and cg17117459) on its promoter. Reduced PARP4 expression was restored by treating cisplatin-sensitive cell lines with a demethylation agent, implicating epigenetic regulation of PARP4 expression by promoter methylation. Depletion of PARP4 expression in cisplatin-resistant cell lines reduced cisplatin chemoresistance and promoted cisplatin-induced DNA fragmentation. The differential mRNA expression and DNA methylation status at specific PARP4 promoter CpG sites (cg18582260 and cg17117459) according to cisplatin responses, were further validated in primary ovarian tumor tissues. The results showed significantly increased PARP4 mRNA expressions and decreased DNA methylation levels at specific PARP4 promoter CpG sites (cg18582260 and cg17117459) in cisplatin-resistant patients. Additionally, the DNA methylation status at cg18582260 CpG sites in ovarian tumor tissues showed fairly clear discrimination between cisplatin-resistant patients and cisplatin-sensitive patients, with high accuracy (area under the curve = 0.86, p = 0.003845). Our findings suggested that the DNA methylation status of PARP4 at the specific promoter site (cg18582260) may be a useful diagnostic biomarker for predicting the response to cisplatin in ovarian cancer patients.

Validation of Myc-Associated Protein X (MAX) regulation in growth hormone secreting and nonfunctional pituitary adenoma.

Many patients with growth hormone-secreting pituitary adenoma (GHPA) fail to achieve biochemical remission, warranting investigation into epigenetic and molecular signatures associated with tumorigenesis and hormonal secretion. Prior work exploring the DNA methylome showed Myc-Associated Protein X (MAX), a transcription factor involved in cell cycle regulation, was differentially methylated between GHPA and nonfunctional pituitary adenoma (NFPA). We aimed to validate the differential DNA methylation and related MAX protein expression profiles between NFPA and GHPA. DNA methylation levels were measured in 52 surgically resected tumors (37 NFPA, 15 GHPA) at ~100,000 known MAX binding sites derived using ChIP-seq analysis from ENCODE. Findings were correlated with MAX protein expression using a constructed tissue microarray (TMA). Gene ontology analysis was performed to explore downstream genetic and signaling pathways regulated by MAX. GHPA had more hypomethylation events across all known MAX binding sites. Of binding sites defined using ChIP-seq analysis, 1,551 sites had significantly different methylation patterns between the two cohorts; 432 occurred near promoter regions potentially regulated by MAX, including promoters of TNF and MMP9. Gene ontology analysis suggested enrichment in genes involved in oxygen response, immune system regulation, and cell proliferation. Thirteen MAX binding sites were within coding regions of genes. GHPA demonstrated significantly increased expression of MAX protein compared to NFPA. GHPA have significantly different DNA methylation and downstream protein expression levels of MAX compared to NFPA. These differences may influence mechanisms involved with cellular proliferation, tumor invasion and hormonal secretion.

Child DNA methylation in a randomised controlled trial of a video-feedback intervention to promote positive parenting and sensitive discipline (VIPP-SD)

IntroductionA major modifiable risk factor for behavioural difficulties is harsh and insensitive parenting, and it has been hypothesised that the biological mechanism by which parenting influences child behaviour is via changes in the child's DNA methylation. We attempted to, in part, address the hypothesis that parenting is associated with child DNA methylation and, in turn, behaviour.MethodsPrimary caregivers of young children with behavioural difficulties (children aged 12–36 months) were randomised to receive a video-feedback Intervention to promote Positive Parenting and Sensitive Discipline (VIPP-SD) (n = 151), or usual care (n = 149). Child buccal samples were collected at a 2-year post-randomisation follow up (children aged 3–5 years, VIPP-SD group n = 106, usual care group n = 117) and were assessed for DNA methylation at the NR3C1, FKBP5 and OXYR genes. Child behaviour was assessed at baseline, post-intervention and 2-years post-randomisation using the Preschool Parental Account of Children's Symptoms (PPACS). We examined group differences in DNA methylation, associations of DNA methylation with behaviour, and sex differences.ResultsFor the NR3C1 and OXYR genes, there were no group differences, sex differences, or associations of DNA methylation with child behaviour, though all non-significant findings were in the hypothesised direction. For FKBP5 DNA methylation, there was a significant interaction between group and sex, such that males in the usual care group had higher DNA methylation than females, but in the intervention group females had higher DNA methylation than males. However, FKBP5 DNA methylation was not associated with behaviour in males or females.DiscussionWe provide the first evidence from a randomised controlled trial focused on improving parenting for sex-specific changes in child DNA methylation at a key gene involved in stress reactivity and psychopathology. This study adds to our understanding of causal mechanisms linking parenting with child behaviour, which is important for developing targeted interventions. A key limitation is that child DNA methylation was only assessed at one time point, so we were unable to assess change in DNA methylation over time. However, we demonstrate that is possible to collect and analyse DNA samples from families with young children receiving parenting interventions in the community, providing impetus for further research on this topic.

Single Cell Transcriptomics of Myeloid Cells in Black Americans with Multiple Sclerosis (P2-3.009)

<h3>Objective:</h3> Identify myeloid cell subsets from peripheral blood and monocyte-derived macrophages associated with MS in Black Americans. <h3>Background:</h3> Prior work in a cohort of patients composed primarily of underrepresented minority populations suggested an association between MS susceptibility and differential DNA methylation of genes that regulate cellular differentiation, proliferation, and invasion. <h3>Design/Methods:</h3> Peripheral blood mononuclear cells (PBMC) were obtained from Black American patients with MS and controls with non-inflammatory neurologic conditions under an IRB-approved protocol (UIC Neuroimmunology Biobank). Single cell libraries for RNAseq were constructed from PBMC (n=9) and monocyte-derived macrophages (MDM, n=6) using the Chromium Next GEM Single Cell 3′ Kit v3.1 and Chromium Next GEM Chip G (10x Genomics). Libraries were sequenced on an Illumina NovaSeq 6000. Expression data was analyzed in R software using the Seurat package (Butler et al., 2018). The SingleR and celldex packages (Aran et al., 2019) were used for identification of human peripheral immune cell types. <h3>Results:</h3> Peripheral blood monocytes from MS patients demonstrated an increase in relative numbers of transitional macrophages. These cells included those from classical CD163+ and non-classical CD16+ monocyte populations. Transcriptomic analysis showed an increase in the expression of a lncRNA (long non-coding RNA), KCNQ1OT1, in these cells. Analysis of MDM revealed expansion of a subset of C1Q1+ macrophages in MS. These cells express high levels of a pattern recognition receptor, FCN1 (ficolin 1). Global transcriptomic analysis revealed increased expression of SNHG5, a small nucleolar RNA host gene, in MS in all macrophage subtypes. <h3>Conclusions:</h3> These results suggest an association between expansion of specific monocyte-macrophage populations in our patient cohort. Increased expression of KCNQ1OT1 in transitional macrophages and SNHG5 in MDM is also consistent with our prior study of epigenetic biomarkers because these non-coding RNA’s are associated with cellular differentiation and proliferation in neoplastic diseases. <b>Disclosure:</b> Mr. Bingen has nothing to disclose. The institution of Dr. Clark has received research support from National Science Foundation. Ms. Wright has nothing to disclose. Fangxiu Xu has nothing to disclose. The institution of Dr. Carrithers has received research support from Biogen.

Intestinal epigenotype of Atlantic salmon (Salmo salar) associates with tenacibaculosis and gut microbiota composition

It remains a challenge to obtain the desired phenotypic traits in aquacultural production of Atlantic salmon, and part of the challenge might come from the effect that host-associated microorganisms have on the fish phenotype. To manipulate the microbiota towards the desired host traits, it is critical to understand the factors that shape it. The bacterial gut microbiota composition can vary greatly among fish, even when reared in the same closed system. While such microbiota differences can be linked to diseases, the molecular effect of disease on host-microbiota interactions and the potential involvement of epigenetic factors remain largely unknown.The aim of this study was to investigate the DNA methylation differences associated with a tenacibaculosis outbreak and microbiota displacement in the gut of Atlantic salmon. Using Whole Genome Bisulfite Sequencing (WGBS) of distal gut tissue from 20 salmon, we compared the genome-wide DNA methylation levels between uninfected individuals and sick fish suffering from tenacibaculosis and microbiota displacement. We discovered >19,000 differentially methylated cytosine sites, often located in differentially methylated regions, and aggregated around genes.The 68 genes connected to the most significant regions had functions related to the ulcerous disease such as epor and slc48a1a but also included prkcda and LOC106590732 whose orthologs are linked to microbiota changes in other species. Although the expression level was not analysed, our epigenetic analysis suggests specific genes potentially involved in host-microbiota interactions and more broadly it highlights the value of considering epigenetic factors in efforts to manipulate the microbiota of farmed fish.

A meta-analysis of epigenome-wide association studies on pregnancy vitamin B12 concentrations and offspring DNA methylation

ABSTRACT Circulating vitamin B12 concentrations during pregnancy are associated with offspring health. Foetal DNA methylation changes could underlie these associations. Within the Pregnancy And Childhood Epigenetics Consortium, we meta-analysed epigenome-wide associations of circulating vitamin B12 concentrations in mothers during pregnancy (n = 2,420) or cord blood (n = 1,029), with cord blood DNA methylation. Maternal and newborn vitamin B12 concentrations were associated with DNA methylation at 109 and 7 CpGs, respectively (False Discovery Rate P-value <0.05). Persistent associations with DNA methylation in the peripheral blood of up to 482 children aged 4–10 y were observed for 40.7% of CpGs associated with maternal vitamin B12 and 57.1% of CpGs associated with newborn vitamin B12. Of the CpGs identified in the maternal meta-analyses, 4.6% were associated with either birth weight or gestational age in a previous work. For the newborn meta-analysis, this was the case for 14.3% of the identified CpGs. Also, of the CpGs identified in the newborn meta-analysis, 14.3% and 28.6%, respectively, were associated with childhood cognitive skills and nonverbal IQ. Of the 109 CpGs associated with maternal vitamin B12, 18.3% were associated with nearby gene expression. In this study, we showed that maternal and newborn vitamin B12 concentrations are associated with DNA methylation at multiple CpGs in offspring blood (P FDR<0.05). Whether this differential DNA methylation underlies associations of vitamin B12 concentrations with child health outcomes, such as birth weight, gestational age, and childhood cognition, should be further examined in future studies.

Role of DNA methylation in the relationship between glioma risk factors and glioma incidence: a two-step Mendelian randomization study

Genetic evidence suggests glioma risk is altered by leukocyte telomere length, allergic disease (asthma, hay fever or eczema), alcohol consumption, childhood obesity, low-density lipoprotein cholesterol (LDLc) and triglyceride levels. DNA methylation (DNAm) variation influences many of these glioma-related traits and is an established feature of glioma. Yet the causal relationship between DNAm variation with both glioma incidence and glioma risk factors is unknown. We applied a two-step Mendelian randomization (MR) approach and several sensitivity analyses (including colocalization and Steiger filtering) to assess the association of DNAm with glioma risk factors and glioma incidence. We used data from a recently published catalogue of germline genetic variants robustly associated with DNAm variation in blood (32,851 participants) and data from a genome-wide association study of glioma risk (12,488 cases and 18,169 controls, sub-divided into 6191 glioblastoma cases and 6305 non-glioblastoma cases). MR evidence indicated that DNAm at 3 CpG sites (cg01561092, cg05926943, cg01584448) in one genomic region (HEATR3) had a putative association with glioma and glioblastoma risk (False discovery rate [FDR] < 0.05). Steiger filtering provided evidence against reverse causation. Colocalization presented evidence against genetic confounding and suggested that differential DNAm at the 3 CpG sites and glioma were driven by the same genetic variant. MR provided little evidence to suggest that DNAm acts as a mediator on the causal pathway between risk factors previously examined and glioma onset. To our knowledge, this is the first study to use MR to appraise the causal link of DNAm with glioma risk factors and glioma onset. Subsequent analyses are required to improve the robustness of our results and rule out horizontal pleiotropy.

The X-factor in ART: does the use of assisted reproductive technologies influence DNA methylation on the X chromosome?

BackgroundAssisted reproductive technologies (ART) may perturb DNA methylation (DNAm) in early embryonic development. Although a handful of epigenome-wide association studies of ART have been published, none have investigated CpGs on the X chromosome. To bridge this knowledge gap, we leveraged one of the largest collections of mother–father–newborn trios of ART and non-ART (natural) conceptions to date to investigate sex-specific DNAm differences on the X chromosome. The discovery cohort consisted of 982 ART and 963 non-ART trios from the Norwegian Mother, Father, and Child Cohort Study (MoBa). To verify our results from the MoBa cohort, we used an external cohort of 149 ART and 58 non-ART neonates from the Australian ‘Clinical review of the Health of adults conceived following Assisted Reproductive Technologies’ (CHART) study. The Illumina EPIC array was used to measure DNAm in both datasets. In the MoBa cohort, we performed a set of X-chromosome-wide association studies (‘XWASs’ hereafter) to search for sex-specific DNAm differences between ART and non-ART newborns. We tested several models to investigate the influence of various confounders, including parental DNAm. We also searched for differentially methylated regions (DMRs) and regions of co-methylation flanking the most significant CpGs. Additionally, we ran an analogous model to our main model on the external CHART dataset.ResultsIn the MoBa cohort, we found more differentially methylated CpGs and DMRs in girls than boys. Most of the associations persisted after controlling for parental DNAm and other confounders. Many of the significant CpGs and DMRs were in gene-promoter regions, and several of the genes linked to these CpGs are expressed in tissues relevant for both ART and sex (testis, placenta, and fallopian tube). We found no support for parental DNAm-dependent features as an explanation for the observed associations in the newborns. The most significant CpG in the boys-only analysis was in UBE2DNL, which is expressed in testes but with unknown function. The most significant CpGs in the girls-only analysis were in EIF2S3 and AMOT. These three loci also displayed differential DNAm in the CHART cohort.ConclusionsGenes that co-localized with the significant CpGs and DMRs associated with ART are implicated in several key biological processes (e.g., neurodevelopment) and disorders (e.g., intellectual disability and autism). These connections are particularly compelling in light of previous findings indicating that neurodevelopmental outcomes differ in ART-conceived children compared to those naturally conceived.