How does one get gene names when using cuffdiff when looking at "gene differential expression testing" results? I am doing some +/- exposure studies with zebrafish embryos and then processing the 50bp single-ended fastq files through Galaxy with particular interest in the cuffdiff readout of differential gene expression. It seems to be working well but my readout is only giving me gene ID's. It would by more efficient if I could get gene names from the output. I suspect I may be able to do this through selection of options from the UCSC TABLE BROWSER when I use its Zv9 danRer7 assembly as a reference genome. I see a place for selected genes but not for all the identified genes. Is this something that could be done directly or is there just a "beyond Galaxy" method of gene gene lists from a list of gene ID's. Thanks, el linney Duke University Medical Center
Hello, If you used RefSeq as your reference annotation GTF from UCSC, the two tables from the UCSC Table browser that you will most likely be initially interested in are refGene and refLink. These can be extracted entirely into Galaxy - as separate data files or together as one file, using the output option "all fields from selected table" or "selected fields from primary and related tables". BED or GTF format will not output the extra fields containing gene symbols, names, descriptions - tip: click on 'describe table schema' to understand a table's contents. Then a tool such as "Join, Subtract and Group -> Join" or "Compare" can be used with the Cuffdiff result's transcript accession and this data to add in the additional labels. If you used a different annotation GTF than RefSeq as your reference, use that track's associated tables from UCSC or Biomart instead. If you need to link genes using overlap - a method similar to the one described in this prior reply could can be used: http://user.list.galaxyproject.org/Genes-from-Regions-td4136535.html Hopefully this helps! Jen Galaxy team On 1/7/13 7:58 AM, Elwood Linney wrote:
How does one get gene names when using cuffdiff when looking at "gene differential expression testing" results?
I am doing some +/- exposure studies with zebrafish embryos and then processing the 50bp single-ended fastq files through Galaxy with particular interest in the cuffdiff readout of differential gene expression. It seems to be working well but my readout is only giving me gene ID's. It would by more efficient if I could get gene names from the output.
I suspect I may be able to do this through selection of options from the UCSC TABLE BROWSER when I use its Zv9 danRer7 assembly as a reference genome. I see a place for selected genes but not for all the identified genes.
Is this something that could be done directly or is there just a "beyond Galaxy" method of gene gene lists from a list of gene ID's.
Thanks,
el linney Duke University Medical Center
___________________________________________________________ The Galaxy User list should be used for the discussion of Galaxy analysis and other features on the public server at usegalaxy.org. Please keep all replies on the list by using "reply all" in your mail client. For discussion of local Galaxy instances and the Galaxy source code, please use the Galaxy Development list:
http://lists.bx.psu.edu/listinfo/galaxy-dev
To manage your subscriptions to this and other Galaxy lists, please use the interface at:
-- Jennifer Hillman-Jackson Galaxy Support and Training http://galaxyproject.org
participants (2)
-
Elwood Linney -
Jennifer Hillman-Jackson