metilene: fast and sensitive calling of differentially methylated regions from bisulfite sequencing data

Mediating effects of BMI on the association between DNA methylation regions and 24-h blood pressure in African Americans

BACKGROUND

DNA methylation is an important epigenetic mechanism that may influence blood pressure (BP) regulation and hypertension risk. Obesity, a major lifestyle factor associated with hypertension, may interact with DNA methylation to affect BP. However, the indirect effect of DNA methylation on 24-h BP measurements mediated by obesity-related phenotypes such as BMI has not been investigated.

METHODS

Causal mediation analysis was applied to examine the mediating role of BMI in the relation between DNA methylation and 24-h BP phenotypes, including SBP, DBP and mean arterial blood pressure (MAP), in 281 African American participants.

RESULTS

Analysis of 38 215 DNA methylation regions, derived from 1 549 368 CpG sites across the genome, identified up to 138 methylation regions that were significantly associated with 24-h BP measurements through BMI mediation. Among them, 38 (19.2%) methylation regions were concurrently associated with SBP, DBP and MAP. Genes associated with BMI-mediated methylation regions are potentially involved in various chronic diseases such as coronary artery disease and renal disease, which are often caused or exacerbated by hypertension. Notably, three genes ( CDH4 , NOTCH1 and COLGALT1 ) showed both direct associations with 24-h BP measurements and indirect associations through BMI after adjusting for age and sex covariates.

CONCLUSION

Our findings suggest that DNA methylation may contribute to the regulation of 24-h BP in African Americans both directly and indirectly through BMI mediation.

Grandmotherhood is associated with reduced OXTR DNA methylation

In mammals, both parental and alloparental care are associated with increased brain oxytocin signaling. Grandmothers are important alloparents in many human families. Based on animal model research showing that peripheral Oxtr methylation is associated with Oxtr expression in the nucleus accumbens, we investigated whether grandmaternal caregiving is associated with lower peripheral OXTR methylation. Results reveal several regions within OXTR where grandmothers have lower DNA methylation compared with non-grandmother controls, and no regions where grandmothers have higher OXTR DNA methylation. Among grandmothers, OXTR methylation was most strongly correlated with the grandmother's assessment of the degree of positive feelings between her and the grandchild, which in turn predicted caregiving engagement. Although there was little evidence that grandmaternal OXTR methylation modulated grandmaternal neural responses to viewing photos of the grandchild within brain regions involved in caregiving motivation, it was negatively correlated with the neural response to an unknown grandchild. Thus, while OT signaling may not be essential for activating grandmaternal brain reward systems in our low-stress experimental context, it may support caregiving motivation towards unrelated children. Future longitudinal research should determine whether the transition to grandmotherhood is associated with a reduction in OXTR methylation.

Reduced representation bisulfite sequencing (RRBS) analysis reveals variation in distribution and levels of DNA methylation in white birch (Betula papyrifera) exposed to nickel

Research in understanding the role of genetics and epigenetics in plant adaptations to environmental stressors such as metals is still in its infancy. The objective of the present study is to assess the effect of nickel on DNA methylation level and distribution in white birch (Betula papyrifera Marshall) using reduced representation bisulfite sequencing (RRBS). The distribution of methylated C sites of each sample revealed that the level of methylation was much higher in CG context varying between 54% and 65%, followed by CHG (24%-31.5%), and then CHH with the methylation rate between 3.3% and 5.2%. The analysis of differentially methylated regions (DMR) revealed that nickel induced both hypermethylation and hypomethylation when compared to water. Detailed analysis showed for the first time that nickel induced a higher level of hypermethylation compared to controls, while potassium triggers a higher level of hypomethylation compared to nickel. Surprisingly, the analysis of the distribution of DMRs revealed that 38%-42% were located in gene bodies, 20%-24% in exon, 19%-20% in intron, 16%-17% in promoters, and 0.03%-0.04% in transcription start site. RRBS was successful in detecting and mapping DMR in plants exposed to nickel.

The chromatin accessibility landscape of mouse oocytes during configuration transition

The transition of chromatin configuration in mammalian oocytes from a non-surrounded nucleolus (NSN) to a surrounded nucleolus (SN) is critical for acquiring the developmental competence. However, the genomic and epigenomic features underlying this process remain poorly understood. In the present study, we first establish the chromatin accessibility landscape of mouse oocytes from NSN to SN stage. Through the integrative analysis of multi-omics, we find that the establishment of DNA methylation in oocytes is independent of the dynamics of chromatin accessibility. In contrast, histone H3K4me3 status is closely associated with the dynamics of accessible regions during configuration transition. Furthermore, by focusing on the actively transcribed genes in NSN and SN oocytes, we discover that chromatin accessibility coupled with histone methylation (H3K4me3 and H3K27me3) participates in the transcriptional control during phase transition. In sum, our data provide a comprehensive resource for probing configuration transition in oocytes, and offer insights into the mechanisms determining chromatin dynamics and oocyte quality.

Paternal high-fat diet affects weight and DNA methylation of their offspring

Obesity poses a public health threat, reaching epidemic proportions. Our hypothesis suggests that some of this epidemic stems from its transmission across generations via paternal epigenetic mechanisms. To investigate this possibility, we focused on examining the paternal transmission of CpG methylation. First-generation male Wistar rats were fed either a high-fat diet (HF) or chow and were mated with females fed chow. We collected sperm from these males. The resulting offspring were raised on a chow diet until day 35, after which they underwent a dietary challenge. Diet-induced obese (DIO) male rats passed on the obesogenic trait to both male and female offspring. We observed significant hypermethylation of the Pomc promoter in the sperm of HF-treated males and in the hypothalamic arcuate nucleus (Arc) of their offspring at weaning. However, these differences in Arc methylation decreased later in life. This hypermethylation is correlated with increased expression of DNMT3B. Further investigating genes in the Arc that might be involved in obesogenic transgenerational transmission, using reduced representation bisulfite sequencing (RRBS) we identified 77 differentially methylated regions (DMRs), highlighting pathways associated with neuronal development. These findings support paternal CpG methylation as a mechanism for transmitting obesogenic traits across generations.

DNA methylation and gene expression profiling reveal potential association of retinol metabolism related genes with hepatocellular carcinoma development

Background

Aberrant DNA methylation patterns play a critical role in the development of hepatocellular carcinoma (HCC). However, the molecular mechanisms associated with these aberrantly methylated genes remain unclear. This study aimed to comprehensively investigate the methylation-driven gene expression alterations in HCC using a multi-omics dataset.

Methods

Whole genome bisulfite sequencing (WGBS) and RNA sequencing (RNA-seq) techniques were used to assess the methylation and gene expression profiles of HCC tissues (HCCs) and normal adjacent tissues (NATs). The candidate genes' potential function was further investigated using single-cell RNA sequencing (scRNA seq) data.

Results

We observed widespread hypomethylation in HCCs compared to NATs. Methylation levels in distinct genomic regions exhibited significant differences between HCCs and NATs. We identified 247,632 differentially methylated regions (DMRs) and 4,926 differentially expressed genes (DEGs) between HCCs and NATs. Integrated analysis of DNA methylation and RNA-seq data identified 987 methylation-driven candidate genes, with 970 showing upregulation and 17 showing downregulation. Four genes involved in the retinol metabolic pathway, namely ADH1A, CYP2A6, CYP2C8, and CYP2C19, were identified as hyper-downregulated genes. Their expression levels could stratify HCCs into three subgroups with distinct survival outcomes, immune cell infiltration, and tumor microenvironments. Validation of these findings in an independent dataset yielded similar outcomes, confirming the high concordance and potential prognostic value of these genes. ScRNA seq data revealed the low expression of these genes in immune cells, emphasizing their role in promoting malignant cell proliferation and migration. In conclusion, this study provides insights into the molecular characteristics of HCC, revealing the involvement of retinol metabolism-related genes in the development and progression of HCC. These findings have implications for HCC diagnosis, prognosis prediction, and the development of therapeutic targets.

Comprehensive Analysis of Methylome and Transcriptome to Identify Potential Genes Regulating Porcine Testis Development

DNA methylation plays a critical role in regulating gene expression during testicular development. However, few studies report on candidate genes related to the DNA methylation regulation of porcine testicular development. This study examined the differentially expressed genes (DEGs) and their methylation levels in testicular tissues from pigs at 60 days of age (60 d) and 180 days of age (180 d) using RNA-Seq and whole genome bisulfite sequencing (WGBS). It was determined that DNA methylation primarily occurs in the cytosine-guanine (CG) context, and the analysis identified 106,282 differentially methylated regions (DMRs) corresponding to 12,385 differentially methylated genes (DMGs). Further integrated analysis of RNA-Seq and WGBS data revealed 1083 DMGs negatively correlated with the expression of DEGs. GO analysis showed that these genes were significantly enriched in spermatogenesis, germ cell development, and spermatid differentiation. The screening of enriched genes revealed that hyper-methylation repressed ADAM30, ADAM3A, DPY19L2, H2BC1, MAK, RPL10L, SPATA16, and YBX2, while hypo-methylation elevated CACNA1I, CADM1, CTNNB1, JAM2, and PAFAH1B3 expression. Additionally, the methylation status of the key genes ADAM3A, ADAM30, YBX2, JAM2, PAFAH1B3, and CTNNB1 was detected by bisulfite sequencing PCR (BSP). This study offers insights into the epigenetic regulation mechanisms underlying porcine testicular development.

