Rem I. Sukernik

Natural selection shaped regional mtDNA variation in humans

Human mtDNA shows striking regional variation, traditionally attributed to genetic drift. However, it is not easy to account for the fact that only two mtDNA lineages (M and N) left Africa to colonize Eurasia and that lineages A, C, D, and G show a 5-fold enrichment from central Asia to Siberia. As an alternative to drift, natural selection might have enriched for certain mtDNA lineages as people migrated north into colder climates. To test this hypothesis we analyzed 104 complete mtDNA sequences from all global regions and lineages. African mtDNA variation did not significantly deviate from the standard neutral model, but European, Asian, and Siberian plus Native American variations did. Analysis of amino acid substitution mutations (nonsynonymous, Ka) versus neutral mutations (synonymous, Ks) (ka/ks) for all 13 mtDNA protein-coding genes revealed that the ATP6 gene had the highest amino acid sequence variation of any human mtDNA gene, even though ATP6 is one of the more conserved mtDNA proteins. Comparison of the ka/ks ratios for each mtDNA gene from the tropical, temperate, and arctic zones revealed that ATP6 was highly variable in the mtDNAs from the arctic zone, cytochrome b was particularly variable in the temperate zone, and cytochrome oxidase I was notably more variable in the tropics. Moreover, multiple amino acid changes found in ATP6, cytochrome b, and cytochrome oxidase I appeared to be functionally significant. From these analyses we conclude that selection may have played a role in shaping human regional mtDNA variation and that one of the selective influences was climate.

Reconstructing Native American population history.

The peopling of the Americas has been the subject of extensive genetic, archaeological and linguistic research; however, central questions remain unresolved. One contentious issue is whether the settlement occurred by means of a single migration or multiple streams of migration from Siberia. The pattern of dispersals within the Americas is also poorly understood. To address these questions at a higher resolution than was previously possible, we assembled data from 52 Native American and 17 Siberian groups genotyped at 364,470 single nucleotide polymorphisms. Here we show that Native Americans descend from at least three streams of Asian gene flow. Most descend entirely from a single ancestral population that we call ‘First American’. However, speakers of Eskimo–Aleut languages from the Arctic inherit almost half their ancestry from a second stream of Asian gene flow, and the Na-Dene-speaking Chipewyan from Canada inherit roughly one-tenth of their ancestry from a third stream. We show that the initial peopling followed a southward expansion facilitated by the coast, with sequential population splits and little gene flow after divergence, especially in South America. A major exception is in Chibchan speakers on both sides of the Panama isthmus, who have ancestry from both North and South America.

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

Here we report the Simons Genome Diversity Project data set: 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 behavioural 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 of other non-Africans.

Adaptations to Climate-Mediated Selective Pressures in Humans

Humans inhabit a remarkably diverse range of environments, and adaptation through natural selection has likely played a central role in the capacity to survive and thrive in extreme climates. Unlike numerous studies that used only population genetic data to search for evidence of selection, here we scan the human genome for selection signals by identifying the SNPs with the strongest correlations between allele frequencies and climate across 61 worldwide populations. We find a striking enrichment of genic and nonsynonymous SNPs relative to non-genic SNPs among those that are strongly correlated with these climate variables. Among the most extreme signals, several overlap with those from GWAS, including SNPs associated with pigmentation and autoimmune diseases. Further, we find an enrichment of strong signals in gene sets related to UV radiation, infection and immunity, and cancer. Our results imply that adaptations to climate shaped the spatial distribution of variation in humans.

Human adaptations to diet, subsistence, and ecoregion are due to subtle shifts in allele frequency

