Question: Feature selection for differential analysis
gravatar for wo_li
24 days ago by
wo_li10 wrote:

Hello all,

I'm doing the single cell RNA-Seq and trying to find out genes showing different expression under two treatments. Thus my first step is to put two datasets from different treatments together and the second is to compare the difference.

When I merge two single cell RNA-seq datasets together, I notice that there are some unique genes from the previous two datasets and they are maintained by Seurat::merge(). I'm wondering whether I should keep these unique genes? Will they affect the downstream analysis of differential analysis?

Thank you in advance!

ADD COMMENTlink modified 23 days ago by kristoffer.vittingseerup2.3k • written 24 days ago by wo_li10

You should not keep the unique genes. If there is only a small number of them then you might be OK just excluding them. If it is a huge proportion of genes then what @kristoffer.vittingseerup suggested below is probably wise.

ADD REPLYlink written 23 days ago by benformatics1.1k

Thank you and agree! It's only a small proportion and their expression are also low in the other dataset.

ADD REPLYlink written 18 days ago by wo_li10
gravatar for kristoffer.vittingseerup
23 days ago by
European Union
kristoffer.vittingseerup2.3k wrote:

There should not be unique genes - then it is because your filtering step(s) have removed them (meaning they are candidate for being highly expressed in one set and not the other). I would suggest merging the data before you start the Seurat analysis so you have a joint QC, filtering etc.

Please also be aware of batch effects - if you ran the single-cell data as two different rounds you might have a large batch effect - then you need to use approaches such as Seurat's data integration.

ADD COMMENTlink written 23 days ago by kristoffer.vittingseerup2.3k

Thank you, Kristoffer! Yes, these genes are filtered first as they have no expression under that treatment. As you suggest, I put them together to do the QC.

ADD REPLYlink written 18 days ago by wo_li10
Please log in to add an answer.


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