Tissue type is a major modifier of the 5-hydroxymethylcytosine content of human genes

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

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

Gene body DNA hydroxymethylation restricts the magnitude of transcriptional changes during aging

DNA hydroxymethylation (5hmC), the most abundant oxidative derivative of DNA methylation, is typically enriched at enhancers and gene bodies of transcriptionally active and tissue-specific genes. Although aberrant genomic 5hmC has been implicated in age-related diseases, its functional role in aging remains unknown. Here, using mouse liver and cerebellum as model organs, we show that 5hmC accumulates in gene bodies associated with tissue-specific function and restricts the magnitude of gene expression changes with age. Mechanistically, 5hmC decreases the binding of splicing associated factors and correlates with age-related alternative splicing events. We found that various age-related contexts, such as prolonged quiescence and senescence, drive the accumulation of 5hmC with age. We provide evidence that this age-related transcriptionally restrictive function is conserved in mouse and human tissues. Our findings reveal that 5hmC regulates tissue-specific function and may play a role in longevity.

Genome-Wide Methylation Analysis Reveals a KCNK3-Prominent Causal Cascade on Hypertension

BACKGROUND

Despite advances in understanding hypertension's genetic structure, how noncoding genetic variants influence it remains unclear. Studying their interaction with DNA methylation is crucial to deciphering this complex disease's genetic mechanisms.

METHODS

We investigated the genetic and epigenetic interplay in hypertension using whole-genome bisulfite sequencing. Methylation profiling in 918 males revealed allele-specific methylation and methylation quantitative trait loci. We engineered rs1275988T/C mutant mice using CRISPR (clustered regularly interspaced short palindromic repeats)/Cas9 (CRISPR-associated protein 9), bred them for homozygosity, and subjected them to a high-salt diet. Telemetry captured their cardiovascular metrics. Protein-DNA interactions were elucidated using DNA pull-downs, mass spectrometry, and Western blots. A wire myograph assessed vascular function, and analysis of the Kcnk3 gene methylation highlighted the mutation's role in hypertension.

RESULTS

We discovered that DNA methylation-associated genetic effects, especially in non-cytosine-phosphate-guanine (non-CpG) island and noncoding distal regulatory regions, significantly contribute to hypertension predisposition. We identified distinct methylation quantitative trait locus patterns in the hypertensive population and observed that the onset of hypertension is influenced by the transmission of genetic effects through the demethylation process. By evidence-driven prioritization and in vivo experiments, we unearthed rs1275988 in a cell type-specific enhancer as a notable hypertension causal variant, intensifying hypertension through the modulation of local DNA methylation and consequential alterations in Kcnk3 gene expression and vascular remodeling. When exposed to a high-salt diet, mice with the rs1275988C/C genotype exhibited exacerbated hypertension and significant vascular remodeling, underscored by increased aortic wall thickness. The C allele of rs1275988 was associated with elevated DNA methylation levels, driving down the expression of the Kcnk3 gene by attenuating Nr2f2 (nuclear receptor subfamily 2 group F member 2) binding at the enhancer locus.

CONCLUSIONS

Our research reveals new insights into the complex interplay between genetic variations and DNA methylation in hypertension. We underscore hypomethylation's potential in hypertension onset and identify rs1275988 as a causal variant in vascular remodeling. This work advances our understanding of hypertension's molecular mechanisms and encourages personalized health care strategies.

Parkinson's disease-associated shifts between DNA methylation and DNA hydroxymethylation in human brain in PD-related genes, including PARK19 (DNAJC6) and PTPRN2 (IA-2β)

Background

The majority of Parkinson's disease (PD) cases are due to a complex interaction between aging, genetics, and environmental factors; epigenetic mechanisms are thought to act as important mediators of these risk factors. While multiple studies to date have explored the role of DNA modifications in PD, few focus on 5-hydroxymethylcytosine (5hmC). Because 5hmC occurs at its highest levels in the brain and is thought to be particularly important in the central nervous system, particularly in the response to neurotoxicants, it is important to explore the potential role of 5hmC in PD. This study expands on our previously published epigenome-wide association study (EWAS) performed on DNA isolated from neuron-enriched nuclei from human postmortem parietal cortex from the Banner Sun Health Research Institute Brain Bank. The study aimed to identify paired changes in 5hmC and 5mC in PD in enriched neuronal nuclei isolated from PD post-mortem parietal cortex and age- and sex-matched controls. We performed oxidative bisulfite (oxBS) conversion and paired it with our previously published bisulfite (BS)-based EWAS on the same samples to identify cytosines with significant shifts between these two related epigenetic marks. Interaction differentially modified cytosines (iDMCs) were identified using our recently published mixed-effects model for co-analyzing βmC and βhmC data.

Results

We identified 1,030 iDMCs with paired changes in 5mC and 5hmC (FDR < 0.05) that map to 695 genes, including PARK19 (DNAJC6), a familial PD gene, and PTPRN2 (IA-2), which has been previously implicated in PD in both epigenetic and mechanistic studies. The majority of iDMC-containing genes have not previously been implicated in PD and were not identified in our previous BS-based EWAS.

Conclusions

These data potentially link epigenetic regulation of the PARK19 and PTPRN2 loci in the pathogenesis of idiopathic PD. In addition, iDMC-containing genes have known functions in synaptic formation and function, cell cycle and senescence, neuroinflammation, and epigenetic regulation. These data suggest that there are significant shifts between 5mC and 5hmC associated with PD in genes relevant to PD pathogenesis that are not captured by analyzing BS-based data alone or by analyzing each mark as a distinct dataset.

Gene body DNA hydroxymethylation restricts the magnitude of transcriptional changes during aging

DNA hydroxymethylation (5hmC), the most abundant oxidative derivative of DNA methylation, is typically enriched at enhancers and gene bodies of transcriptionally active and tissue-specific genes. Although aberrant genomic 5hmC has been implicated in age-related diseases, its functional role in aging remains unknown. Here, using mouse liver and cerebellum as model organs, we show that 5hmC accumulates in gene bodies associated with tissue-specific function and restricts the magnitude of gene expression changes with age. Mechanistically, 5hmC decreases the binding of splicing associated factors and correlates with age-related alternative splicing events. We found that various age-related contexts, such as prolonged quiescence and senescence, drive the accumulation of 5hmC with age. We provide evidence that this age-related transcriptionally restrictive function is conserved in mouse and human tissues. Our findings reveal that 5hmC regulates tissue-specific function and may play a role in longevity.

iPSC-Derived Astrocytes and Neurons Replicate Brain Gene Expression, Epigenetic, Cell Morphology and Connectivity Alterations Found in Autism

Excessive inflammatory reactions and oxidative stress are well-recognized molecular findings in autism and these processes can affect or be affected by the epigenetic landscape. Nonetheless, adequate therapeutics are unavailable, as patient-specific brain molecular markers for individualized therapies remain challenging.

METHODS

We used iPSC-derived neurons and astrocytes of patients with autism vs. controls (5/group) to examine whether they replicate the postmortem brain expression/epigenetic alterations of autism. Additionally, DNA methylation of 10 postmortem brain samples (5/group) was analyzed for genes affected in PSC-derived cells.

RESULTS

We found hyperexpression of TGFB1, TGFB2, IL6 and IFI16 and decreased expression of HAP1, SIRT1, NURR1, RELN, GPX1, EN2, SLC1A2 and SLC1A3 in the astrocytes of patients with autism, along with DNA hypomethylation of TGFB2, IL6, TNFA and EN2 gene promoters and a decrease in HAP1 promoter 5-hydroxymethylation in the astrocytes of patients with autism. In neurons, HAP1 and IL6 expression trended alike. While HAP1 promoter was hypermethylated in neurons, IFI16 and SLC1A3 promoters were hypomethylated and TGFB2 exhibited increased promoter 5-hydroxymethlation. We also found a reduction in neuronal arborization, spine size, growth rate, and migration, but increased astrocyte size and a reduced growth rate in autism. In postmortem brain samples, we found DNA hypomethylation of TGFB2 and IFI16 promoter regions, but DNA hypermethylation of HAP1 and SLC1A2 promoters in autism.

CONCLUSION

Autism-associated expression/epigenetic alterations in iPSC-derived cells replicated those reported in the literature, making them appropriate surrogates to study disease pathogenesis or patient-specific therapeutics.

Adora2A downregulation promotes caffeine neuroprotective effect against LPS-induced neuroinflammation in the hippocampus

Caffeine has been extensively studied in the context of CNS pathologies as many researchers have shown that consuming it reduces pro-inflammatory biomarkers, potentially delaying the progression of neurodegenerative pathologies. Several lines of evidence suggest that adenosine receptors, especially A1 and A2A receptors, are the main targets of its neuroprotective action. We found that caffeine pretreatment 15 min before LPS administration reduced the expression of Il1b in the hippocampus and striatum. The harmful modulation of caffeine-induced inflammatory response involved the downregulation of the expression of A2A receptors, especially in the hippocampus. Caffeine treatment alone promoted the downregulation of the adenosinergic receptor Adora2A; however, this promotion effect was reversed by LPS. Although administering caffeine increased the expression of the enzymes DNA methyltransferases 1 and 3A and decreased the expression of the demethylase enzyme Tet1, this effect was reversed by LPS in the hippocampus of mice that were administered Caffeine + LPS, relative to the basal condition; no significant differences were observed in the methylation status of the promoter regions of adenosine receptors. Finally, the bioinformatics analysis of the expanded network demonstrated the following results: the Adora2B gene connects the extended networks of the adenosine receptors Adora1 and Adora2A; the Mapk3 and Esr1 genes connect the extended Adora1 network; the Mapk4 and Arrb2 genes connect the extended Adora2A network with the extended network of the proinflammatory cytokine Il1β. These results indicated that the anti-inflammatory effects of acute caffeine administration in the hippocampus may be mediated by a complex network of interdependencies between the Adora2B and Adora2A genes.

