Monocyte subpopulation reflects coordinated changes with both the immunotherapy progression and MCC severity.

(A) UMAP embedding of all monocytes colored by monocyte compartment subtype. (B) Violin plots showing the expression of non-classical and classical monocytes, S100A, HLA-DR/DP markers in monocyte subtypes. (C) Several cell phenotypic markers were used to identify monocyte cell subtypes. (D) The proportion changes of two monocyte subtypes at each time point. (E) Dot plot for disease trajectory signature genes in classical monocytes. Color discriminates genes with upregulation (red) or downregulation (blue) expression, and point size represents the number of cells per group expressing the corresponding gene. (F) Dot plot for disease trajectory signature genes in non-classical monocytes. (G) Number of DEGs between RespD376 and ARD614 in two monocyte subtypes. (H) The top 10 GSVA results of non-classical monocytes. Y-axis represents the expression level of the entire pathway. (I) Changes in levels of “antigen processing and presentation pathways” in two monocyte subtypes. (J) Correlation between pathway “primary immunodeficiency” and pathway “antigen processing and presentation” in monocytes.