We have two separate GFF files. One is UTR GFF file, describing the 5' and 3' UTR regions. The other is CDS GFF describing coding sequences. Is there a tool to generate the intron GFF according to these two files? Or, how is the intron GFF file usually produced?
If GFF3, use GenomeTools.
Usage: gt gff3 [option ...] [GFF3_file ...]
Parse, possibly transform, and output GFF3 files.
-addintrons add intron features between existing exon features
default: no
Let me answer your two questions in reverse order.
Or, how is the intron GFF file usually produced?
In my experience, different feature types aren't typically stored in separate files--in other words, you don't have a CDS file, a UTR file, an intron file, etc, you simply have a single file with all the features in it. That doesn't mean your approach is incorrect, it just isn't typical and doesn't provide any immediate benefit (unless of course you are running scripts that have been built to expect it).
Is there a tool to generate the intron GFF according to these two files?
Perhaps, but it shouldn't be too difficult to do with minimal scripting experience. Once you have determined the exon coordinates using the CDS and UTR data, then simply create an intron feature to fill in the space between each adjacent pair of exon features.
Haibao mentioned the very useful GenomeTools utility, but to use that you would still first have to determine the exon coordinates and provide them as input. If you can calculate the exon coordinates from a set of CDS and UTR coordinates, then surely you can calculate intron coordinates from a set of exon coordinates.
You would think that something such basic tasks would be easy to do with a variety of existing tools. The problem is that people use the GFF3 (and GFF3-like) formats so differently as to make it quite difficult to handle every possible case. For example, some people include exons and CDSs but leave out UTRs, some people include CDSs and UTRs but leave out exons, etc. In each case, it is possible to infer the missing features, but it gets complicated when you have to account for every possible case.
I have a similar question. If I have a GFF file with CDS and UTR features how can I find out the exon start and end ?
Thanks! -Abhi
Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
version 2.3.6
Be careful with this approach though. If exon features are not explicitly defined, then gt will not create any intron features. I've made that mistake a few times.
Pretty encouraging. I will try it. Thanks!