I have an issue with my RNA-seq analysis, and was wondering whether any of you know how to solve it? Basically, I want to be able to annotate my analysis with gene symbols. What I have done is to supply Tophat/Bowtie with a gtf file from Ensembl, which contains these gene symbols. However, instead of making use of these gene symbols (the column "gene_name" in the gtf file, for example "DDX11L1"), it seems to use the Ensembl gene ID's instead (the column "gene_id" in the gtf file, for example "ENSG00000223972").
How do I get Tophat/Bowtie to use the "gene_name" column instead? Is this possible?
1 answer
It seems I had been misunderstanding the setup of this for a while now. Precisely because each individual gene_name links to multiple gene_id's (because of for example multiple transcripts of the same gene), this cannot be done directly. It's better to rather annotate the transcripts to gene symbols at a later point during the downstream analysis, I've learnt.
Log in to answer this question.
-
bedmap to annotate bed file with GTF that was converted to BED produces way too many results
written by Tawny
-
9th column of the GTF file is exon_id instead of gene_id, is this normal?
written by Kai_Qi
-
Error with Rseqc infer_experiment - unknown datatype
written by newbie
-
Annotating an assembled GTF file
written by nattzy94 •
-
Problem trying to convert from .gtf to .fasta file
written by nattzy94 •
-
RSEM with ERCC spike-ins formatting issues, ERCC sequences never make it into the reference
written by brian.capaldo •
-
Ensembl GTF format: isn't the tag "transcript_id" mandatory?
written by Pfs
-
Second column in Ensembl GTF: source or biotype?
written by alex.rubinsteyn
-
What's the number at the begining of some lines of transcriptome GTF file means ?
written by jack
-
Tophat / Cufflinks Not Using Gene_Id Or Gene_Name From Gtf File
written by Owen S.
