The RNA methyltransferase NSUN6 suppresses pancreatic cancer development by regulating cell proliferation

Dysregulation of tRNA methylation in cancer: Mechanisms and targeting therapeutic strategies

tRNA is the RNA type that undergoes the most modifications among known RNA, and in recent years, tRNA methylation has emerged as a crucial process in regulating gene translation. Dysregulation of tRNA abundance occurs in cancer cells, along with increased expression and activity of tRNA methyltransferases to raise the level of tRNA modification and stability. This leads to hijacking of translation and synthesis of multiple proteins associated with tumor proliferation, metastasis, invasion, autophagy, chemotherapy resistance, and metabolic reprogramming. In this review, we provide an overview of current research on tRNA methylation in cancer to clarify its involvement in human malignancies and establish a theoretical framework for future therapeutic interventions targeting tRNA methylation processes.

Extracellular vesicle-packaged PIAT from cancer-associated fibroblasts drives neural remodeling by mediating m5C modification in pancreatic cancer mouse models

Perineural invasion (PNI) is a biological characteristic commonly observed in pancreatic cancer. Although PNI plays a key role in pancreatic cancer metastasis, recurrence, and poor postoperative survival, its mechanism is largely unclarified. Clinical sample analysis and endoscopic ultrasonographic elasticity scoring indicated that cancer-associated fibroblasts (CAFs) were closely related to the occurrence of PNI. Furthermore, CAF-derived extracellular vesicles (EVs) were involved in PNI in dorsal root ganglion coculture and mouse sciatic nerve models. Next, we demonstrated that CAFs promoted PNI through extracellular vesicle transmission of PNI-associated transcript (PIAT). Mechanistically, PIAT specifically bound to YBX1 and blocked the YBX1-Nedd4l interaction to inhibit YBX1 ubiquitination and degradation. Furthermore, PIAT enhanced the binding of YBX1 and PNI-associated mRNAs in a 5-methylcytosine (m5C)-dependent manner. Mutation of m5C recognition motifs in YBX1 or m5C sites in downstream target genes reversed PIAT-mediated PNI. Consistent with these findings, analyses using a KPC mouse model demonstrated that the PIAT/YBX1 axis enhanced PNI through m5C modification. Clinical data suggested that the PIAT expression in the serum EVs of patients with pancreatic cancer was associated with the degree of neural invasion and prognosis. Our study revealed the important role of the PIAT/YBX1 signaling axis in the tumor microenvironment (TME) in promoting tumor cell PNI and provided a new target for precise interference with CAFs and RNA methylation in the TME to suppress PNI in pancreatic cancer.

RNA modification gene WDR4 facilitates tumor progression and immunotherapy resistance in breast cancer

INTRODUCTION

Growing interest toward RNA modification in cancer has inspired the exploration of gene sets related to multiple RNA modifications. However, a comprehensive elucidation of the clinical value of various RNA modifications in breast cancer is still lacking.

OBJECTIVES

This study aimed to provide a strategy based on RNA modification-related genes for predicting therapy response and survival outcomes in breast cancer patients.

METHODS

Genes related to thirteen RNA modification patterns were integrated for establishing a nine-gene-containing signature-RMscore. Alterations of tumor immune microenvironment and therapy response featured by different RMscore levels were assessed by bulk transcriptome, single-cell transcriptome and genomics analyses. The biological function of key RMscore-related molecules was investigated by cellular experiments in vitro and in vivo, using flow cytometry, immunohistochemistry and immunofluorescence staining.

RESULTS

This study has raised an effective therapy strategy for breast cancer patients after a well-rounded investigation of RNA modification-related genes. With a great performance of predicting patient prognosis, high levels of the RMscore proposed in this study represented suppressive immune microenvironment and therapy resistance, including adjuvant chemotherapy and PD-L1 blockade treatment. As the key contributor of the RMscore, inhibition of WDR4 impaired breast cancer progression significantly in vitro and in vivo, as well as participated in regulating cell cycle and mTORC1 signaling pathway via m7G modification.

CONCLUSION

Briefly, this study has developed promising and effective tactics to achieve the prediction of survival probabilities and treatment response in breast cancer patients.

Biological function molecular pathways and druggability of DNMT2/TRDMT1

5-methylcytosine (m5C) is among the most common epigenetic modification in DNA and RNA molecules, and plays an important role in the animal development and disease pathogenesis. Interestingly, unlike other m5C DNA methyltransferases (DNMTs), DNMT2/TRDMT1 has the double-substrate specificity and adopts a DNMT-similar catalytic mechanism to methylate RNA. Moreover, it is widely involved in a variety of physiological regulatory processes, such as the gene expression, precise protein synthesis, immune response, and disease occurrence. Thus, comprehending the epigenetic mechanism and function of DNMT2/TRDMT1 will probably provide new strategies to treat some refractory diseases. Here, we discuss recent studies on the spatiotemporal expression pattern and post-translational modifications of DNMT2/TRDMT1, and summarize the research advances in substrate characteristics, catalytic recognition mechanism, DNMT2/TRDMT1-related genes or proteins, pharmacological application, and inhibitor development. This review will shed light on the pharmacological design by targeting DNMT2/TRDMT1 to treat parasitic, viral and oncologic diseases.

The role of RNA methylation in tumor immunity and its potential in immunotherapy

RNA methylation, a prevalent post-transcriptional modification, has garnered considerable attention in research circles. It exerts regulatory control over diverse biological functions by modulating RNA splicing, translation, transport, and stability. Notably, studies have illuminated the substantial impact of RNA methylation on tumor immunity. The primary types of RNA methylation encompass N6-methyladenosine (m6A), 5-methylcytosine (m5C), N1-methyladenosine (m1A), and N7-methylguanosine (m7G), and 3-methylcytidine (m3C). Compelling evidence underscores the involvement of RNA methylation in regulating the tumor microenvironment (TME). By affecting RNA translation and stability through the "writers", "erasers" and "readers", RNA methylation exerts influence over the dysregulation of immune cells and immune factors. Consequently, RNA methylation plays a pivotal role in modulating tumor immunity and mediating various biological behaviors, encompassing proliferation, invasion, metastasis, etc. In this review, we discussed the mechanisms and functions of several RNA methylations, providing a comprehensive overview of their biological roles and underlying mechanisms within the tumor microenvironment and among immunocytes. By exploring how these RNA modifications mediate tumor immune evasion, we also examine their potential applications in immunotherapy. This review aims to provide novel insights and strategies for identifying novel targets in RNA methylation and advancing cancer immunotherapy efficacy.

NSUN6 mediates 5-methylcytosine modification of METTL3 and promotes colon adenocarcinoma progression

Colon adenocarcinoma (COAD) is a common and fatal malignant tumor of digestive system with complex etiology. 5-Methylcytosine (m5C) modification of RNA by the NSUN gene family (NSUN1-NSUN7) and DNMT2 reshape cell biology and regulate tumor development. However, the expression profile, prognostic significance and function of these m5C modifiers in COAD remain largely unclear. By mining multiple integrated tumor databases, we found that NSUN1, NSUN2, NSUN5, and NSUN6 were overexpressed in COAD tumor samples relative to normal samples. Clinically, high expression of NSUN6 was significantly associated with shorter survival (including both disease-free survival and overall survival) in COAD patients. NSUN6 was further confirmed to be upregulated at both tissue and cellular levels of COAD, suggesting that NSUN6 plays a critical role in disease progression. Through comprehensive gene enrichment analysis and cell-based functional validation, it was revealed that NSUN6 promoted the cell cycle progression and cell proliferation of COAD. Mechanistically, NSUN6 upregulates the expression of oncogenic METTL3 and catalyzes its m5C modification in COAD cells. Overexpression of METTL3 significantly relieved the cell cycle inhibition of COAD caused by NSUN6 deficiency. Furthermore, NSUN6 was negatively associated with the abundance of infiltrating immune cells in COAD tumors, such as activated B cells, natural killer cells, effector memory CD8 T cells, and regulatory T cells. Importantly, pan-cancer analysis further uncovered that NSUN6 was dysregulated and heterogeneous in various tumors. Thus our findings extend the role of m5C transferase in COAD and suggest that NSUN6 is a potential biomarker and target for this malignancy.

