Detecting differentially methylated regions with multiple distinct associations

Cohort profile: the Genetics of Glucose regulation in Gestation and Growth (Gen3G) - a prospective prebirth cohort of mother-child pairs in Sherbrooke, Canada, 3-year and 5-year follow-up visits

PURPOSE

Initiated in 2010, the Genetics of Glucose regulation in Gestation and Growth (Gen3G) prospective cohort investigates the pathophysiology of impaired glycaemic regulation in pregnancy and evaluates its impact on both the mothers and her offspring health trajectory. Follow-up visits 3 and 5 years after delivery aimed to investigate pregnancy-related risk factors such as maternal obesity and gestational hyperglycaemia in relation to the mother's metabolic health after pregnancy, and with offspring health outcomes such as risk of obesity and neurodevelopmental problems in early childhood. We also investigated molecular mechanisms involved in the fetal programming of these later health outcomes.

PARTICIPANTS

Of the 1024 women originally recruited in the first trimester of pregnancy, we have targeted the 854 who had complete glucose tolerance test data and the 724 newborns who provided placenta and/or cord blood samples for follow-up recruitment. Of these, 695 mother-child dyads agreed to be contacted for the prospective follow-up visits. 448 and 521 mother-child dyads completed the research visits at 3 and 5 years after delivery respectively.

FINDINGS TO DATE

At both visits, we collected the mother's and child's medical history, lifestyle (using validated questionnaires), sociodemographic status, anthropometric measurements, mother's blood samples, child's saliva samples and growth charts. At the 5-year-old visit, we additionally collected the mother's and child's urine and stool samples and the child's blood samples; we performed a 75 g oral glucose tolerance test in the mothers and assessed the body composition in children using dual-energy X-ray absorptiometry. Using the Gen3G rich longitudinal data set, we have enhanced the understanding of the pathophysiology and characterisation of the heterogeneity of gestational diabetes mellitus, and we have shown that gestational hyperglycaemia and insulin resistance are associated with offspring epigenetics (DNA methylation) variations in the placenta, cord blood and blood at 5 years of age, as well as with offspring anthropometric, metabolic and neurodevelopmental outcomes in early childhood.

FUTURE PLANS

We are currently conducting a prospective follow-up of mothers and their children 12 years after delivery to study how prenatal and early-life metabolic factors may programme childhood adiposity and obesogenic dietary behaviours. This follow-up should be completed by the end of 2026.

Epigenome-wide association study of cerebrospinal fluid-based biomarkers of Alzheimer's disease in cognitively normal individuals

INTRODUCTION

Cerebrospinal fluid (CSF) biomarkers of Alzheimer's disease (AD) are reliable predictors of future AD risk. We investigated whether pre-clinical changes in AD CSF biomarkers are reflected in blood DNA methylation (DNAm) levels in cognitively normal participants.

METHODS

We profiled blood-based DNAm with the EPIC array in participants without a diagnosis of cognitive impairment in the Emory Healthy Brain Study (EHBS; N=495) and ADNI (N=122). Their CSF Aβ42, tTau, and pTau levels were quantified using Elecsys immunoassays. We conducted epigenome-wide association studies to assess associations between DNAm and CSF biomarkers of AD.

RESULTS

In EHBS, no loci were Bonferroni-significant after adjusting for confounding factors. In ADNI, two loci were significant, but they were not replicated in EHBS. There was little agreement between the top loci from EHBS and ADNI.

DISCUSSION

Our study showed little evidence of an association between differential blood-based DNAm and pre-clinical AD CSF biomarkers.

Parental arsenic exposure and tissue-specific DNA methylation in Bangladeshi infants with spina bifida

An emerging hypothesis linking arsenic toxicity involves altered epigenetic mechanisms, such as DNA methylation. In this study, we examined the relationship between parents' arsenic exposure and DNA methylation in tissues obtained from 28 infants with spina bifida from Bangladesh. We analyzed arsenic in parents' toenails using inductively coupled plasma mass spectrometry (ICP-MS). DNA methylation was measured in infants' dural tissue, buccal swabs, and whole blood using the Illumina Infinium MethylationEPIC BeadChip. We performed epigenome-wide association analyses (EWAS) and tested differentially methylated regions (DMRs). In EWAS, DNA methylation at cg24039697 in dural tissue was positively associated (β = 0.59, p = 7.6 × 10-9) with father's toenail arsenic concentrations, adjusting for covariates. We did not identify any CpG sites related to father's arsenic exposure in the other tissues, or any CpG sites related to mother's arsenic exposure. Gene ontology analysis identified many biological pathways of interest, including the Wnt signaling pathways. We identified several DMRs across the tissues related to arsenic exposure that included probes mapping to genes that have previously been identified in studies of neural tube defects. This study emphasizes the potential impact of arsenic exposure in fathers, often understudied in epidemiological studies, on DNA methylation in a unique neurological tissue specific to spina bifida.

