How featureCounts define the gene length
1
1
Entering edit mode
Yiwei Niu ▴ 20
@niuyw-11650
Last seen 2.1 years ago
Beijing, China

Hi, I'm using featureCounts from the Rsubread package. But I have a question about the gene length returned by featrureCounts.
I've read the case study here: http://bioinf.wehi.edu.au/RNAseqCaseStudy/, and the combination of featureCounts and edgeR is wonderful.
We often need to get the R/FPKM values from a RNAseq experiment, and I just used the rpkm() function in edgeR. This function takes the gene length as input, which I got by featureCounts. So I wonder how featureCOunts define the gene length? This will affect the result. I've tried google but got little information.

Any help will be appreciated,thank you very much!

featurecounts rnaseq rsubread • 6.0k views
ADD COMMENT
3
Entering edit mode
@gordon-smyth
Last seen 55 minutes ago
WEHI, Melbourne, Australia

It is the total number of DNA bases that belong to any exon for that gene.

If you type ?featureCounts to get help on featureCounts, then read the description of the output value, it says:

"For each meta-feature, the Length column gives the total length of genomic regions covered by features included in that meta-feature. Note that this length will be less than the sum of lengths of features included in the meta-feature when there are features overlapping with each other."

ADD COMMENT
1
Entering edit mode

Thanks! According the doc of featureCounts, it's the total number of intersection of exons that belong to the gene, not the total number of DNA bases that belong to any exon for that gene.

ADD REPLY
0
Entering edit mode

I have to disagree with you -- the featureCounts documentation does not make any reference to an "intersection of exons". As far I can remember, gene length has always been correctly described in the featureCounts documentation and agrees with my answer. The quote I give from the documentation demonstrates this.

ADD REPLY
0
Entering edit mode

Hi, I'm quite new to RNA seq analysis. Do you know if the length column in the featureCounts output is in kbp?

Thanks!!

ADD REPLY
0
Entering edit mode

Rather than adding a comment to a 5-years old (!) thread, please start a new one.

ADD REPLY
0
Entering edit mode

The length column is exactly as described in the documentation and is literally the number of bases, not kbp.

ADD REPLY

Login before adding your answer.

Traffic: 371 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6