No abstract

Despite extensive research, the exact mechanisms involved in colorectal cancer (CRC) etiology and pathogenesis remain unclear. This study aimed to examine the correlation between tumor-associated alternative splicing (AS) events and tumor immune infiltration (TII) in CRC. We analyzed transcriptome profiling and clinical CRC data from The Cancer Genome Atlas (TCGA) database and lists of AS-related and immune-related signatures from the SpliceSeq and Innate databases, respectively to develop and validate a risk model of differential AS events and subsequently a TII risk model. We then conducted a two-factor survival analysis to study the association between TII and AS risk and evaluated the associations between immune signatures and six types of immune cells based on the TIMER database. Subsequently, we studied the distribution of six types of TII cells in high- and low-risk groups for seven AS events and in total. We obtained the profiles of AS events/genes for 484 patients, which included 473 CRC tumor samples and 41 corresponding normal samples, and detected 22581 AS events in 8122 genes. Exon Skip (ES) (8446) and Mutually Exclusive Exons (ME) (74) exhibited the most and fewest AS events, respectively. We then classified the 433 patients with CRC into low-risk (n = 217) and high-risk (n = 216) groups based on the median risk score in different AS events. Compared with patients with low-risk scores (mortality = 11.8%), patients with high-risk scores were associated with poor overall survival (mortality = 27.6%). The risk score, cancer stage, and pathological stage (T, M, and N) were closely correlated with prognosis in patients with CRC (P < 0.001). We identified 6479 differentially expressed genes from the transcriptome profiles of CRC and intersected 468 differential immune-related signatures. High-AS-risk and high-TII-risk predicted a poor prognosis in CRC. Different AS types were associated with different TII risk characteristics. Alternate Acceptor site (AA) and Alternate Promoter (AP) events directly affected the concentration of CD4T cells, and the level of CD8T cells was closely correlated with Alternate Terminator (AT) and Exon Skip (ES) events. Thus, the concentration of CD4T and CD8T cells in the CRC immune microenvironment was not specifically modulated by AS. However, B cell, dendritic cell, macrophage, and neutrophilic cell levels were strongly correlated with AS events. These results indicate adverse associations between AS event risk levels and immune cell infiltration density. Taken together, our findings show a clear association between tumor-associated alternative splicing and immune cell infiltration events and patient outcome and could form a basis for the identification of novel markers and therapeutic targets for CRC and other cancers in the future.

Objective. This study is aimed at exploring the association between autophagy and tumor immune infiltration (TII) in colorectal cancer (CRC). Methods and Materials. We downloaded the transcriptome profiling and clinical data for CRC from The Cancer Genome Atlas (TCGA) database and obtained the normal colon transcriptome profiling data from Genotype-Tissue Expression Project (GTEx) database. The list of autophagy-related signatures was obtained from the Human Autophagy Database. We isolated the autophagy-related genes from the CRC gene expression matrix and constructed an autophagy-related prognostic (ARP) risk model. Then, we constructed a multiROC curve to validate the prognostic ability of the ARP risk model. CIBERSORT was used to determine the fractions of 22 immune cells in each CRC sample, and the association between these TII cells and CRC clinical variables was further investigated. Finally, we estimated the association of 3 hub-ARP signatures and 20 different types of TII cell distribution. Results. We classified 447 CRC patients into 224 low-risk and 223 high-risk patients using the median ARP risk score. According to the univariate survival test results, except for gender ( P = 0.672 ), age ( P = 0.008 ), cancer stage, and pathological stage T, M, and N were closely correlated with the prognosis of CRC patients ( P < 0.001 ). Multivariate survival analysis results indicate that age and rescore were the only independent prognostic indicators with significant differences ( P < 0.05 ). After merging the immune cell distribution (by CIBERSORT) with the CRC clinical data, the results indicate that activated macrophage M0 cells exhibited the highest clinical response, which included cancer stage and stage T, N, and M. Additionally, six immune cells were closely associated with cancer stage, including regulatory T cells (Tregs), gamma delta T cells, follicular helper T cells, activated memory CD4 T cells, activated NK cells, and resting dendritic cells. Finally, we evaluated the correlation of ARP signatures with TII cell distribution. Compared with the other correlation, NRG1 and plasma cells (↑), risk score and macrophage M1 (↑), NRG1 and dendritic cell activated (↑), CDKN2A and T cell CD4 memory resting (↓), risk score and T cell CD8 (↑), risk score and T cell CD4 memory resting (↓), and DAPK1 and T cell CD4 memory activated (↓) exhibited a stronger association ( P < 0.0001 ). Conclusions. In summary, we explored the correlation between the risk of autophagy and the TII microenvironment in CRC patients. Furthermore, we integrated different CAR signatures with tumor-infiltrating immune cells and found robust associations between different levels of CAR signature expression and immune cell infiltrating density.