Search

Transition of cytosine to thymine in CpG dinucleotides is the most frequent type of mutation in cancer. This increased mutability is commonly attributed to the spontaneous deamination of 5-methylcytosine (5mC), which is normally repaired by the base-excision repair (BER) pathway. However, the contribution of 5mC deamination in the increasing diversity of cancer mutational signatures remains poorly explored. We integrate mutational signatures analysis in a large series of tumor whole genomes with lineage-specific epigenomic data to draw a detailed view of 5mC deamination in cancer. We uncover tumor type-specific patterns of 5mC deamination signatures in CpG and non-CpG contexts. We demonstrate that the BER glycosylase MBD4 preferentially binds to active chromatin and early replicating DNA, which correlates with lower mutational burden in these domains. We validate our findings by modeling BER deficiencies in isogenic cell models. Here, we establish MBD4 as the main actor responsible for 5mC deamination repair in humans., Competing Interests: Competing interests D. Rieke reports advisory agreement with BeiGene and Bayer, honoraria from Bristol Myers Squibb, Bayer and Roche, research support from Seagen, and personal fees from Bayer and Johnson & Johnson, all outside the submitted work. A. Picca reports personal fees from AstraZeneca and Servier, all outside the submitted work. F. Bielle reports funding of research from Abbvie, service agreement for research contracted between his institution and Treefrog Therapeutics as well as Owkin, personal fees from Bristol Myers Squibb and a next-of-kin employed by Bristol Myers Squibb, all outside the submitted work. M.L. Yaspo is COO/CSO and shareholder of Alacris Theranostics without conflict of interest with the submitted work. M. Rodrigues reports non-financial support from AstraZeneca and Merck Sharp and Dohme, grants from Daiichi Sankyo, personal fees from AstraZeneca, Immunocore, Merck Sharp and Dohme and GlaxoSmithKline, all outside the submitted work. M.-H. Stern reports grants from Immunocore and Bionano, and royalties from Myriad Genetics, all outside the submitted work. The remaining authors have no conflict of interest to declare., (© 2024. The Author(s).)

5-Methylcytosine (m 5 C) is closely associated with cancer. However, the role of m 5 C in breast cancer(BC)remains unclear. This study combined single-cell RNA sequencing (scRNA-Seq) and transcriptomics datasets to screen m 5 C regulators associated with BC progression and analyze their clinical values. Firstly, This study elucidates the mechanisms of the m 5 C landscape and the specific roles of m 5 C regulators in BC patients. we found that the dysregulation of m 5 C regulators with m 5 Cscore play the essential role of the carcinogenesis and progression in epithelial cells and myeloid cells of BC at single cell level. External validation was conducted using an independent scRNA-Seq datasets. Then, three distinct m 5 C modification patterns were identified by transcriptomics datasets. Based on the m 5 C differentially expressed regulators, the m 5 Cscore was constructed, and used to divide patients with BC into high and low m 5 Cscore groups. Patients with a high m 5 Cscore had more abundant immune cell infiltration, stronger antitumor immunity, and better prognoses. Finally, Quantitative real-time (PCR) and immunohistochemistry were used for the in vitro experimental validation, which had extensive prognostic value. In this study, we aimed to assess the expression of m 5 C regulators involved in BC and investigate their correlation with the tumor microenvironment, clinicopathological characteristics, and prognosis of BC. The m 5 C regulators could be used to effectively assess the cell specific regulation prognosis of patients with BC and develop more effective immunotherapy strategies., Competing Interests: Declarations Ethics approval and consent to participate The studies involving human participants were reviewed and approved by Ethics Committee of the First Hospital of Shanxi Medical University, the ethics number is K-K0109. The patients provided their written informed consent to participate in this study. Consent for publication The patients involved have obtained ethical approval and written informed consent for the publication of any potentially identifiable images or data included in this article. Competing interests The authors declare no competing interests. Statement All methods were performed in accordance with the relevant guidelines and regulations., (© 2024. The Author(s).)

