Skip to main content

ORIGINAL RESEARCH article

Front. Genet., 28 August 2020
Sec. RNA
This article is part of the Research Topic Perturbation of RNA Binding Protein Regulation in Cancer View all 12 articles

Identification of the Six-RNA-Binding Protein Signature for Prognosis Prediction in Bladder Cancer

\r\nYucai Wu,,,&#x;Yucai Wu1,2,3,4†Yi Liu&#x;Yi Liu5†Anbang He,,,Anbang He1,2,3,4Bao Guan,,,Bao Guan1,2,3,4Shiming He,,,Shiming He1,2,3,4Cuijian Zhang,,,Cuijian Zhang1,2,3,4Zhengjun KangZhengjun Kang5Yanqing Gong,,,*Yanqing Gong1,2,3,4*Xuesong Li,,,*Xuesong Li1,2,3,4*Liqun Zhou,,,Liqun Zhou1,2,3,4
  • 1Department of Urology, Peking University First Hospital, Beijing, China
  • 2Institute of Urology, Peking University, Beijing, China
  • 3National Urological Cancer Center, Beijing, China
  • 4Urogenital Diseases (Male) Molecular Diagnosis and Treatment Center, Peking University, Beijing, China
  • 5Department of Urology, The Fifth Affiliated Hospital of Zhengzhou University, Zhengzhou, China

RNA-binding proteins (RBPs) are a kind of gene regulatory factor that presents a significant biological effect in the initiation and development of various tumors, including bladder cancer (BLCA). However, the RBP-based prognosis signature for BLCA has not been investigated. In this study, we attempted to develop an RBP-based classifier to predict overall survival (OS) for BLCA based on transcriptome analysis. We extracted data of BLCA patients from The Cancer Genome Atlas database (TCGA) and UCSC Xena. Finally, a total of 398 cases without missing clinical data were enrolled and six RBPs (FLNA, HSPG2, AHNAK, FASTKD3, POU5F1, and PCSK9) associated with OS of BLCA were identified through univariate and multivariate Cox regression analysis. Online analyses and immunohistochemistry validated the prognostic value and expression of six RBPs. Risk scores were calculated to divide patients into high-risk and low-risk level, and patients in the high-risk group tended to have a poor prognosis. In addition, the receiver operating characteristic (ROC) curve analysis was performed to assess the prognostic value of RBPs, and the area under the curve (AUC) values were 0.711 and 0.706, respectively, in the training set and validating set. The findings were further validated in an external validation set. Subsequently, the 6-RBP-based signature and pathological stage were used to construct the nomogram to predict the 3- and 5-years OS of BLCA patients. Also, this 6-RBP-based signature was highly related to recurrence-free survival of BLCA. Weighted co-expression network analysis (WGCNA) combined with functional enrichment analysis contributed to study the potential pathways of six RBPs, including keratinocyte differentiation, RHO GTPases activate PNKs, epithelial tube morphogenesis, establishment or maintenance of cell polarity, and so on. In summary, the 6-RBP-based signature holds the potentiality to serve as a novel prognostic predictor of OS for BLCA.

Introduction

Bladder cancer (BLCA) is the 10th most prevalent cancer and the most frequently diagnosed malignancy of the urinary system all over the world (Bray et al., 2018). It has been estimated that there will be 81,440 cases of newly diagnosed BLCA and 17,980 people will die for BLCA in 2020 in the United States (Siegel et al., 2020). Non-metastatic BLCA is separated into non-muscle-invasive bladder cancer (NMIBC) and muscle-invasive bladder cancer (MIBC) and approximately 70% of BLCA patients belong to NMIBC when initially diagnosed (Dobruch et al., 2016). MIBC patients have a more favorable prognosis than those with locally advanced and metastatic BLCA due to the limited effects of surgery on advanced BLCA. In addition, BLCA is the cancer with high recurrence and about half of patients after radical surgery relapse and present with metastases (Alfred Witjes et al., 2017). However, no specific symptoms appeared in the early stage of tumor, which makes it urgent to develop novel biomarkers to predict the survival of BLCA.

RNA-binding proteins (RBPs) are a kind of key factors regulating the process of tumorigenesis, and each step that led to the initiation of malignancy may involve one or more RBPs. Mechanisms of RBPs regulation have been identified in cancer cells, including alternative splicing, polyadenylation, stability, subcellular localization, and translation (Pereira et al., 2017). Post-transcriptional regulation is an essential way of promoting or suppressing oncogenesis. RBPs can interact with other proteins and coding or non-coding RNA to form the ribonucleoprotein complexes. For example, RBPs can interact with microRNAs (miRNAs) (Liang et al., 2020), long non-coding RNAs (lncRNAs) (Jiang et al., 2020), and circular RNAs (circRNAs) (Chen et al., 2019) to affect tumor progression.

