# script for figures 2G

## loading libraries
library(Seurat)
library(ggplot2)
library(ggplotify)
library(ggraph)
library(ggrepel)
library(ggpubr)
library(ComplexHeatmap)
library(gridExtra)
library(RColorBrewer)
library(dittoSeq)
library(openxlsx)
library(reshape2)
library(ggalluvial)
library(dplyr)

# load seurat object

obj <- readRDS("Fig2C_Dataset.rds")
obj <- subset(obj, cells = WhichCells(object = obj,expression = orig.ident %in% c("C04","C02","C03")))
obj <- SetIdent(object = obj, value = "Population")

citokine_dotplot_C04_C02_C03 <- DotPlot(object = obj, 
                                        features = c("KIT","MPL","FLT3","IL6R","IL3RA"),
                                        dot.min = 0.05,
                                        scale.by = "size") + 
  labs(color = "ZscoreExpr", size = "PctExpr")  + 
  scale_color_gradient2(low = "blue", mid = "lightgrey", high = "red") + 
  theme(axis.title = element_blank(), 
        legend.position = "none",           # Places the legend at the bottom
        legend.justification = "left",       # Justifies the legend to the left
        legend.box = "vertical"  
  ) +
  scale_x_discrete(guide = guide_axis(n.dodge = 2))

pdf(file = "Fig2G.pdf", width = 3, height = 4)
citokine_dotplot_C04_C02_C03
dev.off()