Implications of Intratumor Heterogeneity on Consensus Molecular Subtype (CMS) in Colorectal Cancer

Consensus molecular subtyping of metastatic colorectal cancer expands biomarker-directed therapeutic benefit for patients with CMS1 and CMS2 tumors

BACKGROUND

We developed a whole transcriptome sequencing (WTS)-based Consensus Molecular Subtypes (CMS) classifier using FFPE tissue and investigated its prognostic and predictive utility in a large clinico-genomic database of CRC patients (n = 24,939).

METHODS

The classifier was trained against the original CMS datasets using an SVM model and validated in an independent blinded TCGA dataset (88.0% accuracy). Kaplan-Meier estimates of overall survival (OS) and time-on-treatment (TOT) were calculated for each CMS (p < 0.05 considered significant).

RESULTS

CMS2 tumors were enriched on left-side of colon and conferred the longest median OS. In RAS-wildtype mCRC, left-sided tumors and CMS2 classification were associated with longer TOT with anti-EGFR antibodies (cetuximab and panitumumab). When restricting to only CMS2, there was no significant difference in TOT between right- versus left-sided tumors. CMS1 tumors were associated with a longer median TOT with pembrolizumab relative to other CMS groups, even when analyzing only microsatellite stable (MSS) tumors.

DISCUSSION

A WTS-based CMS classifier allowed investigation of a large multi-institutional clinico-genomic mCRC cohort, suggesting anti-EGFR therapy benefit for right-sided RAS-WT CMS2 tumors and immune checkpoint inhibitor benefit for MSS CMS1. Routine CMS classification of CRC provides important treatment associations that should be further investigated.

Gas-Phase Fractionation Data-Independent Acquisition Analysis of 3D Cocultured Spheroid Tumor Model Reveals Altered Translational Processes and Signaling Using Proteomics

Colorectal cancer (CRC) contains considerable heterogeneity; therefore, models of the disease must also reflect the multifarious components. Compared to traditional 2D models, 3D cellular models, such as tumor spheroids, have the utility to determine the drug efficacy of potential therapeutics. Monoculture spheroids are well-known to recapitulate gene expression, cell signaling, and pathophysiological gradients of avascularized tumors. However, they fail to mimic the stromal cell influence present in CRC, which is known to perturb drug efficacy and is associated with metastatic, late-stage colorectal cancer. This study seeks to develop a cocultured spheroid model using carcinoma and noncancerous fibroblast cells. We characterized the proteomic profile of cocultured spheroids in comparison to monocultured spheroids using data-independent acquisition with gas-phase fractionation. Specifically, we determined that proteomic differences related to translation and mTOR signaling are significantly increased in cocultured spheroids compared to monocultured spheroids. Proteins related to fibroblast function, such as exocytosis of coated vesicles and secretion of growth factors, were significantly differentially expressed in the cocultured spheroids. Finally, we compared the proteomic profiles of both the monocultured and cocultured spheroids against a publicly available data set derived from solid CRC tumors. We found that the proteome of the cocultured spheroids more closely resembles that of the patient samples, indicating their potential as tumor mimics.

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.

Single cell analyses of cancer cells identified two regulatorily and functionally distinct categories in differentially expressed genes among tumor subclones

To explore the feature of cancer cells and tumor subclones, we analyzed 101,065 single-cell transcriptomes from 12 colorectal cancer (CRC) patients and 92 single cell genomes from one of these patients. We found cancer cells, endothelial cells and stromal cells in tumor tissue expressed much more genes and had stronger cell-cell interactions than their counterparts in normal tissue. We identified copy number variations (CNVs) in each cancer cell and found correlation between gene copy number and expression level in cancer cells at single cell resolution. Analysis of tumor subclones inferred by CNVs showed accumulation of mutations in each tumor subclone along lineage trajectories. We found differentially expressed genes (DEGs) between tumor subclones had two populations: DEGCNV and DEGreg. DEGCNV, showing high CNV-expression correlation and whose expression differences depend on the differences of CNV level, enriched in housekeeping genes and cell adhesion associated genes. DEGreg, showing low CNV-expression correlation and mainly in low CNV variation regions and regions without CNVs, enriched in cytokine signaling genes. Furthermore, cell-cell communication analyses showed that DEGCNV tends to involve in cell-cell contact while DEGreg tends to involve in secreted signaling, which further support that DEGCNV and DEGreg are two regulatorily and functionally distinct categories.

