Dynamic transcriptomic m5 C and its regulatory role in RNA processing

The covalent nucleotide modifications within plant mRNAs: What we know, how we find them, and what should be done in the future

Although covalent nucleotide modifications were first identified on the bases of transfer RNAs (tRNAs) and ribosomal RNAs (rRNAs), a number of these epitranscriptome marks have also been found to occur on the bases of messenger RNAs (mRNAs). These covalent mRNA features have been demonstrated to have various and significant effects on the processing (e.g. splicing, polyadenylation, etc.) and functionality (e.g. translation, transport, etc.) of these protein-encoding molecules. Here, we focus our attention on the current understanding of the collection of covalent nucleotide modifications known to occur on mRNAs in plants, how they are detected and studied, and the most outstanding future questions of each of these important epitranscriptomic regulatory signals.

Modifications of noncoding RNAs in cancer and their therapeutic implications

In the last 50 years, over 150 various chemical modifications on RNA molecules, including mRNAs, rRNAs, tRNAs, and other noncoding RNAs (ncRNAs), have been identified and characterized. These RNA modifications regulate RNA biogenesis and biological functions and are widely involved in various physiological processes and diseases, including cancer. In recent decades, broad interest has arisen in the epigenetic modification of ncRNAs due to the increased knowledge of the critical roles of ncRNAs in cancer. In this review, we summarize the various modifications of ncRNAs and highlight their roles in cancer initiation and progression. In particular, we discuss the potential of RNA modifications as novel biomarkers and therapeutic targets in cancer.

Aberrant m5C hypermethylation mediates intrinsic resistance to gefitinib through NSUN2/YBX1/QSOX1 axis in EGFR-mutant non-small-cell lung cancer

BACKGROUND

RNA 5-methylcytosine (m⁵C) modification plays critical roles in the pathogenesis of various tumors. However, the function and molecular mechanism of RNA m⁵C modification in tumor drug resistance remain unclear.

METHODS

The correlation between RNA m⁵C methylation, m⁵C writer NOP2/Sun RNA methyltransferase family member 2 (NSUN2) and EGFR-TKIs resistance was determined in non-small-cell lung cancer (NSCLC) cell lines and patient samples. The effects of NSUN2 on EGFR-TKIs resistance were investigated by gain- and loss-of-function assays in vitro and in vivo. RNA-sequencing (RNA-seq), RNA bisulfite sequencing (RNA-BisSeq) and m⁵C methylated RNA immunoprecipitation-qPCR (MeRIP-qPCR) were performed to identify the target gene of NSUN2 involved in EGFR-TKIs resistance. Furthermore, the regulatory mechanism of NSUN2 modulating the target gene expression was investigated by functional rescue and puromycin incorporation assays.

RESULTS

RNA m⁵C hypermethylation and NSUN2 were significantly correlated with intrinsic resistance to EGFR-TKIs. Overexpression of NSUN2 resulted in gefitinib resistance and tumor recurrence, while genetic inhibition of NSUN2 led to tumor regression and overcame intrinsic resistance to gefitinib in vitro and in vivo. Integrated RNA-seq and m⁵C-BisSeq analyses identified quiescin sulfhydryl oxidase 1 (QSOX1) as a potential target of aberrant m⁵C modification. NSUN2 methylated QSOX1 coding sequence region, leading to enhanced QSOX1 translation through m⁵C reader Y-box binding protein 1 (YBX1).

CONCLUSIONS

Our study reveals a critical function of aberrant RNA m⁵C modification via the NSUN2-YBX1-QSOX1 axis in mediating intrinsic resistance to gefitinib in EGFR-mutant NSCLC.

NSUN2-mediated m5 C RNA methylation dictates retinoblastoma progression through promoting PFAS mRNA stability and expression

BACKGROUND

The precise temporal and spatial regulation of N⁵ -methylcytosine (m⁵ C) RNA modification plays essential roles in RNA metabolism, and is necessary for the maintenance of epigenome homeostasis. Howbeit, the mechanism underlying the m⁵ C modification in carcinogenesis remains to be fully addressed.

METHODS

Global and mRNA m⁵ C levels were determined by mRNA isolation and anti-m⁵ C dot blot in both retinoblastoma (RB) cells and clinical samples. Orthotopic intraocular xenografts were established to examine the oncogenic behaviours of RB. Genome-wide multiomics analyses were performed to identify the functional target of NSUN2, including proteomic analysis, transcriptome screening and m⁵ C-methylated RNA immunoprecipitation sequencing (m⁵ C-meRIP-seq). Organoid-based single-cell analysis and gene-correlation analysis were performed to verify the NSUN2/ALYREF/m⁵ C-PFAS oncogenic cascade.

RESULTS

Herein, we report that NSUN2-mediated m⁵ C RNA methylation fuels purine biosynthesis during the oncogenic progression of RB. First, we discovered that global and mRNA m⁵ C levels were significantly enriched in RBs compared to normal retinas. In addition, tumour-specific NSUN2 expression was noted in RB samples and cell lines. Therapeutically, targeted correction of NSUN2 exhibited efficient therapeutic efficacy in RB both in vitro and in vivo. Through multiomics analyses, we subsequently identified phosphoribosylformylglycinamidine synthase (PFAS), a vital enzyme in purine biosynthesis, as a downstream candidate target of NSUN2. The reintroduction of PFAS largely reversed the inhibitory phenotypes in NSUN2-deficient RB cells, indicating that PFAS was a functional downstream target of NSUN2. Mechanistically, we found that the m⁵ C reader protein ALYREF was responsible for the recognition of the m⁵ C modification of PFAS, increasing its expression by enhancing its RNA stability.

CONCLUSIONS

Conclusively, we initially demonstrated that NSUN2 is necessary for oncogenic gene activation in RB, expanding the current understanding of dynamic m⁵ C function during tumour progression. As the NSUN2/ALYREF/m⁵ C-PFAS oncogenic cascade is an important RB trigger, our study suggests that a targeted m⁵ C reprogramming therapeutic strategy may be a novel and efficient anti-tumour therapy approach.

5-methyladenosine regulators play a crucial role in development of chronic hypersensitivity pneumonitis and idiopathic pulmonary fibrosis

5-methyladenosine (m5C) modification regulates gene expression and biological functions in oncologic areas. However, the effect of m5C modification in chronic hypersensitivity pneumonitis (CHP) and idiopathic pulmonary fibrosis (IPF) remains unknown. Expression data for 12 significant m5C regulators were obtained from the interstitial lung disease dataset. Five candidate m5C regulators, namely tet methylcytosine dioxygenase 2, NOP2/Sun RNA methyltransferase 5, Y-box binding protein 1, tRNA aspartic acid methyltransferase 1, and NOP2/Sun RNA methyltransferase 3 were screened using random forest and nomogram models to predict risks of pulmonary fibrosis. Next, we applied the consensus clustering method to stratify the samples with different m5C patterns into two groups (cluster A and B). Finally, we calculated immune cell infiltration scores via single-sample gene set enrichment analysis, then compared immune cell infiltration, related functions as well as the expression of programmed cell death 1 (PD-1, PDCD1) and programmed death protein ligand-1 (PD-L1, CD274) between the two clusters. Principal component analysis of m5C-related scores across the 288 samples revealed that cluster A had higher immune-related expression than B. Notably, T helper cell (Th) 2 type cytokines and Th1 signatures were more abundant in clusters A and B, respectively. Our results suggest that m5C is associated with and plays a crucial role in development of pulmonary fibrosis. These m5C patterns could be potential biomarkers for identification of CHP and IPF, and guide future development of immunotherapy or other new drugs strategies for pulmonary fibrosis.

The lncRNA epigenetics: The significance of m6A and m5C lncRNA modifications in cancer

Most of our transcribed RNAs are represented by non-coding sequences. Long non-coding RNAs (lncRNAs) are transcripts with no or very limited protein coding ability and a length >200nt. They can be epigenetically modified. N6-methyladenosine (m6A), N1-methyladenosine (m1A), 5-methylcytosine (m5C), 7-methylguanosine (m7G) and 2’-O-methylation (Nm) are some of the lncRNAs epigenetic modifications. The epigenetic modifications of RNA are controlled by three classes of enzymes, each playing a role in a specific phase of the modification. These enzymes are defined as “writers”, “readers” and “erasers”. m6A and m5C are the most studied epigenetic modifications in RNA. These modifications alter the structure and properties, thus modulating the functions and interactions of lncRNAs. The aberrant expression of several lncRNAs is linked to the development of a variety of cancers and the epigenetic signatures of m6A- or m5C-related lncRNAs are increasingly recognized as potential biomarkers of prognosis, predictors of disease stage and overall survival. In the present manuscript, the most up to date literature is reviewed with the focus on m6A and m5C modifications of lncRNAs and their significance in cancer.

