Welcome, Guest. Please Login
YaBB - Yet another Bulletin Board
 
  HomeHelpSearchLogin  
 
Comparisons among trees from GABranch (Read 2371 times)
Daryn
YaBB Newbies
*
Offline


Feed your monkey!

Posts: 4
Comparisons among trees from GABranch
May 19th, 2010 at 2:54pm
 
Hi

This may be a silly question, but I’ve been using the GABranch method to look at variation in dN/dS among lineages for different protein domains of a gene. I’d now like to compare the trees to each other to see if there are significant differences among domains. Is there an easy way to do this through Datamonkey or do I need to switch to the HyPhy interface?

Thanks for the help!
Daryn
Back to top
 
 
IP Logged
 
wayne
Global Moderator
*****
Offline


I love YaBB 1G - SP1!

Posts: 57
San Diego, CA
Gender: male
Re: Comparisons among trees from GABranch
Reply #1 - May 19th, 2010 at 3:53pm
 
Hi Daryn, if by trees you mean topologies then the KHTest.bf in ./TemplateBatchFiles in the HyPhy distribution will do this. If by trees you mean the distribution of selection classes among lineages then the answer is not trivial. There is no statistically robust way to do this, but some hackish measures can be attempted. For instance you could do a branch by branch comparison of estimated dNdS values for a full protein alignment versus that of each domain. Clearly, the single result output from GABranch is not sufficient for this and we would need to process the credible set of GABranch models. If you're keen to pursue the hackish approach then please send your job upload id's so i can access the full output files. cheers ./w
Back to top
 

Assistant Project Scientist&&Antiviral Research Center&&Department of Pathology&&University of California, San Diego
WWW WWW  
IP Logged
 
Daryn
YaBB Newbies
*
Offline


Feed your monkey!

Posts: 4
Re: Comparisons among trees from GABranch
Reply #2 - May 20th, 2010 at 12:15pm
 
Sorry, I meant the distribution of selected classes among domains. We’re basically trying to test a long standing hypothesis that there is variation in the strength of purifying selection among domains at this particular protein. I’m definitely up for a hackish approach – thanks! The upload ID for the entire gene is 769413173912666.1 and the upload IDs for the domains are 476578907412493.1; 930161873912609.1; and 99377543463329.1.

In the mean time, I started using HyPhy to compare likelihoods where each domain evolves independently to likelihoods where I constrain the nonsynonymous rate in one domain to that in another domain (doing all pairwise comparisons for both nonsynonymous and synonymous rates). Is that similar to the hackish approach you were thinking of? I know the dN and dS estimates aren’t exactly the same as those provided by GABranch, but this method was all I could think of Smiley

Thanks for the help – I really appreciate it!
Daryn
Back to top
 
 
IP Logged