Integrating bulk RNA-seq and scRNA-seq analyses revealed the function and clinical value of thrombospondins in colon cancer

Background: Acting as mediators in cell-matrix and cell-cell communication, matricellular proteins play a crucial role in cancer progression. Thrombospondins (TSPs), a type of matricellular glycoproteins, are key regulators in cancer biology with multifaceted roles. Although TSPs have been implicated in anti-tumor immunity and epithelial-mesenchymal transition (EMT) in several malignancies, their specific roles to colon cancer remain elusive. Addressing this knowledge gap is essential, as understanding the function of TSPs in colon cancer could identify new therapeutic targets and prognostic markers. Methods: Analyzing 1981 samples from 10 high-throughput datasets, including six bulk RNA-seq, three scRNA-seq, and one spatial transcriptome dataset, our study investigated the prognostic relevance, risk stratification value, immune heterogeneity, and cellular origin of TSPs, as well as their influence on cancer-associated fibroblasts (CAFs). Utilizing survival analysis, unsupervised clustering, and functional enrichment, along with multiple correlation analyses of the tumor-microenvironment (TME) via Gene Set Variation Analysis (GSVA), spatial localization, Monocle2, and CellPhoneDB, we provided insights into the clinical and cellular implications of TSPs. Results: First, we observed significant upregulation of THBS2 and COMP in colon cancer, both of which displayed significant prognostic value. Additionally, we detected a significant positive correlation between TSPs and immune cells, as well as marker genes of EMT. Second, based on TSPs expression, patients were divided into two clusters with distinct prognoses: the high TSPs expression group (TSPs-H) was characterized by pronounced immune and stromal cell infiltration, and notably elevated T-cell exhaustion scores. Subsequently, we found that THBS2 and COMP may be associated with the differentiation of CAFs into pan-iCAFs and pan-dCAFs, which are known for their heightened matrix remodeling activities. Moreover, THBS2 enhanced CAFs communication with vascular endothelial cells and monocyte-macrophages. CAFs expressing THBS2 (THBS2+ CAFs) demonstrated higher scores across multiple signaling pathways, including angiogenic, EMT, Hedgehog, Notch, Wnt, and TGF-β, when compared to THBS2- CAFs. These observations suggest that THBS2 may be associated with stronger pro-carcinogenic activity in CAFs. Conclusions: This study revealed the crucial role of TSPs and the significant correlation between THBS2 and CAFs interactions in colon cancer progression, providing valuable insights for targeting TSPs to mitigate cancer progression.