Mattias Jakobsson

Genotype, haplotype and copy-number variation in worldwide human populations

Genome-wide patterns of variation across individuals provide a powerful source of data for uncovering the history of migration, range expansion, and adaptation of the human species. However, high-resolution surveys of variation in genotype, haplotype and copy number have generally focused on a small number of population groups. Here we report the analysis of high-quality genotypes at 525,910 single-nucleotide polymorphisms (SNPs) and 396 copy-number-variable loci in a worldwide sample of 29 populations. Analysis of SNP genotypes yields strongly supported fine-scale inferences about population structure. Increasing linkage disequilibrium is observed with increasing geographic distance from Africa, as expected under a serial founder effect for the out-of-Africa spread of human populations. New approaches for haplotype analysis produce inferences about population structure that complement results based on unphased SNPs. Despite a difference from SNPs in the frequency spectrum of the copy-number variants (CNVs) detected—including a comparatively large number of CNVs in previously unexamined populations from Oceania and the Americas—the global distribution of CNVs largely accords with population structure analyses for SNP data sets of similar size. Our results produce new inferences about inter-population variation, support the utility of CNVs in human population-genetic research, and serve as a genomic resource for human-genetic studies in diverse worldwide populations.

Clumpak: a program for identifying clustering modes and packaging population structure inferences across K

The identification of the genetic structure of populations from multilocus genotype data has become a central component of modern population‐genetic data analysis. Application of model‐based clustering programs often entails a number of steps, in which the user considers different modelling assumptions, compares results across different predetermined values of the number of assumed clusters (a parameter typically denoted K), examines multiple independent runs for each fixed value of K, and distinguishes among runs belonging to substantially distinct clustering solutions. Here, we present Clumpak (Cluster Markov Packager Across K), a method that automates the postprocessing of results of model‐based population structure analyses. For analysing multiple independent runs at a single K value, Clumpak identifies sets of highly similar runs, separating distinct groups of runs that represent distinct modes in the space of possible solutions. This procedure, which generates a consensus solution for each distinct mode, is performed by the use of a Markov clustering algorithm that relies on a similarity matrix between replicate runs, as computed by the software Clumpp. Next, Clumpak identifies an optimal alignment of inferred clusters across different values of K, extending a similar approach implemented for a fixed K in Clumpp and simplifying the comparison of clustering results across different K values. Clumpak incorporates additional features, such as implementations of methods for choosing K and comparing solutions obtained by different programs, models, or data subsets. Clumpak, available at http://clumpak.tau.ac.il, simplifies the use of model‐based analyses of population structure in population genetics and molecular ecology.

A worldwide survey of haplotype variation and linkage disequilibrium in the human genome

Recent genomic surveys have produced high-resolution haplotype information, but only in a small number of human populations. We report haplotype structure across 12 Mb of DNA sequence in 927 individuals representing 52 populations. The geographic distribution of haplotypes reflects human history, with a loss of haplotype diversity as distance increases from Africa. Although the extent of linkage disequilibrium (LD) varies markedly across populations, considerable sharing of haplotype structure exists, and inferred recombination hotspot locations generally match across groups. The four samples in the International HapMap Project contain the majority of common haplotypes found in most populations: averaging across populations, 83% of common 20-kb haplotypes in a population are also common in the most similar HapMap sample. Consequently, although the portability of tag SNPs based on the HapMap is reduced in low-LD Africans, the HapMap will be helpful for the design of genome-wide association mapping studies in nearly all human populations.

ADZE: a rarefaction approach for counting alleles private to combinations of populations

Motivation: Analysis of the distribution of alleles across populations is a useful tool for examining population diversity and relationships. However, sample sizes often differ across populations, sometimes making it difficult to assess allelic distributions across groups. Results: We introduce a generalized rarefaction approach for counting alleles private to combinations of populations. Our method evaluates the number of alleles found in each of a set of populations but absent in all remaining populations, considering equal-sized subsamples from each population. Applying this method to a worldwide human microsatellite dataset, we observe a high number of alleles private to the combination of African and Oceanian populations. This result supports the possibility of a migration out of Africa into Oceania separate from the migrations responsible for the majority of the ancestry of the modern populations of Asia, and it highlights the utility of our approach to sample size correction in evaluating hypotheses about population history. Availability: We have implemented our method in the computer pro-gram ADZE, which is available for download at http://rosenberglab.bioinformatics.med.umich.edu/adze.html. Contact: szpiechz@umich.edu

