Human population dynamics and Yersinia pestis in ancient northeast Asia

Long shared haplotypes identify the Southern Urals as a primary source for the 10th century Hungarians

During the Hungarian Conquest in the 10th century CE, the early medieval Magyars, a group of mounted warriors from Eastern Europe, settled in the Carpathian Basin. They likely introduced the Hungarian language to this new settlement area, during an event documented by both written sources and archaeological evidence. Previous archaeogenetic research identified the newcomers as migrants from the Eurasian steppe. However, genome-wide ancient DNA from putative source populations has not been available to test alternative theories of their precise source. We generated genome-wide ancient DNA data for 131 individuals from candidate archaeological contexts in the Circum-Uralic region in present-day Russia. Our results tightly link the Magyars to people of the Early Medieval Karayakupovo archaeological horizon on both the European and Asian sides of the southern Urals. Our analyes show that ancestors of the people of the Karayakupovo archaeological horizon were established in the Southern Urals by the Iron Age and that their descendants persisted locally in the Volga-Kama region until at least the 14th century.

Investigating food production-associated DNA methylation changes in paleogenomes: Lack of consistent signals beyond technical noise

The Neolithic transition introduced major diet and lifestyle changes to human populations across continents. Beyond well-documented bioarcheological and genetic effects, whether these changes also had molecular-level epigenetic repercussions in past human populations has been an open question. In fact, methylation signatures can be inferred from UDG-treated ancient DNA through postmortem damage patterns, but with low signal-to-noise ratios; it is thus unclear whether published paleogenomes would provide the necessary resolution to discover systematic effects of lifestyle and diet shifts. To address this we compiled UDG-treated shotgun genomes of 13 pre-Neolithic hunter-gatherers (HGs) and 21 Neolithic farmers (NFs) individuals from West and North Eurasia, published by six different laboratories and with coverage c.1×-58× (median = 9×). We used epiPALEOMIX and a Monte Carlo normalization scheme to estimate methylation levels per genome. Our paleomethylome dataset showed expected genome-wide methylation patterns such as CpG island hypomethylation. However, analyzing the data using various approaches did not yield any systematic signals for subsistence type, genetic sex, or tissue effects. Comparing the HG-NF methylation differences in our dataset with methylation differences between hunter-gatherers versus farmers in modern-day Central Africa also did not yield consistent results. Meanwhile, paleomethylome profiles did cluster strongly by their laboratories of origin. Using larger data volumes, minimizing technical noise and/or using alternative protocols may be necessary for capturing subtle environment-related biological signals from paleomethylomes.

Genomic imputation of ancient Asian populations contrasts local adaptation in pre- and post-agricultural Japan

Early modern humans lived as hunter-gatherers for millennia before agriculture, yet the genetic adaptations of these populations remain a mystery. Here, we investigate selection in the ancient hunter-gatherer-fisher Jomon and contrast pre- and post-agricultural adaptation in the Japanese archipelago. Building on the successful validation of imputation with ancient Asian genomes, we identify selection signatures in the Jomon, particularly robust signals from KITLG variants, which may have influenced dark pigmentation evolution. The Jomon lacks well-known adaptive variants (EDAR, ADH1B, and ALDH2), marking their emergence after the advent of farming in the archipelago. Notably, the EDAR and ADH1B variants were prevalent in the archipelago 1,300 years ago, whereas the ALDH2 variant could have emerged later due to its absence in other ancient genomes. Overall, our study underpins local adaptation unique to the Jomon population, which in turn sheds light on post-farming selection that continues to shape contemporary Asian populations.

Bronze age Northern Eurasian genetics in the context of development of metallurgy and Siberian ancestry

