Development and validation of novel prognostic models based on RNA-binding proteins in breast cancer.
Ontology highlight
ABSTRACT: Objectives
We aimed to construct novel prognostic models based on RNA-binding proteins (RBPs) in breast cancer (BRCA) and explore their roles in this disease and their effects on tumor-infiltrating immune cells (TIICs).Methods
Datasets were downloaded from the Gene Expression Omnibus (GEO) database. Functions and prognostic values of RBPs were systematically investigated using a series of bioinformatics analysis methods. TIICs were assessed using CIBERSORT.Results
Overall, 138 differentially expressed RBPs were identified, of which 86 were upregulated and 52 were downregulated. Of these, 13 RBPs were identified as prognosis-related and adopted to construct an overall survival (OS) model, while 12 RBPs were used for the relapse-free survival (RFS) model. High-risk patients had poorer OS and RFS rates than low-risk patients. The results indicate that the OS and RFS models are good prognostic models with reliable predictive abilities. In addition, the proportions of CD8, CD4 naïve, and CD4 memory resting T cells, as well as resting dendritic cells, were significantly different between the low-risk and high-risk groups in the OS model.Conclusions
OS and RFS signatures can be used as reliable BRCA prognostic biomarkers. This work will help understand the prognostic roles and functions of RBPs in BRCA.
SUBMITTER: Fan W
PROVIDER: S-EPMC9252011 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:BackgroundRNA-binding proteins (RBPs) play a crucial role in regulating RNA turnover and are associated with cancer development. However, little is known about the role of RBPs in esophageal cancer (ESCA). The present study focuses on the association between RBP gene expression and survival in ESCA, addressing the clinical relevance of an RBPs-based prediction model for prognosis.MethodsRNA-sequencing data and clinical information of patients with ESCA were obtained from The Cancer Genome Atlas (TCGA) database. We identified differentially expressed genes in ESCA and intersected them with RBP-encoding genes. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were performed with the identified differentially expressed RBPs. Then, a protein-protein interaction (PPI) network was constructed through the STRING database to determine the hub RBPs. Univariate Cox regression analysis and multivariate Cox regression analysis were applied to construct a novel prognostic model based on RBPs. Based on the R package "Caret", we divided patients into the training set and validation set. The efficacy of the prognostic model was evaluated by the area under the receiver operating characteristic (ROC) curve. A nomogram was developed for the prediction of patient survival outcomes.ResultsA total of 158 ESCA patients from the TCGA database were included in our analysis. We screened out five prognostic RBPs (CLK1, CIRBP, MRPL13, TNRC6A, and TYW3) through univariate and multivariate Cox regression analysis. CLK1, CIRBP, TNRC6A and TYW3 were downregulated in tumor samples, while MRPL13 was upregulated. A prognostic model constructed with these five RBPs in the training data set accurately stratified ESCA patients into high- and low-risk groups. When the same prognostic model was applied to the test data set and entire cohort, the 5-RBP signature remained an independent prognostic factor in multivariate analysis. The areas under the time-dependent ROC curve of the prognostic model for predicting one-year survival in the training data set, test data set, and entire cohort were 0.789, 0.753, and 0.764, respectively, confirming that this model is a good prognostic model. The nomogram based on the five RBPs and clinical variables could improve individualized outcome predictions and highlight the importance of RBPs in the outcomes of patients with ESCA.ConclusionsOur study provides a potential prognostic model for predicting the prognosis of ESCA patients. The prognostic nomogram could improve individualized outcome predictions for patients with ESCA, therefore providing novel insights into future diagnosis and treatment.
| S-EPMC10713301 | biostudies-literature
Project description:Ribonucleic acid-binding proteins (RBPs) are reportedly involved in tumor progression and recurrence; however, the functions and mechanisms of action of RBPs in ovarian serous cystadenocarcinoma (OSC) are not known. To address these issues, gene expression profiles of OSC tissues from The Cancer Genome Atlas (TCGA) and normal tissues from the Genotype-Tissue Expression database were compared in order to identify RBPs that are differentially expressed in OSC. We also analyzed the biological functions of these RBPs and their relationship to clinical outcome. There were 190 RBPs that were differentially expressed between OSC and normal tissues, including 93 that were upregulated and 97 that were downregulated. Five of the RBPs were used to construct a prediction model that was evaluated by univariate and multivariate Cox regression analyses. TCGA data were randomly divided into training and test cohorts, and further categorized into high- and low-risk groups according to risk score in the model. The overall survival (OS) of the high-risk group was shorter than that of the low-risk group (training cohort P = 0.0007596; test cohort P = 0.002219). The area under the receiver operating characteristic curve of the training and test cohorts was 0.701 and 0.638, respectively, demonstrating that the model had good predictive power. A nomogram was established to quantitatively describe the relationship between the five prognostic RBPs and OS in OSC, which can be useful for developing individualized management strategies for patients.
| S-EPMC7593419 | biostudies-literature
Project description:As the most common transcriptional regulators, zinc finer proteins (ZNFs) play vital roles in occurrence and progression of malignant tumors. Whereas, information regarding the roles of ZNFs in soft tissue sarcomas (STS) remains scarce. In this study, a comprehensive bioinformatics analysis investigating roles of ZNFs in STS was performed. Initially, we extracted raw datasets of differentially expressed ZNFs from GSE2719. Using a sequence of bioinformatics methods, we then investigated the prognostic significance, function, and molecular subtype of these differentially expressed ZNFs. In addition, CCK8 and plate clone formation assays were used to explore the effect of ZNF141 on STS cells. A total of 110 differentially expressed ZNFs were identified. Nine ZNFs (HLTF, ZNF292, ZNF141, LDB3, PHF14, ZNF322, PDLIM1, NR3C2, and LIMS2) were selected to establish an overall survival (OS) prediction model, and seven ZNFs (ZIC1, ZNF141, ZHX2, ZNF281, ZNHIT2, NR3C2, and LIMS2) were used to develop a progression-free survival (PFS) prediction model. Compared with patients with low-risk in the TCGA training and testing cohorts, as well as the GEO validation cohorts, patients with high-risk had poorer OS and PFS. Using nomograms constructed with the identified ZNFs predicting OS and PFS, we established a clinically useful model. Four distinct molecular subtypes with different prognostic and immune infiltration characteristics were identified. In vitro experiments showed that ZNF141 promoted the proliferation and viability of STS cells. In conclusion, ZNF-related models are useful as prognostic biomarkers, suggesting their potentials as therapeutic targets in STS. These findings will enable us to develop novel strategies treating STS, which will potentially improve outcomes of patients with STS.
| S-EPMC10188339 | biostudies-literature
Project description:Background:Recent progress in sequencing technologies allows us to explore comprehensive genomic and transcriptomic information to improve the current European LeukemiaNet (ELN) system of acute myeloid leukemia (AML). Methods:We compared the prognostic value of traditional demographic and cytogenetic risk factors, genomic data in the form of somatic aberrations of 25 AML-relevant genes, and whole-transcriptome expression profiling (RNA sequencing) in 267 intensively treated AML patients (Clinseq-AML). Multivariable penalized Cox models (overall survival [OS]) were developed for each data modality (clinical, genomic, transcriptomic), together with an associated prognostic risk score. Results:Of the three data modalities, transcriptomic data provided the best prognostic value, with an integrated area under the curve (iAUC) of a time-dependent receiver operating characteristic (ROC) curve of 0.73. We developed a prognostic risk score (Clinseq-G) from transcriptomic data, which was validated in the independent The Cancer Genome Atlas AML cohort (RNA sequencing, n = 142, iAUC = 0.73, comparing the high-risk group with the low-risk group, hazard ratio [HR]OS = 2.42, 95% confidence interval [CI] = 1.51 to 3.88). Comparison between Clinseq-G and ELN score iAUC estimates indicated strong evidence in favor of the Clinseq-G model (Bayes factor = 26.78). The proposed model remained statistically significant in multivariable analysis including the ELN and other well-known risk factors (HRos = 2.34, 95% CI = 1.30 to 4.22). We further validated the Clinseq-G model in a second independent data set (n = 458, iAUC = 0.66, adjusted HROS = 2.02, 95% CI = 1.33 to 3.08; adjusted HREFS = 2.10, 95% CI = 1.42 to 3.12). Conclusions:Our results indicate that the Clinseq-G prediction model, based on transcriptomic data from RNA sequencing, outperforms traditional clinical parameters and previously reported models based on genomic biomarkers.
| S-EPMC6186516 | biostudies-literature
Project description:Enhancer RNAs (eRNAs) are present specifically in tumors, where they affect the expression of eRNA-regulated genes (ERGs). Owing to this characteristic, ERGs were hypothesized to improve prognosis of overall survival in heterogeneous low-grade and intermediate-grade gliomas. This study aimed to construct and validate an ERG prognostic tool to facilitate clinical management, and offer more effective diagnostic and therapeutic biomarkers for glioma. Survival-related eRNAs were identified, and their ERGs were selected based on eRNA and target gene information. The ERG prognostic model was constructed and validated using internal and external validation cohorts. Finally, biological differences related to the ERG signature were analysed to explore the potential mechanisms influencing survival outcomes. Thirteen ERGs were identified and used to build an ERG risk signature, which included five super-enhancer RNA (seRNA)-regulated genes and five LGG-specific eRNA-regulated genes. The prognostic nomogram established based on combining the ERG score, age, and sex was evaluated by calibration curves, clinical utility, Harrell's concordance index (0.86; 95% CI: 0.83-0.90), and time-dependent receiver operator characteristic curves. We also explored potential immune-related mechanisms that might cause variation in survival. The established prognostic model displayed high validity and robustness. Several immune-related genes regulated by seRNAs or specific eRNAs were identified, indicating that these transcripts or their genes were potential targets for improving immunotherapeutic/therapeutic outcomes. The functions of an important specific eRNA-regulated gene (USP28) were validated in robust vitro experiments. In addition, the ERG risk signature was significantly associated with the immune microenvironment and other immune-related features.
| S-EPMC8921558 | biostudies-literature
Project description:BackgroundMultiple myeloma (MM) is a malignant hematopoietic disease that is usually incurable. RNA-binding proteins (RBPs) are involved in the development of many tumors, but their prognostic significance has not been systematically described in MM. Here, we developed a prognostic signature based on eight RBP-related genes to distinguish MM cohorts with different prognoses.MethodAfter screening the differentially expressed RBPs, univariate Cox regression was performed to evaluate the prognostic relevance of each gene using The Cancer Genome Atlas (TCGA)-Multiple Myeloma Research Foundation (MMRF) dataset. Lasso and stepwise Cox regressions were used to establish a risk prediction model through the training set, and they were validated in three Gene Expression Omnibus (GEO) datasets. We developed a signature based on eight RBP-related genes, which could classify MM patients into high- and low-score groups. The predictive ability was evaluated using bioinformatics methods. Gene ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment, and gene set enrichment analyses were performed to identify potentially significant biological processes (BPs) in MM.ResultThe prognostic signature performed well in the TCGA-MMRF dataset. The signature includes eight hub genes: HNRNPC, RPLP2, SNRPB, EXOSC8, RARS2, MRPS31, ZC3H6, and DROSHA. Kaplan-Meier survival curves showed that the prognosis of the risk status showed significant differences. A nomogram was constructed with age; B2M, LDH, and ALB levels; and risk status as prognostic parameters. Receiver operating characteristic (ROC) curve, C-index, calibration analysis, and decision curve analysis (DCA) showed that the risk module and nomogram performed well in 1, 3, 5, and 7-year overall survival (OS). Functional analysis suggested that the spliceosome pathway may be a major pathway by which RBPs are involved in myeloma development. Moreover, our signature can improve on the R-International Staging System (ISS)/ISS scoring system (especially for stage II), which may have guiding significance for the future.ConclusionWe constructed and verified the 8-RBP signature, which can effectively predict the prognosis of myeloma patients, and suggested that RBPs are promising biomarkers for MM.
| S-EPMC8107400 | biostudies-literature
Project description:PurposeThe function of stanniocalcin-1 (STC-1) in the oncogenesis and progression of tumors has been extensively studied. The purpose of this study was to investigate the relationship between secreted STC-1 and prognosis in patients with breast cancer (BC) and to determine whether STC-1 could be a key prognostic factor in BC.MethodsThe STC-1 level was measured by ELISA and clinical data from 1210 female patients with BC were used to develop and validate nomograms. We then verified the models through the plotting of ROC curves and calibration curves, calculating the C-index, and performing decision curve analyses (DCA).ResultsThe level of STC-1 in the peripheral plasma was significantly correlated with the T stage, N stage, clinical stage, grade, hormone receptors, HER-2 status, and tumor subtype. Cox regression analyses revealed that estrogen receptor(ER) status, N stage, and STC-1 level were risk factors for overall survival (OS), whereas T stage, N stage, and STC-1 level were independent prognostic factors for distant disease-free survival (DDFS) and disease-free survival (DFS). Both the ROC curve and the C-index confirmed the high resolution of these models, while the DCA identified the feasibility of their practical application. In addition, the calibration curves indicated good consistency between the predicted and actual survival rates.ConclusionNomograms were created based on STC-1 levels for 3-, 5-, and 7-year OS, DDFS, and DFS of patients with BC respectively. As a key prognostic factor for BC, peripheral blood STC-1 level can be used clinically as a liquid biopsy indicator.
| S-EPMC10067061 | biostudies-literature
Project description:Immunotherapy has shown excellent therapeutic effects on various malignant tumors; however, to date, immunotherapy for osteosarcoma is still suboptimal. In this study, we performed comprehensive bioinformatic analysis of immune-related genes (IRGs) and tumor-infiltrating immune cells (TIICs). Datasets of differentially expressed IRGs were extracted from the GEO database (GSE16088). The functions and prognostic values of these differentially expressed IRGs were systematically investigated using a series of bioinformatics methods. In addition, CCK8 and plate clone formation assays were used to explore the effect of PGF on osteosarcoma cells, and twenty-nine differentially expressed IRGs were identified, of which 95 were upregulated and 34 were downregulated. Next, PPI was established for Identifying Hub genes and biology networks by Cytoscape. Six IRGs (APLNR, TPM2, PGF, CD86, PROCR, and SEMA4D) were used to develop an overall survival (OS) prediction model, and two IRGs (HLA-B and PGF) were used to develop a relapse-free survival (RFS) prediction model. Compared with the low-risk patients in the training cohort (GSE39058) and TARGET validation cohorts, high-risk patients had poorer OS and RFS. Using these identified IRGs, we used OS and RFS prediction nomograms to generate a clinical utility model. The risk scores of the two prediction models were associated with the infiltration proportions of some TIICs, and the activation of memory CD4 T-cells was associated with OS and RFS. CD86 was associated with CTLA4 and CD28 and influenced the infiltration of different TIICs. In vitro experiments showed that the knockdown of PGF inhibited the proliferation and viability of osteosarcoma cells. In conclusion, these findings help us better understand the prognostic roles of IRGs and TIICs in osteosarcoma, and CD86 and PGF may serve as specific immune targets.
| S-EPMC9019688 | biostudies-literature
Project description:Objective: The stratification of neuroblastoma (NBL) prognosis remains difficult. RNA-based signatures might be able to predict prognosis, but independent cross-platform validation is still rare. Methods: RNA-Seq-based profiles from NBL patients were acquired and then analyzed. The RNA-Seq prognostic index (RPI) and the clinically adjusted RPI (RCPI) were successively established in the training cohort (TARGET-NBL) and then verified in the validation cohort (GSE62564). Survival prediction was assessed using a time-dependent receiver operating characteristic (ROC) curve and area under the ROC curve (AUC). Functional enrichment analysis of the genes was conducted using bioinformatics methods. Results: In the training cohort, 10 gene pairs were eventually integrated into the RPI. In both cohorts, the high-risk group had poor overall survival (OS) (P < 0.001 and P < 0.001, respectively) and favorable event-free survival (EFS) (P = 0.00032 and P = 0.06, respectively). ROC curve analysis also showed that the RPI predicted OS (60 month AUC values of 0.718 and 0.593, respectively) and EFS (60 month AUC values of 0.627 and 0.852, respectively) well in both the training and validation cohorts. Clinicopathological indicators associated with prognosis in the univariate and multivariate regression analyses were identified and added to the RPI to form the RCPI. The RCPI was also used to divide populations into different risk groups, and the high-risk group had poor OS (P < 0.001 and P < 0.001, respectively) and EFS (P < 0.05 and P < 0.05, respectively). Finally, the RCPI had higher accuracy than the RPI for the prediction of OS (60 month AUC values of 0.730 and 0.852, respectively) and EFS (60 month AUC values of 0.663 and 0.763, respectively) in both the training and validation cohorts. Moreover, these differentially expressed genes may be involved in certain NBL-related events. Conclusions: The RCPI could reliably categorize NBL patients based on different risks of death.
| S-EPMC6904333 | biostudies-literature
Project description:Dysregulated expression of RNA-binding proteins (RBPs) is strongly associated with the development and progression of multiple tumors. However, little is known about the role of RBPs in kidney renal clear cell carcinoma (KIRC). In this study, we examined RBP expression profiles using The Cancer Genome Atlas database and identified 133 RBPs that were differentially expressed in KIRC and non-tumor tissues. We then systematically analyzed the potential biological functions of these RBPs and established PPIs. Based on Lasso regression and Cox survival analyses, we constructed a risk model that could independently and accurately predict prognosis based on seven RBPs (NOL12, PABPC1L, RNASE2, RPL22L1, RBM47, OASL, and YBX3). Survival times were shorter in patients with high risk scores for cohorts stratified by different characteristics. Gene set enrichment analysis was also performed to further understand functional differences between high- and low-risk groups. Finally, we developed a clinical nomogram with a concordance index of 0.792 for estimating 3- and 5-year survival probabilities. Our results demonstrate that this risk model could potentially improve individualized diagnostic and therapeutic strategies.
| S-EPMC7803486 | biostudies-literature