PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of nihpaAbout Author manuscriptsSubmit a manuscriptHHS Public Access; Author Manuscript; Accepted for publication in peer reviewed journal;
 
Nature. Author manuscript; available in PMC 2017 March 21.
Published in final edited form as:
PMCID: PMC5161557
NIHMSID: NIHMS798259

The Simons Genome Diversity Project: 300 genomes from 142 diverse populations

Swapan Mallick,1,2,3 Heng Li,2,* Mark Lipson,1,* Iain Mathieson,1,* Melissa Gymrek,2,4,5 Fernando Racimo,6 Mengyao Zhao,1,2,3 Niru Chennagiri,1,2,3 Susanne Nordenfelt,1,2,3 Arti Tandon,1,2 Pontus Skoglund,1,2 Iosif Lazaridis,1,2 Sriram Sankararaman,1,2 Qiaomei Fu,1,2,7 Nadin Rohland,1,2 Gabriel Renaud,8 Yaniv Erlich,9,10,11 Thomas Willems,9,12 Carla Gallo,13 Jeffrey P. Spence,14 Yun S. Song,15,16,17 Giovanni Poletti,13 Francois Balloux,18 George van Driem,19 Peter de Knijff,20 Irene Gallego Romero,21 Aashish R. Jha,22 Doron M. Behar,23 Claudio M. Bravi,24 Cristian Capelli,25 Tor Hervig,26 Andres Moreno-Estrada,27 Olga L. Posukh,28,29 Elena Balanovska,30 Oleg Balanovsky,30,31,32 Sena Karachanak-Yankova,33 Hovhannes Sahakyan,23,24 Draga Toncheva,33 Levon Yepiskoposyan,34 Chris Tyler-Smith,35 Yali Xue,35 M. Syafiq Abdullah,36 Andres Ruiz-Linares,37 Cynthia M. Beall,38 Anna Di Rienzo,22 Choongwon Jeong,22 Elena B. Starikovskaya,39 Ene Metspalu,23,40 Jüri Parik,23 Richard Villems,23,40,41 Brenna M. Henn,42 Ugur Hodoglugil,43 Robert Mahley,44 Antti Sajantila,45 George Stamatoyannopoulos,46 Joseph T. S. Wee,47 Rita Khusainova,48,49 Elza Khusnutdinova,48,49 Sergey Litvinov,23,48,49 George Ayodo,50 David Comas,51 Michael Hammer,52 Toomas Kivisild,23,53 William Klitz,6 Cheryl Winkler,54 Damian Labuda,55 Michael Bamshad,56 Lynn B. Jorde,57 Sarah A. Tishkoff,58 W. Scott Watkins,59 Mait Metspalu,23 Stanislav Dryomov,39,60 Rem Sukernik,39,61 Lalji Singh,62,63 Kumarasamy Thangaraj,62 Svante Pääbo,8 Janet Kelso,8 Nick Patterson,2 and David Reich1,2,3

Abstract

We report the Simons Genome Diversity Project (SGDP) dataset: high quality genomes from 300 individuals from 142 diverse populations. These genomes include at least 5.8 million base pairs that are not present in the human reference genome. Our analysis reveals key features of the landscape of human genome variation, including that the rate of accumulation of mutations has accelerated by about 5% in non-Africans compared to Africans since divergence. We show that the ancestors of some pairs of present-day human populations were substantially separated by 100,000 years ago, well before the archaeologically attested onset of behavioral modernity. We also demonstrate that indigenous Australians, New Guineans and Andamanese do not derive substantial ancestry from an early dispersal of modern humans; instead, their modern human ancestry is consistent with coming from the same source as that in other non-Africans.

To obtain a complete picture of human diversity, it is necessary to sequence the genomes of many individuals from diverse locations. To date, the largest whole-genome sequencing survey, the 1000 Genomes Project, analyzed 26 populations of European, East Asian, South Asian, American, and sub-Saharan African ancestry1. However, this and most other sequencing studies have focused on demographically large populations. Such studies tend to ignore smaller populations that are also important for understanding human diversity. In addition, many of these studies have sequenced genomes to only 4–6-fold coverage. Here, we report the Simons Genome Diversity Project (SGDP): deep genome sequences of 300 individuals from 142 populations chosen to span much of human genetic, linguistic, and cultural variation (Supplementary Data Table 1).

Data set and catalog of novel variants

