Hi,
I need to parse human GTF file for my work. I downloaded it from Ensembl.
Basically I don't know what is the number "1" means at the beginning of some lines ?
also if one gene codes more than one transcript, how can I find it ?
Here is the first few lines of it :
1 havana gene 11869 14409 . + . gene_id "ENSG00000223972"; gene_version "5"; gene_name "DDX11L1"; gene_source "havana"; gene_biotype "transcribed_unprocessed_pseudogene";
1 havana transcript 11869 14409 . + . gene_id "ENSG00000223972"; gene_version "5"; transcript_id "ENST00000456328"; transcript_version "2"; gene_name "DDX11L1"; gene_source "havana"; gene_biotype "transcribed_unprocessed_pseudogene"; transcript_name "DDX11L1-002"; transcript_source "havana"; transcript_biotype "processed_transcript";
1 havana exon 11869 12227 . + . gene_id "ENSG00000223972"; gene_version "5"; transcript_id "ENST00000456328"; transcript_version "2"; exon_number "1"; gene_name "DDX11L1"; gene_source "havana"; gene_biotype "transcribed_unprocessed_pseudogene"; transcript_name "DDX11L1-002"; transcript_source "havana"; transcript_biotype "processed_transcript"; exon_id "ENSE00002234944"; exon_version "1";
1 havana exon 12613 12721 . + . gene_id "ENSG00000223972"; gene_version "5"; transcript_id "ENST00000456328"; transcript_version "2"; exon_number "2"; gene_name "DDX11L1"; gene_source "havana"; gene_biotype "transcribed_unprocessed_pseudogene"; transcript_name "DDX11L1-002"; transcript_source "havana"; transcript_biotype "processed_transcript"; exon_id "ENSE00003582793"; exon_version "1";
1 answer
1 is the chromosome (you'll also see 2, 3, 4, ... X, Y, MT and so on). To find genes that code for multiple transcripts, use a regex to extract the gene_id and the transcript_id. If a given gene_id has more than one transcript_id, then you have your answer.
•
written
by
Devon Ryan
There's no reason to bother parsing a GTF file for that, just use Biomart. Just click on the results on that link and save them to a file, since I already created the query for you.
•
written
by
Devon Ryan
Welcome to bioinformatics - we feel your pain, this is your support group
•
written
by
Istvan Albert
N.B., that link in my last comment had the gene_id rather than the name. Just switch that in the "attributes" section and then get the results of that. Mea culpa!
•
written
by
Devon Ryan
I should have refreshed before replying :)
•
written
by
Devon Ryan
When you extract sequences, you get the results in a fasta file. What it does is place the other attributes you wanted (e.g., transcript ID) in the header for each line and it'll separate multiple attributes with a pipe ("|"). That's convenient enough to parse (really easy in biopython) if you really do want things in columns.
•
written
by
Devon Ryan
I see, but how can I find the transcripts of the a given gene that have different 3' UTR sequnces? for example, if a given gene encode 4 transcripts, and 2 of them might have different 3' UTR sequences, I want to put these two transcripts together which have different 3'UTR sequences. I want to do this for all transcriptom.
Either post-process the biomart output to aggregate sequences by gene and compare them or go back to parsing the GTF file and compare the coordinates of everything past the stop codon (be sure to take strand into account!). I expect that the former is easier.
•
written
by
Devon Ryan
- I'm too busy at the moment to write some code to do this for you. You'll need to do that yourself.
- Jein, to steal a word from German. The only reasonable way for the coordinates to be different is if there's some difference in the underlying sequence. While it's possible that a UTR region is repeated exactly and therefore two different transcripts share the same UTR sequence but derive it from different genomic positions, the odds of that occurring are near 0. It be interesting to check this, though.
•
written
by
Devon Ryan
Log in to answer this question.
-
Gene Specific coverage from WGS data
written by Adyasha •
-
Is there a way to safely remove all underscore from gene id for hg19 from genecode
written by simplitia •
-
Why we have exons that start and end at the same coordinate?!
written by Medhat
-
9th column of the GTF file is exon_id instead of gene_id, is this normal?
written by Kai_Qi
-
Annotating an assembled GTF file
written by nattzy94 •
-
Problem trying to convert from .gtf to .fasta file
written by nattzy94 •
-
Can I get the same annotation data that is in ensembl gene annotation sets through mysql?
written by endrebak
-
RSEM with ERCC spike-ins formatting issues, ERCC sequences never make it into the reference
written by brian.capaldo •
-
Problems uploading custom track to UCSC Genome Browser
written by rs •
-
Choice of Gene Length in TPM Calculations
written by dec986
