Relationships between KRT14, KRT15, and KRT10 in human foreskin (complement to Fig. 2 ). (A and B) UMAP of human foreskin highlighting (A) KRT15 (green) vs. KRT14 levels (red) and (B) KRT15 (green) vs. KRT10 levels (red) across Seurat clusters. (C) sc expression of KRT15 mapped against KRT14 levels in 24,768 keratinocytes in the foreskin dataset. (D) Distribution of KRT14 and KRT15 expression across all keratinocytes in the dataset. (E) Distribution of KRT15, KRT14, and KRT10 expression in cells exhibiting increasing levels of basal score signature genes (composite scores defined in Materials and methods). (F) Analysis of average KRT14 (red, circle), KRT15 (green, square), and KRT10 (blue, triangle) expression level across nine Seurat clusters matching keratinocyte signatures in the foreskin dataset. Clusters are sorted on the x axis from those showing highest basal score to highest differentiation score. Seurat cluster matching the highest level of G2/M composite score is highlighted. (G) Expression of G2/M composite score relative to KRT14 and KRT15 levels across keratinocytes in the foreskin. (H) Comparison of linear pairwise expression between type I and type II keratins KRT14/KRT5 (red, R² = 0.86), KRT10/KRT1 (blue, R² = 0.92), and KRT15/KRT5 (green, R² = 0.66) across all keratinocytes in the foreskin dataset. (H′) Heatmap summarizing Pearson’s (r) correlations between type I and type II keratins across the dataset.