Critical evaluation of the Illumina MethylationEPIC BeadChip microarray for whole-genome DNA methylation profiling

Epigenome-wide association study of long-term psychosocial stress in older adults

Long-term psychosocial stress is strongly associated with negative physical and mental health outcomes, as well as adverse health behaviours; however, little is known about the role that stress plays on the epigenome. One proposed mechanism by which stress affects DNA methylation is through health behaviours. We conducted an epigenome-wide association study (EWAS) of cumulative psychosocial stress (n = 2,689) from the Health and Retirement Study (mean age = 70.4 years), assessing DNA methylation (Illumina Infinium HumanMethylationEPIC Beadchip) at 789,656 CpG sites. For identified CpG sites, we conducted a formal mediation analysis to examine whether smoking, alcohol use, physical activity, and body mass index (BMI) mediate the relationship between stress and DNA methylation. Nine CpG sites were associated with psychosocial stress (all p < 9E-07; FDR q < 0.10). Additionally, health behaviours and/or BMI mediated 9.4% to 21.8% of the relationship between stress and methylation at eight of the nine CpGs. Several of the identified CpGs were in or near genes associated with cardiometabolic traits, psychosocial disorders, inflammation, and smoking. These findings support our hypothesis that psychosocial stress is associated with DNA methylation across the epigenome. Furthermore, specific health behaviours mediate only a modest percentage of this relationship, providing evidence that other mechanisms may link stress and DNA methylation.

Widespread genomic de novo DNA methylation occurs following CD8+ T cell activation and proliferation

This research investigates the intricate dynamics of DNA methylation in the hours following CD8+ T cell activation, during a critical yet understudied temporal window. DNA methylation is an epigenetic modification central to regulation of gene expression and directing immune responses. Our investigation spanned 96-h post-activation and unveils a nuanced tapestry of global and site-specific methylation changes. We identified 15,626 significant differentially methylated CpGs spread across the genome, with the most significant changes occurring within the genes ADAM10, ICA1, and LAPTM5. While many changes had modest effect sizes, approximately 120 CpGs exhibited a log2FC above 1.5, with cell activation and proliferation pathways the most affected. Relatively few of the differentially methylated CpGs occurred along adjacent gene regions. The exceptions were seven differentially methylated gene regions, with the Human T cell Receptor Alpha Joining Genes demonstrating consistent methylation change over a 3kb window. We also investigated whether an inflammatory environment could alter DNA methylation during activation, with proliferating cells exposed to the oxidant glycine chloramine. No substantial differential methylation was observed in this context. The temporal perspective of early activation adds depth to the evolving field of epigenetic immunology, offering insights with implications for therapeutic innovation and expanding our understanding of epigenetic modulation in immune function.

Critical evaluation of the reliability of DNA methylation probes on the Illumina MethylationEPIC v1.0 BeadChip microarrays

DNA methylation (DNAm) plays a crucial role in a number of complex diseases. However, the reliability of DNAm levels measured using Illumina arrays varies across different probes. Previous research primarily assessed probe reliability by comparing duplicate samples between the 450k-450k or 450k-EPIC platforms, with limited investigations on Illumina EPIC v1.0 arrays. We conducted a comprehensive assessment of the EPIC v1.0 array probe reliability using 69 blood DNA samples, each measured twice, generated by the Alzheimer's Disease Neuroimaging Initiative study. We observed higher reliability in probes with average methylation beta values of 0.2 to 0.8, and lower reliability in type I probes or those within the promoter and CpG island regions. Importantly, we found that probe reliability has significant implications in the analyses of Epigenome-wide Association Studies (EWAS). Higher reliability is associated with more consistent effect sizes in different studies, the identification of differentially methylated regions (DMRs) and methylation quantitative trait locus (mQTLs), and significant correlations with downstream gene expression. Moreover, blood DNAm measurements obtained from probes with higher reliability are more likely to show concordance with brain DNAm measurements. Our findings, which provide crucial reliability information for probes on the EPIC v1.0 array, will serve as a valuable resource for future DNAm studies.

A review on trends in development and translation of omics signatures in cancer

The field of cancer genomics and transcriptomics has evolved from targeted profiling to swift sequencing of individual tumor genome and transcriptome. The steady growth in genome, epigenome, and transcriptome datasets on a genome-wide scale has significantly increased our capability in capturing signatures that represent both the intrinsic and extrinsic biological features of tumors. These biological differences can help in precise molecular subtyping of cancer, predicting tumor progression, metastatic potential, and resistance to therapeutic agents. In this review, we summarized the current development of genomic, methylomic, transcriptomic, proteomic and metabolic signatures in the field of cancer research and highlighted their potentials in clinical applications to improve diagnosis, prognosis, and treatment decision in cancer patients.

Epigenetic landscape of 5-hydroxymethylcytosine and associations with gene expression in placenta

5-hydroxymethylcystosine (5hmC), is an intermediate product in the DNA demethylation pathway, but may act as a functional epigenetic modification. We have conducted the largest study of site-specific 5hmC in placenta to date using parallel bisulphite and oxidative bisulphite modification with array-based assessment. Incorporating parallel RNA-sequencing data allowed us to assess associations between 5hmC and gene expression, using expression quantitative trait hydroxymethylation (eQTHM) analysis. We identified ~ 47,000 loci with consistently elevated (systematic) 5hmC proportions. Systematic 5hmC was significantly depleted (p < 0.0001) at CpG islands (CGI), and enriched (p < 0.0001) in 'open sea' regions (CpG >4 kb from CGI). 5hmC was most and least abundant at CpGs in enhancers and active transcription start sites (TSS), respectively (p < 0.05). We identified 499 significant (empirical-p <0.05) eQTHMs within 1 MB of the assayed gene. At most (75.4%) eQTHMs, the proportion of 5hmC was positively correlated with transcript abundance. eQTHMs were significantly enriched among enhancer CpGs and depleted among CpGs in active TSS (p < 0.05 for both). Finally, we identified 107 differentially hydroxymethylated regions (DHMRs, p < 0.05) across 100 genes. Our study provides insight into placental distribution of 5hmC, and sheds light on the functional capacity of this epigenetic modification in placenta.

Associations of Heavy Metals with Cognitive Function: An Epigenome-Wide View of DNA Methylation and Mediation Analysis

OBJECTIVE

Exposure to heavy metals has been reported to be associated with impaired cognitive function, but the underlying mechanisms remain unclear. This pilot study aimed to identify key heavy metal elements associated with cognitive function and further explore the potential mediating role of metal-related DNA methylation.

METHODS

Blood levels of arsenic, cadmium, lead, copper, manganese, and zinc and genome-wide DNA methylations were separately detected in peripheral blood in 155 older adults. Cognitive function was evaluated using the Mini-Mental State Examination (MMSE). Least absolute shrinkage and selection operator penalized regression and Bayesian kernel machine regression were used to identify metals associated with cognitive function. An epigenome-wide association study examined the DNA methylation profile of the identified metal, and mediation analysis investigated its mediating role.

RESULTS

The MMSE scores showed a significant decrease of 1.61 (95% confidence interval [CI]: -2.64, -0.59) with each 1 standard deviation increase in ln-transformed arsenic level; this association was significant in multiple-metal models and dominated the overall negative effect of 6 heavy metal mixture on cognitive function. Seventy-three differentially methylated positions were associated with blood arsenic (p < 1.0 × 10-5). The methylation levels at cg05226051 (annotated to TDRD3) and cg18886932 (annotated to GAL3ST3) mediated 24.8% and 25.5% of the association between blood arsenic and cognitive function, respectively (all p < 0.05).

INTERPRETATION

Blood arsenic levels displayed a negative association with the cognitive function of older adults. This finding shows that arsenic-related DNA methylation alterations are critical partial mediators that may serve as potential biomarkers for further mechanism-related studies. ANN NEUROL 2024;96:87-98.

Identification of a DNA-methylome-based signature for prognosis prediction in driver gene-negative lung adenocarcinoma

"Driver gene-negative" lung adenocarcinoma (LUAD) was of rare treatment options and a poor prognosis. Presently, for them, few biomarkers are available for stratification analysis to make appropriate treatment strategy. This study aimed to develop a DNA-methylome-based signature to realize the precise risk-stratifying. Here, an Illumina MethylationEPIC Beadchip was applied to obtain differentially methylated CpG sites (DMCs). A four-CpG-based signature, named as TLA, was successfully established, whose prognosis-predicting power was well verified in one internal (n = 78) and other external (n = 110) validation cohorts. Patients with high-risk scores had shorter overall survival (OS) in all cohorts [hazard ratio (HR): 11.79, 5.16 and 2.99, respectively]. Additionally, it can effectively divide patients into low-risk and high-risk groups, with significantly different OS in the diverse subgroups stratified by the standard clinical parameters. As an independent prognostic factor, TLA may assist in improving the nomogram's 5-year OS-predicting ability (AUC 0.756, 95% CI:0.695-0.816), superior to TNM alone (AUC 0.644, 95% CI: 0.590-0.698). Additionally, the relationship of TLA-related genes, TAC1, LHX9, and ALX1, with prognosis and tumour invasion made them serve as potential therapy targets for driver gene-negative LUAD.

