Identification of hub genes associated with neutrophils infiltration in colorectal cancer

Colorectal cancer (CRC) is the leading cause of cancer-related mortality in the world. Accumulating evidence indicate that tumour infiltrating immune cells participated in cancer progression. Among them, tumour infiltrating neutrophils (TINs) are reported to play crucial role in various cancers. In this study, we used CIBERSORTx, a digital cytometry tool to evaluate the neutrophils infiltration in CRC based on gene expression data of CRC tissues from GSE39582 data set and The Cancer Genome Atlas data set (TCGA-COAD and TCGA-READ). Weighted gene co-expression network analysis (WGCNA) was conducted in GSE39582 data set to identify hub genes associated with neutrophil infiltration. The association of hub gene and neutrophils was then validated in TCGA cohorts and an independent RJ cohort. Functional analysis was performed to investigate the molecular mechanisms of the interested hub gene. We found that neutrophil infiltration is elevated in CRC tissues, and it is related to a poorer prognosis. A total of 18 gene modules are identified by WGCNA in GSE39582 data set, among which lightcyan module is significantly correlated with neutrophils infiltration. Furthermore, Superoxide Dismutase 2 (SOD2) in lightcyan module was proved to correlated with neutrophils infiltration in various cancer types. In addition, SOD2 expression is highly associated with several chemokines, including CXCL8, a neutrophils-related attractant, and functional analysis revealed that SOD2 is involved in neutrophils recruitment biological process. These results indicate that an 'SOD2-CXCL8-neutrophil recruitment' axis plays a potential role in colorectal cancer progression.