Changes of RNA m6A/m5C Modification Regulatory Molecules in Ferroptosis of T2DM Rat Pancreas

N6-methyladenine (m6A) and 5-methylcytosine (m5C) are two common forms of RNA methylation that play an important role in the epigenetics of type 2 diabetes mellitus (T2DM). One type of cell death, ferroptosis, has been implicated in islet β-cell damage in T2DM. Notably, RNA methylation, an upstream regulatory mechanism of mRNAs, can regulate the expression of ferroptosis signaling molecules, thereby affecting cell proliferation and death. Here, we found that the ferroptosis signaling pathway was activated in pancreas of T2DM rats, followed by significant changes in m6A/m5C modification regulatory molecules. These detection data together with the prediction results that m6A and m5C exist in the mRNAs of ferroptosis molecules, we speculate that m6A and m5C are probably involved in pancreatic cell damage by modifying of ferroptosis signaling molecules. In short, our findings provide a new research idea for future studies on the molecular mechanisms of pancreatic cell damage and point to a new direction for exploring the mechanisms of ferroptosis from the perspective of RNA methylation modification.

Emerging function of main RNA methylation modifications in the immune microenvironment of digestive system tumors

Digestive system tumors have been reported in more than 25% of all cancer cases worldwide, bringing a huge burden on the healthcare system. RNA methylation modification-an important post-transcriptional modification-has become an active research area in gene regulation. It is a dynamic and reversible process involving several enzymes, such as methyltransferases, demethylases, and methylation reader proteins. This review provides insights into the role of three major methylation modifications, namely m6A, m5C, and m1A, in the development of digestive system tumors, specifically in the development of tumor immune microenvironment (TIME) of these malignancies. Abnormal methylation modification affects immunosuppression and antitumor immune response by regulating the recruitment of immune cells and the release of immune factors. Understanding the mechanisms by which RNA methylation regulates digestive system tumors will be helpful in exploring new therapeutic targets.

RNA modification-related genes illuminate prognostic signature and mechanism in esophageal squamous cell carcinoma

Emerging studies have demonstrated the link between RNA modifications and various cancers, while the predictive value and functional mechanisms of RNA modification-related genes (RMGs) in esophageal squamous cell carcinoma (ESCC) remain unclear. Here we established a prognostic signature for ESCC based on five RMGs. The analysis of ESCC clinical samples further verified the prognostic power of the prognostic signature. Moreover, we found that the knockdown of NSUN6 promotes ESCC progression in vitro and in vivo, whereas the overexpression of NSUN6 inhibits the malignant phenotype of ESCC cells. Mechanically, NSUN6 mediated tRNA m5C modifications selectively enhance the translation efficiency of CDH1 mRNA in a codon dependent manner. Rescue assays revealed that E-cadherin is an essential downstream target that mediates NSUN6's function in the regulation of ESCC progression. These findings offer additional insights into the link between ESCC and RMGs, as well as provide potential strategies for ESCC management and therapy.

NSUN2 promotes colorectal cancer progression by enhancing SKIL mRNA stabilization

BACKGROUND

NOP2/Sun domain 2 (NSUN2) is one of the important RNA methyltransferases catalyzing 5-methylcytosine (m5C) formation and participates in many critical bioprocesses. However, the roles and underlying molecular mechanisms of NSUN2-mediated m5C modification in colorectal cancer (CRC) remain unclear.

METHODS

To explore the NSUN2 expression in CRC, fresh tissue samples were collected and Nsun2 knockout mouse was constructed. In vitro and in vivo functional assays were conducted to assess the role of NSUN2. RNA array and bisulfite sequencing were used to investigate the potential targets. The mechanisms of NSUN2 function on SKIL were identified by m5C-methylated-RNA immunoprecipitation and RNA stability assays. Additionally, tissue microarray analysis was conducted and patient-derived tumour xenograft mouse (PDX) models were used to define the potential therapeutic targets.

RESULTS

NSUN2 was highly expressed in CRC and correlated with poor CRC patient survival. Moreover, silencing NSUN2 suppressed CRC tumourigenesis and progression in Nsun2 knockout mouse models. In vitro and in vivo studies suggested that NSUN2 promoted colorectal cancer cell growth. Mechanistically, SKI-like proto-oncogene (SKIL) is positively regulated by NSUN2, and the NSUN2-SKIL axis is clinically relevant to CRC. NSUN2 induced m5C modification of SKIL and stabilized its mRNA, which was mediated by Y-box binding protein 1 (YBX1). Elevated SKIL levels increased transcriptional coactivator with PDZ-binding motif (TAZ) activation.

CONCLUSIONS

Our findings highlight the importance of NSUN2 in the initiation and progression of CRC via m5C-YBX1-dependent stabilization of the SKIL transcript, providing a promising targeted therapeutic strategy for CRC.

The impact of tRNA modifications on translation in cancer: identifying novel therapeutic avenues

Recent advancements have illuminated the critical role of RNA modifications in post-transcriptional regulation, shaping the landscape of gene expression. This review explores how tRNA modifications emerge as critical players, fine-tuning functionalities that not only maintain the fidelity of protein synthesis but also dictate gene expression and translation profiles. Highlighting their dysregulation as a common denominator in various cancers, we systematically investigate the intersection of both cytosolic and mitochondrial tRNA modifications with cancer biology. These modifications impact key processes such as cell proliferation, tumorigenesis, migration, metastasis, bioenergetics and the modulation of the tumor immune microenvironment. The recurrence of altered tRNA modification patterns across different cancer types underscores their significance in cancer development, proposing them as potential biomarkers and as actionable targets to disrupt tumorigenic processes, offering new avenues for precision medicine in the battle against cancer.

The Role of tRNA-Centered Translational Regulatory Mechanisms in Cancer

Cancer is a leading cause of morbidity and mortality worldwide. While numerous factors have been identified as contributing to the development of malignancy, our understanding of the mechanisms involved remains limited. Early cancer detection and the development of effective treatments are therefore critical areas of research. One class of molecules that play a crucial role in the transmission of genetic information are transfer RNAs (tRNAs), which are the most abundant RNA molecules in the human transcriptome. Dysregulated synthesis of tRNAs directly results in translation disorders and diseases, including cancer. Moreover, various types of tRNA modifications and the enzymes responsible for these modifications have been implicated in tumor biology. Furthermore, alterations in tRNA modification can impact tRNA stability, and impaired stability can prompt the cleavage of tRNAs into smaller fragments known as tRNA fragments (tRFs). Initially believed to be random byproducts lacking any physiological function, tRFs have now been redefined as non-coding RNA molecules with distinct roles in regulating RNA stability, translation, target gene expression, and other biological processes. In this review, we present recent findings on translational regulatory models centered around tRNAs in tumors, providing a deeper understanding of tumorigenesis and suggesting new directions for cancer treatment.

Deciphering the vital roles and mechanism of m5C modification on RNA in cancers

