Parbatraj Regmi1†
Fu-Yu Li- 1Department of Biliary Surgery, West China Hospital, Sichuan University, Chengdu, China
- 2Department of Uro Surgery, West China Hospital, Sichuan University, Chengdu, China
- 3Department of Dermatology, West China Hospital, Sichuan University, Chengdu, China
Background: About 90% of liver cancer-related deaths are caused by hepatocellular carcinoma (HCC). N7-methylguanosine (m7G) modification is associated with the biological process and regulation of various diseases. To the best of our knowledge, its role in the pathogenesis and prognosis of HCC has not been thoroughly investigated.
Aim: To identify N7-methylguanosine (m7G) related prognostic biomarkers in HCC. Furthermore, we also studied the association of m7G–related prognostic gene signature with immune infiltration in HCC.
Methods: The TCGA datasets were used as a training and GEO dataset “GSE76427” for validation of the results. Statistical analyses were performed using the R statistical software version 4.1.2.
Results: Functional enrichment analysis identified some pathogenesis related to HCC. We identified 3 m7G-related genes (CDK1, ANO1, and PDGFRA) as prognostic biomarkers for HCC. A risk score was calculated from these 3 prognostic m7G-related genes which showed the high-risk group had a significantly poorer prognosis than the low-risk group in both training and validation datasets. The 3- and 5-years overall survival was predicted better with the risk score than the ideal model in the entire cohort in the predictive nomogram. Furthermore, immune checkpoint genes like CTLA4, HAVCR2, LAG3, and TIGT were expressed significantly higher in the high-risk group and the chemotherapy sensitivity analysis showed that the high-risk groups were responsive to sorafenib treatment.
Conclusion: These 3 m7G genes related signature model can be used as prognostic biomarkers in HCC and a guide for immunotherapy and chemotherapy response. Future clinical study on this biomarker model is required to verify its clinical implications.
Introduction
Liver Cancer is the major cause of cancer-related mortality. Hepatocellular carcinoma (HCC) accounts for approximately 90% of liver cancer-related deaths (Global Burden of Disease Cancer et al., 2017; Bray et al., 2018). Despite the improvement in the surveillance program, most of the patients are diagnosed at an advanced stage with unresectable disease (Bruix et al., 2016; Kudo et al., 2016). In recent years, surgical resection, ablation therapies, transarterial chemoembolization, immunotherapy, and systemic chemotherapeutic agents are used in the treatment of early to advanced staged tumors (Llovet et al., 2021). Furthermore, recent advance in biomedical research has led to the identification of molecular pathogenesis, diagnostic biomarkers, and prognostic biomarkers. However, there are no established biomarkers for predicting the response to therapy which could help the clinicians to guide the appropriate treatment strategy for specific patients (Rizzo and Brandi, 2021; Rizzo and Ricci, 2021).
N7-methylguanosine (m7G), methylation on the N7 position of guanosine is a type of positively charged RNA modification in tRNA, rRNA, or internal mRNA regions (Guy and Phizicky, 2014; Sloan et al., 2017; Malbec et al., 2019). Previous studies have shown that more than 100 genes are associated with methylation in several m7G sites (Song et al., 2020; Xia et al., 2021). N7-methylguanosine (m7G) modification is associated with the biological process and regulation of various diseases (Shaheen et al., 2015; Lin et al., 2018; Ma et al., 2021). Furthermore, few recent studies have shown that methylation and mutation at several m7G sites are associated with the molecular pathogenesis of cancer (Chen et al., 2021; Ma et al., 2021). Therefore identification of the m7G-related genes will help in understanding the molecular pathogenesis of cancer and its role in diagnosis, treatment, and prediction of prognosis. The role of m7G-related genes in the pathogenesis and prognosis of HCC, however, has not been studied in greater detail.
Herein, we perform a bioinformatics analysis to identify the association of several m7G-related genes in the molecular pathogenesis of HCC. Furthermore, we aim to construct a prediction model and nomogram to predict prognosis in patients with HCC, its relation with immune infiltration, and chemotherapy sensitivity.
Materials and Methods
Hepatocellular Carcinoma Datasets
The RNA sequencing data of the HCC patients, normal tissues, and their corresponding clinical information were obtained from the TCGA database on 30 June 2021. The TCGA datasets were downloaded from the UCSC Xane website [UCSC Xena (xenabrowser.net)]. Furthermore, to increase the reliability of our study, we extracted the RNA sequencing data of HCC patients and normal tissues from a GEO dataset “GSE76427” and used it for validation of the results from the TCGA datasets.
Identification of Differentially Expressed m7G-Related Genes in HCC
Lists of the m7G-related genes were extracted from the prior studies (Song et al., 2020; Xia et al., 2021). The expression of these m7G-related genes in the TCGA HCC cohorts and normal tissues were identified using the “limma” package in R.
Functional Enrichment Analysis of Differentially Expressed m7G-Related Genes
The “clusterprofile”, “enrichplot”, “ggplot2”, “GOplot”, and “library (ReactomePA)” packages in R were used to perform Gene Ontology (GO) and Reactome pathway analysis to identify the significant pathways (adj p < 0.05) associated with differentially expressed m7G-related genes.
Development and Validation of PRG Prognostic Model
The “survival” and “glmnet” packages in R were used for the development of a prognostic model. The relevant clinical information of the patients was extracted from the TCGA database and each gene was screened by using the univariate Cox-regression analysis. The least absolute shrinkage and selection operator (LASSO) regression was used to incorporate the m7G-related genes with p-value <0.05 in the univariate analysis. Identified genes from the LASSO analysis were further filtered by the multivariate Cox-regression analysis to establish a risk score. A “forestplot” R package was used to create a forest plot for univariate and multivariate outcomes. The Kaplan-Meier survival curve was constructed for each prognostic gene. Lastly, the identified genes from the multivariate analysis were used to construct the risk score using the formula:
Risk score = (expression level of Gene1 ∗ β1 + expression level of Gene2 ∗ β2 +...+expression level of Genen ∗ βn).
The TCGA HCC cohorts were then divided into high- and low-risk groups based on the median risk score and the corresponding Kaplan-Meier analyses were performed comparing these two subgroups. Finally, the AUC was calculated for 3- and 5-year survival in the ROC curves.
Construction of a Predictive Nomogram
The “survival” “survminer”, and “rms” packages in R were used for the development of the nomogram. Univariate and multivariate cox regression models were used to identify the independent factors affecting the prognosis of HCC patients. Based on the clinicopathologic features of m7G-related genes from multivariate analysis, a cox-proportional hazard (PH) model was used to develop a nomogram to predict the survival probability. The goodness-of-fit test purposed by Schoenfeld was used to test the validity of PH models.
Immune Infiltration and Chemotherapy Sensitivity Analysis
CYBERSORT tool was used to estimate the expression of 22 different immune cells in HCC patients. We calculated the immune, stromal, and ESTIMATE scores in the high and low-risk groups by using “estimate” package in R. Furthermore, the relation of risk-score with immune checkpoints and immune-suppressive cytokines was also assessed. Furthermore, we checked if the risk model could serve as a potential predictor of chemotherapeutic sensitivity to sorafenib. We converted the gene expression TCGA datasets into half inhibitory concentrations (IC50) data matrix of the sorafenib with the “pRRophetic” package and analyzed the IC50 difference between the high-risk and low-risk groups.
Statistical Analyses
Statistical analyses were performed using the R statistical software version 4.1.2. Gene expression levels between the HCC tumor and normal tissue were performed with a single-factor analysis of variance. Categorical variables were compared by using the Pearson chi-square test. Kaplan-Meier analysis was used to compare the survival outcomes. The univariate and multivariate Cox regression analyses were used to evaluate the independent prognostic value of the identified risk model. All the survival outcomes were presented as p-values and hazard ratios (HR) with 95% confidence intervals (CI). Mann-Whitney test was used to compare the immune infiltration and pathway activation between the two groups.
Results
Identification of the Differentially Expressed m7G-Related Genes in HCC
A total of 503 M7G-related genes were extracted from the previous studies. Out of them, 492 genes were expressed in HCC and 18 m7G-related genes were differentially expressed (upregulated: 5 and downregulated: 13) in the HCC than in the normal tissues (Figure 1A). Significantly upregulated genes were CDK1, RBM3, LMNA, COL4A2, and TYMS, while the expression of KBTBD11, BCKDHB, ESR1, IYD, PDGFRA, LDLR, CA2, CCL2, ANO1, SLCO2B1, AGXT, EFHD1, and MTTP were downregulated in the tumor. The heatmap of the DEG is shown in Figure 1B.
FIGURE 1. Identification of the differentially expressed m7G-related genes and functional enrichment analysis. (A) Volcano plot of the differentially expressed m7g-related genes. (B) Heatmap of the differentially expressed m7G-related genes. (C) GO functional enrichment analysis. (D) Reactome pathway analysis.
Functional Enrichment Analyses of the m7G-Related Genes
GO enrichment analysis of the differentially expressed genes showed that the m7G-related genes were associated with several biological processes, cellular components, and molecular functions. The m7G-related genes were significantly associated with the rhythmic process, microvillus, apical part of the cell, mitochondrial matrix, basolateral plasma membrane, actin-based cell projection, basal plasma membrane, and basal part of the cell (Figure 1C). The Reactome pathway analysis showed that significantly expressed m7G-related genes were associated with depolymerization of the nuclear lamina, initiation of the nuclear envelope reformation, G1/S-specific transcription, glyoxylate metabolism and glycine degradation, and ovarian tumor domain proteases (Figure 1D).
Construction of m7G-Related Genes Prognostic Model
The univariate Cox-regression analysis revealed that 6 DEGs (CDK1 (HR: 1.331; 95% CI: 1.119–1.585; p < 0.05), TYMS (HR: 1.176; 95% CI: 1.012–1.366; p < 0.05), AGXT (HR: 0.920; 95% CI: 0.847–1.000; p < 0.05), ANO1 (HR: 0.870; 95% CI: 0.771–0.983; p < 0.05), IYD (HR: 0.740; 95% CI: 0.573–0.957; p < 0.05), and PDGFRA (HR: 1.341; 95% CI: 1.068–1.682; p < 0.05)) were significantly associated with the prognosis of HCC (Figure 2A). Only 5 DEGs (CDK1, AGXT, ANO1, IYD, and PDGFRA) were identified as prognostic genes after LASSO regression analysis (Figures 2B,C). Furthermore, the multivariate Cox-regression analysis identified only 3 DEGs (CDK1, ANO1, and PDGFRA) for the construction of a prognostic model in patients with HCC (Table 1). The patients were then classified into high-risk and low-risk groups based on the median risk score and the Kaplan-Meier survival analysis was performed. The risk score was calculated as: The Risk score: Risk score = (0.251) ∗ CDK1 + (−0.122) ∗ ANO1 + (0.318) ∗ PDGFRA. Kaplan-Meier survival curve showed that the prognosis of patients in the high-risk group was significantly poorer than those in the low-risk group (p < 0.05) (Figure 2D) with AUCs of 0.658 and 0.617 in 3- and 5-year ROC curves respectively (Figure 2E). As the risk score increased, the risk of death increased but the duration of survival decreased (Figures 3A,B). A Heatmap of the expression of three m7G-related genes in HCC is shown in Figure 3C. Furthermore, we identified that the risk score was significantly associated with the patient’s age (>60 years), survival status, tumor grade, and T-stage (Figures 3D–G).
FIGURE 2. Construction of the risk score. (A) Results of the univariate analysis; (B,C) Lasso regression analysis of the genes from univariate analysis. (D) Survival probability based on the high and low-risk group. (E) ROC curve for the 3- and 5-year overall survival.
TABLE 1. The multivariate Cox regression analysis for identification of m7g-related genes associated with prognosis in HCC.
FIGURE 3. Construction of the risk score. (A,B) Distribution of the survival, risk score, and survival of patients with HCC. (C) Heatmap of the expression of three m7G-related genes in HCC. (D) Relation of the risk score and patient’s age. (E) Relation of the risk score with survival status. (F) Relation of the risk score with tumor grade. (G) Relation of the risk score with T-stage.
Predictive Nomogram
A predictive nomogram was constructed based on the clinicopathologic features and the 3 prognostic m7G-related genes. Results of the univariate showed that the T-stage, M-stage, and risk score were significantly associated with prognosis (Figure 4A). The goodness-of-fit test proposed by Schoenfeld was used to check whether the PH assumption was valid. The Schoenfeld individual test p values for risk score age, gender, N-stage, T-stage, M-stage, tumor grade, and risk-score were 0.3689, 0.0605, 0.0896, 0.6864, 0.4132, and 0.8801 respectively. The p-value for the global Schoenfeld test was 0.7701 (Supplementary Table S1) (Supplementary Figure S1). The results indicated that all of the factors were satisfied for the PH assumption. Multivariate analysis revealed that only the risk score and T-stage were the independent factors predicting poor prognosis in patients with HCC (Figure 4B). ROC curve of the clinicopathologic features and risk score revealed that the T-stage (AUC of 0.737) and risk score (AUC of 0.709) can better predict the prognosis than other clinicopathologic factors (Figure 4C). The predictive nomogram showed a better prediction of 3- and 5-years overall survival than the other clinicopathological factors (Concordance index: 0.718) (Figure 4D). The calibration curve for 1-, 3-, and 5-year survival showed that the nomogram can accurately predict survival (Figure 4E).
FIGURE 4. Predictive nomogram. (A) Univariate analysis of the clinicopathologic features and the risk score. (B) Multivariate analysis of the clinicopathologic features and the risk score. (C) ROC curve of the clinicopathologic features and risk score. (D) Nomogram to predict the survival of the HCC patients. (E) Calibration curve for 1-, 3-, and 5-year survival.
Association of m7G-Related Genes With Immune Infiltration in the Tumor Microenvironment and Chemotherapy Sensitivity
We identified that the T cells CD4 memory resting cells were significantly higher in the low-risk group while the T-cells CD4 memory activated cells were higher in the high-risk group (Figure 5A). Furthermore, immune scores and ESTIMATE scores were significantly higher in the high-risk group but there was no significant difference in the stromal score (Figure 5B). Among the several genes involved in immune checkpoints, CTLA4, HAVCR2, LAG3, and TIGT were significantly higher in the high-risk group (Figure 5C). The immune-suppressive cytokines like IL10, IL4, TGFB1, TGFB2, and TGFB3 were also significantly higher in the high-risk group (Figure 5D). We also found that the high-risk group patients were more sensitive to the treatment with sorafenib than the low-risk group (Figure 5E).
FIGURE 5. Immune infiltration analysis and chemotherapy sensitivity. (A) Relation of the immune cells with a risk score. (B) Immune, stromal, and ESTIMATE scores in the high and low-risk groups. (C) Relation of the immune checkpoints with risk-score. (D): Relation of the immune-suppressive cytokines with a risk score. (E) Sorafenib sensitivity in high-risk and low-risk groups.
Validation of the Risk Model Accuracy in the GEO Dataset
The validation of the risk model was performed in the GEO dataset “GSE76427”. The 3 prognostic genes were significantly expressed in the validation dataset (Figure 6A). Similar to the training dataset, as the risk score increased, the risk of death increased but the duration of survival decreased in the validation dataset. (Figures 6B,C). Kaplan-Meier survival curve in the validation dataset showed that the prognosis of patients in the high-risk group was significantly poorer than those in the low-risk group (p < 0.001) (Figure 6D) with AUCs of 0.681 and 0.972 in 3- and 5-year ROC curves respectively (Figure 6E).
FIGURE 6. Validation of the risk score in the GEO dataset. (A) Heatmap of the 3 differentially expressed m7G-related genes in the validation dataset. (B,C) Distribution of the survival, risk score, and survival of patients with HCC. (D) Survival probability based on the high and low-risk group. (E) ROC curve for the 3- and 5-year overall survival.
Discussion
With the advances in high-throughput sequencing technology, more and more studies have published different kinds of RNA modification including the internal mRNA m7G methylation. In clinical research, several previous studies identified the role of N6-methyladenosine (m6A) remodeling in the molecular pathogenesis of several types of carcinogenesis and the role of m6A-related genes in the prognosis of several types of cancers (Chai et al., 2021; Jin et al., 2021; Han et al., 2022). M6A remodeling has been well studied in HCC as some studies have identified m6A related gene signature model for the prediction of HCC (Jin et al., 2021; Wang et al., 2021). Recently, some studies have shown that some of the m7G-related genes like methyltransferase-like protein-1 (METTL1) and WD repeat domain 4 (WDR4) are involved in the proliferation, migration, and invasion of HCC through m7G tRNA modification (Chen et al., 2021; Xia et al., 2021). However, to our knowledge, the construction of the m7G-related genes prognostic model and its association with immunotherapy in patients with HCC has not been studied in detail.
Our study constructed an m7G associated prognostic model and risk score for predicting the prognosis of HCC which will help in understanding its molecular pathogenesis and will further benefit in understanding the diagnosis, treatment, and prognosis. The results of our study identified that CDK1, ANO1, and PGDFRA-based prognostic models can accurately predict the prognosis in patients with HCC. Mutation in m7G methyltransferase is associated with several diseases like Down syndrome, primordial dwarfism, ovarian carcinogenesis, etc., (Pereira et al., 2009; Shaheen et al., 2015; Trimouille et al., 2018; Xia et al., 2021). Furthermore, it is reported that mutation in m7G sites is also associated with neural line differentiation in embryonic stem cells and it also affects cancer cell viability (Barbieri et al., 2017; Lin et al., 2018).
Cyclin-dependent kinase 1 (CDK1) is highly expressed in HCC tissues and is associated with poor prognosis (Yang et al., 2019; Qiang et al., 2021). It may act as a promising diagnostic, therapeutic, and prognostic marker as well as a promising therapeutic target (Zhou et al., 2019; Yin et al., 2021). Anoctamin 1 (ANO1) regulates cell proliferation during the G (1)/S transition of the cell cycle and is associated with tumorigenesis of several gastrointestinal tumors (Stanich et al., 2011; Zhang et al., 2019; Dong et al., 2021). Overexpression of ANO1 is one of the markers of poor prognosis in HCC (Zhang et al., 2019). Platelet-derived growth factor receptor alpha (PGDFRA) mutation is associated with uncontrolled proliferation and cell growth of gastrointestinal and its inhibition blocks the activity of PDGFRA kinase and improves survival (Kovács et al., 2005; Comandone and Boglione, 2015). PDGFRA is associated with the pathogenesis of hepatic fibrosis during chronic liver injury and is also a prognostic biomarker for HCC (Kikuchi et al., 2017; Yu et al., 2017).
Results of the functional enrichment analysis of the m7G-related genes with GO revealed that some of them were associated with the pathogenesis of HCC. Mitochondrial dysregulation is associated with the aggressiveness of HCC (Min et al., 2021). Actin−mediated cellular activity is associated with cell migration and invasion, which may play important role in the development of cancer (Guan et al., 2020). Proteins like PTEN and aquaporin-9 located in the basal plasma and basolateral membranes have shown some anti-tumorigenic effects against HCC in several preclinical studies (Kato et al., 2020; Liao et al., 2020). These molecules may act as a therapeutic target for the treatment of HCC. Reactome pathway analysis showed that m7G- related genes are associated with nuclear envelope dysregulation (depolymerization and reformation). Nuclear envelope dysregulation and laminopathies are associated with mitogen-activated protein kinase (MAPK) signaling and AKT- mammalian target of rapamycin (mTOR) signaling pathways (Choi and Worman, 2014). These signaling pathways are associated with the pathogenesis of HCC (Pang et al., 2020; Luo et al., 2021). G1 to S transition is a crucial step in HCC proliferation and is negatively associated with prognosis (Zeng et al., 2019). Therefore, identification of the m7G modification sites and their host genes could further help in understanding the pathogenesis of several malignancies and identification of target genes and pathways for the development of newer drugs targeting the m7G sites.
Immune infiltration is also associated with carcinogenesis and the prognosis of HCC (Han et al., 2021; Xu et al., 2021). Moreover, immunotherapy with checkpoint inhibitors has strong anti-tumor activity in HCC (Li et al., 2020; Sangro et al., 2021). Activation of the negative co-stimulatory molecules like cytotoxic T-lymphocyte- associated antigen 4 (CTLA4), LAG3, and TIGIT are associated with the progression of HCC (Moreno-Cubero and Larrubia, 2016; Yu et al., 2021). Clinical studies on anti-CTLA-4 antibodies like tremelimumab and ipilimumab are studying the efficacy and safety with some promising results (Sangro et al., 2013). Therefore, immune checkpoint inhibitors targeting the above-mentioned molecules may be one of the novel therapeutic options for the treatment of HCC. LAG3 is also a predictor of tumor response following therapy for HCC (Guo et al., 2021). We found that the immune-suppressive cytokines like IL10, IL4, TGFB1, TGFB2, and TGFB3 were significantly expressed in patients with high-risk scores. Immune-suppressive cytokines are associated with anti-tumor immune response and the cytokines like IL and TGFB are associated with the pathogenesis of HCC (Zhang et al., 2020; Zhang et al., 2022). Sorafenib is the standard treatment for patients with advanced disease with improvement in survival and radiologic progression of the disease (Llovet et al., 2008). We found that the patients with high-risk scores are more sensitive than those with low-risk scores. This suggests that the developed risk model can not only stratify the prognosis but can also help in identifying the appropriate subgroup of patients for chemotherapy. To our knowledge, there is no single established biomarker for predicting the response to immune-related therapy for the management of advanced HCC. The prognostic marker we developed, with further clinical validation in the future may help the clinicians to guide the appropriate treatment strategy for specific patients.
Several previous studies have shown the reliability of using risk scores for the prediction of the prognosis of HCC. Our risk-score model showed that the patients with high-risk scores had significantly worse survival than those with low-risk scores. The risk score was positively correlated with other well-known prognostic factors like T-stage, N-stage, M-stage, and tumor grade. Furthermore, a quantitative nomogram was developed with the inclusion of risk score and several well-known existing prognostic factors of HCC, and its prognostic ability is better than other existing clinicopathologic factors. This further improved the performance and reliability of the risk score. Therefore, we believe that the proposed model can be used for the stratification of prognosis and identify the appropriate therapies. Lastly, the risk score was used to assess the prognosis in a GEO dataset, which showed that the risk score can accurately predict the prognosis. Therefore, we believe with further validation from future clinical trials, this could be one of the novel prognostic models with higher accuracy.
There are certain limitations to our study. First, in vivo and in vitro experiments were not performed to verify the results obtained from our study. Second, the results of this study were not validated in the resected specimens of HCC from our institution. Despite these limitations, we believe the results of our study will provide knowledge about the prognostic significance of m7G-related genes in HCC and will certainly help in guiding future clinical studies.
Conclusion
Our study identified that the 3 m7G-related signature model can be used as a prognostic biomarker in HCC. The functional enrichment analysis of the m7G-related genes showed some of the molecules associated with the pathogenesis and pathways for the progression of HCC. In addition, m7G-related genes were also positively correlated with immune cell infiltration. Future clinical study on this biomarker model is required to verify its clinical implications.
Data Availability Statement
Publicly available databases were analysed in the study. These datasets can be found here: xenabrowser.net and ncbi.nlm.nih.gov/gds.
Author Contributions
PR and Z-QH contributed equally to the data acquisition, statistical analysis, drafted the manuscript, and share the first authorship. TL and AP were involved in editing and reviewing the manuscript. F-YL contributed to the study design and revision of the manuscript.
Funding
This study was funded by the Fundamental Research Funds for the Central Universities (2019SCUH0021), and the Science and Technology Support Project of Sichuan Province (2021YJ0132).
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.2022.918983/full#supplementary-material
References
Barbieri, I., Tzelepis, K., Pandolfini, L., Shi, J., Millán-Zambrano, G., Robson, S. C., et al. (2017). Promoter-bound METTL3 Maintains Myeloid Leukaemia by m6A-Dependent Translation Control. Nature 552, 126–131. doi:10.1038/nature24678
Bray, F., Ferlay, J., Soerjomataram, I., Siegel, R. L., Torre, L. A., and Jemal, A. (2018). Global Cancer Statistics 2018: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. CA A Cancer J. Clin. 68, 394–424. doi:10.3322/caac.21492
Bruix, J., Reig, M., and Sherman, M. (2016). Evidence-Based Diagnosis, Staging, and Treatment of Patients with Hepatocellular Carcinoma. Gastroenterology 150, 835–853. doi:10.1053/j.gastro.2015.12.041
Chai, X.-k., Qi, W., Zou, C.-Y., He, C.-X., Su, M., and Zhao, D.-Q. (2021). Potential Prognostic Value of a Seven m6A-Related LncRNAs Signature and the Correlative Immune Infiltration in Colon Adenocarcinoma. Front. Genet. 12, 774010. doi:10.3389/fgene.2021.774010
Chen, Z., Zhu, W., Zhu, S., Sun, K., Liao, J., Liu, H., et al. (2021). METTL1 Promotes Hepatocarcinogenesis via M7G tRNA Modification-dependent Translation Control. Clin. Transl. Med. 11, e661. doi:10.1002/ctm2.661
Choi, J. C., and Worman, H. J. (2014). Nuclear Envelope Regulation of Signaling Cascades. Adv. Exp. Med. Biol. 773, 187–206. doi:10.1007/978-1-4899-8032-8_9
Comandone, A., and Boglione, A. (2015). The Importance of Mutational Status in Prognosis and Therapy of GIST. Recenti Prog. Med. 106, 17–22. doi:10.1701/1740.18950
Dong, Y., Tian, J., Yan, B., Lv, K., Li, J., and Fu, D. (2021). Liver-Metastasis-Related Genes Are Potential Biomarkers for Predicting the Clinical Outcomes of Patients with Pancreatic Adenocarcinoma. Pathol. Oncol. Res. 27, 1609822. doi:10.3389/pore.2021.1609822
Global Burden of Disease Cancer, C., Fitzmaurice, C., Allen, C., Barber, R. M., Barregard, L., Bhutta, Z. A., et al. (2017). Global, Regional, and National Cancer Incidence, Mortality, Years of Life Lost, Years Lived with Disability, and Disability-Adjusted Life-Years for 32 Cancer Groups, 1990 to 2015: A Systematic Analysis for the Global Burden of Disease Study. JAMA Oncol. 3, 524–548. doi:10.1001/jamaoncol.2016.5688
Guan, X., Guan, X., Dong, C., and Jiao, Z. (2020). Rho GTPases and Related Signaling Complexes in Cell Migration and Invasion. Exp. Cell Res. 388, 111824. doi:10.1016/j.yexcr.2020.111824
Guo, M., Qi, F., Rao, Q., Sun, J., Du, X., Qi, Z., et al. (2021). Serum LAG-3 Predicts Outcome and Treatment Response in Hepatocellular Carcinoma Patients with Transarterial Chemoembolization. Front. Immunol. 12, 754961. doi:10.3389/fimmu.2021.754961
Guy, M. P., and Phizicky, E. M. (2014). Two-subunit Enzymes Involved in Eukaryotic Post-transcriptional tRNA Modification. RNA Biol. 11, 1608–1618. doi:10.1080/15476286.2015.1008360
Han, R., Feng, P., Pang, J., Zou, D., Li, X., Geng, C., et al. (2021). A Novel HCC Prognosis Predictor EEF1E1 Is Related to Immune Infiltration and May Be Involved in EEF1E1/ATM/p53 Signaling. Front. Oncol. 11, 700972. doi:10.3389/fonc.2021.700972
Han, S, Qi, J, Fang, K, Wang, H, Tang, Y, Depei, W, et al. (2022). Characterization of m6A Regulator-Mediated Methylation Modification Patterns and Tumor Microenvironment Infiltration in Acute Myeloid Leukemia. Cancer Med 11(5):1413–1426. doi:10.1002/cam4.4531
Jin, C., Li, R., Deng, T., Li, J., Yang, Y., Li, H., et al. (2021). Identification and Validation of a Prognostic Prediction Model of m6A Regulator-Related LncRNAs in Hepatocellular Carcinoma. Front. Mol. Biosci. 8, 784553. doi:10.3389/fmolb.2021.784553
Kato, T., Yamada, T., Nakamura, H., Igarashi, A., Anders, R. A., Sesaki, H., et al. (2020). The Loss of Nuclear PTEN Increases Tumorigenesis in a Preclinical Mouse Model for Hepatocellular Carcinoma. iScience 23, 101548. doi:10.1016/j.isci.2020.101548
Kikuchi, A., Pradhan-Sundd, T., Singh, S., Nagarajan, S., Loizos, N., and Monga, S. P. (2017). Platelet-Derived Growth Factor Receptor α Contributes to Human Hepatic Stellate Cell Proliferation and Migration. Am. J. Pathology 187, 2273–2287. doi:10.1016/j.ajpath.2017.06.009
Kovács, M., Nagy, P., Pák, G., and Fehér, J. (2005). Gastrointestinal Stromal Tumors (GISTs): Clinical and Pathological Features. Orv. Hetil. 146, 1375–1381.
Kudo, M., Izumi, N., Sakamoto, M., Matsuyama, Y., Ichida, T., Nakashima, O., et al. (2016). Survival Analysis over 28 Years of 173,378 Patients with Hepatocellular Carcinoma in Japan. Liver Cancer 5, 190–197. doi:10.1159/000367775
Li, J., Xing, J., Yang, Y., Liu, J., Wang, W., Xia, Y., et al. (2020). Adjuvant 131I-Metuximab for Hepatocellular Carcinoma after Liver Resection: a Randomised, Controlled, Multicentre, Open-Label, Phase 2 Trial. Lancet Gastroenterology Hepatology 5, 548–560. doi:10.1016/s2468-1253(19)30422-4
Liao, S., Chen, H., Liu, M., Gan, L., Li, C., Zhang, W., et al. (2020). Aquaporin 9 Inhibits Growth and Metastasis of Hepatocellular Carcinoma Cells via Wnt/β-Catenin Pathway. Aging 12, 1527–1544. doi:10.18632/aging.102698
Lin, S., Liu, Q., Lelyveld, V. S., Choe, J., Szostak, J. W., and Gregory, R. I. (2018). Mettl1/Wdr4-Mediated m7G tRNA Methylome Is Required for Normal mRNA Translation and Embryonic Stem Cell Self-Renewal and Differentiation. Mol. Cell 71, 244–255. doi:10.1016/j.molcel.2018.06.001
Llovet, J. M., Ricci, S., Mazzaferro, V., Hilgard, P., Gane, E., Blanc, J.-F., et al. (2008). Sorafenib in Advanced Hepatocellular Carcinoma. N. Engl. J. Med. 359, 378–390. doi:10.1056/nejmoa0708857
Llovet, J. M., De Baere, T., Kulik, L., Haber, P. K., Greten, T. F., Meyer, T., et al. (2021). Locoregional Therapies in the Era of Molecular and Immune Treatments for Hepatocellular Carcinoma. Nat. Rev. Gastroenterol. Hepatol. 18, 293–313. doi:10.1038/s41575-020-00395-0
Luo, Y.-D., Fang, L., Yu, H.-Q., Zhang, J., Lin, X.-T., Liu, X.-Y., et al. (2021). p53 Haploinsufficiency and Increased mTOR Signalling Define a Subset of Aggressive Hepatocellular Carcinoma. J. Hepatology 74, 96–108. doi:10.1016/j.jhep.2020.07.036
Ma, J., Han, H., Huang, Y., Yang, C., Zheng, S., Cai, T., et al. (2021). METTL1/WDR4-mediated m7G tRNA Modifications and m7G Codon Usage Promote mRNA Translation and Lung Cancer Progression. Mol. Ther. 29, 3422–3435. doi:10.1016/j.ymthe.2021.08.005
Malbec, L., Zhang, T., Chen, Y.-S., Zhang, Y., Sun, B.-F., Shi, B.-Y., et al. (2019). Dynamic Methylome of Internal mRNA N7-Methylguanosine and its Regulatory Role in Translation. Cell Res. 29, 927–941. doi:10.1038/s41422-019-0230-z
Min, S., Lee, Y.-K., Hong, J., Park, T. J., Woo, H. G., Kwon, S. M., et al. (2021). MRPS31 Loss Is a Key Driver of Mitochondrial Deregulation and Hepatocellular Carcinoma Aggressiveness. Cell Death Dis. 12, 1076. doi:10.1038/s41419-021-04370-8
Moreno-Cubero, E., and Larrubia, J.-R. (2016). Specific CD8+T Cell Response Immunotherapy for Hepatocellular Carcinoma and Viral Hepatitis. World. J. Gastroe 22, 6469–6483. doi:10.3748/wjg.v22.i28.6469
Pang, Y., Liu, Z., Han, H., Wang, B., Li, W., Mao, C., et al. (2020). Peptide SMIM30 Promotes HCC Development by Inducing SRC/YES1 Membrane Anchoring and MAPK Pathway Activation. J. Hepatology 73, 1155–1169. doi:10.1016/j.jhep.2020.05.028
Pereira, P. L., Magnol, L., Sahún, I., Brault, V., Duchon, A., Prandini, P., et al. (2009). A New Mouse Model for the Trisomy of the Abcg1-U2af1 Region Reveals the Complexity of the Combinatorial Genetic Code of Down Syndrome. Hum. Mol. Genet. 18, 4756–4769. doi:10.1093/hmg/ddp438
Qiang, R., Zhao, Z., Tang, L., Wang, Q., Wang, Y., and Huang, Q. (2021). Identification of 5 Hub Genes Related to the Early Diagnosis, Tumour Stage, and Poor Outcomes of Hepatitis B Virus-Related Hepatocellular Carcinoma by Bioinformatics Analysis. Comput. Math. Methods Med. 2021, 9991255. doi:10.1155/2021/9991255
Rizzo, A., and Brandi, G. (2021). Biochemical Predictors of Response to Immune Checkpoint Inhibitors in Unresectable Hepatocellular Carcinoma. Cancer Treat. Res. Commun. 27, 100328. doi:10.1016/j.ctarc.2021.100328
Rizzo, A., and Ricci, A. D. (2021). Biomarkers for Breast Cancer Immunotherapy: PD-L1, TILs, and Beyond. Expert Opin. Investigational Drugs, 1–7. doi:10.1080/13543784.2022.2008354
Sangro, B., Gomez-Martin, C., de la Mata, M., Iñarrairaegui, M., Garralda, E., Barrera, P., et al. (2013). A Clinical Trial of CTLA-4 Blockade with Tremelimumab in Patients with Hepatocellular Carcinoma and Chronic Hepatitis C. J. Hepatology 59, 81–88. doi:10.1016/j.jhep.2013.02.022
Sangro, B., Sarobe, P., Hervás-Stubbs, S., and Melero, I. (2021). Advances in Immunotherapy for Hepatocellular Carcinoma. Nat. Rev. Gastroenterol. Hepatol. 18, 525–543. doi:10.1038/s41575-021-00438-0
Shaheen, R., Abdel-Salam, G. M. H., Guy, M. P., Alomar, R., Abdel-Hamid, M. S., Afifi, H. H., et al. (2015). Mutation in WDR4 Impairs tRNA m7G46 Methylation and Causes a Distinct Form of Microcephalic Primordial Dwarfism. Genome Biol. 16, 210. doi:10.1186/s13059-015-0779-x
Sloan, K. E., Warda, A. S., Sharma, S., Entian, K.-D., Lafontaine, D. L. J., and Bohnsack, M. T. (2017). Tuning the Ribosome: The Influence of rRNA Modification on Eukaryotic Ribosome Biogenesis and Function. RNA Biol. 14, 1138–1152. doi:10.1080/15476286.2016.1259781
Song, B., Tang, Y., Chen, K., Wei, Z., Rong, R., Lu, Z., et al. (2020). m7GHub: Deciphering the Location, Regulation and Pathogenesis of Internal mRNA N7-Methylguanosine (m7G) Sites in Human. Bioinformatics 36, 3528–3536. doi:10.1093/bioinformatics/btaa178
Stanich, J. E., Gibbons, S. J., Eisenman, S. T., Bardsley, M. R., Rock, J. R., Harfe, B. D., et al. (2011). Ano1 as a Regulator of Proliferation. Am. J. Physiology-Gastrointestinal Liver Physiology 301, G1044–G1051. doi:10.1152/ajpgi.00196.2011
Trimouille, A., Lasseaux, E., Barat, P., Deiller, C., Drunat, S., Rooryck, C., et al. (2018). Further Delineation of the Phenotype Caused by Biallelic Variants in the WDR4 Gene. Clin. Genet. 93, 374–377. doi:10.1111/cge.13074
Wang, Y., Li, N., Tian, D., Zhou, C.-W., Wang, Y.-H., Yang, C., et al. (2021). Analysis of m6A-Related lncRNAs for Prognosis Value and Response to Immune Checkpoint Inhibitors Therapy in Hepatocellular Carcinoma. Cmar 13, 6451–6471. doi:10.2147/cmar.s322179
Xia, P., Zhang, H., Xu, K., Jiang, X., Gao, M., Wang, G., et al. (2021). MYC-targeted WDR4 Promotes Proliferation, Metastasis, and Sorafenib Resistance by Inducing CCNB1 Translation in Hepatocellular Carcinoma. Cell Death Dis. 12, 691. doi:10.1038/s41419-021-03973-5
Xu, Y., He, X., Deng, J., Xiong, L., Li, Y., Zhang, X., et al. (2021). Comprehensive Analysis of the Immune Infiltrates and PD-L1 of m6A RNA Methylation Regulators in Hepatocellular Carcinoma. Front. Cell Dev. Biol. 9, 681745. doi:10.3389/fcell.2021.681745
Yang, W. X., Pan, Y. Y., and You, C. G. (2019). CDK1, CCNB1, CDC20, BUB1, MAD2L1, MCM3, BUB1B, MCM2, and RFC4 May Be Potential Therapeutic Targets for Hepatocellular Carcinoma Using Integrated Bioinformatic Analysis. Biomed. Res. Int. 2019, 1245072. doi:10.1155/2019/1245072
Yin, S., Yang, S., Luo, Y., Lu, J., Hu, G., Wang, K., et al. (2021). Cyclin-dependent Kinase 1 as a Potential Target for Lycorine against Hepatocellular Carcinoma. Biochem. Pharmacol. 193, 114806. doi:10.1016/j.bcp.2021.114806
Yu, J.-H., Kim, J. M., Kim, J. K., Choi, S. J., Lee, K. S., Lee, J.-W., et al. (2017). Platelet-derived Growth Factor Receptor α in Hepatocellular Carcinoma Is a Prognostic Marker Independent of Underlying Liver Cirrhosis. Oncotarget 8, 39534–39546. doi:10.18632/oncotarget.17134
Yu, L., Liu, X., Wang, X., Yan, F., Wang, P., Jiang, Y., et al. (2021). TIGIT+ TIM-3+ NK Cells Are Correlated with NK Cell Exhaustion and Disease Progression in Patients with Hepatitis B Virus Related Hepatocellular Carcinoma. Oncoimmunology 10, 1942673. doi:10.1080/2162402x.2021.1942673
Zeng, G., Wang, J., Huang, Y., Lian, Y., Chen, D., Wei, H., et al. (2019). Overexpressing CCT6A Contributes to Cancer Cell Growth by Affecting the G1-To-S Phase Transition and Predicts A Negative Prognosis in Hepatocellular Carcinoma. Ott 12, 10427–10439. doi:10.2147/ott.s229231
Zhang, R., Ye, J., Huang, H., and Du, X. (2019). Mining Featured Biomarkers Associated with Vascular Invasion in HCC by Bioinformatics Analysis with TCGA RNA Sequencing Data. Biomed. Pharmacother. 118, 109274. doi:10.1016/j.biopha.2019.109274
Zhang, B., Liu, Z., Cao, K., Shan, W., Liu, J., Wen, Q., et al. (2020). Circ‐SPECC1 Modulates TGFβ2 and Autophagy under Oxidative Stress by Sponging miR‐33a to Promote Hepatocellular Carcinoma Tumorigenesis. Cancer Med. 9, 5999–6008. doi:10.1002/cam4.3219
Zhang, M., Zhao, H., and Gao, H. (2022). Interleukin-24 Limits Tumor-Infiltrating T Helper 17 Cell Response in Patients with Hepatitis B Virus-Related Hepatocellular Carcinoma. Viral Immunol. 35 (3), 212–222. doi:10.1089/vim.2021.0174
Keywords: hepatocellular carcinoma, m7G-related genes, prognosis, nomogram, risk-score
Citation: Regmi P, He Z-Q, Lia T, Paudyal A and Li F-Y (2022) N7-Methylguanosine Genes Related Prognostic Biomarker in Hepatocellular Carcinoma. Front. Genet. 13:918983. doi: 10.3389/fgene.2022.918983
Received: 13 April 2022; Accepted: 24 May 2022;
Published: 06 June 2022.
Edited by:
Manoj Kumar Kashyap, Amity University Gurgaon, IndiaReviewed by:
Greesham Tripathi, Amity University Gurgaon, IndiaAshraf A. Tabll, National Research Centre, Egypt
Copyright © 2022 Regmi, He, Lia, Paudyal and Li. 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: Fu-Yu Li, lfy74@hotmail.com
†These authors have contributed equally to this work and share first authorship