Genome-Wide Methylation Profiling of Peripheral T-Cell Lymphomas Identifies TRIP13 as a Critical Driver of Tumor Proliferation and Survival

Cytosine methylation contributes to the regulation of gene expression and normal hematopoiesis in mammals. It is catalyzed by the family of DNA methyltransferases that include DNMT1, DNMT3A, and DNMT3B. Peripheral T-cell lymphomas (PTCLs) represent aggressive mature T-cell malignancies exhibiting a broad spectrum of clinical features with poor prognosis and inadequately understood molecular pathobiology. To better understand the molecular landscape and identify candidate genes involved in disease maintenance, we profiled DNA methylation and gene expression of PTCLs. We found that the methylation patterns in PTCLs are deregulated and heterogeneous but share 767 hypo- and 567 hypermethylated differentially methylated regions (DMRs) along with 231 genes up- and 91 genes downregulated in all samples, suggesting a potential association with tumor development. We further identified 39 hypomethylated promoters associated with increased gene expression in the majority of PTCLs. This putative oncogenic signature included the TRIP13 (thyroid hormone receptor interactor 13) gene whose genetic and pharmacologic inactivation inhibited the proliferation of T-cell lines by inducing G2-M arrest and apoptosis. Our data thus show that human PTCLs have a significant number of recurrent methylation alterations that may affect the expression of genes critical for proliferation whose targeting might be beneficial in anti-lymphoma treatments.

Methylome analysis of endothelial cells suggests new insights on sporadic brain arteriovenous malformation

Arteriovenous malformation of the brain (bAVM) is a vascular phenotype related to brain defective angiogenesis. Involved vessels show impaired expression of vascular differentiation markers resulting in the arteriolar to venule direct shunt. In order to clarify aberrant gene expression occurring in bAVM, here we describe results obtained by methylome analysis performed on endothelial cells (ECs) isolated from bAVM specimens, compared to human cerebral microvascular ECs. Results were validated by quantitative methylation-specific PCR and quantitative realtime-PCR. Differential methylation events occur in genes already linked to bAVM onset, as RBPJ and KRAS. However, among differentially methylated genes, we identified EPHB1 and several other loci involved in EC adhesion as well as in EC/vascular smooth muscle cell (VSMC) crosstalk, suggesting that only endothelial dysfunction might not be sufficient to trigger the bAVM phenotype. Moreover, aberrant methylation pattern was reported for many lncRNA genes targeting transcription factors expressed during neurovascular development. Among these, the YBX1 that was recently shown to target the arteridin coding gene. Finally, in addition to the conventional CpG methylation, we further considered the role of impaired CHG methylation, mainly occurring in brain at embryo stage. We showed as differentially CHG methylated genes are clustered in pathways related to EC homeostasis, as well as to VSMC-EC crosstalk, suggesting as impairment of this interaction plays a prominent role in loss of vascular differentiation, in bAVM phenotype.

Maternal Obesity Induces the Meiotic Defects and Epigenetic Alterations During Fetal Oocyte Development

It has been widely reported that obesity adversely impacts reproductive performance of females. However, the effects of maternal obesity on fetal germ cells remain poorly understood. In the present study, by employing a high-fat diet (HFD)-based mouse model, it is discovered that maternal obesity disrupts the chromosomal synapsis and homologous recombination during fetal oogenesis. Moreover, transcriptomic profiling reveales the potential molecular network controlling this process. Of note, the global hypermethylation of genomic DNA in fetal oocytes from obese mouse is detected. Importantly, time-restricted feeding (TRF) of obese mice not only ameliorate the meiotic defects, but also partly restore the epigenetic remodeling in fetal oocytes. In sum, the evidence are provided showing the deficit fetal oogenesis in obese mother, implicating a mechanism underlying the intergenerational effects of environmental insults. TRF may represent a potentially effective approach for mitigating fertility issues in obese patients.

The genomic mosaic of mitochondrial dysfunction: Decoding nuclear and mitochondrial epigenetic contributions to maternally inherited diabetes and deafness pathogenesis

Aims

Maternally inherited diabetes and deafness (MIDD) is a complex disorder characterized by multiorgan clinical manifestations, including diabetes, hearing loss, and ophthalmic complications. This pilot study aimed to elucidate the intricate interplay between nuclear and mitochondrial genetics, epigenetic modifications, and their potential implications in the pathogenesis of MIDD.

Main methods

A comprehensive genomic approach was employed to analyze a Sicilian family affected by clinically characterized MIDD, negative to the only known causative m.3243 A > G variant, integrating whole-exome sequencing and whole-genome bisulfite sequencing of both nuclear and mitochondrial analyses.

Key findings

Rare and deleterious variants were identified across multiple nuclear genes involved in retinal homeostasis, mitochondrial function, and epigenetic regulation, while complementary mitochondrial DNA analysis revealed a rich tapestry of genetic diversity across genes encoding components of the electron transport chain and ATP synthesis machinery. Epigenetic analyses uncovered significant differentially methylated regions across the genome and within the mitochondrial genome, suggesting a nuanced landscape of epigenetic modulation.

Significance

The integration of genetic and epigenetic data highlighted the potential crosstalk between nuclear and mitochondrial regulation, with specific mtDNA variants influencing methylation patterns and potentially impacting the expression and regulation of mitochondrial genes. This pilot study provides valuable insights into the complex molecular mechanisms underlying MIDD, emphasizing the interplay between nucleus and mitochondrion, tracing the way for future research into targeted therapeutic interventions and personalized approaches for disease management.

Whole-genome bisulfite sequencing data analysis learning module on Google Cloud Platform

This study describes the development of a resource module that is part of a learning platform named 'NIGMS Sandbox for Cloud-based Learning' https://github.com/NIGMS/NIGMS-Sandbox. The overall genesis of the Sandbox is described in the editorial NIGMS Sandbox at the beginning of this Supplement. This module is designed to facilitate interactive learning of whole-genome bisulfite sequencing (WGBS) data analysis utilizing cloud-based tools in Google Cloud Platform, such as Cloud Storage, Vertex AI notebooks and Google Batch. WGBS is a powerful technique that can provide comprehensive insights into DNA methylation patterns at single cytosine resolution, essential for understanding epigenetic regulation across the genome. The designed learning module first provides step-by-step tutorials that guide learners through two main stages of WGBS data analysis, preprocessing and the identification of differentially methylated regions. And then, it provides a streamlined workflow and demonstrates how to effectively use it for large datasets given the power of cloud infrastructure. The integration of these interconnected submodules progressively deepens the user's understanding of the WGBS analysis process along with the use of cloud resources. Through this module, we can enhance the accessibility and adoption of cloud computing in epigenomic research, speeding up the advancements in the related field and beyond. This manuscript describes the development of a resource module that is part of a learning platform named ``NIGMS Sandbox for Cloud-based Learning'' https://github.com/NIGMS/NIGMS-Sandbox. The overall genesis of the Sandbox is described in the editorial NIGMS Sandbox [1] at the beginning of this Supplement. This module delivers learning materials on the analysis of bulk and single-cell ATAC-seq data in an interactive format that uses appropriate cloud resources for data access and analyses.

A multiomic atlas of the aging hippocampus reveals molecular changes in response to environmental enrichment

Aging involves the deterioration of organismal function, leading to the emergence of multiple pathologies. Environmental stimuli, including lifestyle, can influence the trajectory of this process and may be used as tools in the pursuit of healthy aging. To evaluate the role of epigenetic mechanisms in this context, we have generated bulk tissue and single cell multi-omic maps of the male mouse dorsal hippocampus in young and old animals exposed to environmental stimulation in the form of enriched environments. We present a molecular atlas of the aging process, highlighting two distinct axes, related to inflammation and to the dysregulation of mRNA metabolism, at the functional RNA and protein level. Additionally, we report the alteration of heterochromatin domains, including the loss of bivalent chromatin and the uncovering of a heterochromatin-switch phenomenon whereby constitutive heterochromatin loss is partially mitigated through gains in facultative heterochromatin. Notably, we observed the multi-omic reversal of a great number of aging-associated alterations in the context of environmental enrichment, which was particularly linked to glial and oligodendrocyte pathways. In conclusion, our work describes the epigenomic landscape of environmental stimulation in the context of aging and reveals how lifestyle intervention can lead to the multi-layered reversal of aging-associated decline.

Clonal hematopoiesis driven by mutated DNMT3A promotes inflammatory bone loss

Clonal hematopoiesis of indeterminate potential (CHIP) arises from aging-associated acquired mutations in hematopoietic progenitors, which display clonal expansion and produce phenotypically altered leukocytes. We associated CHIP-DNMT3A mutations with a higher prevalence of periodontitis and gingival inflammation among 4,946 community-dwelling adults. To model DNMT3A-driven CHIP, we used mice with the heterozygous loss-of-function mutation R878H, equivalent to the human hotspot mutation R882H. Partial transplantation with Dnmt3aR878H/+ bone marrow (BM) cells resulted in clonal expansion of mutant cells into both myeloid and lymphoid lineages and an elevated abundance of osteoclast precursors in the BM and osteoclastogenic macrophages in the periphery. DNMT3A-driven clonal hematopoiesis in recipient mice promoted naturally occurring periodontitis and aggravated experimentally induced periodontitis and arthritis, associated with enhanced osteoclastogenesis, IL-17-dependent inflammation and neutrophil responses, and impaired regulatory T cell immunosuppressive activity. DNMT3A-driven clonal hematopoiesis and, subsequently, periodontitis were suppressed by rapamycin treatment. DNMT3A-driven CHIP represents a treatable state of maladaptive hematopoiesis promoting inflammatory bone loss.

