RNA 5-methylcytosine modification and its emerging role as an epitranscriptomic mark

Mice with NOP2/sun RNA methyltransferase 5 deficiency die before reaching puberty due to fatal kidney damage

BACKGROUND

NOP2/Sun RNA methyltransferase 5 (NSUN5) is an RNA methyltransferase that has a broad distribution and plays critical roles in various biological processes. However, our knowledge of the biological functions of NSUN5 in mammals is very limited. Therefore, in this study, we investigate the role of NSUN5 in mice.

METHODS

In the present research, we built a mouse model (Nsun5-/-) using the CRISPR/Cas9 system to investigated the specific role of NSUN5.

RESULTS

We observed that Nsun5-/- mice had a reduced body weight compared to wild-type mice. Additionally, their survival rate gradually decreased to 20% after postnatal day (PD) 21. Further examination revealed the Nsun5-/- mice had multiple organ damage, with the most severe damage occurring in the kidneys. Moreover, we observed glycogen deposition and fibrosis, along with a notable shorting of the primary foot processes of glomeruli in Nsun5-/- kidneys. Furthermore, we found that the kidneys of Nsun5-/- mice showed increased expression of the apoptotic signal Caspase-3 and accumulated stronger DNA damage at PD 21.

CONCLUSIONS

In our study, we found that mice lacking NSUN5 died before puberty due to kidney fatal damage caused by DNA damage and cell apoptosis. These results suggest that NSUN5 plays a vital role in preventing the accumulation of DNA damage and cell apoptosis in the kidney.

m6A- and m5C- modified lncRNAs orchestrate the prognosis in cutaneous melanoma and m6A- modified LINC00893 regulates cutaneous melanoma cell metastasis

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.

Non‐m6A RNA modifications in haematological malignancies

Dysregulated RNA modifications, stemming from the aberrant expression and/or malfunction of RNA modification regulators operating through various pathways, play pivotal roles in driving the progression of haematological malignancies. Among RNA modifications, N6-methyladenosine (m6A) RNA modification, the most abundant internal mRNA modification, stands out as the most extensively studied modification. This prominence underscores the crucial role of the layer of epitranscriptomic regulation in controlling haematopoietic cell fate and therefore the development of haematological malignancies. Additionally, other RNA modifications (non-m6A RNA modifications) have gained increasing attention for their essential roles in haematological malignancies. Although the roles of the m6A modification machinery in haematopoietic malignancies have been well reviewed thus far, such reviews are lacking for non-m6A RNA modifications. In this review, we mainly focus on the roles and implications of non-m6A RNA modifications, including N4-acetylcytidine, pseudouridylation, 5-methylcytosine, adenosine to inosine editing, 2'-O-methylation, N1-methyladenosine and N7-methylguanosine in haematopoietic malignancies. We summarise the regulatory enzymes and cellular functions of non-m6A RNA modifications, followed by the discussions of the recent studies on the biological roles and underlying mechanisms of non-m6A RNA modifications in haematological malignancies. We also highlight the potential of therapeutically targeting dysregulated non-m6A modifiers in blood cancer.

NSUN2 mediates distinct pathways to regulate enterovirus 71 replication

Increasing evidences suggest that the methyltransferase NSUN2 catalyzes 5-methylcytosine (m5C) modifications on viral RNAs, which are essential for the replication of various viruses. Despite the function of m5C deposition is well characterized, other potential roles of NSUN2 in regulating viral replication remain largely unknown. In this study, the m5C modified residues catalyzed by NSUN2 on enterovirus 71 (EV71) RNAs were mapped. NSUN2, along with m5C modifications, played multiple roles during the EV71 life cycle. Functional m5C modified nucleotides increased the translational efficiency and stability of EV71 RNAs. Additionally, NSUN2 was found to target the viral protein VP1 for binding and promote its stability by inhibiting the ubiquitination. Furthermore, both viral replication and pathogenicity in mice were largely attenuated when functional m5C residues were mutated. Taken together, this study characterizes distinct pathways mediated by NSUN2 in regulating EV71 replication, and highlights the importance of its catalyzed m5C modifications on EV71 RNAs for the viral replication and pathogenicity.

