Fig. 3

Identification of ligand-receptor pair interaction subtypes in CRC bulk RNA-seq data. A Consensus clustering matrix for the optimal cluster number (k = 3) in the GSE39582 dataset. B, D, F Heat maps of marker genes for three ligand-receptor pair subtypes. C, E, G The consistency between each sample and each subtype was evaluated by the silhouette value. H GO-BP enrichment analysis results of the marker genes of each subtype. I The consistency of the identified subtypes in the three bulk RNA-seq datasets can be found through submap analysis. J Kaplan–Meier survival analysis was used to compare the survival differences of each receptor subtype. **P < 0.01