A Common OXTR Risk Variant Alters Regulation of Gene Expression by DNA Hydroxymethylation in Pregnant Human Myometrium

Postpartum hemorrhage, or excessive bleeding after birth, is a leading cause of maternal morbidity. A major cause of postpartum hemorrhage is uterine atony, tiring of the uterus which leads to ineffective contractions. Uterine contractions depend on oxytocin signaling in the myometrium, which in turn depends on expression of the oxytocin receptor (OXTR). Both genetic and epigenetic factors related to the oxytocin receptor are associated with risk of postpartum hemorrhage, but a mechanism relating these factors to oxytocin receptor activity in myometrium remains unclear. We report a genetic by epigenetic interaction whereby the relationship between DNA hydroxymethylation and OXTR gene expression depends on a common OXTR gene variant (rs53576). We also provide evidence that a similar genetic by epigenetic interaction using blood-derived DNA methylation is associated with relevant clinical outcomes: quantity of oxytocin administration and odds for postpartum hemorrhage. These results provide new avenues for predicting how women will respond to pharmacological agents in the prevention and treatment of postpartum hemorrhage.

Predicting gene expression state and prioritizing putative enhancers using 5hmC signal

BACKGROUND

Like its parent base 5-methylcytosine (5mC), 5-hydroxymethylcytosine (5hmC) is a direct epigenetic modification of cytosines in the context of CpG dinucleotides. 5hmC is the most abundant oxidized form of 5mC, generated through the action of TET dioxygenases at gene bodies of actively-transcribed genes and at active or lineage-specific enhancers. Although such enrichments are reported for 5hmC, to date, predictive models of gene expression state or putative regulatory regions for genes using 5hmC have not been developed.

RESULTS

Here, by using only 5hmC enrichment in genic regions and their vicinity, we develop neural network models that predict gene expression state across 49 cell types. We show that our deep neural network models distinguish high vs low expression state utilizing only 5hmC levels and these predictive models generalize to unseen cell types. Further, in order to leverage 5hmC signal in distal enhancers for expression prediction, we employ an Activity-by-Contact model and also develop a graph convolutional neural network model with both utilizing Hi-C data and 5hmC enrichment to prioritize enhancer-promoter links. These approaches identify known and novel putative enhancers for key genes in multiple immune cell subsets.

CONCLUSIONS

Our work highlights the importance of 5hmC in gene regulation through proximal and distal mechanisms and provides a framework to link it to genome function. With the recent advances in 6-letter DNA sequencing by short and long-read techniques, profiling of 5mC and 5hmC may be done routinely in the near future, hence, providing a broad range of applications for the methods developed here.

From powerhouse to regulator: The role of mitoepigenetics in mitochondrion-related cellular functions and human diseases

Beyond their crucial role in energy production, mitochondria harbor a distinct genome subject to epigenetic regulation akin to that of nuclear DNA. This paper delves into the nascent but rapidly evolving fields of mitoepigenetics and mitoepigenomics, exploring the sophisticated regulatory mechanisms governing mitochondrial DNA (mtDNA). These mechanisms encompass mtDNA methylation, the influence of non-coding RNAs (ncRNAs), and post-translational modifications of mitochondrial proteins. Together, these epigenetic modifications meticulously coordinate mitochondrial gene transcription, replication, and metabolism, thereby calibrating mitochondrial function in response to the dynamic interplay of intracellular needs and environmental stimuli. Notably, the dysregulation of mitoepigenetic pathways is increasingly implicated in mitochondrial dysfunction and a spectrum of human pathologies, including neurodegenerative diseases, cancer, metabolic disorders, and cardiovascular conditions. This comprehensive review synthesizes the current state of knowledge, emphasizing recent breakthroughs and innovations in the field. It discusses the potential of high-resolution mitochondrial epigenome mapping, the diagnostic and prognostic utility of blood or tissue mtDNA epigenetic markers, and the promising horizon of mitochondrial epigenetic drugs. Furthermore, it explores the transformative potential of mitoepigenetics and mitoepigenomics in precision medicine. Exploiting a theragnostic approach to maintaining mitochondrial allostasis, this paper underscores the pivotal role of mitochondrial epigenetics in charting new frontiers in medical science.

DNA Methylation in the Adaptive Response to Exercise

Emerging evidence published over the past decade has highlighted the role of DNA methylation in skeletal muscle function and health, including as an epigenetic transducer of the adaptive response to exercise. In this review, we aim to synthesize the latest findings in this field to highlight: (1) the shifting understanding of the genomic localization of altered DNA methylation in response to acute and chronic aerobic and resistance exercise in skeletal muscle (e.g., promoter, gene bodies, enhancers, intergenic regions, un-annotated regions, and genome-wide methylation); (2) how these global/regional methylation changes relate to transcriptional activity following exercise; and (3) the factors (e.g., individual demographic or genetic features, dietary, training history, exercise parameters, local epigenetic characteristics, circulating hormones) demonstrated to alter both the pattern of DNA methylation after exercise, and the relationship between DNA methylation and gene expression. Finally, we discuss the changes in non-CpG methylation and 5-hydroxymethylation after exercise, as well as the importance of emerging single-cell analyses to future studies-areas of increasing focus in the field of epigenetics. We anticipate that this review will help generate a framework for clinicians and researchers to begin developing and testing exercise interventions designed to generate targeted changes in DNA methylation as part of a personalized exercise regimen.

The Role of Different TET Proteins in Cytosine Demethylation Revealed by Mathematical Modeling

In living cells, some reactions can be conducted by more than one enzyme and sometimes it is difficult to establish which enzyme is responsible. Such is the case with proteins from the TET family, capable of converting 5-methyl-2'-deoxycytidine (5-mdC) in DNA to 5-(hydroxymethyl)-2'-deoxycytidine (5-hmdC) and further to 5-formyl-2'-deoxycytidine (5-fdC) and 5-carboxy-2'-deoxycytidine (5-cadC). The estimation of the efficiency of particular TETs in particular oxidative reactions and different cell types is important but experimentally difficult. Here, we propose an approach with mathematical modeling in which methylation and known deoxycytidine modification pathways are presented by 343 possible model versions with assumed different combinations of TET1, 2, and 3 activities in different pathways. Model parameters were calculated on the basis of 5-mdC, 5-hmdC, 5-fdC, 5-cadC, and 5-hmdU levels experimentally assessed in five human cultured cell lines and previously published. Selection of the model versions that give in simulations the best average fit to experimental data suggested that not all TET proteins participate in all modification reactions and that TET3 activity may be especially important in the reaction of 5-fdC removal.

Associations in cell type-specific hydroxymethylation and transcriptional alterations of pediatric central nervous system tumors

Although intratumoral heterogeneity has been established in pediatric central nervous system tumors, epigenomic alterations at the cell type level have largely remained unresolved. To identify cell type-specific alterations to cytosine modifications in pediatric central nervous system tumors, we utilize a multi-omic approach that integrated bulk DNA cytosine modification data (methylation and hydroxymethylation) with both bulk and single-cell RNA-sequencing data. We demonstrate a large reduction in the scope of significantly differentially modified cytosines in tumors when accounting for tumor cell type composition. In the progenitor-like cell types of tumors, we identify a preponderance differential Cytosine-phosphate-Guanine site hydroxymethylation rather than methylation. Genes with differential hydroxymethylation, like histone deacetylase 4 and insulin-like growth factor 1 receptor, are associated with cell type-specific changes in gene expression in tumors. Our results highlight the importance of epigenomic alterations in the progenitor-like cell types and its role in cell type-specific transcriptional regulation in pediatric central nervous system tumors.

Preservation of 5-Hydroxymethylcytosine Levels in LRIG1 across Genomic DNA and Cell-Free DNA in Glioma Patients

Cell-free DNA (cfDNA) has recently emerged as a promising minimally invasive diagnostic biomarker for various cancers. In this study, our aim was to identify cfDNA biomarkers by investigating genes that displayed significant differences between glioma patients and their corresponding controls. To accomplish this, we utilized publicly available data from the Gene Expression Omnibus, focusing on 5-hydroxymethylcytosine (5hmC) profiles in both cfDNA and genomic DNA (gDNA) from glioma patients and healthy individuals. The intersection of gene lists derived from these comparative analyses unveiled LRIG1 and ZNF703 as the two genes with elevated 5hmC levels in both the cfDNA of glioma patients and gDNA of glioma tissue compared to their respective controls. The gene expression data revealed both genes were upregulated in glioma tissue compared to normal brain tissue. Integration of 5hmC data revealed a strong positive correlation in the glioma tissue group between 5hmC and the gene expression of the LRIG1 gene. Furthermore, exploration using the AmiCa web tool indicated that LRIG1 gene expression was elevated compared to 17 other cancers included in the database, emphasizing its potential as a distinctive biomarker across multiple cancer types.