Stromal Signals Dominate Gene Expression Signature Scores That Aim to Describe Cancer Cell-intrinsic Stemness or Mesenchymality Characteristics

Epithelial-to-mesenchymal transition (EMT) in cancer cells confers migratory abilities, a crucial aspect in the metastasis of tumors that frequently leads to death. In multiple studies, authors proposed gene expression signatures for EMT, stemness, or mesenchymality of tumors based on bulk tumor expression profiling. However, recent studies suggested that noncancerous cells from the microenvironment or macroenvironment heavily influence such signature profiles. Here, we strengthen these findings by investigating 11 published and frequently referenced gene expression signatures that were proposed to describe EMT-related (EMT, mesenchymal, or stemness) characteristics in various cancer types. By analyses of bulk, single-cell, and pseudobulk expression data, we show that the cell type composition of a tumor sample frequently dominates scores of these EMT-related signatures. A comprehensive, integrated analysis of bulk RNA sequencing (RNA-seq) and single-cell RNA-seq data shows that stromal cells, most often fibroblasts, are the main drivers of EMT-related signature scores. We call attention to the risk of false conclusions about tumor properties when interpreting EMT-related signatures, especially in a clinical setting: high patient scores of EMT-related signatures or calls of "stemness subtypes" often result from low cancer cell content in tumor biopsies rather than cancer cell-specific stemness or mesenchymal/EMT characteristics.

SIGNIFICANCE

Cancer self-renewal and migratory abilities are often characterized via gene module expression profiles, also called EMT or stemness gene expression signatures. Using published clinical tumor samples, cancer cell lines, and single cancer cells, we highlight the dominating influence of noncancer cells in low cancer cell content biopsies on their scores. We caution on their application for low cancer cell content clinical cancer samples with the intent to assign such characteristics or subtypes.

Characterizing heterogeneity along EMT and metabolic axes in colorectal cancer reveals underlying consensus molecular subtype-specific trends

Colorectal cancer (CRC) is highly heterogeneous with variable survival outcomes and therapeutic vulnerabilities. A commonly used classification system in CRC is the Consensus Molecular Subtypes (CMS) based on gene expression patterns. However, how these CMS categories connect to axes of phenotypic plasticity and heterogeneity remains unclear. Here, in our analysis of CMS-specific TCGA data and 101 bulk transcriptomic datasets, we found the epithelial phenotype score to be consistently positively correlated with scores of glycolysis, OXPHOS and FAO pathways, while PD-L1 activity scores positively correlated with mesenchymal phenotype scoring, revealing possible interconnections among plasticity axes. Single-cell RNA-sequencing analysis of patient samples revealed that that CMS2 and CMS3 subtype samples were relatively more epithelial as compared to CMS1 and CMS4. CMS1 revealed two subpopulations: one close to CMS4 (more mesenchymal) and the other closer to CMS2 or CMS3 (more epithelial), indicating a partial EMT-like behavior. Consistent observations were made in single-cell analysis of metabolic axes and PD-L1 activity scores. Together, our results quantify the patterns of two functional interconnected axes of phenotypic heterogeneity - EMT and metabolic reprogramming - in a CMS-specific manner in CRC.

Profiling the heterogeneity of colorectal cancer consensus molecular subtypes using spatial transcriptomics

The consensus molecular subtypes (CMS) of colorectal cancer (CRC) is the most widely-used gene expression-based classification and has contributed to a better understanding of disease heterogeneity and prognosis. Nevertheless, CMS intratumoral heterogeneity restricts its clinical application, stressing the necessity of further characterizing the composition and architecture of CRC. Here, we used Spatial Transcriptomics (ST) in combination with single-cell RNA sequencing (scRNA-seq) to decipher the spatially resolved cellular and molecular composition of CRC. In addition to mapping the intratumoral heterogeneity of CMS and their microenvironment, we identified cell communication events in the tumor-stroma interface of CMS2 carcinomas. This includes tumor growth-inhibiting as well as -activating signals, such as the potential regulation of the ETV4 transcriptional activity by DCN or the PLAU-PLAUR ligand-receptor interaction. Our study illustrates the potential of ST to resolve CRC molecular heterogeneity and thereby help advance personalized therapy.

Can immunotherapy reinforce chemotherapy efficacy? a new perspective on colorectal cancer treatment