5-methylcytosine (m5C modification) plays an essential role in tumors, which affects different types of RNA, the expression of downstream target genes, and downstream pathways, thus participating in the tumor process. However, the effect of m5C modification on RNA in tumors and the exact mechanism have not been systematically reviewed. Therefore, we reviewed the status and sites of m5C modification, as well as the expression pattern and biological functions of m5C regulators in tumors, and further summarized the effects and regulation mechanism of m5C modification on messenger RNA (mRNA), ribosomal RNA (rRNA), transfer RNA (tRNA), long non-coding RNA (lncRNA) and other RNA in tumors. Finally, we summed up the interaction network, potential application, and value in clinical diagnosis and treatment of tumors. Taken together, this review benefits revealing the mechanism of m5C modification in tumor progression and provides new strategies for tumor diagnosis and treatment.

NSUN6 Regulates NM23-H1 Expression in an m5C Manner to Affect Epithelial-Mesenchymal Transition in Lung Cancer

PURPOSE

The expression and regulatory mechanism of NSUN6 in lung cancer are still unclear. Our study explored whether NSUN6 mediates progression of lung cancer by affecting NM23-H1 expression in an m5C-dependent manner.

METHODS

qRT-PCR, CCK-8, colony formation, transwell, and Western blot analysis were employed to probe the impact of NSUN6 on lung cancer cell proliferation, migration, and epithelial-mesenchymal transition (EMT). RMVar database was utilized to forecast the downstream genes of NSUN6. The mode of interaction between NSUN6 and NM23-H1 was determined by dot blot, luciferase assay, m5C RIP, and cell function assays. The effect of NSUN6 expression on tumor growth was verified in vivo.

RESULTS

Expression of NSUN6 was reduced in lung cancer cells, and over-expression of NSUN6 restricted the proliferation of lung cancer cells, migration, and EMT. NSUN6 regulated NM23-H1 expression by modifying the 3'-UTR of NM23-H1 mRNA through m5C and inhibited lung cancer cell proliferation, migration, and EMT. In vivo experiments also showed that over-expression of NSUN6 inhibited the occurrence of lung cancer.

CONCLUSION

NSUN6 regulates NM23-H1 expression in an m5C-dependent manner to affect EMT in lung cancer. Thus, NSUN6 may be considered as a potential therapeutic target for lung cancer.

Do RNA modifications contribute to modulation of immune responses in allergic diseases?

RNA modifications have emerged as a fundamental mechanism of post-transcriptional gene regulation, playing vital roles in cellular physiology and the development of various diseases. While the investigation of RNA modifications has seen significant advancements, the exploration of their implication in allergic diseases has been comparatively overlooked. Allergic reactions, including hay fever, asthma, eczema and food allergies, result from hypersensitive immune responses, affecting a considerable population worldwide. Despite the high prevalence, the molecular mechanisms underlying these responses remain partially understood. The potential role of RNA modifications in modulating the hypersensitive immune responses has yet to be fully elucidated. This mini-review seeks to shed light on potential connections between RNA modifications and allergy, highlighting recent findings and potential future research directions. By expanding our understanding of the complex interplay between RNA modifications and allergic responses, we hope to unlock new avenues for allergy diagnosis, prognosis, and therapeutic intervention.

Identification and validation of 5-methylcytosine-associated genes in diffuse large B-cell lymphoma

5-methylcytosine modifications play a significant role in carcinogenesis; however, studies exploring 5-methylcytosine-related genes in diffuse large B-cell lymphoma patients are lacking. In this study, we aimed to understand the potential role and clinical prognostic impact of 5-methylcytosine regulators in diffuse large B-cell lymphoma and identify a prognostic biomarker based on 5-methylcytosine-associated genes. Gene expression profiles and corresponding clinical information of diffuse large B-cell lymphoma patients and normal controls were obtained from The Cancer Genome Atlas, Gene Expression Omnibus, and Genotype-Tissue Expression databases. Diffuse large B-cell lymphoma was divided into three clusters according to the 5-methylcytosine regulators, and differentially expressed genes were screened among the three clusters. Univariate Cox and Lasso-Cox regression analyses were used to screen prognostic genes and construct a prognostic risk model. Kaplan-Meier curve analysis, univariate and multivariate Cox regression analyses, and time-dependent receiver operator characteristic curve analysis were used to evaluate prognostic factors. GSVA was used to enrich potential pathways associated with 5-methylcytosine modification patterns. SsGSEA and CIBERSORT were used to assess immune cell infiltration. Six 5-methylcytosine-related genes (TUBB4A, CD3E, ZNF681, HAP1, IL22RA2, and POSTN) were used to construct a prognostic risk model, which was proved to have a good predictive effect. In addition, univariate and multivariate Cox regression risk scores were independent prognostic factors for diffuse large B-cell lymphoma. Further analysis showed that the 5-methylcytosine risk score was significantly correlated with immune cell infiltration and immune checkpoint of diffuse large B-cell lymphoma. Our study reveals for the first time a potential role for 5-methylcytosine modifications in diffuse large B-cell lymphoma, provides novel insights for future studies on diffuse large B-cell lymphoma, and offers potential prognostic biomarkers and therapeutic targets for patients with diffuse large B-cell lymphoma.

Identification and Validation of Genes Related to RNA Methylation Modification in Diabetic Retinopathy

PURPOSE

To identify and validate the differentially expressed genes related to RNA methylation modification in diabetic retinopathy.

METHODS

The data sets GSE12610 and GSE111465 related to diabetic retinopathy in the Gene Expression Omnibus were selected. The R software package was used to identify differentially expressed genes related to RNA methylation modification in diabetic retinopathy. Protein-protein interaction network was constructed to explore the interactions between proteins and predict proteins. Then, Gene Ontology annotation analysis and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis were used to analyze the potential enrichment pathways and clarify the biological functions of these genes. In addition, the correlation between them and immune cells was visualized, and receiver operating characteristic curves were drawn to evaluate the diagnostic performance of each one of them for diabetic retinopathy. To verify the differentially expressed genes, the mRNA expression of rat retinal vascular endothelial cells cultured in low and high glucose medium separately were detected by RT-qPCR.

RESULTS

The expression of Lrpprc, Nsun4, Nsun6 and Trdmt1 were significantly up-regulated in diabetic retinopathy samples, while the expression of Cbll1, Hnrnpc, Mettl3 and Wtap were significantly down-regulated. Differentially expressed genes were mainly enriched in the RNA-methylation-medication pathways and biological function. The results of immune infiltration analysis proved that eosinophils aggregated more in diabetic group, while T cells follicular helper aggregated more in normal samples. These genes of Cbll1 (AUC = 0.986), Hnrnpc (AUC = 0.819), Lrpprc (AUC = 0.806), Mettl3 (AUC = 0.917), Nsun4 (AUC = 0.819), Nsun6 (AUC = 0.819), Trdmt1 (AUC = 0.972) and Wtap (AUC = 0.972) were respectively used as the diagnostic basis of diabetic retinopathy. According to the RT-qPCR results, the expression of Mettl3 was significantly down-regulated (p < 0.0005) in cells cultured in high glucose, while Trdmt1 (p < 0.05), Nsun4 (p < 0.05) and Nsun6 (p < 0.05) were significantly up-regulated.

CONCLUSION

Differentially expressed genes such as Mettl3, Nsun4, Nsun6, and Trdmt1 should be conducted to explore, and the role of RNA methylation in the process of diabetic retinopathy would be revealed in-depth.

Disease Diagnosis Based on Nucleic Acid Modifications

