HyPhy message board
http://www.hyphy.org/cgi-bin/hyphy_forums/YaBB.pl
Methodology Questions >> How to >> Dn/Ds ratio for each node in a tree?
http://www.hyphy.org/cgi-bin/hyphy_forums/YaBB.pl?num=1166122101

Message started by promateus on Dec 14th, 2006 at 10:48am

Title: Dn/Ds ratio for each node in a tree?
Post by promateus on Dec 14th, 2006 at 10:48am
Hi Sergei, I want to ask if it’s possible to estimate in Hyphy the Dn/Ds ratio for each of nodes in a given tree and if so how to do this?
Assuming that we have separate lentiviruses species (SIV and HIV from human and monkey) and I want to trace the change in the Dn/Ds ratio from the shared common ancestor to the tip of each group of species and subspecies, Is there is a way to do this?  

I hope all the best
promateus

Title: Re: Dn/Ds ratio for each node in a tree?
Post by artpoon on Dec 14th, 2006 at 11:06am
Dear Promateus,

This is a very simple thing to do in HyPhy!  We refer to this exercise as fitting a "local" evolutionary model.  Here's a walkthrough on how to go about obtaining estimates for each node (i.e. branch) in the tree.

1. Open your alignment in the HyPhy GUI by selecting the menu option File > Open > Open Data File..., and locating your file in the browser window that appears.

2. Select the entire alignment by selecting the menu option Edit > Select All.

3. Convert your selection into a 'data partition' by selecting Data > 'Selection -> Partition'

4. A new row will appear in your partition list below your alignment in the data viewer window.  Select the option "Codon" under the column Partition Type, specifying offset and genetic code in the window that appears.

5.  Under the column "Tree", select the tree topology either by choosing a file with a Newick string, loading one from memory, or estimating a tree de novo.

6. Under the column "Model", select a model of codon evolution, e.g. MG94xHKY85_3x4.

7. Under the column "Parameters", the option should read "Local".

8. Generate a likelihood function from these settings using the menu option Likelihood > Build Function.

9. Optimize the LF by selecting the menu option Likelihood > Optimize.

10. After HyPhy has obtained maximum likelihood estimates, a window should appear displaying the parameters of the likelihood function.  Each node has a unique identifier that is either the sequence name or "Node #".  Because you've fit a local model, each node is associated with a dN and dS parameter estimate.

HyPhy message board » Powered by YaBB 2.5.2!
YaBB Forum Software © 2000-2024. All Rights Reserved.