A multi-trait epigenome-wide association study identified DNA methylation signature of inflammation among men with HIV

Inflammation underlies many conditions causing excess morbidity and mortality among people with HIV (PWH). A handful of single-trait epigenome-wide association studies (EWAS) have suggested that inflammation is associated with DNA methylation (DNAm) among PWH. Multi-trait EWAS may further improve statistical power and reveal pathways in common between different inflammatory markers. We conducted single-trait EWAS of three inflammatory markers (soluble CD14, D-dimers and interleukin-6) in the Veterans Aging Cohort Study (n = 920). The study population was all male PWH with an average age of 51 years, and 82.3% self-reported as Black. We then applied two multi-trait EWAS methods-CPASSOC and OmniTest-to combine single-trait EWAS results. CPASSOC and OmniTest identified 189 and 157 inflammation-associated DNAm sites, respectively, of which 112 overlapped. Among the identified sites, 56% were not significant in any single-trait EWAS. Top sites were mapped to inflammation-related genes including IFITM1, PARP9 and STAT1. These genes were significantly enriched in pathways such as "type I interferon signaling" and "immune response to virus." We demonstrate that multi-trait EWAS can improve the discovery of inflammation-associated DNAm sites, genes and pathways. These DNAm sites might hold the key to addressing persistent inflammation in PWH.

Gbdmr: identifying differentially methylated CpG regions in the human genome via generalized beta regressions

BACKGROUND

DNA methylation is a biochemical process in which a methyl group is added to the cytosine-phosphate-guanine (CpG) site on DNA molecules without altering the DNA sequence. Multiple CpG sites in a certain genome region can be differentially methylated across phenotypes. Identifying these differentially methylated CpG regions (DMRs) associated with the phenotypes contributes to disease prediction and precision medicine development.

RESULTS

We propose a novel DMR detection algorithm, gbdmr. In contrast to existing methods under a linear regression framework, gbdmr assumes that DNA methylation levels follow a generalized beta distribution. We compare gbdmr to alternative approaches via simulations and real data analyses, including dmrff, a new DMR detection approach that shows promising performance among competitors, and the traditional EWAS that focuses on single CpG sites. Our simulations demonstrate that gbdmr is superior to the other two when the correlation between neighboring CpG sites is strong, while dmrff shows a higher power when the correlation is weak. We provide an explanation of these phenomena from a theoretical perspective. We further applied the three methods to multiple real DNA methylation datasets. One is from a birth cohort study undertaken on the Isle of Wight, United Kingdom, and the other two are from the Gene Expression Omnibus database repository. Overall, gbdmr identifies more DMR CpGs linked to phenotypes than dmrff, and the simulated results support the findings.

CONCLUSIONS

Gbdmr is an innovative method for detecting DMRs based on generalized beta regression. It demonstrated notable advantages over dmrff and traditional EWAS, particularly when adjacent CpGs exhibited moderate to strong correlations. Our real data analyses and simulated findings highlight the reliability of gbdmr as a robust DMR detection tool. The gbdmr approach is accessible and implemented by R on GitHub: https://github.com/chengzhouwu/gbdmr .

Maternal adverse childhood experiences (ACEs) and DNA methylation of newborns in cord blood

BACKGROUND

Adverse childhood experiences (ACEs) increase the risk of poor health outcomes later in life. Psychosocial stressors may also have intergenerational health effects by which parental ACEs are associated with mental and physical health of children. Epigenetic programming may be one mechanism linking parental ACEs to child health. This study aimed to investigate epigenome-wide associations of maternal preconception ACEs with DNA methylation patterns of children. In the Center for the Health Assessment of Mothers and Children of Salinas study, cord blood DNA methylation was measured using the Illumina HumanMethylation450 BeadChip. Preconception ACEs, which occurred during the mothers' childhoods, were collected using a standard ACE questionnaire including 10 ACE indicators. Maternal ACE exposures were defined in this study as (1) the total number of ACEs; (2) the total number of ACEs categorized as 0, 1-3, and > 4; and (3) individual ACEs. Associations of ACE exposures with differential methylated positions, regions, and CpG modules determined using weighted gene co-expression network analysis were evaluated adjusting for covariates.

RESULTS

