Hello, I have to make a Java program for a college course in which I have to find possible intron retention in a given sample.
I am stuck in the initial part where, given a reference GTF file, I have to parse it and recover all intron regions from in (Making another, pruned, GTF file).
I am not getting how could I find where an intron starts and ends Thanks
1 answer
Hi Alessandro,
Here are a couple relevant links that may be useful:
If you can use UCSC, rather than be forced to use a GTF file: Bed File With Introns Only
A little more detail on a manual solution is here: https://biostar.usegalaxy.org/p/6453/
And some basic information on GTF parsing in Python: http://biopython.org/wiki/GFF_Parsing
It's tough to give an exact solution without seeing the file, but hopefully this is useful. Good luck!
Brent Wilson, PhD | Project Scientist | Cofactor Genomics
4044 Clayton Ave. | St. Louis, MO 63110 | tel. 314.531.4647
Catch the latest from Cofactor on our blog.
Log in to answer this question.
-
Export sequence from specific region from bam file.
written by malkomalktho •
-
Genomic position to genomic feature
written by irfanwustl •
-
How to annotate bam-file using ready annotation gtf
written by Alexandr •
-
Identify reads that span intron-exon junctions in RNA-Seq
written by oriolebaltimore
-
Intron cordinates requires for known exons
written by HK •
-
Retained intron transcripts
written by mforde84
-
Genomic feature graph
written by jb.sequence •
-
Extracting the donor/acceptor sites from mRNA sequence(PERL/ DNA transcription) ?
written by reeki •
-
From Where In The Mzml File Does Toppview Get The Retention Time?
written by Niek De Klein
-
Intron Retention And Alternative 5'3' Splice Site Identification
written by Arun
Are you using a specific GTF file? Can you post a few lines? Latest GTFspec is available at this link.
Hi thank you for the help. But I forgot to mention (I think it's quite important) that my input file contains ALL and ONLY the known exons of a human genome sample. I was thinking on computing from that file for each chromosome and then for each gene and for each transcripts (I saw the same gene can have multiple versions of itself due to splicing events) where is the exons start and end. And then compute the introns as the complementary of this.
How do you think about my algorithm?
Thanks
Please use
ADD REPLY/ADD COMMENTwhen responding to existing posts.Introns are not complementary (not sure what sense you are saying that in). They represent the interval between two exons. e.g.
Exon_1-Intron_1-Exon_2-Intron_2-Exon_3etc. More here: https://en.wikipedia.org/wiki/ExonAlso see this thread for a nice graphic: What'S The Difference Between Cds And Orf?