5-Methyldeoxycytidine (5mC) is a major epigenetic marker that regulates cellular functions in mammals. Endogenous lipid peroxidation can convert 5mC into 3,N4-etheno-5-methylcytosine (ϵ5mC). ϵ5mC is structurally similar to the mutagenic analog 3,N4-ethenocytosine (ϵC), which is repaired by AlkB family enzymes in the direct reversal repair (DRR) pathway and excised by DNA glycosylases in the base excision repair (BER) pathway. However, the repair of ϵ5mC has not been reported. Here, we examined the activities against ϵ5mC by DRR and BER enzymes and TET1-3, enzymes that modify the 5-methyl group in 5mC. We found that the etheno modification of 5mC blocks oxidation by TET1-3. Conversely, three human homologs in the AlkB family, ALKBH2, 3 and FTO were able to repair ϵ5mC to 5mC, which was subsequently modified by TET1 to 5-hydroxymethylcytosine. We also demonstrated that ALKBH2 likely repairs ϵ5mC in MEF cells. Another homolog, ALKBH5, could not repair ϵ5mC. Also, ϵ5mC is not a substrate for BER glycosylases SMUG1, AAG, or TDG. These findings indicate DRR committed by ALKBH2, 3 and FTO could reduce the detrimental effects of ϵ5mC in genetics and epigenetics and may work together with TET enzymes to modulate epigenetic regulations., (© The Author(s) 2024. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Epigenetics plays a pivotal role in regulating gene expression and cellular differentiation. DNA methylation, involving the addition of methyl groups to specific cytosine bases, is a well-known epigenetic modification. The recent discovery of 5-hydroxymethylcytosine (5hmC) has provided new insights into cytosine modifications. 5hmC, derived from the oxidation of 5-methylcytosine (5mC), serves as both an intermediate in demethylation and a stable chemical modification in the genome. In this comprehensive review, we summarize the recent research advancements regarding the functions of 5hmC in development and disease. We discuss its implications in gene expression regulation, cellular differentiation, and its potential role as a diagnostic and prognostic marker in various diseases. Additionally, we highlight the challenges associated with accurately detecting and quantifying 5hmC and present the latest methodologies employed for its detection. Understanding the functional role of 5hmC in epigenetic regulation and further advancing our understanding of gene expression dynamics and cellular processes hold immense promise for the development of novel therapeutic strategies and precision medicine approaches.

Maintenance methylation, of palindromic CpG dinucleotides at DNA replication forks, is crucial for the faithful mitotic inheritance of genomic 5-methylcytosine (5mC) methylation patterns. MBD proteins use two arginine residues to recognize symmetrically-positioned methyl groups in fully-methylated 5mCpG/5mCpG and 5mCpA/TpG dinucleotides. In contrast, C2H2 zinc finger (ZF) proteins recognize CpG and CpA, whether methylated or not, within longer specific sequences in a site- and strand-specific manner. Unmethylated CpG sites, often within CpG island (CGI) promoters, need protection by protein factors to maintain their hypomethylated status. Members of the BEN domain proteins bind CGCG or CACG elements within CGIs to regulate gene expression. Despite their overall structural diversity, MBD, ZF and BEN proteins all use arginine residues to recognize guanine, adopting either a 'straight-on' or 'oblique' conformation. The straight-on conformation accommodates a methyl group in the (5mC/T)pG dinucleotide, while the oblique conformation can clash with the methyl group of 5mC, leading to preferential binding of unmethylated sequences., (© The Author(s) 2024. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Human thymine-DNA glycosylase (TDG) excises T mispaired with G in a CpG context to initiate the base excision repair (BER) pathway. TDG is also involved in epigenetic regulation of gene expression by participating in active DNA demethylation. Here we demonstrate that under extended incubation time the full-length TDG (TDGFL), but neither its isolated catalytic domain (TDGcat) nor methyl-CpG binding domain-containing protein 4 (MBD4) DNA glycosylase, exhibits significant excision activity towards T and C in regular non-damaged DNA duplex in TpG/CpA and CpG/CpG contexts. Time course of the cleavage product accumulation under single-turnover conditions shows that the apparent rate constant for TDGFL-catalysed excision of T from T•A base pairs (0.0014-0.0069 min-1) is 85-330-fold lower than for the excision of T from T•G mispairs (0.47-0.61 min-1). Unexpectedly, TDGFL, but not TDGcat, exhibits prolonged enzyme survival at 37°C when incubated in the presence of equimolar concentrations of a non-specific DNA duplex, suggesting that the disordered N- and C-terminal domains of TDG can interact with DNA and stabilize the overall conformation of the protein. Notably, TDGFL was able to excise 5-hydroxymethylcytosine (5hmC), but not 5-methylcytosine residues from duplex DNA with the efficiency that could be physiologically relevant in post-mitotic cells. Our findings demonstrate that, under the experimental conditions used, TDG catalyses sequence context-dependent removal of T, C and 5hmC residues from regular DNA duplexes. We propose that in vivo the TDG-initiated futile DNA BER may lead to formation of persistent single-strand breaks in non-methylated or hydroxymethylated chromatin regions., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2024 Manapkyzy et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

Background: SLE is a complex autoimmune disease with heterogeneous manifestations and unpredictable outcomes. Early diagnosis is challenging due to non-specific symptoms, and current treatments only manage symptoms. Epigenetic alternations, including 5-Hydroxymethylome (5hmC) modifications, are important contributors to SLE pathogenesis. However, the 5hmC modification status in circulating cell-free DNA (cfDNA) of patients with SLE remains largely unexplored. We investigated the distribution of 5hmC in cfDNA of patients with SLE and healthy controls (HCs), and explored its potential as an SLE diagnosis marker., Methods: We used 5hmC-Seal to generate genome-wide 5hmC profiles of plasma cfDNA and bioinformatics analysis to screen differentially hydroxymethylated regions (DhMRs). In vitro mechanistic exploration was conducted to investigate the regulatory effect of CCCTC-binding factor (CTCF) in 5hmC candidate biomarkers., Results: We found distinct differences in genomic regions and 5hmC modification motif patterns between patients with SLE and HCs, varying with disease progression. Increased 5hmC modification enrichment was detected in SLE. Additionally, we screened 151 genes with hyper-5hmC, which are significantly involved in SLE-related processes, and 5hmC-modified BCL2 , CD83 , ETS1 and GZMB as SLE biomarkers. Our findings suggest that CTCF regulates 5hmC modification of these genes by recruiting TET (ten-eleven translocation) protein, and CTCF knockdown affected the protein expression of these genes in vitro., Conclusions: Our findings demonstrate the increased 5hmC distribution in plasma cfDNA in different disease activity in patients with SLE compared with HCs and relating DhMRs involved in SLE-associated pathways. Furthermore, we identified a panel of SLE relevant biomarkers, and these viewpoints could provide insight into the pathogenesis of SLE., Competing Interests: Competing interests: None declared., (© Author(s) (or their employer(s)) 2024. Re-use permitted under CC BY. Published by BMJ.)

Mutation of tet methylcytosine dioxygenase 2 (encoded by TET2) drives myeloid malignancy initiation and progression 1-3 . TET2 deficiency is known to cause a globally opened chromatin state and activation of genes contributing to aberrant haematopoietic stem cell self-renewal 4,5 . However, the open chromatin observed in TET2-deficient mouse embryonic stem cells, leukaemic cells and haematopoietic stem and progenitor cells 5 is inconsistent with the designated role of DNA 5-methylcytosine oxidation of TET2. Here we show that chromatin-associated retrotransposon RNA 5-methylcytosine (m 5 C) can be recognized by the methyl-CpG-binding-domain protein MBD6, which guides deubiquitination of nearby monoubiquitinated Lys119 of histone H2A (H2AK119ub) to promote an open chromatin state. TET2 oxidizes m 5 C and antagonizes this MBD6-dependent H2AK119ub deubiquitination. TET2 depletion thereby leads to globally decreased H2AK119ub, more open chromatin and increased transcription in stem cells. TET2-mutant human leukaemia becomes dependent on this gene activation pathway, with MBD6 depletion selectively blocking proliferation of TET2-mutant leukaemic cells and largely reversing the haematopoiesis defects caused by Tet2 loss in mouse models. Together, our findings reveal a chromatin regulation pathway by TET2 through retrotransposon RNA m 5 C oxidation and identify the downstream MBD6 protein as a feasible target for developing therapies specific against TET2 mutant malignancies., (© 2024. The Author(s).)

Among over 170 different types of chemical modifications on RNA nucleobases identified so far, RNA methylation is the major type of epitranscriptomic modifications existing on almost all types of RNAs, and has been demonstrated to participate in the entire process of RNA metabolism, including transcription, pre-mRNA alternative splicing and maturation, mRNA nucleus export, mRNA degradation and stabilization, mRNA translation. Attributing to the development of high-throughput detection technologies and the identification of both dynamic regulators and recognition proteins, mechanisms of RNA methylation modification in regulating the normal development of the organism as well as various disease occurrence and developmental abnormalities upon RNA methylation dysregulation have become increasingly clear. Here, we particularly focus on three types of RNA methylations: N 6 -methylcytosine (m 6 A), 5-methylcytosine (m 5 C), and N 7 -methyladenosine (m 7 G). We summarize the elements related to their dynamic installment and removal, specific binding proteins, and the development of high-throughput detection technologies. Then, for a comprehensive understanding of their biological significance, we also overview the latest knowledge on the underlying mechanisms and key roles of these three mRNA methylation modifications in gametogenesis, embryonic development, immune system development, as well as disease and tumor progression., (© 2024. Science China Press.)

Temporal lobe epilepsy (TLE) is a type of focal epilepsy characterized by spontaneous recurrent seizures originating from the hippocampus. The epigenetic reprogramming hypothesis of epileptogenesis suggests that the development of TLE is associated with alterations in gene transcription changes resulting in a hyperexcitable network in TLE. DNA 5-methylcytosine (5-mC) is an epigenetic mechanism that has been associated with chronic epilepsy. However, the contribution of 5-hydroxymethylcytosine (5-hmC), a product of 5-mC demethylation by the Ten-Eleven Translocation (TET) family proteins in chronic TLE is poorly understood. 5-hmC is abundant in the brain and acts as a stable epigenetic mark altering gene expression through several mechanisms. Here, we found that the levels of bulk DNA 5-hmC but not 5-mC were significantly reduced in the hippocampus of human TLE patients and in the kainic acid (KA) TLE rat model. Using 5-hmC hMeDIP-sequencing, we characterized 5-hmC distribution across the genome and found bidirectional regulation of 5-hmC at intergenic regions within gene bodies. We found that hypohydroxymethylated 5-hmC intergenic regions were associated with several epilepsy-related genes, including Gal, SV2, and Kcnj11 and hyperdroxymethylation 5-hmC intergenic regions were associated with Gad65, TLR4, and Bdnf gene expression. Mechanistically, Tet1 knockdown in the hippocampus was sufficient to decrease 5-hmC levels and increase seizure susceptibility following KA administration. In contrast, Tet1 overexpression in the hippocampus resulted in increased 5-hmC levels associated with improved seizure resiliency in response to KA. These findings suggest an important role for 5-hmC as an epigenetic regulator of epilepsy that can be manipulated to influence seizure outcomes., Competing Interests: Declaration of competing interest The authors declare the following financial interests/personal relationships which may be considered as potential competing interests. Farah Lubin reports financial support was provided by National Institute of Neurological Disorders and Stroke. Lawrence Ver Hoef reports was provided by National Institute of Neurological Disorders and Stroke. If there are other authors, they declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 The Authors. Published by Elsevier Inc. All rights reserved.)

DNA demethylases TET2 and TET3 play a fundamental role in thymic invariant natural killer T (iNKT) cell differentiation by mediating DNA demethylation of genes encoding for lineage specifying factors. Paradoxically, differential gene expression analysis revealed that significant number of genes were upregulated upon TET2 and TET3 loss in iNKT cells. This unexpected finding could be potentially explained if loss of TET proteins was reducing the expression of proteins that suppress gene expression. In this study, we discover that TET2 and TET3 synergistically regulate Drosha expression, by generating 5hmC across the gene body and by impacting chromatin accessibility. As DROSHA is involved in microRNA biogenesis, we proceed to investigate the impact of TET2/3 loss on microRNAs in iNKT cells. We report that among the downregulated microRNAs are members of the Let-7 family that downregulate in vivo the expression of the iNKT cell lineage specifying factor PLZF. Our data link TET proteins with microRNA expression and reveal an additional layer of TET mediated regulation of gene expression., Competing Interests: Dr. TÄ is a Director of Data Science at Covera Health. No funding from Covera Health was received for this study. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2024 Gioulbasani, Äijö, Valenzuela, Bettes and Tsagaratou.)