Current and potential roles of RNA modification-mediated autophagy dysregulation in cancer

Autophagy, a cellular lysosomal degradation and survival pathway, supports nutrient recycling and adaptation to metabolic stress and participates in various stages of tumor development, including tumorigenesis, metastasis, and malignant state maintenance. Among the various factors contributing to the dysregulation of autophagy in cancer, RNA modification can regulate autophagy by directly affecting the expression of core autophagy proteins. We propose that autophagy disorder mediated by RNA modification is an important mechanism for cancer development. Therefore, this review mainly discusses the role of RNA modification-mediated autophagy regulation in tumorigenesis. We summarize the molecular basis of autophagy and the core proteins and complexes at different stages of autophagy, especially those involved in cancer development. Moreover, we describe the crosstalk of RNA modification and autophagy and review the recent advances and potential role of the RNA modification/autophagy axis in the development of multiple cancers. Furthermore, the dual role of the RNA modification/autophagy axis in cancer drug resistance is discussed. A comprehensive understanding and extensive exploration of the molecular crosstalk of RNA modifications with autophagy will provide important insights into tumor pathophysiology and provide more options for cancer therapeutic strategies.

Epi-miRNAs: Modern mediators of methylation status in human cancers

Methylation of the fundamental macromolecules, DNA/RNA, and proteins, is remarkably abundant, evolutionarily conserved, and functionally significant in cellular homeostasis and normal tissue/organism development. Disrupted methylation imprinting is strongly linked to loss of the physiological equilibrium and numerous human pathologies, and most importantly to carcinogenesis, tumor heterogeneity, and cancer progression. Mounting recent evidence has documented the active implication of miRNAs in the orchestration of the multicomponent cellular methylation machineries and the deregulation of methylation profile in the epigenetic, epitranscriptomic, and epiproteomic levels during cancer onset and progression. The elucidation of such regulatory networks between the miRNome and the cellular methylation machineries has led to the emergence of a novel subclass of miRNAs, namely "epi-miRNAs" or "epi-miRs." Herein, we have summarized the existing knowledge on the functional role of epi-miRs in the methylation dynamic landscape of human cancers and their clinical utility in modern cancer diagnostics and tailored therapeutics. This article is categorized under: RNA in Disease and Development > RNA in Disease.

m5C-dependent cross-regulation between nuclear reader ALYREF and writer NSUN2 promotes urothelial bladder cancer malignancy through facilitating RABL6/TK1 mRNAs splicing and stabilization

The significance of 5-methylcytosine (m⁵C) methylation in human malignancies has become an increasing focus of investigation. Here, we show that m⁵C regulators including writers, readers and erasers, are predominantly upregulated in urothelial carcinoma of the bladder (UCB) derived from Sun Yat-sen University Cancer Center and The Cancer Genome Atlas cohort. In addition, NOP2/Sun RNA methyltransferase family member 2 (NSUN2) as a methyltransferase and Aly/REF export factor (ALYREF) as a nuclear m⁵C reader, are frequently coexpressed in UCB. By applying patient-derived organoids model and orthotopic xenograft mice model, we demonstrate that ALYREF enhances proliferation and invasion of UCB cells in an m⁵C-dependent manner. Integration of tanscriptome-wide RNA bisulphite sequencing (BisSeq), RNA-sequencing (RNA-seq) and RNA Immunoprecipitation (RIP)-seq analysis revealed that ALYREF specifically binds to hypermethylated m⁵C site in RAB, member RAS oncogene family like 6 (RABL6) and thymidine kinase 1 (TK1) mRNA via its K171 domain. ALYREF controls UCB malignancies through promoting hypermethylated RABL6 and TK1 mRNA for splicing and stabilization. Moreover, ALYREF recognizes hypermethylated m⁵C site of NSUN2, resulting in NSUN2 upregulation in UCB. Clinically, the patients with high coexpression of ALYREF/RABL6/TK1 axis had the poorest overall survival. Our study unveils an m⁵C dependent cross-regulation between nuclear reader ALYREF and m⁵C writer NSUN2 in activation of hypermethylated m⁵C oncogenic RNA through promoting splicing and maintaining stabilization, consequently leading to tumor progression, which provides profound insights into therapeutic strategy for UCB.

Understanding the Epitranscriptome for Avant-Garde Brain Tumour Diagnostics

RNA modifications are diverse, dynamic, and reversible transcript alterations rapidly gaining attention due to their newly defined RNA regulatory roles in cellular pathways and pathogenic mechanisms. The exciting emerging field of 'epitranscriptomics' is predominantly centred on studying the most abundant mRNA modification, N6-methyladenine (m⁶A). The m⁶A mark, similar to many other RNA modifications, is strictly regulated by so-called 'writer', 'reader', and 'eraser' protein species. The abundance of genes coding for the expression of these regulator proteins and m⁶A levels shows great potential as diagnostic and predictive tools across several cancer fields. This review explores our current understanding of RNA modifications in glioma biology and the potential of epitranscriptomics to develop new diagnostic and predictive classification tools that can stratify these highly complex and heterogeneous brain tumours.

Transcriptional and post-transcriptional control of autophagy and adipogenesis by YBX1

Obesity is strongly associated with metabolic diseases, which have become a global health problem. Exploring the underlying mechanism of adipogenesis is crucial for the treatment of excess white fat. Oncogene YBX1 is a multifunctional DNA- and RNA-binding protein that regulates brown adipogenesis. However, the role of YBX1 in white adipogenesis and adipose tissue expansion remains unknown. Here, we showed that YBX1 deficiency inhibited murine and porcine adipocyte differentiation. YBX1 positively regulated adipogenesis through promoting ULK1- and ULK2-mediated autophagy. Mechanistically, we identified YBX1 serves as a 5-methylcytosine (m⁵C)-binding protein directly targeting m⁵C-containing Ulk1 mRNA by using RNA immunoprecipitation. RNA decay assay further proved that YBX1 upregulated ULK1 expression though stabilizing its mRNA. Meanwhile, YBX1 promoted Ulk2 transcription and expression as a transcription factor, thereby enhancing autophagy and adipogenesis. Importantly, YBX1 overexpression in white fat enhanced ULK1/ULK2-mediated autophagy and promoted adipose tissue expansion in mice. Collectively, these findings unveil the post-transcriptional and transcriptional mechanism and functional importance of YBX1 in autophagy and adipogenesis regulation, providing an attractive molecular target for therapies of obesity and metabolic diseases.

MAEL gene contributes to bovine testicular development through the m5C-mediated splicing

Knowledge of RNA molecules regulating testicular development and spermatogenesis in bulls is essential for elite bull selection and an ideal breeding program. Herein, we performed direct RNA sequencing (DRS) to explore the functional characterization of RNA molecules produced in the testicles of 9 healthy Simmental bulls at three testicular development stages (prepuberty, puberty, and postpuberty). We identified 5,043 differentially expressed genes associated with testicular weight. These genes exhibited more alternative splicing at sexual maturity, particularly alternative 3' (A3) and 5' (A5) splice sites usage and exon skipping (SE). The expression of hub genes in testicular developmental stages was also affected by both m6A and m5C RNA modifications. We found m5C-mediated splicing events significantly (p < 0.05) increased MAEL gene expression at the isoform level, likely promoting spermatogenesis. Our findings highlight the complexity of RNA processing and expression as well as the regulation of transcripts involved in testicular development and spermatogenesis.

ALKBH1-mediated m1 A demethylation of METTL3 mRNA promotes the metastasis of colorectal cancer by downregulating SMAD7 expression

