Multiregional transcriptomics identifies congruent consensus subtypes with prognostic value beyond tumor heterogeneity of colorectal cancer.

Intra-tumor heterogeneity compromises the clinical value of transcriptomic classifications of colorectal cancer. We investigated the prognostic effect of transcriptomic heterogeneity and the potential for classifications less vulnerable to heterogeneity in a single-hospital series of 1093 tumor samples from 692 patients, including multiregional samples from 98 primary tumors and 35 primary-metastasis sets. We show that intra-tumor heterogeneity of the consensus molecular subtypes (CMS) is frequent and has poor-prognostic associations independently of tumor microenvironment markers. Multiregional transcriptomics uncover cancer cell-intrinsic and low-heterogeneity signals that recapitulate the intrinsic CMSs proposed by single-cell sequencing. Further subclassification identifies congruent CMSs that explain a larger proportion of variation in patient survival than intra-tumor heterogeneity. Plasticity is indicated by discordant intrinsic phenotypes of matched primary and metastatic tumors. We conclude that multiregional sampling reconciles the prognostic power of tumor classifications from single-cell and bulk transcriptomics in the context of intra-tumor heterogeneity, and phenotypic plasticity challenges the reconciliation of primary and metastatic subtypes.

The prognostic effect of tumor-associated macrophages in stage I-III colorectal cancer depends on T cell infiltration.

BACKGROUND: Tumor-associated macrophages (TAMs) are associated with unfavorable patient prognosis in many cancer types. However, TAMs are a heterogeneous cell population and subsets have been shown to activate tumor-infiltrating T cells and confer a good patient prognosis. Data on the prognostic value of TAMs in colorectal cancer are conflicting. We investigated the prognostic effect of TAMs in relation to tumor-infiltrating T cells in colorectal cancers. METHODS: The TAM markers CD68 and CD163 were analyzed by multiplex fluorescence immunohistochemistry and digital image analysis on tissue microarrays of 1720 primary colorectal cancers. TAM density in the tumor stroma was scored in relation to T cell density (stromal CD3+ and epithelial CD8+ cells) and analyzed in Cox proportional hazards models of 5-year relapse-free survival. Multivariable survival models included clinicopathological factors, MSI status and BRAFV600E mutation status. RESULTS: High TAM density was associated with a favorable 5-year relapse-free survival in a multivariable model of patients with stage I-III tumors (p = 0.004, hazard ratio 0.94, 95% confidence interval 0.90-0.98). However, the prognostic effect was dependent on tumoral T-cell density. High TAM density was associated with a good prognosis in patients who also had high T-cell levels in their tumors, while high TAM density was associated with poorer prognosis in patients with low T-cell levels (pinteraction = 0.0006). This prognostic heterogeneity was found for microsatellite stable tumors separately. CONCLUSIONS: This study supported a phenotypic heterogeneity of TAMs in colorectal cancer, and showed that combined tumor immunophenotyping of multiple immune cell types improved the prediction of patient prognosis.

Evolutionary mode and timing of dissemination of high-grade serous carcinomas.

Dissemination within the peritoneal cavity is a main determinant of poor patient outcomes from high-grade serous carcinomas (HGSCs). The dissemination process is poorly understood from a cancer evolutionary perspective. We reconstructed the evolutionary trajectories across a median of 5 tumor sites and regions from each of 23 patients based on deep whole-exome sequencing. Polyclonal cancer origin was detected in 1 patient. Ovarian tumors had more complex subclonal architectures than other intraperitoneal tumors in each patient, which indicated that tumors developed earlier in the ovaries. Three common modes of dissemination were identified, including monoclonal or polyclonal dissemination of monophyletic (linear) or polyphyletic (branched) subclones. Mutation profiles of initial or disseminated clones varied greatly among cancers, but recurrent mutations were found in 7 cancer-critical genes, including TP53, BRCA1, BRCA2, and DNMT3A, and in the PI3K/AKT1 pathway. Disseminated clones developed late in the evolutionary trajectory models of most cancers, in particular in cancers with DNA damage repair deficiency. Polyclonal dissemination was predicted to occur predominantly as a single and rapid wave, but chemotherapy exposure was associated with higher genomic diversity of disseminated clones. In conclusion, we described three common evolutionary dissemination modes across HGSCs and proposed factors associated with dissemination diversity.

Transcriptomic subtyping of malignant peripheral nerve sheath tumours highlights immune signatures, genomic profiles, patient survival and therapeutic targets.