A major challenge faced by Vibrio cholerae is constant predation by bacteriophage (phage) in aquatic reservoirs and during infection of human hosts. To overcome phage predation, V. cholerae has acquired and/or evolved a myriad of phage defense systems. Although several novel defense systems have been discovered, we hypothesized that more were encoded in V. cholerae given the low diversity of phages that have been isolated, which infect this species. Using a V. cholerae genomic library, we identified a Type IV restriction system consisting of two genes within a 16-kB region of the Vibrio pathogenicity island-2, which we name TgvA and TgvB ( T ype I-embedded g mrSD -like system of V PI-2). We show that both TgvA and TgvB are required for defense against T2, T4, and T6 by targeting glucosylated 5-hydroxymethylcytosine (5hmC). T2 or T4 phages that lose the glucose modifications are resistant to TgvAB defense but exhibit a significant evolutionary tradeoff, becoming susceptible to other Type IV restriction systems that target unglucosylated 5hmC. We also show that the Type I restriction-modification system that embeds the tgvAB genes protects against phage T3, secΦ18, secΦ27, and λ, suggesting that this region is a phage defense island. Our study uncovers a novel Type IV restriction system in V. cholerae , increasing our understanding of the evolution and ecology of V. cholerae, while highlighting the evolutionary interplay between restriction systems and phage genome modification.IMPORTANCEBacteria are constantly being predated by bacteriophage (phage). To counteract this predation, bacteria have evolved a myriad of defense systems. Some of these systems specifically digest infecting phage by recognizing unique base modifications present on the phage DNA. In this study, we discover a Type IV restriction system encoded in V. cholerae, which we name TgvAB, and demonstrate it recognizes and restricts phage that have 5-hydroxymethylcytosine glucosylated DNA. Moreover, the evolution of resistance to TgvAB render phage susceptible to other Type IV restriction systems, demonstrating a significant evolutionary tradeoff. These results enhance our understanding of the evolution of V. cholerae and more broadly how bacteria evade phage predation., Competing Interests: The authors declare no conflict of interest.

Non-enzymatic spontaneous deamination of 5-methylcytosine, producing thymine, is the proposed etiology of cancer mutational signature 1, which is the most predominant signature in all cancers. Here, the proposed mutational process was reconstituted using synthetic DNA and purified proteins. First, single-stranded DNA containing 5-methylcytosine at CpG context was incubated at an elevated temperature to accelerate spontaneous DNA damage. Then, the DNA was treated with uracil DNA glycosylase to remove uracil residues that were formed by deamination of cytosine. The resulting DNA was then used as a template for DNA synthesis by yeast DNA polymerase δ. The DNA products were analyzed by next-generation DNA sequencing, and mutation frequencies were quantified. The observed mutations after this process were exclusively C>T mutations at CpG context, which was very similar to signature 1. When 5-methylcytosine modification and uracil DNA glycosylase were both omitted, C>T mutations were produced on C residues in all sequence contexts, but these mutations were diminished by uracil DNA glycosylase-treatment. These results indicate that the CpG>TpG mutations were produced by the deamination of 5-methylcytosine. Additional mutations, mainly C>G, were introduced by yeast DNA polymerase ζ on the heat-damaged DNA, indicating that G residues of the templates were also damaged. However, the damage on G residues was not converted to mutations with DNA polymerase δ or ε., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2024 Tomohiko Sugiyama. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

Background: N6-methyladenosine (m 6 A) and 5-methylcytosine (m 5 C) play a role in modifying long non-coding RNAs (lncRNAs) implicated in tumorigenesis and progression. This study was performed to evaluate prognostic value of m 6 A- and m 5 C-related lncRNAs and develop an efficient model for prognosis prediction in cervical cancer (CC)., Methods: Using gene expression data of TCGA set, we identified m 6 A- and m 5 C-related lncRNAs. Consensus Clustering Analysis was performed for samples subtyping based on survival-related lncRNAs, followed by analyzing tumor infiltrating immune cells (TIICs). Optimal signature lncRNAs were obtained using lasso Cox regression analysis for constructing a prognostic model and a nomogram to predict prognosis., Results: We built a co-expression network of 23 m 6 A-related genes, 15 m 5 C-related genes, and 62 lncRNAs. Based on 9 m 6 A- and m 5 C-related lncRNAs significantly associated with overall survival (OS) time, two molecular subtypes were obtained, which had significantly different OS time and fractions of TIICs. A prognostic model based on six m 6 A- and m 5 C-related signature lncRNAs was constructed, which could dichotomize patients into two risk subgroups with significantly different OS time. Prognostic power of the model was successfully validated in an independent dataset. We subsequently constructed a nomogram which could accurately predict survival probabilities. Drug sensitivity analysis found preferred chemotherapeutic agents for high and low-risk patients, respectively., Conclusion: Our study reveals that m 6 A- and m 5 C-related lncRNAs are associated with prognosis and immune microenvironment of CC. The m 6 A- and m 5 C-related six-lncRNA signature may be a useful tool for survival stratification in CC and open new avenues for individualized therapies., (© 2024. The Author(s).)

Both DNA 5-methylcytosine (5mC) and RNA N6-methyladenosine (m 6 A) modifications are reported to participate in cellular stress responses including inflammation. Phosphoenolpyruvate carboxykinase 2 (PCK2) is upregulated in Kupffer cells (KCs) to facilitate the proinflammatory phosphorylation signaling cascades upon LPS stimulation, yet the role of 5mC and m 6 A in PCK2 upregulation remain elusive. Here, we report that the significantly augmented PCK2 mRNA and protein levels are associated with global 5mC demethylation coupled with m 6 A hypermethylation in LPS-activated KCs. The suppression of 5mC demethylation or m 6 A hypermethylation significantly alleviates the upregulation of PCK2 and proinflammatory cytokines in LPS-challenged KCs. Further reciprocal tests indicate 5mC demethylation is upstream of m 6 A hypermethylation. Specifically, CpG islands in the promoters of PCK2 and RNA methyltransferase (METTL3 and METTL14) genes are demethylated, while the 3'UTR of PCK2 mRNA is m6A hypermethylated, in LPS-stimulated KCs. These modifications contribute to the transactivation of the PCK2 gene as well as increased PCK2 mRNA stability and protein production via a m 6 A-mediated mechanism with IGF2BP1 as the reader protein. These results indicate that DNA 5mC and RNA m6A collaborate to upregulate PCK2 expression, respectively, at the transcriptional and post-transcriptional levels during KC activation.

Background: Breast tumorigenesis is a complex and multistep process accompanied by both genetic and epigenetic dysregulation. In contrast to the extensive studies on DNA epigenetic modifications 5-hydroxymethylcytosine (5hmC) and 5-methylcytosine (5mC) in malignant breast tumors, their roles in the early phases of breast tumorigenesis remain ambiguous., Results: DNA 5hmC and 5mC exhibited a consistent and significant decrease from usual ductal hyperplasia to atypical ductal hyperplasia and subsequently to ductal carcinoma in situ (DCIS). However, 5hmC showed a modest increase in invasive ductal breast cancer compared to DCIS. Genomic analyses showed that the changes in 5hmC and 5mC levels occurred around the transcription start sites (TSSs), and the modification levels were strongly correlated with gene expression levels. Meanwhile, it was found that differentially hydroxymethylated regions (DhMRs) and differentially methylated regions (DMRs) were overlapped in the early phases and accompanied by the enrichment of active histone marks. In addition, TET2-related DNA demethylation was found to be involved in breast tumorigenesis, and four transcription factor binding sites (TFs: ESR1, FOXA1, GATA3, FOS) were enriched in TET2-related DhMRs/DMRs. Intriguingly, we also identified a certain number of common DhMRs between tumor samples and cell-free DNA (cfDNA)., Conclusions: Our study reveals that dynamic changes in DNA 5hmC and 5mC play a vital role in propelling breast tumorigenesis. Both TFs and active histone marks are involved in TET2-related DNA demethylation. Concurrent changes in 5hmC signals in primary breast tumors and cfDNA may play a promising role in breast cancer screening., (© 2024. The Author(s).)

