Transcriptomic characteristics of circulating myeloid cells in patients with high cfBRAF^V600E. (A) Dot plot showing the average expression level and percentage of cells expressing DEGs involved in the RAS-MAPK-ERK pathway in cell subsets from each mutational group compared with healthy control. Genes with percentage expression >50% in each group are selected for visualization. Genes of RAS-MAPK-ERK pathway originate from GO:0000165, GO:0070372, and GO:0007265. (B) Boxplot showing the scores of the 2 RAS-MAPK-ERK related pathways for each cell using ssGSEA. Differences between groups are compared with Wilcoxon ran sum test with P value adjusted for multiple testing using Benjamini-Hochberg correction. Adjusted P < .05 is considered to be statistically significant. ns, P = 1; *P < .05; **P < .01; ***P < .001; ****P = 0. Cell number of 3 mutational groups of each cell cluster is shown in right column. (C) Correlation network analysis of DEGs (including transcription factors [TFs] and genes involved in RAS-MAPK-ERK related pathways) of CD14⁺ Mo, CD16⁺⁺ Mo, and iMo from high-mutation patients. Pink and brown represent the upregulated genes involved in RAS-MAPK-ERK related pathways and TFs, respectively. Spearman correlation coefficient is calculated using genes with average expression >0.3, and genes with correlation score >0.2 are selected for visualization using cytoscape.