Table S1. Differentially expressed genes of THBS1 obtained based on edgeR package

Table S2.Differentially expressed genes of THBS2 obtained based on edgeR package

Table S3. Differentially expressed genes of THBS3 obtained based on edgeR package

Table S4.Differentially expressed genes of THBS4 obtained based on edgeR package

Table S5. Differentially expressed genes of COMP obtained based on edgeR package

Table S6. Tools and sources

Table S7. Differential expression of TSPs in colon cancer patients with different stages of TCGA-COAD based on wilicox and Kruskal-Wallis analyses

Table S8. Functional pathways are significantly enriched by the TSP family based on GSEA

Table S9. Hallmark pathways score between two clusters based on GSVA

Table S10. Quantification of 28 immune cell infiltrations in patients with two clusters based on ssGSEA

Table S11. Proportion of infiltration of various immune cells in patients with two clusters based on Quantiseq

Table S12. Biomarkers scores between two subtypes based on TIDE

Table S13. Expression levels of T cell exhaustion and immune checkpoint-related genes among clusters

Table S14. Gene Ontology enrichment analysis of five CAFs subtypes

Table S15. Genes closely related to pseudo-time

Table S17. Mean expression levels of interaction partner between various cell types

Table S18. Significance of interaction partner between various cell types

Table S16. Hallmark pathway score between THBS2⁺ and THBS2^- CAFs assessed through GSVA