Cookies on this website
We use cookies to ensure that we give you the best experience on our website. If you click 'Continue' we'll assume that you are happy to receive all cookies and you won't see this message again. Click 'Find out more' for information on how to change your cookie settings.

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.

Original publication

DOI

10.1534/genetics.105.044917

Type

Journal article

Journal

Genetics

Publication Date

03/2006

Volume

172

Pages

1411 - 1425

Addresses

Department of Statistics, University of Oxford, Oxford OX1 3TG, United Kingdom. daniel.wilson@sjc.ox.ac.uk

Keywords

Neisseria meningitidis, Porins, Models, Statistical, Bayes Theorem, Genetics, Population, Recombination, Genetic, Haplotypes, Polymorphism, Genetic, Models, Genetic, Computer Simulation, Selection, Genetic