Thanks Eric, I did look at this tool https://omictools.com/polya-svm-tool but just wasn't sure if this was the way to go. Will give that a try. Thanks again !
Finding polyA site from Ensembl GFF file
Hi, I'm working on a project with pig genome. I'm looking for polyA sites information and downloaded this GFF file ftp://ftp.ncbi.nlm.nih.gov/genomes/all/GCF/000/003/025/GCF_000003025.6_Sscrofa11.1/GCF_000003025.6_Sscrofa11.1_rna.gbff.gz from Ensembl. There are 'three_prime_UTR' regions in this GFF file (field 3). Is there an existing method/tool which can be used directly to get the polyA site information from these 3' UTRs?
I did go through this post How to find polyA sites from gtf/gff?
Any help is much appreciated. Thanks in advance.
1 answer
I'm not aware of any existing tool to find polyA signals directly from GTF. Typically, you need to convert your region coordinates into sequences first. Converting coordinates annotated as 3UTR in your GTF to sequences is a good start.
People have implemented various machine learning algorithms, from simple motif finder to slightly more complicated variations using SVM or HMM. You might find this page useful: https://omictools.com/polyadenylation-prediction-category
You may also find some knockdown RNAseq experiments (CFIM25, etc) targeting at identifying alternative polyadenylation (APA) useful. There are quite a few out there.
Hope this info helps.
Log in to answer this question.
More posts like this
-
Average length of mRNA on NCBI annotation seems too high
written by yh1126 •
-
Some questions about 3'UTR regions from rat 6.0 fasta and gtf files
written by diqixiaoyaoer •
-
NCBI annotated sequence use as input of Roary
written by Manoj •
-
biomartr::getGenome() issue
written by qstefano •
-
Where to find a valid GTF file for stringtie asssembly?
written by lakhujanivijay
-
How to choose correct VCF for Human reference genome?
written by marongiu.luigi
-
what is the corresponding proteome or transcriptome (cDNA fasta in Ensembl) file in NCBI?
written by SMILE
-
How to extract all promoter regions in multi-fasta format from genome using GFF?
written by rimgubaev
-
Receiving low query covers for blast searches on genome files, but not transcript files
written by Jacob •
-
human genome toplevel files
written by biolab
What exactly do you want to do with the information? Is the first base after the 3'UTR not a "good-enough" proxy?
Hi, I want to capture all polyA sites (alternative polyA sites as well) for all transcripts in pig genome. But it doesn't look like the Ensembl GFF file has multiple 'three_prime_UTR' records for a given transcript. I thought it would have them.