Fig. 2

Heterogeneity analysis of tumor cells between KO and WT group. a KEGG result shows the enrichment of signal pathway based on positive top50 expressed gene in KO group. b KEGG result shows the enrichment of signal pathway based on negative top50 expressed gene in KO group. c and d GSEA analysis showed the regulation trend of signal pathway between KO and WT group. e GSVA analysis of KEGG between KO and WT group. f Heatmap of cell interactions in KO group analyzed by CellphoneDB. g The number of ligand–receptor pairs in KO group were shown in the net plot. h Top20 receptor-ligand pair of KO group was shown in heatmap. i Heatmap of cell interactions in WT group. j The number of ligand–receptor pairs in WT group were shown in the net plot. k Top20 receptor-ligand pair of KO group was shown in heatmap. l Top5 different transcription factors between KO and WT group were shown in heatmap. m Different Regulation specificity score (RSS) of the top5 transcription factors were shown between KO and WT group