The genetic prehistory of the New World Arctic

Neandertal ancestry through time: Insights from genomes of ancient and present-day humans

Gene flow from Neandertals has shaped the landscape of genetic and phenotypic variation in modern humans. We identify the location and size of introgressed Neandertal ancestry segments in more than 300 genomes spanning the last 50,000 years. We study how Neandertal ancestry is shared among individuals to infer the time and duration of the Neandertal gene flow. We find the correlation of Neandertal segment locations across individuals and their divergence to sequenced Neandertals, both support a model of single major Neandertal gene flow. Our catalog of introgressed segments through time confirms that most natural selection-positive and negative-on Neandertal ancestry variants occurred immediately after the gene flow, and provides new insights into how the contact with Neandertals shaped human origins and adaptation.

The Allen Ancient DNA Resource (AADR) a curated compendium of ancient human genomes

More than two hundred papers have reported genome-wide data from ancient humans. While the raw data for the vast majority are fully publicly available testifying to the commitment of the paleogenomics community to open data, formats for both raw data and meta-data differ. There is thus a need for uniform curation and a centralized, version-controlled compendium that researchers can download, analyze, and reference. Since 2019, we have been maintaining the Allen Ancient DNA Resource (AADR), which aims to provide an up-to-date, curated version of the world's published ancient human DNA data, represented at more than a million single nucleotide polymorphisms (SNPs) at which almost all ancient individuals have been assayed. The AADR has gone through six public releases at the time of writing and review of this manuscript, and crossed the threshold of >10,000 individuals with published genome-wide ancient DNA data at the end of 2022. This note is intended as a citable descriptor of the AADR.

The genetic history of the Southern Andes from present-day Mapuche ancestry

The southernmost regions of South America harbor some of the earliest evidence of human presence in the Americas. However, connections with the rest of the continent and the contextualization of present-day indigenous ancestries remain poorly resolved. In this study, we analyze the genetic ancestry of one of the largest indigenous groups in South America: the Mapuche. We generate genome-wide data from 64 participants from three Mapuche populations in Southern Chile: Pehuenche, Lafkenche, and Huilliche. Broadly, we describe three main ancestry blocks with a common origin, which characterize the Southern Cone, the Central Andes, and Amazonia. Within the Southern Cone, ancestors of the Mapuche lineages differentiated from those of the Far South during the Middle Holocene and did not experience further migration waves from the north. We find that the deep genetic split between the Central and Southern Andes is followed by instances of gene flow, which may have accompanied the southward spread of cultural traits from the Central Andes, including crops and loanwords from Quechua into Mapudungun (the language of the Mapuche). Finally, we report close genetic relatedness between the three populations analyzed, with the Huilliche characterized additionally by intense recent exchanges with the Far South. Our findings add new perspectives on the genetic (pre)history of South America, from the first settlement through to the present-day indigenous presence. Follow-up fieldwork took these results back to the indigenous communities to contextualize the genetic narrative alongside indigenous knowledge and perspectives.

A paleogenome from a Holocene individual supports genetic continuity in Southeast Alaska

Many specifics of the population histories of the Indigenous peoples of North America remain contentious owing to a dearth of physical evidence. Only few ancient human genomes have been recovered from the Pacific Northwest Coast, a region increasingly supported as a coastal migration route for the initial peopling of the Americas. Here, we report paleogenomic data from the remains of a ∼3,000-year-old female individual from Southeast Alaska, named Tatóok yík yées sháawat (TYYS). Our results demonstrate at least 3,000 years of matrilineal genetic continuity in Southeast Alaska, and that TYYS is most closely related to ancient and present-day northern Pacific Northwest Coast Indigenous Americans. We find no evidence of Paleo-Inuit (represented by Saqqaq) ancestry in present-day or ancient Pacific Northwest peoples. Instead, our analyses suggest the Saqqaq genome harbors Northern Native American ancestry. This study sheds further light on the human population history of the northern Pacific Northwest Coast.

The Impact of Modern Admixture on Archaic Human Ancestry in Human Populations

Admixture, the genetic merging of parental populations resulting in mixed ancestry, has occurred frequently throughout the course of human history. Numerous admixture events have occurred between human populations across the world, which have shaped genetic ancestry in modern humans. For example, populations in the Americas are often mosaics of different ancestries due to recent admixture events as part of European colonization. Admixed individuals also often have introgressed DNA from Neanderthals and Denisovans that may have come from multiple ancestral populations, which may affect how archaic ancestry is distributed across an admixed genome. In this study, we analyzed admixed populations from the Americas to assess whether the proportion and location of admixed segments due to recent admixture impact an individual's archaic ancestry. We identified a positive correlation between non-African ancestry and archaic alleles, as well as a slight increase of Denisovan alleles in Indigenous American segments relative to European segments in admixed genomes. We also identify several genes as candidates for adaptive introgression, based on archaic alleles present at high frequency in admixed American populations but low frequency in East Asian populations. These results provide insights into how recent admixture events between modern humans redistributed archaic ancestry in admixed genomes.

Evolution of the Okvik/Old Bering Sea culture of the Bering Strait as a major transition

Great transitions are thought to embody major shifts in locus of selection, labour diversification and communication systems. Such expectations are relevant for biological and cultural systems as decades of research has demonstrated similar dynamics within the evolution of culture. The evolution of the Neo-Inuit cultural tradition in the Bering Strait provides an ideal context for examination of cultural transitions. The Okvik/Old Bering Sea (Okvik/OBS) culture of Bering Strait is the first representative of the Neo-Inuit tradition. Archaeological evidence drawn for settlement and subsistence data, technological traditions and mortuary contexts suggests that Okvik/OBS fits the definition of a major transition given change in the nature of group membership (from families to political groups with social ranking), task organization (emergent labour specialization) and communication (advent of complex art forms conveying social and ideological information). This permits us to develop a number of implications about the evolutionary process recognizing that transitions may occur on three scales: (1) ephemeral variants, as for example, simple technological entities; (2) integrated systems, spanning modular technology to socio-economic strategies; and (3) simultaneous change across all scales with emergent properties. This article is part of the theme issue 'Human socio-cultural evolution in light of evolutionary transitions'.

Bayesian inference of admixture graphs on Native American and Arctic populations

