Figure S4D

Code
sample.astro <- readRDS("/omics/odcf/analysis/OE0145_projects/idh_gliomas/Figures_Science/revision/IDH_gliomas_book/datasets/integrated_primary_AS.rds")
sample.astro$orig.ident <- stringr::str_remove_all(sample.astro$orig.ident, "OE0145-")

scale.ident <- c("IDH_ACB_AD_785" = "#ae2012",
                 "IDH_ACB_AD_832" = "#ca6702",
                 "IDH_ACB_AD_865" = "#ee9b00",
                 "IDH_NCH2018" = "#e9d8a6",
                 "IDH_NCH2157" = "#457b9d")

# Define set of markers.
markers.suva <- readRDS("/omics/odcf/analysis/OE0145_projects/idh_gliomas/Figures_Science/revision/IDH_gliomas_book/datasets/IDH_gliomas_TB_annotation_kit_with_Suva_programs_and_metaprogram_iterations.rds")

sample.astro <- SCpubr:::compute_enrichment_scores(sample = sample.astro,
                                                   input_gene_list = list("Stemness_program" = markers.suva$Suva_Stemness_Program),
                                                   flavor = "Seurat")
Code
p <- SCpubr::do_BoxPlot(sample = sample.astro[, sample.astro$relabelling == "RA"],
                        feature = "Stemness_program",
                        order = TRUE,
                        group.by = "orig.ident",
                        colors.use = scale.ident,
                        font.size = 20,
                        axis.text.face = "plain") +
     ggplot2::ylab("Stemness enrichment score") +
     ggplot2::xlab("Patient")