Nucleic acid modifications include a wide range of epigenetic and epitranscriptomic factors and impact a wide range of nucleic acids due to their profound influence on biological inheritance, growth, and metabolism. The recently developed methods of mapping and characterizing these modifications have promoted their discovery as well as large-scale studies in eukaryotes, especially in humans. Because of these pioneering strategies, nucleic acid modifications have been shown to have a great impact on human disorders such as cancer. Therefore, whether nucleic acid modifications could become a new type of biomarker remains an open question. In this review, we briefly look back at classical nucleic acid modifications and then focus on the progress made in investigating these modifications as diagnostic biomarkers in clinical therapy and present our perspective on their development prospects.

The integrated single-cell analysis developed an immunogenic cell death signature to predict lung adenocarcinoma prognosis and immunotherapy

BACKGROUND

Research on immunogenic cell death (ICD) in lung adenocarcinoma (LUAD) has been relatively limited. This study aims to create ICD-related signatures for accurate survival prognosis prediction in LUAD patients, addressing the challenge of lacking reliable early prognostic indicators for this type of cancer.

METHODS

Using single-cell RNA sequencing (scRNA-seq) analysis, ICD activity in cells was calculated by AUCell algorithm, divided into high- and low-ICD groups according to median values, and key ICD regulatory genes were identified through differential analysis, and these genes were integrated into TCGA data to construct prognostic signatures using LASSO and COX regression analysis, and multi-dimensional analysis of ICD-related signatures in terms of prognosis, immunotherapy, tumor microenvironment (TME), and mutational landscape.

RESULTS

The constructed signature reveals a pronounced disparity in prognosis between the high- and low-risk groups of LUAD patients. The statistical discrepancies in survival times among LUAD patients from both the TCGA and GEO databases further corroborate this observation. Additionally, heightened levels of immune cell infiltration expression are evidenced in the low-risk group, suggesting a potential benefit from immunotherapeutic interventions for these patients. The expression levels of pivotal risk-associated genes in tissue samples were assessed utilizing qRT-PCR, thereby unveiling PITX3 as a plausible therapeutic target in the context of LUAD.

CONCLUSIONS

Our constructed ICD-related signatures provide help in predicting the prognosis and immunotherapy of LUAD patients, and to some extent guide the clinical treatment of LUAD patients.

Regulation and functions of non-m6A mRNA modifications

Nucleobase modifications are prevalent in eukaryotic mRNA and their discovery has resulted in the emergence of epitranscriptomics as a research field. The most abundant internal (non-cap) mRNA modification is N6-methyladenosine (m6A), the study of which has revolutionized our understanding of post-transcriptional gene regulation. In addition, numerous other mRNA modifications are gaining great attention because of their major roles in RNA metabolism, immunity, development and disease. In this Review, we focus on the regulation and function of non-m6A modifications in eukaryotic mRNA, including pseudouridine (Ψ), N6,2'-O-dimethyladenosine (m6Am), N1-methyladenosine (m1A), inosine, 5-methylcytidine (m5C), N4-acetylcytidine (ac4C), 2'-O-methylated nucleotide (Nm) and internal N7-methylguanosine (m7G). We highlight their regulation, distribution, stoichiometry and known roles in mRNA metabolism, such as mRNA stability, translation, splicing and export. We also discuss their biological consequences in physiological and pathological processes. In addition, we cover research techniques to further study the non-m6A mRNA modifications and discuss their potential future applications.

RNA cytosine methyltransferase NSUN5 promotes protein synthesis and tumorigenic phenotypes in glioblastoma

Glioblastoma (GBM) is the most common and aggressive malignant primary brain tumor in adults. The standard treatment achieves a median overall survival for GBM patients of only 15 months. Hence, novel therapies based on an increased understanding of the mechanistic underpinnings of GBM are desperately needed. In this study, we show that elevated expression of 28S rRNA (cytosine-C(5))-methyltransferase NSUN5, which methylates cytosine 3782 of 28S rRNA in GBM cells, is strongly associated with the poor survival of GBM patients. Moreover, we demonstrate that overexpression of NSUN5 increases protein synthesis in GBM cells. NSUN5 knockdown decreased protein synthesis, cell proliferation, sphere formation, migration, and resistance to temozolomide in GBM cell lines. NSUN5 knockdown also decreased the number and size of GBM neurospheres in vitro. As a corollary, mice harboring U251 tumors wherein NSUN5 was knocked down survived longer than mice harboring control tumors. Taken together, our results suggest that NSUN5 plays a protumorigenic role in GBM by enabling the enhanced protein synthesis requisite for tumor progression. Accordingly, NSUN5 may be a hitherto unappreciated target for the treatment of GBM.

Loss of NSUN6 inhibits osteosarcoma progression by downregulating EEF1A2 expression and activation of Akt/mTOR signaling pathway via m5C methylation

As an important 5-methylcytidine (m5C) methyltransferase, NOP2/Sun RNA methyltransferase family member 6 (NSUN6) has been reported to play an important role in the progression of several diseases. However, the role of NSUN6 in the progression of osteosarcoma (OS) remains unclear. This study aimed to identify the role of NSUN6 in the progression of OS and clarify the potential molecular mechanism. The present study discovered that NSUN6 was upregulated in OS and a higher NSUN6 expression was a strong indicator for poorer prognosis of patients with OS. In addition, the loss of NSUN6 led to reduced proliferation, migration and invasion of OS cells. Through bioinformatics analysis, RNA immunoprecipitation (RIP) and methylated RIP assays, eukaryotic elongation factor 1 α-2 (EEF1A2) was identified and validated as a potential target of NSUN6 in OS. Mechanistically, the expression of EEF1A2 was significantly suppressed following NSUN6 knockdown due to reduced EEF1A2 mRNA stability in an m5C-dependent manner. Meanwhile, NSUN6 deficiency inhibited m5C-dependent activation of Akt/mTOR signaling pathway. In addition, genetic overexpression of EEF1A2 or pharmacological activation of the Akt signaling pathway counteracted the suppressive effects of NSUN6 deficiency on the proliferation, invasion and migration of OS cells. The current findings suggested that NSUN6 may serve as a potential therapeutic target for OS treatment.

RNA epigenetic modifications in ovarian cancer: The changes, chances, and challenges

Ovarian cancer (OC) is the most common female cancer worldwide. Patients with OC have high mortality because of its complex and poorly understood pathogenesis. RNA epigenetic modifications, such as m6 A, m1 A, and m5 C, are closely associated with the occurrence and development of OC. RNA modifications can affect the stability of mRNA transcripts, nuclear export of RNAs, translation efficiency, and decoding accuracy. However, there are few overviews that summarize the link between m6 A RNA modification and OC. Here, we discuss the molecular and cellular functions of different RNA modifications and how their regulation contributes to the pathogenesis of OC. By improving our understanding of the role of RNA modifications in the etiology of OC, we provide new perspectives for their use in OC diagnosis and treatment. This article is categorized under: RNA Processing > RNA Editing and Modification RNA in Disease and Development > RNA in Disease.

RNA modification: mechanisms and therapeutic targets

