Integrated Analysis of Nine Prognostic RNA-Binding Proteins in Soft Tissue Sarcoma

Abstract

RNA-binding proteins (RBPs) have been shown to be dysregulated in cancer transcription and translation, but few studies have investigated their mechanism of action in soft tissue sarcoma (STS). Here, The Cancer Genome Atlas (TCGA) and Genotype-Tissue Expression (GTEx) databases were used to identify differentially expressed RBPs in STS and normal tissues. Through a series of biological information analyses, 329 differentially expressed RBPs were identified. Functional enrichment analysis showed that differentially expressed RBPs were mainly involved in RNA transport, RNA splicing, mRNA monitoring pathways, ribosome biogenesis and translation regulation. Through Cox regression analyses, 9 RBPs (BYSL, IGF2BP3, DNMT3B, TERT, CD3EAP, SRSF12, TLR7, TRIM21 and MEX3A) were all up-regulated in STS as prognosis-related genes, and a prognostic model was established. The model calculated a risk score based on the expression of 9 hub RBPs. The risk score could be used for risk stratification of patients and had a high prognostic value based on the receiver operating characteristic (ROC) curve. We also established a nomogram containing risk scores and 9 key RBPs to predict the 1-year, 3-year, and 5-year survival rates of patients in STS. Afterwards, methylation analysis showed significant changes in the methylation degree of BYSL, CD3EAP and MEX2A. Furthermore, the expression of 9 hub RBPs was closely related to immune infiltration rather than tumor purity. Based on the above studies, these findings may provide new insights into the pathogenesis of STS and will provide candidate biomarkers for the prognosis of STS.

Introduction

Soft tissue sarcoma (STS) accounts for less than 1% of all cancers but is highly heterogeneous in terms of anatomical location, histology, molecular characteristics, and prognosis (1). STS originates from mesenchyme and can occur anywhere in the body, with approximately 50% to 60% of cases happening in the extremities (2). In the past few decades, through the use of multidisciplinary methods to control diseases, including surgery, radiation therapy and systemic therapy, the treatment effect among STS patients has greatly improved (3). However, due to its easy metastasis and recurrence, the prognosis of advanced patients is still poor. At present, the diagnosis of STS mainly depends on ultrasound scanning, imaging examination and tissue biopsy, which have difficulty meeting the clinical requirements (4). To reduce the recurrence rate and mortality, enhance quality of life and improve survival rates in STS patients, early detection, diagnosis, and treatment are essential. Over the years, molecular research has shown excellent results in understanding cellular and molecular mechanisms that drive tumorigenesis and progression (5). Therefore, it is necessary to find genes for the occurrence, development and prognosis of STS to further broaden the research.

RNA binding proteins (RBPs) are involved in various RNA metabolism and multiple biological processes, and play an important role in regulating RNA stability, splicing, modification, localization and protein translation (6, 7). RBPs are mainly involved in posttranscriptional control (8), and their disorders are closely related to the occurrence and development of cancer (9–11). There are numerous human RBPs, but the role of RBPs in cancer is not fully known. Only a few experiments have systematically analyzed the RBPs associated with cancer. For example, RBMS2 plays a role in suppressing cancer in breast cancer and actively regulates the expression of P21 by stabilizing its mRNA (12). The overexpression of ESRP1 in ovarian cancer promotes the transformation of ovarian cancer cells from mesenchymal phenotype to epithelial phenotype (13). In addition, some studies focusing on the comprehensive analysis of RBPs found that DCAF13, EZR, MRPL13, APOBEC3C and EIF4E3 may have a prognostic value in breast cancer (14). It was reported that PUM2 could partly and competitively bind to STARD13 3’UTR with miRNAs to achieve the effect of inhibiting osteosarcoma progression (15). RBM10 was also demonstrated to be tumor suppressor in osteosarcoma (16). However, few studies have investigated the different roles of RBPs in STS. Therefore, based on the above evidence, RBPs may also have a regulatory effect on STS. This provides a new direction for our study of STS.

The methylation and the resulting abnormal expression of target gene play an important role in the tumorigenesis in multiple tumors including STS (17, 18). It was reported that the methylation status was also associated with prognosis in STS (19). Furthermore, immune infiltration has the vital impact on tumor biological behaviors through activating or suppressing innate and adaptive immunity (20, 21). Hence, in this study, we investigated the relationship between STS and RBPs using bioinformatics methods to better investigate the pathological mechanisms of RBP in STS. We also analyzed the signaling pathways, potential functions, prognostic value, gene alterations, methylation levels and immune infiltration of RBPs in STS. These results not only helped us discover new genes related to RBPs for diagnosis but also increased our understanding of STS to further define prognosis and screen new drug targets to improve the treatment of STS.

Materials and Methods

Data Collection and Processing

