Recommendations for the design and analysis of epigenome-wide association studies

An epigenome-wide association study of waist circumference in Chinese monozygotic twins

OBJECTIVES

Central obesity poses significant health risks because it increases susceptibility to multiple chronic diseases. Epigenetic features such as DNA methylation may be associated with specific obesity traits, which could help us understand how genetic and environmental factors interact to influence the development of obesity. This study aims to identify DNA methylation sites associated with the waist circumference (WC) in Northern Han Chinese population, and to elucidate potential causal relationships.

METHODS

A total of 59 pairs of WC discordant monozygotic twins (ΔWC >0) were selected from the Qingdao Twin Registry in China. Generalized estimated equation model was employed to estimate the methylation levels of CpG sites on WC. Causal relationships between methylation and WC were assessed through the examination of family confounding factors using FAmiliaL CONfounding (ICE FALCON). Additionally, the findings of the epigenome-wide analysis were corroborated in the validation stage.

RESULTS

We identified 26 CpG sites with differential methylation reached false discovery rate (FDR) < 0.05 and 22 differentially methylated regions (slk-corrected p < 0.05) strongly linked to WC. These findings provided annotations for 26 genes, with notable emphasis on MMP17, ITGA11, COL23A1, TFPI, A2ML1-AS1, MRGPRE, C2orf82, and NINJ2. ICE FALCON analysis indicated the DNA methylation of ITGA11 and TFPI had a causal effect on WC and vice versa (p < 0.05). Subsequent validation analysis successfully replicated 10 (p < 0.05) out of the 26 identified sites.

CONCLUSIONS

Our research has ascertained an association between specific epigenetic variations and WC in the Northern Han Chinese population. These DNA methylation features can offer fresh insights into the epigenetic regulation of obesity and WC as well as hints to plausible biological mechanisms.

Negative association between higher maternal pre-pregnancy body mass index and breastfeeding outcomes is not mediated by DNA methylation

The benefits of breastfeeding for the health and wellbeing of both infants and mothers are well documented, yet global breastfeeding rates are low. One factor associated with low breast feeding is maternal body mass index (BMI), which is used as a measure of obesity. The negative relationship between maternal obesity and breastfeeding is likely caused by a variety of social, psychological, and physiological factors. Maternal obesity may also have a direct biological association with breastfeeding through changes in maternal DNA methylation. Here, we investigate this potential biological association using data from a UK-based cohort study, the Avon Longitudinal Study of Parents and Children (ALSPAC). We find that pre-pregnancy body mass index (BMI) is associated with lower initiation to breastfeed and shorter breastfeeding duration. We conduct epigenome-wide association studies (EWAS) of pre-pregnancy BMI and breastfeeding outcomes, and run candidate-gene analysis of methylation sites associated with BMI identified via previous meta-EWAS. We find that DNA methylation at cg11453712, annotated to PHTP1, is associated with pre-pregnancy BMI. From our results, neither this association nor those at candidate-gene sites are likely to mediate the link between pre-pregnancy BMI and breastfeeding.

Brain cell-type shifts in Alzheimer's disease, autism, and schizophrenia interrogated using methylomics and genetics

Few neuropsychiatric disorders have replicable biomarkers, prompting high-resolution and large-scale molecular studies. However, we still lack consensus on a more foundational question: whether quantitative shifts in cell types-the functional unit of life-contribute to neuropsychiatric disorders. Leveraging advances in human brain single-cell methylomics, we deconvolve seven major cell types using bulk DNA methylation profiling across 1270 postmortem brains, including from individuals diagnosed with Alzheimer's disease, schizophrenia, and autism. We observe and replicate cell-type compositional shifts for Alzheimer's disease (endothelial cell loss), autism (increased microglia), and schizophrenia (decreased oligodendrocytes), and find age- and sex-related changes. Multiple layers of evidence indicate that endothelial cell loss contributes to Alzheimer's disease, with comparable effect size to APOE genotype among older people. Genome-wide association identified five genetic loci related to cell-type composition, involving plausible genes for the neurovascular unit (P2RX5 and TRPV3) and excitatory neurons (DPY30 and MEMO1). These results implicate specific cell-type shifts in the pathophysiology of neuropsychiatric disorders.

Whole-genome methylation profiling reveals regions associated with painful temporomandibular disorders and active recovery processes

ABSTRACT

Temporomandibular disorders (TMDs), collectively representing one of the most common chronic pain conditions, have a substantial genetic component, but genetic variation alone has not fully explained the heritability of TMD risk. Reasoning that the unexplained heritability may be because of DNA methylation, an epigenetic phenomenon, we measured genome-wide DNA methylation using the Illumina MethylationEPIC platform with blood samples from participants in the Orofacial Pain: Prospective Evaluation and Risk Assessment (OPPERA) study. Associations with chronic TMD used methylation data from 496 chronic painful TMD cases and 452 TMD-free controls. Changes in methylation between enrollment and a 6-month follow-up visit were determined for a separate sample of 62 people with recent-onset painful TMD. More than 750,000 individual CpG sites were examined for association with chronic painful TMD. Six differentially methylated regions were significantly ( P < 5 × 10 -8 ) associated with chronic painful TMD, including loci near genes involved in the regulation of inflammatory and neuronal response. A majority of loci were similarly differentially methylated in acute TMD consistent with observed transience or persistence of symptoms at follow-up. Functional characterization of the identified regions found relationships between methylation at these loci and nearby genetic variation contributing to chronic painful TMD and with gene expression of proximal genes. These findings reveal epigenetic contributions to chronic painful TMD through methylation of the genes FMOD , PM20D1 , ZNF718 , ZFP57 , and RNF39 , following the development of acute painful TMD. Epigenetic regulation of these genes likely contributes to the trajectory of transcriptional events in affected tissues leading to resolution or chronicity of pain.

Neurodevelopmental Programming of Adiposity: Contributions to Obesity Risk

This review analyzes the published evidence regarding maternal factors that influence the developmental programming of long-term adiposity in humans and animals via the central nervous system (CNS). We describe the physiological outcomes of perinatal underfeeding and overfeeding and explore potential mechanisms that may mediate the impact of such exposures on the development of feeding circuits within the CNS-including the influences of metabolic hormones and epigenetic changes. The perinatal environment, reflective of maternal nutritional status, contributes to the programming of offspring adiposity. The in utero and early postnatal periods represent critically sensitive developmental windows during which the hormonal and metabolic milieu affects the maturation of the hypothalamus. Maternal hyperglycemia is associated with increased transfer of glucose to the fetus driving fetal hyperinsulinemia. Elevated fetal insulin causes increased adiposity and consequently higher fetal circulating leptin concentration. Mechanistic studies in animal models indicate important roles of leptin and insulin in central and peripheral programming of adiposity, and suggest that optimal concentrations of these hormones are critical during early life. Additionally, the environmental milieu during development may be conveyed to progeny through epigenetic marks and these can potentially be vertically transmitted to subsequent generations. Thus, nutritional and metabolic/endocrine signals during perinatal development can have lifelong (and possibly multigenerational) impacts on offspring body weight regulation.

Mechanisms of Epigenetic Inheritance in Post-Traumatic Stress Disorder

Post-traumatic stress disorder (PTSD) is a psychiatric disorder that causes debilitating functional impairment in patients. Observations from survivors of traumatic historical events solidify that this disease is not only associated with personal experiences but can also be inherited from familial traumas. Over the past decades, researchers have focused on epigenetic inheritance to understand how responses to adverse experiences can be passed down to future generations. This review aims to present recent findings on epigenetic markers related to PTSD and research in the intergenerational inheritance of trauma. By understanding the information, we hope that epigenetic markers can act as biochemical measurements for future clinical practice.

Trans-ancestry epigenome-wide association meta-analysis of DNA methylation with lifetime cannabis use

Cannabis is widely used worldwide, yet its links to health outcomes are not fully understood. DNA methylation can serve as a mediator to link environmental exposures to health outcomes. We conducted an epigenome-wide association study (EWAS) of peripheral blood-based DNA methylation and lifetime cannabis use (ever vs. never) in a meta-analysis including 9436 participants (7795 European and 1641 African ancestry) from seven cohorts. Accounting for effects of cigarette smoking, our trans-ancestry EWAS meta-analysis revealed four CpG sites significantly associated with lifetime cannabis use at a false discovery rate of 0.05 ( p < 5.85 × 10 - 7 ) : cg22572071 near gene ADGRF1, cg15280358 in ADAM12, cg00813162 in ACTN1, and cg01101459 near LINC01132. Additionally, our EWAS analysis in participants who never smoked cigarettes identified another epigenome-wide significant CpG site, cg14237301 annotated to APOBR. We used a leave-one-out approach to evaluate methylation scores constructed as a weighted sum of the significant CpGs. The best model can explain 3.79% of the variance in lifetime cannabis use. These findings unravel the DNA methylation changes associated with lifetime cannabis use that are independent of cigarette smoking and may serve as a starting point for further research on the mechanisms through which cannabis exposure impacts health outcomes.

