This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Splitting scRNA seq data from Multiome data(scRNA+ scATAC)

Could anyone help me with splitting the scRNA seq data from a Multiome dataset (scATAC + snRNA seq). I'm trying to create a separate anndata from the scRNA seq portion only.

I've tried setting the Default using Signac as RNA and then saving the rds file, but I don't think that works.

scrna-seq signac
• 1,400 views • link
updated by Ram Ram • written by gradstudent • gradstudent

1 answer

You can do split the multiome data in following way-

#Your Multiome data
pbmc
An object of class Seurat 
142657 features across 11909 samples within 2 assays 
Active assay: RNA (36601 features, 0 variable features)
 1 layer present: counts
 1 other assay present: ATAC

#Splitting the scRNA seq data from multiome data
scRNAseq <- CreateSeuratObject(counts = pbmc@assays$RNA)
head(scRNAseq@meta.data)
                      orig.ident nCount_RNA nFeature_RNA
AAACAGCCAAGGAATC-1 SeuratProject       8380         3308
AAACAGCCAATCCCTT-1 SeuratProject       3771         1896
AAACAGCCAATGCGCT-1 SeuratProject       6876         2904
AAACAGCCACACTAAT-1 SeuratProject       1733          846
AAACAGCCACCAACCG-1 SeuratProject       5415         2282
AAACAGCCAGGATAAC-1 SeuratProject       2759         1353

scRNAseq
An object of class Seurat 
36601 features across 11909 samples within 1 assay 
Active assay: RNA (36601 features, 0 variable features)
 1 layer present: counts
• 0 views • link
updated by Ram Ram • written by bk11 • bk11

A small educational note: if an answer was helpful, you should upvote it; if the answer resolved your question, you should mark it as accepted. You can accept more than one answer if they work. This will help future users that might find this post find the right answer.

upvote_bookmark_accept

• 0 views • link
written by Ram Ram

Log in to answer this question.

More posts like this