Epigenetic footprints: Investigating placental DNA methylation in the context of prenatal exposure to phenols and phthalates

BACKGROUND

Endocrine disrupting compounds (EDCs) such as phthalates and phenols can affect placental functioning and fetal health, potentially via epigenetic modifications. We investigated the associations between pregnancy exposure to synthetic phenols and phthalates estimated from repeated urine sampling and genome wide placental DNA methylation.

METHODS

The study is based on 387 women with placental DNA methylation assessed with Infinium MethylationEPIC arrays and with 7 phenols, 13 phthalates, and two non-phthalate plasticizer metabolites measured in pools of urine samples collected twice during pregnancy. We conducted an exploratory analysis on individual CpGs (EWAS) and differentially methylated regions (DMRs) as well as a candidate analysis focusing on 20 previously identified CpGs. Sex-stratified analyses were also performed.

RESULTS

In the exploratory analysis, when both sexes were studied together no association was observed in the EWAS. In the sex-stratified analysis, 114 individual CpGs (68 in males, 46 in females) were differentially methylated, encompassing 74 genes (36 for males and 38 for females). We additionally identified 28 DMRs in the entire cohort, 40 for females and 42 for males. Associations were mostly positive (for DMRs: 93% positive associations in the entire cohort, 60% in the sex-stratified analysis), with the exception of several associations for bisphenols and DINCH metabolites that were negative. Biomarkers associated with most DMRs were parabens, DEHP, and DiNP metabolite concentrations. Some DMRs encompassed imprinted genes including APC (associated with parabens and DiNP metabolites), GNAS (bisphenols), ZIM2;PEG3;MIMT1 (parabens, monoethyl phthalate), and SGCE;PEG10 (parabens, DINCH metabolites). Terms related to adiposity, lipid and glucose metabolism, and cardiovascular function were among the enriched phenotypes associated with differentially methylated CpGs. The candidate analysis identified one CpG mapping to imprinted LGALS8 gene, negatively associated with ethylparaben.

CONCLUSIONS

By combining improved exposure assessment and extensive placental epigenome coverage, we identified several novel genes associated with the exposure, possibly in a sex-specific manner.

Short-Term PM2.5 Exposure and DNA Methylation Changes of Circadian Rhythm Genes: Evidence from Two Experimental Studies

The circadian rhythm regulates many crucial physiological processes, impacting human aging and aging-related outcomes. Observational evidence links circadian rhythm disturbance to PM2.5 exposure, yet the underlying DNA methylation mechanisms remain unclear due to limited PM2.5-dominated experimental settings. Therefore, we investigated the associations between short-term PM2.5 exposure and DNA methylation changes of 1188 CpG candidates across circadian genes among 32 young adults in the FDU study, with the validation in 26 individuals from the PKU study. Further mediation analyses tested whether DNA methylation of circadian genes could mediate the influence of PM2.5 on aging measured by three epigenetic ages: DNAmGrimAge, DunedinPoAm, and the mortality risk score. We identified three CpG sites associated with personal PM2.5 exposure: cg01248361 (CSNK2A2), cg17728065 (RORA), and cg22513396 (PRKAG2). Acute effects of PM2.5 on the three loci could be mediated by several circulating biomarkers, including MDA and EGF, with up to ∼30% of mediated proportions. Three loci further showed varying potentials in mediating the aging acceleration effect of PM2.5. Locus cg17728065 is the key site exhibiting a robust mediating effect (7.54-12.52%) on PM2.5-induced aging acceleration. Our findings demonstrated that PM2.5, even short-term peaks, could leave imprints on human aging via inducing aberrant temporal fluctuation in circadian homeostasis captured by DNA methylation profiles.

Patient-derived organoid biobank identifies epigenetic dysregulation of intestinal epithelial MHC-I as a novel mechanism in severe Crohn's Disease

OBJECTIVE

Epigenetic mechanisms, including DNA methylation (DNAm), have been proposed to play a key role in Crohn's disease (CD) pathogenesis. However, the specific cell types and pathways affected as well as their potential impact on disease phenotype and outcome remain unknown. We set out to investigate the role of intestinal epithelial DNAm in CD pathogenesis.

DESIGN

We generated 312 intestinal epithelial organoids (IEOs) from mucosal biopsies of 168 patients with CD (n=72), UC (n=23) and healthy controls (n=73). We performed genome-wide molecular profiling including DNAm, bulk as well as single-cell RNA sequencing. Organoids were subjected to gene editing and the functional consequences of DNAm changes evaluated using an organoid-lymphocyte coculture and a nucleotide-binding oligomerisation domain, leucine-rich repeat and CARD domain containing 5 (NLRC5) dextran sulphate sodium (DSS) colitis knock-out mouse model.

RESULTS

We identified highly stable, CD-associated loss of DNAm at major histocompatibility complex (MHC) class 1 loci including NLRC5 and cognate gene upregulation. Single-cell RNA sequencing of primary mucosal tissue and IEOs confirmed the role of NLRC5 as transcriptional transactivator in the intestinal epithelium. Increased mucosal MHC-I and NLRC5 expression in adult and paediatric patients with CD was validated in additional cohorts and the functional role of MHC-I highlighted by demonstrating a relative protection from DSS-mediated mucosal inflammation in NLRC5-deficient mice. MHC-I DNAm in IEOs showed a significant correlation with CD disease phenotype and outcomes. Application of machine learning approaches enabled the development of a disease prognostic epigenetic molecular signature.

CONCLUSIONS

Our study has identified epigenetically regulated intestinal epithelial MHC-I as a novel mechanism in CD pathogenesis.

Systemic inflammation, oxidative damage and neurocognition predict telomere length in a transdiagnostic sample stratified by global DNA methylation levels

Epigenetic mechanisms contribute to the maintenance of both type 2 diabetes mellitus (T2DM) and psychiatric disorders. Emerging evidence suggests that molecular pathways and neurocognitive performance regulate epigenetic dynamics in these disorders. The current combined and transdiagnostic study investigated whether inflammatory, oxidative stress, adhesion molecule, neurocognitive and functional performance are significant predictors of telomere dynamics in a sample stratified by global DNA methylation levels. Peripheral blood inflammation, oxidative stress and adhesion molecule biomarkers and neurocognitive function were assessed twice over a 1-year period in 80 individuals, including 16 with schizophrenia (SZ), 16 with bipolar disorder (BD), 16 with major depressive disorder (MDD), 15 with T2DM, and 17 healthy controls (HCs). Leukocyte telomere length (LTL) was measured by qRT-PCR using deoxyribonucleic acid (DNA) extracted from peripheral blood samples. A posteriori, individuals were classified based on their global methylation score (GMS) at baseline into two groups: the below-average methylation (BM) and above-average methylation (AM) groups. Hierarchical and k-means clustering methods, mixed one-way analysis of variance and linear regression analyses were performed. Overall, the BM group showed a significantly higher leukocyte telomere length (LTL) than the AM group at both time points (p = 0.02; η2p = 0.06). Moreover, the BM group had significantly lower levels of tumor necrosis factor alpha (TNF-α) (p = 0.03; η2p = 0.06) and C-reactive protein (CRP) (p = 0.03; η2p = 0.06) than the AM group at the 1-year follow-up. Across all participants, the regression models showed that oxidative stress (reactive oxygen species [ROS]) (p = 0.04) and global cognitive score [GCS] (p = 0.02) were significantly negatively associated with LTL, whereas inflammatory (TNF-α) (p = 0.04), adhesion molecule biomarkers (inter cellular adhesion molecule [ICAM]) (p = 0.009), and intelligence quotient [IQ] (p = 0.03) were significantly positively associated with LTL. Moreover, the model predictive power was increased when tested in both groups separately, explaining 15.8% and 28.1% of the LTL variance at the 1-year follow-up for the AM and BM groups, respectively. Heterogeneous DNA methylation in individuals with T2DM and severe mental disorders seems to support the hypothesis that epigenetic dysregulation occurs in a transdiagnostic manner. Our results may help to elucidate the interplay between epigenetics, molecular processes and neurocognitive function in these disorders. DNA methylation and LTL are potential therapeutic targets for transdiagnostic interventions to decrease the risk of comorbidities.

Transcriptome-wide gene expression outlier analysis pinpoints therapeutic vulnerabilities in colorectal cancer

Multiple strategies are continuously being explored to expand the drug target repertoire in solid tumors. We devised a novel computational workflow for transcriptome-wide gene expression outlier analysis that allows the systematic identification of both overexpression and underexpression events in cancer cells. Here, it was applied to expression values obtained through RNA sequencing in 226 colorectal cancer (CRC) cell lines that were also characterized by whole-exome sequencing and microarray-based DNA methylation profiling. We found cell models displaying an abnormally high or low expression level for 3533 and 965 genes, respectively. Gene expression abnormalities that have been previously associated with clinically relevant features of CRC cell lines were confirmed. Moreover, by integrating multi-omics data, we identified both genetic and epigenetic alternations underlying outlier expression values. Importantly, our atlas of CRC gene expression outliers can guide the discovery of novel drug targets and biomarkers. As a proof of concept, we found that CRC cell lines lacking expression of the MTAP gene are sensitive to treatment with a PRMT5-MTA inhibitor (MRTX1719). Finally, other tumor types may also benefit from this approach.