The je ne sais quoi of 5-methylcytosine in messenger RNA

The potential presence of 5-methylcytosine as a sparse internal modification of mRNA was first raised in 1975, and a first map of the modification was also part of the epitranscriptomics "big bang" in 2012. Since then, the evidence for its presence in mRNA has firmed up, and initial insights have been gained into the molecular function and broader biological relevance of 5-methylcytosine when present in mRNA. Here, we summarize the status quo of the field, outline some of its current challenges, and suggest how to address them in future work.

5-methylcytosine-mediated upregulation of circular RNA 0102913 augments malignant properties of colorectal cancer cells through a microRNA-571/Rac family small GTPase 2 axis

Circular RNAs (circRNAs) are a class of stable non-coding RNAs that have emerged as key regulators in human diseases including cancer. This study investigates the role of circRNA_0102913 (circ_0102913) in malignant behavior of colorectal cancer (CRC) cells and the underpinning mechanisms. By analyzing CRC-related GSE197991, GSE159669, and GSE223001 datasets, we obtained circ_0102913 as an aberrantly upregulated circRNA in CRC. Increased circ_0102913 expression was detected in CRC tissues and cells. By querying multiple bioinformatics systems (circBank, Circular RNA Interactome, TargetScan, miRDIP, miRwalk, and miRDB), we identified microRNA-571 (miR-571) as a target of circ_0102913 and Rac family small GTPase 2 (RAC2) mRNA as a target of miR-571. Biotinylated-RNA pull-down and/or luciferase assays showed that circ_0102913 bound to miR-571 to restore the expression of RAC2 mRNA. Circ_0102913 silencing or miR-571 overexpression repressed proliferation, migration and invasion, and in vivo tumorigenesis abilities of CRC cells. However, the malignant properties of cells were restored by RAC2 overexpression. The increased circ_0102913 expression in CRC cells was attributed to increased 5-methylcytosine (m5C) modification levels. Silencing of NOP2/Sun RNA methyltransferase 5 reduced the m5C level and therefore reduced stability and expression of circ_0102913 expression in CRC cells. In conclusion, this study demonstrates that m5C-mediated upregulation of circ_0102913 augments malignant properties of CRC cells through a miR-571/RAC2 axis.

NSUN2-Mediated m5C Methylation Impairs Endometrial Receptivity

Impaired endometrial decidualization is the primary cause of recurrent implantation failure (RIF). RNA methylation modification, especially NSUN family mediated m5C, is crucial for various physiological events, such as maternal-to-zygotic transition, gametogenesis, embryonic development, organismal lifespan, and cell cycle. However, the regulatory mechanisms between NSUN family mediated m5C modification and RIF remain unknown. We acquired NSUN2 expression data of 15 human endometrium samples at proliferative and secretory stages from reproductive cell atlas. The overall pattern of m5C sites and genes was elucidated through m5C-BS-seq, whereas the overall m5C levels in different groups were revealed by dot blot assay. BrdU and western blotting assays were carried out to evaluate the role of NSUN2 in proliferation and autophagy. The effects of NSUN2-mediated m5C modification on embryo attachment were evaluated by an in vitro model of a confluent monolayer of Ishikawa cells cocultured with BeWo spheroids, and its downstream targets were evaluated by real-time reverse-transcription PCR and western blotting in Ishikawa cells. The molecular mechanism for NSUN2 regulating its downstream targets' expression was determined by Cut&Tag and coimmunoprecipitation assays. NSUN2 was increased in SOX9+ cells and widespread in epithelial cell type at the proliferative stage by previous single-cell RNA sequencing data. NSUN2 overexpression (NSUN2OE) in the Ishikawa cell line elevated m5C levels and promoted cell proliferation and autophagy. NSUN2OE reduced attachment efficiency of BeWo cell spheres. Overexpressed NSUN2 was found to increase STAT1 and MMP14 mRNA expressions by inducing exon skipping. NSUN2 interacted with CLDN4 through m5C modification, and NSUN2OE or NSUN2 knockdown resulted in a similar variation tendency of CLDN4. Overexpression of NSUN2 increased CLDN4 H3K9ac modification by downregulating SIRT4 expression at the protein level, leading to the upregulation of CLDN4 mRNA expression. Our results uncovered a novel intricate regulatory mechanism between NSUN2-mediated m5C and RIF and suggested a potential new therapeutic strategy for RIF.