ChIP-Atlas 3.0: a data-mining suite to explore chromosome architecture together with large-scale regulome data

ChIP-Atlas (https://chip-atlas.org/) presents a suite of data-mining tools for analyzing epigenomic landscapes, powered by the comprehensive integration of over 376 000 public ChIP-seq, ATAC-seq, DNase-seq and Bisulfite-seq experiments from six representative model organisms. To unravel the intricacies of chromatin architecture that mediates the regulome-initiated generation of transcriptional and phenotypic diversity within cells, we report ChIP-Atlas 3.0 that enhances clarity by incorporating additional tracks for genomic and epigenomic features within a newly consolidated 'annotation track' section. The tracks include chromosomal conformation (Hi-C and eQTL datasets), transcriptional regulatory elements (ChromHMM and FANTOM5 enhancers), and genomic variants associated with diseases and phenotypes (GWAS SNPs and ClinVar variants). These annotation tracks are easily accessible alongside other experimental tracks, facilitating better elucidation of chromatin architecture underlying the diversification of transcriptional and phenotypic traits. Furthermore, 'Diff Analysis,' a new online tool, compares the query epigenome data to identify differentially bound, accessible, and methylated regions using ChIP-seq, ATAC-seq and DNase-seq, and Bisulfite-seq datasets, respectively. The integration of annotation tracks and the Diff Analysis tool, coupled with continuous data expansion, renders ChIP-Atlas 3.0 a robust resource for mining the landscape of transcriptional regulatory mechanisms, thereby offering valuable perspectives, particularly for genetic disease research and drug discovery.

Genetic and epigenetic signatures for improved breeding of cultivated blueberry

Blueberry belongs to the Vaccinium genus and is a highly popular fruit crop with significant economic importance. It was not until the early twentieth century that they began to be domesticated through extensive interspecific hybridization. Here, we collected 220 Vaccinium accessions from various geographical locations, including 154 from the United States, 14 from China, eight from Australia, and 29 from Europe and other countries, comprising 164 Vaccinium corymbosum, 15 Vaccinium ashei, 10 lowbush blueberries, seven half-high blueberries, and others. We present the whole-genome variation map of 220 accessions and reconstructed the hundred-year molecular history of interspecific hybridization of blueberry. We focused on the two major blueberry subgroups, the northern highbush blueberry (NHB) and southern highbush blueberry (SHB) and identified candidate genes that contribute to their distinct traits in climate adaptability and fruit quality. Our analysis unveiled the role of gene introgression from Vaccinium darrowii and V. ashei into SHB in driving the differentiation between SHB and NHB, potentially facilitating SHB's adaptation to subtropical environments. Assisted by genome-wide association studies, our analysis suggested VcTBL44 as a pivotal gene regulator governing fruit firmness in SHB. Additionally, we conducted whole-genome bisulfite sequencing on nine NHB and 12 SHB cultivars, and characterized regions that are differentially methylated between the two subgroups. In particular, we discovered that the β-alanine metabolic pathway genes were enriched for DNA methylation changes. Our study provides high-quality genetic and epigenetic variation maps for blueberry, which offer valuable insights and resources for future blueberry breeding.

Single-stranded pre-methylated 5mC adapters uncover the methylation profile of plasma ultrashort Single-stranded cell-free DNA

Whole-genome bisulfite sequencing (BS-Seq) measures cytosine methylation changes at single-base resolution and can be used to profile cell-free DNA (cfDNA). In plasma, ultrashort single-stranded cfDNA (uscfDNA, ∼50 nt) has been identified together with 167 bp double-stranded mononucleosomal cell-free DNA (mncfDNA). However, the methylation profile of uscfDNA has not been described. Conventional BS-Seq workflows may not be helpful because bisulfite conversion degrades larger DNA into smaller fragments, leading to erroneous categorization as uscfDNA. We describe the '5mCAdpBS-Seq' workflow in which pre-methylated 5mC (5-methylcytosine) single-stranded adapters are ligated to heat-denatured cfDNA before bisulfite conversion. This method retains only DNA fragments that are unaltered by bisulfite treatment, resulting in less biased uscfDNA methylation analysis. Using 5mCAdpBS-Seq, uscfDNA had lower levels of DNA methylation (∼15%) compared to mncfDNA and was enriched in promoters and CpG islands. Hypomethylated uscfDNA fragments were enriched in upstream transcription start sites (TSSs), and the intensity of enrichment was correlated with expressed genes of hemopoietic cells. Using tissue-of-origin deconvolution, we inferred that uscfDNA is derived primarily from eosinophils, neutrophils, and monocytes. As proof-of-principle, we show that characteristics of the methylation profile of uscfDNA can distinguish non-small cell lung carcinoma from non-cancer samples. The 5mCAdpBS-Seq workflow is recommended for any cfDNA methylation-based investigations.

DNA methylation patterns in patients with asthenospermia and oligoasthenospermia

BACKGROUND

Spermatogenesis is a highly regulated and complex process in which DNA methylation plays a crucial role. This study aimed to explore the differential methylation profiles in sperm DNA between patients with asthenospermia (AS) and healthy controls (HCs), those with oligoasthenospermia (OAS) and HCs, and patients with AS and those with OAS.

RESULTS

Semen samples and clinical data were collected from five patients with AS, five patients with OAS, and six age-matched HCs. Reduced representation bisulfite sequencing (RRBS) was performed to identify differentially methylated regions (DMRs) in sperm cells among the different types of patients and HCs. A total of 6520, 28,019, and 16,432 DMRs were detected between AS and HC, OAS and HC, and AS and OAS groups, respectively. These DMRs were predominantly located within gene bodies and mapped to 2868, 9296, and 9090 genes in the respective groups. Of note, 12, 9, and 8 DMRs in each group were closely associated with spermatogenesis and male infertility. Furthermore, BDNF, SMARCB1, PIK3CA, and DDX27; RBMX and SPATA17; ASZ1, CDH1, and CHDH were identified as strong differentially methylated candidate genes in each group, respectively. Meanwhile, the GO analysis of DMR-associated genes in the AS vs. HC groups revealed that protein binding, cytoplasm, and transcription (DNA-templated) were the most enriched terms in the biological process (BP), cellular component (CC), and molecular function (MF), respectively. Likewise, in both the OAS vs. HC and AS vs. OAS groups, GO analysis revealed protein binding, nucleus, and transcription (DNA-templated) as the most enriched terms in BP, CC, and MF, respectively. Finally, the KEGG analysis of DMR-annotated genes and these genes at promoters suggested that metabolic pathways were the most significantly associated across all three groups.

CONCLUSIONS

The current study results revealed distinctive sperm DNA methylation patterns in the AS vs. HC and OAS vs. HC groups, particularly between patients with AS and those with OAS. The identification of key genes associated with spermatogenesis and male infertility in addition to the differentially enriched metabolic pathways may contribute to uncovering the potential pathogenesis in different types of abnormal sperm parameters.

Endocrine disruptor-induced epimutagenesis in vitro : Insight into molecular mechanisms

Endocrine disrupting chemicals (EDCs) such as bisphenol S (BPS) are xenobiotic compounds that can disrupt endocrine signaling following exposure due to steric similarities to endogenous hormones within the body. EDCs have been shown to induce disruptions in normal epigenetic programming (epimutations) that accompany dysregulation of normal gene expression patterns that appear to predispose disease states. Most interestingly, the prevalence of epimutations following exposure to many different EDCs often persists over multiple subsequent generations, even with no further exposure to the causative EDC. Many previous studies have described both the direct and prolonged effects of EDC exposure in animal models, but many questions remain about molecular mechanisms by which EDCs initially induce epimutations or contribute to the propagation of EDC-induced epimutations either within the exposed generation or to subsequent generations. Additional questions remain regarding the extent to which there may be differences in cell-type specific susceptibilities to various EDCs, and whether this susceptibility is correlative with expression of relevant hormone receptors and/or the location of relevant hormone response elements (HREs) in the genome. To address these questions, we exposed cultured mouse pluripotent (induced pluripotent stem [iPS]), somatic (Sertoli and granulosa), and germ (primordial germ cell like [PGCLC]) cells to BPS and measured changes in DNA methylation levels at the epigenomic level and gene expression at the transcriptomic level. We found that there was indeed a difference in cell-type specific susceptibility to EDC-induced epimutagenesis and that this susceptibility correlated with differential expression of relevant hormone receptors and, in many cases, tended to generate epimutations near relevant HREs within the genome. Additionally, however, we also found that BPS can induce epimutations in a cell type that does not express relevant receptors and in genomic regions that do not contain relevant HREs, suggesting that both canonical and non-canonical signaling mechanisms can be disrupted by BPS exposure. Most interestingly, we found that when iPS cells were exposed to BPS and then induced to differentiate into PGCLCs, the prevalence of epimutations and differentially expressed genes (DEGs) initially induced in the iPSCs was largely retained in the resulting PGCLCs, however, >90% of the specific epimutations and DEGs were not conserved but were rather replaced by novel epimutations and DEGs following the iPSC to PGCLC transition. These results are consistent with a unique concept that many EDC-induced epimutations may normally be corrected by germline and/or embryonic epigenetic reprogramming but that due to disruption of the underlying chromatin architecture induced by the EDC exposure, many novel epimutations may emerge during the reprogramming process as well. Thus, it appears that following exposure to a disruptive agent such as an EDC, a prevalence of epimutations may transcend epigenetic reprogramming even though most individual epimutations are not conserved during this process.

Genome-wide DNA methylation changes in human spermatogenesis

Sperm production and function require the correct establishment of DNA methylation patterns in the germline. Here, we examined the genome-wide DNA methylation changes during human spermatogenesis and its alterations in disturbed spermatogenesis. We found that spermatogenesis is associated with remodeling of the methylome, comprising a global decline in DNA methylation in primary spermatocytes followed by selective remethylation, resulting in a spermatids/sperm-specific methylome. Hypomethylated regions in spermatids/sperm were enriched in specific transcription factor binding sites for DMRT and SOX family members and spermatid-specific genes. Intriguingly, while SINEs displayed differential methylation throughout spermatogenesis, LINEs appeared to be protected from changes in DNA methylation. In disturbed spermatogenesis, germ cells exhibited considerable DNA methylation changes, which were significantly enriched at transposable elements and genes involved in spermatogenesis. We detected hypomethylation in SVA and L1HS in disturbed spermatogenesis, suggesting an association between the abnormal programming of these regions and failure of germ cells progressing beyond meiosis.

Effects of Developmental Lead and Phthalate Exposures on DNA Methylation in Adult Mouse Blood, Brain, and Liver: A Focus on Genomic Imprinting by Tissue and Sex

BACKGROUND

Maternal exposure to environmental chemicals can cause adverse health effects in offspring. Mounting evidence supports that these effects are influenced, at least in part, by epigenetic modifications. It is unknown whether epigenetic changes in surrogate tissues such as the blood are reflective of similar changes in target tissues such as cortex or liver.

OBJECTIVE

We examined tissue- and sex-specific changes in DNA methylation (DNAm) associated with human-relevant lead (Pb) and di(2-ethylhexyl) phthalate (DEHP) exposure during perinatal development in cerebral cortex, blood, and liver.

METHODS

Female mice were exposed to human relevant doses of either Pb (32  ppm ) via drinking water or DEHP (5 mg / kg-day ) via chow for 2 weeks prior to mating through offspring weaning. Whole genome bisulfite sequencing (WGBS) was utilized to examine DNAm changes in offspring cortex, blood, and liver at 5 months of age. Metilene and methylSig were used to identify differentially methylated regions (DMRs). Annotatr and ChIP-enrich were used for genomic annotations and gene set enrichment tests of DMRs, respectively.

RESULTS

The cortex contained the majority of DMRs associated with Pb (66%) and DEHP (57%) exposure. The cortex also contained the greatest degree of overlap in DMR signatures between sexes (n = 13 and 8 DMRs with Pb and DEHP exposure, respectively) and exposure types (n = 55 and 39 DMRs in males and females, respectively). In all tissues, detected DMRs were preferentially found at genomic regions associated with gene expression regulation (e.g., CpG islands and shores, 5' UTRs, promoters, and exons). An analysis of GO terms associated with DMR-containing genes identified imprinted genes to be impacted by both Pb and DEHP exposure. Of these, Gnas and Grb10 contained DMRs across tissues, sexes, and exposures, with some signatures replicated between target and surrogate tissues. DMRs were enriched in the imprinting control regions (ICRs) of Gnas and Grb10, and we again observed a replication of DMR signatures between blood and target tissues. Specifically, we observed hypermethylation of the Grb10 ICR in both blood and liver of Pb-exposed male animals.

CONCLUSIONS

These data provide preliminary evidence that imprinted genes may be viable candidates in the search for epigenetic biomarkers of toxicant exposure in target tissues. Additional research is needed on allele- and developmental stage-specific effects, as well as whether other imprinted genes provide additional examples of this relationship. https://doi.org/10.1289/EHP14074.

Extraembryonic gut endoderm cells undergo programmed cell death during development

Despite a distinct developmental origin, extraembryonic cells in mice contribute to gut endoderm and converge to transcriptionally resemble their embryonic counterparts. Notably, all extraembryonic progenitors share a non-canonical epigenome, raising several pertinent questions, including whether this landscape is reset to match the embryonic regulation and if extraembryonic cells persist into later development. Here we developed a two-colour lineage-tracing strategy to track and isolate extraembryonic cells over time. We find that extraembryonic gut cells display substantial memory of their developmental origin including retention of the original DNA methylation landscape and resulting transcriptional signatures. Furthermore, we show that extraembryonic gut cells undergo programmed cell death and neighbouring embryonic cells clear their remnants via non-professional phagocytosis. By midgestation, we no longer detect extraembryonic cells in the wild-type gut, whereas they persist and differentiate further in p53-mutant embryos. Our study provides key insights into the molecular and developmental fate of extraembryonic cells inside the embryo.

A comprehensive analysis of minimally differentially methylated regions common to pediatric and adult solid tumors

Cancer is the second most common cause of death in children aged 1-14 years in the United States, with 11,000 new cases and 1200 deaths annually. Pediatric cancers typically have lower mutational burden compared to adult-onset cancers, however, the epigenomes in pediatric cancer are highly altered, with widespread DNA methylation changes. The rarity of pediatric cancers poses a significant challenge to developing cancer-type specific biomarkers for diagnosis, prognosis, or treatment monitoring. In the current study, we explored the potential of a DNA methylation profile common across various pediatric cancers. To do this, we conducted whole genome bisulfite sequencing (WGBS) on 31 recurrent pediatric tumor tissues, 13 normal tissues, and 20 plasma cell-free (cf)DNA samples, representing 11 different pediatric cancer types. We defined minimal focal regions that were differentially methylated across samples in the multiple cancer types which we termed minimally differentially methylated regions (mDMRs). These methylation changes were also observed in 506 pediatric and 5691 adult cancer samples accessed from publicly available databases, and in 44 pediatric cancer samples we analyzed using a targeted hybridization probe capture assay. Finally, we found that these methylation changes were detectable in cfDNA and could serve as potential cfDNA methylation biomarkers for early detection or minimal residual disease.

Wemics: A Single-Base Resolution Methylation Quantification Method for Enhanced Prediction of Epigenetic Regulation

DNA methylation, an epigenetic mechanism that alters gene expression without changing DNA sequence, is essential for organism development and key biological processes like genomic imprinting and X-chromosome inactivation. Despite tremendous efforts in DNA methylation research, accurate quantification of cytosine methylation remains a challenge. Here, a single-base methylation quantification approach is introduced by weighting methylation of consecutive CpG sites (Wemics) in genomic regions. Wemics quantification of DNA methylation better predicts its regulatory impact on gene transcription and identifies differentially methylated regions (DMRs) with more biological relevance. Most Wemics-quantified DMRs in lung cancer are epigenetically conserved and recurrently occurred in other primary cancers from The Cancer Genome Atlas (TCGA), and their aberrant alterations can serve as promising pan-cancer diagnostic markers. It is further revealed that these detected DMRs are enriched in transcription factor (TF) binding motifs, and methylation of these TF binding motifs and TF expression synergistically regulate target gene expression. Using Wemics on epigenomic-transcriptomic data from the large lung cancer cohort, a dozen novel genes with oncogenic potential are discovered that are upregulated by hypomethylation but overlooked by other quantification methods. These findings increase the understanding of the epigenetic mechanism by which DNA methylation regulates gene expression.

Exploring the potential of Paris polyphylla var. yunnanensis pollen manipulation in modifying seed dormancy

Paris polyphylla var. yunnanensis, a well-known Chinese medicinal herb, shows a unique physiological trait characterized by the cyclic opening and closing of its anthers after pollen maturation. The aim of this study was to explore the implications of this phenomenon on breeding. RNA sequencing coupled with methylation sequencing was used to scrutinize and compare gene expression profiles and methylation alterations in pollen and seeds during anther opening and closing, along with cold exposure. Genes enriched within Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were examined to identify gene clusters susceptible to temperature-related methylation changes in both pollen and seeds. Four pollen treatment models, namely, normal control, "pollen protected from low temperatures," "pollen from just-opened anther," and "pollen from close-blocked anther," were used to produce corresponding seeds via artificial pollination. Subsequently, qRT-PCR was used to validate modifications in the expression patterns of marker genes in pollinated seeds under diverse treatment scenarios. Genes exhibiting significant differences in expression between anthers and normal tissues, along with gene regions linked to methylation variations attributed to low-temperature-treated pollen and seeds, were identified through transcriptomic analysis. Convergence was observed in three signaling pathways: oxidative phosphorylation (ko00190), plant hormone signal transduction (Ko04075), and zeatin biosynthesis (ko00908). Notably, gene clusters prone to temperature-induced methylation changes, such as NADH-ubiquinone oxidoreductase chain 5, plasma membrane ATPase 4, cytochrome c oxidase subunit 2, cis-zeatin O-glucosyltransferase, ABSCISIC ACID-INSENSITIVE 5-like protein 4, and indole-3-acetic acid-amido synthetase (IAAS), were identified. Evaluation using various pollen pollination models revealed altered expression patterns of five dormancy-regulating marker genes: IAAS, sucrose synthase (SUS), gibberellin 2-oxidase (GA2ox), ABA INSENSITIVE 2 (ABI2), and auxin-repressed protein (ARP), in seeds pollinated with pollen from close-blocked anthers, cold-protected pollen, and pollen from freshly opened anthers. The close-blocked anther treatment led to significantly upregulated expression of IAAS, SUS, GA2ox, and ABI2, whereas ARP expression decreased markedly, indicating a propensity toward prolonged seed dormancy. Conversely, in the low-temperature-protected anther model, SUS, ARP, GA2ox, and IAAS exhibited reduced expression levels, whereas the expression of ABI2 was upregulated, overall facilitating seed germination.

A platform for rapid patient-derived cutaneous neurofibroma organoid establishment and screening

Localized cutaneous neurofibromas (cNFs) are benign tumors that arise in the dermis of patients affected by neurofibromatosis type 1 syndrome. cNFs are benign lesions: they do not undergo malignant transformation or metastasize. Nevertheless, they can cover a significant proportion of the body, with some individuals developing hundreds to thousands of lesions. cNFs can cause pain, itching, and disfigurement resulting in substantial socio-emotional repercussions. Currently, surgery and laser desiccation are the sole treatment options but may result in scarring and potential regrowth from incomplete removal. To identify effective systemic therapies, we introduce an approach to establish and screen cNF organoids. We optimized conditions to support the ex vivo growth of genomically diverse cNFs. Patient-derived cNF organoids closely recapitulate cellular and molecular features of parental tumors as measured by immunohistopathology, methylation, RNA sequencing, and flow cytometry. Our cNF organoid platform enables rapid screening of hundreds of compounds in a patient- and tumor-specific manner.

Epigenomic and transcriptomic persistence of heat stress memory in strawberry (Fragaria vesca)

BACKGROUND

In plants, epigenetic stress memory has so far been found to be largely transient. Here, we wanted to assess the heritability of heat stress-induced epigenetic and transcriptomic changes following woodland strawberry (Fragaria vesca) reproduction. Strawberry is an ideal model to study epigenetic inheritance because it presents two modes of reproduction: sexual (self-pollinated plants) and asexual (clonally propagated plants named daughter plants). Taking advantage of this model, we investigated whether heat stress-induced DNA methylation changes can be transmitted via asexual reproduction.

RESULTS

Our genome-wide study provides evidence for stress memory acquisition and maintenance in F. vesca. We found that specific DNA methylation marks or epimutations are stably transmitted over at least three asexual generations. Some of the epimutations were associated with transcriptional changes after heat stress.

CONCLUSION

Our findings show that the strawberry methylome and transcriptome respond with a high level of flexibility to heat stress. Notably, independent plants acquired the same epimutations and those were inherited by their asexual progenies. Overall, the asexual progenies can retain some information in the genome of past stresses encountered by their progenitors. This molecular memory, also documented at the transcriptional level, might be involved in functional plasticity and stress adaptation. Finally, these findings may contribute to novel breeding approaches for climate-ready plants.

Novel genome-wide DNA methylation profiling reveals distinct epigenetic landscape, prognostic model and cellular composition of early-stage lung adenocarcinoma

BACKGROUND

Lung adenocarcinoma (LUAD) has been a leading cause of cancer-related mortality worldwide. Early intervention can significantly improve prognosis. DNA methylation could occur in the early stage of tumor. Comprehensive understanding the epigenetic landscape of early-stage LUAD is crucial in understanding tumorigenesis.

METHODS

Enzymatic methyl sequencing (EM-seq) was performed on 23 tumors and paired normal tissue to reveal distinct epigenetic landscape, for compared with The Cancer Genome Atlas (TCGA) 450K methylation microarray data. Then, an integrative analysis was performed combined with TCGA LUAD RNA-seq data to identify significant differential methylated and expressed genes. Subsequently, the prognostic risk model was constructed and cellular composition was analyzed.

RESULTS

Methylome analysis of EM-seq comparing tumor and normal tissues identified 25 million cytosine-phosphate-guanine (CpG) sites and 30,187 differentially methylated regions (DMR) with a greater number of untraditional types. EM-seq identified a significantly higher number of CpG sites and DMRs compared to the 450K microarray. By integrating the differentially methylated genes (DMGs) with LUAD-related differentially expressed genes (DEGs) from the TCGA database, we constructed prognostic model based on six differentially methylated-expressed genes (MEGs) and verified our prognostic model in GSE13213 and GSE42127 dataset. Finally, cell deconvolution based on the in-house EM-seq methylation profile was used to estimate cellular composition of early-stage LUAD.

CONCLUSIONS

This study firstly delves into novel pattern of epigenomic DNA methylation and provides a multidimensional analysis of the role of DNA methylation revealed by EM-seq in early-stage LUAD, providing distinctive insights into its potential epigenetic mechanisms.

Nonlinear DNA methylation trajectories in aging male mice

Although DNA methylation data yields highly accurate age predictors, little is known about the dynamics of this quintessential epigenomic biomarker during lifespan. To narrow the gap, we investigate the methylation trajectories of male mouse colon at five different time points of aging. Our study indicates the existence of sudden hypermethylation events at specific stages of life. Precisely, we identify two epigenomic switches during early-to-midlife (3-9 months) and mid-to-late-life (15-24 months) transitions, separating the rodents' life into three stages. These nonlinear methylation dynamics predominantly affect genes associated with the nervous system and enrich in bivalently marked chromatin regions. Based on groups of nonlinearly modified loci, we construct a clock-like classifier STageR (STage of aging estimatoR) that accurately predicts murine epigenetic stage. We demonstrate the universality of our clock in an independent mouse cohort and with publicly available datasets.

Increased DNA methylation contributes to the early ripening of pear fruits during domestication and improvement

BACKGROUND

DNA methylation is an essential epigenetic modification. However, its contribution to trait changes and diversity in the domestication of perennial fruit trees remains unknown.

RESULTS

Here, we investigate the variation in DNA methylation during pear domestication and improvement using whole-genome bisulfite sequencing in 41 pear accessions. Contrary to the significant decrease during rice domestication, we detect a global increase in DNA methylation during pear domestication and improvement. We find this specific increase in pear is significantly correlated with the downregulation of Demeter-like1 (DML1, encoding DNA demethylase) due to human selection. We identify a total of 5591 differentially methylated regions (DMRs). Methylation in the CG and CHG contexts undergoes co-evolution during pear domestication and improvement. DMRs have higher genetic diversity than selection sweep regions, especially in the introns. Approximately 97% of DMRs are not associated with any SNPs, and these DMRs are associated with starch and sucrose metabolism and phenylpropanoid biosynthesis. We also perform correlation analysis between DNA methylation and gene expression. We find genes close to the hypermethylated DMRs that are significantly associated with fruit ripening. We further verify the function of a hyper-DMR-associated gene, CAMTA2, and demonstrate that overexpression of CAMTA2 in tomato and pear callus inhibits fruit ripening.

CONCLUSIONS

Our study describes a specific pattern of DNA methylation in the domestication and improvement of a perennial pear tree and suggests that increased DNA methylation plays an essential role in the early ripening of pear fruits.

Blood leukocytes as a non-invasive diagnostic tool for thyroid nodules: a prospective cohort study

BACKGROUND

Thyroid nodule (TN) patients in China are subject to overdiagnosis and overtreatment. The implementation of existing technologies such as thyroid ultrasonography has indeed contributed to the improved diagnostic accuracy of TNs. However, a significant issue persists, where many patients undergo unnecessary biopsies, and patients with malignant thyroid nodules (MTNs) are advised to undergo surgery therapy.

METHODS

This study included a total of 293 patients diagnosed with TNs. Differential methylation haplotype blocks (MHBs) in blood leukocytes between MTNs and benign thyroid nodules (BTNs) were detected using reduced representation bisulfite sequencing (RRBS). Subsequently, an artificial intelligence blood leukocyte DNA methylation (BLDM) model was designed to optimize the management and treatment of patients with TNs for more effective outcomes.

RESULTS

The DNA methylation profiles of peripheral blood leukocytes exhibited distinctions between MTNs and BTNs. The BLDM model we developed for diagnosing TNs achieved an area under the curve (AUC) of 0.858 in the validation cohort and 0.863 in the independent test cohort. Its specificity reached 90.91% and 88.68% in the validation and independent test cohorts, respectively, outperforming the specificity of ultrasonography (43.64% in the validation cohort and 47.17% in the independent test cohort), albeit with a slightly lower sensitivity (83.33% in the validation cohort and 82.86% in the independent test cohort) compared to ultrasonography (97.62% in the validation cohort and 100.00% in the independent test cohort). The BLDM model could correctly identify 89.83% patients whose nodules were suspected malignant by ultrasonography but finally histological benign. In micronodules, the model displayed higher specificity (93.33% in the validation cohort and 92.00% in the independent test cohort) and accuracy (88.24% in the validation cohort and 87.50% in the independent test cohort) for diagnosing TNs. This performance surpassed the specificity and accuracy observed with ultrasonography. A TN diagnostic and treatment framework that prioritizes patients is provided, with fine-needle aspiration (FNA) biopsy performed only on patients with indications of MTNs in both BLDM and ultrasonography results, thus avoiding unnecessary biopsies.

CONCLUSIONS

This is the first study to demonstrate the potential of non-invasive blood leukocytes in diagnosing TNs, thereby making TN diagnosis and treatment more efficient in China.

FinaleMe: Predicting DNA methylation by the fragmentation patterns of plasma cell-free DNA

Analysis of DNA methylation in cell-free DNA reveals clinically relevant biomarkers but requires specialized protocols such as whole-genome bisulfite sequencing. Meanwhile, millions of cell-free DNA samples are being profiled by whole-genome sequencing. Here, we develop FinaleMe, a non-homogeneous Hidden Markov Model, to predict DNA methylation of cell-free DNA and, therefore, tissues-of-origin, directly from plasma whole-genome sequencing. We validate the performance with 80 pairs of deep and shallow-coverage whole-genome sequencing and whole-genome bisulfite sequencing data.

HBCR_DMR: A Hybrid Method Based on Beta-Binomial Bayesian Hierarchical Model and Combination of Ranking Method to Detect Differential Methylation Regions in Bisulfite Sequencing Data

DNA methylation is a key epigenetic modification involved in gene regulation, contributing to both physiological and pathological conditions. For a more profound comprehension, it is essential to conduct a precise comparison of DNA methylation patterns between sample groups that represent distinct statuses. Analysis of differentially methylated regions (DMRs) using computational approaches can help uncover the precise relationships between these phenomena. This paper describes a hybrid model that combines the beta-binomial Bayesian hierarchical model with a combination of ranking methods known as HBCR_DMR. During the initial phase, we model the actual methylation proportions of the CpG sites (CpGs) within the replicates. This modeling is achieved through beta-binomial distribution, with parameters set by a group mean and a dispersion parameter. During the second stage, we establish the selection of distinguishing CpG sites based on their methylation status, employing multiple ranking techniques. Finally, we combine the ranking lists of differentially methylated CpG sites through a voting system. Our analyses, encompassing simulations and real data, reveal outstanding performance metrics, including a sensitivity of 0.72, specificity of 0.89, and an F1 score of 0.76, yielding an overall accuracy of 0.82 and an AUC of 0.94. These findings underscore HBCR_DMR's robust capacity to distinguish methylated regions, confirming its utility as a valuable tool for DNA methylation analysis.

Heritable changes of epialleles near genes in maize can be triggered in the absence of CHH methylation

Trans-chromosomal interactions resulting in changes in DNA methylation during hybridization have been observed in several plant species. However, little is known about the causes or consequences of these interactions. Here, we compared DNA methylomes of F1 hybrids that are mutant for a small RNA biogenesis gene, Mop1 (Mediator of paramutation1), with that of their parents, wild-type siblings, and backcrossed progeny in maize (Zea mays). Our data show that hybridization triggers global changes in both trans-chromosomal methylation (TCM) and trans-chromosomal demethylation (TCdM), most of which involved changes in CHH methylation. In more than 60% of these TCM differentially methylated regions (DMRs) in which small RNAs are available, no significant changes in the quantity of small RNAs were observed. Methylation at the CHH TCM DMRs was largely lost in the mop1 mutant, although the effects of this mutant varied depending on the location of these DMRs. Interestingly, an increase in CHH at TCM DMRs was associated with enhanced expression of a subset of highly expressed genes and suppressed expression of a small number of lowly expressed genes. Examination of the methylation levels in backcrossed plants demonstrates that both TCM and TCdM can be maintained in the subsequent generation, but that TCdM is more stable than TCM. Surprisingly, although increased CHH methylation in most TCM DMRs in F1 plants required Mop1, initiation of a new epigenetic state of these DMRs did not require a functional copy of this gene, suggesting that initiation of these changes is independent of RNA-directed DNA methylation.

MOF-mediated acetylation of UHRF1 enhances UHRF1 E3 ligase activity to facilitate DNA methylation maintenance

The multi-domain protein UHRF1 (ubiquitin-like, containing PHD and RING finger domains, 1) recruits DNMT1 for DNA methylation maintenance during DNA replication. Here, we show that MOF (males absent on the first) acetylates UHRF1 at K670 in the pre-RING linker region, whereas HDAC1 deacetylates UHRF1 at the same site. We also identify that K667 and K668 can also be acetylated by MOF when K670 is mutated. The MOF/HDAC1-mediated acetylation in UHRF1 is cell-cycle regulated and peaks at G1/S phase, in line with the function of UHRF1 in recruiting DNMT1 to maintain DNA methylation. In addition, UHRF1 acetylation significantly enhances its E3 ligase activity. Abolishing UHRF1 acetylation at these sites attenuates UHRF1-mediated H3 ubiquitination, which in turn impairs DNMT1 recruitment and DNA methylation. Taken together, these findings identify MOF as an acetyltransferase for UHRF1 and define a mechanism underlying the regulation of DNA methylation maintenance through MOF-mediated UHRF1 acetylation.

Methylation changes of liver DNA during the formation of gallstones

Aim: To explore the overall methylation changes in liver tissues during the formation of gallstones, as well as the key pathways and genes involved in the process. Methods: Reduced-representation bisulfite sequencing and RNA sequencing were conducted on the liver tissues of mice with gallstones and control normal mice. Results: A total of 8705 differentially methylated regions in CpG and 1410 differentially expressed genes were identified. The joint analysis indicated that aberrant DNA methylation may be associated with dysregulated gene expression in key pathways such as cholesterol metabolism and bile secretion. Conclusion: We propose for the first time that methylation changes in some key pathway genes in liver tissue may be involved in the formation of gallstones.

BSXplorer: analytical framework for exploratory analysis of BS-seq data

BACKGROUND

Bisulfite sequencing detects and quantifies DNA methylation patterns, contributing to our understanding of gene expression regulation, genome stability maintenance, conservation of epigenetic mechanisms across divergent taxa, epigenetic inheritance and, eventually, phenotypic variation. Graphical representation of methylation data is crucial in exploring epigenetic regulation on a genome-wide scale in both plants and animals. This is especially relevant for non-model organisms with poorly annotated genomes and/or organisms where genome sequences are not yet assembled on chromosome level. Despite being a technology of choice to profile DNA methylation for many years now there are surprisingly few lightweight and robust standalone tools available for efficient graphical analysis of data in non-model systems. This significantly limits evolutionary studies and agrigenomics research. BSXplorer is a tool specifically developed to fill this gap and assist researchers in explorative data analysis and in visualising and interpreting bisulfite sequencing data more easily.

RESULTS

BSXplorer provides in-depth graphical analysis of sequencing data encompassing (a) profiling of methylation levels in metagenes or in user-defined regions using line plots and heatmaps, generation of summary statistics charts, (b) enabling comparative analyses of methylation patterns across experimental samples, methylation contexts and species, and (c) identification of modules sharing similar methylation signatures at functional genomic elements. The tool processes methylation data quickly and offers API and CLI capabilities, along with the ability to create high-quality figures suitable for publication.

CONCLUSIONS

BSXplorer facilitates efficient methylation data mining, contrasting and visualization, making it an easy-to-use package that is highly useful for epigenetic research.

The impact of parental and developmental stress on DNA methylation in the avian hypothalamic-pituitary-adrenal axis

The hypothalamic-pituitary-adrenal (HPA) axis coordinates an organism's response to environmental stress. The responsiveness and sensitivity of an offspring's stress response may be shaped not only by stressors encountered in their early post-natal environment but also by stressors in their parent's environment. Yet, few studies have considered how stressors encountered in both of these early life environments may function together to impact the developing HPA axis. Here, we manipulated stressors in the parental and post-natal environments in a population of house sparrows (Passer domesticus) to assess their impact on changes in DNA methylation (and corresponding gene expression) in a suite of genes within the HPA axis. We found that nestlings that experienced early life stress across both life-history periods had higher DNA methylation in a critical HPA axis gene, the glucocorticoid receptor (NR3C1). In addition, we found that the life-history stage when stress was encountered impacted some genes (HSD11B1, NR3C1 and NR3C2) differently. We also found evidence for the mitigation of parental stress by post-natal stress (in HSD11B1 and NR3C2). Finally, by assessing DNA methylation in both the brain and blood, we were able to evaluate cross-tissue patterns. While some differentially methylated regions were tissue-specific, we found cross-tissue changes in NR3C2 and NR3C1, suggesting that blood is a suitable tissue for assessing DNA methylation as a biomarker of early life stress. Our results provide a crucial first step in understanding the mechanisms by which early life stress in different life-history periods contributes to changes in the epigenome of the HPA axis.

Long-read sequencing unveils novel somatic variants and methylation patterns in the genetic information system of early lung cancer

Lung cancer poses a global health challenge, necessitating advanced diagnostics for improved outcomes. Intensive efforts are ongoing to pinpoint early detection biomarkers, such as genomic variations and DNA methylation, to elevate diagnostic precision. We conducted long-read sequencing on cancerous and adjacent non-cancerous tissues from a patient with lung adenocarcinoma. We identified somatic structural variations (SVs) specific to lung cancer by integrating data from various SV calling methods and differentially methylated regions (DMRs) that were distinct between these two tissue samples, revealing a unique methylation pattern associated with lung cancer. This study discovered over 40,000 somatic SVs and over 180,000 DMRs linked to lung cancer. We identified approximately 700 genes of significant relevance through comprehensive analysis, including genes intricately associated with many lung cancers, such as NOTCH1, SMOC2, CSMD2, and others. Furthermore, we observed that somatic SVs and DMRs were substantially enriched in several pathways, such as axon guidance signaling pathways, which suggests a comprehensive multi-omics impact on lung cancer progression across various biological investigation levels. These datasets can potentially serve as biomarkers for early lung cancer detection and may hold significant value in clinical diagnosis and treatment applications.

Mechanisms of action and resistance in histone methylation-targeted therapy

Epigenomes enable the rectification of disordered cancer gene expression, thereby providing new targets for pharmacological interventions. The clinical utility of targeting histone H3 lysine trimethylation (H3K27me3) as an epigenetic hallmark has been demonstrated1-7. However, in actual therapeutic settings, the mechanism by which H3K27me3-targeting therapies exert their effects and the response of tumour cells remain unclear. Here we show the potency and mechanisms of action and resistance of the EZH1-EZH2 dual inhibitor valemetostat in clinical trials of patients with adult T cell leukaemia/lymphoma. Administration of valemetostat reduced tumour size and demonstrated durable clinical response in aggressive lymphomas with multiple genetic mutations. Integrative single-cell analyses showed that valemetostat abolishes the highly condensed chromatin structure formed by the plastic H3K27me3 and neutralizes multiple gene loci, including tumour suppressor genes. Nevertheless, subsequent long-term treatment encounters the emergence of resistant clones with reconstructed aggregate chromatin that closely resemble the pre-dose state. Acquired mutations at the PRC2-compound interface result in the propagation of clones with increased H3K27me3 expression. In patients free of PRC2 mutations, TET2 mutation or elevated DNMT3A expression causes similar chromatin recondensation through de novo DNA methylation in the H3K27me3-associated regions. We identified subpopulations with distinct metabolic and gene translation characteristics implicated in primary susceptibility until the acquisition of the heritable (epi)mutations. Targeting epigenetic drivers and chromatin homeostasis may provide opportunities for further sustained epigenetic cancer therapies.

Cell-free DNA methylation reveals cell-specific tissue injury and correlates with disease severity and patient outcomes in COVID-19

BACKGROUND

The recently identified methylation patterns specific to cell type allows the tracing of cell death dynamics at the cellular level in health and diseases. This study used COVID-19 as a disease model to investigate the efficacy of cell-specific cell-free DNA (cfDNA) methylation markers in reflecting or predicting disease severity or outcome.

METHODS

Whole genome methylation sequencing of cfDNA was performed for 20 healthy individuals, 20 cases with non-hospitalized COVID-19 and 12 cases with severe COVID-19 admitted to intensive care unit (ICU). Differentially methylated regions (DMRs) and gene ontology pathway enrichment analyses were performed to explore the locus-specific methylation difference between cohorts. The proportion of cfDNA derived from lung and immune cells to a given sample (i.e. tissue fraction) at cell-type resolution was estimated using a novel algorithm, which reflects lung injuries and immune response in COVID-19 patients and was further used to evaluate clinical severity and patient outcome.

RESULTS

COVID‑19 patients had globally reduced cfDNA methylation level compared with healthy controls. Compared with non-hospitalized COVID-19 patients, the cfDNA methylation pattern was significantly altered in severe patients with the identification of 11,156 DMRs, which were mainly enriched in pathways related to immune response. Markedly elevated levels of cfDNA derived from lung and more specifically alveolar epithelial cells, bronchial epithelial cells, and lung endothelial cells were observed in COVID-19 patients compared with healthy controls. Compared with non-hospitalized patients or healthy controls, severe COVID-19 had significantly higher cfDNA derived from B cells, T cells and granulocytes and lower cfDNA from natural killer cells. Moreover, cfDNA derived from alveolar epithelial cells had the optimal performance to differentiate COVID-19 with different severities, lung injury levels, SOFA scores and in-hospital deaths, with the area under the receiver operating characteristic curve of 0.958, 0.941, 0.919 and 0.955, respectively.

CONCLUSION

Severe COVID-19 has a distinct cfDNA methylation signature compared with non-hospitalized COVID-19 and healthy controls. Cell type-specific cfDNA methylation signature enables the tracing of COVID-19 related cell deaths in lung and immune cells at cell-type resolution, which is correlated with clinical severities and outcomes, and has extensive application prospects to evaluate tissue injuries in diseases with multi-organ dysfunction.

DNA methylation landscape reveals GNAS as a decitabine-responsive marker in patients with acute myeloid leukemia

BACKGROUND

The demethylation agent decitabine (DAC) is a pivotal non-intensive alternative treatment for acute myeloid leukemia (AML). However, patient responses to DAC are highly variable, and predictive biomarkers are warranted. Herein, the DNA methylation landscape of patients treated with a DAC-based combination regimen was compared with that of patients treated with standard chemotherapy to develop a molecular approach for predicting clinical response to DAC.

METHODS

Twenty-five non-M3 AML patients were enrolled and subjected to DNA methylation sequencing and profiling to identify differentially methylated regions (DMRs) and genes of interest. Moreover, the effects of a DAC-based regimen on apoptosis and gene expression were explored using Kasumi-1 and K562 cells.

RESULTS

Overall, we identified 541 DMRs that were specifically responsive to DAC, among which 172 DMRs showed hypomethylation patterns upon treatment and were aligned with the promoter regions of 182 genes. In particular, GNAS was identified as a critical DAC-responsive gene, with in vitro GNAS downregulation leading to reduced cell apoptosis induced by DAC and cytarabine combo treatment.

CONCLUSIONS

We found that GNAS is a DAC-sensitive gene in AML and may serve as a prognostic biomarker to assess the responsiveness of patients with AML to DAC-based therapy.

Genome-wide methylation profiling of Peripheral T-cell lymphomas identifies TRIP13 as a critical driver of tumor proliferation and survival

Cytosine methylation of genomic DNA contributes to the regulation of gene expression and is involved in normal development including hematopoiesis in mammals. It is catalyzed by the family of DNA methyltransferases (DNMTs) that include DNMT1, DNMT3A, and DNMT3B. Peripheral T-cell lymphomas (PTCLs) represent a diverse group of aggressive mature T-cell malignancies accounting for approximately 10-15% of non-Hodgkin lymphoma cases in the US. PTCLs exhibit a broad spectrum of clinical, histological, and immunophenotypic features with poor prognosis and inadequately understood molecular pathobiology. To better understand the molecular landscape and identify candidate genes involved in disease maintenance, we used high-resolution Whole Genome Bisulfite Sequencing (WGBS) and RNA-seq to profile DNA methylation and gene expression of PTCLs and normal T-cells. We found that the methylation patterns in PTCLs are deregulated and heterogeneous but share 767 hypo- and 567 hypermethylated differentially methylated regions (DMRs) along with 231 genes up- and 91 genes downregulated in all samples suggesting a potential association with tumor development. We further identified 39 hypomethylated promoters associated with increased gene expression in the majority of PTCLs. This putative oncogenic signature included the TRIP13 (thyroid hormone receptor interactor 13) gene whose both genetic and pharmacologic inactivation, inhibited cellular growth of PTCL cell lines by inducing G2-M arrest accompanied by apoptosis suggesting that such an approach might be beneficial in human lymphoma treatment. Altogether we show that human PTCLs are characterized by a large number of recurrent methylation alterations, and demonstrated that TRIP13 is critical for PTCL maintenance in vitro.

Epigenetic variation in early and late flowering plants of the rubber-producing Russian dandelion Taraxacum koksaghyz provides insights into the regulation of flowering time

The Russian dandelion (Taraxacum koksaghyz) grows in temperate zones and produces large amounts of poly(cis-1,4-isoprene) in its roots, making it an attractive alternative source of natural rubber. Most T. koksaghyz plants require vernalization to trigger flower development, whereas early flowering varieties that have lost their vernalization dependence are more suitable for breeding and domestication. To provide insight into the regulation of flowering time in T. koksaghyz, we induced epigenetic variation by in vitro cultivation and applied epigenomic and transcriptomic analysis to the resulting early flowering plants and late flowering controls, allowing us to identify differences in methylation patterns and gene expression that correlated with flowering. This led to the identification of candidate genes homologous to vernalization and photoperiodism response genes in other plants, as well as epigenetic modifications that may contribute to the control of flower development. Some of the candidate genes were homologous to known floral regulators, including those that directly or indirectly regulate the major flowering control gene FT. Our atlas of genes can be used as a starting point to investigate mechanisms that control flowering time in T. koksaghyz in greater detail and to develop new breeding varieties that are more suited to domestication.

Source cell-type epigenetic memory persists in induced pluripotent cells but is lost in subsequently derived germline cells

Introduction: Retention of source cell-type epigenetic memory may mitigate the potential for induced pluripotent stem cells (iPSCs) to fully achieve transitions in cell fate in vitro. While this may not preclude the use of iPSC-derived somatic cell types for therapeutic applications, it becomes a major concern impacting the potential use of iPSC-derived germline cell types for reproductive applications. The transition from a source somatic cell type to iPSCs and then on to germ-cell like cells (GCLCs) recapitulates two major epigenetic reprogramming events that normally occur during development in vivo-embryonic reprogramming in the epiblast and germline reprogramming in primordial germ cells (PGCs). We examined the extent of epigenetic and transcriptomic memory persisting first during the transition from differentiated source cell types to iPSCs, and then during the transition from iPSCs to PGC-like cells (PGCLCs). Methods: We derived iPSCs from four differentiated mouse cell types including two somatic and two germ cell types and tested the extent to which each resulting iPSC line resembled a) a validated ES cell reference line, and b) their respective source cell types, on the basis of genome-wide gene expression and DNA methylation patterns. We then induced each iPSC line to form PGCLCs, and assessed epigenomic and transcriptomic memory in each compared to endogenous PGCs/M-prospermatogonia. Results: In each iPSC line, we found residual gene expression and epigenetic programming patterns characteristic of the corresponding source differentiated cell type from which each was derived. However, upon deriving PGCLCs, we found very little evidence of lingering epigenetic or transcriptomic memory of the original source cell type. Discussion: This result indicates that derivation of iPSCs and then GCLCs from differentiated source cell types in vitro recapitulates the two-phase epigenetic reprogramming that normally occurs in vivo, and that, to a significant extent, germline cell types derived in vitro from pluripotent cells accurately recapitulate epigenetic programming and gene expression patterns corresponding to equivalent endogenous germ cell types, suggesting that they have the potential to form the basis of in vitro gametogenesis as a useful therapeutic strategy for treatment of infertility.

