Originally published as Genetics Published Articles Ahead of Print on April 15, 2007.
Genetics, Vol. 176, 1741-1757, July 2007, Copyright © 2007
doi:10.1534/genetics.106.066233
Estimating the Number of Ancestral Lineages Using a Maximum-Likelihood Method Based on Rejection Sampling
Michael G. B. Blum1 and
Noah A. Rosenberg
Department of Human Genetics, University of Michigan, Ann Arbor, Michigan 48109
1 Corresponding author: Laboratoire TIMC, Faculté de Médecine, 38706 La Tronche, France.
E-mail: michael.blum{at}imag.fr
Estimating the number of ancestral lineages of a sample of DNA sequences at time t in the past can be viewed as a variation on the problem of estimating the time to the most recent common ancestor. To estimate the number of ancestral lineages, we develop a maximum-likelihood approach that takes advantage of a prior model of population demography, in addition to the molecular data summarized by the pattern of polymorphic sites. The method relies on a rejection sampling algorithm that is introduced for simulating conditional coalescent trees given a fixed number of ancestral lineages at time t. Computer simulations show that the number of ancestral lineages can be estimated accurately, provided that the number of mutations that occurred since time t is sufficiently large. The method is applied to 986 present-day human sequences located in hypervariable region 1 of the mitochondrion to estimate the number of ancestral lineages of modern humans at the time of potential admixture with the Neanderthal population. Our estimates support a view that the proportion of the modern population consisting of Neanderthal contributions must be relatively small, less than
5%, if the admixture happened as recently as 30,000 years ago.
Related articles in Genetics:
- ISSUE HIGHLIGHTS
Genetics 2007 176: NP.
[Full Text]
Copyright © 2007 by the Genetics Society of America.