1 / 22

Rooting Phylogenetic Trees with Non-reversible Substitution Models

Rooting Phylogenetic Trees with Non-reversible Substitution Models. Von Bing Yap* and Terry Speed § *Statistics and Applied Probability, National University of Singapore § Statistics, University of California Reference: BMC Evolutionary Biology 5:1 (2005). Molecular Phylogenetics.

leigh
Download Presentation

Rooting Phylogenetic Trees with Non-reversible Substitution Models

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. Rooting Phylogenetic Trees with Non-reversible Substitution Models Von Bing Yap* and Terry Speed§ *Statistics and Applied Probability, National University of Singapore §Statistics, University of California Reference: BMC Evolutionary Biology 5:1 (2005)

  2. Molecular Phylogenetics • From alignments to trees. • Many methods: parsimony, distance, stochastic models.

  3. Reversible Models • Almost all substitution models are reversible: for example, Pr(anc=A, des=C) = Pr(anc=C, des=A). • Rooted trees that give the same unrooted tree are indistinguishable.

  4. Stationary Models • Character states have the same frequencies everywhere on the tree. • Root can be identified (Yang 1994, Huelsenbeck et al. 2001).

  5. Nonstationary Models • Yang and Roberts (1985) • Galtier and Gouy (1998)

  6. SUBSTITUTION MODELS NON-STATIONARY STATIONARY REVERSIBLE

  7. The Simplest NSTA Model • Parameters: rooted tree topology θ: root base frequency Q: rate matrix (calibrated) branch lengths No relationship between θ and Q.

  8. Specialisations • If θ is the equilibrium distribution of Q, get STA. • If in addition, Q satisfies the detailed balance conditions, get REV.

  9. Probability of alignment • Felsenstein’s algorithm can be used to compute the probability of one site. • Multiplying across sites gives probability of alignment.

  10. Tree Inference • Fix a rooted tree. • Find the most likely parameter values. • The maximum likelihood is the support of the tree. • Choose tree with highest support.

  11. Site Heterogeneity • Codon positions, secondary structure. • Deterministic or random relative rates can be accommodated in the model. • Two deterministic models: codon position, and codon position + fast/slow.

  12. Two deterministic models • codon: 3 fixed unknown rates, corresponding to codon positions, with weighted average 1. • codonsite: get two classes of amino acids (fast/slow) from CLUSTAL alignment output. Coupled with codon positions, get 6 unknown rates with weighted average 1.

  13. Test Data Sets • A: human, chimp, gorilla • B: human, mouse, rat • C: human, chimp, gorilla, orangutan • D: human, chimp, mouse, rat • E: human, mouse, chicken, frog • 13 mitochondrial protein-coding genes

  14. Method • Unrooted tree is assumed known. • For each rooted tree consistent with the unrooted tree, its support is the maximum loglikelihood upon finding the MLE of the process parameters and branch lengths.

  15. Method (continued) • Three processes: REV, STA, NSTA • Three site models: novar (no variations), codon (3 classes), codonsite (6 classes).

  16. Method (continued) • Two outcomes (a) number of genes for which the correct rooted tree is the most likely (b) does the model get the right rooted tree when the loglikelihoods are summed over genes?

  17. Number of successes

  18. Combined genes: Does it get the right tree?

  19. Discussion (1) • In general, NSTA fits much better than STA, which fits much better than REV, by the likelihood ratio test criterion. • Not only does NSTA get the right tree more often than STA, it is also more discriminative: the best tree has much larger support compared to the other trees.

  20. Discussion (2) • The codon+site model of site variation is very crude, and this may explain why the performance is worse than codon model. • Need to use better methods. Also need to compare with random model, like discrete gamma.

  21. Discussion (3) • The NSTA only has 3 more parameters than STA, and 6 more than REV, so the extra computation is not heavy. • Also, since it is possible to identify the root, perhaps NSTA should be used routinely.

  22. Discussion (4) • Constraint on NSTA: base compositions of sequences that are equally distant from the root are the same. This may not hold. • Software freely available upon request. Email stayapvb@nus.edu.sg

More Related