Thanks to visit codestin.com
Credit goes to github.com

Skip to content

Adding SCENIC AUC matrix is not working #39

@saeedfc

Description

@saeedfc

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)?

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions