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:As a key mechanism, alternative splicing (AS) plays a role in the cancer initiation and development. However, in papillary thyroid cancer (PTC), data for the comprehensive AS event profile and its clinical implications are lacking. Herein, a genome-wide AS event profiling using RNA-Seq data and its correlation with matched clinical information was performed using a 389 PTC patient cohort from the project of The Cancer Genome Atlas (TCGA). We identified 1,925 cancer-associated AS events (CASEs) by comparing paired tumors and neighboring healthy tissues. Parent genes with CASEs remarkably enriched in the pathways were linked with carcinogenesis, such as P53, KRAS, IL6-JAK-STAT3, apoptosis, and MYC signaling. The regulatory networks of AS implied an obvious correlation between the expression of splicing factor and CASE. We identified eight CASEs as predictors for overall survival (OS) and disease-free survival (DFS). The established risk score model based on DFS-associated CASEs successfully predicted the prognosis of PTC patients. From the unsupervised clustering analysis results, it is found that different clusters based on AS correlated with prognosis, molecular features, and immune characteristics. Taken together, the comprehensive genome-wide AS landscape analysis in PTC showed new AS events linked with tumorigenesis and prognosis, which provide new insights for clinical monitoring and therapy for PTC.
| S-EPMC8544310 | 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: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:BackgroundAberrant alternative splicing (AS) events could be viewed as prognostic indicators in a large number of malignancies. This study aims to identify prognostic AS events, illuminate the function of the splicing variants biomarkers and provide reliable evidence for formulating public health strategies for gastric cancer (GC) surveillance.MethodsRNA-Seq data, clinical information and percent spliced in (PSI) values were available in The cancer genome atlas (TCGA) and TCGA SpliceSeq data portal. A three-step regression method was conducted to identify prognostic AS events and construct multi-AS-based signatures. The associations between prognostic AS events and splicing factors were also investigated.ResultsWe identified a total of 1,318 survival-related AS events in GC, parent genes of which were implicated in numerous oncogenic pathways. The final prognostic signatures stratified by seven types of AS events or not stratified performed well in risk prediction for GC patients. Moreover, five signatures based on AA, AD, AT, ES and RI events function as independent prognostic indicators after multivariate adjustment of other clinical variables. Splicing network also showed marked correlation between the expression of splicing factors and PSI value of AS events in GC patients.ConclusionOur findings provide a landscape of AS events and regulatory network in GC, indicating that AS events might serve as prognostic biomarkers and therapeutic targets for GC.
| S-EPMC9131537 | 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: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