Nitric oxide inhibits ten-eleven translocation DNA demethylases to regulate 5mC and 5hmC across the genome

DNA methylation at cytosine bases of eukaryotic DNA (5-methylcytosine, 5mC) is a heritable epigenetic mark that can regulate gene expression in health and disease. Enzymes that metabolize 5mC have been well-characterized, yet the discovery of endogenously produced signaling molecules that regulate DNA methyl-modifying machinery have not been described. Herein, we report that the free radical signaling molecule nitric oxide (NO) can directly inhibit the Fe(II)/2-OG-dependent DNA demethylases ten-eleven translocation (TET) and human AlkB homolog 2 (ALKBH2). Physiologic NO concentrations reversibly inhibited TET and ALKBH2 demethylase activity by binding to the mononuclear non-heme iron atom which formed a dinitrosyliron complex (DNIC) preventing cosubstrates (2-OG and O2) from binding. In cancer cells treated with exogenous NO, or cells endogenously synthesizing NO, there was a global increase in 5mC and 5-hydroxymethylcytosine (5hmC) in DNA, the substrates for TET, that could not be attributed to increased DNA methyltransferase activity. 5mC was also elevated in NO-producing cell-line-derived mouse xenograft and patient-derived xenograft tumors. Genome-wide DNA methylome analysis of cells chronically treated with NO (10 days) demonstrated enrichment of 5mC and 5hmC at gene-regulatory loci which correlated to changes in the expression of NO-regulated tumor-associated genes. Regulation of DNA methylation is distinctly different from canonical NO signaling and represents a novel epigenetic role for NO.

Single-Molecule DNA Methylation Reveals Unique Epigenetic Identity Profiles of T Helper Cells

Both identity and plasticity of CD4 T helper (Th) cells are regulated in part by epigenetic mechanisms. However, a method that reliably and readily profiles DNA base modifications is still needed to finely study Th cell differentiation. Cytosine methylation in CpG context (5mCpG) and cytosine hydroxymethylation (5hmCpG) are DNA modifications that identify stable cell phenotypes, but their potential to characterize intermediate cell transitions has not yet been evaluated. To assess transition states in Th cells, we developed a method to profile Th cell identity using Cas9-targeted single-molecule nanopore sequencing. Targeting as few as 10 selected genomic loci, we were able to distinguish major in vitro polarized murine T cell subtypes, as well as intermediate phenotypes, by their native DNA 5mCpG patterns. Moreover, by using off-target sequences, we were able to infer transcription factor activities relevant to each cell subtype. Detection of 5mCpG and 5hmCpG was validated on intestinal Th17 cells escaping transforming growth factor β control, using single-molecule adaptive sampling. A total of 21 differentially methylated regions mapping to the 10-gene panel were identified in pathogenic Th17 cells relative to their nonpathogenic counterpart. Hence, our data highlight the potential to exploit native DNA methylation profiling to study physiological and pathological transition states of Th cells.

5-Hydroxymethylcytosine: the many faces of the sixth base of mammalian DNA

Epigenetic phenomena play a central role in cell regulatory processes and are important factors for understanding complex human disease. One of the best understood epigenetic mechanisms is DNA methylation. In the mammalian genome, cytosines (C) in CpG dinucleotides were long known to undergo methylation at the 5-position of the pyrimidine ring (mC). Later it was found that mC can be oxidized to 5-hydroxymethylcytosine (hmC) or even further to 5-formylcytosine (fC) and to 5-carboxylcytosine (caC) by the action of 2-oxoglutarate-dependent dioxygenases of the TET family. These findings unveiled a long elusive mechanism of active DNA demethylation and bolstered a wave of studies in the area of epigenetic regulation in mammals. This review is dedicated to critical assessment of recent data on biochemical and chemical aspects of the formation and conversion of hmC in DNA, analytical techniques used for detection and mapping of this nucleobase in mammalian genomes as well as epigenetic roles of hmC in DNA replication, transcription, cell differentiation and human disease.

Whole-Genome Bisulfite Sequencing Protocol for the Analysis of Genome-Wide DNA Methylation and Hydroxymethylation Patterns at Single-Nucleotide Resolution

The analysis of genome-wide epigenomic alterations including DNA methylation and hydroxymethylation has become a subject of intensive research for many biological and clinical questions. DNA methylation analysis bears the particular promise to supplement or replace biochemical and imaging-based tests for the next generation of personalized medicine. Whole-genome bisulfite sequencing (WGBS) using next-generation sequencing technologies is currently considered the gold standard for a comprehensive and quantitative analysis of DNA methylation throughout the genome. However, bisulfite conversion does not allow distinguishing between cytosine methylation and hydroxymethylation requiring an additional chemical or enzymatic step to identify hydroxymethylated cytosines. Here, we provide a detailed protocol based on a commercial kit for the preparation of sequencing libraries for the comprehensive whole-genome analysis of DNA methylation and/or hydroxymethylation. The protocol is based on the construction of sequencing libraries from limited amounts of input DNA by ligation of methylated adaptors to the fragmented DNA prior to bisulfite conversion. For analyses requiring a quantitative distinction between 5-methylcytosine and 5-hydroxymethylcytosines levels, an oxidation step is included in the same workflow to perform oxidative bisulfite sequencing (OxBs-Seq). In this case, two sequencing libraries will be generated and sequenced: a classic methylome following bisulfite conversion and analyzing modified cytosines (not distinguishing between methylated and hydroxymethylated cytosines) and a methylome analyzing only methylated cytosines, respectively. Hydroxymethylation levels are deduced from the differences between the two reactions. We also provide a step-by-step description of the data analysis using publicly available bioinformatic tools. The described protocol has been successfully applied to different human and plant samples and yields robust and reproducible results.

MicroRNA biogenesis machinery activation and lncRNA and REST overexpression as neuroprotective responses to fight inflammation in the hippocampus

Brain Long non-coding RNA (lncRNA) and microRNAs (miRs) play essential roles in the regulation of several important biological processes, including neuronal activity, cognitive processes, neurogenesis, angiogenesis, and neuroinflammation. In this context, the transcriptional repressor, RE1 silencing transcription factor (Rest), acts regulating the expression of neuronal genes as well as of lncRNAs and multiple miRNAs in the central nervous system. Nevertheless, its role in neuroinflammation was less explored. Here, we demonstrate, using an in vivo model of neuroinflammation induced by i.p. injection of LPS (0.33 mg/kg), that neuroinflammation increases gene expression of pro-inflammatory cytokines concomitant with the native and truncated forms of Rest and of non-coding RNAs. Additionally, the increased expression of enzymes Drosha ribonuclease III) (Drosha), Exportin 5 (Xpo5) and Endoribonuclease dicer (Dicer), associated with high expression of neuroprotective miRs 22 and 132 are indicative that the activation of biogenesis of miRs in the hippocampal region is a Central Nervous System (CNS) protective mechanism for the deleterious effects of neuroinflammation. Our results indicate that positive regulation of Rest gene expression in the hippocampal region by neuroinflammation correlates directly with the expression of miRs 22 and 132 and inversely with miR 335. In parallel, the confirmation of the possible alignment between the lncRNAs with miR 335 by bioinformatics corroborates with the sponge effect of Hottip and Hotair hybridizing and inhibiting the pro-inflammatory action of miR 335. This suggests the existence of a possible correlation between the activation of miR biogenesis machinery with increased expression of the transcription factor Rest, contributing to neuroprotection.

Steadfast Toll Like Receptor 4 (TLR4) 5-Hydroxymethylcytosine Levels in Cell-Free DNA: A Promising Consistency Marker for Colorectal Cancer Patients

Cell-free DNA (cfDNA) from patient blood is emerging as a noninvasive diagnostic avenue for various cancers. We aimed to identify reliable biomarkers in cfDNA by investigating genes exhibiting significant differences between colorectal cancer and control samples. Our objective was to identify genes that showed a positive difference between cancer and control samples. To achieve this, we conducted an in silico analysis to identify genes that exhibit no significant variation in methylation between genomic DNA (gDNA) and cfDNA. We collected experimental data from publicly available repositories, which included 5-hydroxymethylcytosine (5hmC) profiles of gDNA and cfDNA samples from both cancer patients and healthy individuals. By comparing and overlapping these two groups, we identified 187 genes of interest, of which 53 genes had a positive difference among colon cancer patients and healthy individuals. Next, we performed an ANOVA test on these genes, resulting in the identification of 12 genes that showed statistically significant higher levels of 5hmC in cfDNA and gDNA from cancer patients compared to healthy individuals. Additionally, we compared the 5hmC status of these genes between cfDNA and gDNA from cancer patients. Interestingly, we found that the 5hmC of the toll like receptor 4 (TLR4) gene was not statistically different between cfDNA and gDNA from cancer patients, indicating consistency between cfDNA and gDNA. These findings have important implications, not only for experimental validation but also for the development of more sensitive and robust noninvasive methods to improve diagnostic, prognostic, and treatment options for colon cancer.