The Eurasian Bronze Age (BA) has been described as a period of substantial human migrations, the emergence of pastoralism, horse domestication, and development of metallurgy. This study focuses on two north Eurasian sites sharing Siberian genetic ancestry. One of the sites, Rostovka, is associated with the Seima-Turbino (ST) phenomenon (~2200-1900 BCE) that is characterized by elaborate metallurgical objects found throughout Northern Eurasia. The genetic profiles of Rostovka individuals vary widely along the forest-tundra Siberian genetic cline represented by many modern Uralic-speaking populations, and the genetic heterogeneity observed is consistent with the current understanding of the ST being a transcultural phenomenon. Individuals from the second site, Bolshoy Oleni Ostrov in Kola, in comparison form a tighter cluster on the Siberian ancestry cline. We further explore this Siberian ancestry profile and assess the role of the ST phenomenon and other contemporaneous BA cultures in the spread of Uralic languages and Siberian ancestry.

Reconstructing the Genetic Relationship between Ancient and Present-Day Siberian Populations

Human populations across a vast area in northern Eurasia, from Fennoscandia to Chukotka, share a distinct genetic component often referred to as the Siberian ancestry. Most enriched in present-day Samoyedic-speaking populations such as Nganasans, its origins and history still remain elusive despite the growing list of ancient and present-day genomes from Siberia. Here, we reanalyze published ancient and present-day Siberian genomes focusing on the Baikal and Yakutia, resolving key questions regarding their genetic history. First, we show a long-term presence of a unique genetic profile in southern Siberia, up to 6,000 yr ago, which distinctly shares a deep ancestral connection with Native Americans. Second, we provide plausible historical models tracing genetic changes in West Baikal and Yakutia in fine resolution. Third, the Middle Neolithic individual from Yakutia, belonging to the Belkachi culture, serves as the best source so far available for the spread of the Siberian ancestry into Fennoscandia and Greenland. These findings shed light on the genetic legacy of the Siberian ancestry and provide insights into the complex interplay between different populations in northern Eurasia throughout history.

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.

Going local with ancient DNA: A review of human histories from regional perspectives

Ancient DNA (aDNA) has added a wealth of information about our species' history, including insights on genetic origins, migrations and gene flow, genetic admixture, and health and disease. Much early work has focused on continental-level questions, leaving many regional questions, especially those relevant to the Global South, comparatively underexplored. A few success stories of aDNA studies from smaller laboratories involve more local aspects of human histories and health in the Americas, Africa, Asia, and Oceania. In this Review, we cover some of these contributions by synthesizing finer-scale questions of importance to the archaeogenetics field, as well as to Indigenous and Descendant communities. We further highlight the potential of aDNA to uncover past histories in regions where colonialism has neglected the oral histories of oppressed peoples.

Mitogenomics of the Koryaks and Evens of the northern coast of the Sea of Okhotsk

Due to the geographical proximity of the northern coast of the Sea of Okhotsk and Kamchatka Peninsula to the Beringia, the indigenous populations of these territories are of great interest for elucidating the human settlement history of northern Asia and America. Meanwhile, there is a clear shortage of genetic studies of the indigenous populations of the northern coast of the Sea of Okhotsk. Here, in order to examine their fine-scale matrilineal genetic structure, ancestry and relationships with neighboring populations, we analyzed 203 complete mitogenomes (174 of which are new) from population samples of the Koryaks and Evens of the northern coast of the Sea of Okhotsk and the Chukchi of the extreme northeast Asia. The patterns observed underscore the reduced level of genetic diversity found in the Koryak, Even, and Chukchi populations, which, along with the high degree of interpopulation differentiation, may be the result of genetic drift. Our phylogeographic analysis reveals common Paleo-Asiatic ancestry for 51.1% of the Koryaks and 17.8% of the Evens. About third of the mitogenomes found in the Koryaks and Evens might be considered as ethno-specific, as these are virtually absent elsewhere in North, Central and East Asia. Coalescence ages of most of these lineages coincide well with the emergence and development of the Tokarev and Old Koryak archaeological cultures associated with the formation of the Koryaks, as well as with the period of separation and split of the North Tungusic groups migrated northwards from the Lake Baikal or the Amur River area.

