This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Constructing the ancestral state of amino acid

Hello Everyone,

I have the results from codeml analysis where I have identified few sites that are under positive selection in my target species. I am now looking for the change in amino acid in the nearest species and the ancestral amino acid. I have around 20 such genes where I have identified positive selection in orthologs of these genes. I can easily identify the amino acid in nearest species from the MSA of these genes. I am wondering if it is fair enough to identify the ancestral amino acid same way- looking at MSA for amino acid in the distant species retained in the orthology cluster. or there is any other way to construct the ancestral state of these sites?

I hope I am clear in asking the question.

Thanks

paml codeml

2 answers

FastML http://fastml.tau.ac.il/ | http://nar.oxfordjournals.org/content/40/W1/W580.long is a web server for ancestral sequence reconstruction, using the sequence MSA plus phylogenetic tree as input. I don't know if that answers your question, but I have seen papers that use both PAML and FastML and compare the results.

You could try from the DNA sequence using the Ortheus tool. (http://www.ncbi.nlm.nih.gov/pmc/articles/PMC2577868/). It allows to reconstruct the ancestral state of DNA sequences aligned.

Log in to answer this question.