BACKGROUND: Malignant peripheral nerve sheath tumour (MPNST) is an aggressive orphan disease commonly affecting adolescents or young adults. Current knowledge of molecular tumour biology has been insufficient for development of rational treatment strategies. We aimed to discover molecular subtypes of potential clinical relevance. METHODS: Fresh frozen samples of MPNSTs (n = 94) and benign neurofibromas (n = 28) from 115 patients in a European multicentre study were analysed by DNA copy number and/or transcriptomic profiling. Unsupervised transcriptomic subtyping was performed and the subtypes characterized for genomic aberrations, clinicopathological associations and patient survival. FINDINGS: MPNSTs were classified into two transcriptomic subtypes defined primarily by immune signatures and proliferative processes. "Immune active" MPNSTs (44%) had sustained immune signals relative to neurofibromas, were more frequently low-grade (P = 0.01) and had favourable prognostic associations in a multivariable model of disease-specific survival with clinicopathological factors (hazard ratio 0.25, P = 0.003). "Immune deficient" MPNSTs were more aggressive and characterized by proliferative signatures, high genomic complexity, aberrant TP53 and PRC2 loss, as well as high relative expression of several potential actionable targets (EGFR, ERBB2, EZH2, KIF11, PLK1, RRM2). Integrated gene-wise analyses suggested a DNA copy number-basis for proliferative transcriptomic signatures in particular, and the tumour copy number burden further stratified the transcriptomic subtypes according to patient prognosis (P < 0.01). INTERPRETATION: Approximately half of MPNSTs belong to an "immune deficient" transcriptomic subtype associated with an aggressive disease course, PRC2 loss and expression of several potential therapeutic targets, providing a rationale for molecularly-guided intervention trials. FUNDING: Research grants from non-profit organizations, as stated in the Acknowledgements.

Spatial analysis and CD25-expression identify regulatory T cells as predictors of a poor prognosis in colorectal cancer.

Regulatory T cells (Tregs) are a heterogeneous cell population that can either suppress or stimulate immune responses. Tumor-infiltrating Tregs are associated with an adverse outcome from most cancer types, but have generally been found to be associated with a good prognosis in colorectal cancer (CRC). We investigated the prognostic heterogeneity of Tregs in CRC by co-expression patterns and spatial analyses with diverse T cell markers, using multiplex fluorescence immunohistochemistry and digital image analysis in two consecutive series of primary CRCs (total n = 1720). Treg infiltration in tumors, scored as FOXP3+ or CD4+/CD25+/FOXP3+ (triple-positive) cells, was strongly correlated to the overall amount of CD3+ and CD8+ T cells, and consequently associated with a favorable 5-year relapse-free survival rate among patients with stage I-III CRC who underwent complete tumor resection. However, high relative expression of the activation marker CD25 in triple-positive Tregs was independently associated with an adverse outcome in a multivariable model incorporating clinicopathological and known molecular prognostic markers (hazard ratio = 1.35, p = 0.028). Furthermore, spatial marker analysis based on Voronoi diagrams and permutation testing of cellular neighborhoods revealed a statistically significant proximity between Tregs and CD8+-cells in 18% of patients, and this was independently associated with a poor survival (multivariable hazard ratio = 1.36, p = 0.017). These results show prognostic heterogeneity of different Treg populations in primary CRC, and highlight the importance of multi-marker and spatial analyses for accurate immunophenotyping of tumors in relation to patient outcome.

Expressed prognostic biomarkers for primary prostate cancer independent of multifocality and transcriptome heterogeneity.

The majority of prostate cancer patients are diagnosed with multiple primary malignant foci. The distinct foci are exceptionally heterogeneous with regard to DNA mutations, but whether this is recapitulated at the transcriptome level remains unknown. In this study, inter- and intrafocal heterogeneity has been assessed by whole-transcriptome sequencing of 87 tissue samples from 23 patients with localized prostate cancer treated with radical prostatectomy. From each patient, multiple samples were taken from one or more malignant foci, in addition to one sample from benign prostate tissue. Transcriptomic profiles of different malignant foci from the same patient showed a similar level of heterogeneity as tumors from different patients. This applies to expression of genes, fusion genes, and somatic mutations. Within-patient pair-wise analyses identified expression patterns linked to ETS status and extraprostatic extension. A set of 62 genes were found with low intrapatient heterogeneity and high interpatient heterogeneity, retaining stable expression profiles across foci within the same patient. Among these, 16 genes are associated with biochemical recurrence in a separately published study and are therefore nominated as biomarkers with prognostic value regardless of which malignant focus is sampled. In conclusion, an extensive heterogeneity in multifocal prostate cancer is confirmed at the gene expression level. Diagnostic biomarkers were identified for ETS positive samples and samples from extraprostatic extensions. Finally, prognostic biomarkers independent of multifocal heterogeneity were found.

