Charting a dynamic DNA methylation landscape of the human genome

Implementation of the Methyl-Seq platform to identify tissue- and sex-specific DNA methylation differences in the rat epigenome

Epigenomic annotations for the rat lag far behind those of human and mouse, despite the rat's immense utility in pharmacological and behavioral studies and the need to understand their epigenetic mechanisms. We have designed a targeted-enrichment method followed by next-generation sequencing (Methyl-Seq) to identify DNA methylation (DNAm) signatures across the rat genome. The design reflected an attempt to create a more comprehensive investigation of the rat epigenome, as it included promoters, CpG islands, and island shores of all RefSeq genes. In this study, we implemented the rat Methyl-Seq platform and tested its ability to distinguish differentially methylated regions (DMRs) among three different tissue types, three distinct brain regions, and, in the hippocampus, between males and females. These comparisons yielded DNAm differences of differing magnitudes, many of which were independently validated by bisulfite pyrosequencing, including autosomal regions that were predicted to show the least degree of difference in DNAm between males and females. Quantitative reverse transcription PCR revealed that most genes associated with the DMRs showed tissue-, brain region-, and sex-specific differences in expression. In particular, we found evidence for sex-specific DNAm and expression differences at Tubb6, Lrrn2, Tex26, and Sox5l1, all of which play important roles in neurodevelopment and have been implicated in studies examining sex differences. Our results demonstrate the utility of the rat Methyl-Seq platform and suggest the presence of DNAm differences between the male and female hippocampus. The rat Methyl-Seq has the potential to provide epigenomic insights into pharmacological and behavioral studies performed in the rat.

Methylation and transcriptomic profiling reveals short term and long term regulatory responses in polarized macrophages

Macrophage plasticity allows the adoption of distinct functional states in response to environmental cues. While unique transcriptomic profiles define these states, focusing solely on transcription neglects potential long-term effects. The investigation of epigenetic changes can be used to understand how temporary stimuli can result in lasting effects. Epigenetic alterations play an important role in the pathophysiology of macrophages, including their trained innate immunity, enabling faster and more efficient inflammatory responses upon subsequent encounters to the same pathogen or insult. In this study, we used a multi-omics approach to elucidate the interplay between gene expression and DNA-methylation, to explore the potential long-term effects of diverse polarizing environments on macrophage activity. We identified a common core set of genes that are differentially methylated regardless of exposure type, indicating a potential common fundamental mechanism for adaptation to various stimuli. Functional analysis revealed that processes requiring rapid responses displayed transcriptomic regulation, whereas functions critical for long-term adaptations exhibited co-regulation at both transcriptomic and epigenetic levels. Our study uncovers a novel set of genes linked to the long-term effects of macrophage polarization. This discovery underscores the potential of epigenetics in elucidating how macrophages establish long-term memory and influence health outcomes.

Diet Quality and Epigenetic Aging in the Women's Health Initiative

BACKGROUND

Higher diet quality scores are associated with a lower risk for many chronic diseases and all-cause mortality; however, it is unclear if diet quality is associated with aging biology.

OBJECTIVE

This study aimed to examine the association between diet quality and a measure of biological aging known as epigenetic aging.

DESIGN

A cross-sectional data analysis was used to examine the association between three diet quality scores based on self-reported food frequency questionnaire data and five measures of epigenetic aging based on DNA methylation (DNAm) data from peripheral blood.

PARTICIPANTS/SETTING

This study included 4,500 postmenopausal women recruited from multiple sites across the United States (1993-98), aged 50 to 79 years, with food frequency questionnaire and DNAm data available from the Women's Health Initiative baseline visit.

MAIN OUTCOME MEASURES

Five established epigenetic aging measures were generated from HumanMethylation450 Beadchip DNAm data, including AgeAccelHannum, AgeAccelHorvath, AgeAccelPheno, AgeAccelGrim, and DunedinPACE.

STATISTICAL ANALYSES PERFORMED

Linear mixed models were used to test for associations between three diet quality scores (Healthy Eating Index, Dietary Approaches to Stop Hypertension, and alternate Mediterranean diet scores) and epigenetic aging measures, adjusted for age, race and ethnicity, education, tobacco smoking, physical activity, Women's Health Initiative substudy from which DNAm data were obtained, and DNAm-based estimates of leukocyte proportions.

RESULTS

Healthy Eating Index, Dietary Approaches to Stop Hypertension, and alternate Mediterranean diet scores were all inversely associated with AgeAccelPheno, AgeAccelGrim, and DunedinPACE (P < 0.05), with the largest effects with DunedinPACE. A one standard deviation increment in diet quality scores was associated with a decrement (β ± SE) in DunedinPACE z score of -0.097 ± 0.014 (P = 9.70 x 10-13) for Healthy Eating Index, -0.107 ± 0.014 (P = 1.53 x 10-14) for Dietary Approaches to Stop Hypertension, and -0.068 ± 0.013 (P = 2.31 x 10-07) for the alternate Mediterranean diet.

CONCLUSIONS

In postmenopausal women, diet quality scores were inversely associated with DNAm-based measures of biological aging, particularly DunedinPACE.

DNA methylation of transposons pattern aging differences across a diverse cohort of dogs from the Dog Aging Project

Within a species, larger individuals often have shorter lives and higher rates of age-related disease. Despite this well-known link, we still know little about underlying age-related epigenetic differences, which could help us better understand inter-individual variation in aging and the etiology, onset, and progression of age-associated disease. Dogs exhibit this negative correlation between size, health, and longevity and thus represent an excellent system in which to test the underlying mechanisms. Here, we quantified genome-wide DNA methylation in a cohort of 864 dogs in the Dog Aging Project. Age strongly patterned the dog epigenome, with the majority (66% of age-associated loci) of regions associating age-related loss of methylation. These age effects were non-randomly distributed in the genome and differed depending on genomic context. We found the LINE1 (long interspersed elements) class of TEs (transposable elements) were the most frequently hypomethylated with age (FDR < 0.05, 40% of all LINE1 regions). This LINE1 pattern differed in magnitude across breeds of different sizes- the largest dogs lost 0.26% more LINE1 methylation per year than the smallest dogs. This suggests that epigenetic regulation of TEs, particularly LINE1s, may contribute to accelerated age and disease phenotypes within a species. Since our study focused on the methylome of immune cells, we looked at LINE1 methylation changes in golden retrievers, a breed highly susceptible to hematopoietic cancers, and found they have accelerated age-related LINE1 hypomethylation compared to other breeds. We also found many of the LINE1s hypomethylated with age are located on the X chromosome and are, when considering X chromosome inactivation, counter-intuitively more methylated in males. These results have revealed the demethylation of LINE1 transposons as a potential driver of inter-species, demographic-dependent aging variation.

A comprehensive lifestyle index and its associations with DNA methylation and type 2 diabetes among Ghanaian adults: the rodam study

BACKGROUND

A series of modifiable lifestyle factors, such as diet quality, physical activity, alcohol intake, and smoking, may drive the rising burden of type 2 diabetes (T2DM) among sub-Saharan Africans globally. It is unclear whether epigenetic changes play a mediatory role in the associations between these lifestyle factors and T2DM. We assessed the associations between a comprehensive lifestyle index, DNA methylation and T2DM among Ghanaian adults.

METHODS

We used whole-blood Illumina 450 k DNA methylation data from 713 Ghanaians from the Research on Obesity and Diabetes among African Migrants (RODAM) study. We constructed a comprehensive lifestyle index based on established cut-offs for diet quality, physical activity, alcohol intake, and smoking status. In the T2DM-free discovery cohort (n = 457), linear models were fitted to identify differentially methylated positions (DMPs) and differentially methylated regions (DMRs) associated with the lifestyle index after adjustment for age, sex, body mass index (BMI), and technical covariates. Associations between the identified DMPs and the primary outcome (T2DM), as well as secondary outcomes (fasting blood glucose (FBG) and HbA1c), were determined via logistic and linear regression models, respectively.

RESULTS

In the present study population (mean age: 52 ± 10 years; male: 42.6%), the comprehensive lifestyle index showed a significant association with one DMP annotated to an intergenic region on chromosome 7 (false discovery rate (FDR) = 0.024). Others were annotated to ADCY7, SMARCE1, AHRR, LOXL2, and PTBP1 genes. One DMR was identified and annotated to the GFPT2 gene (familywise error rate (FWER) from bumphunter bootstrap = 0.036). None of the DMPs showed significant associations with T2DM; directions of effect were positive for the DMP in the AHRR and inverse for all the other DMPs. Higher methylation of the ADCY7 DMP was associated with higher FBG (p = 0.024); LOXL2 DMP was associated with lower FBG (p = 0.023) and HbA1c (p = 0.049); and PTBP1 DMP was associated with lower HbA1c (p = 0.002).

CONCLUSIONS