We sequenced the samples to an average coverage of 43-fold (range 34–83 fold) at Illumina Ltd.; almost all samples (278) were prepared using the same PCR-free library preparation2. We aligned reads to the human reference genome hs37d5/hg19 using BWA-MEM (BWA-0.7.12)3 (Supplementary Information section 1). We genotyped each sample separately using the Genome Analysis Toolkit (GATK)4, with a modification to eliminate bias toward genotypes matching the reference (Supplementary Information section 1). We developed a filtering procedure that generates a sample-specific mask. At “filter level 1” which we recommend for most analyses, we retain an average of 2.13 Gb of sequence per sample and identify 34.4 million single nucleotide polymorphisms (SNPs) and 2.1 million insertion/deletion polymorphisms (indels) (Supplementary Information section 2). We have made the GATK-processed data available in a file small enough to download by FTP, along with software to analyze these data (Supplementary Information section 3). The SGDP dataset highlights the incompleteness of current catalogs of human variation, with the fraction of heterozygous positions not discovered by the 1000 Genomes Project being 11% in the KhoeSan and 5% in New Guineans and Australians (Extended Data Fig. 1; Supplementary Data Table 1). We used FermiKit5 to map short reads against each other, store the assemblies in a compressed form that retains all the information required for polymorphism discovery and analysis, and identified SNPs by comparing against the human reference. We find that FermiKit has comparable sensitivity and specificity to GATK for SNP discovery and genotyping, and is more accurate for indels (Supplementary Information section 4). FermiKit also identified 5.8 Mb of contigs that are present in the SGDP but absent in the human reference genome presumably because they are deleted there; these contigs which we have made publicly available can be used as “decoys” to improve read mapping (Supplementary Information section 5). Finally, we called copy number variants6 and used lobSTR7,8 to genotype 1.6 million short tandem repeats (STRs) (Supplementary Information section 6). The high quality of the STR genotypes (r2=0.92 to capillary sequencing calls) is evident from their accurate reconstruction of population relationships, even for difficult-to-genotype mononucleotide repeats (Extended Data Fig. 2).

The structure of human genetic diversity

To obtain an overview of population relationships, we carried out ADMIXTURE9 (Extended Data Fig. 3) and principal component analysis10 (Extended Data Fig. 4a). We also built neighbor-joining trees based on pairwise divergence per nucleotide (Fig. 1a) and FST (Extended Data Fig. 4b) whose topologies are consistent with previous findings that the deepest splits among human populations are among Africans. We computed heterozygosity – the proportion of diallelic genotypes per base pair – and recapitulate previous findings that the highest genetic diversity is found in sub-Saharan Africa and that there is a much lower ratio of X-to-autosome diversity in non-Africans than in Africans (Fig. 1b)11. A surprise is that African “Pygmy” hunter-gatherers have reduced X-to-autosome diversity ratios relative to all other sub-Saharan Africans. This pattern remains even after we remove the third of chromosome X known to be subject to the strongest natural selection, suggesting that the finding is driven by demographic history rather than by natural selection (Supplementary Information section 7). It has been suggested that the reduced X-to-autosome heterozygosity ratio in non-Africans is due to ongoing male-driven admixture11,12. Male non-Pygmy admixture into Pygmies is well-documented13,14, so this process could explain these findings.

Figure 1
Genetic variation in the SGDP

Comparisons of ancient to present-day human genomes have shown that all non-Africans today possess Neanderthal ancestry15 with more in eastern non-Africans16,17, and that Australo-Melanesians and to a lesser extent other eastern non-Africans possess Denisovan ancestry1820. However, these studies only analyzed genomes from a handful of populations. We computed statistics informative about Neanderthal and Denisovan ancestry and provide a fine-scale view of these ancestry distributions worldwide (Fig. 1c,d; Supp. Data Table 1; Supplementary Information section 8). We do not detect any population with a higher proportion of Neanderthal ancestry than is present in East Asians. However, we do find suggestive evidence of an excess of Denisovan ancestry in some South Asians compared to other Eurasians. This signal may not have been detected before because earlier surveys of archaic introgression largely excluded South Asians (Fig. 1d; Supp. Data Table 1).

The time course of human population separation

We studied demographic history by leveraging the fact that variation across the genome in divergent sites per base pair can be used to reconstruct population size changes and separations. We used the Pairwise Sequential Markovian Coalescent (PSMC)21 to reconstruct population size changes ,and the multiple sequentially Markovian coalescent22 (MSMC) to study the time course of population separations. We infer that the population ancestral to all present day humans began to develop substructure at least two hundred thousand years ago (kya), which is most apparent when comparing the ancestors of some present-day African hunter-gatherers (southern African KhoeSan and central African Mbuti Pygmies) and other populations (Fig. 2a). However, it is also clear that this substructure developed slowly, as all pairs of present-day populations including African hunter-gatherer share a substantial subset of their ancestors as recently as a hundred thousand years ago2326. Quoting the time at which MSMC infers that more than 50% (25–75%) of lineages for a pair of populations are descended from the same ancestral population, we estimate that non-Africans separated substantially from KhoeSan 131 (82–173) kya and almost as anciently from the Mbuti around 112 (67–171) kya. Within Africa (Fig. 2a–b), we infer that the Yoruba separated substantially from the KhoeSan 87 (58–120) kya; from the Mbuti 56 (32–85) kya; and from the Dinka 19 (9–25) kya. We estimate a relatively rapid 21 (21–36) kya separation of northern and southern KhoeSan24,27 potentially reflecting isolation since the last glacial maximum; and 38 (27–44) kya separation between western (Biaka) and eastern (Mbuti) Pygmies, confirming very old substructure between these two central African hunter gatherer groups28. Outside Africa, the most ancient structure dates to around 50 kya (Fig. 2c) during or shortly after the deepest part of the shared non-African bottleneck 40–60 kya, consistent with the archaeological evidence of the dispersal of modern humans into Eurasia during this period. We are not confident about the estimates of the date of separation of Australians, New Guineans and Andamanese from other populations because we find that these inferences change depending on the computational method we use for phasing, likely due to these populations not being represented in the 1000 Genomes haploid genome reference panel (Supplementary Information section 9). We caution that the date estimates also do not take into account uncertainty about the true value of the human mutation rate, which could plausibly be 30% higher or lower than the point estimate we use29.

