Systematic Analysis Identifies a Specific RNA-Binding Protein-Related Gene Model for Prognostication and Risk-Adjustment in HBV-Related Hepatocellular Carcinoma

Objective: Increasing evidence shows that dysregulated RNA binding proteins (RBPs) modulate the progression of several malignancies. Nevertheless, their clinical implications of RBPs in HBV-related hepatocellular carcinoma (HCC) remain largely undefined. Here, this study systematically analyzed the associations of RBPs with HBV-related HCC prognosis.

Methods: Based on differentially expressed RBPs between HBV-related HCC and control specimens, prognosis-related RBPs were screened by univariate analyses. A LASSO model was then created. Kaplan-Meier curves, ROCs, multivariate analyses, subgroup analyses and external verification were separately applied to assess the efficacy of this model in predicting prognosis and recurrence of patients. A nomogram was created by incorporating the model and clinical indicators, which was verified by ROCs, calibration curves and decision curve analyses. By CIBERSORT algorithm, the association between the risk score and immune cell infiltrations was evaluated.

Results: Totally, 54 RBPs were distinctly correlated to prognosis of HBV-related HCC. An 11-RBP model was created, containing POLR2L, MRPS12, DYNLL1, ZFP36, PPIH, RARS, SRP14, DDX41, EIF2B4, and NOL12. This risk score sensitively and accurately predicted one-, three- and five-year overall survival, disease-free survival, and progression-free interval. Compared to other clinical parameters, this risk score had the best predictive efficacy. Also, the clinical generalizability of the model was externally verified in the GSE14520 dataset. The nomogram may predict patients’ survival probabilities. Also, the risk score was related to the components in the immune microenvironment.

Conclusion: Collectively, RBPs may act as critical elements in the malignant progression of HBV-related HCC and possess potential implications on prognostication and therapy decision.

Introduction

Hepatocellular carcinoma (HCC) is the most prevalent type of liver cancer and represents a common malignant neoplasm globally (Liu et al., 2019). Because of the high risk of recurrence and metastasis, the 5-year survival probabilities of advanced HCC are still undesirable. The epidemiology of HCC is affected by underlying liver diseases especially hepatitis B virus (HBV) (Chabrolles et al., 2020). It has been estimated that HBV infection is responsible for 50% of HCC worldwide (Li et al., 2019). HBV that integrates into cancer-relevant genes may drive hepatocarcinogenesis (Nakagawa et al., 2019). Nevertheless, the mechanism by which HBV infection contributes to HCC remain deficiently expounded (Torresi et al., 2019). The extensive decrease in HCC cases demands a broader range of universal HBV vaccination application and efficient therapy of HBV-relevant chronic hepatitis, which has a long way to go (Sagnelli et al., 2020).

RNA-binding protein (RBP) represents a critical mediator of cancer phenotype (Müller-McNicoll and Neugebauer, 2013). RBP acts dynamically and multifunctionally on multiple levels of post-transcriptional gene expression, such as mRNA splicing, stability and translation (Hou et al., 2019). Over 1,500 human RBPs have been discovered, which possess 600 structure-specific RNA-binding domains (Schneider et al., 2019). Most of them are characterized by evolutionary conservatism and ubiquitous expression to maintain cellular homeostasis (Chabrolles et al., 2020). Genetic and proteomic data highlight that alterations in RBP expression display profound implications on HCC (Lin et al., 2019). For instance, RBP YTHDF2 facilitates cancer stem cell phenotype and metastasis in HCC through regulating OCT4 N⁶-methyladenosine methylation (Zhang et al., 2020). RBP RPS3 leads to hepatocarcinogenesis through up-regulating SIRT1 at a post-transcriptional level (Zhao et al., 2019). RBP RBM3 induces proliferation of HCC cells via regulating circRNA-SCD production (Dong et al., 2019). Despite this, the functions of most RBPs in HBV-related HCC remain still unclear.

This study systematically dissected the prognosis-related RBPs of HBV-related HCC, and established and externally verified an RBP model for predicting prognosis and recurrence by applying least absolute shrinkage selection operator (LASSO). This model might serve as a potential prognostic stratification tool and offer several therapeutic targets for HBV-related HCC.

Materials and Methods

Data Acquisition

RNA-seq transcriptome data and complete clinical information of 374 HCC and 50 control specimens were retrieved from the Cancer Genome Atlas (TCGA) database¹. Among them, 108 HBV-related HCC were extracted from TCGA dataset, which were employed as the discovery set. Expression profiling and clinical features of 224 HBV-related HCC samples were obtained from the GSE14520 dataset in the Gene Expression Omnibus (GEO;²) repository. This GSE14520 dataset was based on the GPL571 and GPL3921 platforms. This dataset was used as the validation set. Table 1 listed the clinical characteristics of HBV-related HCC patients. Each probe was transformed to the corresponding gene symbol. If multiple probes matched the same gene symbol, the average value was determined as the expression value of this gene. Based on published articles, a list of 1,542 RBPs was collected (Supplementary Table 1).