Initial assessment of BLCA has been explored in recent times. In clinical practice, lncRNAs, miRNAs, and clinicopathological factors including TNM stage and lymph node status have been gradually used to assess BLCA prognosis. Recently increasing researches have revealed that RBPs were associated with the prognosis of patients (Busa et al., 2007; Vo et al., 2012). Therefore, we aimed to identify a number of RBPs as potential biomarkers based on transcriptome analysis to predict the outcome of BLCA. We constructed a 6-RBP-based classifier for OS by using the multivariable Cox regression, which could optimize the predictivity of the current TNM staging system. Patients with gene sequencing data from the GSE13507 database were adopted as the external validation. In addition, this 6-RBP-based classifier was also highly relevant to recurrence-free survival (RFS) in BLCA. Our results demonstrated that the 6-RBP-based classifier could be used as reliable prognostic predictors of BLCA survival and recurrence.

Materials and Methods

Data Acquisition

The TCGA database was used to obtain transcriptome profiling data of tumor and normal tissues. Then, 19 normal samples and 411 BLCA samples were obtained. The matrix of mRNA expression was extracted separately by annotations using Gencode (GENCODE v 26) GTF file. Those genes whose expression was “zero” in 90% of BLCA patients were eliminated. Clinical data were downloaded from the UCSC Xena website1. To analyze the correlation of gene expression signatures with the OS of BLCA patients, we filtered out samples without clinical survival information; then, we selected a total of 398 patients and these patients were divided into training (n = 265) and validating set (n = 133) randomly at a 2:1 ratio for further analysis. Microarray study and its clinical information (GSE13507) in Gene Expression Omnibus (GEO) database2 (n = 165) were extracted, profiled by the Illumina human-6 v2.0 expression BeadChip platform. A total of 1348 genes coding for RBPs including those with high confidence for RNA binding and those annotated as RNA binding in Ensembl were summarized from the published literature (Baltz et al., 2012; Castello et al., 2012; Kwon et al., 2013; Cunningham et al., 2015).

Analysis of Differentially Expressed Genes

We used the R package edgeR to obtain differentially expressed genes (DEGs), where | log2 fold change (FC)| >1, P <0.05, and false discovery rate (FDR) <0.05 were used as the cutoffs. Then, we filtered the DEGs coding for RNA-binding proteins (DERBPs). R package “heatmap” was performed to display the selected six DERBPs.

Data Processing and Risk Score Calculation

The DERBPs were subjected to univariable Cox regression analysis to select DERBPs that were associated with OS of BLCA patients. We selected those DERBPs with P <0.001 into multivariable Cox regression to obtain the coefficients. Then, six DERBPs significantly correlated with OS were identified to build up the prediction model weighted by their coefficients. A risk score formula for OS was constructed, and each patient had been assigned a risk score according to this risk score formula that was a linear combination of the expression levels of the significant DERBPs weighted by their respective Cox regression coefficients. Then, we divided patients into low-risk group and high-risk group according to the median risk score.

Weighted Co-expression Network Analysis (WGCNA)

Considering our risk score model was built based on the expression of six RBPs, we constructed a weight co-expression network of risk score gene with DEGs in BLCA to explore the biological function by the R package “WGCNA” (Langfelder and Horvath, 2008). We selected 3 as the soft thresholding power to produce a scale-free network and the enrolled genes were hierarchically clustered into 16 modules, where the red module was found to be most relevant to risk score.

Pathway Enrichment Analysis

In order to explore the potential functions of the 6-RBP signature, genes in the red module were picked up for enrichment analysis. Pathway enrichment was conducted using an online web tool “Metascape3.” The significance threshold of FDR for significantly enriched biological processes and pathways was set at 0.05.

Statistical Analysis

We use Chi-squared test or Fisher’s exact test to measure the difference between training and validating sets and the relationship between clinical data and risk score. Both univariable and multivariable Cox regression analysis were performed using the R package “survival.” The Kaplan–Meier survival curve was drawn to demonstrate the relationship between DERBPs and OS or RFS. The log-rank test was constructed to test the significance of the difference between the two groups. ROC analysis was performed to measure prognostic accuracy. All statistical tests were two-sided, and P <0.05 was considered statistically significant. All analyses were performed in SPSS version 25.0 (SPSS Inc., Chicago, IL, United States) or R version 3.5.24 with the following packages: “edgeR” (Robinson et al., 2010), “glmnet,” “gplot,” and “survivalROC.”

Results

Data Source and Processing

As shown in Figure 1, we obtained 19 normal samples and 411 BLCA samples from TCGA database, and 4456 DEGs with | Log2FC| >1 and FDR <0.05 were identified using edgeR (Figure 2A). A total of 1348 genes coding known or predicted RBPs were matched with the 4456 DEGs and then 109 RBPs remained. After that, univariate Cox regression was performed to choose factors to predict prognostic of patients and 12 DERBPs with P <0.001 were retained for further analysis. Clinical characters of BLCA patients were downloaded from the UCSC database, and these cases were randomly divided into training set (n = 265) and validating set (n = 133) at a 2:1 ratio. There were no significant differences in age, gender, pathological stage, histologic grade, and diagnosis subtype between two sets (Table 1). Then, we identified six DERBPs, which were strongly associated with OS of BLCA by multivariate Cox regression analysis in the training set (Figure 2B), and the detailed information of these RBPs including FLNA (Filamin A), HSPG2 (Heparan Sulfate Proteoglycan 2), AHNAK (AHNAK Nucleoprotein), FASTKD3 (FAST Kinase Domains 3), POU5F1 (POU Class 5 Homeobox 1), and PCSK9 (Proprotein Convertase Subtilisin/Kexin Type 9) were listed in Table 2. Among these genes, higher expression of HSPG2, AHNAK, and PCSK9 was associated with decreased survival. On the contrary, higher expression of FLNA, FASTKD3, and POU5F1 was related to increased survival.