TET3 gene rs828867 G>A polymorphism reduces neuroblastoma risk in Chinese children

Objective

Neuroblastoma (NB) is a prevalent pediatric tumor originating from primordial neural crest cells. As one of the latest epigenetics investigations focuses, RNA 5-methylcytosine (m5C) is closely related to cancer risk. TET methylcytosine dioxygenase 3 (TET3) is a demethylase for m5C modification. Whether there is an association between TET3 gene polymorphisms and neuroblastoma risk remains unclear.

Methods

We conducted an epidemiological study in 402 patients and 473 controls to evaluate the relationship between TET3 gene SNPs (rs7560668 T > C, rs828867 G > A, and rs6546891 A > G) and NB susceptibility.

Results

Our results showed that rs828867 G > A significantly reduced NB risk in Chinese children [GA vs. GG, adjusted odds ratio (OR) = 0.72, 95% confidence interval (CI) = 0.52-0.98, P=0.040; GA/AA vs. GG, adjusted OR = 0.74, 95% CI = 0.55-0.998, P=0.048]. Individuals with 2-3 risk genotypes had a significantly higher NB risk than those with 0-1 risk genotypes (adjusted OR = 1.40, 95% CI = 1.04-1.88, P=0.027). The stratified analysis showed that the rs828867 G > A associated with decreased NB risk is remarkable among children aged >18 months (adjusted OR = 0.67, 95% CI = 0.46-0.96, P=0.029) and patients at clinical III + IV stages (adjusted OR = 0.67, 95% CI = 0.45-0.98, P=0.040). Compared with the 0-1 risk genotype, the concurrence of 2-3 risk genotypes significantly increased NB risk in the following subgroups: children aged >18 months and patients at clinical III + IV stages. GTEx analysis suggested that rs828867 G > A was significantly associated with RP11-287D1.4 and POLE4 mRNA expression.

Conclusions

Overall, our results revealed that rs828867 G > A in the TET3 gene is significantly associated with predisposition to NB.

Tet-dependent 5-hydroxymethyl-Cytosine modification of mRNA regulates axon guidance genes in Drosophila

Modifications of mRNA, especially methylation of adenosine, have recently drawn much attention. The much rarer modification, 5-hydroxymethylation of cytosine (5hmC), is not well understood and is the subject of this study. Vertebrate Tet proteins are 5-methylcytosine (5mC) hydroxylases and catalyze the transition of 5mC to 5hmC in DNA. These enzymes have recently been shown to have the same function in messenger RNAs in both vertebrates and in Drosophila. The Tet gene is essential in Drosophila as Tet knock-out animals do not reach adulthood. We describe the identification of Tet-target genes in the embryo and larval brain by mapping one, Tet DNA-binding sites throughout the genome and two, the Tet-dependent 5hmrC modifications transcriptome-wide. 5hmrC modifications are distributed along the entire transcript, while Tet DNA-binding sites are preferentially located at the promoter where they overlap with histone H3K4me3 peaks. The identified mRNAs are preferentially involved in neuron and axon development and Tet knock-out led to a reduction of 5hmrC marks on specific mRNAs. Among the Tet-target genes were the robo2 receptor and its slit ligand that function in axon guidance in Drosophila and in vertebrates. Tet knock-out embryos show overlapping phenotypes with robo2 and both Robo2 and Slit protein levels were markedly reduced in Tet KO larval brains. Our results establish a role for Tet-dependent 5hmrC in facilitating the translation of modified mRNAs primarily in cells of the nervous system.

