- THIS ARTICLE
- Full Text
- Full Text (PDF)
-
All Versions of this Article:
genetics.103.025692v1
170/1/419 most recent - Alert me when this article is cited
- Alert me if a correction is posted
- SERVICES
- Similar articles in this journal
- Similar articles in PubMed
- Alert me to new issues of the journal
- Download to citation manager
- Reprints & Permissions
- CITING ARTICLES
- Citing Articles via HighWire
- Citing Articles via Google Scholar
- GOOGLE SCHOLAR
- Articles by Suchard, M. A.
- Search for Related Content
- PUBMED
- PubMed Citation
- Articles by Suchard, M. A.
Originally published as Genetics Published Articles Ahead of Print on March 21, 2005.
Genetics, Vol. 170, 419-431, May 2005, Copyright © 2005
doi:10.1534/genetics.103.025692
Stochastic Models for Horizontal Gene Transfer
Taking a Random Walk Through Tree Space
Marc A. Suchard1
Department of Biomathematics, David Geffen School of Medicine, University of California, Los Angeles, California 90095-1766
1 Address for correspondence: Department of Biomathematics, David Geffen School of Medicine, UCLA, 650 Charles Young Dr., Box 951766, Los Angeles, CA 90095-1766.
E-mail: msuchard{at}ucla.edu
Horizontal gene transfer (HGT) plays a critical role in evolution across all domains of life with important biological and medical implications. I propose a simple class of stochastic models to examine HGT using multiple orthologous gene alignments. The models function in a hierarchical phylogenetic framework. The top level of the hierarchy is based on a random walk process in "tree space" that allows for the development of a joint probabilistic distribution over multiple gene trees and an unknown, but estimable species tree. I consider two general forms of random walks. The first form is derived from the subtree prune and regraft (SPR) operator that mirrors the observed effects that HGT has on inferred trees. The second form is based on walks over complete graphs and offers numerically tractable solutions for an increasing number of taxa. The bottom level of the hierarchy utilizes standard phylogenetic models to reconstruct gene trees given multiple gene alignments conditional on the random walk process. I develop a well-mixing Markov chain Monte Carlo algorithm to fit the models in a Bayesian framework. I demonstrate the flexibility of these stochastic models to test competing ideas about HGT by examining the complexity hypothesis. Using 144 orthologous gene alignments from six prokaryotes previously collected and analyzed, Bayesian model selection finds support for (1) the SPR model over the alternative form, (2) the 16S rRNA reconstruction as the most likely species tree, and (3) increased HGT of operational genes compared to informational genes.
This article has been cited by other articles:
![]() |
S. Linz, A. Radtke, and A. von Haeseler A Likelihood Framework to Measure Horizontal Gene Transfer Mol. Biol. Evol., June 1, 2007; 24(6): 1312 - 1319. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. Ane, B. Larget, D. A. Baum, S. D. Smith, and A. Rokas Bayesian Estimation of Concordance among Gene Trees Mol. Biol. Evol., February 1, 2007; 24(2): 412 - 426. [Abstract] [Full Text] [PDF] |
||||
