The consensus molecular subtype (CMS) classification of colorectal cancer is based on bulk transcriptomics. The underlying epithelial cell diversity remains unclear. We analyzed 373,058 single-cell transcriptomes from 63 patients, focusing on 49,155 epithelial cells. We identified a pervasive genetic and transcriptomic dichotomy of malignant cells, based on distinct gene expression, DNA copy number and gene regulatory network. We recapitulated these subtypes in bulk transcriptomes from 3,614 patients. The two intrinsic subtypes, iCMS2 and iCMS3, refine CMS. iCMS3 comprises microsatellite unstable (MSI-H) cancers and one-third of microsatellite-stable (MSS) tumors. iCMS3 MSS cancers are transcriptomically more similar to MSI-H cancers than to other MSS cancers. CMS4 cancers had either iCMS2 or iCMS3 epithelium; the latter had the worst prognosis. We defined the intrinsic epithelial axis of colorectal cancer and propose a refined ‘IMF’ classification with five subtypes, combining intrinsic epithelial subtype (I), microsatellite instability status (M) and fibrosis (F).

Chemo-resistance is one of the major causes of cancer-related deaths. Here we used single-cell transcriptomics to investigate divergent modes of chemo-resistance in tumor cells. We observed that higher degree of phenotypic intra-tumor heterogeneity (ITH) favors selection of pre-existing drug-resistant cells, whereas phenotypically homogeneous cells engage covert epigenetic mechanisms to trans-differentiate under drug-selection. This adaptation was driven by selection-induced gain of H3K27ac marks on bivalently poised resistance-associated chromatin, and therefore not expressed in the treatment-naïve setting. Mechanistic interrogation of this phenomenon revealed that drug-induced adaptation was acquired upon the loss of stem factor SOX2, and a concomitant gain of SOX9. Strikingly we observed an enrichment of SOX9 at drug-induced H3K27ac sites, suggesting that tumor evolution could be driven by stem cell-switch-mediated epigenetic plasticity. Importantly, JQ1 mediated inhibition of BRD4 could reverse drug-induced adaptation. These results provide mechanistic insights into the modes of therapy-induced cellular plasticity and underscore the use of epigenetic inhibitors in targeting tumor evolution.

Genomics-driven cancer therapeutics has gained prominence in personalized cancer treatment. However, its utility in indications lacking biomarker-driven treatment strategies remains limited. Here we present a “phenotype-driven precision-oncology” approach, based on the notion that biological response to perturbations, chemical or genetic, in ex vivo patient-individualized models can serve as predictive biomarkers for therapeutic response in the clinic. We generated a library of “screenable” patient-derived primary cultures (PDCs) for head and neck squamous cell carcinomas that reproducibly predicted treatment response in matched patient-derived-xenograft models. Importantly, PDCs could guide clinical practice and predict tumour progression in two n = 1 co-clinical trials. Comprehensive “-omics” interrogation of PDCs derived from one of these models revealed YAP1 as a putative biomarker for treatment response and survival in ~24% of oral squamous cell carcinoma. We envision that scaling of the proposed PDC approach could uncover biomarkers for therapeutic stratification and guide real-time therapeutic decisions in the future.

While understanding molecular heterogeneity across patients underpins precision oncology, there is increasing appreciation for taking intra-tumor heterogeneity into account. Based on large-scale analysis of cancer omics datasets, we highlight the importance of intra-tumor transcriptomic heterogeneity (ITTH) for predicting clinical outcomes. Leveraging single-cell RNA-seq (scRNA-seq) with a recommender system (CaDRReS-Sc), we show that heterogeneous gene-expression signatures can predict drug response with high accuracy (80%). Using patient-proximal cell lines, we established the validity of CaDRReS-Sc’s monotherapy (Pearson r>0.6) and combinatorial predictions targeting clone-specific vulnerabilities (>10% improvement). Applying CaDRReS-Sc to rapidly expanding scRNA-seq compendiums can serve as in silico screen to accelerate drug-repurposing studies. Availability: https://github.com/CSB5/CaDRReS-Sc.