FIGURE 1
www.frontiersin.org

Figure 1. Study flowchart showing the process of constructing the 6-RBP-based signature to predict overall survival (OS) of bladder cancer (BLCA).

FIGURE 2
www.frontiersin.org

Figure 2. Screening for differentially expressed genes in bladder cancer (BLCA). (A) Volcano plot of differentially expressed genes in the TCGA-BLCA cohort. The red dot represented the upregulated gene, while the green dot represented the downregulated gene. (B) Forest plot exhibited the hazard ratios (HRs) with 95% confidence interval (95% CI) of prognostic RBPs in BLCA on the basis of the multivariate Cox regression result.

TABLE 1
www.frontiersin.org

Table 1. Clinical features of BLCA patients in the training and validating sets.

TABLE 2
www.frontiersin.org

Table 2. Six RBPs significantly associated with the OS of BLCA patients in the training set.

Validation the Prognostic Value and Expression of Six RBPs

To further explore the prognostic value of six RBPs in BLCA, the Kaplan–Meier plotter was used to determine the relationship between six RBPs and OS. Five of the six RBPs (AHNAK, HSPG2, PCSK9, POU5F1, and FASTKD3) were identified. Results of log-rank test demonstrated that the high expressions of AHNAK, HSPG2, and PCSK9 were associated with the low OS, while the high expression of POU5F1 and FASTKD3 was associated with the high OS of BLCA patients (Figure 3). To further validate the expression of these RBPs in BLCA, we analyzed immunohistochemistry data from the Human Protein Atlas (HPA) database5 to show that FLNA, FASTKD3, and POU5F1 were significantly decreased in BLCA compared with normal urinary bladder tissue (Figures 4A–C). Besides, the staining level of HSPG2 was increased in BLCA (Figure 4E). However, the staining level of AHNAK was relatively reduced in normal urinary bladder tissue and the result of PCSK9 protein expression was not detectable (Figure 4D). These results showed that expression of each of the six RBPs was related to prognosis of BLCA.

FIGURE 3
www.frontiersin.org

Figure 3. Kaplan–Meier analysis for overall survival (OS) of bladder cancer (BLCA) patients based on the mRNA expression of AHNAK (A), HSPG2 (B), PCSK9 (C), POU5F1 (D), and FASTKD3 (E).

FIGURE 4
www.frontiersin.org

Figure 4. Validation of the expression of RNA-binding proteins (RBPs) in bladder cancer (BLCA) and normal tissues in the Human Protein Atlas (HPA) database. (A) FLNA, (B) FASTKD3, (C) POU5F1, (D) AHNAK, and (E) HSPG2.

Development and Validation a 6-RBP-Based Classifier to Predict OS of BLCA

To assess the ability of the 6-RBP-based model predicting survival of BLCA, we created a risk score according to the expression of six RBPs as follows: Risk score = (3.50 expression value of HSPG2) + (1.74 expression value of AHNAK) + (0.35 expression value of PCSK9) - (1.49 expression value of FLNA) - (1.45 expression value of FASTKD3) - (0.28 expression value of POU5F1). Then, we calculated risk score according to this formula and cases were divided into high-risk and low-risk group based on the cutoff of median risk score (Figure 5A). The mortality was higher in the high-risk group than that in the low-risk group [HR: 2.274 (95% CI: 1.562–3.312), p < 0.001]. Moreover, HSPG2, AHNAK, and PCSK9 were highly expressed in the high-risk group, while FLNA, FASTKD3, and POU5F1 were highly expressed in the low-risk group. Results in the validating set were consistent with findings described above (Figure 5B). Kaplan–Meier curves revealed that patients in the high-risk group had shorter OS than those in the low-risk group (p < 0.001) in the training set (Figure 5C), and this result was further confirmed in the validating set (p < 0.001) (Figure 5D).

FIGURE 5
www.frontiersin.org

Figure 5. Construction of the 6-RBP-based risk model of bladder cancer (BLCA). (A) The 6-RBP-based risk score distribution, living status of BLCA patients, and heatmap of six RBP expression profiles in the training set. (B) The 6-RBP-based risk score distribution, living status of BLCA patients, and heatmap of six RBP expression profiles in the validating set. (C,D) Kaplan–Meier analysis for overall survival (OS) of BLCA patients based on the risk stratification in the training set (C) and validating set (D).

We also extracted BLCA samples from the GSE31507 database (n = 165) to validate the ability of the 6-RBP-based classifier predicting OS of BLCA. Results were compatible with that in the training and validating set derived from TCGA database. Patients in the high-risk group had shorter OS than those in the low-risk group (p = 0.014) (Figure 6A). We also found that our 6-RBP-based classifier performed well in predicting 5-year OS of BLCA (AUCs = 0.721) (Figure 6B). The mortality was higher in the high-risk group than that in the low-risk group (Figure 6C).

