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

Usual dietary intake and change in DNA methylation over years: EWAS in KORA FF4 and KORA fit

Introduction

Changes in DNA methylation can increase or suppress the expression of health-relevant genes. We investigated for the first time the relationship between habitual food consumption and changes in DNA methylation.

Methods

The German KORA FF4 and KORA Fit studies were used to study the change in methylation over a median follow-up of 4 years. Only subjects participating in both surveys and with available dietary and methylation data were included in the analysis (n = 465). DNA methylation was measured using the Infinium MethylationEPIC BeadChip (Illumina), resulting in 735,527 shared CpGs across both studies. Generalized estimating equation models with an interaction term of exposure and time point were used to analyze the association of 34 food groups, folic acid, and two dietary patterns with changes in DNA methylation over time.

Results

The results were corrected for genomic inflation. Significant interaction terms indicate different effects between both time points. We observed only a few significant associations between food intake and change in DNA methylation, except for cream and spirit consumption. The annotated genes include CLN3, PROM1, DLEU7, TLL2, and UGT1A10.

Discussion

We identified weak associations between food consumption and DNA methylation change. The differential results for cream and spirits, both consumed in low quantities, require replication in independent studies.

Alterations in leukocyte DNA methylome are associated to immunosuppression in severe clinical phenotypes of septic patients

Introduction

Sepsis patients experience a complex interplay of host pro- and anti-inflammatory processes which compromise the clinical outcome. Despite considering the latest clinical and scientific research, our comprehension of the immunosuppressive events in septic episodes remains incomplete. Additionally, a lack of data exists regarding the role of epigenetics in modulating immunosuppression, subsequently impacting patient survival.

Methods

To advance the current understanding of the mechanisms underlying immunosuppression, in this study we explored the dynamics of DNA methylation using the Infinium Methylation EPIC v1.0 BeadChip Kit in leukocytes from patients suffering from sepsis, septic shock, and critically ill patients as controls, within the first 24 h after admission in the Intensive Care Unit of a tertiary hospital.

Results and discussion

Employing two distinct analysis approaches (DMRcate and mCSEA) in comparing septic shock and critically ill patients, we identified 1,256 differentially methylated regions (DMRs) intricately linked to critical immune system pathways. The examination of the top 100 differentially methylated positions (DMPs) between septic shock and critically ill patients facilitated a clear demarcation among the three patient groups. Notably, the top 6,657 DMPs exhibited associations with organ dysfunction and lactate levels. Among the individual genes displaying significant differential methylation, IL10, TREM1, IL1B, and TNFAIP8 emerged with the most pronounced methylation alterations across the diverse patient groups when subjected to DNA bisulfite pyrosequencing analysis. These findings underscore the dynamic nature of DNA methylation profiles, highlighting the most pronounced alterations in patients with septic shock, and revealing their close association with the disease.

Identifying epigenetic aging moderators using the epigenetic pacemaker

Epigenetic clocks are DNA methylation-based chronological age prediction models that are commonly employed to study age-related biology. The difference between the predicted and observed age is often interpreted as a form of biological age acceleration, and many studies have measured the impact of environmental and disease-associated factors on epigenetic age. Most epigenetic clocks are fit using approaches that minimize the error between the predicted and observed chronological age, and as a result, they may not accurately model the impact of factors that moderate the relationship between the actual and epigenetic age. Here, we compare epigenetic clocks that are constructed using penalized regression methods to an evolutionary framework of epigenetic aging with the epigenetic pacemaker (EPM), which directly models DNA methylation as a function of a time-dependent epigenetic state. In simulations, we show that the value of the epigenetic state is impacted by factors such as age, sex, and cell-type composition. Next, in a dataset aggregated from previous studies, we show that the epigenetic state is also moderated by sex and the cell type. Finally, we demonstrate that the epigenetic state is also moderated by toxins in a study on polybrominated biphenyl exposure. Thus, we find that the pacemaker provides a robust framework for the study of factors that impact epigenetic age acceleration and that the effect of these factors may be obscured in traditional clocks based on linear regression models.

Multi-omic and functional analysis for classification and treatment of sarcomas with FUS-TFCP2 or EWSR1-TFCP2 fusions

Linking clinical multi-omics with mechanistic studies may improve the understanding of rare cancers. We leverage two precision oncology programs to investigate rhabdomyosarcoma with FUS/EWSR1-TFCP2 fusions, an orphan malignancy without effective therapies. All tumors exhibit outlier ALK expression, partly accompanied by intragenic deletions and aberrant splicing resulting in ALK variants that are oncogenic and sensitive to ALK inhibitors. Additionally, recurrent CKDN2A/MTAP co-deletions provide a rationale for PRMT5-targeted therapies. Functional studies show that FUS-TFCP2 blocks myogenic differentiation, induces transcription of ALK and truncated TERT, and inhibits DNA repair. Unlike other fusion-driven sarcomas, TFCP2-rearranged tumors exhibit genomic instability and signs of defective homologous recombination. DNA methylation profiling demonstrates a close relationship with undifferentiated sarcomas. In two patients, sarcoma was preceded by benign lesions carrying FUS-TFCP2, indicating stepwise sarcomagenesis. This study illustrates the potential of linking precision oncology with preclinical research to gain insight into the classification, pathogenesis, and therapeutic vulnerabilities of rare cancers.

Using saliva epigenetic data to develop and validate a multivariable predictor of esophageal cancer status

Background: Salivary epigenetic biomarkers may detect esophageal cancer. Methods: A total of 256 saliva samples from esophageal adenocarcinoma patients and matched volunteers were analyzed with Illumina EPIC methylation arrays. Three datasets were created, using 64% for discovery, 16% for testing and 20% for validation. Modules of gene-based methylation probes were created using weighted gene coexpression network analysis. Module significance to disease and gene importance to module were determined and a random forest classifier generated using best-scoring gene-related epigenetic probes. A cost-sensitive wrapper algorithm maximized cancer diagnosis. Results: Using age, sex and seven probes, esophageal adenocarcinoma was detected with area under the curve of 0.72 in discovery, 0.73 in testing and 0.75 in validation datasets. Cancer sensitivity was 88% with specificity of 31%. Conclusion: We have demonstrated a potentially clinically viable classifier of esophageal cancer based on saliva methylation.

Association of psychological distress and DNA methylation: A 5-year longitudinal population-based twin study

AIM

To identify the psychological distress (PD)-associated 5'-cytosine-phosphate-guanine-3' sites (CpGs), and investigate the temporal relationship between dynamic changes in DNA methylation (DNAm) and PD.

METHODS

This study included 1084 twins from the Chinese National Twin Register (CNTR). The CNTR conducted epidemiological investigations and blood withdrawal twice in 2013 and 2018. These included twins were used to perform epigenome-wide association studies (EWASs) and to validate the previously reported PD-associated CpGs selected from previous EWASs in PubMed, Embase, and the EWAS catalog. Next, a cross-lagged study was performed to examine the temporality between changes in DNAm and PD in 308 twins who completed both 2013 and 2018 surveys.

RESULTS

The EWAS analysis of our study identified 25 CpGs. In the validation analysis, 741 CpGs from 29 previous EWASs on PD were selected for validation, and 101 CpGs were validated to be significant at a false discovery rate <0.05. The cross-lagged analysis found a unidirectional path from PD to DNAm at 14 CpGs, while no sites showed significance from DNAm to PD.

CONCLUSIONS

This study identified and validated PD-related CpGs in a Chinese twin population, and suggested that PD may be the cause of changes in DNAm over time. The findings provide new insights into the molecular mechanisms underlying PD pathophysiology.

Characterization of HOXB13 expression patterns in localized and metastatic castration-resistant prostate cancer

HOXB13 is a key lineage homeobox transcription factor that plays a critical role in the differentiation of the prostate gland. Several studies have suggested that HOXB13 alterations may be involved in prostate cancer development and progression. Despite its potential biological relevance, little is known about the expression of HOXB13 across the disease spectrum of prostate cancer. To this end, we validated a HOXB13 antibody using genetic controls and investigated HOXB13 protein expression in murine and human developing prostates, localized prostate cancers, and metastatic castration-resistant prostate cancers. We observed that HOXB13 expression increases during later stages of murine prostate development. All localized prostate cancers showed HOXB13 protein expression. Interestingly, lower HOXB13 expression levels were observed in higher-grade tumors, although no significant association between HOXB13 expression and recurrence or disease-specific survival was found. In advanced metastatic prostate cancers, HOXB13 expression was retained in the majority of tumors. While we observed lower levels of HOXB13 protein and mRNA levels in tumors with evidence of lineage plasticity, 84% of androgen receptor-negative castration-resistant prostate cancers and neuroendocrine prostate cancers (NEPCs) retained detectable levels of HOXB13. Notably, the reduced expression observed in NEPCs was associated with a gain of HOXB13 gene body CpG methylation. In comparison to the commonly used prostate lineage marker NKX3.1, HOXB13 showed greater sensitivity in detecting advanced metastatic prostate cancers. Additionally, in a cohort of 837 patients, 383 with prostatic and 454 with non-prostatic tumors, we found that HOXB13 immunohistochemistry had a 97% sensitivity and 99% specificity for prostatic origin. Taken together, our studies provide valuable insight into the expression pattern of HOXB13 during prostate development and cancer progression. Furthermore, our findings support the utility of HOXB13 as a diagnostic biomarker for prostate cancer, particularly to confirm the prostatic origin of advanced metastatic castration-resistant tumors. © 2023 The Pathological Society of Great Britain and Ireland.