As one of the main threats to human life (the fourth most dangerous and prevalent cancer), colorectal cancer affects many people yearly, decreases patients' quality of life, and causes irreparable financial and social damages. In addition, this type of cancer can metastasize and involve the liver in advanced stages. However, current treatments can't completely eradicate this disease. Chemotherapy and subsequent surgery can be mentioned among the current main treatments for this disease. Chemotherapy has many side effects, and regarding the treatment of this type of tumor, chemotherapy can lead to liver damage, such as steatohepatitis, steatosis, and sinus damage. These damages can eventually lead to liver failure and loss of its functions. Therefore, it seems that other treatments can be used in addition to chemotherapy to increase its efficiency and reduce its side effects. Biological therapies and immunotherapy are one of the leading suggestions for combined treatment. Antibodies (immune checkpoint blockers) and cell therapy (DC and CAR-T cells) are among the immune system-based treatments used to treat tumors. Immunotherapy targets various aspects of the tumor that may lead to 1) the recruitment of immune cells, 2) increasing the immunogenicity of tumor cells, and 3) leading to the elimination of inhibitory mechanisms established by the tumor. Therefore, immunotherapy can be used as a complementary treatment along with chemotherapy. This review will discuss different chemotherapy and immunotherapy methods for colorectal cancer. Then we will talk about the studies that have dealt with combined treatment.

Computational pathology to improve biomarker testing in breast cancer: how close are we?

The recent advancements in breast cancer precision medicine have highlighted the urgency for the precise and reproducible characterization of clinically actionable biomarkers. Despite numerous standardization efforts, biomarker testing by conventional methodologies is challenged by several issues such as high inter-observer variabilities, the spatial heterogeneity of biomarkers expression, and technological heterogeneity. In this respect, artificial intelligence-based digital pathology approaches are being increasingly recognized as promising methods for biomarker testing and subsequently improved clinical management. Here, we provide an overview on the most recent advances for artificial intelligence-assisted biomarkers testing in breast cancer, with a particular focus on tumor-infiltrating lymphocytes, programmed death-ligand 1, phosphatidylinositol-3 kinase catalytic alpha, and estrogen receptor 1. Challenges and solutions for this integrative analysis in pathology laboratories are also provided.

Long Non-Coding RNAs in Colorectal Cancer: Navigating the Intersections of Immunity, Intercellular Communication, and Therapeutic Potential

This comprehensive review elucidates the intricate roles of long non-coding RNAs (lncRNAs) within the colorectal cancer (CRC) microenvironment, intersecting the domains of immunity, intercellular communication, and therapeutic potential. lncRNAs, which are significantly involved in the pathogenesis of CRC, immune evasion, and the treatment response to CRC, have crucial implications in inflammation and serve as promising candidates for novel therapeutic strategies and biomarkers. This review scrutinizes the interaction of lncRNAs with the Consensus Molecular Subtypes (CMSs) of CRC, their complex interplay with the tumor stroma affecting immunity and inflammation, and their conveyance via extracellular vesicles, particularly exosomes. Furthermore, we delve into the intricate relationship between lncRNAs and other non-coding RNAs, including microRNAs and circular RNAs, in mediating cell-to-cell communication within the CRC microenvironment. Lastly, we propose potential strategies to manipulate lncRNAs to enhance anti-tumor immunity, thereby underlining the significance of lncRNAs in devising innovative therapeutic interventions in CRC.

Identification of collateral sensitivity and evolutionary landscape of chemotherapy-induced drug resistance using cellular barcoding technology

Background: One of the most significant challenges impeding cancer treatment effectiveness is drug resistance. Combining evolutionary understanding with drug resistance can pave the way for the identification of second-line drug options that can overcome drug resistance. Although capecitabine and irinotecan are commonly used therapeutic agents in the treatment of CRC patients, resistance to these agents is common. The underlying clonal dynamics of resistance to these agents using high-resolution barcode technology and identification of effective second-line drugs in this context remain unclear. Methods and materials: Caco-2 and HT-29 cell lines were barcoded, and then capecitabine and irinotecan resistant derivatives of these cell lines were established. The frequencies of barcodes from resistant cell lines and harvested medium, longitudinally, were determined. Collateral drug sensitivity testing was carried out on resistant Caco-2 and HT-29 cell lines using single agents or drug combinations. The SyngeryFinder tool was used to analyse drug combination testing. Results: In Caco-2 and HT-29 cell lines, barcode frequency measurements revealed clonal dynamics of capecitabine and irinotecan formed by both pre-existing and de novo barcodes, indicating the presence of polyclonal drug resistance. The temporal dynamics of clonal evolution in Caco-2 and HT-29 cell lines were demonstrated by longitudinal analysis of pre-existing and de novo barcodes from harvested medium. In Caco-2 and HT-29 cell lines, collateral drug sensitivity revealed a number of drugs that were effective alone and in combination. Conclusion: The use of barcoding technology reveals the clonal dynamics of chemotherapy-induced drug resistance not only from harvested cell populations, but also from longitudinal sampling throughout the course of clonal evolution. Second-line drugs that sensitize drug-resistant CRC cell lines are identified through collateral drug testing.

