This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Getting the overlap between two GTF files

Hello,

I have two GTF files which contain the information of transcripts, and I want to get the overlap of transcripts between the two GTF file. Can anyone give me some advice?

Thanks!

rna-seq gtf
• 3,748 views • link
updated by Ram Ram • written by feather-W • feather-W

Ok, thanks for all of your help very much! I will try it.

• 0 views • link
written by feather-W • feather-W

4 answers

Maybe somebody knows something I don't, but I feel like bedtools should be able to handle your gtf files directly.

If it's not done already, you can filter for transcripts only of each file, then use bedtools:

awk '$3 == "transcript"' file1.gtf > file1.txOnly.gtf
awk '$3 == "transcript"' file2.gtf > file2.txOnly.gtf

bedtools intersect -u -a file1.txOnly.gtf -b file2.txOnly.gtf > file1_tx_overlapping_file2_tx.gtf

depending on your exact goals for what overlaps you want reported you can switch the -a and -b files or adjust options.

• 0 views • link
written by rfran010 • rfran010 
bedtools intersect \
    -a  <(awk '/^[^#]/ {printf("%s\t%d\t%s\t%s\n",$1,int($4)-1,$5,$0);}' file1.gtf  | sort -t $'\t' -k1,1 -k2,2n ) \
    -b  <(awk '/^[^#]/ {printf("%s\t%d\t%s\t%s\n",$1,int($4)-1,$5,$0);}' file2.gtf  | sort -t $'\t' -k1,1 -k2,2n )
• 0 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

Something in AGAT should work: https://agat.readthedocs.io/en/latest/?badge=latest

• 0 views • link
written by GenoMax GenoMax

The script agat_sp_compare_two_annotations.pl or agat_sp_sensitivity_specificity.pl
But if you are interested in the ranges then bedtools bedops or awk will be your friend

• 0 views • link
written by Juke34 Juke34

It would be helpful if recommendations for toolkits included actual solutions.

• 0 views • link
written by Alex Reynolds Alex Reynolds

Since I had not recommended a specific script in AGAT I had posted this as a comment so OP can check the documentation on their own. If Juke34 (author of AGAT) wants to post a stand-alone answer we can delete my original comment.

• 0 views • link
written by GenoMax GenoMax

Using bedops --intersect and gtf2bed will get their common genomic space:

bedops --intersect <(gtf2bed < transcripts.gtf) <(gtf2bed < transcripts.gtf) > answer.bed

If you want to know what transcripts overlap other transcripts, specifically, you could use bedmap --echo --echo-map:

bedmap --echo --echo-map <(gtf2bed < transcripts.gtf) > answer.bed

More information at: https://bedops.readthedocs.io/

• 0 views • link
written by Alex Reynolds Alex Reynolds

Log in to answer this question.

More posts like this