Serum DNA methylome of the colorectal cancer serrated pathway enables non-invasive detection

The clinical relevance of the colorectal cancer serrated pathway is evident, but the screening of serrated lesions remains challenging. We aimed to characterize the serum methylome of the serrated pathway and to evaluate circulating cell-free DNA (cfDNA) methylomes as a potential source of biomarkers for the non-invasive detection of serrated lesions. We collected serum samples from individuals with serrated adenocarcinoma (SAC), traditional serrated adenomas, sessile serrated lesions, hyperplastic polyps and individuals with no colorectal findings. First, we quantified cfDNA methylation with the MethylationEPIC array. Then, we compared the methylation profiles with tissue and serum datasets. Finally, we evaluated the utility of serum cfDNA methylation biomarkers. We identified a differential methylation profile able to distinguish high-risk serrated lesions from no serrated neoplasia, showing concordance with tissue methylation from SAC and sessile serrated lesions. Serum methylation profiles are pathway-specific, clearly separating serrated lesions from conventional adenomas. The combination of ninjurin 2 (NINJ2) and glutamate-rich 1 (ERICH1) methylation discriminated high-risk serrated lesions and SAC with 91.4% sensitivity (64.4% specificity), while zinc finger protein 718 (ZNF718) methylation reported 100% sensitivity for the detection of SAC (96% specificity). This is the first study exploring the serum methylome of serrated lesions. Differential methylation of cfDNA can be used for the non-invasive detection of colorectal serrated lesions.

Blood-based DNA methylation in advanced Nasopharyngeal Carcinoma exhibited distinct CpG methylation signature

The TNM staging system is currently used to detect cancer stages. Regardless, a small proportion of cancer patients recur even after therapy, suggesting more specific molecular tools are required to justify the stage-specific detection and prompt cancer diagnosis. Thus, we aimed to explore the blood-based DNA methylation signature of metastatic nasopharyngeal carcinoma (NPC) to establish a holistic methylation biomarker panel. For the identification of methylation signature, the EPIC BeadChip-based array was performed. Comparative analysis for identifying unique probes, validation, and functional studies was investigated by analyzing GEO and TCGA datasets. We observed 4093 differentially methylated probes (DMPs), 1232 hydroxymethylated probes, and 25 CpG islands. Gene expression study revealed both upregulated and downregulated genes. Correlation analysis suggested a positive (with a positive r, p ≤ 0.05) and negative (with a negative r, p ≤ 0.05) association with different cancers. TFBS analysis exhibited the binding site for many TFs. Furthermore, gene enrichment analysis indicated the involvement of those identified genes in biological pathways. However, blood-based DNA methylation data uncovered a distinct DNA methylation pattern, which might have an additive role in NPC progression by altering the TFs binding. Moreover, based on tissue-specificity, a variation of correlation between methylation and gene expression was noted in different cancers.

Methylomic, Proteomic, and Metabolomic Correlates of Traffic-Related Air Pollution in the Context of Cardiorespiratory Health: A Systematic Review, Pathway Analysis, and Network Analysis

A growing body of literature has attempted to characterize how traffic-related air pollution (TRAP) affects molecular and subclinical biological processes in ways that could lead to cardiorespiratory disease. To provide a streamlined synthesis of what is known about the multiple mechanisms through which TRAP could lead to cardiorespiratory pathology, we conducted a systematic review of the epidemiological literature relating TRAP exposure to methylomic, proteomic, and metabolomic biomarkers in adult populations. Using the 139 papers that met our inclusion criteria, we identified the omic biomarkers significantly associated with short- or long-term TRAP and used these biomarkers to conduct pathway and network analyses. We considered the evidence for TRAP-related associations with biological pathways involving lipid metabolism, cellular energy production, amino acid metabolism, inflammation and immunity, coagulation, endothelial function, and oxidative stress. Our analysis suggests that an integrated multi-omics approach may provide critical new insights into the ways TRAP could lead to adverse clinical outcomes. We advocate for efforts to build a more unified approach for characterizing the dynamic and complex biological processes linking TRAP exposure and subclinical and clinical disease and highlight contemporary challenges and opportunities associated with such efforts.

World Trade Center Exposure, DNA Methylation Changes, and Cancer: A Review of Current Evidence

Introduction: Known carcinogens in the dust and fumes from the destruction of the World Trade Center (WTC) towers on 9 November 2001 included metals, asbestos, and organic pollutants, which have been shown to modify epigenetic status. Epigenome-wide association analyses (EWAS) using uniform (Illumina) methodology have identified novel epigenetic profiles of WTC exposure. Methods: We reviewed all published data, comparing differentially methylated gene profiles identified in the prior EWAS studies of WTC exposure. This included DNA methylation changes in blood-derived DNA from cases of cancer-free "Survivors" and those with breast cancer, as well as tissue-derived DNA from "Responders" with prostate cancer. Emerging molecular pathways related to the observed DNA methylation changes in WTC-exposed groups were explored and summarized. Results: WTC dust exposure appears to be associated with DNA methylation changes across the genome. Notably, WTC dust exposure appears to be associated with increased global DNA methylation; direct dysregulation of cancer genes and pathways, including inflammation and immune system dysregulation; and endocrine system disruption, as well as disruption of cholesterol homeostasis and lipid metabolism. Conclusion: WTC dust exposure appears to be associated with biologically meaningful DNA methylation changes, with implications for carcinogenesis and development of other chronic diseases.

Retroelement-Age Clocks: Epigenetic Age Captured by Human Endogenous Retrovirus and LINE-1 DNA methylation states

Human endogenous retroviruses (HERVs), the remnants of ancient viral infections embedded within the human genome, and long interspersed nuclear elements 1 (LINE-1), a class of autonomous retrotransposons, are silenced by host epigenetic mechanisms including DNA methylation. The resurrection of particular retroelements has been linked to biological aging. Whether the DNA methylation states of locus specific HERVs and LINEs can be used as a biomarker of chronological age in humans remains unclear. We show that highly predictive epigenetic clocks of chronological age can be constructed from retroelement DNA methylation states in the immune system, across human tissues, and pan-mammalian species. We found retroelement epigenetic clocks were reversed during transient epigenetic reprogramming, accelerated in people living with HIV-1, responsive to antiretroviral therapy, and accurate in estimating long-term culture ages of human brain organoids. Our findings support the hypothesis of epigenetic dysregulation of retroelements as a potential contributor to the biological hallmarks of aging.

Epigenetics of the far northern Yakutian population

BACKGROUND

Yakuts are one of the indigenous populations of the subarctic and arctic territories of Siberia characterized by a continental subarctic climate with severe winters, with the regular January average temperature in the regional capital city of Yakutsk dipping below - 40 °C. The epigenetic mechanisms of adaptation to such ecologies and environments and, in particular, epigenetic age acceleration in the local population have not been studied before.

RESULTS

This work reports the first epigenetic study of the Yakutian population using whole-blood DNA methylation data, supplemented with the comparison to the residents of Central Russia. Gene set enrichment analysis revealed, among others, geographic region-specific differentially methylated regions associated with adaptation to climatic conditions (water consumption, digestive system regulation), aging processes (actin filament activity, cell fate), and both of them (channel activity, regulation of steroid and corticosteroid hormone secretion). Further, it is demonstrated that the epigenetic age acceleration of the Yakutian representatives is significantly higher than that of Central Russia counterparts. For both geographic regions, we showed that epigenetically males age faster than females, whereas no significant sex differences were found between the regions.

CONCLUSIONS

We performed the first study of the epigenetic data of the Yakutia cohort, paying special attention to region-specific features, aging processes, age acceleration, and sex specificity.

Integration of datasets for individual prediction of DNA methylation-based biomarkers

BACKGROUND

Epigenetic scores (EpiScores) can provide biomarkers of lifestyle and disease risk. Projecting new datasets onto a reference panel is challenging due to separation of technical and biological variation with array data. Normalisation can standardise data distributions but may also remove population-level biological variation.

RESULTS

We compare two birth cohorts (Lothian Birth Cohorts of 1921 and 1936 - nLBC1921 = 387 and nLBC1936 = 498) with blood-based DNA methylation assessed at the same chronological age (79 years) and processed in the same lab but in different years and experimental batches. We examine the effect of 16 normalisation methods on a novel BMI EpiScore (trained in an external cohort, n = 18,413), and Horvath's pan-tissue DNA methylation age, when the cohorts are normalised separately and together. The BMI EpiScore explains a maximum variance of R2=24.5% in BMI in LBC1936 (SWAN normalisation). Although there are cross-cohort R2 differences, the normalisation method makes a minimal difference to within-cohort estimates. Conversely, a range of absolute differences are seen for individual-level EpiScore estimates for BMI and age when cohorts are normalised separately versus together. While within-array methods result in identical EpiScores whether a cohort is normalised on its own or together with the second dataset, a range of differences is observed for between-array methods.

CONCLUSIONS

