Genome wide inherited modifications of the tomato epigenome by trans-activated bacterial CG methyltransferase.

BACKGROUND: Epigenetic variation is mediated by epigenetic marks such as DNA methylation occurring in all cytosine contexts in plants. CG methylation plays a critical role in silencing transposable elements and regulating gene expression. The establishment of CG methylation occurs via the RNA-directed DNA methylation pathway and CG methylation maintenance relies on METHYLTRANSFERASE1, the homologue of the mammalian DNMT1. PURPOSE: Here, we examined the capacity to stably alter the tomato genome methylome by a bacterial CG-specific M.SssI methyltransferase expressed through the LhG4/pOP transactivation system. RESULTS: Methylome analysis of M.SssI expressing plants revealed that their euchromatic genome regions are specifically hypermethylated in the CG context, and so are most of their genes. However, changes in gene expression were observed only with a set of genes exhibiting a greater susceptibility to CG hypermethylation near their transcription start site. Unlike gene rich genomic regions, our analysis revealed that heterochromatic regions are slightly hypomethylated at CGs only. Notably, some M.SssI-induced hypermethylation persisted even without the methylase or transgenes, indicating inheritable epigenetic modification. CONCLUSION: Collectively our findings suggest that heterologous expression of M.SssI can create new inherited epigenetic variations and changes in the methylation profiles on a genome wide scale. This open avenues for the conception of epigenetic recombinant inbred line populations with the potential to unveil agriculturally valuable tomato epialleles.

Fusobacterium nucleatum infection modulates the transcriptome and epigenome of HCT116 colorectal cancer cells in an oxygen-dependent manner.

Fusobacterium nucleatum, a gram-negative oral bacterium, has been consistently validated as a strong contributor to the progression of several types of cancer, including colorectal (CRC) and pancreatic cancer. While previous in vitro studies have shown that intracellular F. nucleatum enhances malignant phenotypes such as cell migration, the dependence of this regulation on features of the tumor microenvironment (TME) such as oxygen levels are wholly uncharacterized. Here we examine the influence of hypoxia in facilitating F. nucleatum invasion and its effects on host responses focusing on changes in the global epigenome and transcriptome. Using a multiomic approach, we analyze epigenomic alterations of H3K27ac and global transcriptomic alterations sustained within a hypoxia and normoxia conditioned CRC cell line HCT116 at 24 h following initial infection with F. nucleatum. Our findings reveal that intracellular F. nucleatum activates signaling pathways and biological processes in host cells similar to those induced upon hypoxia conditioning in the absence of infection. Furthermore, we show that a hypoxic TME favors F. nucleatum invasion and persistence and therefore infection under hypoxia may amplify malignant transformation by exacerbating the effects induced by hypoxia alone. These results motivate future studies to investigate host-microbe interactions in tumor tissue relevant conditions that more accurately define parameters for targeted cancer therapies.

Kernel-based testing for single-cell differential analysis.

Single-cell technologies offer insights into molecular feature distributions, but comparing them poses challenges. We propose a kernel-testing framework for non-linear cell-wise distribution comparison, analyzing gene expression and epigenomic modifications. Our method allows feature-wise and global transcriptome/epigenome comparisons, revealing cell population heterogeneities. Using a classifier based on embedding variability, we identify transitions in cell states, overcoming limitations of traditional single-cell analysis. Applied to single-cell ChIP-Seq data, our approach identifies untreated breast cancer cells with an epigenomic profile resembling persister cells. This demonstrates the effectiveness of kernel testing in uncovering subtle population variations that might be missed by other methods.

Multimodal analysis of cfDNA methylomes for early detecting esophageal squamous cell carcinoma and precancerous lesions.

Detecting early-stage esophageal squamous cell carcinoma (ESCC) and precancerous lesions is critical for improving survival. Here, we conduct whole-genome bisulfite sequencing (WGBS) on 460 cfDNA samples from patients with non-metastatic ESCC or precancerous lesions and matched healthy controls. We develop an expanded multimodal analysis (EMMA) framework to simultaneously identify cfDNA methylation, copy number variants (CNVs), and fragmentation markers in cfDNA WGBS data. cfDNA methylation markers are the earliest and most sensitive, detectable in 70% of ESCCs and 50% of precancerous lesions, and associated with molecular subtypes and tumor microenvironments. CNVs and fragmentation features show high specificity but are linked to late-stage disease. EMMA significantly improves detection rates, increasing AUCs from 0.90 to 0.99, and detects 87% of ESCCs and 62% of precancerous lesions with >95% specificity in validation cohorts. Our findings demonstrate the potential of multimodal analysis of cfDNA methylome for early detection and monitoring of molecular characteristics in ESCC.