Emerging Role of Environmental Epitranscriptomics and RNA Modifications in Parkinson's Disease

Environmental risk factors and gene-environment interactions play a critical role in Parkinson's disease (PD). However, the relatively large contribution of environmental risk factors in the overwhelming majority of PD cases has been widely neglected in the field. A "PD prevention agenda" proposed in this journal laid out a set of research priorities focused on preventing PD through modification of environmental risk factors. This agenda includes a call for preclinical studies to employ new high-throughput methods for analyzing transcriptomics and epigenomics to provide a deeper understanding of the effects of exposures linked to PD. Here, we focus on epitranscriptomics as a novel area of research with the potential to add to our understanding of the interplay between genes and environmental exposures in PD. Both epigenetics and epitranscriptomics have been recognized as potential mediators of the complex relationship between genes, environment, and disease. Multiple studies have identified epigenetic alterations, such as DNA methylation, associated with PD and PD-related exposures in human studies and preclinical models. In addition, recent technological advancements have made it possible to study epitranscriptomic RNA modifications, such as RNA N6-methyladenosine (m6A), and a handful of recent studies have begun to explore epitranscriptomics in PD-relevant exposure models. Continued exploration of epitranscriptomic mechanisms in environmentally relevant PD models offers the opportunity to identify biomarkers, pre-degenerative changes that precede symptom onset, and potential mitigation strategies for disease prevention and treatment.

Identification and validation of novel prognostic signatures based on m5C methylation patterns and tumor EMT profiles in head and neck squamous cell carcinoma

The role of 5-methylcytosine (m5C) in tumor initiation and progression has been increasingly recognized. However, the precise association between the regulation of m5C and the progression, metastasis, and prognosis of head and neck squamous cell carcinoma (HNSCC) has not yet been fully explored. Data from 545 HNSCC patients obtained from The Cancer Genome Atlas (TCGA) database were analyzed. Unsupervised cluster analysis was conducted using the expression levels of m5C regulatory genes. Additionally, gene set variation analysis (GSVA), single-sample gene set enrichment analysis (ssGSEA), and Cox regression analysis were utilized. Quantitative reverse transcription polymerase chain reaction (RT-qPCR), colony formation assay, transwell experiments and western blots were performed in the HNSCC cell line UM-SCC-17B to assess the expression and functional role of one of the novel signatures, CNFN. Significant expression differences were found in m5C regulatory genes between tumor and normal tissues in HNSCC. Two distinct m5C modification patterns, characterized by substantial prognostic differences, were identified. Cluster-2, which exhibited a strong association with epithelial-mesenchymal transition (EMT), was found to be associated with a poorer prognosis. Based on the m5C clusters and EMT status, differentially expressed genes (DEGs) were identified. Using DEGs, an 8-gene signature (CAMK2N1, WNT7A, F2RL1, AREG, DEFB1, CNFN, TGFBI, and CAV1) was established to develop a prognostic model. The performance of this signature was validated in both the training and external validation datasets, demonstrating its promising efficacy. Furthermore, additional investigations using RT-qPCR on clinical specimens and experimental assays in cell lines provided compelling evidence suggesting that CNFN, one of the genes in the signature, could play a role in HNSCC progression and metastasis through the EMT pathway. This study highlighted the role of m5C in HNSCC progression and metastasis. The relationship between m5C and EMT has been elucidated for the first time. A robust prognostic model was developed for accurately predicting HNSCC patients' survival outcomes. Potential molecular mechanisms underlying these associations have been illuminated through this research.

The regulation of antiviral innate immunity through non-m6A RNA modifications