Descent, marriage, and residence practices of a 3,800-year-old pastoral community in Central Eurasia

Our understanding of prehistoric societal organization at the family level is still limited. Here, we generated genome data from 32 individuals from an approximately 3,800-y-old burial mound attributed to the Bronze Age Srubnaya-Alakul cultural tradition at the site of Nepluyevsky, located in the Southern Ural region of Central Eurasia. We found that life expectancy was generally very low, with adult males living on average 8 y longer than females. A total of 35 first-degree, 40 second-degree, and 48 third-degree biological relationships connected 23 of the studied individuals, allowing us to propose a family tree spanning three generations with six brothers at its center. The oldest of these brothers had eight children with two women and the most children overall, whereas the other relationships were monogamous. Notably, related female children above the age of five were completely absent from the site, and adult females were more genetically diverse than males. These results suggest that biological relationships between male siblings played a structural role in society and that descent group membership was based on patrilineality. Women originated from a larger mating network and moved to join the men, with whom they were buried. Finally, the oldest brother likely held a higher social position, which was expressed in terms of fertility.

Modeling of African population history using f-statistics is biased when applying all previously proposed SNP ascertainment schemes

f-statistics have emerged as a first line of analysis for making inferences about demographic history from genome-wide data. Not only are they guaranteed to allow robust tests of the fits of proposed models of population history to data when analyzing full genome sequencing data-that is, all single nucleotide polymorphisms (SNPs) in the individuals being analyzed-but they are also guaranteed to allow robust tests of models for SNPs ascertained as polymorphic in a population that is an outgroup in a phylogenetic sense to all groups being analyzed. True "outgroup ascertainment" is in practice impossible in humans because our species has arisen from a substructured ancestral population that does not descend from a homogeneous ancestral population going back many hundreds of thousands of years into the past. However, initial studies suggested that non-outgroup-ascertainment schemes might produce robust enough results using f-statistics, and that motivated widespread fitting of models to data using non-outgroup-ascertained SNP panels such as the "Affymetrix Human Origins array" which has been genotyped on thousands of modern individuals from hundreds of populations, or the "1240k" in-solution enrichment reagent which has been the source of about 70% of published genome-wide data for ancient humans. In this study, we show that while analyses of population history using such panels work well for studies of relationships among non-African populations and one African outgroup, when co-modeling more than one sub-Saharan African and/or archaic human groups (Neanderthals and Denisovans), fitting of f-statistics to such SNP sets is expected to frequently lead to false rejection of true demographic histories, and failure to reject incorrect models. Analyzing panels of SNPs polymorphic in archaic humans, which has been suggested as a solution for the ascertainment problem, has limited statistical power and retains important biases. However, by carrying out simulations of diverse demographic histories, we show that bias in inferences based on f-statistics can be minimized by ascertaining on variants common in a union of diverse African groups; such ascertainment retains high statistical power while allowing co-analysis of archaic and modern groups.

The role of Beringia in human adaptation to Arctic conditions based on results of genomic studies of modern and ancient populations

The results of studies in Quaternary geology, archeology, paleoanthropology and human genetics demonstrate that the ancestors of Native Americans arrived in mid-latitude North America mainly along the Pacific Northwest Coast, but had previously inhabited the Arctic and during the last glacial maximum were in a refugium in Beringia, a land bridge connecting Eurasia and North America. The gene pool of Native Americans is represented by unique haplogroups of mitochondrial DNA and the Y chromosome, the evolutionary age of which ranges from 13 to 22 thousand years. The results of a paleogenomic analysis also show that during the last glacial maximum Beringia was populated by human groups that had arisen as a result of interaction between the most ancient Upper Paleolithic populations of Northern Eurasia and newcomer groups from East Asia. Approximately 20 thousand years ago the Beringian populations began to form, and the duration of their existence in relative isolation is estimated at about 5 thousand years. Thus, the adaptation of the Beringians to the Arctic conditions could have taken several millennia. The adaptation of Amerindian ancestors to high latitudes and cold climates is supported by genomic data showing that adaptive genetic variants in Native Americans are associated with various metabolic pathways: melanin production processes in the skin, hair and eyes, the functioning of the cardiovascular system, energy metabolism and immune response characteristics. Meanwhile, the analysis of the existing hypotheses about the selection of some genetic variants in the Beringian ancestors of the Amerindians in connection with adaptation to the Arctic conditions (for example, in the FADS, ACTN3, EDAR genes) shows the ambiguity of the testing results, which may be due to the loss of some traces of the "Beringian" adaptation in the gene pools of modern Native Americans. The most optimal strategy for further research seems to be the search for adaptive variant.

