Hello,
I am doing the single cell RNA-seq data analyzing now. I have a treated sample and a non-treated sample. To compare the differences, I have combined these two samples using the official integrating tutorial. My cell populations vary a lot between two samples, so after clustering, in some clusters, there are only very few cells in one sample but lots of cells in other sample.
When I do the Differential expression testing using FindMarkers between these two samples in one cluster, the top genes showed in the results only express in several cells in my many-cell samples, but I want to see the genes specifically broadly expressed in my many-cell samples.
I think the problem is the logic of FindMarkers is: Only the genes expressed in both samples can be compared. But some of my genes may not even emerge in one of my sample. So, are there any way to identify the genes only expressed in one sample in my combined file?
Thank you so much for your help!
1 answer
I think the problem is the logic of FindMarkers is: Only the genes expressed in both samples can be compared
Where are you getting this information from? This is not true as far as I'm aware. If you want to ensure any genes found are more broadly expressed in one of your populations, you can try increasing min.pct parameter:
min.pct: only test genes that are detected in a minimum fraction of min.pct cells in either of the two populations. Meant to speed up the function by not testing genes that are very infrequently expressed. Default is 0.1
Log in to answer this question.