Latest
Open
RNA-Seq
ChIP-Seq
SNP
Assembly
Tutorials
Tools
Jobs
Forum
Planet
All »
Home
Log In
Question: How To Calculate Substitution Rate Of Noncoding Rnas
1
gravatar for lazyplant
7.6 years ago by
lazyplant • 20
lazyplant • 20 wrote:

hi, everyone I am working on calculate the substitution rate of noncoding rnas now. if i have a multiple sequence alignment from three species, the topology of tree is ((1,2),3), how can i use baseml and the REV substitution model to get my desire result? i have tried the example 'brown.nuc' and 'brown.trees' in paml4.7 package, but i cannot make out which number is subtitution rate of that sequences. i have read a lot of papers, but the formula is too difficult to understand, can you give me suggestion?any reply will be appreciated!
• 3.5k views
ADD COMMENT • link •
modified 17 months ago by Biostar ♦♦ 20 • written 7.6 years ago by lazyplant • 20
1
gravatar for Charles Warden
7.0 years ago by
Charles Warden ♦ 8.0k
Duarte, CA
Charles Warden ♦ 8.0k wrote:

So, it's been a while since I did this type of analysis, but there are two levels of modeling:

1) Modeling substitution rates between two samples (to try and take the probability of multiple substitutions at the same site into consideration). This wikipedia page talks about the models for doing this:

http://en.wikipedia.org/wiki/Models_of_DNA_evolution

2) Using an outgroup (sample 3 in your case) to try and estimate the ancestral state (or, perhaps more precisely, the rate of divergence from an ancestral state). In this particular paper (looking at substitutions at intergenic and synonymous sites), this was done with a relative rate test. See links below for paper and description of strategy:

http://www.plosgenetics.org/article/info%3Adoi%2F10.1371%2Fjournal.pgen.0020163

http://en.wikipedia.org/wiki/Relative_rate_test

From a practical standpoint, I don't remember what specific program outputs look like. However, I know MEGA is typically one of the most popular packages for evolutionary genomics analysis, and I think it has fairly comprehensive documentation:

http://www.megasoftware.net/

If you can get a hold of this textbook, it might also be helpful to you:

http://www.amazon.com/Fundamentals-Molecular-Evolution-Dan-Graur/dp/0878932666
ADD COMMENT • link modified 7.0 years ago • written 7.0 years ago by Charles Warden ♦ 8.0k
Please log in to add an answer.

Similar posts • Search »

Content
Help
Access
Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 2254 users visited in the last hour
_