Genetic Evidence for Long-Term Population Decline in a Savannah-Dwelling Primate: Inferences from a Hierarchical Bayesian Model

Jay F. Storz*, Mark A. Beaumont{dagger} and Susan C. Alberts{ddagger}

*Department of Ecology and Evolutionary Biology, University of Arizona;
{dagger}School of Animal and Microbial Sciences, University of Reading, Whiteknights, Reading, United Kingdom;
{ddagger}Department of Biology, Duke University, Durham, North Carolina


    Abstract
 TOP
 Abstract
 Introduction
 Methods
 Results
 Discussion
 Acknowledgements
 References
 
The purpose of this study was to test for evidence that savannah baboons (Papio cynocephalus) underwent a population expansion in concert with a hypothesized expansion of African human and chimpanzee populations during the late Pleistocene. The rationale is that any type of environmental event sufficient to cause simultaneous population expansions in African humans and chimpanzees would also be expected to affect other codistributed mammals. To test for genetic evidence of population expansion or contraction, we performed a coalescent analysis of multilocus microsatellite data using a hierarchical Bayesian model. Markov chain Monte Carlo (MCMC) simulations were used to estimate the posterior probability density of demographic and genealogical parameters. The model was designed to allow interlocus variation in mutational and demographic parameters, which made it possible to detect aberrant patterns of variation at individual loci that could result from heterogeneity in mutational dynamics or from the effects of selection at linked sites. Results of the MCMC simulations were consistent with zero variance in demographic parameters among loci, but there was evidence for a 10- to 20-fold difference in mutation rate between the most slowly and most rapidly evolving loci. Results of the model provided strong evidence that savannah baboons have undergone a long-term historical decline in population size. The mode of the highest posterior density for the joint distribution of current and ancestral population size indicated a roughly eightfold contraction over the past 1,000 to 250,000 years. These results indicate that savannah baboons apparently did not share a common demographic history with other codistributed primate species.


    Introduction
 TOP
 Abstract
 Introduction
 Methods
 Results
 Discussion
 Acknowledgements
 References
 
The occurrence of past vicariant events can often be inferred from geographically congruent patterns of phyletic subdivision within populations of codistributed species (e.g., Patton and da Silva 1998Citation ). Similarly, past environmental changes can be inferred from patterns of demographic history (i.e., population expansion or contraction) among codistributed species with similar or dissimilar ecological requirements (Storz and Beaumont 2002Citation ). Thus, information from multiple species about the direction, magnitude, and timing of past changes in population size can provide a basis for testing hypotheses about regional zoogeographic history.

Fortunately, genetic data can now be used to test such hypotheses, thanks to recent advances in theoretical population genetics and computational statistics (Emerson, Paradis, and Thébaud 2001Citation ). In particular, coalescent-based modeling has provided a powerful new means of estimating demographic and mutational parameters from patterns of DNA variation in contemporary populations (Hudson 1990Citation ; Donnelly and Tavaré 1995Citation ; Nordborg 2001Citation ). A typical approach is to infer parameters by comparing summary statistics calculated from an empirical data set with a distribution generated by Monte Carlo simulations of the coalescent process. When using sequence data to test for evidence of population expansion or contraction, one of the most commonly used summary statistics is Tajima's (1989a)Citation D, which measures whether the observed frequency distribution of segregating sites is consistent with a neutral-equilibrium model. But as pointed out by Hey and Harris (1999)Citation , Tajima's D and other related statistics generally provide inadequate summaries of the data for the purpose of testing complex parameter-rich models of demographic history (such as those considered in theories of human origins; Relethford 2001Citation ). For such models, likelihood analysis of the raw sample configuration generally permits a much broader scope of inference (Donnelly and Tavaré 1995Citation ; Stephens 2001Citation ). Specifically, coalescent simulations can be used to calculate the likelihood of the observed sample configuration under the stationary distribution of a specific demographic model, thereby permitting full use of the data in a maximum likelihood or Bayesian framework.

Inferences about the biogeographic history of East Africa during the late Quaternary are of special interest because of the implications for the evolution of anatomically modern humans and other hominoid primates. Archaeological evidence indicates that the expansion of the African human population coincided with the development and spread of Upper Paleolithic technology during the late Pleistocene, suggesting a causal role for cultural innovation in driving the expansion (Klein 1989Citation ). An alternative hypothesis is that the human population expansion in Africa was primarily attributable to climatically induced environmental change at or before the glacial-interglacial transition that marked the beginning of oxygen isotope stage three (~60,000 years BP; Rogers and Jorde 1995Citation ; Ambrose 1998Citation ). Such a broad-scale environmental change would be expected to affect other codistributed mammals (Rogers and Jorde 1995Citation ; Ambrose 1998Citation ; Harpending and Rogers 2000Citation ). As stated by Ruvolo (1996Citation , p. 211), "Any significant environmental or climatological event producing dramatic genetic effects within humans should have affected chimpanzees and gorillas, as well as a diverse array of other species."

Genetic inferences about the timing of the human population expansion remain highly controversial (Harpending and Rogers 2000Citation ; Wall and Przeworski 2000Citation ). Much of the debate centers on the effects of nonequilibrium demographic processes and the potentially confounding influence of selection at linked loci (Hey 1997Citation ; Fay and Wu 1999Citation ; Hey and Harris 1999Citation ; Wakeley 1999Citation ; Wall and Przeworski 2000Citation ). But available genetic data generally present a more consistent picture of late Pleistocene human history when attention is restricted to Africa alone. Data based on mtDNA variation generally suggest that the expansion of the African human population started ~30,000–130,000 years ago (Rogers and Harpending 1992Citation ; Hasegawa, Di Rienzo, and Wilson 1993Citation ; Sherry et al. 1994Citation ; Rogers 1995Citation , 1996Citation ; Rogers and Jorde 1995Citation ; Watson et al. 1997Citation ; Harpending et al. 1998Citation ; Excoffier and Schneider 1999Citation ; Harpending and Rogers 2000Citation ; Ingman et al. 2000Citation ). Similarly, mtDNA data for chimpanzees (Pan troglodytes) suggest that this species underwent a demographic expansion in equatorial Africa during this same approximate time interval (~20,000–60,000 years ago; Goldberg and Ruvolo 1997Citation ; Gagneux et al. 1999Citation ). The apparent similarity in the mismatch distributions of homologous mtDNA sequences from chimpanzees and African humans has been interpreted as evidence for a simultaneous response to a common environmental event in the late Pleistocene (Rogers and Jorde 1995Citation ; Ambrose 1998Citation ; Harpending and Rogers 2000Citation ). The purpose of this study was to evaluate this hypothesis by analyzing data from another codistributed African primate, the savannah baboon (Papio cynocephalus). We used a coalescent-based analysis of multilocus microsatellite data to test for evidence of population expansion or contraction. Specifically, we used a hierarchical Bayesian model to test for evidence that savannah baboons underwent a population expansion in concert with the hypothesized expansion of chimpanzees and African humans.


    Methods
 TOP
 Abstract
 Introduction
 Methods
 Results
 Discussion
 Acknowledgements
 References
 
