De novo transcriptomic subtyping of colorectal cancer liver metastases in the context of tumor heterogeneity.

BACKGROUND: Gene expression-based subtyping has the potential to form a new paradigm for stratified treatment of colorectal cancer. However, current frameworks are based on the transcriptomic profiles of primary tumors, and metastatic heterogeneity is a challenge. Here we aimed to develop a de novo metastasis-oriented framework. METHODS: In total, 829 transcriptomic profiles from patients with colorectal cancer were analyzed, including primary tumors, liver metastases, and non-malignant liver samples. High-resolution microarray gene expression profiling was performed of 283 liver metastases from 171 patients treated by hepatic resection, including multiregional and/or multi-metastatic samples from each of 47 patients. A single randomly selected liver metastasis sample from each patient was used for unsupervised subtype discovery by nonnegative matrix factorization, and a random forest prediction model was trained to classify multi-metastatic samples, as well as liver metastases from two independent series of 308 additional patients. RESULTS: Initial comparisons with non-malignant liver samples and primary colorectal tumors showed a highly variable degree of influence from the liver microenvironment in metastases, which contributed to inter-metastatic transcriptomic heterogeneity, but did not define subtype distinctions. The de novo liver metastasis subtype (LMS) framework recapitulated the main distinction between epithelial-like and mesenchymal-like tumors, with a strong immune and stromal component only in the latter. We also identified biologically distinct epithelial-like subtypes originating from different progenitor cell types. LMS1 metastases had several transcriptomic features of cancer aggressiveness, including secretory progenitor cell origin, oncogenic addictions, and microsatellite instability in a microsatellite stable background, as well as frequent RAS/TP53 co-mutations. The poor-prognostic association of LMS1 metastases was independent of mutation status, clinicopathological variables, and current subtyping frameworks (consensus molecular subtypes and colorectal cancer intrinsic subtypes). LMS1 was also the least heterogeneous subtype in comparisons of multiple metastases per patient, and tumor heterogeneity did not confound the prognostic value of LMS1. CONCLUSIONS: We report the first large study of multi-metastatic gene expression profiling of colorectal cancer. The new metastasis-oriented subtyping framework showed potential for clinically relevant transcriptomic classification in the context of metastatic heterogeneity, and an LMS1 mini-classifier was constructed to facilitate prognostic stratification and further clinical testing.

Metastatic heterogeneity of the consensus molecular subtypes of colorectal cancer.

Gene expression-based subtypes of colorectal cancer have clinical relevance, but the representativeness of primary tumors and the consensus molecular subtypes (CMS) for metastatic cancers is not well known. We investigated the metastatic heterogeneity of CMS. The best approach to subtype translation was delineated by comparisons of transcriptomic profiles from 317 primary tumors and 295 liver metastases, including multi-metastatic samples from 45 patients and 14 primary-metastasis sets. Associations were validated in an external data set (n = 618). Projection of metastases onto principal components of primary tumors showed that metastases were depleted of CMS1-immune/CMS3-metabolic signals, enriched for CMS4-mesenchymal/stromal signals, and heavily influenced by the microenvironment. The tailored CMS classifier (available in an updated version of the R package CMScaller) therefore implemented an approach to regress out the liver tissue background. The majority of classified metastases were either CMS2 or CMS4. Nonetheless, subtype switching and inter-metastatic CMS heterogeneity were frequent and increased with sampling intensity. Poor-prognostic value of CMS1/3 metastases was consistent in the context of intra-patient tumor heterogeneity.

Genomic and prognostic heterogeneity among RAS/BRAFV600E /TP53 co-mutated resectable colorectal liver metastases.

Hepatic resection is potentially curative for patients with colorectal liver metastases, but the treatment benefit varies. KRAS/NRAS (RAS)/TP53 co-mutations are associated with a poor prognosis after resection, but there is large variation in patient outcome within the mutation groups, and genetic testing is currently not used to evaluate benefit from surgery. We have investigated the potential for improved prognostic stratification by combined biomarker analysis with DNA copy number aberrations (CNAs), and taking tumor heterogeneity into account. We determined the mutation status of RAS, BRAFV600 , and TP53 in 441 liver lesions from 171 patients treated by partial hepatectomy for metastatic colorectal cancer. CNAs were profiled in 232 tumors from 67 of the patients. Mutations and high-level amplifications of cancer-critical genes, the latter including ERBB2 and EGFR, were predominantly homogeneous within patients. RAS/BRAFV600E and TP53 co-mutations were associated with a poor patient outcome (hazard ratio, HR, 3.9, 95% confidence interval, CI, 1.3-11.1, P = 0.012) in multivariable analyses with clinicopathological variables. The genome-wide CNA burden and intrapatient intermetastatic CNA heterogeneity varied within the mutation groups, and the CNA burden had prognostic associations in univariable analysis. Combined prognostic analyses of RAS/BRAFV600E /TP53 mutations and CNAs, either as a high CNA burden or high intermetastatic CNA heterogeneity, identified patients with a particularly poor outcome (co-mutation/high CNA burden: HR 2.7, 95% CI 1.2-5.9, P = 0.013; co-mutation/high CNA heterogeneity: HR 2.5, 95% CI 1.1-5.6, P = 0.022). In conclusion, DNA copy number profiling identified genomic and prognostic heterogeneity among patients with resectable colorectal liver metastases with co-mutated RAS/BRAFV600E /TP53.

Patient-Derived Organoids from Multiple Colorectal Cancer Liver Metastases Reveal Moderate Intra-patient Pharmacotranscriptomic Heterogeneity.

PURPOSE: Molecular tumor heterogeneity may have important implications for the efficacy of targeted therapies in metastatic cancers. Inter-metastatic heterogeneity of sensitivity to anticancer agents has not been well explored in colorectal cancer. EXPERIMENTAL DESIGN: We established a platform for ex vivo pharmacogenomic profiling of patient-derived organoids (PDO) from resected colorectal cancer liver metastases. Drug sensitivity testing (n = 40 clinically relevant agents) and gene expression profiling were performed on 39 metastases from 22 patients. RESULTS: Three drug-response clusters were identified among the colorectal cancer metastases, based primarily on sensitivities to EGFR and/or MDM2 inhibition, and corresponding with RAS mutations and TP53 activity. Potentially effective therapies, including off-label use of drugs approved for other cancer types, could be nominated for eighteen patients (82%). Antimetabolites and targeted agents lacking a decisive genomic marker had stronger differential activity than most approved chemotherapies. We found limited intra-patient drug sensitivity heterogeneity between PDOs from multiple (2-5) liver metastases from each of ten patients. This was recapitulated at the gene expression level, with a highly proportional degree of transcriptomic and pharmacological variation. One PDO with a multi-drug resistance profile, including resistance to EGFR inhibition in a RAS-mutant background, showed sensitivity to MEK plus mTOR/AKT inhibition, corresponding with low-level PTEN expression. CONCLUSIONS: Intra-patient inter-metastatic pharmacological heterogeneity was not pronounced and ex vivo drug screening may identify novel treatment options for metastatic colorectal cancer. Variation in drug sensitivities was reflected at the transcriptomic level, suggesting potential to develop gene expression-based predictive signatures to guide experimental therapies.