How to deal with single nucleus data in Seurat
Entering edit mode
2.1 years ago

Dear all, I am by no means a bioinformatician, but because of the pandemic I have been looking for parts of my work I can do from home, so please forgive me if I ask extremely basic questions.

I obtained a data set of single nuclei from the Drosophila antenna, and my boss suggested using the Seurat pipeline to analyse the data. I did the Seurat tutorials but none of them are specific for nuclei.

Does anyone have experience how to deal with this kind of data or can point me into the direction of some more information? For example in the QC step mitochondrial genes are excluded but my nuclei data should not have any of those to begin with.

Thank you for your advice, any help is appreciated.

RNA-Seq seurat single nucleus RNA sequencing • 1.9k views
Entering edit mode

You can proceed as normal with single nuclei data.

Entering edit mode
21 months ago
Pratik ▴ 850

Hi Kristina,

Rpolicastro is right.

There should be no difference in either using single cell RNA seq data or single nuclear RNA seq data with Seurat.

A good place to start is definitely:

Many of the steps are modifiable, based on your judgement. If you know that your data doesn't have any mitochondrial gene contamination, than you could definitely just skip that step.

pbmc[[""]] <- PercentageFeatureSet(pbmc, pattern = "^MT-")
VlnPlot(pbmc, features = c("nFeature_RNA", "nCount_RNA", ""), ncol = 3)

or the following which I use cause because my expression matrix lists the mitochondrial genes as MT... rather than MT-...

pbmc[[""]] <- PercentageFeatureSet(pbmc, pattern = "^MT")
VlnPlot(pbmc, features = c("nFeature_RNA", "nCount_RNA", ""), ncol = 3)

just shows the mitochondrial genes.

The subset step where you do:

pbmc <- subset(pbmc, subset = nFeature_RNA > 200 & nFeature_RNA < 2500 & < 5)

You can certainly just remove the & < 5 (MT gene filtering) to just pbmc <- subset(pbmc, subset = nFeature_RNA > 200 & nFeature_RNA < 2500). Furthermore, as you may know the filtering numbers are also up to you to decide where you want to draw the line in QC, based on your analysis of the previous plots generated.

Hope this helps!

Feel free to ask any other questions. I think I've got a decent grasp on Seurat. : )

Sincerely, Pratik


Login before adding your answer.

Traffic: 2025 users visited in the last hour
Help About
Access RSS

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6