Multi-label transcriptional classification of colorectal cancer reflects tumor cell population heterogeneity

BACKGROUND

Transcriptional classification has been used to stratify colorectal cancer (CRC) into molecular subtypes with distinct biological and clinical features. However, it is not clear whether such subtypes represent discrete, mutually exclusive entities or molecular/phenotypic states with potential overlap. Therefore, we focused on the CRC Intrinsic Subtype (CRIS) classifier and evaluated whether assigning multiple CRIS subtypes to the same sample provides additional clinically and biologically relevant information.

METHODS

A multi-label version of the CRIS classifier (multiCRIS) was applied to newly generated RNA-seq profiles from 606 CRC patient-derived xenografts (PDXs), together with human CRC bulk and single-cell RNA-seq datasets. Biological and clinical associations of single- and multi-label CRIS were compared. Finally, a machine learning-based multi-label CRIS predictor (ML2CRIS) was developed for single-sample classification.

RESULTS

Surprisingly, about half of the CRC cases could be significantly assigned to more than one CRIS subtype. Single-cell RNA-seq analysis revealed that multiple CRIS membership can be a consequence of the concomitant presence of cells of different CRIS class or, less frequently, of cells with hybrid phenotype. Multi-label assignments were found to improve prediction of CRC prognosis and response to treatment. Finally, the ML2CRIS classifier was validated for retaining the same biological and clinical associations also in the context of single-sample classification.

CONCLUSIONS

These results show that CRIS subtypes retain their biological and clinical features even when concomitantly assigned to the same CRC sample. This approach could be potentially extended to other cancer types and classification systems.

Mimicking Tumor Cell Heterogeneity of Colorectal Cancer in a Patient-derived Organoid-Fibroblast Model

BACKGROUND & AIMS

Patient-derived organoid cancer models are generated from epithelial tumor cells and reflect tumor characteristics. However, they lack the complexity of the tumor microenvironment, which is a key driver of tumorigenesis and therapy response. Here, we developed a colorectal cancer organoid model that incorporates matched epithelial cells and stromal fibroblasts.

METHODS

Primary fibroblasts and tumor cells were isolated from colorectal cancer specimens. Fibroblasts were characterized for their proteome, secretome, and gene expression signatures. Fibroblast/organoid co-cultures were analyzed by immunohistochemistry and compared with their tissue of origin, as well as on gene expression levels compared with standard organoid models. Bioinformatics deconvolution was used to calculate cellular proportions of cell subsets in organoids based on single-cell RNA sequencing data.

RESULTS

Normal primary fibroblasts, isolated from tumor adjacent tissue, and cancer associated fibroblasts retained their molecular characteristics in vitro, including higher motility of cancer associated compared with normal fibroblasts. Importantly, both cancer-associated fibroblasts and normal fibroblasts supported cancer cell proliferation in 3D co-cultures, without the addition of classical niche factors. Organoids grown together with fibroblasts displayed a larger cellular heterogeneity of tumor cells compared with mono-cultures and closely resembled the in vivo tumor morphology. Additionally, we observed a mutual crosstalk between tumor cells and fibroblasts in the co-cultures. This was manifested by considerably deregulated pathways such as cell-cell communication and extracellular matrix remodeling in the organoids. Thrombospondin-1 was identified as a critical factor for fibroblast invasiveness.

CONCLUSION

We developed a physiological tumor/stroma model, which will be vital as a personalized tumor model to study disease mechanisms and therapy response in colorectal cancer.

Improved Delineation of Colorectal Cancer Molecular Subtypes and Functional Profiles with a 62-Gene Panel

