Bayes Estimation of Species Divergence Times and Ancestral Population Sizes Using DNA Sequences From Multiple Loci

Corresponding author: Ziheng Yang, Darwin Bldg., Gower St., London WC1E 6BT, England., z.yang{at}ucl.ac.uk (E-mail)

Communicating editor: N. TAKAHATA

The effective population sizes of ancestral as well as modern species are important parameters in models of population genetics and human evolution. The commonly used method for estimating ancestral population sizes, based on counting mismatches between the species tree and the inferred gene trees, is highly biased as it ignores uncertainties in gene tree reconstruction. In this article, we develop a Bayes method for simultaneous estimation of the species divergence times and current and ancestral population sizes. The method uses DNA sequence data from multiple loci and extracts information about conflicts among gene tree topologies and coalescent times to estimate ancestral population sizes. The topology of the species tree is assumed known. A Markov chain Monte Carlo algorithm is implemented to integrate over uncertain gene trees and branch lengths (or coalescence times) at each locus as well as species divergence times. The method can handle any species tree and allows different numbers of sequences at different loci. We apply the method to published noncoding DNA sequences from the human and the great apes. There are strong correlations between posterior estimates of speciation times and ancestral population sizes. With the use of an informative prior for the human-chimpanzee divergence date, the population size of the common ancestor of the two species is estimated to be 20,000, with a 95% credibility interval (8000, 40,000). Our estimates, however, are affected by model assumptions as well as data quality. We suggest that reliable estimates have yet to await more data and more realistic models.

This article has been cited by other articles:

				R. Zhou, K. Zeng, W. Wu, X. Chen, Z. Yang, S. Shi, and C.-I Wu Population Genetics of Speciation in Nonmodel Organisms: I. Ancestral Polymorphism in Mangroves Mol. Biol. Evol., December 1, 2007; 24(12): 2746 - 2754. [Abstract] [Full Text] [PDF]

				S. C. Leman, M. K. Uyenoyama, M. Lavine, and Y. Chen The evolutionary forest algorithm Bioinformatics, August 1, 2007; 23(15): 1962 - 1968. [Abstract] [Full Text] [PDF]

				S. V. Edwards, L. Liu, and D. K. Pearl High-resolution species trees without concatenation PNAS, April 3, 2007; 104(14): 5936 - 5941. [Abstract] [Full Text] [PDF]

				J. F. Baines and B. Harr Reduced X-Linked Diversity in Derived Populations of House Mice Genetics, April 1, 2007; 175(4): 1911 - 1921. [Abstract] [Full Text] [PDF]

				J. Hey and R. Nielsen Integration within the Felsenstein equation for improved Markov chain Monte Carlo methods in population genetics PNAS, February 20, 2007; 104(8): 2785 - 2790. [Abstract] [Full Text] [PDF]

				J. Charlesworth and A. Eyre-Walker The Rate of Adaptive Evolution in Enteric Bacteria Mol. Biol. Evol., July 1, 2006; 23(7): 1348 - 1356. [Abstract] [Full Text] [PDF]

				M. Slatkin and J. L. Pollack The Concordance of Gene Trees and Species Trees at Two Linked Loci Genetics, March 1, 2006; 172(3): 1979 - 1984. [Abstract] [Full Text] [PDF]

				Z. Yang and B. Rannala Bayesian Estimation of Species Divergence Times Under a Molecular Clock Using Multiple Fossil Calibrations with Soft Bounds Mol. Biol. Evol., January 1, 2006; 23(1): 212 - 226. [Abstract] [Full Text] [PDF]

				H. E. Walsh, I. L. Jones, and V. L. Friesen A Test of Founder Effect Speciation Using Multiple Loci in the Auklets (Aethia spp.) Genetics, December 1, 2005; 171(4): 1885 - 1894. [Abstract] [Full Text] [PDF]

				L. Y. Yampolsky, F. A. Kondrashov, and A. S. Kondrashov Distribution of the strength of selection against amino acid replacements in human proteins Hum. Mol. Genet., November 1, 2005; 14(21): 3191 - 3201. [Abstract] [Full Text] [PDF]

				S. C. Leman, Y. Chen, J. E. Stajich, M. A. F. Noor, and M. K. Uyenoyama Likelihoods From Summary Statistics: Recent Divergence Between Species Genetics, November 1, 2005; 171(3): 1419 - 1436. [Abstract] [Full Text] [PDF]

				P. D. Keightley, G. V. Kryukov, S. Sunyaev, D. L. Halligan, and D. J. Gaffney Evolutionary constraints in conserved nongenic sequences of mammals Genome Res., October 1, 2005; 15(10): 1373 - 1378. [Abstract] [Full Text] [PDF]

				T.-K. Seo, H. Kishino, and J. L. Thorne Incorporating gene-specific variation when inferring and evaluating optimal evolutionary tree topologies from multilocus sequence data PNAS, March 22, 2005; 102(12): 4436 - 4441. [Abstract] [Full Text] [PDF]

				J. Hey and R. Nielsen Multilocus Methods for Estimating Population Sizes, Migration Rates and Divergence Time, With Applications to the Divergence of Drosophila pseudoobscura and D. persimilis Genetics, June 1, 2004; 167(2): 747 - 760. [Abstract] [Full Text] [PDF]