Question

Seurat DefaultAssay integrated vs RNA

0

Entering edit mode

4.3 years ago

Payal ▴ 160

Hi,

Whats the difference between the following in Seurat objects?

DefaultAssay(combined.data) <- "integrated"

&

DefaultAssay(combined.data) <- "RNA"

Thanks

RNA-Seq sequencing R • 6.8k views

ADD COMMENT • link 4.3 years ago by Payal ▴ 160

0

Entering edit mode

Can I ask I followup question -

I need to first find markers for clusters for integrated dataset and the do Differential gene expression per cluster per condition. For that should I choose DefaultAssay(combined.data) <- "integrated" or DefaultAssay(combined.data) <- "RNA". How will the results differ?

ADD REPLY • link 4.3 years ago by Payal ▴ 160

0

Entering edit mode

Please use ADD COMMENT/ADD REPLY when responding to existing posts to keep threads logically organized. SUBMIT ANSWER is for new answers to original question

ADD REPLY • link 4.3 years ago by GenoMax 147k

1

Entering edit mode

In this vignette it is suggested DefaultAssay(combined.data) <- "RNA" to find markers.

António

ADD REPLY • link 4.3 years ago by antonioggsousa 3.2k

score 3 · Answer 1 · 2020-08-20

The RNA assay contains the raw counts, and if you use their older count normalization method (not SCTransform), the normalized and scaled counts. The integration assay is created after normalization and integration, as detailed in their integration vignette. The counts here are slightly adjusted so that cells that are (probably) similar between samples and technologies will cluster closer together with PCA and UMAP dimension reduction.

score 1 · Answer 2 · 2020-08-20

Hi,

You're setting your default assay as "integrated" or "RNA". A Seurat object has many data sets, usually the "RNA" data slot contains multiple data sets: https://github.com/satijalab/seurat/wiki/Assay

I do not remember quite well, but I believe that If you performed integration, apart from the "RNA", you'll have also "integrated" slot, with possible multiple data sets, like gene expressions tables normalized, transformed etc.

When you do - DefaultAssay(combined.data) <- "integrated" - you're setting the default assay "integrated". This means when you call other Seurat functions they will use this "integrated" slot data. If you want to use the "RNA" instead you do: DefaultAssay(combined.data) <- "RNA"

I hope this answers your question.

António