RNA modifications are dynamic and reversible chemical modifications on substrate RNA that are regulated by specific modifying enzymes. They play important roles in the regulation of many biological processes in various diseases, such as the development of cancer and other diseases. With the help of advanced sequencing technologies, the role of RNA modifications has caught increasing attention in human diseases in scientific research. In this review, we briefly summarized the basic mechanisms of several common RNA modifications, including m6A, m5C, m1A, m7G, Ψ, A-to-I editing and ac4C. Importantly, we discussed their potential functions in human diseases, including cancer, neurological disorders, cardiovascular diseases, metabolic diseases, genetic and developmental diseases, as well as immune disorders. Through the "writing-erasing-reading" mechanisms, RNA modifications regulate the stability, translation, and localization of pivotal disease-related mRNAs to manipulate disease development. Moreover, we also highlighted in this review all currently available RNA-modifier-targeting small molecular inhibitors or activators, most of which are designed against m6A-related enzymes, such as METTL3, FTO and ALKBH5. This review provides clues for potential clinical therapy as well as future study directions in the RNA modification field. More in-depth studies on RNA modifications, their roles in human diseases and further development of their inhibitors or activators are needed for a thorough understanding of epitranscriptomics as well as diagnosis, treatment, and prognosis of human diseases.

Revealing the role of regulatory T cells in the tumor microenvironment of lung adenocarcinoma: a novel prognostic and immunotherapeutic signature

Background

Regulatory T cells (Tregs), are a key class of cell types in the immune system. In the tumor microenvironment (TME), the presence of Tregs has important implications for immune response and tumor development. Relatively little is known about the role of Tregs in lung adenocarcinoma (LUAD).

Methods

Tregs were identified using but single-cell RNA sequencing (scRNA-seq) analysis and interactions between Tregs and other cells in the TME were investigated. Next, we used multiple bulk RNA-seq datasets to construct risk models based on marker genes of Tregs and explored differences in prognosis, mutational landscape, immune cell infiltration and immunotherapy between high- and low-risk groups, and finally, qRT-PCR and cell function experiments were performed to validate the model genes.

Results

The cellchat analysis showed that MIF-(CD74+CXCR4) pairs play a key role in the interaction of Tregs with other cell subpopulations, and the Tregs-associated signatures (TRAS) could well classify multiple LUAD cohorts into high- and low-risk groups. Immunotherapy may offer greater potential benefits to the low-risk group, as indicated by their superior survival, increased infiltration of immune cells, and heightened expression of immune checkpoints. Finally, the experiment verified that the model genes LTB and PTTG1 were relatively highly expressed in cancer tissues, while PTPRC was relatively highly expressed in paracancerous tissues. Colony Formation assay confirmed that knockdown of PTTG1 reduced the proliferation ability of LUAD cells.

Conclusion

TRAS were constructed using scRNA-seq and bulk RNA-seq to distinguish patient risk subgroups, which may provide assistance in the clinical management of LUAD patients.

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.

Emerging roles of tRNA in cancer

Transfer RNAs (tRNAs) play pivotal roles in the transmission of genetic information, and abnormality of tRNAs directly leads to translation disorders and causes diseases, including cancer. The complex modifications enable tRNA to execute its delicate biological function. Alteration of appropriate modifications may affect the stability of tRNA, impair its ability to carry amino acids, and disrupt the pairing between anticodons and codons. Studies confirmed that dysregulation of tRNA modifications plays an important role in carcinogenesis. Furthermore, when the stability of tRNA is impaired, tRNAs are cleaved into small tRNA fragments (tRFs) by specific RNases. Though tRFs have been found to play vital regulatory roles in tumorigenesis, its formation process is far from clear. Understanding improper tRNA modifications and abnormal formation of tRFs in cancer is conducive to uncovering the role of metabolic process of tRNA under pathological conditions, which may open up new avenues for cancer prevention and treatment.

Identification of m6A/m5C/m1A-associated LncRNAs for prognostic assessment and immunotherapy in pancreatic cancer

Methylation of RNA plays an important role in cancer. Classical forms of such modifications include N6-methyladenine (m6A), 5-methylcytosine (m5C), and N1-methyladenine (m1A). Methylation-regulated long non-coding (lnc) RNAs are involved in various biological processes, such as tumor proliferation, apoptosis, immune escape, invasion, and metastasis. Therefore, we performed an analysis of transcriptomic and clinical data of pancreatic cancer samples in The Cancer Genome Atlas (TCGA). Using the co-expression method, we summarized 44 m6A/m5C/m1A-related genes and obtained 218 methylation-associated lncRNAs. Next, with COX regression, we screened 39 lncRNAs that are strongly associated with prognosis and found that their expression differed significantly between normal tissues and pancreatic cancer samples (P < 0.001). We then used the least absolute shrinkage and selection operator (LASSO) to construct a risk model comprising seven lncRNAs. In validation set, the nomogram generated by combining clinical characteristics accurately predicted the survival probability of pancreatic cancer patients at 1, 2, and 3 years after diagnosis (AUC = 0.652, 0.686, and 0.740, respectively). Tumor microenvironment analysis showed that the high-risk group had significantly more resting memory CD4 T cells, M0 macrophages, and activated dendritic cells and fewer naïve B cells, plasma cells, and CD8 T cells than the low-risk group (both P < 0.05). Most immune-checkpoint genes were significantly different between the high- and low-risk groups (P < 0.05). The Tumor Immune Dysfunction and Exclusion score showed that high-risk patients benefited more from treatment with immune checkpoint inhibitors (P < 0.001). Overall survival was also lower in high-risk patients with more tumor mutations than in low-risk patients with fewer mutations (P < 0.001). Finally, we explored the sensitivity of the high- and low-risk groups to seven candidate drugs. Our findings indicated that m6A/m5C/m1A-associated lncRNAs are potentially useful biomarkers for the early diagnosis and estimating the prognosis of, and ascertaining the responses to immunotherapy in, patients with pancreatic cancer.

DNMT2/TRDMT1 gene knockout compromises doxorubicin-induced unfolded protein response and sensitizes cancer cells to ER stress-induced apoptosis

The acidic, hypoxic and nutrient-deprived tumor microenvironment may induce endoplasmic reticulum (ER) stress and the unfolded protein response (UPR) may exert an important cytoprotective role by promoting folding of newly synthesized proteins and cancer cell survival. The lack of DNMT2/TRDMT1 methyltransferase-mediated C38 tRNA methylation compromises translational fidelity that may result in the accumulation of misfolded and aggregated proteins leading to proteotoxic stress-related cell death. In the present study, DNMT2/TRDMT1 gene knockout-mediated effects were investigated during doxorubicin (DOX)-induced ER stress and PERK-, IRE1- and ATF6-orchestrated UPR in four genetically different cellular models of cancer (breast and cervical cancer, osteosarcoma and glioblastoma cells). Upon DOX stimulation, DNMT2/TRDMT1 gene knockout impaired PERK activation and modulated NSUN and 5-methylcytosine RNA-based responses and microRNA profiles. The lack of DNMT2/TRDMT1 gene in DOX-treated four cancer cell lines resulted in decreased levels of four microRNAs, namely, miR-23a-3p, miR-93-5p, miR-125a-5p and miR-191-5p involved in the regulation of several pathways such as ubiquitin-mediated proteolysis, amino acid degradation and translational misregulation in cancer. We conclude that DNMT2/TRDMT1 gene knockout, at least in selected cellular cancer models, affects adaptive responses associated with protein homeostasis networks that during prolonged ER stress may result in increased sensitivity to apoptotic cell death.

PLAGL1 is associated with prognosis and cell proliferation in pancreatic adenocarcinoma

BACKGROUND

Emerging evidence has shown the crucial roles of pleomorphic adenoma gene (PLAG) family genes in multiple cancers. However, their functions and mechanisms in pancreatic adenocarcinoma (PAAD) remain poorly understood.

METHODS

We analyzed the expression levels of PLAG family genes in both The Cancer Genome Atlas (TCGA) database and a Gene Expression Omnibus (GEO) database, and confirmed the results in our three independent cohorts of 382 PAAD tissues and 362 adjacent nontumor pancreatic tissues. Integrated analyses were carried out to explore the function, mechanism and prognostic value of the selected PLAG family gene in PAAD patients.

