1 / 30

GRAPPA: Large-scale whole genome phylogenies based upon gene order evolution

GRAPPA: Large-scale whole genome phylogenies based upon gene order evolution. Tandy Warnow, UT-Austin Department of Computer Sciences Institute for Cellular and Molecular Biology Program in Evolution, Ecology, and Behavior Center for Computational Biology and Bioinformatics. A. C. A. D.

danica
Télécharger la présentation

GRAPPA: Large-scale whole genome phylogenies based upon gene order evolution

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. GRAPPA: Large-scale whole genome phylogenies based upon gene order evolution Tandy Warnow, UT-Austin Department of Computer Sciences Institute for Cellular and Molecular Biology Program in Evolution, Ecology, and Behavior Center for Computational Biology and Bioinformatics

  2. A C A D X E Y B E Z W C F B D F Whole-Genome Phylogenetics

  3. Genomes As Signed Permutations 1 –5 3 4 -2 -6or6 2 -4 –3 5 –1 etc.

  4. 1 2 3 –8 –7 –6 –5 -4 9 10 1 2 3 9 -8 –7 –6 –5 –4 10 1 2 3 9 4 5 6 7 8 10 Genomes Evolve by Rearrangements 1 2 3 4 5 6 7 8 9 10 • Inversion (Reversal) • Transposition • Inverted Transposition

  5. Genome Rearrangement Has A Huge State Space • DNA sequences : 4 states per site • Signed circular genomes with n genes: states, 1 site • Circular genomes (1 site) • with 37 genes: states • with 120 genes: states

  6. Why use gene orders? • “Rare genomic changes”: huge state space and relative infrequency of events (compared to site substitutions) could make the inference of deep evolution easier, or more accurate. • Our research shows this is true, but accurate analysis of gene order data is computationally very intensive!

  7. Phylogeny reconstruction from gene orders • Distance-based reconstruction: estimate pairwise distances, and apply methods like Neighbor-Joining or Weighbor • “Maximum Parsimony”: find tree with the minimum length (inversions, transpositions, or other edit distances) • Maximum Likelihood: find tree and parameters of evolution most likely to generate the observed data

  8. This talk • Distance-based methods: we show how to estimate genomic distances appropriately for the Generalized Nadeau-Taylor model • Parsimony-style methods: we can find very good solutions to NP-hard problems (inversion and breakpoint phylogeny) quite quickly • Validation of these approaches on real and simulated data

  9. Distance-based methods

  10. Genomic Distance Estimators • Standard: • Breakpoint distance • (Minimum) Inversion distance • Our estimators: We attempt to estimate the actual number of events (the “true evolutionary distance”): • EDE [Moret et al, ISMB’01] • Approx-IEBP [Wang and Warnow, STOC’01] • Exact-IEBP [Wang, WABI’01]

  11. Breakpoint Distance • Breakpoint distance=5 1 2 3 4 5 6 7 8 9 10 1 –3 –2 4 5 9 6 7 8 10

  12. Minimum Inversion Distance • Inversion distance=3 1 2 3 4 5 6 7 8 9 10 1 2 3 –8 –7 –6 –5 –4 9 10 1 8 –3 –2 –7 –6 –5 –4 9 10 1 8 –3 7 2 –6 –5 –4 9 10

  13. Measured Distance vs. Actual Number of Events Breakpoint Distance Inversion Distance 120 genes, inversion-only evolution

  14. Generalized Nadeau-Taylor Model • Three types of events: • Inversions • Transpositions • Inverted Transpositions • Events of the same type are equiprobable • Probability of the three types have fixed ratio: Inv : Trp : Inv.Trp = (1-a-b):a:b

  15. Estimating True Evolutionary Distances for Genomes Given fixed probabilities for each type of event, we estimate the expected breakpoint distance after k random events, or the expected inversion distance after k random events. Inverting these functions gives us a better estimate of true evolutionary distances. • Approx-IEBP [Wang and Warnow 2001] • Exact-IEBP [Wang 2001] • EDE [Moret et al, ISMB 2001] (inversion based)

  16. Estimating True Evolutionary Distances for Genomes (cont.) Estimating the expected Inversion distance: EDE [Moret, Wang, Warnow, Wyman 2001] • Closed-form formula based upon an empirical estimation of the expected inversion distance after k random events (based upon 120 genes and inversion only, but robust to errors in the model) . • Polynomial time, fastest of the three.

  17. Absolute Difference • 120 genes • Inversion only evolution • (Similar relative • performance under • other models)

  18. Accuracy of Neighbor Joining Using Distance Estimators • 120 genes • All three event types equiprobable • 10, 20, 40, 80, and 160 genomes • Similar relative • performance under • other models

  19. Summary of Distance-based Reconstruction Methods • Statistically-based estimation of genomic distances improves NJ analyses. • Our IEBP estimators assume knowledge of the probabilities of each type of event, but are robust to model violations. • EDE is based upon an inversion-only evolutionary model, but is robust. • Best performing method: Weighbor(EDE); second best is NJ(EDE); both are robust to model violations. • Worst performing is NJ(BP). • Accuracy is very good, except when very close to saturation.

  20. A A D D B B 3 3 Total length = 18 6 C C E F 4 2 Maximum Parsimony on Rearranged Genomes (MPRG) • The leaves are rearranged genomes. • Find the tree that minimizes the total number of rearrangement events

  21. Optimization problems for gene order phylogeny • Breakpoint phylogeny: find the phylogeny which minimizes the total number of breakpoints (NP-hard, even to find the median of three genomes) • Inversion phylogeny: find the phylogeny which minimizes the sum of inversion distances on the edges (NP-hard, even to find the median of three genomes)

  22. Inversion and Breakpoint phylogenies • When the data are close to saturated, even Weighbor(EDE) analyses are insufficiently accurate. In these cases, our initial investigations suggest that the inversion and breakpoint phylogeny approaches may be superior. • Problem: finding the best trees is enormously hard, since even the “point estimation” problem is hard (worse than estimating branch lengths in ML). Local optimum MP score Global optimum Phylogenetic trees

  23. GRAPPA (Genome Rearrangement Analysis under Parsimony and other Phylogenetic Algorithms) http://www.cs.unm.edu/~moret/GRAPPA/ • Heuristics for NP-hard optimization problems • Fast polynomial time distance-based methods • Contributors: U. New Mexico,U. Texas at Austin, Universitá di Bologna, Italy • Freely available in source code at this site. • Project leader: Bernard Moret (UNM) (moret@cs.unm.edu)

  24. Benchmark gene order dataset: Campanulaceae • 12 genomes + 1 outgroup (Tobacco), 105 gene segments • NP-hard optimization problems: breakpoint and inversion phylogenies (techniques score every tree) 1997: BPAnalysis (Blanchette and Sankoff): 200 years (est.)

  25. Benchmark gene order dataset: Campanulaceae • 12 genomes + 1 outgroup (Tobacco), 105 gene segments • NP-hard optimization problems: breakpoint and inversion phylogenies (techniques score every tree) 1997: BPAnalysis (Blanchette and Sankoff): 200 years (est.) 2000: Using GRAPPA v1.1 on the 512-processor Los Lobos Supercluster machine: 2 minutes (200,000-fold speedup per processor)

  26. Benchmark gene order dataset: Campanulaceae • 12 genomes + 1 outgroup (Tobacco), 105 gene segments • NP-hard optimization problems: breakpoint and inversion phylogenies (techniques score every tree) 1997: BPAnalysis (Blanchette and Sankoff): 200 years (est.) 2000: Using GRAPPA v1.1 on the 512-processor Los Lobos Supercluster machine: 2 minutes (200,000-fold speedup per processor) 2003: Using latest version of GRAPPA: 2 minutes on a single processor (1-billion-fold speedup per processor)

  27. Summary • Weighbor(EDE) and NJ(EDE) are highly accurate polynomial time distance-based reconstructions, except when datasets are close to saturated • GRAPPA (inversion phylogeny or breakpoint phylogeny) produces highly accurate estimates of trees, and even of ancestral gene orders, acceptably fast

  28. Limitations and ongoing research • Current methods limited to single chromosomes with equal gene content (or very small amounts of deletions and duplications) -- we are working on developing reliable techniques for genomes with unequal gene content

  29. Acknowledgements • Funding: The David and Lucile Packard Foundation, and The National Science Foundation. • Collaborators: Bernard Moret (UNM), David Bader (UNM), Bob Jansen (UT), Linda Raubeson (CWU) • Students: Li-San Wang (now postdoc of Junhyong Kim at Penn), Jijun Tang (UNM), and others at UNM

  30. Phylolab, U. Texas Please visit us at http://www.cs.utexas.edu/users/phylo/

More Related