Association of biological aging with frailty and post-transplant outcomes among adults with cirrhosis

Frailty is classically associated with advanced age but is also an important predictor of clinical outcomes in comparatively young adults with cirrhosis. We examined the association of biological aging with frailty and post-transplant outcomes in a pilot of adults with cirrhosis undergoing liver transplantation (LT). Frailty was measured via the Liver Frailty Index (LFI). The primary epigenetic clock DNA methylation (DNAm) PhenoAge was calculated from banked peripheral blood mononuclear cells; we secondarily explored two first-generation clocks (Hannum; Horvath) and two additional second-generation clocks (GrimAge; GrimAge2). Twelve adults were included: seven frail (LFI ≥ 4.4, mean age 55 years) and five robust (LFI < 3.2, mean age 55 years). Mean PhenoAge age acceleration (AgeAccel) was + 2.5 years (P = 0.23) for frail versus robust subjects. Mean PhenoAge AgeAccel was + 2.7 years (P = 0.19) for subjects who were readmitted or died within 30 days of discharge post-LT versus those without this outcome. When compared with first-generation clocks, the second-generation clocks demonstrated greater average AgeAccel for subjects with frailty or poor post-LT outcomes. Measuring biological age using DNAm-derived epigenetic clocks is feasible in adults undergoing LT. While frail and robust subjects had the same average chronological age, average biological age as measured by second-generation epigenetic clocks tended to be accelerated among those who were frail or experienced a poor post-LT outcome. These results suggest that frailty in these relatively young subjects with cirrhosis may involve similar aging mechanisms as frailty classically observed in chronologically older adults and warrant validation in a larger cohort.

Childhood adversity, accelerated GrimAge, and associated health consequences

Childhood adversity is linked to psychological, behavioral, and physical health problems, including obesity and cardiometabolic disease. Epigenetic alterations are one pathway through which the effects of early life stress and adversity might persist into adulthood. Epigenetic mechanisms have also been proposed to explain why cardiometabolic health can vary greatly between individuals with similar Body Mass Index (BMIs). We evaluated two independent cross-sectional cohorts of adults without known medical illness, one of which explicitly recruited individuals with early life stress (ELS) and control participants (n = 195), and the other a general community sample (n = 477). In these cohorts, we examine associations between childhood adversity, epigenetic aging, and metabolic health. Childhood adversity was associated with increased GrimAge Acceleration (GAA) in both cohorts, both utilizing a dichotomous yes/no classification (both p < 0.01) as well as a continuous measure using the Childhood Trauma Questionnaire (CTQ) (both p < 0.05). Further investigation demonstrated that CTQ subscales for physical and sexual abuse (both p < 0.05) were associated with increased GAA in both cohorts, whereas physical and emotional neglect were not. In both cohorts, higher CTQ was also associated with higher BMI and increased insulin resistance (both p < 0.05). Finally, we demonstrate a moderating effect of BMI on the relationship between GAA and insulin resistance where GAA correlated with insulin resistance specifically at higher BMIs. These results, which were largely replicated between two independent cohorts, suggest that interactions between epigenetics, obesity, and metabolic health may be important mechanisms through which childhood adversity contributes to long-term physical and metabolic health effects.

Benchmarking of methods for DNA methylome deconvolution

Defining the number and abundance of different cell types in tissues is important for understanding disease mechanisms as well as for diagnostic and prognostic purposes. Typically, this is achieved by immunohistological analyses, cell sorting, or single-cell RNA-sequencing. Alternatively, cell-specific DNA methylome information can be leveraged to deconvolve cell fractions from a bulk DNA mixture. However, comprehensive benchmarking of deconvolution methods and modalities was not yet performed. Here we evaluate 16 deconvolution algorithms, developed either specifically for DNA methylome data or more generically. We assess the performance of these algorithms, and the effect of normalization methods, while modeling variables that impact deconvolution performance, including cell abundance, cell type similarity, reference panel size, method for methylome profiling (array or sequencing), and technical variation. We observe differences in algorithm performance depending on each these variables, emphasizing the need for tailoring deconvolution analyses. The complexity of the reference, marker selection method, number of marker loci and, for sequencing-based assays, sequencing depth have a marked influence on performance. By developing handles to select the optimal analysis configuration, we provide a valuable source of information for studies aiming to deconvolve array- or sequencing-based methylation data.

DNA methylation signatures of youth-onset type 2 diabetes and exposure to maternal diabetes

OBJECTIVE

Youth-onset type 2 diabetes (T2D) is physiologically distinct from adult-onset, but it is not clear how the two diseases differ at a molecular level. In utero exposure to maternal type 2 diabetes (T2D) is known to be a specific risk factor for youth-onset T2D. DNA methylation (DNAm) changes associated with T2D but which differ between youth- and adult-onset might delineate the impacts of T2D development at different ages and could also determine the contribution of exposure to in utero diabetes.

METHODS

We performed an epigenome-wide analysis of DNAm on whole blood from 218 youth with T2D and 77 normoglycemic controls from the iCARE (improving renal Complications in Adolescents with type 2 diabetes through REsearch) cohort. Associations were tested using multiple linear regression models while adjusting for maternal diabetes, sex, age, BMI, smoking status, second-hand smoking exposure, cell-type proportions and genetic ancestry.

RESULTS

We identified 3830 differentially methylated sites associated with youth T2D onset, of which 3794 were moderately (adjusted p-value < 0.05 and effect size estimate > 0.01) associated and 36 were strongly (adjusted p-value < 0.05 and effect size estimate > 0.05) associated. A total of 3725 of these sites were not previously reported in the EWAS Atlas as associated with T2D, adult obesity or youth obesity. Moreover, three CpGs associated with youth-onset T2D in the PFKFB3 gene were also associated with maternal T2D exposure (FDR < 0.05 and effect size > 0.01). This is the first study to link PFKFB3 and T2D in youth.

CONCLUSION

Our findings support that T2D in youth has different impacts on DNAm than adult-onset, and suggests that changes in DNAm could provide an important link between in utero exposure to maternal diabetes and the onset of T2D.

Placental PFAS concentrations are associated with perturbations of placental DNA methylation at loci with important roles on cardiometabolic health

The placenta is crucial for fetal development, is affected by PFAS toxicity, and evidence is accumulating that gestational PFAS perturb the epigenetic activity of the placenta. Gestational PFAS exposure is can adversely affect offspring, yet individual and cumulative impacts of PFAS on the placental epigenome remain underexplored. Here, we conducted an epigenome-wide association study (EWAS) to examine the relationships between placental PFAS levels and DNA methylation in a cohort of mother-infant dyads in Arkansas. We measured 17 PFAS in human placental tissues and quantified placental DNA methylation levels via the Illumina EPIC Microarray. We tested for differential DNA methylation with individual PFAS, and with mixtures of multiple PFAS. Our results demonstrated that numerous epigenetic loci were perturbed by PFAS, with PFHxS exhibiting the most abundant effects. Mixture analyses suggested cumulative effects of PFOA and PFOS, while PFHxS may act more independently. We additionally explored whether sex-specific effects may be present and concluded that future large studies should explicitly test for sex-specific effects. The genes that are annotated to our PFAS-associated epigenetic loci are primarily involved in growth processes and cardiometabolic health, while some genes are involved in neurodevelopment. These findings shed light on how prenatal PFAS exposures affect birth outcomes and children's health, emphasizing the importance of understanding PFAS mechanisms in the in-utero environment.

Genetic variation and pesticide exposure influence blood DNA methylation signatures in females with early-stage Parkinson's disease

Although sex, genetics, and exposures can individually influence risk for sporadic Parkinson's disease (PD), the joint contributions of these factors to the epigenetic etiology of PD have not been comprehensively assessed. Here, we profiled sex-stratified genome-wide blood DNAm patterns, SNP genotype, and pesticide exposure in agricultural workers (71 early-stage PD cases, 147 controls) and explored replication in three independent samples of varying demographics (n = 218, 222, and 872). Using a region-based approach, we found more associations of blood DNAm with PD in females (69 regions) than in males (2 regions, Δβadj| ≥0.03, padj ≤ 0.05). For 48 regions in females, models including genotype or genotype and pesticide exposure substantially improved in explaining interindividual variation in DNAm (padj ≤ 0.05), and accounting for these variables decreased the estimated effect of PD on DNAm. The results suggested that genotype, and to a lesser degree, genotype-exposure interactions contributed to variation in PD-associated DNAm. Our findings should be further explored in larger study populations and in experimental systems, preferably with precise measures of exposure.

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.

Genome-wide DNA methylation-analysis of blastic plasmacytoid dendritic cell neoplasm identifies distinct molecular features