Developmental exposures to common environmental contaminants, DEHP and lead, alter adult brain and blood hydroxymethylation in mice

Introduction: The developing epigenome changes rapidly, potentially making it more sensitive to toxicant exposures. DNA modifications, including methylation and hydroxymethylation, are important parts of the epigenome that may be affected by environmental exposures. However, most studies do not differentiate between these two DNA modifications, possibly masking significant effects. Methods: To investigate the relationship between DNA hydroxymethylation and developmental exposure to common contaminants, a collaborative, NIEHS-sponsored consortium, TaRGET II, initiated longitudinal mouse studies of developmental exposure to human-relevant levels of the phthalate plasticizer di(2-ethylhexyl) phthalate (DEHP), and the metal lead (Pb). Exposures to 25 mg DEHP/kg of food (approximately 5 mg DEHP/kg body weight) or 32 ppm Pb-acetate in drinking water were administered to nulliparous adult female mice. Exposure began 2 weeks before breeding and continued throughout pregnancy and lactation, until offspring were 21 days old. At 5 months, perinatally exposed offspring blood and cortex tissue were collected, for a total of 25 male mice and 17 female mice (n = 5-7 per tissue and exposure). DNA was extracted and hydroxymethylation was measured using hydroxymethylated DNA immunoprecipitation sequencing (hMeDIP-seq). Differential peak and pathway analysis was conducted comparing across exposure groups, tissue types, and animal sex, using an FDR cutoff of 0.15. Results: DEHP-exposed females had two genomic regions with lower hydroxymethylation in blood and no differences in cortex hydroxymethylation. For DEHP-exposed males, ten regions in blood (six higher and four lower) and 246 regions (242 higher and four lower) and four pathways in cortex were identified. Pb-exposed females had no statistically significant differences in blood or cortex hydroxymethylation compared to controls. Pb-exposed males, however, had 385 regions (all higher) and six pathways altered in cortex, but no differential hydroxymethylation was identified in blood. Discussion: Overall, perinatal exposure to human-relevant levels of two common toxicants showed differences in adult DNA hydroxymethylation that was specific to sex, exposure type, and tissue, but male cortex was most susceptible to hydroxymethylation differences by exposure. Future assessments should focus on understanding if these findings indicate potential biomarkers of exposure or are related to functional long-term health effects.

Ten-Eleven Translocation 1 and 2 Enzymes Affect Human Skin Fibroblasts in an Age-Related Manner

Ten-eleven translocation (TET) enzymes catalyze the oxidation of 5-methylcytosine (5mC), first to 5-hydroxymethylcytosine (5hmC), then to 5-formylcytosine (5fC), and finally to 5-carboxycytosine (5caC). Evidence suggests that changes in TET expression may impact cell function and the phenotype of aging. Proliferation, apoptosis, markers of autophagy and double-strand DNA break repair, and the expression of Fibulin 5 were assessed by flow cytometry in TET1 and TET2-overexpressing fibroblasts isolated from sun-unexposed skin of young (23-35 years) and age-advanced (75-94 years) individuals. In cells derived from young individuals, TET1 overexpression resulted in the inhibition of proliferation and apoptosis by 37% (p = 0.03) and 24% (p = 0.05), respectively, while the overexpression of TET2 caused a decrease in proliferation by 46% (p = 0.01). Notably, in cells obtained from age-advanced individuals, TETs exhibited different effects. Specifically, TET1 inhibited proliferation and expression of autophagy marker Beclin 1 by 45% (p = 0.05) and 28% (p = 0.048), respectively, while increasing the level of γH2AX, a marker of double-strand DNA breaks necessary for initiating the repair process, by 19% (p = 0.04). TET2 inhibited proliferation by 64% (p = 0.053) and increased the level of γH2AX and Fibulin 5 by 46% (p = 0.007) and 29% (p = 0.04), respectively. These patterns of TET1 and TET2 effects suggest their involvement in regulating various fibroblast functions and that some of their biological actions depend on the donor's age.

Mediation effects of DNA methylation and hydroxymethylation on birth outcomes after prenatal per- and polyfluoroalkyl substances (PFAS) exposure in the Michigan mother-infant Pairs cohort

BACKGROUND

Per- and polyfluoroalkyl substances (PFAS) are chemicals that are resistant to degradation and ubiquitous in our environments. PFAS may impact the developing epigenome, but current human evidence is limited to assessments of total DNA methylation. We assessed associations between first trimester PFAS exposures with newborn DNA methylation, including 5-methylcytosine (5-mC) and 5-hydroxymethylcytosine (5-hmC). DNA methylation mediation of associations between PFAS and birth outcomes were explored in the Michigan Mother Infant Pairs cohort. Nine PFAS were measured in maternal first trimester blood. Seven were highly detected and included for analysis: PFHxS, PFOA, PFOS, PFNA, PFDA, PFUnDA, and MeFOSAA. Bisulfite-converted cord blood DNA (n = 141) and oxidative-bisulfite-converted cord blood (n = 70) were assayed on Illumina MethylationEPIC BeadChips to measure total DNA methylation (5-mC + 5-hmC) and 5-mC/5-hmC. Correcting for multiple comparisons, beta regressions were used to assess associations between levels of PFAS and total methylation, 5-mC, or 5-hmC. Nonlinear mediation analyses were used to assess the epigenetic meditation effect between PFAS and birth outcomes.

RESULTS

PFAS was significantly associated with total methylation (q < 0.05: PFHxS-12 sites; PFOS-19 sites; PFOA-2 sites; PFNA-3 sites; PFDA-4 sites). In 72 female infants and 69 male infants, there were sex-specific associations between five PFAS and DNA methylation. 5-mC and 5-hmC were each significantly associated with thousands of sites for PFHxS, PFOS, PFNA, PFDA, PFUnDA, and MeFOSAA (q < 0.05). Clusters of 5-mC and 5-hmC sites were significant mediators between PFNA and PFUnDA and decreased gestational age (q < 0.05).

CONCLUSIONS

This study demonstrates the mediation role of specific types of DNA methylation on the relationship between PFAS exposure and birth outcomes. These results suggest that 5-mC and 5-hmC may be more sensitive to the developmental impacts of PFAS than total DNA methylation.

A View on Uterine Leiomyoma Genesis through the Prism of Genetic, Epigenetic and Cellular Heterogeneity

Uterine leiomyomas (ULs), frequent benign tumours of the female reproductive tract, are associated with a range of symptoms and significant morbidity. Despite extensive research, there is no consensus on essential points of UL initiation and development. The main reason for this is a pronounced inter- and intratumoral heterogeneity resulting from diverse and complicated mechanisms underlying UL pathobiology. In this review, we comprehensively analyse risk and protective factors for UL development, UL cellular composition, hormonal and paracrine signalling, epigenetic regulation and genetic abnormalities. We conclude the need to carefully update the concept of UL genesis in light of the current data. Staying within the framework of the existing hypotheses, we introduce a possible timeline for UL development and the associated key events—from potential prerequisites to the beginning of UL formation and the onset of driver and passenger changes.

EBS-seq: enrichment-based method for accurate analysis of 5-hydroxymethylcytosine at single-base resolution

BACKGROUND

A growing body of research has emphasized 5-hydroxymethylcytosine (5hmC) as an important epigenetic mark. High-resolution methods to detect 5hmC require high sequencing depth and are therefore expensive. Many studies have used enrichment-based methods to detect 5hmC; however, conventional enrichment-based methods have limited resolution. To overcome these limitations, we developed EBS-seq, a cost-efficient method for 5hmC detection with single-base resolution that combines the advantages of high-resolution methods and enrichment-based methods.

RESULTS

EBS-seq uses selective labeling of 5hmC, deamination of cytosine and 5-methylcytosine, pull-down of labeled 5hmC, and C-to-T conversion during DNA amplification. Using this method, we profiled 5hmC in HEK293T cells and two colorectal cancer samples. Compared with conventional enrichment-based 5hmC detection, EBS-seq improved 5hmC signals by localizing them at single-base resolution. Furthermore, EBS-seq was able to determine 5hmC levels in CpG-dense regions where distortion of signals can occur, such as CpG islands and CpG shores. Comparing EBS-seq and conventional high-resolution 5hmC detection by ACE-seq, we showed that EBS-seq is more effective at finding 5hmC sites. Using EBS-seq, we found strong associations between gene expression and gene-body 5hmC content in both HEK293T cells and colorectal cancer samples.

CONCLUSIONS

EBS-seq is a reliable and cost-efficient method for 5hmC detection because it simultaneously enriches 5hmC-containing DNA fragments and localizes 5hmC signals at single-base resolution. This method is a promising choice for 5hmC detection in challenging clinical samples with low 5hmC levels, such as cancer tissues.

Benchmarking of Nanopore R10.4 and R9.4.1 flow cells in single-cell whole-genome amplification and whole-genome shotgun sequencing