FIGURE 6
www.frontiersin.org

Figure 6. Further validation of the 6-RBP-based risk model in the GSE13507 dataset. (A) Kaplan–Meier analysis for overall survival (OS) of bladder cancer (BLCA) patients based on the risk stratification. (B) Receiver operating characteristic (ROC) analysis for OS prediction including 1–, 3–, and 5–years of BLCA patients. (C) The living status of BLCA patients.

Correlation Between RBPs Classifier and Clinicopathologic Characteristics

As shown in Table 3, clinicopathologic characteristics age (p = 0.032), pathological stage (p < 0.001), histologic grade (p = 0.011), and diagnosis subtype (p = 0.002) showed significant differences between the high-risk group and low-risk group in the training set. Only age (p = 0.046), pathological stage (p = 0.002), and histologic grade (p = 0.007) displayed distinct differences in the validating set. Patients with high pathological stage and histologic grade were prone to get a high-risk score in two sets (Figure 7).

TABLE 3
www.frontiersin.org

Table 3. Clinicopathological characteristics of the 6-marker-based classifier with OS in the training set and validating set.

FIGURE 7
www.frontiersin.org

Figure 7. The risk score was associated with histological grade and pathologic stage of bladder cancer (BLCA). (A,B) Boxplot of risk score in patients with different histological grade in the training set (A) and validating set (B). (C,D) Boxplot of risk score in patients with different pathologic stage in the training set (C) and validating set (D).

Prognostic Value of RBPs Classifier for Assessing Overall Survival

The 6-RBP-based classifier, age, and pathological stage were significantly related to the OS in the univariate Cox regression analysis. After the multivariate Cox regression analysis of the abovementioned factors, the 6-RBP-based classifier model and pathological stage were retained to be dependable factors for OS in the training set. Except for age, similar results were observed in the validating set (Table 4). Our result showed that the 6-RBP-based classifier was an independent prognostic factor for OS in BLCA in two sets.

TABLE 4
www.frontiersin.org

Table 4. Univariate and multivariate Cox regression analysis of the 6-marker-based classifier with OS in the training set and validating set.

In order to evaluate the capabilities of the 6-RBP-based signature to predict OS of BLCA, we plotted the ROC curves and AUC was also calculated in both cohorts. AUC values in the training set and validating set were 0.711 and 0.706, respectively, which showed that the RBP-based classifier model had an obvious better predictive accuracy compared with the TNM staging (0.670 and 0.674, respectively, in the training set and validating set) (Figures 8A,B). In consideration of the role of TNM staging in clinical practice, we combined the RBP-based model and TNM staging to predict OS of BLCA. AUC values of this joint prediction model in the training and validating set were 0.753 and 0.728, indicating that this model was more accurate than models enrolled in RBPs or TNM staging solely. Subsequently, the 6-RBP-based risk score and pathological stage were used to construct the nomogram to predict the 3- and 5-year OS of BLCA patients (Figure 8C).

FIGURE 8
www.frontiersin.org

Figure 8. Estimate the predictive ability of RBP-based classifier to predict survival in bladder cancer (BLCA) patients. (A,B) The receiver operating characteristic (ROC) for TNM stage only, the RBP-based risk score (BLCARBP6), and the RBP-based risk score combined with TNM stage for overall survival (OS) in the training set (A) and validating set (B). (C) Nomogram to predict the 3- and 5-years OS. (D) Kaplan–Meier analysis for recurrence-free survival (RFS) of BLCA patients based on the risk stratification.

Prognostic Value of the RBP-Based Classifier for Assessing RFS

We further explored whether this 6-RBP-based classifier was related to RFS of BLCA. As shown in Table 5, univariate and multivariate Cox regression analysis was conducted to identify prognostic factors for RFS. Finally, outcomes of univariate and multivariate analysis indicated that pathological stage and the 6-RBP-based classifier were independent risk factors of RFS in BLCA patients. Survival analysis revealed that the RFS of patients in the high-risk group were considerably shorter than that in the low-risk group (Figure 8D). Our results demonstrated that this 6-RBP-based classifier could also be used as a reliable prognostic predictor of BLCA recurrence.

TABLE 5
www.frontiersin.org

Table 5. Univariate and multivariate Cox regression analysis of the 6-marker-based classifier with RFS in TCGA database.

Pathway Enrichment Analysis of DERBPs

To explore the biological function of the 6-RBP signature in BLCA, the WGCNA method was performed to choose genes associated with risk score. Sixteen modules were identified by hierarchical clustering analysis, and the red module was significantly associated risk score (Figure 9A). In addition, genes in the red module served to perform pathway enrichment analysis. Results of enrichment analysis showed that genes were mostly enriched in keratinocyte differentiation, and RHO GTPases activate PNKs, epithelial tube morphogenesis, establishment or maintenance of cell polarity, and so on, suggesting that these pathways were correlated with the disease progression of BLCA with high-risk score (Figures 9B,C).

