Fig. 4

Cellular heterogeneity associated with LACC Sub1-6. A and G Uniform manifold approximation and projection (UMAP) plot showing subclusters of ECs (A) or myeloid cells (G). Associations of the individual subclusters with Sub1-6 are indicated in parenthesis. The stacked bar graph shows percentages of ECs belonging to individual subclusters. B Marker genes upregulated predominantly in each subcluster. Numbers in parenthesis represents the number of marker genes in the corresponding subcluster. The color bar represents the gradient of log₂-FC with respect to their mean values. C, H, and K Dot plot showing the enrichment p-value of marker genes for each subcluster of ECs (C), myeloid cells (H), or CAFs (K) in the indicated mRNA and protein signatures for Sub1-6 (see legend in bottom right). The associated subtype for each signature is denoted in parenthesis. Dot size represents the fraction of overlapping marker genes as shown in legend (bottom left). D, I, and L Signature scores of marker genes for each subcluster of ECs (D), myeloid cells (I), or CAFs (L) in Sub1-6. n = 33, 17, 27, 28, 18, and 23 in Sub1-6, respectively. E, J, and M Cellular pathways enriched by marker genes for each subcluster of ECs (E), myeloid cells (J), or CAFs (M). The color bar represents the gradient of –log₁₀(enrichment p-value). F UMAP plot showing ADC ECs with transferred labels of the subclusters of SCC ECs. The stacked bar graph shows percentages of ADC ECs having individual transferred labels