Human populations use a variety of subsistence strategies to exploit an exceptionally broad range of ecoregions and dietary components. These aspects of human environments have changed dramatically during human evolution, giving rise to new selective pressures. To understand the genetic basis of human adaptations, we combine population genetics data with ecological information to detect variants that increased in frequency in response to new selective pressures. Our approach detects SNPs that show concordant differences in allele frequencies across populations with respect to specific aspects of the environment. Genic and especially nonsynonymous SNPs are overrepresented among those most strongly correlated with environmental variables. This provides genome-wide evidence for selection due to changes in ecoregion, diet, and subsistence. We find particularly strong signals associated with polar ecoregions, with foraging, and with a diet rich in roots and tubers. Interestingly, several of the strongest signals overlap with those implicated in energy metabolism phenotypes from genome-wide association studies, including SNPs influencing glucose levels and susceptibility to type 2 diabetes. Furthermore, several pathways, including those of starch and sucrose metabolism, are enriched for strong signals of adaptations to a diet rich in roots and tubers, whereas signals associated with polar ecoregions are overrepresented in genes associated with energy metabolism pathways.

Mitochondrial DNA Diversity in Indigenous Populations of the Southern Extent of Siberia, and the Origins of Native American Haplogroups

In search of the ancestors of Native American mitochondrial DNA (mtDNA) haplogroups, we analyzed the mtDNA of 531 individuals from nine indigenous populations in Siberia. All mtDNAs were subjected to high‐resolution RFLP analysis, sequencing of the control‐region hypervariable segment I (HVS‐I), and surveyed for additional polymorphic markers in the coding region. Furthermore, the mtDNAs selected according to haplogroup/subhaplogroup status were completely sequenced. Phylogenetic analyses of the resulting data, combined with those from previously published Siberian arctic and sub‐arctic populations, revealed that remnants of the ancient Siberian gene pool are still evident in Siberian populations, suggesting that the founding haplotypes of the Native American A‐D branches originated in different parts of Siberia. Thus, lineage A complete sequences revealed in the Mansi of the Lower Ob and the Ket of the Lower Yenisei belong to A1, suggesting that A1 mtDNAs occasionally found in the remnants of hunting‐gathering populations of northwestern and northern Siberia belonged to a common gene pool of the Siberian progenitors of Paleoindians. Moreover, lineage B1, which is the most closely related to the American B2, occurred in the Tubalar and Tuvan inhabiting the territory between the upper reaches of the Ob River in the west, to the Upper Yenisei region in the east. Finally, the sequence variants of haplogroups C and D, which are most similar to Native American C1 and D1, were detected in the Ulchi of the Lower Amur. Overall, our data suggest that the immediate ancestors of the Siberian/Beringian migrants who gave rise to ancient (pre‐Clovis) Paleoindians have a common origin with aboriginal people of the area now designated the Altai‐Sayan Upland, as well as the Lower Amur/Sea of Okhotsk region.

The Dual Origin and Siberian Affinities of Native American Y Chromosomes

The Y chromosomes of 549 individuals from Siberia and the Americas were analyzed for 12 biallelic markers, which defined 15 haplogroups. The addition of four microsatellite markers increased the number of haplotypes to 111. The major Native American founding lineage, haplogroup M3, accounted for 66% of male Y chromosomes and was defined by the biallelic markers M89, M9, M45, and M3. The founder haplotype also harbored the microsatellite alleles DYS19 (10 repeats), DYS388 (11 repeats), DYS390 (11 repeats), and DYS391 (10 repeats). In Siberia, the M3 haplogroup was confined to the Chukotka peninsula, adjacent to Alaska. The second major group of Native American Y chromosomes, haplogroup M45, accounted for about one-quarter of male lineages. M45 was subdivided by the biallelic marker M173 and by the four microsatellite loci alleles into two major subdivisions: M45a, which is found throughout the Americas, and M45b, which incorporates the M173 variant and is concentrated in North and Central America. In Siberia, M45a haplotypes, including the direct ancestor of haplogroup M3, are concentrated in Middle Siberia, whereas M45b haplotypes are found in the Lower Amur River and Sea of Okhotsk regions of eastern Siberia. Among the remaining 5% of Native American Y chromosomes is haplogroup RPS4Y-T, found in North America. In Siberia, this haplogroup, along with haplogroup M45b, is concentrated in the Lower Amur River/Sea of Okhotsk region. These data suggest that Native American male lineages were derived from two major Siberian migrations. The first migration originated in southern Middle Siberia with the founding haplotype M45a (10-11-11-10). In Beringia, this gave rise to the predominant Native American lineage, M3 (10-11-11-10), which crossed into the New World. A later migration came from the Lower Amur/Sea of Okhkotsk region, bringing haplogroup RPS4Y-T and subhaplogroup M45b, with its associated M173 variant. This migration event contributed to the modern genetic pool of the Na-Dene and Amerinds of North and Central America.