Origins and genetic legacy of Neolithic farmers and hunter-gatherers in Europe.

Farmer Displaced European Hunters Our understanding of prehistoric demography and human evolution has been improved by analysis of ancient DNA. Skoglund et al. (p. 466) describe the retrieval and analysis of genomic DNA from ancient (~5000-year-old) northern European Neolithic individuals within modern Sweden. These include three hunter-gatherers from the Pitted Ware Culture horizon and one farmer ascribed to the Mid-Neolithic North-Central TRB culture. The hunter-gatherers displayed a distinct genetic signature, similar to that of extant northern Europeans, whereas the farmers genetic signature more closely resembled southern Europeans, suggesting migration and admixture during the spread of farming. Genomic analysis of ancient Scandinavians reveals that agricultural expansion was driven by long-range population movement. The farming way of life originated in the Near East some 11,000 years ago and had reached most of the European continent 5000 years later. However, the impact of the agricultural revolution on demography and patterns of genomic variation in Europe remains unknown. We obtained 249 million base pairs of genomic DNA from ~5000-year-old remains of three hunter-gatherers and one farmer excavated in Scandinavia and find that the farmer is genetically most similar to extant southern Europeans, contrasting sharply to the hunter-gatherers, whose distinct genetic signature is most similar to that of extant northern Europeans. Our results suggest that migration from southern Europe catalyzed the spread of agriculture and that admixture in the wake of this expansion eventually shaped the genomic landscape of modern-day Europe.

Upper Palaeolithic Siberian genome reveals dual ancestry of Native Americans

The origins of the First Americans remain contentious. Although Native Americans seem to be genetically most closely related to east Asians, there is no consensus with regard to which specific Old World populations they are closest to. Here we sequence the draft genome of an approximately 24,000-year-old individual (MA-1), from Mal’ta in south-central Siberia, to an average depth of 1×. To our knowledge this is the oldest anatomically modern human genome reported to date. The MA-1 mitochondrial genome belongs to haplogroup U, which has also been found at high frequency among Upper Palaeolithic and Mesolithic European hunter-gatherers, and the Y chromosome of MA-1 is basal to modern-day western Eurasians and near the root of most Native American lineages. Similarly, we find autosomal evidence that MA-1 is basal to modern-day western Eurasians and genetically closely related to modern-day Native Americans, with no close affinity to east Asians. This suggests that populations related to contemporary western Eurasians had a more north-easterly distribution 24,000 years ago than commonly thought. Furthermore, we estimate that 14 to 38% of Native American ancestry may originate through gene flow from this ancient population. This is likely to have occurred after the divergence of Native American ancestors from east Asian ancestors, but before the diversification of Native American populations in the New World. Gene flow from the MA-1 lineage into Native American ancestors could explain why several crania from the First Americans have been reported as bearing morphological characteristics that do not resemble those of east Asians. Sequencing of another south-central Siberian, Afontova Gora-2 dating to approximately 17,000 years ago, revealed similar autosomal genetic signatures as MA-1, suggesting that the region was continuously occupied by humans throughout the Last Glacial Maximum. Our findings reveal that western Eurasian genetic signatures in modern-day Native Americans derive not only from post-Columbian admixture, as commonly thought, but also from a mixed ancestry of the First Americans.

The genome of a Late Pleistocene human from a Clovis burial site in western Montana