TABLE 1

Table 1. Clinical characteristics of HBV-related HCC patients.

Differential Expression Analysis

Differences in expression levels of RBPs were analyzed between 108 HBV-related HCC and 50 control samples via limma package based on RNA-seq transcriptome data (Ritchie et al., 2015). Under the threshold of | fold-change| > 2 and adjusted p < 0.001, up- and down-regulated RBPs were screened for HBV-related HCC.

Functional Annotation Analysis

Gene ontology (GO) enrichment analysis primarily contains three categories: biological process (BP), cellular component (CC) and molecular function (MF). Furthermore, Kyoto Encyclopedia of Genes and Genomes (KEGG) can provide signaling pathways involved in RBPs. Here, GO, and KEGG enrichment analyses of differentially expressed RBPs were carried out via clusterProfiler package (Ashburner et al., 2000). Terms with false discovery rate (FDR) < 0.05 were significantly enriched by above RBPs. Functional association between HBV-related RBPs was analyzed through the STRING database³ (Szklarczyk et al., 2017). A protein-protein interaction (PPI) network was conducted utilizing Cytoscape software (Doncheva et al., 2019).

Determining Candidate Prognosis-Related RBPs

Univariate Cox regression analysis was applied for analyzing the associations between overall survival (OS) of HBV-related HCC patients and differentially expressed RBPs utilizing survival package. Prognosis-related RBPs were determined with the threshold of p < 0.05. Then, key prognosis-related RBPs were screened through LASSO regression analysis with glmnet package (Engebretsen and Bohlin, 2019). Normalized regression coefficients of key prognosis-related RBPs were calculated based on multivariate regression analysis.

Establishment of an RBP-Related Prognostic Model

A risk score was established on the basis of the key prognosis-related RBPs for HBV-related HCC patients in the TCGA dataset. This study calculated the risk score of each patient, according to the formula: risk score = expression of RBP1 × regression coefficient of RBP1 + expression of RBP2 × regression coefficient of RBP2 + … + expression of RBPn × regression coefficient of RBPn. According to the median risk score, HBV-related HCC patients were assigned into high- and low-risk groups. Hierarchical clustering analysis was applied for showing the associations between expression patterns of above RBPs and clinical characteristics (stage, grade, gender, and age). Kaplan-Meier curve analyses were presented for investigating the differences in OS, disease-free survival (DFS) and progression-free interval (PFI) between two groups using Survival package. P values were determined with log-rank tests. Time-dependent receiver operating characteristic (ROC) curves under one-, three- and five-year OS, DFS, and PFI were generated by SurvivalROC package (Heagerty et al., 2000). Then, the area under the curve (AUC) was calculated to evaluate the predictive usefulness of the risk score.

Prognostic Model Verification

The associations between risk score, age, gender, grade, stage, and prognosis were evaluated utilizing univariate cox regression analysis. Prognostic factors with p < 0.05 were incorporated for multivariate cox regression analysis. Independent prognostic factors were then identified for HBV-related HCC. Time-independent ROCs of risk score, age, gender, grade, and stage were separately constructed and the predictive power was compared. In published literature, Fang and Chen (2020) proposed a two-m⁶A RNA methylation regulator prognostic model (HNRNPA2B1 and RBM15) for HBV-related HCC prognosis. By ROCs, the predictive efficacy was compared with our prognostic model. For evaluating the clinical generalizability of the model, this model was verified in the GSE14520 dataset. With the same formula, the risk scores of patients were calculated in this dataset. Based on the median risk score, OS differences between high- and low-risk groups were assessed by Kaplan-Meier curves and log-rank tests. The predictive performance was verified by ROCs.

Subgroup Analysis

HBV-related HCC subjects in TCGA dataset were stratified into subgroups according to clinical characteristics, as follows: age ≥ 65 and < 65 subgroups, female and male subgroups, grade 1–2 and 3–4 subgroups and stage I-II and stage III-IV subgroups. In each subgroup, Kaplan-Meier curves of OS were conducted between high- and low-risk patients. Survival differences were estimated with log-rank tests.

Nomogram Construction

A nomogram by incorporating gender, age, grade, stage, and risk score was created for predicting one-, three-, and five-year survival probabilities of HBV-related HCC patients in the TCGA dataset. The accuracy in predicting prognosis was evaluated by ROC curves, calibration curves and decision curve analysis (DCA) (Vickers and Elkin, 2006). ROC curves were conducted for evaluating the predictive performance of this nomogram for one-, three- and five-year OS. By calibration curves, discrepancy between nomogram-estimated and actual one-, three-, and five-year survival duration was analyzed. DCA was applied for quantifying the clinical practical use with survival outcomes of a decision considered.