Figure 2
Cross-coalescence rates and effective population sizes for selected population pairs

Early modern human dispersals contributed little to non-African populations

There is intense debate about whether present-day Australians, New Guineans and Asian “Negrito” populations are descended from the same source population as mainland Eurasians, or whether they also derive some ancestry from an early, independent dispersal of modern humans into Asia3032. To explore this scenario rigorously, we fit an admixture graph33—a phylogenetic tree incorporating mixture events—to the allele frequency correlations among Neanderthals, Denisovans, Upper Paleolithic Europeans, East Asians, New Guineans, Australians, and Andamanese. We obtain a good fit to the data if we include known Neanderthal and Denisovan introgression and model all modern human ancestry in New Guineans, Australians and Andamanese as part of an eastern clade together with mainland East Asians (Supplementary Information section 11; Fig. 3). Furthermore, when we manually introduce a deeply diverging modern human lineage contributing ancestry to Australians, New Guineans, and Andamanese (or when we repeat the analysis in a model without Andamanese), no position or proportion of the deep lineage improves the fit. If this putative source population branched off the main lineage leading to non-Africans more than about 10–20 ky prior to the separation of European and East Asian ancestors, we obtain an upper bound of a few percent for the possible contribution to Australians and New Guineans (Fig. 3 inset; Supplementary Information section 11). These results are at odds with an inference of substantial early dispersal ancestry in a previous analysis of an Australian genome32; however, that study used a less complete model that, notably, did not include the known Denisovan admixture into Australo-Melanesians18. The findings for Australians are also unlikely to be due to some unusual feature of the individuals we sequenced, as when we compared three different Australian samples for which there is published genome-wide data, they are all consistent with descending from a common homogeneous population since separation from New Guineans (Supplementary Information section 10). These results are not in conflict with skeletal and archaeological evidence of an early modern human presence outside of Africa30,34, as early migrations could have occurred but not contributed substantially to present-day populations. The possibility of populations that once flourished but did not contribute substantially to living groups is especially plausible now that ancient DNA from the ~45 kya Ust’-Ishim29 and the ~40 kya Oase 1 individuals35 has documented directly their existence.

Figure 3
Present-day populations have negligible ancestry from an early dispersal of modern humans out of Africa

More mutation accumulation in non-Africans than in Africans

The SGDP data provide an opportunity to compare the rates at which mutations have accumulated across populations. We restricted our analyses to samples for which our genotypes are likely to be most reliable (this included restricting to samples which were all processed in the same way), and we used the highest level of filtering (“level 9”) (Supplementary Information section 7). We pooled samples by region to increase power, and for all pairs of regions, computed the expected number of positions where, if we picked a random chromosome from both, region A would mismatch chimpanzee and region B would be identical to chimpanzee (or vice versa). If the rate of accumulation of mutation has been the same since the two populations diverged, these numbers are expected to be equal36. However, when we compute the ratio of mutations on one lineage or the other since separation, we find a subtle (average of 0.5%) but significant excess of mutations in non-Africans relative to sub-Saharan Africans (3.3<|Z|<9.4 standard errors from zero; Extended Data Table 1). Because any difference must reflect events since non-African / African population divergence which is a less than a tenth of average genetic divergence (Fig. 2a), this implies a greater difference in mutation accumulation rates since population divergence (~5%). We were concerned that these results might be biased by the fact that the human genome reference sequence is more closely related to non-Africans than to Africans, or by higher levels of heterozygosity in Africans, as both these issues could make detection of divergent sites in Africans more difficult. However, we replicated the findings after remapping to chimpanzee, which is equally distant to all present populations, and after restricting analyses to the X chromosome in males (males only have a single X chromosome, and so this procedure avoids bias due to different error rates in detecting heterozygous genotypes in populations with different rates of heterozygosity) (Extended Data Fig. 5). These observations are most likely to be explained by acceleration in the rate of mutation accumulation in non-Africans, since the same signal appears in comparisons to sub-Saharan Africans related in different ways to non-Africans (Extended Data Table 1). It is known that the rate of CCT>CTT mutations differs across human populations. However, this particular mutation class was found to be enriched relative to Africans in Europeans but not in East Asians, and thus cannot explain our signal37. One of several possible explanations for these findings is a decrease in the generation interval in non-Africans compared to Africans since separation38.

