Question: Differential expression of scRNA-seq after Seurat
0
gravatar for lshepard
17 months ago by
lshepard220
United States
lshepard220 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 • 2.9k views
ADD COMMENTlink modified 16 months ago by igor7.6k • written 17 months ago by lshepard220
1
gravatar for igor
16 months ago by
igor7.6k
United States
igor7.6k 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 16 months ago • written 16 months ago by igor7.6k

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 16 months ago • written 16 months ago by lshepard220
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 16 months ago by igor7.6k

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 16 months ago by lshepard220
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: 1881 users visited in the last hour