Dataset Information

Identification of tumour immune microenvironment-related alternative splicing events for the prognostication of pancreatic adenocarcinoma

ABSTRACT:

Purpose

Pancreatic adenocarcinoma (PAAD) is characterized by low antitumour immune cell infiltration in an immunosuppressive microenvironment. This study aimed to systematically explore the impact on prognostic alternative splicing events (ASs) of tumour immune microenvironment (TIME) in PAAD.

Methods

The ESTIMATE algorithm was implemented to compute the stromal/immune-related scores of each PAAD patient, followed by Kaplan–Meier (KM) survival analysis of patients with different scores grouped by X-tile software. TIME-related differentially expressed ASs (DEASs) were determined and evaluated through functional annotation analysis. In addition, Cox analyses were implemented to construct a TIME-related signature and an AS clinical nomogram. Moreover, comprehensive analyses, including gene set enrichment analysis (GSEA), immune infiltration, immune checkpoint gene expression, and tumour mutation were performed between the two risk groups to understand the potential mechanisms. Finally, Cytoscape was implemented to illuminate the AS-splicing factor (SF) regulatory network.

Results

A total of 437 TIME-related DEASs significantly related to PAAD tumorigenesis and the formation of the TIME were identified. Additionally, a robust TIME-related prognostic signature based on seven DEASs was generated, and an AS clinical nomogram combining the signature and four clinical predictors also exhibited prominent discrimination by ROC (0.762 ~ 0.804) and calibration curves. More importantly, the fractions of CD8 T cells, regulatory T cells and activated memory CD4 T cells were lower, and the expression of four immune checkpoints—PD-L1, CD47, CD276, and PVR—was obviously higher in high-risk patients. Finally, functional analysis and tumour mutations revealed that aberrant immune signatures and activated carcinogenic pathways in high-risk patients may be the cause of the poor prognosis.

Conclusion

We extracted a list of DEASs associated with the TIME through the ESTIMATE algorithm and constructed a prognostic signature on the basis of seven DEASs to predict the prognosis of PAAD patients, which may guide advanced decision-making for personalized precision intervention.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12885-021-08962-7.

SUBMITTER: Chen B

PROVIDER: S-EPMC8590242 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:BackgroundAlternative splicing (AS) is a gene regulatory mechanism that drives protein diversity. Dysregulation of AS is thought to play an essential role in cancer initiation and development. This study aimed to construct a prognostic signature based on AS and explore the role in the tumor immune microenvironment (TIME) in lung adenocarcinoma.MethodsWe analyzed transcriptome profiling and clinical lung adenocarcinoma data from The Cancer Genome Atlas (TCGA) database and lists of AS-related and immune-related signatures from the SpliceSeq. Prognosis-related AS events were analyzed by univariate Cox regression analysis. Gene set enrichment analyses (GSEA) were performed for functional annotation. Prognostic signatures were identified and validated using univariate and multivariate Cox regression, LASSO regression, Kaplan-Meier survival analyses, and proportional hazards model. The context of TIME in lung adenocarcinoma was also analyzed. Gene and protein expression data of Cyclin-Dependent Kinase Inhibitor 2A (CDKN2A) were obtained from ONCOMINE and Human Protein Atlas. Splicing factor (SF) regulatory networks were visualized.ResultsA total of 19,054 survival-related AS events in lung adenocarcinoma were screened in 1,323 genes. Exon skip (ES) and mutually exclusive exons (ME) exhibited the most and fewest AS events, respectively. Based on AS subtypes, eight AS prognostic signatures were constructed. Patients with high-risk scores were associated with poor overall survival. A nomogram with good validity in prognostic prediction was generated. AUCs of risk scores at 1, 2, and 3 years were 0.775, 0.736, and 0.759, respectively. Furthermore, the prognostic signatures were significantly correlated with TIME diversity and immune checkpoint inhibitor (ICI)-related genes. Low-risk patients had a higher StromalScore, ImmuneScore, and ESTIMATEScore. AS-based risk score signature was positively associated with CD8+ T cells. CDKN2A was also found to be a prognostic factor in lung adenocarcinoma. Finally, potential functions of SFs were determined by regulatory networks.ConclusionTaken together, our findings show a clear association between AS and immune cell infiltration events and patient outcome, which could provide a basis for the identification of novel markers and therapeutic targets for lung adenocarcinoma. SF networks provide information of regulatory mechanisms.

Project description:BackgroundHepatocellular carcinoma (HCC) is one of the most common malignant tumors in the world, and its 5-year survival rate is less than 20%, despite various treatments being available. Increasing evidence indicates that alternative splicing (AS) plays a nonnegligible role in the formation and development of the tumor microenvironment (TME). However, the comprehensive analysis of the impact on prognostic AS events on immune-related perspectives in HCC is lacking but urgently needed.MethodsThe transcriptional data and clinical information of HCC patients were downloaded from TCGA (The Cancer Genome Atlas) database for calculating immune and stromal scores by ESTIMATE algorithm. We then divided patients into high/low score groups and explored their prognostic significance using Kaplan-Meier curves. Based on stromal and immune scores, differentially expressed AS events (DEASs) were screened and evaluated with functional enrichment analysis. Additionally, a risk score model was established by applying univariate and multivariate Cox regression analyses. Finally, gene set variation analysis (GSVA) was adopted to explore differences in biological behaviors between the high- and low-risk subgroups.ResultsA total of 370 HCC patients with complete and qualified corresponding data were included in the subsequent analysis. According to the results of ESTIMATE analysis, we observed that the high immune/stromal score group had a longer survival probability, which was significantly correlated with prognosis in HCC patients. In addition, 467 stromal/immune score-related DEASs were identified, and enrichment analysis revealed that DEASs were significantly enriched in pathways related to HCC tumorigenesis and the immune microenvironment. More importantly, the final prognostic signature containing 16 DEASs showed powerful predictive ability. Finally, GSVA demonstrated that activation of carcinogenic pathways and immune-related pathways in the high-risk group may lead to poor prognosis.ConclusionsCollectively, these outcomes revealed prognostic AS events related to carcinogenesis and the immune microenvironment, which may yield new directions for HCC immunotherapy.