Identification of BIK as an unfavorable prognostic marker and novel therapeutic target in microsatellite stable colorectal cancer harboring KRAS mutations.

ABSTRACT

KRAS mutations lead to persistent activation of multiple downstream effectors that drive the cancer phenotype. Approximately 30%-50% of colorectal cancer (CRC) patients harbor KRAS mutations, which confer more aggressive tumor biology and shorter overall survival (OS), especially in microsatellite stable (MSS) metastatic CRC. Given that KRAS mutant protein has been proven difficult to target directly, identifying genes that function closely with KRAS and targeting these genes seems to be a promising therapeutic strategy for KRAS-mutated MSS CRC. Here, KRAS function-sensitive genes were identified by assessing the correlation between gene dependency scores from CRISPR knockout screens and KRAS mRNA expression in KRAS-mutated MSS CRC cell lines in the Cancer Cell Line Encyclopedia (CCLE) database. If the correlation coefficient was ≥ 0.6, the gene was considered a KRAS function-sensitive gene. Then KRAS function-sensitive genes related to prognosis were screened out in The Cancer Genome Atlas (TCGA) cohort, and the prognostic value was validated in the Gene Expression Omnibus (GEO) cohort. Single-sample gene set enrichment analysis (ssGSEA) was performed to investigate the potential mechanisms. PockDrug-Server was used to predict the druggability of candidate genes. The results showed that in 20 KRAS-mutated MSS CRC cell lines, 13 genes were identified as KRAS function-sensitive genes. Of these 13 genes, only BIK expression was significantly associated with progression-free survival (PFS) and OS, and the BIK-high patients had significantly poorer PFS (HR=3.18, P=0.020) and OS (HR=4.74, P=0.030) than the BIK-low patients. Multivariate Cox regression analysis revealed high BIK expression as an independent predictor for poorer prognosis in KRAS-mutated MSS CRC. The prognostic value of BIK was also successfully validated in a GEO cohort. The results of ssGSEA showed that the BIK-high group was more prone to strong metastasis activity than the BIK-low group. Pocket druggability prediction analysis presented that BIK had three druggable pockets, and their druggability scores were above 0.8. These findings suggested that BIK is a promising prognostic marker and therapeutic target in KRAS-mutated MSS CRC.