Admixture graphs are mathematical structures that describe the ancestry of populations in terms of divergence and merging (admixing) of ancestral populations as a graph. An admixture graph consists of a graph topology, branch lengths, and admixture proportions. The branch lengths and admixture proportions can be estimated using numerous numerical optimization methods, but inferring the topology involves a combinatorial search for which no polynomial algorithm is known. In this paper, we present a reversible jump MCMC algorithm for sampling high-probability admixture graphs and show that this approach works well both as a heuristic search for a single best-fitting graph and for summarizing shared features extracted from posterior samples of graphs. We apply the method to 11 Native American and Siberian populations and exploit the shared structure of high-probability graphs to characterize the relationship between Saqqaq, Inuit, Koryaks, and Athabascans. Our analyses show that the Saqqaq is not a good proxy for the previously identified gene flow from Arctic people into the Na-Dene speaking Athabascans.

Bidirectional dispersals during the peopling of the North American Arctic

It is thought that Paleo-Inuit were the first people that settled the American Arctic about 5000 BP (before the present) from a migration that crossed Beringia from Northeast Asia. It is theorized that this group initially migrated to the North Slopes of Alaska and subsequently expanded eastward, eventually reaching Greenland. A second circumpolar dispersal of Neo-Inuit from the North Slopes associated with the Thule-Inuk culture has been postulated to have extended eastward around 800 BP, totally replacing the original Paleo-Inuit without admixing. Although generally accepted, this migration scenario is incompatible with previously reported indications of east to west gene flow across the American Arctic. Here we report on the Y-chromosome haplogroup and Y-STR diversity of the four circumpolar populations of the Tuva Republic (N = 24), Northeast Siberia (N = 9), Bethel, Alaska (N = 40), and Barrow, Alaska (N = 31). Four haplogroup lineages (Q-NWT01, Q-M3, Q-M346, and Q-M120) were detected, Q-NWT01 and Q-M3 being the most abundant at 11.11 and 66.67% in Northeast Siberia, 32.50 and 65.00% in Bethel, and 67.74 and 32.26% in Barrow, respectively. The same samples genotyped for Y-chromosome SNPs were typed for 17 Y-STYR loci using the AmpFlSTR Yfiler system. Age estimates and diversity values for the Q-NWT01 and Q-M3 mutations suggest extensive movement of male individuals along the entire longitudinal stretch of the American circumpolar region. Throughout the entire region, Q-M3 exhibits a west to east decreasing gradient in age and diversity while Q-NWT01 indicates the opposite with older TMRCA and higher diversity values running from east to west with the most recent estimates in Canada and Alaska. The high age and diversity values in Greenland are congruent with an origin of the Q-NWT01 mutation in the east of the circumpolar range about 2000-3000 ya. This scenario is incompatible with a complete biological replacement starting about 700 BP of Paleo-Inuit like the Dorset by the Thule-Inuit (Neo-Inuit), as is currently thought, and more parsimonious with gene flow carrying the NWT01 mutation from a pre-Thule population to the ancestors of the present-day Inuit.

The impact of modern admixture on archaic human ancestry in human populations

Admixture, the genetic merging of parental populations resulting in mixed ancestry, has occurred frequently throughout the course of human history. Numerous admixture events have occurred between human populations across the world, as well as introgression between humans and archaic humans, Neanderthals and Denisovans. One example are genomes from populations in the Americas, as these are often mosaics of different ancestries due to recent admixture events as part of European colonization. In this study, we analyzed admixed populations from the Americas to assess whether the proportion and location of admixed segments due to recent admixture impact an individual’s archaic ancestry. We identified a positive correlation between non-African ancestry and archaic alleles, as well as a slight enrichment of Denisovan alleles in Indigenous American segments relative to European segments in admixed genomes. We also identify several genes as candidates for adaptive introgression, based on archaic alleles present at high frequency in admixed American populations but low frequency in East Asian populations. These results provide insights into how recent admixture events between modern humans redistributed archaic ancestry in admixed genomes.

A novel splice-affecting HNF1A variant with large population impact on diabetes in Greenland

Background

The genetic disease architecture of Inuit includes a large number of common high-impact variants. Identification of such variants contributes to our understanding of the genetic aetiology of diseases and improves global equity in genomic personalised medicine. We aimed to identify and characterise novel variants in genes associated with Maturity Onset Diabetes of the Young (MODY) in the Greenlandic population.

Methods

Using combined data from Greenlandic population cohorts of 4497 individuals, including 448 whole genome sequenced individuals, we screened 14 known MODY genes for previously identified and novel variants. We functionally characterised an identified novel variant and assessed its association with diabetes prevalence and cardiometabolic traits and population impact.

Findings

We identified a novel variant in the known MODY gene HNF1A with an allele frequency of 1.9% in the Greenlandic Inuit and absent elsewhere. Functional assays indicate that it prevents normal splicing of the gene. The variant caused lower 30-min insulin (β = -232 pmol/L, β_SD = -0.695, P = 4.43 × 10^-4) and higher 30-min glucose (β = 1.20 mmol/L, β_SD = 0.441, P = 0.0271) during an oral glucose tolerance test. Furthermore, the variant was associated with type 2 diabetes (OR 4.35, P = 7.24 × 10^-6) and HbA1c (β = 0.113 HbA1c%, β_SD = 0.205, P = 7.84 × 10^-3). The variant explained 2.5% of diabetes variance in Greenland.

Interpretation

The reported variant has the largest population impact of any previously reported variant within a MODY gene. Together with the recessive TBC1D4 variant, we show that close to 1 in 5 cases of diabetes (18%) in Greenland are associated with high-impact genetic variants compared to 1-3% in large populations.

Funding

Novo Nordisk Foundation, Independent Research Fund Denmark, and Karen Elise Jensen's Foundation.

Life-threatening viral disease in a novel form of autosomal recessive IFNAR2 deficiency in the Arctic

Type I interferons (IFN-I) play a critical role in human antiviral immunity, as demonstrated by the exceptionally rare deleterious variants of IFNAR1 or IFNAR2. We investigated five children from Greenland, Canada, and Alaska presenting with viral diseases, including life-threatening COVID-19 or influenza, in addition to meningoencephalitis and/or hemophagocytic lymphohistiocytosis following live-attenuated viral vaccination. The affected individuals bore the same homozygous IFNAR2 c.157T>C, p.Ser53Pro missense variant. Although absent from reference databases, p.Ser53Pro occurred with a minor allele frequency of 0.034 in their Inuit ancestry. The serine to proline substitution prevented cell surface expression of IFNAR2 protein, small amounts of which persisted intracellularly in an aberrantly glycosylated state. Cells exclusively expressing the p.Ser53Pro variant lacked responses to recombinant IFN-I and displayed heightened vulnerability to multiple viruses in vitro-a phenotype rescued by wild-type IFNAR2 complementation. This novel form of autosomal recessive IFNAR2 deficiency reinforces the essential role of IFN-I in viral immunity. Further studies are warranted to assess the need for population screening.