Increased sensitivity to SMAC mimetic LCL161 identified by longitudinal ex vivo pharmacogenomics of recurrent, KRAS mutated rectal cancer liver metastases.

Tumor heterogeneity is a primary cause of treatment failure. However, changes in drug sensitivity over time are not well mapped in cancer. Patient-derived organoids (PDOs) may predict clinical drug responses ex vivo and offer an opportunity to evaluate novel treatment strategies in a personalized fashion. Here we have evaluated spatio-temporal functional and molecular dynamics of five PDO models established after hepatic re-resections and neoadjuvant combination chemotherapies in a patient with microsatellite stable and KRAS mutated metastatic rectal cancer. Histopathological differentiation phenotypes of the PDOs corresponded with the liver metastases, and ex vivo drug sensitivities generally reflected clinical responses and selection pressure, assessed in comparison to a reference data set of PDOs from metastatic colorectal cancers. PDOs from the initial versus the two recurrent metastatic settings showed heterogeneous cell morphologies, protein marker expression, and drug sensitivities. Exploratory analyses of a drug screen library of 33 investigational anticancer agents showed the strongest ex vivo sensitivity to the SMAC mimetic LCL161 in PDOs of recurrent disease compared to those of the initial metastasis. Functional analyses confirmed target inhibition and apoptosis induction in the LCL161 sensitive PDOs from the recurrent metastases. Gene expression analyses indicated an association between LCL161 sensitivity and tumor necrosis factor alpha signaling and RIPK1 gene expression. In conclusion, LCL161 was identified as a possible experimental therapy of a metastatic rectal cancer that relapsed after hepatic resection and standard systemic treatment.

The expressed mutational landscape of microsatellite stable colorectal cancers.

BACKGROUND: Colorectal cancer is the 2nd leading cause of cancer-related deaths with few patients benefiting from biomarker-guided therapy. Mutation expression is essential for accurate interpretation of mutations as biomarkers, but surprisingly, little has been done to analyze somatic cancer mutations on the expression level. We report a large-scale analysis of allele-specific mutation expression. METHODS: Whole-exome and total RNA sequencing was performed on 137 samples from 121 microsatellite stable colorectal cancers, including multiregional samples of primary and metastatic tumors from 4 patients. Data were integrated with allele-specific resolution. Results were validated in an independent set of 241 colon cancers. Therapeutic associations were explored by pharmacogenomic profiling of 15 cell lines or patient-derived organoids. RESULTS: The median proportion of expressed mutations per tumor was 34%. Cancer-critical mutations had the highest expression frequency (gene-wise mean of 58%), independent of frequent allelic imbalance. Systematic deviation from the general pattern of expression levels according to allelic frequencies was detected, including preferential expression of mutated alleles dependent on the mutation type and target gene. Translational relevance was suggested by correlations of KRAS/NRAS or TP53 mutation expression levels with downstream oncogenic signatures (p < 0.03), overall survival among patients with stage II and III cancer (KRAS/NRAS: hazard ratio 6.1, p = 0.0070), and targeted drug sensitivity. The latter was demonstrated for EGFR and MDM2 inhibition in pre-clinical models. CONCLUSIONS: Only a subset of mutations in microsatellite stable colorectal cancers were expressed, and the "expressed mutation dose" may provide an opportunity for more fine-tuned biomarker interpretations.

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.

Treatment outcomes and prognostic factors after chemoradiotherapy for anal cancer.