Genetic population structure of the Xiongnu Empire at imperial and local scales

The Xiongnu established the first nomadic imperial power, controlling the Eastern Eurasian steppe from ca. 200 BCE to 100 CE. Recent archaeogenetic studies identified extreme levels of genetic diversity across the empire, corroborating historical records of the Xiongnu Empire being multiethnic. However, it has remained unknown how this diversity was structured at the local community level or by sociopolitical status. To address this, we investigated aristocratic and local elite cemeteries at the western frontier of the empire. Analyzing genome-wide data from 18 individuals, we show that genetic diversity within these communities was comparable to the empire as a whole, and that high diversity was also observed within extended families. Genetic heterogeneity was highest among the lowest-status individuals, implying diverse origins, while higher-status individuals harbored less genetic diversity, suggesting that elite status and power was concentrated within specific subsets of the broader Xiongnu population.

The impacts of bronze age in the gene pool of Chinese: Insights from phylogeographics of Y-chromosomal haplogroup N1a2a-F1101

Objectives: Previous studies of archaeology and history suggested that the rise and prosperity of Bronze Age culture in East Asia had made essential contribution to the formation of early state and civilization in this region. However, the impacts in perspective of genetics remain ambiguous. Previous genetic researches indicated the Y-chromosome Q1a1a-M120 and N1a2a-F1101 may be the two most important paternal lineages among the Bronze Age people in ancient northwest China. Here, we investigated the 9,000-years history of haplogroup N1a2a-F1101 with revised phylogenetic tree and spatial autocorrelation analysis.

Materials and Methods: In this study, 229 sequences of N1a2a-F1101 were analyzed. We developed a highly-revised phylogenetic tree with age estimates for N1a2a-F1101. In addition, we also explored the geographical distribution of sub-lineages of N1a2a-F1101, and spatial autocorrelation analysis was conducted for each sub-branch.

Results: The initial differentiation location of N1a2a-F1101 and its most closely related branch, N1a2b-P43, a major lineage of Uralic-speaking populations in northern Eurasia, is likely the west part of northeast China. After ~4 thousand years of bottleneck effect period, haplgroup N1a2a-F1101 experienced continuous expansion during the Chalcolithic age (~ 4.5 kya to 4 kya) and Bronze age (~ 4 kya to 2.5 kya) in northern China. Ancient DNA evidence supported that this haplogroup is the lineage of ruling family of Zhou Dynasty (~ 3 kya-2.2 kya) of ancient China.

Discussion: In general, we proposed that the Bronze Age people in the border area between the eastern Eurasian steppe and northern China not only played a key role in promoting the early state and civilization of China, but also left significant traces in the gene pool of Chinese people.

Middle Holocene Siberian genomes reveal highly connected gene pools throughout North Asia

