This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Determining the correct reading frames for transcripts for dn/ds calculation

Hi everyone,

I am trying to calculate genome-wide dn/ds ratio between two closely related rat species using PAML. I have genome co-ordinates for transcript orthologs for the two species. However, I am not able to figure out how to decide the reading frame for translating the sequences in to amino acid sequences. Currently, I am translating into reading frame 1 for all genes on forward strand and into reading frame 6 for all on reverse strand. I know this is gross approximation. Can anyone give me any ideas on how I could tackle this issue? Where can I find information on transcript specific reading frame? I am ignoring the problem of alternate splicing for now.

Without the information on reading frames, I am ending up with biased dn/ds ration. Please help me. Thank you.

dn/ds ka/ks orf comparative genomics translation
• 2,707 views • link
written by shreyasibiswas88 • shreyasibiswas88

You just want the sequences of all the transcripts first?

• 0 views • link
written by Joe Joe

I have the nucleotide sequences of all transcripts. I can't decide which reading frame to translate each one into. Should I select the longest ORF for each transcript?

• 0 views • link
written by shreyasibiswas88 • shreyasibiswas88

It is typical to take the longest ORF. But here is a relevant study.

• 0 views • link
written by microfuge microfuge

0 answers

No answers yet.

Log in to answer this question.

More posts like this