Colorectal cancer (CRC) is one of the most common malignancies, and the main cause of death from CRC is tumor metastasis. m¹ A RNA modification plays critical role in many biological processes. However, the role of m¹ A modification in CRC remains unclear. Here, we find that the m¹ A demethylase alkB homolog 1, histone H2A dioxygenase (ALKBH1) is overexpressed in CRC and is associated with metastasis and poor prognosis. Upregulation of ALKBH1 expression promotes CRC metastasis in vitro and in vivo. Mechanistically, knockdown of ALKBH1 results in a decrease in methyltransferase 3, N6-adenosine-methyltransferase complex catalytic subunit (METTL3) expression, probably due to m¹ A modification of METTL3 mRNA, followed by m⁶ A demethylation of SMAD family member 7 (SMAD7) mRNA. In addition, downregulation of SMAD7 establishes an aggressive phenotype. More importantly, the cell migration and invasion defects caused by ALKBH1 depletion or METTL3 depletion are significantly reversed by SMAD7 silencing. Considering these results collectively, we propose that ALKBH1 promotes CRC metastasis by destabilizing SMAD7 through METTL3.

Interplay of four types of RNA modification writers revealed distinct tumor microenvironment and biological characteristics in pancreatic cancer

Background

Pancreatic cancer (PC) is one of the most lethal malignancies and carries a dismal mortality and morbidity. Four types of RNA modification (namely m6A, m1A, APA and A-to-I) could be catalyzed by distinct enzymatic compounds ("writers"), mediating numerous epigenetic events in carcinogenesis and immunomodulation. We aim to investigate the interplay mechanism of these writers in immunogenomic features and molecular biological characteristics in PC.

Methods

We first accessed the specific expression pattern and transcriptional variation of 26 RNA modification writers in The Cancer Genome Atlas (TCGA) dataset. Unsupervised consensus clustering was performed to divide patients into two RNA modification clusters. Then, based on the differentially expressed genes (DEGs) among two clusters, RNA modification score (WM_Score) model was established to determine RNA modification-based subtypes and was validated in International Cancer Genome Consortium (ICGC) dataset. What's more, we manifested the unique status of WM_Score in transcriptional and post-transcriptional regulation, molecular biological characteristics, targeted therapies and immunogenomic patterns.

Results

We documented the tight-knit correlations between transcriptional expression and variation of RNA modification writers. We classified patients into two distinct RNA modification patterns (WM_Score_high and _low), The WM_Score_high subgroup was correlated with worse prognosis, Th2/Th17 cell polarization and oncogenic pathways (e.g. EMT, TGF-β, and mTORC1 signaling pathways), whereas the WM_Score_low subgroup associated with favorable survival rate and Th1 cell trend. WM_Score model also proved robust predictive power in interpreting transcriptional and post-transcriptional events. Additionally, the potential targeted compounds with related pathways for the WM_Score model were further identified.

Conclusions

Our research unfolds a novel horizon on the interplay network of four RNA modifications in PC. This WM_Score model demonstrated powerful predictive capacity in epigenetic, immunological and biological landscape, providing a theoretical basis for future clinical judgments of PC.

RNA Epigenetics in Chronic Lung Diseases

Chronic lung diseases are highly prevalent worldwide and cause significant mortality. Lung cancer is the end stage of many chronic lung diseases. RNA epigenetics can dynamically modulate gene expression and decide cell fate. Recently, studies have confirmed that RNA epigenetics plays a crucial role in the developing of chronic lung diseases. Further exploration of the underlying mechanisms of RNA epigenetics in chronic lung diseases, including lung cancer, may lead to a better understanding of the diseases and promote the development of new biomarkers and therapeutic strategies. This article reviews basic information on RNA modifications, including N⁶ methylation of adenosine (m⁶A), N¹ methylation of adenosine (m¹A), N⁷-methylguanosine (m⁷G), 5-methylcytosine (m⁵C), 2'O-methylation (2'-O-Me or Nm), pseudouridine (5-ribosyl uracil or Ψ), and adenosine to inosine RNA editing (A-to-I editing). We then show how they relate to different types of lung disease. This paper hopes to summarize the mechanisms of RNA modification in chronic lung disease and finds a new way to develop early diagnosis and treatment of chronic lung disease.

RNA-methyltransferase Nsun5 controls the maternal-to-zygotic transition by regulating maternal mRNA stability

BACKGROUND

RNA modification-induced ovarian dysgenesis appears to be necessary for ovary development. However, how m⁵ C (5-methylcytosine)-coordinating modificatory transcripts are dynamically regulated during oogenesis, and ovarian development is unknown. The purpose of this study was to determine whether NOP2/Sun RNA methyltransferase 5 (Nsun5) deletion leads to suppression of ovarian function and arrest of embryonic development. The regulation of mRNA decay and stability by m⁵ C modification is essential at multiple stages during the maternal-to-zygotic (MZT) transition.

METHODS

Mouse ovaries and oocytes with Nsun5^KO and the KGN cell line were subjected to m⁵ C identification, alternative splicing analysis and protein expression. BS-m⁵ C-seq, real-time polymerase chain reaction, Western blot, immunofluorescence and actinomycin D treatment assays were used. In particular, BS-m⁵ C-seq revealed a dynamic pattern of m⁵ C sites and genes in the ovaries between Nsun5^KO and WT mice at the 2-month and 6-month stages. Diverse bioinformatic tools were employed to identify target genes for Nsun5.

RESULTS

Here, a maternal mRNA stability study showed that deletion of the m⁵ C methyltransferase Nsun5 obstructs follicular development and ovarian function, which leads directly to inhibition of embryogenesis and embryo development. Dynamic analysis of m⁵ C revealed that the level of m⁵ C decreased in a time-dependent manner after Nsun5 knockout. Regarding the molecular mechanism, we found that Nsun5 deficiency caused a m⁵ C decline in the exon and 3'UTR regions that influenced the translation efficiency of Mitotic arrest deficient 2 like 2 (MAD2L2) and Growth differentiation factor 9 (GDF9) in the ovary. Mechanistic investigation of alternative splicing indicated that Nsun5^KO triggers aberrant events in the exon region of Brd8.

CONCLUSIONS

Nsun5 loss arrests follicular genesis and development in ovarian aging, indicating that Nsun5/m⁵ C-regulated maternal mRNA stabilization is essential for MZT transition.

RNA Modifications Meet Tumors

RNA modifications occur through the whole process of gene expression regulation, including transcription, translation, and post-translational processes. They are closely associated with gene expression, RNA stability, and cell cycle. RNA modifications in tumor cells play a vital role in tumor development and metastasis, changes in the tumor microenvironment, drug resistance in tumors, construction of tumor cell-cell "internet", etc. Several types of RNA modifications have been identified to date and have various effects on the biological characteristics of different tumors. In this review, we discussed the function of RNA modifications, including N 6-methyladenine (m6A), 5-methylcytosine (m5C), N 7-methyladenosine (m7G), N 1-methyladenosine (m1A), pseudouridine (Ψ), and adenosine-to-inosine (A-to-I), in the microenvironment and therapy of solid and liquid tumors.

The role of RNA modification in the generation of acquired drug resistance in glioma

Glioma is the most common malignant tumor in the central nervous system. The clinical treatment strategy is mainly surgery combined with concurrent temozolomide chemotherapy, but patients can develop drug resistance during treatment, which severely limits its therapeutic efficacy. Epigenetic regulation at the RNA level is plastic and adaptable, and it can induce a variety of tumor responses to drugs. The regulators of RNA modification include methyltransferases, demethylases, and methylation binding proteins; these are also considered to play an important role in the development, prognosis, and therapeutic response of gliomas, which provides a basis for finding new targets of epigenetic drugs and resetting the sensitivity of tumor cells to temozolomide. This review discusses the relationship between the development of adaptive drug resistance and RNA modification in glioma and summarizes the progress of several major RNA modification strategies in this field, especially RNA m6A modification, m5C modification, and adenosine-to-inosine editing.

Emerging Role of RNA m5C Modification in Cardiovascular Diseases

Epitranscriptomics is the emerging field of research that comprises the study of epigenetics changes in RNAs. Progressing development in the field of epigenetics has helped to manage and comprehend human diseases. RNA methylation regulates all aspects of RNA functions, which are involved in the pathogenesis of human diseases. Interestingly, RNA m5C methylation is significantly linked to various types of human disease, including cardiovascular diseases (CVD). The m5C methylation is controlled by m5C regulatory proteins, which act as methyltransferase, demethyltransferase, and RNA-binding protein. Dysregulated expression in m5C regulatory proteins is significantly associated with cardiovascular disease, and these regulatory proteins have crucial roles in biological and cellular functions. This review is mainly focused on the role of RNA m5C modification in CVD and mitochondrial dysfunction. Thus, m5C will contribute to discovering the new diagnostic marker and therapeutic target for CVD.

Detection technologies for RNA modifications