Metastable epialleles in humans

First identified in isogenic mice, metastable epialleles (MEs) are loci where the extent of DNA methylation (DNAm) is variable between individuals but correlates across tissues derived from different germ layers within a given individual. This property, termed systemic interindividual variation (SIV), is attributed to stochastic methylation establishment before germ layer differentiation. Evidence suggests that some putative human MEs are sensitive to environmental exposures in early development. In this review we introduce key concepts pertaining to human MEs, describe methods used to identify MEs in humans, and review their genomic features. We also highlight studies linking DNAm at putative human MEs to early environmental exposures and postnatal (including disease) phenotypes.

A systematic review of neuroimaging epigenetic research: calling for an increased focus on development

Epigenetic mechanisms, such as DNA methylation (DNAm), have gained increasing attention as potential biomarkers and mechanisms underlying risk for neurodevelopmental, psychiatric and other brain-based disorders. Yet, surprisingly little is known about the extent to which DNAm is linked to individual differences in the brain itself, and how these associations may unfold across development - a time of life when many of these disorders emerge. Here, we systematically review evidence from the nascent field of Neuroimaging Epigenetics, combining structural or functional neuroimaging measures with DNAm, and the extent to which the developmental period (birth to adolescence) is represented in these studies. We identified 111 articles published between 2011-2021, out of which only a minority (21%) included samples under 18 years of age. Most studies were cross-sectional (85%), employed a candidate-gene approach (67%), and examined DNAm-brain associations in the context of health and behavioral outcomes (75%). Nearly half incorporated genetic data, and a fourth investigated environmental influences. Overall, studies support a link between peripheral DNAm and brain imaging measures, but there is little consistency in specific findings and it remains unclear whether DNAm markers present a cause, correlate or consequence of brain alterations. Overall, there is large heterogeneity in sample characteristics, peripheral tissue and brain outcome examined as well as the methods used. Sample sizes were generally low to moderate (median nall = 98, ndevelopmental = 80), and attempts at replication or meta-analysis were rare. Based on the strengths and weaknesses of existing studies, we propose three recommendations on how advance the field of Neuroimaging Epigenetics. We advocate for: (1) a greater focus on developmentally oriented research (i.e. pre-birth to adolescence); (2) the analysis of large, prospective, pediatric cohorts with repeated measures of DNAm and imaging to assess directionality; and (3) collaborative, interdisciplinary science to identify robust signals, triangulate findings and enhance translational potential.

DNA methylation analysis is used to identify novel genetic loci associated with circulating fibrinogen levels in blood

BACKGROUND

Fibrinogen plays an essential role in blood coagulation and inflammation. Circulating fibrinogen levels may be determined based on interindividual differences in DNA methylation at cytosine-phosphate-guanine (CpG) sites and vice versa.

OBJECTIVES

To perform an EWAS to examine an association between blood DNA methylation levels and circulating fibrinogen levels to better understand its biological and pathophysiological actions.

METHODS

We performed an epigenome-wide association study of circulating fibrinogen levels in 18 037 White, Black, American Indian, and Hispanic participants, representing 14 studies from the Cohorts for Heart and Aging Research in Genomic Epidemiology consortium. Circulating leukocyte DNA methylation was measured using the Illumina 450K array in 12 904 participants and using the EPIC array in 5133 participants. In each study, an epigenome-wide association study of fibrinogen was performed using linear mixed models adjusted for potential confounders. Study-specific results were combined using array-specific meta-analysis, followed by cross-replication of epigenome-wide significant associations. We compared models with and without CRP adjustment to examine the role of inflammation.

RESULTS

We identified 208 and 87 significant CpG sites associated with fibrinogen levels from the 450K (p < 1.03 × 10-7) and EPIC arrays (p < 5.78 × 10-8), respectively. There were 78 associations from the 450K array that replicated in the EPIC array and 26 vice versa. After accounting for overlapping sites, there were 83 replicated CpG sites located in 61 loci, of which only 4 have been previously reported for fibrinogen. The examples of genes located near these CpG sites were SOCS3 and AIM2, which are involved in inflammatory pathways. The associations of all 83 replicated CpG sites were attenuated after CRP adjustment, although many remained significant.

CONCLUSION

We identified 83 CpG sites associated with circulating fibrinogen levels. These associations are partially driven by inflammatory pathways shared by both fibrinogen and CRP.

Examining the immunoepigenetic-gut microbiome axis in the context of self-esteem among Native Hawaiians and other Pacific Islanders

Introduction: Native Hawaiian and other Pacific Islander (NHPI) populations experience higher rates of immunometabolic diseases compared to other racial-ethnic groups in Hawaii. As annual NHPI mortality rates for suicide and type 2 diabetes mellitus (T2DM) exceed those of the state as a whole, understanding the social and biological mechanisms underlying these disparities are urgently needed to enable preventive strategies. Methods: A community-based approach was used to investigate the immunoepigenetic-gut microbiome axis in an NHPI-enriched cohort of Oahu residents (N = 68). Self-esteem (SE) data was collected using a modified Rosenberg self-esteem (SE) assessment as a proxy measure for mental wellbeing in consideration for cultural competency. T2DM status was evaluated using point-of-care A1c (%) tests. Stool samples were collected for 16s-based metagenomic sequencing analyses. Plasma from blood samples were isolated by density-gradient centrifugation. Peripheral blood mononuclear cells (PBMCs) were collected from the same samples and enriched for monocytes using negative selection techniques. Flow-cytometry was used for immunoprofiling assays. Monocyte DNA was extracted for Illumina EPIC array-based methylation analysis. Results: Compared to individuals with normal SE (NSE), those with low SE (LSE) exhibited significantly higher plasma concentrations (pg/ml) of proinflammatory cytokines IL-8 (p = 0.051) and TNF-α (p = 0.011). Metagenomic analysis revealed that the relative abundance (%) of specific gut bacteria significantly differed between SE groups - some of which directly correlated with SE scores. Gene ontology analysis revealed that 104 significantly differentially methylated loci (DML) between SE groups were preferentially located at genes involved in immunometabolic processes. Horvath clock analyses indicated epigenetic age (Epi-Age) deceleration in individuals with LSE and acceleration in individuals with NSE (p = 0.042), yet was not reproduced by other clocks. Discussion: These data reveal novel differences in the immunoepigenetic-gut microbiome axis with respect to SE, warranting further investigation into its relationship to brain activity and mental health in NHPI. Unexpected results from Epi-Age analyses warrant further investigation into the relationship between biological age and disparate health outcomes among the NHPI population. The modifiable component of epigenetic processes and the gut microbiome makes this axis an attractive target for potential therapeutics, biomarker discovery, and novel prevention strategies.

Calculating detection limits and uncertainty of reference-based deconvolution of whole-blood DNA methylation data

DNA methylation (DNAm)-based cell mixture deconvolution (CMD) has become a quintessential part of epigenome-wide association studies where DNAm is profiled in heterogeneous tissue types. Despite being introduced over a decade ago, detection limits, which represent the smallest fraction of a cell type in a mixed biospecimen that can be reliably detected, have yet to be determined in the context of DNAm-based CMD. Moreover, there has been little attention given to approaches for quantifying the uncertainty associated with DNAm-based CMD. Here, analytical frameworks for determining both cell-specific limits of detection and quantification of uncertainty associated with DNAm-based CMD are described. This work may contribute to improved rigor, reproducibility and replicability of epigenome-wide association studies involving CMD.

Brain-derived neurotrophic factor (BDNF) epigenomic modifications and brain-related phenotypes in humans: A systematic review

Epigenomic modifications of the brain-derived neurotrophic factor (BDNF) gene have been postulated to underlie the pathogenesis of neurodevelopmental, psychiatric, and neurological conditions. This systematic review summarizes current evidence investigating the association of BDNF epigenomic modifications (DNA methylation, non-coding RNA, histone modifications) with brain-related phenotypes in humans. A novel contribution is our creation of an open access web-based application, the BDNF DNA Methylation Map, to interactively visualize specific positions of CpG sites investigated across all studies for which relevant data were available. Our literature search of four databases through September 27, 2021 returned 1701 articles, of which 153 met inclusion criteria. Our review revealed exceptional heterogeneity in methodological approaches, hindering the identification of clear patterns of robust and/or replicated results. We summarize key findings and provide recommendations for future epigenomic research. The existing literature appears to remain in its infancy and requires additional rigorous research to fulfill its potential to explain BDNF-linked risk for brain-related conditions and improve our understanding of the molecular mechanisms underlying their pathogenesis.

