Sierra dataframe has 0 length

[top334]

Hi everyone:

I met a problem when I run the DUTest as below, so I find my peaks.seurat@tools$Sierra doesn't have any thing except the rownames.

res.table = DUTest(peaks.seurat,

•                population.1 = "W-P", 
  

•                population.2 = "L-P",
  

•                exp.thresh = 0.1, 
  

•                #feature.type = c("UTR3", "exon")
  

•                )
  

[1] "0 expressed peaks in feature types UTR3, exon" [1] "0 genes detected with multiple peak sites expressed" [1] "0 individual peak sites to test" Error in [.data.frame(Tool(apa.seurat.object, "Sierra"), , c("Gene_name", : 选择了未定义的列 And my annotation_peak.txt as below:

I have no idea why it lost the information of Sierra dataframe when combine the seruat object and annotated peaks. Can anyone help me?

[rj-patrick]

Hi @top334, Could you please give some more information about the steps you took to create your Seurat object from the Sierra output? Without knowing anything else, I'd say the most likely cause for your issue is a conflict in the names of the cell barcodes in your Seurat object vs the Sierra output. Cheers, Ralph