- Department of Pathology, School of Basic Medical Sciences, Fudan University, Shanghai, China
The dysregulation of RNA binding proteins (RBPs) is closely related to tumorigenesis and development. However, the role of RBPs in Colon adenocarcinoma (COAD) is still poorly understood. We downloaded COAD’s RNASeq data from the Cancer Genome Atlas (TCGA) database, screened the differently expressed RBPs in normal tissues and tumor, and constructed a protein interaction network. COAD patients were randomly divided into a training set (N = 315) and a testing set (N = 132). In the training set, univariate Cox analysis identified 12 RBPs significantly related to the prognosis of COAD. By multivariate COX analysis, we constructed a prognostic model composed of five RBPs (CELF4, LRRFIP2, NOP14, PPARGC1A, ZNF385A) based on the lowest Akaike information criterion. Each COAD patient was scored according to the model formula. Further analysis showed that compared with the low-risk group, the overall survival rate (OS) of patients in the high-risk group was significantly lower. The area under the curve of the time-dependent receiver operator characteristic (ROC) curve was 0.722 in the training group and 0.738 in the test group, which confirmed a good prediction feature. In addition, a nomogram was constructed based on clinicopathological characteristics and risk scores. C-index and calibration curve proved the accuracy in predicting the 1-, 3-, and 5-year survival rates of COAD patients. In short, we constructed a superior prognostic and diagnostic signature composed of five RBPs, which indicates new possibilities for individualized treatment of COAD patients.
Introduction
Colorectal cancer is the third most common cancer and the third most common cause of cancer-related death (1). By 2030, global new cases of colorectal cancer infection are expected to exceed 2.2 million, and the death toll will reach 1.1 million (2). Colon adenocarcinoma (COAD) is the most common type of colorectal cancer (3) and primarily occurs in the intestinal mucosa. COAD usually grows into the intestinal lumen and spreads to adjacent organs. It is a highly aggressive malignant tumor with a high mortality and recurrence rate (4, 5). Although the clinical treatments for colon cancer, including surgical techniques, radiotherapy and chemotherapy, have been improved, the prognosis of patients is still poor (6, 7). Genes were affected by a variety of regulatory mechanisms, including but not being limited to DNA methylation, histone deacetylation and miRNA expression, thereby promoting the occurrence and development of tumors (8, 9). Therefore, discovery of new regulatory factors and therapeutic targets for COAD is imperative to improve our understanding of cancer occurrence and disease progression.
Post-transcriptional gene regulation (PTGR) is necessary in order to maintain cell metabolism and coordinate the maturation, transportation, stabilization and degradation of various types of RNA. RNA-binding proteins (RBPs) play a key role in the processes of PTGR (10). They interact with target gene mRNAs and manipulate the processing of these mRNAs to determine cell behaviors (11). In fact, each of these events is regulated by the formation of ribonucleoprotein (RNP) complex with different core RBPs (12). Recent studies have revealed that RBPs are associated with neurodegenerative diseases, muscle atrophy, diabetes, as well as different cancers and developmental disorders (13–16). However, the role of RBPs in cancer continues to be poorly understood.
RBPs play a major role in the development of colon cancer. RBM3, a member of RBPs, is up-regulated in a phase-dependent manner, and its overexpression can induce oncogenic transformation (17). Increased expression of CELF1 leads to growth arrest of intestinal epithelial cells in G1 phase, while its silence promotes cell proliferation (18). Overexpression of LIN28B is associated with an aggressive phenotype, worsened survival rate and increased recurrence of tumor (19, 20). However, these studies are still far from enough to explore the panorama of RBPs in COAD. Therefore, we obtained the expression data and patient data of COAD patients from TCGA, explored the functions of RBPs through a series of bioinformatics and statistical analysis, and constructed a prognostic signature composed of five RBPs to predict the prognosis of COAD patients. In addition, we constructed a nomogram, hoping to be able to provide valuable insights for the individualized treatment of COAD patients.
Materials and Methods
Data Processing
We obtained RNA sequencing data of 41 normal colon tissue samples and 473 tumor samples and corresponding data of 447 patients from the official website of the Cancer Genome Atlas Database (https://portal.gdc.cancer.gov/). M stage was removed in multivariate analysis because of too much missing information. Moreover, we obtained 1542 RBPs from these articles (10, 21). We used limma (22) package to analyze the difference of RBPs, | logFC | ≥ 0.5, FDR pvalue < 0.05 was the cut-off value.
KEGG Pathway and GO Enrichment Analysis
We used clusterProfiler (23) package to do GO enrichment and Kyoto Genome Encyclopedia (KEGG) analysis method to comprehensively analyze the biological functions of differentially expressed RBPs. GO analysis terms include cell component (CC), molecular function (MF) and biological process (BP). FDR P value<0.05 as a filter condition.
PPI Network Construction
Differentially expressed RBPs were submitted to the string database to determine the information on protein-protein interactions (https://string-db.org/). Using Cytoscape V3.8.0 software was used for visualization and the most relevant sub network modules were obtained by using the molecular complex detection (MCODE) (24) plug-in. P ≤0.05 represents significant difference.
Construction and Verification of Prediction Model
Univariate Cox analysis was performed in the training set to screen the differences of RBPs related to overall survival (OS) in COAD patients, and then based on multivariate Cox regression analysis of the lowest Akaike information criterion (AIC) value, a proportional risk regression model of 5 RBPs prognosis signature were obtained. The risk score was calculated based on the expression of these 5 genes. The risk score formula was as follows: represent the regression coefficient and the expression value of each prognosis related RBPs, respectively. According to the median value, COAD patients were divided into high and low risk groups. The Kaplan-Meier survival curve was used to compare survival differences. Principal component analysis (PCA) was used to visualize gene expression patterns. The receiver-operating characteristic (ROC) curves were performed evaluate the prognostic accuracy of the model. In addition, the testing set was applied to confirm the prediction ability of the prediction model.
Establishment and Validation of Nomogram
In the training set, we constructed a nomogram to predict the survival of COAD patients in 1, 3, 5 years by combining the clinicopathological characteristics of age, gender, stage, T stage and N stage, as well as the risk score obtained from prognostic signature. In addition, the concordance index (C-index) was used to evaluate the descriptive and predictive capabilities of nomograms. The calibration curve of nomogram was performed to test whether the predicted survival rate was consistent with the actual survival rate.
Verification of Express Level and Prognostic Significance
The Wilcoxon test was used to detect the expression of 5 RBPs at the transcriptional level in COAD patients. Kaplan Meier survival curve proves the prognostic valued of RBPs.
Statistical Analysis
In this study, Strawberry Perl for windows (Version5.18.2) was used for data processing, and R (3.6.2) was performed for data analysis. P < 0.05 was considered statistically significant.
Result
The flow of this research was presented in Figure 1. A total of 447 COAD patients from the TCGA-COAD cohort were enrolled and were randomly divided into training set (N = 315) and testing set (N = 132). Table 1 summarized the detailed clinical characteristics of these patients.
Identification of Different RBPs
We analyzed RNA-seq data from 473 COAD samples and 41 normal colon tissues from the TCGA database and obtained 1,542 RBPs gene expression data. We followed (FDR P value < 0.05, | logFC | > 0.5) as the screening criteria to obtain 321 up-regulated genes and 151 down-regulated genes (Figures 2A, B).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/627504/fonc-11-627504-HTML/image_m/fonc-11-627504-g002.jpg)
Figure 2 The differentially expressed RBPs in COAD. (A) Heat map. Columns represent tissue samples, and rows represent differentially expressed RBP. N, normal; T, tumor. (B) Volcano plot. The selection criteria were |logFC| > 0.5 and FDR p value < 0.05.
Functional Enrichment Analysis of Differential RBPs
In order to explore the biological processes and functions of these differentially expressed genes, we performed gene ontology (GO) and Kyoto encyclopedia of genes and genomes (KEGG) analysis on 321 up-regulated genes and 151 down-regulated genes, respectively. The result showed that up-regulated RBPs were significantly enriched in the biological process (BP) ribosome-related processes, including ribosome biogenesis, ribonucleoprotein complex biogenesis, and rRNA metabolic process. Down-regulated RBPs were mainly enriched in RNA cleavage, including RNA fragmentation, mRNA metabolic process, and RNA splicing. In terms of cell components (CC), the up-regulated differentially expressed RBPs were mainly ribosomal components, such as 90S preribosome, small-subunit processome, and ribonucleoprotein granule. Down-regulation of RBPs mainly existed in spliceosomal complex and catalytic step 2 spliceosome. However, molecular function analysis showed that up-regulated RBPs increased significantly nuclear activity, ribonuclease activity and endonuclease activity, while down-regulated RBPs significantly decreased translation reporter activity, ribonuclease activity and mRNA 3 ‘- UTR binding (Figures 3A, B). Furthermore, KEGG analysis showed similar results. Up-regulated RBPs were significantly enriched in ribosome and RNA degradation (Figure 3C), while down-regulated PBRs were enriched in RNA transport and spliceosome (Figure 3D).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/627504/fonc-11-627504-HTML/image_m/fonc-11-627504-g003.jpg)
Figure 3 KEGG pathway and GO enrichment analysis of aberrantly expressed RBPs. (A) Enrichment of up-regulated RBPs in GO. (B) Enrichment of down-regulated RBPs in GO. (C) Enrichment of up-regulated RBPs in KEGG. (D) Enrichment of down-regulated RBPs in KEGG.
Protein-Protein Interaction Network Construction
We further studied the role of differentially expressed RBPs in COAD. Through the STRING database, we obtained the PPI network with required interaction score greater than 0.9, which includes 464 nodes and 2,288 edges. Cytoscape software was used for visualization (Figure 4A). The MCODE tool was applied to identify possible key modules. Figures 4B and D showed the three most important modules. The first important module was shown in Figure 4B. Ribosome biogenesis, rRNA metabolic process, and rRNA processing were significantly enriched in module 1 with 63 RBPs (Table S1).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/627504/fonc-11-627504-HTML/image_m/fonc-11-627504-g004.jpg)
Figure 4 Protein-protein interaction network and modules analysis. (A) Protein-protein interaction network of differentially expressed RBPs. (B) Module 1 in PPI network. (C) Module 2 in PPI network. (D) Module 3 in PPI network. Blue circles: down-regulation; red circles: up-regulation.
Construction and Evaluation of Prognostic Signature
Then we constructed a prognostic model in the training set. Univariate Cox regression analysis showed that 12 RBPs (NOP14, POP1, EIF2AK3, PPARGC1A, ZNF385A, RP9, NSUN5, G3BP2, PNLDC1, CELF4, LRRFIP2, CAPRIN2) gene expression levels were significantly correlated with the overall survival rate (Figure 5A; P < 0.05). Hazard ratios (HRs) identify risk-related genes (HR > 1) and protective genes (HR < 1). Subsequently, multivariate Cox regression analysis was performed on the candidate RBPs-related genes to evaluate their roles as independent prognostic factors for patient survival. Based on the lowest Akaike information criterion (AIC) value, we identified 5 RBPs (CELF4, LRRFIP2, NOP14, PPARGC1A, ZNF385A) as prognostic signature (Figure 5B).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/627504/fonc-11-627504-HTML/image_m/fonc-11-627504-g005.jpg)
Figure 5 Establishment of RBPs-related prognosis model in training set. (A) Univariate Cox regression analysis for identification of prognosis related RBPs. (B) Multivariate Cox regression analysis for constructing model.
The risk score of each patient was determined according to the prognostic signature model, and COAD patients were divided into high and low-risk groups in the light of the median value of risk value. Kaplan Meier survival curve showed that the survival rate of COAD patients in high-risk group was significantly lower than that in low-risk group (Figure 6A). Receiver operating characteristic (ROC) curve was performed to evaluate the prognosis of the model (Figure 6B). The area under the curve (AUC) was 0.722 which indicated good prediction ability. Principal component analysis showed that patients with high and low risk had different distribution patterns (Figure 6C). We ranked patients from low to high according to risk value, and the scatter plot demonstrated the survival rate of high-risk patients was much lower than that of low-risk patients (Figures 6D, E). The heat map displayed the difference of 5 RBPs genes expressions in high and low risk groups (Figure 6F).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/627504/fonc-11-627504-HTML/image_m/fonc-11-627504-g006.jpg)
Figure 6 Validation of the RBPs-related prognostic signature in training set. (A) Kaplan-Meier survival curve showed that the survival time of the high-risk group was significantly shorter than that of the low-risk score group. (B) ROC curve analysis of risk score as an independent prognostic factor. (C) PCA showed that five RBPs are distributed differently between high-risk groups and low-risk groups. (D) As the risk score increases, the distribution of high and low risk COAD patients. (E) The scatter plot displayed the correlation between the survival status of high- and low-risk COAD patients. (F) The heat map showed the differences in the expression of five RBPs between high and low risk patients.
Moreover, testing set was applied for external verification. As shown in Figures 7A–C, the survival rate of patients in the high-risk group was lower than that in the low-risk group. The area under ROC was 0.738, and PCA analysis showed different distribution. These data proved that RBPs-related genes can accurately predict the survival of COAD patients. Similarly, we analyzed the risk characteristics of the training set (Figures 7D–F). As we assumed, the risk signatures of the five RBPs-related genes were robust.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/627504/fonc-11-627504-HTML/image_m/fonc-11-627504-g007.jpg)
Figure 7 Evaluation of the RBPs-related prognostic signature in the testing set. (A) Survival curve for low- and high-risk subgroups. (B) ROC curves for forecasting OS based on risk score. (C) Distribution pattern of PCA in different risk subgroups. (D) Risk score distribution of COAD patients. (E) Survival status of COAD patients. (F) Expression heat map.
RBPs-Related Gene Prognostic Signature Is an Independent Prognostic Factor
Next, we performed univariate and multivariate Cox regression analysis based on clinicopathological characteristics (age, gender, stage, T stage, N stage) and risk scores of COAD patients, to determine whether RBPs-related gene prognostic characteristics are independent predictors for COAD patients. Univariate analysis showed that age (P < 0.018), stage (P < 0.001), T stage (P < 0.001), N stage (P < 0.001) and risk score (P < 0.001) were significantly correlated with OS in both training and testing sets (Figures 8A, C). Multivariate analysis showed that only age and risk score had significant correlation with OS in COAD patients (P < 0.05; Figures 8B, D). These data suggest that RBPs-related gene prognostic signatures are independent factors affecting the prognosis of COAD patients.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/627504/fonc-11-627504-HTML/image_m/fonc-11-627504-g008.jpg)
Figure 8 Risk score based on RBPs-related prognostic signature was an independent predictor for COAD patients. (A) Univariate Cox regression analysis in training set. (B) Multivariate Cox regression analysis in training set. (C) Univariate Cox regression analysis in testing set. (D) Multivariate Cox regression analysis in testing set.
Construction of a Nomogram
The nomogram integrates multiple prognostic factors to evaluate the survival probability of an individual at a specific time point and displays it graphically (25). To further predict the survival of COAD patients, we constructed a nomogram consisting of clinicopathological features (age, gender, stage, T stage, N stage) and risk score (Figure 9A). Nomography predicted the 1-, 3-, 5-year survival rate of COAD patients. The calibration curve showed that the actual patient survival was consistent with the predicted value (Figures 9B–D). The concordance index (C index) of nomogram was 0.770, which proves the accurate prediction performance of the nomogram. These results indicate that the nomogram with risk score can accurately predict the 1-, 3-, 5-year survival rate of patients, and provide valuable insights for individualized treatment of COAD patients.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/627504/fonc-11-627504-HTML/image_m/fonc-11-627504-g009.jpg)
Figure 9 Construction and verification of nomogram. (A) The prognostic nomogram constructed based on the risk score of RBPs and clinicopathological parameters predicted the survival rate of COAD patients at 1, 3, and 5 years. (B–D) Calibration curves showed the concordance between predicted and observed 1-, 3-, and 5-year survival rates.
Validation the Prognostic Value and Expression of RBPs
We further explored the prognostic value and expression value of five RBPs in COAD patients. The Kaplan-Meier survival curve showed that the expressions of CELF4, NOP14, PPARGC1A, and ZNF385A were correlated with the survival of COAD patients (all P < 0.05; Figures 10A–D), while the remaining LRRFIP2 was not significantly correlated with survival (P = 0.072; Figure 10E). We then analyzed whether there were differences in the expressions of these genes between normal and cancer tissues. Surprisingly, NOP14 was up-regulated in tumors and the other four (CELF4, PPARGC1A, ZFF385A, LRRFIP2) were down-regulated in tumors (Figure 11).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/627504/fonc-11-627504-HTML/image_m/fonc-11-627504-g010.jpg)
Figure 10 Kaplan Meier survival curve analysis verified the prognostic value of RBPs in COAD. (A) CELF4; (B) ZNF385A; (C) NOP14; (D) PPARGC1A; (E) LRRFIP2.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/627504/fonc-11-627504-HTML/image_m/fonc-11-627504-g011.jpg)
Figure 11 Verification of the transcription expression of RBPs in COAD and normal colon tissues. (A) CELF4; (B) LRRFIP2; (C) NOP14; (D) PPARGC1A; (E) ZNF385A.
Conclusion and Discussion
RNA-binding proteins (RBPs) play an important role in cellular homeostasis by controlling gene expression at the post-transcriptional level (26). RBPs are dysregulated in several cancer types, affecting the expressions and functions of both tumor proteins and tumor-suppressor proteins (27). In this study, we obtained the original RNASeq data of COAD patients from TCGA and identified 472 differentially expressed RBPs genes. Then we performed enrichment analysis on up-regulated genes and down-regulated genes. The results showed that the up-regulated RBPs were significantly enriched in the ribosome-related processes, including ribosome biosynthesis, ribonucleoprotein complex biosynthesis, and rRNA metabolism. Down-regulated RBPs were mainly enriched in RNA cleavage, including RNA fragmentation, mRNA metabolic process, and RNA splicing. Relevant studies have proved that the regulation of RNA processing and RNA metabolism is related to the development of COAD (28, 29). In addition, we constructed a PPI network to comprehensively display the correlation between differentially expressed RBPs proteins in COAD.
RBPs are a group of genes that regulate the growth, development and survival of cancer cells, and are closely related to cancer progression (30, 31). For instance, IGF2BP3 may contributes to lung tumorigenesis by regulating the alternative splicing of PKM (32), HuR promotes the progression of head and neck squamous cell carcinoma and bladder cancer (33, 34). Therefore, RBPs are potential biomarkers that most likely predict cancer risk and survival outcomes. In this study, we systematically analyzed the prediction accuracy of the prognostic signature of RBPs in COAD using bioinformatics and statistical tools. We randomly divided COAD patients into training set (N = 315) and testing set (N = 132) at 7:3. The multivariate COX regression analysis identified 5 genes as the optimal items for constructing a prognostic model based on the lowest AIC value. We verified the accuracy of RBPs in predicting the prognosis of COAD patients both in the training set and testing set. Furthermore, multivariate COX analysis confirmed that the risk score based on the prognostic signature of RBPs was an independent predictor of COAD patients (P < 0.05). Nomography is widely used to predict the survival rate of cancer patients (35). We constructed and evaluated a nomogram composed of clinicopathological characteristics (age, gender, stage, T stage, N stage) and risk scores. These results suggest that the present RBPs prognostic signature accurately predicts the survival of COAD patients.
Recently, a Cox model based on seven RBPs was reported to predict the survival of patients with COAD, two of which (PPARGC1A and LRRFIP2) were also present in our model (36). Among the five RBPs in our model, the other three genes (CELF4, NOP14, ZNF385A) were not previously reported to be associated with COAD. We used the Kaplan-Meier survival curve to explore the prognostic value of these five genes. CELF4, LRRFIP2, NOP14 and ZNF385A were associated with OS in COAD patients (P < 0.05), while PPARGC1A was not significantly correlated with the overall survival of COAD patients (P = 0.072). However, since we applied the average gene expression value instead of the best P value as the cutoff criteria to divide the high and low expression groups, we cannot completely deny the prognostic value of PPARGC1A.
RBPs regulate the expressions of genes required for many aspects of cancer behaviors including its sensitivity to chemotherapy in post-transcriptional network (37). In the identified five RBPs, the expressions of CELF4, LRRFIP2, ZNF385A, and PPARGC1A were down-regulated, while the expression of NOP14 was up-regulated (all P < 0.001). As we know, the main treatment of colon cancer is surgical resection, supplemented by chemotherapy and radiotherapy. Chemoresistance in colorectal cancer is urged to be conquered. CELF4 (CUGBP, ELAV-like family member 4) is one of six mammalian CELF proteins that function in mRNA metabolism. Current studies have shown that CELF4 genetic variation contributes to chemotherapy-related cardiac dysfunction (38, 39). Whether CELF4 contributes to chemoresistance in COAD patients is not reported yet. LRRFIP2 (Leucine-rich repeat flightless-interacting protein 2) is a signal regulator that interacts with Toll-like receptor (TLR) adaptor protein MyD88 to regulate NF-κB activity (40). NF-κB activates the expression of c-MYC, ICAM-1, and VEGFA, and promotes the growth and proliferation of colon cancer cells (41). This suggests that targeting LRRFIP2 may be effective for COAD treatment. ZNF385A (Zinc finger protein 385A) acts as a transcription factor which modulates the activation of PAK-2p34 by proteasome-mediated degradation. ZNF385A is related to cognitive decline in one’s later years (42), however, its relationship with cancer is rarely mentioned. Since we found that the high level of ZNF385A indicates the poor prognosis of COAD patients, it may be a potential target for further COAD research. NOP14 (NOP14 Nucleolar Protein) plays a role in pre-18s rRNA processing and small ribosomal subunit assembly. It is associated with a variety of cancer progression including pancreatic cancer (43), melanoma (44) and breast cancer (45). Although NOP14 lacks detailed research in the field of COAD, we speculate that it has the ability to promote tumor progression based on the literature. Thus, we screened out five RBPs which have a huge excavation potential and provide probable targets for COAD treatment.
In conclusion, we constructed a prognostic signature of RBPs which can accurately predict the survival outcome of COAD patients. We combined the prognostic signature and other clinicopathological characteristics to establish and verify a prognostic nomogram. Our data indicates that the identified five RBPs are the potential prognostic and diagnostic biomarkers, which provide a valuable reference for practitioners’ clinical decision and enable individualization of COAD treatment.
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.
Author Contributions
CY performed the statistical analyses and drafted the manuscript. KC supervised the statistical analyses, participated in data analysis and interpretation, and provided critical feedback. KC and CY provided data for fine-mapping. All authors contributed to the article and approved the submitted version.
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.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2021.627504/full#supplementary-material
References
1. Siegel RL, Miller KD, Jemal A. Cancer statistics, 2018. CA Cancer J Clin (2018) 68(1):7–30. doi: 10.3322/caac.21442
2. Arnold M, Sierra MS, Laversanne M, Soerjomataram I, Jemal A, Bray F. Global patterns and trends in colorectal cancer incidence and mortality. Gut (2017) 66(4):683–91. doi: 10.1136/gutjnl-2015-310912
3. Mutch MG. Molecular profiling and risk stratification of adenocarcinoma of the colon. J Surg Oncol (2007) 96(8):693–703. doi: 10.1002/jso.20915
4. McDermott FT, Hughes ES, Pihl E, Milne BJ, Price AB. Prognosis in relation to symptom duration in colon cancer. Br J Surg (1981) 68(12):846–9. doi: 10.1002/bjs.1800681206
5. Labianca R, Nordlinger B, Beretta GD, Mosconi S, Mandala M, Cervantes A, et al. Early colon cancer: ESMO Clinical Practice Guidelines for diagnosis, treatment and follow-up. Ann Oncol (2013) 24(Suppl 6):vi64–72. doi: 10.1093/annonc/mdt354
6. Bresalier RS, Kopetz S, Brenner DE. Blood-based tests for colorectal cancer screening: do they threaten the survival of the FIT test? Dig Dis Sci (2015) 60(3):664–71. doi: 10.1007/s10620-015-3575-2
7. Siegel R, DeSantis C, Virgo K, Stein K, Mariotto A, Smith T, et al. Cancer treatment and survivorship statistics, 2012. CA Cancer J Clin (2012) 62(4):220–41. doi: 10.3322/caac.21149
8. Perri F, Longo F, Giuliano M, Sabbatino F, Favia G, Ionna F, et al. Epigenetic control of gene expression: Potential implications for cancer treatment. Crit Rev Oncol Hematol (2017) 111:166–72. doi: 10.1016/j.critrevonc.2017.01.020
9. Della Vittoria Scarpati G, Calura E, Di Marino M, Romualdi C, Beltrame L, Malapelle U, et al. Analysis of differential miRNA expression in primary tumor and stroma of colorectal cancer patients. BioMed Res Int (2014) 2014:840921. doi: 10.1155/2014/840921
10. Gerstberger S, Hafner M, Tuschl T. A census of human RNA-binding proteins. Nat Rev Genet (2014) 15(12):829–45. doi: 10.1038/nrg3813
11. Mittal N, Scherrer T, Gerber AP, Janga SC. Interplay between posttranscriptional and posttranslational interactions of RNA-binding proteins. J Mol Biol (2011) 409(3):466–79. doi: 10.1016/j.jmb.2011.03.064
12. Masuda K, Kuwano Y. Diverse roles of RNA-binding proteins in cancer traits and their implications in gastrointestinal cancers. Wiley Interdiscip Rev RNA (2019) 10(3):e1520. doi: 10.1002/wrna.1520
13. De Conti L, Baralle M, Buratti E. Neurodegeneration and RNA-binding proteins. Wiley Interdiscip Rev RNA (2017) 8(2). doi: 10.1002/wrna.1394
14. Yao P, Fox PL. Aminoacyl-tRNA synthetases in medicine and disease. EMBO Mol Med (2013) 5(3):332–43. doi: 10.1002/emmm.201100626
15. Nutter CA, Kuyumcu-Martinez MN. Emerging roles of RNA-binding proteins in diabetes and their therapeutic potential in diabetic complications. Wiley Interdiscip Rev RNA (2018) 9(2). doi: 10.1002/wrna.1459
16. Gerstberger S, Hafner M, Ascano M, Tuschl T. Evolutionary conservation and expression of human RNA-binding proteins and their role in human genetic disease. Adv Exp Med Biol (2014) 825:1–55. doi: 10.1007/978-1-4939-1221-6_1
17. Sureban SM, Ramalingam S, Natarajan G, May R, Subramaniam D, Bishnupuri KS, et al. Translation regulatory factor RBM3 is a proto-oncogene that prevents mitotic catastrophe. Oncogene (2008) 27(33):4544–56. doi: 10.1038/onc.2008.97
18. Liu L, Ouyang M, Rao JN, Zou T, Xiao L, Chung HK, et al. Competition between RNA-binding proteins CELF1 and HuR modulates MYC translation and intestinal epithelium renewal. Mol Biol Cell (2015) 26(10):1797–810. doi: 10.1091/mbc.E14-11-1500
19. Madison BB, Liu Q, Zhong X, Hahn CM, Lin N, Emmett MJ, et al. LIN28B promotes growth and tumorigenesis of the intestinal epithelium via Let-7. Genes Dev (2013) 27(20):2233–45. doi: 10.1101/gad.224659.113
20. Tu HC, Schwitalla S, Qian Z, LaPier GS, Yermalovich A, Ku YC, et al. LIN28 cooperates with WNT signaling to drive invasive intestinal and colorectal adenocarcinoma in mice and humans. Genes Dev (2015) 29(10):1074–86. doi: 10.1101/gad.256693.114
21. Li W, Gao LN, Song PP, You CG. Development and validation of a RNA binding protein-associated prognostic model for lung adenocarcinoma. Aging (Albany NY) (2020) 12(4):3558–73. doi: 10.18632/aging.102828
22. Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, et al. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res (2015) 43(7):e47. doi: 10.1093/nar/gkv007
23. Yu G, Wang LG, Han Y, He QY. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS (2012) 16(5):284–7. doi: 10.1089/omi.2011.0118
24. Bader GD, Hogue CW. An automated method for finding molecular complexes in large protein interaction networks. BMC Bioinf (2003) 4:2. doi: 10.1186/1471-2105-4-2
25. Di Trapani E, Sanchez-Salas R, Gandaglia G, Rocchini L, Moschini M, Lizee D, et al. A nomogram predicting the cancer-specific mortality in patients eligible for radical cystectomy evaluating clinical data and neoadjuvant cisplatinum-based chemotherapy. World J Urol (2016) 34(2):207–13. doi: 10.1007/s00345-015-1640-2
26. Mittal N, Roy N, Babu MM, Janga SC. Dissecting the expression dynamics of RNA-binding proteins in posttranscriptional regulatory networks. Proc Natl Acad Sci USA (2009) 106(48):20300–5. doi: 10.1073/pnas.0906940106
27. Pereira B, Billaud M, Almeida R. 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
28. Chen L, Luo C, Shen L, Liu Y, Wang Q, Zhang C, et al. SRSF1 Prevents DNA Damage and Promotes Tumorigenesis through Regulation of DBF4B Pre-mRNA Splicing. Cell Rep (2017) 21(12):3406–13. doi: 10.1016/j.celrep.2017.11.091
29. Zhou X, Li X, Cheng Y, Wu W, Xie Z, Xi Q, et al. BCLAF1 and its splicing regulator SRSF10 regulate the tumorigenic potential of colon cancer cells. Nat Commun (2014) 5:4581. doi: 10.1038/ncomms5581
30. Wu Y, Liu Y, He A, Guan B, He S, Zhang C, et al. Identification of the Six-RNA-Binding Protein Signature for Prognosis Prediction in Bladder Cancer. Front Genet (2020) 11:992. doi: 10.3389/fgene.2020.00992
31. Kim MY, Hur J, Jeong S. Emerging roles of RNA and RNA-binding protein network in cancer cells. BMB Rep (2009) 42(3):125–30. doi: 10.5483/bmbrep.2009.42.3.125
32. Xueqing H, Jun Z, Yueqiang J, Xin L, Liya H, Yuanyuan F, et al. IGF2BP3 May Contributes to Lung Tumorigenesis by Regulating the Alternative Splicing of PKM. Front Bioeng Biotechnol (2020) 8:679. doi: 10.3389/fbioe.2020.00679
33. Xu CZ, Jiang C, Wu Q, Liu L, Yan X, Shi R. A Feed-Forward Regulatory Loop between HuR and the Long Noncoding RNA HOTAIR Promotes Head and Neck Squamous Cell Carcinoma Progression and Metastasis. Cell Physiol Biochem (2016) 40(5):1039–51. doi: 10.1159/000453160
34. Yu D, Zhang C, Gui J. RNA-binding protein HuR promotes bladder cancer progression by competitively binding to the long noncoding HOTAIR with miR-1. Onco Targets Ther (2017) 10:2609–19. doi: 10.2147/OTT.S132728
35. Shariat SF, Margulis V, Lotan Y, Montorsi F, Karakiewicz PI. Nomograms for bladder cancer. Eur Urol (2008) 54(1):41–53. doi: 10.1016/j.eururo.2008.01.004
36. Wang Y, Chen Y, Xiao S, Fu K. Integrated Analysis of the Functions and Prognostic Values of RNA-Binding Proteins in Colorectal Cancer. Front Cell Dev Biol (2020) 8:595605. doi: 10.3389/fcell.2020.595605
37. Blagden S, Abdel Mouti M, Chettle J. Ancient and modern: hints of a core post-transcriptional network driving chemotherapy resistance in ovarian cancer. Wiley Interdiscip Rev RNA (2018) 9(1). doi: 10.1002/wrna.1432
38. Linschoten M, Teske AJ, Cramer MJ, van der Wall E, Asselbergs FW. Chemotherapy-Related Cardiac Dysfunction: A Systematic Review of Genetic Variants Modulating Individual Risk. Circ Genom Precis Med (2018) 11(1):e001753. doi: 10.1161/CIRCGEN.117.001753
39. Wang X, Sun CL, Quinones-Lombrana A, Singh P, Landier W, Hageman L, et al. CELF4 Variant and Anthracycline-Related Cardiomyopathy: A Children’s Oncology Group Genome-Wide Association Study. J Clin Oncol (2016) 34(8):863–70. doi: 10.1200/JCO.2015.63.4550
40. Gunawardena HP, Huang Y, Kenjale R, Wang H, Xie L, Chen X. Unambiguous characterization of site-specific phosphorylation of leucine-rich repeat Fli-I-interacting protein 2 (LRRFIP2) in Toll-like receptor 4 (TLR4)-mediated signaling. J Biol Chem (2011) 286(13):10897–910. doi: 10.1074/jbc.M110.168179
41. Li X, Wei R, Wang M, Ma L, Zhang Z, Chen L, et al. MGP Promotes Colon Cancer Proliferation by Activating the NF-kappaB Pathway through Upregulation of the Calcium Signaling Pathway. Mol Ther Oncolytics (2020) 17:371–83. doi: 10.1016/j.omto.2020.04.005
42. Yu L, Dawe RJ, Boyle PA, Gaiteri C, Yang J, Buchman AS, et al. Association Between Brain Gene Expression, DNA Methylation, and Alteration of Ex Vivo Magnetic Resonance Imaging Transverse Relaxation in Late-Life Cognitive Decline. JAMA Neurol (2017) 74(12):1473–80. doi: 10.1001/jamaneurol.2017.2807
43. Du YX, Li ZZ, Liu CZ, Liu ZW. [Expression of nucleolar protein 14 and CD31 in pancreatic cancer mouse model and its correlation with tumor progression]. Zhonghua Wai Ke Za Zhi (2017) 55(6):463–7. doi: 10.3760/cma.j.issn.0529-5815.2017.06.013
44. Li J, Zhao R, Fang R, Wang J. [miR-122-5p inhibits the proliferation of melanoma cells by targeting NOP14]. Nan Fang Yi Ke Da Xue Xue Bao (2018) 38(11):1360–5. doi: 10.12122/j.issn.1673-4254.2018.11.14
Keywords: colon adenocarcinoma (COAD), RNA binding proteins (RBPs), prognostic, signature, nomograph
Citation: Chang K, Yuan C and Liu X (2021) A New RBPs-Related Signature Predicts the Prognosis of Colon Adenocarcinoma Patients. Front. Oncol. 11:627504. doi: 10.3389/fonc.2021.627504
Received: 09 November 2020; Accepted: 27 January 2021;
Published: 09 March 2021.
Edited by:
Alessandro Ottaiano, Istituto Nazionale Tumori Fondazione G. Pascale (IRCCS), ItalyReviewed by:
Angela Lombardi, University of Campania Luigi Vanvitelli, ItalyFrancesco Perri, Istituto Nazionale Tumori Fondazione G. Pascale (IRCCS), Italy
Copyright © 2021 Chang, Yuan and Liu. 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: Xueguang Liu, glxg69@shmu.edu.cn
†These authors have contributed equally to this work