Distinct genome-wide DNA methylation and gene expression signatures in classical monocytes from African American patients with systemic sclerosis

BACKGROUND

Systemic sclerosis (SSc) is a multisystem autoimmune disorder that has an unclear etiology and disproportionately affects women and African Americans. Despite this, African Americans are dramatically underrepresented in SSc research. Additionally, monocytes show heightened activation in SSc and in African Americans relative to European Americans. In this study, we sought to investigate DNA methylation and gene expression patterns in classical monocytes in a health disparity population.

METHODS

Classical monocytes (CD14+ + CD16-) were FACS-isolated from 34 self-reported African American women. Samples from 12 SSc patients and 12 healthy controls were hybridized on MethylationEPIC BeadChip array, while RNA-seq was performed on 16 SSc patients and 18 healthy controls. Analyses were computed to identify differentially methylated CpGs (DMCs), differentially expressed genes (DEGs), and CpGs associated with changes in gene expression (eQTM analysis).

RESULTS

We observed modest DNA methylation and gene expression differences between cases and controls. The genes harboring the top DMCs, the top DEGs, as well as the top eQTM loci were enriched for metabolic processes. Genes involved in immune processes and pathways showed a weak upregulation in the transcriptomic analysis. While many genes were newly identified, several other have been previously reported as differentially methylated or expressed in different blood cells from patients with SSc, supporting for their potential dysregulation in SSc.

CONCLUSIONS

While contrasting with results found in other blood cell types in largely European-descent groups, the results of this study support that variation in DNA methylation and gene expression exists among different cell types and individuals of different genetic, clinical, social, and environmental backgrounds. This finding supports the importance of including diverse, well-characterized patients to understand the different roles of DNA methylation and gene expression variability in the dysregulation of classical monocytes in diverse populations, which might help explaining the health disparities.

Systematic and benchmarking studies of pipelines for mammal WGBS data in the novel NGS platform

BACKGROUND

Whole genome bisulfite sequencing (WGBS), possesses the aptitude to dissect methylation status at the nucleotide-level resolution of 5-methylcytosine (5-mC) on a genome-wide scale. It is a powerful technique for epigenome in various cell types, and tissues. As a recently established next-generation sequencing (NGS) platform, GenoLab M is a promising alternative platform. However, its comprehensive evaluation for WGBS has not been reported. We sequenced two bisulfite-converted mammal DNA in this research using our GenoLab M and NovaSeq 6000, respectively. Then, we systematically compared those data via four widely used WGBS tools (BSMAP, Bismark, BatMeth2, BS-Seeker2) and a new bisulfite-seq tool (BSBolt). We interrogated their computational time, genome depth and coverage, and evaluated their percentage of methylated Cs.

RESULT

Here, benchmarking a combination of pre- and post-processing methods, we found that trimming improved the performance of mapping efficiency in eight datasets. The data from two platforms uncovered ~ 80% of CpG sites genome-wide in the human cell line. Those data sequenced by GenoLab M achieved a far lower proportion of duplicates (~ 5.5%). Among pipelines, BSMAP provided an intriguing representation of 5-mC distribution at CpG sites with 5-mC levels > ~ 78% in datasets from human cell lines, especially in the GenoLab M. BSMAP performed more advantages in running time, uniquely mapped reads percentages, genomic coverage, and quantitative accuracy. Finally, compared with the previous methylation pattern of human cell line and mouse tissue, we confirmed that the data from GenoLab M performed similar consistency and accuracy in methylation levels of CpG sites with that from NovaSeq 6000.

CONCLUSION

Together we confirmed that GenoLab M was a qualified NGS platform for WGBS with high performance. Our results showed that BSMAP was the suitable pipeline that allowed for WGBS studies on the GenoLab M platform.

Using epigenomics to understand cellular responses to environmental influences in diseases

It is a generally accepted model that environmental influences can exert their effects, at least in part, by changing the molecular regulators of transcription that are described as epigenetic. As there is biochemical evidence that some epigenetic regulators of transcription can maintain their states long term and through cell division, an epigenetic model encompasses the idea of maintenance of the effect of an exposure long after it is no longer present. The evidence supporting this model is mostly from the observation of alterations of molecular regulators of transcription following exposures. With the understanding that the interpretation of these associations is more complex than originally recognised, this model may be oversimplistic; therefore, adopting novel perspectives and experimental approaches when examining how environmental exposures are linked to phenotypes may prove worthwhile. In this review, we have chosen to use the example of nonalcoholic fatty liver disease (NAFLD), a common, complex human disease with strong environmental and genetic influences. We describe how epigenomic approaches combined with emerging functional genetic and single-cell genomic techniques are poised to generate new insights into the pathogenesis of environmentally influenced human disease phenotypes exemplified by NAFLD.

Methylome-wide association study of different responses to risperidone in schizophrenia

Background: Accumulating evidence shows that DNA methylation plays a role in antipsychotic response. However, the mechanisms by which DNA methylation changes are associated with antipsychotic responses remain largely unknown. Methods: We performed a methylome-wide association study (MWAS) to evaluate the association between DNA methylation and the response to risperidone in schizophrenia. Genomic DNA methylation patterns were assessed using the Agilent Human DNA Methylation Microarray. Results: We identified numerous differentially methylated positions (DMPs) and regions (DMRs) associated with antipsychotic response. CYP46A1, SPATS2, and ATP6V1E1 had the most significant DMPs, with p values of 2.50 × 10^-6, 3.53 × 10^-6, and 5.71 × 10^-6, respectively. The top-ranked DMR was located on chromosome 7, corresponding to the PTPRN2 gene with a Šidák-corrected p-value of 9.04 × 10^-13. Additionally, a significant enrichment of synaptic function and neurotransmitters was found in the differentially methylated genes after gene ontology and pathway analysis. Conclusion: The identified DMP- and DMR-overlapping genes associated with antipsychotic response are related to synaptic function and neurotransmitters. These findings may improve understanding of the mechanisms underlying antipsychotic response and guide the choice of antipsychotic in schizophrenia.

Childhood Trauma and Epigenetics: State of the Science and Future

PURPOSE OF REVIEW

There is a great deal of interest regarding the biological embedding of childhood trauma and social exposures through epigenetic mechanisms, including DNA methylation (DNAm), but a comprehensive understanding has been hindered by issues of limited reproducibility between studies. This review presents a summary of the literature on childhood trauma and DNAm, highlights issues in the field, and proposes some potential solutions.

RECENT FINDINGS

Investigations of the associations between DNAm and childhood trauma are commonly performed using candidate gene approaches, specifically involving genes related to neurological and stress pathways. Childhood trauma is defined in a wide range of ways in several societal contexts. However, although variations in DNAm are frequently found in stress-related genes, unsupervised epigenome-wide association studies (EWAS) have shown limited reproducibility both between studies and in relating these changes to exposures. The reproducibility of childhood trauma DNAm studies, and the field of social epigenetics in general, may be improved by increasing sample sizes, standardizing variables, making use of effect size thresholds, collecting longitudinal and intervention samples, appropriately accounting for known confounding factors, and applying causal analysis wherever possible, such as "two-step epigenetic Mendelian randomization."

Epigenetic impact of a 1-week intensive multimodal group program for adolescents with multiple adverse childhood experiences

Adverse childhood experiences (ACEs, i.e., abuse, neglect, household dysfunction) represent a potential risk factor for a wide range of long-lasting diseases and shorter life expectancy. We recently described a 1-week residential group program, based on mindfulness training, artistic expression and EMDR group therapy, that significantly reduced PTSD-related symptoms and increased attention/awareness-related outcomes in adolescent girls with multiple ACEs in a randomized controlled study. Since epigenetic mechanisms (i.e., DNA methylation) have been associated with the long-lasting effects of ACEs, the present report extends these prior findings by exploring genome-wide DNA methylation changes following the program. Saliva samples from all participants (n = 44) were collected and genomic DNA was extracted prior (T1) and following (T2) the intervention. Genome-wide DNA methylation analysis using the MethylationEPIC beadchip array (Illumina) revealed 49 differentially methylated loci (DML; p value < 0.001; methylation change > 10%) that were annotated to genes with roles in biological processes linked to early childhood adversity (i.e., neural, immune, and endocrine pathways, cancer and cardiovascular disease). DNA sequences flanking these DML showed significant enrichment of transcription factor binding sites involved in inflammation, cancer, cardiovascular disease, and brain development. Methylation changes in SIRT5 and TRAPPC2L genes showed associations with changes in trauma-related psychological measures. Results presented here suggest that this multimodal group program for adolescents with multiple victimization modulates the DNA methylome at sites of potential relevance for health and behavioral disorders associated with ACEs.

