Inference of Gain and Loss Events from Phyletic Patterns Using Stochastic Mapping and Maximum Parsimony—A Simulation Study
Open Access
- 1 January 2011
- journal article
- research article
- Published by Oxford University Press (OUP) in Genome Biology and Evolution
- Vol. 3, 1265-1275
- https://doi.org/10.1093/gbe/evr101
Abstract
Bacterial evolution is characterized by frequent gain and loss events of gene families. These events can be inferred from phyletic pattern data—a compact representation of gene family repertoire across multiple genomes. The maximum parsimony paradigm is a classical and prevalent approach for the detection of gene family gains and losses mapped on specific branches. We and others have previously developed probabilistic models that aim to account for the gain and loss stochastic dynamics. These models are a critical component of a methodology termed stochastic mapping, in which probabilities and expectations of gain and loss events are estimated for each branch of an underlying phylogenetic tree. In this work, we present a phyletic pattern simulator in which the gain and loss dynamics are assumed to follow a continuous-time Markov chain along the tree. Various models and options are implemented to make the simulation software useful for a large number of studies in which binary (presence/absence) data are analyzed. Using this simulation software, we compared the ability of the maximum parsimony and the stochastic mapping approaches to accurately detect gain and loss events along the tree. Our simulations cover a large array of evolutionary scenarios in terms of the propensities for gene family gains and losses and the variability of these propensities among gene families. Although in all simulation schemes, both methods obtain relatively low levels of false positive rates, stochastic mapping outperforms maximum parsimony in terms of true positive rates. We further studied the factors that influence the performance of both methods. We find, for example, that the accuracy of maximum parsimony inference is substantially reduced when the goal is to map gain and loss events along internal branches of the phylogenetic tree. Furthermore, the accuracy of stochastic mapping is reduced with smaller data sets (limited number of gene families) due to unreliable estimation of branch lengths. Our simulator and simulation results are additionally relevant for the analysis of other types of binary-coded data, such as the existence of homologues restriction sites, gaps, and introns, to name a few. Both the simulation software and the inference methodology are freely available at a user-friendly server: http://gloome.tau.ac.il/.Keywords
This publication has 60 references indexed in Scilit:
- Networks of Gene Sharing among 329 Proteobacterial Genomes Reveal Differences in Lateral Gene Transfer Frequency at Different Phylogenetic DepthsMolecular Biology and Evolution, 2010
- Inference and Characterization of Horizontally Transferred Gene Families Using Stochastic MappingMolecular Biology and Evolution, 2009
- Genomics of bacteria and archaea: the emerging dynamic view of the prokaryotic worldNucleic Acids Research, 2008
- A likelihood framework to analyse phyletic patternsPhilosophical Transactions Of The Royal Society B-Biological Sciences, 2008
- Taxon Sampling Affects Inferences of Macroevolutionary Processes from Phylogenetic TreesSystematic Biology, 2008
- Coevolution of gene families in prokaryotesGenome Research, 2008
- Counting labeled transitions in continuous-time Markov models of evolutionJournal of Mathematical Biology, 2007
- Three distinct modes of intron dynamics in the evolution of eukaryotesGenome Research, 2007
- Horizontal gene transfer, genome innovation and evolutionNature Reviews Microbiology, 2005
- Language-tree divergence times support the Anatolian theory of Indo-European originNature, 2003