To date, more than 170 chemical modifications have been characterized in RNA, providing a new layer of gene expression regulation termed the 'epitranscriptome'. RNA modification detection methods and tools advance the functional studies of the epitranscriptome. According to the detection throughput and principles, existing RNA modification detection technologies can be categorized into four classes, including quantification methods, locus-specific detection methods, next-generation sequencing-based detection technologies and nanopore direct RNA sequencing-based technologies. In this review, we summarize the current knowledge about these RNA modification detection technologies and discuss the challenges for the existing detection tools, providing information for a comprehensive understanding of the epitranscriptome.

New insights into the epitranscriptomic control of pluripotent stem cell fate

Each cell in the human body has a distinguishable fate. Pluripotent stem cells are challenged with a myriad of lineage differentiation options. Defects are more likely to be fatal to stem cells than to somatic cells due to the broad impact of the former on early development. Hence, a detailed understanding of the mechanisms that determine the fate of stem cells is needed. The mechanisms by which human pluripotent stem cells, although not fully equipped with complex chromatin structures or epigenetic regulatory mechanisms, accurately control gene expression and are important to the stem cell field. In this review, we examine the events driving pluripotent stem cell fate and the underlying changes in gene expression during early development. In addition, we highlight the role played by the epitranscriptome in the regulation of gene expression that is necessary for each fate-related event.

The role of RNA modification in hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is a highly mortal type of primary liver cancer. Abnormal epigenetic modifications are present in HCC, and RNA modification is dynamic and reversible and is a key post-transcriptional regulator. With the in-depth study of post-transcriptional modifications, RNA modifications are aberrantly expressed in human cancers. Moreover, the regulators of RNA modifications can be used as potential targets for cancer therapy. In RNA modifications, N6-methyladenosine (m6A), N7-methylguanosine (m7G), and 5-methylcytosine (m5C) and their regulators have important regulatory roles in HCC progression and represent potential novel biomarkers for the confirmation of diagnosis and treatment of HCC. This review focuses on RNA modifications in HCC and the roles and mechanisms of m6A, m7G, m5C, N1-methyladenosine (m1A), N3-methylcytosine (m3C), and pseudouridine (ψ) on its development and maintenance. The potential therapeutic strategies of RNA modifications are elaborated for HCC.

An m5C methylation regulator-associated signature predicts prognosis and therapy response in pancreatic cancer

Pancreatic ductal adenocarcinoma (PDAC) is the most aggressive digestive malignancy due to frequent late-stage diagnosis, rapid progression and resistance to therapy. With increasing PDAC incidence worldwide, there is an urgent need for new prognostic biomarkers and therapy targets. Recently, RNA methylation has emerged as a new tumorigenic mechanism in different cancers. 5-methylcytosine (m5C) is one of the most frequent RNA modifications and occurs on a variety of RNA species including mRNA, thereby regulating gene expression. Here we investigated the prognostic role of m5C-regulator-associated transcriptional signature in PDAC. We evaluated m5C-regulator status and expression in 239 PDAC samples from publicly available datasets. We used unsupervised consensus clustering analyses to classify PDACs based on m5C-regulator expression. From the resulting signature of differentially expressed genes (DEGs), we selected prognosis-relevant DEGs to stratify patients and build a scoring signature (m5C-score) through LASSO and multivariate Cox regression analyses. The m5C-score represented a highly significant independent prognostic marker. A high m5C-score correlated with poor prognosis in different PDAC cohorts, and was associated with the squamous/basal subtype as well as activated cancer-related pathways including Ras, MAPK and PI3K pathways. Furthermore, the m5C-score correlated with sensitivity to pathway-specific inhibitors of PARP, EGFR, AKT, HER2 and mTOR. Tumors with high m5C-score were characterized by overall immune exclusion, low CD8⁺ T-cell infiltration, and higher PD-L1 expression. Overall, the m5C-score represented a robust predictor of prognosis and therapy response in PDAC, which was associated with unfavorable molecular subtypes and immune microenvironment.

Post-Transcriptional Modifications of RNA as Regulators of Apoptosis in Glioblastoma

This review is devoted to changes in the post-transcriptional maturation of RNA in human glioblastoma cells, which leads to disruption of the normal course of apoptosis in them. The review thoroughly highlights the latest information on both post-transcriptional modifications of certain regulatory RNAs, associated with the process of apoptosis, presents data on the features of apoptosis in glioblastoma cells, and shows the relationship between regulatory RNAs and the apoptosis in tumor cells. In conclusion, potential target candidates are presented that are necessary for the development of new drugs for the treatment of glioblastoma.

mRNA m5C inhibits adipogenesis and promotes myogenesis by respectively facilitating YBX2 and SMO mRNA export in ALYREF-m5C manner

Although 5-methylcytosine (m5C) has been identified as a novel and abundant mRNA modification and associated with energy metabolism, its regulation function in adipose tissue and skeletal muscle is still limited. This study aimed at investigating the effect of mRNA m5C on adipogenesis and myogenesis using Jinhua pigs (J), Yorkshire pigs (Y) and their hybrids Yorkshire-Jinhua pigs (YJ). We found that Y grow faster than J and YJ, while fatness-related characteristics observed in Y were lower than those of J and YJ. Besides, total mRNA m5C levels and expression rates of NSUN2 were higher both in backfat layer (BL) and longissimus dorsi muscle (LDM) of Y compared to J and YJ, suggesting that higher mRNA m5C levels positively correlate with lower fat and higher muscle mass. RNA bisulfite sequencing profiling of m5C revealed tissue-specific and dynamic features in pigs. Functionally, hyper-methylated m5C-containing genes were enriched in pathways linked to impaired adipogenesis and enhanced myogenesis. In in vitro, m5C inhibited lipid accumulation and promoted myogenic differentiation. Furthermore, YBX2 and SMO were identified as m5C targets. Mechanistically, YBX2 and SMO mRNAs with m5C modification were recognized and exported into the cytoplasm from the nucleus by ALYREF, thus leading to increased YBX2 and SMO protein expression and thereby inhibiting adipogenesis and promoting myogenesis, respectively. Our work uncovered the critical role of mRNA m5C in regulating adipogenesis and myogenesis via ALYREF-m5C-YBX2 and ALYREF-m5C-SMO manners, providing a potential therapeutic target in the prevention and treatment of obesity, skeletal muscle dysfunction and metabolic disorder diseases.

The Epitranscriptome in miRNAs: Crosstalk, Detection, and Function in Cancer

The epitranscriptome encompasses all post-transcriptional modifications that occur on RNAs. These modifications can alter the function and regulation of their RNA targets, which, if dysregulated, result in various diseases and cancers. As with other RNAs, miRNAs are highly modified by epitranscriptomic modifications such as m⁶A methylation, 2′-O-methylation, m⁵C methylation, m⁷G methylation, polyuridine, and A-to-I editing. miRNAs are a class of small non-coding RNAs that regulates gene expression at the post-transcriptional level. miRNAs have gathered high clinical interest due to their role in disease, development, and cancer progression. Epitranscriptomic modifications alter the targeting, regulation, and biogenesis of miRNAs, increasing the complexity of miRNA regulation. In addition, emerging studies have revealed crosstalk between these modifications. In this review, we will summarize the epitranscriptomic modifications—focusing on those relevant to miRNAs—examine the recent crosstalk between these modifications, and give a perspective on how this crosstalk expands the complexity of miRNA biology.

The Role of RNA Modification in HIV-1 Infection

RNA plays an important role in biology, and more than 170 RNA modifications have been identified so far. Post-transcriptional modification of RNA in cells plays a crucial role in the regulation of its stability, transport, processing, and gene expression. So far, the research on RNA modification and the exact role of its enzymes is becoming more and more comprehensive. Human immunodeficiency virus 1 (HIV-1) is an RNA virus and the causative agent of acquired immunodeficiency syndrome (AIDS), which is one of the most devastating viral pandemics in history. More and more studies have shown that HIV has RNA modifications and regulation of its gene expression during infection and replication. This review focuses on several RNA modifications and their regulatory roles as well as the roles that different RNA modifications play during HIV-1 infection, in order to find new approaches for the development of anti-HIV-1 therapeutics.

Cross-Talk of Multiple Types of RNA Modification Regulators Uncovers the Tumor Microenvironment and Immune Infiltrates in Soft Tissue Sarcoma

Background