No evidence for intervention-associated DNA methylation changes in monocytes of patients with posttraumatic stress disorder

DNA methylation patterns can be responsive to environmental influences. This observation has sparked interest in the potential for psychological interventions to influence epigenetic processes. Recent studies have observed correlations between DNA methylation changes and therapy outcome. However, most did not control for changes in cell composition. This study had two aims: first, we sought to replicate therapy-associated changes in DNA methylation of commonly assessed candidate genes in isolated monocytes from 60 female patients with post-traumatic stress disorder (PTSD). Our second, exploratory goal was to identify novel genomic regions with substantial pre-to-post intervention DNA methylation changes by performing whole-genome bisulfite sequencing (WGBS) in two patients with PTSD. Equivalence testing and Bayesian analyses provided evidence against physiologically meaningful intervention-associated DNA methylation changes in monocytes of PTSD patients in commonly investigated target genes (NR3C1, FKBP5, SLC6A4, OXTR). Furthermore, WGBS yielded only a limited set of candidate regions with suggestive evidence of differential DNA methylation pre- to post-therapy. These differential DNA methylation patterns did not prove replicable when investigated in the entire cohort. We conclude that there is no evidence for major, recurrent intervention-associated DNA methylation changes in the investigated genes in monocytes of patients with PTSD.

Social Factors, Epigenomics and Lupus in African American Women (SELA) Study: protocol for an observational mechanistic study examining the interplay of multiple individual and social factors on lupus outcomes in a health disparity population

Introduction

Despite the disproportional impact of SLE on historically marginalised communities, the individual and sociocultural factors underlying these health disparities remain elusive. We report the design and methods for a study aimed at identifying epigenetic biomarkers associated with racism and resiliency that affect gene function and thereby influence SLE in a health disparity population.

Methods and analysis

The Social Factors, Epigenomics and Lupus in African American Women (SELA) Study is a cross-sectional, case–control study. A total of 600 self-reported African American women will be invited to participate. All participants will respond to questionnaires that capture detailed sociodemographic and medical history, validated measures of racial discrimination, social support, as well as disease activity and damage for cases. Participants who wish will receive their genetic ancestry estimates and be involved in research. Blood samples are required to provide peripheral blood mononuclear cell counts, DNA and RNA. The primary goals of SELA are to identify variation in DNA methylation (DNAm) associated with self-reported exposure to racial discrimination and social support, to evaluate whether social DNAm sites affect gene expression, to identify the synergistic effects of social factors on DNAm changes on SLE and to develop a social factors-DNAm predictive model for disease outcomes. This study is conducted in cooperation with the Sea Island Families Project Citizen Advisory Committee.

Discussion and dissemination

SELA will respond to the pressing need to clarify the interplay and regulatory mechanism by which various positive and negative social exposures influence SLE. Results will be published and shared with patients and the community. Knowledge of the biological impact of social exposures on SLE, as informed by the results of this study, can be leveraged by advocacy efforts to develop psychosocial interventions that prevent or mitigate risk exposures, and services or interventions that promote positive exposures. Implementation of such interventions is paramount to the closure of the health disparities gap.

LARGE-SCALE INTEGRATION OF DNA METHYLATION AND GENE EXPRESSION ARRAY PLATFORMS IDENTIFIES BOTH cis AND trans RELATIONSHIPS

Although epigenome-wide association studies (EWAS) have been successful in identifying DNA methylation (DNAm) patterns associated with disease states, any further characterization of etiologic mechanisms underlying disease remains elusive. This knowledge gap does not originate from a lack of DNAm–trait associations, but rather stems from study design issues that affect the interpretability of EWAS results. Despite known limitations in predicting the function of a particular CpG site, most EWAS maintain the broad assumption that altered DNAm results in a concomitant change of transcription at the most proximal gene. This study integrated DNAm and gene expression (GE) measurements in two cohorts, the Adolescent and Young Adult Twin Study (AYATS) and the Pregnancy, Race, Environment, Genes (PREG) study, to improve the understanding of epigenomic regulatory mechanisms. CpG sites associated with GE in cis were enriched in areas of transcription factor binding and areas of intermediate-to-low CpG density. CpG sites associated with trans GE were also enriched in areas of known regulatory significance, including enhancer regions. These results highlight issues with restricting DNAm-transcript annotations to small genomic intervals and question the validity of assuming a cis DNAm–GE pathway. Based on these findings, the interpretation of EWAS results is limited in studies without multi-omic support and further research should identify genomic regions in which GE-associated DNAm is overrepresented. An in-depth characterization of GE-associated CpG sites could improve predictions of the downstream functional impact of altered DNAm and inform best practices for interpreting DNAm–trait associations generated by EWAS.

Network Approaches for Charting the Transcriptomic and Epigenetic Landscape of the Developmental Origins of Health and Disease

The early developmental phase is of critical importance for human health and disease later in life. To decipher the molecular mechanisms at play, current biomedical research is increasingly relying on large quantities of diverse omics data. The integration and interpretation of the different datasets pose a critical challenge towards the holistic understanding of the complex biological processes that are involved in early development. In this review, we outline the major transcriptomic and epigenetic processes and the respective datasets that are most relevant for studying the periconceptional period. We cover both basic data processing and analysis steps, as well as more advanced data integration methods. A particular focus is given to network-based methods. Finally, we review the medical applications of such integrative analyses.

Association of Glyphosate Exposure with Blood DNA Methylation in a Cross-Sectional Study of Postmenopausal Women

BACKGROUND

Glyphosate is the most commonly used herbicide in the world and is purported to have a variety of health effects, including endocrine disruption and an elevated risk of several types of cancer. Blood DNA methylation has been shown to be associated with many other environmental exposures, but to our knowledge, no studies to date have examined the association between blood DNA methylation and glyphosate exposure.

OBJECTIVE

We conducted an epigenome-wide association study to identify DNA methylation loci associated with urinary glyphosate and its metabolite aminomethylphosphonic acid (AMPA) levels. Secondary goals were to determine the association of epigenetic age acceleration with glyphosate and AMPA and develop blood DNA methylation indices to predict urinary glyphosate and AMPA levels.

METHODS

For 392 postmenopausal women, white blood cell DNA methylation was measured using the Illumina Infinium MethylationEPIC BeadChip array. Glyphosate and AMPA were measured in two urine samples per participant using liquid chromatography-tandem mass spectrometry. Methylation differences at the probe and regional level associated with glyphosate and AMPA levels were assessed using a resampling-based approach. Probes and regions that had an false discovery rate q < 0.1 in ≥ 90 % of 1,000 subsamples of the study population were considered differentially methylated. Differentially methylated sites from the probe-specific analysis were combined into a methylation index. Epigenetic age acceleration from three epigenetic clocks and an epigenetic measure of pace of aging were examined for associations with glyphosate and AMPA.

RESULTS

We identified 24 CpG sites whose methylation level was associated with urinary glyphosate concentration and two associated with AMPA. Four regions, within the promoters of the MSH4, KCNA6, ABAT, and NDUFAF2/ERCC8 genes, were associated with glyphosate levels, along with an association between ESR1 promoter hypomethylation and AMPA. The methylation index accurately predicted glyphosate levels in an internal validation cohort. AMPA, but not glyphosate, was associated with greater epigenetic age acceleration.

DISCUSSION

Glyphosate and AMPA exposure were associated with DNA methylation differences that could promote the development of cancer and other diseases. Further studies are warranted to replicate our results, determine the functional impact of glyphosate- and AMPA-associated differential DNA methylation, and further explore whether DNA methylation could serve as a biomarker of glyphosate exposure. https://doi.org/10.1289/EHP10174.

A Novel Framework for the Identification of Reference DNA Methylation Libraries for Reference-Based Deconvolution of Cellular Mixtures

