Integrating bulk RNA-seq and scRNA-seq analyses revealed the function and clinical value of thrombospondins in colon cancer
Jing Li,Ying Tang,Fei Long,Luyao Tian,Ao Tang,LiHui Ding,Juan Chen,Mingwei Liu
DOI: https://doi.org/10.1016/j.csbj.2024.05.021
IF: 6.155
2024-05-19
Computational and Structural Biotechnology Journal
Abstract:Highlights • Exploring the role of thrombospondins with multi-dimensional data, including bulk RNA-seq, scRNA-seq, spatial transcriptomics, cell line drug sensitivity, Whole Genome Sequencing, and immunohistochemical images. • We discovered that pan-dCAFs are one of the major sources of the TSP family in the colon cancer microenvironment. This CAFs subtype possesses high matrix remodeling activity. • Based on TSPs expression, patients can be divided into two clusters with distinct prognoses and immune microenvironment. • We found that THBS2 may be associated with CAFs differentiation and enhanced pro-carcinogenic activity. Moreover, we observed that the communication intensity between THBS2 + CAFs and other cells was stronger than that in THBS2 - CAFs. • A web server presenting the main findings of this study is available at http://tsp.liumwei.org/Col.Cancer . This includes all cohorts and tools used, results, and the appendix. 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 1,981 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. Graphical abstract Download : Download high-res image (213KB) Download : Download full-size image
biochemistry & molecular biology