- 1Department of Oncology, Tongji Hospital, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China
- 2Department of Epidemiology, Center for Global Health, School of Public Health, Nanjing Medical University, Nanjing, China
- 3Department of Medical Oncology, National Cancer Center, National Clinical Research Center for Cancer, Cancer Hospital, Chinese Academy of Medical Sciences and Peking Union Medical College, Beijing, China
RNA modification plays important role in the occurrence and development of hepatocellular carcinoma. The best characterized RNA modification is m6A, while other kinds of RNA modifications have not been fully investigated in hepatocellular carcinoma (HCC). In the current study, we investigated the roles of one hundred RNA modification regulators belonging to eight different types of cancer-related RNA modifications in HCC. Expression analysis revealed that nearly 90% RNA regulators exhibited significantly higher expression in tumors than normal tissues. By consensus clustering, we identified two clusters with distinct biological characteristics, immune microenvironment, and prognostic pattern. An RNA modification score (RMScore) was constructed and stratified patients into high- and low-risk group, which showed significantly different prognosis. Moreover, a nomogram including clinicopathologic features and the RMScore could well predict the survival in HCC patients. This study indicated the important role of eight types of RNA modification in HCC and develop a RMScore, which will be a new method to forecast the prognosis of HCC patients.
Introduction
Hepatocellular carcinoma (HCC) is one of the leading cause of cancer-related deaths worldwide, with few effective therapeutic options (Sung et al., 2021). Most patients have already progressed to the middle and late stages when they are diagnosed, leading to a poor prognosis of HCC patients (Sia et al., 2017). The complicated etiologies and complex pathogenesis lead to the high heterogeneity of HCC, making the efficiency of treatment very low (Park et al., 2010; Seehawer et al., 2018). Identifying novel treatment and prognostic targets for HCC is of urgent need.
RNA modification is a critical posttranscriptional regulators of gene expression program (Nombela et al., 2021). With the rapid development of molecular and sequencing technologies, RNA modification has become a hot research topic in recent years. Emerging reports confirm that dysregulation of RNA modification gives rise to a variety of human diseases, particularly hepatocellular carcinoma (Xu et al., 2021). Various types of cancer-related RNA modifications can be found in eukaryotes, including N6-methyladenosine (m6A), N1-methyladenosine (m1A), and 5-methylcytosine (m5C), 2′-O-methylation (Nm), m7G N7-methylguanosine (m7G), pseudouridylation (Ψ), adenosine-to-inosine (A-to-I), 5-methoxycarbonylmethyl-2-thiouridine (mcm5s2U) (Boccaletto et al., 2018; Frye et al., 2018). Among them, m6A is the most investigated and best characterized one (Shi et al., 2019; Wei and He, 2021). For example, KIAA1429 was demonstrated to promote HCC invasion and migration by altering the methylation of m6A in ID2 and GATA3 mRNA (Cheng et al., 2019; Lan et al., 2019). ALKBH5 suppresses proliferation and invasion of HCC via m6A-guided epigenetic inhibition of LYPD1 (Chen et al., 2020). METTL3 promotes HCC progression via post-transcriptional silencing of SOCS2 (Chen et al., 2018). Other methylations are also reported to be related to HCC. For example, aberrant NSUN2-mediated m5C modification of H19 lncRNA is proved to be associated with poor differentiation of HCC (Sun et al., 2020). TRMT6/TRMT61A-mediated m1A methylation is required for self-renewal of liver cancer stem cells and tumorigenesis (Wang et al., 2021a). However, these studies mostly focused on the functions of only one RNA modification regulator. The crosstalk between different regulators remains unknown. Recently, the interactions between different regulators within one type of RNA modification have been investigated. For example, Fang et al. reported a two m6A gene-based signature (HNRNPA2B1 and RBM15) could predict the prognosis of HBV-related HCC (Fang and Chen, 2020). A risk signature consisted of four m1A regulators (TRMT6, TRMT61A, TRMT10C, and YTHDF1) was remarkably associated with HCC patient prognosis (Shi et al., 2020). A risk signature comprised of six m6A-related genes was reported to be able to predict the prognosis of HCC(19). However, these studies ignored the relationship and crosstalk between different types of RNA regulators. It is very necessary to systematically characterize the relationship and crosstalk of other kinds of RNA modification regulators in HCC.
In the current research, we profiled one hundred RNA regulators from eight types of RNA modifications in HCC for the first time. We clustered HCC patients into two RNA modification-related subtypes with different prognosis and cancer hallmarks. An RNA modification score was constructed to serve as an independent prognosis predictor for HCC patients. These results would provide insight into the prognosis prediction and therapeutic management of HCC in the future.
Materials and methods
Data collection
The gene mutation, expression and clinical information of HCC patients was downloaded from the TCGA and ICGC database.
There were huge number of RNA regulators identified in eukaryotes. However, we only choose RNA regulators which were proved to play critical roles in HCC and other cancers (Delaunay and Frye, 2019; Barbieri and Kouzarides, 2020; Han et al., 2021). As a result, a total of one hundred RNA regulators which belong to eight types of RNA modifications (m6A, m1A, m5C, Nm, m7G, Ψ, A-to-I, and mcm5s2U) were included in our study. They were listed in Supplementary Table S1.
Expression and mutation analysis of RNA modification regulators
Unpaired t-test were used to compare the expression of 100 RNA regulators between tumor and adjacent normal tissues. The “Maftools” R package was used to analyze the mutation and interaction of 100 RNA modification genes. The protein-protein interaction between 100 RNA modification regulators were analyzed by the STRING database and visualized by using Cytoscape software (Szklarczyk et al., 2019).
Identification of RNA modification regulator-related subgroups
To figure out the relationship between the expression of 100 RNA modification regulator and HCC subtypes, consensus clustering analysis was conducted in TCGA-LIHC patients by using the “ConsensusClusterPlus” R package.
Pathway enrichment analysis
To reveal the biological function enriched in different group of patients, the Gene Set Variation Analysis (GSVA) was performed by utilizing the R package “GSVA”. The annotated gene set “h.all.v7.4.symbols.gmt” was used for GSVA analysis. The enrichment scores of immune cells in sample from different groups were determined with single-sample gene set enrichment analysis (ssGSEA).
Construction and validation of the RNA modification score (RMScore)
To construct the RMScore, we first identified 1,427 differentially expressed genes (DEGs) in cluster 1 and cluster 2 groups. We divided TCGA-LIHC into training set and validation set. In the training set, univariate cox regression analysis revealed that 266 DEGs were identified to be significantly related to the prognosis of HCC patients (p < 0.001). Then LASSO regression analysis was conducted to select survival-associated DEGs. Multivariate cox regression analysis was conducted to construct the final prognostic model. The RMScore was calculated by the following equation: RMScore = Exp1* Coe 1 + Exp2* Coe 2 + … Exp N* Coe N.
The validation dataset of TCGA-LIHC, whole TCGA-LIHC dataset, and ICGC dataset were utilized as external validation dataset. The RMScore of each sample was calculated by the formula presented above. Patients were divided into low RMScore and high RMScore groups by the median value of the RMScore. The “Survival” R package was conducted to analysis the survival of two groups. The “SurvivalROC” R package was used to plot the ROC curves.
Chemotherapy response prediction
To predict the sensitivity to chemotherapeutic drugs, we use the R package “pRRophetic” to estimate the half-maximal inhibitory concentration (IC50) in individual patients. The R package “ggplot2” was utilized to generate the plots.
Construction of a nomogram
The stage and RMScore, both of which were identified to be an independent prognostic value by multivariate Cox proportional hazards analysis, were included to construct the nomogram. The calibration curve was plotted to observe the prediction probabilities of the nomogram against the observed rates. The prognostic ability of the nomogram was evaluated by receiver operating (ROC) curves and decision curve analysis (DCA) curves.
Results
Dysregulation and mutation of eight types of RNA regulators in HCC
We choose a total of one hundred RNA modification regulators which were reported to play roles in at least one kind of cancer to perform following analysis (Supplementary Table S1).
In order to characterize the expression profile of 100 RNA regulators, unpaired t-test was conducted to compare the gene expression between tumor and normal tissues. Notably, nearly 90% RNA regulators exhibited significantly higher expression in tumors than normal tissues (Figure 1A). Among the eight types of RNA modification regulators, all regulators in mcm5s2U, Ψ and A-to-I showed elevated expression in tumor tissues. The dysregulation of these RNA regulators may play key roles in the development of HCC.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1127301/fgene-14-1127301-HTML/image_m/fgene-14-1127301-g001.jpg)
FIGURE 1. The dysregulation and mutation landscape of 100 RNA modification regulators included in our study in HCC. (A) The gene expression level of 100 RNA modification regulators in normal and tumor tissues in TCGA-LIHC. (B) Genomic mutational landscape of RNA modification regulators with mutation rates ≥1% in the TCGA-LIHC cohort. (C) Correlation of the mutation between genes with mutation rates ≥1%. (D) The protein-protein interactions between one hundred RNA regulators analyzed by the STRING database. (p-value < 0.05 and correlation coefficient ≥0.6). (ns, not significant; *, p < 0.05: **, p < 0.01; ***, p < 0.001).
The mutation profiles of 100 hundred regulators in HCC were depicted together. As a results, 94 of 371 samples (about 25.34%) showed RNA modification regulator mutations. M5C eraser TET1 and A-to-I writer ADARB2 had a mutation rate of 2%. About 40% RNA regulators had a mutation rate of 1% (Figure 1B). Genetic interaction analysis showed that mutation of m5C eraser TET3 was positively correlated with mutations of TRMT10B, ADAT2, PRRC2A, YTHDC1 and TET2. Mutation of m6A reader HNRNPA2B1 was positively correlated with mutations of WTAP, ELAVL1, MRM2, HNRNPC, HNRNPA1, YTHDC1, PRRC2A. Mutation of m6A writer ZC3H13 was positively correlated with mutations of TRMT10B, TET3, TRMT13, PRRC2A, YTHDC1, ADAR, YTHDC2 and ZCCHC4, while negatively correlated with mutations of METTL1 and CTU2 (Figure 1C).
To reveal the interactive relationships between proteins encoded by 100 RNA regulators, we constructed a PPI network using STRING and Cytoscape software. Correlations with adjusted p-value <0.05 and |r| ≥ 0.6 were shown in Figure 1D. We can see that these RNA regulators, especially regulators of m6A and m5C, showed close correlations with each other.
Results above collectively demonstrated that there were huge interaction and crosstalk between these 100 RNA regulators. They may form a complicated network to synergistically mediate the development and progression of HCC.
Identification of different RNA modification patterns in HCC
In order to excavate the RNA modification patterns of HCC, consensus clustering analysis were conducted in 374 patients in TCGA-LIHC dataset. Patients in TCGA-LIHC dataset were separated into diverse clusters (K = 2–9) according to the expression of 100 RNA modification genes via an unsupervised consensus clustering analysis. The cumulative distribution function (CDF) curve and the area under the CDF curve was plotted, which indicated that the cohort was well distributed when K = 2 (Figures 2A–C). Kaplan-Meier survival analysis showed that patients in cluster 1 had a significantly better prognosis than patients in cluster 2 (p < 00.001) (Figure 2D). GSVA and ssGSEA enrichment analysis were performed to analyze biological process in two clusters. As shown in Figure 2E, most signaling pathways including DNA repair, MYC targets, G2M checkpoint, E2F targets were significantly enriched in Cluster 2. These may lead to the high degree of malignancy and poor prognosis of tumors in Cluster 2. Signaling pathways such as oxidative phosphorylation, reactive oxygen species pathway, cholesterol homeostasis and fatty acid metabolism were enriched in Cluster 1. There were more NK cells, eosinophils, neutrophils, effector memory CD8 T cells in Cluster 1, while more activated CD4 T cells, T follicular helper cells, type 2 T helper cell and effector memory CD4 T cells were accumulated in Cluster 2. This result implied two clusters had different immune infiltrating patterns.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1127301/fgene-14-1127301-HTML/image_m/fgene-14-1127301-g002.jpg)
FIGURE 2. Identification of consensus clusters by eight types of RNA modification regulators in LIHC. (A) Cumulative distribution function (CDF) for k = 2–9. (B) The area under the CDF curve k = 2–9. (C) Consensus clustering matrix for k = 2. (D) Kaplan–Meier survival curve of two clusters. (E) The heatmap of the different cancer hallmarks enriched in two clusters. (F) The difference of the immune cells infiltrated in tumor microenvironments between two clusters.
Results above collectively suggested that 100 RNA modification regulators defined two RNA modification patterns in HCC, which exhibited different malignancy, immune microenvironment, and prognosis.
Establishment and verification of RNA modification score in HCC
To further investigate the potential genetic alterations between two RNA modification patterns, differential expression analyses were carried out between two clusters (Figure 3A). A total of 1,472 differentially expressed genes (DEGs) were identified (Supplementary Table S2). Next, we divided the TCGA-LIHC queue into a training cohort and a verification cohort. Among the 1,427 DEGs, 266 genes were associated with the prognosis of HCC patients in training cohort by univariate COX regression analysis (Supplementary Table S3). Then, the LASSO regression analysis was conducted to further narrow the survival-related DEGs (Supplementary Table S4). Finally, six genes were identified by stepwise multivariate regression analysis and subsequently used to construct the RNA modification score (RMScore). The six genes identified were CYP2C9, CBX2, NDRG1, EPS8L3, FAM83D and MYCN. The RMScore = 0.1034*Exp of CBX2 + 0.0087*Exp of NDRG1 + 0.0392*Exp of EPS8L3 + 0.0517*Exp of FAM83D + 0.0596*Exp of MYCN−0.0031*Exp of CYP2C9.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1127301/fgene-14-1127301-HTML/image_m/fgene-14-1127301-g003.jpg)
FIGURE 3. Development and validation of RNA modification related signature for HCC. (A) The workflow of the construction of the RMScore. (B) Survival curves stratified by the RMScore in the TCGA training cohort. (C) Receiver operating characteristic (ROC) curves of RMScore for predicting survival in the TGGA training cohort. (D) Survival curves stratified by the RMScore in the TCGA validation cohort. (E) ROC curves of RMScore for predicting survival in the TGGA validation cohort. (F) Survival curves stratified by the RMScore in the ICGC cohort. (G) ROC curves of RMScore for predicting survival in the ICGC cohort.
By Kaplan-Meier analysis, we found that patients in low-risk group had a clear survival advantage over the high-risk group in training set (p < 0.001) (Figure 3B). The same survival advantage was also observed in the validation cohort (p = 0.033) and the ICGC cohort (p < 0.001) (Figures 3D, F). The time-dependent ROC curves were utilized in assessing the performance of RMScore, which showed the RMScore had good sensitivity and specificity in prognosis prediction (Figures 3C, E, G). These results suggested that the RMScore had good prediction ability for the survival of HCC patients.
Independent prognostic value of RMScore
The RMScore together with clinicopathological factors, such as age, gender, grade, and stage were included for univariate and multivariate Cox regression analysis. In the univariate analysis, the RMScore was significantly related to the prognosis of HCC patients (hazard ratio [HR], 1.053; 95% confidence interval [CI], 1.037–1.069; p < 0.001) (Figure 4A). Multivariate analysis also demonstrated that the RMScore was an independent risk of the prognosis of HCC patients (hazard ratio [HR], 1.045; 95% confidence interval [CI], 1.028–1.062; p < 0.001) (Figure 4B). In the heatmap incorporating gene expression and clinicopathological traits, the majority of patients in Cluster 1 were gathered in low RMScore group, which demonstrated that the RMScore well represented the traits of patients in Cluster 1. Among the six gene which constructing the RMScore, only CYP2C9 was upregulated in the Cluster 1 as well as low RMScore group. The other five genes (CBX2, NDRG1, EPSBL3, FAM83D and MYCN) were upregulated in Cluster 2 as well as high RMScore group (Figure 4C). We explored the pathway enrichment associated with the RMScore to identify the internal mechanisms involved. Similar to the pathway enrichment results between Cluster 1 and Cluster 2, the high RMScore group was mainly enriched in signaling pathways related to DNA repair, MYC targets, G2M checkpoint, and E2F targets, while the low RMScore group was enriched in oxidative phosphorylation, reactive oxygen species pathway, cholesterol homeostasis and fatty acid metabolism (Figure 4D).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1127301/fgene-14-1127301-HTML/image_m/fgene-14-1127301-g004.jpg)
FIGURE 4. (A) Univariate Cox regression analysis and (B) multivariate Cox regression analysis verifying independent prognostic value of RMScore for the prognosis of HCC patients. (C) Heat map for the association of gene and clinicopathologic features with prognosis of HCC patients. (D) Visualization of biological processes analyzed by GSVA in two RMScore subgroups.
Potential predictive biomarker for chemotherapy and target therapy
We investigated the association between the RMScore and the sensitivity to chemotherapeutics in HCC patients. As a result, significant differences in the estimated IC50 between the two risk groups were observed for common chemotherapy or targeted drugs including paclitaxel (Figure 5A), sorafenib (Figure 5B), 5-flurouracil (Figure 5C), doxorubicin (Figure 5D), gemcitabine (Figure 5E) and erlotinib (Figure 5F). Lower RMScore was related to higher IC50 among paclitaxel, sorafenib, 5-flurouracil, gemcitabine, suggesting that patients with low RMScore were less sensitive to these treatments. However, patients with low RMScore had a lower IC50 to erlotinib, suggesting that they may be more sensitive to erlotinib than high RMScore patients.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1127301/fgene-14-1127301-HTML/image_m/fgene-14-1127301-g005.jpg)
FIGURE 5. Estimated IC50 for paclitaxel (A), sorafenib (B), 5-flurouracil (C), doxorubicin (D), gemcitabine (E) and erlotinib (F) between high- and low-RMScore group.
Development of a predictive nomogram
To facilitate the clinical applicability and availability of the RMScore, a predictive nomogram for 1-, 3-, 5-year OS combined with tumor stage and RMScore was developed (Figure 6A). The calibration curves showed that actual overall survival of patients was almost consistent with the predicted overall survival of patients (Figure 6B). The DCA plot suggested that the nomogram had a superior efficiency than RMScore in predicting patient survival outcomes (Figure 6C). The AUCs of the nomogram in the 1-, 3-, and 5-year ROC curves were 0.750, 0.743 and 0.735, respectively (Figure 6D).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1127301/fgene-14-1127301-HTML/image_m/fgene-14-1127301-g006.jpg)
FIGURE 6. Construction of a nomogram by clinicopathological factors and RMScore. (A) Prognostic nomogram based on RMScore and tumor stage. (B) The 1-, 3-, and 5-year calibration plots of the nomogram. (C) DCA curves of the nomogram and other clinicopathological parameters. (D) The 1-, 3-, and 5-year ROC curves of the nomogram.
Discussion
The current interventions for HCC are not satisfactory, and more precise prognostic indicator and promising strategies need to be explored. Recent studies point out that RNA modifiers play vital roles in HCC cell proliferation, exacerbation, and metastasis (Zhu et al., 2021). However, these studies are mainly limited to several RNA modification regulators. The function of other kinds of RNA modification regulators and the interaction with each other remain to be elucidated.
Here, we systematically analyzed the expression, mutation, and the prognostic value of a total of one hundred RNA modifiers from eight types of RNA modifications in HCC. We found that most of the RNA modifiers were highly expressed in the tumor tissues compared to normal tissues. At the genetic level, nearly half of genes showed a least one kind of mutation. Most of the mutations were positively interacted with each other. These results suggested that eight types of RNA modifications might play important role in HCC.
Two RNA modification patterns with distinct prognosis were identified in HCC patients.
Then, we constructed RMScore to calculate the modification status of HCC patient. Patients with low RMScore showed better prognosis than patients with high RMScore. Similar to the pathway enrichment results in Cluster 1, pathway related to the cell cycle and proliferation, such as “DNA repair,” “G2M checkpoint,” “E2F targets” were all significantly enriched in the low RMScore group. All these results implied that the RMScore could well represent the modification status of the HCC patients.
There were some studies which documented different RNA modification patterns in various type of tumors, including HCC. For example, Qi et al. identified two distinct m6A modification patterns, one of which showed higher expression of m6A modification regulators and poor prognosis (Qi et al., 2020). Similarly, Li et al. reported two distinct clusters with different prognosis and clinical features based on the expression of 45 m6A/m5C/m1A regulated genes (Li et al., 2022). In addition, RNA modification patterns, especially m6A modification-related patterns, are also reported to have prognostic predictive value in lung cancer (Liu et al., 2020; Zhang et al., 2021), gastric cancer (Zhang et al., 2020), prostate cancer (Liu et al., 2021), kidney renal clear cell carcinoma (Li et al., 2021) and pancreatic adenocarcinoma (Wang et al., 2021b). In the current study, a total of one hundred cancer-related RNA regulators from eight-types of RNA modifications were included. We identified two subgroups with distinct cancer hallmarks, tumor microenvironments, and overall survival. These results showed the importance of other types of RNA modifications in HCC.
Finding optimal strategy to select patients most likely to benefit from available chemotherapeutic regimens can facilitate therapeutic decision-making. Recently, RNA modification was proved to be associated with sensitivity to anticancer therapies (Lan et al., 2021). For example, m6A writer METTL3 was documented to regulate sorafenib resistance though FOXO3-mediated autophagy in HCC (Lin et al., 2020). M7G writer WDR4 was demonstrated to play a role in sorafenib resistance by inducing CCNB1 translation in HCC. Herein, we used GDSC database to assess the differences in the sensitivity to various chemotherapeutics between the two RMScore groups. In our study, patients with high RMScore had a lower IC50, implying to be more sensitive to chemotherapy drugs including paclitaxel, sorafenib, 5-flurouracil, doxorubicin, mitomycin C and gemcitabine. The RMScore could be served as a good measure to predict the drug efficacy for clinicians.
Several limitations for our study should be recognized. Firstly, there were several cancer-related RNA modification regulators, such as m6Am regulators PCIF1 and METTL4, were not included in our analysis. This may lead to subtle deviations in our results, but does not affect the main conclusions. Secondly, due to the limited sample size, we still need more patients to verify the reliability of our conclusions. Thirdly, the potential mechanism of the role of genes contained in the RMScore was not based on experimental evidence, and required further in vivo and in vitro experiment. Last but not the least, as most of HCC patients were Caucasian (most were Hepatitis C virus-related HCC), whether the RMScore is also applicable to non-Caucasian is unknown. However, our study was the first study to give an insight into the mutation and expression of eight types of RNA modifications in HCC. This study paves the way for further investigations of the roles of eight types of RNA modifications in the HCC in the future.
Data availability statement
Publicly available datasets were analyzed in this study. This data can be found here: https://portal.gdc.cancer.gov; https://dcc.icgc.org.
Author contributions
Conception and design: WQ and JZ; Acquisition of data: CJ; Analysis of data: WQ and CJ; Writing the manuscript: WQ and JZ; All authors read and approved the final manuscript.
Funding
This work was supported by the National Natural Science Foundation of China for Young Scholars (Grant No.82203972).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2023.1127301/full#supplementary-material
References
Barbieri, I., and Kouzarides, T. (2020). Role of RNA modifications in cancer. Nat. Rev. Cancer 20 (6), 303–322. doi:10.1038/s41568-020-0253-2
Boccaletto, P., Machnicka, M. A., Purta, E., Piatkowski, P., Baginski, B., Wirecki, T. K., et al. (2018). Modomics: A database of RNA modification pathways. 2017 update. Nucleic Acids Res. 46 (D1), D303–D307. doi:10.1093/nar/gkx1030
Chen, M., Wei, L., Law, C. T., Tsang, F. H., Shen, J., Cheng, C. L., et al. (2018). RNA N6-methyladenosine methyltransferase-like 3 promotes liver cancer progression through YTHDF2-dependent posttranscriptional silencing of SOCS2. Hepatology 67 (6), 2254–2270. doi:10.1002/hep.29683
Chen, Y., Zhao, Y., Chen, J., Peng, C., Zhang, Y., Tong, R., et al. (2020). ALKBH5 suppresses malignancy of hepatocellular carcinoma via m(6)A-guided epigenetic inhibition of LYPD1. Mol. Cancer 19 (1), 123. doi:10.1186/s12943-020-01239-w
Cheng, X., Li, M., Rao, X., Zhang, W., Li, X., Wang, L., et al. (2019). KIAA1429 regulates the migration and invasion of hepatocellular carcinoma by altering m6A modification of ID2 mRNA. Onco Targets Ther. 12, 3421–3428. doi:10.2147/OTT.S180954
Delaunay, S., and Frye, M. (2019). RNA modifications regulating cell fate in cancer. Nat. Cell Biol. 21 (5), 552–559. doi:10.1038/s41556-019-0319-0
Fang, Q., and Chen, H. (2020). The significance of m6A RNA methylation regulators in predicting the prognosis and clinical course of HBV-related hepatocellular carcinoma. Mol. Med. 26(1), 60. doi:10.1186/s10020-020-00185-z
Frye, M., Harada, B. T., Behm, M., and He, C. (2018). RNA modifications modulate gene expression during development. Science 361 (6409), 1346–1349. doi:10.1126/science.aau1646
Han, X., Wang, M., Zhao, Y. L., Yang, Y., and Yang, Y. G. (2021). RNA methylations in human cancers. Semin. Cancer Biol. 75, 97–115. doi:10.1016/j.semcancer.2020.11.007
Lan, Q., Liu, P. Y., Bell, J. L., Wang, J. Y., Huttelmaier, S., Zhang, X. D., et al. (2021). The emerging roles of RNA m(6)A methylation and demethylation as critical regulators of tumorigenesis, drug sensitivity, and resistance. Cancer Res. 81 (13), 3431, doi:10.1158/0008-5472.CAN-20-4107
Lan, T., Li, H., Zhang, D., Xu, L., Liu, H., Hao, X., et al. (2019). KIAA1429 contributes to liver cancer progression through N6-methyladenosine-dependent post-transcriptional modification of GATA3. Mol. Cancer 18 (1), 186. doi:10.1186/s12943-019-1106-z
Li, D., Li, K., Zhang, W., Yang, K. W., Mu, D. A., Jiang, G. J., et al. (2022). The m6A/m5C/m1A regulated gene signature predicts the prognosis and correlates with the immune status of hepatocellular carcinoma. Front. Immunol. 13, 918140. doi:10.3389/fimmu.2022.918140
Li, H., Hu, J., Yu, A., Othmane, B., Guo, T., Liu, J., et al. (2021). RNA modification of N6-methyladenosine predicts immune phenotypes and therapeutic opportunities in kidney renal clear cell carcinoma. Front. Oncol. 11, 642159, doi:10.3389/fonc.2021.642159
Lin, Z., Niu, Y., Wan, A., Chen, D., Liang, H., Chen, X., et al. (2020). RNA m(6) A methylation regulates sorafenib resistance in liver cancer through FOXO3-mediated autophagy. EMBO J. 39 (12), e103181. doi:10.15252/embj.2019103181
Liu, Y., Guo, X., Zhao, M., Ao, H., Leng, X., Liu, M., et al. (2020). Contributions and prognostic values of m(6) A RNA methylation regulators in non-small-cell lung cancer. J. Cell Physiol. 235 (9), 6043–6057. doi:10.1002/jcp.29531
Liu, Z., Zhong, J., Zeng, J., Duan, X., Lu, J., Sun, X., et al. (2021). Characterization of the m6A-associated tumor immune microenvironment in prostate cancer to aid immunotherapy. Front. Immunol. 12, 735170. doi:10.3389/fimmu.2021.735170
Nombela, P., Miguel-Lopez, B., and Blanco, S. (2021). The role of m(6)A, m(5)C and Psi RNA modifications in cancer: Novel therapeutic opportunities. Mol. Cancer 20 (1), 18. doi:10.1186/s12943-020-01263-w
Park, E. J., Lee, J. H., Yu, G. Y., He, G., Ali, S. R., Holzer, R. G., et al. (2010). Dietary and genetic obesity promote liver inflammation and tumorigenesis by enhancing IL-6 and TNF expression. Cell 140 (2), 197–208. doi:10.1016/j.cell.2009.12.052
Qi, L. W., Jia, J. H., Jiang, C. H., and Hu, J. M. (2020). Contributions and prognostic values of N6-methyladenosine RNA methylation regulators in hepatocellular carcinoma. Front. Genet. 11, 614566. doi:10.3389/fgene.2020.614566
Seehawer, M., Heinzmann, F., D'Artista, L., Harbig, J., Roux, P. F., Hoenicke, L., et al. (2018). Necroptosis microenvironment directs lineage commitment in liver cancer. Nature 562 (7725), 69–75. doi:10.1038/s41586-018-0519-y
Shi, H., Wei, J., and He, C. (2019). Where, when, and how: Context-dependent functions of RNA methylation writers, readers, and erasers. Mol. Cell 74 (4), 640–650. doi:10.1016/j.molcel.2019.04.025
Shi, Q., Xue, C., Yuan, X., He, Y., and Yu, Z. (2020). Gene signatures and prognostic values of m1A-related regulatory genes in hepatocellular carcinoma. Sci. Rep. 10 (1), 15083. doi:10.1038/s41598-020-72178-1
Sia, D., Villanueva, A., Friedman, S. L., and Llovet, J. M. (2017). Liver cancer cell of origin, molecular class, and effects on patient prognosis. Gastroenterology 152 (4), 745–761. doi:10.1053/j.gastro.2016.11.048
Sun, Z., Xue, S., Zhang, M., Xu, H., Hu, X., Chen, S., et al. (2020). Aberrant NSUN2-mediated m(5)C modification of H19 lncRNA is associated with poor differentiation of hepatocellular carcinoma. Oncogene 39 (45), 6906–6919. doi:10.1038/s41388-020-01475-w
Sung, H., Ferlay, J., Siegel, R. L., Laversanne, M., Soerjomataram, I., Jemal, A., et al. (2021). Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J. Clin. 71, 209. doi:10.3322/caac.21660
Szklarczyk, D., Gable, A. L., Lyon, D., Junge, A., Wyder, S., Huerta-Cepas, J., et al. (2019). STRING v11: Protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 47 (D1), D607–D613. doi:10.1093/nar/gky1131
Wang, L., Zhang, S., Li, H., Xu, Y., Wu, Q., Shen, J., et al. (2021). Quantification of m6A RNA methylation modulators pattern was a potential biomarker for prognosis and associated with tumor immune microenvironment of pancreatic adenocarcinoma. BMC Cancer 21 (1), 876. doi:10.1186/s12885-021-08550-9
Wang, Y., Wang, J., Li, X., Xiong, X., Wang, J., Zhou, Z., et al. (2021). N(1)-methyladenosine methylation in tRNA drives liver tumourigenesis by regulating cholesterol metabolism. Nat. Commun. 12(1), 6314, doi:10.1038/s41467-021-26718-6
Wei, J., and He, C. (2021). Chromatin and transcriptional regulation by reversible RNA methylation. Curr. Opin. Cell Biol. 70, 109–115. doi:10.1016/j.ceb.2020.11.005
Xu, Y., Zhang, M., Zhang, Q., Yu, X., Sun, Z., He, Y., et al. (2021). Role of main RNA methylation in hepatocellular carcinoma: N6-Methyladenosine, 5-methylcytosine, and N1-methyladenosine. Front. Cell Dev. Biol. 9, 767668. doi:10.3389/fcell.2021.767668
Zhang, B., Wu, Q., Li, B., Wang, D., Wang, L., and Zhou, Y. L. (2020). m(6 A regulator-mediated methylation modification patterns and tumor microenvironment infiltration characterization in gastric cancer. Mol. Cancer 19 (1), 53. doi:10.1186/s12943-020-01170-0
Zhang, H., Hu, J., Liu, A., Qu, H., Jiang, F., Wang, C., et al. (2021). An N6-methyladenosine-related gene set variation score as a prognostic tool for lung adenocarcinoma. Front. Cell Dev. Biol. 9, 651575. doi:10.3389/fcell.2021.651575
Keywords: RNA modification, hepatocellular carcinoma, prognosis, nomogram, therapeutic responses
Citation: Qin W, Jin C and Zou J (2023) Analysis of eight types of RNA modification regulators and their correlation with the prognosis in hepatocellular carcinoma. Front. Genet. 14:1127301. doi: 10.3389/fgene.2023.1127301
Received: 19 December 2022; Accepted: 20 February 2023;
Published: 16 March 2023.
Edited by:
Iben Lyskjær, Aarhus University, DenmarkReviewed by:
Jiangbo Wei, The University of Chicago, United StatesPeng Du, Peking University, China
Mark Jon Olsen, Midwestern University, United States
Copyright © 2023 Qin, Jin and Zou. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jun Zou, 263390138@qq.com; Chen Jin, jinchen0223@126.com