In this explorative epigenome-wide association study among Ghanaians, we identified one DMP and DMR associated with a comprehensive lifestyle index not previously associated with individual lifestyle factors. Based on our findings, we infer that lifestyle factors in combination, affect DNA methylation, thereby influencing the risk of T2DM among Ghanaian adults living in different contexts.

The role of DNA methylation in chondrogenesis of human iPSCs as a stable marker of cartilage quality

BACKGROUND

Lack of insight into factors that determine purity and quality of human iPSC (hiPSC)-derived neo-cartilage precludes applications of this powerful technology toward regenerative solutions in the clinical setting. Here, we set out to generate methylome-wide landscapes of hiPSC-derived neo-cartilages from different tissues-of-origin and integrated transcriptome-wide data to identify dissimilarities in set points of methylation with associated transcription and the respective pathways in which these genes act.

METHODS

We applied in vitro chondrogenesis using hiPSCs generated from two different tissue sources: skin fibroblasts and articular cartilage. Upon differentiation toward chondrocytes, these are referred to as hFiCs and hCiC, respectively. Genome-wide DNA methylation and RNA sequencing datasets were generated of the hiPSC-derived neo-cartilages, and the epigenetically regulated transcriptome was compared to that of neo-cartilage deposited by human primary articular cartilage (hPAC).

RESULTS

Methylome-wide landscapes of neo-cartilages of hiPSCs reprogrammed from two different somatic tissues were 85% similar to that of hPACs. By integration of transcriptome-wide data, differences in transcriptionally active CpGs between hCiC relative to hPAC were prioritized. Among the CpG-gene pairs lower expressed in hCiCs relative to hPACs, we identified genes such as MGP, GDF5, and CHAD enriched in closely related pathways and involved in cartilage development that likely mark phenotypic differences in chondrocyte states. Vice versa, among the CpG-gene pairs higher expressed, we identified genes such as KIF1A or NKX2-2 enriched in neurogenic pathways and likely reflecting off target differentiation.

CONCLUSIONS

We did not find significant variation between the neo-cartilages derived from hiPSCs of different tissue sources, suggesting that application of a robust differentiation protocol such as we applied here is more important as compared to the epigenetic memory of the cells of origin. Results of our study could be further exploited to improve quality, purity, and maturity of hiPSC-derived neo-cartilage matrix, ultimately to realize introduction of sustainable, hiPSC-derived neo-cartilage implantation into clinical practice.

Improved detection of methylation in ancient DNA

Reconstructing premortem DNA methylation levels in ancient DNA has led to breakthrough studies such as the prediction of anatomical features of the Denisovan. These studies rely on computationally inferring methylation levels from damage signals in naturally deaminated cytosines, which requires expensive high-coverage genomes. Here, we test two methods for direct methylation measurement developed for modern DNA based on either bisulfite or enzymatic methylation treatments. Bisulfite treatment shows the least reduction in DNA yields as well as the least biases during methylation conversion, demonstrating that this method can be successfully applied to ancient DNA.

Methylation patterns at the adjacent CpG sites within enhancers are a part of cell identity

BACKGROUND

It is generally accepted that methylation status of CpG sites spaced up to 50 bp apart is correlated, and accumulation of locally disordered methylation at adjacent CpG sites is involved in neoplastic transformation, acting in similar way as stochastic accumulation of mutations.

RESULTS

We used EPIC microarray data from 596 samples, representing 12 healthy tissue and cell types, as well as 572 blood cancer specimens to analyze methylation status of adjacent CpG sites across human genome, and subsequently validated our findings with NGS and Sanger sequencing. Our analysis showed that there is a subset of the adjacent CpG sites in human genome, with cytosine at one CpG site methylated and the other devoid of methyl group. These loci map to enhancers that are targeted by families of transcription factors involved in cell differentiation. Moreover, our results suggest that the methylation at these loci differ between alleles within a cell, what allows for remarkable level of heterogeneity of methylation patterns. However, different types of specialized cells acquire only one specific and stable pattern of methylation at each of these loci and that pattern is to a large extent lost during neoplastic transformation.

CONCLUSIONS

We identified a substantial number of adjacent CpG loci in human genome that display remarkably stable and cell type specific methylation pattern. The methylation pattern at these loci appears to reflect different methylation of alleles in cells. Furthermore, we showed that changes of methylation status at those loci are likely to be involved in regulation of the activity of enhancers and contribute to neoplastic transformation.

Mitochondrial DNA content and methylation in sperm of patients with asthenozoospermia

PURPOSE

The aim of the current study was to investigate the mtDNA methylation levels and mtDNA copy numbers in the sperm of patients with asthenozoospermia and compare them to those observed in controls with normozoospermia.

METHODS

Pyrosequencing analysis of the methylation levels of the mitochondrial D-loop and MT-CO1/chr1:631,907-632083/chrX:26,471,887-126,472,063 (hereinafter referred to as "MT-CO1-AVG") region and quantitative PCR analysis of the mtDNA copy number were performed on sperm from 30 patients with asthenozoospermia and 30 controls with normozoospermia.

RESULTS

Compared with those of controls with normozoospermia, the methylation levels of D-loop and MT-CO1-AVG regions and mtDNA copy number were significantly higher in patients with asthenozoospermia. The methylation level of the D-loop region in patients with asthenozoospermia and controls with normozoospermia and that of MT-CO1-AVG region in patients with asthenozoospermia showed a decreasing tendency with increasing total sperm motility. A significant inverse correlation between the mtDNA copy number and total sperm motility was observed in patients with asthenozoospermia but not in controls with normozoospermia. In patients with asthenozoospermia, but not in controls with normozoospermia, we observed a significant inverse correlation between D-loop methylation levels and mtDNA copy number, while no significant correlation was observed between MT-CO1-AVG methylation levels and mtDNA copy number.

CONCLUSION

These results reveal the occurrence of mtDNA methylation in human sperm and altered D-loop and MT-CO1-AVG methylation levels in patients with asthenozoospermia. Additional research is needed to determine the function of these features in the etiology and course of asthenozoospermia.

Fast matrix completion in epigenetic methylation studies with informative covariates

DNA methylation is an important epigenetic mark that modulates gene expression through the inhibition of transcriptional proteins binding to DNA. As in many other omics experiments, the issue of missing values is an important one, and appropriate imputation techniques are important in avoiding an unnecessary sample size reduction as well as to optimally leverage the information collected. We consider the case where relatively few samples are processed via an expensive high-density whole genome bisulfite sequencing (WGBS) strategy and a larger number of samples is processed using more affordable low-density, array-based technologies. In such cases, one can impute the low-coverage (array-based) methylation data using the high-density information provided by the WGBS samples. In this paper, we propose an efficient Linear Model of Coregionalisation with informative Covariates (LMCC) to predict missing values based on observed values and covariates. Our model assumes that at each site, the methylation vector of all samples is linked to the set of fixed factors (covariates) and a set of latent factors. Furthermore, we exploit the functional nature of the data and the spatial correlation across sites by assuming some Gaussian processes on the fixed and latent coefficient vectors, respectively. Our simulations show that the use of covariates can significantly improve the accuracy of imputed values, especially in cases where missing data contain some relevant information about the explanatory variable. We also showed that our proposed model is particularly efficient when the number of columns is much greater than the number of rows-which is usually the case in methylation data analysis. Finally, we apply and compare our proposed method with alternative approaches on two real methylation datasets, showing how covariates such as cell type, tissue type or age can enhance the accuracy of imputed values.

NF-κB and TET2 promote macrophage reprogramming in hypoxia that overrides the immunosuppressive effects of the tumor microenvironment

Macrophages orchestrate tissue homeostasis and immunity. In the tumor microenvironment (TME), macrophage presence is largely associated with poor prognosis because of their reprogramming into immunosuppressive cells. We investigated the effects of hypoxia, a TME-associated feature, on the functional, epigenetic, and transcriptional reprogramming of macrophages and found that hypoxia boosts their immunogenicity. Hypoxic inflammatory macrophages are characterized by a cluster of proinflammatory genes undergoing ten-eleven translocation-mediated DNA demethylation and overexpression. These genes are regulated by NF-κB, while HIF1α dominates the transcriptional reprogramming, demonstrated through ChIP-seq and pharmacological inhibition. In bladder and ovarian carcinomas, hypoxic inflammatory macrophages are enriched in immune-infiltrated tumors, correlating with better patient prognoses. Coculture assays and cell-cell communication analyses support that hypoxic-activated macrophages enhance T cell-mediated responses. The NF-κB-associated hypomethylation signature is displayed by a subset of hypoxic inflammatory macrophages, isolated from ovarian tumors. Our results challenge paradigms regarding the effects of hypoxia on macrophages and highlight actionable target cells to modulate anticancer immune responses.

A machine learning-based method for feature reduction of methylation data for the classification of cancer tissue origin

BACKGROUND

Genome DNA methylation profiling is a promising yet costly method for cancer classification, involving substantial data. We developed an ensemble learning model to identify cancer types using methylation profiles from a limited number of CpG sites.