Data on maternal ACEs and cord blood DNA methylation were available for 196 mother/newborn pairs. One differential methylated position was associated with maternal experience of emotional abuse (cg05486260/FAM135B gene; q value < 0.05). Five differential methylated regions were significantly associated with the total number of ACEs, and 36 unique differential methylated regions were associated with individual ACEs (Šidák p value < 0.05). Fifteen CpG modules were significantly correlated with the total number of ACEs or individual ACEs, of which 8 remained significant in fully adjusted models (p value < 0.05). Significant modules were enriched for pathways related to neurological and immune development and function.

CONCLUSIONS

Maternal ACEs prior to conception were associated with cord blood DNA methylation of offspring at birth. Although there was limited overlap between differential methylated regions and CpGs in modules associated with ACE exposures, statistically significant regions and networks were related to genes involved in neurological and immune function. Findings may provide insights to pathways linking psychosocial stressors to health. Further research is needed to understand the relationship between changes in DNA methylation and child health.

Epigenome-wide association studies of allergic disease and the environment

The epigenome is at the intersection of the environment, genotype, and cellular response. DNA methylation of cytosine nucleotides, the most studied epigenetic modification, has been systematically evaluated in human studies by using untargeted epigenome-wide association studies (EWASs) and shown to be both sensitive to environmental exposures and associated with allergic diseases. In this narrative review, we summarize findings from key EWASs previously conducted on this topic; interpret results from recent studies; and discuss the strengths, challenges, and opportunities regarding epigenetics research on the environment-allergy relationship. The majority of these EWASs have systematically investigated select environmental exposures during the prenatal and early childhood periods and allergy-associated epigenetic changes in leukocyte-isolated DNA and more recently in nasal cells. Overall, many studies have found consistent DNA methylation associations across cohorts for certain exposures, such as smoking (eg, aryl hydrocarbon receptor repressor gene [AHRR] gene), and allergic diseases (eg, EPX gene). We recommend the integration of both environmental exposures and allergy or asthma within long-term prospective designs to strengthen causality as well as biomarker development. Future studies should collect paired target tissues to examine compartment-specific epigenetic responses, incorporate genetic influences in DNA methylation (methylation quantitative trait locus), replicate findings across diverse populations, and carefully interpret epigenetic signatures from bulk, target tissue or isolated cells.

Maternal glycemia in pregnancy is longitudinally associated with blood DNAm variation at the FSD1L gene from birth to 5 years of age

BACKGROUND

In utero exposure to maternal hyperglycemia has been associated with an increased risk for the development of chronic diseases in later life. These predispositions may be programmed by fetal DNA methylation (DNAm) changes that persist postnatally. However, although some studies have associated fetal exposure to gestational hyperglycemia with DNAm variations at birth, and metabolic phenotypes in childhood, no study has yet examined how maternal hyperglycemia during pregnancy may be associated with offspring DNAm from birth to five years of age.

HYPOTHESIS

Maternal hyperglycemia is associated with variation in offspring DNAm from birth to 5 years of age.

METHODS

We estimated maternal hyperglycemia using the area under the curve for glucose (AUCglu) following an oral glucose tolerance test conducted at 24-30 weeks of pregnancy. We quantified DNAm levels in cord blood (n = 440) and peripheral blood at five years of age (n = 293) using the Infinium MethylationEPIC BeadChip (Illumina). Our total sample included 539 unique dyads (mother-child) with 194 dyads having DNAm at both time-points. We first regressed DNAm M-values against the cell types and child age for each time-point separately to account for the difference by time of measurement for these variables. We then used a random intercept model from the linear mixed model (LMM) framework to assess the longitudinal association between maternal AUCglu and the repeated measures of residuals of DNAm. We adjusted for the following covariates as fixed effects in the random intercept model: maternal age, gravidity, smoking status, child sex, maternal body mass index (BMI) (measured at first trimester of pregnancy), and a binary variable for time-point.

RESULTS

In utero exposure to higher maternal AUCglu was associated with lower offspring blood DNAm levels at cg00967989 located in FSD1L gene (β = - 0.0267, P = 2.13 × 10-8) in adjusted linear regression mixed models. Our study also reports other CpG sites for which DNAm levels were suggestively associated (P < 1.0 × 10-5) with in utero exposure to gestational hyperglycemia. Two of these (cg12140144 and cg07946633) were found in the promotor region of PRDM16 gene (β: - 0.0251, P = 4.37 × 10-07 and β: - 0.0206, P = 2.24 × 10-06, respectively).

CONCLUSION

Maternal hyperglycemia is associated with offspring DNAm longitudinally assessed from birth to 5 years of age.

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