Blastic plasmacytoid dendritic cell neoplasm (BPDCN) constitutes a rare and aggressive malignancy originating from plasmacytoid dendritic cells (pDCs) with a primarily cutaneous tropism followed by dissemination to the bone marrow and other organs. We conducted a genome-wide analysis of the tumor methylome in an extended cohort of 45 BPDCN patients supplemented by WES and RNA-seq as well as ATAC-seq on selected cases. We determined the BPDCN DNA methylation profile and observed a dramatic loss of DNA methylation during malignant transformation from early and mature DCs towards BPDCN. DNA methylation profiles further differentiate between BPDCN, AML, CMML, and T-ALL exhibiting the most striking global demethylation, mitotic stress, and merely localized DNA hypermethylation in BPDCN resulting in pronounced inactivation of tumor suppressor genes by comparison. DNA methylation-based analysis of the tumor microenvironment by MethylCIBERSORT yielded two, prognostically relevant clusters (IC1 and IC2) with specific cellular composition and mutational spectra. Further, the transcriptional subgroups of BPDCN (C1 and C2) differ by DNA methylation signatures in interleukin/inflammatory signaling genes but also by higher transcription factor activity of JAK-STAT and NFkB signaling in C2 in contrast to an EZH2 dependence in C1-BPDCN. Our integrative characterization of BPDCN offers novel molecular insights and potential diagnostic applications.

Integrated Epigenetic and Transcriptomic Analysis Identifies Interleukin 17 DNA Methylation Signature of Malignant Peripheral Nerve Sheath Tumor Progression and Metastasis

PURPOSE

Sarcomas are a complex group of highly aggressive and metastatic tumors with over 100 distinct subtypes. Because of their diversity and rarity, it is challenging to generate multisarcoma signatures that are predictive of patient outcomes.

MATERIALS AND METHODS

Here, we identify a DNA methylation signature for progression and metastasis of numerous sarcoma subtypes using multiple epigenetic and genomic patient data sets. Malignant Peripheral Nerve Sheath Tumors (MPNSTs) are highly metastatic sarcomas with frequent loss of the histone methyltransferase, PRC2. Loss of PRC2 is associated with MPNST metastasis and plays a critical noncanonical role in DNA methylation.

RESULTS

We found that over 900 5'-C-phosphate-G-3' (CpGs) were hypermethylated in MPNSTs with PRC2 loss. Furthermore, we identified eight differentially methylated CpGs in the IL17D/RD family that correlate with the progression and metastasis of MPNSTs in two independent patient data sets. Similar trends were identified in other sarcoma subtypes, including osteosarcoma, rhabdomyosarcoma, and synovial sarcoma. Analysis of scRNAseq data sets determined that IL17D/RD expression occurs in both the tumor cells and the surrounding stromal populations.

CONCLUSION

These results might have broad implications for the clinical management and surveillance of sarcoma.

Differences in DNA methylation status explain phenotypic variability in patients with 5p- syndrome

Cri Du Chat syndrome, or 5p- syndrome, is characterized by a terminal or interstitial deletion on the short arm of chromosome 5 that causes variable clinical manifestations, including high-pitched cry in newborns, delayed growth, and global development. Different cytogenomic rearrangements, family history, and environmental factors may hinder the genotype-phenotype association. Thus, the phenotypic variability of this syndrome may not be limited only to variations in gene structure, such as deletions and duplications. It is possible that other mechanisms related to the activation or inactivation of promoters and/or exons of actively transcribed genes, such as DNA methylation are involved. Therefore, we studied the genome-wide methylation status profile of peripheral blood samples from fifteen patients with Cri du Chat Syndrome and nine control samples through the array method to look for Differentially Methylated Regions. We found that Differentially Methylated Regions outside the 5p region are mainly associated with regulating gene transcription, splicing, and chromatin remodeling. Most biological pathways are related to transcription, histone and chromatin binding, spliceosome and ribosomal complex, and RNA processing. Our results suggest that changes in the 5p region can cause an imbalance in other chromosomal regions capable of affecting gene modulation and thus explain the phenotypic differences in patients with 5p-.

Methylation profiles at birth linked to early childhood obesity

Childhood obesity represents a significant global health concern and identifying its risk factors is crucial for developing intervention programs. Many "omics" factors associated with the risk of developing obesity have been identified, including genomic, microbiomic, and epigenomic factors. Here, using a sample of 48 infants, we investigated how the methylation profiles in cord blood and placenta at birth were associated with weight outcomes (specifically, conditional weight gain, body mass index, and weight-for-length ratio) at age six months. We characterized genome-wide DNA methylation profiles using the Illumina Infinium MethylationEpic chip, and incorporated information on child and maternal health, and various environmental factors into the analysis. We used regression analysis to identify genes with methylation profiles most predictive of infant weight outcomes, finding a total of 23 relevant genes in cord blood and 10 in placenta. Notably, in cord blood, the methylation profiles of three genes (PLIN4, UBE2F, and PPP1R16B) were associated with all three weight outcomes, which are also associated with weight outcomes in an independent cohort suggesting a strong relationship with weight trajectories in the first six months after birth. Additionally, we developed a Methylation Risk Score (MRS) that could be used to identify children most at risk for developing childhood obesity. While many of the genes identified by our analysis have been associated with weight-related traits (e.g., glucose metabolism, BMI, or hip-to-waist ratio) in previous genome-wide association and variant studies, our analysis implicated several others, whose involvement in the obesity phenotype should be evaluated in future functional investigations.

Developmental epigenomic effects of maternal financial problems

Early-life adversity as neglect or low socioeconomic status is associated with negative physical/mental health outcomes and plays an important role in health trajectories through life. The early-life environment has been shown to be encoded as changes in epigenetic markers that are retained for many years.We investigated the effect of maternal major financial problems (MFP) and material deprivation (MD) on their children's epigenome in the Avon Longitudinal Study of Parents and Children (ALSPAC) cohort. Epigenetic aging, measured with epigenetic clocks, was weakly accelerated with increased MFP. In subsequent EWAS, MFP, and MD showed strong, independent programing effects on children's genomes. MFP in the period from birth to age seven was associated with genome-wide epigenetic modifications on children's genome visible at age 7 and partially remaining at age 15.These results support the hypothesis that physiological processes at least partially explain associations between early-life adversity and health problems later in life. Both maternal stressors (MFP/MD) had similar effects on biological pathways, providing preliminary evidence for the mechanisms underlying the effects of low socioeconomic status in early life and disease outcomes later in life. Understanding these associations is essential to explain disease susceptibility, overall life trajectories and the transition from health to disease.

Methylation-driven mechanisms of allergic rhinitis during pollen and non-pollen seasons using integrated bioinformatics analysis

Background

Allergic rhinitis (AR) is a widespread allergic airway disease that results from a complex interplay between genetic and environmental factors and affects approximately 10%-40% of the global population. Pollen is a common allergen, and exposure to pollen can cause epigenetic changes. However, the mechanism underlying pollen-induced DNA methylation changes and their potential effects on the allergic march are still unclear. The purpose of this study was to explore the methylation-driven mechanisms of AR during the pollen and non-pollen seasons using bioinformatics analysis and to investigate their relationship with asthma.

Methods

We downloaded DNA methylation and gene expression data from the GEO database (GSE50387: GSE50222, GSE50101) and identified differentially methylated positions (DMPs) and differentially expressed genes (DEGs) during the pollen and non-pollen seasons using the CHAMP and limma packages. Through correlation analysis, we identified methylation-driven genes and performed pathway enrichment analysis to annotate their functions. We incorporated external data on AR combined with asthma (GSE101720) for analysis to identify key CpGs that promote the transformation of AR to asthma. We also utilized external data on olive pollen allergy (GSE54522) for analysis to validate the methylation-driven genes. Weighted correlation network analysis (WGCNA) was employed to identify gene modules significantly correlated with pollen allergy. We extracted genes related to the key methylation-driven gene ZNF667-AS1 from the significant module and performed pathway intelligent clustering using KOBAS-i. We also utilized gene set enrichment analysis to explore the potential function of ZNF667-AS1.

Results

We identified 20 and 24 CpG-Gene pairings during the pollen and non-pollen seasons. After incorporating external data from GSE101720, we found that ZNF667-AS1 is a key gene that may facilitate the transformation of AR into asthma during the pollen season. This finding was further validated in another external dataset, GSE54522, which is associated with pollen allergy. WGCNA identified 17 modules, among which the blue module showed significant correlation with allergies. ZNF667-AS1 was located in the blue module. We performed pathway analysis on the genes correlated with ZNF667-AS1 extracted from the blue module and identified a prominent cluster of pathways in the KOBAS-i results, including Toll-like receptor (TLR) family, MyD88, MAPK, and oxidative stress. Gene set enrichment analysis around cg05508084 (paired with ZNF667-AS1) also indicated its potential involvement in initiating and modulating allergic inflammation from the perspective of TLR and MAPK signaling.

Conclusion

We identified methylation-driven genes and their related pathways during the pollen and non-pollen seasons in patients with AR and identified key CpGs that promote the transformation of AR into asthma due to pollen exposure. This study provides new insights into the underlying molecular mechanisms of the transformation of AR to asthma.

Ultrasensitive amplification-free quantification of a methyl CpG-rich cancer biomarker by single-molecule kinetic fingerprinting