METHODS

Analyzing methylation data from 890 samples across 10 cancer types from the TCGA database, we utilized ANOVA and Gain Ratio to select the most significant CpG sites, then employed Gradient Boosting to reduce these to just 100 sites.

RESULTS

This approach maintained high accuracy across multiple machine learning models, with classification accuracy rates between 87.7% and 93.5% for methods including Extreme Gradient Boosting, CatBoost, and Random Forest. This method effectively minimizes the number of features needed without losing performance, helping to classify primary organs and uncover subgroups within specific cancers like breast and lung.

CONCLUSIONS

Using a gradient boosting feature selector shows potential for streamlining methylation-based cancer classification.

Comparative DNA methylation reveals epigenetic adaptation to high altitude in snub-nosed monkeys

DNA methylation plays a crucial role in environmental adaptations. Here, using whole-genome bisulfite sequencing, we generated comprehensive genome-wide DNA methylation profiles for the high-altitude Yunnan snub-nosed monkey ( Rhinopithecus bieti) and the closely related golden snub-nosed monkey ( R. roxellana). Our findings indicated a slight increase in overall DNA methylation levels in golden snub-nosed monkeys compared to Yunnan snub-nosed monkeys, suggesting a higher prevalence of hypermethylated genomic regions in the former. Comparative genomic methylation analysis demonstrated that genes associated with differentially methylated regions were involved in membrane fusion, vesicular formation and trafficking, hemoglobin function, cell cycle regulation, and neuronal differentiation. These results suggest that the high-altitude-related epigenetic modifications are extensive, involving a complete adaptation process from the inhibition of single Ca 2+ channel proteins to multiple proteins collaboratively enhancing vesicular function or inhibiting cell differentiation and proliferation. Functional assays demonstrated that overexpression or down-regulation of candidate genes, such as SNX10, TIMELESS, and CACYBP, influenced cell viability under stress conditions. Overall, this research suggests that comparing DNA methylation across closely related species can identify novel candidate genomic regions and genes associated with local adaptations, thereby deepening our understanding of the mechanisms underlying environmental adaptations.

Mining nucleic acid "omics" to boost liquid biopsy in cancer

Treatments for cancer patients are becoming increasingly complex, and there is a growing desire from clinicians and patients for biomarkers that can account for this complexity to support informed decisions about clinical care. To achieve precision medicine, the new generation of biomarkers must reflect the spatial and temporal heterogeneity of cancer biology both between patients and within an individual patient. Mining the different layers of 'omics in a multi-modal way from a minimally invasive, easily repeatable, liquid biopsy has increasing potential in a range of clinical applications, and for improving our understanding of treatment response and resistance. Here, we detail the recent developments and methods allowing exploration of genomic, epigenomic, transcriptomic, and fragmentomic layers of 'omics from liquid biopsy, and their integration in a range of applications. We also consider the specific challenges that are posed by the clinical implementation of multi-omic liquid biopsies.

Cost-effective solutions for high-throughput enzymatic DNA methylation sequencing

Characterizing DNA methylation patterns is important for addressing key questions in evolutionary biology, geroscience, and medical genomics. While costs are decreasing, whole-genome DNA methylation profiling remains prohibitively expensive for most population-scale studies, creating a need for cost-effective, reduced representation approaches (i.e., assays that rely on microarrays, enzyme digests, or sequence capture to target a subset of the genome). Most common whole genome and reduced representation techniques rely on bisulfite conversion, which can damage DNA resulting in DNA loss and sequencing biases. Enzymatic methyl sequencing (EM-seq) was recently proposed to overcome these issues, but thorough benchmarking of EM-seq combined with cost-effective, reduced representation strategies has not yet been performed. To do so, we optimized Targeted Methylation Sequencing protocol (TMS)-which profiles ∼4 million CpG sites-for miniaturization, flexibility, and multispecies use at a cost of ∼$80. First, we tested modifications to increase throughput and reduce cost, including increasing multiplexing, decreasing DNA input, and using enzymatic rather than mechanical fragmentation to prepare DNA. Second, we compared our optimized TMS protocol to commonly used techniques, specifically the Infinium MethylationEPIC BeadChip (n=55 paired samples) and whole genome bisulfite sequencing (n=6 paired samples). In both cases, we found strong agreement between technologies (R² = 0.97 and 0.99, respectively). Third, we tested the optimized TMS protocol in three non-human primate species (rhesus macaques, geladas, and capuchins). We captured a high percentage (mean=77.1%) of targeted CpG sites and produced methylation level estimates that agreed with those generated from reduced representation bisulfite sequencing (R² = 0.98). Finally, we applied our protocol to profile age-associated DNA methylation variation in two subsistence-level populations-the Tsimane of lowland Bolivia and the Orang Asli of Peninsular Malaysia-and found age-methylation patterns that were strikingly similar to those reported in high income cohorts, despite known differences in age-health relationships between lifestyle contexts. Altogether, our optimized TMS protocol will enable cost-effective, population-scale studies of genome-wide DNA methylation levels across human and non-human primate species.

DNA methylation and histone modifications associated with antipsychotic treatment: a systematic review

Antipsychotic medications are essential when treating schizophrenia spectrum and other psychotic disorders, but the efficacy and tolerability of these medications vary from person to person. This interindividual variation is likely mediated, at least in part, by epigenomic processes that have yet to be fully elucidated. Herein, we systematically identified and evaluated 65 studies that examine the influence of antipsychotic drugs on epigenomic changes, including global methylation (9 studies), genome-wide methylation (22 studies), candidate gene methylation (16 studies), and histone modification (18 studies). Our evaluation revealed that haloperidol was consistently associated with increased global hypermethylation, which corroborates with genome-wide analyses, mostly performed by methylation arrays. In contrast, clozapine seems to promote hypomethylation across the epigenome. Candidate-gene methylation studies reveal varying effects post-antipsychotic therapy. Some genes like Glra1 and Drd2 are frequently found to undergo hypermethylation, whereas other genes such as SLC6A4, DUSP6, and DTNBP1 are more likely to exhibit hypomethylation in promoter regions. In examining histone modifications, the literature suggests that clozapine changes histone methylation patterns in the prefrontal cortex, particularly elevating H3K4me3 at the Gad1 gene and affecting the transcription of genes like mGlu2 by modifying histone acetylation and interacting with HDAC2 enzymes. Risperidone and quetiapine, however, exhibit distinct impacts on histone marks across different brain regions and cell types, with risperidone reducing H3K27ac in the striatum and quetiapine modifying global H3K9me2 levels in the prefrontal cortex, suggesting antipsychotics demonstrate selective influence on histone modifications, which demonstrates a complex and targeted mode of action. While this review summarizes current knowledge, the intricate dynamics between antipsychotics and epigenetics clearly warrant more exhaustive exploration with the potential to redefine our understanding and treatment of psychiatric conditions. By deciphering the epigenetic changes associated with drug treatment and therapeutic outcomes, we can move closer to personalized medicine in psychiatry.

Cognitive Dysfunction and Exercise: From Epigenetic to Genetic Molecular Mechanisms

Maintaining good health is crucial, and exercise plays a vital role in achieving this goal. It offers a range of positive benefits for cognitive function, regardless of age. However, as our population ages and life expectancy increases, cognitive impairment has become a prevalent issue, often coexisting with age-related neurodegenerative conditions. This can result in devastating consequences such as memory loss, difficulty speaking, and confusion, greatly hindering one's ability to lead an ordinary life. In addition, the decrease in mental capacity has a significant effect on an individual's physical and emotional well-being, greatly reducing their overall level of contentment and causing a significant financial burden for communities. While most current approaches aim to slow the decline of cognition, exercise offers a non-pharmacological, safe, and accessible solution. Its effects on cognition are intricate and involve changes in the brain's neural plasticity, mitochondrial stability, and energy metabolism. Moreover, exercise triggers the release of cytokines, playing a significant role in the body-brain connection and its impact on cognition. Additionally, exercise can influence gene expression through epigenetic mechanisms, leading to lasting improvements in brain function and behavior. Herein, we summarized various genetic and epigenetic mechanisms that can be modulated by exercise in cognitive dysfunction.

NMR metabolomics-guided DNA methylation mortality predictors

BACKGROUND

1H-NMR metabolomics and DNA methylation in blood are widely known biomarkers predicting age-related physiological decline and mortality yet exert mutually independent mortality and frailty signals.

METHODS

Leveraging multi-omics data in four Dutch population studies (N = 5238, ∼40% of which male) we investigated whether the mortality signal captured by 1H-NMR metabolomics could guide the construction of DNA methylation-based mortality predictors.

FINDINGS