Third-generation sequencing can be used in human cancer genomics and epigenomic research. Oxford Nanopore Technologies (ONT) recently released R10.4 flow cell, which claimed an improved read accuracy compared to R9.4.1 flow cell. To evaluate the benefits and defects of R10.4 flow cell for cancer cell profiling on MinION devices, we used the human non-small-cell lung-carcinoma cell line HCC78 to construct libraries for both single-cell whole-genome amplification (scWGA) and whole-genome shotgun sequencing. The R10.4 and R9.4.1 reads were benchmarked in terms of read accuracy, variant detection, modification calling, genome recovery rate and compared with the next generation sequencing (NGS) reads. The results highlighted that the R10.4 outperforms R9.4.1 reads, achieving a higher modal read accuracy of over 99.1%, superior variation detection, lower false-discovery rate (FDR) in methylation calling, and comparable genome recovery rate. To achieve high yields scWGA sequencing in the ONT platform as NGS, we recommended multiple displacement amplification with a modified T7 endonuclease Ⅰ cutting procedure as a promising method. In addition, we provided a possible solution to filter the likely false positive sites among the whole genome region with R10.4 by using scWGA sequencing result as a negative control. Our study is the first benchmark of whole genome single-cell sequencing using ONT R10.4 and R9.4.1 MinION flow cells by clarifying the capacity of genomic and epigenomic profiling within a single flow cell. A promising method for scWGA sequencing together with the methylation calling results can benefit researchers who work on cancer cell genomic and epigenomic profiling using third-generation sequencing.

Hydroxymethylation alterations in progenitor-like cell types of pediatric central nervous system tumors are associated with cell type-specific transcriptional changes

Although intratumoral heterogeneity has been established in pediatric central nervous system tumors, epigenomic alterations at the cell type level have largely remained unresolved. To identify cell type-specific alterations to cytosine modifications in pediatric central nervous system tumors we utilized a multi-omic approach that integrated bulk DNA cytosine modification data (methylation and hydroxymethylation) with both bulk and single-cell RNA-sequencing data. We demonstrate a large reduction in the scope of significantly differentially modified cytosines in tumors when accounting for tumor cell type composition. In the progenitor-like cell types of tumors, we identified a preponderance differential CpG hydroxymethylation rather than methylation. Genes with differential hydroxymethylation, like HDAC4 and IGF1R, were associated with cell type-specific changes in gene expression in tumors. Our results highlight the importance of epigenomic alterations in the progenitor-like cell types and its role in cell type-specific transcriptional regulation in pediatric CNS tumors.

Genome-Wide Mapping Implicates 5-Hydroxymethylcytosines in Diabetes Mellitus and Alzheimer's Disease

BACKGROUND

Diabetes mellitus (DM) is a recognized risk factor for dementia. Because DM is a potentially modifiable condition, greater understanding of the mechanisms linking DM to the clinical expression of Alzheimer's disease dementia may provide insights into much needed dementia therapeutics.

OBJECTIVE

In this feasibility study, we investigated DM as a dementia risk factor by examining genome-wide distributions of the epigenetic DNA modification 5-hydroxymethylcytosine (5hmC).

METHODS

We obtained biologic samples from the Rush Memory and Aging Project and used the highly sensitive 5hmC-Seal technique to perform genome-wide profiling of 5hmC in circulating cell-free DNA (cfDNA) from antemortem serum samples and in genomic DNA from postmortem prefrontal cortex brain tissue from 80 individuals across four groups: Alzheimer's disease neuropathologically defined (AD), DM clinically defined, AD with DM, and individuals with neither disease (controls).

RESULTS

Distinct 5hmC signatures and biological pathways were enriched in persons with both AD and DM versus AD alone, DM alone, or controls, including genes inhibited by EGFR signaling in oligodendroglia and those activated by constitutive RHOA. We also demonstrate the potential diagnostic value of 5hmC profiling in circulating cfDNA. Specifically, an 11-gene weighted model distinguished AD from non-AD/non-DM controls (AUC = 91.8%; 95% CI, 82.9-100.0%), while a 4-gene model distinguished DM-associated AD from AD alone (AUC = 87.9%; 95% CI, 77.5-98.3%).

CONCLUSION

We demonstrate in this small sample, the feasibility of detecting and characterizing 5hmC in DM-associated AD and of using 5hmC information contained in circulating cfDNA to detect AD in high-risk individuals, such as those with diabetes.

Preconception exposure to malathion and glucose homeostasis in rats: Effects on dams during pregnancy and post-term periods, and on their progeny

Understanding the individual and global impact of pesticides on human physiology and the different stages of life is still a challenge in environmental health. We analyzed here whether administration of the organophosphate insecticide malathion before pregnancy could affect glucose homeostasis during pregnancy and, in addition, generate possible later consequences in mothers and offspring. For this, adult Wistar rats were allocated into two groups and were treated daily (intragastric) with malathion (14 or 140 mg/kg, body mass (bm)) for 21-25 days. Corn oil was used as vehicle in the Control group. Subgroups were defined based on the absence (nulliparous) or presence (pregnant) of a copulatory plug. Pregnant rats were followed by an additional period of 2 months after the term (post-term), without continuing malathion treatment. Fetuses and adult offspring of males and females were also evaluated. We ran an additional experimental design with rats exposed to malathion before pregnancy at a dose of 0.1 mg/kg bm. Malathion exposure resulted in glucose intolerance in the mothers during pregnancy and post-term period, regardless of the exposure dose. This was accompanied by increased visceral adipose tissue mass, dyslipidemia, unchanged pancreatic β-cell mass, and varying insulin responses to glucose in vivo. The number of total newborns and birthweight was not affected by malathion exposure. Adult offspring from both sexes also became glucose-intolerant, regardless of the pesticide dose their dams were exposed to. This alteration could be associated with changes at the epigenomic level, as reduced hepatic mRNA content of DNA methylases and demethylases was found. We demonstrated that periconceptional exposure to malathion with doses aiming to mimic from work environment to indirect contamination predisposes progenitors and offspring rats to glucose intolerance. Thus, we conclude that subchronic exposure to malathion is a risk factor for gestational diabetes and prediabetes later in life.

The 5-Hydroxymethylcytosine Landscape of Prostate Cancer

Analysis of DNA methylation is a valuable tool to understand disease progression and is increasingly being used to create diagnostic and prognostic clinical biomarkers. While conversion of cytosine to 5-methylcytosine (5mC) commonly results in transcriptional repression, further conversion to 5-hydroxymethylcytosine (5hmC) is associated with transcriptional activation. Here we perform the first study integrating whole-genome 5hmC with DNA, 5mC, and transcriptome sequencing in clinical samples of benign, localized, and advanced prostate cancer. 5hmC is shown to mark activation of cancer drivers and downstream targets. Furthermore, 5hmC sequencing revealed profoundly altered cell states throughout the disease course, characterized by increased proliferation, oncogenic signaling, dedifferentiation, and lineage plasticity to neuroendocrine and gastrointestinal lineages. Finally, 5hmC sequencing of cell-free DNA from patients with metastatic disease proved useful as a prognostic biomarker able to identify an aggressive subtype of prostate cancer using the genes TOP2A and EZH2, previously only detectable by transcriptomic analysis of solid tumor biopsies. Overall, these findings reveal that 5hmC marks epigenomic activation in prostate cancer and identify hallmarks of prostate cancer progression with potential as biomarkers of aggressive disease.

SIGNIFICANCE

In prostate cancer, 5-hydroxymethylcytosine delineates oncogene activation and stage-specific cell states and can be analyzed in liquid biopsies to detect cancer phenotypes. See related article by Wu and Attard, p. 3880.

Placental DNA Methylation in pregnancies complicated by maternal diabetes and/or obesity: State of the Art and research gaps

SUMMARYMaternal diabetes and/or obesity in pregnancy are undoubtedly associated with later disease-risk in the offspring. The placenta, interposed between the mother and the fetus, is a potential mediator of this risk through epigenetic mechanisms, including DNA methylation. In recent years, multiple studies have identified differentially methylated CpG sites in the placental tissue DNA in pregnancies complicated by diabetes and obesity. We reviewed all published original research relevant to this topic and analyzed our findings with the focus of identifying overlaps, contradictions and gaps. Most studies focused on the association of gestational diabetes and/or hyperglycemia in pregnancy and DNA methylation in placental tissue at term. We identified overlaps in results related to specific candidate genes, but also observed a large research gap of pregnancies affected by type 1 diabetes. Other unanswered questions relate to analysis of specific placental cell types and the timing of DNA methylation change in response to diabetes and obesity during pregnancy. Maternal metabolism is altered already in the first trimester involving structural and functional changes in the placenta, but studies into its effects on placental DNA methylation during this period are lacking and urgently needed. Fetal sex is also an important determinant of pregnancy outcome, but only few studies have taken this into account. Collectively, we provide a reference work for researchers working in this large and evolving field. Based on the results of the literature review, we formulate suggestions for future focus of placental DNA methylation studies in pregnancies complicated by diabetes and obesity.

Tight correlation of 5-hydroxymethylcytosine expression with the scarring damage of discoid lupus erythematosus

OBJECTIVE

Cutaneous lupus erythematosus (CLE) is a heterogenous skin disease. The two most common subtypes are discoid LE (DLE) characterized by scarring skin damage and acute CLE (ACLE) presenting with transiently reversible skin lesions. It remains unknown what causes the difference of skin lesions. Studies have shown the existence of tissue-specific 5-Hydroxymethylcytosine (5 hmC)-modified regions in human tissues, which may affect the tissue-related diseases. Here, we aim to assess the expression of 5 hmc in DLE and ACLE lesions and explore the relationship of 5 hmc with scarring damage in DLE.