Normalisation methods returning similar EpiScores, whether cohorts are analysed separately or together, will minimise technical variation when projecting new data onto a reference panel. These methods are important for cases where raw data is unavailable and joint normalisation of cohorts is computationally expensive.

Waterpipe and cigarette epigenome analysis reveals markers implicated in addiction and smoking type inference

Waterpipe smoking is frequent in the Middle East and Africa with emerging prevalence worldwide. The epigenome acts as a molecular sensor to exposures and a crucial driver in several diseases. With the widespread use of waterpipe smoking, it is timely to investigate its epigenomic markers and their role in addiction, as a central player in disease prevention and therapeutic strategies. DNA methylome-wide profiling was performed on an exposure-rich population from the Middle East, constituting of 216 blood samples split equally between never, cigarette-only and waterpipe-only smokers. Waterpipe smokers showed predominantly distinct epigenetic markers from cigarette smokers, even though both smoking forms are tobacco-based. Moreover, each smoking form could be accurately (∼90 %) inferred from the DNA methylome using machine learning. Top markers showed dose-response relationship with extent of smoking and were validated using independent technologies and additional samples (total N = 284). Smoking markers were enriched in regulatory regions and several biological pathways, primarily addiction. The epigenetically altered genes were not associated with genetic etiology of tobacco use, and the methylation levels of addiction genes, in particular, were more likely to reverse after smoking cessation. In contrast, other epigenetic markers continued to feature smoking exposure after cessation, which may explain long-term health effects observed in former smokers. This study reports, for the first time, blood epigenome-wide markers of waterpipe smokers and reveals new markers of cigarette smoking, with implications in mechanisms of addiction and the capacity to discriminate between different smoking types. These markers may offer actionable targets to reverse the epigenetic memory of addiction and can guide future prevention strategies for tobacco smoking as the most preventable cause of illnesses worldwide.

Linking Prenatal Environmental Exposures to Lifetime Health with Epigenome-Wide Association Studies: State-of-the-Science Review and Future Recommendations

BACKGROUND

The prenatal environment influences lifetime health; epigenetic mechanisms likely predominate. In 2016, the first international consortium paper on cigarette smoking during pregnancy and offspring DNA methylation identified extensive, reproducible exposure signals. This finding raised expectations for epigenome-wide association studies (EWAS) of other exposures.

OBJECTIVE

We review the current state-of-the-science for DNA methylation associations across prenatal exposures in humans and provide future recommendations.

METHODS

We reviewed 134 prenatal environmental EWAS of DNA methylation in newborns, focusing on 51 epidemiological studies with meta-analysis or replication testing. Exposures spanned cigarette smoking, alcohol consumption, air pollution, dietary factors, psychosocial stress, metals, other chemicals, and other exogenous factors. Of the reproducible DNA methylation signatures, we examined implementation as exposure biomarkers.

RESULTS

Only 19 (14%) of these prenatal EWAS were conducted in cohorts of 1,000 or more individuals, reflecting the still early stage of the field. To date, the largest perinatal EWAS sample size was 6,685 participants. For comparison, the most recent genome-wide association study for birth weight included more than 300,000 individuals. Replication, at some level, was successful with exposures to cigarette smoking, folate, dietary glycemic index, particulate matter with aerodynamic diameter < 10 μ m and < 2.5 μ m , nitrogen dioxide, mercury, cadmium, arsenic, electronic waste, PFAS, and DDT. Reproducible effects of a more limited set of prenatal exposures (smoking, folate) enabled robust methylation biomarker creation.

DISCUSSION

Current evidence demonstrates the scientific premise for reproducible DNA methylation exposure signatures. Better powered EWAS could identify signatures across many exposures and enable comprehensive biomarker development. Whether methylation biomarkers of exposures themselves cause health effects remains unclear. We expect that larger EWAS with enhanced coverage of epigenome and exposome, along with improved single-cell technologies and evolving methods for integrative multi-omics analyses and causal inference, will expand mechanistic understanding of causal links between environmental exposures, the epigenome, and health outcomes throughout the life course. https://doi.org/10.1289/EHP12956.

Biological stability of DNA methylation measurements over varying intervals of time and in the presence of acute stress

Identifying factors that influence the stability of DNA methylation measurements across biological replicates is of critical importance in basic and clinical research. Using a within-person between-group experimental design (n = 31, number of observations = 192), we report the stability of biological replicates over a variety of unique temporal scenarios, both in the absence and presence of acute psychosocial stress, and between individuals who have experienced early life adversity (ELA) and non-exposed individuals. We found that varying time intervals, acute stress, and ELA exposure influenced the stability of repeated DNA methylation measurements. In the absence of acute stress, probes were less stable as time passed; however, stress exerted a stabilizing influence on probes over longer time intervals. Compared to non-exposed individuals, ELA-exposed individuals had significantly lower probe stability immediately following acute stress. Additionally, we found that across all scenarios, probes used in most epigenetic-based algorithms for estimating epigenetic age or immune cell proportions had average or below-average stability, except for the Principal Component and DunedinPACE epigenetic ageing clocks, which were enriched for more stable probes. Finally, using highly stable probes in the absence of stress, we identified multiple probes that were hypomethylated in the presence of acute stress, regardless of ELA status. Two hypomethylated probes are located near the transcription start site of the glutathione-disulfide reductase gene (GSR), which has previously been shown to be an integral part of the stress response to environmental toxins. We discuss implications for future studies concerning the reliability and reproducibility of DNA methylation measurements.Abbreviations: DNAm - DNA methylation, CpG - 5'-cytosine-phosphate-guanine-3,' ICC - Interclass correlation coefficient, ELA - Early-life adversity, PBMCs - Peripheral blood mononuclear cells, mQTL - Methylation quantitative trait loci, TSS - Transcription start site, GSR - Glutathione-disulfide reductase gene, TSST - Trier social stress test, PC - Principal component.

A novel principal component based method for identifying differentially methylated regions in Illumina Infinium MethylationEPIC BeadChip data

Differentially methylated regions (DMRs) are genomic regions with methylation patterns across multiple CpG sites that are associated with a phenotype. In this study, we proposed a Principal Component (PC) based DMR analysis method for use with data generated using the Illumina Infinium MethylationEPIC BeadChip (EPIC) array. We obtained methylation residuals by regressing the M-values of CpGs within a region on covariates, extracted PCs of the residuals, and then combined association information across PCs to obtain regional significance. Simulation-based genome-wide false positive (GFP) rates and true positive rates were estimated under a variety of conditions before determining the final version of our method, which we have named DMRPC. Then, DMRPC and another DMR method, coMethDMR, were used to perform epigenome-wide analyses of several phenotypes known to have multiple associated methylation loci (age, sex, and smoking) in a discovery and a replication cohort. Among regions that were analysed by both methods, DMRPC identified 50% more genome-wide significant age-associated DMRs than coMethDMR. The replication rate for the loci that were identified by only DMRPC was higher than the rate for those that were identified by only coMethDMR (90% for DMRPC vs. 76% for coMethDMR). Furthermore, DMRPC identified replicable associations in regions of moderate between-CpG correlation which are typically not analysed by coMethDMR. For the analyses of sex and smoking, the advantage of DMRPC was less clear. In conclusion, DMRPC is a new powerful DMR discovery tool that retains power in genomic regions with moderate correlation across CpGs.

Validation of the new EPIC DNA methylation microarray (900K EPIC v2) for high-throughput profiling of the human DNA methylome

DNA methylation, one of the best characterized epigenetic marks in the human genome, plays a pivotal role in gene transcription regulation and other biological processes in humans. On top of that, the DNA methylome undergoes profound changes in cancer and other disorders. However, large-scale and population-based studies are limited by high costs and the need for considerable expertise in data analysis for whole-genome bisulphite-sequencing methodologies. Following the success of the EPIC DNA methylation microarray, the newly developed Infinium HumanMethylationEPIC version 2.0 (900K EPIC v2) is now available. This new array contains more than 900,000 CpG probes covering the human genome and excluding masked probes from the previous version. The 900K EPIC v2 microarray adds more than 200,000 probes covering extra DNA cis-regulatory regions such as enhancers, super-enhancers and CTCF binding regions. Herein, we have technically and biologically validated the new methylation array to show its high reproducibility and consistency among technical replicates and with DNA extracted from FFPE tissue. In addition, we have hybridized primary normal and tumoural tissues and cancer cell lines from different sources and tested the robustness of the 900K EPIC v2 microarray when analysing the different DNA methylation profiles. The validation highlights the improvements offered by the new array and demonstrates the versatility of this updated tool for characterizing the DNA methylome in human health and disease.

Gaps in current methods to detect polymorphic CpGs from Illumina Infinium human methylation microarrays and exploring their potential impact in multi-EWAS analyses