The most well-studied epigenetic marker in humans is the 5-methyl modification of cytosine in DNA, which has great potential as a disease biomarker in liquid biopsies of cell-free DNA. Currently, quantification of DNA methylation relies heavily on bisulfite conversion followed by PCR amplification and NGS or microarray analysis. PCR is subject to potential bias in differential amplification of bisulfite-converted methylated versus unmethylated sequences. Here, we combine bisulfite conversion with single-molecule kinetic fingerprinting to develop an amplification-free assay for DNA methylation at the branched-chain amino acid transaminase 1 (BCAT1) promoter. Our assay selectively responds to methylated sequences with a limit of detection below 1 fM and a specificity of 99.9999%. Evaluating complex genomic DNA matrices, we reliably distinguish 2-5% DNA methylation at the BCAT1 promoter in whole blood DNA from completely unmethylated whole-genome amplified DNA. Taken together, these results demonstrate the feasibility and sensitivity of our amplification-free, single-molecule quantification approach to improve the early detection of methylated cancer DNA biomarkers.

The association of cigarette smoking with DNA methylation and gene expression in human tissue samples

Cigarette smoking adversely affects many aspects of human health, and epigenetic responses to smoking may reflect mechanisms that mediate or defend against these effects. Prior studies of smoking and DNA methylation (DNAm), typically measured in leukocytes, have identified numerous smoking-associated regions (e.g., AHRR). To identify smoking-associated DNAm features in typically inaccessible tissues, we generated array-based DNAm data for 916 tissue samples from the GTEx (Genotype-Tissue Expression) project representing 9 tissue types (lung, colon, ovary, prostate, blood, breast, testis, kidney, and muscle). We identified 6,350 smoking-associated CpGs in lung tissue (n = 212) and 2,735 in colon tissue (n = 210), most not reported previously. For all 7 other tissue types (sample sizes 38-153), no clear associations were observed (false discovery rate 0.05), but some tissues showed enrichment for smoking-associated CpGs reported previously. For 1,646 loci (in lung) and 22 (in colon), smoking was associated with both DNAm and local gene expression. For loci detected in both lung and colon (e.g., AHRR, CYP1B1, CYP1A1), top CpGs often differed between tissues, but similar clusters of hyper- or hypomethylated CpGs were observed, with hypomethylation at regulatory elements corresponding to increased expression. For lung tissue, 17 hallmark gene sets were enriched for smoking-associated CpGs, including xenobiotic- and cancer-related gene sets. At least four smoking-associated regions in lung were impacted by lung methylation quantitative trait loci (QTLs) that co-localize with genome-wide association study (GWAS) signals for lung function (FEV1/FVC), suggesting epigenetic alterations can mediate the effects of smoking on lung health. Our multi-tissue approach has identified smoking-associated regions in disease-relevant tissues, including effects that are shared across tissue types.

The omics era: a nexus of untapped potential for Mendelian chromatinopathies

The OMICs cascade describes the hierarchical flow of information through biological systems. The epigenome sits at the apex of the cascade, thereby regulating the RNA and protein expression of the human genome and governs cellular identity and function. Genes that regulate the epigenome, termed epigenes, orchestrate complex biological signaling programs that drive human development. The broad expression patterns of epigenes during human development mean that pathogenic germline mutations in epigenes can lead to clinically significant multi-system malformations, developmental delay, intellectual disabilities, and stem cell dysfunction. In this review, we refer to germline developmental disorders caused by epigene mutation as "chromatinopathies". We curated the largest number of human chromatinopathies to date and our expanded approach more than doubled the number of established chromatinopathies to 179 disorders caused by 148 epigenes. Our study revealed that 20.6% (148/720) of epigenes cause at least one chromatinopathy. In this review, we highlight key examples in which OMICs approaches have been applied to chromatinopathy patient biospecimens to identify underlying disease pathogenesis. The rapidly evolving OMICs technologies that couple molecular biology with high-throughput sequencing or proteomics allow us to dissect out the causal mechanisms driving temporal-, cellular-, and tissue-specific expression. Using the full repertoire of data generated by the OMICs cascade to study chromatinopathies will provide invaluable insight into the developmental impact of these epigenes and point toward future precision targets for these rare disorders.

Analytical challenges in omics research on asthma and allergy: A National Institute of Allergy and Infectious Diseases workshop

Studies of asthma and allergy are generating increasing volumes of omics data for analysis and interpretation. The National Institute of Allergy and Infectious Diseases (NIAID) assembled a workshop comprising investigators studying asthma and allergic diseases using omics approaches, omics investigators from outside the field, and NIAID medical and scientific officers to discuss the following areas in asthma and allergy research: genomics, epigenomics, transcriptomics, microbiomics, metabolomics, proteomics, lipidomics, integrative omics, systems biology, and causal inference. Current states of the art, present challenges, novel and emerging strategies, and priorities for progress were presented and discussed for each area. This workshop report summarizes the major points and conclusions from this NIAID workshop. As a group, the investigators underscored the imperatives for rigorous analytic frameworks, integration of different omics data types, cross-disciplinary interaction, strategies for overcoming current limitations, and the overarching goal to improve scientific understanding and care of asthma and allergic diseases.

Detailed DNA methylation characterisation of phyllodes tumours identifies a signature of malignancy and distinguishes phyllodes from metaplastic breast carcinoma

Phyllodes tumours (PTs) are rare fibroepithelial lesions of the breast that are classified as benign, borderline, or malignant. As little is known about the molecular underpinnings of PTs, current diagnosis relies on histological examination. However, accurate classification is often difficult, particularly for distinguishing borderline from malignant PTs. Furthermore, PTs can be misdiagnosed as other tumour types with shared histological features, such as fibroadenoma and metaplastic breast cancers. As DNA methylation is a recognised hallmark of many cancers, we hypothesised that DNA methylation could provide novel biomarkers for diagnosis and tumour stratification in PTs, whilst also allowing insight into the molecular aetiology of this otherwise understudied tumour. We generated whole-genome methylation data using the Illumina EPIC microarray in a novel PT cohort (n = 33) and curated methylation microarray data from published datasets including PTs and other potentially histopathologically similar tumours (total n = 817 samples). Analyses revealed that PTs have a unique methylome compared to normal breast tissue and to potentially histopathologically similar tumours (metaplastic breast cancer, fibroadenoma and sarcomas), with PT-specific methylation changes enriched in gene sets involved in KRAS signalling and epithelial-mesenchymal transition. Next, we identified 53 differentially methylated regions (DMRs) (false discovery rate < 0.05) that specifically delineated malignant from non-malignant PTs. The top DMR in both discovery and validation cohorts was hypermethylation at the HSD17B8 CpG island promoter. Matched PT single-cell expression data showed that HSD17B8 had minimal expression in fibroblast (putative tumour) cells. Finally, we created a methylation classifier to distinguish PTs from metaplastic breast cancer samples, where we revealed a likely misdiagnosis for two TCGA metaplastic breast cancer samples. In conclusion, DNA methylation alterations are associated with PT histopathology and hold the potential to improve our understanding of PT molecular aetiology, diagnostics, and risk stratification. © 2024 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

The Sociodemographic and Lifestyle Correlates of Epigenetic Aging in a Nationally Representative U.S. Study of Younger Adults

Importance

Epigenetic clocks represent molecular evidence of disease risk and aging processes and have been used to identify how social and lifestyle characteristics are associated with accelerated biological aging. However, most of this research is based on older adult samples who already have measurable chronic disease.

Objective

To investigate whether and how sociodemographic and lifestyle characteristics are related to biological aging in a younger adult sample across a wide array of epigenetic clock measures.

Design

Nationally representative prospective cohort study.

Setting

United States (U.S.).

Participants

Data come from the National Longitudinal Study of Adolescent to Adult Health, a national cohort of adolescents in grades 7-12 in U.S. in 1994 followed for 25 years over five interview waves. Our analytic sample includes participants followed-up through Wave V in 2016-18 who provided blood samples for DNA methylation (DNAm) testing (n=4237) at Wave V.

Exposure

Sociodemographic (sex, race/ethnicity, immigrant status, socioeconomic status, geographic location) and lifestyle (obesity status, exercise, tobacco, and alcohol use) characteristics.

Main Outcome

Biological aging assessed from blood DNAm using 16 epigenetic clocks when the cohort was aged 33-44 in Wave V.

Results

While there is considerable variation in the mean and distribution of epigenetic clock estimates and in the correlations among the clocks, we found sociodemographic and lifestyle factors are more often associated with biological aging in clocks trained to predict current or dynamic phenotypes (e.g., PhenoAge, GrimAge and DunedinPACE) as opposed to clocks trained to predict chronological age alone (e.g., Horvath). Consistent and strong associations of faster biological aging were found for those with lower levels of education and income, and those with severe obesity, no weekly exercise, and tobacco use.

Conclusions and Relevance

Our study found important social and lifestyle factors associated with biological aging in a nationally representative cohort of younger-aged adults. These findings indicate that molecular processes underlying disease risk can be identified in adults entering midlife before disease is manifest and represent useful targets for interventions to reduce social inequalities in heathy aging and longevity.

Key Points

Question: Are epigenetic clocks, measures of biological aging developed mainly on older-adult samples, meaningful for younger adults and associated with sociodemographic and lifestyle characteristics in expected patterns found in prior aging research?Findings: Sociodemographic and lifestyle factors were associated with biological aging in clocks trained to predict morbidity and mortality showing accelerated aging among those with lower levels of education and income, and those with severe obesity, no weekly exercise, and tobacco use.Meaning: Age-related molecular processes can be identified in younger-aged adults before disease manifests and represent potential interventions to reduce social inequalities in heathy aging and longevity.