Gene Set Enrichment Analyses (GSEA)

HBV-related HCC specimens were stratified into high- and low-risk groups. By Gene Set Enrichment Analyses (GSEA)⁴, potential mechanisms of this prognosis-related model were elucidated (Subramanian et al., 2005). GSEA was carried out for finding enriched KEGG pathways, with KEGG gene set “c2.cp.kegg.v7.0.symbols.gmt” as a reference. Pathways with nominal p < 0.05 were distinctly enriched.

Estimation of Immune Cell Infiltrations and HLA Expression

By applying CIBERSORT algorithm⁵, the proportions of 22 immune cells in HBV-related HCC and control specimens were quantified based on their expression profiling (Newman et al., 2015). The proportions of all immune cells in each specimen were equal to 1. The LM22 signature was employed as a reference set. The permutations were set as 1,000. Samples with p < 0.05 were screened for further analyses. The correlations between risk score and infiltrations of immune cells and HLA family expression were assessed with Spearson correlation analysis.

Results

Expression and Functions of RBPs in HBV-Related HCC

This study collected 1,542 RBPs and analyzed their expression in 108 HBV-related HCC and 50 control tissues. | Fold-change| > 2 and adjusted p < 0.001 were set as the screening thresholds. As a result, 340 RBPs were up-regulated in HBV-related HCC than control specimens (Figures 1A,B). The first 20 up-regulated RBPs were shown in Table 2. Meanwhile, there were six down-regulated RBPs (GSPT2, PAIP2B, ANG, AZGP1, CPEB3, and ZFP36) in HBV-related HCC compared to controls (Figures 1A,B and Table 2). These RBPs were primarily enriched in mRNA metabolic biological processes such as nuclear-transcribed mRNA catabolic process, RNA catabolic process, nuclear-transcribed mRNA catabolic process, mRNA catabolic process and ncRNA processing (Figure 1C). Our KEGG analysis demonstrated that above RBPs were distinctly related to key pathways including ribosome, spliceosome, RNA transport, mRNA surveillance pathway, ribosome biogenesis in eukaryotes and RNA degradation, confirming their roles on post-transcriptional gene regulation (Figure 1D). A PPI network was conducted for revealing the interactions between these HCC-related RBPs (Figure 1E).

FIGURE 1

Figure 1. Expression patterns and functions of RBPs in HBV-related HCC tissues. (A) Volcano plots of up- and down-regulated RBPs in HBV-related HCC than control specimens. Red dots represented up-regulated RBPs; blue dots represented down-regulated RBPs; gray dots represented no differentially expressed RBPs. (B) Heatmap visualizing expression values of up- (red) and down-regulated (blue) RBPs in HBV-related HCC and control specimens. (C) The first ten biological processes (BPs), cellular components (CCs) and molecular functions (MFs) involved in differentially expressed RBPs. (D) KEGG pathways enriched by above RBPs. The longer the band, the greater the number of enriched genes. The more the color tended to red, the smaller adjusted p-value. (E) A PPI network based on the differentially expressed RBPs. Red indicated up-regulated RBPs and green indicated down-regulated RBPs.

TABLE 2

Table 2. The first 20 up-regulated RBPs and 6 down-regulated RBPs in HBV-related HCC.

Establishing a RBP Model for Predicting HBV-Related HCC Patients’ Prognosis and Recurrence

By univariate analyses, we investigated which dysregulated RBPs were in relation to HBV-related HCC patients’ survival outcomes. With the cutoff of p < 0.05, 54 RBPs were distinctly correlated to prognosis (Table 3). These prognosis-related RBPs were utilized for LASSO analysis. As a result, 10 key RBPs were screened for constructing a prognostic model (Figures 2A,B). The regression coefficients of above RBPs were as follows: POLR2L = 0.453206443450835; MRP S12 = 0.00334358058970182; DYNLL1 = 0.0556454917799665; ZFP36 = 0.215686692231002; PPIH = 0.300574508736105; RARS = 1.0223374037773; SRP14 = −1.08473362012663; DDX 41 = 0.00122492284215762; EIF2B4 = 0.348820298750318; NOL12 = 0.155516694418799. The risk score of each subject from TCGA dataset was determined according to expressions and regression coefficients of RBPs. Heatmap depicted the correlations between expression patterns of POLR2L, MRPS12, DYNLL1, ZFP36, PPIH, RARS, SRP14, DDX41, EIF2B4, and NOL12 and clinical characteristics of HBV-related HCC patients (Figure 2C). Patients from the TCGA dataset were assigned into high- and low-risk groups. As a result, high-risk subjects were indicative of unfavorable OS (p = 3.734e-06; Figure 2D), DFS (p = 6.495e-02; Figure 2E) and PFI (p = 2.135e-02; Figure 2F). The ROCs were plotted for evaluation of the predictive performance. The AUCs under one-, three-, and five-year OS were separately 0.900, 0.945 and 0.886, demonstrating that this model could be precisely predictive of OS probabilities (Figure 2G). Meanwhile, the AUCs under one-, three-, and five-year DFS were 0.686, 0.729 and 0.668 (Figure 2H) as well as the AUCs under one-, three-, and five-year PFI (Figure 2I) were 0.673, 0.784, and 0.667, which were indicative that this model possessed the well performance on predicting HBV-related HCC recurrence and progression.

