This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How is it possible to have more than 25000 genes in counts file?

Hello everyone,

I have a question about RNA-seq for Homo sapiens. Humans have about 25000 genes which are protein-coding. But in my recent analysis, there are more than 58000 genes with their Ensembl Gene IDs inside my counts file. How is it possible? Did I miss something while processing the fastq's and it caused an error? Thanks...

rna-seq rna-seq sequencing
• 1,704 views • link
updated by Friederike Friederike • written by Batu Batu

2 answers

Not all genes are protein coding.

• 1 views • link
written by Devon Ryan Devon Ryan

and 58,000 is about right for the full gene set including non-coding.

• 1 views • link
written by i.sudbery i.sudbery

How is it possible to have more than 25000 genes in counts file.

If your annotation file (most likely a GTF file?) that you used for generating the counts contained more than 25K gene IDs.

For more details than you might ever want to know about annotation caveats, I recommend reading "A comprehensive evaluation of ensembl, RefSeq, and UCSC annotations in the context of RNA-seq read mapping and gene quantification."

• 1 views • link
written by Friederike Friederike

Log in to answer this question.

More posts like this