Originally published as Genetics Published Articles Ahead of Print on December 30, 2005.

Genetics, Vol. 172, 1411-1425, March 2006, Copyright © 2006
doi:10.1534/genetics.105.044917

Estimating Diversifying Selection and Functional Constraint in the Presence of Recombination

Department of Statistics, University of Oxford, Oxford OX1 3TG, United Kingdom

1 Corresponding author: Department of Statistics, 1 South Parks Rd., Oxford OX1 3TG, United Kingdom.
E-mail: daniel.wilson{at}sjc.ox.ac.uk

Models of molecular evolution that incorporate the ratio of nonsynonymous to synonymous polymorphism (dN/dS ratio) as a parameter can be used to identify sites that are under diversifying selection or functional constraint in a sample of gene sequences. However, when there has been recombination in the evolutionary history of the sequences, reconstructing a single phylogenetic tree is not appropriate, and inference based on a single tree can give misleading results. In the presence of high levels of recombination, the identification of sites experiencing diversifying selection can suffer from a false-positive rate as high as 90%. We present a model that uses a population genetics approximation to the coalescent with recombination and use reversible-jump MCMC to perform Bayesian inference on both the dN/dS ratio and the recombination rate, allowing each to vary along the sequence. We demonstrate that the method has the power to detect variation in the dN/dS ratio and the recombination rate and does not suffer from a high false-positive rate. We use the method to analyze the porB gene of Neisseria meningitidis and verify the inferences using prior sensitivity analysis and model criticism techniques.




This article has been cited by other articles:


Home page
Mol Biol EvolHome page
T. Miyake, N. Takebayashi, and D. E. Wolf
Possible Diversifying Selection in the Imprinted Gene, MEDEA, in Arabidopsis
Mol. Biol. Evol., April 1, 2009; 26(4): 843 - 857.
[Abstract] [Full Text] [PDF]


Home page
Mol Biol EvolHome page
M. Anisimova and C. Kosiol
Investigating Protein-Coding Sequence Evolution with Probabilistic Codon Substitution Models
Mol. Biol. Evol., February 1, 2009; 26(2): 255 - 271.
[Abstract] [Full Text] [PDF]


Home page
Brief BioinformHome page
W. Delport, K. Scheffler, and C. Seoighe
Models of coding sequence evolution
Brief Bioinform, January 1, 2009; 10(1): 97 - 109.
[Abstract] [Full Text] [PDF]


Home page
Proc R Soc BHome page
J. A Andres, L. S Maroja, and R. G Harrison
Searching for candidate speciation genes using a proteomic approach: seminal proteins in field crickets
Proc R Soc B, September 7, 2008; 275(1646): 1975 - 1983.
[Abstract] [Full Text] [PDF]


Home page
GeneticsHome page
J. Vieira, E. Teles, R. A. M. Santos, and C. P. Vieira
Recombination at Prunus S-Locus Region SLFL1 Gene
Genetics, September 1, 2008; 180(1): 483 - 491.
[Abstract] [Full Text] [PDF]


Home page
Mol Biol EvolHome page
A. M. Gunasekera, T. Wickramarachchi, D. E. Neafsey, I. Ganguli, L. Perera, P. H. Premaratne, D. Hartl, S. M. Handunnetti, P. V. Udagama-Randeniya, and D. F. Wirth
Genetic Diversity and Selection at the Plasmodium vivax Apical Membrane Antigen-1 (PvAMA-1) Locus in a Sri Lankan Population
Mol. Biol. Evol., April 1, 2007; 24(4): 939 - 947.
[Abstract] [Full Text] [PDF]


Home page
GeneticsHome page
C. T. T. Edwards, E. C. Holmes, O. G. Pybus, D. J. Wilson, R. P. Viscidi, E. J. Abrams, R. E. Phillips, and A. J. Drummond
Evolution of the Human Immunodeficiency Virus Envelope Gene Is Dominated by Purifying Selection
Genetics, November 1, 2006; 174(3): 1441 - 1453.
[Abstract] [Full Text] [PDF]