Reference-based deconvolution methods use reference libraries of cell-specific DNA methylation (DNAm) measurements as a means toward deconvoluting cell proportions in heterogeneous biospecimens (e.g., whole-blood). As the accuracy of such methods depends highly on the CpG loci comprising the reference library, recent research efforts have focused on the selection of libraries to optimize deconvolution accuracy. While existing approaches for library selection work extremely well, the best performing approaches require a training data set consisting of both DNAm profiles over a heterogeneous cell population and gold-standard measurements of cell composition (e.g., flow cytometry) in the same samples. Here, we present a framework for reference library selection without a training dataset (RESET) and benchmark it against the Legacy method (minfi:pickCompProbes), where libraries are constructed based on a pre-specified number of cell-specific differentially methylated loci (DML). RESET uses a modified version of the Dispersion Separability Criteria (DSC) for comparing different libraries and has four main steps: 1) identify a candidate set of cell-specific DMLs, 2) randomly sample DMLs from the candidate set, 3) compute the Modified DSC of the selected DMLs, and 4) update the selection probabilities of DMLs based on their contribution to the Modified DSC. Steps 2–4 are repeated many times and the library with the largest Modified DSC is selected for subsequent reference-based deconvolution. We evaluated RESET using several publicly available datasets consisting of whole-blood DNAm measurements with corresponding measurements of cell composition. We computed the RMSE and R² between the predicted cell proportions and their measured values. RESET outperformed the Legacy approach in selecting libraries that improve the accuracy of deconvolution estimates. Additionally, reference libraries constructed using RESET resulted in cellular composition estimates that explained more variation in DNAm as compared to the Legacy approach when evaluated in the context of epigenome-wide association studies (EWAS) of several publicly available data sets. This finding has implications for the statistical power of EWAS. RESET combats potential challenges associated with existing approaches for reference library assembly and thus, may serve as a viable strategy for library construction in the absence of a training data set.

Risks and rewards of big-data in epigenomics research: an interview with Melanie Ehrlich

Melanie Ehrlich, PhD, is a professor in the Tulane Cancer Center, the Tulane Center for Medical Bioinformatics and Genomics and the Hayward Human Genetics Program at Tulane Medical School, New Orleans, LA. She obtained her PhD in molecular biology in 1971 from the State University of New York at Stony Brook and completed postdoctoral research at Albert Einstein College of Medicine in 1972. She has been working on various aspects of epigenetics, starting with DNA methylation, since 1973. Her group made many first findings about DNA methylation (see below). For example, in 1982 and 1983, in collaboration with Charles Gehrke at the University of Missouri, she was the first to report tissue-specific and cancer-specific differences in overall DNA methylation in humans. In 1985, Xian-Yang Zhang and Richard Wang in her lab discovered a class of human DNA sequences specifically hypomethylated in sperm. In 1998, her group was the first to describe extensive losses of DNA methylation in pericentromeric and centromeric DNA repeats in human cancer. Her lab's many publications on the prevalence of both DNA hypermethylation and hypomethylation in the same cancers brought needed balance to our understanding of the epigenetics of cancer and to its clinical implications [1]. Besides working on cancer epigenetics, her research group has helped elucidate cytogenetic and gene expression abnormalities in the immunodeficiency, centromeric and facial anomalies (ICF) syndrome, a rare recessive disease often caused by mutations in DNMT3B. Her group also studied the epigenetics and transcriptomics of facioscapulohumeral muscular dystrophy (FSHD), whose disease locus is a tandem 3.3-kb repeat at subtelomeric 4q (that happens to be hypomethylated in ICF DNA [2]). Her study of FSHD has taken her in the direction of muscle (skeletal muscle, heart and aorta) epigenetics [3-6]. Recently, she has led research that applies epigenetics much more rigorously than usual to the evaluation of genetic variants from genome-wide association studies (GWAS) of osteoporosis and obesity. In continued collaboration with Sriharsa Pradhan at New England Biolabs and Michelle Lacey at Tulane University, she has compared 5-hydroxymethylcytosine and 5-methylcytosine clustering in various human tissues [7] and is studying myoblast methylomes that they generated by a new high-resolution enzymatic technique (enzymatic methyl-seq).

Periconception and Prenatal Exposure to Maternal Perceived Stress and Cord Blood DNA Methylation

BACKGROUND

Maternal prenatal stress is associated with physiologic and adverse mental health outcomes in the offspring, but the underlying biologic mechanisms are unknown. We examined the associations of maternal perceived stress, including preconception exposure, with DNA methylation (DNAm) alterations in the cord blood buffy coats of 358 singleton infants.

METHODS

Maternal perceived stress was measured prior to and throughout pregnancy in a cohort of women enrolled in Effects of Aspirin in Gestation and Reproduction Trial (EAGeR) trial. Perceived stress assessments based on a standardized Likert-scale were obtained in periconception (~2 months preconception and 2-8 weeks of gestation) and pregnancy (8-36 weeks of gestation). Cumulative perceived stress was estimated by calculating the predicted area under the curve of stress reported prior to and during pregnancy. DNAm was measured by the Infinium MethylationEPIC BeadChip. Multivariable robust linear regression was used to assess associations of perceived stress with individual CpG probes.

RESULTS

Based on a 0 to 3 scale, average reported preconception and early pregnancy stress were 0.76 (0.60) and 0.67 (0.50), respectively. Average mid- to late-pregnancy stress, based on a 0 to 10 scale, was 4.9 (1.6). Neither periconception nor pregnancy perceived stress were associated with individual CpG sites in neonatal cord blood (all false discovery rate [FDR] >5%).

CONCLUSION

No effects of maternal perceived stress exposure on array-wide cord blood neonatal methylation differences were found.

"The study of expression levels of DNA methylation regulators in patients affected with congenital heart defects (CHDs)"

BACKGROUND

Congenial heart defects (CHDs) have multifactorial etiology with complex interplay of genetic and environmental factors. Environmental impact can have epigenetic mechanism of CHD development. Many studies have reported the causal association between CHD and distinct DNA methylation profile which is one of the key epigenetic events, which has vital role in normal embryonic development. The products of DNMT1, DNMT3A, DNMT3B, and MBD2 are important regulators of DNA methylation process. Changes in the expression of these genes are implicated in congenital structural cardiac defects. Hence, in this proof-of-concept study, we have compared the expression levels of these genes in the blood samples of healthy controls and CHD cases while investigating the etiology of CHD.

METHODS

In this study with 48 CHD cases and 47 healthy controls, total RNA was isolated from the whole blood samples using TRI reagent. Quantitative RT PCR (qRT-PCR) was used to analyze the mRNA levels of DNMT1, DNMT3A, DNMT3B, and MBD2. The expression levels have been analyzed by relative quantification.

RESULTS

We observed that DNMT3B (fold change = -2.563; p = .0018) and DNMT3A (fold change = -2.169; p = .05) were significantly downregulated in CHD patients, whereas the expression of DNMT1 and MBD2 was not significantly different between cases and controls.

CONCLUSIONS

Lower expression of de novo methyltransferases, namely, DNMT3B and DNMT3A in CHD cases, may be an important contributor to the mechanism of CHD pathogenesis. Further studies with age-matched controls and analysis of global DNA methylation profile are required to investigate the proposed causal association.

A history of childhood maltreatment is associated with altered DNA methylation levels of DNA methyltransferase 1 in maternal but not neonatal mononuclear immune cells

Childhood maltreatment (CM) is associated with alterations in DNA methylation (DNAm) especially in stress response genes. Due to the higher risk of overall health complications of individuals with a parental history of CM, intergenerational transmission of CM-associated DNAm changes has been investigated but remains unclear. In this study, we investigated if different severities of CM have any influence on the DNAm of DNA methyltransferase 1 (DNMT1), an important enzyme of the DNAm machinery, in immune and buccal cells of mother-newborn dyads. DNAm was assessed by mass spectrometry using immune cell DNA from mothers (N = 117) and their newborns (N = 113), and buccal cell DNA of mother-newborn dyads (N = 68 each). Mothers with a history of CM had lower mean methylation of DNMT1 in immune cells compared to the mothers without a CM history. CM status only influenced maternal DNMT1 gene expression when at least moderate CM was reported. Buccal cell DNAm was not associated with CM status. Maternal history of CM was not linked to any alterations in DNMT1 mean DNAm in any of the cell types studied in newborns. We conclude that the CM-associated alterations in DNMT1 DNAm might point to allostatic load and can be physiologically relevant, especially in individuals with more severe CM experiences, resulting in an activated DNA methylation machinery that might influence stress response genes. Our lack of significant findings in buccal cells shows the tissue-specific effects of CM on DNAm. In our sample with low to moderate maternal CM history, there was no intergenerational transmission of DNMT1 DNAm in newborns.

