Question: Re: Cuffdiff Using The Combined Gtf Files From Cuffcompare
gravatar for Jennifer Hillman Jackson
8.0 years ago by
United States
Jennifer Hillman Jackson25k 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 "". 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 COMMENTlink modified 8.0 years ago by David Matthews630 • written 8.0 years ago by Jennifer Hillman Jackson25k
gravatar for David Matthews
8.0 years ago by
United Kingdom
David Matthews630 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
ADD COMMENTlink written 8.0 years ago by David Matthews630
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 180 users visited in the last hour