Soft-tissue sarcoma (STS) represents a rare and diverse cohort of solid tumors, and encompasses over 100 various histologic and molecular subtypes. In recent years, RNA modifications including m6A, m5C, m1A, and m7G have been demonstrated to regulate immune response and tumorigenesis. Nevertheless, the cross-talk among these RNA modification regulators and related effects upon the tumor microenvironment (TME), immune infiltrates, and immunotherapy in STS remain poorly understood.

Methods

In this study, we comprehensively investigated transcriptional and genetic alterations of 32 RNA modification regulators in STS patients from The Cancer Genome Atlas (TCGA) cohort and validated them in the Gene Expression Omnibus (GEO) cohort. Single-cell transcriptomes were introduced to identify regulators within specific cell types, with own sequencing data and RT-qPCR conducted for biological validation. Distinct regulator clusters and regulator gene subtypes were identified by using unsupervised consensus clustering analysis. We further built the regulator score model based on the prognostic regulator-related differentially expressed genes (DEGs), which could be used to quantitatively assess the risk for individual STS patients. The clinical and biological characteristics of different regulator score groups were further examined.

Results

A total of 455 patients with STS were included in this analysis. The network of 32 RNA modification regulators demonstrated significant correlations within multiple different RNA modification types. Distinct regulator clusters and regulator gene subtypes were characterized by markedly different prognoses and TME landscapes. The low regulator score group in the TCGA-SARC cohort was characterized by poor prognosis. The robustness of the scoring model was further confirmed by the external validation in GSE30929 and GSE17674. The regulator score was negatively correlated with CD4+ T cell, Th2 cell, and Treg cell recruitment and most immunotherapy-predicted pathways, and was also associated with immunotherapy efficacy.

Conclusions

Overall, our study is the first to demonstrate the cross-talk of RNA modification regulators and the potential roles in TME and immune infiltrates in STS. The individualized assessment based on the regulator score model could facilitate and optimize personalized treatment.

RNA 5-methylcytosine status is associated with DNMT2/TRDMT1 nuclear localization in osteosarcoma cell lines

•

Selected phenotypic features of three osteosarcoma (OS) cell lines were evaluated.

•

Redox disequilibrium promoted sustained AKT and ERK1/2 activation.

•

Redox imbalance modulated cell death pathways in OS cells.

•

Nuclear levels of TRDMT1 methyltransferase were associated with RNA methylation.

•

A novel marker for predicting therapy response in OS patients is proposed.

Osteosarcoma (OS) is a pediatric malignant bone tumor with unsatisfying improvements in survival rates due to limited understanding of OS biology and potentially druggable targets. The present study aims to better characterize osteosarcoma U-2 OS, SaOS-2, and MG-63 cell lines that are commonly used as in vitro models of OS. We focused on evaluating the differences in cell death pathways, redox equilibrium, the activity of proliferation-related signaling pathways, DNA damage response, telomere maintenance, DNMT2/TRDMT1-based responses and RNA 5-methylcytosine status. SaOS-2 cells were characterized by higher levels of superoxide and nitric oxide that promoted AKT and ERK1/2 activation thus modulating cell death pathways. OS cell lines also differed in the levels and localization of DNA repair regulator DNMT2/TRDMT1. SaOS-2 cells possessed the lowest levels of total, cytoplasmic and nuclear DNMT2/TRDMT1, whereas in MG-63 cells, the highest levels of nuclear DNMT2/TRDMT1 were associated with the most pronounced status of RNA 5-methylcytosine. In silico analysis revealed potential phosphorylation sites at DNMT2/TRDMT1 that may be related to the regulation of DNMT2/TRDMT1 localization. We postulate that redox homeostasis, proliferation-related pathways and DNMT2/TRDMT1-based effects can be modulated as a part of anti-osteosarcoma strategy reflecting diverse phenotypic features of OS cells.

Comprehensive Analysis of m5C Methylation Regulatory Genes and Tumor Microenvironment in Prostate Cancer

Background

5-Methylcytidine (m5C) methylation is an emerging epigenetic modification in recent years, which is associated with the development and progression of various cancers. However, the prognostic value of m5C regulatory genes and the correlation between m5C methylation and the tumor microenvironment (TME) in prostate cancer remain unknown.

Methods

In the current study, the genetic and transcriptional alterations and prognostic value of m5C regulatory genes were investigated in The Cancer Genome Atlas and Gene Expression Omnibus datasets. Then, an m5C prognostic model was established by LASSO Cox regression analysis. Gene set variation analyses (GSVA), gene set enrichment analysis (GSEA), clinical relevance, and TME analyses were conducted to explain the biological functions and quantify the TME scores between high-risk and low-risk subgroups. m5C regulatory gene clusters and m5C immune subtypes were identified using consensus unsupervised clustering analysis. The Cell-type Identification By Estimating Relative Subsets of RNA Transcripts algorithm was used to calculate the contents of immune cells.

Results

TET3 was upregulated at transcriptional levels in PCa compared with normal tissues, and a high TET3 expression was associated with poor prognosis. An m5C prognostic model consisting of 3 genes (NSUN2, TET3, and YBX1) was developed and a nomogram was constructed for improving the clinical applicability of the model. Functional analysis revealed the enrichment of pathways and the biological processes associated with RNA regulation and immune function. Significant differences were also found in the expression levels of m5C regulatory genes, TME scores, and immune cell infiltration levels between different risk subgroups. We identified two distinct m5C gene clusters and found their correlation with patient prognosis and immune cell infiltration characteristics. Naive B cells, CD8+ T cells, M1 macrophages and M2 macrophages were obtained and 2 m5C immune subtypes were identified. CTLA4, NSUN6, TET1, and TET3 were differentially expressed between immune subtypes. The expression of CTLA4 was found to be correlated with the degree of immune cell infiltration.

Conclusions

Our comprehensive analysis of m5C regulatory genes in PCa demonstrated their potential roles in the prognosis, clinical features, and TME. These findings may improve our understanding of m5C regulatory genes in the tumor biology of PCa.

A Cross-Tissue Investigation of Molecular Targets and Physiological Functions of Nsun6 Using Knockout Mice

The 5-methylcytosine (m5C) modification on an mRNA molecule is deposited by Nsun2 and its paralog Nsun6. While the physiological functions of Nsun2 have been carefully studied using gene knockout (KO) mice, the physiological functions of Nsun6 remain elusive. In this study, we generated an Nsun6-KO mouse strain, which exhibited no apparent phenotype in both the development and adult stages as compared to wild-type mice. Taking advantage of this mouse strain, we identified 80 high-confident Nsun6-dependent m5C sites by mRNA bisulfite sequencing in five different tissues and systematically analyzed the transcriptomic phenotypes of Nsun6-KO tissues by mRNA sequencing. Our data indicated that Nsun6 is not required for the homeostasis of these organs under laboratory housing conditions, but its loss may affect immune response in the spleen and oxidoreductive reaction in the liver under certain conditions. Additionally, we further investigated T-cell-dependent B cell activation in KO mice and found that Nsun6 is not essential for the germinal center B cell formation but is associated with the formation of antibody-secreting plasma cells. Finally, we found that Nsun6-mediated m5C modification does not have any evident influence on the stability of Nsun6 target mRNAs, suggesting that Nsun6-KO-induced phenotypes may be associated with other functions of the m5C modification or Nsun6 protein.

The Methylation Game: Epigenetic and Epitranscriptomic Dynamics of 5-Methylcytosine

DNA and RNA methylation dynamics have been linked to a variety of cellular processes such as development, differentiation, and the maintenance of genome integrity. The correct deposition and removal of methylated cytosine and its oxidized analogues is pivotal for cellular homeostasis, rapid responses to exogenous stimuli, and regulated gene expression. Uncoordinated expression of DNA/RNA methyltransferases and demethylase enzymes has been linked to genome instability and consequently to cancer progression. Furthermore, accumulating evidence indicates that post-transcriptional DNA/RNA modifications are important features in DNA/RNA function, regulating the timely recruitment of modification-specific reader proteins. Understanding the biological processes that lead to tumorigenesis or somatic reprogramming has attracted a lot of attention from the scientific community. This work has revealed extensive crosstalk between epigenetic and epitranscriptomic pathways, adding a new layer of complexity to our understanding of cellular programming and responses to environmental cues. One of the key modifications, m5C, has been identified as a contributor to regulation of the DNA damage response (DDR). However, the various mechanisms of dynamic m5C deposition and removal, and the role m5C plays within the cell, remains to be fully understood.

5-methylcytosine RNA methyltransferases and their potential roles in cancer