We trained DNA methylation-based surrogates for 64 metabolomic analytes and found that analytes marking inflammation, fluid balance, or HDL/VLDL metabolism could be accurately reconstructed using DNA-methylation assays. Interestingly, a previously reported multi-analyte score indicating mortality risk (MetaboHealth) could also be accurately reconstructed. Sixteen of our derived surrogates, including the MetaboHealth surrogate, showed significant associations with mortality, independent of relevant covariates.

INTERPRETATION

The addition of our metabolic analyte-derived surrogates to the well-established epigenetic clock GrimAge demonstrates that our surrogates potentially represent valuable mortality signal.

FUNDING

BBMRI-NL, X-omics, VOILA, Medical Delta, NWO, ERC.

Map of epigenetic age acceleration: A worldwide analysis

We present a systematic analysis of epigenetic age acceleration based on by far the largest collection of publicly available DNA methylation data for healthy samples (93 datasets, 23 K samples), focusing on the geographic (25 countries) and ethnic (31 ethnicities) aspects around the world. We employed the most popular epigenetic tools for assessing age acceleration and examined their quality metrics and ability to extrapolate to epigenetic data from different tissue types and age ranges different from the training data of these models. In most cases, the models proved to be inconsistent with each other and showed different signs of age acceleration, with the PhenoAge model tending to systematically underestimate and different versions of the GrimAge model tending to systematically overestimate the age prediction of healthy subjects. Referring to data availability and consistency, most countries and populations are still not represented in GEO, moreover, different datasets use different criteria for determining healthy controls. Because of this, it is difficult to fully isolate the contribution of "geography/environment", "ethnicity" and "healthiness" to epigenetic age acceleration. Among the explored metrics, only the DunedinPACE, which measures aging rate, appears to adequately reflect the standard of living and socioeconomic indicators in countries, although it has a limited application to blood methylation data only. Invariably, by epigenetic age acceleration, males age faster than females in most of the studied countries and populations.

Non-nucleoside inhibitors of DNMT1 and DNMT3 for targeted cancer therapy

DNA methylation can deactivate tumor suppressor genes thus causing cancers. Two DNA methylation inhibitors have been approved by the Food and Drug Administration (FDA) and have entered clinical use. However, these inhibitors are nucleoside analogues that can be incorporated into DNA or RNA and induce significant side effects. DNMT1 and DNMT3 are key enzymes involved in DNA methylation. In the acute myeloid leukemia model, a non-nucleoside DNMT1-specific inhibitor has shown lower toxicity and improved pharmacokinetics compared to traditional nucleoside drugs. DNMT3 is also implicated in certain specific cancers. Thus, developing non-nucleoside inhibitors for DNMT1 or DNMT3 can help in understanding their roles in carcinogenesis and provide targeted treatment options in certain cancers. Although no non-nucleoside inhibitors have yet entered clinical trials, in this review, we focus on DNMT1 or DNMT3 selective inhibitors. For DNMT1 selective inhibitors, we have compiled information on the repurposed drugs, derivative compounds and selective inhibitors identified through virtual screening. Additionally, we have outlined potential targets for DNMT1, including protein-protein complex, RNA mimics and aptamers. Compared to DNMT1, research on DNMT3-specific inhibitors has been less extensive. In this context, our exploration has identified a limited number of molecular inhibitors, and we have proposed specific long non-coding RNAs (lncRNAs) as potential contributors to the selective inhibition of DNMT3. This collective effort aims to offer valuable insights into the development of non-nucleoside inhibitors that selectively target DNMT1 or DNMT3.

Biomarkers of Cellular Senescence and Aging: Current State-of-the-Art, Challenges and Future Perspectives

Population aging has increased the global prevalence of aging-related diseases, including cancer, sarcopenia, neurological disease, arthritis, and heart disease. Understanding aging, a fundamental biological process, has led to breakthroughs in several fields. Cellular senescence, evinced by flattened cell bodies, vacuole formation, and cytoplasmic granules, ubiquitously plays crucial roles in tissue remodeling, embryogenesis, and wound repair as well as in cancer therapy and aging. The lack of universal biomarkers for detecting and quantifying senescent cells, in vitro and in vivo, constitutes a major limitation. The applications and limitations of major senescence biomarkers, including senescence-associated β-galactosidase staining, telomere shortening, cell-cycle arrest, DNA methylation, and senescence-associated secreted phenotypes are discussed. Furthermore, explore senotherapeutic approaches for aging-associated diseases and cancer. In addition to the conventional biomarkers, this review highlighted the in vitro, in vivo, and disease models used for aging studies. Further, technologies from the current decade including multi-omics and computational methods used in the fields of senescence and aging are also discussed in this review. Understanding aging-associated biological processes by using cellular senescence biomarkers can enable therapeutic innovation and interventions to improve the quality of life of older adults.

Neuron and astrocyte specific 5mC and 5hmC signatures of BDNF's receptor, TrkB

Brain derived neurotrophic factor (BDNF) is the most studied trophic factor in the central nervous system (CNS), and its role in the maturation of neurons, including synapse development and maintenance has been investigated intensely for over three decades. The primary receptor for BDNF is the tropomyosin receptor kinase B (TrkB), which is broadly expressed as two primary isoforms in the brain; the full length TrkB (TrkB.FL) receptor, expressed mainly in neurons and the truncated TrkB (TrkB.T1) receptor. We recently demonstrated that TrkB.T1 is predominately expressed in astrocytes, and appears critical for astrocyte morphological maturation. Given the critical role of BDNF/TrkB pathway in healthy brain development and mature CNS function, we aimed to identify molecular underpinnings of cell-type specific expression of each TrkB isoform. Using Nanopore sequencing which enables direct, long read sequencing of native DNA, we profiled DNA methylation patterns of the entire TrkB gene, Ntrk2, in both neurons and astrocytes. Here, we identified robust differences in cell-type specific isoform expression associated with significantly different methylation patterns of the Ntrk2 gene in each cell type. Notably, astrocytes demonstrated lower 5mC methylation, and higher 5hmC across the entire gene when compared to neurons, including differentially methylated sites (DMSs) found in regions flanking the unique TrkB.T1 protein coding sequence (CDS). These data suggest DNA methylation patterns may provide instruction for isoform specific TrkB expression across unique CNS cell types.

Temporally discordant chromatin accessibility and DNA demethylation define short and long-term enhancer regulation during cell fate specification

Epigenetic mechanisms govern the transcriptional activity of lineage-specifying enhancers; but recent work challenges the dogma that joint chromatin accessibility and DNA demethylation are prerequisites for transcription. To understand this paradox, we established a highly-resolved timeline of DNA demethylation, chromatin accessibility, and transcription factor occupancy during neural progenitor cell differentiation. We show thousands of enhancers undergo rapid, transient accessibility changes associated with distinct periods of transcription factor expression. However, most DNA methylation changes are unidirectional and delayed relative to chromatin dynamics, creating transiently discordant epigenetic states. Genome-wide detection of 5-hydroxymethylcytosine further revealed active demethylation begins ahead of chromatin and transcription factor activity, while enhancer hypomethylation persists long after these activities have dissipated. We demonstrate that these timepoint specific methylation states predict past, present and future chromatin accessibility using machine learning models. Thus, chromatin and DNA methylation collaborate on different timescales to mediate short and long-term enhancer regulation during cell fate specification.

Emerging advances in defining the molecular and therapeutic landscape of small-cell lung cancer

Small-cell lung cancer (SCLC) has traditionally been considered a recalcitrant cancer with a dismal prognosis, with only modest advances in therapeutic strategies over the past several decades. Comprehensive genomic assessments of SCLC have revealed that most of these tumours harbour deletions of the tumour-suppressor genes TP53 and RB1 but, in contrast to non-small-cell lung cancer, have failed to identify targetable alterations. The expression status of four transcription factors with key roles in SCLC pathogenesis defines distinct molecular subtypes of the disease, potentially enabling specific therapeutic approaches. Overexpression and amplification of MYC paralogues also affect the biology and therapeutic vulnerabilities of SCLC. Several other attractive targets have emerged in the past few years, including inhibitors of DNA-damage-response pathways, epigenetic modifiers, antibody-drug conjugates and chimeric antigen receptor T cells. However, the rapid development of therapeutic resistance and lack of biomarkers for effective selection of patients with SCLC are ongoing challenges. Emerging single-cell RNA sequencing data are providing insights into the plasticity and intratumoural and intertumoural heterogeneity of SCLC that might be associated with therapeutic resistance. In this Review, we provide a comprehensive overview of the latest advances in genomic and transcriptomic characterization of SCLC with a particular focus on opportunities for translation into new therapeutic approaches to improve patient outcomes.

The correlation between CpG methylation and gene expression is driven by sequence variants

