# Source code for producing figure 2I

library(ggplot2)
library(dplyr)
library(RColorBrewer)
library(UCell)
library(scGate)
library(viridis)
library(stringr)   
library(patchwork)
library(dittoSeq)
library(gridExtra)
library(Seurat)
library(openxlsx)
library(ggalluvial)
library(gridExtra)

obj <- readRDS(file = "Fig2H_dataset.rds")

umap_wpre <- dittoSeq::dittoDimPlot(object = obj, reduction.use = "umap.harmony.pt.tp",order = "increasing",
                                    var = "WPRE_condition", opacity = 0.8, do.label = F,labels.size = 4,
                                    labels.highlight = F, labels.repel = T, size = 1) + 
  scale_color_manual(values = c("#A0A0A0B3","#FF3333B3")) + 
  xlab("UMAP_1") +
  ylab(label = "UMAP_2") +
  theme(legend.position = "none",
        legend.title = element_blank(),
        legend.key.size = unit(0.001, "in"),
        legend.key.spacing = unit(0,"in"),
        plot.title = element_blank(),
        axis.text = element_text(size = 16)
  )

pdf(file = "Fig2I.pdf", width = 4, height = 4)
umap_wpre
dev.off()