Effects of maternal methyl donor intake during pregnancy on ileum methylation and function in an intrauterine growth restriction pig model

BACKGROUND

Intrauterine growth retardation (IUGR) affects intestinal growth, morphology, and function, which leads to poor growth performance and high mortality. The present study explored whether maternal dietary methyl donor (MET) supplementation alleviates IUGR and enhances offspring's growth performance by improving intestinal growth, function, and DNA methylation of the ileum in a porcine IUGR model.

METHODS

Forty multiparous sows were allocated to the control or MET diet groups from mating until delivery. After farrowing, 8 pairs of IUGR and normal birth weight piglets from 8 litters were selected for sampling before suckling colostrum.

RESULTS

The results showed that maternal MET supplementation tended to decrease the IUGR incidence and increased the average weaning weight of piglets. Moreover, maternal MET supplementation significantly reduced the plasma concentrations of isoleucine, cysteine, urea, and total amino acids in sows and newborn piglets. It also increased lactase and sucrase activity in the jejunum of newborn piglets. MET addition resulted in lower ileal methionine synthase activity and increased betaine homocysteine S-methyltransferase activity in the ileum of newborn piglets. DNA methylation analysis of the ileum showed that MET supplementation increased the methylation level of DNA CpG sites in the ileum of newborn piglets. Down-regulated differentially methylated genes were enriched in folic acid binding, insulin receptor signaling pathway, and endothelial cell proliferation. In contrast, up-regulated methylated genes were enriched in growth hormone receptor signaling pathway and nitric oxide biosynthetic process.

