Fig. 2

Major contributors to intratumor heterogeneity. A UMAP plots of malignant epithelial cell reclustering. Each cell is labelled with tissue origin, clinical typing and first clustering result. B Scatter plot demonstrating GSEA results for metabolic pathways weighted by principal component analysis. C Scatter plot demonstrating GSEA enrichment results for hallmark gene sets. D Ranking genes of integrated malignant tumour cells according to principal components. The plot shows GSEA enrichment results for OXPHOS, glycolytic and hypoxic gene sets. E, F Density scatter plot demonstrating glycolysis, OXPHOS activity and response to hypoxia in malignant epithelial cells based on single-cell and bulk data. The colour of the dots indicates the local density. G Pearson correlation between glycolysis, OXPHOS activity and response to hypoxia in other cell types of TME. The height of the bar reflects the correlation coefficient, and the number of * reflects P-value. The * represents P ≤ 0.05, ** represents P ≤ 0.01, *** represents P ≤ 0.001, and ns represents P > 0.05. H Enrichment of upregulated genes in GO from low activity cells in the OXPHOS, glycolytic and hypoxic pathways