METHODS

84 CLE samples were included in the study. We evaluated the skin damage score and reviewed the histopathologic sections. Immunohistochemical staining was performed to detect the expression of 5 hmc in the appendage and periappendageal inflammatory cells. The 5 hmc expression in periappendageal lymphocytic cells was investigated by multi-spectrum immunohistochemistry staining.

RESULTS

Scarring/atrophy was the most significant damage in differentiating the DLE from ACLE. Perifollicular inflammatory infiltration was present in all patients with DLE scarring alopecia (DLESA). The 5 hmc expression in the appendage and periappendageal inflammatory cells was significantxly increased in DLESA than ACLE. Similar expression pattern was seen in the staining of IFN-alpha/beta Receptor (IFNAR). The expression of 5 hmc in the appendage was positively correlated with that in the periappendageal inflammatory cells. There was an increased 5 hmc expression in lymphocytes cluster around hair follicle consisting of CD4⁺ cells, CD8⁺ cells, and CD19⁺ cells in DLESA lesions.

CONCLUSION

These data demonstrate a close association of the expression pattern of 5 hmc with the histopathological characteristic distribution, and with the type I interferons (IFNs) signals in DLESA, supporting the importance of 5 hmc in the amplification of appendage damage and periappendageal inflammation, thereby offering a novel insight into the scarring damage of DLE and the heterogeneity of CLE skin lesions.

Distinct cytosine modification profiles define epithelial-to-mesenchymal cell-state transitions

Background: Epithelial-to-mesenchymal transition (EMT) is an early step in the invasion-metastasis cascade, involving progression through intermediate cell states. Due to challenges with isolating intermediate cell states, genome-wide cytosine modifications that define transition are not completely understood. Methods: The authors measured multiple DNA cytosine modification marks and chromatin accessibility across clonal populations residing in specific EMT states. Results: Clones exhibiting more intermediate EMT phenotypes demonstrated increased 5-hydroxymethylcytosine and decreased 5-methylcytosine. Open chromatin regions containing increased 5-hydroxymethylcytosine CpG loci were enriched in EMT transcription factor motifs and were associated with Rho GTPases. Conclusion: The results indicate the importance of both distinct and shared epigenetic profiles associated with EMT processes that may be targeted to prevent EMT progression.

Regulation of 5-Hydroxymethylcytosine by TET2 Contributes to Squamous Cell Carcinoma Tumorigenesis

DNA methylation is a key regulatory event controlling a variety of physiological processes and can have dramatic effects on gene transcription. Methylated cytosine (5-methylcytosine) can be oxidized by the TET family of enzymes to 5-hydroxymethylcytosine (5-hmC), a key intermediate in the demethylation cycle, and 5-hmC levels are reduced in malignancies such as acute myeloid leukemia and melanoma. We constructed a tissue microarray of human cutaneous squamous cell carcinoma tumors and found a global reduction in 5-hmC levels compared with that in the adjacent skin. Using a murine K14-CreER system, we have found that loss of Tet2 promotes carcinogen-induced squamous cell carcinoma and cooperates with loss of Tp53 to drive spontaneous squamous cell carcinoma tumors in epithelial tissues. Analysis of changes in 5-hmC and gene expression after loss of Tet2 in the epidermis revealed focal alterations in 5-hmC levels and an increase in hair follicle transient amplifying cell genes along with a reduction in epidermal differentiation genes. These results show a role for TET2 in epidermal lineage specification, consistent with reported roles for TET enzymes in controlling lineage commitment in hematopoietic stem cells and embryonic stem cells and establishing TET2 as a bone fide tumor suppressor in squamous cell carcinoma.

5-Hydroxymethylcytosine (5hmC) at or near cancer mutation hot spots as potential targets for early cancer detection

Objective

Universal noninvasive genomic screening to detect cancer and/or fetal DNA in plasma at all stages of development is highly warranted. Since 5-hydroxymethylcytosine (5hmC) emerged as an intermediate metabolite in active DNA demethylation, there have been increasing efforts to elucidate its function as a stable modification of the genome. In the current study, we demonstrate that discrete 5hmC sites within 80 bp hotspot regions exist in a greater proportion of cancer versus normal cells.

Result

5hmC was detected in 16 of 17 known hotspots having C to T or G to A mutations. The results show the presence of two characteristically distinct 5hmC groups: Tier 1 Group with 3 to eightfold more 5hmCs detected in tumor-cells than in normal-cell derived DNA (as observed in 6 of 11 CpG sites). Tier 2 group with equal allele frequency of 5hmC among normal and tumor-cell derived DNA at 5 CpG hotspot sites as well as 5 non-CpG hotspots. Thus, detection and quantification of the Tier 1 group of 5hmC sites or its prevalence at or near cancer mutation hot spots in cells may enable early detection, screening and potentially prediction of the likelihood of cancer occurrence or the severity of the cancer.

DNA Hydroxymethylation in Smoking-Associated Cancers

5-hydroxymethylcytosine (5-hmC) was first detected in mammalian DNA five decades ago. However, it did not take center stage in the field of epigenetics until 2009, when ten-eleven translocation 1 (TET1) was found to oxidize 5-methylcytosine to 5-hmC, thus offering a long-awaited mechanism for active DNA demethylation. Since then, a remarkable body of research has implicated DNA hydroxymethylation in pluripotency, differentiation, neural system development, aging, and pathogenesis of numerous diseases, especially cancer. Here, we focus on DNA hydroxymethylation in smoking-associated carcinogenesis to highlight the diagnostic, therapeutic, and prognostic potentials of this epigenetic mark. We describe the significance of 5-hmC in DNA demethylation, the importance of substrates and cofactors in TET-mediated DNA hydroxymethylation, the regulation of TETs and related genes (isocitrate dehydrogenases, fumarate hydratase, and succinate dehydrogenase), the cell-type dependency and genomic distribution of 5-hmC, and the functional role of 5-hmC in the epigenetic regulation of transcription. We showcase examples of studies on three major smoking-associated cancers, including lung, bladder, and colorectal cancers, to summarize the current state of knowledge, outstanding questions, and future direction in the field.

Epigenomic and enhancer dysregulation in uterine leiomyomas

BACKGROUND

Uterine leiomyomas, also known as uterine fibroids or myomas, are the most common benign gynecological tumors and are found in women of reproductive and postmenopausal age. There is an exceptionally high prevalence of this tumor in women by the age of 50 years. Black women are particularly affected, with an increased incidence, earlier age of onset, larger and faster growing fibroids and greater severity of symptoms as compared to White women. Although advances in identifying genetic and environmental factors to delineate these fibroids have already been made, only recently has the role of epigenomics in the pathogenesis of this disease been considered.

OBJECTIVE AND RATIONALE

Over recent years, studies have identified multiple epigenomic aberrations that may contribute to leiomyoma development and growth. This review will focus on the most recent discoveries in three categories of epigenomic changes found in uterine fibroids, namely aberrant DNA methylation, histone tail modifications and histone variant exchange, and their translation into altered target gene architecture and transcriptional outcome. The findings demonstrating how the altered 3D shape of the enhancer can regulate gene expression from millions of base pairs away will be discussed. Additionally, translational implications of these discoveries and potential roadblocks in leiomyoma treatment will be addressed.

SEARCH METHODS

A comprehensive PubMed search was performed to identify published articles containing keywords relevant to the focus of the review, such as: uterine leiomyoma, uterine fibroids, epigenetic alterations, epigenomics, stem cells, chromatin modifications, extracellular matrix [ECM] organization, DNA methylation, enhancer, histone post-translational modifications and dysregulated gene expression. Articles until September 2021 were explored and evaluated to identify relevant updates in the field. Most of the articles focused on in the discussion were published between 2015 and 2021, although some key discoveries made before 2015 were included for background information and foundational purposes. We apologize to the authors whose work was not included because of space restrictions or inadvertent omission.

OUTCOMES

Chemical alterations to the DNA structure and of nucleosomal histones, without changing the underlying DNA sequence, have now been implicated in the phenotypic manifestation of uterine leiomyomas. Genome-wide DNA methylation analysis has revealed subsets of either suppressed or overexpressed genes accompanied by aberrant promoter methylation. Furthermore, differential promoter access resulting from altered 3D chromatin structure and histone modifications plays a role in regulating transcription of key genes thought to be involved in leiomyoma etiology. The dysregulated genes function in tumor suppression, apoptosis, angiogenesis, ECM formation, a variety of cancer-related signaling pathways and stem cell differentiation. Aberrant DNA methylation or histone modification is also observed in altering enhancer architecture, which leads to changes in enhancer-promoter contact strength, producing novel explanations for the overexpression of high mobility group AT-hook 2 and gene dysregulation found in mediator complex subunit 12 mutant fibroids. While many molecular mechanisms and epigenomic features have been investigated, the basis for the racial disparity observed among those in the Black population remains unclear.

WIDER IMPLICATIONS

A comprehensive understanding of the exact pathogenesis of uterine leiomyoma is lacking and requires attention as it can provide clues for prevention and viable non-surgical treatment. These findings will widen our knowledge of the role epigenomics plays in the mechanisms related to uterine leiomyoma development and highlight novel approaches for the prevention and identification of epigenome targets for long-term non-invasive treatment options of this significantly common disease.

