Estimating absolute rates of synonymous and nonsynonymous nucleotide substitution in order to characterize natural selection and date species divergences.

The rate of molecular evolution can vary among lineages. Sources of this variation have differential effects on synonymous and nonsynonymous substitution rates. Changes in effective population size or patterns of natural selection will mainly alter nonsynonymous substitution rates. Changes in generation length or mutation rates are likely to have an impact on both synonymous and nonsynonymous substitution rates. By comparing changes in synonymous and nonsynonymous rates, the relative contributions of the driving forces of evolution can be better characterized. Here, we introduce a procedure for estimating the chronological rates of synonymous and nonsynonymous substitutions on the branches of an evolutionary tree. Because the widely used ratio of nonsynonymous and synonymous rates is not designed to detect simultaneous increases or simultaneous decreases in synonymous and nonsynonymous rates, the estimation of these rates rather than their ratio can improve characterization of the evolutionary process. With our Bayesian approach, we analyze cytochrome oxidase subunit I evolution in primates and infer that nonsynonymous rates have a greater tendency to change over time than do synonymous rates. Our analysis of these data also suggests that rates have been positively correlated.

[1]  J. Thompson,et al.  CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. , 1994, Nucleic acids research.

[2]  M. Nei,et al.  Positive Darwinian selection after gene duplication in primate ribonuclease genes. , 1998, Proceedings of the National Academy of Sciences of the United States of America.

[3]  N. Goldman,et al.  A codon-based model of nucleotide substitution for protein-coding DNA sequences. , 1994, Molecular biology and evolution.

[4]  A. Lapedes,et al.  Timing the ancestor of the HIV-1 pandemic strains. , 2000, Science.

[5]  Michael J. Sanderson,et al.  R8s: Inferring Absolute Rates of Molecular Evolution, Divergence times in the Absence of a Molecular Clock , 2003, Bioinform..

[6]  R. Adkins,et al.  Evolution of eutherian cytochrome c oxidase subunit II: heterogeneous rates of protein evolution and altered interaction with cytochrome c. , 1996, Molecular biology and evolution.

[7]  M. Sanderson.  Estimating absolute rates of molecular evolution and divergence times: a penalized likelihood approach. , 2002, Molecular biology and evolution.

[8]  T. Ohta,et al.  Theoretical study of near neutrality. I. Heterozygosity and rate of mutant substitution. , 1990, Genetics.

[9]  W. Li,et al.  Evidence for higher rates of nucleotide substitution in rodents than in man. , 1985, Proceedings of the National Academy of Sciences of the United States of America.

[10]  T. Ohta,et al.  Synonymous and nonsynonymous substitutions in mammalian genes and the nearly neutral theory. , 1995, Journal of molecular evolution.

[11]  B S Weir,et al.  Testing for equality of evolutionary rates. , 1992, Genetics.

[12]  Z. Yang,et al.  Estimation of primate speciation dates using local molecular clocks. , 2000, Molecular biology and evolution.

[13]  L. Chao,et al.  THE MOLECULAR CLOCK AND THE RELATIONSHIP BETWEEN POPULATION SIZE AND GENERATION TIME , 1993, Evolution; international journal of organic evolution.

[14]  M. Kimura,et al.  The neutral theory of molecular evolution. , 1983, Scientific American.

[15]  M. I. Lomax,et al.  Molecular Evolution of Cytochrome c Oxidase Subunit IV: Evidence for Positive Selection in Simian Primates , 1997, Journal of Molecular Evolution.

[16]  S. Tavaré,et al.  Using the fossil record to estimate the age of the last common ancestor of extant primates , 2002, Nature.

[17]  G. Glazko,et al.  Estimation of divergence times for major lineages of primate species. , 2003, Molecular biology and evolution.

[18]  J. Huelsenbeck,et al.  A compound poisson process for relaxing the molecular clock. , 2000, Genetics.

[19]  David B. Dunson,et al.  Bayesian Data Analysis , 2010 .

[20]  Robert V. Hogg,et al.  Introduction to Mathematical Statistics. , 1966 .

[21]  Hirohisa Kishino,et al.  Divergence time and evolutionary rate estimation with multilocus data. , 2002, Systematic biology.

[22]  J. Felsenstein.  Evolutionary trees from DNA sequences: A maximum likelihood approach , 2005, Journal of Molecular Evolution.

[23]  Ronald M. Adkins,et al.  Evolution of the primate cytochrome c oxidase subunit II gene , 1994, Journal of Molecular Evolution.