Genomes from Verteba cave suggest diversity within the Trypillians in Ukraine

The transition to agriculture occurred relatively late in Eastern Europe, leading researchers to debate whether it was a gradual, interactive process or a colonisation event. In the forest and forest-steppe regions of Ukraine, farming appeared during the fifth millennium BCE, associated with the Cucuteni-Trypillia cultural complex (CTCC, ~ 5000–3000 BCE). Across Europe, the Neolithisation process was highly variable across space and over time. Here, we investigate the population dynamics of early agriculturalists from the eastern forest-steppe region based on the analyses of 20 ancient genomes from the site of Verteba Cave (3935–825 cal BCE). Results reveal that the CTCC individuals’ ancestry is related to both western hunter-gatherers and Near Eastern farmers, has no local ancestry associated with Ukrainian Neolithic hunter-gatherers and has steppe ancestry. An Early Bronze Age individual has an ancestry profile related to the Yamnaya expansions but with 20% of ancestry related to the other Trypillian individuals, which suggests admixture between the Trypillians and the incoming populations carrying steppe-related ancestry. A Late Bronze Age individual dated to 980–825 cal BCE has a genetic profile indicating affinity to Beaker-related populations, detected close to 1000 years after the end of the Bell Beaker phenomenon during the third millennium BCE.

Loss of Sucrase-Isomaltase Function Increases Acetate Levels and Improves Metabolic Health in Greenlandic Cohorts

BACKGROUND & AIMS

The sucrase-isomaltase (SI) c.273_274delAG loss-of-function variant is common in Arctic populations and causes congenital sucrase-isomaltase deficiency, which is an inability to break down and absorb sucrose and isomaltose. Children with this condition experience gastrointestinal symptoms when dietary sucrose is introduced. We aimed to describe the health of adults with sucrase-isomaltase deficiency.

METHODS

The association between c.273_274delAG and phenotypes related to metabolic health was assessed in 2 cohorts of Greenlandic adults (n = 4922 and n = 1629). A sucrase-isomaltase knockout (Sis-KO) mouse model was used to further elucidate the findings.

RESULTS

Homozygous carriers of the variant had a markedly healthier metabolic profile than the remaining population, including lower body mass index (β [standard error], -2.0 [0.5] kg/m²; P = 3.1 × 10^-5), body weight (-4.8 [1.4] kg; P = 5.1 × 10^-4), fat percentage (-3.3% [1.0%]; P = 3.7 × 10^-4), fasting triglyceride (-0.27 [0.07] mmol/L; P = 2.3 × 10^-6), and remnant cholesterol (-0.11 [0.03] mmol/L; P = 4.2 × 10^-5). Further analyses suggested that this was likely mediated partly by higher circulating levels of acetate observed in homozygous carriers (β [standard error], 0.056 [0.002] mmol/L; P = 2.1 × 10^-26), and partly by reduced sucrose uptake, but not lower caloric intake. These findings were verified in Sis-KO mice, which, compared with wild-type mice, were leaner on a sucrose-containing diet, despite similar caloric intake, had significantly higher plasma acetate levels in response to a sucrose gavage, and had lower plasma glucose level in response to a sucrose-tolerance test.

CONCLUSIONS

These results suggest that sucrase-isomaltase constitutes a promising drug target for improvement of metabolic health, and that the health benefits are mediated by reduced dietary sucrose uptake and possibly also by higher levels of circulating acetate.

A unified genealogy of modern and ancient genomes

The sequencing of modern and ancient genomes from around the world has revolutionized our understanding of human history and evolution. However, the problem of how best to characterize ancestral relationships from the totality of human genomic variation remains unsolved. Here, we address this challenge with nonparametric methods that enable us to infer a unified genealogy of modern and ancient humans. This compact representation of multiple datasets explores the challenges of missing and erroneous data and uses ancient samples to constrain and date relationships. We demonstrate the power of the method to recover relationships between individuals and populations as well as to identify descendants of ancient samples. Finally, we introduce a simple nonparametric estimator of the geographical location of ancestors that recapitulates key events in human history.

Variation of Genomic Sites Associated with Severe Covid-19 Across Populations: Global and National Patterns

Background

Information about the distribution of clinically significant genetic markers in different populations may be helpful in elaborating personalized approaches to the clinical management of COVID-19 in the absence of consensus guidelines.

Aim

Analyze frequencies and distribution patterns of two markers associated with severe COVID-19 (rs11385942 and rs657152) and look for potential correlations between these markers and deaths from COVID-19 among populations in Russia and across the world.

Methods

We genotyped 1883 samples from 91 ethnic groups pooled into 28 populations representing Russia and its neighbor states. We also compiled a dataset on 32 populations from other regions using genotypes extracted or imputed from the available databases. Geographic maps showing the frequency distribution of the analyzed markers were constructed using the obtained data.

Results

The cartographic analysis revealed that rs11385942 distribution follows the West Eurasian pattern: the marker is frequent among the populations of Europe, West Asia and South Asia but rare or absent in all other parts of the globe. Notably, the transition from high to low rs11385942 frequencies across Eurasia is not abrupt but follows the clinal variation pattern instead. The distribution of rs657152 is more homogeneous. The analysis of correlations between the frequencies of the studied markers and the epidemiological characteristics of COVID-19 in a population revealed that higher frequencies of both risk alleles correlated positively with mortality from this disease. For rs657152, the correlation was especially strong (r = 0.59, p = 0.02). These reasonable correlations were observed for the "Russian" dataset only: no such correlations were established for the "world" dataset. This could be attributed to the differences in methodology used to collect COVID-19 statistics in different countries.

Conclusion

Our findings suggest that genetic differences between populations make a small yet tangible contribution to the heterogeneity of the pandemic worldwide.

The influence of a lost society, the Sadlermiut, on the environment in the Canadian Arctic

High latitude freshwater ecosystems are sentinels of human activity and environmental change. The lakes and ponds that characterize Arctic landscapes have a low resilience to buffer variability in climate, especially with increasing global anthropogenic stressors in recent decades. Here, we show that a small freshwater pond in proximity of the archaeological site “Native Point” on Southampton Island (Nunavut, Arctic Canada) is a highly sensitive environmental recorder. The sediment analyses allowed for pinpointing the first arrival of Sadlermiut culture at Native Point to ~ 1250 CE, followed by a dietary shift likely in response to the onset of cooling in the region ~ 1400 CE. The influence of the Sadlermiut on the environment persisted long after the last of their population perished in 1903. Presently, the pond remains a distorted ecosystem that has experienced fundamental shifts in the benthic invertebrate assemblages and accumulated anthropogenic metals in the sediment. Our multi-proxy paleolimnological investigation using geochemical and biological indicators emphasizes that direct and indirect anthropogenic impacts have long-term environmental implications on high latitude ecosystems.