Clovis, with its distinctive biface, blade and osseous technologies, is the oldest widespread archaeological complex defined in North America, dating from 11,100 to 10,700 14C years before present (bp) (13,000 to 12,600 calendar years bp). Nearly 50 years of archaeological research point to the Clovis complex as having developed south of the North American ice sheets from an ancestral technology. However, both the origins and the genetic legacy of the people who manufactured Clovis tools remain under debate. It is generally believed that these people ultimately derived from Asia and were directly related to contemporary Native Americans. An alternative, Solutrean, hypothesis posits that the Clovis predecessors emigrated from southwestern Europe during the Last Glacial Maximum. Here we report the genome sequence of a male infant (Anzick-1) recovered from the Anzick burial site in western Montana. The human bones date to 10,705 ± 35 14C years bp (approximately 12,707–12,556 calendar years bp) and were directly associated with Clovis tools. We sequenced the genome to an average depth of 14.4× and show that the gene flow from the Siberian Upper Palaeolithic Mal’ta population into Native American ancestors is also shared by the Anzick-1 individual and thus happened before 12,600 years bp. We also show that the Anzick-1 individual is more closely related to all indigenous American populations than to any other group. Our data are compatible with the hypothesis that Anzick-1 belonged to a population directly ancestral to many contemporary Native Americans. Finally, we find evidence of a deep divergence in Native American populations that predates the Anzick-1 individual.

Genomic Variation in Seven Khoe-San Groups Reveals Adaptation and Complex African History

African Origins Humans originated in Africa and then spread across the globe. The high genetic diversity found in sub-Saharan Africa is consistent with this view, but the relationships among and within African populations have been less well explored. Schlebusch et al. (p. 374, published online 20 September) genotyped 220 individuals from 11 populations representing groups from Southern Africa to determine their relationships and history. The data suggest that modern-day human populations arose from a complex blend of mixing between groups and genetic stratification. Cutting-edge genomic approaches test hypotheses about the roots of human history in southern African indigenous populations. The history of click-speaking Khoe-San, and African populations in general, remains poorly understood. We genotyped ∼2.3 million single-nucleotide polymorphisms in 220 southern Africans and found that the Khoe-San diverged from other populations ≥100,000 years ago, but population structure within the Khoe-San dated back to about 35,000 years ago. Genetic variation in various sub-Saharan populations did not localize the origin of modern humans to a single geographic region within Africa; instead, it indicated a history of admixture and stratification. We found evidence of adaptation targeting muscle function and immune response; potential adaptive introgression of protection from ultraviolet light; and selection predating modern human diversification, involving skeletal and neurological development. These new findings illustrate the importance of African genomic diversity in understanding human evolutionary history.

Demographic history of european populations of Arabidopsis thaliana.

The model plant species Arabidopsis thaliana is successful at colonizing land that has recently undergone human-mediated disturbance. To investigate the prehistoric spread of A. thaliana, we applied approximate Bayesian computation and explicit spatial modeling to 76 European accessions sequenced at 876 nuclear loci. We find evidence that a major migration wave occurred from east to west, affecting most of the sampled individuals. The longitudinal gradient appears to result from the plant having spread in Europe from the east ∼10,000 years ago, with a rate of westward spread of ∼0.9 km/year. This wave-of-advance model is consistent with a natural colonization from an eastern glacial refugium that overwhelmed ancient western lineages. However, the speed and time frame of the model also suggest that the migration of A. thaliana into Europe may have accompanied the spread of agriculture during the Neolithic transition.

Genomic evidence for the Pleistocene and recent population history of Native Americans