In recent years, 5-methylcytosine (m⁵C) RNA modification has emerged as a key player in regulating RNA metabolism and function through coding as well as non-coding RNAs. Accumulating evidence has shown that m⁵C modulates the stability, translation, transcription, nuclear export, and cleavage of RNAs to mediate cell proliferation, differentiation, apoptosis, stress responses, and other biological functions. In humans, m⁵C RNA modification is catalyzed by the NOL1/NOP2/sun (NSUN) family and DNA methyltransferase 2 (DNMT2). These RNA modifiers regulate the expression of multiple oncogenes such as fizzy-related-1, forkhead box protein C2, Grb associated-binding protein 2, and TEA domain transcription factor 1, facilitating the pathogenesis and progression of cancers. Furthermore, the aberrant expression of methyltransferases have been identified in various cancers and used to predict the prognosis of patients. In this review, we present a comprehensive overview of m⁵C RNA methyltransferases. We specifically highlight the potential mechanism of action of m⁵C in cancer. Finally, we discuss the prospect of m⁵C-relative studies.

Roles of m5C RNA Modification Patterns in Biochemical Recurrence and Tumor Microenvironment Characterization of Prostate Adenocarcinoma

Background

Prostate cancer is the second most common cancer with a high risk of biochemical recurrence (BCR) among men. Recently, 5-methylcytosine (m5C) modification has attracted more attention as a new layer of RNA post-transcriptional regulation. Hence, we aimed at investigating the potential roles of m5C modification regulators in the BCR of prostate adenocarcinoma (PRAD).

Methods

CNV data, mutation annotation data, mRNA expression profiles, and clinical data were downloaded from TCGA and GEO databases. Kaplan-Meier curves analysis, log-rank test, univariate and multivariate Cox regression, and time-dependent ROC curves analysis were performed to evaluate the prognostic factors. Principal components analysis (PCA) was applied to validate the distinction between subgroups. Gene set variation analysis (GSVA) was used to investigate the underlying pathways associated with m5C modification patterns. Single sample gene set enrichment analysis (ssGSEA) was utilized to assess the infiltration of distinct immune cells. Tumor Immune Dysfunction and Exclusion (TIDE) prediction was carried out to assess the potential response to immune checkpoint blockade (ICB) therapy. The m5C modification signature was constructed via LASSO Cox's proportional hazards regression method.

Results

After comprehensively analyzing various types of data from TCGA dataset, and exploring the differential expression and prognostic value of each m5C regulator, we identified m5C modification patterns based on 17 m5C regulators. Two patterns presented a significant difference in the risk of BCR, the tumor microenvironment (TME), and immunotherapy response in PRAD. We found that TET2, which was highly expressed in adjacent normal tissues compared to tumor tissues, was closely associated with many infiltrating immune cells. The m5C modification signature was constructed for the clinical application. Risk score calculated by m5C signature was associated with T stage, N stage, Gleason score, and the possibility of BCR (HR, 4.197; 95% CI, 3.016-5.842; p < 0.001). A higher risk score also represented the possibility of immunotherapy response. Finally, the potential roles of m5C modification signature were validated in the testing dataset.

Conclusions

Our study revealed the potential roles of m5C modification in the PRAD BCR and TME diversity, which may provide new insight into the field of prostate cancer in future research.

Comprehensive Analysis of 5-Methylcytosine (m5C) Regulators and the Immune Microenvironment in Pancreatic Adenocarcinoma to Aid Immunotherapy

Background

Pancreatic adenocarcinoma (PAAD) is one of the most malignant cancers and has a poor prognosis. As a critical RNA modification, 5-methylcytosine (m⁵C) has been reported to regulate tumor progression, including PAAD progression. However, a comprehensive analysis of m⁵C regulators in PAAD is lacking.

Methods

In the present study, PAAD datasets were obtained from the Gene Expression Omnibus (GEO), The Cancer Genome Atlas (TCGA), International Cancer Genome Consortium (ICGC), and ArrayExpress databases. The expression pattern of m⁵C regulators were analyzed and patients were divided into different m⁵C clusters according to consensus clustering based on m⁵C regulators. Additionally, m⁵C differentially expressed genes (DEGs) were determined using Limma package. Based on m⁵C DEGs, patients were divided into m⁵C gene clusters. Moreover, m⁵C gene signatures were derived from m⁵C DEGs and a quantitative indicator, the m⁵C score, was developed from the m⁵C gene signatures.

Results

Our study showed that m⁵C regulators were differentially expressed in patients with PAAD. The m⁵C clusters and gene clusters based on m⁵C regulators and m⁵C DEGs were related to immune cell infiltration, immune-related genes and patient survival status, indicating that m⁵C modification play a central role in regulating PAAD development partly by modulating immune microenvironment. Additionally, a quantitative indicator, the m⁵C score, was also developed and was related to a series of immune-related indicators. Moreover, the m⁵C score precisely predicted the immunotherapy response and prognosis of patients with PAAD.

Conclusion

In summary, we confirmed that m⁵C regulators regulate PAAD development by modulating the immune microenvironment. In addition, a quantitative indicator, the m⁵C score, was developed to predict immunotherapy response and prognosis and assisted in identifying PAAD patients suitable for tailored immunotherapy strategies.

m5CRegpred: Epitranscriptome Target Prediction of 5-Methylcytosine (m5C) Regulators Based on Sequencing Features

5-methylcytosine (m5C) is a common post-transcriptional modification observed in a variety of RNAs. m5C has been demonstrated to be important in a variety of biological processes, including RNA structural stability and metabolism. Driven by the importance of m5C modification, many projects focused on the m5C sites prediction were reported before. To better understand the upstream and downstream regulation of m5C, we present a bioinformatics framework, m5CRegpred, to predict the substrate of m5C writer NSUN2 and m5C readers YBX1 and ALYREF for the first time. After features comparison, window lengths selection and algorism comparison on the mature mRNA model, our model achieved AUROC scores 0.869, 0.724 and 0.889 for NSUN2, YBX1 and ALYREF, respectively in an independent test. Our work suggests the substrate of m5C regulators can be distinguished and may help the research of m5C regulators in a special condition, such as substrates prediction of hyper- or hypo-expressed m5C regulators in human disease.

Chemical biology and medicinal chemistry of RNA methyltransferases

RNA methyltransferases (MTases) are ubiquitous enzymes whose hitherto low profile in medicinal chemistry, contrasts with the surging interest in RNA methylation, the arguably most important aspect of the new field of epitranscriptomics. As MTases become validated as drug targets in all major fields of biomedicine, the development of small molecule compounds as tools and inhibitors is picking up considerable momentum, in academia as well as in biotech. Here we discuss the development of small molecules for two related aspects of chemical biology. Firstly, derivates of the ubiquitous cofactor S-adenosyl-l-methionine (SAM) are being developed as bioconjugation tools for targeted transfer of functional groups and labels to increasingly visible targets. Secondly, SAM-derived compounds are being investigated for their ability to act as inhibitors of RNA MTases. Drug development is moving from derivatives of cosubstrates towards higher generation compounds that may address allosteric sites in addition to the catalytic centre. Progress in assay development and screening techniques from medicinal chemistry have led to recent breakthroughs, e.g. in addressing human enzymes targeted for their role in cancer. Spurred by the current pandemic, new inhibitors against coronaviral MTases have emerged at a spectacular rate, including a repurposed drug which is now in clinical trial.

Biological roles of RNA m5C modification and its implications in Cancer immunotherapy

Epigenetics including DNA and RNA modifications have always been the hotspot field of life sciences in the post-genome era. Since the first mapping of N6-methyladenosine (m⁶A) and the discovery of its widespread presence in mRNA, there are at least 160-170 RNA modifications have been discovered. These methylations occur in different RNA types, and their distribution is species-specific. 5-methylcytosine (m⁵C) has been found in mRNA, rRNA and tRNA of representative organisms from all kinds of species. As reversible epigenetic modifications, m⁵C modifications of RNA affect the fate of the modified RNA molecules and play important roles in various biological processes including RNA stability control, protein synthesis, and transcriptional regulation. Furthermore, accumulative evidence also implicates the role of RNA m⁵C in tumorigenesis. Here, we review the latest progresses in the biological roles of m⁵C modifications and how it is regulated by corresponding “writers”, “readers” and “erasers” proteins, as well as the potential molecular mechanism in tumorigenesis and cancer immunotherapy.

