Uncovering the Association Between m5C Regulator-Mediated Methylation Modification Patterns and Tumour Microenvironment Infiltration Characteristics in Hepatocellular Carcinoma

RNA m5C methylation: a potential modulator of innate immune pathways in hepatocellular carcinoma

RNA 5-methylcytosine (m5C) methylation plays a crucial role in hepatocellular carcinoma (HCC). As reported, aberrant m5C methylation is closely associated with the progression, therapeutic efficacy, and prognosis of HCC. The innate immune system functions as the primary defense mechanism in the body against pathogenic infections and tumors since it can activate innate immune pathways through pattern recognition receptors to exert anti-infection and anti-tumor effects. Recently, m5C methylation has been demonstrated to affect the activation of innate immune pathways including TLR, cGAS-STING, and RIG-I pathways by modulating RNA function, unveiling new mechanisms underlying the regulation of innate immune responses by tumor cells. However, research on m5C methylation and its interplay with innate immune pathways is still in its infancy. Therefore, this review details the biological significance of RNA m5C methylation in HCC and discusses its potential regulatory relationship with TLR, cGAS-STING, and RIG-I pathways, thereby providing fresh insights into the role of RNA methylation in the innate immune mechanisms and treatment of HCC.

RNA methylation-related inhibitors: Biological basis and therapeutic potential for cancer therapy

RNA methylation is widespread in nature. Abnormal expression of proteins associated with RNA methylation is strongly associated with a number of human diseases including cancer. Increasing evidence suggests that targeting RNA methylation holds promise for cancer treatment. This review specifically describes several common RNA modifications, such as the relatively well-studied N6-methyladenosine, as well as 5-methylcytosine and pseudouridine (Ψ). The regulatory factors involved in these modifications and their roles in RNA are also comprehensively discussed. We summarise the diverse regulatory functions of these modifications across different types of RNAs. Furthermore, we elucidate the structural characteristics of these modifications along with the development of specific inhibitors targeting them. Additionally, recent advancements in small molecule inhibitors targeting RNA modifications are presented to underscore their immense potential and clinical significance in enhancing therapeutic efficacy against cancer. KEY POINTS: In this paper, several important types of RNA modifications and their related regulatory factors are systematically summarised. Several regulatory factors related to RNA modification types were associated with cancer progression, and their relationships with cancer cell migration, invasion, drug resistance and immune environment were summarised. In this paper, the inhibitors targeting different regulators that have been proposed in recent studies are summarised in detail, which is of great significance for the development of RNA modification regulators and cancer treatment in the future.

The impact of epitranscriptomic modifications on liver disease

RNA modifications have emerged as important mechanisms of gene regulation. Developmental, metabolic, and cell cycle regulatory processes are all affected by epitranscriptomic modifications, which control gene expression in a dynamic manner. The hepatic tissue is highly metabolically active and has an impressive ability to regenerate after injury. Cell proliferation, differentiation, and metabolism, which are all essential to the liver response to injury and regeneration, are regulated via RNA modification. Two such modifications, N6-methyladenosine (m6A)and 5-methylcytosine (m5C), have been identified as prognostic disease markers and potential therapeutic targets for liver diseases. Here, we describe progress in understanding the role of RNA modifications in liver biology and disease and discuss specific areas where unexpected results could lead to improved future understanding.

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

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

RNA m5C methylation orchestrates BLCA progression via macrophage reprogramming

Recently, epigenetics showed essential roles in tumour microenvironment (TME) and immunotherapy response, however, the functions of RNA 5-methylcytosine (m5C) modification in TME remains unknown. According to 13 m5C regulators, we evaluated 412 BLCA patients from The Cancer Genome Atlas (TCGA) database. The m5C score was constructed by unsupervised clustering analysis and principal component analysis (PCA) algorithms. Gene set variation analysis (GSVA), ESTIMATE algorithm, and immunohistochemical (IHC) staining were performed. Macrophage chemotaxis assay was used to assess the M2 macrophages. Among the 412 patients, the frequency of mutation was 13%. m5C regulators was expressed significantly in BLCA tissue compared with normal tissue. Then, two m5C methylation modification patterns were identified with dissimilar TME cell infiltration patterns. The C1 alteration pattern in the m5C cluster was connected with better survival. In addition, we found that NSUN6 was highly correlated with recruitment of macrophages via bioinformatics and IHC. Further experiment validated that NSUN6 promoted HDAC10 expression by mediating m5C methylation, inhibited the transcription of macrophage-associated chemokines and thus inhibited the recruitment of M2 macrophages. The m5C score constructed by m5C modification pattern showed that high m5C score group had a better prognosis. This study uncovered the significant roles of m5C modifications in modulating the TME and indicated that NSUN6 could inhibit the recruitment of M2 macrophages via m5C methylation, which provided novel insight into epigenetic regulation of TME and clinical suggestions for immunotherapeutic strategies.

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.