Three-dimensional assessments are necessary to determine the true, spatially-resolved composition of tissues

Methods for spatially resolved cellular profiling using thinly cut sections have enabled in-depth quantitative tissue mapping to study inter-sample and intra-sample differences in normal human anatomy and disease onset and progression. These methods often profile extremely limited regions, which may impact the evaluation of heterogeneity due to tissue sub-sampling. Here, we applied CODA, a deep learning-based tissue mapping platform, to reconstruct the three-dimensional (3D) microanatomy of grossly normal and cancer-containing human pancreas biospecimens obtained from individuals who underwent pancreatic resection. To compare inter- and intra-sample heterogeneity, we assessed bulk and spatially resolved tissue composition in a cohort of two-dimensional (2D) whole slide images (WSIs) and a cohort of thick slabs of pancreas tissue that were digitally reconstructed in 3D from serial sections. To demonstrate the marked under sampling of 2D assessments, we simulated the number of WSIs and tissue microarrays (TMAs) necessary to represent the compositional heterogeneity of 3D data within 10% error to reveal that tens of WSIs and hundreds of TMA cores are sometimes needed. We show that spatial correlation of different pancreatic structures decay significantly within a span of microns, demonstrating that 2D histological sections may not be representative of their neighboring tissues. In sum, we demonstrate that 3D assessments are necessary to accurately assess tissue composition in normal and abnormal specimens and in order to accurately determine neoplastic content. These results emphasize the importance of intra-sample heterogeneity in tissue mapping efforts.

Epigenome-wide Association Study Shows Differential DNA Methylation of MDC1, KLF9, and CUTA in Autoimmune Thyroid Disease

CONTEXT

Autoimmune thyroid disease (AITD) includes Graves disease (GD) and Hashimoto disease (HD), which often run in the same family. AITD etiology is incompletely understood: Genetic factors may account for up to 75% of phenotypic variance, whereas epigenetic effects (including DNA methylation [DNAm]) may contribute to the remaining variance (eg, why some individuals develop GD and others HD).

OBJECTIVE

This work aimed to identify differentially methylated positions (DMPs) and differentially methylated regions (DMRs) comparing GD to HD.

METHODS

Whole-blood DNAm was measured across the genome using the Infinium MethylationEPIC array in 32 Australian patients with GD and 30 with HD (discovery cohort) and 32 Danish patients with GD and 32 with HD (replication cohort). Linear mixed models were used to test for differences in quantile-normalized β values of DNAm between GD and HD and data were later meta-analyzed. Comb-p software was used to identify DMRs.

RESULTS

We identified epigenome-wide significant differences (P < 9E-8) and replicated (P < .05) 2 DMPs between GD and HD (cg06315208 within MDC1 and cg00049440 within KLF9). We identified and replicated a DMR within CUTA (5 CpGs at 6p21.32). We also identified 64 DMPs and 137 DMRs in the meta-analysis.

CONCLUSION

Our study reveals differences in DNAm between GD and HD, which may help explain why some people develop GD and others HD and provide a link to environmental risk factors. Additional research is needed to advance understanding of the role of DNAm in AITD and investigate its prognostic and therapeutic potential.

Toward a more holistic approach to the study of exposures and child outcomes

Aim: The current work was designed to demonstrate the application of the exposome framework in examining associations between exposures and children's long-term neurodevelopmental and behavioral outcomes. Methods: Longitudinal data were collected from birth through age 6 from 402 preterm infants. Three statistical methods were utilized to demonstrate the exposome framework: exposome-wide association study, cumulative exposure and machine learning models, with and without epigenetic data. Results: Each statistical approach answered a distinct research question regarding the impact of exposures on longitudinal child outcomes. Findings highlight associations between exposures, epigenetics and executive function. Conclusion: Findings demonstrate how an exposome-based approach can be utilized to understand relationships between internal (e.g., DNA methylation) and external (e.g., prenatal risk) exposures and long-term developmental outcomes in preterm children.

Epigenetic associations in HPA axis genes related to bronchopulmonary dysplasia and antenatal steroids

BACKGROUND

Bronchopulmonary dysplasia (BPD), a common morbidity among very preterm infants, is associated with chronic disease and neurodevelopmental impairments. A hypothesized mechanism for these outcomes lies in altered glucocorticoid (GC) activity. We hypothesized that BPD and its treatments may result in epigenetic differences in the hypothalamic-pituitary-adrenal (HPA) axis, which is modulated by GC, and could be ascertained using an established GC risk score and DNA methylation (DNAm) of HPA axis genes.

METHODS

DNAm was quantified from buccal tissue (ECHO-NOVI) and from neonatal blood spots (ELGAN ECHO) via the EPIC microarray. Prenatal maternal characteristics, pregnancy complication, and neonatal medical complication data were collected from medical record review and maternal interviews.

RESULTS

The GC score was not associated with steroid exposure or BPD. However, six HPA genes involved in stress response regulation demonstrated differential methylation with antenatal steroid exposure; two CpGs within FKBP5 and POMC were differentially methylated with BPD severity. These findings were sex-specific in both cohorts; males had greater magnitude of differential methylation within these genes.

CONCLUSIONS

These findings suggest that BPD severity and antenatal steroids are associated with DNAm at some HPA genes in very preterm infants and the effects appear to be sex-, tissue-, and age-specific.

IMPACT

This study addresses bronchopulmonary dysplasia (BPD), an important health outcome among preterm neonates, and interrogates a commonly studied pathway, the hypothalamic-pituitary-adrenal (HPA) axis. The combination of BPD, the HPA axis, and epigenetic markers has not been previously reported. In this study, we found that BPD itself was not associated with epigenetic responses in the HPA axis in infants born very preterm; however, antenatal treatment with steroids was associated with epigenetic responses.

Proteogenomic analysis of enriched HGSOC tumor epithelium identifies prognostic signatures and therapeutic vulnerabilities

We performed a deep proteogenomic analysis of bulk tumor and laser microdissection enriched tumor cell populations from high-grade serous ovarian cancer (HGSOC) tissue specimens spanning a broad spectrum of purity. We identified patients with longer progression-free survival had increased immune-related signatures and validated proteins correlating with tumor-infiltrating lymphocytes in 65 tumors from an independent cohort of HGSOC patients, as well as with overall survival in an additional 126 HGSOC patient cohort. We identified that homologous recombination deficient (HRD) tumors are enriched in pathways associated with metabolism and oxidative phosphorylation that we validated in independent patient cohorts. We further identified that polycomb complex protein BMI-1 is elevated in HR proficient (HRP) tumors, that elevated BMI-1 correlates with poor overall survival in HRP but not HRD HGSOC patients, and that HRP HGSOC cells are uniquely sensitive to BMI-1 inhibition.

Epigenome-wide association study of total nicotine equivalents in multiethnic current smokers from three prospective cohorts

The impact of tobacco exposure on health varies by race and ethnicity and is closely tied to internal nicotine dose, a marker of carcinogen uptake. DNA methylation is strongly responsive to smoking status and may mediate health effects, but study of associations with internal dose is limited. We performed a blood leukocyte epigenome-wide association study (EWAS) of urinary total nicotine equivalents (TNEs; a measure of nicotine uptake) and DNA methylation measured using the MethylationEPIC v1.0 BeadChip (EPIC) in six racial and ethnic groups across three cohort studies. In the Multiethnic Cohort Study (discovery, n = 1994), TNEs were associated with differential methylation at 408 CpG sites across >250 genomic regions (p < 9 × 10-8). The top significant sites were annotated to AHRR, F2RL3, RARA, GPR15, PRSS23, and 2q37.1, all of which had decreasing methylation with increasing TNEs. We identified 45 novel CpG sites, of which 42 were unique to the EPIC array and eight annotated to genes not previously linked with smoking-related DNA methylation. The most significant signal in a novel gene was cg03748458 in MIR383;SGCZ. Fifty-one of the 408 discovery sites were validated in the Singapore Chinese Health Study (n = 340) and the Southern Community Cohort Study (n = 394) (Bonferroni corrected p < 1.23 × 10-4). Significant heterogeneity by race and ethnicity was detected for CpG sites in MYO1G and CYTH1. Furthermore, TNEs significantly mediated the association between cigarettes per day and DNA methylation at 15 sites (average 22.5%-44.3% proportion mediated). Our multiethnic study highlights the transethnic and ethnic-specific methylation associations with internal nicotine dose, a strong predictor of smoking-related morbidities.

Characterisation and reproducibility of the HumanMethylationEPIC v2.0 BeadChip for DNA methylation profiling

BACKGROUND

The Illumina family of Infinium Methylation BeadChip microarrays has been widely used over the last 15 years for genome-wide DNA methylation profiling, including large-scale and population-based studies, due to their ease of use and cost effectiveness. Succeeding the popular HumanMethylationEPIC BeadChip (EPICv1), the recently released Infinium MethylationEPIC v2.0 BeadChip (EPICv2) claims to extend genomic coverage to more than 935,000 CpG sites. Here, we comprehensively characterise the reproducibility, reliability and annotation of the EPICv2 array, based on bioinformatic analysis of both manifest data and new EPICv2 data from diverse biological samples.

RESULTS

