Question: Tree input for PAML ancestral state reconstruction
gravatar for memory_donk
5.2 years ago by
memory_donk310 wrote:

Hi Biostars,

I have a set of about 6,000 groups of orthologous proteins. Each orthologous group has a representative protein from anywhere between 7 and 16 species (1 per species). I'm trying to do ancestral state reconstructions for each of these groups based on their well-established phylogenies, and then to compare each amino acid position in a species of interest to the reconstructed sequences at key ancestral nodes. 

I've been trying to use PAML (codeml) for this purpose, following the method used in this blog post

My problem is that when using clock = 0 (no molecular clock), PAML requires an unrooted tree. According to PAML's manual " a rooted tree has a bifurcation at the root, while an unrooted tree has a trifurcation or multifurcation at the root."

This is problem is that the mammalian tree that I have (composed of eutherians and marsupials) is a well-established bifurcating tree. How can I do an ancestral state reconstruction with PAML which requires a multifurcating tree when the true phylogeny is bifurcating? Adding an outgroup like platypus or chicken would still be a rooted tree since they're both outgroups and a polytomy of chicken/platypus, eutherians and marsupials would be false. I'm sure I've just deeply misunderstood something along the way. Any help would be really greatly appreciated! 

software error gene • 3.1k views
ADD COMMENTlink modified 5.2 years ago by Brice Sarver3.6k • written 5.2 years ago by memory_donk310
gravatar for Brice Sarver
5.2 years ago by
Brice Sarver3.6k
United States
Brice Sarver3.6k wrote:

The short answer is to unroot the tree. This can be accomplished easily in R using unroot() in ape.

It is possible for a tree that is normally bifurcating when rooted to have a node that has an order greater than two when unrooted. Consider a three-taxa case. There is one possible unrooted tree (shaped like a 'Y') and three possible rooted trees (whether rooted using the first, second, or third taxon) that are bifurcating. Remember that the placement of the root is a hypothesis.

ADD COMMENTlink modified 14 months ago by Ram32k • written 5.2 years ago by Brice Sarver3.6k

I can confirm the above. 'unroot' will work well for this.

On Brice's suggestion I have been using 'ape' for a while know. It's easily the most useful phylogenetics packages for R out there.

ADD REPLYlink modified 14 months ago by Ram32k • written 5.2 years ago by confusedious420

Thanks Brice! Codeml is happy now.

ADD REPLYlink modified 14 months ago by Ram32k • written 5.2 years ago by memory_donk310
Please log in to add an answer.


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