Vulnerability of the North Water ecosystem to climate change

High Arctic ecosystems and Indigenous livelihoods are tightly linked and exposed to climate change, yet assessing their sensitivity requires a long-term perspective. Here, we assess the vulnerability of the North Water polynya, a unique seaice ecosystem that sustains the world’s northernmost Inuit communities and several keystone Arctic species. We reconstruct mid-to-late Holocene changes in sea ice, marine primary production, and little auk colony dynamics through multi-proxy analysis of marine and lake sediment cores. Our results suggest a productive ecosystem by 4400–4200 cal yrs b2k coincident with the arrival of the first humans in Greenland. Climate forcing during the late Holocene, leading to periods of polynya instability and marine productivity decline, is strikingly coeval with the human abandonment of Greenland from c. 2200–1200 cal yrs b2k. Our long-term perspective highlights the future decline of the North Water ecosystem, due to climate warming and changing sea-ice conditions, as an important climate change risk.

The North Water polynya is a unique but vulnerable ecosystem, home to Indigenous people and Arctic keystone species. New palaeoecological records from Greenland suggest human abandonment c. 2200–1200 cal yrs BP occurred during climate-forced polynya instability, foreshadowing future ecosystem declines.

Uniparental genetic markers in Native Americans: A summary of all available data from ancient and contemporary populations

OBJECTIVES

The aim of this study was to create a comprehensive summary of available mtDNA and Y-chromosome data for Native Americans from North, Central, and South America, including both modern and ancient DNA. To illustrate the usefulness of this dataset we present a broad picture of the genetic variation for both markers across the Americas.

METHODS

We searched PubMed, ResearchGate, Google Scholar for studies about mtDNA or Y-chromosome variation in Native American populations, including geographic, linguistic, ecological (ecoregion), archeological and chronological information. We used AMOVA to estimate the genetic structure associated with language and ecoregion grouping and Mantel tests to evaluate the correlation between genetic and geographic distances.

RESULTS

Genetic data were obtained from 321 primary sources, including 22,569 individuals from 298 contemporary populations, and 3628 individuals from 202 archeological populations. MtDNA lineages of probable non-Amerindian origin were rare, in contrast with Y-chromosome lineages. Mantel tests showed a statistically significant correlation for the whole continent considering mtDNA but not the Y-chromosome. Genetic structure between groups was always stronger for mtDNA than for the Y-chromosome.

CONCLUSIONS

This study summarizes decades of research conducted in Native American populations for both mtDNA and the Y-chromosome. Continental or sub-continental patterns of variation reveal that most of the genetic variation occurs within populations rather than among linguistic or ecoregional groups, and that isolation by distance is barely detectable in most population sets. The genetic structure among groups was always larger for mtDNA than for the Y-chromosome, suggesting between-sex differences in gene flow.

Genomic Insights Into the Admixture History of Mongolic- and Tungusic-Speaking Populations From Southwestern East Asia

As a major part of the modern Trans-Eurasian or Altaic language family, most of the Mongolic and Tungusic languages were mainly spoken in northern China, Mongolia, and southern Siberia, but some were also found in southern China. Previous genetic surveys only focused on the dissection of genetic structure of northern Altaic-speaking populations; however, the ancestral origin and genomic diversification of Mongolic and Tungusic-speaking populations from southwestern East Asia remain poorly understood because of the paucity of high-density sampling and genome-wide data. Here, we generated genome-wide data at nearly 700,000 single-nucleotide polymorphisms (SNPs) in 26 Mongolians and 55 Manchus collected from Guizhou province in southwestern China. We applied principal component analysis (PCA), ADMIXTURE, f statistics, qpWave/qpAdm analysis, qpGraph, TreeMix, Fst, and ALDER to infer the fine-scale population genetic structure and admixture history. We found significant genetic differentiation between northern and southern Mongolic and Tungusic speakers, as one specific genetic cline of Manchu and Mongolian was identified in Guizhou province. Further results from ADMIXTURE and f statistics showed that the studied Guizhou Mongolians and Manchus had a strong genetic affinity with southern East Asians, especially for inland southern East Asians. The qpAdm-based estimates of ancestry admixture proportion demonstrated that Guizhou Mongolians and Manchus people could be modeled as the admixtures of one northern ancestry related to northern Tungusic/Mongolic speakers or Yellow River farmers and one southern ancestry associated with Austronesian, Tai-Kadai, and Austroasiatic speakers. The qpGraph-based phylogeny and neighbor-joining tree further confirmed that Guizhou Manchus and Mongolians derived approximately half of the ancestry from their northern ancestors and the other half from southern Indigenous East Asians. The estimated admixture time ranged from 600 to 1,000 years ago, which further confirmed the admixture events were mediated via the Mongolians Empire expansion during the formation of the Yuan dynasty.

Helicobacter pylori's historical journey through Siberia and the Americas

The peopling of Siberia and the Americas is intriguing for archaeologists, linguists, and human geneticists, but despite significant recent developments, many details remain controversial. Here, we provide insights based on genetic diversity within Helicobacter pylori, a bacterium that infects 50% of all humans. H. pylori strains were collected from across eastern Eurasia and the Americas. Sequence analyses indicated that Siberia contains both anciently diverged and recently admixed bacteria, supporting both human persistence over the last glacial maximum and more recent human recolonization. We inferred a single migration across the Bering land bridge, accompanied by a dramatic reduction in effective population size, followed by bidirectional Holocene gene flow between Asia and the Americas.

The gastric bacterium Helicobacter pylori shares a coevolutionary history with humans that predates the out-of-Africa diaspora, and the geographical specificities of H. pylori populations reflect multiple well-known human migrations. We extensively sampled H. pylori from 16 ethnically diverse human populations across Siberia to help resolve whether ancient northern Eurasian populations persisted at high latitudes through the last glacial maximum and the relationships between present-day Siberians and Native Americans. A total of 556 strains were cultivated and genotyped by multilocus sequence typing, and 54 representative draft genomes were sequenced. The genetic diversity across Eurasia and the Americas was structured into three populations: hpAsia2, hpEastAsia, and hpNorthAsia. hpNorthAsia is closely related to the subpopulation hspIndigenousAmericas from Native Americans. Siberian bacteria were structured into five other subpopulations, two of which evolved through a divergence from hpAsia2 and hpNorthAsia, while three originated though Holocene admixture. The presence of both anciently diverged and recently admixed strains across Siberia support both Pleistocene persistence and Holocene recolonization. We also show that hspIndigenousAmericas is endemic in human populations across northern Eurasia. The evolutionary history of hspIndigenousAmericas was reconstructed using approximate Bayesian computation, which showed that it colonized the New World in a single migration event associated with a severe demographic bottleneck followed by low levels of recent admixture across the Bering Strait.

