Aberrant expression and prognostic value of CXCL11 in colon cancer
As shown in Figure 1a, the Kaplan–Meier survival curve revealed that patients with high CXCL11 mRNA expression had better prognosis (p=0.0053) in TCGA. Next, TISIDB was utilized to examine the association between CXCL11 and clinical outcome across all TCGA tumors and we found that patients with high expression of CXCL11 mRNA had good prognosis especially in COAD (colon cancer) (Figure 1b). Furthermore, we explored the differential expression of CXCL11 between tumor and adjacent normal tissues across all TCGA tumors in TIMER and found that CXCL11 was significantly up-regulated in tumor than adjacent normal tissues across several tumors (Figure 1c). In addition, we also validated the up-regulated expression of CXCL11 in tumor compared to adjacent normal tissues (Figure 2a and 2b) by immunochemistry staining in the Yijishan Hospital cohort (YJSHC).
Prognostic value of CXCL11+ cells and predictive value of CXCL11+ cells for response to ACT
Patients with high abundance of intratumoral CXCL11+ cells had better OS (overall survival, p=0.001) within the YJSHC, as shown in Figure 2c. Meanwhile, in Figure 2d, it could be observed that the effect of intratumoral CXCL11+cells infiltration kept significant after confounders were adjusted (HR= 0.35; 95% CI 0.17-0.72; p= 0.004).
Association of CXCL11 with tumor immune microenvironment
Firstly, the differential tumor immune microenvironment in different groups of CXCL11 mRNA expression in TCGA was investigated by CIBERSORT. High CXCL11 mRNA expression group had higher proportion of anti-tumor immune cells, such as: CD8+ T cells (CD8T, p<0.001), activated natural killer cells (NKa, p<0.05), as shown in Figure 3a. On the contrary, high CXCL11 mRNA expression group had lower proportion of pro-tumor immune cells, such as: M0 macrophages (M0, p<0.05), resting natural killer cells (NKr, p<0.001), monocytes (Mono, p<0.05).
Next, TISIDB[28] was employed to validate which kinds of intratumoral lymphocytes might be regulated by CXCL11 across all TCGA tumors (supplementary Figure S1a). Surprisingly, we found that CXCL11 positively correlates with activated CD8+ T cells (Act CD8, r=0.55, p<0.001; Figure 3b), natural killer T cells (NKT, r=0.438, p<0.001; Figure 3c), activated dendritic cells (Act DC, r=0.488, p<0.001; Figure 3d) in colon cancer.
Genetic alteration and enriched biological process in the high and low CXCL11 metagene expression groups
Firstly, we performed GSEA and identified that biological pathways enriched in high-CXCL11 tumor samples were immune-activated process, which consisted of chemokine signaling pathway, cytokine-cytokine receptor interaction, antigen processing and presentation, T cell receptor signaling pathway,natural killer cell mediated cytotoxicity, intestinal immune network for IgA production, interferon alpha response,interferon gamma response , as shown in Figure 3e-3l.
We then analyzed the genetic alteration and found that differentially expressed genes in tumors revealed several immune-activated genes in high CXCL11 mRNA expression group. The volcano plot was shown in Figure 4a. As shown in Figure 4b and 4c, CXCL11 positively correlated with a gene set associated with DC, NK and T recruiting genes: CCL4 (r=0.543, p<0.001), CCL5 (r=0.59, p<0.001), CXCL9 (r=0.717, p<0.001), CXCL10 (r=0.821, p<0.001).
Cytotoxic genes [32-34] including IFNG, GZMA, GZMB, GZMK, GZMM and PRF1 were positively correlated with CXCL11 mRNA expression (Spearman’s q=0.49, 0.55, 0.22, 0.37, 0.21 and 0.37; p<0.001, p <0.001, p<0.001, p=0.045, and p=0.004; respectively), as shown in Figure 4d. Interestingly, CXCL11 positively correlated with several immunosuppressive molecules [34] including of PDCD1,PD-L1, CTLA4 (Spearman’s q=0.35, 0.56, and 0.44; p<0.001, p<0.001 and p<0.001; respectively), as shown in Figure 4e. Surprisingly, we also found the up-regulated expression of CXCL11 in tumor compared to adjacent normal tissue (Figure 4f and 4h), and the positive correlation between the expression of CXCL11 and PD-L1 (r=0.66, r=0.78; p<0.001, p<0.001; respectively; Figure 4g and 4i) in colon cancer single-cell RNA-seq dataset GSE146771 and GSE132465.
Validating the association of CXCL11+ cells with TILs and PD-L1
Further analysis in YJSHC was conducted to validate the findings in TCGA. As was expected, tumor with high abundance of CXCL11+ cells infiltration tended to have high abundance of intratumoral CD8+ T cells and CD56+ NK cells infiltration. In addition, high abundance of CXCL11+ cells infiltration was associated with high abundance of intratumoral PD-L1+ cells. Accordingly, the expression of PD-L1, CD8A, CD56 was higher in CXCL11-high group than CXCL11-low group (Figure 5b-5d) and the positive correlation between the expression of CXCL11 and PD-L1, CD8A, CD56 (r=0.62, p<0.001; r=0.34, p<0.001; r=0.32, p<0.001; respectively; Figure 5e-5g) was significant within YJSHC.