TABLE 3

Table 3. Prognosis-related RBPs in HBV-related HCC patients by univariate cox regression analysis.

FIGURE 2

Figure 2. Establishing an RBP prognostic model for HBV-related HCC patients in the TCGA dataset. (A) Relationships of λ and partial likelihood deviance. The vertical slash represents the optimal λ value. (B) LASSO regression coefficients of RBPs in this model. (C) Heatmap for the correlations between expression values of RBPs and clinical features. Kaplan-Meier curves of (D) OS, (E) DFS and (F) PFI between high- and low-risk patients. The ROCs of (G) OS, (H) DFS, and (I) PFI based on the risk score.

The RBP Model Displays Independent and Well Predictive Power for HBV-Related HCC Patients

By univariate analyses, we investigated the associations between survival outcomes and risk score and clinical parameters in the TCGA dataset. In Figure 3A, risk score (p < 0.001) and stage (p = 0.009) were both risk factors of HBV-related HCC. Following multivariate analyses, the risk score was independently predictive of survival outcomes (p < 0.00l; Figure 3B). In comparison to other clinical parameters, this risk score exhibited the highest AUC of OS (0.940). This demonstrated that the risk score possessed more excellent predictive performance than other parameters (Figure 3C). Compared to the published prognostic model (HNRNPA2B1 and RBM15), higher AUC value was investigated in our model (Figure 3D). We further evaluated the clinical generalizability of this model. In the GSE14520 dataset, high-risk scores were also indicative of unfavorable survival outcomes (p = 2.999e-03; Figure 3E) and the AUC value was 0.656 (Figure 3F). Collectively, this model displayed the independent and well power in predicting prognosis.

FIGURE 3

Figure 3. Verification of the independency and accuracy of the prognostic model. (A) Univariate and (B) multivariate cox regression analyses of risk score, age, gender, grade, and stage with HBV-related HCC patients in the TCGA dataset. (C) The ROCs for comparing the AUCs of the risk score with age, gender, grade, and stage. (D) The ROCs for comparing the AUCs of the risk score with the model constructed by Fang et al. (E) Survival differences between high- and low-risk patients in the GSE14520 dataset. P value was estimated with log-rank test. (F) The ROC for evaluation of the predictive performance of the risk score.

Subgroup Analysis of This Prognostic RBP Model in HBV-Related HCC

For investigating the predictive sensitivity of this model in HBV-related HCC prognosis, survival analyses were carried out in different subgroups. Our data demonstrated that high-risk scores were distinctly predictive of poorer survival outcomes than low-risk score in age ≥ 65 (p = 0.067; Figure 4A) and < 65 (p < 0.001; Figure 4B) subgroups, female (p = 0.042; Figure 4C) and male (p < 0.001; Figure 4D) subgroups, grade 1–2 (p = 0.027; Figure 4E) and grade 3–4 (p < 0.001; Figure 4F) subgroups and stage I-II (p < 0.001; Figure 4G) and stage III-IV (p = 0.044; Figure 4H) subgroups.

FIGURE 4

Figure 4. Subgroup analysis for assessing the sensibility of the risk score to predict HBV-related HCC patients’ prognosis. Kaplan-Meier curves for survival outcomes of high- and low-risk patients in (A) age ≥ 65 and (B) < 65 subgroups, (C) female and (D) male subgroups, (E) grade 1–2 (G1-2) and (F) grade 3–4 (G3-4) subgroups, (G) stage I–II and (H) stage III–IV subgroups. P values were estimated with log-rank tests.

Constructing a Prognostic Nomogram for HBV-Related HCC

To facilitate personalized treatment, we constructed a nomogram by incorporating gender, age, grade, stage, and risk score in TCGA dataset. This nomogram was utilized for estimating one-, three-, and five-year survival probabilities (Figure 5A). ROC curves demonstrated the well performance on predicting one-, three- and five-year (Figure 5B). As shown in our calibration plots, nomogram-estimated one-, three- and five-year survival probabilities were close to actual survival consequences (Figures 5C–E). Meanwhile, DCA showed that this nomogram exhibited the best net benefit for one-, three- and five-year survival duration (Figures 5F–H). Hence, the nomogram model could assist clinical management and decision.

FIGURE 5