Diet and physical activity in Greenland: genetic interactions and associations with obesity and diabetes

The Inuit in Greenland have gone through dramatic lifestyle changes during the last half century. More time is spent being sedentary and imported foods replaces traditional foods like seal and whale. The population has also experienced a rapid growth in obesity and metabolic disturbances and diabetes is today common despite being almost unknown few decades ago. In this paper, we describe and discuss the role of lifestyle changes and genetics for Inuit metabolic health. Novelty: Cardiometabolic disease risk has increased in Greenland. Lifestyle changes and possibly gene-lifestyle interactions play a role.

Peopling of the Americas as inferred from ancient genomics

In less than a decade, analyses of ancient genomes have transformed our understanding of the Indigenous peopling and population history of the Americas. These studies have shown that this history, which began in the late Pleistocene epoch and continued episodically into the Holocene epoch, was far more complex than previously thought. It is now evident that the initial dispersal involved the movement from northeast Asia of distinct and previously unknown populations, including some for whom there are no currently known descendants. The first peoples, once south of the continental ice sheets, spread widely, expanded rapidly and branched into multiple populations. Their descendants-over the next fifteen millennia-experienced varying degrees of isolation, admixture, continuity and replacement, and their genomes help to illuminate the relationships among major subgroups of Native American populations. Notably, all ancient individuals in the Americas, save for later-arriving Arctic peoples, are more closely related to contemporary Indigenous American individuals than to any other population elsewhere, which challenges the claim-which is based on anatomical evidence-that there was an early, non-Native American population in the Americas. Here we review the patterns revealed by ancient genomics that help to shed light on the past peoples who created the archaeological landscape, and together lead to deeper insights into the population and cultural history of the Americas.

The genetic history of Greenlandic-European contact

The Inuit ancestors of the Greenlandic people arrived in Greenland close to 1,000 years ago.¹ Since then, Europeans from many different countries have been present in Greenland. Consequently, the present-day Greenlandic population has ~25% of its genetic ancestry from Europe.² In this study, we investigated to what extent different European countries have contributed to this genetic ancestry. We combined dense SNP chip data from 3,972 Greenlanders and 8,275 Europeans from 14 countries and inferred the ancestry contribution from each of these 14 countries using haplotype-based methods. Due to the rapid increase in population size in Greenland over the past ~100 years, we hypothesized that earlier European interactions, such as pre-colonial Dutch whalers and early German and Danish-Norwegian missionaries, as well as the later Danish colonists and post-colonial immigrants, all contributed European genetic ancestry. However, we found that the European ancestry is almost entirely Danish and that a substantial fraction is from admixture that took place within the last few generations.

The Greenlandic Inuit have had extensive historical contact with Europeans, and the present-day Greenlandic population has substantial amounts of European ancestry. Waples et al. use genetic data to investigate the origin of this ancestry. They show that much of it is Danish and find little evidence of it being from pre-colonial European contact.

The genomic history of the Aegean palatial civilizations

The Cycladic, the Minoan, and the Helladic (Mycenaean) cultures define the Bronze Age (BA) of Greece. Urbanism, complex social structures, craft and agricultural specialization, and the earliest forms of writing characterize this iconic period. We sequenced six Early to Middle BA whole genomes, along with 11 mitochondrial genomes, sampled from the three BA cultures of the Aegean Sea. The Early BA (EBA) genomes are homogeneous and derive most of their ancestry from Neolithic Aegeans, contrary to earlier hypotheses that the Neolithic-EBA cultural transition was due to massive population turnover. EBA Aegeans were shaped by relatively small-scale migration from East of the Aegean, as evidenced by the Caucasus-related ancestry also detected in Anatolians. In contrast, Middle BA (MBA) individuals of northern Greece differ from EBA populations in showing ∼50% Pontic-Caspian Steppe-related ancestry, dated at ca. 2,600-2,000 BCE. Such gene flow events during the MBA contributed toward shaping present-day Greek genomes.

Predicting sample success for large-scale ancient DNA studies on marine mammals

In recent years, nonhuman ancient DNA studies have begun to focus on larger sample sizes and whole genomes, offering the potential to reveal exciting and hitherto unknown answers to ongoing biological and archaeological questions. However, one major limitation to such studies is the substantial financial and time investments still required during sample screening, due to uncertainty regarding successful sample selection. This study investigates the effect of a wide range of sample properties including latitude, sample age, skeletal element, collagen preservation, and context on endogenous content and DNA damage profiles for 317 ancient and historic pinniped samples collected from across the North Atlantic and surrounding regions. Using generalised linear and mixed-effect models, we found that a range of factors affected DNA preservation within each of the species under consideration. The most important findings were that endogenous content varied significantly within species according to context, the type of skeletal element, the collagen content and collection year. There also appears to be an effect of the sample's geographic origin, with samples from the Arctic generally showing higher endogenous content and lower damage rates. Both latitude and sample age were found to have significant relationships with damage levels, but only for walrus samples. Sex, ontogenetic age and extraction material preparation were not found to have any significant relationship with DNA preservation. Overall, skeletal element and sample context were found to be the most influential factors and should therefore be considered when selecting samples for large-scale ancient genome studies.

Population History and Gene Divergence in Native Mexicans Inferred from 76 Human Exomes

Native American genetic variation remains underrepresented in most catalogs of human genome sequencing data. Previous genotyping efforts have revealed that Mexico’s Indigenous population is highly differentiated and substructured, thus potentially harboring higher proportions of private genetic variants of functional and biomedical relevance. Here we have targeted the coding fraction of the genome and characterized its full site frequency spectrum by sequencing 76 exomes from five Indigenous populations across Mexico. Using diffusion approximations, we modeled the demographic history of Indigenous populations from Mexico with northern and southern ethnic groups splitting 7.2 KYA and subsequently diverging locally 6.5 and 5.7 KYA, respectively. Selection scans for positive selection revealed BCL2L13 and KBTBD8 genes as potential candidates for adaptive evolution in Rarámuris and Triquis, respectively. BCL2L13 is highly expressed in skeletal muscle and could be related to physical endurance, a well-known phenotype of the northern Mexico Rarámuri. The KBTBD8 gene has been associated with idiopathic short stature and we found it to be highly differentiated in Triqui, a southern Indigenous group from Oaxaca whose height is extremely low compared to other Native populations.

