Abstract
Background: This study was to explore differential RNA splicing patterns and elucidate the function of the splice variants served as prognostic biomarkers in colorectal cancer (CRC). Methods: Genome-wide profiling of prognostic alternative splicing (AS) events using RNA-seq data from The Cancer Genome Atlas (TCGA) program was conducted to evaluate the roles of seven AS patterns in 330 colorectal cancer cohort. The prognostic predictors models were assessed by integrated Cox proportional hazards regression. Based on the correlations between survival associated AS events and splicing factors, splicing networks were built. Results: A total of 2,158 survival associated AS events in CRC were identified. Interestingly, most of these top 20 survival associated AS events were adverse prognostic factors. The prognostic models were built by each type of splicing patterns, performing well for risk stratification in CRC patients. The area under curve (AUC) of receiver operating characteristic (ROC) for the combined prognostic predictors model could reach 0.963. Splicing network also suggested distinguished correlation between the expression of splicing factors and AS events in CRC patients. Conclusion: The ideal prognostic predictors model for risk stratification in CRC patients was constructed by differential splicing patterns of 13 genes. Our findings enriched knowledge about differential RNA splicing patterns and the regulation of splicing, providing generous biomarker candidates and potential targets for the treatment of CRC.