FIGURE 9
www.frontiersin.org

Figure 9. Pathway enrichment analysis of the RBPs signature. (A) Weighted co-expression network analysis (WGCNA) plot was performed to cluster genes associated with the 6-RBP-based risk score, and the data of the risk score were added to construct the eigengene adjacency heatmap. (B) Pathways associated with the 6-RBP-based signature were enriched by Metascape. (C) The histogram of the top 10 enriched pathways associated with the 6-RBP-based risk score. The abscissa was the value of -Log10P, and the longitudinal axis denotes different enrichment pathways, sorted by the value of -Log10P.

Discussion

RBPs were regarded acting as amplifiers of oncogenic driver mutations. Increasing literature demonstrated that RBPs were of vital importance in the initiation, development, and recurrence of multi-malignancies. For instance, LIN28B overexpression promoted tumorigenesis and metastasis of colon cancer via repressing the level of let-7 microRNAs (King et al., 2011). RBPs also played a critical role in the initiation of BLCA. It was reported that HuR was upregulated in the BLCA tissue compared with adjacent normal tissue and it could promote BLCA progression by competitively binding to the long non-coding HOX transcript antisense RNA with miR-1 (Yu et al., 2017). RBPs also exhibit potentiality as novel biomarkers. Pancreatic ductal adenocarcinoma patients with higher levels of ESRP1 showed longer survival than those with low ESRP1 expression (Ueda et al., 2014). In conclusion, RBPs can regulate the biology of cancer and hold potential as novel biomarkers.

Diverse models to predict the outcome of BLCA have been created, including miRNA-based signature (Hofbauer et al., 2018; Yin et al., 2019), clinical character-based nomogram (Zhang et al., 2019), and lncRNA-based model (Wang et al., 2020). Parts of them performed well in predicting the survival of BLCA. However, no RBP-based classifier for predicting survival in BLCA has been established. RBPs are a subset molecule of regulating progression and development of malignancies. Considering the limited capabilities of single RBP for prediction prognosis, we constructed a predicted model based on mRNA expression of six RBPs by univariate and multivariate Cox regression analysis. Patients were divided into two categories based on the median-risk score and high-risk patients have shorter survival than low-risk patients, suggesting that our model displayed powerful strength to forecast the OS for BLCA. Clinicopathologic characteristics including age, pathological stage, and histologic grade were associated with the expression of six RBPs and patients with high stage or grade tended to get high-risk score. Our 6-RBP-based classifier with AUCs being 0.711 and 0.706 in the training and validating sets presented a strong ability to predict OS of BLCA. In addition, TNM stage is the most commonly used index to assess the treatment way and outcome of BLCA. We also combined the 6-RBP-based signature together with TNM stage to assess the prognosis and the AUCs showed that this model was more accurate than the 6-RBP-based model. Furthermore, we explored the efficiency of the 6-RBP-based classifier in predicting RFS of BLCA. Results showed that pathological stage and the 6-RBP-based classifier were independent factors of RFS in BLCA patients and patients in the low-risk group showed a significantly longer RFS than those in the high-risk group.

Six prognosis-related RBPs including FLNA, HSPG2, AHNAK, FASTKD3, POU5F1, and PCSK9 were selected to build a classifier. It has been confirmed that FLNA was downregulated in BLCA tissue and overexpression of FLNA repressed migration, invasion, and migration of BLCA by regulating autophagy (Wang et al., 2018). AHNAK downregulated in BLCA performed accurately on discriminating between benign urothelial lesion and bladder urothelial carcinoma using voided-urine liquid-based cytology (Lee et al., 2018). In addition, low expression of POU5F1 was associated with shorter cancer-related survival and might be a novel biomarker for BLCA (Chang et al., 2008). However, functions of HSPG2, FASTKD3, and PCSK9 in BLCA have not been explored yet. These confirmed or predicted prognosis-related RBPs support evidence that our model is capable of assessing the outcome of BLCA.

In order to explore the biological function of the 6-RBP signature, we performed WGCNA and pathway enrichment analysis. Results showed that those genes relevant to risk score were mainly enriched in keratinocyte differentiation, and RHO GTPases activate PNKs, epithelial tube morphogenesis, establishment or maintenance of cell polarity, etc. Interestingly, we compared pathways predicted with annotations of these six RBPs in GeneCards6 and found that HSPG2 acted as an anti-angiogenic and anti-tumor peptide that inhibited endothelial cell migration and collagen-induced endothelial tube morphogenesis. Angiogenesis is thought to be a critical procedure of promoting BLCA progression and associated with poor survival (Bochner et al., 1995; Roudnicky et al., 2017). Among these prognostic RBPs, HSPG2 demonstrated an anti-angiogenesis effect via binding to α2β1 integrin and interacting with VEGFR2 at the surface of endothelial cells (Woodall et al., 2008; Goyal et al., 2012; Poluzzi et al., 2016). A previous study has demonstrated that the breakdown of cell polarity programs could promote the occurrence of aggressive, invasive tumors (VanderVorst et al., 2018). In BLCA, BMP4 could induce monocyte/macrophage polarization toward M2 phenotype macrophages, which promoted the progression of BLCA (Martinez et al., 2017). FLNA also regulated cell polarity by interacting with FilGAP, a Rac-specific GTPase-activating protein (Nakamura et al., 2009), but this mechanism had not been elucidated in BLCA. Exact mechanisms of these RBPs remain largely unknown and more research is required to investigate their roles in BLCA.

