Fig. 4
Intra- and intertumoral heterogeneity of cancer-associated fibroblasts. (A) Tumor phylogenetic tree constructed by hierarchical clustering using the 150 branch genes. (B) UMAP plot showing the major cell types in dataset GSE154600. (C) Bar plot showing the origins of cell types in three subtypes of branching evolution. (D) UMAP plot showing the subtypes of cancer-associated fibroblasts (CAFs). (E) Bar plot showing the origins of CAFs in the three evolutionary subtypes. (F) WGCNA results showing the gene modules in distinct CAF subtypes. Columns represent cell types. The colors from blue to red indicate low to high correlation between the gene module and cell subtypes (Pearson correlation test). (G) GO enrichment analysis of hub genes of the BR3 enrichment subtype (F_CXCL12). (H) Number of significant ligand-receptor pairs between CAF and epithelial subtypes. The edge width is proportional to the indicated number of ligand-receptor pairs. EPI_1, epithelial subtype with high expression of MMP7 and ELF3; EPI_3, epithelial subtype with high expression of HES1 and CD24. (I) Dot plot showing the ligand-receptor pairs between CAFs and epithelial cells. Rows represent the ligand receptor (L-R) pairs, and columns represent cell subset–cell subset pairs. The color gradient from black/blue to red indicates the mean values of the L–R pairs from low to high, and the circle size indicates the significance of the pairs. P-values were calculated via a permutation test using CellChat.