CONCLUSIONS

Maternal MET supplementation can reduce the incidence of IUGR and increase the weaning litter weight of piglets, which may be associated with better intestinal function and methylation status.

Rapid and accurate remethylation of DNA in Dnmt3a-deficient hematopoietic cells with restoration of DNMT3A activity

Here, we characterize the DNA methylation phenotypes of bone marrow cells from mice with hematopoietic deficiency of Dnmt3a or Dnmt3b (or both enzymes) or expressing the dominant-negative Dnmt3aR878H mutation [R882H in humans; the most common DNMT3A mutation found in acute myeloid leukemia (AML)]. Using these cells as substrates, we defined DNA remethylation after overexpressing wild-type (WT) DNMT3A1, DNMT3B1, DNMT3B3 (an inactive splice isoform of DNMT3B), or DNMT3L (a catalytically inactive "chaperone" for DNMT3A and DNMT3B in early embryogenesis). Overexpression of DNMT3A for 2 weeks reverses the hypomethylation phenotype of Dnmt3a-deficient cells or cells expressing the R878H mutation. Overexpression of DNMT3L (which is minimally expressed in AML cells) also corrects the hypomethylation phenotype of Dnmt3aR878H/+ marrow, probably by augmenting the activity of WT DNMT3A encoded by the residual WT allele. DNMT3L reactivation may represent a previously unidentified approach for restoring DNMT3A activity in hematopoietic cells with reduced DNMT3A function.