BACKGROUND

Assisted 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.

RESULTS

In 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.

CONCLUSIONS

Genes 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.

DNA methylation in cocaine use disorder-An epigenome-wide approach in the human prefrontal cortex

BACKGROUND

Cocaine use disorder (CUD) is characterized by a loss of control over cocaine intake and is associated with structural, functional, and molecular alterations in the human brain. At the molecular level, epigenetic alterations are hypothesized to contribute to the higher-level functional and structural brain changes observed in CUD. Most evidence of cocaine-associated epigenetic changes comes from animal studies while only a few studies have been performed using human tissue.

METHODS

We investigated epigenome-wide DNA methylation (DNAm) signatures of CUD in human post-mortem brain tissue of Brodmann area 9 (BA9). A total of N = 42 BA9 brain samples were obtained from N = 21 individuals with CUD and N = 21 individuals without a CUD diagnosis. We performed an epigenome-wide association study (EWAS) and analyzed CUD-associated differentially methylated regions (DMRs). To assess the functional role of CUD-associated differential methylation, we performed Gene Ontology (GO) enrichment analyses and characterized co-methylation networks using a weighted correlation network analysis. We further investigated epigenetic age in CUD using epigenetic clocks for the assessment of biological age.

RESULTS

While no cytosine-phosphate-guanine (CpG) site was associated with CUD at epigenome-wide significance in BA9, we detected a total of 20 CUD-associated DMRs. After annotation of DMRs to genes, we identified Neuropeptide FF Receptor 2 (NPFFR2) and Kalirin RhoGEF Kinase (KALRN) for which a previous role in the behavioral response to cocaine in rodents is known. Three of the four identified CUD-associated co-methylation modules were functionally related to neurotransmission and neuroplasticity. Protein-protein interaction (PPI) networks derived from module hub genes revealed several addiction-related genes as highly connected nodes such as Calcium Voltage-Gated Channel Subunit Alpha1 C (CACNA1C), Nuclear Receptor Subfamily 3 Group C Member 1 (NR3C1), and Jun Proto-Oncogene, AP-1 Transcription Factor Subunit (JUN). In BA9, we observed a trend toward epigenetic age acceleration (EAA) in individuals with CUD remaining stable even after adjustment for covariates.

CONCLUSION

Results from our study highlight that CUD is associated with epigenome-wide differences in DNAm levels in BA9 particularly related to synaptic signaling and neuroplasticity. This supports findings from previous studies that report on the strong impact of cocaine on neurocircuits in the human prefrontal cortex (PFC). Further studies are needed to follow up on the role of epigenetic alterations in CUD focusing on the integration of epigenetic signatures with transcriptomic and proteomic data.

EWAS of post-COVID-19 patients shows methylation differences in the immune-response associated gene, IFI44L, three months after COVID-19 infection

Although substantial progress has been made in managing COVID-19, it is still difficult to predict a patient’s prognosis. We explored the epigenetic signatures of COVID-19 in peripheral blood using data from an ongoing prospective observational study of COVID-19 called the Norwegian Corona Cohort Study. A series of EWASs were performed to compare the DNA methylation profiles between COVID-19 cases and controls three months post-infection. We also investigated differences associated with severity and long-COVID. Three CpGs—cg22399236, cg03607951, and cg09829636—were significantly hypomethylated (FDR < 0.05) in COVID-19 positive individuals. cg03607951 is located in IFI44L which is involved in innate response to viral infection and several systemic autoimmune diseases. cg09829636 is located in ANKRD9, a gene implicated in a wide variety of cellular processes, including the degradation of IMPDH2. The link between ANKRD9 and IMPDH2 is striking given that IMPDHs are considered therapeutic targets for COVID-19. Furthermore, gene ontology analyses revealed pathways involved in response to viruses. The lack of significant differences associated with severity and long-COVID may be real or reflect limitations in sample size. Our findings support the involvement of interferon responsive genes in the pathophysiology of COVID-19 and indicate a possible link to systemic autoimmune diseases.

CLINICOPATHOLOGIC AND MOLECULAR ANALYSIS OF A SINGLE BCOR-CCNB3+ UNDIFFERENTIATED SARCOMA OF THE KIDNEY CONFERS SIGNIFICANT EPIGENETIC ALTERATIONS