Figure 5. Establishment of the nomogram for HBV-related HCC patients’ prognosis. (A) The nomogram that contained gender, age, grade, stage, and risk score for estimating one-, three- and five-year survival probabilities. (B) The ROC for evaluation of the predictive performance of this nomogram. (C–E) The calibration curves for investigating the discrepancy between nomogram-estimated and actual one-, three- and five-year survival duration. (F–H) The DCA for calculating the clinical net benefit of the nomogram, clinical factor, and prognostic factor in comparison to all or none strategies.

Activated Pathways in High-Risk HBV-Related HCC

For observing potential pathways involved in unfavorable survival outcomes, GSEA was carried out. We found that endocytosis (Figure 6A), RNA degradation (Figure 6B), spliceosome (Figure 6C) and ubiquitin-mediated proteolysis (Figure 6D) were distinctly activated in high-risk HBV-related HCC specimens.

FIGURE 6

Figure 6. Activated KEGG pathways in high-risk HBV-related HCC samples from the TCGA dataset by GSEA. (A) Endocytosis; (B) RNA degradation; (C) spliceosome; (D) ubiquitin-mediated proteolysis.

The Risk Score Is Associated With Immune Microenvironment of HBV-Related HCC

CIBERSORT algorithm was applied for inferring the proportions of 22 immune cells in HBV-related HCC tissues, including B cells naïve, B cells memory, plasma cells, T cells CD8, T cells CD4 memory resting, T cells CD4 memory activated, T cells follicular helper, T cells regulatory (Tregs), T cells gamma delta, NK cells resting, NK cells activated, monocytes, macrophages M0, macrophages M1, macrophages M2, dendritic cells resting, dendritic cells activated, mast cells resting, mast cells activated, eosinophils and neutrophils (Figure 7A). There was the heterogeneity in the immune microenvironment among subjects. The close crosstalk between immune cells was found, as shown in Figure 7B. Also, the risk score was associated with mast cells resting, NK cells resting, neutrophils, mast cells activated, macrophages M0, Tregs and eosinophils. Moreover, the risk score displayed the significant correlations to HLA expression in HBV-related HCC specimens (Figure 7C). The data indicated that the risk score was related to the immune microenvironment of HBV-related HCC.

FIGURE 7

Figure 7. The risk score is related to the immune microenvironment of HBV-related HCC in the TCGA dataset by CIBERSORT. (A) The proportions of immune cells in HBV-related HCC samples. (B) Heatmap of the correlations between the risk score and immune cells. (C) Heatmap visualizing the correlations between the risk score and HLA expression. Correlation coefficients were marked in the box.

Discussion

In this study, an RBP-related gene model was created, which could robustly predict prognosis and recurrence of HBV-related HCC individuals. High risk scores were indicative of undesirable survival outcomes. Our data confirmed RBPs as critical elements in the malignant progression of HBV-related HCC. The gene model acted as a key clinical implication in prognostication and therapy decision.

Alterations in RBP expression may lead to carcinogenesis. Here, we identified 340 dysregulated RBPs in HBV-related HCC. These RBPs were distinctly related to mRNA metabolic biological processes such as nuclear-transcribed mRNA catabolic process, RNA catabolic process, nuclear-transcribed mRNA catabolic process, mRNA catabolic process and ncRNA processing as well as key pathways including ribosome, spliceosome, RNA transport, mRNA surveillance pathway, ribosome biogenesis in eukaryotes and RNA degradation. Thus, these RBPs acted as key regulators on post-transcriptional gene expression.

By LASSO analysis, we created an RBP prognostic model, which contained POLR2L, MRPS12, DYNLL1, ZFP36, PPIH, RARS, SRP14, DDX41, EIF2B4, and NOL12. Following external verification, this model possessed higher accuracy and sensitivity on prognostication in comparison to other clinical parameters. The biological implications of above RBPs in this model have been reported in previous research. Liu et al. (2020) found that POLR2L displayed a correlation to survival duration and alternative splicing in lung squamous cell carcinoma patients. MRPS12 functioned as an oncogene and a prognostic candidate in ovarian carcinoma (Qiu et al., 2021). DYNLL1 hypomethylation and upregulation was characterized by stage- and grade-dependent manners and correlated to unfavorable survival outcomes in HCC (Berkel and Cacan, 2020). ZFP36 down-regulation was detected in HCC tissues and served as a tumor suppressor (Kröhler et al., 2019). PPIH was highly expressed in stomach adenocarcinoma and down-regulated PPIH suppressed cellular migratory and invasive behaviors (Li et al., 2021). Also, PPIH up-regulation contributed to undesirable survival outcomes. DDX41 displayed a correlation to tumor stage and grade in HCC (Qi et al., 2020). NOL12 was in relation to kidney renal clear cell carcinoma prognosis (Xiang Y. et al., 2020). Nevertheless, biological roles and clinical implications of these RBPs may require in-depth exploration in HBV-related HCC.

