Re: Cuffdiff Using The Combined Gtf Files From Cuffcompare

Question: Re: Cuffdiff Using The Combined Gtf Files From Cuffcompare

8.0 years ago by

United States

Jennifer Hillman Jackson ♦ 25k wrote:

Hello David, You are correct about the tools, so the problem is most likely with the original GTF file. If gene_id is not assigned there correctly, then the data will not be sorted by gene_id. Although GTF format is consistent (mostly!) between sources, the actual content can vary. One example is from UCSC - the GTF format from the Table browser will have the transcript name assigned to both the gene_id and the transcript_id tags in the attributes field (f9). Post processing to extract gene name from the track and swapping it into the GTF file's gene_id attribute tag would be a necessary pre-processing step before using the downstream tools with functionality that would use the attribute. The good news is that you should be able to use Galaxy's Text Manipulation tools to do whatever file processing you need to do, from whatever input source you are using, once you have the data content loaded into your history. Create->save->use a workflow so that you only have to work out tedious file conversions step-by-step one time. If you need more help, please let us know and share your history: Options -> Share or Publish -> Share with a user "jen@bx.pus.edu". Thanks, Jen Galaxy team ps. It is best to send data questions to galaxy-use mail list to help the community learn from each other. I am going to forward this answer there now, since this question has come up a few times recently after the addition of the new tools.

galaxy • 854 views

ADD COMMENT • link •

modified 8.0 years ago by David Matthews • 630 • written 8.0 years ago by Jennifer Hillman Jackson ♦ 25k

8.0 years ago by

David Matthews • 630

United Kingdom

David Matthews • 630 wrote:

Hi, Sorry about not replying via the user group mail list - lack of concentration on my part. I also emailed Adam Roberts about this problem of multiple identical gene names and his reply is below: So it seems as though it is probably/possibly a problem of when are two transcripts really the same gene (!)? Not always a straightforward question. However, for simple summing of expression to get a simple answer I used the "group" command. I take the isoforms.fpkm file and group on gene name and sum the fpkm values - that gives me a straightforward list of expressed genes with a single number attached. However, I now have two more problems (!!). I have a time course of samples T0, T8 and T24 and so I compare T0 with T8 and T8 with T24. However, the expression values for T8 are not the same - the T0/T8 values are slightly higher for T8 than the T8/T24 samples. I am using the latest ensembl.gtf (I fixed the chr issue) and I made my own female_hg19 file (the data is from HeLa cells). Again, I've emailed Adam about this to see what he thinks. I ran the analysis twice to (hopefully) make sure I did not set it up wrong. So I have installed the latest cufflinks here on my mac and I'm running it again but with all three at the same time - it will take a day or two to complete. My second and hopefully final problem is that the p_id and cds files from cuffdiff are empty - can't figure this one out. I've tried looking on seq answers but with no avail although others have seen the same thing. The ensembl.gtf file is mentioned in the cufflinks manual so I assume its the best one to use but maybe I did something unintended when I added the "chr" to the chromosome numbers. I added you as a user so you can see for yourself, any suggestions gratefully received. Best Wishes, David. P.S. Just love the site - excellent stuff - bored everyone to tears here about how great it is! __________________________________ Dr David A. Matthews Senior Lecturer in Virology Room E49 Department of Cellular and Molecular Medicine, School of Medical Sciences University Walk, University of Bristol Bristol. BS8 1TD U.K. Tel. +44 117 3312058 D.A.Matthews@bristol.ac.uk

ADD COMMENT • link written 8.0 years ago by David Matthews • 630

Similar posts • Search »