Inter-patient and intra-tumour heterogeneity (ITH) have prompted the need for a more personalised approach to cancer therapy. Although patient-derived xenograft (PDX) models can generate drug response specific to patients, they are not sustainable in terms of cost and time and have limited scalability. Tumour Organ-on-Chip (OoC) models are in vitro alternatives that can recapitulate some aspects of the 3D tumour microenvironment and can be scaled up for drug screening. While many tumour OoC systems have been developed to date, there have been limited validation studies to ascertain whether drug responses obtained from tumour OoCs are comparable to those predicted from patient-derived xenograft (PDX) models. In this study, we established a multiplexed tumour OoC device, that consists of an 8 × 4 array (32-plex) of culture chamber coupled to a concentration gradient generator. The device enabled perfusion culture of primary PDX-derived tumour spheroids to obtain dose-dependent response of 5 distinct standard-of-care (SOC) chemotherapeutic drugs for 3 colorectal cancer (CRC) patients. The in vitro efficacies of the chemotherapeutic drugs were rank-ordered for individual patients and compared to the in vivo efficacy obtained from matched PDX models. We show that quantitative correlation analysis between the drug efficacies predicted via the microfluidic perfusion culture is predictive of response in animal PDX models. This is a first study showing a comparative framework to quantitatively correlate the drug response predictions made by a microfluidic tumour organ-on-chip (OoC) model with that of PDX animal models.

The Eph family of receptor tyrosine kinases is crucial for assembly and maintenance of healthy tissues. Dysfunction in Eph signaling is causally associated with cancer progression. In breast cancer cells, dysregulated Eph signaling has been linked to alterations in receptor clustering abilities. Here, we implemented a single-cell assay and a scoring scheme to systematically probe the spatial organization of activated EphA receptors in multiple carcinoma cells. We show that cancer cells retain EphA clustering phenotype over several generations, and the degree of clustering reported for migration potential both at population and single-cell levels. Finally, using patient-derived cancer lines, we probed the evolution of EphA signalling in cell populations that underwent metastatic transformation and acquisition of drug resistance. Taken together, our scalable approach provides a reliable scoring scheme for EphA clustering that is consistent over multiple carcinomas and can assay heterogeneity of cancer cell populations in a cost- and time-effective manner.

SummaryWhile understanding heterogeneity in molecular signatures across patients underpins precision oncology, there is increasing appreciation for taking intra-tumor heterogeneity into account. Single-cell RNA-seq (scRNA-seq) technologies have facilitated investigations into the role of intra-tumor transcriptomic heterogeneity (ITTH) in tumor biology and evolution, but their application to in silico models of drug response has not been explored. Based on large-scale analysis of cancer omics datasets, we highlight the utility of ITTH for predicting clinical outcomes. We then show that heterogeneous gene expression signatures obtained from scRNA-seq data can be accurately analyzed (80%) in a recommender system framework (CaDRReS-Sc) for in silico drug response prediction. Patient-derived cell lines capturing transcriptomic heterogeneity from primary and metastatic tumors were used as in vitro proxies for validating monotherapy predictions (Pearson r>0.6), as well as optimal drug combinations to target different subclonal populations (>10% improvement). Applying CaDRReS-Sc to the increasing number of publicly available tumor scRNA-seq datasets can serve as an in silico screen for further in vitro and in vivo drug repurposing studies.Graphical abstractHighlightsLarge-scale analysis to establish the impact of transcriptomic heterogeneity within tumors on clinical outcomesCalibrated recommender system for drug response prediction based on single-cell RNA-seq data (CaDRReS-Sc)Prediction of drug response in patient-derived cell lines with transcriptomic heterogeneityIn silico identification of drug combinations that work based on clonal vulnerabilities

Cancer treatment is gradually evolving from the classical use of nonspecific cytotoxic drugs targeting generic mechanisms of cell growth and proliferation. Instead, new “patient‐specific treatment paradigms” that are based on an individual patient's tumor‐specific molecular features are emerging, and these include “druggable” genomic alterations such as oncogenic driver mutations, downstream activities of cancer‐signaling pathways, and the expression of specific genes involved in tumorigenesis and cancer progression. This evolving landscape of making evidence‐based treatment decisions forms the foundation of precision oncology, which aims to deliver “the right drug, to the right patient and at the right time”. The long‐term vision for this approach is to maximize the treatment efficacy while minimizing exposure to ineffective therapy and reducing co‐morbidity‐related side effects. Successful clinical translation and implementation of this vision have the potential to revolutionize treatment paradigms from predominantly reactive, to more evidence‐based, proactive and predictive care. In this article, we review the past and current approaches in precision oncology, and describe their remarkable power and limitations. We also speculate on the evolution of newly emerging methodologies of the future that can be used to address some of the key challenges associated with the existing paradigms. This article is categorized under: Cancer > Genetics/Genomics/Epigenetics Cancer > Molecular and Cellular Physiology Cancer > Computational Models