Hi, all
Recently, I am working with ensembl GTF annotation files, and try to detect the useful exon I wanted.
I am confused with ensembl Exon ID. For example, the three exons (see below) are belongs to Gene ENSG00000000003 and have the same start site and end site.
chrX protein_coding exon 99890555 99890743 . - . gene_id "ENSG00000000003"; transcript_id "ENST00000373020"; exon_number "2"; gene_name "TSPAN6"; gene_biotype "protein_coding"; transcript_name "TSPAN6-001"; exon_id "ENSE00003662440";
chrX processed_transcript exon 99890555 99890743 . - . gene_id "ENSG00000000003"; transcript_id "ENST00000496771"; exon_number "2"; gene_name "TSPAN6"; gene_biotype "protein_coding"; transcript_name "TSPAN6-003"; exon_id "ENSE00003512331";
chrX processed_transcript exon 99890555 99890743 . - . gene_id "ENSG00000000003"; transcript_id "ENST00000494424"; exon_number "3"; gene_name "TSPAN6"; gene_biotype "protein_coding"; transcript_name "TSPAN6-002"; exon_id "ENSE00003512331";
My questions:
- Why the first exon (ENSE00003662440) and last two exons (ENSE00003512331) are annotated with different Exon ID?
- Could anybody explain the method of Exon ID annotation? (I don't find any document on ensembl site about the Exon annotation)
Thanks
1 answer
I'd guess part of the explanation is that the transcripts come from different sources and, despite what your GTF file states, have different biotypes.
Take a look at the region here. Transcripts TSPAN6-002 and TSPAN6-003 are Havana transcripts of type "processed transcript". TSPAN6-001 is an Ensembl/Havana merge transcript of type "known protein coding". So the 2 former exons are considered "the same exon"; the latter exon has the same coordinates but a different source so is considered a different exon.
There is some (not detailed) information about the annotation process here. Also note that your data appear to come from genome build GRCh37 and things are somewhat different in the latest build. There are now 5 transcripts of 3 types and consequently 3 exon IDs.
Log in to answer this question.
-
Why we have exons that start and end at the same coordinate?!
written by Medhat
-
How to add "transcript" feature to a gtf file?
written by pomodoro_sinensis •
-
Gene feature information missing in Stringtie merged assembly
written by waqaskhokhar999
-
Can I get the same annotation data that is in ensembl gene annotation sets through mysql?
written by endrebak
-
modify gff file
written by manishbiotechie •
-
How can I count read that falls into exon-intron and exon-exon junction?
written by biplab
-
What is agi and bgi in the Ensembl gtf file
written by lakhujanivijay
-
How to had CDS feature with cufflinks (cuffmerge) using custom GTF?
written by A. Domingues
-
What's the number at the begining of some lines of transcriptome GTF file means ?
written by jack
-
Ensembl And Ucsc Gtfs To Match
written by dfernan
I'm also a bit confused about this GTF files. What does it mean "exon version"? What means the first "1" on every column?
Thanks in advance!
first "1" on every column is name of the chromosome ( chromosome 1 ). exon_version: The stable identifier version for this exon.
you can find gtf format detail from: ftp://ftp.ensembl.org/pub/release-81/gtf/homo_sapiens/README