Global diversity, population stratification, and selection of human copy-number variation.

Duplications and deletions in the human genome Duplications and deletions can lead to variation in copy number for genes and genomic loci among humans. Such variants can reveal evolutionary patterns and have implications for human health. Sudmant et al. examined copy-number variation across 236 individual genomes from 125 human populations. Deletions were under more selection, whereas duplications showed more population-specific structure. Interestingly, Oceanic populations retain large duplications postulated to have originated in an ancient Denisovan lineage. Science, this issue 10.1126/science.aab3761 Copy-number variation reveals how selection affects the human genome across the globe. INTRODUCTION Most studies of human genetic variation have focused on single-nucleotide variants (SNVs). However, copy-number variants (CNVs) affect more base pairs of DNA among humans, and yet our understanding of CNV diversity among human populations is limited. RATIONALE We aimed to understand the pattern, selection, and diversity of copy-number variation by analyzing deeply sequenced genomes representing the diversity of all humans. We compared the selective constraints of deletions versus duplications to understand population stratification in the context of the ancestral human genome and to assess differences in CNV load between African and non-African populations. RESULTS We sequenced 236 individual genomes from 125 distinct human populations and identified 14,467 autosomal CNVs and 545 X-linked CNVs with a sequence read-depth approach. Deletions exhibit stronger selective pressure and are better phylogenetic markers of population relationships than duplication polymorphisms. We identified 1036 population-stratified copy-number–variable regions, 295 of which intersect coding regions and 199 of which exhibit extreme signatures of differentiation. Duplicated loci were 1.8-fold more likely to be stratified than deletions but were poorly correlated with flanking genetic diversity. Among these, we highlight a duplication polymorphism restricted to modern Oceanic populations yet also present in the genome of the archaic Denisova hominin. This 225–kilo–base pair (kbp) duplication includes two microRNA genes and is almost fixed among human Papuan-Bougainville genomes. The data allowed us to reconstruct the ancestral human genome and create a more accurate evolutionary framework for the gain and loss of sequences during human evolution. We identified 571 loci that segregate in the human population and another 2026 loci of fixed-copy 2 in all human genomes but absent from the reference genome. The total deletion and duplication load between African and non-African population groups showed no difference after we account for ancestral sequences missing from the human reference. However, we did observe that the relative number of base pairs affected by CNVs compared to single-nucleotide polymorphisms is higher among non-Africans than Africans. CONCLUSION Deletions, duplications, and CNVs have shaped, to different extents, the genetic diversity of human populations by the combined forces of mutation, selection, and demography. Figure Global human CNV diversity and archaic introgression of a chromosome 16 duplication. (Left) The geographic coordinates of populations sampled are indicated on a world map (colored dots). The pie charts show the continental population allele frequency of a single ~225-kbp duplication polymorphism found exclusively among Oceanic populations and an archaic Denisova. (Right) The ancestral structure of this duplication locus (1) and the Denisova duplication structure (2) are shown in relation to their position on chromosome 16. We estimate that the duplication emerged ~440 thousand years ago (ka) in the Denisova and then introgressed into ancestral Papuan populations ~40 ka. In order to explore the diversity and selective signatures of duplication and deletion human copy-number variants (CNVs), we sequenced 236 individuals from 125 distinct human populations. We observed that duplications exhibit fundamentally different population genetic and selective signatures than deletions and are more likely to be stratified between human populations. Through reconstruction of the ancestral human genome, we identify megabases of DNA lost in different human lineages and pinpoint large duplications that introgressed from the extinct Denisova lineage now found at high frequency exclusively in Oceanic populations. We find that the proportion of CNV base pairs to single-nucleotide–variant base pairs is greater among non-Africans than it is among African populations, but we conclude that this difference is likely due to unique aspects of non-African population history as opposed to differences in CNV load.