Founder variants and population genomes-Toward precision medicine

Human migration and community specific cultural practices have contributed to founder events and enrichment of the variants associated with genetic diseases. While many founder events in isolated populations have remained uncharacterized, the application of genomics in clinical settings as well as for population scale studies in the recent years have provided an unprecedented push towards identification of founder variants associated with human health and disease. The discovery and characterization of founder variants could have far reaching implications not only in understanding the history or genealogy of the disease, but also in implementing evidence based policies and genetic testing frameworks. This further enables precise diagnosis and prevention in an attempt towards precision medicine. This review provides an overview of founder variants along with methods and resources cataloging them. We have also discussed the public health implications and examples of prevalent disease associated founder variants in specific populations.

Children and innovation: play, play objects and object play in cultural evolution

Cultural evolutionary theory conceptualises culture as an information-transmission system whose dynamics take on evolutionary properties. Within this framework, however, innovation has been likened to random mutations, reducing its occurrence to chance or fortuitous transmission error. In introducing the special collection on children and innovation, we here place object play and play objects - especially functional miniatures - from carefully chosen archaeological contexts in a niche construction perspective. Given that play, including object play, is ubiquitous in human societies, we suggest that plaything construction, provisioning and use have, over evolutionary timescales, paid substantial selective dividends via ontogenetic niche modification. Combining findings from cognitive science, ethology and ethnography with insights into hominin early developmental life-history, we show how play objects and object play probably had decisive roles in the emergence of innovative capabilities. Importantly, we argue that closer attention to play objects can go some way towards addressing changes in innovation rates that occurred throughout human biocultural evolution and why innovations are observable within certain technological domains but not others.

AuthentiCT: a model of ancient DNA damage to estimate the proportion of present-day DNA contamination

Contamination from present-day DNA is a fundamental issue when studying ancient DNA from historical or archaeological material, and quantifying the amount of contamination is essential for downstream analyses. We present AuthentiCT, a command-line tool to estimate the proportion of present-day DNA contamination in ancient DNA datasets generated from single-stranded DNA libraries. The prediction is based solely on the patterns of post-mortem damage observed on ancient DNA sequences. The method has the power to quantify contamination from as few as 10,000 mapped sequences, making it particularly useful for analysing specimens that are poorly preserved or for which little data is available.

Ancient Jomon genome sequence analysis sheds light on migration patterns of early East Asian populations

Anatomically modern humans reached East Asia more than 40,000 years ago. However, key questions still remain unanswered with regard to the route(s) and the number of wave(s) in the dispersal into East Eurasia. Ancient genomes at the edge of the region may elucidate a more detailed picture of the peopling of East Eurasia. Here, we analyze the whole-genome sequence of a 2,500-year-old individual (IK002) from the main-island of Japan that is  characterized with a typical Jomon culture. The phylogenetic analyses support multiple waves of migration, with IK002 forming a basal lineage to the East and Northeast Asian genomes examined, likely representing some of the earliest-wave migrants who went north from Southeast Asia to East Asia. Furthermore, IK002 shows strong genetic affinity with the indigenous Taiwan aborigines, which may support a coastal route of the Jomon-ancestry migration. This study highlights the power of ancient genomics to provide new insights into the complex history of human migration into East Eurasia.

Takashi Gakuhari, Shigeki Nakagome et al. report the genomic analysis on a 2.5 kya individual from the ancient Jomon culture in present-day Japan. Phylogenetic analysis with comparison to other Eurasian sequences suggests early migration patterns in Asia and provides insight into the genetic affinities between peoples of the region.

Robust genome-wide ancestry inference for heterogeneous datasets: illustrated using the 1,000 genome project with 3D facial images

Estimates of individual-level genomic ancestry are routinely used in human genetics, and related fields. The analysis of population structure and genomic ancestry can yield insights in terms of modern and ancient populations, allowing us to address questions regarding admixture, and the numbers and identities of the parental source populations. Unrecognized population structure is also an important confounder to correct for in genome-wide association studies. However, it remains challenging to work with heterogeneous datasets from multiple studies collected by different laboratories with diverse genotyping and imputation protocols. This work presents a new approach and an accompanying open-source toolbox that facilitates a robust integrative analysis for population structure and genomic ancestry estimates for heterogeneous datasets. We show robustness against individual outliers and different protocols for the projection of new samples into a reference ancestry space, and the ability to reveal and adjust for population structure in a simulated case–control admixed population. Given that visually evident and easily recognizable patterns of human facial characteristics co-vary with genomic ancestry, and based on the integration of three different sources of genome data, we generate average 3D faces to illustrate genomic ancestry variations within the 1,000 Genome project and for eight ancient-DNA profiles, respectively.

The history behind the mosaic of the Americas

Focusing on literature published in 2018-2020, we review inferences about: (i) how ancient DNA is contributing to clarify the peopling of the Americas and the dispersal of its first inhabitants, (ii) how the interplay between environmental diversity and culture has influenced the genetic structure and adaptation of Andean and Amazon populations, (iii) how genetics has contributed to our understanding of the Pre-Columbian Tupi expansion in Eastern South America, (iv) the subcontinental origins and dynamics of Post-Columbian admixture in the Americas, and finally, (v) episodes of adaptive natural selection in the American continent, particularly in the high altitudes of the Andes.

Arctic-adapted dogs emerged at the Pleistocene-Holocene transition

Although sled dogs are one of the most specialized groups of dogs, their origin and evolution has received much less attention than many other dog groups. We applied a genomic approach to investigate their spatiotemporal emergence by sequencing the genomes of 10 modern Greenland sled dogs, an ~9500-year-old Siberian dog associated with archaeological evidence for sled technology, and an ~33,000-year-old Siberian wolf. We found noteworthy genetic similarity between the ancient dog and modern sled dogs. We detected gene flow from Pleistocene Siberian wolves, but not modern American wolves, to present-day sled dogs. The results indicate that the major ancestry of modern sled dogs traces back to Siberia, where sled dog-specific haplotypes of genes that potentially relate to Arctic adaptation were established by 9500 years ago.

Genetic study of the Arctic CPT1A variant suggests that its effect on fatty acid levels is modulated by traditional Inuit diet

