Identification of key genes involved in tumor immune cell infiltration and cetuximab resistance in colorectal cancer

BACKGROUND: The anti-epidermal growth factor receptor (EGFR) antibody introduces adaptable variations to the transcriptome and triggers tumor immune infiltration, resulting in colorectal cancer (CRC) treatment resistance. We intended to identify genes that play essential roles in cetuximab resistance and tumor immune cell infiltration. METHODS: A cetuximab-resistant CACO2 cellular model was established, and its transcriptome variations were detected by microarray. Meanwhile, public data from the Gene Expression Omnibus and The Cancer Genome Atlas (TCGA) database were downloaded. Integrated bioinformatics analysis was applied to detect differentially expressed genes (DEGs) between the cetuximab-resistant and the cetuximab-sensitive groups. Then, we investigated correlations between DEGs and immune cell infiltration. The DEGs from bioinformatics analysis were further validated in vitro and in clinical samples. RESULTS: We identified 732 upregulated and 1259 downregulated DEGs in the induced cellular model. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analyses, along with Gene Set Enrichment Analysis and Gene Set Variation Analysis, indicated the functions of the DEGs. Together with GSE59857 and GSE5841, 12 common DEGs (SATB-2, AKR1B10, ADH1A, ADH1C, MYB, ATP10B, CDX-2, FAR2, EPHB2, SLC26A3, ORP-1, VAV3) were identified and their predictive values of cetuximab treatment were validated in GSE56386. In online Genomics of Drug Sensitivity in Cancer (GDSC) database, nine of twelve DEGs were recognized in the protein-protein (PPI) network. Based on the transcriptome profiles of CRC samples in TCGA and using Tumor Immune Estimation Resource Version 2.0, we bioinformatically determined that SATB-2, ORP-1, MYB, and CDX-2 expressions were associated with intensive infiltration of B cell, CD4(+) T cell, CD8(+) T cell and macrophage, which was then validated the correlation in clinical samples by immunohistochemistry. We found that SATB-2, ORP-1, MYB, and CDX-2 were downregulated in vitro with cetuximab treatment. Clinically, patients with advanced CRC and high ORP-1 expression exhibited a longer progression-free survival time when they were treated with anti-EGFR therapy than those with low ORP-1 expression. CONCLUSIONS: SATB-2, ORP-1, MYB, and CDX-2 were related to cetuximab sensitivity as well as enhanced tumor immune cell infiltration in patients with CRC.