Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
. 2010 Apr 22;5(4):e10284.
doi: 10.1371/journal.pone.0010284.

Formulating a historical and demographic model of recent human evolution based on resequencing data from noncoding regions

Guillaume Laval  1 , Etienne Patin, Luis B Barreiro, Lluís Quintana-Murci
Affiliations

Formulating a historical and demographic model of recent human evolution based on resequencing data from noncoding regions

Guillaume Laval et al. PLoS One. .

Abstract

Background: Estimating the historical and demographic parameters that characterize modern human populations is a fundamental part of reconstructing the recent history of our species. In addition, the development of a model of human evolution that can best explain neutral genetic diversity is required to identify confidently regions of the human genome that have been targeted by natural selection.

Methodology/principal findings: We have resequenced 20 independent noncoding autosomal regions dispersed throughout the genome in 213 individuals from different continental populations, corresponding to a total of approximately 6 Mb of diploid resequencing data. We used these data to explore and co-estimate an extensive range of historical and demographic parameters with a statistical framework that combines the evaluation of multiple models of human evolution via a best-fit approach, followed by an Approximate Bayesian Computation (ABC) analysis. From a methodological standpoint, evaluating the accuracy of the parameter co-estimation allowed us to identify the most accurate set of statistics to be used for the estimation of each of the different historical and demographic parameters characterizing recent human evolution.

Conclusions/significance: Our results support a model in which modern humans left Africa through a single major dispersal event occurring approximately 60,000 years ago, corresponding to a drastic reduction of approximately 5 times the effective population size of the ancestral African population of approximately 13,800 individuals. Subsequently, the ancestors of modern Europeans and East Asians diverged much later, approximately 22,500 years ago, from the population of ancestral migrants. This late diversification of Eurasians after the African exodus points to the occurrence of a long maturation phase in which the ancestral Eurasian population was not yet diversified.

PubMed Disclaimer

Conflict of interest statement

Competing Interests: The authors have declared that no competing interests exist.

Figures