Several recent studies have found signs of recent selection on the carnitine palmitoyl-transferase 1A (CPT1A) gene in the ancestors of Arctic populations likely as a result of their traditional diet. CPT1A is involved in fatty acid transportation and is known to affect circulating fatty acid profiles in Inuit as does the unique traditional diet rich in marine animals. We aimed to assess which fatty acids may have driven the selection of rs80356779, a c.1436C>T (p.(Pro479Leu)) variant in CPT1A, by analyzing a potential interaction between the variant and traditional Inuit diet. We included 3005 genome-wide genotyped individuals living in Greenland, who had blood cell membrane fatty acid levels measured. Consumption of 25 traditional food items was expressed as percentage of total energy intake. We tested for CPT1A × traditional diet interaction while taking relatedness and admixture into account. Increasing intakes of traditional diet was estimated to attenuate the effect of 479L on 20:3 omega-6 levels (p = 0.000399), but increase the effect of the variant on 22:5 omega-3 levels (p = 0.000963). The 479L effect on 22:5 omega-3 more than doubled in individuals with a high intake of traditional diet (90% percentile) compared with individuals with a low intake (10% percentile). Similar results were found when assessing interactions with marine foods. Our results suggest that the association between traditional diet and blood cell fatty acid composition is affected by the CPT1A genotype, or other variants in linkage disequilibrium, and support the hypothesis that omega-3 fatty acids may have been important for adaptation to the Arctic diet.

Muskox status, recent variation, and uncertain future

Muskoxen (Ovibos moschatus) are an integral component of Arctic biodiversity. Given low genetic diversity, their ability to respond to future and rapid Arctic change is unknown, although paleontological history demonstrates adaptability within limits. We discuss status and limitations of current monitoring, and summarize circumpolar status and recent variations, delineating all 55 endemic or translocated populations. Acknowledging uncertainties, global abundance is ca 170 000 muskoxen. Not all populations are thriving. Six populations are in decline, and as recently as the turn of the century, one of these was the largest population in the world, equaling ca 41% of today's total abundance. Climate, diseases, and anthropogenic changes are likely the principal drivers of muskox population change and result in multiple stressors that vary temporally and spatially. Impacts to muskoxen are precipitated by habitat loss/degradation, altered vegetation and species associations, pollution, and harvest. Which elements are relevant for a specific population will vary, as will their cumulative interactions. Our summaries highlight the importance of harmonizing existing data, intensifying long-term monitoring efforts including demographics and health assessments, standardizing and implementing monitoring protocols, and increasing stakeholder engagement/contributions.

Genomic landscape of the signals of positive natural selection in populations of Northern Eurasia: A view from Northern Russia

Natural selection of beneficial genetic variants played a critical role in human adaptation to a wide range of environmental conditions. Northern Eurasia, despite its severe climate, is home to lots of ethnically diverse populations. The genetic variants associated with the survival of these populations have hardly been analyzed. We searched for the genomic signatures of positive selection in (1) the genome-wide microarray data of 432 people from eight different northern Russian populations and (2) the whole-genome sequences of 250 people from Northern Eurasia from a public repository through testing the extended haplotype homozigosity (EHH) and direct comparison of allele frequency, respectively. The 20 loci with the strongest selection signals were characterized in detail. Among the top EHH hits were the NRG3 and NBEA genes, which are involved in the development and functioning of the neural system, the PTPRM gene, which mediates cell-cell interactions and adhesion, and a region on chromosome 4 (chr4:28.7-28.9 Mb) that contained several loci affiliated with different classes of non-coding RNAs (RN7SL101P, MIR4275, MESTP3, and LINC02364). NBEA and the region on chromosome 4 were novel selection targets that were identified for the first time in Western Siberian populations. Cross-population comparisons of EHH profiles suggested a particular role for the chr4:28.7-28.9 Mb region in the local adaptation of Western Siberians. The strongest selection signal identified in Siberian sequenced genomes was formed by six SNPs on chromosome 11 (chr11:124.9-125.2 Mb). This region included well-known genes SLC37A2 and PKNOX2. SLC37A2 is most-highly expressed in the gut. Its expression is regulated by vitamin D, which is often deficient in northern regions. The PKNOX2 gene is a transcription factor of the homeobox family that is expressed in the brain and many other tissues. This gene is associated with alcohol addiction, which is widespread in many Northern Eurasian populations.

Molecular analysis of an ancient Thule population at Nuvuk, Point Barrow, Alaska

OBJECTIVES

The North American archaeological record supports a Holocene origin of Arctic Indigenous peoples. Although the Paleo-Inuit were present for millennia, archaeological and genetic studies suggest that modern peoples descend from a second, more recent tradition known as the Neo-Inuit. Origins of the Neo-Inuit and their relations to the earlier and later Indigenous peoples are an area of active study. Here, we genetically analyze the maternal lineages present at Nuvuk, once the northernmost community in Alaska and located in a region identified as a possible origin point of the Neo-Inuit Thule. The cemetery at Nuvuk contains human remains representing a nearly one thousand year uninterrupted occupation from early Thule to post-contact Iñupiat.

MATERIALS AND METHODS

We selected 44 individuals from Nuvuk with calibrated dates between 981 AD and 1885 AD for molecular analysis. We amplified and sequenced the hypervariable segment I of the mitogenome. We compared the Nuvuk data with previously published sequences from 68 modern and ancient communities from across Asia and North America. Phylogeographic analyses suggest possible scenarios of Holocene Arctic and sub-Arctic population movements.

RESULTS

We successfully retrieved sequence data from 39 individuals. Haplogroup frequencies in Nuvuk were typed as 66.7% A2b1, 25.6% A2a, and 7.7% D4b1a2a1a. These results suggest that the population at Nuvuk was closest to the ancient Thule and modern Inuit of Canada, and to the Siberian Naukan people. We confirm that haplogroups A2a, A2b1, D2a, and D4b1a2a1a appear at high frequency in Arctic and sub-Arctic populations of North America and Chukotka. Sister clades D2b and D4b1a2a1b are present in Asian and Eastern European populations.

DISCUSSION

The ancient mitochondrial sequences from Nuvuk confirm the link between the North Slope and the Thule who later spread east, and the maternal discontinuity between the Neo-Inuit and Paleo-Inuit. We suggest haplogroups A2a, A2b, and D4b1a2a1a are linked to the ancestors of the Thule in eastern Beringia, whereas the D2 and D4b1a2a1 clades appear to have Asian Holocene origins. Further Siberian and Alaskan genomes are necessary to clarify these population migrations beyond a simple two-wave scenario of Neo-Inuit and Paleo-Inuit.

Specialized sledge dogs accompanied Inuit dispersal across the North American Arctic

