ABSTRACT
BACKGROUND: Recent advancements in high-throughput genomics and targeted therapies have provided tremendous potential to identify and therapeutically target distinct mutations associated with cancers. However, to date the majority of targeted therapies are used to treat all functional mutations within the same gene, regardless of affected codon or phenotype. RESULTS: In this study, we developed a functional genomic analysis workflow with a unique isogenic cell line panel bearing two distinct hotspot PIK3CA mutations, E545K and H1047R, to accurately identify targetable differences between mutations within the same gene. We performed RNA-seq and ATAC-seq and identified distinct transcriptomic and epigenomic differences associated with each PIK3CA hotspot mutation. We used this data to curate a select CRISPR knock out screen to identify mutation-specific gene pathway vulnerabilities. These data revealed AREG as a E545K-preferential target that was further validated through in vitro analysis and publicly available patient databases. CONCLUSIONS: Using our multi-modal genomics framework, we discover distinct differences in genomic regulation between PIK3CA hotspot mutations, suggesting the PIK3CA mutations have different regulatory effects on the function and downstream signaling of the PI3K complex. Our results demonstrate the potential to rapidly uncover mutation specific molecular targets, specifically AREG and a proximal gene regulatory region, that may provide clinically relevant therapeutic targets. The methods outlined provide investigators with an integrative strategy to identify mutation-specific targets for the treatment of other oncogenic mutations in an isogenic system.
Subject(s)
Breast Neoplasms , Class I Phosphatidylinositol 3-Kinases , Genomics , Mutation , Class I Phosphatidylinositol 3-Kinases/genetics , Humans , Breast Neoplasms/genetics , Genomics/methods , Cell Line, Tumor , Female , Gene Expression Regulation, NeoplasticABSTRACT
Formalin-fixed paraffin-embedded (FFPE) samples represent the cornerstone of tissue-based analysis in precision medicine. Targeted next-generation sequencing panels are routinely used to analyze a limited number of genes to guide treatment decision-making for advanced-stage patients. The number and complexity of genetic alterations to be investigated are rapidly growing; in several instances, a comprehensive genomic profiling analysis is needed. The poor quality of genetic material extracted from FFPE samples may impact the feasibility/reliability of sequencing data. We sampled 9 colorectal cancers to allow 4 parallel fixations: (1) neutral buffered formalin (NBF), (2) acid-deprived formalin fixation (ADF), (3) precooled ADF (coldADF), and (4) glyoxal acid free (GAF). DNA extraction, fragmentation analysis, and sequencing by 2 large next-generation sequencing panels (OCAv3 and TSO500) followed. We comprehensively analyzed library and sequencing quality controls and the quality of sequencing results. Libraries from coldADF samples showed significantly longer reads than the others with both panels. ADF-derived and coldADF-derived libraries showed the lowest level of noise and the highest levels of uniformity with the OCAv3 panel, followed by GAF and NBF samples. The data uniformity was confirmed by the TSO500 results, which also highlighted the best performance in terms of the total region sequenced for the ADF and coldADF samples. NBF samples had a significantly smaller region sequenced and displayed a significantly lower number of evaluable microsatellite loci and a significant increase in single-nucleotide variations compared with other protocols. Mutational signature 1 (aging and FFPE artifact related) showed the highest (37%) and lowest (17%) values in the NBF and coldADF samples, respectively. Most of the identified genetic alterations were shared by all samples in each lesion. Five genes showed a different mutational status across samples and/or panels: 4 discordant results involved NBF samples. In conclusion, acid-deprived fixatives (GAF and ADF) guarantee the highest DNA preservation/sequencing performance, thus allowing more complex molecular profiling of tissue samples.
Subject(s)
Artifacts , DNA , Humans , Tissue Fixation/methods , Reproducibility of Results , DNA/genetics , DNA/analysis , Formaldehyde , Genomics , Paraffin Embedding , High-Throughput Nucleotide SequencingABSTRACT
BACKGROUND & AIMS: Only limited therapeutic options are currently available for hepatocellular carcinoma (HCC), making the development of effective alternatives essential. Based on the recent finding that systemic or local hypothyroidism is associated with HCC development in humans and rodents, we investigated whether the thyroid hormone triiodothyronine (T3) could inhibit the progression of HCCs. METHODS: Different rat and mouse models of hepatocarcinogenesis were investigated. The effect of T3 on tumorigenesis and metabolism/differentiation was evaluated by transcriptomic analysis, quantitative reverse transcription PCR, immunohistochemistry, and enzymatic assay. RESULTS: A short treatment with T3 caused a shift in the global expression profile of the most aggressive preneoplastic nodules towards that of normal liver. This genomic reprogramming preceded the disappearance of nodules and involved reprogramming of metabolic genes, as well as pro-differentiating transcription factors, including Kruppel-like factor 9, a target of the thyroid hormone receptor ß (TRß). Treatment of HCC-bearing rats with T3 strongly reduced the number and burden of HCCs. Reactivation of a local T3/TRß axis, a switch from Warburg to oxidative metabolism and loss of markers of poorly differentiated hepatocytes accompanied the reduced burden of HCC. This effect persisted 1 month after T3 withdrawal, suggesting a long-lasting effect of the hormone. The antitumorigenic effect of T3 was further supported by its inhibitory activity on cell growth and the tumorigenic ability of human HCC cell lines. CONCLUSIONS: Collectively, these findings suggest that reactivation of the T3/TRß axis induces differentiation of neoplastic cells towards a more benign phenotype and that T3 or its analogs, particularly agonists of TRß, could be useful tools in HCC therapy. LAY SUMMARY: Hepatocellular carcinoma (HCC) represents an important challenge for global health. Recent findings showed that systemic or local hypothyroidism is associated with HCC development. In rat models, we showed that administration of the thyroid hormone T3 impaired HCC progression, even when given at late stages. This is relevant from a translational point of view as HCC is often diagnosed at an advanced stage when it is no longer amenable to curative treatments. Thyroid hormones and/or thyromimetics could be useful for the treatment of patients with HCC.
Subject(s)
Anticarcinogenic Agents/administration & dosage , Carcinogenesis/drug effects , Carcinogenesis/metabolism , Carcinoma, Hepatocellular/drug therapy , Carcinoma, Hepatocellular/metabolism , Cell Differentiation/drug effects , Disease Progression , Liver Neoplasms/drug therapy , Liver Neoplasms/metabolism , Triiodothyronine/administration & dosage , Aged , Animals , Carcinoma, Hepatocellular/pathology , Disease Models, Animal , Female , Gene Expression Profiling , Hep G2 Cells , Humans , Liver Neoplasms/pathology , Male , Mice , Mice, Nude , Middle Aged , Rats , Rats, Inbred F344 , Rats, Wistar , Thyroid Hormone Receptors beta/metabolism , Transcriptome , Triiodothyronine/metabolismABSTRACT
Background: Recent advancements in high-throughput genomics and targeted therapies have provided tremendous potential to identify and therapeutically target distinct mutations associated with cancers. However, to date the majority of targeted therapies are used to treat all functional mutations within the same gene, regardless of affected codon or phenotype. Results: In this study, we developed a functional genomic analysis workflow with a unique isogenic cell line panel bearing two distinct hotspot PIK3CA mutations, E545K and H1047R, to accurately identify targetable differences between mutations within the same gene. We performed RNA-seq and ATAC-seq and identified distinct transcriptomic and epigenomic differences associated with each PIK3CA hotspot mutation. We used this data to curate a select CRISPR knock out screen to identify mutation-specific gene pathway vulnerabilities. These data revealed AREG as a E545K-preferential target that was further validated through in vitro analysis and publicly available patient databases. Conclusions: Using our multi-modal genomics framework, we discover distinct differences in genomic regulation between PIK3CA hotspot mutations, suggesting the PIK3CA mutations have different regulatory effects on the function and downstream signaling of the PI3K complex. Our results demonstrate the potential to rapidly uncover mutation specific molecular targets, specifically AREG and a proximal gene regulatory region, that may provide clinically relevant therapeutic targets. The methods outlined provide investigators with an integrative strategy to identify mutation-specific targets for the treatment of other oncogenic mutations in an isogenic system.
ABSTRACT
The synthetic lethal effect observed with the use of PARP inhibitors (PARPi) with tumors characterized by the loss of key players in the homologous recombination (HR) pathway, commonly referred to as "BRCAness", is maintaining high interest in oncology. While BRCAness is a well-established feature in breast, ovarian, prostate, and pancreatic carcinomas, our recent findings indicate that up to 15% of colorectal cancers (CRC) also harbor defects in the HR pathway, presenting promising opportunities for innovative therapeutic strategies in CRC patients. We developed a new tool called HRDirect, which builds upon the HRDetect algorithm and is able to predict HR deficiency (HRD) from reference-free tumor samples. We validated HRDirect using matched breast cancer and CRC patient samples. Subsequently, we assessed its efficacy in predicting response to the PARP inhibitor olaparib by comparing it with two other commercial assays: AmoyDx HRD by Amoy Diagnostics and the TruSight Oncology 500 HRD (TSO500-HRD) panel by Illumina NGS technology. While all three approaches successfully identified the most PARPi-sensitive CRC models, HRDirect demonstrated superior precision in distinguishing resistant models compared to AmoyDX and TSO500-HRD, which exhibited overlapping scores between sensitive and resistant cells. Furthermore, we propose integrating HRDirect scoring with ATM and RAD51C immunohistochemical analysis as part of our "composite biomarker approach" to enhance the identification of HRD tumors, with an immediate translational and clinical impact for CRC personalized treatment.
ABSTRACT
Extrachromosomal circular DNAs (eccDNAs) have emerged as important intra-cellular mobile genetic elements that affect gene copy number and exert in trans regulatory roles within the cell nucleus. Here, we describe scCircle-seq, a method for profiling eccDNAs and unraveling their diversity and complexity in single cells. We implement and validate scCircle-seq in normal and cancer cell lines, demonstrating that most eccDNAs vary largely between cells and are stochastically inherited during cell division, although their genomic landscape is cell type-specific and can be used to accurately cluster cells of the same origin. eccDNAs are preferentially produced from chromatin regions enriched in H3K9me3 and H3K27me3 histone marks and are induced during replication stress conditions. Concomitant sequencing of eccDNAs and RNA from the same cell uncovers the absence of correlation between eccDNA copy number and gene expression levels, except for a few oncogenes, including MYC, contained within a large eccDNA in colorectal cancer cells. Lastly, we apply scCircle-seq to one prostate cancer and two breast cancer specimens, revealing cancer-specific eccDNA landscapes and a higher propensity of eccDNAs to form in amplified genomic regions. scCircle-seq is a scalable tool that can be used to dissect the complexity of eccDNAs across different cell and tissue types, and further expands the potential of eccDNAs for cancer diagnostics.
Subject(s)
DNA, Circular , DNA , Male , Humans , DNA, Circular/genetics , Chromosomes , Cell Line , OncogenesABSTRACT
PURPOSE: Trastuzumab is an HER2-specific agent approved as the gold-standard therapy for advanced HER2-positive (HER2+) gastric cancer (GC), but the high rate and rapid appearance of resistance limit its clinical efficacy, resulting in the need to identify new vulnerabilities. Defining the drivers influencing HER2+ cancer stem cell (CSC) maintenance/survival could represent a clinically useful strategy to counteract tumor growth and therapy resistance. Accumulating evidence show that targeting crucial metabolic hubs, as the fatty acid synthase (FASN), may be clinically relevant. METHODS: FASN protein and transcript expression were examined by WB and FACS and by qRT-PCR and GEP analyses, respectively, in trastuzumab-sensitive and trastuzumab-resistant HER2+ GC cell lines cultured in adherent (2D) or gastrosphere promoting (3D) conditions. Molecular data were analyzed in silico in public HER2+ GC datasets. The effectiveness of the FASN inhibitor TVB3166 to overcome anti-HER2 therapy resistance was tested in vitro in gastrospheres forming efficiency bioassays and in vivo in mice bearing trastuzumab-resistant GC cells. RESULTS: We compared the transcriptome profiles of HER2+ GC cells cultured in 2D versus 3D conditions finding a significant enrichment of FASN in 3D cultures. FASN upregulation significantly correlated with high stemness score and poor prognosis in HER2+ GC cases. TVB3166 treatment significantly decreased GCSCs in all cell targets. HER2 and FASN cotargeting significantly decreased the capability to form gastrospheres versus monotherapy and reduced the in vivo growth of trastuzumab-resistant GC cells. CONCLUSION: Our findings indicate that cotargeting HER2 and FASN increase the benefit of anti-HER2 therapy representing a new opportunity for metabolically combating trastuzumab-resistant HER2+ GC.
Subject(s)
Receptor, ErbB-2 , Stomach Neoplasms , Animals , Mice , Receptor, ErbB-2/metabolism , Stomach Neoplasms/pathology , Trastuzumab/pharmacology , Fatty Acid Synthases/metabolism , Fatty Acid Synthases/therapeutic use , Cell Line, TumorABSTRACT
PURPOSE: In JACOB trial, pertuzumab added to trastuzumab-chemotherapy did not significantly improve survival of patients with HER2-positive metastatic gastric cancer, despite 3.3 months increase versus placebo. HER2 copy-number variation (CNV) and AMNESIA panel encompassing primary resistance alterations (KRAS/PIK3CA/MET mutations, KRAS/EGFR/MET amplifications) may improve patients' selection for HER2 inhibition. EXPERIMENTAL DESIGN: In a post hoc analysis of JACOB on 327 samples successfully sequenced by next-generation sequencing (NGS; Oncomine Focus DNA), HER2 CNV, HER2 expression by IHC, and AMNESIA were correlated with overall response rate (ORR), progression-free survival (PFS), and overall survival (OS) by univariable/multivariable models. RESULTS: Median HER2 CNV was 4.7 (interquartile range, 2.2-16.9). HER2 CNV-high versus low using the median as cutoff was associated with longer median PFS (10.5 vs. 6.4 months; HR = 0.48; 95% confidence interval: 0.38-0.62; P < 0.001) and OS (20.3 vs. 13.0 months; HR = 0.54; 0.42-0.72; P < 0.001). Combining HER2 CNV and IHC improved discriminative ability, with better outcomes restricted to HER2-high/HER2 3+ subgroup. AMNESIA positivity was found in 51 (16%), with unadjusted HR = 1.35 (0.98-1.86) for PFS; 1.43 (1.00-2.03) for OS.In multivariable models, only HER2 CNV status remained significant for PFS (P < 0.001) and OS (P = 0.004). Higher ORR was significantly associated with IHC 3+ [61% vs. 34% in 2+; OR = 3.11 (1.89-5.17)] and HER2-high [59% vs. 43% in HER2-low; OR = 1.84 (1.16-2.94)], with highest OR in the top CNV quartile. These biomarkers were not associated with treatment effect of pertuzumab. CONCLUSIONS: HER2 CNV-high assessed by NGS may be associated with better ORR, PFS, and OS in a JACOB subgroup, especially if combined with HER2 3+. The negative prognostic role of AMNESIA requires further clinical validation.
Subject(s)
Breast Neoplasms , Stomach Neoplasms , Humans , Female , Trastuzumab/therapeutic use , Trastuzumab/metabolism , Receptor, ErbB-2/metabolism , Stomach Neoplasms/drug therapy , Stomach Neoplasms/genetics , Stomach Neoplasms/pathology , DNA Copy Number Variations , Proto-Oncogene Proteins p21(ras)/genetics , Antineoplastic Combined Chemotherapy Protocols/adverse effects , Breast Neoplasms/drug therapyABSTRACT
BACKGROUND: Long-Interspersed Nuclear Element (L1) retrotransposons are silenced in healthy tissues but unrepressed in cancer. Even if L1 reactivation has been associated with reduced overall survival in breast cancer (BC) patients, a comprehensive correlation with clinicopathological features is still missing. METHODS: Using quantitative, reverse-transcription PCR, we assessed L1 mRNA expression in 12 BC cells, 210 BC patients and in 47 normal mammary tissues. L1 expression was then correlated with molecular and clinicopathological data. RESULTS: We identified a tumor-exclusive expression of L1s, absent in normal mammary cells and tissues. A positive correlation between L1 expression and tumor dedifferentiation, lymph-node involvement and increased immune infiltration was detected. Molecular subtyping highlighted an enrichment of L1s in basal-like cells and cancers. By exploring disease-free survival, we identified L1 overexpression as an independent biomarker for patients with a high risk of recurrence in hormone-receptor-negative BCs. CONCLUSIONS: Overall, L1 reactivation identified BCs with aggressive features and patients with a worse clinical fate.
Subject(s)
Retroelements , Triple Negative Breast Neoplasms , B7-H1 Antigen/metabolism , Hormones , Humans , Neoplasm Recurrence, Local/genetics , Neoplasm Recurrence, Local/pathology , Triple Negative Breast Neoplasms/pathologyABSTRACT
BACKGROUND: The "HER2-low" nomenclature identifies breast carcinomas (BCs) displaying a HER2 score of 1+/2+ in immunohistochemistry and lacking ERBB2 amplification. Whether HER2-low BCs (HLBCs) constitute a distinct entity is debated. METHODS: We performed DNA and RNA high-throughput analysis on 99 HLBC samples (n = 34 cases with HER2 score 1+/HLBC-1, n = 15 cases with HER2 score 2+ and ERBB2 not amplified/HLBC-2N, and n = 50 cases with score 2+ and ERBB2 copy number in the equivocal range/HLBC-2E). We compared the mutation rates with data from 1317 samples in the Memorial Sloan-Kettering Cancer Center (MSKCC) BC cohort and gene expression data with those from an internal cohort of HER2-negative and HER2-positive BCs. RESULTS: The most represented mutations affected PIK3CA (31/99, 31%), GATA3 (18/99, 18%), TP53 (17/99, 17%), and ERBB2 (8/99, 8%, private to HLBC-2E). Tumor mutational burden was significantly higher in HLBC-1 compared to HLBC-2E/N (P = 0.04). Comparison of mutation spectra revealed that HLBCs were different from both HER2-negative and HER2-positive BCs, with HLBC-1 resembling more HER2-negative tumors and HLBC-2 mutationally related to HER2-addicted tumors. Potentially actionable alterations (annotated by using OncoKB/ESCAT classes) affected 52 patients. Intra-group gene expression revealed overlapping features between HLBC-1 and control HER2-negative BCs, whereas the HLBC-2E tumors showed the highest diversity overall. The RNA-based class discovery analysis unveiled four subsets of tumors with (i) lymphocyte activation, (ii) unique enrichment in HER2-related features, (iii) stromal remodeling alterations, and (iv) actionability of PIK3CA mutations (LAURA classification). CONCLUSIONS: HLBCs harbor distinct genomic features when compared with HER2-positive and HER2-negative BCs; however, differences across IHC classes were also unveiled thus dissecting the full picture of heterogeneity across HER2-low disease. The HLBC-2E category harbors most distinctive features, whereas HLBC-1 seems superimposable to HER2-negative disease. Further studies are needed to ascertain whether the four genomic-driver classes of the LAURA classification hold prognostic and/or predictive implications.
Subject(s)
Breast Neoplasms , Transcriptome , Biomarkers, Tumor/genetics , Breast Neoplasms/pathology , Female , Gene Expression Profiling , Genomics , Humans , Mutation , RNAABSTRACT
In recent years, cytoreductive surgery (CRS) and hyperthermic intraperitoneal chemotherapy (HIPEC) have substantially improved the clinical outcome of pseudomyxoma peritonei (PMP) originating from mucinous appendiceal cancer. However, current histopathological grading of appendiceal PMP frequently fails in predicting disease outcome. We recently observed that the integration of cancer cell transcriptional traits with those of cancer-associated fibroblasts (CAFs) improves prognostic prediction for tumors of the large intestine. We therefore generated global expression profiles on a consecutive series of 24 PMP patients treated with CRS plus HIPEC. Multiple lesions were profiled for nine patients. We then used expression data to stratify the samples by a previously published "high-risk appendiceal cancer" (HRAC) signature and by a CAF signature that we previously developed for colorectal cancer, or by a combination of both. The prognostic value of the HRAC signature was confirmed in our cohort and further improved by integration of the CAF signature. Classification of cases profiled for multiple lesions revealed the existence of outlier samples and highlighted the need of profiling multiple PMP lesions to select representative samples for optimal performances. The integrated predictor was subsequently validated in an independent PMP cohort. These results provide new insights into PMP biology, revealing a previously unrecognized prognostic role of the stromal component and supporting integration of standard pathological grade with the HRAC and CAF transcriptional signatures to better predict disease outcome.
ABSTRACT
PURPOSE: Overexpression of miR-100 in stem cells derived from basal-like breast cancers causes loss of stemness, induction of luminal breast cancer markers and response to endocrine therapy. We, therefore, explored miR-100 as a novel biomarker in patients with luminal breast cancer. METHODS: miR-100 expression was studied in 90 patients with oestrogen-receptor-positive/human-epidermal growth factor receptor 2-negative breast cancer enrolled in a prospective study of endocrine therapy given either preoperatively, or for the treatment of de novo metastatic disease. Response was defined as a Ki67 ≤2.7% after 21±3 days of treatment. The prognostic role of miR-100 expression was evaluated in the Molecular Taxonomy of Breast Cancer International Consortium (METABRIC) and The Cancer Genome Atlas (TCGA) breast cancer datasets. Additionally, we explored the correlation between miR-100 and the expression its targets reported as being associated with endocrine resistance. Finally, we evaluated whether a signature based on miR-100 and its target genes could predict the luminal A molecular subtype. RESULTS: Baseline miR-100 was significantly anticorrelated with baseline and post-treatment Ki67 (p<0.001 and 0.004, respectively), and independently associated with response to treatment (OR 3.329, p=0.047). In the METABRIC dataset, high expression of miR-100 identified women with luminal A tumours treated with adjuvant endocrine therapy with improved overall survival (HR 0.55, p<0.001). miR-100 was negatively correlated with PLK1, FOXA1, mTOR and IGF1R expression, potentially explaining its prognostic effect. Finally, a miR-100-based signature developed in patients enrolled in the prospective study outperformed Ki67 alone in predicting the luminal A phenotype. CONCLUSIONS: Our findings suggest that miR-100 should be further explored as a biomarker in patients with luminal breast cancer.
Subject(s)
Breast Neoplasms , MicroRNAs , Breast Neoplasms/drug therapy , Breast Neoplasms/genetics , Female , Hepatocyte Nuclear Factor 3-alpha , Humans , MicroRNAs/genetics , Prognosis , Prospective StudiesABSTRACT
BACKGROUND: Clustered protocadherins (PCDHs) map in tandem at human chromosome 5q31 and comprise three multi-genes clusters: α-, ß- and γ-PCDH. The expression of this cluster consists of a complex mechanism involving DNA hub formation through DNA-CCTC binding factor (CTCF) interaction. Methylation alterations can affect this interaction, leading to transcriptional dysregulation. In cancer, clustered PCDHs undergo a mechanism of long-range epigenetic silencing by hypermethylation. RESULTS: In this study, we detected frequent methylation alterations at CpG islands associated to these clustered PCDHs in all the solid tumours analysed (colorectal, gastric and biliary tract cancers, pilocytic astrocytoma), but not hematologic neoplasms such as chronic lymphocytic leukemia. Importantly, several altered CpG islands were associated with CTCF binding sites. Interestingly, our analysis revealed a hypomethylation event in pilocytic astrocytoma, suggesting that in neuronal tissue, where PCDHs are highly expressed, these genes become hypomethylated in this type of cancer. On the other hand, in tissues where PCDHs are lowly expressed, these CpG islands are targeted by DNA methylation. In fact, PCDH-associated CpG islands resulted hypermethylated in gastrointestinal tumours. CONCLUSIONS: Our study highlighted a strong alteration of the clustered PCDHs methylation pattern in the analysed solid cancers and suggested these methylation aberrations in the CpG islands associated with PCDH genes as powerful diagnostic biomarkers.
Subject(s)
Cadherins/genetics , DNA Methylation , High-Throughput Nucleotide Sequencing/methods , Neoplasms/genetics , CpG Islands , Epigenesis, Genetic , Gene Expression Regulation, Neoplastic , Humans , Multigene Family , Promoter Regions, Genetic , Sequence Analysis, DNAABSTRACT
Recent studies identified a poor-prognosis stem/serrated/mesenchymal (SSM) transcriptional subtype of colorectal cancer (CRC). We noted that genes upregulated in this subtype are also prominently expressed by stromal cells, suggesting that SSM transcripts could derive from stromal rather than epithelial cancer cells. To test this hypothesis, we analyzed CRC expression data from patient-derived xenografts, where mouse stroma supports human cancer cells. Species-specific expression analysis showed that the mRNA levels of SSM genes were mostly due to stromal expression. Transcriptional signatures built to specifically report the abundance of cancer-associated fibroblasts (CAFs), leukocytes or endothelial cells all had significantly higher expression in human CRC samples of the SSM subtype. High expression of the CAF signature was associated with poor prognosis in untreated CRC, and joint high expression of the stromal signatures predicted resistance to radiotherapy in rectal cancer. These data show that the distinctive transcriptional and clinical features of the SSM subtype can be ascribed to its particularly abundant stromal component.