Immunological significance of prognostic alternative splicing signature in hepatocellular carcinoma.
Ontology highlight
ABSTRACT: Background
Hepatocellular 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.Methods
Univariate 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.Results
A 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).Conclusion
Herein, 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.
SUBMITTER: Xu Q
PROVIDER: S-EPMC8017877 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
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:Alternative splicing (AS) is a pervasive and vital mechanism involved in the progression of cancer by expanding genomic encoding capacity and increasing protein complexity. However, the systematic analysis of AS in hepatocellular carcinoma (HCC) is lacking and urgently required. In the present study, genome‑wide AS events with corresponding clinical information were profiled in 290 patients with HCC from the Cancer Genome Atlas and SpliceSeq software. Functional enrichment analyses revealed the pivotal biological process of AS regulation. Univariate Cox regression analyses were performed, followed by stepwise forward multivariate analysis to develop the prognostic signatures. Spearman's correlation analyses were also used to construct potential regulatory network between the AS events and aberrant splicing factors. A total of 34,163 AS events were detected, among which 1,805 AS events from 1,314 parent genes were significantly associated with the overall survival (OS) of patients with HCC, and their parent genes serve crucial roles in HCC‑related oncogenic processes, including the p53 signaling pathway, AMPK signaling pathway and HIF‑1 signaling pathway. A prognostic AS signature was established that was found to be an independent prognostic factor for OS in stratified cohorts, harboring a noteworthy ability to distinguish between the distinct prognoses of patients with HCC (high‑risk vs. low‑risk, 827 vs. 3,125 days, P<2e‑16). Time‑dependent receiver‑-operator characteristic curves confirmed its robustness and clinical efficacy, with the area under the curves maintained >0.9 for short‑term and long‑term prognosis prediction. The splicing correlation network suggested a trend in the interactions between splicing factors and prognostic AS events, further revealing the underlying mechanism of AS in the oncogenesis of HCC. In conclusion, the present study provides a comprehensive portrait of global splicing alterations involved in the progression and HCC in addition to valuable prognostic factors for patients, which may represent as underappreciated hallmark and provide novel clues of therapeutic targets in HCC.
| S-EPMC6826324 | biostudies-literature
Project description:BackgroundBladder cancer (BLCA) is the most common malignant tumor in the genitourinary system, and the complex tumor microenvironment (TME) of BLCA is the main factor in its difficult treatment. Accumulated evidence supports that alternative splicing (AS) events frequently occur in cancer and are closely related to the TME. Therefore, there is an urgent need to comprehensively analyze the prognostic value of AS events in BLCA.MethodThe clinical, transcriptome and AS data of BLCA were downloaded from the Cancer Genome Atlas database, and a Cox proportional hazard regression model and LASSO regression were used to establish a prognostic signature. Then, the prognostic value of the signature was verified by clinical survival status, clinicopathologic features, tumor immune microenvironment (TIME), and immune checkpoint. Next, we screened the AS-related genes with the largest expression differences between tumor and normal samples by gene differential expression analysis. Finally, the regulatory network of AS-splicing factors (SFs) was established to unravel the potential regulatory mechanism of AS events in BLCA.ResultsA BLCA prognostic signature related to seven AS events was constructed, and the prognostic value of the signature was also verified from multiple perspectives. Moreover, there was significant abnormal expression of PTGER3, a gene implicated in AS events, the expression of which was associated with the survival, clinicopathological features, TIME, and immunotherapy of BLCA, suggesting that it has potential clinical application value. Furthermore, the AS-SF regulatory network indicated that splicing factors (PRPF39, LUC7L, HSPA8 and DDX21) might be potential biomarkers of BLCA.ConclusionsOur study revealed the potential role of AS events in the prognosis, TIME and immunotherapy of BLCA and yielded new insights into the molecular mechanisms of and personalized immunotherapy for BLCA.
| S-EPMC8873598 | biostudies-literature
Project description:BackgroundRenal cell carcinoma (RCC) accounts for 90% of renal cancers, of which clear cell carcinoma (ccRCC) is the most usual histological type. The process of alternative splicing (AS) contributes to protein diversity, and the dysregulation of protein diversity may have a great influence on tumorigenesis. We developed a prognostic signature and comprehensively analyzed the role of tumor immune microenvironment (TIME) and immune checkpoint blocking (ICB) treatment in ccRCC.MethodsTo identify prognosis-related AS events, univariate Cox regression was used and functional annotation was performed using gene set enrichment analysis (GSEA). In this study, prognostic signatures were developed based on multivariate Cox, univariate Cox, and LASSO regression models. Moreover, to assess the prognostic value, the proportional hazards model, Kruskal-Wallis analysis, and ROC curves were used. To obtain a better understanding of TIME in ccRCC, the ESTIMATE R package, single sample gene set enrichment analysis (ssGSEA) algorithm, CIBERSORT method, and the tumor immune estimation resource (TIMER) were applied. The database was searched to verify the expression of C4OF19 in tumor and normal samples. Regulatory networks for AS-splicing factors (SFs) were visualized using Cytoscape 3.9.1.ResultsThere were 9,347 AS cases associated with the survival of ccRCC patients screened. A total of eight AS prognostic signatures were developed with stable prognostic predictive accuracy based on splicing subtypes. In addition, a qualitative prognostic nomogram was developed, and the prognostic prediction showed high effectiveness. In addition, we found that the combined signature was significantly associated with the diversity of TIME and ICB treatment-related genes. C4ORF19 might become an important prognostic factor for ccRCC. Finally, the AS-SF regulatory network was established to clearly reveal the potential function of SFs.ConclusionWe found novel and robust indicators (i.e., risk signature, prognostic nomogram, etc.) for the prognostic prediction of ccRCC. A new and reliable prognostic nomogram was established to quantitatively predict the clinical outcome. The AS-SF networks could provide a new way for the study of potential regulatory mechanisms, and the important roles of AS events in the context of TIME and immunotherapy efficiency were exhibited. C4ORF19 was found to be a vital gene in TIME and ICB treatment.
| S-EPMC10824562 | biostudies-literature
Project description:ObjectivesBreast cancer (BC) currently has the highest incidence rate. Epigenetic regulation could alter gene expression and is closely related to BC initiation. This study aimed to develop an alternative splicing (AS)-based prognostic signature and clarify its relevance to the tumor immune microenvironment (TIME) status and immunotherapy of BC.MethodsCox regression analysis was conducted to screen for prognosis-related AS events. Thereafter, LASSO with Cox regression analyses was designed to construct a prognostic signature model. Kaplan-Meier survival analysis, receiver operating characteristic curves, and proportional hazard model were then utilized to confirm the prognostic value. Multiple methods were employed to reveal the context of TIME in BC. QPCR, western blotting and immunofluorescence microscopy were carried out to detect myc-associated zinc finger protein (MAZ) expression in different cell lines, and BC and paracancerous tissues.ResultsA total of 1,787 prognosis-related AS events were screened. Eight AS prognostic signatures were constructed with robust predictive accuracy based on the splicing subtypes. Furthermore, the establishment of a quantitative prognostic nomogram and consolidated signature was significantly correlated with TIME diversity and immune checkpoint blockade-related genes. MAZ was detected to be upregulated in BC. Finally, a newly constructed splicing regulatory network model revealed the potential functions of splicing factors.ConclusionsAS-splicing factor networks may enable a new approach to investigating potential regulatory mechanisms. Moreover, pivotal players in AS events with regards to TIME and efficiency of immunotherapy were uncovered and could facilitate clinical decision-making and individual determination of BC prognosis.
| S-EPMC9274553 | biostudies-literature
Project description:Alternative splicing events (ASEs) play a role in cancer development and progression. We investigated whether ASEs are prognostic for overall survival (OS) in hepatocellular carcinoma (HCC). RNA sequencing data was obtained for 343 patients included in The Cancer Genome Atlas. Matched splicing event data for these patients was then obtained from the TCGASpliceSeq database, which includes data for seven types of ASEs. Univariate and multivariate Cox regression analysis demonstrated that 3,814 OS-associated splicing events (OS-SEs) were correlated with OS. Prognostic indices were developed based on the most significant OS-SEs. The prognostic index based on all seven types of ASEs (PI-ALL) demonstrated superior efficacy in predicting OS of HCC patients at 2,000 days compared to those based on single ASE types. Patients were stratified into two risk groups (high and low) based on the median prognostic index. Kaplan-Meier survival analysis demonstrated that PI-ALL had the greatest capacity to distinguish between patients with favorable vs. poor outcomes. Finally, univariate Cox regression analysis demonstrated that the expression of 23 splicing factors was correlated with OS-SEs in the HCC cohort. Our data indicate that a prognostic index based on ASEs is prognostic for OS in HCC.
| S-EPMC6660027 | biostudies-literature
Project description:Objectives: Bladder cancer is the fourth most common malignancy in men in the United States. Aberrant alternative splicing (AS) events are involved in the carcinogenesis, but the association between AS and bladder cancer remains unclear. This study aimed to construct an AS-based prognostic signature and elucidate the role of the tumor immune microenvironment (TIME) and the response to immunotherapy and chemotherapy in bladder cancer. Methods: Univariate Cox regression analysis was performed to detect prognosis-related AS events. The least absolute shrinkage and selection operator (LASSO) and multivariate Cox analyses were employed to build prognostic signatures. Kaplan-Meier survival analysis, multivariate Cox regression analysis, and receiver operating characteristic (ROC) curves were conducted to validate the prognostic signatures. Then, the Estimation of Stromal and Immune cells in MAlignant Tumor tissues using Expression data (ESTIMATE) and tumor immune estimation resource (TIMER) databases were searched and the single-sample gene set enrichment analysis (ssGSEA) algorithm and CIBERSORT method were performed to uncover the context of TIME in bladder cancer. The Tumor Immune Dysfunction and Exclusion (TIDE) web tool and pRRophetic algorithm were used to predict the response to immunotherapy and chemotherapy. Finally, we constructed a correlation network between splicing factors (SFs) and survival-related AS events. Results: A total of 4684 AS events were significantly associated with overall survival in patients with bladder cancer. Eight prognostic signatures of bladder cancer were established, and a clinical survival prediction model was built. In addition, the consolidated prognostic signature was closely related to immune infiltration and the response to immunotherapy and chemotherapy. Furthermore, the correlation identified EIF3A, DDX21, SDE2, TNPO1, and RNF40 as hub SFs, and function analysis found ubiquitin-mediated proteolysis is correlated most significantly with survival-associated AS events. Conclusion: Our findings highlight the prognostic value of AS for patients with bladder cancer and reveal pivotal players of AS events in the context of TIME and the response to immunotherapy and chemotherapy, which may be important for patient management and treatment.
| S-EPMC9083046 | biostudies-literature
Project description:Background: Hepatocellular carcinoma (HCC) is a tumor with high morbidity and high mortality worldwide. DNA methylation, one of the most common epigenetic changes, might serve a vital regulatory role in cancer. Methods: To identify categories based on DNA methylation data, consensus clustering was employed. The risk signature was yielded by systematic bioinformatics analyses based on the remarkably methylated CpG sites of cluster 1. Kaplan-Meier analysis, variable regression analysis, and ROC curve analysis were further conducted to validate the prognosis predictive ability of risk signature. Gene set enrichment analysis (GSEA) was performed for functional annotation. To uncover the context of tumor immune microenvironment (TIME) of HCC, we employed the ssGSEA algorithm and CIBERSORT method and performed TIMER database exploration and single-cell RNA sequencing analysis. Additionally, quantitative real-time polymerase chain reaction was employed to determine the LRRC41 expression and preliminarily explore the latent role of LRRC41 in prognostic prediction. Finally, mutation data were analyzed by employing the "maftools" package to delineate the tumor mutation burden (TMB). Results: HCC samples were assigned into seven subtypes with different overall survival and methylation levels based on 5'-cytosine-phosphate-guanine-3' (CpG) sites. The risk prognostic signature including two candidate genes (LRRC41 and KIAA1429) exhibited robust prognostic predictive accuracy, which was validated in the external testing cohort. Then, the risk score was significantly correlated with the TIME and immune checkpoint blockade (ICB)-related genes. Besides, a prognostic nomogram based on the risk score and clinical stage presented powerful prognostic ability. Additionally, LRRC41 with prognostic value was corroborated to be closely associated with TIME characterization in both expression and methylation levels. Subsequently, the correlation regulatory network uncovered the potential targets of LRRC41 and KIAA1429. Finally, the methylation level of KIAA1429 was correlated with gene mutation status. Conclusion: In summary, this is the first to identify HCC samples into distinct clusters according to DNA methylation and yield the CpG-based prognostic signature and quantitative nomogram to precisely predict prognosis. And the pivotal player of DNA methylation of genes in the TIME and TMB status was explored, contributing to clinical decision-making and personalized prognosis monitoring of HCC.
| S-EPMC8187884 | 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:Accumulating evidence indicates an unexpected role of aberrant splicing in hepatocellular carcinoma (HCC) that has been seriously neglected in previous studies. There is a need for a detailed analysis of alternative splicing (AS) and its underlying biological and clinical relevance in HCC. In this study, clinical information and corresponding RNA sequencing data of HCC patients were obtained from The Cancer Genome Atlas. Percent spliced in (PSI) values and transcriptional splicing patterns of genes were determined from the original RNA sequencing data using SpliceSeq. Then, based on the PSI values of AS events in different patients, a series of bioinformatics methods was used to identify differentially expressed AS events (DEAS), determine potential regulatory relationships, and investigate the correlation between DEAS and the patients' clinicopathological features. Finally, 25,934 AS events originating from 8,795 genes were screened with high reliability; 263 of these AS events were identified as DEAS. The parent genes of these DEAS formed an intricate network with roles in the regulation of cancer-related pathway and liver metabolism. In HCC, 36 splicing factors were involved in the dysregulation of part DEAS, 100 DEAS events were correlated with overall survival, and 71 DEAS events were correlated with disease-free survival. Stratifying HCC patients according to DEAS resulted in four clusters with different survival patterns. Significant variations in AS occurred during HCC initiation and maintenance; these are likely to be vital both for biological processes and in prognosis. The HCC-related AS events identified here and the splicing networks constructed will be valuable in deciphering the underlying role of AS in HCC.
| S-EPMC7432180 | biostudies-literature