Calcium handling genes are regulated by promoter DNA methylation in colorectal cancer cells

Calcium signaling regulates various cellular processes, including proliferation and cell death. DNA methylation of gene promoters is an epigenetic modification that facilitates transcriptional suppression. Disruption of calcium homeostasis and DNA methylation in cancer are each linked to tumor development and progression. However, the possible connection between these two processes has not been thoroughly studied. Therefore, we measured the expression of six gene families involved in calcium regulation (ATP2A, ITPR, ORAI, RyR, STIM, and TRPC) in a colorectal cancer cell model, HCT116, with either genetic (Double Knock-out/DKO) or pharmacological (5-aza-2'-deoxycytidine/DAC) inhibition of DNA methyltransferases. Fourteen of the 20 examined calcium handling genes were expressed at higher levels in DKO cells as compared to HCT116. Expression of five genes was increased in HCT116 cells treated with DAC, three matching DKO. Due to a unique expression pattern of the three ATP2A genes in our model, encoding the Sarcoplasmic/Endoplasmic Reticulum Calcium ATPase (SERCA) pumps, we chose to evaluate the methylation status of these genes, protein expression, and potential associated physiological effects, using the SERCA inhibitor thapsigarin (TG). We observed an expected pattern of promoter methylation coinciding with reduced expression and vice versa. This differential mRNA expression was associated with altered SERCA3 protein expression and cytosolic calcium levels with TG exposure. As a result, DKO cells displayed less TG-induced cytotoxicity, as compared to HCT116 cells. Overall, it is likely that at least several calcium regulatory genes are transcriptionally regulated by DNA methylation, and this may play a role in tumorigenesis through altering apoptosis in cancer.

Dysregulated TET Family Genes and Aberrant 5mC Oxidation in Breast Cancer: Causes and Consequences

Simple Summary

Both genetic and epigenetic mechanisms contribute to the pathogenesis of breast cancer. Since Tahiliani et al. identified TET1 as the first methyl-cytosine dioxygenase in 2009, accumulating evidence has shown that aberrant 5mC oxidation and dysregulated TET family genes are associated with diseases, including breast cancer. In this review we provide an overview on the diagnosis and prognosis values of aberrant 5mC oxidation in breast cancer and emphasize the causes and consequences of such epigenetic alterations.

Abstract

DNA methylation (5-methylcytosine, 5mC) was once viewed as a stable epigenetic modification until Rao and colleagues identified Ten-eleven translocation 1 (TET1) as the first 5mC dioxygenase in 2009. TET family genes (including TET1, TET2, and TET3) encode proteins that can catalyze 5mC oxidation and consequently modulate DNA methylation, not only regulating embryonic development and cellular differentiation, but also playing critical roles in various physiological and pathophysiological processes. Soon after the discovery of TET family 5mC dioxygenases, aberrant 5mC oxidation and dysregulation of TET family genes have been reported in breast cancer as well as other malignancies. The impacts of aberrant 5mC oxidation and dysregulated TET family genes on the different aspects (so-called cancer hallmarks) of breast cancer have also been extensively investigated in the past decade. In this review, we summarize current understanding of the causes and consequences of aberrant 5mC oxidation in the pathogenesis of breast cancer. The challenges and future perspectives of this field are also discussed.

Oxytocin prevents cue-induced reinstatement of oxycodone seeking: Involvement of DNA methylation in the hippocampus

Oxycodone is one of the most commonly used analgesics in the clinic. However, long-term use can contribute to drug dependence. Accumulating evidence of changes in DNA methylation after opioid relapse has provided insight into mechanisms underlying drug-associated memory. The neuropeptide oxytocin is reported to be a potential treatment for addiction. The present study sought to identify changes in global and synaptic gene methylation after cue-induced reinstatement of oxycodone conditioned place preference (CPP) and the effect of oxytocin. We analyzed hippocampal mRNA of synaptic genes and also synaptic density in response to oxycodone CPP. We determined the mRNA levels of DNA methyltransferases (Dnmts) and ten-eleven translocations (Tets), observed global 5-methylcytosine (5-mC) and 5-hydroxymethylcytosine (5-hmC) levels, and measured DNA methylation status of four synaptic genes implicated in learning and memory (Arc, Dlg1, Dlg4, and Syn1). Both synaptic density and the transcription of 15 hippocampal synaptic genes significantly increased following cue-induced reinstatement of oxycodone CPP. Oxycodone relapse was also related to markedly decreased 5-mC levels and decreased transcription of Dnmt1, Dnmt3a, and Dnmt3b; in contrast, 5-hmC levels and the transcription of Tet1 and Tet3 were increased. Oxycodone exposure induced DNA hypomethylation at the exons of the Arc, Dlg1, Dlg4, and Syn1 genes. Intracerebroventricular (ICV) administration of oxytocin (2.5 μg/μl) specifically blocked oxycodone relapse, possibly by inhibition of Arc, Dlg1, Dlg4, and Syn1 hypomethylation in oxycodone-treated rats. Together, these data indicate the occurrence of epigenetic changes in the hippocampus following oxycodone relapse and the potential role of oxytocin in oxycodone addiction.

Cell-free DNA 5-hydroxymethylcytosine profiles of long non-coding RNA genes enable early detection and progression monitoring of human cancers

BACKGROUND

5-Hydroxymethylcytosine (5hmC) is a significant DNA epigenetic modification. However, the 5hmC modification alterations in genomic regions encoding long non-coding RNA (lncRNA) and their clinical significance remain poorly characterized.

RESULTS

A three-phase discovery-modeling-validation study was conducted to explore the potential of the plasma-derived 5hmC modification level in genomic regions encoding lncRNAs as a superior alternative biomarker for cancer diagnosis and surveillance. Genome-wide 5hmC profiles in the plasma circulating cell-free DNA of 1632 cancer and 1379 non-cancerous control samples from different cancer types and multiple centers were repurposed and characterized. A large number of altered 5hmC modifications were distributed at genomic regions encoding lncRNAs in cancerous compared with healthy subjects. Furthermore, most 5hmC-modified lncRNA genes were cancer-specific, with only a relatively small number of 5hmC-modified lncRNA genes shared by various cancer types. A 5hmC-LncRNA diagnostic score (5hLD-score) comprising 39 tissue-shared 5hmC-modified lncRNA gene markers was developed using elastic net regularization. The 5hLD-score was able to accurately distinguish tumors from healthy controls with an area under the curve (AUC) of 0.963 [95% confidence interval (CI) 0.940-0.985] and 0.912 (95% CI 0.837-0.987) in the training and internal validation cohorts, respectively. Results from three independent validations confirmed the robustness and stability of the 5hLD-score with an AUC of 0.851 (95% CI 0.786-0.916) in Zhang's non-small cell lung cancer cohort, AUC of 0.887 (95% CI 0.852-0.922) in Tian's esophageal cancer cohort, and AUC of 0.768 (95% CI 0.746-0.790) in Cai's hepatocellular carcinoma cohort. In addition, a significant association was identified between the 5hLD-score and the progression from hepatitis to liver cancer. Finally, lncRNA genes modified by tissue-specific 5hmC alteration were again found to be capable of identifying the origin and location of tumors.

CONCLUSION

The present study will contribute to the ongoing effort to understand the transcriptional programs of lncRNA genes, as well as facilitate the development of novel invasive genomic tools for early cancer detection and surveillance.

Ascorbate Inhibits Proliferation and Promotes Myeloid Differentiation in TP53-Mutant Leukemia

Loss-of-function mutations in the DNA demethylase TET2 are associated with the dysregulation of hematopoietic stem cell differentiation and arise in approximately 10% of de novo acute myeloid leukemia (AML). TET2 mutations coexist with other mutations in AML, including TP53 mutations, which can indicate a particularly poor prognosis. Ascorbate can function as an epigenetic therapeutic in pathological contexts involving heterozygous TET2 mutations by restoring TET2 activity. How this response is affected when myeloid leukemia cells harbor mutations in both TET2 and TP53 is unknown. Therefore, we examined the effects of ascorbate on the SKM-1 AML cell line that has mutated TET2 and TP53. Sustained treatment with ascorbate inhibited proliferation and promoted the differentiation of these cells. Furthermore, ascorbate treatment significantly increased 5-hydroxymethylcytosine, suggesting increased TET activity as the likely mechanism. We also investigated whether ascorbate affected the cytotoxicity of Prima-1^Met, a drug that reactivates some p53 mutants and is currently in clinical trials for AML. We found that the addition of ascorbate had a minimal effect on Prima-1^Met–induced cytotoxicity, with small increases or decreases in cytotoxicity being observed depending on the timing of treatment. Collectively, these data suggest that ascorbate could exert a beneficial anti-proliferative effect on AML cells harboring both TET2 and TP53 mutations whilst not interfering with targeted cytotoxic therapies such as Prima-1^Met.

Epigenetic Deregulation of the Histone Methyltransferase KMT5B Contributes to Malignant Transformation in Glioblastoma

Glioblastoma multiforme (GBM) is the most common and aggressive type of brain tumor in adulthood. Epigenetic mechanisms are known to play a key role in GBM although the involvement of histone methyltransferase KMT5B and its mark H4K20me2 has remained largely unexplored. The present study shows that DNA hypermethylation and loss of DNA hydroxymethylation is associated with KMT5B downregulation and genome-wide reduction of H4K20me2 levels in a set of human GBM samples and cell lines as compared with non-tumoral specimens. Ectopic overexpression of KMT5B induced tumor suppressor-like features in vitro and in a mouse tumor xenograft model, as well as changes in the expression of several glioblastoma-related genes. H4K20me2 enrichment was found immediately upstream of the promoter regions of a subset of deregulated genes, thus suggesting a possible role for KMT5B in GBM through the epigenetic modulation of key target cancer genes.

Ten-eleven translocation 1 mediated-DNA hydroxymethylation is required for myelination and remyelination in the mouse brain

Ten-eleven translocation (TET) proteins, the dioxygenase for DNA hydroxymethylation, are important players in nervous system development and diseases. However, their role in myelination and remyelination after injury remains elusive. Here, we identify a genome-wide and locus-specific DNA hydroxymethylation landscape shift during differentiation of oligodendrocyte-progenitor cells (OPC). Ablation of Tet1 results in stage-dependent defects in oligodendrocyte (OL) development and myelination in the mouse brain. The mice lacking Tet1 in the oligodendrocyte lineage develop behavioral deficiency. We also show that TET1 is required for remyelination in adulthood. Transcriptomic, genomic occupancy, and 5-hydroxymethylcytosine (5hmC) profiling reveal a critical TET1-regulated epigenetic program for oligodendrocyte differentiation that includes genes associated with myelination, cell division, and calcium transport. Tet1-deficient OPCs exhibit reduced calcium activity, increasing calcium activity rescues the differentiation defects in vitro. Deletion of a TET1-5hmC target gene, Itpr2, impairs the onset of OPC differentiation. Together, our results suggest that stage-specific TET1-mediated epigenetic programming and intracellular signaling are important for proper myelination and remyelination in mice.

Epigenetic mapping of the somatotropic axis in Nile tilapia reveals differential DNA hydroxymethylation marks associated with growth

In vertebrates, the somatotropic axis comprising the pituitary gland, liver and muscle plays a major role in myogenesis. Its output in terms of muscle growth is highly affected by nutritional and environmental cues, and thus likely epigenetically regulated. Hydroxymethylation is emerging as a DNA modification that modulates gene expression but a holistic characterization of the hydroxymethylome of the somatotropic axis has not been investigated to date. Using reduced representation 5-hydroxymethylcytosine profiling we demonstrate tissue-specific localization of 5-hydroxymethylcytosines at single nucleotide resolution. Their abundance within gene bodies and promoters of several growth-related genes supports their pertinent role in gene regulation. We propose that cytosine hydroxymethylation may contribute to the phenotypic plasticity of growth through epigenetic regulation of the somatotropic axis.

Genome-wide detection of cytosine methylation by single molecule real-time sequencing

Single molecule real-time (SMRT) sequencing theoretically offers the opportunity to directly assess certain base modifications of native DNA molecules without any prior chemical/enzymatic conversions and PCR amplification, using kinetic signals of a DNA polymerase. However, the kinetic signal changes caused by 5mC modification are extremely subtle. Hence, the robust genome-wide measurement of 5mC modification has not been achieved. We enhanced 5mC detection using SMRT sequencing by holistically analyzing kinetic signals of a DNA polymerase and sequence context for every base within a measurement window. We employed a convolutional neural network to train a methylation classification model, leading to genome-wide 5mC detection. The sensitivity and specificity reached 90% and 94%, with a 99% correlation of overall methylation level with bisulfite sequencing.

5-Methylcytosine (5mC) is an important type of epigenetic modification. Bisulfite sequencing (BS-seq) has limitations, such as severe DNA degradation. Using single molecule real-time sequencing, we developed a methodology to directly examine 5mC. This approach holistically examined kinetic signals of a DNA polymerase (including interpulse duration and pulse width) and sequence context for every nucleotide within a measurement window, termed the holistic kinetic (HK) model. The measurement window of each analyzed double-stranded DNA molecule comprised 21 nucleotides with a cytosine in a CpG site in the center. We used amplified DNA (unmethylated) and M.SssI-treated DNA (methylated) (M.SssI being a CpG methyltransferase) to train a convolutional neural network. The area under the curve for differentiating methylation states using such samples was up to 0.97. The sensitivity and specificity for genome-wide 5mC detection at single-base resolution reached 90% and 94%, respectively. The HK model was then tested on human–mouse hybrid fragments in which each member of the hybrid had a different methylation status. The model was also tested on human genomic DNA molecules extracted from various biological samples, such as buffy coat, placental, and tumoral tissues. The overall methylation levels deduced by the HK model were well correlated with those by BS-seq (r = 0.99; P < 0.0001) and allowed the measurement of allele-specific methylation patterns in imprinted genes. Taken together, this methodology has provided a system for simultaneous genome-wide genetic and epigenetic analyses.

Prenatal Lead (Pb) Exposure and Peripheral Blood DNA Methylation (5mC) and Hydroxymethylation (5hmC) in Mexican Adolescents from the ELEMENT Birth Cohort

BACKGROUND

Gestational lead (Pb) exposure can adversely affect offspring health through multiple mechanisms, including epigenomic alterations via DNA methylation (5mC) and hydroxymethylation (5hmC), an intermediate in oxidative demethylation. Most current methods do not distinguish between 5mC and 5hmC, limiting insights into their individual roles.

OBJECTIVE

Our study sought to identify the association of trimester-specific (T1, T2, T3) prenatal Pb exposure with 5mC and 5hmC levels at multiple cytosine-phosphate-guanine sites within gene regions previously associated with prenatal Pb (HCN2, NINJ2, RAB5A, TPPP) in whole blood leukocytes of children ages 11-18 years of age.

METHODS

Participants from the Early Life Exposure in Mexico to Environmental Toxicants (ELEMENT) birth cohorts were selected (n=144) for pyrosequencing analysis following oxidative or standard sodium bisulfite treatment. This workflow directly quantifies total methylation (5mC+5hmC) and 5mC only; 5hmC is estimated by subtraction.

RESULTS

Participants were 51% male, and mean maternal blood lead levels (BLL) were 6.43±5.16μg/dL in Trimester 1 (T1), 5.66±5.21μg/dL in Trimester 2 (T2), and 5.86±4.34μg/dL in Trimester 3 (T3). In addition, 5hmC levels were calculated for HCN2 (mean±standard deviation(SD), 2.08±4.18%), NINJ2 (G/C: 2.01±5.95; GG: 0.90±3.97), RAB5A (0.66±0.80%), and TPPP (1.11±6.67%). Furthermore, 5mC levels were measured in HCN2 (81.3±9.63%), NINJ2 (heterozygotes: 38.6±7.39%; GG homozygotes: 67.3±9.83%), RAB5A (1.41±1.21%), and TPPP (92.5±8.03%). Several significant associations between BLLs and 5mC/5hmC were identified: T1 BLLs with 5mC in HCN2 (β=-0.37, p=0.03) and 5hmC in NINJ2 (β=0.49, p=0.003); T2 BLLs with 5mC in HCN2 (β=0.37, p=0.03) and 5hmC in NINJ2 (β=0.27, p=0.008); and T3 BLLs with 5mC in HCN2 (β=0.50, p=0.01) and NINJ2 (β=-0.35, p=0.004) and 5hmC in NINJ2 (β=0.45, p<0.001). NINJ2 5mC was negatively correlated with gene expression (Pearson r=-0.5, p-value=0.005), whereas 5hmC was positively correlated (r=0.4, p-value=0.04).

DISCUSSION

These findings suggest there is variable 5hmC in human whole blood and that prenatal Pb exposure is associated with gene-specific 5mC and 5hmC levels at adolescence, providing evidence to consider 5hmC as a regulatory mechanism that is responsive to environmental exposures. https://doi.org/10.1289/EHP8507.

Ten-eleven translocation protein 1 modulates medulloblastoma progression

BACKGROUND

Medulloblastoma (MB) is the most common malignant pediatric brain tumor that originates in the cerebellum and brainstem. Frequent somatic mutations and deregulated expression of epigenetic regulators in MB highlight the substantial role of epigenetic alterations. 5-hydroxymethylcytosine (5hmC) is a highly abundant cytosine modification in the developing cerebellum and is regulated by ten-eleven translocation (TET) enzymes.

RESULTS

We investigate the alterations of 5hmC and TET enzymes in MB and their significance to cerebellar cancer formation. We show total abundance of 5hmC is reduced in MB, but identify significant enrichment of MB-specific 5hmC marks at regulatory regions of genes implicated in stem-like properties and Nanog-binding motifs. While TET1 and TET2 levels are high in MBs, only knockout of Tet1 in the smoothened (SmoA1) mouse model attenuates uncontrolled proliferation, leading to a favorable prognosis. The pharmacological Tet1 inhibition reduces cell viability and platelet-derived growth factor signaling pathway-associated genes.

CONCLUSIONS

These results together suggest a potential key role of 5hmC and indicate an oncogenic nature for TET1 in MB tumorigenesis, suggesting it as a potential therapeutic target for MBs.