DNA methylation (DNAm) epigenome-wide association studies (EWAS) have been performed on diverse ethnicities to discover novel biomarkers associated with various diseases, such as cancers, autoimmune diseases, and neurological disorders. However, genetic polymorphisms can influence DNAm levels resulting in methylation quantitative trait loci (meQTL). These can be either direct effects, by altering the sequence of the methylation (CpG) site itself, or, in the case of array-based measures, indirectly altering the detection probe-binding site interaction. Given that genetic variant frequencies associated with meQTL can differ between population groups, these have the potential to confound EWAS observations, particularly in multi-ethnic populations. In this study, we analysed publicly available DNA methylation profiles (450K array), consisting of 1342 individuals from 6 distinct ancestral groups. We investigate two distinct tools (GapHunter and MethylToSNP) specifically designed to identify CpG sites that may be influenced by genetic variation. Results from this aggregated trans-ancestral epigenome-wide dataset suggest that both tools fail to consistently identify not only rarer (MAF < 0.05) genetic variant effects but also more than half of sites predicted to be associated with variants with much higher allele frequencies (MAF >0.2). In addition, there is a relatively low concordance in the detection of polymorphic CpGs between GapHunter and MethylToSNP. Screening of CpG site associations from EWAS using either of these tools is unlikely to be a robust or comprehensive means of identifying all genetic variant confounding effects.

Differential methylation pattern in pubertal girls associated with biochemical premature adrenarche

Biochemical premature adrenarche is defined by elevated serum DHEAS [≥40 μg/dL] before age 8 y in girls. This condition is receiving more attention due to its association with obesity, hyperinsulinemia, dyslipidemia, and polycystic ovary syndrome. Nevertheless, the link between early androgen excess and these risk factors remains unknown. Epigenetic modifications, and specifically DNA methylation, have been associated with the initiation and progression of numerous disorders, including obesity and insulin resistance. The aim of this study was to determine if prepubertal androgen exposure is associated with a different methylation profile in pubertal girls. Eighty-six healthy girls were studied. At age 7 y, anthropometric measurements were begun and DHEAS levels were determined. Girls were classified into Low DHEAS (LD) [<42 μg/dL] and High DHEAS (HD) [≥42 μg/dL] groups. At Tanner stages 2 and 4 a DNA methylation microarray was performed to identify differentially methylated CpG positions (DMPs) between HD and LD groups. We observed a differential methylation pattern between pubertal girls with and without biochemical PA. Moreover, a set of DNA methylation markers, selected by the LASSO method, successfully distinguished between HD and LD girls regardless of Tanner stage. Additionally, a subset of these markers were significantly associated with glucose-related measures such as insulin level, HOMA-IR, and glycaemia. This pilot study provides evidence consistent with the hypothesis that high DHEAS concentration, or its hormonally active metabolites, may induce a unique blood methylation signature in pubertal girls, and that this methylation pattern is associated with altered glucose metabolism.

Comparison of global DNA methylation analysis by whole genome bisulfite sequencing and the Infinium Mouse Methylation BeadChip using fresh and fresh-frozen mouse epidermis

Until recently, studying the murine methylome was restricted to sequencing-based methods. In this study we compared the global DNA methylation levels of hairless mouse epidermis using the recently released Infinium Mouse Methylation BeadChip from Illumina and whole genome bisulphite sequencing (WGBS). We also studied the effect of sample storage conditions by using fresh and fresh-frozen epidermis. The DNA methylation levels of 123,851 CpG sites covered by both the BeadChip and WGBS were compared. DNA methylation levels obtained with WGBS and the BeadChip were strongly correlated (Pearson correlation r = 0.984). We applied a threshold of 15 reads for the WGBS methylation analysis. Even at a threshold of 10 reads, we observed no substantial difference in DNA methylation levels compared with that obtained with the BeadChip. The DNA methylation levels from the fresh and the fresh-frozen samples were strongly correlated when analysed with both the BeadChip (r = 0.999) and WGBS (r = 0.994). We conclude that the two methods of analysis generally work equally well for studies of DNA methylation of mouse epidermis and find that fresh and fresh-frozen epidermis can generally be used equally well. The choice of method will depend on the specific study's aims and the available resources in the laboratory.

An epigenome-wide analysis of sex hormone levels and DNA methylation in male blood samples

Endogenous sex hormones and DNA methylation both play important roles in various diseases. However, their interplay is largely unknown. A deeper understanding of their interrelationships could provide new insights into the pathology of disease development. We, therefore, investigated associations between circulating sex hormones, sex hormone binding globulin (SHBG), and DNA methylation in blood, using samples from 77 men (65 with repeated samples), from the population-based Northern Sweden Health and Disease Study (NSHDS). DNA methylation was measured in buffy coat using the Infinium Methylation EPIC BeadChip (Illumina). Sex hormone (oestradiol, oestrone, testosterone, androstenedione, dehydroepiandrosterone, and progesterone) and SHBG concentrations were measured in plasma using a high-performance liquid chromatography tandem mass spectrometry (LC/MS-MS) method and an enzyme-linked immunoassay, respectively. Associations between sex hormones, SHBG, and DNA methylation were estimated using both linear regression and mixed-effects models. Additionally, we used the comb-p method to identify differentially methylated regions based on nearby P values. We identified one novel CpG site (cg14319657), at which DNA methylation was associated with dehydroepiandrosterone, surpassing a genome-wide significance level. In addition, more than 40 differentially methylated regions were associated with levels of sex hormones and SHBG and several of these mapped to genes involved in hormone-related diseases. Our findings support a relationship between circulating sex hormones and DNA methylation and suggest that further investigation is warranted, both for validation, further exploration and to gain a deeper understanding of the mechanisms and potential consequences for health and disease.

Prenatal exposures to endocrine disrupting chemicals: The role of multi-omics in understanding toxicity

Endocrine disrupting chemicals (EDCs) are a diverse group of toxicants detected in populations globally. Prenatal EDC exposures impact birth and childhood outcomes. EDCs work through persistent changes at the molecular, cellular, and organ level. Molecular and biochemical signals or 'omics' can be measured at various functional levels - including the epigenome, transcriptome, proteome, metabolome, and the microbiome. In this narrative review, we introduce each omics and give examples of associations with prenatal EDC exposures. There is substantial research on epigenomic modifications in offspring exposed to EDCs during gestation, and a growing number of studies evaluating the transcriptome, proteome, metabolome, or microbiome in response to these exposures. Multi-omics, integrating data across omics layers, may improve understanding of disrupted function pathways related to early life exposures. We highlight several data integration methods to consider in multi-omics studies. Information from multi-omics can improve understanding of the biological processes and mechanisms underlying prenatal EDC toxicity.

Multifaceted immune dysregulation characterizes individuals at-risk for rheumatoid arthritis

Molecular markers of autoimmunity, such as antibodies to citrullinated protein antigens (ACPA), are detectable prior to inflammatory arthritis (IA) in rheumatoid arthritis (RA) and may define a state that is 'at-risk' for future RA. Here we present a cross-sectional comparative analysis among three groups that include ACPA positive individuals without IA (At-Risk), ACPA negative individuals and individuals with early, ACPA positive clinical RA (Early RA). Differential methylation analysis among the groups identifies non-specific dysregulation in peripheral B, memory and naïve T cells in At-Risk participants, with more specific immunological pathway abnormalities in Early RA. Tetramer studies show increased abundance of T cells recognizing citrullinated (cit) epitopes in At-Risk participants, including expansion of T cells reactive to citrullinated cartilage intermediate layer protein I (cit-CILP); these T cells have Th1, Th17, and T stem cell memory-like phenotypes. Antibody-antigen array analyses show that antibodies targeting cit-clusterin, cit-fibrinogen and cit-histone H4 are elevated in At-Risk and Early RA participants, with the highest levels of antibodies detected in those with Early RA. These findings indicate that an ACPA positive at-risk state is associated with multifaceted immune dysregulation that may represent a potential opportunity for targeted intervention.

Nepalese indoor cookstove smoke extracts alter human airway epithelial gene expression, DNA methylation and hydroxymethylation

Household air pollution caused by inefficient cooking practices causes 4 million deaths a year worldwide. In Nepal, 86% of the rural population use solid fuels for cooking. Over 25% of premature deaths associated with air pollution are respiratory in nature. Here we aimed to identify molecular signatures of different cookstove and fuel type exposures in human airway epithelial cells, to understand the mechanisms mediating cook stove smoke induced lung disease. Primary human airway epithelial cells in submerged culture were exposed to traditional cook stove (TCS), improved cook stove (ICS) and liquefied petroleum gas (LPG) stove smoke extracts. Changes to gene expression, DNA methylation and hydroxymethylation were measured by bulk RNA sequencing and HumanMethylationEPIC BeadChip following oxidative bisulphite conversion, respectively. TCS smoke extract alone reproducibly caused changes in the expression of 52 genes enriched for oxidative stress pathways. TCS, ICS and LPG smoke extract exposures were associated with distinct changes to DNA methylation and hydroxymethylation. A subset of TCS induced genes were associated with differentially methylated and/or hydroxymethylated CpGs sites, and enriched for the ferroptosis pathway and the upstream regulator NFE2L2. DNA methylation and hydroxymethylation changes not associated with a concurrent change in gene expression, were linked to biological processes and molecular pathways important to airway health, including neutrophil function, transforming growth factor beta signalling, GTPase activity, and cell junction organisation. Our data identified differential impacts of TCS, ICS and LPG cook stove smoke on the human airway epithelium transcriptome, DNA methylome and hydroxymethylome and provide further insight into the association between indoor air pollution exposure and chronic lung disease mechanisms.

Associations between the New DNA-Methylation-Based Telomere Length Estimator, the Mediterranean Diet and Genetics in a Spanish Population at High Cardiovascular Risk