Since its establishment in 2015, the transcriptomics-based consensus molecular subtype (CMS) classification has unified our understanding of colorectal cancer. Each of the four CMS exhibited distinctive high-level molecular signatures that correlated well with prognosis and treatment response. Nonetheless, many key aspects of colorectal cancer progression and intra-subtype heterogeneity remain unresolved. This is partly because the bulk transcriptomic data used to define CMS contain substantial interference from non-tumor cells. Here, we propose a concise panel of 62 genes that not only accurately recapitulates all key characteristics of the four original CMS but also identifies three additional subpopulations with unique molecular signatures. Validation on independent cohorts confirms that the new CMS4 intra-subtypes coincide with single-cell-derived intrinsic subtypes and that the panel consists of many immune cell-type markers that can capture the status of tumor microenvironment. Furthermore, a 2D embedding of CMS structure based on the proposed gene panel provides a high-resolution view of the functional pathways and cell-type markers that underlie each CMS intra-subtype and the continuous progression from CMS2 to CMS4 subtypes. Our gene panel and 2D visualization refined the delineation of colorectal cancer subtypes and could aid further discovery of molecular mechanisms in colorectal cancer.

IMPLICATIONS

: Well-selected gene panel and representation can capture both the continuum of cancer cell states and tumor microenvironment status.

Genetic and microenvironmental intra-tumor heterogeneity impacts colorectal cancer evolution and metastatic development

Colorectal cancer (CRC) is a highly diverse disease, where different genomic instability pathways shape genetic clonal diversity and tumor microenvironment. Although intra-tumor heterogeneity has been characterized in primary tumors, its origin and consequences in CRC outcome is not fully understood. Therefore, we assessed intra- and inter-tumor heterogeneity of a prospective cohort of 136 CRC samples. We demonstrate that CRC diversity is forged by asynchronous forms of molecular alterations, where mutational and chromosomal instability collectively boost CRC genetic and microenvironment intra-tumor heterogeneity. We were able to depict predictor signatures of cancer-related genes that can foresee heterogeneity levels across the different tumor consensus molecular subtypes (CMS) and primary tumor location. Finally, we show that high genetic and microenvironment heterogeneity are associated with lower metastatic potential, whereas late-emerging copy number variations favor metastasis development and polyclonal seeding. This study provides an exhaustive portrait of the interplay between genetic and microenvironment intra-tumor heterogeneity across CMS subtypes, depicting molecular events with predictive value of CRC progression and metastasis development.

Refining colorectal cancer classification and clinical stratification through a single-cell atlas

BACKGROUND

Colorectal cancer (CRC) consensus molecular subtypes (CMS) have different immunological, stromal cell, and clinicopathological characteristics. Single-cell characterization of CMS subtype tumor microenvironments is required to elucidate mechanisms of tumor and stroma cell contributions to pathogenesis which may advance subtype-specific therapeutic development. We interrogate racially diverse human CRC samples and analyze multiple independent external cohorts for a total of 487,829 single cells enabling high-resolution depiction of the cellular diversity and heterogeneity within the tumor and microenvironmental cells.

RESULTS

Tumor cells recapitulate individual CMS subgroups yet exhibit significant intratumoral CMS heterogeneity. Both CMS1 microsatellite instability (MSI-H) CRCs and microsatellite stable (MSS) CRC demonstrate similar pathway activations at the tumor epithelial level. However, CD8+ cytotoxic T cell phenotype infiltration in MSI-H CRCs may explain why these tumors respond to immune checkpoint inhibitors. Cellular transcriptomic profiles in CRC exist in a tumor immune stromal continuum in contrast to discrete subtypes proposed by studies utilizing bulk transcriptomics. We note a dichotomy in tumor microenvironments across CMS subgroups exists by which patients with high cancer-associated fibroblasts (CAFs) and C1Q+TAM content exhibit poor outcomes, providing a higher level of personalization and precision than would distinct subtypes. Additionally, we discover CAF subtypes known to be associated with immunotherapy resistance.

CONCLUSIONS

Distinct CAFs and C1Q+ TAMs are sufficient to explain CMS predictive ability and a simpler signature based on these cellular phenotypes could stratify CRC patient prognosis with greater precision. Therapeutically targeting specific CAF subtypes and C1Q + TAMs may promote immunotherapy responses in CRC patients.

Artificial intelligence, molecular subtyping, biomarkers, and precision oncology

A targeted cancer therapy is only useful if there is a way to accurately identify the tumors that are susceptible to that therapy. Thus rapid expansion in the number of available targeted cancer treatments has been accompanied by a robust effort to subdivide the traditional histological and anatomical tumor classifications into molecularly defined subtypes. This review highlights the history of the paired evolution of targeted therapies and biomarkers, reviews currently used methods for subtype identification, and discusses challenges to the implementation of precision oncology as well as possible solutions.