Question: Differential expression of scRNA-seq after Seurat
0
gravatar for lshepard
2.8 years ago by
lshepard420
United States
lshepard420 wrote:

Hi,

I am currently planning a scRNA-seq analysis and I would like to have some feedback from this community. As this is a Drop-seq run, I will be using the Seurat package to identify clusters of cells. However, the program does not perform differential expression among experimental conditions (only among populations of cells).

From what I gathered, bulk RNA-seq methods have been shown to perform well (DESeq2, EdgeR), along with scRNA methods like MAST. Recently, Seurat has included DESeq2 and MAST as part of the tests of expression for cell clusters. I am leaning towards DESeq2 for my approach as I like how the design can be constructed, however, I wanted to see if this sounds reasonable among users in here. As a further question, Seurat performs global normalization in the data, which confuses me a bit in the need to re-normalize in DESeq2 (although I think it is necessary as it normalizes for sequencing depth - thoughts?).

Thanks for the input.

rna-seq • 6.1k views
ADD COMMENTlink modified 2.8 years ago by igor11k • written 2.8 years ago by lshepard420
1
gravatar for igor
2.8 years ago by
igor11k
United States
igor11k wrote:

However, the program does not perform differential expression among experimental conditions (only among populations of cells).

You can set cell identities to anything. Then, differential expression can be performed between any of the identities.

You can run DESeq2 or MAST as part of Seurat's FindMarkers() function. That way, you don't have to worry about moving the data between different packages manually.

Check the Seurat differential expression vignette: http://satijalab.org/seurat/de_vignette.html

ADD COMMENTlink modified 2.8 years ago • written 2.8 years ago by igor11k

I edited my reply because your answer made me think of the following:

Since you need to keep everything within a single object, I could merge the objects (control vs treatment) and follow this up with FindMarkers correct?

ADD REPLYlink modified 2.8 years ago • written 2.8 years ago by lshepard420
2

I didn't realize your samples were in separate objects.

You can merge multiple Seurat objects into a single Seurat object. All the annotations in the meta.data slot should be preserved.

ADD REPLYlink written 2.8 years ago by igor11k

Thanks, I don't know why I was thinking it was better to keep the treatment groups as separate objects. It makes more sense now to merge them. I appreciate the input!

ADD REPLYlink written 2.8 years ago by lshepard420
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 829 users visited in the last hour