The peopling history of North Asia remains largely unexplored due to the limited number of ancient genomes analyzed from this region. Here, we report genome-wide data of ten individuals dated to as early as 7,500 years before present from three regions in North Asia, namely Altai-Sayan, Russian Far East, and the Kamchatka Peninsula. Our analysis reveals a previously undescribed Middle Holocene Siberian gene pool in Neolithic Altai-Sayan hunter-gatherers as a genetic mixture between paleo-Siberian and ancient North Eurasian (ANE) ancestries. This distinctive gene pool represents an optimal source for the inferred ANE-related population that contributed to Bronze Age groups from North and Inner Asia, such as Lake Baikal hunter-gatherers, Okunevo-associated pastoralists, and possibly Tarim Basin populations. We find the presence of ancient Northeast Asian (ANA) ancestry-initially described in Neolithic groups from the Russian Far East-in another Neolithic Altai-Sayan individual associated with different cultural features, revealing the spread of ANA ancestry ∼1,500 km further to the west than previously observed. In the Russian Far East, we identify 7,000-year-old individuals that carry Jomon-associated ancestry indicating genetic links with hunter-gatherers in the Japanese archipelago. We also report multiple phases of Native American-related gene flow into northeastern Asia over the past 5,000 years, reaching the Kamchatka Peninsula and central Siberia. Our findings highlight largely interconnected population dynamics throughout North Asia from the Early Holocene onward.

Modeling of African population history using f -statistics can be highly biased and is not addressed by previously suggested SNP ascertainment schemes

f -statistics have emerged as a first line of analysis for making inferences about demographic history from genome-wide data. These statistics can provide strong evidence for either admixture or cladality, which can be robust to substantial rates of errors or missing data. f -statistics are guaranteed to be unbiased under "SNP ascertainment" (analyzing non-randomly chosen subsets of single nucleotide polymorphisms) only if it relies on a population that is an outgroup for all groups analyzed. However, ascertainment on a true outgroup that is not co-analyzed with other populations is often impractical and uncommon in the literature. In this study focused on practical rather than theoretical aspects of SNP ascertainment, we show that many non-outgroup ascertainment schemes lead to false rejection of true demographic histories, as well as to failure to reject incorrect models. But the bias introduced by common ascertainments such as the 1240K panel is mostly limited to situations when more than one sub-Saharan African and/or archaic human groups (Neanderthals and Denisovans) or non-human outgroups are co-modelled, for example, f 4 -statistics involving one non-African group, two African groups, and one archaic group. Analyzing panels of SNPs polymorphic in archaic humans, which has been suggested as a solution for the ascertainment problem, cannot fix all these problems since for some classes of f -statistics it is not a clean outgroup ascertainment, and in other cases it demonstrates relatively low power to reject incorrect demographic models since it provides a relatively small number of variants common in anatomically modern humans. And due to the paucity of high-coverage archaic genomes, archaic individuals used for ascertainment often act as sole representatives of the respective groups in an analysis, and we show that this approach is highly problematic. By carrying out large numbers of simulations of diverse demographic histories, we find that bias in inferences based on f -statistics introduced by non-outgroup ascertainment can be minimized if the derived allele frequency spectrum in the population used for ascertainment approaches the spectrum that existed at the root of all groups being co-analyzed. Ascertaining on sites with variants common in a diverse group of African individuals provides a good approximation to such a set of SNPs, addressing the great majority of biases and also retaining high statistical power for studying population history. Such a "pan-African" ascertainment, although not completely problem-free, allows unbiased exploration of demographic models for the widest set of archaic and modern human populations, as compared to the other ascertainment schemes we explored.

Beringia and the peopling of the Western Hemisphere

Did Beringian environments represent an ecological barrier to humans until less than 15 000 years ago or was access to the Americas controlled by the spatial-temporal distribution of North American ice sheets? Beringian environments varied with respect to climate and biota, especially in the two major areas of exposed continental shelf. The East Siberian Arctic Shelf ('Great Arctic Plain' (GAP)) supported a dry steppe-tundra biome inhabited by a diverse large-mammal community, while the southern Bering-Chukchi Platform ('Bering Land Bridge' (BLB)) supported mesic tundra and probably a lower large-mammal biomass. A human population with west Eurasian roots occupied the GAP before the Last Glacial Maximum (LGM) and may have accessed mid-latitude North America via an interior ice-free corridor. Re-opening of the corridor less than 14 000 years ago indicates that the primary ancestors of living First Peoples, who already had spread widely in the Americas at this time, probably dispersed from the NW Pacific coast. A genetic 'arctic signal' in non-arctic First Peoples suggests that their parent population inhabited the GAP during the LGM, before their split from the former. We infer a shift from GAP terrestrial to a subarctic maritime economy on the southern BLB coast before dispersal in the Americas from the NW Pacific coast.