Biological aging is a relevant risk factor for chronic diseases, and several indicators for measuring this factor have been proposed, with telomere length (TL) among the most studied. Oxidative stress may regulate telomere shortening, which is implicated in the increased risk. Using a novel estimator for TL, we examined whether adherence to the Mediterranean diet (MedDiet), a highly antioxidant-rich dietary pattern, is associated with longer TL. We determined TL using DNA methylation algorithms (DNAmTL) in 414 subjects at high cardiovascular risk from Spain. Adherence to the MedDiet was assessed by a validated score, and genetic variants in candidate genes and at the genome-wide level were analyzed. We observed several significant associations (p < 0.05) between DNAmTL and candidate genes (TERT, TERF2, RTEL1, and DCAF4), contributing to the validity of DNAmTL as a biomarker in this population. Higher adherence to the MedDiet was associated with lower odds of having a shorter TL in the whole sample (OR = 0.93; 95% CI: 0.85-0.99; p = 0.049 after fully multivariate adjustment). Nevertheless, this association was stronger in women than in men. Likewise, in women, we observed a direct association between adherence to the MedDiet score and DNAmTL as a continuous variable (beta = 0.015; SE: 0.005; p = 0.003), indicating that a one-point increase in adherence was related to an average increase of 0.015 ± 0.005 kb in TL. Upon examination of specific dietary items within the global score, we found that fruits, fish, "sofrito", and whole grains exhibited the strongest associations in women. The novel score combining these items was significantly associated in the whole population. In the genome-wide association study (GWAS), we identified ten polymorphisms at the suggestive level of significance (p < 1 × 10-5) for DNAmTL (intergenics, in the IQSEC1, NCAPG2, and ABI3BP genes) and detected some gene-MedDiet modulations on DNAmTL. As this is the first study analyzing the DNAmTL estimator, genetics, and modulation by the MedDiet, more studies are needed to confirm these findings.

Associations Between Epigenetic Age Acceleration and microRNA Expression Among U.S. Firefighters

Epigenetic changes may be biomarkers of health. Epigenetic age acceleration (EAA), the discrepancy between epigenetic age measured via epigenetic clocks and chronological age, is associated with morbidity and mortality. However, the intersection of epigenetic clocks with microRNAs (miRNAs) and corresponding miRNA-based health implications have not been evaluated. We analyzed DNA methylation and miRNA profiles from blood sampled among 332 individuals enrolled across 2 U.S.-based firefighter occupational studies (2015-2018 and 2018-2020). We considered 7 measures of EAA in leukocytes (PhenoAge, GrimAge, Horvath, skin-blood, and Hannum epigenetic clocks, and extrinsic and intrinsic epigenetic age acceleration). We identified miRNAs associated with EAA using individual linear regression models, adjusted for sex, race/ethnicity, chronological age, and cell type estimates, and investigated downstream effects of associated miRNAs with miRNA enrichment analyses and genomic annotations. On average, participants were 38 years old, 88% male, and 75% non-Hispanic white. We identified 183 of 798 miRNAs associated with EAA (FDR q < 0.05); 126 with PhenoAge, 59 with GrimAge, 1 with Horvath, and 1 with the skin-blood clock. Among miRNAs associated with Horvath and GrimAge, there were 61 significantly enriched disease annotations including age-related metabolic and cardiovascular conditions and several cancers. Enriched pathways included those related to proteins and protein modification. We identified miRNAs associated with EAA of multiple epigenetic clocks. PhenoAge had more associations with individual miRNAs, but GrimAge and Horvath had greater implications for miRNA-associated pathways. Understanding the relationship between these epigenetic markers could contribute to our understanding of the molecular underpinnings of aging and aging-related diseases.

Stability of DNA-Methylation Profiles of Biological Aging in Children and Adolescents

Background and Objectives

Methylation profile scores (MPSs) index biological aging and aging-related disease in adults and are cross-sectionally associated with social determinants of health in childhood. MPSs thus provide an opportunity to trace how aging-related biology responds to environmental changes in early life. Information regarding the stability of MPSs in early life is currently lacking.

Method

We use longitudinal data from children and adolescents ages 8-18 (N = 428, M age = 12.15 years) from the Texas Twin Project. Participants contributed two waves of salivary DNA-methylation data (mean lag = 3.94 years), which were used to construct four MPSs reflecting multi-system physiological decline and mortality risk (PhenoAgeAccel and GrimAgeAccel), pace of biological aging (DunedinPACE), and cognitive function (Epigenetic-g). Furthermore, we exploit variation among participants in whether they were exposed to the COVID-19 pandemic during the course of study participation, in order to test how a historical period characterized by environmental disruption might affect children's aging-related MPSs.

Results

All MPSs showed moderate longitudinal stability (test-retest rs = 0.42, 0.44, 0.46, 0.51 for PhenoAgeAccel, GrimAgeAccel, and Epigenetic-g, and DunedinPACE, respectively). No differences in the stability of MPSs were apparent between those whose second assessment took place after the onset of the COVID-19 pandemic vs. those for whom both assessments took place prior to the pandemic.

Conclusions

Aging-related DNA-methylation patterns are less stable in childhood than has been previously observed in adulthood. Further developmental research on the methylome is necessary to understand which environmental perturbations in childhood impact trajectories of biological aging and when children are most sensitive to those impacts.

Genome-wide association study of placental weight identifies distinct and shared genetic influences between placental and fetal growth

A well-functioning placenta is essential for fetal and maternal health throughout pregnancy. Using placental weight as a proxy for placental growth, we report genome-wide association analyses in the fetal (n = 65,405), maternal (n = 61,228) and paternal (n = 52,392) genomes, yielding 40 independent association signals. Twenty-six signals are classified as fetal, four maternal and three fetal and maternal. A maternal parent-of-origin effect is seen near KCNQ1. Genetic correlation and colocalization analyses reveal overlap with birth weight genetics, but 12 loci are classified as predominantly or only affecting placental weight, with connections to placental development and morphology, and transport of antibodies and amino acids. Mendelian randomization analyses indicate that fetal genetically mediated higher placental weight is causally associated with preeclampsia risk and shorter gestational duration. Moreover, these analyses support the role of fetal insulin in regulating placental weight, providing a key link between fetal and placental growth.

Epigenome-Wide Meta-analysis Reveals Associations Between Dietary Glycemic Index and Glycemic Load and DNA Methylation in Children and Adolescents of Different Body Sizes

OBJECTIVE

Dietary glycemic index (GI) and glycemic load (GL) are associated with cardiometabolic health in children and adolescents, with potential distinct effects in people with increased BMI. DNA methylation (DNAm) may mediate these effects. Thus, we conducted meta-analyses of epigenome-wide association studies (EWAS) between dietary GI and GL and blood DNAm of children and adolescents.

RESEARCH DESIGN AND METHODS

We calculated dietary GI and GL and performed EWAS in children and adolescents (age range: 4.5-17 years) from six cohorts (N = 1,187). We performed stratified analyses of participants with normal weight (n = 801) or overweight or obesity (n = 386). We performed look-ups for the identified cytosine-phosphate-guanine (CpG) sites (false discovery rate [FDR] <0.05) with tissue-specific gene expression of 832 blood and 223 subcutaneous adipose tissue samples from children and adolescents.

RESULTS

Dietary GL was positively associated with DNAm of cg20274553 (FDR <0.05), annotated to WDR27. Several CpGs were identified in the normal-weight (GI: 85; GL: 17) and overweight or obese (GI: 136; GL: 298; FDR <0.05) strata, and none overlapped between strata. In participants with overweight or obesity, identified CpGs were related to RNA expression of genes associated with impaired metabolism (e.g., FRAT1, CSF3).

CONCLUSIONS

We identified 537 associations between dietary GI and GL and blood DNAm, mainly in children and adolescents with overweight or obesity. High-GI and/or -GL diets may influence epigenetic gene regulation and thereby promote metabolic derangements in young people with increased BMI.

Genome-wide neonatal epigenetic changes associated with maternal exposure to the COVID-19 pandemic

BACKGROUND

During gestation, stressors to the fetus, including viral exposure or maternal psychological distress, can fundamentally alter the neonatal epigenome, and may be associated with long-term impaired developmental outcomes. The impact of in utero exposure to the COVID-19 pandemic on the newborn epigenome has yet to be described.

METHODS

This study aimed to determine whether there are unique epigenetic signatures in newborns who experienced otherwise healthy pregnancies that occurred during the COVID-19 pandemic (Project RESCUE). The pre-pandemic control and pandemic cohorts (Project RESCUE) included in this study are part of a prospective observational and longitudinal cohort study that evaluates the impact of elevated prenatal maternal stress during the COVID-19 pandemic on early childhood neurodevelopment. Using buccal swabs collected at birth, differential DNA methylation analysis was performed using the Infinium MethylationEPIC arrays and linear regression analysis. Pathway analysis and gene ontology enrichment were performed on resultant gene lists.

RESULTS

Widespread differential methylation was found between neonates exposed in utero to the pandemic and pre-pandemic neonates. In contrast, there were no apparent epigenetic differences associated with maternal COVID-19 infection during pregnancy. Differential methylation was observed among genomic sites that underpin important neurological pathways that have been previously reported in the literature to be differentially methylated because of prenatal stress, such as NR3C1.

CONCLUSIONS