[24]  Andrew Rambaut,et al.  Estimating the rate of molecular evolution: incorporating non-contemporaneous sequences into maximum likelihood phylogenies , 2000, Bioinform..

[25]  M. Nei,et al.  Simple methods for estimating the numbers of synonymous and nonsynonymous nucleotide substitutions. , 1986, Molecular biology and evolution.

[26]  Stéphane Aris-Brosou,et al.  Effects of models of rate evolution on estimation of divergence dates with special reference to the metazoan 18S ribosomal RNA phylogeny. , 2002, Systematic biology.

[27]  Jack Porter,et al.  Effciency of Covariance Matrix Estimators for Maximum Likelihood Estimation , 2002 .

[28]  Michael J. Sanderson,et al.  A Nonparametric Approach to Estimating Divergence Times in the Absence of Rate Constancy , 1997 .

[29]  Z. Yang,et al.  Estimating synonymous and nonsynonymous substitution rates under realistic evolutionary models. , 2000, Molecular biology and evolution.

[30]  C. Wiuf,et al.  A codon-based model designed to describe lentiviral evolution. , 1998, Molecular biology and evolution.

[31]  W. Bruno,et al.  Performance of a divergence time estimation method under a probabilistic model of rate evolution. , 2001, Molecular biology and evolution.

[32]  Simon Easteal,et al.  Evolutionary Rate Acceleration of Cytochrome c Oxidase Subunit I in Simian Primates , 2000, Journal of Molecular Evolution.

[33]  T. Ohta.  Population size and rate of evolution , 1972, Journal of Molecular Evolution.

[34]  T. Andrews,et al.  Accelerated Evolution of Cytochrome b in Simian Primates: Adaptive Evolution in Concert with Other Mitochondrial Proteins? , 1998, Journal of Molecular Evolution.

[35]  A. Purvis.  A composite estimate of primate phylogeny. , 1995, Philosophical transactions of the Royal Society of London. Series B, Biological sciences.

[36]  S. Muse,et al.  A likelihood approach for comparing synonymous and nonsynonymous nucleotide substitution rates, with application to the chloroplast genome. , 1994, Molecular biology and evolution.

[37]  T. Ohta.  Very slightly deleterious mutations and the molecular clock , 2005, Journal of Molecular Evolution.

[38]  Z. Yang,et al.  Likelihood ratio tests for detecting positive selection and application to primate lysozyme evolution. , 1998, Molecular biology and evolution.

[39]  M. Goodman,et al.  Molecular evolution of cytochrome c oxidase subunit I in primates: is there coevolution between mitochondrial and nuclear genomes? , 2000, Molecular phylogenetics and evolution.

[40]  O. Ryder,et al.  Human and ape molecular clocks and constraints on paleontological hypotheses. , 2001, The Journal of heredity.

[41]  Ziheng Yang,et al.  Codon-substitution models to detect adaptive evolution that account for heterogeneous selective pressures among site classes. , 2002, Molecular biology and evolution.

[42]  Ziheng Yang,et al.  Statistical methods for detecting molecular adaptation , 2000, Trends in Ecology & Evolution.

[43]  D. Wildman,et al.  Episodic positive selection in ape cytochrome c oxidase subunit IV. , 2002, Molecular biology and evolution.

[44]  M. Goodman,et al.  Molecular evolution of the COX7A gene family in primates. , 1999, Molecular biology and evolution.

[45]  D. Wildman,et al.  Molecular evolution of aerobic energy metabolism in primates. , 2001, Molecular phylogenetics and evolution.

[46]  R. Nielsen,et al.  Codon-substitution models for detecting molecular adaptation at individual sites along specific lineages. , 2002, Molecular biology and evolution.

[47]  H. Kishino,et al.  Estimating the rate of evolution of the rate of molecular evolution. , 1998, Molecular biology and evolution.

[48]  A. Rodrigo,et al.  Reconstructing genealogies of serial samples under the assumption of a molecular clock using serial-sample UPGMA. , 2000, Molecular biology and evolution.

[49]  W. K. Hastings,et al.  Monte Carlo Sampling Methods Using Markov Chains and Their Applications , 1970 .

[50]  M. Hasegawa,et al.  Interordinal relationships and timescale of eutherian evolution as inferred from mitochondrial genome data. , 2000, Gene.

[51]  Y. Ina,et al.  New methods for estimating the numbers of synonymous and nonsynonymous substitutions , 1995, Journal of Molecular Evolution.