Find corresponding branch from gene tree in species tree
0
0
Entering edit mode
7 weeks ago

I have constructed a gene tree for a single gene family across various species. The internal nodes of this tree represent gene duplication events and also speciation events. The leaves are individual genes of the gene families.

Say if I also have the corresponding species tree. How can I find the branch in the species tree that corresponds to any given branch in the gene tree? For example, say there is a duplication node in the gene tree I am particularly interested in, and I would like to find out where (on which branch) in the species tree this duplication event took place. How can I do that? And do the gene tree and species tree need to be completely reconciled in order to do this?

phylogeny family-tree gene species phylogenetics • 348 views
ADD COMMENT
0
Entering edit mode

check Notung: link

ADD REPLY
0
Entering edit mode

Following your suggestion I've been playing around with Notung and seems like its got a lot of great features. I ran:

java -jar notung -g GENE_TREE.nhx -s SPECIES_TREE.nh --reconcile

And got the file GENE_TREE.nhx.reconciled. I also tried the --savepng option to have a look at the result. My original GENE_TREE.nhx has the duplication of interest labelled as [&&NHX:D=Y:retro=Y], where retro=Y is used to identify this node of interest. However I don't find retro=Y in the output file GENE_TREE.nhx.reconciled, nor have I seen a way to annotate the corresponding branch in the png. Any Idea how I can do this?

ADD REPLY
0
Entering edit mode

I use the Notung's graphical user interface to display the duplication events and use a different software to highlight the nodes in the tree in which that event occurs.

ADD REPLY
0
Entering edit mode

This is for a pipeline where I will process hundreds of trees, so I need a programmatic solution :)

ADD REPLY
0
Entering edit mode

Well, good luck then ;)

ADD REPLY

Login before adding your answer.

Traffic: 974 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6