No evidence for species-wide sweeps since the origin of anatomically modern humans

We finally used the SGDP dataset to address the hypothesis that the widespread appearance of modern human behavior in the archaeological record after ~50 kya was driven by one or a few changes in neurological genes that swept through the population shortly before this time39. We first applied the 3P-CLR method40 to search for locations in the genome with low allele frequency differentiation between KhoeSan and other modern humans, combined with high differentiation between modern and archaic (Neanderthal and Denisovan) humans, as might be expected from a selective sweep in the ancestors of all modern humans (Supplementary Information section 12) (Extended Data Figure 6). We found no strong outlier signals, although a caveat is that our scan has imperfect power and we could not apply it to filtered sections of the genome. We also applied the PSMC method21 to estimate the average time since the most recent common ancestor (TMRCA) of individuals’ two chromosomes in the genomic regions within the largest 3P-CLR peaks (38 peaks corresponding the top 0.1%). In none of the regions did we find that the great majority of all pairs of modern humans are inferred to share a common ancestor <100 kya, as would be expected for a sweep just prior to ~50 kya years ago (Supplementary Data Table 2).

As a second approach to scanning for species-wide selective sweeps, we applied the PSMC to infer TMRCA for SGDP samples across the entire genome. This analysis found no regions where the great majority of pairs of human genomes are inferred to share a common ancestor <100 kya (the largest fraction seen anywhere in the genome is 68%; Extended Data Fig. 7).

Taken together, these results do not rule out the possibility that genetic changes played a meaningful role in the changes in human behaviors after 50 kya; for example, changing selection can produce shifts in the frequencies of pre-existing mutations to bring a population to a new and advantageous set-point for a phenotype as occurred in the case of height differences between northern and southern Europeans41. For polygenic selection, however, genetics is not the creative force, but instead is responding to selection pressures imposed by new environmental conditions or lifestyles. Thus, our results provide evidence against a model in which one or a few mutations were responsible for the rapid developments in human behavior in the last 50 kya. Instead, changes in lifestyles due to cultural innovation or exposure to new environments are likely to have been the ultimate driving forces behind the rapid transformations in human behavior that became evident after 50 ky42,43.

Extended Data

Extended Data Figure 1

An external file that holds a picture, illustration, etc.
Object name is nihms798259f4.jpg
Heatmap of fraction of heterozygous sites missed in the 1000 Genomes Project

For each sample, we examine all heterozygous sites passing filter level 1, and compute the fraction included as known polymorphisms in the 1000 Genomes Project.

Extended Data Figure 2

An external file that holds a picture, illustration, etc.
Object name is nihms798259f5.jpg
Worldwide variation in human short tandem repeats

A: Mean STR length is reported as the average of the length difference (in base pairs) from the GRCh37 reference for each genotype. Bubble area scales with the number of calls compared at each point. B: and C: show the first two principal components after performing principal component analysis on tetranucleotide and homopolymer genotypes, respectively. Colors represent the region of origin of each sample. D: Pairwise FST values between populations computed using only SNPs vs. using combined SNP+STR loci. E: Block jackknife standard errors for the SNP vs. SNP+STR FST analysis. The red dashed lines give the best-fit line, described by the formula in red. The black dashed line denotes the diagonal.

Extended Data Figure 3

An external file that holds a picture, illustration, etc.
Object name is nihms798259f6.jpg
ADMIXTURE analysis

We carried out unsupervised ADMIXTURE 1.239,44 analysis over the 300 SGDP individuals in 20 replicates with randomly chosen initial seeds, varying the number of ancestral populations between K=2 and K=12 and using default 5-fold cross-validation (--cv flag). We used genotypes of at least filter level 1, and restricted analysis to sites where at least two individuals carried the variant allele (as singleton variants are non-informative for population clustering). After further filtering sites with at least 99% completeness and performing linkage-disequilibrium based pruning in PLINK 1.945,46 with parameters (--indep-pairwise 1000 100 0.2), a total of 482,515 single nucleotide polymorphisms remained. This figure shows the highest likelihood replicate for each value of K. We found that log likelihood monotonically increases with K, while the value K=5 minimizes cross-validation error (not shown). The solution at K=5 corresponds to major continental groups (Sub-Saharan Africans, Oceanians, East Asians, Native Americans, and West Eurasians), but we show the full range of K here as they illustrate finer-scale population structure that may be useful to users of the data.