We find a high degree of reproducibility with EPICv1, evidenced by comparable sensitivity and precision from empirical cross-platform comparison incorporating whole genome bisulphite sequencing (WGBS), and high correlation between technical sample replicates, including between samples with DNA input levels below the manufacturer's recommendation. We provide a full assessment of probe content, evaluating genomic distribution and changes from previous array versions. We characterise EPICv2's new feature of replicated probes and provide recommendations as to the superior probes. In silico analysis of probe sequences demonstrates that probe cross-hybridisation remains a significant problem in EPICv2. By mapping the off-target sites at single nucleotide resolution and comparing with WGBS we show empirical evidence for preferential off-target binding.

CONCLUSIONS

Overall, we find EPICv2 a worthy successor to the previous Infinium methylation microarrays, however some technical issues remain. To support optimal EPICv2 data analysis we provide an expanded version of the EPICv2 manifest to aid researchers in understanding probe design, data processing, choosing appropriate probes for analysis and for integration with methylation datasets from previous versions of the Infinium Methylation BeadChip.

Epigenetic variation in neonatal tissues in infants conceived using capacitation-in vitro maturation vs. in vitro fertilization

OBJECTIVE

To investigate alterations of the global DNA methylation profile in placenta, cord blood, and neonatal buccal smears in infants conceived using in vitro maturation (IVM) with a prematuration step (capacitation-IVM [CAPA-IVM]) vs. in vitro fertilization (IVF).

DESIGN

Analysis of data from the offspring of participants in a randomized controlled trial.

SETTING

Private clinic.

PATIENTS

Forty-six women with polycystic ovary syndrome and/or high antral follicle count and their offspring (58 newborns).

INTERVENTION(S)

Women with polycystic ovary syndrome and/or a high antral follicle count participating in the clinical trial were randomized to undergo CAPA-IVM or conventional IVF.

MAIN OUTCOME MEASURE(S)

At delivery, biological samples including cord blood, placental tissue, and a neonatal buccal smear were collected. Genome-wide DNA methylation was determined using the Illumina Infinium MethylationEPIC BeadChip. Variability in methylation was also considered, and mean variances for the two treatment categories were compared.

RESULTS

In neonatal buccal smears, there were no significant differences between the CAPA-IVM and conventional IVF groups on the basis of the CpG probe after linear regression analysis using a significant cut-off of false-discovery rate <0.05 and |Δβ|≥0.05. In cord blood, only one CpG site showed a significant gain of methylation in the CAPA-IVM group. In the placenta, CAPA-IVM was significantly associated with changes in methylation at five CpG sites. Significantly more variable DNA methylation was found in five probes in the placenta, 54 in cord blood, and two in buccal smears after IVM of oocytes. In cord blood samples, 20 CpG sites had more variable methylation in the conventional IVF vs. IVM group. Isolated CpG sites showing differences in methylation in cord blood were not associated with changes in gene expression of the overlapping genes.

CONCLUSION(S)

Capacitation-IVM appeared to be associated with only a small amount of epigenetic variation in cord blood, placental tissue, and neonate buccal smears.

CLINICAL TRIAL REGISTRATION NUMBER

NCT03405701 (www.

CLINICALTRIALS

gov).

The astrocyte-produced growth factor HB-EGF limits autoimmune CNS pathology

Central nervous system (CNS)-resident cells such as microglia, oligodendrocytes and astrocytes are gaining increasing attention in respect to their contribution to CNS pathologies including multiple sclerosis (MS). Several studies have demonstrated the involvement of pro-inflammatory glial subsets in the pathogenesis and propagation of inflammatory events in MS and its animal models. However, it has only recently become clear that the underlying heterogeneity of astrocytes and microglia can not only drive inflammation, but also lead to its resolution through direct and indirect mechanisms. Failure of these tissue-protective mechanisms may potentiate disease and increase the risk of conversion to progressive stages of MS, for which currently available therapies are limited. Using proteomic analyses of cerebrospinal fluid specimens from patients with MS in combination with experimental studies, we here identify Heparin-binding EGF-like growth factor (HB-EGF) as a central mediator of tissue-protective and anti-inflammatory effects important for the recovery from acute inflammatory lesions in CNS autoimmunity. Hypoxic conditions drive the rapid upregulation of HB-EGF by astrocytes during early CNS inflammation, while pro-inflammatory conditions suppress trophic HB-EGF signaling through epigenetic modifications. Finally, we demonstrate both anti-inflammatory and tissue-protective effects of HB-EGF in a broad variety of cell types in vitro and use intranasal administration of HB-EGF in acute and post-acute stages of autoimmune neuroinflammation to attenuate disease in a preclinical mouse model of MS. Altogether, we identify astrocyte-derived HB-EGF and its epigenetic regulation as a modulator of autoimmune CNS inflammation and potential therapeutic target in MS.

Aberrant PRDM2 methylation as an early event in serrated lesions destined to evolve into microsatellite-instable colorectal cancers

Up to 30% of colorectal cancers (CRCs) develop from sessile serrated lesions (SSLs). Within the serrated neoplasia pathway, at least two principally distinct oncogenetic routes exist generating microsatellite-stable and microsatellite-instable CRCs, respectively. Aberrant DNA methylation (DNAm) is found early in the serrated pathway and might play a role in both oncogenetic routes. We studied a cohort of 23 SSLs with a small focus (<10 mm) of dysplasia or cancer, 10 of which were MLH1 deficient and 13 MLH1 proficient. By comparing, for each SSL, the methylation status of (1) the region of dysplasia or cancer (SSL-D), (2) the nondysplastic SSL (SSL), and (3) adjacent normal mucosa, differentially methylated probes (DMPs) and regions (DMRs) were assessed both genome-wide as well as in a tumor-suppressor gene-focused approach. By comparing DNAm of MLH1-deficient SSL-Ds with their corresponding SSLs, we identified five DMRs, including those annotating for PRDM2 and, not unexpectedly, MLH1. PRDM2 gene promotor methylation was associated with MLH1 expression status, as it was largely hypermethylated in MLH1-deficient SSL-Ds and hypomethylated in MLH1-proficient SSL-Ds. Significantly increased DNAm levels of PRDM2 and MLH1, in particular at 'critical' MLH1 probe sites, were to some extent already visible in SSLs as compared to normal mucosa (p = 0.02, p = 0.01, p < 0.0001, respectively). No DMRs, nor DMPs, were identified for SSLs destined to evolve into MLH1-proficient SSL-Ds. Our data indicate that, within both arms of the serrated CRC pathway, the majority of the epigenetic alterations are introduced early during SSL formation. Promoter hypermethylation of PRDM2 and MLH1 on the other hand specifically initiates in SSLs destined to transform into MLH1-deficient CRCs suggesting that the fate of SSLs may not necessarily result from a stochastic process but possibly is already imprinted and predisposed.

DNA methylation profiling reveals potential biomarkers of β-lactams induced fatal anaphylactic shock

Anaphylaxis is a serious reaction of systemic hypersensitivity with that rapid onset and sudden death. Drug hypersensitivity, particularly induced by β-lactams, is one of the most frequent causes of anaphylaxis in adults. But identification of anaphylactic shock, in forensic sciences recently, is difficult, because it mainly depends on nonspecific characteristic morphological changes, as well as exclusion and circumstantial evidence. Here, we detected DNA methylation signatures of β-lactams-induced fatal anaphylactic shock with the Illumina Infinium Human Methylation EPIC BeadChip, to screen potential forensic biomarkers and reveal the molecular mechanisms of drug-induced anaphylaxis with fatal shock and sudden death. Our results indicated that DNA methylation was associated with β-lactams-induced fatal anaphylactic shock, in which the hypomethylation played a vital role. We found that 1459 differentially methylated positions (DMPs) were mainly involved in β-lactams-induced fatal anaphylactic shock by regulating MAPK and other signaling pathways. 18 DNA methylation signatures that could separate β-lactams-induced anaphylactic shock from healthy individuals were identified. The altered methylation of DMPs can affect the transcription of corresponding genes and promote β-lactams-induced fatal anaphylactic shock. The results suggest that DNA methylation can detect forensic identification markers of drug-induced anaphylaxis with fatal shock and sudden death, and it is an effective method for the forensic diagnosis.

The potential of epigenetic therapy to target the 3D epigenome in endocrine-resistant breast cancer

Three-dimensional (3D) epigenome remodeling is an important mechanism of gene deregulation in cancer. However, its potential as a target to counteract therapy resistance remains largely unaddressed. Here, we show that epigenetic therapy with decitabine (5-Aza-mC) suppresses tumor growth in xenograft models of pre-clinical metastatic estrogen receptor positive (ER+) breast tumor. Decitabine-induced genome-wide DNA hypomethylation results in large-scale 3D epigenome deregulation, including de-compaction of higher-order chromatin structure and loss of boundary insulation of topologically associated domains. Significant DNA hypomethylation associates with ectopic activation of ER-enhancers, gain in ER binding, creation of new 3D enhancer-promoter interactions and concordant up-regulation of ER-mediated transcription pathways. Importantly, long-term withdrawal of epigenetic therapy partially restores methylation at ER-enhancer elements, resulting in a loss of ectopic 3D enhancer-promoter interactions and associated gene repression. Our study illustrates the potential of epigenetic therapy to target ER+ endocrine-resistant breast cancer by DNA methylation-dependent rewiring of 3D chromatin interactions, which are associated with the suppression of tumor growth.