RNA 5-Methylcytosine Regulators Contribute to Metabolism Heterogeneity and Predict Prognosis in Ovarian Cancer

5-Methylcytosine (m⁵C) is an abundant and highly conserved modification in RNAs. The dysregulation of RNA m⁵C methylation has been reported in cancers, but the regulatory network in ovarian cancer of RNA m⁵C methylation-related genes and its implication in metabolic regulation remain largely unexplored. In this study, RNA-sequencing data and clinical information of 374 ovarian cancer patients were downloaded from The Cancer Genome Atlas database, and a total of 14 RNA m⁵C regulators were included. Through unsupervised consensus clustering, two clusters with different m⁵C modification patterns were identified with distinct survivals. According to enrichment analyses, glycosaminoglycan and collagen metabolism–related pathways were specifically activated in cluster 1, whereas fatty acid metabolism–related pathways were enriched in cluster 2, which had better overall survival (OS). Besides the metabolism heterogeneity, the higher sensitivity to platinum and paclitaxel in cluster 2 can further explain the improved OS. Ultimately, a least absolute shrinkage and selection operator prediction model formed by ALYREF, NOP2, and TET2 toward OS was constructed. In conclusion, distinct m⁵C modification pattern exhibited metabolism heterogeneity, different chemotherapy sensitivity, and consequently survival difference, providing evidence for risk stratification.

Distinct Roles of m5C RNA Methyltransferase NSUN2 in Major Gynecologic Cancers

RNA methylation has recently emerged as an important category of epigenetic modifications, which plays diverse physiopathological roles in various cancers. Recent studies have confirmed the presence of 5-methylcytosine (m5C) modification on mammalian mRNAs, mainly modified by NOP2/Sun RNA methyltransferase family member 2 (NSUN2), but little is known about the underlying functions of m5C. Gynecologic cancers are malignancies starting from women’s reproductive organs. The prevalence of gynecologic cancers leads to a massive economic burden and public health concern. In this study, we investigated the potential biological functions of NSUN2 in common gynecologic cancers including cervical cancer, ovarian cancer, and endometrial cancer. Remarkably, distinct scenarios were found. The levels of NSUN2 did not show alteration in endometrial cancer, and in ovarian cancer, depletion of upregulated NSUN2 did not reduce carcinogenesis in cancer cells, suggesting that the upregulated NSUN2 might be an incidental effect. On the contrary, NSUN2 played a role in tumorigenesis of cervical cancer; depletion of upregulated NSUN2 notably inhibited migration and invasion of cancer cells, and only wild-type but not catalytically inactive NSUN2 rescued these malignant phenotypes of cancer cells. Mechanistically, NSUN2 promoted migration and invasion by leading to m5C methylation on keratin 13 (KRT13) transcripts, and methylated KRT13 transcripts would be recognized and stabilized by an m5C reader, Y-box binding protein 1 (YBX1). Collectively, these results not only displayed the nature of diversity among human malignancies, but also demonstrated a novel NSUN2-dependent m5C-YBX1-KRT13 oncogenic regulatory pathway.

Exploration of Potential Roles of m5C-Related Regulators in Colon Adenocarcinoma Prognosis

Objectives: The purpose of this study was to investigate the role of 13 m⁵C-related regulators in colon adenocarcinoma (COAD) and determine their prognostic value.

Methods: Gene expression and clinicopathological data were obtained from The Cancer Genome Atlas (TCGA) datasets. The expression of m⁵C-related regulators was analyzed with clinicopathological characteristics and alterations within m⁵C-related regulators. Subsequently, different subtypes of patients with COAD were identified. Then, the prognostic value of m⁵C-related regulators in COAD was confirmed via univariate Cox regression and least absolute shrinkage and selection operator (LASSO) Cox regression analyses. The prognostic value of risk scores was evaluated using the Kaplan-Meier method, receiver operating characteristic (ROC) curve. The correlation between the two m⁵C-related regulators, risk score, and clinicopathological characteristics were explored. Additionally, Gene Set Enrichment Analysis (GSEA), Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, and Gene Ontology (GO) analysis were performed for biological functional analysis. Finally, the expression level of two m⁵C-related regulators in clinical samples and cell lines was detected by quantitative reverse transcription-polymerase chain reaction and through the Human Protein Atlas database.

Results: m⁵C-related regulators were found to be differentially expressed in COAD with different clinicopathological features. We observed a high alteration frequency in these genes, which were significantly correlated with their mRNA expression levels. Two clusters with different prognostic features were identified. Based on two independent prognostic m⁵C-related regulators (NSUN6 and ALYREF), a risk signature with good predictive significance was constructed. Univariate and multivariate Cox regression analyses suggested that the risk score was an independent prognostic factor. Furthermore, this risk signature could serve as a prognostic indicator for overall survival in subgroups of patients with different clinical characteristics. Biological processes and pathways associated with cancer, immune response, and RNA processing were identified.

Conclusion: We revealed the genetic signatures and prognostic values of m⁵C-related regulators in COAD. Together, this has improved our understanding of m⁵C RNA modification and provided novel insights to identify predictive biomarkers and develop molecular targeted therapy for COAD.

m5C-Atlas: a comprehensive database for decoding and annotating the 5-methylcytosine (m5C) epitranscriptome

5-Methylcytosine (m5C) is one of the most prevalent covalent modifications on RNA. It is known to regulate a broad variety of RNA functions, including nuclear export, RNA stability and translation. Here, we present m5C-Atlas, a database for comprehensive collection and annotation of RNA 5-methylcytosine. The database contains 166 540 m5C sites in 13 species identified from 5 base-resolution epitranscriptome profiling technologies. Moreover, condition-specific methylation levels are quantified from 351 RNA bisulfite sequencing samples gathered from 22 different studies via an integrative pipeline. The database also presents several novel features, such as the evolutionary conservation of a m5C locus, its association with SNPs, and any relevance to RNA secondary structure. All m5C-atlas data are accessible through a user-friendly interface, in which the m5C epitranscriptomes can be freely explored, shared, and annotated with putative post-transcriptional mechanisms (e.g. RBP intermolecular interaction with RNA, microRNA interaction and splicing sites). Together, these resources offer unprecedented opportunities for exploring m5C epitranscriptomes. The m5C-Atlas database is freely accessible at https://www.xjtlu.edu.cn/biologicalsciences/m5c-atlas.

Ovarian cancer subtypes based on the regulatory genes of RNA modifications: Novel prediction model of prognosis

BACKGROUND

Ovarian cancer (OC) is a female reproductive system tumor. RNA modifications play key roles in gene expression regulation. The growing evidence demonstrates that RNA methylation is critical for various biological functions, and that its dysregulation is related to the progression of cancer in human.

METHOD

OC samples were classified into different subtypes (Clusters 1 and 2) based on various RNA-modification regulatory genes (RRGs) in the process of RNA modifications (m1A, m6A, m6Am, m5C, m7G, ac4C, m3C, and Ψ) by nonnegative matrix factorization method (NMF). Based on differently expressed RRGs (DERRGs) between clusters, a pathologically specific RNA-modification regulatory gene signature was constructed with Lasso regression. Kaplan-Meier analysis and receiver operating characteristic (ROC) curves were used to evaluate the prognostic ability of the identified model. The correlations of clinicopathological features, immune subtypes, immune scores, immune cells, and tumor mutation burden (TMB) were also estimated between different NMF clusters and riskscore groups.

RESULTS

In this study, 59 RRGs in the process of RNA modifications (m1A, m6A, m6Am, m5C, m7G, ac4C, m3C, and Ψ) were obtained from TCGA database. These RRGs were interactional, and sample clusters based on these regulators were significantly correlated with survival rate, clinical characteristics (involving survival status and pathologic stage), drug sensibility, and immune microenvironment. Furthermore, Lasso regression based on these 21 DERRGs between clusters 1 and 2 constructed a four-DERRG signature (ALYREF, ZC3H13, WTAP, and METTL1). Based on this signature, 307 OC patients were classified into high- and low-risk groups based on median value of riskscores from lasso regression. This identified signature was significantly associated with overall survival, radiation therapy, age, clinical stage, cancer status, and immune cells (involving CD4+ memory resting T cells, plasma cells, and Macrophages M1) of ovarian cancer patients. Further, GSEA revealed that multiple biological behaviors were significantly enriched in different groups.

CONCLUSIONS

