Read Counts Difference From Different Tools
0
1
Entering edit mode
11.7 years ago
camelbbs ▴ 710

Hi, I got reads counts from rnaseq from different tools, I am not sure why the cufflinks gives so different counts result. Anyone have experience on that?

    ACTB    40264.3    38071.2    26648.7    41528.2    39006.5    40361.8    47723.2    25179.8   cufflinks
    ACTB    2889    1102    1793    2264    2747    1836    624    2316   htseq
    ACTB    2084    875    1362    1719    1919    1345    475    1722    partek

Thanks, Cam

rnaseq • 3.3k views
ADD COMMENT
0
Entering edit mode

Check this post Deseq, Edger And Cuffdiff - Different Result, might give some hints :)

ADD REPLY
0
Entering edit mode

Is that the correct link? I do not see that phrase, nor anything which answers this question.

ADD REPLY
0
Entering edit mode

now corrected !!

ADD REPLY
0
Entering edit mode

I think it was about differential expression, and that the problems were eventually with cuffdiff, not cufflinks. I think Cam's problem is at the read count estimation stage.

ADD REPLY
0
Entering edit mode

You can check these links:

How did you run cufflinks? Its estimates are orders of magnitudes higher than the other tools.

ADD REPLY
0
Entering edit mode

Thanks, I run cuffdiff like this: "cuffdiff -p 8 --uper-quartile-norm -u -b hg19.fa -o output merged.gtf s1.bam s2.bam s3.bam ..."

ADD REPLY
0
Entering edit mode

Cam, could you try to run cuffdiff without the quartile normalizatioin?

ADD REPLY
0
Entering edit mode

Hi massyah, I tried with no --uper-quartile-norm option and got the same result... So I think the problem might not be here..

ADD REPLY
0
Entering edit mode

Could it be that cufflinks returns FPKM and htseq and partek return the other two read counts? I just read the htseq man page here http://www-huber.embl.de/users/anders/HTSeq/doc/count.html#count (I never used the software so ... beware!) and I think I understood that the software actually count the reads mapping, and not FPKM. I don't know what partek does. You can track read counts in cuffdiff output files such as "genes.count_tracking", but I never didi it.

ADD REPLY

Login before adding your answer.

Traffic: 1386 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