Epigenetically upregulated NSUN2 confers ferroptosis resistance in endometrial cancer via m5C modification of SLC7A11 mRNA

Endometrial cancer (EC) is a prevalent gynecological malignancy worldwide, and 5-methylcytosine (m5C) modification of mRNA is a crucial epigenetic modification associated with the development and occurrence of several cancers. However, the precise function of m5C modification in EC remains elusive. This study aimed to investigate the expression and clinical significance of the primary m5C modification writer, NSUN2, in EC. Our findings indicated that NSUN2 exhibited a substantial up-regulation in EC as a result of an epigenetic augmentation in H3K4me3 levels within the promoter region, which was triggered by the down-regulation of KDM5A. Moreover, gain- and loss-of-function experiments revealed the role of NSUN2 in enhancing m5C modification of mRNA, thereby promoting EC cell proliferation. RNA bisulfite sequencing and transcriptomic sequencing were employed to elucidate the involvement of NSUN2 in the regulation of ferroptosis. Subsequent in vitro experiments confirmed that the knockdown of NSUN2 significantly up-regulated the levels of lipid peroxides and lipid ROS in EC cells, thereby augmenting the susceptibility of EC to ferroptosis. Mechanistically, NSUN2 stimulated the m5C modification of SLC7A11 mRNA, and the m5C reader YBX1 exhibited direct recognition and binding to the m5C sites on SLC7A11 mRNA via its internal cold shock domain (CSD), leading to an increase in SLC7A11 mRNA stability and elevated levels of SLC7A11. Additionally, rescue experiments showed that NSUN2 functioned as a suppressor of ferroptosis, which was dependent on SLC7A11. Overall, targeting the NSUN2/SLC7A11 axis inhibited tumor growth by increasing lipid peroxidation and ferroptosis of EC cells both in vitro and in vivo. Therefore, our study provides new insight into the role of NSUN2, suggesting that NSUN2 may serve as a prognostic biomarker and therapeutic target in patients with EC.

