This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Seurat DefaultAssay integrated vs RNA

Hi,

Whats the difference between the following in Seurat objects?

DefaultAssay(combined.data) <- "integrated"

&

DefaultAssay(combined.data) <- "RNA"

Thanks

rna-seq sequencing r

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?

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

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

António

2 answers

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.

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

Log in to answer this question.