Extended Data Figure 4

An external file that holds a picture, illustration, etc.
Object name is nihms798259f7.jpg
Principal component analysis and neighbor joining tree

A: Principal component analysis. B: Neighbor-joining tree based on FST values for all populations with at least two samples.

Extended Data Figure 5

An external file that holds a picture, illustration, etc.
Object name is nihms798259f8.jpg
Fewer accumulated mutations in Africans than in non-Africans confirmed by mapping to chimpanzee

We compute a statistic D(Population A, Population B, Chimp), measuring the difference in the rate of matching to chimpanzee in Population A compared to Population B. The evidence of mismatching to chimpanzee is seen when we restrict to the male X chromosome to eliminate possible effects due to differences in heterozygosity across populations, and map to the chimpanzee genome which is phylogenetically symmetrically related to all present-day humans. We find that in 78 randomly chosen Population A = African and Population B = non-African pairs of males, transversion substitutions show no consistent skew from zero, but transition substitutions do.

Extended Data Figure 6

An external file that holds a picture, illustration, etc.
Object name is nihms798259f9.jpg
3P-CLR scan for positive selection

The red line denotes the 99.9% quantile cutoff. The genes in the top 5 regions are labeled. A: Scan for selection on the San terminal branch. B: Scan for selection on the non-San terminal branch. C: Scan for selection on the ancestral modern human branch.

Extended Data Figure 7

An external file that holds a picture, illustration, etc.
Object name is nihms798259f10.jpg
Scan for genomic locations where the great majority of present-day humans share a recent common ancestor

We carried out PSMC analysis on 40 pairs of haploid genomes chosen to sample some of the most deeply divergent present-day human lineages. We recorded the time since the most recent common ancestor (TMRCA) at each position, and rescaled to obtain an estimate of absolute time (Supplementary Information section 12). A: Distribution across the genome of the fraction of TMRCAs below specified date cutoffs. For the 100 kya cutoff, the maximum fraction observed anywhere in the genome is 68%. B: Distribution across the genome of the date T at which specified fractions of sample pairs are inferred to have a TMRCA less than T. C: Percentile points of the cumulative distribution function of B.

Extended Data Table 1

Fewer accumulated mutations in Africans than in non-Africans.

All
autosomes
All X
chromosome
Lowest B
quintile
Highest B
quintile
Population APopulation BD×100ZD×100ZD×100ZD×100Z
KhoesanOceania−0.35−8.2−0.70−2.7−0.68−6.4−0.14−1.7
AfricaAmerica−0.33−9.4−0.73−2.8−0.65−7.3−0.18−2.6
KhoesanWestEurasia−0.30−7.5−0.68−3.1−0.63−6.3−0.17−2.1
AfricaOceania−0.29−8.5−0.66−3.2−0.55−6.6−0.07−1.0
AfricaWestEurasia−0.25−8.5−0.66−3.1−0.49−6.4−0.11−1.8
KhoesanSouthAsia−0.24−6.0−0.56−2.7−0.61−6.3−0.11−1.4
AfricaEastAsia−0.20−6.6−0.65−2.5−0.42−5.2−0.10−1.5
AfricaCentralAsiaSiberia−0.20−6.2−0.55−2.2−0.48−6.3−0.05−0.7
PygmyWestEurasia−0.19−4.8−0.46−1.4−0.43−4.6−0.04−0.5
AfricaSouthAsia−0.18−6.4−0.50−2.0−0.46−6.3−0.03−0.5
CentralAsiaSiberiaOceania−0.13−3.9−0.15−0.6−0.09−1.1−0.03−0.4
PygmySouthAsia−0.13−3.3−0.38−1.1−0.38−4.20.020.2
EastAsiaOceania−0.13−4.10.000.0−0.17−2.10.040.6
KhoesanPygmy−0.10−2.6−0.14−0.4−0.16−1.6−0.12−1.5
SouthAsiaWestEurasia−0.08−4.3−0.20−1.2−0.05−1.0−0.10−2.7
CentralAsiaSiberiaWestEurasia−0.06−2.2−0.16−0.8−0.01−0.2−0.09−1.6
EastAsiaWestEurasia−0.06−2.1−0.00−0.0−0.08−1.0−0.02−0.3
CentralAsiaSiberiaEastAsia−0.00−0.2−0.18−1.10.071.2−0.08−1.8
AfricaPygmy−0.00−0.1−0.06−0.20.030.4−0.06−0.8
EastAsiaSouthAsia0.020.70.221.7−0.04−0.70.081.7
CentralAsiaSiberiaSouthAsia0.020.70.050.30.020.4−0.00−0.0
AmericaOceania0.030.90.110.40.101.10.131.7
OceaniaWestEurasia0.082.3−0.03−0.10.101.1−0.04−0.6
AfricaKhoesan0.102.90.170.70.232.60.071.0
AmericaWestEurasia0.113.60.110.40.192.20.081.3
CentralAsiaSiberiaPygmy0.143.40.320.90.434.5−0.04−0.4
OceaniaSouthAsia0.144.80.220.90.131.70.040.7
EastAsiaPygmy0.153.60.491.40.373.90.040.5
AmericaEastAsia0.185.90.090.30.283.60.111.8
AmericaCentralAsiaSiberia0.186.20.341.70.232.90.183.1
AmericaSouthAsia0.186.40.341.50.223.00.183.1
OceaniaPygmy0.245.40.461.30.454.60.020.2
CentralAsiaSiberiaKhoesan0.256.00.572.90.646.30.091.1
EastAsiaKhoesan0.256.20.683.20.595.90.141.7
AmericaPygmy0.265.90.581.60.585.70.091.0
AmericaKhoesan0.378.70.763.30.777.30.222.5

