Hello everyone,
I used DESeq2 to do differentially gene expression analysis for 6 RNA-seq samples per species.
tximport was used to input salmon mapping results to DESeq2. I measured the expression value (txi.salmon$counts) for the differentially expressed genes per sample per species, and found that it is lower than that not suggested as differentially expressed genes. For example, the average expression value for the differentially expressed genes is 100, while the value is 300 for the gene not suggested as differentially expressed.
I have tested 4 species, and found same problem. Could anybody explain or give suggestion?
Thank you
Chen Lingyun
Hi Michael,
Thanks for your great answer. It really helps me.
Is there any evidence/opinion shows that genes with overall low expression level across samples are more likely to be incorrectly identified as differentially expressed genes? In other words, genes with high expression level are less likely to be identified as differentially expressed genes?
Thank you very much.
Best wishes,
Chen
"Is there any evidence/opinion shows that genes with overall low expression level across samples are more likely to be incorrectly identified as differentially expressed genes? "
I don't think so. The methods take into account of the amount of precision from low and high counts genes. I'd recommend you look at plotCounts() for individual genes.
"In other words, genes with high expression level are less likely to be identified as differentially expressed genes?"
I think that the true DE genes can be skewed low, or highly expressed, or scatter uniformly. This is more about the biology.
Thank you so much. I will use plotCounts to see what happen