Systematic Analysis of Alternative Splicing Landscape in Pancreatic Adenocarcinoma Reveals Regulatory Network Associated with Tumorigenesis and Immune Response.

BACKGROUND Pancreatic ductal adenocarcinoma (PDAC) is one of the most aggressive gastrointestinal tumors and has an extremely high mortality rate. Recent studies indicate that alternative splicing (AS), a common post-transcriptional process, has important roles in tumor biological behaviors and may provide novel immunotherapeutic targets. This study systematically analyzes AS profiles in PDAC and reveals their potential regulatory effects on cancer immune response. MATERIAL AND METHODS AS event, RNA sequencing, and splicing factor (SF) data were extracted from SpliceSeq, The Cancer Genome Atlas, and SpliceAid2, respectively. Overall survival (OS)-associated AS events and SFs were identified with univariate analysis. The LASSO method and multivariate Cox regression analysis were used to construct predictive signatures for the prediction of patient prognosis. The proportions of immune cells within PDAC samples were evaluated using the CIBERSORT algorithm. The correlations among AS events, SFs, and immune cell proportions were calculated using Spearman correlation analysis. Consensus clustering and immune classification were performed on the PDAC cohort. RESULTS A total of 4812 OS-related AS events from 3341 parent genes were identified, and 8 AS-based predictive models were constructed for PDAC. An OS-related SF-AS regulatory network was constructed. The AS events regulated by ELAVL4 exhibited strong correlations with CD8 T cells and regulatory T cells. In addition, AS-based clusters demonstrated distinct OS outcomes and immune features. CONCLUSIONS AS-based predictive models with high accuracy were constructed to facilitate prognosis prediction and treatment of PDAC. An SF-AS regulatory network was constructed, revealing the potential relationships among SF, AS, and immune response.