This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Tophat Gene Models Merging

Am looking to use multiple human gene models for tophat, like ensembl, refseq and custom gtf files. Can one simply concatenate all these gtf (from the same genome build) and use the final file to guide tophat for rna-seq data?

These would obviously contain exact duplicates and overlaps. One can remove the exact duplicates if that helps tophat performance and efficiency, but overlaps become tricky...

thanks in advance..

tophat rna-seq
• 3,999 views • link
updated by Wen.Huang Wen.Huang • written by Bioinfosm Bioinfosm

1 answer

cufflinks has a utility called "cuffmerge" that does perhaps what you want.

• 1 views • link
written by Wen.Huang Wen.Huang

THanks Wen! But I think cuffmerge is at the latter end of analysis to merge the cufflinks results and annotate. What am looking for is the front end merging of different annotation sources like refseq, ensembl, etc. so I can have all the possible models to guide tophat analysis. If its not necessary to merge and no loss of efficiency, I can simply use a concatenate of all those annotations

• 0 views • link
written by Bioinfosm Bioinfosm

cuffmerge DOES give you a union of all possible models. it takes ANY gtf files. what is better, it filters out redundant ones. but if you don't care about duplicates, then yes concatenation is the simplest way.

• 0 views • link
written by Wen.Huang Wen.Huang

Log in to answer this question.

More posts like this