Needed BiocManager::install("DelayedMatrixStats") to run the step "Run Slingshot"
reference_sce <- slingshot(reference_sce, clusterLabels = "seurat_clusters", start.clus = "2", end.clus = "1", reducedDim = "UMAP")
#Rename the Pseudotime column to work with GeneSwitches colData(reference_sce)$Pseudotime <- reference_sce$slingPseudotime_1