Welcome, Guest. Please Login
YaBB - Yet another Bulletin Board
 
  HomeHelpSearchLogin  
 
Relative rates with mid-point rooting? (Read 2767 times)
lichgen
YaBB Newbies
*
Offline


Curious HyPhy user

Posts: 2
Relative rates with mid-point rooting?
Feb 14th, 2012 at 11:39am
 
Hi,
HyPhy is a great tool for quickly calculating (pair-wise) relative rates in protein alignments (Standard Analysis -> PairwiseRelativeRate.bf). This method however requires to define an outgroup. Is there any way to calculate these rates without outgroups, but by using mid-point rooting instead?
Thanks
Back to top
 
 
IP Logged
 
Sergei
YaBB Administrator
*****
Offline


Datamonkeys are forever...

Posts: 1658
UCSD
Gender: male
Re: Relative rates with mid-point rooting?
Reply #1 - Feb 14th, 2012 at 11:51am
 
Hi lichgen,

I am not sure you can easily extend the test to use a mid-point root. Interpretation issues aside (under time reversible models, mid-point rooting is effectively uninformative), the standard relative rate test requires three sequences, one of which is designated as the outgroup. How does one define this sequence based on mid-point rooting?

It is possible to inquire whether path-lengths in the tree from the mid-point to any two taxa are of different lengths; this test is not the same as the relative rate test.

What is the biological question you are trying to answer?

Sergei
Back to top
 

Associate Professor
Division of Infectious Diseases
Division of Biomedical Informatics
School of Medicine
University of California San Diego
WWW WWW  
IP Logged
 
lichgen
YaBB Newbies
*
Offline


Curious HyPhy user

Posts: 2
Re: Relative rates with mid-point rooting?
Reply #2 - Feb 14th, 2012 at 2:16pm
 
Hi Sergei,
I am trying to identify phylogenetically old genes in a symbiotic fungal species, by looking at reduced evolutionary rates, compared to homologues in a defined group of related fungal taxa. One of these reference taxa is more distant and used as outgroup. I considered a gene of my fungal species with a reduced rate, if it constantly showed the lowest rate in all triplet comparisons within the HyPhy rate test (I processed several thousand alignments, with exactly one sequence per species).  That worked well and let me identify a small number of genes.

The biological idea behind this approach is that symbiosis was acquired early in evolution, and some (not all) key genes haven't much evolved since then, in order to maintain the symbiotic lifestyle.

In a second approach I am now working with MCL clustered gene families, instead of well conserved homologue groups. These gene families are naturally more diverse, with usually more than one, non monophyletic outgroup sequence, and can therefore not be rooted. My idea was now to identify different path lengths (as you described it above) for each two species in an unrooted tree instead.

Since I have several thousand alignments it would be very helpful if a test existed, which could perform this mid-point path length calculation, without the need of building a tree first.

Olaf 
Back to top
 
 
IP Logged