Our results showed that endocytosis, RNA degradation, spliceosome and ubiquitin-mediated proteolysis pathways were markedly activated in high-risk HBV-related HCC specimens, indicating that these pathways might modulate HCC progress and metastasis. RBPs may regulate the stability of mRNAs encoding immune-related proteins, thereby ornamenting the immune microenvironment (Käfer et al., 2019). For instance, RBP ZFP36, as an inflammatory regulator, restrained T cell activation and anti-viral immunity (Moore et al., 2018). RBP YTHDF3 restrained interferon-dependent anti-viral response through increasing FOXO3 translation (Zhang et al., 2019). Suppressing YTHDF1 in dendritic cells induced durable neoantigen-specific immunity and enhanced the efficacy of anti-PD-L1 therapy (Han et al., 2019). PCBP1 served as an intracellular immune checkpoint toward maintaining T cell functions (Ansa-Addo et al., 2020). Recently, the roles of RBPs have been investigated thoroughly in HCC immune microenvironment of HCC. Xiang J et al. (2020) reported that RBP SIRT7 enhanced the efficacies of anti-PD-L1 therapy through MEF2D in HCC cells. Here, our data demonstrated the close associations of the risk score with immune cells in HCC tissues such as mast cells resting, NK cells resting, neutrophils, mast cells activated, macrophages M0, Tregs and eosinophils. More studies will be carried out for verifying the roles of the risk score on ornamenting immune microenvironment in HCC.

Conclusion

Collectively, this study established an RBP model for predicting OS and recurrence of HBV-related HCC individuals. Following external verification, this model possessed the well predictive efficacy and acted as a robust and specific prognostic indicator. Thus, our findings might assist guide clinical decision and personalized therapies.

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.

Author Contributions

QM conceived and designed the study. ML, ZL, and JW conducted most of the experiments and data analysis, and wrote the manuscript. HL, HG, SL, and MJ participated in collecting data and helped to draft the manuscript. All authors reviewed and approved the manuscript.

Funding

This work was funded by National Science and Technology Major Project Specialized for Infectious Diseases Prevention and Treatment (2017ZX10203201-006 and 2017ZX10202201-0 04-010).

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.2021.707305/full#supplementary-material

Supplementary Table 1 | A list of 1,542 RBPs.

Abbreviations

HCC, hepatocellular carcinoma; HBV, hepatitis B virus; RBP, RNA-binding protein; LASSO, least absolute shrinkage selection operator; TCGA, The Cancer Genome Atlas; GEO, Gene Expression Omnibus; BP, biological process; CC, cellular component; MF, molecular function; KEGG, Kyoto Encyclopedia of Genes and Genomes; FDR, false discovery rate; OS, overall survival; DFS, disease-free survival; PFI, progression-free interval; ROC, receiver operating characteristic; AUC, area under the curve; DCA, decision curve analysis; GSEA, Gene Set Enrichment Analyses.

Footnotes

1. ^ https://portal.gdc.cancer.gov/
2. ^ https://www.ncbi.nlm.nih.gov/gds/
3. ^ http://string-db.org/
4. ^ http://software.broadinstitute.org/gsea/index.jsp
5. ^ http://cibersort.stanford.edu/

References

Ansa-Addo, E. A., Huang, H. C., Riesenberg, B., Iamsawat, S., Borucki, D., Nelson, M. H., et al. (2020). RNA binding protein PCBP1 is an intracellular immune checkpoint for shaping T cell responses in cancer immunity. Sci. Adv. 6:eaaz3865. doi: 10.1126/sciadv.aaz3865

PubMed Abstract | CrossRef Full Text | Google Scholar

Ashburner, M., Ball, C. A., Blake, J. A., Botstein, D., Butler, H., Cherry, J. M., et al. (2000). Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat. Genet. 25, 25–29. doi: 10.1038/75556

PubMed Abstract | CrossRef Full Text | Google Scholar

Berkel, C., and Cacan, E. (2020). DYNLL1 is hypomethylated and upregulated in a tumor stage- and grade-dependent manner and associated with increased mortality in hepatocellular carcinoma. Exp. Mol. Pathol. 117:104567. doi: 10.1016/j.yexmp.2020.104567

PubMed Abstract | CrossRef Full Text | Google Scholar

Chabrolles, H., Auclair, H., Vegna, S., Lahlali, T., Pons, C., Michelet, M., et al. (2020). Hepatitis B virus Core protein nuclear interactome identifies SRSF10 as a host RNA-binding protein restricting HBV RNA production. PLoS Pathog. 16:e1008593. doi: 10.1371/journal.ppat.1008593

PubMed Abstract | CrossRef Full Text | Google Scholar

Doncheva, N. T., Morris, J. H., Gorodkin, J., and Jensen, L. J. (2019). Cytoscape StringApp: network analysis and visualization of proteomics data. J. Proteome Res. 18, 623–632. doi: 10.1021/acs.jproteome.8b00702

PubMed Abstract | CrossRef Full Text | Google Scholar