Figure 1
Figure 1. Minor allele and derived allele frequency spectra.
(A) Minor allele frequency (MAF) and (B) derived allele frequency (DAF) spectra computed by merging the 20 non coding autosomal DNA sequences. The expected MAF and DAF spectra (grey bars) were obtained assuming constant population sizes (Material and Methods). To focus on low frequency bins, the MAF spectrum display values lower than 35 counts in each continental population. To show the derived alleles that are fixed in each continental population, we arbitrarily removed intermediate bins in the DAF spectrum.
Figure 2
Figure 2. Sequenced-based summary statistics in Africans, Europeans and East-Asians.
Biplots of Tajima's D and Fu and Li's F* computed for each genomic region separately, in Africans (A), Europeans (B) and East-Asians (C). Significant Tajima's D values (P<0.05) are indicated in blue, in green for Fu and Li's F* only, and in red for both. Biplots of the number of haplotypes (K) and polymorphisms (S) computed for each genomic region separately in Africans (D), Europeans (E) and East-Asians (F). Significant K values (P<0.05) are indicated in blue, in green for S, and in red for both. The grey dots indicate the expected values of each genomic region simulated assuming a constant population size model (simulation procedure and significance of each region are described in the Materials and Methods section).
Figure 3
Figure 3. Model and parameter best-fitted estimations.
(A) Simulations considering different levels of replacement of archaic hominids by modern humans. We performed 8 sets of 105 simulations: one set for a replacement rate δ = 0, one for δ = 1, 3 sets for 0≤δ≤0.01, 0≤δ≤0.1 and 0≤δ≤0.5, and 3 sets for δ≥0.5, δ≥0.9 and δ≥0.99. For each of the 8 sets, we considered three models of ancestral migration (represented by black arrows): a residual ancestral migration rate (m0∼10−10), an ancestral migration rate with the same range (10−6 to 4×10−3) as m the current migration rate (represented by gray arrows), and an ancestral migration twice higher than m. Among the 24 models tested, the model assuming a complete replacement rate of archaic hominids (δ = 1) and a residual ancestral migration (m0∼10−10) exhibited the significantly highest ψ1 except when compared with the model assuming an almost complete replacement rate of archaic hominids (δ≥0.99). This best-fitted range of parameters (δ≥0.99 and m0∼10−10), indicated by the yellow/orange/white area (A), was therefore used to simulate the African expansion (B) and the non African bottleneck (C). We performed three sets of 105 simulations for the onset tA: 0≤tA≤25 Kyears, 25≤tA≤50 Kyears and 50≤tA≤75 Kyears. For each of the three sets, we considered 5 models of growth rate αA parameters; αA = 0, 0≤αA≤0.005, 0.005≤αA≤0.01, 0.01≤αA≤0.015 and 0.015≤αA≤0.02. Among the 15 models tested, the best-fitted ranges of parameters (ψ1 significantly higher than ψ1 of the constant size model αA = 0, P<0.01) are indicated by the yellow/orange/white area (B). Likewise, we performed 5 sets of 105 simulations assuming bottlenecks intensities βOoA, starting at the time of the out-of-Africa exodus (TOoA) and ending at the independent Neolithic expansions in Europe and east-Asia: βOoA = 1, 1≤βOoA≤2, 2≤βOoA≤20, 20≤βOoA≤40 and 40≤βOoA≤60. The best-fitted range of parameter (ψ1 significant higher than ψ1 of the constant size model βOoA = 1, P<0.01), indicated by the yellow/orange/white area (C), was obtained with the set of priors 2≤βOoA≤20. The distributions used are specified in Table S4.
Figure 4
Figure 4. Models of recent African origin involving different dispersal scenarios.
(A) General RAOEB model best fitting the data, with parameter ranges given in Table 2. This model assumes a single out-of-Africa dispersal followed by the European and East-Asian split. (B) RAOEB model involving two independent, concomitant dispersals out of Africa, each giving rise to Europeans and East-Asians. (C) RAOEB model involving two independent dispersals out of Africa occurring at different times, the earlier giving rise to Europeans. (D) RAOEB model involving two independent dispersals out of Africa occurring at different times, the earlier giving rise to East Asians. For models B–D, the ranges of parameters are the same as those given in Table 2. The alternative dispersal model B (two independent dispersals at the same time) was performed using a split of the two non Africans populations concomitant with the time of out-of-Africa exodus (TOoA) simulated with the same prior reported in Table 2. The two alternative dispersal models C and D (two independent dispersals at different times) were simulated using times for the first out-of-Africa exodus drawn from the first half of the prior distribution of TOoA (Table 2), while times for the second out-of-Africa exodus were drawn from the second half of the prior distribution of TOoA. (E) Posterior probability estimated for the 4 possible dispersal models represented in A, B, C, and D.
Figure 5
Figure 5. Approximate posterior distributions of historical and demographic parameters.
This figure gives the estimated ABC posterior distributions of the historical and demographic parameters (Table 3) using the RAOEB model (Figure 4A) with best-fitted priors (Table 2). Black lines represent the prior distributions and grey bars the posterior distributions. The times were translated into years using a generation time equal to 25 years. The posterior distributions of the parameters where the estimations were not validated by means of the accuracy evaluation procedure are not presented (i.e. NA and δ).
See this image and copyright information in PMC

Similar articles

See all similar articles

Cited by

See all "Cited by" articles

References

    1. Cavalli-Sforza LL, Menozzi P, Piazza A. Princeton: Princeton Univ. Press; 1994. The History and Geography of Human Genes.
    1. Excoffier L. Human demographic history: refining the recent African origin model. Curr Opin Genet Dev. 2002;12:675–682. - PubMed
    1. Garrigan D, Hammer MF. Reconstructing human origins in the genomic era. Nat Rev Genet. 2006;7:669–680. - PubMed
    1. Harpending HC, Batzer MA, Gurven M, Jorde LB, Rogers AR, et al. Genetic traces of ancient demography. Proc Natl Acad Sci USA. 1998;95:1961–1967. - PMC - PubMed
    1. Stringer CB. The emergence of modern humans. Sci Am. 1990;263:98–104. - PubMed

Publication types

Substances