Pathway level subtyping identifies a slow-cycling biological phenotype associated with poor clinical outcomes in colorectal cancer.

Molecular stratification using gene-level transcriptional data has identified subtypes with distinctive genotypic and phenotypic traits, as exemplified by the consensus molecular subtypes (CMS) in colorectal cancer (CRC). Here, rather than gene-level data, we make use of gene ontology and biological activation state information for initial molecular class discovery. In doing so, we defined three pathway-derived subtypes (PDS) in CRC: PDS1 tumors, which are canonical/LGR5+ stem-rich, highly proliferative and display good prognosis; PDS2 tumors, which are regenerative/ANXA1+ stem-rich, with elevated stromal and immune tumor microenvironmental lineages; and PDS3 tumors, which represent a previously overlooked slow-cycling subset of tumors within CMS2 with reduced stem populations and increased differentiated lineages, particularly enterocytes and enteroendocrine cells, yet display the worst prognosis in locally advanced disease. These PDS3 phenotypic traits are evident across numerous bulk and single-cell datasets, and demark a series of subtle biological states that are currently under-represented in pre-clinical models and are not identified using existing subtyping classifiers.

KRAS allelic imbalance drives tumour initiation yet suppresses metastasis in colorectal cancer in vivo.

Oncogenic KRAS mutations are well-described functionally and are known to drive tumorigenesis. Recent reports describe a significant prevalence of KRAS allelic imbalances or gene dosage changes in human cancers, including loss of the wild-type allele in KRAS mutant cancers. However, the role of wild-type KRAS in tumorigenesis and therapeutic response remains elusive. We report an in vivo murine model of colorectal cancer featuring deletion of wild-type Kras in the context of oncogenic Kras. Deletion of wild-type Kras exacerbates oncogenic KRAS signalling through MAPK and thus drives tumour initiation. Absence of wild-type Kras potentiates the oncogenic effect of KRASG12D, while incidentally inducing sensitivity to inhibition of MEK1/2. Importantly, loss of the wild-type allele in aggressive models of KRASG12D-driven CRC significantly alters tumour progression, and suppresses metastasis through modulation of the immune microenvironment. This study highlights the critical role for wild-type Kras upon tumour initiation, progression and therapeutic response in Kras mutant CRC.