This is a test version of Biostars. For the public version, visit https://www.biostars.org.
how to fetch intron coordinates from transcriptome.gtf file ? is there any direct tool ?

Hello

Is there any direct tool to fetch intron coordinates from transcriptome.gtf file ?

Is there any direct tool available ?

Thank you

intron
• 3,694 views • link
updated by Ram Ram • written by isha.lily20 • isha.lily20

It’s simple to parse. Simply take the coordinates annotated as .`gene` transcript and those annotated as exon. Take the complement, so gene transcript minus exon, that’s the introns.

If you google the question there are multiple threads on that here at biostars already.

• 0 views • link
written by ATpoint ATpoint

If you want the exact intron ranges you would need to modify the above by looking per transcript. Doing it per gene would give you regions contained within the introns of all transcripts, but not proper introns per say.

• 0 views • link
written by rpolicastro rpolicastro

Indeed, good catch!

• 0 views • link
written by ATpoint ATpoint

i have done intergenic minus exon for intron coordinates, that is wrong ?

• 0 views • link
written by isha.lily20 • isha.lily20

In a sense, yes it is wrong because it doesn't capture biological complexity of alternative splicing. If you have simple gene models where each gene has a single annotated transcript, then it would suffice, but only because the gene models do not capture the real biology well.

Imagine these hypothetical transcripts for the same gene:

T1: oo     t1e1]--------[        t1e2]----------[  t1e3  ooo>
T2: ooo    t2e1]--------[    t2e2]-----[ t2e3]--[ t2e4 oo>

 oo: UTR's or leader sequences
 [  ]: internal exon boundaries
 --: intronic sequence

If you look at genes in a genome browser you will find even more complex cases. Therefore, it becomes very clear that introns depend on the actual transcript.

• 0 views • link
written by Michael Dondrup Michael Dondrup

1 answer

If GFF/GTF output format is ok you can use agat_sp_add_introns.pl from AGAT

• 0 views • link
written by Juke34 Juke34

Log in to answer this question.

More posts like this