Background: By regulating the functions of multiple RNAs, 5-methylcytosine (m 5 C) RNA methylation, particularly mediated by NOP2, is involved in tumorigenesis and developments. However, the specific functions and potential mechanisms of m 5 C, especially involving NOP2, in clear-cell renal cell carcinoma (ccRCC), remain unclear. Methods: NOP2 expression in cell lines and patient tissues was detected using western blotting, quantitative real-time polymerase chain reaction (RT-qPCR), and immunohistochemistry. The biological effects of NOP2 on ccRCC cells were investigated through a series of in vitro and in vivo experiments. To explore the potential regulatory mechanisms by which NOP2 affects ccRCC progression, m 5 C bisulfite sequencing, RNA-sequencing, RNA immunoprecipitation and methylated RNA immunoprecipitation (RIP/MeRIP) RT-qPCR assay, luciferase reporter assay, RNA stability assay, and bioinformatic analysis were performed. Results: NOP2 expression was significantly upregulated in ccRCC tissues and was associated with poor prognosis. Moreover, loss-of-function and gain-of-function assays demonstrated that NOP2 altered ccRCC cell proliferation, migration, and invasion. Mechanistically, NOP2 stimulated m 5 C modification of apolipoprotein L1 (APOL1) mRNA, and m 5 C reader YBX1 stabilized APOL1 mRNA through recognizing and binding to m 5 C site in the 3'-untranslated regions. Silencing APOL1 expression inhibited ccRCC cell proliferation in vitro and tumor formation in vivo . Furthermore, NOP2/APOL1 affected ccRCC progression via the PI3K-Akt signaling pathway. Conclusion: NOP2 functions as an oncogene in ccRCC by promoting tumor progression through the m 5 C-dependent stabilization of APOL1, which in turn regulates the PI3K-Akt signaling pathway, suggesting a potential therapeutic target for ccRCC., Competing Interests: Competing Interests: The authors have declared that no competing interest exists., (© The author(s).)

RNA modifications play an important role in actively controlling recently created formation in cellular regulation mechanisms, which link them to gene expression and protein. The RNA modifications have numerous alterations, presenting broad glimpses of RNA's operations and character. The modification process by the TET enzyme oxidation is the crucial change associated with cytosine hydroxymethylation. The effect of CR is an alteration in specific biochemical ways of the organism, such as gene expression and epigenetic alterations. Traditional laboratory systems that identify 5-hydroxymethylcytosine (5hmC) samples are expensive and time-consuming compared to other methods. To address this challenge, the paper proposed XGB5hmC, a machine learning algorithm based on a robust gradient boosting algorithm (XGBoost), with different residue based formulation methods to identify 5hmC samples. Their results were amalgamated, and six different frequency residue based encoding features were fused to form a hybrid vector in order to enhance model discrimination capabilities. In addition, the proposed model incorporates SHAP (Shapley Additive Explanations) based feature selection to demonstrate model interpretability by highlighting the high contributory features. Among the applied machine learning algorithms, the XGBoost ensemble model using the tenfold cross-validation test achieved improved results than existing state-of-the-art models. Our model reported an accuracy of 89.97%, sensitivity of 87.78%, specificity of 94.45%, F1-score of 0.8934%, and MCC of 0.8764%. This study highlights the potential to provide valuable insights for enhancing medical assessment and treatment protocols, representing a significant advancement in RNA modification analysis., (© 2024. The Author(s).)

In the era of precision oncology, identifying abnormal genetic and epigenetic alterations has transformed the way cancer is diagnosed, managed, and treated. 5-hydroxymethylcytosine (5hmC) is an emerging epigenetic modification formed through the oxidation of 5-methylcytosine (5mC) by ten-eleven translocase (TET) enzymes. DNA hydroxymethylation exhibits tissue- and cancer-specific patterns and is essential in DNA demethylation and gene regulation. Recent advancements in 5hmC detection methods and the discovery of 5hmC in cell-free DNA (cfDNA) have highlighted the potential for cell-free 5hmC as a cancer biomarker. This review explores the current and emerging techniques and applications of DNA hydroxymethylation in cancer, particularly in the context of cfDNA.

Motivation: 5-Hydroxymethylcytosine (5hmC), a crucial epigenetic mark with a significant role in regulating tissue-specific gene expression, is essential for understanding the dynamic functions of the human genome. Despite its importance, predicting 5hmC modification across the genome remains a challenging task, especially when considering the complex interplay between DNA sequences and various epigenetic factors such as histone modifications and chromatin accessibility., Results: Using tissue-specific 5hmC sequencing data, we introduce Deep5hmC, a multimodal deep learning framework that integrates both the DNA sequence and epigenetic features such as histone modification and chromatin accessibility to predict genome-wide 5hmC modification. The multimodal design of Deep5hmC demonstrates remarkable improvement in predicting both qualitative and quantitative 5hmC modification compared to unimodal versions of Deep5hmC and state-of-the-art machine learning methods. This improvement is demonstrated through benchmarking on a comprehensive set of 5hmC sequencing data collected at four developmental stages during forebrain organoid development and across 17 human tissues. Compared to DeepSEA and random forest, Deep5hmC achieves close to 4% and 17% improvement of Area Under the Receiver Operating Characteristic (AUROC) across four forebrain developmental stages, and 6% and 27% across 17 human tissues for predicting binary 5hmC modification sites; and 8% and 22% improvement of Spearman correlation coefficient across four forebrain developmental stages, and 17% and 30% across 17 human tissues for predicting continuous 5hmC modification. Notably, Deep5hmC showcases its practical utility by accurately predicting gene expression and identifying differentially hydroxymethylated regions (DhMRs) in a case-control study of Alzheimer's disease (AD). Deep5hmC significantly improves our understanding of tissue-specific gene regulation and facilitates the development of new biomarkers for complex diseases., Availability and Implementation: Deep5hmC is available via https://github.com/lichen-lab/Deep5hmC., (© The Author(s) 2024. Published by Oxford University Press.)

Stunting is the result of chronic malnutrition due to the lack of micronutrient-based methyl donors required for epigenetic programming during the first 1000 days of life. Methylation studies using bisulfite conversion from blood DNA are invasive and may not be practical for large-scale epidemiological investigation or nutrition intervention programs. Buccal epithelial methylation may reflect early germline methylation. Therefore, buccal cells can serve as convenient sample sources to collect biomarkers associated with the risk of stunting. This study aims to describe the feasibility of nanopore adaptive sampling in detecting DNA methylation from children's buccal DNA. We used adaptive sampling of Oxford Nanopore Technology on barcoded samples to describe differential methylation associated with malnutrition. Overall, the level of 5-methylcytosine (5mC) was lower in stunted children than in normal children. We also found differentially methylated regions at the MIR6724 and RNA45SN1 gene loci on chromosome 21, which was higher in stunted children than in normal children. We described and detected differential DNA methylation in the locus previously not known to be associated with stunting. Interestingly, this locus on chromosome 21 has been implicated in the stunted phenotype of Down syndrome.

TET1/2/3 dioxygenases iteratively demethylate 5-methylcytosine, beginning with the formation of 5-hydroxymethylcytosine (5hmC). The post-mitotic brain maintains higher levels of 5hmC than most peripheral tissues, and TET1 ablation studies have underscored the critical role of TET1 in brain physiology. However, deletion of Tet1 precludes the disentangling of the catalytic and non-catalytic functions of TET1. Here, we dissect these functions of TET1 by comparing adult cortex of Tet1 wildtype ( Tet1 WT), a novel Tet1 catalytically dead mutant ( Tet1 HxD), and Tet1 knockout ( Tet1 KO) mice. Using DNA methylation array, we uncover that Tet1 HxD and KO mutations perturb the methylation status of distinct subsets of CpG sites. Gene ontology (GO) analysis on specific differential 5hmC regions indicates that TET1's catalytic activity is linked to neuronal-specific functions. RNA-Seq further shows that Tet1 mutations predominantly impact the genes that are associated with alternative splicing. Lastly, we performed High-performance Liquid Chromatography Mass-Spectrometry lipidomics on WT and mutant cortices and uncover accumulation of lysophospholipids lysophosphatidylethanolamine and lysophosphatidylcholine in Tet1 HxD cortex. In summary, we show that Tet1 HxD does not completely phenocopy Tet1 KO, providing evidence that TET1 modulates distinct cortical functions through its catalytic and non-catalytic roles.

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.