Conclusion

In conclusion, we identified six RBPs associated with prognosis of BLCA and constructed the 6-RBP-based classifier to help clinical decision, while optimizing the predictive ability of the current TNM staging system. This study takes the initiative report that the RBP-based classifier could predict the prognosis in human BLCA. Nevertheless, large-scale, multi-center, and prospective studies are necessary to confirm our results before the 6-RBP-based signature can be applied in the clinic.

Data Availability Statement

All datasets generated for this study are included in the article/supplementary material.

Ethics Statement

This study is a secondary data analysis, and the raw data come from public databases. The institutional and/or national research ethics committee has approved the data collection and management process.

Author Contributions

YW and YL: design, analysis, and interpretation of data, drafting of the manuscript, and critical revision of the manuscript. AH and BG: acquisition of data and statistical analysis. SH, CZ, and ZK: research direction. YG, XL, and LZ: critical revision of the manuscript for important intellectual content, administrative support, obtaining funding, and supervision. All authors contributed to the article and approved the submitted version.

Funding

This study was funded by the National Natural Science Foundation of China (Nos. 81772703 and 81972380) and the Natural Science Foundation of Beijing (7172219).

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.

Abbreviations

AUCs, Area under the curves; BLCA, Bladder cancer; circRNAs, Circular RNAs; DEGs, Differently expressed genes; DERBPs, DEGs coding for RNA-binding proteins; FD, Fold change; FDR, False discovery rate; GEO, Gene Expression Omnibus; HPA, Human Protein Atlas; HR, Hazard ratio; lncRNAs, Long non-coding RNAs; MIBC, Muscle-invasive bladder cancer; miRNAs, MicroRNAs; NMIBC, Non-muscle-invasive bladder cancer; OS, Overall survival; RBPs, RNA-binding proteins; RFS, Recurrence-free survival; ROC, Receiver operating characteristic; TCGA, The Cancer Genome Atlas; WGCNA, Weighted co-expression network analysis.

Footnotes

  1. ^ https://xena.ucsc.edu/
  2. ^ http://www.ncbi.nlm.nih.gov/geo/
  3. ^ http://metascape.org/
  4. ^ http://www.r-project.org/
  5. ^ https://www.proteinatlas.org/
  6. ^ https://www.genecards.org/

References

Alfred Witjes, J., Lebret, T., Comperat, E. M., Cowan, N. C., De Santis, M., Bruins, H. M., et al. (2017). Updated 2016 EAU Guidelines on muscle-invasive and metastatic bladder Cancer. Eur. Urol. 71, 462–475. doi: 10.1016/j.eururo.2016.06.020

PubMed Abstract | CrossRef Full Text | Google Scholar

Baltz, A. G., Munschauer, M., Schwanhausser, B., Vasile, A., Murakawa, Y., Schueler, M., et al. (2012). The mRNA-bound proteome and its global occupancy profile on protein-coding transcripts. Mol. Cell 46, 674–690. doi: 10.1016/j.molcel.2012.05.021

PubMed Abstract | CrossRef Full Text | Google Scholar

Bochner, B. H., Cote, R. J., Weidner, N., Groshen, S., Chen, S. C., Skinner, D. G., et al. (1995). Angiogenesis in bladder cancer: relationship between microvessel density and tumor prognosis. J. Natl. Cancer Inst. 87, 1603–1612. doi: 10.1093/jnci/87.21.1603

PubMed Abstract | CrossRef Full Text | Google Scholar

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 Cancer J. Clin. 68, 394–424. doi: 10.3322/caac.21492

PubMed Abstract | CrossRef Full Text | Google Scholar

Busa, R., Paronetto, M. P., Farini, D., Pierantozzi, E., Botti, F., Angelini, D. F., et al. (2007). The RNA-binding protein Sam68 contributes to proliferation and survival of human prostate cancer cells. Oncogene 26, 4372–4382. doi: 10.1038/sj.onc.1210224

PubMed Abstract | CrossRef Full Text | Google Scholar

Castello, A., Fischer, B., Eichelbaum, K., Horos, R., Beckmann, B. M., Strein, C., et al. (2012). Insights into RNA biology from an atlas of mammalian mRNA-binding proteins. Cell 149, 1393–1406. doi: 10.1016/j.cell.2012.04.031

PubMed Abstract | CrossRef Full Text | Google Scholar

Chang, C. C., Shieh, G. S., Wu, P., Lin, C. C., Shiau, A. L., and Wu, C. L. (2008). Oct-3/4 expression reflects tumor progression and regulates motility of bladder cancer cells. Cancer Res. 68, 6281–6291. doi: 10.1158/0008-5472.CAN-08-0094

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, Y., Yang, F., Fang, E., Xiao, W., Mei, H., Li, H., et al. (2019). Circular RNA circAGO2 drives cancer progression through facilitating HuR-repressed functions of AGO2-miRNA complexes. Cell Death. Differ. 26, 1346–1364. doi: 10.1038/s41418-018-0220-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Cunningham, F., Amode, M. R., Barrell, D., Beal, K., Billis, K., Brent, S., et al. (2015). Ensembl 2015. Nucleic Acids Res. 43, D662–D669. doi: 10.1093/nar/gku1010