Gene promoter and enhancer sequences are bound by transcription factors and are depleted of methylated CpG sites (cytosines preceding guanines in DNA). The absence of methylated CpGs in these sequences typically correlates with increased gene expression, indicating a regulatory role for methylation. We used nanopore sequencing to determine haplotype-specific methylation rates of 15.3 million CpG units in 7,179 whole-blood genomes. We identified 189,178 methylation depleted sequences where three or more proximal CpGs were unmethylated on at least one haplotype. A total of 77,789 methylation depleted sequences (~41%) associated with 80,503 cis-acting sequence variants, which we termed allele-specific methylation quantitative trait loci (ASM-QTLs). RNA sequencing of 896 samples from the same blood draws used to perform nanopore sequencing showed that the ASM-QTL, that is, DNA sequence variability, drives most of the correlation found between gene expression and CpG methylation. ASM-QTLs were enriched 40.2-fold (95% confidence interval 32.2, 49.9) among sequence variants associating with hematological traits, demonstrating that ASM-QTLs are important functional units in the noncoding genome.

A varying-coefficient model for the analysis of methylation sequencing data

DNA methylation is an important epigenetic modification involved in gene regulation. Advances in the next generation sequencing technology have enabled the retrieval of DNA methylation information at single-base-resolution. However, due to the sequencing process and the limited amount of isolated DNA, DNA-methylation-data are often noisy and sparse, which complicates the identification of differentially methylated regions (DMRs), especially when few replicates are available. We present a varying-coefficient model for detecting DMRs by using single-base-resolved methylation information. The model simultaneously smooths the methylation profiles and allows detection of DMRs, while accounting for additional covariates. The proposed model takes into account possible overdispersion by using a beta-binomial distribution. The overdispersion itself can be modeled as a function of the genomic region and explanatory variables. We illustrate the properties of the proposed model by applying it to two real-life case studies.

Long-term storage does not affect the DNA methylation profiles of vitrified-warmed human embryos

With the widespread application of embryo cryopreservation in assisted reproductive techniques, it is necessary to assess the safety of long-term cryopreservation of human embryos and it is unclear whether storage time has an impact on the DNA methylation profiles of human embryos. Nine women who received IVF treatment were recruited for this study. The retrieved eight-cell human embryos were classified into three groups including fresh embryos, cryopreserved embryos stored for 3 years, and cryopreserved embryos stored for 8 years. Single-cell whole-genome bisulfite sequencing (scWGBS) was conducted. The genome-wide methylation pattern of the fresh and two cryopreserved groups were similar. In addition, the methylation level in different genomic regions showed comparable patterns and no significant differences were observed in the methylation level of imprinted genes among the three groups. A total of 587 differentially methylated regions (DMRs) in the 3-year group and 540 DMRs in the 8-year group were identified comparing to fresh group. However, they were not enriched in promoters and had a similar genome-wide distributions, suggesting that these DMRs may not contribute to the changes in corresponding gene expressions. Our study illustrated that long-term cryopreservation will not affect the DNA methylation profiles of human eight-cell embryos at single-cell level.

Beyond the base pairs: comparative genome-wide DNA methylation profiling across sequencing technologies

Deoxyribonucleic acid (DNA) methylation plays a key role in gene regulation and is critical for development and human disease. Techniques such as whole-genome bisulfite sequencing (WGBS) and reduced representation bisulfite sequencing (RRBS) allow DNA methylation analysis at the genome scale, with Illumina NovaSeq 6000 and MGI Tech DNBSEQ-T7 being popular due to their efficiency and affordability. However, detailed comparative studies of their performance are not available. In this study, we constructed 60 WGBS and RRBS libraries for two platforms using different types of clinical samples and generated approximately 2.8 terabases of sequencing data. We systematically compared quality control metrics, genomic coverage, CpG methylation levels, intra- and interplatform correlations, and performance in detecting differentially methylated positions. Our results revealed that the DNBSEQ platform exhibited better raw read quality, although base quality recalibration indicated potential overestimation of base quality. The DNBSEQ platform also showed lower sequencing depth and less coverage uniformity in GC-rich regions than did the NovaSeq platform and tended to enrich methylated regions. Overall, both platforms demonstrated robust intra- and interplatform reproducibility for RRBS and WGBS, with NovaSeq performing better for WGBS, highlighting the importance of considering these factors when selecting a platform for bisulfite sequencing.

Epigenetic modifications in Parkinson's disease: A critical review

Parkinson's Disease (PD) is a progressive neurodegenerative disorder expected to increase by over 50% by 2030 due to increasing life expectancy. The disease's hallmarks include slow movement, tremors, and postural instability. Impaired protein processing is a major factor in the pathophysiology of PD, leading to the buildup of aberrant protein aggregates, particularly misfolded α-synuclein, also known as Lewy bodies. These Lewy bodies lead to inflammation and further death of dopaminergic neurons, leading to imbalances in excitatory and inhibitory neurotransmitters, causing excessive uncontrollable movements called dyskinesias. It was previously suggested that a complex interplay involving hereditary and environmental variables causes the specific death of neurons in PD; however, the exact mechanism of the association involving the two primary modifiers is yet unknown. An increasing amount of research points to the involvement of epigenetics in the onset and course of several neurological conditions, such as PD. DNA methylation, post-modifications of histones, and non-coding RNAs are the primary examples of epigenetic alterations, that is defined as alterations to the expression of genes and functioning without modifications in DNA sequence. Epigenetic modifications play a significant role in the development of PD, with genes such as Parkin, PTEN-induced kinase 1 (PINK1), DJ1, Leucine-Rich Repeat Kinase 2 (LRRK2), and alpha-synuclein associated with the disease. The aberrant epigenetic changes implicated in the pathophysiology of PD and their impact on the design of novel therapeutic approaches are the primary focus of this review.

Decoding the Epigenetic Landscape: Insights into 5mC and 5hmC Patterns in Mouse Cortical Cell Types

