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

Abstract

Background: 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).

Patients and methods: 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]).

Results: 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).

Conclusions: 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.

Keywords: CMS; colorectal cancer; immune; microsatellite instability; prognosis; stromal.

Figure 1.

Adjusted disease-free survival (DFS) Kaplan−Meier curves stratified by CytoLym and CAF groups, estimated with multivariable Cox proportional hazard model controlling for pT, pN, age, sex, primary tumor location, MSI status and CMS4 scores:chemotherapy untreated (A), chemotherapy treated (B), stage II (C) and stage III (D) CRC cohorts.

Figure 2.

Relative proportion of explained variation in disease-free survival (DFS) of the full multivariable model in chemotherapy untreated (A), chemotherapy treated (B), stage II (C) and stage III (D) CRC cohorts.