RESULTS

By analyzing the TCGA and GEO databases, PLAGL1 was identified to be downregulated in PAAD tissues, and its decreasing levels of both mRNA and protein were verified in our three independent PAAD cohorts. PLAGL1 expression was inversely correlated with clinicopathological factors including the Ki67⁺ cell rate and pathologic stage. Further GSEA of the TCGA-PAAD cohort demonstrated that multiple signaling pathways implicated in cell proliferation were enriched in the lower PLAGL1 expressing PAAD group. Moreover, we demonstrated that PLAGL1 expression was obviously negatively associated with patients' overall survival outcome in both the TCGA-PAAD cohort and our verification cohorts. Additionally, through MTS and BrdU assays, we further demonstrated in vitro that PLAGL1 had the impact of preventing the proliferation of pancreatic cancer cells.

CONCLUSIONS

Our present study suggested that downregulated PLAGL1 might act as a biomarker in predicts poor prognosis and one of important factors in increasing cell proliferation in PAAD. This study provides us with a novel prognostic marker and therapeutic strategy for PAAD, which deserves further study.

Genetic and epigenetic defects of the RNA modification machinery in cancer

Cancer was initially considered to be an exclusively genetic disease, but an interplay of dysregulated genetic and epigenetic mechanisms is now known to contribute to the cancer phenotype. More recently, chemical modifications of RNA molecules - the so-called epitranscriptome - have been found to regulate various aspects of RNA function and homeostasis. Specific enzymes, known as RNA-modifying proteins (RMPs), are responsible for depositing, removing, and reading chemical modifications in RNA. Intensive investigations in the epitranscriptomic field in recent years, in conjunction with great technological advances, have revealed the critical role of RNA modifications in regulating numerous cellular pathways. Furthermore, growing evidence has revealed that RNA modification machinery is often altered in human cancers, highlighting the enormous potential of RMPs as pharmacological targets or diagnostic markers.

5-Methylcytosine (m5C) modification in peripheral blood immune cells is a novel non-invasive biomarker for colorectal cancer diagnosis

Current non-invasive tumor biomarkers failed to accurately identify patients with colorectal cancer (CRC), delaying CRC diagnosis and thus leading to poor prognosis. Dysregulation of 5-Methylcytosine (m⁵C) RNA has gradually been reported in various cancers, but their role in tumor diagnosis is rarely mentioned. Our study aimed to determine the role of m⁵C methylation modification in blood immune cells for the diagnosis of CRC. Peripheral blood samples were obtained from a total of 83 healthy controls and 196 CRC patients. We observed that m⁵C RNA contents in blood immune cells of CRC patients were markedly enhanced in both training set and validation set. Moreover, levels of m⁵C increased with CRC progression and metastasis but reduced after treatment. Compared with common blood tumor biomarkers, m⁵C levels in peripheral blood immune cells had superior discrimination and reclassification performance in diagnosing CRC. Besides, bioinformatics and qRT-PCR analysis identified increased expression of m⁵C-modified regulators NSUN5 and YBX1 in CRC patients’ blood. A series of animal models and cell co-culture models further demonstrated that CRC tumor cells could increase immune cells’ m⁵C levels and m⁵C-modified regulators. Monocyte was the predominant m⁵C-modified immune cell type in CRC patients’ blood by Gene set variation analysis (GSVA). Taken together, m⁵C methylation modification in peripheral blood immune cells was a promising biomarker for non-invasive diagnosis of CRC.

Prognostic Value and Immune Landscapes of m5C-Related lncRNAs in Lung Squamous Cell Carcinoma

5-methylcytosine (m5C) modification is involved in tumor progression. However, the lncRNAs associated with m5C in lung squamous cell carcinoma (LUSC) have not been elucidated. The Cancer Genome Atlas database was used to get the open-accessed transcriptional profiling and clinical information of LUSC patients. All the statistical analyses were performed based on R software v 4.0.0 and SPSS13.0. First, there were 614 m5C-related lncRNAs identified under the criterion of |R|>0.4 and p < 0.001 with m5C genes. Next, a prognosis model based on ERICD, AL021068.1, LINC01341, AC254562.3, and AP002360.1 was established, which showed good prediction efficiency in both the training and validation cohorts. Next, a nomogram plot was established by combining the risk score and clinical features for a better application in clinical settings. Pathway enrichment analysis showed that the pathways of angiogenesis, TGF-β signaling, IL6-JAK-STAT3 signaling, protein secretion, androgen response, interferon-α response, and unfolded protein response were significantly enriched in the high-risk patients. Immune infiltration analysis showed that the risk score was positively correlated with neutrophils, resting CD4+ memory T cells, and M2 macrophages, yet negatively correlated with follicular helper T cells, CD8+ T cells, and activated NK cells. Moreover, we found that high-risk patients might be more sensitive to immunotherapy, imatinib, yet resistant to erlotinib, gefitinib, and vinorelbine. In summary, our prognosis model is an effective tool that could robustly predict LUSC patient prognosis, which had the potential for clinical guidance.

Genetic variants in RNA m5 C modification genes associated with survival and chemotherapy efficacy of colorectal cancer

BACKGROUND

Colorectal cancer is one of the most common malignant digestive tract tumors with a poor prognosis. RNA 5-methylcytosine (m⁵ C) is an important posttranscriptional widespread modification involved in many biological processes. However, the association between genetic variations of m⁵ C modification genes and the prognostic value of colorectal cancer remains unclear.

METHODS

We investigated the association between candidate single nucleotide polymorphisms (SNPs) in 13 m⁵ C modification genes and colorectal cancer overall survival (OS) after chemotherapy by the Cox regression model. The combined effect of selected SNPs on OS, progression-free survival (PFS), and disease control rate (DCR) was assessed by the number of risk alleles (NRA). The GTEx and TCGA database were used to perform expression qualitative trait locus (eQTL) analysis.

RESULTS

We identified that two SNPs in YBX1 were associated with OS after chemotherapy (HR = 1.43, p = 0.001 for rs10890208; HR = 1.36, p = 0.025 for rs3862218). A striking dose-response effect between NRA and OS after chemotherapy was found (p_trend  = 0.002). The DCR of patients receiving oxaliplatin chemotherapy in the 3-4 NRA group was markedly reduced in comparison to that in the 0-2 NRA group (OR = 1.49, p = 0.036). Moreover, YBX1 mRNA expression was significantly overexpressed in tumor tissues (p < 0.05) in the TCGA database, and eQTL analysis demonstrated that the two SNPs were associated with YBX1 (p = 0.003 for rs10890208 and p = 0.024 for rs3862218).

CONCLUSION

Our study indicates that genetic variants in m⁵ C modification genes may mediate changes in YBX1 mRNA levels and affect the chemotherapeutic efficacy of colorectal cancer patients.

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.

Role of main RNA modifications in cancer: N6-methyladenosine, 5-methylcytosine, and pseudouridine

Cancer is one of the major diseases threatening human life and health worldwide. Epigenetic modification refers to heritable changes in the genetic material without any changes in the nucleic acid sequence and results in heritable phenotypic changes. Epigenetic modifications regulate many biological processes, such as growth, aging, and various diseases, including cancer. With the advancement of next-generation sequencing technology, the role of RNA modifications in cancer progression has become increasingly prominent and is a hot spot in scientific research. This review studied several common RNA modifications, such as N⁶-methyladenosine, 5-methylcytosine, and pseudouridine. The deposition and roles of these modifications in coding and noncoding RNAs are summarized in detail. Based on the RNA modification background, this review summarized the expression, function, and underlying molecular mechanism of these modifications and their regulators in cancer and further discussed the role of some existing small-molecule inhibitors. More in-depth studies on RNA modification and cancer are needed to broaden the understanding of epigenetics and cancer diagnosis, treatment, and prognosis.

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.