[Explainable artificial intelligence in pathology]

With the advancements in precision medicine, the demands on pathological diagnostics have increased, requiring standardized, quantitative, and integrated assessments of histomorphological and molecular pathological data. Great hopes are placed in artificial intelligence (AI) methods, which have demonstrated the ability to analyze complex clinical, histological, and molecular data for disease classification, biomarker quantification, and prognosis estimation. This paper provides an overview of the latest developments in pathology AI, discusses the limitations, particularly concerning the black box character of AI, and describes solutions to make decision processes more transparent using methods of so-called explainable AI (XAI).

Genome-wide DNA methylation analysis of cannabis use disorder in a veteran cohort enriched for posttraumatic stress disorder

Cannabis use has been increasing over the past decade, not only in the general US population, but particularly among military veterans. With this rise in use has come a concomitant increase in cannabis use disorder (CUD) among veterans. Here, we performed an epigenome-wide association study for lifetime CUD in an Iraq/Afghanistan era veteran cohort enriched for posttraumatic stress disorder (PTSD) comprising 2,310 total subjects (1,109 non-Hispanic black and 1,201 non-Hispanic white). We also investigated CUD interactions with current PTSD status and examined potential indirect effects of DNA methylation (DNAm) on the relationship between CUD and psychiatric diagnoses. Four CpGs were associated with lifetime CUD, even after controlling for the effects of current smoking (AHRR cg05575921, LINC00299 cg23079012, VWA7 cg22112841, and FAM70A cg08760398). Importantly, cg05575921, a CpG strongly linked to smoking, remained associated with lifetime CUD even when restricting the analysis to veterans who reported never smoking cigarettes. Moreover, CUD interacted with current PTSD to affect cg05575921 and cg23079012 such that those with both CUD and PTSD displayed significantly lower DNAm compared to the other groups. Finally, we provide preliminary evidence that AHRR cg05575921 helps explain the association between CUD and any psychiatric diagnoses, specifically mood disorders.

Former smoking associated with epigenetic modifications in human granulosa cells among women undergoing assisted reproduction

Smoking exposure during adulthood can disrupt oocyte development in women, contributing to infertility and possibly adverse birth outcomes. Some of these effects may be reflected in epigenome profiles in granulosa cells (GCs) in human follicular fluid. We compared the epigenetic modifications throughout the genome in GCs from women who were former (N = 15) versus never smokers (N = 44) undergoing assisted reproductive technologies (ART). This study included 59 women undergoing ART. Smoking history including time since quitting was determined by questionnaire. GCs were collected during oocyte retrieval and DNA methylation (DNAm) levels were profiled using the Infinium MethylationEPIC BeadChip. We performed an epigenome-wide association study with robust linear models, regressing DNAm level at individual loci on smoking status, adjusting for age, ovarian stimulation protocol, and three surrogate variables. We performed differentially methylated regions (DMRs) analysis and over-representation analysis of the identified CpGs and corresponding gene set. 81 CpGs were differentially methylated among former smokers compared to never smokers (FDR < 0.05). We identified 2 significant DMRs (KCNQ1 and RHBDD2). The former smoking-associated genes were enriched in oxytocin signaling, adrenergic signaling in cardiomyocytes, platelet activation, axon guidance, and chemokine signaling pathway. These epigenetic variations have been associated with inflammatory responses, reproductive outcomes, cancer development, neurodevelopmental disorder, and cardiometabolic health. Secondarily, we examined the relationships between time since quitting and DNAm at significant CpGs. We observed three CpGs in negative associations with the length of quitting smoking (p < 0.05), which were cg04254052 (KCNIP1), cg22875371 (OGDHL), and cg27289628 (LOC148145), while one in positive association, which was cg13487862 (PLXNB1). As a pilot study, we demonstrated epigenetic modifications associated with former smoking in GCs. The study is informative to potential biological pathways underlying the documented association between smoking and female infertility and biomarker discovery for smoking-associated reproductive outcomes.

Epigenomic profiling of the infrapatellar fat pad in osteoarthritis

Osteoarthritis is a prevalent, complex disease of the joints, and affects multiple intra-articular tissues. Here, we have examined genome-wide DNA methylation profiles of primary infrapatellar fat pad and matched blood samples from 70 osteoarthritis patients undergoing total knee replacement surgery. Comparing the DNA methylation profiles between these tissues reveal widespread epigenetic differences. We produce the first genome-wide methylation quantitative trait locus (mQTL) map of fat pad, and make the resource available to the wider community. Using two-sample Mendelian randomization and colocalization analyses, we resolve osteoarthritis GWAS signals and provide insights into the molecular mechanisms underpinning disease aetiopathology. Our findings provide the first view of the epigenetic landscape of infrapatellar fat pad primary tissue in osteoarthritis.

Epigenome-wide association study identifies neonatal DNA methylation associated with two-year attention problems in children born very preterm

Prior research has identified epigenetic predictors of attention problems in school-aged children but has not yet investigated these in young children, or children at elevated risk of attention problems due to preterm birth. The current study evaluated epigenome-wide associations between neonatal DNA methylation and attention problems at age 2 years in children born very preterm. Participants included 441 children from the Neonatal Neurobehavior and Outcomes in Very Preterm Infants (NOVI) Study, a multi-site study of infants born < 30 weeks gestational age. DNA methylation was measured from buccal swabs collected at NICU discharge using the Illumina MethylationEPIC Bead Array. Attention problems were assessed at 2 years of adjusted age using the attention problems subscale of the Child Behavior Checklist (CBCL). After adjustment for multiple testing, DNA methylation at 33 CpG sites was associated with child attention problems. Differentially methylated CpG sites were located in genes previously linked to physical and mental health, including several genes associated with ADHD in prior epigenome-wide and genome-wide association studies. Several CpG sites were located in genes previously linked to exposure to prenatal risk factors in the NOVI sample. Neonatal epigenetics measured at NICU discharge could be useful in identifying preterm children at risk for long-term attention problems and related psychiatric disorders, who could benefit from early prevention and intervention efforts.

Reconstructing DNA methylation maps of ancient populations

Studying premortem DNA methylation from ancient DNA (aDNA) provides a proxy for ancient gene activity patterns, and hence valuable information on evolutionary changes in gene regulation. Due to statistical limitations, current methods to reconstruct aDNA methylation maps are constrained to high-coverage shotgun samples, which comprise a small minority of available ancient samples. Most samples are sequenced using in-situ hybridization capture sequencing which targets a predefined set of genomic positions. Here, we develop methods to reconstruct aDNA methylation maps of samples that were not sequenced using high-coverage shotgun sequencing, by way of pooling together individuals to obtain a DNA methylation map that is characteristic of a population. We show that the resulting DNA methylation maps capture meaningful biological information and allow for the detection of differential methylation across populations. We offer guidelines on how to carry out comparative studies involving ancient populations, and how to control the rate of falsely discovered differentially methylated regions. The ability to reconstruct DNA methylation maps of past populations allows for the development of a whole new frontier in paleoepigenetic research, tracing DNA methylation changes throughout human history, using data from thousands of ancient samples.

Comparing feature selection and machine learning approaches for predicting CYP2D6 methylation from genetic variation

Introduction

Pharmacogenetics currently supports clinical decision-making on the basis of a limited number of variants in a few genes and may benefit paediatric prescribing where there is a need for more precise dosing. Integrating genomic information such as methylation into pharmacogenetic models holds the potential to improve their accuracy and consequently prescribing decisions. Cytochrome P450 2D6 (CYP2D6) is a highly polymorphic gene conventionally associated with the metabolism of commonly used drugs and endogenous substrates. We thus sought to predict epigenetic loci from single nucleotide polymorphisms (SNPs) related to CYP2D6 in children from the GUSTO cohort.

Methods

Buffy coat DNA methylation was quantified using the Illumina Infinium Methylation EPIC beadchip. CpG sites associated with CYP2D6 were used as outcome variables in Linear Regression, Elastic Net and XGBoost models. We compared feature selection of SNPs from GWAS mQTLs, GTEx eQTLs and SNPs within 2 MB of the CYP2D6 gene and the impact of adding demographic data. The samples were split into training (75%) sets and test (25%) sets for validation. In Elastic Net model and XGBoost models, optimal hyperparameter search was done using 10-fold cross validation. Root Mean Square Error and R-squared values were obtained to investigate each models' performance. When GWAS was performed to determine SNPs associated with CpG sites, a total of 15 SNPs were identified where several SNPs appeared to influence multiple CpG sites.

Results

Overall, Elastic Net models of genetic features appeared to perform marginally better than heritability estimates and substantially better than Linear Regression and XGBoost models. The addition of nongenetic features appeared to improve performance for some but not all feature sets and probes. The best feature set and Machine Learning (ML) approach differed substantially between CpG sites and a number of top variables were identified for each model.

Discussion

The development of SNP-based prediction models for CYP2D6 CpG methylation in Singaporean children of varying ethnicities in this study has clinical application. With further validation, they may add to the set of tools available to improve precision medicine and pharmacogenetics-based dosing.