Dong, W., Dai, Z. H., Liu, F. C., Guo, X. G., Ge, C. M., Ding, J., et al. (2019). The RNA-binding protein RBM3 promotes cell proliferation in hepatocellular carcinoma by regulating circular RNA SCD-circRNA 2 production. EBioMedicine 45, 155–167. doi: 10.1016/j.ebiom.2019.06.030

PubMed Abstract | CrossRef Full Text | Google Scholar

Engebretsen, S., and Bohlin, J. (2019). Statistical predictions with glmnet. Clin. Epigenetics 11:123. doi: 10.1186/s13148-019-0730-1

PubMed Abstract | CrossRef Full Text | Google Scholar

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:60. doi: 10.1186/s10020-020-00185-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Han, D., Liu, J., Chen, C., Dong, L., Liu, Y., Chang, R., et al. (2019). Anti-tumour immunity controlled through mRNA m(6)A methylation and YTHDF1 in dendritic cells. Nature 566, 270–274. doi: 10.1038/s41586-019-0916-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Heagerty, P. J., Lumley, T., and Pepe, M. S. (2000). Time-dependent ROC curves for censored survival data and a diagnostic marker. Biometrics 56, 337–344. doi: 10.1111/j.0006-341x.2000.00337.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Hou, J., Zhang, H., Liu, J., Zhao, Z., Wang, J., Lu, Z., et al. (2019). YTHDF2 reduction fuels inflammation and vascular abnormalization in hepatocellular carcinoma. Mol. Cancer 18:163. doi: 10.1186/s12943-019-1082-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Käfer, R., Schmidtke, L., Schrick, K., Montermann, E., Bros, M., Kleinert, H., et al. (2019). The RNA-Binding Protein KSRP Modulates Cytokine Expression of CD4(+) T Cells. J. Immunol. Res. 2019:4726532. doi: 10.1155/2019/4726532

PubMed Abstract | CrossRef Full Text | Google Scholar

Kröhler, T., Kessler, S. M., Hosseini, K., List, M., Barghash, A., Patial, S., et al. (2019). The mRNA-binding Protein TTP/ZFP36 in Hepatocarcinogenesis and Hepatocellular Carcinoma. Cancers 11:1754. doi: 10.3390/cancers11111754

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, A., Wu, J., Zhai, A., Qian, J., Wang, X., Qaria, M. A., et al. (2019). HBV triggers APOBEC2 expression through miR-122 regulation and affects the proliferation of liver cancer cells. Int. J. Oncol. 55, 1137–1148. doi: 10.3892/ijo.2019.4870

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, J., Zhou, W., Wei, J., Xiao, X., An, T., Wu, W., et al. (2021). Prognostic Value and Biological Functions of RNA Binding Proteins in Stomach Adenocarcinoma. Onco. Targets Ther. 14, 1689–1705. doi: 10.2147/ott.S297973

PubMed Abstract | CrossRef Full Text | Google Scholar

Lin, Y., Liang, R., Qiu, Y., Lv, Y., Zhang, J., Qin, G., et al. (2019). Expression and gene regulation network of RBM8A in hepatocellular carcinoma based on data mining. Aging 11, 423–447. doi: 10.18632/aging.101749

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, Y., Jia, W., Li, J., Zhu, H., and Yu, J. (2020). Identification of Survival-Associated Alternative Splicing Signatures in Lung Squamous Cell Carcinoma. Front. Oncol. 10:587343. doi: 10.3389/fonc.2020.587343

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, Z., Jiang, Y., Yuan, H., Fang, Q., Cai, N., Suo, C., et al. (2019). The trends in incidence of primary liver cancer caused by specific etiologies: results from the Global Burden of Disease Study 2016 and implications for liver cancer prevention. J. Hepatol. 70, 674–683. doi: 10.1016/j.jhep.2018.12.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Moore, M. J., Blachere, N. E., Fak, J. J., Park, C. Y., Sawicka, K., Parveen, S., et al. (2018). ZFP36 RNA-binding proteins restrain T cell activation and anti-viral immunity. Elife 7:e33057. doi: 10.7554/eLife.33057

PubMed Abstract | CrossRef Full Text | Google Scholar

Müller-McNicoll, M., and Neugebauer, K. M. (2013). How cells get the message: dynamic assembly and function of mRNA-protein complexes. Nat. Rev. Genet. 14, 275–287. doi: 10.1038/nrg3434

PubMed Abstract | CrossRef Full Text | Google Scholar

Nakagawa, H., Fujita, M., and Fujimoto, A. (2019). Genome sequencing analysis of liver cancer for precision medicine. Semin. Cancer Biol. 55, 120–127. doi: 10.1016/j.semcancer.2018.03.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Newman, A. M., Liu, C. L., Green, M. R., Gentles, A. J., Feng, W., Xu, Y., et al. (2015). Robust enumeration of cell subsets from tissue expression profiles. Nat. Methods 12, 453–457. doi: 10.1038/nmeth.3337