Developing an m5C regulator-mediated RNA methylation modification signature to predict prognosis and immunotherapy efficacy in rectal cancer

Background

Currently, a very small number of patients with colorectal cancer (CRC) respond to immune checkpoint inhibitor (ICI) treatment. Therefore, there is an urgent need to investigate effective biomarkers to determine the responsiveness to ICI treatment. Recently, aberrant 5-methylcytosine (m⁵C) RNA modification has emerged as a key player in the pathogenesis of cancer. Thus, we aimed to explore the predictive signature based on m⁵C regulator-related genes for characterizing the immune landscapes and predicting the prognosis and response to therapies.

Methods

The Cancer Genome Atlas (TCGA) cohort was used as the training set, while GEO data sets, real-time quantitative PCR (RT-qPCR) analysis from paired frozen tissues, and immunohistochemistry (IHC) data from tissue microarray (TMA) were used for validation. We constructed a novel signature based on three m⁵C regulator-related genes in patients with rectal adenocarcinoma (READ) using a least absolute shrinkage and selection operator (LASSO)-Cox regression and unsupervised consensus clustering analyses. Additionally, we correlated the three-gene signature risk model with the tumor immune microenvironment, immunotherapy efficiency, and potential applicable drugs.

Results

The m⁵C methylation-based signature was an independent prognostic factor, where low-risk patients showed a stronger immunoreactivity phenotype and a superior response to ICI therapy. Conversely, the high-risk patients had enriched pathways of cancer hallmarks and presented immune-suppressive state, which demonstrated that they are more insensitive to immunotherapy. Additionally, the signature markedly correlated with drug susceptibility.

Conclusions

We developed a reliable m⁵C regulator-based risk model to predict the prognosis, clarify the molecular and tumor microenvironment status, and identify patients who would benefit from immunotherapy or chemotherapy. Our study could provide vital guidance to improve prognostic stratification and optimize personalized therapeutic strategies for patients with rectal cancer.

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.

Scoring model based on the signature of non-m6A-related neoantigen-coding lncRNAs assists in immune microenvironment analysis and TCR-neoantigen pair selection in gliomas

Background

Small peptides encoded by long non-coding RNAs (lncRNAs) have attracted attention for their various functions. Recent studies indicate that these small peptides participate in immune responses and antigen presentation. However, the significance of RNA modifications remains unclear.

Methods

Thirteen non-m6A-related neoantigen-coding lncRNAs were selected for analysis from the TransLnc database. Next, a neoantigen activation score (NAS) model was established based on the characteristics of the lncRNAs. Machine learning was employed to expand the model to two additional RNA-seq and two single-cell sequencing datasets for further validation. The DLpTCR algorithm was used to predict T cell receptor (TCR)-peptide binding probability.

Results

The non-m6A-related NAS model predicted patients’ overall survival outcomes more precisely than the m6A-related NAS model. Furthermore, the non-m6A-related NAS was positively correlated with tumor cells’ evolutionary level, immune infiltration, and antigen presentation. However, high NAS gliomas also showed more PD-L1 expression and high mutation frequencies of T-cell positive regulators. Interestingly, results of intercellular communication analysis suggest that T cell-high neoplastic cell interaction is weaker in both of the NAS groups which might arise from decreased IFNGR1 expression. Moreover, we identified unique TCR-peptide pairs present in all glioma samples based on peptides encoded by the 13 selected lncRNAs. And increased levels of neoantigen-active TCR patterns were found in high NAS gliomas.

Conclusions

Our work suggests that non-m6A-related neoantigen-coding lncRNAs play an essential role in glioma progression and that screened TCR clonotypes might provide potential avenues for chimeric antigen receptor T cell (CAR-T) therapy for gliomas.

Prognostic stratification based on m5C regulators acts as a novel biomarker for immunotherapy in hepatocellular carcinoma

Background

Immunotherapy is a promising anti-cancer strategy in hepatocellular carcinoma (HCC). However, a limited number of patients can benefit from it. There are currently no reliable biomarkers available to find the potential beneficiaries. Methylcytosine (m⁵C) is crucial in HCC, but its role in forecasting clinical responses to immunotherapy has not been fully clarified.

Methods

In this study, we analyzed 371 HCC patients from The Cancer Genome Atlas (TCGA) database and investigated the expression of 18 m⁵C regulators. We selected 6 differentially expressed genes (DEGs) to construct a prognostic risk model as well as 2 m⁵C-related diagnostic models.