The post-transcriptional RNA modifications impact the dynamic regulation of gene expression in diverse biological and physiological processes. Host RNA modifications play an indispensable role in regulating innate immune responses against virus infection in mammals. Meanwhile, the viral RNAs can be deposited with RNA modifications to interfere with the host immune responses. The N6-methyladenosine (m6A) has boosted the recent emergence of RNA epigenetics, due to its high abundance and a transcriptome-wide widespread distribution in mammalian cells, proven to impact antiviral innate immunity. However, the other types of RNA modifications are also involved in regulating antiviral responses, and the functional roles of these non-m6A RNA modifications have not been comprehensively summarized. In this Review, we conclude the regulatory roles of 2'-O-methylation (Nm), 5-methylcytidine (m5C), adenosine-inosine editing (A-to-I editing), pseudouridine (Ψ), N1-methyladenosine (m1A), N7-methylguanosine (m7G), N6,2'-O-dimethyladenosine (m6Am), and N4-acetylcytidine (ac4C) in antiviral innate immunity. We provide a systematic introduction to the biogenesis and functions of these non-m6A RNA modifications in viral RNA, host RNA, and during virus-host interactions, emphasizing the biological functions of RNA modification regulators in antiviral responses. Furthermore, we discussed the recent research progress in the development of antiviral drugs through non-m6A RNA modifications. Collectively, this Review conveys knowledge and inspiration to researchers in multiple disciplines, highlighting the challenges and future directions in RNA epitranscriptome, immunology, and virology.

Tet-dependent 5-hydroxymethyl-Cytosine modification of mRNA regulates axon guidance genes in Drosophila

Modifications of mRNA, especially methylation of adenosine, have recently drawn much attention. The much rarer modification, 5-hydroxymethylation of cytosine (5hmC), is not well understood and is the subject of this study. Vertebrate Tet proteins are 5-methylcytosine (5mC) hydroxylases and catalyze the transition of 5mC to 5hmC in DNA. These enzymes have recently been shown to have the same function in messenger RNAs in both vertebrates and in Drosophila. The Tet gene is essential in Drosophila as Tet knock-out animals do not reach adulthood. We describe the identification of Tet-target genes in the embryo and larval brain by mapping one, Tet DNA-binding sites throughout the genome and two, the Tet-dependent 5hmrC modifications transcriptome-wide. 5hmrC modifications are distributed along the entire transcript, while Tet DNA-binding sites are preferentially located at the promoter where they overlap with histone H3K4me3 peaks. The identified mRNAs are preferentially involved in neuron and axon development and Tet knock-out led to a reduction of 5hmrC marks on specific mRNAs. Among the Tet-target genes were the robo2 receptor and its slit ligand that function in axon guidance in Drosophila and in vertebrates. Tet knock-out embryos show overlapping phenotypes with robo2 and both Robo2 and Slit protein levels were markedly reduced in Tet KO larval brains. Our results establish a role for Tet-dependent 5hmrC in facilitating the translation of modified mRNAs primarily in cells of the nervous system.

Vital roles of m5C RNA modification in cancer and immune cell biology

RNA modification plays an important role in epigenetics at the posttranscriptional level, and 5-methylcytosine (m5C) has attracted increasing attention in recent years due to the improvement in RNA m5C site detection methods. By influencing transcription, transportation and translation, m5C modification of mRNA, tRNA, rRNA, lncRNA and other RNAs has been proven to affect gene expression and metabolism and is associated with a wide range of diseases, including malignant cancers. RNA m5C modifications also substantially impact the tumor microenvironment (TME) by targeting different groups of immune cells, including B cells, T cells, macrophages, granulocytes, NK cells, dendritic cells and mast cells. Alterations in immune cell expression, infiltration and activation are highly linked to tumor malignancy and patient prognosis. This review provides a novel and holistic examination of m5C-mediated cancer development by examining the exact mechanisms underlying the oncogenicity of m5C RNA modification and summarizing the biological effects of m5C RNA modification on tumor cells as well as immune cells. Understanding methylation-related tumorigenesis can provide useful insights for the diagnosis as well as the treatment of cancer.

Epigenetic inactivation of the 5-methylcytosine RNA methyltransferase NSUN7 is associated with clinical outcome and therapeutic vulnerability in liver cancer

BACKGROUND

RNA modifications are important regulators of transcript activity and an increasingly emerging body of data suggests that the epitranscriptome and its associated enzymes are altered in human tumors.