Undifferentiated soft tissue sarcomas (UDSTSs) are a group of mesenchymal tumors that remain a diagnostic challenge because of their morphologic heterogeneity and unclear histologic origin (Peters et al., Mod Pathol28: 575 [2015]). In this case report, we present the first multiomics molecular signature for a BCOR–CCNB3 sarcoma (BCS) that includes mutation analysis, gene expression, DNA methylation, and micro RNA (miRNA) expression. We identify a paucity of additional mutations in this tumor and detail that there is significant dysregulation of gene expression of epigenetic remodeling agents including key members of the PRC, Sin3A/3b, NuRD, and NcoR/SMRT complexes and the DNA methyltransferases DNMT1, DNMT3a, and DNMT3b. This is accompanied by significant DNA methylation changes and dysregulation of multiple miRNAs with known links to tumorigenesis. This study significantly increases our understanding of the BCOR effects on fusion-positive undifferentiated sarcomas at both the genomic and epigenomic level and suggests that as better-tailored and more refined treatment algorithms continue to evolve, epigenetic modifying agents should be further evaluated for their efficacy against these tumors.

Prenatal metal exposure, cord blood DNA methylation and persistence in childhood: an epigenome-wide association study of 12 metals

BACKGROUND

Prenatal exposure to essential and non-essential metals impacts birth and child health, including fetal growth and neurodevelopment. DNA methylation (DNAm) may be involved in pathways linking prenatal metal exposure and health. In the Project Viva cohort, we analyzed the extent to which metals (As, Ba, Cd, Cr, Cs, Cu, Hg, Mg, Mn, Pb, Se, and Zn) measured in maternal erythrocytes were associated with differentially methylated positions (DMPs) and regions (DMRs) in cord blood and tested if associations persisted in blood collected in mid-childhood. We measured metal concentrations in first-trimester maternal erythrocytes, and DNAm in cord blood (N = 361) and mid-childhood blood (N = 333, 6-10 years) with the Illumina HumanMethylation450 BeadChip. For each metal individually, we tested for DMPs using linear models (considered significant at FDR < 0.05), and for DMRs using comb-p (Sidak p < 0.05). Covariates included biologically relevant variables and estimated cell-type composition. We also performed sex-stratified analyses.

RESULTS

Pb was associated with decreased methylation of cg20608990 (CASP8) (FDR = 0.04), and Mn was associated with increased methylation of cg02042823 (A2BP1) in cord blood (FDR = 9.73 × 10-6). Both associations remained significant but attenuated in blood DNAm collected at mid-childhood (p < 0.01). Two and nine Mn-associated DMPs were identified in male and female infants, respectively (FDR < 0.05), with two and six persisting in mid-childhood (p < 0.05). All metals except Ba and Pb were associated with ≥ 1 DMR among all infants (Sidak p < 0.05). Overlapping DMRs annotated to genes in the human leukocyte antigen (HLA) region were identified for Cr, Cs, Cu, Hg, Mg, and Mn.

CONCLUSIONS

Prenatal metal exposure is associated with DNAm, including DMRs annotated to genes involved in neurodevelopment. Future research is needed to determine if DNAm partially explains the relationship between prenatal metal exposures and health outcomes.

Associations between infant sex and DNA methylation across umbilical cord blood, artery, and placenta samples

DNA methylation (DNAm) is vulnerable to dysregulation by environmental exposures during epigenetic reprogramming that occurs in embryogenesis. Sexual dimorphism in environmentally induced DNAm dysregulation has been identified and therefore it is important to understand sex-specific DNAm patterns. DNAm at several autosomal sites has been consistently associated with sex in cord blood and placental foetal tissues. However, there is limited research comparing sex-specific DNAm across tissues, particularly differentially methylated regions (DMRs). This study leverages DNAm data measured using the Illumina HumanMethylation450 BeadChip in cord blood (N = 179), placenta (N = 229), and umbilical artery samples (N = 229) in the PRogramming of Intergenerational Stress Mechanisms (PRISM) cohort to identify autosomal DMRs and differentially methylated positions (DMPs). A replication analyses was conducted in an independent cohort (GEO Accession GSE129841). We identified 183, 257, and 419 DMRs and 2119, 2281, and 3405 DMPs (p_Bonferroni < 0.05) in cord blood, placenta, and artery samples, respectively. Thirty-nine DMRs overlapped in all three tissues, overlapping with genes involved in spermatogenesis (NKAPL, PIWIL2 and AURKC) and X-inactivation (LRIF1). In replication analysis, 85% of DMRs overlapped with those identified in PRISM. Overall, DMRs and DMPs had higher methylation levels among females in cord blood and artery samples, but higher methylation levels among males in placenta samples. Further research is necessary to understand biological mechanisms that contribute to differences in sex-specific DNAm signatures across tissues, as well as to determine if sexual dimorphism in the epigenome impacts response to environmental stressors.