Machado-Joseph disease (MJD) is an autosomal dominant spinocerebellar ataxia (SCA) caused by a polyglutamine expansion in the ataxin-3 protein, which initiates a cascade of pathogenic events, including transcriptional dysregulation. Genotype-phenotype correlations in MJD are incomplete, suggesting an influence of additional factors, such as epigenetic modifications, underlying the MJD pathogenesis. DNA methylation is known to impact the pathophysiology of neurodegenerative disorders through gene expression regulation and increased methylation has been reported for other SCAs. In this work we aimed to analyse global methylation in MJD carriers. Global 5-mC levels were quantified in blood samples of 33 MJD mutation carriers (patients and preclinical subjects) and 33 healthy controls, matched by age, sex, and smoking status. For a subset of 16 MJD subjects, a pilot follow-up analysis with two time points was also conducted. No differences were found in median global 5-mC levels between MJD mutation carriers and controls and no correlations between methylation levels and clinical or genetic variables were detected. Also, no alterations in global 5-mC levels were observed over time. Our findings do not support an increase in global blood methylation levels associated with MJD.

As the most common and aggressive primary malignant brain tumor, glioblastoma is still lacking a satisfactory curative approach. The standard management consisting of gross total resection followed by radiotherapy and chemotherapy with temozolomide only prolongs patients' life moderately. In recent years, many therapeutics have failed to give a breakthrough in GBM treatment. In the search for new treatment solutions, we became interested in the repurposing of existing medicines, which have established safety profiles. We focused on the possible implementation of well-known drugs, metformin, and arginine. Metformin is widely used in diabetes treatment, but arginine is mainly a cardiovascular protective drug. We evaluated the effects of metformin and arginine on total DNA methylation, as well as the oxidative stress evoked by treatment with those agents. In glioblastoma cell lines, a decrease in 5-methylcytosine contents was observed with increasing drug concentration. When combined with temozolomide, both guanidines parallelly increased DNA methylation and decreased 8-oxo-deoxyguanosine contents. These effects can be explained by specific interactions of the guanidine group with m 5 CpG dinucleotide. We showed that metformin and arginine act on the epigenetic level, influencing the foreground and potent DNA regulatory mechanisms. Therefore, they can be used separately or in combination with temozolomide, in various stages of disease, depending on desired treatment effects.

Enzymatic modification of DNA nucleobases can coordinate gene expression, nuclease protection, or mutagenesis. We recently discovered a clade of phage-specific cytosine methyltransferase (MT) and 5-methylpyrimidine dioxygenase (5mYOX) enzymes that produce 5-hydroxymethylcytosine (5hmC) as a precursor for enzymatic hypermodifications on viral genomes. Here, we identify phage MT- and 5mYOX-associated glycosyltransferases (GTs) that catalyze linkage of diverse sugars to 5hmC nucleobase substrates. Metavirome mining revealed thousands of biosynthetic gene clusters containing enzymes with predicted roles in cytosine sugar hypermodification. We developed a platform for high-throughput screening of GT-containing pathways, relying on the Escherichia coli metabolome as a substrate pool. We successfully reconstituted several pathways and isolated diverse sugar modifications appended to cytosine, including mono-, di-, or tri-saccharides comprised of hexoses, N-acetylhexosamines, or heptose. These findings expand our knowledge of hypermodifications on nucleic acids and the origins of corresponding sugar-installing enzymes., Competing Interests: Declaration of interests J.D.P., S.R.L., K.H.O., and L.S. are current employees of New England Biolabs, a commercial producer and supplier of molecular biology reagents. This affiliation does not impact the adherence to ethical standards of the scientific community, experimental impartiality of the authors, or availability of the data presented. All authors are named inventors on a patent application related to the work described in this manuscript., (Copyright © 2024 The Author(s). Published by Elsevier Inc. All rights reserved.)