PubMed Abstract | CrossRef Full Text | Google Scholar

Dobruch, J., Daneshmand, S., Fisch, M., Lotan, Y., Noon, A. P., Resnick, M. J., et al. (2016). Gender and bladder Cancer: a collaborative review of etiology. biology, and outcomes. Eur. Urol. 69, 300–310. doi: 10.1016/j.eururo.2015.08.037

PubMed Abstract | CrossRef Full Text | Google Scholar

Goyal, A., Poluzzi, C., Willis, C. D., Smythies, J., Shellard, A., Neill, T., et al. (2012). Endorepellin affects angiogenesis by antagonizing diverse vascular endothelial growth factor receptor 2 (VEGFR2)-evoked signaling pathways: transcriptional repression of hypoxia-inducible factor 1alpha and VEGFA and concurrent inhibition of nuclear factor of activated T cell 1 (NFAT1) activation. J. Biol. Chem. 287, 43543–43556. doi: 10.1074/jbc.M112.401786

PubMed Abstract | CrossRef Full Text | Google Scholar

Hofbauer, S. L., de Martino, M., Lucca, I., Haitel, A., Susani, M., Shariat, S. F., et al. (2018). A urinary microRNA (miR) signature for diagnosis of bladder cancer. Urol Oncol.36, 531.e1–531.e8. doi: 10.1016/j.urolonc.2018.09.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Jiang, D., Zhang, Y., Yang, L., Lu, W., Mai, L., Guo, H., et al. (2020). Long noncoding RNA HCG22 suppresses proliferation and metastasis of bladder cancer cells by regulation of PTBP1. J. Cell. Physiol. 235, 1711–1722. doi: 10.1002/jcp.29090

PubMed Abstract | CrossRef Full Text | Google Scholar

King, C. E., Cuatrecasas, M., Castells, A., Sepulveda, A. R., Lee, J. S., and Rustgi, A. K. (2011). LIN28B promotes colon cancer progression and metastasis. Cancer Res. 71, 4260–4268. doi: 10.1158/0008-5472.CAN-10-4637

PubMed Abstract | CrossRef Full Text | Google Scholar

Kwon, S. C., Yi, H., Eichelbaum, K., Fohr, S., Fischer, B., You, K. T., et al. (2013). The RNA-binding protein repertoire of embryonic stem cells. Nat. Struct. Mol. Biol. 20, 1122–1130. doi: 10.1038/nsmb.2638

PubMed Abstract | CrossRef Full Text | Google Scholar

Langfelder, P., and Horvath, S. (2008). WGCNA: an R package for weighted correlation network analysis. BMC Bioinformatics 9:559. doi: 10.1186/1471-2105-9-559

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, H., Kim, K., Woo, J., Park, J., Kim, H., Lee, K. E., et al. (2018). Quantitative proteomic analysis identifies AHNAK (Neuroblast Differentiation-associated Protein AHNAK) as a novel candidate biomarker for bladder urothelial carcinoma diagnosis by liquid-based cytology. Mol. Cell. Proteomics 17, 1788–1802. doi: 10.1074/mcp.RA118.000562

PubMed Abstract | CrossRef Full Text | Google Scholar

Liang, G., Meng, W., Huang, X., Zhu, W., Yin, C., Wang, C., et al. (2020). miR-196b-5p-mediated downregulation of TSPAN12 and GATA6 promotes tumor progression in non-small cell lung cancer. Proc. Natl. Acad. Sci. U.S.A. 117, 4347–4357. doi: 10.1073/pnas.1917531117

PubMed Abstract | CrossRef Full Text | Google Scholar

Martinez, V. G., Rubio, C., Martinez-Fernandez, M., Segovia, C., Lopez-Calderon, F., Garin, M. I., et al. (2017). BMP4 Induces M2 macrophage polarization and favors tumor progression in Bladder Cancer. Clin. Cancer Res 23, 7388–7399. doi: 10.1158/1078-0432.CCR-17-1004

PubMed Abstract | CrossRef Full Text | Google Scholar

Nakamura, F., Heikkinen, O., Pentikainen, O. T., Osborn, T. M., Kasza, K. E., Weitz, D. A., et al. (2009). Molecular basis of filamin A-FilGAP interaction and its impairment in congenital disorders associated with filamin A mutations. PLoS One 4:e4928. doi: 10.1371/journal.pone.0004928

PubMed Abstract | CrossRef Full Text | Google Scholar