Sampling and Genotyping
Whole blood was obtained from a total of 100 baboons that were sampled from a contiguous portion of the species range in East Africa (Altmann et al. 1996Citation ; Storz, Ramakrishnan, and Alberts 2002Citation ). Genomic DNA was isolated from blood samples according to the protocol described by Bruford et al. (1992)Citation or by using QIAamp extraction columns (Qiagen, Valencia, Calif.). The genetic analysis was based on a total of 14 autosomal microsatellites (10 di- and 4 tetranucleotide repeat loci). Baboon microsatellites were amplified using primers designed from flanking sequences of homologous loci in the human genome. Locus names and chromosomal assignments (in parentheses) are as follows: D7S503 (3), D13S159 (17), D6S311 (4), D6S271 (4), D16S402 (20), D4S431 (5), D2S141 (13), D16S420 (20), D11S925 (14), D17S791 (16), D5S1457 (6), D10S611 (9), D14S306 (7), and D4S243 (5). Chromosomal assignments of microsatellite loci (following karyotypic nomenclature of Cambefort, Mounie, and Colombies 1976Citation ) were based on linkage map data from Rogers et al. (2000)Citation . Although the markers were originally selected in humans on the basis of polymorphism levels, there is no correlation in levels of allelic diversity at homologous microsatellite loci between baboons and humans (Morin et al. 1998Citation ). Thus, there is no reason to expect that ascertainment bias has affected the observed frequency spectra (Beaumont 1999Citation ; Wakeley et al. 2001Citation ).

PCR protocols were reported previously (Storz, Ramakrishnan, and Alberts 2002Citation ). Allele sizes were quantified using an ABI 3700 automated sequencer and analyzed using GENESCAN software (PE Biosystems, Foster City, Calif.).

The Model
We consider the demographic history of a closed population that increases or decreases exponentially from an initial size N1 to the current size N0 over a time interval xa. The population size at any time x is given by


where time is increasing into the past. The data consist of frequency counts of chromosomes scored for microsatellite allele length at a number of unlinked loci. The loci are assumed to have independent genealogies and are evolving according to a strict single-step mutation model, with mutation rate µ. The distribution of allelic frequencies is assumed to arise from a genealogical history, G, consisting of the genealogy and the sequence of mutation laid down on the genealogy. Given the demographic parameters and mutation rate, results of coalescent theory can be used to calculate the probability of obtaining any genealogical history consistent with the allelic-frequency data, as described in Beaumont (1999)Citation .

We infer the parameters using the hierarchical Bayesian model of Storz and Beaumont (2002)Citation . In this model the mutational and demographic parameters are allowed to vary among loci. Although there are empirically justified reasons to expect interlocus variation in mutation rate (e.g., Chakraborty et al. 1997Citation ; Di Rienzo et al. 1998Citation ; Deka et al. 1999Citation ; Ellegren 2000aCitation ; Schlötterer 2000Citation ; Renwick et al. 2001Citation ), there would not seem to be any theoretically justified reason to expect interlocus variation in demographic parameters. This is because all genes have been transmitted through the same organismal pedigree and have therefore experienced the same demographic history. But the joint effects of linkage and selection can produce variation in effective population size across different regions of the genome (Hill and Robertson 1966Citation ; Barton 1995Citation ). Thus, allowing demographic parameters to vary makes it possible to detect aberrant patterns of variation at individual loci that may arise from genetic hitchhiking (Maynard-Smith and Haigh 1974Citation ) or background selection (Charlesworth B., Morgan, and Charlesworth D. 1993Citation ), either of which could make the demographic history of a single locus appear different from that of other unlinked loci. This approach should be viewed as a heuristic approximation to a full analysis in which the locus-specific effects of selection are modeled explicitly. The robustness of this modeling approach stems from the fact that likelihoods from separate loci are not combined in a multiplicative fashion, as in the model described by Beaumont (1999)Citation . Thus, aberrant loci will not strongly influence the results.

The hierarchical model is built as follows. There are separate parameters for each of k loci, giving a set of parameters {Phi} = {N0i, N1i, xai, µi}i=1,k. The parameters for each locus are assumed to be drawn from log-normal distributions, each with means (on a log scale) of = {MN0, MN1, Mµ, Mxa} and standard deviations (SDs) of = {VN0, VN1, Vµ, Vxa}. These distributions were then used to calculate the probability density p({Phi}|, ;th). Prior distributions for the means were assumed to be specified by normal distributions with means {alpha}N0, {alpha}N1, {alpha}µ, and {alpha}xa and SDs {sigma}N0, {sigma}N1, {sigma}µ, and {sigma}xa. Priors for the SDs were taken to be normal distributions truncated at zero with means ßN0, ßN1, ßµ, and ßxa and SDs {tau}N0, {tau}N1, {tau}µ, and {tau}xa. This set of parameters is designated as . The use of log-normal priors was motivated by the observation that posterior distributions for demographic and mutational parameters are generally skewed. It is therefore preferable to model the parameter distributions using log-transformed priors on a log scale (Beaumont 1999Citation ; Storz and Beaumont 2002Citation ) or skewed priors on a natural scale (Wilson and Balding 1998Citation ; Pritchard et al. 1999Citation ).

The above distributions were then used to calculate the probability density p(, ;th|). We used Metropolis-Hastings simulation, as described in Beaumont (1999)Citation and Storz and Beaumont (2002)Citation , to estimate the joint posterior distribution of the parameters and the genealogical history:


where D denotes the observed sample configuration (for a thorough discussion of this type of analysis, see Stephens and Donnelly 2000Citation and Stephens 2001Citation ). In the Metropolis-Hastings algorithm, parameters are updated by choosing trial values and deciding whether or not they should be accepted. The simulated sequence of parameter values tends to the joint posterior distribution of the parameters as the sequence length becomes large. The goal of the analysis is to report marginal posterior distributions for and , which can easily be obtained from the simulation output by analyzing the values of these parameters independently of the others.

For the hierarchical model, the acceptance probability was calculated as:


