SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
multiple lines per single gene in gene_exp.diff output of cuffdiff reut Bioinformatics 2 03-09-2015 04:16 AM
Cufflinks (Cuffdiff) v0.9.3 multiple samples keebs42 Bioinformatics 4 11-06-2013 02:20 AM
how to get gene name from accession ID_in batch ahli1981 Bioinformatics 0 12-19-2011 12:32 PM
multiple FPKM problem for single gene in gene_exp.diff after running cuffdiff ngs RNA Sequencing 4 03-30-2011 01:55 PM
Multiple gene entries in cuffdiff output ae_ucla RNA Sequencing 1 10-22-2010 03:54 PM

Reply
 
Thread Tools
Old 12-13-2011, 06:59 AM   #1
apadr007
Member
 
Location: washington DC

Join Date: Oct 2011
Posts: 21
Default Cuffdiff multiple accession numbers

Hello,

I wanted to know why the Cuffdiff "gene.exp.diff" and "cds.diff" files were producing multiple accession numbers in the same row. For example:


XLOC_000122 - - XLOC_000122 AGAP004916,AGAP004918 TSS334,TSS335



As you can tell, this is a single row with two "AGAP" accession numbers, separated by a comma. Why is this happening?

Thanks.
apadr007 is offline   Reply With Quote
Old 12-13-2011, 10:29 AM   #2
NM_010117
Member
 
Location: US

Join Date: Oct 2010
Posts: 14
Default

This is most likely due to it reporting the different isoforms that are part of that gene. Most mammalian genes have at least two isoforms, if not three or four (or more!). That is why you are also getting two different transcription start sites (TSS): TSS334 and TSS335. One is for AGAP004916 and the other is the transcription start site for AGAP004918.

Last edited by NM_010117; 12-13-2011 at 10:30 AM. Reason: grammar
NM_010117 is offline   Reply With Quote
Old 12-13-2011, 12:40 PM   #3
apadr007
Member
 
Location: washington DC

Join Date: Oct 2011
Posts: 21
Default

This makes absolutely no sense to me. Cuffdiff generates an isoform file so it should simply report isoforms of a particular gene there, not in a file it calls "gene.exp.diff".

I guess it may do with the way isoforms are clumped up into genes. Maybe the algorithm misses some that are very different and therefore calls it a new gene?
apadr007 is offline   Reply With Quote
Reply

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:35 AM.


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