Genetic admixture and language shift in the medieval Volga-Oka interfluve

The Volga-Oka interfluve in northwestern Russia has an intriguing history of population influx and language shift during the Common Era. Today, most inhabitants of the region speak Russian, but until medieval times, northwestern Russia was inhabited by Uralic-speaking peoples.^1,2,3 A gradual shift to Slavic languages started in the second half of the first millennium with the expansion of Slavic tribes, which led to the foundation of the Kievan Rus' state in the late 9^th century CE. The medieval Rus' was multicultural and multilingual-historical records suggest that its northern regions comprised Slavic and Uralic peoples ruled by Scandinavian settlers.^4,5,6 In the 10^th-11^th centuries, the introduction of Christianity and Cyrillic literature raised the prestige status of Slavic, driving a language shift from Uralic to Slavic.³ This eventually led to the disappearance of the Uralic languages from northwestern Russia. Here, we study a 1,500-year time transect of 30 ancient genomes and stable isotope values from the Suzdal region in the Volga-Oka interfluve. We describe a previously unsampled local Iron Age population and a gradual genetic turnover in the following centuries. Our time transect captures the population shift associated with the spread of Slavic languages and illustrates the ethnically mixed state of medieval Suzdal principality, eventually leading to the formation of the admixed but fully Slavic-speaking population that inhabits the area today. We also observe genetic outliers that highlight the importance of the Suzdal region in medieval times as a hub of long-reaching contacts via trade and warfare.

Genomic Analysis of 18th-Century Kazakh Individuals and Their Oral Microbiome

The Asian Central Steppe, consisting of current-day Kazakhstan and Russia, has acted as a highway for major migrations throughout history. Therefore, describing the genetic composition of past populations in Central Asia holds value to understanding human mobility in this pivotal region. In this study, we analyse paleogenomic data generated from five humans from Kuygenzhar, Kazakhstan. These individuals date to the early to mid-18th century, shortly after the Kazakh Khanate was founded, a union of nomadic tribes of Mongol Golden Horde and Turkic origins. Genomic analysis identifies that these individuals are admixed with varying proportions of East Asian ancestry, indicating a recent admixture event from East Asia. The high amounts of DNA from the anaerobic Gram-negative bacteria Tannerella forsythia, a periodontal pathogen, recovered from their teeth suggest they may have suffered from periodontitis disease. Genomic analysis of this bacterium identified recently evolved virulence and glycosylation genes including the presence of antibiotic resistance genes predating the antibiotic era. This study provides an integrated analysis of individuals with a diet mostly based on meat (mainly horse and lamb), milk, and dairy products and their oral microbiome.

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.

The deep population history of northern East Asia from the Late Pleistocene to the Holocene

Northern East Asia was inhabited by modern humans as early as 40 thousand years ago (ka), as demonstrated by the Tianyuan individual. Using genome-wide data obtained from 25 individuals dated to 33.6-3.4 ka from the Amur region, we show that Tianyuan-related ancestry was widespread in northern East Asia before the Last Glacial Maximum (LGM). At the close of the LGM stadial, the earliest northern East Asian appeared in the Amur region, and this population is basal to ancient northern East Asians. Human populations in the Amur region have maintained genetic continuity from 14 ka, and these early inhabitants represent the closest East Asian source known for Ancient Paleo-Siberians. We also observed that EDAR V370A was likely to have been elevated to high frequency after the LGM, suggesting the possible timing for its selection. This study provides a deep look into the population dynamics of northern East Asia.