PubMed Abstract | CrossRef Full Text | Google Scholar

Qi, Z., Yan, F., Chen, D., Xing, W., Li, Q., Zeng, W., et al. (2020). Identification of prognostic biomarkers and correlations with immune infiltrates among cGAS-STING in hepatocellular carcinoma. Biosci. Rep. 40:BSR20202603. doi: 10.1042/bsr20202603

PubMed Abstract | CrossRef Full Text | Google Scholar

Qiu, X., Guo, D., Du, J., Bai, Y., and Wang, F. (2021). A novel biomarker, MRPS12 functions as a potential oncogene in ovarian cancer and is a promising prognostic candidate. Medicine 100:e24898. doi: 10.1097/md.0000000000024898

PubMed Abstract | CrossRef Full Text | Google Scholar

Ritchie, M. E., Phipson, B., Wu, D., Hu, Y., Law, C. W., Shi, W., et al. (2015). limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 43:e47. doi: 10.1093/nar/gkv007

PubMed Abstract | CrossRef Full Text | Google Scholar

Sagnelli, E., Macera, M., Russo, A., Coppola, N., and Sagnelli, C. (2020). Epidemiological and etiological variations in hepatocellular carcinoma. Infection 48, 7–17. doi: 10.1007/s15010-019-01345-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Schneider, T., Hung, L. H., Aziz, M., Wilmen, A., Thaum, S., Wagner, J., et al. (2019). Combinatorial recognition of clustered RNA elements by the multidomain RNA-binding protein IMP3. Nat. Commun. 10:2266. doi: 10.1038/s41467-019-09769-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Subramanian, A., Tamayo, P., Mootha, V. K., Mukherjee, S., Ebert, B. L., Gillette, M. A., et al. (2005). Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc. Natl. Acad. Sci. U. S. A. 102, 15545–15550. doi: 10.1073/pnas.0506580102

PubMed Abstract | CrossRef Full Text | Google Scholar

Szklarczyk, D., Morris, J. H., Cook, H., Kuhn, M., Wyder, S., Simonovic, M., et al. (2017). The STRING database in 2017: quality-controlled protein-protein association networks, made broadly accessible. Nucleic Acids Res. 45, D362–D368. doi: 10.1093/nar/gkw937

PubMed Abstract | CrossRef Full Text | Google Scholar

Torresi, J., Tran, B. M., Christiansen, D., Earnest-Silveira, L., Schwab, R. H. M., and Vincan, E. (2019). HBV-related hepatocarcinogenesis: the role of signalling pathways and innovative ex vivo research models. BMC Cancer 19:707. doi: 10.1186/s12885-019-5916-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Vickers, A. J., and Elkin, E. B. (2006). Decision curve analysis: a novel method for evaluating prediction models. Med. Decis. Making 26, 565–574. doi: 10.1177/0272989x06295361

PubMed Abstract | CrossRef Full Text | Google Scholar

Xiang, J., Zhang, N., Sun, H., Su, L., Zhang, C., Xu, H., et al. (2020). Disruption of SIRT7 Increases the Efficacy of Checkpoint Inhibitor via MEF2D Regulation of Programmed Cell Death 1 Ligand 1 in Hepatocellular Carcinoma Cells. Gastroenterology 158, 664–678.e24. doi: 10.1053/j.gastro.2019.10.025

PubMed Abstract | CrossRef Full Text | Google Scholar

Xiang, Y., Zhou, S., Hao, J., Zhong, C., Ma, Q., Sun, Z., et al. (2020). Development and validation of a prognostic model for kidney renal clear cell carcinoma based on RNA binding protein expression. Aging 12, 25356–25372. doi: 10.18632/aging.104137

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, C., Huang, S., Zhuang, H., Ruan, S., Zhou, Z., Huang, K., et al. (2020). YTHDF2 promotes the liver cancer stem cell phenotype and cancer metastasis by regulating OCT4 expression via m6A RNA methylation. Oncogene 39, 4507–4518. doi: 10.1038/s41388-020-1303-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Y., Wang, X., Zhang, X., Wang, J., Ma, Y., Zhang, L., et al. (2019). RNA-binding protein YTHDF3 suppresses interferon-dependent antiviral responses by promoting FOXO3 translation. Proc. Natl. Acad. Sci. U. S. A. 116, 976–981. doi: 10.1073/pnas.1812536116

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, L., Cao, J., Hu, K., Wang, P., Li, G., He, X., et al. (2019). RNA-binding protein RPS3 contributes to hepatocarcinogenesis by post-transcriptionally up-regulating SIRT1. Nucleic Acids Res. 47, 2011–2028. doi: 10.1093/nar/gky1209

PubMed Abstract | CrossRef Full Text | Google Scholar