Genetic history of Native Americans Several theories have been put forth as to the origin and timing of when Native American ancestors entered the Americas. To clarify this controversy, Raghavan et al. examined the genomic variation among ancient and modern individuals from Asia and the Americas. There is no evidence for multiple waves of entry or recurrent gene flow with Asians in northern populations. The earliest migrations occurred no earlier than 23,000 years ago from Siberian ancestors. Amerindians and Athabascans originated from a single population, splitting approximately 13,000 years ago. Science, this issue 10.1126/science.aab3884 Genetic variation within ancient and extant Native American populations informs on their migration into the Americas. INTRODUCTION The consensus view on the peopling of the Americas is that ancestors of modern Native Americans entered the Americas from Siberia via the Bering Land Bridge and that this occurred at least ~14.6 thousand years ago (ka). However, the number and timing of migrations into the Americas remain controversial, with conflicting interpretations based on anatomical and genetic evidence. RATIONALE In this study, we address four major unresolved issues regarding the Pleistocene and recent population history of Native Americans: (i) the timing of their divergence from their ancestral group, (ii) the number of migrations into the Americas, (iii) whether there was ~15,000 years of isolation of ancestral Native Americans in Beringia (Beringian Incubation Model), and (iv) whether there was post-Pleistocene survival of relict populations in the Americas related to Australo-Melanesians, as suggested by apparent differences in cranial morphologies between some early (“Paleoamerican”) remains and those of more recent Native Americans. We generated 31 high-coverage modern genomes from the Americas, Siberia, and Oceania; 23 ancient genomic sequences from the Americas dating between ~0.2 and 6 ka; and SNP chip genotype data from 79 present-day individuals belonging to 28 populations from the Americas and Siberia. The above data sets were analyzed together with published modern and ancient genomic data from worldwide populations, after masking some present-day Native Americans for recent European admixture. RESULTS Using three different methods, we determined the divergence time for all Native Americans (Athabascans and Amerindians) from their Siberian ancestors to be ~20 ka, and no earlier than ~23 ka. Furthermore, we dated the divergence between Athabascans (northern Native American branch, together with northern North American Amerindians) and southern North Americans and South and Central Americans (southern Native American branch) to be ~13 ka. Similar divergence times from East Asian populations and a divergence time between the two branches that is close in age to the earliest well-established archaeological sites in the Americas suggest that the split between the branches occurred within the Americas. We additionally found that several sequenced Holocene individuals from the Americas are related to present-day populations from the same geographical regions, implying genetic continuity of ancient and modern populations in some parts of the Americas over at least the past 8500 years. Moreover, our results suggest that there has been gene flow between some Native Americans from both North and South America and groups related to East Asians and Australo-Melanesians, the latter possibly through an East Asian route that might have included ancestors of modern Aleutian Islanders. Last, using both genomic and morphometric analyses, we found that historical Native American groups such as the Pericúes and Fuego-Patagonians were not “relicts” of Paleoamericans, and hence, our results do not support an early migration of populations directly related to Australo-Melanesians into the Americas. CONCLUSION Our results provide an upper bound of ~23 ka on the initial divergence of ancestral Native Americans from their East Asian ancestors, followed by a short isolation period of no more than ~8000 years, and subsequent entrance and spread across the Americas. The data presented are consistent with a single-migration model for all Native Americans, with later gene flow from sources related to East Asians and, indirectly, Australo-Melanesians. The single wave diversified ~13 ka, likely within the Americas, giving rise to the northern and southern branches of present-day Native Americans. Population history of present-day Native Americans. The ancestors of all Native Americans entered the Americas as a single migration wave from Siberia (purple) no earlier than ~23 ka, separate from the Inuit (green), and diversified into “northern” and “southern” Native American branches ~13 ka. There is evidence of post-divergence gene flow between some Native Americans and groups related to East Asians/Inuit and Australo-Melanesians (yellow). How and when the Americas were populated remains contentious. Using ancient and modern genome-wide data, we found that the ancestors of all present-day Native Americans, including Athabascans and Amerindians, entered the Americas as a single migration wave from Siberia no earlier than 23 thousand years ago (ka) and after no more than an 8000-year isolation period in Beringia. After their arrival to the Americas, ancestral Native Americans diversified into two basal genetic branches around 13 ka, one that is now dispersed across North and South America and the other restricted to North America. Subsequent gene flow resulted in some Native Americans sharing ancestry with present-day East Asians (including Siberians) and, more distantly, Australo-Melanesians. Putative “Paleoamerican” relict populations, including the historical Mexican Pericúes and South American Fuego-Patagonians, are not directly related to modern Australo-Melanesians as suggested by the Paleoamerican Model.