Pereira, B., Billaud, M., and Almeida, R. (2017). RNA-binding proteins in cancer: old players and new actors. Trends Cancer 3, 506–528. doi: 10.1016/j.trecan.2017.05.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Poluzzi, C., Iozzo, R. V., and Schaefer, L. (2016). Endostatin and endorepellin: a common route of action for similar angiostatic cancer avengers. Adv. Drug Deliv. Rev. 97, 156–173. doi: 10.1016/j.addr.2015.10.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Robinson, M. D., McCarthy, D. J., and Smyth, G. K. (2010). edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 26, 139–140. doi: 10.1093/bioinformatics/btp616

PubMed Abstract | CrossRef Full Text | Google Scholar

Roudnicky, F., Dieterich, L. C., Poyet, C., Buser, L., Wild, P., Tang, D., et al. (2017). High expression of insulin receptor on tumour-associated blood vessels in invasive bladder cancer predicts poor overall and progression-free survival. J. Pathol. 242, 193–205. doi: 10.1002/path.4892

PubMed Abstract | CrossRef Full Text | Google Scholar

Siegel, R. L., Miller, K. D., and Jemal, A. (2020). Cancer statistics, 2020. CA J. Cancer Clin. 70, 7–30. doi: 10.3322/caac.21590

PubMed Abstract | CrossRef Full Text | Google Scholar

Ueda, J., Matsuda, Y., Yamahatsu, K., Uchida, E., Naito, Z., Korc, M., et al. (2014). Epithelial splicing regulatory protein 1 is a favorable prognostic factor in pancreatic cancer that attenuates pancreatic metastases. Oncogene 33, 4485–4495. doi: 10.1038/onc.2013.392

PubMed Abstract | CrossRef Full Text | Google Scholar

VanderVorst, K., Hatakeyama, J., Berg, A., Lee, H., and Carraway, K. L. III (2018). Cellular and molecular mechanisms underlying planar cell polarity pathway contributions to cancer malignancy. Semin. Cell Dev. Biol. 81, 78–87. doi: 10.1016/j.semcdb.2017.09.026

PubMed Abstract | CrossRef Full Text | Google Scholar

Vo, D. T., Subramaniam, D., Remke, M., Burton, T. L., Uren, P. J., Gelfond, J. A., et al. (2012). The RNA-binding protein Musashi1 affects medulloblastoma growth via a network of cancer-related genes and is an indicator of poor prognosis. Am. J. Pathol. 181, 1762–1772. doi: 10.1016/j.ajpath.2012.07.031

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Y., Du, L., Yang, X., Li, J., Li, P., Zhao, Y., et al. (2020). A nomogram combining long non-coding RNA expression profiles and clinical factors predicts survival in patients with bladder cancer. Aging 12, 2857–2879. doi: 10.18632/aging.102782

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Z., Li, C., Jiang, M., Chen, J., Yang, M., and Pu, J. (2018). Filamin A (FLNA) regulates autophagy of bladder carcinoma cell and affects its proliferation, invasion and metastasis. Int. Urol. Nephrol. 50, 263–273. doi: 10.1007/s11255-017-1772-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Woodall, B. P., Nystrom, A., Iozzo, R. A., Eble, J. A., Niland, S., Krieg, T., et al. (2008). Integrin alpha2beta1 is the required receptor for endorepellin angiostatic activity. J. Biol. Chem. 283, 2335–2343. doi: 10.1074/jbc.M708364200

PubMed Abstract | CrossRef Full Text | Google Scholar

Yin, X. H., Jin, Y. H., Cao, Y., Wong, Y., Weng, H., Sun, C., et al. (2019). Development of a 21-miRNA signature associated with the prognosis of patients with bladder Cancer. Front. Oncol. 9:729. doi: 10.3389/fonc.2019.00729

PubMed Abstract | CrossRef Full Text | Google Scholar

Yu, D., Zhang, C., and Gui, J. (2017). RNA-binding protein HuR promotes bladder cancer progression by competitively binding to the long noncoding HOTAIR with miR-1. Onco Targets Ther. 10, 2609–2619. doi: 10.2147/OTT.S132728

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Y., Hong, Y. K., Zhuang, D. W., He, X. J., and Lin, M. E. (2019). Bladder cancer survival nomogram: development and validation of a prediction tool, using the SEER and TCGA databases. Medicine 98:e17725. doi: 10.1097/MD.0000000000017725

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: bladder cancer, RNA-binding proteins, overall survival, recurrence-free survival, prognosis

Citation: Wu Y, Liu Y, He A, Guan B, He S, Zhang C, Kang Z, Gong Y, Li X and Zhou L (2020) Identification of the Six-RNA-Binding Protein Signature for Prognosis Prediction in Bladder Cancer. Front. Genet. 11:992. doi: 10.3389/fgene.2020.00992

Received: 11 May 2020; Accepted: 05 August 2020;
Published: 28 August 2020.

Edited by:

Yongsheng Kevin Li, Harbin Medical University, China

Reviewed by:

Yunyan Gu, Harbin Medical University, China
Zhaoying Yang, Jilin University, China

Copyright © 2020 Wu, Liu, He, Guan, He, Zhang, Kang, Gong, Li and Zhou. 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: Yanqing Gong, yqgong@bjmu.edu.cn; Xuesong Li, pineneedle@sina.com

These authors have contributed equally to this work

Disclaimer: 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.