This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Cutoff of reads in htseq

Hello,

To my understanding the standard number of cutoff for falsely expressed genes in trancriptome is 30 reads.

Am I right? and what is the reason for a this specific cutoff?

Thank you

rna-seq
• 3,437 views • link
updated by Ram Ram • written by olga.lubinsky • olga.lubinsky

Can you explain more ? As the heading of question, htseq, there is no min or max cutoff to calculate no. Of read/fragments mapped per genomic location. If you are talking about the differential expression, the cutoff will be based on statistical test but not on any standard number.

If you can provide any source for your understanding, that would be good.

• 0 views • link
written by GouthamAtla GouthamAtla

2 answers

There is no standard cut-off and there could never be one. The best practice is to use independent filtering.

• 0 views • link
written by Devon Ryan Devon Ryan
there is no standard threshold. you can use the one that is more suitable for your data. in deseq2 you will find an implementation of the independent filtering, suggested by Devon. it basically sets the threshold in a way to get the highest number of differentially expressed calls. check it afterwards, since it could be too high (or too low), depending on your expectations.
• 0 views • link
written by Martombo Martombo

I'm not dealing with differentially expressed gene but with the transcriptome itself.

I have the reads of each gene , how do I decide which number of reads is too low to consider as a signficantly expressed gene?

• 0 views • link
updated by Ram Ram • written by olga.lubinsky • olga.lubinsky

There's no standard for that as well (there never could be). The best you can do is plot the FPKM/RPKMs and see if you have an obviously bimodal distribution. If not, you're left with using the top X% as expressed.

• 0 views • link
written by Devon Ryan Devon Ryan

Hello Devon, can i ask you after the observing the bimodal distribution, which steps do you suggest for filtering, im dealing with tpm. Thanks in advance

• 0 views • link
written by lessismore lessismore

Log in to answer this question.

More posts like this