Domestic dogs have been central to life in the North American Arctic for millennia. The ancestors of the Inuit were the first to introduce the widespread usage of dog sledge transportation technology to the Americas, but whether the Inuit adopted local Palaeo-Inuit dogs or introduced a new dog population to the region remains unknown. To test these hypotheses, we generated mitochondrial DNA and geometric morphometric data of skull and dental elements from a total of 922 North American Arctic dogs and wolves spanning over 4500 years. Our analyses revealed that dogs from Inuit sites dating from 2000 BP possess morphological and genetic signatures that distinguish them from earlier Palaeo-Inuit dogs, and identified a novel mitochondrial clade in eastern Siberia and Alaska. The genetic legacy of these Inuit dogs survives today in modern Arctic sledge dogs despite phenotypic differences between archaeological and modern Arctic dogs. Together, our data reveal that Inuit dogs derive from a secondary pre-contact migration of dogs distinct from Palaeo-Inuit dogs, and probably aided the Inuit expansion across the North American Arctic beginning around 1000 BP.

Arctic Archaeology and Climate Change

An enduring debate in the field of Arctic archaeology has been the extent to which climate change impacted cultural developments in the past. Long-term culture change across the circumpolar Arctic was often highly dynamic, with episodes of rapid migration, regional abandonment, and—in some cases—the disappearance or wholesale replacement of entire cultural traditions. By the 1960s, researchers were exploring the possibility that warming episodes had positive effects on cold-adapted premodern peoples in the Arctic by ( a) reducing the extent of sea ice, ( b) expanding the size and range of marine mammal populations, and ( c) opening new waterways and hunting areas for marine-adapted human groups. Although monocausal climatic arguments for change are now regarded as overly simplistic, the growing threat of contemporary Arctic warming to Indigenous livelihoods has given wider relevance to research into long-term culture–climate interactions. With their capacity to examine deeper cultural responses to climate change, archaeologists are in a unique position to generate human-scale climate adaptation insights that may inform future planning and mitigation efforts. The exceptionally well-preserved cultural and paleo-ecological sequences of the Arctic make it one of the best-suited regions on Earth to address such problems. Ironically, while archaeologists employ an exciting and highly promising new generation of methods and approaches to examine long-term fragility and resilience in Arctic social-ecological systems, many of these frozen paleo-societal archives are fast disappearing due to anthropogenic warming.

New middle chapter in the story of human evolution

Analyzing genomic data from ancient humans illuminates South Asian ancestry

Genetic architecture and adaptations of Nunavik Inuit

The Canadian Inuit have a distinct population background that may entail particular implications for the health of its individuals. However, the number of genetic studies examining this Inuit population is limited, and much remains to be discovered in regard to its genetic characteristics. In this study, we generated whole-exome sequences and genomewide genotypes for 170 Nunavik Inuit, a small and isolated founder population of Canadian Arctic indigenous people. Our study revealed the genetic background of Nunavik Inuit to be distinct from any known present-day population. The majority of Nunavik Inuit show little evidence of gene flow from European or present-day Native American peoples, and Inuit living around Hudson Bay are genetically distinct from those around Ungava Bay. We also inferred that Nunavik Inuit have a small effective population size of 3,000 and likely split from Greenlandic Inuit ∼10.5 kya. Nunavik Inuit went through a bottleneck at approximately the same time and might have admixed with a population related to the Paleo-Eskimos. Our study highlights population-specific genomic signatures in coding regions that show adaptations unique to Nunavik Inuit, particularly in pathways involving fatty acid metabolism and cellular adhesion (CPNE7, ICAM5, STAT2, and RAF1). Subsequent analyses in selection footprints and the risk of intracranial aneurysms (IAs) in Nunavik Inuit revealed an exonic variant under weak negative selection to be significantly associated with IA (rs77470587; P = 4.6 × 10^-8).

Predicting the loss of organic archaeological deposits at a regional scale in Greenland

Across the Arctic, microbial degradation is actively destroying irreplaceable cultural and environmental records that have been preserved within archaeological deposits for millennia. Because it is not possible to survey the many sites in this remote part of the world, new methods are urgently needed to detect and assess the potential degradation. Here, we investigate organic deposits at seven archaeological sites located along the dominating west-east climatic gradient in West Greenland. We show that, regardless of age, depositional history and environmental conditions, all organic deposits are highly vulnerable to degradation. A state-of-the-art model that simulates the effect of future climate change on degradation indicates that 30–70% of the archaeological fraction of organic carbon (OC) could disappear within the next 80 years. This range reflects the variation within the climatic gradient and the future climate scenario applied (RCP 4.5 and RCP 8.5). All archaeological deposits are expected to experience a substantial loss, but the most rapid degradation seems to occur in the continental inland areas of the region, dominated by dry and warm summers. This suggests that organic remains from the Norse Viking Age settlers are especially under threat in the coming years.

The population history of northeastern Siberia since the Pleistocene

Northeastern Siberia has been inhabited by humans for more than 40,000 years but its deep population history remains poorly understood. Here we investigate the late Pleistocene population history of northeastern Siberia through analyses of 34 newly recovered ancient genomes that date to between 31,000 and 600 years ago. We document complex population dynamics during this period, including at least three major migration events: an initial peopling by a previously unknown Palaeolithic population of 'Ancient North Siberians' who are distantly related to early West Eurasian hunter-gatherers; the arrival of East Asian-related peoples, which gave rise to 'Ancient Palaeo-Siberians' who are closely related to contemporary communities from far-northeastern Siberia (such as the Koryaks), as well as Native Americans; and a Holocene migration of other East Asian-related peoples, who we name 'Neo-Siberians', and from whom many contemporary Siberians are descended. Each of these population expansions largely replaced the earlier inhabitants, and ultimately generated the mosaic genetic make-up of contemporary peoples who inhabit a vast area across northern Eurasia and the Americas.

Palaeo-Eskimo genetic ancestry and the peopling of Chukotka and North America

Much of the American Arctic was first settled 5,000 years ago, by groups of people known as Palaeo-Eskimos. They were subsequently joined and largely displaced around 1,000 years ago by ancestors of the present-day Inuit and Yup'ik1-3. The genetic relationship between Palaeo-Eskimos and Native American, Inuit, Yup'ik and Aleut populations remains uncertain4-6. Here we present genomic data for 48 ancient individuals from Chukotka, East Siberia, the Aleutian Islands, Alaska, and the Canadian Arctic. We co-analyse these data with data from present-day Alaskan Iñupiat and West Siberian populations and published genomes. Using methods based on rare-allele and haplotype sharing, as well as established techniques4,7-9, we show that Palaeo-Eskimo-related ancestry is ubiquitous among people who speak Na-Dene and Eskimo-Aleut languages. We develop a comprehensive model for the Holocene peopling events of Chukotka and North America, and show that Na-Dene-speaking peoples, people of the Aleutian Islands, and Yup'ik and Inuit across the Arctic region all share ancestry from a single Palaeo-Eskimo-related Siberian source.