Cross-Talk Between m6A- and m5C-Related lncRNAs to Construct a Novel Signature and Predict the Immune Landscape of Colorectal Cancer Patients

Background

N6-methyladenosine (m⁶A) and 5-methylcytosine (m⁵C) can modify long non-coding RNAs (lncRNAs), thereby affecting tumorigenesis and tumor progression. However, there is a lack of knowledge regarding the potential roles and cross-talk of m⁶A- and m⁵C-related lncRNAs in the tumor microenvironment (TME) and their effect on prognosis.

Methods

We systematically evaluated the expression patterns of m⁶A- and m⁵C-related lncRNAs in 1358 colorectal cancer (CRC) samples from four datasets. Consensus clustering was conducted to identify molecular subtypes of CRC, and the clinical significance, TME, tumor-infiltrating immune cells (TIICs), and immune checkpoints in the different molecular subtypes were analyzed. Finally, we established a m⁶A- and m⁵C-related lncRNA signature and a prognostic nomogram.

Results

We identified 141 m⁶A- and m⁵C-related lncRNAs by co-expression analysis, among which 23 lncRNAs were significantly associated with the overall survival (OS) of CRC patients. Two distinct molecular subtypes (cluster A and cluster B) were identified, and these two distinct molecular subtypes could predict clinicopathological features, prognosis, TME stromal activity, TIICs, immune checkpoints. Next, a m⁶A- and m⁵C-related lncRNA signature for predicting OS was constructed, and its predictive capability in CRC patients was validated. We then constructed a highly accurate nomogram for improving the clinical applicability of the signature. Analyses of clinicopathological features, prognosis, TIICs, cancer stem cell (CSC), and drug response revealed significant differences between two risk groups. In addition, we found that patients with a low-risk score exhibited enhanced response to anti-PD-1/L1 immunotherapy. Functional enrichment analysis showed that these lncRNAs related to the high-risk group were involved in the development and progression of CRC.

Conclusions

We conducted a comprehensive analysis of m⁶A- and m⁵C-related lncRNAs in CRC and revealed their potential functions in predicting tumor-immune-stromal microenvironment, clinicopathological features, and prognosis, and determined their role in immunotherapy. These findings may improve our understanding of the cross-talk between m⁶A- and m⁵C-related lncRNAs in CRC and pave a new road for prognosis assessment and more effective immunotherapy strategies.

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.

Roles of RNA Modifications in Diverse Cellular Functions

Chemical modifications of RNA molecules regulate both RNA metabolism and fate. The deposition and function of these modifications are mediated by the actions of writer, reader, and eraser proteins. At the cellular level, RNA modifications regulate several cellular processes including cell death, proliferation, senescence, differentiation, migration, metabolism, autophagy, the DNA damage response, and liquid-liquid phase separation. Emerging evidence demonstrates that RNA modifications play active roles in the physiology and etiology of multiple diseases due to their pervasive roles in cellular functions. Here, we will summarize recent advances in the regulatory and functional role of RNA modifications in these cellular functions, emphasizing the context-specific roles of RNA modifications in mammalian systems. As m⁶A is the best studied RNA modification in biological processes, this review will summarize the emerging advances on the diverse roles of m⁶A in cellular functions. In addition, we will also provide an overview for the cellular functions of other RNA modifications, including m⁵C and m¹A. Furthermore, we will also discuss the roles of RNA modifications within the context of disease etiologies and highlight recent advances in the development of therapeutics that target RNA modifications. Elucidating these context-specific functions will increase our understanding of how these modifications become dysregulated during disease pathogenesis and may provide new opportunities for improving disease prevention and therapy by targeting these pathways.

Identification of a Methylation-Regulating Genes Prognostic Signature to Predict the Prognosis and Aid Immunotherapy of Clear Cell Renal Cell Carcinoma

Methylation is one of the most extensive modifications of biological macromolecules and affects cell-fate determination, development, aging, and cancer. Several methylation modifications, including 5-methylcytosine and N6-methyladenosine, play an essential role in many cancers. However, little is known about the relationship between methylation and the prognosis of clear cell renal cell carcinoma (ccRCC). Here, we established a methylation-regulating genes prognostic signature (MRGPS) to predict the prognoses of ccRCC patients. We obtained ccRCC samples from The Cancer Genome Atlas and identified methylation-regulatingd genes (MRGs) from the Gene Set Enrichment Analysis database. We also determined differentially expressed genes (DEGs) and performed cluster analysis to identify candidate genes. Subsequently, we established and validated an MRGPS to predict the overall survival of ccRCC patients. This was also verified in 15 ccRCC samples collected from the Fujian Provincial Hospital via quantitative real-time transcription (qRT-PCR). While 95 MRGs were differentially expressed (DEGs1) between tumor and normal tissues, 17 MRGs were differentially expressed (DEGs2) between cluster 1 and 2. Notably, 13 genes common among DEGs1 and DEGs2 were identified as hub genes. In fact, we established three genes (NOP2, NSUN6, and TET2) to be an MRGPS based on their multivariate Cox regression analysis coefficients (p &lt; 0.05). A receiver operating characteristic curve analysis confirmed this MRGPS to have a good prognostic performance. Moreover, the MRGPS was associated with characteristics of the tumor immune microenvironment and responses to inhibitor checkpoint inhibitors. Data from “IMvigor 210” demonstrated that patients with a low MRGPS would benefit more from atelozumab (p &lt; 0.05). Furthermore, a multivariate analysis revealed that MRGPS was an independent risk factor associated with ccRCC prognosis (p &lt; 0.05). Notably, a nomogram constructed by combining with clinical characteristics (age, grade, stage, and MRGPS risk score) to predict the overall survival of a ccRCC patient had a favorable predictive value. Eventually, our qRT-PCR results showed that tumor tissues had higher NOP2 and NSUN6 expression levels and lower TET2 expression than normal tissues of ccRCC samples. While the proposed MRGPS comprising NOP2, NSUN6, and TET2 can be an alternative prognostic biomarker for ccRCC patients, it is a promising index for personalized ICI treatments against ccRCC.

An Overview of Epigenetic Methylation in Pancreatic Cancer Progression

Over the past decades, the aberrant epigenetic modification, apart from genetic alteration, has emerged as dispensable events mediating the transformation of pancreatic cancer (PC). However, the understanding of molecular mechanisms of methylation modifications, the most abundant epigenetic modifications, remains superficial. In this review, we focused on the mechanistic insights of DNA, histone, and RNA methylation that regulate the progression of PC. The methylation regulators including writer, eraser and reader participate in the modification of gene expression associated with cell proliferation, invasion and apoptosis. Some of recent clinical trials on methylation drug targeting were also discussed. Understanding the novel regulatory mechanisms in the methylation modification may offer alternative opportunities to improve therapeutic efficacy to fight against this dismal disease.

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.

tRNA modifications and their potential roles in pancreatic cancer

