Hi, I have the gtf file for mouse, and I wanted to create bed files for the 5' and 3' UTRs using this gtf file ( I know there is method using the UCSC Table website, but I want to generate using the GTF file itself).
As far as I know, information about the UTR regions is not explicitly stated in the gtf file. So, what would be a good strategy to extract these bed files from the gtf file?
Thanks!
•
written
by
c_u
2 answers
Another option.
python extract_transcript_regions.py -i genes.gtf -o genes --gtf
now convert this blockbed (bed12) to bed6
cat genes_3utr.bed | bed12ToBed6 -i stdin -n > genes_3utr_bed6.bed
cat genes_5utr.bed | bed12ToBed6 -i stdin -n > genes_5utr_bed6.bed
The 5'UTR comprises the beginning of exon1 up to the base right upstream of the start codon. In contrast, the 3' UTR is the base right downstream of the stop codon to the end of the last exon. One would therefore need the GTF and a fasta (or BSgenome in case you want to use R) of the genome to get the required information. There is a Git repository that apparently has a script to do exactly what you want , but I never tested it towards accuracy. Maybe give it a careful try.
In case you want to put code together yourself, I would do for the 5'UTR:
- extract the first exon of every gene
- get the nucleotide sequence of it
- get coordinates of the start codon ATG
- set start(exon1) to start(ATG)-1 as 5'UTR
For the 3' UTR the same, but search for the last exon and the stops (TGA, TAA, TAG) and define 3'UTR as end(stopCodon)+1 to end(lastExon). Be sure to be strand-specific.
Log in to answer this question.
-
Some questions about 3'UTR regions from rat 6.0 fasta and gtf files
written by diqixiaoyaoer •
-
Extract data from GTF file
written by bioinformaticsnewbie •
-
How to extract exon, intron and UTR regions from gene?
written by ganeshram •
-
Genes with promoter and enhancer regions as GTF
written by hkarakurt
-
get CDS cordinates for genes, excluding UTR
written by J.F.Jiang
-
BedTools errors on RefSeq genes but not UCSC genes
written by bradbenjamin526 •
-
Download genome annotations from UCSC's MySQL database
written by James Ashmore
-
Determining the 5'UTR of Saccharomyces cerevisiae?
written by dashdp123 •
-
5"UTR and 3"UTR of Nile Tilapia Genome (Oreochromis niloticus)
written by wanziyi89 •
-
How Can We Find The Info For 3'Utr And 5'Utr In Gencode Gtf File?
written by J.F.Jiang
GTFs downloaded from Gencode have UTR annotations. The designation of 3' and 5' can be derived from the UTR coordinates and their corresponding strand.