Relative contribution of clinicopathological variables, genomic markers, transcriptomic subtyping and microenvironment features for outcome prediction in stage II/III colorectal cancer

It remains unknown to what extent consensus molecular subtype (CMS) groups and immune-stromal infiltration patterns improve our ability to predict outcomes over tumor–node–metastasis (TNM) staging and microsatellite instability (MSI) status in early-stage colorectal cancer (CRC). We carried out a comprehensive retrospective biomarker analysis of prognostic markers in adjuvant chemotherapy-untreated (N = 1656) and treated (N = 980), stage II (N = 1799) and III (N = 837) CRCs. We defined CMS scores and estimated CD8+ cytotoxic lymphocytes (CytoLym) and cancer-associated fibroblasts (CAF) infiltration scores from bulk tumor tissue transcriptomes (CMSclassifier and MCPcounter R packages); constructed a stratified multivariable Cox model for disease-free survival (DFS); and calculated the relative proportion of explained variation by each marker (clinicopathological [ClinPath], genomics [Gen: MSI, BRAF and KRAS mutations], CMS scores [CMS] and microenvironment cells [MicroCells: CytoLym+CAF]). In multivariable models, only ClinPath and MicroCells remained significant prognostic factors, with both CytoLym and CAF infiltration scores improving survival prediction beyond other markers. The explained variation for DFS models of ClinPath, MicroCells, Gen markers and CMS4 scores was 77%, 14%, 5.3% and 3.7%, respectively, in stage II; and 55.9%, 35.1%, 4.1% and 0.9%, respectively, in stage III. Patients whose tumors were CytoLym high/CAF low had better DFS than other strata [HR=0.71 (0.6–0.9); P = 0.004]. Microsatellite stable tumors had the strongest signal for improved outcomes with CytoLym high scores (interaction P = 0.04) and the poor prognosis linked to high CAF scores was limited to stage III disease (interaction P = 0.04). Our results confirm that tumor microenvironment infiltration patterns represent potent determinants of the risk for distant dissemination in early-stage CRC. Multivariable models suggest that the prognostic value of MSI and CMS groups is largely explained by CytoLym and CAF infiltration patterns.