Maternal haemoglobin levels in pregnancy and child DNA methylation: a study in the pregnancy and childhood epigenetics consortium

Altered maternal haemoglobin levels during pregnancy are associated with pre-clinical and clinical conditions affecting the fetus. Evidence from animal models suggests that these associations may be partially explained by differential DNA methylation in the newborn with possible long-term consequences. To test this in humans, we meta-analyzed the epigenome-wide associations of maternal haemoglobin levels during pregnancy with offspring DNA methylation in 3,967 newborn cord blood and 1,534 children and 1,962 adolescent whole-blood samples derived from 10 cohorts. DNA methylation was measured using Illumina Infinium Methylation 450K or MethylationEPIC arrays covering 450,000 and 850,000 methylation sites, respectively. There was no statistical support for the association of maternal haemoglobin levels with offspring DNA methylation either at individual methylation sites or clustered in regions. For most participants, maternal haemoglobin levels were within the normal range in the current study, whereas adverse perinatal outcomes often arise at the extremes. Thus, this study does not rule out the possibility that associations with offspring DNA methylation might be seen in studies with more extreme maternal haemoglobin levels.

DNA Methylation Profiles of Purified Cell Types in Bronchoalveolar Lavage: Applications for Mixed Cell Paediatric Pulmonary Studies

In epigenome-wide association studies analysing DNA methylation from samples containing multiple cell types, it is essential to adjust the analysis for cell type composition. One well established strategy for achieving this is reference-based cell type deconvolution, which relies on knowledge of the DNA methylation profiles of purified constituent cell types. These are then used to estimate the cell type proportions of each sample, which can then be incorporated to adjust the association analysis. Bronchoalveolar lavage is commonly used to sample the lung in clinical practice and contains a mixture of different cell types that can vary in proportion across samples, affecting the overall methylation profile. A current barrier to the use of bronchoalveolar lavage in DNA methylation-based research is the lack of reference DNA methylation profiles for each of the constituent cell types, thus making reference-based cell composition estimation difficult. Herein, we use bronchoalveolar lavage samples collected from children with cystic fibrosis to define DNA methylation profiles for the four most common and clinically relevant cell types: alveolar macrophages, granulocytes, lymphocytes and alveolar epithelial cells. We then demonstrate the use of these methylation profiles in conjunction with an established reference-based methylation deconvolution method to estimate the cell type composition of two different tissue types; a publicly available dataset derived from artificial blood-based cell mixtures and further bronchoalveolar lavage samples. The reference DNA methylation profiles developed in this work can be used for future reference-based cell type composition estimation of bronchoalveolar lavage. This will facilitate the use of this tissue in studies examining the role of DNA methylation in lung health and disease.

Comparing the epigenetic landscape in myonuclei purified with a PCM1 antibody from a fast/glycolytic and a slow/oxidative muscle

Muscle cells have different phenotypes adapted to different usage, and can be grossly divided into fast/glycolytic and slow/oxidative types. While most muscles contain a mixture of such fiber types, we aimed at providing a genome-wide analysis of the epigenetic landscape by ChIP-Seq in two muscle extremes, the fast/glycolytic extensor digitorum longus (EDL) and slow/oxidative soleus muscles. Muscle is a heterogeneous tissue where up to 60% of the nuclei can be of a different origin. Since cellular homogeneity is critical in epigenome-wide association studies we developed a new method for purifying skeletal muscle nuclei from whole tissue, based on the nuclear envelope protein Pericentriolar material 1 (PCM1) being a specific marker for myonuclei. Using antibody labelling and a magnetic-assisted sorting approach, we were able to sort out myonuclei with 95% purity in muscles from mice, rats and humans. The sorting eliminated influence from the other cell types in the tissue and improved the myo-specific signal. A genome-wide comparison of the epigenetic landscape in EDL and soleus reflected the differences in the functional properties of the two muscles, and revealed distinct regulatory programs involving distal enhancers, including a glycolytic super-enhancer in the EDL. The two muscles were also regulated by different sets of transcription factors; e.g. in soleus, binding sites for MEF2C, NFATC2 and PPARA were enriched, while in EDL MYOD1 and SIX1 binding sites were found to be overrepresented. In addition, more novel transcription factors for muscle regulation such as members of the MAF family, ZFX and ZBTB14 were identified.

Complex tissues like skeletal muscle contain a variety of cells which confound the analysis of each cell type when based on homogenates, thus only about half of the cell nuclei in muscles reside inside the muscle cells. We here describe a labelling and sorting technique that allowed us to study the epigenetic landscape in purified muscle cell nuclei leaving the other cell types out. Differences between a fast/glycolytic and a slow/oxidative muscle were studied. While all skeletal muscle fibers have a similar make up and basic function, they differ in their physiology and the way they are used. Thus, some fibers are fast contracting but fatigable, and are used for short lasting explosive tasks such as sprinting. Other fibers are slow and are used for more prolonged tasks such as standing or long distance running. Since fiber type correlate with metabolic profile these features can also be related to metabolic diseases. We here show that the epigenetic landscape differed in gene loci corresponding to the differences in functional properties, and revealed that the two types are enriched in different gene regulatory networks. Exercise can alter muscle phenotype, and the epigenetic landscape might be related to how plastic different properties are.

Integrating genetic and social factors to understand health disparities in lupus

PURPOSE OF REVIEW

Both social and genetic factors are associated with health outcomes in systemic lupus erythematosus (SLE), thus playing a role in its health disparities. Despite the growing list of social and genetic factors associated with SLE outcomes, studies integrating sociocultural and individual determinants of health to understand health disparities in SLE are lacking. We review the contributions of different social and genetic factors to the disparities in SLE, and propose a socioecological model to integrate and examine the complex interactions between individual and social factors in SLE outcomes.

RECENT FINDINGS

Multiple studies collecting comprehensive social data and biospecimens from diverse populations are underway, which will contribute to the elucidation of the interplay and underlying mechanisms by which positive and negative social determinants of health influence epigenomic variation, and how the resulting biological changes may contribute to the lupus health disparities.

SUMMARY

There is growing awareness of the need to integrate genomic and health disparities research to understand how social exposures affect disease outcomes. Understanding the contributions of these factors to the SLE health disparity will inform the development of interventions to eliminate risk exposures and close the health disparity gap.

Host methylation predicts SARS-CoV-2 infection and clinical outcome

BACKGROUND

Since the onset of the SARS-CoV-2 pandemic, most clinical testing has focused on RT-PCR1. Host epigenome manipulation post coronavirus infection2-4 suggests that DNA methylation signatures may differentiate patients with SARS-CoV-2 infection from uninfected individuals, and help predict COVID-19 disease severity, even at initial presentation.

METHODS

We customized Illumina's Infinium MethylationEPIC array to enhance immune response detection and profiled peripheral blood samples from 164 COVID-19 patients with longitudinal measurements of disease severity and 296 patient controls.

RESULTS

Epigenome-wide association analysis revealed 13,033 genome-wide significant methylation sites for case-vs-control status. Genes and pathways involved in interferon signaling and viral response were significantly enriched among differentially methylated sites. We observe highly significant associations at genes previously reported in genetic association studies (e.g. IRF7, OAS1). Using machine learning techniques, models built using sparse regression yielded highly predictive findings: cross-validated best fit AUC was 93.6% for case-vs-control status, and 79.1%, 80.8%, and 84.4% for hospitalization, ICU admission, and progression to death, respectively.

CONCLUSIONS

In summary, the strong COVID-19-specific epigenetic signature in peripheral blood driven by key immune-related pathways related to infection status, disease severity, and clinical deterioration provides insights useful for diagnosis and prognosis of patients with viral infections.

Maternal caffeine intake and DNA methylation in newborn cord blood

BACKGROUND

Epigenetic mechanisms may underlie associations between maternal caffeine consumption and adverse childhood metabolic outcomes. However, limited studies have examined neonate DNA methylation (DNAm) patterns in the context of preconception or prenatal exposure to caffeine metabolites.

OBJECTIVES

We examined preconception and pregnancy caffeine exposure with DNAm alterations in neonate cord blood (n = 378).

METHODS

In a secondary analysis of the Effects of Aspirin in Gestation and Reproduction Trial (EAGeR), we measured maternal caffeine, paraxanthine, and theobromine concentrations from stored serum collected preconception (on average 2 months before pregnancy) and at 8 weeks of gestation. In parallel, self-reported caffeinated beverage intake was captured via administration of questionnaires and daily diaries. We profiled DNAm from the cord blood buffy coat of singletons using the MethylationEPIC BeadChip. We assessed associations of maternal caffeine exposure and methylation β values using multivariable robust linear regression. A false discovery rate (FDR) correction was applied using the Benjamini-Hochberg method.