Since the breakthrough discovery of N6-methyladenosine (m⁶A), the field of RNA epitranscriptomics has attracted increasing interest in the biological sciences. Transfer RNAs (tRNAs) are extensively modified, and various modifications play a crucial role in the formation and stability of tRNA, which is universally required for accurate and efficient functioning of tRNA. Abnormal tRNA modification can lead to tRNA degradation or specific cleavage of tRNA into fragmented derivatives, thus affecting the translation process and frequently accompanying a variety of human diseases. Increasing evidence suggests that tRNA modification pathways are also misregulated in human cancers. In this review, we summarize tRNA modifications and their biological functions, describe the type and frequency of tRNA modification alterations in cancer, and highlight variations in tRNA-modifying enzymes and the multiple functions that they regulate in different types of cancers. Furthermore, the current implications and the potential role of tRNA modifications in the progression of pancreatic cancer are discussed. Collectively, this review describes recent advances in tRNA modification in cancers and its potential significance in pancreatic cancer. Further study of the mechanism of tRNA modifications in pancreatic cancer may provide possibilities for therapies targeting enzymes responsible for regulating tRNA modifications in pancreatic cancer.

Prognostic Risk Model and Tumor Immune Environment Modulation of m5C-Related LncRNAs in Pancreatic Ductal Adenocarcinoma

RNA methylation modification is a key process in epigenetics that regulates posttranscriptional gene expression. With advances in next-generation sequencing technology, 5-methylcytosine (m5C) modification has also been found in multiple RNAs. Long non-coding RNAs (lncRNAs) were proved to have a key role in cancer progression and closely related to the tumor immune microenvironment. Thus, based on the PDAC patients' clinical information and genetic transcriptome data from the TCGA database, we performed a detailed bioinformatic analysis to establish a m5C-related lncRNA prognostic risk model for PDAC patients and discovered the relationship between the risk model and PDAC immune microenvironment. Pearson correlation coefficient analysis was applied to conduct a m5C regulatory gene and m5C-related lncRNA co-expression network. Expression of m5C-related lncRNAs screened by univariate regression analysis with prognostic value showed a significant difference between pancreatic cancer and normal tissues. The least absolute shrinkage and selection operator (LASSO) Cox regression method was applied to determine an 8-m5C-related lncRNA prognostic risk model. We used principal component analysis to indicate that the risk model could distinguish all the samples clearly. The clinical nomogram also accurately predicted 1-, 1.5-, 2-, and 3-year survival time among PDAC patients. Additionally, this risk model was validated in the entire group and sub-test groups using KM analysis and ROC analysis. Combined with the clinical characteristics, the risk score was found to be an independent factor for predicting the survival of PDAC patients. Furthermore, the association between the risk model and tumor immune microenvironment was evaluated via the ESTIMATE R package and CIBERSORT method. Consequently, the results indicated that immune cells were associated with m5C-related lncRNA risk model scores and had different distribution in the high- and low-risk groups. Based on all these analyses, the m5C-related lncRNA risk model could be a reliable prognostic tool and therapeutic target for PDAC patients.

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.

RNA methylation and cancer treatment

To this date, over 100 different types of RNA modification have been identified. Methylation of different RNA species has emerged as a critical regulator of transcript expression. RNA methylation and its related downstream signaling pathways are involved in plethora biological processes, including cell differentiation, sex determination and stress response, and others. It is catalyzed by the RNA methyltransferases, is demethylated by the demethylases (FTO and ALKBH5) and read by methylation binding protein (YTHDF1 and IGF2BP1). Increasing evidence indicates that this process closely connected to cancer cell proliferation, cellular stress, metastasis, immune response. And RNA methylation related protein has been becoming a promising targets of cancer therapy. This review outlines the relationship between different types of RNA methylation and cancer, and some FTO inhibitors in cancer treatment.

The role of RNA m5C modification in cancer metastasis

Epigenetic modification plays a crucial regulatory role in the biological processes of eukaryotic cells. The recent characterization of DNA and RNA methylation is still ongoing. Tumor metastasis has long been an unconquerable feature in the fight against cancer. As an inevitable component of the epigenetic regulatory network, 5-methylcytosine is associated with multifarious cellular processes and systemic diseases, including cell migration and cancer metastasis. Recently, gratifying progress has been achieved in determining the molecular interactions between m⁵C writers (DNMTs and NSUNs), demethylases (TETs), readers (YTHDF2, ALYREF and YBX1) and RNAs. However, the underlying mechanism of RNA m⁵C methylation in cell mobility and metastasis remains unclear. The functions of m⁵C writers and readers are believed to regulate gene expression at the post-transcription level and are involved in cellular metabolism and movement. In this review, we emphatically summarize the recent updates on m⁵C components and related regulatory networks. The content will be focused on writers and readers of the RNA m⁵C modification and potential mechanisms in diseases. We will discuss relevant upstream and downstream interacting molecules and their associations with cell migration and metastasis.

m5C-Related Signatures for Predicting Prognosis in Cutaneous Melanoma with Machine Learning

Background

Cutaneous melanoma (CM) is one of the most life-threatening primary skin cancers and is prone to distant metastases. A widespread presence of posttranscriptional modification of RNA, 5-methylcytosine (m5C), has been observed in human cancers. However, the potential mechanism of the tumorigenesis and prognosis in CM by dysregulated m5C-related regulators is obscure.

Methods

We use comprehensive bioinformatics analyses to explore the expression of m5C regulators in CM, the prognostic implications of the m5C regulators, the frequency of the copy number variant (CNV), and somatic mutations in m5C regulators. Additionally, the CM patients were divided into three clusters for better predicting clinical features and outcomes via consensus clustering of m5C regulators. Then, the risk score was established via Lasso Cox regression analysis. Next, the prognosis value and clinical characteristics of m5C-related signatures were further explored. Then, machine learning was used to recognize the outstanding m5C regulators to risk score. Finally, the expression level and clinical value of USUN6 were analyzed via the tissue microarray (TMA) cohort.

Results

We found that m5C regulators were dysregulated in CM, with a high frequency of somatic mutations and CNV alterations of the m5C regulatory gene in CM. Furthermore, 16 m5C-related proteins interacted with each other frequently, and we divided CM patients into three clusters to better predicting clinical features and outcomes. Then, five m5C regulators were selected as a risk score based on the LASSO model. The XGBoost algorithm recognized that NOP2 and NSUN6 were the most significant risk score contributors. Immunohistochemistry has verified that low expression of USUN6 was closely correlated with CM progression.

Conclusion

The m5C-related signatures can be used as new prognostic biomarkers and therapeutic targets for CM, and NSUN6 might play a vital role in tumorigenesis and malignant progression.

The lack of functional DNMT2/TRDMT1 gene modulates cancer cell responses during drug-induced senescence

Cellular senescence may be a side effect of chemotherapy and other anti-cancer treatments that may promote inflammation and paracrine secondary senescence in healthy tissues. DNMT2/TRDMT1 methyltransferase is implicated in the regulation of cellular lifespan and DNA damage response (DDR). In the present study, the responses to senescence inducing concentrations of doxorubicin and etoposide in different cancer cells with DNMT2/TRDMT1 gene knockout were evaluated, namely changes in the cell cycle, apoptosis, autophagy, interleukin levels, genetic stability and DDR, and 5-mC and NSUN1-6 levels. Moreover, the effect of azacytidine post-treatment was considered. Diverse responses were revealed that was based on type of cancer cells (breast and cervical cancer, osteosarcoma and glioblastoma cells) and anti-cancer drugs. DNMT2/TRDMT1 gene knockout in drug-treated glioblastoma cells resulted in decreased number of apoptotic and senescent cells, IL-8 levels and autophagy, and increased number of necrotic cells, DNA damage and affected DDR compared to drug-treated glioblastoma cells with unmodified levels of DNMT2/TRDMT1. We suggest that DNMT2/TRDMT1 gene knockout in selected experimental settings may potentiate some adverse effects associated with chemotherapy-induced senescence.