Abstract
KRAS mutation is the most frequent oncogenic aberration in colorectal cancer (CRC). The molecular mechanism and clinical implications of KRAS mutation in CRC remain unclear and show high heterogeneity within these tumors. Methods: We harnessed the multi-omics data (genomic, transcriptomic, proteomic, and phosphoproteomic etc.) of KRAS-mutant CRC tumors and performed unsupervised clustering to identify proteomics-based subgroups and molecular characterization. Results: In-depth analysis of the tumor microenvironment by single-cell transcriptomic revealed the cellular landscape of KRAS-mutant CRC tumors and identified the specific cell subsets with KRAS mutation. Integrated multi-omics analyses separated the KRAS-mutant tumors into two distinct molecular subtypes, termed KRAS-M1 (KM1) and KRAS-M2 (KM2). The two subtypes had a similar distribution of mutated residues in KRAS (G12D/V/C etc.) but were characterized by distinct features in terms of prognosis, genetic alterations, microenvironment dysregulation, biological phenotype, and potential therapeutic approaches. Proteogenomic analyses revealed that the EMT, TGF-β and angiogenesis pathways were enriched in the KM2 subtype and that the KM2 subtype was associated with the mesenchymal phenotype-related CMS4 subtype, which indicated stromal invasion and worse prognosis. The KM1 subtype was characterized predominantly by activation of the cell cycle, E2F and RNA transcription and was associated with the chromosomal instability (CIN)-related ProS-E proteomic subtype, which suggested cyclin-dependent features and better survival outcomes. Moreover, drug sensitivity analyses based on three compound databases revealed subgroup-specific agents for KM1 and KM2 tumors. Conclusions: This study clarifies the molecular heterogeneity of KRAS-mutant CRC and reveals new biological subtypes and therapeutic possibilities for these tumors.