Hi Everyone,

I'm trying to find orthologs and lineage specific paralogs between two species. I tried using the ensembl homology pipeline but both my species are not on the database. Therefore, I tried to write my own similar pipeline. So far I've accomplished the following:

1. Blast all for every gene in both genomes
2. Filtering of blast results based of evalue and alignment length
3. Single Linkage clustering with MCL to form gene families
4. For each gene family, I did a protein alignment with PRANK and built a tree with Treebest, which also takes in the species trees and tries to build a gene tree accordingly.

My question deals with parsing these gene trees. I want to use these gene trees to find paralogs and orthologs between my two species but I'm not sure how to parse all of the topologies of these trees and how to determine paralogous or orthologous relationships.

Are there any programs that can take in gene trees and output a list of paralogs and orthologs?

Thanks.

LC
The <a href="http://etetoolkit.org">ETE toolkit</a> is indeed capable of doing that. You can use a species overlap algorithm to detect duplication and speciation events, or reconcile your gene tree with the expected species trees (<a href="http://www.plosone.org/article/info%3Adoi%2F10.1371%2Fjournal.pone.0004357">this paper</a> includes a comparison of both methods).

Briefly, you will need to load your gene tree as a PhyloTree object, and then call any of the tree.get_descendant_evol_events or tree.reconcile methods. Both methods will process your tree, label the nodes as speciation or duplication and return a list of speciation and duplication events. Then you can visualize your tree with `tree.show()` to see the predictions or process the list of events. If you used `tree.reconcile`, a reconciled-tree PhyloTree object will also be returned, including the inferred lost branches (also visible con tree.show). There are some examples in the ETE tutorial showing how to get orthology/paralogy prediction based on gene trees: http://pythonhosted.org/ete2/tutorial/tutorial_phylogeny.html#detecting-evolutionary-events