DNA methyltransferases are drug targets for myelodysplastic syndrome (MDS), chronic myelomonocytic leukemia (CMML), acute myelogenous leukemia (AML) and possibly β-hemoglobinopathies. We characterize the interaction of nucleoside analogues in DNA with a prokaryotic CpG-specific DNA methyltransferase (M.MpeI) as a model for mammalian DNMT1 methyltransferases. We tested DNA containing 5-hydroxymethylcytosine (5hmC), 5-hydroxycytosine (5OHC), 5-methyl-2-pyrimidinone (in the ribosylated form known as 5-methylzebularine, 5mZ), 5,6-dihydro-5-azacytosine (dhaC), 5-fluorocytosine (5FC), 5-chlorocytosine (5ClC), 5-bromocytosine (5BrC) and 5-iodocytosine (5IC). Covalent complex formation was by far most efficient for 5FC. Non-covalent complexes were most abundant for dhaC and 5mZ. Surprisingly, we observed methylation of 5IC and 5BrC, and to a lesser extent 5ClC and 5FC, in the presence, but not the absence of small molecule thiol nucleophiles. For 5IC and 5BrC, we demonstrated by mass spectrometry that the reactions were due to methyltransferase driven dehalogenation, followed by methylation. Crystal structures of M.MpeI-DNA complexes capture the 'in' conformation of the active site loop for analogues with small or rotatable (5mZ) 5-substituents and its 'out' form for bulky 5-substituents. Since very similar 'in' and 'out' loop conformations were also observed for DNMT1, it is likely that our conclusions generalize to other DNA methyltransferases., (© The Author(s) 2024. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Background: Obstructive sleep apnoea (OSA) is a sleep-disordered breathing characterized by intermittent hypoxia (IH) that may cause cognitive dysfunction. However, the impact of IH on molecular processes involved in cognitive function remains unclear., Methods: C57BL / 6 J mice were exposed to either normoxia (control) or IH for 6 weeks. DNA hydroxymethylation was quantified by hydroxymethylated DNA immunoprecipitation (hMeDIP) sequencing. ten-eleven translocation 1 (Tet1) was knocked down by lentivirus. Specifically, cognitive function was assessed by behavioral experiments, pathological features were assessed by HE staining, the hippocampal DNA hydroxymethylation was examined by DNA dot blot and immunohistochemical staining, while the Wnt signaling pathway and its downstream effects were studied using qRT-PCR, immunofluorescence staining, and Luminex liquid suspension chip analysis., Results: IH mice showed pathological changes and cognitive dysfunction in the hippocampus. Compared with the control group, IH mice exhibited global DNA hydroxylmethylation in the hippocampus, and the expression of three hydroxylmethylases increased significantly. The Wnt signaling pathway was activated, and the mRNA and 5hmC levels of Wnt3a, Ccnd2, and Prickle2 were significantly up-regulated. Further caused downstream neurogenesis abnormalities and neuroinflammatory activation, manifested as increased expression of IBA1 (a marker of microglia), GFAP (a marker of astrocytes), and DCX (a marker of immature neurons), as well as a range of inflammatory cytokines (e.g. TNFa, IL3, IL9, and IL17A). After Tet1 knocked down, the above indicators return to normal., Conclusion: Activation of Wnt signaling pathway by hippocampal Tet1 is associated with cognitive dysfunction induced by IH., (© 2024. The Author(s).)

While the significance of N6-methyladenosine (m 6 A) in viral regulation has been extensively studied, the functions of 5-methylcytosine (m 5 C) modification in viral biology remain largely unexplored. In this study, we demonstrate that m 5 C is more abundant than m 6 A in severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and provide a comprehensive profile of the m 5 C landscape of SARS-CoV-2 RNA. Knockout of NSUN2 reduces m 5 C levels in SARS-CoV-2 virion RNA and enhances viral replication. Nsun2 deficiency mice exhibited higher viral burden and more severe lung tissue damages. Combined RNA-Bis-seq and m 5 C-MeRIP-seq identified the NSUN2-dependent m 5 C-methylated cytosines across the positive-sense genomic RNA of SARS-CoV-2, and the mutations of these cytosines enhance RNA stability. The progeny SARS-CoV-2 virions from Nsun2 deficiency mice with low levels of m 5 C modification exhibited a stronger replication ability. Overall, our findings uncover the vital role played by NSUN2-mediated m 5 C modification during SARS-CoV-2 replication and propose a host antiviral strategy via epitranscriptomic addition of m 5 C methylation to SARS-CoV-2 RNA.

Tamoxifen, a selective estrogen receptor modulator (SERM), exhibits dual agonist or antagonist effects contingent upon its binding to either G-protein-coupled estrogen receptor (GPER) or estrogen nuclear receptor (ESR). Estrogen signaling plays a pivotal role in initiating epigenetic alterations and regulating estrogen-responsive genes in breast cancer. Employing three distinct breast cancer cell lines-MCF-7 (ESR+; GPER+), MDA-MB-231 (ESR-; GPER-), and SkBr3 (ESR-; GPER+)-this study subjected them to treatment with two tamoxifen derivatives: 4-hydroxytamoxifen (4-HT) and endoxifen (Endox). Through 2D high-performance liquid chromatography with tandem mass spectrometry detection (HPLC-MS/MS), varying levels of 5-methylcytosine (5-mC) were found, with MCF-7 displaying the highest levels. Furthermore, TET3 mRNA expression levels varied among the cell lines, with MCF-7 exhibiting the lowest expression. Notably, treatment with 4-HT induced significant changes in TET3 expression across all cell lines, with the most pronounced increase seen in MCF-7 and the least in MDA-MB-231. These findings underscore the influence of tamoxifen derivatives on DNA methylation patterns, particularly through modulating TET3 expression, which appears to be contingent on the presence of estrogen receptors. This study highlights the potential of targeting epigenetic modifications for personalized anti-cancer therapy, offering a novel avenue to improve treatment outcomes.

Epigenetic modifications, such as 5-methylcytosine (5mC), can sometimes be transmitted between generations, provoking speculation that epigenetic changes could play a role in adaptation and evolution. Here, we use experimental evolution to investigate how 5mC levels evolve in populations of biparental insect (Nicrophorus vespilloides) derived from a wild source population and maintained independently under different regimes of parental care in the lab. We show that 5mC levels in the transcribed regions of genes (gene bodies) diverge between populations that have been exposed to different levels of care for 30 generations. These changes in 5mC do not reflect changes in the levels of gene expression. However, the accumulation of 5mC within genes between populations is associated with reduced variability in gene expression within populations. Our results suggest that evolved change in 5mC could contribute to phenotypic evolution by influencing variability in gene expression in invertebrates., (© 2024. The Author(s).)

Inheritance of 5-methylcytosine from one cell generation to the next by DNA methyltransferase 1 (DNMT1) plays a key role in regulating cellular identity. While recent work has shown that the activity of DNMT1 is imprecise, it remains unclear how the fidelity of DNMT1 is tuned in different genomic and cell state contexts. Here we describe Dyad-seq, a method to quantify the genome-wide methylation status of cytosines at the resolution of individual CpG dinucleotides to find that the fidelity of DNMT1-mediated maintenance methylation is related to the local density of DNA methylation and the landscape of histone modifications. To gain deeper insights into methylation/demethylation turnover dynamics, we first extended Dyad-seq to quantify all combinations of 5-methylcytosine and 5-hydroxymethylcytosine at individual CpG dyads. Next, to understand how cell state transitions impact maintenance methylation, we scaled the method down to jointly profile genome-wide methylation levels, maintenance methylation fidelity and the transcriptome from single cells (scDyad&T-seq). Using scDyad&T-seq, we demonstrate that, while distinct cell states can substantially impact the activity of the maintenance methylation machinery, locally there exists an intrinsic relationship between DNA methylation density, histone modifications and DNMT1-mediated maintenance methylation fidelity that is independent of cell state., (© 2024. The Author(s), under exclusive licence to Springer Nature America, Inc.)

Background: DNA methylation in the form of 5-methylcytosine (5mC) is the most abundant base modification in animals. However, 5mC levels vary widely across taxa. While vertebrate genomes are hypermethylated, in most invertebrates, 5mC concentrates on constantly and highly transcribed genes (gene body methylation; GbM) and, in some species, on transposable elements (TEs), a pattern known as "mosaic". Yet, the role and developmental dynamics of 5mC and how these explain interspecies differences in DNA methylation patterns remain poorly understood, especially in Spiralia, a large clade of invertebrates comprising nearly half of the animal phyla., Results: Here, we generate base-resolution methylomes for three species with distinct genomic features and phylogenetic positions in Annelida, a major spiralian phylum. All possible 5mC patterns occur in annelids, from typical invertebrate intermediate levels in a mosaic distribution to hypermethylation and methylation loss. GbM is common to annelids with 5mC, and methylation differences across species are explained by taxon-specific transcriptional dynamics or the presence of intronic TEs. Notably, the link between GbM and transcription decays during development, alongside a gradual and global, age-dependent demethylation in adult stages. Additionally, reducing 5mC levels with cytidine analogs during early development impairs normal embryogenesis and reactivates TEs in the annelid Owenia fusiformis., Conclusions: Our study indicates that global epigenetic erosion during development and aging is an ancestral feature of bilateral animals. However, the tight link between transcription and gene body methylation is likely more important in early embryonic stages, and 5mC-mediated TE silencing probably emerged convergently across animal lineages., (© 2024. The Author(s).)

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., (© 2024. This is a U.S. Government work and not under copyright protection in the US; foreign copyright protection may apply.)

Multiomics require concerted recording of independent information, ideally from a single experiment. In this study, we introduce RIMS-seq2, a high-throughput technique to simultaneously sequence genomes and overlay methylation information while requiring only a small modification of the experimental protocol for high-throughput DNA sequencing to include a controlled deamination step. Importantly, the rate of deamination of 5-methylcytosine is negligible and thus does not interfere with standard DNA sequencing and data processing. Thus, RIMS-seq2 libraries from whole- or targeted-genome sequencing show the same germline variation calling accuracy and sensitivity compared with standard DNA-seq. Additionally, regional methylation levels provide an accurate map of the human methylome., (© 2024 Yan et al.; Published by Cold Spring Harbor Laboratory Press.)

5-Methylcytosine (5mC) is a widespread silencing mechanism that controls genomic parasites. In eukaryotes, 5mC has gained complex roles in gene regulation beyond parasite control, yet 5mC has also been lost in many lineages. The causes for 5mC retention and its genomic consequences are still poorly understood. Here, we show that the protist closely related to animals Amoebidium appalachense features both transposon and gene body methylation, a pattern reminiscent of invertebrates and plants. Unexpectedly, hypermethylated genomic regions in Amoebidium derive from viral insertions, including hundreds of endogenized giant viruses, contributing 14% of the proteome. Using a combination of inhibitors and genomic assays, we demonstrate that 5mC silences these giant virus insertions. Moreover, alternative Amoebidium isolates show polymorphic giant virus insertions, highlighting a dynamic process of infection, endogenization, and purging. Our results indicate that 5mC is critical for the controlled coexistence of newly acquired viral DNA into eukaryotic genomes, making Amoebidium a unique model to understand the hybrid origins of eukaryotic DNA.

Maps of the RNA modification 5-methylcytosine (m 5 C) often diverge markedly not only because of differences in detection methods, data depand analysis pipelines but also biological factors. We re-analysed bisulfite RNA sequencing datasets from five human cell lines and seven tissues using a coherent m 5 C site calling pipeline. With the resulting union list of 6,393 m 5 C sites, we studied site distribution, enzymology, interaction with RNA-binding proteins and molecular function. We confirmed tRNA:m 5 C methyltransferases NSUN2 and NSUN6 as the main mRNA m 5 C "writers," but further showed that the rRNA:m 5 C methyltransferase NSUN5 can also modify mRNA. Each enzyme recognises mRNA features that strongly resemble their canonical substrates. By analysing proximity between mRNA m 5 C sites and footprints of RNA-binding proteins, we identified new candidates for functional interactions, including the RNA helicases DDX3X, involved in mRNA translation, and UPF1, an mRNA decay factor. We found that lack of NSUN2 in HeLa cells affected both steady-state levels of, and UPF1-binding to, target mRNAs. Our studies emphasise the emerging diversity of m 5 C writers and readers and their effect on mRNA function., (© 2024 Guarnacci et al.)

Background: Radioresistance is the leading cause of death in advanced cervical cancer (CC). Dysregulation of RNA modification has recently emerged as a regulatory mechanism in radiation and drug resistance. We aimed to explore the biological function and clinical significance of 5-methylcytosine (m 5 C) in cervical cancer radiosensitivity., Methods: The abundance of RNA modification in radiotherapy-resistant and sensitive CC specimens was quantified by liquid chromatography-tandem mass spectrometry. The essential RNA modification-related genes involved in CC radiosensitivity were screened via RNA sequencing. The effect of NSUN6 on radiosensitivity was verified in CC cell lines, cell-derived xenograft (CDX), and 3D bioprinted patient-derived organoid (PDO). The mechanisms of NSUN6 in regulating CC radiosensitivity were investigated by integrative m 5 C sequencing, mRNA sequencing, and RNA immunoprecipitation., Results: We found a higher abundance of m 5 C modification in resistant CC samples, and NSUN6 was the essential m 5 C-regulating gene concerning radiosensitivity. NSUN6 overexpression was clinically correlated with radioresistance and poor prognosis in cervical cancer. Functionally, higher NSUN6 expression was associated with radioresistance in the 3D PDO model of cervical cancer. Moreover, silencing NSUN6 increased CC radiosensitivity in vivo and in vitro. Mechanistically, NDRG1 was one of the downstream target genes of NSUN6 identified by integrated m 5 C-seq, mRNA-seq, and functional validation. NSUN6 promoted the m 5 C modification of NDRG1 mRNA, and the m 5 C reader ALYREF bound explicitly to the m 5 C-labeled NDRG1 mRNA and enhanced NDRG1 mRNA stability. NDRG1 overexpression promoted homologous recombination-mediated DNA repair, which in turn led to radioresistance in cervical cancer., Conclusions: Aberrant m 5 C hypermethylation and NSUN6 overexpression drive resistance to radiotherapy in cervical cancer. Elevated NSUN6 expression promotes radioresistance in cervical cancer by activating the NSUN6/ALYREF-m 5 C-NDRG1 pathway. The low expression of NSUN6 in cervical cancer indicates sensitivity to radiotherapy and a better prognosis., (© 2024. The Author(s).)

DNA methylation plays an important role in various biological processes, including cell differentiation, ageing, and cancer development. The most important methylation in mammals is 5-methylcytosine mostly occurring in the context of CpG dinucleotides. Sequencing methods such as whole-genome bisulfite sequencing successfully detect 5-methylcytosine DNA modifications. However, they suffer from the serious drawbacks of short read lengths and might introduce an amplification bias. Here we present Rockfish, a deep learning algorithm that significantly improves read-level 5-methylcytosine detection by using Nanopore sequencing. Rockfish is compared with other methods based on Nanopore sequencing on R9.4.1 and R10.4.1 datasets. There is an increase in the single-base accuracy and the F1 measure of up to 5 percentage points on R.9.4.1 datasets, and up to 0.82 percentage points on R10.4.1 datasets. Moreover, Rockfish shows a high correlation with whole-genome bisulfite sequencing, requires lower read depth, and achieves higher confidence in biologically important regions such as CpG-rich promoters while being computationally efficient. Its superior performance in human and mouse samples highlights its versatility for studying 5-methylcytosine methylation across varied organisms and diseases. Finally, its adaptable architecture ensures compatibility with new versions of pores and chemistry as well as modification types., (© 2024. The Author(s).)

Background: As the most important modifications on the RNA level, N6-methyladenosine (m6A-) and 5-methylcytosine (m5C-) modification could have a direct influence on the RNAs. Long non-coding RNAs (lncRNAs) could also be modified by methylcytosine modification. Compared with mRNAs, the function of lncRNAs could be more potent to some extent in biological processes like tumorigenesis. Until now, rare reports have been done associated with cutaneous melanoma. Herein, we wonder if the m6A- and m5C- modified lncRNAs could influence the immune landscape and prognosis in melanoma, and we also want to find some lncRNAs which could directly affect the malignant behaviors of melanoma., Methods: Systematically, we explored the expression pattern of m6A- and m5C- modified lncRNAs in melanoma from datasets including UCSC Xena and NCBI GEO, and the prognostic lncRNAs were selected. Then, according to the expression pattern of lncRNAs, melanoma samples from these datasets were divided into several subtypes. Prognostic model, nomogram survival model, drug sensitivity, GO, and KEGG pathway analysis were performed. Furthermore, among several selected lncRNAs, we identified one lncRNA named LINC00893 and investigated its expression pattern and its biological function in melanoma cell lines., Results: We identified 27 m6A- and m5C- related lncRNAs which were significantly associated with survival, and we made a subtype analysis of melanoma samples based on these 27 lncRNAs. Among the two subtypes, we found differences of immune cells infiltration between these two subtypes. Then, LASSO algorithm was used to screen the optimized lncRNAs combination including ZNF252P-AS1, MIAT, FAM13A-AS1, LINC-PINT, LINC00893, AGAP2-AS1, OIP5-AS1, and SEMA6A-AS1. We also found that there was a significant correlation between the different risk groups predicted based on RS model and the actual prognosis. The nomogram survival model based on independent survival prognostic factors was also constructed. Besides, sensitivity to chemotherapeutic agents, GO and KEGG analysis were performed. In different risk groups, a total of 14 drug molecules with different distributions were obtained, which included AZD6482, AZD7762, AZD8055, camptothecin, dasatinib, erlotinib, gefitinib, gemcitabine, GSK269962A, nilotinib, rapamycin, and sorafenib. A total of 55 significantly related biological processes and 17 KEGG signaling pathways were screened. At last, we noticed that LINC00893 had a relatively lower expression in melanoma tissue and cell lines compared with adjacent tissues and epidermal melanocyte, and down-regulation of LINC00893 could promote the malignant behavior of melanoma cells in A875 and MV3. In these two melanoma cell lines, down-regulation of m6A-related molecules like YTHDF3 and METTL3 could promote the expression of LINC00893., Conclusion: We made an analysis of m6A- and m5C- related lncRNAs in melanoma samples and a prediction of these lncRNAs' role in prognosis, tumor microenvironment, immune infiltration, and clinicopathological features. We also found that LINC00893, which is potentially regulated by m6A modification, could serve as a tumor-suppressor in melanoma and play an inhibitory role in melanoma metastasis., (© 2024 The Author(s). Skin Research and Technology published by John Wiley & Sons Ltd.)

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., (© The Author(s) 2024. Published by Oxford University Press on behalf of Nucleic Acids Research.)

The incidence of papillary thyroid carcinoma (PTC) has increased significantly in recent years, and for patients with metastatic and recurrent PTC, the options for treatment currently available are insufficient. To date, the exact molecular mechanism underlying PTC is still not fully understood. 5-Methylcytosine (m5C) RNA methylation is associated with the prognosis of a variety of tumors. However, the molecular mechanisms and biomarkers associated with m5C in the diagnosis, treatment, and prognosis of this disease have not been fully elucidated. Ten m5C regulators with significantly different expression levels were included in this study. Immune infiltration analysis revealed significant negative correlations between most of these regulators and regulatory T cells. TRDMT1, NSUN5, and NSUN6 had high weights and strong correlations in the protein-protein interaction network. Using gene ontology, Kyoto Encyclopedia of Genes and Genomes, and gene set enrichment analysis, 1489 differentially expressed genes were screened from The Cancer Genome Atlas messenger RNA matrix, indicating that these differentially expressed genes were significantly enriched in various pathways and functions related to cancers. Four m5C regulators, NSUN2, NSUN4, NSUN6, and DNMT3B, were screened as prognostic markers by least absolute shrinkage and selection operator regression analysis, and NSUN2 and NSUN6 were identified as risk factors for poor prognosis. We found that the prognostic prediction model constructed using the m5C regulators NSUN2, NSUN4, NSUN6, and DNMT3B showed good prognostic prediction ability and diagnostic ability. This model was applied to predict the survival probability of patients with PTC, the prediction ability of 5-year survival was the best. The multi-factor prognostic prediction model combined with the tumor node metastasis stage and risk score grouping showed better prognostic predictive power., Competing Interests: The authors have no conflicts of interest to disclose., (Copyright © 2024 the Author(s). Published by Wolters Kluwer Health, Inc.)

Coordination of neuronal differentiation with expansion of the neuroepithelial/neural progenitor cell (NEPC/NPC) pool is essential in early brain development. Our in vitro and in vivo studies identify independent and opposing roles for two neural-specific and differentially expressed non-coding RNAs derived from the same locus: the evolutionarily conserved lncRNA Rncr3 and the embedded microRNA miR124a-1. Rncr3 regulates NEPC/NPC proliferation and controls the biogenesis of miR124a, which determines neuronal differentiation. Rncr3 conserved exons 2/3 are cytosine methylated and bound by methyl-CpG binding protein MeCP2, which restricts expression of miR124a embedded in exon 4 to prevent premature neuronal differentiation, and to orchestrate proper brain growth. MeCP2 directly binds cytosine-methylated Rncr3 through previously unrecognized lysine residues and suppresses miR124a processing by recruiting PTBP1 to block access of DROSHA-DGCR8. Thus, miRNA processing is controlled by lncRNA m 5 C methylation along with the defined m 5 C epitranscriptomic RNA reader protein MeCP2 to coordinate brain development., (© 2024. The Author(s).)

Despite extensive research on 5-methylcytosine (5mC) in relation to smoking, there has been limited exploration into the interaction between smoking and 5-hydroxymethylcytosine (5hmC). In this study, total DNA methylation (5mC+5hmC), true DNA methylation (5mC) and hydroxymethylation (5hmC) levels were profiled utilizing conventional bisulphite (BS) and oxidative bisulphite (oxBS) treatment, measured with the Illumina Infinium Methylation EPIC BeadChip. An epigenome-wide association study (EWAS) of 5mC+5hmC methylation revealed a total of 38,575 differentially methylated positions (DMPs) and 2023 differentially methylated regions (DMRs) associated with current smoking, along with 82 DMPs and 76 DMRs associated with former smoking (FDR-adjusted p < 0.05). Additionally, a focused examination of 5mC identified 33 DMPs linked to current smoking and 1 DMP associated with former smoking (FDR-adjusted p < 0.05). In the 5hmC category, eight DMPs related to current smoking and two DMPs tied to former smoking were identified, each meeting a suggestive threshold ( p < 1 × 10 -5 ). The substantial number of recognized DMPs, including 5mC+5hmC (7069/38,575, 2/82), 5mC (0/33, 1/1), and 5hmC (2/8, 0/2), have not been previously reported. Our findings corroborated previously established methylation positions and revealed novel candidates linked to tobacco smoking. Moreover, the identification of hydroxymethylated CpG sites with suggestive links provides avenues for future research.

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., (© 2024. The Author(s).)

Nanopore sequencing platforms combined with supervised machine learning (ML) have been effective at detecting base modifications in DNA such as 5-methylcytosine (5mC) and N6-methyladenine (6mA). These ML-based nanopore callers have typically been trained on data that span all modifications on all possible DNA [Formula: see text]-mer backgrounds-a complete training dataset. However, as nanopore technology is pushed to more and more epigenetic modifications, such complete training data will not be feasible to obtain. Nanopore calling has historically been performed with hidden Markov models (HMMs) that cannot make successful calls for [Formula: see text]-mer contexts not seen during training because of their independent emission distributions. However, deep neural networks (DNNs), which share parameters across contexts, are increasingly being used as callers, often outperforming their HMM cousins. It stands to reason that a DNN approach should be able to better generalize to unseen [Formula: see text]-mer contexts. Indeed, herein we demonstrate that a common DNN approach (DeepSignal) outperforms a common HMM approach (Nanopolish) in the incomplete data setting. Furthermore, we propose a novel hybrid HMM-DNN approach, amortized-HMM, that outperforms both the pure HMM and DNN approaches on 5mC calling when the training data are incomplete. This type of approach is expected to be useful for calling other base modifications such as 5-hydroxymethylcytosine and for the simultaneous calling of different modifications, settings in which complete training data are not likely to be available.

Obesity exacerbates tissue degeneration and compromises the integrity and reparative potential of mesenchymal stem/stromal cells (MSCs), but the underlying mechanisms have not been sufficiently elucidated. Mitochondria modulate the viability, plasticity, proliferative capacity, and differentiation potential of MSCs. We hypothesized that alterations in the 5-hydroxymethylcytosine (5hmC) profile of mitochondria-related genes may mediate obesity-driven dysfunction of human adipose-derived MSCs. MSCs were harvested from abdominal subcutaneous fat of obese and age/sex-matched non-obese subjects (n = 5 each). The 5hmC profile and expression of nuclear-encoded mitochondrial genes were examined by hydroxymethylated DNA immunoprecipitation sequencing (h MeDIP-seq) and mRNA-seq, respectively. MSC mitochondrial structure (electron microscopy) and function, metabolomics, proliferation, and neurogenic differentiation were evaluated in vitro, before and after epigenetic modulation. hMeDIP-seq identified 99 peaks of hyper-hydroxymethylation and 150 peaks of hypo-hydroxymethylation in nuclear-encoded mitochondrial genes from Obese- versus Non-obese-MSCs. Integrated hMeDIP-seq/mRNA-seq analysis identified a select group of overlapping (altered levels of both 5hmC and mRNA) nuclear-encoded mitochondrial genes involved in ATP production, redox activity, cell proliferation, migration, fatty acid metabolism, and neuronal development. Furthermore, Obese-MSCs exhibited decreased mitochondrial matrix density, membrane potential, and levels of fatty acid metabolites, increased superoxide production, and impaired neuronal differentiation, which improved with epigenetic modulation. Obesity elicits epigenetic changes in mitochondria-related genes in human adipose-derived MSCs, accompanied by structural and functional changes in their mitochondria and impaired fatty acid metabolism and neurogenic differentiation capacity. These observations may assist in developing novel therapies to preserve the potential of MSCs for tissue repair and regeneration in obese individuals., (© 2024. The Author(s).)

Oxidative modification of 5-methylcytosine (5mC) by ten-eleven translocation (TET) DNA dioxygenases generates 5-hydroxymethylcytosine (5hmC), the most abundant form of oxidized 5mC. Existing single-cell bisulfite sequencing methods cannot resolve 5mC and 5hmC, leaving the cell-type-specific regulatory mechanisms of TET and 5hmC largely unknown. Here, we present joint single-nucleus (hydroxy)methylcytosine sequencing (Joint-snhmC-seq), a scalable and quantitative approach that simultaneously profiles 5hmC and true 5mC in single cells by harnessing differential deaminase activity of APOBEC3A toward 5mC and chemically protected 5hmC. Joint-snhmC-seq profiling of single nuclei from mouse brains reveals an unprecedented level of epigenetic heterogeneity of both 5hmC and true 5mC at single-cell resolution. We show that cell-type-specific profiles of 5hmC or true 5mC improve multimodal single-cell data integration, enable accurate identification of neuronal subtypes and uncover context-specific regulatory effects on cell-type-specific genes by TET enzymes., (© 2023. The Author(s), under exclusive licence to Springer Nature America, Inc.)

Varicose vein disease (VVD) is a common health problem worldwide. Microfibril-associated protein 5 (MFAP5) is one of the potential key players in its pathogenesis. Our previous microarray analysis revealed the cg06256735 and cg15815843 loci in the regulatory regions of the MFAP5 gene as hypomethylated in varicose veins which correlated with its up-regulation. The aim of this work was to validate preliminary microarray data, estimate the level of 5-hydroxymethylcytosine (5hmC) at these loci, and determine the methylation status of one of them in different layers of the venous wall. For this, methyl- and hydroxymethyl-sensitive restriction techniques were used followed by real-time PCR and droplet digital PCR, correspondingly, as well as bisulfite pyrosequencing of +/- oxidized DNA. Our microarray data on hypomethylation at the cg06256735 and cg15815843 loci in whole varicose vein segments were confirmed and it was also demonstrated that the level of 5hmC at these loci is increased in VVD. Specifically, among other layers of the venous wall, tunica (t.) intima is the main contributor to hypomethylation at the cg06256735 locus in varicose veins. Thus, it was shown that hypomethylation at the cg06256735 and cg15815843 loci takes place in VVD, with evidence to suggest that it happens through their active demethylation leading to up-regulation of the MFAP5 gene, and t. intima is most involved in this biochemical process., (© 2024 The Author(s).)

RNA 5-methylcytosine (m 5 C) methylation plays a crucial role in hepatocellular carcinoma (HCC). As reported, aberrant m 5 C methylation is closely associated with the progression, therapeutic efficacy, and prognosis of HCC. The innate immune system functions as the primary defense mechanism in the body against pathogenic infections and tumors since it can activate innate immune pathways through pattern recognition receptors to exert anti-infection and anti-tumor effects. Recently, m 5 C methylation has been demonstrated to affect the activation of innate immune pathways including TLR, cGAS-STING, and RIG-I pathways by modulating RNA function, unveiling new mechanisms underlying the regulation of innate immune responses by tumor cells. However, research on m 5 C methylation and its interplay with innate immune pathways is still in its infancy. Therefore, this review details the biological significance of RNA m 5 C methylation in HCC and discusses its potential regulatory relationship with TLR, cGAS-STING, and RIG-I pathways, thereby providing fresh insights into the role of RNA methylation in the innate immune mechanisms and treatment of HCC., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2024 Meng, Jiangtao, Haisong, Mei, Long and Shanfeng.)