DNA methylation in the wild: epigenetic transgenerational inheritance can mediate adaptation in clones of wild strawberry (Fragaria vesca)

Due to the accelerating climate change, it is crucial to understand how plants adapt to rapid environmental changes. Such adaptation may be mediated by epigenetic mechanisms like DNA methylation, which could heritably alter phenotypes without changing the DNA sequence, especially across clonal generations. However, we are still missing robust evidence of the adaptive potential of DNA methylation in wild clonal populations. Here, we studied genetic, epigenetic and transcriptomic variation of Fragaria vesca, a predominantly clonally reproducing herb. We examined samples from 21 natural populations across three climatically distinct geographic regions, as well as clones of the same individuals grown in a common garden. We found that epigenetic variation was partly associated with climate of origin, particularly in non-CG contexts. Importantly, a large proportion of this variation was heritable across clonal generations. Additionally, a subset of these epigenetic changes affected the expression of genes mainly involved in plant growth and responses to pathogen and abiotic stress. These findings highlight the potential influence of epigenetic changes on phenotypic traits. Our findings indicate that variation in DNA methylation, which can be environmentally inducible and heritable, may enable clonal plant populations to adjust to their environmental conditions even in the absence of genetic adaptation.

Exertional heat stroke causes long-term skeletal muscle epigenetic reprogramming, altered gene expression, and impaired satellite cell function in mice