RESULTS

In preconception, the majority of women reported consuming 1 or fewer servings/day of caffeine on average, and caffeine and paraxanthine metabolite levels were 88 and 36 µmol/L, respectively. Preconception serum caffeine metabolites were not associated with individual cytosine-guanine (CpG) sites (FDR >5%), though pregnancy theobromine was associated with DNAm at cg09460369 near RAB2A (β = 0.028; SE = 0.005; FDR P = 0.012). Preconception self-reported caffeinated beverage intake compared to no intake was associated with DNAm at cg09002832 near GLIS3 (β = -0.013; SE = 0.002; FDR P = 0.036). No associations with self-reported intake during pregnancy were found.

CONCLUSIONS

Few effects of maternal caffeine exposure on neonate methylation differences in leukocytes were identified in this population with relatively low caffeine consumption.

Differential Methylation and Transcriptome Integration Analysis Identified Differential Methylation Annotation Genes and Functional Research Related to Hair Follicle Development in Sheep

Hair follicle growth and development are a complex and long-term physiological process, which is regulated by a variety of physical factors and signal pathways. Increasing the understanding of the epigenetic regulation and function of candidate genes related to hair follicle development will help to better understand the molecular regulatory mechanisms of hair follicle development. In this study, the methylated DNA immunoprecipitation sequencing (MeDIP-seq) was used to obtain the genome-wide methylation map of the hair follicular development of Super Merino sheep in six stages (fetal skin tissue at 65d, 85d, 105d, 135d, 7d, and 30d after birth). Combined with the results of previous RNA-sequencing, 65 genes were screened out that were both differential methylation and differential expression, including EDN1, LAMC2, NR1D1, RORB, MyOZ3, and WNT2 gene. Differential methylation genes were enriched in Wnt, TNF, TGF-beta, and other signaling pathways related to hair follicle development. The bisulfite sequencing PCR results and MeDIP-seq were basically consistent, indicating that the sequencing results were accurate. As a key gene in the Wnt signaling pathway, both differential methylation and expression gene identified by MeDIP-seq and RNA-seq, further exploration of the function of WNT2 gene revealed that the DNA methylation of exon 5 (CpG11 site) promoted the expression of WNT2 gene. The overexpression vector of lentivirus pLEX-MCS-WNT2 was constructed, and WNT2 gene effectively promoted the proliferation of sheep skin fibroblasts. The results showed that WNT2 gene could promote the growth and development of skin and hair follicles. The results of this study will provide a theoretical basis for further research on sheep hair follicle development and gene regulation mechanisms.

Ten Years of EWAS

Epigenome-wide association study (EWAS) has been applied to analyze DNA methylation variation in complex diseases for a decade, and epigenome as a research target has gradually become a hot topic of current studies. The DNA methylation microarrays, next-generation, and third-generation sequencing technologies have prepared a high-quality platform for EWAS. Here, the progress of EWAS research is reviewed, its contributions to clinical applications, and mainly describe the achievements of four typical diseases. Finally, the challenges encountered by EWAS and make bold predictions for its future development are presented.

Epigenetics Is Implicated in the Basis of Gender Incongruence: An Epigenome-Wide Association Analysis

Introduction

The main objective was to carry out a global DNA methylation analysis in a population with gender incongruence before gender-affirming hormone treatment (GAHT), in comparison to a cisgender population.

Methods

A global CpG (cytosine-phosphate-guanine) methylation analysis was performed on blood from 16 transgender people before GAHT vs. 16 cisgender people using the Illumina© Infinium Human Methylation 850k BeadChip, after bisulfite conversion. Changes in the DNA methylome in cisgender vs. transgender populations were analyzed with the Partek^® Genomics Suite program by a 2-way ANOVA test comparing populations by group and their sex assigned at birth.

Results

The principal components analysis (PCA) showed that both populations (cis and trans) differ in the degree of global CpG methylation prior to GAHT. The 2-way ANOVA test showed 71,515 CpGs that passed the criterion FDR p < 0.05. Subsequently, in male assigned at birth population we found 87 CpGs that passed both criteria (FDR p < 0.05; fold change ≥ ± 2) of which 22 were located in islands. The most significant CpGs were related to genes: WDR45B, SLC6A20, NHLH1, PLEKHA5, UBALD1, SLC37A1, ARL6IP1, GRASP, and NCOA6. Regarding the female assigned at birth populations, we found 2 CpGs that passed both criteria (FDR p < 0.05; fold change ≥ ± 2), but none were located in islands. One of these CpGs, related to the MPPED2 gene, is shared by both, trans men and trans women. The enrichment analysis showed that these genes are involved in functions such as negative regulation of gene expression (GO:0010629), central nervous system development (GO:0007417), brain development (GO:0007420), ribonucleotide binding (GO:0032553), and RNA binding (GO:0003723), among others.

Strengths and Limitations

It is the first time that a global CpG methylation analysis has been carried out in a population with gender incongruence before GAHT. A prospective study before/during GAHT would provide a better understanding of the influence of epigenetics in this process.

Conclusion

The main finding of this study is that the cis and trans populations have different global CpG methylation profiles prior to GAHT. Therefore, our results suggest that epigenetics may be involved in the etiology of gender incongruence.

Origins of human disease: the chrono-epigenetic perspective

Epigenetics has enriched human disease studies by adding new interpretations to disease features that cannot be explained by genetic and environmental factors. However, identifying causal mechanisms of epigenetic origin has been challenging. New opportunities have risen from recent findings in intra-individual and cyclical epigenetic variation, which includes circadian epigenetic oscillations. Cytosine modifications display deterministic temporal rhythms, which may drive ageing and complex disease. Temporality in the epigenome, or the 'chrono' dimension, may help the integration of epigenetic, environmental and genetic disease studies, and reconcile several disparities stemming from the arbitrarily delimited research fields. The ultimate goal of chrono-epigenetics is to predict disease risk, age of onset and disease dynamics from within individual-specific temporal dynamics of epigenomes.

DNA methylation of Vesicular Glutamate Transporters in the mesocorticolimbic brain following early-life stress and adult ethanol exposure-an explorative study

DNA methylation and gene expression can be altered by early life stress (ELS) and/or ethanol consumption. The present study aimed to investigate whether DNA methylation of the Vesicular Glutamate Transporters (Vglut)1-3 is related to previously observed Vglut1-3 transcriptional differences in the ventral tegmental area (VTA), nucleus accumbens (Acb), dorsal striatum (dStr) and medial prefrontal cortex (mPFC) of adult rats exposed to ELS, modelled by maternal separation, and voluntary ethanol consumption. Targeted next-generation bisulfite sequencing was performed to identify the methylation levels on 61 5′-cytosine-phosphate-guanosine-3′ sites (CpGs) in potential regulatory regions of Vglut1, 53 for Vglut2, and 51 for Vglut3. In the VTA, ELS in ethanol-drinking rats was associated with Vglut1-2 CpG-specific hypomethylation, whereas bidirectional Vglut2 methylation differences at single CpGs were associated with ELS alone. Exposure to both ELS and ethanol, in the Acb, was associated with lower promoter and higher intronic Vglut3 methylation; and in the dStr, with higher and lower methylation in 26% and 43% of the analyzed Vglut1 CpGs, respectively. In the mPFC, lower Vglut2 methylation was observed upon exposure to ELS or ethanol. The present findings suggest Vglut1-3 CpG-specific methylation signatures of ELS and ethanol drinking, underlying previously reported Vglut1-3 transcriptional differences in the mesocorticolimbic brain.

A multi-ethnic epigenome-wide association study of leukocyte DNA methylation and blood lipids

Here we examine the association between DNA methylation in circulating leukocytes and blood lipids in a multi-ethnic sample of 16,265 subjects. We identify 148, 35, and 4 novel associations among Europeans, African Americans, and Hispanics, respectively, and an additional 186 novel associations through a trans-ethnic meta-analysis. We observe a high concordance in the direction of effects across racial/ethnic groups, a high correlation of effect sizes between high-density lipoprotein and triglycerides, a modest overlap of associations with epigenome-wide association studies of other cardio-metabolic traits, and a largely non-overlap with lipid loci identified to date through genome-wide association studies. Thirty CpGs reached significance in at least 2 racial/ethnic groups including 7 that showed association with the expression of an annotated gene. CpGs annotated to CPT1A showed evidence of being influenced by triglycerides levels. DNA methylation levels of circulating leukocytes show robust and consistent association with blood lipid levels across multiple racial/ethnic groups.

