aertslab/SCopeLoomR

Adding SCENIC AUC matrix is not working

Opened this issue · 0 comments

Hi,
Thanks you for Scope and ScopeLoomR.
I tried the following code to add the AUC_Matrix. The loom file works, but I cannot visualize any regulon activity in scope.

library(Seurat)
library(SCopeLoomR)
exp_mtx <- GetAssayData(obj, assay = "RNA", slot = "data")
exp_mtx <- as.matrix(exp_mtx)
## Read in the auc matrix obtained ferom pySCENIC pipeline
scenic_mtx <- read.csv(paste0("../REGA_Lungs/output/", auc_files[1]), row.names = 1)
## Remove the appended '...' character at the end of regulon names
colnames(scenic_mtx) <- gsub(pattern = "\\.\\.\\.$", replacement = "", colnames(scenic_mtx))
## Prepend "TF_" to regulon names
colnames(scenic_mtx)  <- paste0("TF_", colnames(scenic_mtx) )


build_loom(
  file.name=file.name,
  dgem=exp_mtx,
  title=obj_name[i], # A namr for your data
  genome="Mouse", # Just for user information, not used internally
  default.embedding=umap_mtx,
  default.embedding.name="Seurat_UMAP"
)


loom <- open_loom(file.name, mode = 'r+')

add_regulons_AUC(
    loom,
    t(scenic_mtx),
    column.attr.name="MotifRegulonsAUC") 

close_loom(loom= loom)

instead add_regulon_AUC, I also tried the function add_scenic_regulons_auc_matrix(loom=loom, regulons.AUC = t(scenic_mtx))

t(scenic_mtx) is a matrix with same cells as columns and regulons(TFs) as rows.
Code runs without any error. I can open the file in Scope. I can see the gene expression. But when I type transcription factor names, I cannot see anything other than genes ( I mean cannot see regulons). As you might notice, the regulons all start with "TF_", so in theory when i try to type "TF.." in the field in scope, it should pop down with regulon names. Also, I cannot see anything if I go to the 'regulon' tab on scope.

I also tried without adding "TF_" in auc matrix regulon names. In such cases, If I type a transcription factor name, I should in theory see the gene in 'gene' category and 'regulon' category.

I cannot figure out, what could be going wrong.

Is there an alternative to add AUC matrix as an additional assay (I tried with add_matrix function. But didn't work)?