Temporal dynamics of the multi-omic response to endurance exercise training.

Regular exercise promotes whole-body health and prevents disease, but the underlying molecular mechanisms are incompletely understood1-3. Here, the Molecular Transducers of Physical Activity Consortium4 profiled the temporal transcriptome, proteome, metabolome, lipidome, phosphoproteome, acetylproteome, ubiquitylproteome, epigenome and immunome in whole blood, plasma and 18 solid tissues in male and female Rattus norvegicus over eight weeks of endurance exercise training. The resulting data compendium encompasses 9,466 assays across 19 tissues, 25 molecular platforms and 4 training time points. Thousands of shared and tissue-specific molecular alterations were identified, with sex differences found in multiple tissues. Temporal multi-omic and multi-tissue analyses revealed expansive biological insights into the adaptive responses to endurance training, including widespread regulation of immune, metabolic, stress response and mitochondrial pathways. Many changes were relevant to human health, including non-alcoholic fatty liver disease, inflammatory bowel disease, cardiovascular health and tissue injury and recovery. The data and analyses presented in this study will serve as valuable resources for understanding and exploring the multi-tissue molecular effects of endurance training and are provided in a public repository ( https://motrpac-data.org/ ).

Blood-based epigenome-wide analyses of chronic low-grade inflammation across diverse population cohorts.

Chronic inflammation is a hallmark of age-related disease states. The effectiveness of inflammatory proteins including C-reactive protein (CRP) in assessing long-term inflammation is hindered by their phasic nature. DNA methylation (DNAm) signatures of CRP may act as more reliable markers of chronic inflammation. We show that inter-individual differences in DNAm capture 50% of the variance in circulating CRP (N = 17,936, Generation Scotland). We develop a series of DNAm predictors of CRP using state-of-the-art algorithms. An elastic-net-regression-based predictor outperformed competing methods and explained 18% of phenotypic variance in the Lothian Birth Cohort of 1936 (LBC1936) cohort, doubling that of existing DNAm predictors. DNAm predictors performed comparably in four additional test cohorts (Avon Longitudinal Study of Parents and Children, Health for Life in Singapore, Southall and Brent Revisited, and LBC1921), including for individuals of diverse genetic ancestry and different age groups. The best-performing predictor surpassed assay-measured CRP and a genetic score in its associations with 26 health outcomes. Our findings forge new avenues for assessing chronic low-grade inflammation in diverse populations.

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.

[Postnatal epigenome-mediated aging control and global trends].

The epigenome can adequately regulate the on/off states of genes in response to external environmental factors and stress. In recent years, it has been observed that the epigenome, which is modulated through DNA methylation, histone modifications, and chromatin remodeling, changes with age. Alterations in the epigenome lead to the loss of cell-specific epigenome/identity, which in turn triggers a decline in tissue function. In mammals, postnatal epigenomic variations are not only caused by metabolic diseases, such as diabetes or DNA damage, but also by social stress and infectious diseases. Unlike Genome-Wide Association Studies (GWAS), dynamically changing epigenomes, along with their cellular roles, need to be established as objective biomarkers in conjunction with various biological signals, such as walking speed, brain waves, and clinical data. The biological age/aging clock, determined by methylated DNA, has attracted attention, and calorie restriction not only slows the progression of aging, but also seems to suppress it. However, as indicated by gene expression analysis in aging mice, aging is not a linear model, but is represented by nonlinear dynamic changes. Consequently, the development of experimental models and analytical methods that enhance temporal resolution through time-series analysis, tailored to spatial resolution, such as cell distribution and organ specificity, is progressing. Moreover, in recent years, in addition to anti-aging efforts targeting epigenomic variations, global attention has increasingly focused on research and development aimed at rejuvenating treatments, thus leading to the birth of many biotech companies. Aging Hallmarks such as inflammation, stem cells, metabolism, genomic instability, and autophagy, interact closely with the epigenome. Various postnatal and reversible epigenomic controls of aging, including Yamanaka factors (OKSM and OSK), are now entering a new phase. In the future, the development of aging control using diverse modalities, such as mRNA, artificial peptides, and genome editing, is expected, along with an improved molecular understanding of aging and identification of useful biomarkers.

Measuring cancer driving force of chromosomal aberrations through multi-layer Boolean implication networks.

Multi-layer Complex networks are commonly used for modeling and analysing biological entities. This paper presents the advantage of using COMBO (Combining Multi Bio Omics) to suggest a new role of the chromosomal aberration as a cancer driver factor. Exploiting the heterogeneous multi-layer networks, COMBO integrates gene expression and DNA-methylation data in order to identify complex bilateral relationships between transcriptome and epigenome. We evaluated the multi-layer networks generated by COMBO on different TCGA cancer datasets (COAD, BLCA, BRCA, CESC, STAD) focusing on the effect of a specific chromosomal numerical aberration, broad gain in chromosome 20, on different cancer histotypes. In addition, the effect of chromosome 8q amplification was tested in the same TCGA cancer dataset. The results demonstrate the ability of COMBO to identify the chromosome 20 amplification cancer driver force in the different TCGA Pan Cancer project datasets.

Genome-wide profiling of DNA methylome and transcriptome reveals epigenetic regulation of Urechis unicinctus response to sulfide stress.

Sulfide is a well-known environmental pollutant that can have detrimental effects on most organisms. However, few metazoans living in sulfide-rich environments have developed mechanisms to tolerate and adapt to sulfide stress. Epigenetic mechanisms, including DNA methylation, have been shown to play a vital role in environmental stress adaptation. Nevertheless, the precise function of DNA methylation in biological sulfide adaptation remains unclear. Urechis unicinctus, a benthic organism inhabiting sulfide-rich intertidal environments, is an ideal model organism for studying adaptation to sulfide environments. In this study, we conducted a comprehensive analysis of the DNA methylome and transcriptome of U. unicinctus after exposure to 50 µM sulfide. The results revealed dynamic changes in the DNA methylation (5-methylcytosine) landscape in response to sulfide stress, with U. unicinctus exhibiting elevated DNA methylation levels following stress exposure. Integrating differentially expressed genes (DEGs) and differentially methylated regions (DMRs), we identified a crucial role of gene body methylation in predicting gene expression. Furthermore, using a DNA methyltransferase inhibitor, we validated the involvement of DNA methylation in the sulfide stress response and the gene regulatory network influenced by DNA methylation. The results indicated that by modulating DNA methylation levels during sulfide stress, the expression of glutathione S-transferase, glutamyl aminopeptidase, and cytochrome c oxidase could be up-regulated, thereby facilitating the metabolism and detoxification of exogenous sulfides. Moreover, DNA methylation was found to regulate and enhance the oxidative phosphorylation pathway, including NADH dehydrogenase, isocitrate dehydrogenase, and ATP synthase. Additionally, DNA methylation influenced the regulation of Cytochrome P450 and macrophage migration inhibitory factor, both of which are closely associated with oxidative stress and stress resistance. Our findings not only emphasize the role of DNA methylation in sulfide adaptation but also provide novel insights into the potential mechanisms through which marine organisms adapt to environmental changes.

DNA Methylation in Ctenophores.

Epigenomic regulation and dynamic DNA methylation, in particular, are widespread mechanisms orchestrating the genome operation across time and species. Whole-genome bisulfite sequencing (WGBS) is currently the only method for unbiasedly capturing the presence of 5-methylcytosine (5-mC) DNA methylation patterns across an entire genome with single-nucleotide resolution. Bisulfite treatment converts unmethylated cytosines to uracils but leaves methylated cytosines intact, thereby creating a map of all methylated cytosines across a genome also known as a methylome. These epigenomic patterns of DNA methylation have been found to regulate gene expression and influence gene evolution rates between species. While protocols have been optimized for vertebrate methylome production, little adaptation has been done for invertebrates. Creating a methylome reference allows comparisons to be made between rates of transcription and epigenomic patterning in animals. Here we present a method of library construction for bisulfite sequencing optimized for non-bilateral metazoans such as the ctenophore, Mnemiopsis leidyi. We have improved upon our previously published method by including spike-in genomic DNA controls to measure methylation conversion rates. By pooling two bisulfite conversion reactions from the same individual, we also produced sequencing libraries that yielded a higher percentage of sequenced reads uniquely mapping to the reference genome. We successfully detected 5-mC in whole-animal methylomes at CpG, CHG, and CHH sites and visualized datasets using circos diagrams. The proof-of-concept tests were performed both under control conditions and following injury tests with changes in methylation patterns of genes encoding innexins, toxins and neuropeptides. Our approach can be easily adapted to produce epigenomes from other fragile marine animals.

Integrative analysis in head and neck cancer reveals distinct role of miRNome and methylome as tumour epigenetic drivers.

Head and neck cancer is the sixth most common malignancy worldwide, with the relatively low 5-year survival rate, mainly because it is diagnosed at a late stage. Infection with HPV is a well known aetiology, which affects the nature of these cancers and patients' survival. Besides, it is considered that the main driving force for this type of cancer could be epigenetics. In this study we aimed to find potential epigenetic biomarkers, by integrating miRNome, methylome, and transcriptome analyses. From the fresh head and neck cancer tissue samples, we chose a group for miRNome, methylome and transcriptome profiling, in comparison to adequate control samples. Bioinformatics analyses are performed in R v4.2.2. Count normalisation and group differential expression for mRNA and the previously obtained miRNA count data was performed with DESeq2 v1.36. Gene set enrichment analysis was performed and visualised using gProfiler2 v0.2.1 Identification of miRNA targets was performed by querying in miRTarBase using multiMiR v1.18.0. Annotation of CpG sites merging into islands was obtained from RnBeads.hg19 v1.28.0. package. For the integrative analysis we performed kmeans clustering using stats v4.2.2 package, using 8-12 clusters and nstart 100. We found that transcriptome analysis divides samples into cancers and controls clusters, with no relation to HPV status or cancer anatomical location. Differentially expressed genes (n = 2781) were predominantly associated with signalling pathways of tumour progression. We identified a cluster of genes under the control of the transcription factor E2F that are significantly underexpressed in cancer tissue, as well as T cell immunity genes and genes related to regulation of transcription. Among overexpressed genes in tumours we found those that belong to cell cycle regulation and vasculature. A small number of genes were found significantly differentially expressed in HPV-positive versus HPV-negative tumours (for example NEFH, ZFR2, TAF7L, ZNF541, and TYMS). In this comprehensive study on an overlapping set of samples where the integration of miRNome, methylome and transcriptome analysis were performed for head and neck cancer, we demonstrated that the majority of genes were associated exclusively with miRNome or methylome and, to a lesser extent, under the control of both epigenetic mechanisms.

Genome- and epigenome-wide association studies identify susceptibility of CpG sites and regions for metabolic syndrome in a Korean population.

BACKGROUND: While multiple studies have investigated the relationship between metabolic syndrome (MetS) and its related traits (fasting glucose, triglyceride, HDL cholesterol, blood pressure, waist circumference) and DNA methylation, our understanding of the epigenetic mechanisms in MetS remains limited. Therefore, we performed an epigenome-wide meta-analysis of blood DNA methylation to identify differentially methylated probes (DMPs) and differentially methylated regions (DMRs) associated with MetS and its components using two independent cohorts comprising a total of 2,334 participants. We also investigated the specific genetic effects on DNA methylation, identified methylation quantitative trait loci (meQTLs) through genome-wide association studies and further utilized Mendelian randomization (MR) to assess how these meQTLs subsequently influence MetS status. RESULTS: We identified 40 DMPs and 27 DMRs that are significantly associated with MetS. In addition, we identified many novel DMPs and DMRs underlying inflammatory and steroid hormonal processes. The most significant associations were observed in 3 DMPs (cg19693031, cg26974062, cg02988288) and a DMR (chr1:145440444-145441553) at the TXNIP, which are involved in lipid metabolism. These CpG sites were identified as coregulators of DNA methylation in MetS, TG and FAG levels. We identified a total of 144 cis-meQTLs, out of which only 13 were found to be associated with DMPs for MetS. Among these, we confirmed the identified causal mediators of genetic effects at CpG sites cg01881899 at ABCG1 and cg00021659 at the TANK genes for MetS. CONCLUSIONS: This study observed whether specific CpGs and methylated regions act independently or are influenced by genetic effects for MetS and its components in the Korean population. These associations between the identified DNA methylation and MetS, along with its individual components, may serve as promising targets for the development of preventive interventions for MetS.

The Interplay between Endocrine-Disrupting Chemicals and the Epigenome towards Metabolic Dysfunction-Associated Steatotic Liver Disease: A Comprehensive Review.

Metabolic dysfunction-associated steatotic liver disease (MASLD), described as the most prominent cause of chronic liver disease worldwide, has emerged as a significant public health issue, posing a considerable challenge for most countries. Endocrine-disrupting chemicals (EDCs), commonly found in daily use items and foods, are able to interfere with nuclear receptors (NRs) and disturb hormonal signaling and mitochondrial function, leading, among other metabolic disorders, to MASLD. EDCs have also been proposed to cause transgenerationally inherited alterations leading to increased disease susceptibility. In this review, we are focusing on the most prominent linking pathways between EDCs and MASLD, their role in the induction of epigenetic transgenerational inheritance of the disease as well as up-to-date practices aimed at reducing their impact.

From morphology to methylome: epigenetic studies of Müllerian mesonephric-like adenocarcinoma reveal similarities to cervical mesonephric adenocarcinoma†.

Mesonephric adenocarcinomas (MAs) and mesonephric-like adenocarcinomas (MLAs) are rare, aggressive neoplasms that arise in the gynecologic tract and show overlapping morphologic, immunohistochemical, and molecular features. While MAs occur in the cervix and are thought to arise from mesonephric remnants, MLAs occur in the endometrium and ovary and are believed to originate from transdifferentiation of Müllerian lesions. Both MAs and MLAs show a variety of architectural patterns, exhibit frequent expression of GATA3 by immunohistochemistry, and harbor KRAS mutations. In a recent article published in The Journal of Pathology, Kommoss and colleagues used DNA methylation profiling to extend these similarities and showed that MLAs and MAs cluster together based on their epigenetic signatures and are epigenetically distinct from other Müllerian adenocarcinomas. They also showed that MLAs and MAs harbor a high number of global copy number alterations. This study provides evidence that MLAs more closely resemble MAs than Müllerian carcinomas on an epigenetic level. As a result, the authors argue that MLA should be renamed 'mesonephric-type adenocarcinoma.' Further research is needed to establish the relationship between these two entities, their etiology, and pathogenesis. © 2024 The Pathological Society of Great Britain and Ireland.

Cortico-striatal differences in the epigenome in attention-deficit/ hyperactivity disorder.

While epigenetic modifications have been implicated in ADHD through studies of peripheral tissue, to date there has been no examination of the epigenome of the brain in the disorder. To address this gap, we mapped the methylome of the caudate nucleus and anterior cingulate cortex in post-mortem tissue from fifty-eight individuals with or without ADHD. While no single probe showed adjusted significance in differential methylation, several differentially methylated regions emerged. These regions implicated genes involved in developmental processes including neurogenesis and the differentiation of oligodendrocytes and glial cells. We demonstrate a significant association between differentially methylated genes in the caudate and genes implicated by GWAS not only in ADHD but also in autistic spectrum, obsessive compulsive and bipolar affective disorders through GWAS. Using transcriptomic data available on the same subjects, we found modest correlations between the methylation and expression of genes. In conclusion, this study of the cortico-striatal methylome points to gene and gene pathways involved in neurodevelopment, consistent with studies of common and rare genetic variation, as well as the post-mortem transcriptome in ADHD.

Alterations in Blood Methylome as Potential Epigenetic Biomarker in Sporadic Parkinson's Disease.

OBJECTIVE: To characterize DNA methylation (DNAm) differences between sporadic Parkinson's disease (PD) and healthy control (HC) individuals enrolled in the Parkinson's Progression Markers Initiative (PPMI). METHODS: Using whole blood, we characterized longitudinal differences in DNAm between sporadic PD patients (n = 196) and HCs (n = 86) enrolled in PPMI. RNA sequencing (RNAseq) was used to conduct gene expression analyses for genes mapped to differentially methylated cytosine-guanine sites (CpGs). RESULTS: At the time of patient enrollment, 5,178 CpGs were differentially methylated (2,683 hypermethylated and 2,495 hypomethylated) in PD compared to HC. Of these, 579 CpGs underwent significant methylation changes over 3 years. Several differentially methylated CpGs were found near the cytochrome P450 family 2 subfamily E member 1 (CYP2E1) gene. Additionally, multiple hypermethylated CpGs were associated with the N-myc downregulated gene family member 4 (NDRG4) gene. RNA-Seq analyses showed 75 differentially expressed genes in PD patients compared to controls. An integrative analysis of both differentially methylated sites and differentially expressed genes revealed 20 genes that exhibited hypomethylation concomitant with overexpression. Additionally, 1 gene, cathepsin H (CTSH), displayed hypermethylation that was associated with its decreased expression. INTERPRETATION: We provide initial evidence of alterations in DNAm in blood of PD patients that may serve as potential epigenetic biomarker of disease. To evaluate the significance of these changes throughout the progression of PD, additional profiling at longer intervals and during the prodromal stages of disease will be necessary. ANN NEUROL 2024;95:1162-1172.

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.