The present study reveals potential associations between exposure to the COVID-19 pandemic during pregnancy and subsequent changes in the newborn epigenome. While this finding warrants further investigation, it is a point that should be considered in any study assessing newborn DNA methylation studies obtained during this period, even in otherwise healthy pregnancies.

DNA methylation of insulin signaling pathways is associated with HOMA2-IR in primary myoblasts from older adults

BACKGROUND

While ageing is associated with increased insulin resistance (IR), the molecular mechanisms underlying increased IR in the muscle, the primary organ for glucose clearance, have yet to be elucidated in older individuals. As epigenetic processes are suggested to contribute to the development of ageing-associated diseases, we investigated whether differential DNA methylation was associated with IR in human primary muscle stem cells (myoblasts) from community-dwelling older individuals.

METHODS

We measured DNA methylation (Infinium HumanMethylationEPIC BeadChip) in myoblast cultures from vastus lateralis biopsies (119 males/females, mean age 78.24 years) from the Hertfordshire Sarcopenia Study extension (HSSe) and examined differentially methylated cytosine phosphate guanine (CpG) sites (dmCpG), regions (DMRs) and gene pathways associated with HOMA2-IR, an index for the assessment of insulin resistance, and levels of glycated hemoglobin HbA1c.

RESULTS

Thirty-eight dmCpGs (false discovery rate (FDR) < 0.05) were associated with HOMA2-IR, with dmCpGs enriched in genes linked with JNK, AMPK and insulin signaling. The methylation signal associated with HOMA2-IR was attenuated after the addition of either BMI (6 dmCpGs), appendicular lean mass index (ALMi) (7 dmCpGs), grip strength (15 dmCpGs) or gait speed (23 dmCpGs) as covariates in the model. There were 8 DMRs (Stouffer < 0.05) associated with HOMA2-IR, including DMRs within T-box transcription factor (TBX1) and nuclear receptor subfamily-2 group F member-2 (NR2F2); the DMRs within TBX1 and NR2F2 remained associated with HOMA2-IR after adjustment for BMI, ALMi, grip strength or gait speed. Forty-nine dmCpGs and 21 DMRs were associated with HbA1c, with cg13451048, located within exoribonuclease family member 3 (ERI3) associated with both HOMA2-IR and HbA1c. HOMA2-IR and HbA1c were not associated with accelerated epigenetic ageing.

CONCLUSIONS

These findings suggest that insulin resistance is associated with differential DNA methylation in human primary myoblasts with both muscle mass and body composition making a significant contribution to the methylation changes associated with IR.

Open access-enabled evaluation of epigenetic age acceleration in colorectal cancer and development of a classifier with diagnostic potential

Aberrant DNA methylation (DNAm) is known to be associated with the aetiology of cancer, including colorectal cancer (CRC). In the past, the availability of open access data has been the main driver of innovative method development and research training. However, this is increasingly being eroded by the move to controlled access, particularly of medical data, including cancer DNAm data. To rejuvenate this valuable tradition, we leveraged DNAm data from 1,845 samples (535 CRC tumours, 522 normal colon tissues adjacent to tumours, 72 colorectal adenomas, and 716 normal colon tissues from healthy individuals) from 14 open access studies deposited in NCBI GEO and ArrayExpress. We calculated each sample's epigenetic age (EA) using eleven epigenetic clock models and derived the corresponding epigenetic age acceleration (EAA). For EA, we observed that most first- and second-generation epigenetic clocks reflect the chronological age in normal tissues adjacent to tumours and healthy individuals [e.g., Horvath (r = 0.77 and 0.79), Zhang elastic net (EN) (r = 0.70 and 0.73)] unlike the epigenetic mitotic clocks (EpiTOC, HypoClock, MiAge) (r < 0.3). For EAA, we used PhenoAge, Wu, and the above mitotic clocks and found them to have distinct distributions in different tissue types, particularly between normal colon tissues adjacent to tumours and cancerous tumours, as well as between normal colon tissues adjacent to tumours and normal colon tissue from healthy individuals. Finally, we harnessed these associations to develop a classifier using elastic net regression (with lasso and ridge regularisations) that predicts CRC diagnosis based on a patient's sex and EAAs calculated from histologically normal controls (i.e., normal colon tissues adjacent to tumours and normal colon tissue from healthy individuals). The classifier demonstrated good diagnostic potential with ROC-AUC = 0.886, which suggests that an EAA-based classifier trained on relevant data could become a tool to support diagnostic/prognostic decisions in CRC for clinical professionals. Our study also reemphasises the importance of open access clinical data for method development and training of young scientists. Obtaining the required approvals for controlled access data would not have been possible in the timeframe of this study.

Critical evaluation of the reliability of DNA methylation probes on the Illumina MethylationEPIC 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 arrays. We conducted a comprehensive assessment of the EPIC array probe reliability using 138 duplicated blood DNAm samples generated by the Alzheimer's Disease Neuroimaging Initiative study. We introduced a novel statistical measure, the modified intraclass correlation, to better account for the disagreement in duplicate measurements. 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 reliable information for probes on the EPIC array, will serve as a valuable resource for future DNAm studies.

Critical evaluation of the reliability of DNA methylation probes on the Illumina MethylationEPIC 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 arrays. We conducted a comprehensive assessment of the EPIC array probe reliability using 138 duplicated blood DNAm samples generated by the Alzheimer's Disease Neuroimaging Initiative study. We introduced a novel statistical measure, the modified intraclass correlation, to better account for the disagreement in duplicate measurements. 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 reliable information for probes on the EPIC array, will serve as a valuable resource for future DNAm studies.

Diagnostic Utility of Genome-wide DNA Methylation Analysis in Genetically Unsolved Developmental and Epileptic Encephalopathies and Refinement of a CHD2 Episignature

Sequence-based genetic testing currently identifies causative genetic variants in ∼50% of individuals with developmental and epileptic encephalopathies (DEEs). Aberrant changes in DNA methylation are implicated in various neurodevelopmental disorders but remain unstudied in DEEs. Rare epigenetic variations ("epivariants") can drive disease by modulating gene expression at single loci, whereas genome-wide DNA methylation changes can result in distinct "episignature" biomarkers for monogenic disorders in a growing number of rare diseases. Here, we interrogate the diagnostic utility of genome-wide DNA methylation array analysis on peripheral blood samples from 516 individuals with genetically unsolved DEEs who had previously undergone extensive genetic testing. We identified rare differentially methylated regions (DMRs) and explanatory episignatures to discover causative and candidate genetic etiologies in 10 individuals. We then used long-read sequencing to identify DNA variants underlying rare DMRs, including one balanced translocation, three CG-rich repeat expansions, and two copy number variants. We also identify pathogenic sequence variants associated with episignatures; some had been missed by previous exome sequencing. Although most DEE genes lack known episignatures, the increase in diagnostic yield for DNA methylation analysis in DEEs is comparable to the added yield of genome sequencing. Finally, we refine an episignature for CHD2 using an 850K methylation array which was further refined at higher CpG resolution using bisulfite sequencing to investigate potential insights into CHD2 pathophysiology. Our study demonstrates the diagnostic yield of genome-wide DNA methylation analysis to identify causal and candidate genetic causes as ∼2% (10/516) for unsolved DEE cases.

Epigenetic and Metabolomic Biomarkers for Biological Age: A Comparative Analysis of Mortality and Frailty Risk

Biological age captures a person's age-related risk of unfavorable outcomes using biophysiological information. Multivariate biological age measures include frailty scores and molecular biomarkers. These measures are often studied in isolation, but here we present a large-scale study comparing them. In 2 prospective cohorts (n = 3 222), we compared epigenetic (DNAm Horvath, DNAm Hannum, DNAm Lin, DNAm epiTOC, DNAm PhenoAge, DNAm DunedinPoAm, DNAm GrimAge, and DNAm Zhang) and metabolomic-based (MetaboAge and MetaboHealth) biomarkers in reflection of biological age, as represented by 5 frailty measures and overall mortality. Biomarkers trained on outcomes with biophysiological and/or mortality information outperformed age-trained biomarkers in frailty reflection and mortality prediction. DNAm GrimAge and MetaboHealth, trained on mortality, showed the strongest association with these outcomes. The associations of DNAm GrimAge and MetaboHealth with frailty and mortality were independent of each other and of the frailty score mimicking clinical geriatric assessment. Epigenetic, metabolomic, and clinical biological age markers seem to capture different aspects of aging. These findings suggest that mortality-trained molecular markers may provide novel phenotype reflecting biological age and strengthen current clinical geriatric health and well-being assessment.

Epigenomic signature of major congenital heart defects in newborns with Down syndrome

BACKGROUND

Congenital heart defects (CHDs) affect approximately half of individuals with Down syndrome (DS), but the molecular reasons for incomplete penetrance are unknown. Previous studies have largely focused on identifying genetic risk factors associated with CHDs in individuals with DS, but comprehensive studies of the contribution of epigenetic marks are lacking. We aimed to identify and characterize DNA methylation differences from newborn dried blood spots (NDBS) of DS individuals with major CHDs compared to DS individuals without CHDs.

METHODS