We compute a statistic D(Population A, Population B, Chimp), measuring the difference in the rate of matching to chimpanzee in Population A compared to Population B. For all the autosomes, we observe highly significant signals (3.3<|Z|<9.4) of excess mismatching to chimpanzee in non-Africans compared to Africans, using a standard error from a Block Jackknife. We highlight |D|>0.002 in blue, and |Z|>3 in yellow. The deviations from zero are greatest in subsets of the genome where the time since two populations split comprises a relatively larger fraction of the total genetic divergence time between the populations; this is the direction expected from a mutation accumulation change since divergence. Compared to all the autosomes as a baseline, a least squares fit indicate that the deviations are 2.2-times higher on chromosome X, 2.0 times higher in the quintile of lowest B-statistic (closest to functionally important regions), and 0.43 times as high in the quintile of lowest B-statistic (furthest from functional regions).

Supplementary Material

supp_datatable1

supp_datatable2

supp_info

Acknowledgments

We thank the volunteers who donated samples. We thank Hélène Blanche, Nicole Boivin, Howard Cann (deceased), Evan Eichler, Hank Greely, Michael Petraglia, Kay Prüfer, Alan Rogers, Matthias Steinrücken, Udo Stenzel and Peter Sudmant for comments, critiques, discussions, or advice on assembling samples. The sequencing was funded by the Simons Foundation (SFARI 280376) and the U.S. National Science Foundation (BCS-1032255). I.M. was supported by a Long Term Fellowship grant LT001095/2014 from the Human Frontier Science program. P.S. was supported by the Wenner-Gren foundation and the Swedish Research Council (VR grant 2014-453). T.W and M.G were supported by an NIJ grant 2014-DN-BX-K089. Y.E. was supported by a Career Award at the Scientific Interface from the Burroughs Wellcome Fund and by NIJ grant 2014-DN-BX-K089. T.K. was supported by ERC Starting Investigator grant FP7 - 261213. R.S. received support from Russian Foundation for Basic Research (#15-04-02543). S.D. received support from Russian Foundation for Basic Research (#16-34-00599). R.K., E.K., and S.L. were supported by the Russian Foundation for Basic Research (11-04-00725-a). E.B. was supported by the Russian Scientific Fund (16-14-10354) and by the Russian Foundation for Basic Research (16-06-00303). O.B. was supported by the Russian Foundation for Basic Research (16-06-00364). D.M.B., H.S., E.M., R.V., and M.M. were supported by Institutional Research Funding from the Estonian Research Council IUT24-1 and by the European Regional Development Fund (European Union) through the Centre of Excellence in Genomics to Estonian Biocentre and University of Tartu. D.C. was supported by the Spanish MINECO grant CGL-44351-P. L.B.J. and W.S.W were supported by NIH grant GM59290. S.A.T. was supported by NIH grants 5DP1ES022577 05, 1R01DK104339-01, and 1R01GM113657-01. C.T.S. was supported by The Wellcome Trust grant 098051. C.M.B. was supported by NSF grants 0924726 and 1153911. K.T was supported by CSIR Network Project grant (GENESIS: BSC0121). J.P.S. and Y.S.S. were supported in part by an NIH grant R01-GM094402, and a Packard Fellowship for Science and Engineering. G.R. J.K and S.P. were funded by the Max Planck Society. N.P. and D.R. were supported by NIH grant GM100233 and D.R. is a Howard Hughes Medical Institute investigator.

Footnotes

Author information and data access: Raw data for 277 genomes are available through the EBI European Nucleotide Archive under accession numbers PRJEB9586 and ERP010710. For the remaining 23 genomes, the informed consent documentation is not consistent with fully public data release; data for these genomes can be accessed through a password protected link by researchers who send S.M. and D.R. a signed letter containing the following language: “With regards to the non-public samples from the Simons Genome Diversity Project, I agree that: (a) I will not secondarily distribute the data to anyone, (b) I will not post it publicly, (c) I will make no attempt to connect the genetic data to personal identifiers for the samples, (d) I will not use the data for any commercial purposes.” The identifiers for these 23 genomes are BR_Kashmiri_Pandit-1, BR_Kharia-1, BR_Kurumba-1, BR_Mala-1, BR_Onge-1, BR_Onge-2, S_Igbo-1, S_Igbo-2, S_Kongo-2, S_Lemande-1, S_Lemande-2, S_Chipewyan-1, S_Chipewyan-2, S_Cree-1, S_Cree-2, S_Nahua-1, S_Nahua-2, T_Sherpa-2, T_Tibetan-1, T_Tibetan-2 and T_Sherpa-1, and are designated by code “Y” in the seventh column of Supplementary Data Table 1. Compact versions of the SGDP dataset and software for accessing it are available at (http://genetics.med.harvard.edu/reichlab/Reich_Lab/Datasets.html).

The authors declare competing financial interests. Ugur Hodoglugil is employed by NextBio, a division of Illumina Ltd.

Readers are welcome to comment on the online version of the paper.

Author contributions: S.M., Y.E., Y.S.S., S.P., J.K., N.P. and D.R. supervised the study. S.N., N.R., C.G., G.P., F.B., G.D., I.G.R., A.R.J., P.D., D.M.B., C.M.B., C.C., T.H., A.M.-E., O.L.P., E.B., O.B., S.K.-Y., H.S., D.T., L.Y., C.T.-S., Y.X., M.S.A., A.R.-L., C.B., A.D.R., C.J., E.B.S., E.M., J.P., R.V., B.M.H., U.H., R.W.M., A.S., G.S., J.T.S.W., R.K., E.K., S.L., G.A., D.C., M.H., T.K., W.K., C.W., D.L., M.B., L.B.J., S.A.T., W.S.W., M.M., S.D., R.S., L.S., K.T. and D.R. assembled samples. S.M., H.L., M.L., I.M., M.G., F.R., J.P.S. M.Z., N.C., A.T., P.S., I.L., S.S., Q.F., G.R., Y.S., N.P. and D.R. performed analyses. S.M., H.L., M.L., I.M., M.G., F.R., M.Z., N.P. and D.R. wrote the manuscript with help from all co-authors.

References

1. Genomes Project C, et al. An integrated map of genetic variation from 1,092 human genomes. Nature. 2012;491:56–65. [PMC free article] [PubMed]
3. Li H, Durbin R. Fast and accurate long-read alignment with Burrows-Wheeler transform. Bioinformatics. 2010;26:589–595. [PMC free article] [PubMed]
4. McKenna A, et al. The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome research. 2010;20:1297–1303. [PubMed]
6. Sudmant PH, et al. Global diversity, population stratification, and selection of human copy number variation. Science. 2015 [PMC free article] [PubMed]
7. Gymrek M, Erlich Y. Profiling short tandem repeats from short reads. Methods Mol Biol. 2013;1038:113–135. [PubMed]
8. Gymrek M, Golan D, Rosset S, Erlich Y. lobSTR: A short tandem repeat profiler for personal genomes. Genome research. 2012;22:1154–1162. [PubMed]
9. Alexander DH, Novembre J, Lange K. Fast model-based estimation of ancestry in unrelated individuals. Genome research. 2009;19:1655–1664. [PubMed]
10. Patterson N, Price AL, Reich D. Population structure and eigenanalysis. PLoS genetics. 2006;2:e190. [PubMed]
11. Keinan A, Mullikin JC, Patterson N, Reich D. Accelerated genetic drift on chromosome X during the human dispersal out of Africa. Nature genetics. 2009;41:66–70. [PMC free article] [PubMed]
12. Keinan A, Reich D. Can a sex-biased human demography account for the reduced effective population size of chromosome X in non-Africans? Molecular biology and evolution. 2010;27:2312–2321. [PMC free article] [PubMed]
13. Verdu P, et al. Sociocultural behavior, sex-biased admixture, and effective population sizes in Central African Pygmies and non-Pygmies. Molecular biology and evolution. 2013;30:918–937. [PMC free article] [PubMed]
14. Joiris DV. THE FRAMEWORK OF CENTRAL AFRICAN HUNTER-GATHERERS AND NEIGHBOURING SOCIETIES. African Study Monographs. 2003;(Suppl.28):57–79. x.
15. Green RE, et al. A draft sequence of the Neandertal genome. Science. 2010;328:710–722. [PMC free article] [PubMed]
16. Meyer M, et al. A High-Coverage Genome Sequence from an Archaic Denisovan Individual. Science. 2012 [PMC free article] [PubMed]
17. Wall JD, et al. Higher levels of neanderthal ancestry in East Asians than in Europeans. Genetics. 2013;194:199–209. [PubMed]
18. Reich D, et al. Genetic history of an archaic hominin group from Denisova Cave in Siberia. Nature. 2010;468:1053–1060. [PMC free article] [PubMed]
19. Prufer K, et al. The complete genome sequence of a Neanderthal from the Altai Mountains. Nature. 2014;505:43–49. [PMC free article] [PubMed]
20. Skoglund P, Jakobsson M. Archaic human ancestry in East Asia. Proceedings of the National Academy of Sciences of the United States of America. 2011;108:18301–18306. [PubMed]
21. Li H, Durbin R. Inference of human population history from individual whole-genome sequences. Nature. 2011;475:493–496. [PMC free article] [PubMed]
22. Schiffels S, Durbin R. Inferring human population size and separation history from multiple genome sequences. Nature genetics. 2014;46:919–925. [PMC free article] [PubMed]
23. Gronau I, Hubisz MJ, Gulko B, Danko CG, Siepel A. Bayesian inference of ancient human demography from individual genome sequences. Nature genetics. 2011;43:1031–1034. [PMC free article] [PubMed]
24. Schlebusch CM, et al. Genomic variation in seven Khoe-San groups reveals adaptation and complex African history. Science. 2012;338:374–379. [PubMed]
25. Veeramah KR, et al. An early divergence of KhoeSan ancestors from those of other modern humans is supported by an ABC-based analysis of autosomal resequencing data. Molecular biology and evolution. 2012;29:617–630. [PMC free article] [PubMed]
26. Labuda D, Zietkiewicz E, Yotova V. Archaic lineages in the history of modern humans. Genetics. 2000;156:799–808. [PubMed]
27. Pickrell JK, et al. The genetic prehistory of southern Africa. Nature communications. 2012;3:1143. [PMC free article] [PubMed]
28. Patin E, et al. Inferring the demographic history of African farmers and pygmy hunter-gatherers using a multilocus resequencing data set. PLoS genetics. 2009;5:e1000448. [PMC free article] [PubMed]
29. Fu Q, et al. Genome sequence of a 45,000-year-old modern human from western Siberia. Nature. 2014;514:445–449. [PMC free article] [PubMed]
30. Groucutt HS, et al. Rethinking the dispersal of Homo sapiens out of Africa. Evol Anthropol. 2015;24:149–164. [PubMed]
31. Reyes-Centeno H, Hubbe M, Hanihara T, Stringer C, Harvati K. Testing modern human out-of-Africa dispersal models and implications for modern human origins. Journal of human evolution. 2015 [PubMed]
32. Rasmussen M, et al. An Aboriginal Australian genome reveals separate human dispersals into Asia. Science. 2011;334:94–98. [PMC free article] [PubMed]
33. Patterson NJ, et al. Ancient Admixture in Human History. Genetics. 2012
34. Liu W, et al. The earliest unequivocally modern humans in southern China. Nature. 2015;526:696–699. [PubMed]
35. Fu Q, et al. An early modern human from Romania with a recent Neanderthal ancestor. Nature. 2015 [PMC free article] [PubMed]
36. Do R, et al. No evidence that selection has been less effective at removing deleterious mutations in Europeans than in Africans. Nature genetics. 2015;47:126–131. [PMC free article] [PubMed]
37. Harris K. Evidence for recent, population-specific evolution of the human mutation rate. Proceedings of the National Academy of Sciences of the United States of America. 2015;112:3439–3444. [PubMed]
38. Segurel L, Wyman MJ, Przeworski M. Determinants of mutation rate variation in the human germline. Annu Rev Genomics Hum Genet. 2014;15:47–70. [PubMed]
39. Klein RG, Edgar B. The dawn of human culture. Wiley; 2002.
40. Racimo F. Testing for Ancient Selection Using Cross-population Allele Frequency Differentiation. Genetics. 2015 [PubMed]
41. Turchin MC, et al. Evidence of widespread selection on standing variation in Europe at height-associated SNPs. Nature genetics. 2012;44:1015–1019. [PMC free article] [PubMed]
42. McBrearty S, Brooks AS. The revolution that wasn't: a new interpretation of the origin of modern human behavior. Journal of human evolution. 2000;39:453–563. [PubMed]
43. Renfrew C. Prehistory : the making of the human mind. Modern library paperback edn. Modern Library; 2009.
44. Alexander DH, Lange K. Enhancements to the ADMIXTURE algorithm for individual ancestry estimation. BMC bioinformatics. 2011;12:246. [PMC free article] [PubMed]
45. Chang C, et al. Second-generation PLINK: rising to the challenge of larger and richer datasets. GigaScience. 2015;4:7. [PMC free article] [PubMed]
46. Purcell S, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. American journal of human genetics. 2007;81:559–575. [PubMed]