Comprehensive Analysis of Prognostic Alternative Splicing Signature Reveals Recurrence Predictor for Papillary Thyroid Cancer
Ontology highlight
ABSTRACT: Background
Alternative splicing (AS) plays a key role in the diversity of proteins and is closely associated with tumorigenicity. The aim of this study was to systemically analyze RNA alternative splicing (AS) and identify its prognostic value for papillary thyroid cancer (PTC). Methods
AS percent-splice-in (PSI) data of 430 patients with PTC were downloaded from the TCGA SpliceSeq database. We successfully identified recurrence-free survival (RFS)-associated AS events through univariate Cox regression, LASSO regression and multivariate regression and then constructed different types of prognostic prediction models. Gene function enrichment analysis revealed the relevant signaling pathways involved in RFS-related AS events. Simultaneously, a regulatory network diagram of AS and splicing factors (SFs) was established. Results
We identified 1397 RFS-related AS events which could be used as the potential prognostic biomarkers for PTC. Based on these RFS-related AS events, we constructed a ten-AS event prognostic prediction signature that could distinguish high-and low-risk patients and was highly capable of predicting PTC patient prognosis. ROC curve analysis revealed the excellent predictive ability of the ten-AS events model, with an area under the curve (AUC) value of 0.889; the highest prediction intensity for one-year RFS was 0.923, indicating that the model could be used as a prognostic biomarker for PTC. In addition, the nomogram constructed by the risk score of the ten-AS model also showed high predictive efficiency for the prognosis of PTC patients. Finally, the constructed SF-AS network diagram revealed the regulatory role of SFs in PTC. Conclusion
Through the limited analysis, AS events could be regarded as reliable prognostic biomarkers for PTC. The splicing correlation network also provided new insight into the potential molecular mechanisms of PTC.
SUBMITTER: Liu M
PROVIDER: S-EPMC8548661 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:AimCervical cancer is a common malignant carcinoma of the gynecological tract with high morbidity and mortality. Therefore, it is crucial to elucidate the pathogenesis, prevention, diagnosis and prognosis of cervical cancer by searching for the involved key genes.MethodIn this study, the alternative splicing (AS) events of 253 patients with cervical cancer were analyzed, and 41,766 AS events were detected in 9961 genes. Univariate analysis was performed to screen prognostic AS events. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis was used to identify the pathways in which these AS events were involved.ResultsWe found that exon skip (ES) is the main AS event in patients with cervical cancer. There was pronounced consistency between the genes involved in overall survival and those involved in recurrence. At the same time, we found that a gene may exhibit several different types of AS events, and these different AS events may be related to prognosis. Four characteristic genes, HSPA14, SDHAF2, CAMKK2 and TM9SF1, that can be used as prognostic markers for cervical cancer were selected.ConclusionThe importance of AS events in the development of cervical cancer and prediction of prognosis was revealed by a large amount of data at the whole genome level, which may provide a potential target for cervical cancer treatment. We also provide a new method for exploring the pathogenesis of cervical cancer to determine clinical treatment and prognosis more accurately.
| S-EPMC6865056 | biostudies-literature
Project description:BackgroundAlternative splicing (AS) is a molecular event that drives protein diversity through the generation of multiple mRNA isoforms. Growing evidence demonstrates that dysregulation of AS is associated with tumorigenesis. However, an integrated analysis in identifying the AS biomarkers attributed to esophageal carcinoma (ESCA) is largely unexplored.MethodsAS percent-splice-in (PSI) data were obtained from the TCGA SpliceSeq database. Univariate and multivariate Cox regression analysis was successively performed to identify the overall survival (OS)-associated AS events, followed by the construction of AS predictor through different splicing patterns. Then, a nomogram that combines the final AS predictor and clinicopathological characteristics was established. Finally, a splicing regulatory network was created according to the correlation between the AS events and the splicing factors (SF).ResultsWe identified a total of 2389 AS events with the potential to be used as prognostic markers that are associated with the OS of ESCA patients. Based on splicing patterns, we then built eight AS predictors that are highly capable in distinguishing high- and low-risk patients, and in predicting ESCA prognosis. Notably, the area under curve (AUC) value for the exon skip (ES) prognostic predictor was shown to reach a score of 0.885, indicating that ES has the highest prediction strength in predicting ESCA prognosis. In addition, a nomogram that comprises the pathological stage and risk group was shown to be highly efficient in predicting the survival possibility of ESCA patients. Lastly, the splicing correlation network analysis revealed the opposite roles of splicing factors (SFs) in ESCA.ConclusionIn this study, the AS events may provide reliable biomarkers for the prognosis of ESCA. The splicing correlation networks could provide new insights in the identification of potential regulatory mechanisms during the ESCA development.
| S-EPMC7387693 | biostudies-literature
Project description:BackgroundAlternative splicing (AS) is a key factor in protein-coding gene diversity, and is associated with the development and progression of malignant tumours. However, the role of AS in cervical cancer is unclear.MethodsThe AS data for cervical squamous cell carcinoma and endocervical adenocarcinoma (CESC) were downloaded from The Cancer Genome Atlas (TCGA) SpliceSeq website. Few prognostic AS events were identified through univariate Cox analysis. We further identified the prognostic prediction models of the seven subtypes of AS events and assessed their predictive power. We constructed a clinical prediction model through global analysis of prognostic AS events and established a nomogram using the risk score calculated from the prognostic model and relevant clinical information. Unsupervised cluster analysis was used to explore the relationship between prognostic AS events in the model and clinical features.ResultsA total of 2860 prognostic AS events in cervical cancer were identified. The best predictive effect was shown by a single alternate acceptor subtype with an area under the curve of 0.96. Our clinical prognostic model included a nine-AS event signature, and the c-index of the predicted nomogram model was 0.764. SNRPA and CCDC12 were hub genes for prognosis-associated splicing factors. Unsupervised cluster analysis through the nine prognostic AS events revealed three clusters with different survival patterns.ConclusionsAS events affect the prognosis and biological progression of cervical cancer. The identified prognostic AS events and splicing regulatory networks can increase our understanding of the underlying mechanisms of cervical cancer, providing new therapeutic strategies.
| S-EPMC7282181 | biostudies-literature
Project description:Background: Pancreatic cancer (PC) is an aggressive cancer with worse survival in the world. Emerging evidence suggested that the imbalance of alternative splicing (AS) is a hallmark of cancer and indicated poor prognosis of patients. Genes-derived splicing events can produce neoepitopes for immunotherapy. However, the profound study of splicing profiling in PC is still elusive. We aimed to identification of novel prognostic signature across a comprehensive splicing landscape and reveal their relationship with tumor-infiltrating immune cells in pancreatic cancer microenvironment. Methods: Based on integrated analysis of splicing profiling and clinical data, differentially splicing events were filtered out. Then, stepwise Cox regression analysis was applied to identify survival-related splicing events and construct prognostic signature. Functional enrichment analysis was performed to explore biology function. Kaplan-Meier curves and receiver operating characteristic (ROC) curves were performed to validate the predictive effect of predictive signature. We also verified the clinical value of prognostic signature under the influence of different clinical parameters. For deeper analysis, we evaluated the correlation between prognostic signature and infiltrating immune cells by CIBERSORT. Results: According to systematic analyzing, a final six splicing events were identified and validated the good prognostic capability in entire TCGA dataset, validation set 1 and validation set 2 by Kaplan-Meier curves (P < 0.0001). The area under the curve (AUC) of ROC curves were also confirmed the high predictive efficiency of the prognostic signature in these three cohorts (AUC = 0.857, 0.895 and 0.788). In order to validate whether prognostic signature highlights a correlation between AS and immune contexture, CIBERSORT was performed to analyze the proportion of tumor-infiltrating immune cells in PC. Based on prognostic signature, we identified survival-related immune cells including CD8 T cells (P = 0.0111), activated CD4 memory T cells (P = 0.0329) and resting mast cells (P = 0.0352). Conclusion: In conclusion, our study contribute to provide a promising prognostic signature based on six splicing events and revealed prognosis-related immune cells which indeed represented novel tumor drivers and provide potential targets for personalized therapeutic.
| S-EPMC7545682 | biostudies-literature
Project description:BackgroundThe incidence of thyroid cancer (THCA) continues to increase in recent decades. Accumulating evidence showed that the unbalanced alternative splicing (AS) promotes the occurrence of cancers and leads to poor prognosis of patients. However, the research on alternative splicing events in THCA is lacking, and its underlying mechanism is not fully understood. This study identifies a novel prognostic signature based on AS events to reveal the relationship of AS with tumor immune microenvironment.MethodsBased on the AS data, transcriptional data, and clinical information, the differentially expressed alternative splicings (DEASs) were screened out. Least absolute shrinkage and selection operator (LASSO) regression and multi-Cox regression analyses were employed to identify prognostic results related to AS events and establish a prognostic signature. The predictive ability of the signature was assessed by Kaplan-Meier (K-M) survival curve, risk plots, and receiver operating characteristic (ROC) curves. Furthermore, correlations between tumor-infiltrating immune cells, immune checkpoints, immune score and prognostic signature were analyzed.ResultsAccording to the LASSO regression analysis, a total of five AS events were selected to construct the signature. K-M survival curve showed that the higher the risk score, the worse the OS of the patients. Risk plots further confirmed this result. ROC curves indicated the high predictive efficiency of the prognostic signature. As for tumor immune microenvironment, patients in the high-risk group had a higher proportion of immune cells, including plasma cell, CD8+ T cell, macrophages (M0 and M2), and activated dendritic cell. Immune checkpoint proteins, such as PDCD1LG2, HAVCR2, CD274, etc., were significantly higher in the high-risk group. We also found that the ESTIMATE score, stromal score, and immune score were lower in the high-risk group, while the result of tumor purity was the opposite.ConclusionsCollectively, a prognostic signature consisting of five AS events in THCA was established. Furthermore, there was an inextricable correlation between immune cell infiltration, immune checkpoint proteins, and AS events. This study will provide a basis for THCA immunotherapy in the future.
| S-EPMC8558422 | biostudies-literature
Project description:Previous studies have demonstrated the role of abnormal alternative splicing (AS) in tumor progression. This study examines the prognostic index (PI) of alternative splices (ASs) in patients with hepatocellular carcinoma (HCC). The clinical features and splicing events of patients with HCC were downloaded from The Cancer Genome Atlas (TCGA). Differentially expressed AS (DEAS) were compared between HCC and adjacent normal samples. Univariate Cox regression analysis was used to determine changes in DEAS associated with overall survival (OS). A PI was generated from OS-associated DEASs using Kaplan-Meier curves, receiver operating characteristic (ROC) curves, multivariate Cox regression, and cluster analysis. Then, the correlation between DEASs and splicing factors was assessed, followed by functional and pathway enrichment analysis. We identified 34 163 ASs of 8985 genes in HCC, and 153 OS-ASs were identified using univariate Cox regression analysis. Low- and high-PI groups were determined based on the median "PI-ALL" value according to significantly different survival (P = 2.2e - 16). The ROC curve of all PI (PI-ALL) had an area under the curve (AUC) of 0.993 for survival status in patients with HCC. A potential regulatory network associated with prognosis of patients with HCC was established. Enrichment analysis also resulted in the identification of several pathways potentially associated with carcinogenesis and progression of HCC. Four clusters were identified that were associated with clinical features and prognosis. Our study generated comprehensive profiles of ASs in HCC. The interaction network and functional connections were used to elucidate the underlying mechanisms of AS in HCC.
| S-EPMC7064038 | biostudies-literature
Project description:BackgroundIncreasing evidence indicated a close relationship between aberrant splicing variants and carcinoma, whereas comprehensive analysis of prognostic alternative splicing (AS) profiling in breast cancer (BRCA) is lacking and largely unknown.MethodsRNA-seq data and corresponding clinical information of BRCA patients were obtained and integrated from The Cancer Genome Atlas (TCGA). Then SpliceSeq software was used to assess seven AS types and calculate the Percent Spliced In (PSI) value. Univariate followed by stepwise multivariate Cox regression analyses identified survival associated AS events and constructed the AS signature, which were further sent for enrichment analysis, respectively. Besides, the splicing correlation network was constructed. Additionally, nomogram incorporating AS signature and clinicopathological characteristics was developed and its efficacy was evaluated with respect to discrimination, calibration and clinical utility.ResultsA total of 45,421 AS events were detected, among which 3071 events were found associated with overall survival (OS) after strict filtering. Parent genes of these prognostic events were involved in BRCA-related processes including NF-kappaB and HIF-1 signaling pathway. Besides, the final prognostic signature built with 20 AS events performed well with an area under the curve (AUC) of receiver operating characteristic (ROC) curve up to 0.957 for 5 years. And gene set enrichment analysis (GSEA) also confirmed the candidate 20 AS events contributed to progression of BRCA. Moreover, the nomogram that incorporated 20-AS-event-based classifier, age, pathological stage and Her-2 status showed good calibration and moderate discrimination, with C-index of 0.883 (95% CI, 0.844-0.921). Decision curve analysis (DCA) confirmed more benefit was added to survival prediction with our nomogram, especially in 5 or 8 years with threshold probability up to 80%. Finally, splicing correlation network revealed an obvious regulatory pattern of prognostic splicing factors (SF) in BRCA.ConclusionThis study provided a systematic portrait of survival-associated AS events involved in BRCA and further presented a AS-clinicopathological nomogram, which could be conveniently used to assist the individualized prediction of long-term survival probability for BRCA patients. And a series of bioinformatic analysis provided a promising perspective for further uncovering the underlying mechanisms of AS events and validating therapeutic targets for BRCA.
| S-EPMC6448481 | biostudies-literature
Project description:Background: Head and neck squamous cell carcinoma (HNSC) is a prevalent and heterogeneous malignancy with poor prognosis and high mortality rates. There is significant evidence of alternative splicing (AS) contributing to tumor development, suggesting its potential in predicting prognosis and therapeutic efficacy. This study aims to establish an AS-based prognostic signature in HNSC patients. Methods: The expression profiles and clinical information of 486 HNSC patients were downloaded from the TCGA database, and the AS data were downloaded from the TCGA SpliceSeq database. The survival-associated AS events were identified by conducting a Cox regression analysis and utilized to develop a prognostic signature by fitting into a LASSO-regularized Cox regression model. Survival analysis, univariate and multivariate Cox regression analysis, and receiver operating characteristic (ROC) curve analysis were performed to evaluate the signature and an independent cohort was used for validation. The immune cell function and infiltration were analyzed by CIBERSORT and the ssGSEA algorithm. Results: Univariate Cox regression analysis identified 2726 survival-associated AS events from 1714 genes. The correlation network reported DDX39B, PRPF39, and ARGLU1 as key splicing factors (SF) regulating these AS events. Eight survival-associated AS events were selected and validated by LASSO regression to develop a prognostic signature. It was confirmed that this signature could predict HNSC outcomes independent of other variables via multivariate Cox regression analysis. The risk score AUC was more than 0.75 for 3 years, highlighting the signature's prediction capability. Immune infiltration analysis reported different immune cell distributions between the two risk groups. The immune cell content was higher in the high-risk group than in the low-risk group. The correlation analysis revealed a significant correlation between risk score, immune cell subsets, and immune checkpoint expression. Conclusion: The prognostic signature developed from survival-associated AS events could predict the prognosis of HNSC patients and their clinical response to immunotherapy. However, this signature requires further research and validation in larger cohort studies.
| S-EPMC9633855 | biostudies-literature
Project description:Purpose:Papillary thyroid carcinoma (PTC), the most frequent type of malignant thyroid tumor, lacks novel and reliable biomarkers of patients' prognosis. In the current study, we mined The Cancer Genome Atlas (TCGA) to develop lncRNA signature of PTC. Patients and methods:The intersection of PTC lncRNAs was obtained from the TCGA database using integrative computational method. By the univariate and multivariate Cox analysis, key lncRNAs were identified to construct the prognostic model. Then, all patients were divided into the high-risk group and low-risk group to perform the Kaplan-Meier (K-M) survival curves and time-dependent receiver operating characteristic (ROC) curve, estimating the prognostic power of the prognostic model. Functional enrichment analysis was also performed. Finally, we verified the results of the TCGA analysis by the Gene Expression Omnibus (GEO) databases and quantitative real-time PCR (qRT-PCR). Results:After the comprehensive analysis, a three-lncRNA signature (PRSS3P2, KRTAP5-AS1 and PWAR5) was obtained. Interestingly, patients with low-risk scores tended to gain obviously longer survival time, and the area under the time-dependent ROC curve was 0.739. Furthermore, gene ontology (GO) and pathway analysis revealed the tumorigenic and prognostic function of the three lncRNAs. We also found three potential transcription factors to help understand the mechanisms of the PTC-specific lncRNAs. Finally, the GEO databases and qRT-PCR validation were consistent with our TCGA bioinformatics results. Conclusion:We built a three-lncRNA signature by mining the TCGA database, which could effectively predict the prognosis of PTC.
| S-EPMC6183593 | biostudies-literature
Project description:BackgroundHepatocellular carcinoma (HCC) ranks the sixth prevalent tumors with high mortality globally. Alternative splicing (AS) drives protein diversity, the imbalance of which might act an important factor in tumorigenesis. This study aimed to construct of AS-based prognostic signature and elucidate the role in tumor immune microenvironment (TIME) and immunotherapy in HCC.MethodsUnivariate Cox regression analysis was performed to determine the prognosis-related AS events and gene set enrichment analysis (GSEA) was employed for functional annotation, followed by the development of prognostic signatures using univariate Cox, LASSO and multivariate Cox regression. K-M survival analysis, proportional hazards model, and ROC curves were conducted to validate prognostic value. ESTIMATE R package, ssGSEA algorithm and CIBERSORT method and TIMER database exploration were performed to uncover the context of TIME in HCC. Quantitative real-time polymerase chain reaction was implemented to detect ZDHHC16 mRNA expression. Cytoscape software 3.8.0 were employed to visualize AS-splicing factors (SFs) regulatory networks.ResultsA total of 3294 AS events associated with survival of HCC patients were screened. Based on splicing subtypes, eight AS prognostic signature with robust prognostic predictive accuracy were constructed. Furthermore, quantitative prognostic nomogram was developed and exhibited robust validity in prognostic prediction. Besides, the consolidated signature was significantly correlated with TIME diversity and ICB-related genes. ZDHHC16 presented promising prospect as prognostic factor in HCC. Finally, the splicing regulatory network uncovered the potential functions of splicing factors (SFs).ConclusionHerein, exploration of AS patterns may provide novel and robust indicators (i.e., risk signature, prognostic nomogram, etc.,) for prognostic prediction of HCC. The AS-SF networks could open up new approach for investigation of potential regulatory mechanisms. And pivotal players of AS events in context of TIME and immunotherapy efficiency were revealed, contributing to clinical decision-making and personalized prognosis monitoring of HCC.
| S-EPMC8017877 | biostudies-literature