Methylation of CYP1A1 and VKORC1 promoter associated with stable dosage of warfarin in Chinese patients

Objective

To investigate the association between DNA methylation and the stable warfarin dose through genome-wide DNA methylation analysis and pyrosequencing assay.

Method

This study included 161 patients and genome-wide DNA methylation analysis was used to screen potential warfarin dose-associated CpGs through Illumina Infinium HumanMethylation 450 K BeadChip; then, the pyrosequencing assay was used to further validate the association between the stable warfarin dose and alterations in the methylation of the screened CpGs. GenomeStudio Software and R were used to analyze the differentially methylated CpGs.

Results

The methylation levels of CpGs surrounding the xenobiotic response element (XRE) within the CYP1A1 promoter, differed significantly between the different dose groups (P < 0.05), and these CpGs presented a positive correlation (r> 0, P < 0.05) with an increase in the stable dose of warfarin. At the VKORC1 promoter, two CpGs methylation levels were significantly different between the differential dose groups (P < 0.05), and one CpG (Chr16: 31106793) presented a significant negative correlation (r <  0, P <  0.05) among different dose (low, medium, and high) groups.

Conclusion

This is a novel report of the methylation levels of six CpGs surrounding the XRE within the CYP1A1 promoter and one differential CpG at the VKORC1 promoter associated with stable warfarin dosage; these methylation levels might be applied as molecular signatures for warfarin.

Association of mammographic density with blood DNA methylation

Background: Altered DNA methylation may be an intermediate phenotype between breast cancer risk factors and disease. Mammographic density is a strong risk factor for breast cancer. However, no studies to date have identified an epigenetic signature of mammographic density. We performed an epigenome-wide association study of mammographic density.

Methods: White blood cell DNA methylation was measured for 385 postmenopausal women using the Illumina Infinium MethylationEPIC BeadChip array. Differential methylation was assessed using genome-wide, probe-level, and regional analyses. We implemented a resampling-based approach to improve the stability of our findings.

Results: On average, women with elevated mammographic density exhibited DNA hypermethylation within CpG islands and gene promoters compared to women with lower mammographic density. We identified 250 CpG sites for which DNA methylation was significantly associated with mammographic density. The top sites were located within genes associated with cancer, including HDLBP, TGFB2, CCT4, and PAX8, and were more likely to be located in regulatory regions of the genome. We also identified differential DNA methylation in 37 regions, including within the promoters of PAX8 and PF4, a gene involved in the regulation of angiogenesis. Overall, our results paint a picture of epigenetic dysregulation associated with mammographic density.

Conclusion: Mammographic density is associated with differential DNA methylation throughout the genome, including within genes associated with cancer. Our results suggest the potential involvement of several genes in the biological mechanisms behind differences in breast density between women. Further studies are warranted to explore these potential mechanisms and potential links to breast cancer risk.

Cord blood DNA methylation modifications in infants are associated with white matter microstructure in the context of prenatal maternal depression and anxiety

Maternal and environmental factors influence brain networks and architecture via both physiological pathways and epigenetic modifications. In particular, prenatal maternal depression and anxiety symptoms appear to impact infant white matter (WM) microstructure, leading us to investigate whether epigenetic modifications (i.e., DNA methylation) contribute to these WM differences. To determine if infants of women with depression and anxiety symptoms exhibit epigenetic modifications linked to neurodevelopmental changes, 52 umbilical cord bloods (CBs) were profiled. We observed 219 differentially methylated genomic positions (DMPs; FDR p < 0.05) in CB that were associated with magnetic resonance imaging measures of WM microstructure at 1 month of age and in regions previously described to be related to maternal depression and anxiety symptoms. Genomic characterization of these associated DMPs revealed 143 unique genes with significant relationships to processes involved in neurodevelopment, GTPase activity, or the canonical Wnt signaling pathway. Separate regression models for female (n = 24) and male (n = 28) infants found 142 associated DMPs in females and 116 associated DMPs in males (nominal p value < 0.001, R > 0.5), which were annotated to 98 and 81 genes, respectively. Together, these findings suggest that umbilical CB DNA methylation levels at birth are associated with 1-month WM microstructure.

Neonatal DNA methylation and childhood low prosocial behavior: An epigenome-wide association meta-analysis

Low prosocial behavior in childhood has been consistently linked to later psychopathology, with evidence supporting the influence of both genetic and environmental factors on its development. Although neonatal DNA methylation (DNAm) has been found to prospectively associate with a range of psychological traits in childhood, its potential role in prosocial development has yet to be investigated. This study investigated prospective associations between cord blood DNAm at birth and low prosocial behavior within and across four longitudinal birth cohorts from the Pregnancy And Childhood Epigenetics (PACE) Consortium. We examined (a) developmental trajectories of "chronic-low" versus "typical" prosocial behavior across childhood in a case-control design (N = 2,095), and (b) continuous "low prosocial" scores at comparable cross-cohort time-points (N = 2,121). Meta-analyses were performed to examine differentially methylated positions and regions. At the cohort-specific level, three CpGs were found to associate with chronic low prosocial behavior; however, none of these associations was replicated in another cohort. Meta-analysis revealed no epigenome-wide significant CpGs or regions. Overall, we found no evidence for associations between DNAm patterns at birth and low prosocial behavior across childhood. Findings highlight the importance of employing multi-cohort approaches to replicate epigenetic associations and reduce the risk of false positive discoveries.

The role of DNA methylation in syndromic and non-syndromic congenital heart disease

Congenital heart disease (CHD) is a common structural birth defect worldwide, and defects typically occur in the walls and valves of the heart or enlarged blood vessels. Chromosomal abnormalities and genetic mutations only account for a small portion of the pathogenic mechanisms of CHD, and the etiology of most cases remains unknown. The role of epigenetics in various diseases, including CHD, has attracted increased attention. The contributions of DNA methylation, one of the most important epigenetic modifications, to CHD have not been illuminated. Increasing evidence suggests that aberrant DNA methylation is related to CHD. Here, we briefly introduce DNA methylation and CHD and then review the DNA methylation profiles during cardiac development and in CHD, abnormalities in maternal genome-wide DNA methylation patterns are also described. Whole genome methylation profile and important differentially methylated genes identified in recent years are summarized and clustered according to the sample type and methodologies. Finally, we discuss the novel technology for and prospects of CHD-related DNA methylation.

Human methylome variation across Infinium 450K data on the Gene Expression Omnibus

While DNA methylation (DNAm) is the most-studied epigenetic mark, few recent studies probe the breadth of publicly available DNAm array samples. We collectively analyzed 35 360 Illumina Infinium HumanMethylation450K DNAm array samples published on the Gene Expression Omnibus. We learned a controlled vocabulary of sample labels by applying regular expressions to metadata and used existing models to predict various sample properties including epigenetic age. We found approximately two-thirds of samples were from blood, one-quarter were from brain and one-third were from cancer patients. About 19% of samples failed at least one of Illumina's 17 prescribed quality assessments; signal distributions across samples suggest modifying manufacturer-recommended thresholds for failure would make these assessments more informative. We further analyzed DNAm variances in seven tissues (adipose, nasal, blood, brain, buccal, sperm and liver) and characterized specific probes distinguishing them. Finally, we compiled DNAm array data and metadata, including our learned and predicted sample labels, into database files accessible via the recountmethylation R/Bioconductor companion package. Its vignettes walk the user through some analyses contained in this paper.

Influence of degree of DNA degradation in formalin-fixed and paraffin-embedded tissue samples on accuracy of genome-wide DNA methylation analysis

Aim: Depending upon the degree of DNA degradation of formalin-fixed and paraffin-embedded tissue samples, accuracy of measurement by Infinium MethylationEPIC BeadChip assay (Illumina, CA, USA) was assessed. Materials & methods: DNA quality of six formalin-fixed and paraffin-embedded lung tissue samples with different formalin fixation periods was assessed by Illumina quality control, DNA copy number and DNA integrity number value. Infinium data from restored bisulfite treated DNA were compared with datum from a fresh-frozen sample. Results: The correlation coefficient decreased from 0.993 to 0.970 depending upon DNA degradation, even if the Illumina quality control was met. Exclusion of specific probes improved the correlation regardless of tissue. Conclusion: Poor DNA quality can be assessed as an amplifiable DNA copy number and DNA integrity number value. Probe filtering has the potential to improve assay accuracy.