This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Can you suggest tools or script to get canonical sequences from a cds fasta file?

I have a CDS fasta file from NCBI, and there are around 11k duplicated transcripts. I would like to only get the canonical sequences (would be the longest transcript, I suppose?), so there won't be duplicated coding sequences in my downstream analyses. Can you please suggest tools I can easily use?

Thanks!

genome cds fasta
• 1,276 views • link
updated by Ram Ram • written by Yuto • Yuto

The longest transcript is not necessary canonical, as: encoding the most abundant protein/most conserved among species etc. The 11k seems to be to small for human/mammalian transcripts. Anyway, ENSEMBL GTF files should have canonical tag (no idea if this is valid for annotations of more exotic species) so you can get the canonical selecting these and extracting sequences from fasta using i.e. bedtools.

• 0 views • link
written by Darked89 Darked89

1 answer

Which organism is it and for what purpose do you need the data? For human and mouse you can use the CCDS database. The meaning of "canonical" may also vary by organism. In the end it may be rather a subjective choice (for human it is "what's annotated as canonical in the database") and depends on your purpose. I would prefer experimentally validated or manually annotated over the longest sequence.

• 0 views • link
written by Michael Dondrup Michael Dondrup

Log in to answer this question.

More posts like this