We used the Illumina EPIC array and whole-genome bisulfite sequencing (WGBS) to quantitate DNA methylation for 86 NDBS samples from the California Biobank Program: (1) 45 DS-CHD (27 female, 18 male) and (2) 41 DS non-CHD (27 female, 14 male). We analyzed global CpG methylation and identified differentially methylated regions (DMRs) in DS-CHD versus DS non-CHD comparisons (both sex-combined and sex-stratified) corrected for sex, age of blood collection, and cell-type proportions. CHD DMRs were analyzed for enrichment in CpG and genic contexts, chromatin states, and histone modifications by genomic coordinates and for gene ontology enrichment by gene mapping. DMRs were also tested in a replication dataset and compared to methylation levels in DS versus typical development (TD) WGBS NDBS samples.

RESULTS

We found global CpG hypomethylation in DS-CHD males compared to DS non-CHD males, which was attributable to elevated levels of nucleated red blood cells and not seen in females. At a regional level, we identified 58, 341, and 3938 CHD-associated DMRs in the Sex Combined, Females Only, and Males Only groups, respectively, and used machine learning algorithms to select 19 Males Only loci that could distinguish CHD from non-CHD. DMRs in all comparisons were enriched for gene exons, CpG islands, and bivalent chromatin and mapped to genes enriched for terms related to cardiac and immune functions. Lastly, a greater percentage of CHD-associated DMRs than background regions were differentially methylated in DS versus TD samples.

CONCLUSIONS

A sex-specific signature of DNA methylation was detected in NDBS of DS-CHD compared to DS non-CHD individuals. This supports the hypothesis that epigenetics can reflect the variability of phenotypes in DS, particularly CHDs.

Epigenomic response to albuterol treatment in asthma-relevant airway epithelial cells

BACKGROUND

Albuterol is the first-line asthma medication used in diverse populations. Although DNA methylation (DNAm) is an epigenetic mechanism involved in asthma and bronchodilator drug response (BDR), no study has assessed whether albuterol could induce changes in the airway epithelial methylome. We aimed to characterize albuterol-induced DNAm changes in airway epithelial cells, and assess potential functional consequences and the influence of genetic variation and asthma-related clinical variables.

RESULTS

We followed a discovery and validation study design to characterize albuterol-induced DNAm changes in paired airway epithelial cultures stimulated in vitro with albuterol. In the discovery phase, an epigenome-wide association study using paired nasal epithelial cultures from Puerto Rican children (n = 97) identified 22 CpGs genome-wide associated with repeated-use albuterol treatment (p < 9 × 10-8). Albuterol predominantly induced a hypomethylation effect on CpGs captured by the EPIC array across the genome (probability of hypomethylation: 76%, p value = 3.3 × 10-5). DNAm changes on the CpGs cg23032799 (CREB3L1), cg00483640 (MYLK4-LINC01600), and cg05673431 (KSR1) were validated in nasal epithelia from 10 independent donors (false discovery rate [FDR] < 0.05). The effect on the CpG cg23032799 (CREB3L1) was cross-tissue validated in bronchial epithelial cells at nominal level (p = 0.030). DNAm changes in these three CpGs were shown to be influenced by three independent genetic variants (FDR < 0.05). In silico analyses showed these polymorphisms regulated gene expression of nearby genes in lungs and/or fibroblasts including KSR1 and LINC01600 (6.30 × 10-14 ≤ p ≤ 6.60 × 10-5). Additionally, hypomethylation at the CpGs cg10290200 (FLNC) and cg05673431 (KSR1) was associated with increased gene expression of the genes where they are located (FDR < 0.05). Furthermore, while the epigenetic effect of albuterol was independent of the asthma status, severity, and use of medication, BDR was nominally associated with the effect on the CpG cg23032799 (CREB3L1) (p = 0.004). Gene-set enrichment analyses revealed that epigenomic modifications of albuterol could participate in asthma-relevant processes (e.g., IL-2, TNF-α, and NF-κB signaling pathways). Finally, nine differentially methylated regions were associated with albuterol treatment, including CREB3L1, MYLK4, and KSR1 (adjusted p value < 0.05).

CONCLUSIONS

This study revealed evidence of epigenetic modifications induced by albuterol in the mucociliary airway epithelium. The epigenomic response induced by albuterol might have potential clinical implications by affecting biological pathways relevant to asthma.

Methylomic, proteomic, and metabolomic correlates of traffic-related air pollution: A systematic review, pathway analysis, and network analysis relating traffic-related air pollution to subclinical and clinical cardiorespiratory outcomes

A growing body of literature has attempted to characterize how traffic-related air pollution (TRAP) affects molecular and subclinical biological processes in ways that could lead to cardiorespiratory disease. To provide a streamlined synthesis of what is known about the multiple mechanisms through which TRAP could lead cardiorespiratory pathology, we conducted a systematic review of the epidemiological literature relating TRAP exposure to methylomic, proteomic, and metabolomic biomarkers in adult populations. Using the 139 papers that met our inclusion criteria, we identified the omic biomarkers significantly associated with short- or long-term TRAP and used these biomarkers to conduct pathway and network analyses. We considered the evidence for TRAP-related associations with biological pathways involving lipid metabolism, cellular energy production, amino acid metabolism, inflammation and immunity, coagulation, endothelial function, and oxidative stress. Our analysis suggests that an integrated multi-omics approach may provide critical new insights into the ways TRAP could lead to adverse clinical outcomes. We advocate for efforts to build a more unified approach for characterizing the dynamic and complex biological processes linking TRAP exposure and subclinical and clinical disease, and highlight contemporary challenges and opportunities associated with such efforts.

Placental DNA methylation analysis of selective fetal growth restriction in monochorionic twins reveals aberrant methylated CYP11A1 gene for fetal growth restriction

Fetal growth restriction (FGR) is associated with increased susceptibility to perinatal morbidity and mortality. Evidence suggests that epigenetic changes play critical roles in the regulation of fetal growth. We sought to present a comprehensive analysis of the associations between placental DNA methylation and selective fetal growth restriction (sFGR), which is a severe complication of monochorionic twin pregnancies, characterized by one fetus experiencing restricted growth. Genome-wide methylation analysis was performed on 24 placental samples obtained from 12 monochorionic twins with sFGR (Cohort 1) using Illumina Infinium MethylationEPIC BeadChip. Integrative analysis of our EPIC data and two previous placental methylation studies of sFGR (a total of 30 placental samples from 15 sFGR twins) was used to identify convincing differential promoter methylation. Validation analysis was performed on the placentas from 15 sFGR twins (30 placental samples), 15 FGR singletons, and 14 control singletons (Cohort 2) using pyrosequencing, quantitative real-time polymerase chain reaction, western blot, and immunohistochemistry (IHC). A globe shift toward hypomethylation was identified in the placentas of growth-restricted fetuses compared with the placentas of normal fetuses in monochorionic twins, including 5625 hypomethylated CpGs and 452 hypermethylated CpGs, especially in the regions of CpG islands, gene-body and promoters. The analysis of pathways revealed dysregulation primarily in steroid hormone biosynthesis, metabolism, cell adhesion, signaling transduction, and immune response. Integrative analysis revealed a differentially methylated promoter region in the CYP11A1 gene, encoding a rate-limiting enzyme of steroidogenesis converting cholesterol to pregnenolone. The CYP11A1 gene was validated to have hypomethylation and higher mRNA expression in sFGR twins and FGR singletons. In conclusion, our findings suggested that the changes in placental DNA methylation pattern in sFGR may have functional implications for differentially methylated genes and regulatory regions. The study provides reliable evidence for identifying abnormally methylated CYP11A1 gene in the placenta of sFGR.

Associations between blood leukocyte DNA methylation and sustained attention in mid-to-late childhood

Aims: To identify associations between DNA methylation (DNAm) across the epigenome and symptoms related to attention-deficit/hyperactivity disorder in a population of Hispanic children. Materials & methods: Among 517 participants in the ELEMENT study aged 9-18 years, we conducted an epigenome-wide association study examining associations between blood leukocyte DNAm and performance on the Conners' continuous performance test (CPT3). Results: DNAm at loci in or near ZNF814, ELF4 and OR6K6 and functional enrichment for gene pathways pertaining to ferroptosis, inflammation, immune response and neurotransmission were significantly related to CPT3 scores. Conclusion: DNAm was associated with CPT3 performance. Further analysis is warranted to understand how these genes and enriched pathways contribute to attention-deficit/hyperactivity disorder.

Getting the chronological age out of DNA: using insights of age-dependent DNA methylation for forensic DNA applications

BACKGROUND

DNA analysis for forensic investigations has a long tradition with important developments and optimizations since its first application. Traditionally, short tandem repeats analysis has been the most powerful method for the identification of individuals. However, in addition, epigenetic changes, i.e., DNA methylation, came into focus of forensic DNA research. Chronological age prediction is one promising application to allow for narrowing the pool of possible individuals who caused a trace, as well as to support the identification of unknown bodies and for age verification of living individuals.

OBJECTIVE

This review aims to provide an overview of the current knowledge, possibilities, and (current) limitations about DNA methylation-based chronological age prediction with emphasis on forensic application.

METHODS

The development, implementation and application of age prediction tools requires a deep understanding about the biological background, the analysis methods, the age-dependent DNA methylation markers, as well as the mathematical models for age prediction and their evaluation. Furthermore, additional influences can have an impact. Therefore, the literature was evaluated in respect to these diverse topics.