where the prime denotes a trial updated value chosen from an updating schedule given in Storz and Beaumont (2002)Citation . The last term in the expression is the Hastings correction, which ensures that there is no bias in the updating procedure toward any particular region of parameter space. The details of the Hastings correction are given in Beaumont (1999)Citation and Storz and Beaumont (2002)Citation . For P >= 1, the trial value was accepted, otherwise it was accepted with probability P.

Specification of the Priors
Simulations were based on the following assumptions. Prior means for the current and ancestral population sizes ({alpha}N0 and {alpha}N1) were set equal to a log10-transformed value of 5.9. This value was obtained by estimating the total area of the distributional range of savannah baboons (~3.0 x 106 km2), the population density (0.80 baboons/km2; Samuels and Altmann 1991Citation ), and the ratio of effective population size to census number (0.33; Storz, Ramakrishnan, and Alberts 2001Citation , 2002Citation ). We allowed for a large variation in this mean, {sigma}N0 = {sigma}N1 = 3, to allow for uncertainty in these estimates. The prior mean time since the population started changing in size, {alpha}xa, was taken to be 60,000 years BP, to match the hypothesized date of population expansion for African humans and chimpanzees (Rogers and Jorde 1995Citation ; Ambrose 1998Citation ; Harpending and Rogers 2000Citation ), giving {alpha}xa = 4.8. Although savannah baboons, chimpanzees, and archaic humans surely had somewhat different ecological requirements, the types of environmental catastrophe envisioned by previous authors (Rogers and Jorde 1995Citation ; Ambrose 1998Citation ; Harpending and Rogers 2000Citation ) are expected to have caused population bottlenecks (and subsequent expansions) in all three species. The SD of the mean around the hypothesized time of expansion ({sigma}xa = 2) was chosen such that more recent and more ancient dates also had support under the prior. We assumed a point prior of g = 10.15 years per generation (=average of male and female generation times; Storz, Ramakrishnan, and Alberts 2002Citation ). We chose priors for the mean mutation rate that supported values between 10-3 and 10-4: {alpha}µ = -3.3 and {sigma}µ = 0.25. This gives good support for a mean mutation rate among loci of 5 x 10-4, a value widely assumed in demographic models (e.g., Goldstein et al. 1995Citation ), but which also allows for the somewhat higher rates indicated by pedigree-based estimates for autosomal microsatellites in humans (Brinkmann et al. 1998Citation ; Ellegren 2000bCitation ).

The priors specified above are for the means of the demographic parameters and mutation rate, = {MN0, MN1, Mµ, Mxa}. To model the variation in parameters among loci, = {VN0, VN1, Vµ, Vxa}, we assumed in all cases that zero variation was most likely, so that ßN0 = ßN1 = ßµ, = ßxa = 0. In the case of the demographic parameters, we assumed that aberrant loci would be sufficiently rare that there would be relatively little variation among loci. Consequently, we chose {tau}N0 = {tau}N1 = {tau}xa = 0.05. By contrast, it is reasonable to expect substantial variation in the mutation rate among loci, so we chose a vague prior {tau}µ = 1, which allows the SD in the log10 mutation rate among loci to vary from 0 to ~2.

