This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Using Gencode annotation with cufflinks

Hello,

I am experimenting with the latest Gencode human annotation releases (v20 and v21) with my RNA-seq pipeline. After filtering the annotation file for known and level 1 and level 2 transcripts, I find that cufflinks reports multiple FPKM values for some ENSG genes in the gene.fpkm_tracking file. A manual inspection suggests that the multiple FPKM values correspond to distinct transcripts for these genes, but these values should ideally get rolled up into the gene level value (as it does for a majority of the other multi-transcript genes).

Has anyone else noticed this, and is there a secret sauce to working with the Gencode annotation GTF files with cufflinks?

rna-seq cufflinks
• 3,087 views • link
updated by Ram Ram • written by drollix • drollix

1 answer

If you experience problem with V20, try this on your GTF and check your results (I had problem saying that duplicate GeneID, in your case cufflinks did not throw it as duplicate but instead it has two different expression values).

sed -i 's/\"\ tag\ \"PAR\"\;/\"\;\ tag\ \"PAR\"\;/g' gencode.v20.annotation.gtf
• 0 views • link
updated by Ram Ram • written by EagleEye EagleEye

Thanks, I'll give it a try, though the genes with duplicated FPKM's were not on chromosome Y. I also didn't find any pattern (yet) about the transcript structures of the multi-FPKM genes. Some had disjoint transcripts while some had completely overlapping ones.

• 0 views • link
updated by Ram Ram • written by drollix • drollix

Log in to answer this question.

More posts like this