This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Output FindMarkers()

Hi all,

Would you please share why I got peaks instead of genes when using this command:

compare_0 <- FindMarkers(merged_seurat, ident.1 = 'WT_0', ident.2 = 'MT_0')
head(compare_0)

enter image description here

Thank you so much!

seurat
• 1,888 views • link
written by Chris • Chris

1 answer

Hey, it looks like you are doing scATAC-seq data analysis. Here you are checking for Differentially accessible peaks between cell types.

• 0 views • link
written by bk11 • bk11

enter image description hereHi bk11, I am doing single cell multiome RNA + ATAC but I only read RNA assay into RStudio to analyze first.

• 0 views • link
written by Chris • Chris

You need to change the default assays before running FindMarkers then-

#To find markers in your RNA data slot
DefaultAssay(merged_seurat) <- "RNA"
compare_0 <- FindMarkers(merged_seurat, ident.1 = 'WT_0', ident.2 = 'MT_0')

#To find markers in your ATAC data slot
DefaultAssay(merged_seurat) <- "ATAC"
compare_0 <- FindMarkers(merged_seurat, ident.1 = 'WT_0', ident.2 = 'MT_0')
• 0 views • link
written by bk11 • bk11

Unfortunately, it doesn't work as you see the only assay I have in this case is RNA.

• 0 views • link
written by Chris • Chris

Then your upstream steps of FindMarkers are off. I would make sure I read the correct data.

• 0 views • link
written by bk11 • bk11

I used Readmtx() to read 3 files: matrix, feature and barcode. The file filtered_feature_bc_matrix.h5 will have ATAC seq assay but I didn't run.

• 0 views • link
written by Chris • Chris

Log in to answer this question.

More posts like this