METHODS

Combining data mining and conventional experimental procedures, NSUN7 methylation and expression status was assessed in liver cancer cell lines and primary tumors. Loss-of-function and transfection-mediated recovery experiments coupled with RNA bisulfite sequencing and proteomics determined the activity of NSUN7 in downstream targets and drug sensitivity.

RESULTS

In this study, the initial screening for genetic and epigenetic defects of 5-methylcytosine RNA methyltransferases in transformed cell lines, identified that the NOL1/NOP2/Sun domain family member 7 (NSUN7) undergoes promoter CpG island hypermethylation-associated with transcriptional silencing in a cancer-specific manner. NSUN7 epigenetic inactivation was common in liver malignant cells and we coupled bisulfite conversion of cellular RNA with next-generation sequencing (bsRNA-seq) to find the RNA targets of this poorly characterized putative RNA methyltransferase. Using knock-out and restoration-of-function models, we observed that the mRNA of the coiled-coil domain containing 9B (CCDC9B) gene required NSUN7-mediated methylation for transcript stability. Most importantly, proteomic analyses determined that CCDC9B loss impaired protein levels of its partner, the MYC-regulator Influenza Virus NS1A Binding Protein (IVNS1ABP), creating sensitivity to bromodomain inhibitors in liver cancer cells exhibiting NSUN7 epigenetic silencing. The DNA methylation-associated loss of NSUN7 was also observed in primary liver tumors where it was associated with poor overall survival. Interestingly, NSUN7 unmethylated status was enriched in the immune active subclass of liver tumors.

CONCLUSION

The 5-methylcytosine RNA methyltransferase NSUN7 undergoes epigenetic inactivation in liver cancer that prevents correct mRNA methylation. Furthermore, NSUN7 DNA methylation-associated silencing is associated with clinical outcome and distinct therapeutic vulnerability.

Epitranscriptomics in the development, functions, and disorders of cancer stem cells

Biomolecular modifications play an important role in the development of life, and previous studies have investigated the role of DNA and proteins. In the last decade, with the development of sequencing technology, the veil of epitranscriptomics has been gradually lifted. Transcriptomics focuses on RNA modifications that affect gene expression at the transcriptional level. With further research, scientists have found that changes in RNA modification proteins are closely linked to cancer tumorigenesis, progression, metastasis, and drug resistance. Cancer stem cells (CSCs) are considered powerful drivers of tumorigenesis and key factors for therapeutic resistance. In this article, we focus on describing RNA modifications associated with CSCs and summarize the associated research progress. The aim of this review is to identify new directions for cancer diagnosis and targeted therapy.

NSUN2 promotes osteosarcoma progression by enhancing the stability of FABP5 mRNA via m5C methylation

5-methylcytosine (m⁵C) modification, which is mainly induced by the RNA methyltransferase NSUN2 (NOP2/Sun domain family, member 2), is an important chemical posttranscriptional modification in mRNA and has been proven to play important roles in the progression of many cancers. However, the functions and underlying molecular mechanisms of NSUN2-mediated m⁵C in osteosarcoma (OS) remain unclear. In this study, we found NSUN2 was highly expressed in OS tissues and cells. We also discovered that higher expression of NSUN2 predicted poorer prognosis of OS patients. Our study showed that NSUN2 could promote the progression of OS cells. Moreover, we employed RNA sequencing, RNA immunoprecipitation (RIP), and methylated RIP to screen and validate the candidate targets of NSUN2 and identified FABP5 as the target. We observed that NSUN2 stabilized FABP5 mRNA by inducing m⁵C modification and further promoted fatty acid metabolism in OS cells. Moreover, both knocking down the expression of FABP5 and adding fatty acid oxidation inhibitor could counterbalance the promoting effect of NSUN2 on the progression of OS. Our study confirms that NSUN2 can up-regulate the expression of FABP5 by improving the stability of FABP5 mRNA via m⁵C, so as to promote fatty acid metabolism in OS cells, and finally plays the role in promoting the progression of OS. Our findings suggest that NSUN2 is a promising prognostic marker for OS patients and may serve as a potential therapeutic target for OS treatment. A schematic illustration was proposed to summarize our findings.