Corresponding gene expression profiles and clinical data were obtained from The Cancer Genome Atlas database (TCGA) and Genotype-Tissue Expression (GTEx) database using the University of California at Santa Cruz (UCSC, https://genome.ucsc.edu) Genome Website. RNA-sequencing data (FPKM values) of STS in TCGA-SARC cohort (https://xenabrowser.net/datapages/?cohort=GDC%20TCGA%20Sarcoma%20(SARC)&removeHub=https%3A%2F%2Fxena.treehouse.gi.ucsc.edu%3A443) and normal fat tissue samples in GTEx database (https://xenabrowser.net/datapages/?cohort=GTEX&removeHub=https%3A%2F%2Fxena.treehouse.gi.ucsc.edu%3A443) were included in this study. After downloading gene expression datasets, ENSEMBLE identifiers were matched to official gene symbols. We created a unified processing before the two sets of data were merged for analysis. Because the obtained STS gene expression profiles have been transformed by log2(x+1), the normal fat tissue gene expression data were also manipulated by using log2(x+1) transformation. Then, two gene files were merged and standardized, and the expression matrix of RBPs was extracted for differential analysis. Ultimately, differentially expressed RBPs were identified based on |logfold change (FC)| > 0.5 and false discovery rate (FDR) < 0.05, using Wilcoxon test in R software. Wilcoxon test is one of the most common nonparametric test methods (22).

Functional Enrichment Analyses

Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were performed to find the biological functions of these differentially expressed RBPs. GO categories were comprised of three terms: molecular function (MF), biological process (BP), and cellular component (CC). KEGG analysis was used to find the key pathways that these genes were involved in STS. All enrichment analyses were performed based on background genes from GO or KEGG database using the “clusterProfiler” package (23) in R software. Both P and FDR values less than 0.05 were the threshold for the identification of significant GO terms and KEGG pathways.

Protein–Protein Interaction (PPI) Network Construction and Module Screening

Differentially expressed RBPs were uploaded to the STRING database (https://srting-db.org/cgi/input.pl) (24) to obtain a protein-protein interaction (PPI) network. Subsequently, Cytoscape3.7.0 software was used to study PPI network. We used Centiscape 2.2 to obtain centrality indexes of nodes of network. Molecular Complex Detection (MCODE) (25) was used to find and visualize the top 3 clusters of the PPI network. Selection criteria were as follows: degree cutoff = 2, node score cutoff = 0.2, K - Core = 2. After that, the shared genes between the PPI network and the differentially expressed RBPs were extracted for subsequent model construction.

Prognostic Model Construction

All samples were randomly divided into 2 groups (training set and test set) with the ratio of 1:1. The training group was used to construct a prognostic model, while the test group was used to verify the result of the training group. Since there was little clinical information in the training and testing groups, we also validated the prognostic model in the whole population. First, univariate Cox regression analysis of all samples was performed to screen differentially expressed RBPs (P < 0.01 was considered statistically significant). Then, multivariate Cox regression analysis was performed to find the optimal genes to construct a prognostic model in the training group. These optimal genes were called as hub genes. Risk score for STS patients was calculated based on the following formula: risk score = where i represented each gene, n represented the total number of hub genes in the multivariate Cox regression analysis, βi represented gene coefficient value and Expi represented the expression value of each gene. According to the median risk score, STS patients were divided into a high-risk group and a low-risk group. The Kaplan-Meier survival curve was drawn to compare the survival difference between the two groups in the training set. In addition, we used the “survivalROC” package to draw a receiver operation curve (ROC) (26) to evaluate the predictive ability of the prognostic model. To further validate our results, the prognostic model was also used to calculate the risk score for each patient in the test group and the whole cohort. Similarly, Kaplan-Meier survival curve and the ROC curve were plotted to evaluate the predictive ability of the model in both the validation and whole cohorts. P < 0.05 was considered statistically significant for survival differences between the high and low-risk groups.

Identification and Analysis of Independent Prognostic Factors

Univariate and multivariate Cox regression analyses were used to evaluate the effect of multiple clinical factors (age, sex, metastasis, radiotherapy, and risk score) on the prognosis of STS patients to examine the accuracy of the STS prognostic model and determine independent prognostic factors (P < 0.05 in both analyses). Finally, we used the “rms” package in the R software to construct a nomogram to accurately predict each patient’s survival rate.

Alterations of Hub RBPs

Rates of hub RBP changes in STS were carried out by querying the online cBioPortal (http://cbioportal.org/), which provides complex tumor genomics (27). Through this website, we explored the impact of hub RBP changes on the survival of STS.

Methylation Analyses

Human Disease Methylation Database Version 2.0 (DiseaseMeth2.0 (http://biobigdata.hrbmu.edu.cn/diseasemeth/) was used to compare the methylation degree of hub RBPs between STS and normal tissue. This database collects and analyzes DNA methylation data in various human diseases (28, 29). In addition, we also used MEXPRESS (http://mexpress.be) (30) to study the relationship between hub RBPs and their methylation status.

Immune Infiltration Analysis

There is evidence that tumor cell immune infiltration is closely related to prognosis (31). For example, in LUSC, T follicular helper cells are associated with a good prognosis, and an increase in neutrophils indicates a poor prognosis (32). Tumor immune estimation resource (TIMER) is a comprehensive resource for systematic analysis of various malignant tumors. TIMER contains molecular characterization of 32 types of cancer and 6 types of immune cells (33). Based on the TIMER database, we further evaluated the relationship between immune cell types (CD4+T cells, CD8+T cells, B cells, dendritic cells, macrophages, neutrophils and tumor purity) and 9 hub RBPs in STS.

Potential Interactions of Hub RBPs

RNA Interactome Database (RNAInter) (http://www.rna-society.org/rnainter/ or http://www.rna-society.org/raid/.) (34) is a platform to investigate RNA associated interactions, which involve many physiological and pathological processes. We used RNAInter to explore the potential targets of 9 hub RBPs and RBPs-associated interactions to facilitate the understanding of 9 hub RBPs and further research of RBPs on STS.

Results

The Identification of RBPs in STS Patients

In this study, 263 tumor samples and 517 normal tissue samples were downloaded from TCGA and GTEx databases, respectively. Only 259 patients had survival data. There were seven histological types in 259 patients and specific information was presented in Supplementary File 1. A total of 1492 RBPs were extracted for the identification of differentially expressed RBPs between the STS group and normal fat tissue group (Supplementary File 2). According to the screening criteria (FDR < 0.05, |logFC|>0.5), 329 RBPs were obtained in this study, including 191 upregulated and 138 downregulated (Figure 1). The detailed information of 329 RBPs was summarized in Supplementary File 3.

Figure 1

GO Terms and KEGG Pathway Analysis of the Differentially Expressed RBPs

To elucidate the potential biological functions and mechanisms related to the differentially expressed RBPs, GO and KEGG functional analyses were performed. In GO analysis, downregulated genes were enriched in RNA splicing, nucleic acid transport, RNA localization, cytoplasmic ribonucleoprotein granule, nuclear speck, and ribonuclease activity (Figures 2A, B), and upregulated genes were enriched in ncRNA processing, RNA modification, tRNA metabolic process, ribonucleoprotein granule, mRNA 3’−UTR binding, and methyltransferase activity (Figures 2C, D). In addition, KEGG analysis of these genes showed that they were mainly involved in the RNA transport, mRNA surveillance pathway, spliceosome and ribosome biogenesis in eukaryotes (Figures 2E–H).

Figure 2

Protein-Protein Interaction (PPI) Network Construction and Key Module Selection

To explore the interaction relationship between the differentially expressed RBPs, we constructed a PPI network containing 2110 edges in the STRING database (Figure 3A, Supplementary File 4). These results were further analyzed and visualized by Cytoscape, and we obtained a PPI network with 302 nodes and 2075 edges (Figure 3B). Supplementary File 5 showed centrality indexes of nodes in PPI network. To explore the hub RBPs in STS, the PPI network was analyzed by MCODE in Cytoscape, and three important modules were selected (Figures 3C–E and Supplementary File 6). Module 1 contained 21 points and 209 edges, Module 2 contained 27 points and 190 edges, and module 3 contained 44 points and 236 edges (Figures 3C–E and Supplementary File 6). The genes from three modules were visualized in a network (Figure 3F). Finally, 272 shared genes (Supplementary File 3) between the PPI network and 329 differentially expressed RBPs were extracted for subsequent analyses.

Figure 3

Prognosis-Related RBP Selection

By combining gene expression and survival information, 259 patients were used for subsequent model construction. The 259 patients were randomly divided into two groups: training set (N = 131) and test set (N = 128). To select RBPs associated with overall survival (OS), univariate Cox regression analysis was performed on 272 RBPs in 259 patients. The results showed that 23 differentially expressed RBPs were significantly correlated with OS (P < 0.01, Figure 4). Then, the above 23 RBPs were used for further multivariate Cox regression analysis, and finally 9 RBPs highly correlated with STS survival were obtained to establish a prognostic model in the training group (Figure 5).

Figure 4

Figure 5

Construction and Evaluation of a Prognostic Model Associated With RBPs

We obtained 9 RBPs to calculate a risk score and construct a prognostic model. The patient Risk Assessment formula was as follows (Table 1):

According to the median risk score, 131 patients in the training group from TCGA database were divided into the high-risk group and the low-risk group. Survival analysis showed that high-risk patients had shorter OS (P = 6.404e-06, Figure 6A). Next, the area under the ROC curve (AUC) was 0.813, indicating that this model had high predictive ability in the training group (Figure 6D). Similarly, the risk score could also predict the survival rate (P = 5.679e-03 and P = 1.545e-07, respectively; Figures 6B, C) of patients and had an accurate predictive ability (AUC = 0.668 and AUC = 0.726, respectively; Figures 6E, F) in the test group and the whole cohort, respectively. These results suggested that the risk score can predict patient survival well, and the model had good sensitivity and specificity. Moreover, the expression of hub RBPs was also significantly different between high- and low-risk groups (Figures 6G–I).

Figure 6

Identification of Independent Prognostic Factors and Construction of the Nomogram

As shown in Figure 7A, the univariate Cox regression analysis showed that factors including age, metastasis and risk score affected the survival of STS patients (P < 0.05). Subsequently, multivariate Cox regression analysis further demonstrated that metastases and risk scores were independent prognostic factors (P < 0.05, Figure 7B). Based on the above regression analyses, the risk score results further confirmed the accuracy of the prognostic model. Finally, we constructed a nomogram to predict 1-year, 3-year, and 5-year survival rates of patients (Figure 7C).

Figure 7

Relationship Between Hub RBP Alterations and OS

Here, cBioPortal was used to study the hub RBPs in 261 STS samples. Figures 8A, B showed that alterations of hub RBPs were found in 57% (148/261) of the patients. Alterations in these RBPs led to a poor prognosis (P < 0.05, Figure 8C).

Figure 8

Methylation Analysis of Key RBPs

In this study, we used the DiseaseMeth version 2.0 to further analyze the methylation of these 9 differentially expressed RBPs. Figures 9A–C showed significant differences in methylation of BYSL (P = 1.803e-03), CD3EAP (P = 5.989e-05) and MEX3A (P = 2.520e-02). Furthermore, using MEXPRESS analysis, we observed that multiple methylation sites in the RBP sequence were negatively correlated with their own expression levels (P < 0.05, Supplementary Figure 1).

Figure 9

Hub RBP Immune Infiltration Analysis

Given that immune cells are involved in the composition of the tumor microenvironment, which is of great significance to the prognosis of the tumor, we studied the potential connection between 9 hub RBPs and immune infiltration (purity, B cell, CD8+T cell, CD4+T cell, macrophage, neutrophil, dendritic cell) in the STS. Figure 10A showed that TLR7 was positively correlated with B cells (partial. cor = 0.336, P = 1.02e-07), CD8+ T cells (partial. cor = 0.328, P = 1.98e-07), CD4+ T cells (partial. cor = 0.5, P = 1.49e-16), macrophages (partial. cor = 0.622, P = 1.56e-26), neutrophils (partial. cor = 0.571, P = 2.28e-22) and dendritic cells (partial. cor = 0.564, P = 1.18e-21). Similarly, TRIM21 and IGF2BP3 were also positively correlated with most immune cells (Figures 10I, E). BYSL and MEX3A were negatively associated with most immune cells (Figures 10H, D). Figures 10A–C indicated that tumor purity was only significantly correlated with TLR7 (cor = -0.456, P = 5.8e-14), TERT (cor=-0.185, P = 3.66e-03) and SRSF12 (cor = 0.28, P = 8.31e-06). In addition, we found that DNMT3B (partial. cor = -0.179, P = 5.63e-03) and CD3EAP (partial. cor = -0.245, P = 1.33e-04) were only negatively associated with CD4+ T cells (Figures 10F, G).

Figure 10

Identification of Potential Interactions of Hub RBPs

We displayed potential interactions of RBPs using “Exact Search” function of RNAInter. Supplementary Figure 2 showed the top 100 potential interactions of RBPs. BYSL, CD3EAP had more interactions with transcription factors and RBPs. DNMT3B, IGF2BP3, MEX3A, SRSF12 and ILR7 had more interactions with microRNA and RBPs. TERT and TRIM21 had more interactions with microRNA, transcription factors and RBPs.

Discussion

Currently, accurate diagnosis and prediction of biological behavior have become major challenges because of the STS rarity and complexity. In the course of tumor development, molecular and cellular components are considered to be potential prognostic factors (35). Therefore, finding new treatment targets for STS and elucidating their mechanisms can improve patient prognosis. Recently, bioinformatics methods, such as machine learning and weighted gene co-expression network analysis, have been used in the identification of diagnostic or prognostic biomarkers, as well as molecular subtypes to better understand STS (36, 37). Some of biomarkers have also been validated by basic experiments. Nevertheless, further research is needed on STS. In addition, there were numerous studies about RBPs in cancers (10, 38, 39), but few studies of RBPs in STS. In this study, TCGA and GTEx databases were used to identify potential prognostic genes associated with RBPs to explore the pathways and functions of the RBPs involved in STS. First, we obtained 329 differentially expressed RBPs. Among them, CDC5L (40) and IGF2BP1 (41) have been indicated to be biomarkers of osteosarcoma and affect the prognosis of osteosarcoma. To elucidate the underlying biomolecular mechanism of differentially expressed RBPs, we used GO and KEGG enrichment analyses. Then, the PPI network was constructed by STRING and visualized by Cytoscape. Finally, through univariate and multivariate Cox regression analyses, we obtained 9 RBPs and a prognostic model with high prognostic value. To further explore the role of these 9 RBPs in STS, we further performed the alternation, methylation degree, immune infiltration analyses and identified potential interactions of RBPs. These findings may help to identify RBPs associated with STS and be used for the diagnosis and treatment of STS.

At present, the research of RBPs mainly focuses on posttranscriptional events, which are extremely related to the occurrence of tumor. In this study, GO biological process analysis showed that the differentially expressed RBPs mainly focused on RNA splicing, ribonucleoprotein granule, cytoplasmic ribonucleoprotein granule and nuclease activity. Similar to previous studies, these pathways affected the progression of human diseases (11, 42, 43). RNA splicing is a very important biological process in the human body, and splicing dysregulation is involved in the pathogenesis of various cancers (44). This mechanism has been reported to affect the prognosis and biological behavior of STS (45). RNA modification is a dynamic regulatory process, mainly involved in biological processes such as cell differentiation, and its dysregulation is one of the causes of cancer (46). Translation disorder is common in cancer (47). Both carcinogens and tumor suppressor factors affect the mechanism of translation, making translation abnormalities common in cancer. In addition, KEGG analysis revealed mRNA monitoring pathways, ribosomal biogenesis, and ribosomal lysis pathways that have been reported to be associated with cancer. Previous studies have suggested that ribosomal proteins affect tumor development by regulating the P53 pathway and mRNA translation (48, 49). Overall, GO and KEGG analyses suggested that these RBPs were closely related to the progress of STS through the abovementioned pathways.

To identify the key RBPs, 272 shared genes were extracted from the PPI network and 329 differentially expressed RBPs for further analysis. Then, 9 key RBPs (BYSL, CD3EP, DNMT3B, IGF2BP3, MAX3A, SRSF12, TERT, TLR7 and TRIM21) were obtained using univariate and multivariate Cox regression analyses. They have been reported to participate in the progression and play a key role in the prognosis of different cancers, such as lung cancer, hepatocellular carcinoma, breast cancer, rhabdomyosarcoma, glioma, etc. BYSL is a nucleolar protein involved in the biogenesis of ribosomes through 18s rRNA processing in mammals to affect cell proliferation (50). It is positively expressed in hepatocellular carcinoma (51), prostate cancer (52) and ovarian cancer (53). CD3EAP (also known as ASE1, CAST, ERCC1 antisense) is a component of the RNA polymerase I complex, which is involved in ribosome biosynthesis and promotes cell proliferation similar to BYSL (54). CD3EAP has been reported to be involved in the occurrence and metastasis of lung cancer tissues (55). Although BYSL and CD3EAP play important roles in human cancers, their specific mechanisms in STS remain unclear and need further exploration. DNMT3B, a DNA methyltransferase responsible for de novo methylation during human development (56), was significantly elevated in osteosarcoma and rhabdomyosarcoma (57, 58). IGF2BP3 is an embryonic protein that belongs to the mRNA binding protein family and is considered an oncogene (59). It is re-expressed in some cancers and promotes the invasion and migration of tumor cells. In Ewing sarcoma, IGF2BP3 serves as an effective indicator of poor prognosis and predicts the recurrence of patients (60). TERT is a telomerase reverse transcriptase that participates in various human activities through telomere elongation and has enzyme activity in tumors (61). Past studies have confirmed that TERT is upregulated and altered in most cancer patients (62). Among STS, patients with TERT alternation have a shorter lifespan (63). TLR7 is a member of the Toll-like receptor family and plays a key role in the proliferation of tumor cells. It has been shown to be highly expressed in pancreatic cancer (64) and lung cancer (65). TRIM21 is a Fc receptor that can bind IgG, IgM and IgA (66). It has different expressions in different cancers, such as low expression in diffuse large B-cell lymphoma and breast cancer (67, 68), and high expression in human glioma and nasopharyngeal carcinoma cells (69, 70). Overexpression of TRIM21 makes osteosarcoma more resistant to various stresses and promotes its proliferation (71). In our study, TRIM21 expression was one of the reasons for poor STS prognosis. However, there are few studies on these 9 RBPs. In particular, it is difficult to find disease-related studies on MEX3A and SRSF12. Furthermore, we found that 9 RBPs had close interactions with transcription factors, microRNA and other RBPs. Previous research demonstrated that RBPs had similar patterns of somatic copy number analysis to transcription factors across 15 human cancer (72), which could support these findings. Therefore, these results might promote further research of RBPs in STS.

In addition, we compared the prognostic ability of these nine RBPs and clinical factors for STS. Univariate and multivariate Cox regression analyses showed that tumor metastasis and risk scores were independent prognostic factors. To accurately predict patient survival, we constructed a nomogram containing 9 RBPs. Previous nomograms have been developed to predict the OS of patients with STS based on tumor size, tumor grade, histological subtype and complete surgical resection (73). However, there are few nomograms for studying genes. Our study integrated the risk score and expression of 9 RBPs in the nomogram to further accurately predict the survival of patients with STS.

Subsequently, we further analyzed the changes of hub RBPs in STS and found that the alterations of key RBPs in STS patients are related to the shorter survival time, which further suggested that these 9 RBPs had high prognostic ability. In addition, we found that multiple methylation sites of upregulated genes were negatively correlated with their own expression in STS. These data showed that methylation changes may cause abnormal gene expression. The relationship between methylation of these RBPs and STS has rarely been studied before, and thus our study may provide a new direction for the treatment of STS.

Studies have shown that immune infiltration is an important factor affecting the prognosis of STS (74). We discussed the relationship between these 9 hub RBPs and immune infiltration. In this study, we found that most of these 9 hub RBPs have no correlation with tumor purity. Among them, TLR7 and TRIM21 were positively correlated with most immune cell infiltration, and MEX3A and SRSF12 were negatively correlated with most immune cell infiltration. DNMT3B was negatively correlated with CD4+ T cells, while TERT was positively correlated with CD4+ T cells. Some studies showed that TLR7 interacts with miR-25-3p, thereby stimulating the secretion of IL6 by macrophages and promoting the growth and spread of liposarcoma (75). Activation of TRIM21 expression induced infiltration of most immune cells except macrophages (76). IGF2BP3 was associated with immune cell infiltration in patients with liver cancer and affected prognosis (77). Activation of CD3EAP may cause T cells to produce IL2 to play an important role in non-small cell lung cancer (78). Therefore, the abnormal expression of these RBPs may alter the immune infiltrating components in STS, thereby achieving the purpose of treatment and improving the survival of patients.

Overall, based on the above series of analyses, the expression, participation mechanism and prognostic effects of the 9 hub RBPs were studied. These RBPs can be used as new biomarkers to predict the prognosis of STS and provide new ideas for the diagnosis and treatment of STS in the future. As far as we know, this was the first report to establish a prognostic model for RBPs related to STS. In short, our results will greatly facilitate the development of new diagnostic and therapeutic strategies and provide directions for future research.

References

• 1

  JemalASiegelRXuJWardE. Cancer statistics, 2010. CA Cancer J Clin (2010) 60(5):277–300. doi: 10.3322/caac.20073

  • CrossRef
  • Google Scholar

• 2

  von MehrenMRandallRLBenjaminRSBolesSBuiMMCasperESet al. Soft tissue sarcoma, version 2.2014. J Natl Compr Canc Netw (2014) 12(4):473–83. doi: 10.6004/jnccn.2014.0053

  • CrossRef
  • Google Scholar

• 3

  BourcierKLe CesneATselikasLAdamJMirOHonoreCet al. Basic Knowledge in Soft Tissue Sarcoma. Cardiovasc Intervent Radiol (2019) 42(9):1255–61. doi: 10.1007/s00270-019-02259-w

  • CrossRef
  • Google Scholar

• 4

  DangoorASeddonBGerrandCGrimerRWhelanJJudsonI. UK guidelines for the management of soft tissue sarcomas. Clin Sarcoma Res (2016) 6:20. doi: 10.1186/s13569-016-0060-4

  • CrossRef
  • Google Scholar

• 5

  KaliaM. Biomarkers for personalized oncology: recent advances and future challenges. Metabolism (2015) 64(3 Suppl 1):S16–21. doi: 10.1016/j.metabol.2014.10.027

  • CrossRef
  • Google Scholar

• 6

  HentzeMWCastelloASchwarzlTPreissT. A brave new world of RNA-binding proteins. Nat Rev Mol Cell Biol (2018) 19(5):327–41. doi: 10.1038/nrm.2017.130

  • CrossRef
  • Google Scholar

• 7

  PerronGJandaghiPSolankiSSafisamghabadiMStorozCKarimzadehMet al. A General Framework for Interrogation of mRNA Stability Programs Identifies RNA-Binding Proteins that Govern Cancer Transcriptomes. Cell Rep (2018) 23(6):1639–50. doi: 10.1016/j.celrep.2018.04.031

  • CrossRef
  • Google Scholar

• 8

  CalabrettaSRichardS. Emerging Roles of Disordered Sequences in RNA-Binding Proteins. Trends Biochem Sci (2015) 40(11):662–72. doi: 10.1016/j.tibs.2015.08.012

  • CrossRef
  • Google Scholar

• 9

  NeelamrajuYGonzalez-PerezABhat-NakshatriPNakshatriHJangaSC. Mutational landscape of RNA-binding proteins in human cancers. RNA Biol (2018) 15(1):115–29. doi: 10.1080/15476286.2017.1391436

  • CrossRef
  • Google Scholar

• 10

  PereiraBBillaudMAlmeidaR. RNA-Binding Proteins in Cancer: Old Players and New Actors. Trends Cancer (2017) 3(7):506–28. doi: 10.1016/j.trecan.2017.05.003

  • CrossRef
  • Google Scholar

• 11

  WurthL. Versatility of RNA-Binding Proteins in Cancer. Comp Funct Genomics (2012) 2012:178525. doi: 10.1155/2012/178525

  • CrossRef
  • Google Scholar

• 12

  SunXHuYWuJShiLZhuLXiPWet al. RBMS2 inhibits the proliferation by stabilizing P21 mRNA in breast cancer. J Exp Clin Cancer Res (2018) 37(1):298. doi: 10.1186/s13046-018-0968-z

  • CrossRef
  • Google Scholar

• 13

  JeongHMHanJLeeSHParkHJLeeHJChoiJSet al. ESRP1 is overexpressed in ovarian cancer and promotes switching from mesenchymal to epithelial phenotype in ovarian cancer cells. Oncogenesis (2017) 6(10):e389. doi: 10.1038/oncsis.2017.87

  • CrossRef
  • Google Scholar

• 14

  WangKLiLFuLYuanYDaiHZhuTet al. Integrated Bioinformatics Analysis the Function of RNA Binding Proteins (RBPs) and Their Prognostic Value in Breast Cancer. Front Pharmacol (2019) 10:140. doi: 10.3389/fphar.2019.00140

  • CrossRef
  • Google Scholar

• 15

  HuRZhuXChenCXuRLiYXuW. RNA-binding protein PUM2 suppresses osteosarcoma progression via partly and competitively binding to STARD13 3’UTR with miRNAs. Cell Prolif (2018) 51(6):e12508. doi: 10.1111/cpr.12508

  • CrossRef
  • Google Scholar

• 16

  HanLPWangCPHanSL. Overexpression of RBM10 induces osteosarcoma cell apoptosis and inhibits cell proliferation and migration. Med Sci (Paris) (2018) 34:Focus issue F1:81–6. doi: 10.1051/medsci/201834f114

  • CrossRef
  • Google Scholar

• 17

  EstellerM. Epigenetics in cancer. New Engl J Med (2008) 358(11):1148–59. doi: 10.1056/NEJMra072067

  • CrossRef
  • Google Scholar

• 18

  RennerMWolfTMeyerHHartmannWPenzelRUlrichAet al. Integrative DNA methylation and gene expression analysis in high-grade soft tissue sarcomas. Genome Biol (2013) 14(12):r137. doi: 10.1186/gb-2013-14-12-r137

  • CrossRef
  • Google Scholar

• 19

  SeidelCBartelFRastetterMBluemkeKWurlPTaubertHet al. Alterations of cancer-related genes in soft tissue sarcomas: hypermethylation of RASSF1A is frequently detected in leiomyosarcoma and associated with poor prognosis in sarcoma. Int J Cancer (2005) 114(3):442–7. doi: 10.1002/ijc.20707

  • CrossRef
  • Google Scholar

• 20

  GajewskiTFSchreiberHFuYX. Innate and adaptive immune cells in the tumor microenvironment. Nat Immunol (2013) 14(10):1014–22. doi: 10.1038/ni.2703

  • CrossRef
  • Google Scholar

• 21

  PetitprezFde ReynièsAKeungEZChenTWSunCMCalderaroJet al. B cells are associated with survival and immunotherapy response in sarcoma. Nature (2020) 577(7791):556–60. doi: 10.1038/s41586-019-1906-8

  • CrossRef
  • Google Scholar

• 22

  TudorGKochGG. Review of nonparametric methods for the analysis of crossover studies. Stat Methods Med Res (1994) 3(4):345–81. doi: 10.1177/096228029400300404

  • CrossRef
  • Google Scholar

• 23

  YuGWangLGHanYHeQY. clusterProfiler: an R package for comparing biological themes among gene clusters. Omics (2012) 16(5):284–7. doi: 10.1089/omi.2011.0118

  • CrossRef
  • Google Scholar

• 24

  SzklarczykDFranceschiniAWyderSForslundKHellerDHuerta-CepasJet al. STRING v10: protein-protein interaction networks, integrated over the tree of life. Nucleic Acids Res (2015) 43(Database issue):D447–52. doi: 10.1093/nar/gku1003

  • CrossRef
  • Google Scholar

• 25

  BaderGDHogueCW. An automated method for finding molecular complexes in large protein interaction networks. BMC Bioinf (2003) 4:2. doi: 10.1186/1471-2105-4-2

  • CrossRef
  • Google Scholar

• 26

  HeagertyPJZhengY. Survival model predictive accuracy and ROC curves. Biometrics (2005) 61(1):92–105. doi: 10.1111/j.0006-341X.2005.030814.x

  • CrossRef
  • Google Scholar

• 27

  GaoJAksoyBADogrusozUDresdnerGGrossBSumerSOet al. Integrative analysis of complex cancer genomics and clinical profiles using the cBioPortal. Sci Signal (2013) 6(269):pl1. doi: 10.1126/scisignal.2004088

  • CrossRef
  • Google Scholar

• 28

  XiongYWeiYGuYZhangSLyuJZhangBet al. DiseaseMeth version 2.0: a major expansion and update of the human disease methylation database. Nucleic Acids Res (2017) 45(D1):D888–d95. doi: 10.1093/nar/gkw1123

  • CrossRef
  • Google Scholar

• 29

  LvJLiuHSuJWuXLiuHLiBet al. DiseaseMeth: a human disease methylation database. Nucleic Acids Res (2012) 40(Database issue):D1030–5. doi: 10.1093/nar/gkr1169

  • CrossRef
  • Google Scholar

• 30

  KochAJeschkeJVan CriekingeWvan EngelandMDe MeyerT. MEXPRESS update 2019. Nucleic Acids Res (2019) 47(W1):W561–w5. doi: 10.1093/nar/gkz445

  • CrossRef
  • Google Scholar

• 31

  BindeaGMlecnikBTosoliniMKirilovskyAWaldnerMObenaufACet al. Spatiotemporal dynamics of intratumoral immune cells reveal the immune landscape in human cancer. Immunity (2013) 39(4):782–95. doi: 10.1016/j.immuni.2013.10.003

  • CrossRef
  • Google Scholar

• 32

  LiuXWuSYangYZhaoMZhuGHouZ. The prognostic landscape of tumor-infiltrating immune cell and immunomodulators in lung cancer. BioMed Pharmacother (2017) 95:55–61. doi: 10.1016/j.biopha.2017.08.003

  • CrossRef
  • Google Scholar

• 33

  LiTFanJWangBTraughNChenQLiuJSet al. TIMER: A Web Server for Comprehensive Analysis of Tumor-Infiltrating Immune Cells. Cancer Res (2017) 77(21):e108–e10. doi: 10.1158/0008-5472.Can-17-0307

  • CrossRef
  • Google Scholar

• 34

  LinYLiuTCuiTWangZZhangYTanPet al. RNAInter in 2020: RNA interactome repository with increased coverage and annotation. Nucleic Acids Res (2020) 48(D1):D189–D97. doi: 10.1093/nar/gkz804

  • CrossRef
  • Google Scholar

• 35

  WestphalMLamszusK. Circulating biomarkers for gliomas. Nat Rev Neurol (2015) 11(10):556–66. doi: 10.1038/nrneurol.2015.171

  • CrossRef
  • Google Scholar

• 36

  vanIDGPSzuhaiKBriaire-de BruijnIHKostineMKuijjerMLBovéeJ. Machine learning analysis of gene expression data reveals novel diagnostic and prognostic biomarkers and identifies therapeutic targets for soft tissue sarcomas. PloS Comput Biol (2019) 15(2):e1006826. doi: 10.1371/journal.pcbi.1006826

  • CrossRef
  • Google Scholar

• 37

  ZhuZJinZDengYWeiLYuanXZhangMet al. Co-expression Network Analysis Identifies Four Hub Genes Associated With Prognosis in Soft Tissue Sarcoma. Front Genet (2019) 10:37: doi: 10.3389/fgene.2019.00037

  • CrossRef
  • Google Scholar

• 38

  CooperTAWanLDreyfussG. RNA and disease. Cell (2009) 136(4):777–93. doi: 10.1016/j.cell.2009.02.011

  • CrossRef
  • Google Scholar

• 39

  KapeliKMartinezFJYeoGW. Genetic mutations in RNA-binding proteins and their roles in ALS. Hum Genet (2017) 136(9):1193–214. doi: 10.1007/s00439-017-1830-7

  • CrossRef
  • Google Scholar

• 40

  MartinJWChilton-MacNeillSKotiMvan WijnenAJSquireJAZielenskaM. Digital expression profiling identifies RUNX2, CDC5L, MDM2, RECQL4, and CDK4 as potential predictive biomarkers for neo-adjuvant chemotherapy response in paediatric osteosarcoma. PloS One (2014) 9(5):e95843. doi: 10.1371/journal.pone.0095843

  • CrossRef
  • Google Scholar

• 41

  WangLAiretiAAihaitiALiK. Expression of microRNA-150 and its Target Gene IGF2BP1 in Human Osteosarcoma and their Clinical Implications. Pathol Oncol Res (2019) 25(2):527–33. doi: 10.1007/s12253-018-0454-0

  • CrossRef
  • Google Scholar

• 42

  WangZTangWYuanJQiangBHanWPengX. Integrated Analysis of RNA-Binding Proteins in Glioma. Cancers (Basel) (2020) 12(4). doi: 10.3390/cancers12040892

  • CrossRef
  • Google Scholar

• 43

  MontesMSanfordBLComiskeyDFChandlerDS. RNA Splicing and Disease: Animal Models to Therapies. Trends Genet (2019) 35(1):68–87. doi: 10.1016/j.tig.2018.10.002

  • CrossRef
  • Google Scholar

• 44

  ChenJWeissWA. Alternative splicing in cancer: implications for biology and therapy. Oncogene (2015) 34(1):1–14. doi: 10.1038/onc.2013.570

  • CrossRef
  • Google Scholar

• 45

  YangXHuangWTHeRQMaJLinPXieZCet al. Determining the prognostic significance of alternative splicing events in soft tissue sarcoma using data from The Cancer Genome Atlas. J Transl Med (2019) 17(1):283. doi: 10.1186/s12967-019-2029-6

  • CrossRef
  • Google Scholar

• 46

  DelaunaySFryeM. RNA modifications regulating cell fate in cancer. Nat Cell Biol (2019) 21(5):552–9. doi: 10.1038/s41556-019-0319-0

  • CrossRef
  • Google Scholar

• 47

  BhatMRobichaudNHuleaLSonenbergNPelletierJTopisirovicI. Targeting the translation machinery in cancer. Nat Rev Drug Discovery (2015) 14(4):261–78. doi: 10.1038/nrd4505

  • CrossRef
  • Google Scholar

• 48

  GoudarziKMLindströmMS. Role of ribosomal protein mutations in tumor development (Review). Int J Oncol (2016) 48(4):1313–24. doi: 10.3892/ijo.2016.3387

  • CrossRef
  • Google Scholar

• 49

  GenuthNRBarnaM. The Discovery of Ribosome Heterogeneity and Its Implications for Gene Regulation and Organismal Life. Mol Cell (2018) 71(3):364–74. doi: 10.1016/j.molcel.2018.07.018

  • CrossRef
  • Google Scholar

• 50

  MiyoshiMOkajimaTMatsudaTFukudaMNNadanoD. Bystin in human cancer cells: intracellular localization and function in ribosome biogenesis. Biochem J (2007) 404(3):373–81. doi: 10.1042/bj20061597

  • CrossRef
  • Google Scholar

• 51

  WangHXiaoWZhouQChenYYangSShengJet al. Bystin-like protein is upregulated in hepatocellular carcinoma and required for nucleologenesis in cancer cell proliferation. Cell Res (2009) 19(10):1150–64. doi: 10.1038/cr.2009.99

  • CrossRef
  • Google Scholar

• 52

  AyalaGEDaiHLiRIttmannMThompsonTCRowleyDet al. Bystin in perineural invasion of prostate cancer. Prostate (2006) 66(3):266–72. doi: 10.1002/pros.20323

  • CrossRef
  • Google Scholar

• 53

  GodoyHMhawech-FaucegliaPBeckAMiliottoAMillerALeleSet al. Developmentally restricted differentiation antigens are targets for immunotherapy in epithelial ovarian carcinoma. Int J Gynecol Pathol (2013) 32(6):536–40. doi: 10.1097/PGP.0b013e318275a550

  • CrossRef
  • Google Scholar

• 54

  WhiteheadCMWinkfeinRJFritzlerMJRattnerJB. ASE-1: a novel protein of the fibrillar centres of the nucleolus and nucleolus organizer region of mitotic chromosomes. Chromosoma (1997) 106(8):493–502. doi: 10.1007/s004120050271

  • CrossRef
  • Google Scholar

• 55

  YuTXuePCuiSZhangLZhangGXiaoMet al. Rs3212986 polymorphism, a possible biomarker to predict smoking-related lung cancer, alters DNA repair capacity via regulating ERCC1 expression. Cancer Med (2018) 7(12):6317–30. doi: 10.1002/cam4.1842

  • CrossRef
  • Google Scholar

• 56

  GagliardiMStrazzulloMMatarazzoMR. DNMT3B Functions: Novel Insights From Human Disease. Front Cell Dev Biol (2018) 6:140. doi: 10.3389/fcell.2018.00140

  • CrossRef
  • Google Scholar

• 57

  GongHLTaoYMaoXZSongDYYouDNiJD. MicroRNA-29a suppresses the invasion and migration of osteosarcoma cells by regulating the SOCS1/NF-κB signalling pathway through negatively targeting DNMT3B. Int J Mol Med (2019) 44(4):1219–32. doi: 10.3892/ijmm.2019.4287

  • CrossRef
  • Google Scholar

• 58

  MegiorniFCameroSCeccarelliSMcDowellHPMannarinoOMaramponFet al. DNMT3B in vitro knocking-down is able to reverse embryonal rhabdomyosarcoma cell phenotype through inhibition of proliferation and induction of myogenic differentiation. Oncotarget (2016) 7(48):79342–56. doi: 10.18632/oncotarget.12688

  • CrossRef
  • Google Scholar

• 59

  LedererMBleyNSchleiferCHüttelmaierS. The role of the oncofetal IGF2 mRNA-binding protein 3 (IGF2BP3) in cancer. Semin Cancer Biol (2014) 29:3–12. doi: 10.1016/j.semcancer.2014.07.006

  • CrossRef
  • Google Scholar

• 60

  MancarellaCPaselloMVenturaSGrilliACalzolariLToracchioLet al. Insulin-Like Growth Factor 2 mRNA-Binding Protein 3 is a Novel Post-Transcriptional Regulator of Ewing Sarcoma Malignancy. Clin Cancer Res (2018) 24(15):3704–16. doi: 10.1158/1078-0432.Ccr-17-2602

  • CrossRef
  • Google Scholar

• 61

  AkıncılarSCKhattarEBoonPLUnalBFullwoodMJTergaonkarV. Long-Range Chromatin Interactions Drive Mutant TERT Promoter Activation. Cancer Discovery (2016) 6(11):1276–91. doi: 10.1158/2159-8290.Cd-16-0177

  • CrossRef
  • Google Scholar

• 62

  BarthelFPWeiWTangMMartinez-LedesmaEHuXAminSBet al. Systematic analysis of telomere length and somatic alterations in 31 cancer types. Nat Genet (2017) 49(3):349–57. doi: 10.1038/ng.3781

  • CrossRef
  • Google Scholar

• 63

  KillelaPJReitmanZJJiaoYBettegowdaCAgrawalNDiazLAJr.et al. TERT promoter mutations occur frequently in gliomas and a subset of tumors derived from cells with low rates of self-renewal. Proc Natl Acad Sci USA (2013) 110(15):6021–6. doi: 10.1073/pnas.1303607110

  • CrossRef
  • Google Scholar

• 64

  MichaelisKANorgardMAZhuXLevasseurPRSivagnanamSLiudahlSMet al. The TLR7/8 agonist R848 remodels tumor and host responses to promote survival in pancreatic cancer. Nat Commun (2019) 10(1):4682. doi: 10.1038/s41467-019-12657-w

  • CrossRef
  • Google Scholar

• 65

  Cherfils-ViciniJPlatonovaSGillardMLauransLValidirePCaliandroRet al. Triggering of TLR7 and TLR8 expressed by human lung cancer cells induces cell survival and chemoresistance. J Clin Invest (2010) 120(4):1285–97. doi: 10.1172/jci36551

  • CrossRef
  • Google Scholar

• 66

  RhodesDAIsenbergDA. TRIM21 and the Function of Antibodies inside Cells. Trends Immunol (2017) 38(12):916–26. doi: 10.1016/j.it.2017.07.005

  • CrossRef
  • Google Scholar

• 67

  BraunerSZhouWBacklinCGreenTMFolkersenLIvanchenkoMet al. Reduced expression of TRIM21/Ro52 predicts poor prognosis in diffuse large B-cell lymphoma patients with and without rheumatic disease. J Intern Med (2015) 278(3):323–32. doi: 10.1111/joim.12375

  • CrossRef
  • Google Scholar

• 68

  JinYZhangYLiBZhangJDongZHuXet al. TRIM21 mediates ubiquitination of Snail and modulates epithelial to mesenchymal transition in breast cancer cells. Int J Biol Macromol (2019) 124:846–53. doi: 10.1016/j.ijbiomac.2018.11.269

  • CrossRef
  • Google Scholar

• 69

  ZhaoZWangYYunDHuangQMengDLiQet al. TRIM21 overexpression promotes tumor progression by regulating cell proliferation, cell migration and cell senescence in human glioma. Am J Cancer Res (2020) 10(1):114–30.

  • Google Scholar

• 70

  ZhangPLiXHeQZhangLSongKYangXet al. TRIM21-SERPINB5 aids GMPS repression to protect nasopharyngeal carcinoma cells from radiation-induced apoptosis. J BioMed Sci (2020) 27(1):30. doi: 10.1186/s12929-020-0625-7

  • CrossRef
  • Google Scholar

• 71

  ZengQZLiuWTLuJLLiuXHZhangYFLiuLXet al. YWHAZ Binds to TRIM21 but Is Not Involved in TRIM21-stimulated Osteosarcoma Cell Proliferation. BioMed Environ Sci (2018) 31(3):186–96. doi: 10.3967/bes2018.024

  • CrossRef
  • Google Scholar

• 72

  WangZLLiBLuoYXLinQLiuSRZhangXQet al. Comprehensive Genomic Characterization of RNA-Binding Proteins across Human Cancers. Cell Rep (2018) 22(1):286–98. doi: 10.1016/j.celrep.2017.12.035

  • CrossRef
  • Google Scholar

• 73

  CallegaroDMiceliRMarianiLRautCPGronchiA. Soft tissue sarcoma nomograms and their incorporation into practice. Cancer (2017) 123(15):2802–20. doi: 10.1002/cncr.30721

  • CrossRef
  • Google Scholar

• 74

  Cancer Genome Atlas Research Network. Comprehensive and Integrated Genomic Characterization of Adult Soft Tissue Sarcomas. Cell (2017) 171(4):950–65.e28. doi: 10.1016/j.cell.2017.10.014

  • CrossRef
  • Google Scholar

• 75

  CasadeiLCaloreFCreightonCJGuesciniMBatteKIwenofuOHet al. Exosome-Derived miR-25-3p and miR-92a-3p Stimulate Liposarcoma Progression. Cancer Res (2017) 77(14):3846–56. doi: 10.1158/0008-5472.Can-16-2984

  • CrossRef
  • Google Scholar

• 76

  BottermannMFossSvan TienenLMVaysburdMCruickshankJO’ConnellKet al. TRIM21 mediates antibody inhibition of adenovirus-based gene delivery and vaccination. Proc Natl Acad Sci USA (2018) 115(41):10440–5. doi: 10.1073/pnas.1806314115

  • CrossRef
  • Google Scholar

• 77

  WangDLiuJLiuSLiW. Identification of Crucial Genes Associated With Immune Cell Infiltration in Hepatocellular Carcinoma by Weighted Gene Co-expression Network Analysis. Front Genet (2020) 11:342. doi: 10.3389/fgene.2020.00342

  • CrossRef
  • Google Scholar

• 78

  FengXFangXXiaLRenYLiXQuanXet al. Association of PPP1R13L and CD3EAP polymorphisms with risk and survival of non-small cell lung cancer in Chinese non-smoking females. Oncotarget (2017) 8(45):78749–56. doi: 10.18632/oncotarget.20224

  • CrossRef
  • Google Scholar