Mitochondrial Genome Diversity in Arctic Siberians, with Particular Reference to the Evolutionary History of Beringia and Pleistocenic Peopling of the Americas

Through extended survey of mitochondrial DNA (mtDNA) diversity in the Nganasan, Yukaghir, Chuvantsi, Chukchi, Siberian Eskimos, and Commander Aleuts, we filled important gaps in previously unidentified internal sequence variation within haplogroups A, C, and D, three of five (A-D and X) canonical mtDNA lineages that defined Pleistocenic extension from the Old to the New World. Overall, 515 mtDNA samples were analyzed via high-resolution SNP analysis and then complete sequencing of the 84 mtDNAs. A comparison of the data thus obtained with published complete sequences has resulted in the most parsimonious phylogenetic structure of mtDNA evolution in Siberia-Beringia. Our data suggest that although the latest inhabitants of Beringia are well genetically reflected in the Chukchi-, Eskimo-Aleut-, and Na-Dene-speaking Indians, the direct ancestors of the Paleosiberian-speaking Yukaghir are primarily drawn from the southern belt of Siberia when environmental conditions changed, permitting recolonization the high arctic since early Postglacial. This study further confirms that (1) Alaska seems to be the ancestral homeland of haplogroup A2 originating in situ approximately 16.0 thousand years ago (kya), (2) an additional founding lineage for Native American D, termed here D10, arose approximately 17.0 kya in what is now the Russian Far East and eventually spread northward along the North Pacific Rim. The maintenance of two refugial sources, in the Altai-Sayan and mid-lower Amur, during the last glacial maximum appears to be at odds with the interpretation of limited founding mtDNA lineages populating the Americas as a single migration.

Analysis of Mitochondrial DNA Diversity in the Aleuts of the Commander Islands and Its Implications for the Genetic History of Beringia

The Aleuts are aboriginal inhabitants of the Aleutian archipelago, including Bering and Copper (Medny) Islands of the Commanders, and seem to be the survivors of the inhabitants of the southern belt of the Bering Land Bridge that connected Chukotka/Kamchatka and Alaska during the end of the Ice Age. Thirty mtDNA samples collected in the Commanders, as well as seven mtDNA samples from Sireniki Eskimos in Chukotka who belong to the Beringian-specific subhaplogroup D2, were studied through complete sequencing. This analysis has provided evidence that all 37 of these mtDNAs are closely related, since they share the founding haplotype for subhaplogroup D2. We also demonstrated that, unlike the Eskimos and Na-Dene, the Aleuts of the Commanders were founded by a single lineage of haplogroup D2, which had acquired the novel transversion mutation 8910A. The phylogeny of haplogroup D complete sequences showed that (1) the D2 root sequence type originated among the latest inhabitants of Beringia and (2) the Aleut 8910A sublineage of D2 is a part of larger radiation of rooted D2, which gave rise to D2a (Na-Dene), D2b (Aleut), and D2c (Eskimo) sublineages. The geographic specificity and remarkable intrinsic diversity of D2 lineages support the refugial hypothesis, which assumes that the founding population of Eskimo-Aleut originated in Beringan/southwestern Alaskan refugia during the early postglacial period, rather than having reached the shores of Alaska as the result of recent wave of migration from interior Siberia.