BACKGROUND: Squamous cell carcinoma of the anus (SCCA) is a rare malignancy with rising incidence, associated with human papilloma virus (HPV). Chemoradiotherapy (CRT) is the preferred treatment. The purpose was to investigate treatment failure, survival and prognostic factors after CRT. MATERIAL AND METHODS: In this prospective observational study from a large regional centre, 141 patients were included from 2013 to 2017, and 132 were eligible for analysis. The main inclusion criteria were SCCA, planned radiotherapy, and performance status (ECOG) ≤2. Patient characteristics, disease stage, treatment, and treatment response were prospectively registered. Disease-free survival (DFS), overall survival (OS), and locoregional treatment failure after CRT were analysed. Hazard ratios (HRs) were estimated with Cox`s proportional hazards model. RESULTS: Median follow-up was 54 (range 6-71) months. Eighteen patients (14%) had treatment failures after CRT; of these 10 (8%) had residual tumour, and 8 (6%) relapse as first failure. The first treatment failure was locoregional (11 patients), distant (5 patients), and both (2 patients). Salvage abdomino-perineal resection was performed in 10 patients, 2 had resections of metastases, and 3 both. DFS was 85% at 3 years and 78% at 5 years. OS was 93% at 3 years and 86% at 5 years. In analyses adjusted for age and gender, HPV negative tumours (HR 2.5, p = 0.024), N3 disease (HR 2.6, p = 0.024), and tumour size ≥4 cm (HR 2.4, p = 0.038) were negative prognostic factors for DFS. CONCLUSION: State-of-the-art chemoradiotherapy for SCCA resulted in excellent outcomes, and improved survival compared with previous national data, with <15% treatment failures and a 3-year DFS of >80%.

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.

Multiregional assessment of CIMP in primary colorectal cancers: Phenotype concordance but marker variability.

Intratumor heterogeneity of colorectal cancers (CRCs) is manifested both at the genomic and epigenomic levels. Early genetic aberrations in carcinogenesis are clonal and present throughout the tumors, but less is known about the heterogeneity of the epigenetic CpG island methylator phenotype (CIMP). CIMP characterizes a subgroup of CRCs thought to originate from specific precursor lesions, and it is defined by widespread DNA methylation within promoter regions. In this work, we investigated CIMP in two to four multiregional samples from 30 primary tumors (n = 86 samples) using the consensus Weisenberger gene panel (CACNA1G, IGF2, NEUROG1, RUNX3 and SOCS1). Twenty-nine of 30 tumors (97%) showed concordant CIMP status in all samples, and percent methylated reference (PMR) values of all five markers had higher intertumor than intratumor variation (P value = 1.5e-09). However, a third of the CIMP+ tumors exhibited discrepancies in methylation status in at least one of the five gene markers. To conclude, CIMP status was consistent within primary CRCs, and it is likely a clonal phenotype. However, spatial discordances of the individual genes suggest that large-scale analysis of multiregional samples could be of interest for identifying CIMP markers that are robust to intratumor heterogeneity.

Prediction of relapse-free survival according to adjuvant chemotherapy and regulator of chromosome condensation 2 (RCC2) expression in colorectal cancer.

BACKGROUND: There is a need for improved selection of patients for adjuvant chemotherapy after resection of non-metastatic colorectal cancer (CRC). Regulator of chromosome condensation 2 (RCC2) is a potential prognostic biomarker. We report on the establishment of a robust protocol for RCC2 expression analysis and prognostic tumour biomarker evaluation in patients who did and did not receive adjuvant chemotherapy. MATERIALS AND METHODS: RCC2 was analysed in 2916 primary CRCs from the QUASAR2 randomised trial and two single-hospital Norwegian series. A new protocol using fluorescent antibody staining and digital image analysis was optimised. Biomarker value for 5-year relapse-free survival was analysed in relation to tumour stage, adjuvant chemotherapy and the molecular markers microsatellite instability, KRAS/BRAFV600E/TP53 mutations and CDX2 expression. RESULTS: Low RCC2 expression was scored in 41% of 2696 evaluable samples. Among patients with stage I-III CRC who had not received adjuvant chemotherapy, low RCC2 expression was an independent marker of inferior 5-year relapse-free survival in multivariable Cox models including clinicopathological factors and molecular markers (HR 1.45, 95% CI 1.09 to 1.94, p=0.012, N=521). RCC2 was not prognostic in patients who had received adjuvant chemotherapy, neither in QUASAR2 nor the pooled Norwegian series. The interaction between RCC2 and adjuvant chemotherapy for prediction of patient outcome was significant in stage III, and strongest among patients with microsatellite stable tumours (pinteraction=0.028). CONCLUSIONS: Low expression of RCC2 is a biomarker for poor prognosis in patients with stage I-III CRC and seems to be a predictive biomarker for effect of adjuvant chemotherapy.

Molecular correlates of sensitivity to PARP inhibition beyond homologous recombination deficiency in pre-clinical models of colorectal cancer point to wild-type TP53 activity.

BACKGROUND: PARP inhibitors are active in various tumour types beyond BRCA-mutant cancers, but their activity and molecular correlates in colorectal cancer (CRC) are not well studied. METHODS: Mutations and genome-wide mutational patterns associated with homologous recombination deficiency (HRD) were investigated in 255 primary CRCs with whole-exome sequencing and/or DNA copy number data. Efficacy of five PARP inhibitors and their molecular correlates were evaluated in 93 CRC cell lines partly annotated with mutational-, DNA copy number-, and/or gene expression profiles. Post-treatment gene expression profiling and specific protein expression analyses were performed in two pairs of PARP inhibitor sensitive and resistant cell lines. FINDINGS: A subset of microsatellite stable (MSS) CRCs had truncating mutations in homologous recombination-related genes, but these were not associated with genomic signatures of HRD. Eight CRC cell lines (9%) were sensitive to PARP inhibition, but sensitivity was not predicted by HRD-related genomic and transcriptomic signatures. In contrast, drug sensitivity in MSS cell lines was strongly associated with TP53 wild-type status (odds ratio 15.7, p = 0.023) and TP53-related expression signatures. Increased downstream TP53 activity was among the primary response mechanisms, and TP53 inhibition antagonized the effect of PARP inhibitors. Wild-type TP53-mediated suppression of RAD51 was identified as a possible mechanism of action for sensitivity to PARP inhibition. INTERPRETATION: PARP inhibitors are active in a subset of CRC cell lines and preserved TP53 function may increase the likelihood of response.

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.

Improved prognostication of glioblastoma beyond molecular subtyping by transcriptional profiling of the tumor microenvironment.

Glioblastoma (GBM), the most aggressive form of brain cancer, is characterized by a high level of molecular heterogeneity, and infiltration by various immune and stromal cell populations. Important advances have been made in deciphering the microenvironment of GBMs, but its association with existing molecular subtypes and its potential prognostic role remain elusive. We have investigated the abundance of infiltrating immune and stromal cells in silico, from gene expression profiles. Two cohorts, including in-house normal brain and glioma samples (n = 70) and a large sample set from TCGA (n = 393), were combined into a single exploratory dataset. A third independent cohort (n = 124) was used for validation. Tumors were clustered based on their microenvironment infiltration profiles, and associations with known GBM molecular subtypes and patient outcome were tested a posteriori in a multivariable setting. We identified a subset of GBM samples with significantly higher abundances of most immune and stromal cell populations. This subset showed increased expression of both immune suppressor and immune effector genes compared to other GBMs and was enriched for the mesenchymal molecular subtype. Survival analyses suggested that tumor microenvironment infiltration pattern was an independent prognostic factor for GBM patients. Among all, patients with the mesenchymal subtype with low immune and stromal infiltration had the poorest survival. By combining molecular subtyping with gene expression measures of tumor infiltration, the present work contributes with improving prognostic models in GBM.

High Concordance and Negative Prognostic Impact of RAS/BRAF/PIK3CA Mutations in Multiple Resected Colorectal Liver Metastases.

BACKGROUND: The prevalence and clinical implications of genetic heterogeneity in patients with multiple colorectal liver metastases remain largely unknown. In a prospective series of patients undergoing resection of colorectal liver metastases, the aim was to investigate the inter-metastatic and primary-to-metastatic heterogeneity of mutations in KRAS, NRAS, BRAF, and PIK3CA and their prognostic impact. PATIENTS AND METHODS: We analyzed the mutation status among 372 liver metastases and 78 primary tumors from 106 patients by methods used in clinical routine testing, by Sanger sequencing, by next-generation sequencing (NGS), and/or by droplet digital polymerase chain reaction. The 3-year cancer-specific survival (CSS) was analyzed using the Kaplan-Meier method. RESULTS: Although Sanger sequencing indicated inter-metastatic mutation heterogeneity in 14 of 97 patients (14%), almost all cases were refuted by high-sensitive NGS. Also, heterogeneity among metastatic deposits was concluded only for PIK3CA in 2 patients. Similarly, primary-to-metastatic heterogeneity was indicated in 8 of 78 patients (10%) using Sanger sequencing but for only 2 patients after NGS, showing the emergence of 1 KRAS and 1 PIK3CA mutation in the metastatic lesions. KRAS mutations were present in 53 of 106 patients (50%) and were associated with poorer 3-year CSS after liver resection (37% vs. 61% for KRAS wild-type; P = .004). Poor prognostic associations were found also for the combination of KRAS/NRAS/BRAF mutations compared with triple wild-type (P = .002). CONCLUSION: Intra-patient mutation heterogeneity was virtually undetected, both between the primary tumor and the liver metastases and among the metastatic deposits. KRAS mutations separately, and KRAS/NRAS/BRAF mutations combined, were associated with poor patient survival after partial liver resection.