The DNA modifications, 5-methylcytosine (5mC) and 5-hydroxymethylcytosine (5hmC), represent powerful epigenetic regulators of temporal and spatial gene expression. Yet, how the cooperation of these genome-wide, epigenetic marks determine unique transcriptional signatures across different brain cell populations is unclear. Here we applied Nanopore sequencing of native DNA to obtain a complete, genome-wide, single-base resolution atlas of 5mC and 5hmC modifications in neurons, astrocytes and microglia in the mouse cortex (99% genome coverage, 40 million CpG sites). In tandem with RNA sequencing, analysis of 5mC and 5hmC patterns across cell types reveals astrocytes drive uniquely high brain 5hmC levels and support two decades of research regarding methylation patterns, gene expression and alternative splicing, benchmarking this resource. As such, we provide the most comprehensive DNA methylation data in mouse brain as an interactive, online tool (NAM-Me, https://olsenlab.shinyapps.io/NAMME/) to serve as a resource dataset for those interested in the methylome landscape.

Transposable elements-mediated recruitment of KDM1A epigenetically silences HNF4A expression to promote hepatocellular carcinoma

Transposable elements (TEs) contribute to gene expression regulation by acting as cis-regulatory elements that attract transcription factors and epigenetic regulators. This research aims to explore the functional and clinical implications of transposable element-related molecular events in hepatocellular carcinoma, focusing on the mechanism through which liver-specific accessible TEs (liver-TEs) regulate adjacent gene expression. Our findings reveal that the expression of HNF4A is inversely regulated by proximate liver-TEs, which facilitates liver cancer cell proliferation. Mechanistically, liver-TEs are predominantly occupied by the histone demethylase, KDM1A. KDM1A negatively influences the methylation of histone H3 Lys4 (H3K4) of liver-TEs, resulting in the epigenetic silencing of HNF4A expression. The suppression of HNF4A mediated by KDM1A promotes liver cancer cell proliferation. In conclusion, this study uncovers a liver-TE/KDM1A/HNF4A regulatory axis that promotes liver cancer growth and highlights KDM1A as a promising therapeutic target. Our findings provide insight into the transposable element-related molecular mechanisms underlying liver cancer progression.

CpG Island Definition and Methylation Mapping of the T2T-YAO Genome

Precisely defining and mapping all cytosine (C) positions and their clusters, known as CpG islands (CGIs), as well as their methylation status, are pivotal for genome-wide epigenetic studies, especially when population-centric reference genomes are ready for timely application. Here, we first align the two high-quality reference genomes, T2T-YAO and T2T-CHM13, from different ethnic backgrounds in a base-by-base fashion and compute their genome-wide density-defined and position-defined CGIs. Second, by mapping some representative genome-wide methylation data from selected organs onto the two genomes, we find that there are about 4.7%-5.8% sequence divergency of variable categories depending on quality cutoffs. Genes among the divergent sequences are mostly associated with neurological functions. Moreover, CGIs associated with the divergent sequences are significantly different with respect to CpG density and observed CpG/expected CpG (O/E) ratio between the two genomes. Finally, we find that the T2T-YAO genome not only has a greater CpG coverage than that of the T2T-CHM13 genome when whole-genome bisulfite sequencing (WGBS) data from the European and American populations are mapped to each reference, but also shows more hyper-methylated CpG sites as compared to the T2T-CHM13 genome. Our study suggests that future genome-wide epigenetic studies of the Chinese populations rely on both acquisition of high-quality methylation data and subsequent precision CGI mapping based on the Chinese T2T reference.

A dynamic microRNA profile that tracks a chemotherapy resistance phenotype in osteosarcoma. Implications for novel therapeutics

Osteosarcoma is a rare primary bone tumor for which no significant therapeutic advancement has been made since the late 1980s despite ongoing efforts. Overall, the five-year survival rate remains about 65%, and is much lower in patients with tumors unresponsive to methotrexate, doxorubicin, and cisplatin therapy. Genetic studies have not revealed actionable drug targets, but our group, and others, have reported that epigenomic biomarkers, including regulatory RNAs, may be useful prognostic tools for osteosarcoma. We tested if microRNA (miRNA) transcriptional patterns mark the transition from a chemotherapy sensitive to resistant tumor phenotype. Small RNA sequencing was performed using 14 patient matched pre-chemotherapy biopsy and post-chemotherapy resection high-grade osteosarcoma frozen tumor samples. Independently, small RNA sequencing was performed using 14 patient matched biopsy and resection samples from untreated tumors. Separately, miRNA specific Illumina DASL arrays were used to assay an independent cohort of 65 pre-chemotherapy biopsy and 26 patient matched post-chemotherapy resection formalin fixed paraffin embedded (FFPE) tumor samples. mRNA specific Illumina DASL arrays were used to profile 37 pre-chemotherapy biopsy and five post-chemotherapy resection FFPE samples, all of which were also used for Illumina DASL miRNA profiling. The National Cancer Institute Therapeutically Applicable Research to Generate Effective Treatments dataset, including PCR based miRNA profiling and RNA-seq data for 86 and 93 pre-chemotherapy tumor samples, respectively, was also used. Paired differential expression testing revealed a profile of 17 miRNAs with significantly different transcriptional levels following chemotherapy. Genes targeted by the miRNAs were differentially expressed following chemotherapy, suggesting the miRNAs may regulate transcriptional networks. Finally, an in vitro pharmacogenomic screen using miRNAs and their target transcripts predicted response to a set of candidate small molecule therapeutics which potentially reverse the chemotherapy resistance phenotype and synergize with chemotherapy in otherwise treatment resistant tumors. Importantly, these novel therapeutic targets are distinct from targets identified by a similar pharmacogenomic analysis of previously published prognostic miRNA profiles from pre chemotherapy biopsy specimens.

The roles of epigenetic regulation in graft-versus-host disease

Allogeneic hematopoietic stem cell transplantation (aHSCT) is utilized as a potential curative treatment for various hematologic malignancies. However, graft-versus-host disease (GVHD) post-aHSCT is a severe complication that significantly impacts patients' quality of life and overall survival, becoming a major cause of non-relapse mortality. In recent years, the association between epigenetics and GVHD has garnered increasing attention. Epigenetics focuses on studying mechanisms that affect gene expression without altering DNA sequences, primarily including DNA methylation, histone modifications, non-coding RNAs (ncRNAs) regulation, and RNA modifications. This review summarizes the role of epigenetic regulation in the pathogenesis of GVHD, with a focus on DNA methylation, histone modifications, ncRNA, RNA modifications and their involvement and applications in the occurrence and development of GVHD. It also highlights advancements in relevant diagnostic markers and drugs, aiming to provide new insights for the clinical diagnosis and treatment of GVHD.

Disruption of maternal vascular remodeling by a fetal endoretrovirus-derived gene in preeclampsia

BACKGROUND

Preeclampsia, one of the most lethal pregnancy-related diseases, is associated with the disruption of uterine spiral artery remodeling during placentation. However, the early molecular events leading to preeclampsia remain unknown.

RESULTS

By analyzing placentas from preeclampsia, non-preeclampsia, and twin pregnancies with selective intrauterine growth restriction, we show that the pathogenesis of preeclampsia is attributed to immature trophoblast and maldeveloped endothelial cells. Delayed epigenetic reprogramming during early extraembryonic tissue development leads to generation of excessive immature trophoblast cells. We find reduction of de novo DNA methylation in these trophoblast cells results in selective overexpression of maternally imprinted genes, including the endoretrovirus-derived gene PEG10 (paternally expressed gene 10). PEG10 forms virus-like particles, which are transferred from the trophoblast to the closely proximate endothelial cells. In normal pregnancy, only a low amount of PEG10 is transferred to maternal cells; however, in preeclampsia, excessive PEG10 disrupts maternal vascular development by inhibiting TGF-beta signaling.

CONCLUSIONS

Our study reveals the intricate epigenetic mechanisms that regulate trans-generational genetic conflict and ultimately ensure proper maternal-fetal interface formation.

DNA methylation patterns of FKBP5 regulatory regions in brain and blood of humanized mice and humans

Humanized mouse models can be used to explore human gene regulatory elements (REs), which frequently lie in non-coding and less conserved genomic regions. Epigenetic modifications of gene REs, also in the context of gene x environment interactions, have not yet been explored in humanized mouse models. We applied high-accuracy measurement of DNA methylation (DNAm) via targeted bisulfite sequencing (HAM-TBS) to investigate DNAm in three tissues/brain regions (blood, prefrontal cortex and hippocampus) of mice carrying the human FK506-binding protein 5 (FKBP5) gene, an important candidate gene associated with stress-related psychiatric disorders. We explored DNAm in three functional intronic glucocorticoid-responsive elements (at introns 2, 5, and 7) of FKBP5 at baseline, in cases of differing genotype (rs1360780 single nucleotide polymorphism), and following application of the synthetic glucocorticoid dexamethasone. We compared DNAm patterns in the humanized mouse (N = 58) to those in human peripheral blood (N = 447 and N = 89) and human postmortem brain prefrontal cortex (N = 86). Overall, DNAm patterns in the humanized mouse model seem to recapitulate DNAm patterns observed in human tissue. At baseline, this was to a higher extent in brain tissue. The animal model also recapitulated effects of dexamethasone on DNAm, especially in peripheral blood and to a lesser extent effects of genotype on DNAm. The humanized mouse model could thus assist in reverse translation of human findings in psychiatry that involve genetic and epigenetic regulation in non-coding elements.

Enhanced resolution profiling in twins reveals differential methylation signatures of type 2 diabetes with links to its complications

BACKGROUND

Type 2 diabetes (T2D) susceptibility is influenced by genetic and environmental factors. Previous findings suggest DNA methylation as a potential mechanism in T2D pathogenesis and progression.

METHODS

We profiled DNA methylation in 248 blood samples from participants of European ancestry from 7 twin cohorts using a methylation sequencing platform targeting regulatory genomic regions encompassing 2,048,698 CpG sites.

FINDINGS

We find and replicate 3 previously unreported T2D differentially methylated CpG positions (T2D-DMPs) at FDR 5% in RGL3, NGB and OTX2, and 20 signals at FDR 25%, of which 14 replicated. Integrating genetic variation and T2D-discordant monozygotic twin analyses, we identify both genetic-based and genetic-independent T2D-DMPs. The signals annotate to genes with established GWAS and EWAS links to T2D and its complications, including blood pressure (RGL3) and eye disease (OTX2).

INTERPRETATION

The results help to improve our understanding of T2D disease pathogenesis and progression and may provide biomarkers for its complications.

FUNDING

Funding acknowledgements for each cohort can be found in the Supplementary Note.

Understanding the role of DNA methylation in colorectal cancer: Mechanisms, detection, and clinical significance

Colorectal cancer (CRC) is one of the deadliest malignancies worldwide, ranking third in incidence and second in mortality. Remarkably, early stage localized CRC has a 5-year survival rate of over 90%; in stark contrast, the corresponding 5-year survival rate for metastatic CRC (mCRC) is only 14%. Compounding this problem is the staggering lack of effective therapeutic strategies. Beyond genetic mutations, which have been identified as critical instigators of CRC initiation and progression, the importance of epigenetic modifications, particularly DNA methylation (DNAm), cannot be underestimated, given that DNAm can be used for diagnosis, treatment monitoring and prognostic evaluation. This review addresses the intricate mechanisms governing aberrant DNAm in CRC and its profound impact on critical oncogenic pathways. In addition, a comprehensive review of the various techniques used to detect DNAm alterations in CRC is provided, along with an exploration of the clinical utility of cancer-specific DNAm alterations.

Sex-Specific Deflection of Age-Related DNA Methylation and Gene Expression in Mouse Heart by Perinatal Toxicant Exposures

Background

Global and site-specific changes in DNA methylation and gene expression are associated with cardiovascular aging and disease, but how toxicant exposures during early development influence the normal trajectory of these age-related molecular changes, and whether there are sex differences, has not yet been investigated.

Objectives

We used an established mouse model of developmental exposures to investigate the effects of perinatal exposure to either lead (Pb) or diethylhexyl phthalate (DEHP), two ubiquitous environmental contaminants strongly associated with CVD, on age-related cardiac DNA methylation and gene expression.

Methods

Dams were randomly assigned to receive human physiologically relevant levels of Pb (32 ppm in water), DEHP (25 mg/kg chow), or control water and chow. Exposures started two weeks prior to mating and continued until weaning at postnatal day 21 (3 weeks of age). Approximately one male and one female offspring per litter were followed to 3 weeks, 5 months, or 10 months of age, at which time whole hearts were collected (n ≥ 5 per sex per exposure). Enhanced reduced representation bisulfite sequencing (ERRBS) was used to assess the cardiac DNA methylome at 3 weeks and 10 months, and RNA-seq was conducted at all 3 time points. MethylSig and edgeR were used to identify age-related differentially methylated regions (DMRs) and differentially expressed genes (DEGs), respectively, within each sex and exposure group. Cell type deconvolution of bulk RNA-seq data was conducted using the MuSiC algorithm and publicly available single cell RNA-seq data.

Results

Thousands of DMRs and hundreds of DEGs were identified in control, DEHP, and Pb-exposed hearts across time between 3 weeks and 10 months of age. A closer look at the genes and pathways showing differential DNA methylation revealed that the majority were unique to each sex and exposure group. Overall, pathways governing development and differentiation were most frequently altered with age in all conditions. A small number of genes in each group showed significant changes in DNA methylation and gene expression with age, including several that were altered by both toxicants but were unchanged in control. We also observed subtle, but significant changes in the proportion of several cell types due to age, sex, and developmental exposure.

Discussion

Together these data show that perinatal Pb or DEHP exposures deflect normal age-related gene expression, DNA methylation programs, and cellular composition across the life course, long after cessation of exposure, and highlight potential biomarkers of developmental toxicant exposures. Further studies are needed to investigate how these epigenetic and transcriptional changes impact cardiovascular health across the life course.

Pregnancy is linked to faster epigenetic aging in young women

A central prediction of evolutionary theory is that energy invested into reproduction comes at the expense of somatic maintenance and repair, accelerating biological aging. Supporting this prediction are findings that high fertility among women predicts shorter lifespan and poorer health later in life. However, biological aging is thought to begin before age-related health declines, limiting the applicability of morbidity and mortality for studying the aging process earlier in life. Here, we examine the relationship between reproductive history and biological aging in a sample of young (20 to 22yo) men and women from the Cebu Longitudinal Health and Nutrition Survey, located in the Philippines (n = 1,735). We quantify biological aging using six measures, collectively known as epigenetic clocks, reflecting various facets of cellular aging, health, and mortality risk. In a subset of women, we test whether longitudinal changes in gravidity between young and early-middle adulthood (25 to 31yo) are associated with changes in epigenetic aging during that time. Cross-sectionally, gravidity was associated with all six measures of accelerated epigenetic aging in women (n = 825). Furthermore, longitudinal increases in gravidity were linked to accelerated epigenetic aging in two epigenetic clocks (n = 331). In contrast, the number of pregnancies a man reported fathering was not associated with epigenetic aging among same-aged cohort men (n = 910). These effects were robust to socioecological, environmental, and immunological factors, consistent with the hypothesis that pregnancy accelerates biological aging and that these effects can be detected in young women in a high-fertility context.

Epigenomic insights into common human disease pathology

The epigenome-the chemical modifications and chromatin-related packaging of the genome-enables the same genetic template to be activated or repressed in different cellular settings. This multi-layered mechanism facilitates cell-type specific function by setting the local sequence and 3D interactive activity level. Gene transcription is further modulated through the interplay with transcription factors and co-regulators. The human body requires this epigenomic apparatus to be precisely installed throughout development and then adequately maintained during the lifespan. The causal role of the epigenome in human pathology, beyond imprinting disorders and specific tumour suppressor genes, was further brought into the spotlight by large-scale sequencing projects identifying that mutations in epigenomic machinery genes could be critical drivers in both cancer and developmental disorders. Abrogation of this cellular mechanism is providing new molecular insights into pathogenesis. However, deciphering the full breadth and implications of these epigenomic changes remains challenging. Knowledge is accruing regarding disease mechanisms and clinical biomarkers, through pathogenically relevant and surrogate tissue analyses, respectively. Advances include consortia generated cell-type specific reference epigenomes, high-throughput DNA methylome association studies, as well as insights into ageing-related diseases from biological 'clocks' constructed by machine learning algorithms. Also, 3rd-generation sequencing is beginning to disentangle the complexity of genetic and DNA modification haplotypes. Cell-free DNA methylation as a cancer biomarker has clear clinical utility and further potential to assess organ damage across many disorders. Finally, molecular understanding of disease aetiology brings with it the opportunity for exact therapeutic alteration of the epigenome through CRISPR-activation or inhibition.

Tissue informative cell-free DNA methylation sites in amyotrophic lateral sclerosis

Cell-free DNA (cfDNA) is increasingly recognized as a promising biomarker candidate for disease monitoring. However, its utility in neurodegenerative diseases, like amyotrophic lateral sclerosis (ALS), remains underexplored. Existing biomarker discovery approaches are tailored to a specific disease context or are too expensive to be clinically practical. Here, we address these challenges through a new approach combining advances in molecular and computational technologies. First, we develop statistical tools to select tissue-informative DNA methylation sites relevant to a disease process of interest. We then employ a capture protocol to select these sites and perform targeted methylation sequencing. Multi-modal information about the DNA methylation patterns are then utilized in machine learning algorithms trained to predict disease status and disease progression. We applied our method to two independent cohorts of ALS patients and controls (n=192). Overall, we found that the targeted sites accurately predicted ALS status and replicated between cohorts. Additionally, we identified epigenetic features associated with ALS phenotypes, including disease severity. These findings highlight the potential of cfDNA as a non-invasive biomarker for ALS.

Genomic origin, fragmentomics, and transcriptional properties of long cell-free DNA molecules in human plasma

Recent studies have revealed an unexplored population of long cell-free DNA (cfDNA) molecules in human plasma using long-read sequencing technologies. However, the biological properties of long cfDNA molecules (>500 bp) remain largely unknown. To this end, we have investigated the origins of long cfDNA molecules from different genomic elements. Analysis of plasma cfDNA using long-read sequencing reveals an uneven distribution of long molecules from across the genome. Long cfDNA molecules show overrepresentation in euchromatic regions of the genome, in sharp contrast to short DNA molecules. We observe a stronger relationship between the abundance of long molecules and mRNA gene expression levels, compared with short molecules (Pearson's r = 0.71 vs. -0.14). Moreover, long and short molecules show distinct fragmentation patterns surrounding CpG sites. Leveraging the cleavage preferences surrounding CpG sites, the combined cleavage ratios of long and short molecules can differentiate patients with hepatocellular carcinoma (HCC) from non-HCC subjects (AUC = 0.87). We also investigated knockout mice in which selected nuclease genes had been inactivated in comparison with wild-type mice. The proportion of long molecules originating from transcription start sites are lower in Dffb-deficient mice but higher in Dnase1l3-deficient mice compared with that of wild-type mice. This work thus provides new insights into the biological properties and potential clinical applications of long cfDNA molecules.

DNA methylation haplotype block signatures responding to Staphylococcus aureus subclinical mastitis and association with production and health traits

BACKGROUND

DNA methylation has been documented to play vital roles in diseases and biological processes. In bovine, little is known about the regulatory roles of DNA methylation alterations on production and health traits, including mastitis.

RESULTS

Here, we employed whole-genome DNA methylation sequencing to profile the DNA methylation patterns of milk somatic cells from sixteen cows with naturally occurring Staphylococcus aureus (S. aureus) subclinical mastitis and ten healthy control cows. We observed abundant DNA methylation alterations, including 3,356,456 differentially methylated cytosines and 153,783 differential methylation haplotype blocks (dMHBs). The DNA methylation in regulatory regions, including promoters, first exons and first introns, showed global significant negative correlations with gene expression status. We identified 6435 dMHBs located in the regulatory regions of differentially expressed genes and significantly correlated with their corresponding genes, revealing their potential effects on transcriptional activities. Genes harboring DNA methylation alterations were significantly enriched in multiple immune- and disease-related pathways, suggesting the involvement of DNA methylation in regulating host responses to S. aureus subclinical mastitis. In addition, we found nine discriminant signatures (differentiates cows with S. aureus subclinical mastitis from healthy cows) representing the majority of the DNA methylation variations related to S. aureus subclinical mastitis. Validation of seven dMHBs in 200 cows indicated significant associations with mammary gland health (SCC and SCS) and milk production performance (milk yield).

CONCLUSIONS

In conclusion, our findings revealed abundant DNA methylation alterations in milk somatic cells that may be involved in regulating mammary gland defense against S. aureus infection. Particularly noteworthy is the identification of seven dMHBs showing significant associations with mammary gland health, underscoring their potential as promising epigenetic biomarkers. Overall, our findings on DNA methylation alterations offer novel insights into the regulatory mechanisms of bovine subclinical mastitis, providing further avenues for the development of effective control measures.

DNA methylation may partly explain psychotropic drug-induced metabolic side effects: results from a prospective 1-month observational study

BACKGROUND

Metabolic side effects of psychotropic medications are a major drawback to patients' successful treatment. Using an epigenome-wide approach, we aimed to investigate DNA methylation changes occurring secondary to psychotropic treatment and evaluate associations between 1-month metabolic changes and both baseline and 1-month changes in DNA methylation levels. Seventy-nine patients starting a weight gain inducing psychotropic treatment were selected from the PsyMetab study cohort. Epigenome-wide DNA methylation was measured at baseline and after 1 month of treatment, using the Illumina Methylation EPIC BeadChip.

RESULTS

A global methylation increase was noted after the first month of treatment, which was more pronounced (p < 2.2 × 10-16) in patients whose weight remained stable (< 2.5% weight increase). Epigenome-wide significant methylation changes (p < 9 × 10-8) were observed at 52 loci in the whole cohort. When restricting the analysis to patients who underwent important early weight gain (≥ 5% weight increase), one locus (cg12209987) showed a significant increase in methylation levels (p = 3.8 × 10-8), which was also associated with increased weight gain in the whole cohort (p = 0.004). Epigenome-wide association analyses failed to identify a significant link between metabolic changes and methylation data. Nevertheless, among the strongest associations, a potential causal effect of the baseline methylation level of cg11622362 on glycemia was revealed by a two-sample Mendelian randomization analysis (n = 3841 for instrument-exposure association; n = 314,916 for instrument-outcome association).

CONCLUSION

These findings provide new insights into the mechanisms of psychotropic drug-induced weight gain, revealing important epigenetic alterations upon treatment, some of which may play a mediatory role.

UHRF1 ubiquitin ligase activity supports the maintenance of low-density CpG methylation

The RING E3 ubiquitin ligase UHRF1 is an established cofactor for DNA methylation inheritance. Nucleosomal engagement through histone and DNA interactions directs UHRF1 ubiquitin ligase activity toward lysines on histone H3 tails, creating binding sites for DNMT1 through ubiquitin interacting motifs (UIM1 and UIM2). Here, we profile contributions of UHRF1 and DNMT1 to genome-wide DNA methylation inheritance and dissect specific roles for ubiquitin signaling in this process. We reveal DNA methylation maintenance at low-density CpGs is vulnerable to disruption of UHRF1 ubiquitin ligase activity and DNMT1 ubiquitin reading activity through UIM1. Hypomethylation of low-density CpGs in this manner induces formation of partially methylated domains (PMD), a methylation signature observed across human cancers. Furthermore, disrupting DNMT1 UIM2 function abolishes DNA methylation maintenance. Collectively, we show DNMT1-dependent DNA methylation inheritance is a ubiquitin-regulated process and suggest a disrupted UHRF1-DNMT1 ubiquitin signaling axis contributes to the development of PMDs in human cancers.

Increased CpG methylation at the CDH1 locus in inflamed ileal mucosa of patients with Crohn disease

BACKGROUND

E-cadherin, a major actor of cell adhesion in the intestinal barrier, is encoded by the CDH1 gene associated with susceptibility to Crohn Disease (CD) and colorectal cancer. Since epigenetic mechanisms are suspected to contribute to the multifactorial pathogenesis of CD, we studied CpG methylation at the CDH1 locus. The methylation of the CpG island (CGI) and of the 1st enhancer, two critical regulatory positions, was quantified in surgical specimens of inflamed ileal mucosa and in peripheral blood mononuclear cells (PBMC) of 21 CD patients. Sixteen patients operated on for a non-inflammatory bowel disease, although not normal controls, provided a macroscopically normal ileal mucosa and PBMC for comparison.

RESULTS

In ileal mucosa, 19/21 (90%) CD patients vs 8/16 control patients (50%) (p < 0.01) had a methylated CDH1 promoter CGI. In PBMC, CD patients with methylated CGI were 11/21 (52%) vs 7/16 controls (44%), respectively. Methylation in the 1st enhancer of CDH1 was also higher in the CD group for each of the studied CpGs and for their average value (45 ± 17% in CD patients vs 36 ± 17% in controls; p < 0.001). Again, methylation was comparable in PBMC. Methylation of CGI and 1st enhancer were not correlated in mucosa or PBMC.

CONCLUSIONS

Methylation of several CpGs at the CDH1 locus was increased in the inflamed ileal mucosa, not in the PBMC, of CD patients, suggesting the association of CDH1 methylation with ileal inflammation. Longitudinal studies will explore if this increased methylation is a risk marker for colorectal cancer.

Identification of differentially methylated regions associated with both liver fibrosis and hepatocellular carcinoma

BACKGROUND

Liver fibrosis is a major risk factor for hepatocellular carcinoma (HCC). We have previously reported that differentially methylated regions (DMRs) are correlated with the fibrosis stages of metabolic dysfunction-associated steatotic liver disease (MASLD). In this study, the methylation levels of those DMRs in liver fibrosis and subsequent HCC were examined.

METHODS

The methylation levels of DMRs were investigated using alcoholic cirrhosis and HCC (GSE60753). The data of hepatitis C virus-infected cirrhosis and HCC (GSE60753), and two datasets (GSE56588 and GSE89852) were used for replication analyses. The transcriptional analyses were performed using GSE114564, GSE94660, and GSE142530.

RESULTS

Hypomethylated DMR and increased transcriptional level of zinc finger and BTB domain containing 38 (ZBTB38) were observed in HCC. Hypermethylated DMRs, and increased transcriptional levels of forkhead box K1 (FOXK1) and zinc finger CCCH-type containing 3 (ZC3H3) were observed in HCC. The methylation levels of DMR of kazrin, periplakin interacting protein (KAZN) and its expression levels were gradually decreased as cirrhosis progressed to HCC.

CONCLUSIONS

Changes in the methylation and transcriptional levels of ZBTB38, ZC3H3, FOXK1, and KAZN are important for the development of fibrosis and HCC; and are therefore potential therapeutic targets and diagnostic tools for cirrhosis and HCC.

Differential methylation of linoleic acid pathway genes is associated with PTSD symptoms - a longitudinal study with Burundian soldiers returning from a war zone

Soldiers may be exposed to traumatic stress during combat deployment and thus are at risk for developing posttraumatic stress disorder (PTSD). Genetic and epigenetic evidence suggests that PTSD is linked to forming stress-related memories. In the current study, we investigated post-deployment associations of PTSD symptoms with differential DNA methylation in a sample of Burundian soldiers returning from the African Union Mission in Somalia's war zone. We used a matched longitudinal study design to explore epigenetic changes associated with PTSD symptoms in N = 191 participants. PTSD symptoms and saliva samples were collected at 1-3 (t1) and 9-14 months (t2) after the return of the soldiers to their home base. Individuals with either worsening or improving PTSD symptoms were matched for age, stressful, traumatic and self-perpetrated events prior to the post-assessment, traumatic and violent experiences between the post- and the follow-up assessment, and violence experienced during childhood. A mixed model analysis was conducted to identify top nominally significantly differentially methylated genes, which were then used to perform a gene enrichment analysis. The linoleic acid metabolism pathway was significantly associated with post-deployment PTSD symptoms, after accounting for multiple comparisons. Linoleic acid has been linked to memory and immune related processes in previous research. Our findings suggest that differential methylation of linoleic acid pathway genes is associated with PTSD and thus may merit closer inspection as a possible mediator of resilience.

Changes in DNA methylation associated with a specific mode of delivery: a pilot study

Background

The mode of delivery represents an epigenetic factor with potential to affect further development of the individual by multiple mechanisms. DNA methylation may be one of them, representing a major epigenetic mechanism involving direct chemical modification of the individual's DNA. This pilot study aims to examine whether a specific mode of delivery induces changes of DNA methylation by comparing the umbilical cord blood and peripheral blood of the newborns.

Methods

Blood samples from infants born by vaginal delivery and caesarean section were analysed to prepare the Methylseq library according to NEBNext enzymatic Methyl-seq Methylation Library Preparation Kit with further generation of target-enriched DNA libraries using the Twist Human Methylome Panel. DNA methylation status was determined using Illumina next-generation sequencing (NGS).

Results

We identified 168 differentially methylated regions in umbilical cord blood samples and 157 regions in peripheral blood samples. These were associated with 59 common biological, metabolic and signalling pathways for umbilical cord and peripheral blood samples.

Conclusion

Caesarean section is likely to represent an important epigenetic factor with the potential to induce changes in the genome that could play an important role in development of a broad spectrum of disorders. Our results could contribute to the elucidation of how epigenetic factors, such as a specific mode of delivery, could have adverse impact on health of an individual later in their life.