OC patients were classified into two subtypes per these RRGs. This study identified four-DERRG signature (ALYREF, ZC3H13, WTAP, and METTL1) in OC, which was an independent prognostic model for patient stratification, prognostic evaluation, and prediction of response to immunotherapy in ovarian cancer by classifying OC patients into high- and low-risk groups.

RNA 5-methylcytosine regulates YBX2-dependent liquid-liquid phase separation

5-Methylcytosine (m5C) is one of the most prevalent internal modifications of messenger RNA (mRNA) in higher eukaryotes. Here we report that Y box protein 2 (YBX2) serves as a novel mammalian m5C binding protein to undergo liquid-liquid phase separation (LLPS) both in vivo and in vitro, and this YBX2-dependent LLPS is enhanced by m5C marked RNA. Furthermore, the crystal structure assay revealed that W100, as a distinct m5C binding site of YBX2, is critical in mediating YBX2 phase separation. Our study resolved the relationship between RNA m5C and phase separation, providing a clue for a new regulatory layer of epigenetics.

RNA nucleotide methylation: 2021 update

Among RNA modifications, transfer of methylgroups from the typical cofactor S-adenosyl-l-methionine by methyltransferases (MTases) to RNA is by far the most common reaction. Since our last review about a decade ago, the field has witnessed the re-emergence of mRNA methylation as an important mechanism in gene regulation. Attention has then spread to many other RNA species; all being included into the newly coined concept of the "epitranscriptome." The focus moved from prokaryotes and single cell eukaryotes as model organisms to higher eukaryotes, in particular to mammals. The perception of the field has dramatically changed over the past decade. A previous lack of phenotypes in knockouts in single cell organisms has been replaced by the apparition of MTases in numerous disease models and clinical investigations. Major driving forces of the field include methylation mapping techniques, as well as the characterization of the various MTases, termed "writers." The latter term has spilled over from DNA modification in the neighboring epigenetics field, along with the designations "readers," applied to mediators of biological effects upon specific binding to a methylated RNA. Furthermore "eraser" enzymes effect the newly discovered oxidative removal of methylgroups. A sense of reversibility and dynamics has replaced the older perception of RNA modification as a concrete-cast, irreversible part of RNA maturation. A related concept concerns incompletely methylated residues, which, through permutation of each site, lead to inhomogeneous populations of numerous modivariants. This review recapitulates the major developments of the past decade outlined above, and attempts a prediction of upcoming trends. This article is categorized under: RNA Processing > RNA Editing and Modification.

Emerging Role of Epitranscriptomics in Diabetes Mellitus and Its Complications

Diabetes mellitus (DM) and its related complications are among the leading causes of disability and mortality worldwide. Substantial studies have explored epigenetic regulation that is involved in the modifications of DNA and proteins, but RNA modifications in diabetes are still poorly investigated. In recent years, posttranscriptional epigenetic modification of RNA (the so-called 'epitranscriptome') has emerged as an interesting field of research. Numerous modifications, mainly N⁶ -methyladenosine (m⁶A), have been identified in nearly all types of RNAs and have been demonstrated to have an indispensable effect in a variety of human diseases, such as cancer, obesity, and diabetes. Therefore, it is particularly important to understand the molecular basis of RNA modifications, which might provide a new perspective for the pathogenesis of diabetes mellitus and the discovery of new therapeutic targets. In this review, we aim to summarize the recent progress in the epitranscriptomics involved in diabetes and diabetes-related complications. We hope to provide some insights for enriching the understanding of the epitranscriptomic regulatory mechanisms of this disease as well as the development of novel therapeutic targets for future clinical benefit.

The RNA Methylation Modification 5-Methylcytosine Impacts Immunity Characteristics, Prognosis and Progression of Oral Squamous Cell Carcinoma by Bioinformatics Analysis

Disorders pertaining to 5-methylcytosine (m5C) modifications are involved in the pathological process of many diseases. However, the effect of m5C on the tumorigenesis and progression of oral squamous cell carcinoma (OSCC) remains unclear. In this study, we integrated the genomic and clinical data of 558 OSCC samples to comprehensively evaluate m5C modification patterns. Based on 16 m5C methylation regulators, two m5C modification clusters were identified with distinct tumor immune microenvironment (TIME) characteristics and prognosis in OSCC. We then performed weighted gene co-expression network analysis (WGCNA) to identify m5C modification cluster-related modules. Genes in the selected module were chosen to construct the m5Cscore scoring system for evaluating m5C modification pattern in individual OSCC patients. Patients with a high m5Cscore had higher immune, stromal, and ESTIMATE scores; lower tumor purity score; lower immune activity; and higher tumor mutational burden. The overall survival rate and progression-free survival rate were markedly worse and the tumor recurrence rate was higher in OSCC patients with a high m5Cscore. Furthermore, patients with oral leukoplakia who also had a high m5Cscore had a higher risk of deterioration to OSCC. This study demonstrated that m5C modification patterns might affect the TIME in OSCC. m5Cscore may provide a new approach for predicting the prognosis and progression of OSCC.

The Prognostic Value and Immune Landscapes of a m6A/m5C/m1A-Related LncRNAs Signature in Head and Neck Squamous Cell Carcinoma

Background: N6-methyladenosine (m6A), 5-methylcytosine (m5C) and N1-methyladenosine (m1A) are the main RNA methylation modifications involved in the progression of cancer. However, it is still unclear whether m6A/m5C/m1A-related long non-coding RNAs (lncRNAs) affect the prognosis of head and neck squamous cell carcinoma (HNSCC). Methods: We summarized 52 m6A/m5C/m1A-related genes, downloaded 44 normal samples and 501 HNSCC tumor samples with RNA-seq data and clinical information from The Cancer Genome Atlas (TCGA) database, and then searched for m6A/m5C/m1A-related genes co-expressed lncRNAs. We adopt the least absolute shrinkage and selection operator (LASSO) Cox regression to obtain m6A/m5C/m1A-related lncRNAs to construct a prognostic signature of HNSCC. Results: This prognostic signature is based on six m6A/m5C/m1A-related lncRNAs (AL035587.1, AC009121.3, AF131215.5, FMR1-IT1, AC106820.5, PTOV1-AS2). It was found that the high-risk subgroup has worse overall survival (OS) than the low-risk subgroup. Moreover, the results showed that most immune checkpoint genes were significantly different between the two risk groups (p < 0.05). Immunity microenvironment analysis showed that the contents of NK cell resting, macrophages M2, and neutrophils in samples of low-risk group were significantly lower than those of high-risk group (p < 0.05), while the contents of B cells navie, plasma cells, and T cells regulatory (Tregs) were on the contrary (p < 0.05). In addition, patients with high tumor mutational burden (TMB) had the worse overall survival than those with low tumor mutational burden. Conclusion: Our study elucidated how m6A/m5C/m1A-related lncRNAs are related to the prognosis, immune microenvironment, and TMB of HNSCC. In the future, these m6A/m5C/m1A-related lncRNAs may become a new choice for immunotherapy of HNSCC.

Comprehensive Analysis of the Prognostic Value and Immune Infiltrates of the Three-m5C Signature in Colon Carcinoma

Background

The 5-methylcytosine (m5C) is one of the important forms of RNA post modification, and its regulatory mechanism in tumors has received increasing attention. However, its potential role in colorectal cancer remains unclear.

Materials and Methods

Here, we systematically investigated the genetic variation and prognostic value of the 14 m5c RNA methylation regulators in colon cancer. The prognostic risk score was constructed using three m5C regulators, which was verified in the GSE17536 (N=177), GSE41258 (N=248) and GSE38832 (N=122) datasets.

Results

The risk score developed from the three-m5C signature represents an independent prognostic factor, which can accurately predict the prognosis of patients with colon cancer in multiple datasets. The cytokine–cytokine receptor interaction and chemokine signaling pathway were significantly enriched in the low-risk score group. Further analysis showed that the three-m5C signature was related to tumor immune microenvironment (TIME), affecting the abundance of tumor-infiltrating immune cells. Especially, patients with low risk score had higher immune score than those with high risk score. In addition, gene set enrichment analysis (GSEA) confirmed that all three regulatory factors are associated with the MAPK/p38 signaling pathway.

Conclusion

In conclusion, our study illustrates that the three-m5C signature may be involved in the regulation of colon cancer immune microenvironment in synergy with the MAPK signaling pathway. Therefore, further studying the three-m5C signature regulatory mechanisms might provide promising targets for improving the responsiveness of colon cancer to immunotherapy.