CONCLUSION

The numerous research efforts in recent years have led to a rapid change in our understanding of the application of DNA methylation for chronological age prediction, which is now on the way to implementation and validation. Knowledge of the various aspects leads to a better understanding and allows a more informed interpretation of DNAm quantification results, as well as the obtained results by the age prediction tools.

Cell-free DNA Methylation as a Predictive Biomarker of Response to Neoadjuvant Chemotherapy for Patients with Muscle-invasive Bladder Cancer in SWOG S1314

BACKGROUND

Neoadjuvant chemotherapy (NAC) is the standard of care in muscle-invasive bladder cancer (MIBC). However, treatment is intense, and the overall benefit is small, necessitating effective biomarkers to identify patients who will benefit most.

OBJECTIVE

To characterize cell-free DNA (cfDNA) methylation in patients receiving NAC in SWOG S1314, a prospective cooperative group trial, and to correlate the methylation signatures with pathologic response at radical cystectomy.

DESIGN, SETTING, AND PARTICIPANTS

SWOG S1314 is a prospective cooperative group trial for patients with MIBC (cT2-T4aN0M0, ≥5 mm of viable tumor), with a primary objective of evaluating the coexpression extrapolation (COXEN) gene expression signature as a predictor of NAC response, defined as achieving pT0N0 or ≤pT1N0 at radical cystectomy. For the current exploratory analysis, blood samples were collected prospectively from 72 patients in S1314 before and during NAC, and plasma cfDNA methylation was measured using the Infinium MethylationEPIC BeadChip array.

INTERVENTION

No additional interventions besides plasma collection.

OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS

Differential methylation between pathologic responders (≤pT1N0) and nonresponders was analyzed, and a classifier predictive of treatment response was generated using the Random Forest machine learning algorithm.

RESULTS AND LIMITATIONS

Using prechemotherapy plasma cfDNA, we developed a methylation-based response score (mR-score) predictive of pathologic response. Plasma samples collected after the first cycle of NAC yielded mR-scores with similar predictive ability. Furthermore, we used cfDNA methylation data to calculate the circulating bladder DNA fraction, which had a modest but independent predictive ability for treatment response. In a model combining mR-score and circulating bladder DNA fraction, we correctly predicted pathologic response in 79% of patients based on their plasma collected at baseline and after one cycle of chemotherapy. Limitations of this study included a limited sample size and relatively low circulating bladder DNA levels.

CONCLUSIONS

Our study provides the proof of concept that cfDNA methylation can be used to generate classifiers of NAC response in bladder cancer patients.

PATIENT SUMMARY

In this exploratory analysis of S1314, we demonstrated that cell-free DNA methylation can be profiled to generate biomarker signatures associated with neoadjuvant chemotherapy response. With validation in additional cohorts, this minimally invasive approach may be used to predict chemotherapy response in locally advanced bladder cancer and perhaps also in metastatic disease.

Parenting Practices May Buffer the Impact of Adversity on Epigenetic Age Acceleration Among Young Children With Developmental Delays

This study examined whether children exposed to adversity would exhibit lower epigenetic age acceleration in the context of improved parenting. Children with developmental delays and externalizing behavior problems (N = 62; Mage = 36.26 months; 70.97% boys, 29.03% girls; 71% Latinx, 22.6% Black) were drawn from a larger randomized controlled trial (RCT), which randomized them to receive Internet-delivered parent-child interaction therapy (iPCIT; n = 30) or community referrals as usual (RAU; n = 32). Epigenetic age acceleration was estimated with the pediatric buccal epigenetic clock, using saliva. Adversity was assessed using parent, family, and neighborhood-level cumulative-risk indicators. Adversity interacted with Time 2 (T2) observations of positive and negative-parenting practices to predict epigenetic age acceleration 1.5 years later, regardless of treatment assignment. Children exposed to more adversity displayed lower epigenetic age acceleration when parents evidenced increased positive (b = -0.15, p = .001) and decreased negative (b = -0.12, p = .01) parenting practices.

Alpha-ketoglutarate supplementation and BiologicaL agE in middle-aged adults (ABLE)-intervention study protocol

Targeting molecular processes of aging will enable people to live healthier and longer lives by preventing age-related diseases. Geroprotectors are compounds with the potential to increase healthspan and lifespan. Even though many of them have been tested in animal models, the translation to humans is limited. Alpha-Ketoglutarate (AKG) has been studied widely in model animals, but there are few studies testing its geroprotective properties in humans. ABLE is a double blinded placebo-controlled randomized trial (RCT) of 1 g sustained release Ca-AKG versus placebo for 6 months of intervention and 3 months follow up including 120 40-60-year-old healthy individuals with a higher DNA methylation age compared to their chronological age. The primary outcome is the decrease in DNA methylation age from baseline to the end of the intervention. A total of 120 participants will be randomized to receive either sustained release Ca-AKG or placebo. Secondary outcomes include changes in the inflammatory and metabolic parameters in blood, handgrip strength and leg extension strength, arterial stiffness, skin autofluorescence, and aerobic capacity from baseline to 3 months, 6 months, and 9 months. This study will recruit middle-aged participants with an older DNA methylation age compared to their chronological age, and test whether supplementation with Ca-AKG can reduce DNA methylation age. This study is unique in its inclusion of biologically older participants.

Fresh and frozen cardiac tissue are comparable in DNA methylation array β-values, but formalin-fixed, paraffin-embedded tissue may overestimate DNA methylation levels

Untreated fresh cardiac tissue is the optimal tissue material for investigating DNA methylation patterns of cardiac biology and diseases. However, fresh tissue is difficult to obtain. Therefore, tissue stored as frozen or formalin-fixed, paraffin-embedded (FFPE) is widely used for DNA methylation studies. It is unknown whether storage conditions alter the DNA methylation in cardiac tissue. In this study, we compared the DNA methylation patterns of fresh, frozen, and FFPE cardiac tissue to investigate if the storage method affected the DNA methylation results. We used the Infinium MethylationEPIC assay to obtain genome-wide methylation levels in fresh, frozen, and FFPE tissues from nine individuals. We found that the DNA methylation levels of 21.4% of the examined CpG sites were overestimated in the FFPE samples compared to that of fresh and frozen tissue, whereas 5.7% were underestimated. Duplicate analyses of the DNA methylation patterns showed high reproducibility (precision) for frozen and FFPE tissues. In conclusion, we found that frozen and FFPE tissues gave reproducible DNA methylation results and that frozen and fresh tissues gave similar results.

Comprehensive Evaluation of The Infinium Human MethylationEPIC v2 BeadChip

Infinium Methylation BeadChips are widely used to profile DNA cytosine modifications in large cohort studies for reasons of cost-effectiveness, accurate quantification, and user-friendly data analysis in characterizing these canonical epigenetic marks. In this work, we conducted a comprehensive evaluation of the updated Infinium MethylationEPIC v2 BeadChip (EPICv2). Our evaluation revealed that EPICv2 offers significant improvements over its predecessors, including expanded enhancer coverage, applicability to diverse ancestry groups, support for low-input DNA down to one nanogram, coverage of existing epigenetic clocks, cell type deconvolution panels, and human trait associations, while maintaining accuracy and reproducibility. Using EPICv2, we were able to identify epigenome and sequence signatures in cell line models of DNMT and SETD2 loss and/or hypomorphism. Furthermore, we provided probe-wise evaluation and annotation to facilitate the use of new features on this array for studying the interplay between somatic mutations and epigenetic landscape in cancer genomics. In conclusion, EPICv2 provides researchers with a valuable tool for studying epigenetic modifications and their role in development and disease.

Epigenome-wide DNA methylation profiling in comparison between pathological and physiological hypertrophy of human cardiomyocytes

Background: Physiological and pathological stimuli result in distinct forms of cardiac hypertrophy, but the molecular regulation comparing the two, especially at the DNA methylation level, is not well understood. Methods: We conducted an in vitro study using human cardiomyocytes exposed to angiotensin II (AngII) and insulin-like growth factor 1 (IGF-1) to mimic pathologically and physiologically hypertrophic heart models, respectively. Whole genome DNA methylation patterns were profiled by the Infinium human MethylationEPIC platform with >850 K DNA methylation loci. Two external datasets were used for comparisons and qRT-PCR was performed for examining expression of associated genes of those identified DNA methylation loci. Results: We detected 194 loci that are significantly differentially methylated after AngII treatment, and 206 significant loci after IGF-1 treatment. Mapping the significant loci to genes, we identified 158 genes corresponding to AngII treatment and 175 genes to IGF-1 treatment. Using the gene-set enrichment analysis, the PI3K-Akt signaling pathway was identified to be significantly enriched for both AngII and IGF-1 treatment. The Hippo signaling pathway was enriched after IGF-1 treatment, but not for AngII treatment. CDK6 and RPTOR are components of the PI3K-Akt pathway but have different DNA methylation patterns in response to AngII and IGF-1. qRT-PCR confirmed the different gene expressions of CDK6 and PRTOR. Conclusion: Our study is pioneering in profiling epigenome DNA methylation changes in adult human cardiomyocytes under distinct stress conditions: pathological (AngII) and physiological (IGF-1). The identified DNA methylation loci, genes, and pathways might have the potential to distinguish between pathological and physiological cardiac hypertrophy.