SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
cufflinks FPKM >>> Cuffdiff FPKM peromhc Bioinformatics 6 10-17-2012 01:07 PM
the gene fpkm and isoform fpkm are zero with Cufflink ??? fulxie RNA Sequencing 5 05-10-2012 12:41 AM
Weird Illumina library ashchin Illumina/Solexa 7 02-17-2012 11:34 AM
VariantFiltration weird error francois.sabot Bioinformatics 0 12-11-2011 11:26 PM
weird 454 reads eoh001 454 Pyrosequencing 2 09-21-2011 06:25 PM

Reply
 
Thread Tools
Old 10-10-2012, 09:53 AM   #1
bioliyezhang
Member
 
Location: Boston, MA

Join Date: Mar 2011
Posts: 19
Question weird fpkm value

Hi, all:

when I checked my paired-end RNAseq data, I saw some weird result.

http://www.freeimagehosting.net/o9akt

basically two genome features have same coordinate(as highlighted in red), however, their fpkm value is different. I wonder anyone have some ideas why this occur?

By the way,
I am using Tophat v2.0.4, Cufflinks v2.0.2 and cummeRbund v2.0.0.

Thanks.

Best,
Liye
bioliyezhang is offline   Reply With Quote
Old 10-11-2012, 02:18 AM   #2
blanco
Member
 
Location: Iceland

Join Date: Apr 2012
Posts: 28
Default

Hi,
I suppose there could be a couple of reasons for this. It could be two different transcripts from the same gene. Or it could be antisense transcription.

If you look at the column gene_short_name and look up those names you will indeed see that the second one is an antisense transcript. What I find a bit strange is when I look up DDX11L10 it is listed on chromosome 16, not chromosome 1. The gene on chromosome 1 at this locus is called DDX11L1.
blanco is offline   Reply With Quote
Old 10-11-2012, 01:04 PM   #3
bioliyezhang
Member
 
Location: Boston, MA

Join Date: Mar 2011
Posts: 19
Default

Hi, Blanco:

Thanks for checking that out.
Yes, the annotation file I used have some potential problem(DDX11L1 gene have one annotation at chr1 and chr16 as well). But I think I know what is the reason.

I check the annotation for DDX11L1 of chr1, which spans 11869 to 13670.
And RP11-34P13.3 spans 29554 to 30039.
And Has-mir-1302-2 starts from 30366 to 30503.

Therefore, I think the locus represents a transcript contig defined by cufflinks, which may include multiple transcripts. And each can have an unique ID. Therefore, both XLOC0001 and XLOC_00002 is within contig chr1:11868-31109.

Best,
Liye
bioliyezhang is offline   Reply With Quote
Reply

Tags
fpkm, rnaseq

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 03:36 AM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO