Using the nucleotide substitution rate matrix to detect horizontal gene transfer.
BACKGROUND
Horizontal gene transfer (HGT) has allowed bacteria to evolve many new capabilities. Because transferred genes perform many medically important functions, such as conferring antibiotic resistance, improved detection of horizontally transferred genes from sequence data would be an important advance. Existing sequence-based methods for detecting HGT focus on changes in nucleotide composition or on differences between gene and genome phylogenies; these methods have high error rates.


RESULTS
First, we introduce a new class of methods for detecting HGT based on the changes in nucleotide substitution rates that occur when a gene is transferred to a new organism. Our new methods discriminate simulated HGT events with an error rate up to 10 times lower than does GC content. Use of models that are not time-reversible is crucial for detecting HGT. Second, we show that using combinations of multiple predictors of HGT offers substantial improvements over using any single predictor, yielding as much as a factor of 18 improvement in performance (a maximum reduction in error rate from 38% to about 3%). Multiple predictors were combined by using the random forests machine learning algorithm to identify optimal classifiers that separate HGT from non-HGT trees.


CONCLUSIONS
The new class of HGT-detection methods introduced here combines advantages of phylogenetic and compositional HGT-detection techniques. These new techniques offer order-of-magnitude improvements over compositional methods because they are better able to discriminate HGT from non-HGT trees under a wide range of simulated conditions. We also found that combining multiple measures of HGT is essential for detecting a wide range of HGT events. These novel indicators of horizontal transfer will be widely useful in detecting HGT events linked to the evolution of important bacterial traits, such as antibiotic resistance and pathogenicity.
Hamady
Micah
M
Department of Chemistry and Biochemistry, University of Colorado, Boulder, CO 80309, USA. hamady@colorado.edu <hamady@colorado.edu>
Betterton
M D
MD
Knight
Rob
R
D019295:Computational Biology / Q000379:methods*; D003198:Computer Simulation; D022761:Gene Transfer, Horizontal; D015894:Genome, Human; D006801:Humans; D008390:Markov Chains; D008957:Models, Genetic; D008962:Models, Theoretical; D009711:Nucleotides / Q000737:chemistry*; D010802:Phylogeny; D015203:Reproducibility of Results; D012984:Software
D009711:Nucleotides