Implementation
To obtain convergence in a reasonable amount of time with 14 loci, a subset of 100 chromosomes was sampled for each locus. The chain was run for 4 x 109 steps, recording parameter values every 2 x 105 steps to give 20,000 draws from the posterior distribution. Five independent chains were run for each analysis (see Results). The output was checked for convergence using the Gelman-Rubin statistic calculated in CODA (Best, Cowles, and Vines 1995Citation ), as implemented in R (Ihaka and Gentleman 1996Citation ; http://www.r-project.org/). Posterior densities from individual runs were also examined to check for overall consistency in shape. The last half of each run was then combined to produce an overall set of 50,000 points. Each run took ca. 70 h on a 700 MHz Pentium (a small "Beowulf" cluster of machines was used to run chains in parallel). Density estimation was carried out using the program Locfit (Loader 1996Citation ), as implemented in R. This was used to estimate modes and highest probability density (HPD) limits (Gelman et al. 1995Citation ). The HPD limits specify points of equal probability density enclosing a region (or possibly disjoint regions in multimodal data), with probability equal to some specified value. For example, in a bivariate plot, the 90% HPD limit is represented as a contour enclosing a region (or regions) within which the parameter values are expected to lie with 90% probability. The strength of evidence of population decline versus population growth was assessed using Bayes factors, as described in Beaumont (1999)Citation . The Bayes factor is a ratio where the numerator is the posterior probability of one model (e.g., population contraction) divided by its prior probability and the denominator is the posterior probability of another model (e.g., population expansion) divided by its prior probability. The Bayes factor gives the relative likelihood of the two different models. Bayes factors greater than exp(2) are generally considered significant. The ratio of posterior probabilities can be estimated from the simulated chain by counting the proportion of iterations in which the population has expanded and then dividing this by the proportion of iterations in which it has contracted.


    Results
 TOP
 Abstract
 Introduction
 Methods
 Results
 Discussion
 Acknowledgements
 References
 
For the total sample of savannah baboons, the number of alleles per locus averaged 7.7 (range = 3–11) and expected heterozygosity averaged 0.678 (range = 0.271–0.830). No departures from Hardy-Weinberg equilibrium were observed at any loci (Storz, Ramakrishnan, and Alberts 2002Citation ). A significant degree of genetic subdivision would have been reflected by a heterozyogote deficit. No such Wahlund effect was detected in our population sample as the average inbreeding coefficient for all 14 loci (f = 0.011) was not significantly different from zero.

Variation in Parameters Among Loci
There was no evidence for substantial variation in the demographic parameters among loci, but there was evidence for variation in mutation rate. The five independent runs of the Markov chain Monte Carlo (MCMC) showed good convergence of the posterior distributions of = {VN0, VN1, Vµ, Vxa}, the SDs of the parameters across all 14 loci. Neither the Gelman-Rubin statistic nor its 0.975 quantile was greater than 1.02 for any parameter. The posterior distributions for VN0, VN1, and Vxa were very similar to the priors. The posterior density at zero was always high, and the estimated modes were at or near zero. Although we cannot rule out some variation in the demographic parameters, the results are consistent with zero variation. By contrast, there was strong evidence for interlocus variation in mutation rate (fig. 1 ). The estimated HPD limits of the prior distribution were 0–1.3 with a mode at zero. The estimated HPD limits of the posterior distribution of Vµ were 0.10–0.59 with a mode at 0.31, indicating that there is little support for Vµ = 0. The mode of 0.31 was consistent with a 10- to 20-fold difference in mutation rate between the most slowly and most rapidly evolving loci.



View larger version (14K):
[in this window]
[in a new window]
 
Fig. 1.—The posterior distribution (on a log10 scale) of the SD of mutation rate among loci (Vµ). The prior is shown as a dotted line

 
Convergence of MCMC Simulations
The five independent simulations show reasonable convergence for the means of the parameter values, although the degree of convergence for MN0 and Mxa was less than ideal. The Gelman-Rubin statistics (with 0.975 quantile in parentheses) for each parameter were MN0, 1.15 (1.39); MN1, 1.00 (1.01); Mµ, 1.00 (1.00); and Mxa, 1.11 (1.28). As a rule of thumb, values of the Gelman-Rubin statistic less than 1.2 indicate reasonable convergence (Gelman et al. 1995Citation ). The reason for poor convergence for MN0 and Mxa is that there was a thick tail for low values of MN0 and Mxa, and this was sampled to different degrees among the independent simulations. Differences in the degree of convergence corresponding to different values of the Gelman-Rubin statistic are illustrated in figure 2 .



View larger version (14K):
[in this window]
[in a new window]
 
Fig. 2.—Posterior distributions for current population size (MN0; top) and ancestral population size (MN1; bottom) estimated from five independent MCMC simulations

 
Evidence for Historical Population Decline
The data indicate that P. cynocephalus has undergone a pronounced population decline (fig. 3 ). The 90% HPD limits for the joint distribution of the ancestral and current population size lie well above the line of equality. The Bayes factor for contraction versus expansion was 9,999 (in only 5 out of 50,000 points was the current size greater than the ancestral size). The mode of the joint distribution was MN0 = 3.36, MN1 = 4.24, indicating an eightfold contraction in population size. The summaries of the marginal posterior distributions for each parameter are given in table 1 . Marginally, the ancestral effective population size was estimated to be between 5,000 and 60,000, and the current effective population size was estimated to be between 200 and 10,000.



View larger version (17K):
[in this window]
[in a new window]
 
Fig. 3.—Joint posterior distribution (on a log10 scale) of the mean current population size (MN0) and the mean ancestral population size (MN1) of the savannah baboon, P. cynocephalus. The 90%, 50%, and 10% HPD limits are plotted as continuous lines. The HPD limits for the prior are plotted as dotted lines. The diagonal line corresponds to MN0 = MN1

 

View this table:
[in this window]
[in a new window]
 
Table 1 Posterior and Prior Distributions of the Means of the Parameter Values Used in the Model

 
Timing of the Population Decline
The data were not very informative about the possible time scale of population decline. The posterior distribution of the time since the population started to decline (conditional on having declined rather than expanded) is plotted in figure 4 , and the summary of this distribution is given in table 1 . The distribution is skewed with a long tail to the left. The estimated mode for the time since the population started to decline was around 30,000 years, and it is unlikely to have occurred more than 250,000 years ago or more recently than ~1,000 years ago. It can be seen from figure 5 that the joint posterior distribution of MN0 and Mxa was highly correlated, and short times since the population has declined imply that the current effective size is quite small. If we had used more informative priors for the population size, the tail in Mxa would be reduced, supporting a more ancient date for the onset of the population decline. An additional caveat is that the information on times comes from the prior chosen for Mµ, and if higher (or lower) values had been chosen, the times would be shorter (or longer) in direct proportion. It can be seen from table 1 that there is negligible information on Mµ, and the posterior distribution is almost identical to the prior.



View larger version (12K):
[in this window]
[in a new window]
 
Fig. 4.—The posterior distribution (on a log10 scale) of the mean time since the population started to change in size (Mxa). The prior is shown as a dashed line

 


View larger version (16K):
[in this window]
[in a new window]
 
Fig. 5.—Joint posterior distribution (on a log10 scale) of the mean current population size (MN0) and the mean time since the population started to change in size (Mxa). The 90%, 50%, and 10% HPD limits are plotted as continuous lines. The HPD limits for the prior are plotted as dotted lines

 

    Discussion
 TOP
 Abstract
 Introduction
 Methods
 Results
 Discussion
 Acknowledgements
 References
 
Results of the coalescent analysis of microsatellite variation provide strong evidence that savannah baboons have undergone a historical decline in population size. The data suggest that the onset of the decline occurred in the late Pleistocene or early Holocene. A previous analysis of microsatellite allelic frequencies in savannah baboons revealed no evidence for a recent genetic bottleneck (Storz, Ramakrishnan, and Alberts 2002Citation ). The test used to detect a bottleneck was based on theory developed for neutral loci that relates the expected heterozygosity at mutation-drift equilibrium to the observed number of alleles in the population sample (Nei, Maruyama, and Chakraborty 1975Citation ; Watterson 1984Citation ; Maruyama and Fuerst 1985Citation ; Cornuet and Luikart 1996Citation ). Because allelic diversity is reduced more rapidly than heterozygosity during a bottleneck (Nei, Maruyama, and Chakraborty 1975Citation ), an abrupt reduction in population size is expected to produce a transient excess of heterozygosity, relative to the expected value at mutation-drift equilibrium. The fact that this test revealed no significant departure from mutation-drift equilibrium suggests that the savannah baboon population has remained relatively stable in recent history. This supports the results of the likelihood-based coalescent analysis, which indicated that the onset of population decline occurred well over 1,000 years ago. It thus appears that the size reduction of the baboon population is attributable to a long-term decline and not to a recent bottleneck.

The evidence for a relatively ancient onset of population decline also suggests that the inferred distribution of coalescence times cannot be attributed to a recent change in the degree of population subdivision (Wakeley 1999Citation ; Wakeley et al. 2001Citation ). A more recent onset of population decline would have been more difficult to distinguish from the effects of subdivision because shorter coalescence times could simply reflect rates of drift within partially isolated local demes rather than in the population as a whole.

Disentangling the Effects of Demography and Selection
Using genetic data to make inferences about demographic history is complicated by the potentially confounding influence of selection at linked sites (Nachman et al. 1998Citation ; Przeworski, Hudson, and Di Rienzo 2000Citation ). Departures from expectations of the standard neutral-equilibrium model can result from the direct or indirect effects of selection or from changes in the size or spatial contiguity of the population under study. For example, negative values of Tajima's D indicate an excess of low-frequency variants relative to null expectations and may be caused by population expansion or positive selection at linked sites. Positive values of D indicate an excess of intermediate-frequency variants relative to null expectations and may be caused by population contraction, population subdivision, or balancing selection at linked sites (Tajima 1989a,Citation 1989bCitation ). With regard to marker-based inferences about historical demography, the confounding influence of selection may be especially problematic in studies based on variation in mtDNA and the nonrecombining portion of the Y-chromosome. Because of the absence (or near absence; Awadalla, Eyre-Walker, and Maynard-Smith 1999Citation ) of recombination, positive selection at linked sites or purifying selection against mildly deleterious mutations could have produced the excess of low-frequency variants observed in mtDNA surveys of human populations (Excoffier 1990Citation ; Nachman et al. 1996Citation ; Wise, Sraml, and Easteal 1998Citation ).

Global surveys of human genetic variation have revealed profound differences in the frequency spectra of mitochondrial and nuclear sequence variation (Hey 1997Citation ). It seems likely that the observed discrepancy is attributable to the effects of selection on one or both classes of sequence variation. Whatever the underlying cause, the disparate patterns of variation have led to conflicting inferences about human demographic history. Harpending and Rogers (2000)Citation argued that mtDNA variation preserves a faithful record of human demographic history, whereas the frequency spectrum of nuclear sequence variation often reflects a history of balancing selection that has acted independently on multiple loci. According to Harpending and Rogers (2000)Citation , the pervasive effects of balancing selection may explain why nuclear sequence data are often characterized by an excess of intermediate-frequency variants relative to neutral-equilibrium expectations (see Przeworski, Hudson, and Di Rienzo 2000Citation ; Nachman 2001Citation ). By contrast, Wall and Przeworski (2000)Citation suggested that the consensus pattern of sequence variation across multiple, unlinked nuclear loci may often provide the most reliable basis for inferring aspects of human demographic history and that the frequency spectrum of mtDNA variation likely reflects a history of selective sweeps or purifying selection on mildly deleterious mutations. According to Wall and Przeworski (2000)Citation , the effects of either positive or negative selection may explain why mtDNA data typically exhibit an excess of low-frequency variants relative to neutral-equilibrium expectations (see also Hey 1997Citation ).

Although any given nuclear locus is potentially subject to the direct or indirect effects of selection, the potential for selection to confound marker-based inferences about demographic history can be minimized by combining data from multiple unlinked regions of the genome. Comparison among unlinked markers can be used to distinguish between the effects of selection (which are locus specific, at least to a degree determined by local recombination rate; Maynard-Smith and Haigh 1974Citation ; Kaplan, Hudson, and Langley 1989Citation ) and the effects of demographic processes (which are expected to have relatively uniform effects across the genome). One of the chief virtues of using microsatellite markers for the estimation of neutral parameters is that a large number of allelic length polymorphisms can be assayed across disparate regions of the genome. Moreover, data for humans indicate that recombination rates exhibit no significant association with levels of microsatellite polymorphism (Payseur and Nachman 2000Citation ).

But even if information is combined from many microsatellite loci, the effect of selection on a single locus can potentially bias the results of a multilocus analysis. In a model where demographic parameters are not allowed to vary among loci, aberrant patterns of variation at loci under selection can have a strong influence if likelihoods for individual loci are combined in a multiplicative fashion. In our hierarchical model, the likelihoods are combined in a more additive fashion. Thus, results will not be unduly affected by loci that deviate from the baseline level of neutral variation. One of the chief virtues of the hierarchical model is that it provides a means of identifying aberrant loci that could potentially bias the estimation of neutral parameters (cf. Beaumont and Nichols 1996Citation ; Vitalis, Dawson, and Boursot 2001Citation ).

Our results provide strong evidence for interlocus variation in some parameters. Given our prior assumption that heterogeneity in mutation rate is the most likely cause of such variation, we conclude that there is no evidence for substantial interlocus variation in demographic parameters. The posterior distributions most strongly support SDs close to zero and are generally indistinguishable from the priors. This suggests that there is insufficient power to detect differences among loci, most likely because each individual locus conveys relatively little information about demographic history.

Disentangling the Effects of Demography and Mutation
Compared with single-nucleotide polymorphism, standing levels of microsatellite variation generally will be less strongly affected by the interaction of linkage and selection because loci with high mutation rates are expected to undergo a more rapid return to equilibrium after a selective sweep (Wiehe 1998Citation ). But uncertainty about the mode of mutation presents more of a problem for microsatellite analysis. Convergent or parallel mutations that produce allelic size homoplasy can potentially obscure the signal of past demographic events (Di Rienzo et al. 1998Citation ; Gonser et al. 2000Citation ). As an example, consider the genealogical history of a sample under a model of population expansion. Looking backward in time, coalescence rates will be lowest near the branch tips (closer to the present when the population is large) and will be highest near the root of the tree (closer to the start date of the expansion when the population was small; Donnelly and Tavaré 1995Citation ). At any given locus, the inferred relationship between two randomly sampled alleles will depend on the depth of the tree (i.e., the number of generations since the most recent common ancestor). Because mutations are expected to occur independently along the branches leading from the root (i.e., common ancestor) to the sampled alleles, branch lengths should be proportional to the number of mutational steps separating the alleles. Thus, inferences about a past history of population expansion or contraction may be sensitive to the assumptions of the mutation model (Donnelly 1999Citation ).

The relative frequency of multistep mutations remains a subject of considerable debate (Di Rienzo et al. 1994Citation , 1998Citation ; Estoup and Cornuet 1999Citation ; Ellegren 2000bCitation ; Renwick et al. 2001Citation ). Pedigree data from humans suggest that >90%–95% of microsatellite mutations involve single-step changes in allele size, and the remainder typically involve a small number of steps (Brinkmann et al. 1998Citation ; Kayser et al. 2000Citation ). Moreover, analyses of Y-chromosome microsatellites based on the stepwise mutation model have provided strong evidence for demographic expansion of the human population, corroborating inferences from other markers (Pritchard et al. 1999Citation ; Thomson et al. 2000Citation ). If multistep mutational changes are infrequent relative to single-step changes, as suggested by the available pedigree evidence, then only a small number of loci should be affected. Our choice of priors for the mutational parameters allowed for a considerable degree of interlocus variation in µ. In a similar manner, the model accounts for some degree of variation in the mode of mutation (occasional multistep changes) by allowing interlocus variation in demographic parameters. Overall, the hierarchical model will tend to give less weight to loci that exhibit aberrant patterns of variation, thereby ensuring that our conclusions are not overly sensitive to violations of the strict stepwise mutation model.

Discordant Patterns of Demographic History Among Codistributed African Primates
The main objective of this analysis was to test for evidence that savannah baboons underwent a population expansion in concert with the hypothesized expansion of African human and chimpanzee populations in the late Pleistocene (Harpending and Rogers 2000Citation ). The rationale is that any type of environmental event sufficient to cause simultaneous population expansions in African humans and chimpanzees would also be expected to affect other codistributed mammals (Rogers and Jorde 1995Citation ; Ruvolo 1996Citation ; Goldberg and Ruvolo 1997Citation ; Ambrose 1998Citation ; Harpending and Rogers 2000Citation ). Although the microsatellite data indicate that savannah baboons actually underwent a pronounced population contraction during the late Pleistocene, our results do not necessarily refute the hypothesis that the expansions of African humans and chimpanzees were driven by a single environmental event. Several authors (Rogers and Jorde 1995Citation ; Ambrose 1998Citation ; Harpending and Rogers 2000Citation ) have suggested that the population expansions occurred as a result of recovery from some form of environmental catastrophe, such as the volcanic winter that followed the eruption of Mt. Toba (which occurred ~73,500 years BP; Rampino and Self 1992Citation , 1993Citation ). If the expansions were preceded by a population bottleneck (for example, due to a climatically induced range contraction into one or more refugia), the nature of the subsequent recovery would depend on whether the populations expanded from a single refugium (as suggested for African humans and chimpanzees; Ambrose 1998Citation ; Harpending and Rogers 2000Citation ) or several. As stated by Harpending and Rogers (2000Citation , p. 370): "... a few species that showed a different pattern would not refute the hypothesis of a Pleistocene bottleneck, but a few showing the same pattern would greatly strengthen it." Data from a more diverse array of codistributed mammals should eventually provide a clearer picture of the zoogeographic history of East Africa and the possible role of environmental catastrophes in the evolution of anatomically modern humans and other hominoid primates.


    Acknowledgements
 TOP
 Abstract
 Introduction
 Methods
 Results
 Discussion
 Acknowledgements
 References
 
We thank S. Palumbi, B. Payseur, and two anonymous reviewers for comments on the manuscript. This work was supported by NSF grant IBN-0002342 to S.C.A. and by NSF grant IBN-9985910 and Chicago Zoological Society funds to J. Altmann. The office of the President of Kenya and the Kenya Wildlife Service provided permission to work in Kenya. S.C.A. thanks R. Leakey, N. Kio, D. Western, J. Else, M. Ishakia, N. Rotich, J. Mwenda, C. Mlay, and the staff at Amboseli National Park for cooperation and assistance and J. Altmann, R. S. Mututua, S. Sayialel, and J. K. Warutere for the collection of blood samples. The members of the pastoralist communities of Amboseli and Longido and the Institute for Primate Research in Nairobi provided assistance and local sponsorship. J.F.S. acknowledges support from an NIH Postdoctoral Fellowship (F32 HL68487-01) and a Fellowship in Computational Molecular Biology from the Alfred P. Sloan Foundation and U.S. Department of Energy.


    Footnotes
 
Stephen Palumbi, Reviewing Editor

Keywords: baboon Bayesian inference coalescent demographic history Markov chain Monte Carlo microsatellite DNA Papio cynocephalus Back

Author for correspondence and reprints: Jay F. Storz, Department of Ecology and Evolutionary Biology, University of Arizona, Biosciences West, Tucson, Arizona 85721. storz{at}email.arizona.edu Back


    References
 TOP
 Abstract
 Introduction
 Methods
 Results
 Discussion
 Acknowledgements
 References
 

    Altmann J., S. C. Alberts, S. A. Haines, et al. (13 co-authors) 1996 Behavior predicts genetic structure in a wild primate group Proc. Natl. Acad. Sci. USA 93:5797-5801[Abstract/Free Full Text]

    Ambrose S. H., 1998 Late Pleistocene human population bottlenecks, volcanic winter, and differentiation of modern humans J. Hum. Evol 34:623-651[ISI][Medline]

    Awadalla P., A. Eyre-Walker, J. Maynard-Smith, 1999 Linkage disequilibrium and recombination in hominoid mitochondrial DNA Science 286:2524-2525[Abstract/Free Full Text]

    Barton N. H., 1995 Linkage and the limits to natural selection Genetics 140:821-841[Abstract/Free Full Text]

    Beaumont M. A., 1999 Detecting population expansion and decline using microsatellites Genetics 153:2013-2029[Abstract/Free Full Text]

    Beaumont M. A., R. A. Nichols, 1996 Evaluating loci for use in the genetic analysis of population structure Proc. R. Soc. Lond. B 263:1619-1626[ISI]

    Best N. G., M. K. Cowles, S. K. Vines, 1995 CODA manual version 0.30 MRC Biostatistics Unit, Cambridge, U.K

    Brinkman B., M. Klintschar, F. Neuhuber, J. Huhne, B. Rolf, 1998 Mutation rate in human microsatellites: influence of the structure and length of the tandem repeat Am. J. Hum. Genet 62:1408-1415[ISI][Medline]

    Bruford M. W., O. Hanotte, J. F. Y. Brookfield, T. Burke, 1992 Multi- and single-locus DNA fingerprinting Pp. 225–269 in A. R. Hoelzel, ed. Molecular analysis of populations: a practical approach. IRL Press, Oxford, U.K

    Cambefort Y., C. Mounie, P. Colombies, 1976 Chromosome banding patterns of Papio papio Ann. Genet 19:5-9[ISI][Medline]

    Chakraborty R., M. Kimmel, D. N. Stivers, L. J. Davison, R. Deka, 1997 Relative mutation rates at di-, tri-, and tetranucleotide microsatellite loci Proc. Natl. Acad. Sci. USA 94:1041-1046[Abstract/Free Full Text]

    Charlesworth B., M. T. Morgan, D. Charlesworth, 1993 The effect of deleterious mutations on neutral molecular variation Genetics 134:1289-1303[Abstract/Free Full Text]

    Cornuet J.-M., G. Luikart, 1996 Description and power analysis of two tests for detecting recent population bottlenecks from allele frequency data Genetics 144:2001-2014[Abstract/Free Full Text]

    Deka R., S. Guangyun, D. Smelser, et al. (6 co-authors) 1999 Rate and directionality of mutations and effects of allele size constraints at anonymous, gene-associated, and disease-causing trinucleotide loci Mol. Biol. Evol 16:1166-1177.[Abstract]

    Di Rienzo A., P. Donnelly, C. Toomajian, et al. (8 co-authors) 1998 Heterogeneity of microsatellite mutations within and between loci, and implications for human demographic histories Genetics 148:1269-1284.[Abstract/Free Full Text]

    Di Rienzo A., A. C. Peterson, J. C. Garza, A. M. Valdes, M. Slatkin, N. B. Freimer, 1994 Mutational processes of simple-sequence repeat loci in human populations Proc. Natl. Acad. Sci. USA 91:3166-3170[Abstract]

    Donnelly P., 1999 The coalescent and microsatellite variability Pp. 116–128 in D. B. Goldstein and C. Schlötterer, eds. Microsatellites: evolution and applications. Oxford University Press, Oxford, U.K

    Donnelly P., S. Tavaré, 1995 Coalescents and genealogical structure under neutrality Annu. Rev. Genet 29:410-421

    Ellegren H., 2000a. Heterogeneous mutation processes in human microsatellite DNA sequences Nat. Genet 24:400-402[ISI][Medline]

    ———. 2000b. Microsatellite mutations in the germline: implications for evolutionary inference Trends Genet 16:551-558[ISI][Medline]

    Emerson B. C., E. Paradis, C. Thébaud, 2001 Revealing the demographic histories of species using DNA sequences Trends. Ecol. Evol 16:707-716[ISI]

    Estoup A., J.-M. Cornuet, 1999 Microsatellite evolution: inferences from population data Pp. 49–65 in D. B. Goldstein and C. Schlötterer, eds. Microsatellites: evolution and applications. Oxford University Press, Oxford, U.K

    Excoffier L., 1990 Evolution of human mitochondrial DNA: evidence for departure from a pure neutral model of populations at equilibrium J. Mol. Evol 30:125-139[ISI][Medline]

    Excoffier L., S. Schneider, 1999 Why hunter-gatherer populations do not show signs of Pleistocene demographic expansions Proc. Natl. Acad. Sci. USA 96:10597-10602[Abstract/Free Full Text]

    Fay J. C., C.-I. Wu, 1999 A human population bottleneck can account for the discordance between patterns of mitochondrial and nuclear DNA variation Mol. Biol. Evol 16:1003-1005[Free Full Text]

    Gagneux P., C. Wills, U. Gerloff, D. Tautz, P. A. Morin, C. Boesch, B. Fruth, G. Hohmann, O. A. Ryder, D. S. Woodruff, 1999 Mitochondrial sequences show diverse evolutionary histories of African hominoids Proc. Natl. Acad. Sci. USA 96:5077-5082[Abstract/Free Full Text]

    Gelman A., J. B. Carlin, H. S. Stern, D. B. Rubin, 1995 Bayesian data analysis Chapman and Hall, London, U.K

    Goldberg T. L., M. Ruvolo, 1997 The geographic apportionment of mitochondrial genetic diversity in East African chimpanzees, Pan troglodytes schweinfurthii Mol. Biol. Evol 14:976-984[Abstract]

    Goldstein D. B., A. Ruiz Linares, L. L. Cavalli-Sforza, M. W. Feldman, 1995 Genetic absolute dating based on microsatellites and the origin of modern humans Proc. Natl. Acad. Sci. USA 92:6723-6727[Abstract]

    Gonser R., P. Donnelly, G. Nicholson, A. Di Rienzo, 2000 Microsatellite mutations and inferences about human demography Genetics 154:1793-1807[Abstract/Free Full Text]

    Harpending H. C., M. A. Batzer, M. Gurven, L. B. Jorde, A. R. Rogers, S. T. Sherry, 1998 Genetic traces of ancient demography Proc. Natl. Acad. Sci. USA 95:1961-1967[Abstract/Free Full Text]

    Harpending H. C., A. Rogers, 2000 Genetic perspectives on human origins and differentiation Annu. Rev. Genomics Hum. Genet 1:361-385[ISI][Medline]

    Hasegawa M., A. Di Rienzo, A. C. Wilson, 1993 Toward a more accurate time scale for the human mitochondrial DNA tree J. Mol. Evol 37:347-354[ISI][Medline]

    Hey J., 1997 Mitochondrial and nuclear genes present conflicting portraits of human origins Mol. Biol. Evol 14:166-172[Abstract]

    Hey J., E. Harris, 1999 Population bottlenecks and patterns of human polymorphism Mol. Biol. Evol 16:1423-1426[Free Full Text]

    Hill W. G., A. Robertson, 1966 The effect of linkage on the limits to artificial selection Genet. Res 8:269-294[ISI][Medline]

    Hudson R. R., 1990 Gene genealogies and the coalescent process Oxf. Surv. Evol. Biol 7:1-44

    Ihaka R., R. Gentleman, 1996 R: a language for data analysis and graphics J. Comput. Graph. Stat 5:299-314

    Ingman M., H. Kaessmann, S. Pääbo, U. Gyllensten, 2000 Mitochondrial genome variation and the origin of modern humans Nature 408:708-713[ISI][Medline]

    Kaplan N. L., R. R. Hudson, C. H. Langley, 1989 The "hitchhiking effect" revisited Genetics 123:887-899[Abstract/Free Full Text]

    Kayser M., L. Roewer, M. Hedman, et al. (14 co-authors) 2000 Characteristics and frequency of germline mutations at microsatellite loci from the human Y chromosome, as revealed by direct observation in father/son pairs Am. J. Hum. Genet 66:1580-1588[ISI][Medline]

    Klein R. G., 1989 The human career University of Chicago Press, Chicago

    Loader C. R., 1996 Local likelihood density estimation Ann. Stat 24:1602-1618[ISI]

    Maruyama T., P. A. Fuerst, 1985 Population bottlenecks and nonequilibrium models in population genetics. II. Number of alleles in a small population that was formed by a recent bottleneck Genetics 111:675-689[Abstract/Free Full Text]

    Maynard-Smith J., J. Haigh, 1974 The hitchhiking effect of a favorable gene Genet. Res 23:23-35[ISI][Medline]

    Morin P. A., P. Mahboubi, S. Wedel, J. Rogers, 1998 Rapid screening and comparison of human microsatellite markers in baboons: allele size is conserved, but allele number is not Genomics 53:12-20[ISI][Medline]

    Nachman M. W., 2001 Single nucleotide polymorphisms and recombination rate in humans Trends Genet 17:481-485[ISI][Medline]

    Nachman M. W., V. L. Bauer, S. L. Crowell, C. F. Aquadro, 1998 DNA variability and recombination rates at X-linked loci in humans Genetics 150:1133-1141[Abstract/Free Full Text]

    Nachman M. W., W. M. Brown, M. Stoneking, C. F. Aquadro, 1996 Nonneutral mitochondrial DNA variation in humans and chimpanzees Genetics 142:953-963[Abstract/Free Full Text]

    Nei M., T. Maruyama, R. Chakraborty, 1975 The bottleneck effect and genetic variability in populations Evolution 29:1-10[ISI]

    Nordborg M., 2001 Coalescent theory Pp. 179–212 in D. J. Balding, M. Bishop, and C. Cannings, eds. Handbook of statistical genetics. John Wiley and Sons, New York

    Patton J. L., M. N. F. da Silva, 1998 Rivers, refuges, and ridges: the geography of speciation of Amazonian mammals Pp. 202–213 in D. J. Howard and S. H. Berlocher, eds. Endless forms: species and speciation. Oxford University Press, Oxford, U.K

    Payseur B., M. W. Nachman, 2000 Microsatellite variation and recombination rate in the human genome Genetics 156:1285-1298[Abstract/Free Full Text]

    Pritchard J. K., M. T. Seielstad, A. Perez-Lezaun, M. W. Feldman, 1999 Population growth of human Y chromosomes: a study of Y chromosome microsatellites Mol. Biol. Evol 16:1791-1798[Abstract/Free Full Text]

    Przeworski M., R. R. Hudson, A. Di Rienzo, 2000 Adjusting the focus on human variation Trends Genet 16:296-302[ISI][Medline]

    Rampino M. R., S. Self, 1992 Volcanic winter and accelerated glaciation following the Toba super-eruption Nature 359:50-52[ISI]

    ———. 1993 Bottleneck in human evolution and the Toba eruption Science 262:1955.[Medline]

    Relethford J. H., 2001 Genetic history of the human species Pp. 813–846 in D. J. Balding, M. Bishop, and C. Cannings, eds. Handbook of statistical genetics. John Wiley and Sons, New York

    Renwick A., L. Davison, H. Spratt, J. P. King, M. Kimmel, 2001 DNA dinucleotide evolution in humans: fitting theory to facts Genetics 159:737-747[Abstract/Free Full Text]

    Rogers A. R., 1995 Genetic evidence for a Pleistocene population explosion Evolution 49:608-615[ISI]

    ———. 1996 Mitochondrial mismatch analysis is insensitive to the mutational process Mol. Biol. Evol 13:895-902[Abstract/Free Full Text]

    Rogers A. R., H. C. Harpending, 1992 Population growth makes waves in the distribution of pairwise genetic differences Mol. Biol. Evol 9:552-569[Abstract]

    Rogers A. R., L. B. Jorde, 1995 Genetic evidence on modern human origins Hum. Biol 67:1-36[ISI][Medline]

    Rogers J., M. C. Mahaney, S. M. Witte, et al. (17 co-authors) 2000 A genetic linkage map of the baboon (Papio hamadryas) genome based on human microsatellite polymorphisms Genomics 67:237-247[ISI][Medline]

    Ruvolo M., 1996 A new approach to studying modern human origins: hypothesis testing with coalescence time distributions Mol. Phylogenet. Evol 5:202-219[ISI][Medline]

    Samuels A., J. Altmann, 1991 Baboons of the Amboseli basin: demographic stability and change Int. J. Primatol 12:1-19[ISI]

    Schlötterer C., 2000 Evolutionary dynamics of microsatellite DNA Chromosoma 109:365-371[ISI][Medline]

    Sherry S., A. R. Rogers, H. C. Harpending, H. Soodyall, T. Jenkins, M. Stoneking, 1994 Mismatch distributions of mtDNA reveal recent human population expansions Hum. Biol 66:761-775[ISI][Medline]

    Stephens M., 2001 Inference under the coalescent Pp. 213–238 in D. J. Balding, M. Bishop, and C. Cannings, eds. Handbook of statistical genetics. John Wiley and Sons, New York

    Stephens M., P. Donnelly, 2000 Inference in molecular population genetics J. Roy. Stat. Soc. B 62:605-635[ISI]

    Storz J. F., M. A. Beaumont, 2002 Testing for genetic evidence of population expansion and contraction: an empirical analysis of microsatellite DNA variation using a hierarchical Bayesian model Evolution 56:154-166[ISI][Medline]

    Storz J. F., U. Ramakrishnan, S. C. Alberts, 2001 Determinants of effective population size for loci with different modes of inheritance J. Hered 92:497-502[Abstract/Free Full Text]

    ———. 2002 Genetic effective size of a wild primate population: influence of current and historical demography Evolution 56:817-829[ISI][Medline]

    Tajima F., 1989a. Statistical method for testing the neutral mutation hypothesis by DNA polymorphism Genetics 123:585-595[Abstract/Free Full Text]

    ———. 1989b. The effect of change in population size on DNA polymorphism Genetics 123:597-601[Abstract/Free Full Text]

    Thomson R., J. K. Pritchard, P. D. Shen, P. J. Oefner, M. W. Feldman, 2000 Recent common ancestry of human Y chromosomes: evidence from DNA sequence data Proc. Natl. Acad. Sci. USA 97:7360-7365[Abstract/Free Full Text]

    Vitalis R., K. Dawson, P. Boursot, 2001 Interpretation of variation across marker loci as evidence of selection Genetics 158:1811-1823[Abstract/Free Full Text]

    Wakeley J., 1999 Nonequilibrium migration in human history Genetics 153:1863-1871[Abstract/Free Full Text]

    Wakeley J., R. Nielsen, S. N. Liu-Cordero, K. Ardlie, 2001 The discovery of single-nucleotide polymorphisms—and inferences about human demographic history Am. J. Hum. Genet 69:1332-1347[ISI][Medline]

    Wall J. D., M. Przeworski, 2000 When did the human population size start increasing? Genetics 155:1865-1874[Abstract/Free Full Text]

    Watson E., P. Forster, M. Richards, N.-J. Bandelt, 1997 Mitochondrial footprints of human expansions in Africa Am. J. Hum. Genet 61:691-704[ISI][Medline]

    Watterson G. A., 1984 Allele frequencies after a bottleneck Theor. Popul. Biol 26:387-407[ISI]

    Wiehe T., 1998 The effect of selective sweeps on the variance of the allele distribution of a linked multiallele locus: hitchhiking of microsatellites Theor. Popul. Biol 53:272-283[ISI][Medline]

    Wilson I. J., D. J. Balding, 1998 Genealogical inference from microsatellite data Genetics 150:499-510[Abstract/Free Full Text]

    Wise C. A., M. Sraml, S. Easteal, 1998 Departure from neutrality at the mitochondrial NADH dehydrogenase subunit 2 gene in humans, but not in chimpanzees Genetics 148:409-421[Abstract/Free Full Text]

Accepted for publication July 17, 2002.