The effect of exertional heat stroke (EHS) exposure on skeletal muscles is incompletely understood. Muscle weakness is an early symptom of EHS but is not considered a major target of multiorgan injury. Previously, in a preclinical mouse model of EHS, we observed the vulnerability of limb muscles to a second EHS exposure, suggesting hidden processes contributing to declines in muscle resilience. Here, we evaluated the possible molecular origins of EHS-induced declines in muscle resilience. Female C57BL/6 mice [total n = 56; 28/condition, i.e., EHS and exercise control (EXC)] underwent forced wheel running at 37.5°C/40% relative humidity until symptom limitation (unconsciousness). EXC mice exercised identically at room temperature (22-23°C). After 1 mo of recovery, the following were assessed: 1) specific force and caffeine-induced contracture in soleus (SOL) and extensor digitorum longus (EDL) muscles; 2) transcriptome and DNA methylome responses in gastrocnemius (GAST); and 3) primary satellite cell function (proliferation and differentiation). There were no differences in specific force in either SOL or EDL from EXC. Only EHS solei exhibited lower caffeine sensitivity. EHS GAST exhibited higher RNA expression of genes encoding structural proteins of slow fibers, heat shock proteins, and myogenesis. A total of ∼2,500 differentially methylated regions of DNA that could potentially affect many cell functions were identified. Primary satellite cells exhibited suppressed proliferation rates but normal differentiation responses. Results demonstrate long-term changes in skeletal muscles 1 mo after EHS that could contribute to declines in muscle resilience. Skeletal muscle may join other, more recognized tissues considered vulnerable to long-term effects of EHS.NEW & NOTEWORTHY Exertional heat stroke (EHS) in mice induces long-term molecular and functional changes in limb muscle that could reflect a loss of "resilience" to further stress. The phenotype was characterized by altered caffeine sensitivity and suppressed satellite cell proliferative potential. This was accompanied by changes in gene expression and DNA methylation consistent with ongoing muscle remodeling and stress adaptation. We propose that EHS may induce a prolonged vulnerability of skeletal muscle to further stress or injury.