Tet-dependent 5-hydroxymethyl-Cytosine modification of mRNA regulates the axon guidance genes robo2 and slit in Drosophila

Modifications of mRNA, especially methylation of adenosine, have recently drawn much attention. The much rarer modification, 5-hydroxymethylation of cytosine (5hmC), is not well understood and is the subject of this study. Vertebrate Tet proteins are 5-methylcytosine (5mC) hydroxylases enzymes catalyzing the transition of 5mC to 5hmC in DNA and have recently been shown to have the same function in messenger RNAs in both vertebrates and in Drosophila. The Tet gene is essential in Drosophila because Tet knock-out animals do not reach adulthood. We describe the identification of Tet-target genes in the embryo and larval brain by determining Tet DNA-binding sites throughout the genome and by mapping the Tet-dependent 5hmrC modifications transcriptome-wide. 5hmrC-modified sites can be found along the entire transcript and are preferentially located at the promoter where they overlap with histone H3K4me3 peaks. The identified mRNAs are frequently involved in neuron and axon development and Tet knock-out led to a reduction of 5hmrC marks on specific mRNAs. Among the Tet-target genes were the robo2 receptor and its slit ligand that function in axon guidance in Drosophila and in vertebrates. Tet knock-out embryos show overlapping phenotypes with robo2 and are sensitized to reduced levels of slit. Both Robo2 and Slit protein levels were markedly reduced in Tet KO larval brains. Our results establish a role for Tet-dependent 5hmrC in facilitating the translation of modified mRNAs, primarily in developing nerve cells.

Genetic variants in m5C modification core genes are associated with the risk of Chinese pediatric acute lymphoblastic leukemia: A five-center case-control study

Objective

To explore the functions of the polymorphisms in 5-methylcytosine (m5C) modification-related coding genes on the susceptibility of pediatric acute lymphoblastic leukemia (ALL).

Methods

Case-control study and multinomial logistic regression analysis were performed to construct models to evaluate the susceptibility of pediatric ALL. The relationship between five functional SNPs in m5C modification-coding genes and pediatric ALL risk was analyzed. Genotyping of 808 cases and 1,340 healthy samples from South China was identified using a TaqMan assay; odds ratios (ORs) and 95% confidence intervals (CIs) were calculated to estimate the relationship between the five selected SNPs and pediatric ALL susceptibility.

Results

Among the five analyzed SNPs, NOL1 rs3764909 and NSUN4 rs10252 variants significantly increased the susceptibility of pediatric ALL, while NSUN3 rs7653521, NSUN5 rs1880948, and NSUN6 rs3740102 variants were not associated with the risk of ALL. Stratification analyses demonstrated that NOL1 rs3764909 C>A exhibited a significant association with increased pediatric ALL risk in subgroups of common B ALL, pre-B ALL, T-cell ALL, low and middle risk, other gene fusion types, non-gene fusion, hypodiploid, normal diploid, primitive lymphocytes in marrow < 5% on week 12, and minimal residual disease (MRD) <0.01% on week 12 after induced therapy; NSUN4 rs10252 G>A was related to increased risk of ALL children in subgroups of age ≥ 120 months, normal white blood cell (WBC) number, middle risk, non-gene fusion, MRD ≥ 0.01 on days 15-19, and primitive lymphocytes in marrow < 5% on day 33 after induced therapy. Compared with the reference haplotype CAGTA, children who harbored haplotypes CCGTG and ACATA were remarkably related to increased ALL susceptibility. rs3764909 and rs10252 varieties of alleles were not associated with MRD levels after the selected chemotherapeutics.

Conclusions

In conclusion, NOL1 rs3764909 and NSUN4 rs10252 variants were enhanced by pediatric ALL risk and were suggested to be potential biomarkers for pediatric ALL.