Results

The 1-, 3-, and 5-year area under the curve (AUC) of m⁵C scores for the overall survival (OS) was 0.781/0.762/0.711, indicating the m⁵C score system had an ideal distinction of prognostic prediction for HCC. The survival analysis showed that patients with high-risk scores present a worse prognosis than the patients with low-risk scores (p< 0.0001). We got consistent results in 6 public cohorts and validated them in Xiangya real-world cohort by quantitative real-time PCR and immunohistochemical (IHC) assays. The high-m⁵C score group was predicted to be in an immune evasion state and showed low sensitivity to immunotherapy, but high sensitivity to chemotherapy and potential targeted drugs and agents, such as sepantronium bromide (YM-155), axitinib, vinblastine and docetaxel. Meanwhile, we also constructed two diagnostic models to distinguish HCC tumors from normal liver tissues or liver cirrhosis.

Conclusion

In conclusion, our study helps to early screen HCC patients and select patients who can benefit from immunotherapy. Step forwardly, for the less likely beneficiaries, this study provides them with new potential targeted drugs and agents for choice to improve their prognosis.

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.

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.

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.

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 (m6A) 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 (m5C) has been found in mRNA, rRNA and tRNA of representative organisms from all kinds of species. As reversible epigenetic modifications, m5C 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 m5C in tumorigenesis. Here, we review the latest progresses in the biological roles of m5C 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.

Identification of the Expression Patterns and Potential Prognostic Role of 5-Methylcytosine Regulators in Hepatocellular Carcinoma

Background: Hepatocellular carcinoma (HCC) is the most common primary liver cancer with a poor prognosis. 5-methylcytosine (m5C) modification plays a nonnegligible role in tumor pathogenesis and progression. However, little is known about the role of m5C regulators in HCC.

Methods: Based on 9 m5C regulators, the m5C modification patterns of HCC samples extracted from public databases were systematically evaluated and correlated with tumor immune and prognosis characteristics. An integrated model called the “m5Cscore” was constructed using principal component analysis, and its prognostic value was evaluated.

Results: Almost all m5C regulators were differentially expressed between HCC and normal tissues. Through unsupervised clustering, three different m5Cclusters were ultimately uncovered; these clusters were characterized by differences in prognosis, immune cell infiltration, and pathway signatures. The m5Cscore was constructed to quantify the m5C modifications of individual patients. Subsequent analysis revealed that the m5Cscore was an independent prognostic factor of HCC and could be a novel indicator to predict the prognosis of HCC.

Conclusion: This study comprehensively explored and systematically profiled the features of m5C modification in HCC. m5C modification patterns play a crucial role in the tumor immune microenvironment (TIME) and prognosis of HCC. The m5Cscore provides a more holistic understanding of m5C modification in HCC and provides a practical tool for predicting the prognosis of HCC. This study will help clinicians identify effective indicators of HCC to improve the poor prognosis of this disease.

Analysis of Multi-Layer RNA Modification Patterns for the Characterization of Tumor Immune Microenvironment in Hepatocellular Carcinoma

Background: RNA modifications have emerged as important posttranscriptional changes in multiple tumor cellular processes and tumorigenesis, including hepatocellular carcinoma (HCC). However, the potential roles and the interaction between regulators of RNA modifications and the tumor microenvironment (TME) are unclear in HCC. Methods: The gene expression profiles of 26 RNA modification "writers" were investigated in the TCGA cohort. The unsupervised clustering approach was used to class these RNA modification regulators. The characteristics of immune cell infiltration from TME for each cluster was tested by the CIBERSORT method. Additionally, we established a scoring model to evaluate the RNA modification characteristics of individual tumors. The associations between the scoring model and genetic as well as clinical characteristics, drug sensitivity, and response to immunotherapy were also analyzed. Results: We mapped the somatic mutations and somatic copy number variation of the RNA modification regulators. The expression of all selected regulators was detected, and two modification patterns were identified that featured distinct immune cell infiltration characteristics. Subsequently, we developed a score model (termed as WM-Score model). Furthermore, the survival analysis showed that the WM-Score value was associated with HCC patient prognosis. The results of the ROC curves analysis and multivariate analysis all confirmed that the WM-Score value was strongly associated with anti-cancer drug resistance and therapeutic efficacy of immunotherapy, thus could be used as an independent risk factor in HCC. Conclusion: Our research identified two RNA modification patterns characterized by distinct TME, and the WM-Score model was developed that might serve as reliable prognostic and immunotherapeutic effect predictor of HCC.