Bronze and Iron Age population movements underlie Xinjiang population history

Urbanization and genetic homogenization in the medieval Low Countries revealed through a ten-century paleogenomic study of the city of Sint-Truiden

BACKGROUND

Processes shaping the formation of the present-day population structure in highly urbanized Northern Europe are still poorly understood. Gaps remain in our understanding of when and how currently observable regional differences emerged and what impact city growth, migration, and disease pandemics during and after the Middle Ages had on these processes.

RESULTS

We perform low-coverage sequencing of the genomes of 338 individuals spanning the eighth to the eighteenth centuries in the city of Sint-Truiden in Flanders, in the northern part of Belgium. The early/high medieval Sint-Truiden population was more heterogeneous, having received migrants from Scotland or Ireland, and displayed less genetic relatedness than observed today between individuals in present-day Flanders. We find differences in gene variants associated with high vitamin D blood levels between individuals with Gaulish or Germanic ancestry. Although we find evidence of a Yersinia pestis infection in 5 of the 58 late medieval burials, we were unable to detect a major population-scale impact of the second plague pandemic on genetic diversity or on the elevated differentiation of immunity genes.

CONCLUSIONS

This study reveals that the genetic homogenization process in a medieval city population in the Low Countries was protracted for centuries. Over time, the Sint-Truiden population became more similar to the current population of the surrounding Limburg province, likely as a result of reduced long-distance migration after the high medieval period, and the continuous process of local admixture of Germanic and Gaulish ancestries which formed the genetic cline observable today in the Low Countries.

Ancient DNA indicates 3,000 years of genetic continuity in the Northern Iranian Plateau, from the Copper Age to the Sassanid Empire

In this study, we present new ancient DNA data from prehistoric and historic populations of the Iranian Plateau. By analysing 50 samples from nine archaeological sites across Iran, we report 23 newly sequenced mitogenomes and 13 nuclear genomes, spanning 4700 BCE to 1300 CE. We integrate an extensive reference sample set of previously published ancient DNA datasets from Western and South-Central Asia, enhancing our understanding of genetic continuity and diversity within ancient Iranian populations. A new Early Chalcolithic sample, predating all other Chalcolithic genomes from Iran, demonstrates mostly Early Neolithic Iranian genetic ancestry. This finding reflects long-term cultural and biological continuity in and around the Zagros area, alongside evidence of some western genetic influence. Our sample selection prioritizes northern Iran, with a particular focus on the Achaemenid, Parthian, and Sassanid periods (355 BCE-460 CE). The genetic profiles of historical samples from this region position them as intermediates on an east-west genetic cline across the Persian Plateau. They also exhibit strong connections to local and South-Central Asian Bronze Age populations, underscoring enduring genetic connections across these regions. Diachronic analyses of uniparental lineages on the Iranian Plateau further highlight population stability from prehistoric to modern times.

Demographics of genetic admixture and expansion

We explore the dynamics of genetic admixture and expansion, as well as language assimilation, through mathematical-demographic modeling. Our primary goal is to address the population-genetics 'paradox' wherein autosomes and allosomes present markedly different, if not contradictory, pictures of past migrations. We demonstrate that this paradox may find a purely demographic explanation, as single-sex and two-sex reproduction models exhibit markedly distinct dynamics. We illustrate that the three processes-allosomal expansion, autosomal admixture, and language assimilation-occur at significantly different modes, potentially explaining the varied outcomes of these processes upon the completion of ethnogenetic transitions. Our research offers valuable insights into the intricate interplay between demography, genetics, and social organization, providing implications for historical scenarios and enhancing our understanding of the long-term consequences of migration and social cohesion.

Ancient genomes shed light on the genetic history of the Iron Age to historical central Xinjiang, northwest China

BACKGROUND

The genetic profile of the population in Xinjiang, northwest China, has been shaped by interregional movement and admixture since the Bronze Age. However, the detailed and intraregional population history of Xinjiang, especially central Xinjiang, has been unsolved due to uneven sample distribution.

RESULTS

Here, we reported the ancient genomes from 8 individuals between the Iron Age and the historical period in central Xinjiang. We observed an east-west admixed ancestry profile and a degree of genetic continuity between the Iron Age and historical central Xinjiang individuals. Furthermore, these central Xinjiang individuals harboured ancestry related to ancient farmers of the Yellow River. We also identified a temporal change of the Yellow River farmers-related ancestry in central Xinjiang, showing an increase the Yelllow River affinity from Iron Age to Historical Era.

CONCLUSIONS

The finding indicated that the genetic structure of the central Xinjiang population since the Iron Age could have resulted from immigration from northern China, which was attributed to geopolitical factors. Hence, our results indicated that the geopolitical change with the deepening of Central Plains' management has influenced the genetic profile of central Xinjiang.

Population genomics advances in frontier ethnic minorities in China

China, with its large geographic span, possesses rich genetic diversity across vast frontier regions in addition to the Han Chinese majority. Importantly, demographic events and various natural and cultural environments in Chinese frontier regions have shaped the genomic diversity of ethnic minorities via local adaptations. Thus, insights into the genetic diversity and adaptive evolution of these under-represented ethnic groups are crucial for understanding evolutionary scenarios and biomedical implications in East Asian populations. Here, we focus on ethnic minorities in Chinese frontier regions and review research advances regarding genomic diversity, genetic structure, population history, genetic admixture, and local adaptation. We first provide an overview of the extensive genetic diversity across populations in different Chinese frontier regions. Next, we summarize research progress regarding genetic ancestry, demographic history, the adaptive process, and the archaic identification of multiple ethnic minorities in different Chinese frontier regions. Finally, we discuss the gaps and opportunities in genomic studies of Chinese populations and the need for a more comprehensive understanding of genomic diversity and the evolution of populations of East Asian ancestry in the post-genomic era.

Pilot work of the 10K Chinese People Genomic Diversity Project along the Silk Road suggests a complex east-west admixture landscape and biological adaptations

Genomic sources from China are underrepresented in the population-specific reference database. We performed whole-genome sequencing or genome-wide genotyping on 1,207 individuals from four linguistically diverse groups (1,081 Sinitic, 56 Mongolic, 40 Turkic, and 30 Tibeto-Burman people) living in North China included in the 10K Chinese People Genomic Diversity Project (10K_CPGDP) to characterize the genetic architecture and adaptative history of ethnic groups in the Silk Road Region of China. We observed a population split between Northwest Chinese minorities (NWCMs) and Han Chinese since the Upper Paleolithic and later Neolithic genetic differentiation within NWCMs. The observed population substructures among ethnically/linguistically diverse NWCMs suggested that differentiated admixture events contributed to the differences in their genomic and phenotypic diversity. We estimated that the Dongxiang, Tibetan, and Yugur people inherited more than 10% of the Western Eurasian ancestry, which is much greater than that of the Salar and Tu people (<7%), while Han neighbors showed less West Eurasian ancestry (∼1%-3%). Male-biased admixture introduced Western Eurasian ancestry in the Dongxiang, Tibetan, and Yugur populations. We found that the eastern-western admixture in NWCMs occurred ∼800-1,100 years ago, coinciding with intensive economic and cultural exchanges during the Tang and Song dynasties. Additionally, we identified the signatures of natural selection associated with cardiovascular system diseases or lipid metabolism and developmental/neurogenetic disorders. Moreover, the EPAS1 gene showed relatively high population branch statistic values in NWCMs. The well-fitted demographical models presented the vast landscape of complex admixture processes of the Silk Road people, and the newly reported functionally important variations suggested the importance of including ethnolinguistically diverse populations in Chinese genetic studies for uncovering the genetic basis of complex traits/diseases.

Ancient DNA analysis of elite nomadic warrior from Chinge-Tey I funerary commemorative complex in the "Valley of the Kings", Tuva

BACKGROUND

In the Ist millennium BC bearers of the Scythian-type nomadic cultures inhabited the steppes of Eurasia, from Northern China to the Carpathians. According to archaeological data, the origin of nomadic life style and economy can be traced to the eastern part of this steppe "corridor", primarily to the territory of the present-day Republic of Tuva in Russia. Here, in the Turan-Uyuk Basin, also known as the "Valley of the Kings", some of the earliest known Scythian-type archaeological sites called Arzhan-1, Arzhan-2, Chinge-Tey I, Tunnug 1 were studied. Each of them is a large-scale funerary commemorative complex with burials of tribal nomadic leaders, surrounded by graves of supposed members of their families or associates. All these people belonged to the societies which are associated with the earliest nomadic cultures in Asia. Representatives of similar cultures will later be known and described as the Scythians/the Saka in Assyrian, Achaemenid, and Greek sources. Arzhan 2 and Chinge-Tey I elite level sites as well as ordinary pastoralist burials of the early-Scythian period in Tuva are attributed to the Aldy-Bel archaeological culture of the Early Iron Age (8th- 6th century BC). Taking the first step to shed light on the genetic origin of Aldy-Bel elites, we carried out a comparative genome-wide analysis of an elite level person buried in grave 9 at Chinge-Tey I (7th- 6th centuries BC) and two published earlier genomes of individuals, whose burials (graves 14 and 22) accompanied the 'royal couple' (grave 5) at Arzhan-2. This study aims also at checking a hypothesis of genetic kinship between human individuals buried in the large-scale burial complexes of the "Valley of the Kings" and brings up the issue of possible dynastic connections of local elites, buried under different kurgans of the valley.

RESULTS

First, ancient DNA analysis of an elite nomadic warrior from Chinge-Tey I has been carried out, thus a third wide-genome dataset for Aldy-Bel culture- one of the earliest nomadic cultures in Asia, is presented in this study. Second, we undertook a comparative analysis of genome-wide data of three mentioned Aldy-Bel culture representatives and individuals of the other Bronze and Early Iron Age population groups of Asia to estimate their possible genetic connections. Then, kinship analysis was undertaken for these three Aldy-Bel culture individuals. Finally, mitochondrial and Y-chromosome haplogroups of Chinge-Tey princely person were compared to those of other Aldy-Bel culture representatives and to individuals of subsequent Scythian-type Uyuk-Sagly culture in Tuva.

CONCLUSION

(1) Generating the third wide-genome of the enabled us to undertake its comparison with two other genomes of Aldy-Bel culture representatives (Arzhan-2, graves 14 and 22) and with other Bronze and Early Iron Age population groups in Asia to trace the origin and genetic connection of Aldy-Bel population, representing one of the earliest Scythian-type nomadic group. (2) The results obtained show that the princely individual from Chinge-Tey I and two 'king's associates' from Arzhan-2 were genetically close to nomads of simultaneous Tasmola culture in Eastern and Central Kazakhstan and pastoralists buried in the Early Iron Age cemeteries of present-day Xinjiang (first of all, Abusanteer archaeological site). Aldy-Bel culture representatives appeared also close to individuals of the Middle Bronze Age Okunevo culture in the Minusinsk Basin. Besides, Aldy-Bel pastoralists turned out genetically close to nomads of the subsequent Uyuk-Sagly culture in Mongolia (5th - 3rd centuries BC). (3) Ancient DNA kinship analyses, undertaken for three Aldy-Bel culture individuals pointed out to the absence of their tribe kinship. (4) On the other hand, Chinge-Tey warrior's mitochondrial haplogroup G was previously described in two (graves 14 and 5) individuals from Arzhan-2, including a female individual from the "royal" tomb 5. This result provided a possibility of maternal kinship among this so called 'queen' from Arzhan-2 and the princely person from Chinge-Tey I. This possibility supported a hypothesis of their family ties suggested on archaeological materials. Y-chromosome haplogroup Q1b1, revealed for the princely person, was widely distributed among local people of Aldy-Bel and subsequent Uyuk-Sagly cultures.

Ancient genomes reveal trans-Eurasian connections between the European Huns and the Xiongnu Empire

The Huns appeared in Europe in the 370s, establishing an Empire that reshaped West Eurasian history. Yet until today their origins remain a matter of extensive debate. Traditional theories link them to the Xiongnu, the founders of the first nomadic empire of the Mongolian steppe. The Xiongnu empire dissolved, however, ~300 y before the Huns appeared in Europe, and there is little archaeological and historical evidence of Huns in the steppe during this time gap. Furthermore, despite the rich 5th to 6th centuries current era (CE) archaeological record of the Carpathian Basin, the cultural elements of connections with the steppe are limited to few findings and even fewer solitary eastern-type burials. In this study, we coanalyze archaeological evidence with 35 newly sequenced and published genomic data for a total of 370 individuals-from 5th to 6th century CE contexts in the Carpathian Basin including 10 Hun-period eastern-type burials, 2nd to 5th century sites across Central Asia and 2nd c. before current era (BCE) to 1st c. CE Xiongnu period sites across the Mongolian steppe. We find no evidence for the presence of a large eastern/steppe descent community among the Hun- and post-Hun-period Carpathian Basin population. We also observe a high genetic diversity among the eastern-type burials that recapitulates the variability observed across the Eurasian Steppe. This suggests a mixed origin of the incoming steppe conquerors. Nevertheless, long-shared genomic tracts provide compelling evidence of genetic lineages directly connecting some individuals of the highest Xiongnu-period elite with 5th to 6th century CE Carpathian Basin individuals, showing that some European Huns descended from them.

The genetic origin of the Indo-Europeans

The Yamnaya archaeological complex appeared around 3300 BC across the steppes north of the Black and Caspian Seas, and by 3000 BC it reached its maximal extent, ranging from Hungary in the west to Kazakhstan in the east. To localize Yamnaya origins among the preceding Eneolithic people, we assembled ancient DNA from 435 individuals, demonstrating three genetic clines. A Caucasus-lower Volga (CLV) cline suffused with Caucasus hunter-gatherer1 ancestry extended between a Caucasus Neolithic southern end and a northern end at Berezhnovka along the lower Volga river. Bidirectional gene flow created intermediate populations, such as the north Caucasus Maikop people, and those at Remontnoye on the steppe. The Volga cline was formed as CLV people mixed with upriver populations of Eastern hunter-gatherer2 ancestry, creating hypervariable groups, including one at Khvalynsk. The Dnipro cline was formed when CLV people moved west, mixing with people with Ukraine Neolithic hunter-gatherer ancestry3 along the Dnipro and Don rivers to establish Serednii Stih groups, from whom Yamnaya ancestors formed around 4000 BC and grew rapidly after 3750-3350 BC. The CLV people contributed around four-fifths of the ancestry of the Yamnaya and, entering Anatolia, probably from the east, at least one-tenth of the ancestry of Bronze Age central Anatolians, who spoke Hittite4,5. We therefore propose that the final unity of the speakers of 'proto-Indo-Anatolian', the language ancestral to both Anatolian and Indo-European people, occurred in CLV people some time between 4400 BC and 4000 BC.

East and West admixture in eastern China of Tang Dynasty inferred from ancient human genomes

The ancestry composition and Sinicisation process of the descendants of the immigrants of Hu people living in ancient China are largely unknown due to the lack of genetic evidence. Tang Dynasty people in Fudamen cemetery () excavated from Shandong province in eastern China are believed to be related to the descendants of Hu people, as some of the individuals with the surnames An (). The genetic origin of the Fudamen population requires genetic clarification using ancient DNA data. Here we successfully obtain genome-wide SNP data for 17 Tang dynasty individuals from Fudamen cemetery. Based on autosomal data, although all Fudamen individuals show high levels of middle reaches of Yellow River-related ancestry as previously published historical period Shandong populations, 2 Fudamen individuals require ~5% Western Eurasian/Central Asian-related ancestry to describe their ancestry composition best. To the best of our knowledge, it is the first evidence of such ancestry in Shandong, the eastern part of today's China. Moreover, the admixture pattern is also reflected by the presence of both west and east Eurasian-specific mtDNA and Y chromosomal haplogroups in Fudamen people. The estimated admixture time is also consistent with periods when Sogdians and other non-Han populations were active in ancient China. These genomic findings suggest that intermarriage with Han Chinese involved the Sinicization process of the Hu people.

Directional Selection and Evolution of Polygenic Traits in Eastern Eurasia: Insights from Ancient DNA

This study explores directional selection on physical and psychosocial phenotypes in Eastern Eurasian populations, utilizing a dataset of 1245 ancient genomes. By analyzing polygenic scores (PGS) for traits including height, educational attainment (EA), IQ, autism, schizophrenia, and others, we observed significant temporal trends spanning the Holocene era. The results suggest positive selection for cognitive-related traits such as IQ, EA and autism spectrum disorder (ASD), alongside negative selection for anxiety and depression. The results for height were mixed and showed nonlinear relationships with Years Before Present (BP). These trends were partially mediated by genetic components linked to distinct ancestral populations. Regression models incorporating admixture, geography, and temporal variables were used to account for biases in population composition over time. Latitude showed a positive effect on ASD PGS, EA and height, while it had a negative effect on skin pigmentation scores. Additionally, latitude exhibited significant nonlinear effects on multiple phenotypes. The observed patterns highlight the influence of climate-mediated selection pressures on trait evolution. Spline regression revealed that several polygenic scores had nonlinear relationships with years BP. The findings provide evidence for complex evolutionary dynamics, with distinct selective pressures shaping phenotypic diversity across different timescales and environments.

YHSeqY3000 panel captures all founding lineages in the Chinese paternal genomic diversity database

BACKGROUND

The advancements in second-/third-generation sequencing technologies, alongside computational innovations, have significantly enhanced our understanding of the genomic structure of Y-chromosomes and their unique phylogenetic characteristics. These researches, despite the challenges posed by the lack of population-scale genomic databases, have the potential to revolutionize our approach to high-resolution, population-specific Y-chromosome panels and databases for anthropological and forensic applications.

OBJECTIVES

This study aimed to develop the highest-resolution Y-targeted sequencing panel, utilizing time-stamped, core phylogenetic informative mutations identified from high-coverage sequences in the YanHuang cohort. This panel is intended to provide a new tool for forensic complex pedigree search and paternal biogeographical ancestry inference, as well as explore the general patterns of the fine-scale paternal evolutionary history of ethnolinguistically diverse Chinese populations.

RESULTS

The sequencing performance of the East Asian-specific Y-chromosomal panel, including 2999-core SNP variants, was found to be robust and reliable. The YHSeqY3000 panel was designed to capture the genetic diversity of Chinese paternal lineages from 3500 years ago, identifying 408 terminal lineages in 2097 individuals across 41 genetically and geographically distinct populations. We identified a fine-scale paternal substructure that was correlating with ancient population migrations and expansions. New evidence was provided for extensive gene flow events between minority ethnic groups and Han Chinese people, based on the integrative Chinese Paternal Genomic Diversity Database.

CONCLUSIONS

This work successfully integrated Y-chromosome-related basic genomic science with forensic and anthropological translational applications, emphasizing the necessity of comprehensively characterizing Y-chromosome genomic diversity from genomically under-representative populations. This is particularly important in the second phase of our population-specific medical or anthropological genomic cohorts, where dense sampling strategies are employed.

Traces of Bronze Age globalization in East Asia: Insights from a revised phylogeography of the Y-chromosome haplogroup Q1a1a-M120

OBJECTIVE

In this study, we aim to explore the genetic imprint of Bronze Age globalization in East Asia from a phylogeographic perspective by examining the Y-chromosome haplogroup Q1a1a-M120, and to identify key demographic processes involved in the formation of early China and the ancient Huaxia people.

METHODS

Over the past few decades, we have collected the sequences of 347 Y chromosomes from the haplogroup Q1a1a-M120. These sequences were utilized to analyze and reconstruct a highly revised phylogenetic tree with age estimates. And we analyzed the geographical distribution and spatial autocorrelation of nine major sub-branches of Q1a1a-M120. Finally, we observed the expansion of Q1a1a-M120 from the beginning of the Bronze Age in East Asia, along with the continuous dissemination of its sub-lineages among East Asian populations.

RESULTS

We suggest that certain sub-lineages played a significant role in the formation of states and early civilizations in China, as well as in the development of the ancient Huaxia people, who are the direct ancestors of the Han population. Overall, we propose that haplogroup Q-M120 played a role in the introduction of Bronze Age culture to the central region of East Asia. Therefore, it is haplogroup Q-M120, rather than the Western Eurasian paternal lineage, that expanded and contributed to the gene pool of the East Asian population.

CONCLUSION

In summary, the globalization of the Bronze Age led to large-scale population replacement and admixture across various regions of Eurasia; our findings highlight the unique demographic processes that occurred in East Asia during this period.

Genetic Polymorphism of Y-Chromosome in Turkmen Population from Turkmenistan

This study investigates the Y-chromosome genetic diversity of the Turkmen population in Turkmenistan, analyzing 23 Y-STR loci for the first time in a sample of 100 individuals. Combined with comparative data from Turkmen populations in Afghanistan, Iran, Iraq, Russia, and Uzbekistan, this analysis offers insights into the genetic structure and relationships among Turkmen populations across regions across Central Asia and the Near East. High haplotype diversity in the Turkmen of Turkmenistan is shaped by founder effects (lineage expansions) from distinct haplogroups, with haplogroups Q and R1a predominating. Subhaplogroups Q1a and Q1b identified in Turkmenistan trace back to ancient Y-chromosome lineages from the Bronze Age. Comparative analyses, including genetic distance (RST), median-joining network, and multidimensional scaling (MDS), highlight the genetic proximity of the Turkmen in Turkmenistan to those in Afghanistan and Iran, while Iraqi Turkmen display unique characteristics, aligning with Near Eastern populations. This study underscores the Central Asian genetic affinity across most Turkmen populations. It demonstrates the value of deep-sequencing Y-chromosome data in tracing the patrilineal history of Central Asia for future studies. These findings contribute to a more comprehensive understanding of Turkmen genetic ancestry and add new data to the ongoing study of Central Asian population genetics.

Ancient genomes from the Tang Dynasty capital reveal the genetic legacy of trans-Eurasian communication at the eastern end of Silk Road

BACKGROUND

Ancient Chang'an in the Tang Dynasty (618-907 AD) was one of the world's largest and most populated cities and acted as the eastern end of the world-famous Silk Road. However, little is known about the genetics of Chang'an people and whether the Western Regions-related gene flows have been prevalent in this cosmopolitan city.

RESULTS

Here, we present seven genomes from Xingfulindai (XFLD) sites dating to the Tang Dynasty in Chang'an. We observed that four of seven XFLD individuals (XFLD_1) were genetically homogenous with the Late Neolithic Wadian, Pingliangtai, and Haojiatai populations from the middle reaches of the Yellow River Basin (YR_LN), with no genetic influence from the Western Eurasian or other non-Yellow River-related lineages. The remaining three XFLD individuals were a mixture of YR_LN-related ancestry and ~ 3-15% Western Eurasian-related ancestry. Mixtures of XFLD_1 and Western Eurasian-related ancestry drove the main gradient of genetic variation in northern and central Shaanxi Province today.

CONCLUSIONS

Our study underlined the widespread distribution of the YR_LN-related ancestry alongside the Silk Road within the territory of China during the historical era and provided direct evidence of trans-Eurasian communication in Chang'an from a genetic perspective.

Bronze Age cheese reveals human-Lactobacillus interactions over evolutionary history

Despite the long history of consumption of fermented dairy, little is known about how the fermented microbes were utilized and evolved over human history. Here, by retrieving ancient DNA of Bronze Age kefir cheese (∼3,500 years ago) from the Xiaohe cemetery, we explored past human-microbial interactions. Although it was previously suggested that kefir was spread from the Northern Caucasus to Europe and other regions, we found an additional spreading route of kefir from Xinjiang to inland East Asia. Over evolutionary history, the East Asian strains gained multiple gene clusters with defensive roles against environmental stressors, which can be a result of the adaptation of Lactobacillus strains to various environmental niches and human selection. Overall, our results highlight the role of past human activities in shaping the evolution of human-related microbes, and such insights can, in turn, provide a better understanding of past human behaviors.

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.

Improving data archiving practices in ancient genomics

Ancient DNA is producing a rich record of past genetic diversity in humans and other species. However, unless the primary data is appropriately archived, its long-term value will not be fully realised. I surveyed publicly archived data from 42 recent ancient genomics studies. Half of the studies archived incomplete datasets, preventing accurate replication and representing a loss of data of potential future use. No studies met all criteria that could be considered best practice. Based on these results, I make six recommendations for data producers: (1) archive all sequencing reads, not just those that aligned to a reference genome, (2) archive read alignments too, but as secondary analysis files, (3) provide correct experiment metadata on samples, libraries and sequencing runs, (4) provide informative sample metadata, (5) archive data from low-coverage and negative experiments, and (6) document archiving choices in papers, and peer review these. Given the reliance on destructive sampling of finite material, ancient genomics studies have a particularly strong responsibility to ensure the longevity and reusability of generated data.

Multiple Human Population Movements and Cultural Dispersal Events Shaped the Landscape of Chinese Paternal Heritage

Large-scale genomic projects and ancient DNA innovations have ushered in a new paradigm for exploring human evolutionary history. However, the genetic legacy of spatiotemporally diverse ancient Eurasians within Chinese paternal lineages remains unresolved. Here, we report an integrated Y-chromosome genomic database encompassing 15,563 individuals from both modern and ancient Eurasians, including 919 newly reported individuals, to investigate the Chinese paternal genomic diversity. The high-resolution, time-stamped phylogeny reveals multiple diversification events and extensive expansions in the early and middle Neolithic. We identify four major ancient population movements, each associated with technological innovations that have shaped the Chinese paternal landscape. First, the expansion of early East Asians and millet farmers from the Yellow River Basin predominantly carrying O2/D subclades significantly influenced the formation of the Sino-Tibetan people and facilitated the permanent settlement of the Tibetan Plateau. Second, the dispersal of rice farmers from the Yangtze River Valley carrying O1 and certain O2 sublineages reshapes the genetic makeup of southern Han Chinese, as well as the Tai-Kadai, Austronesian, Hmong-Mien, and Austroasiatic people. Third, the Neolithic Siberian Q/C paternal lineages originated and proliferated among hunter-gatherers on the Mongolian Plateau and the Amur River Basin, leaving a significant imprint on the gene pools of northern China. Fourth, the J/G/R paternal lineages derived from western Eurasia, which were initially spread by Yamnaya-related steppe pastoralists, maintain their presence primarily in northwestern China. Overall, our research provides comprehensive genetic evidence elucidating the significant impact of interactions with culturally distinct ancient Eurasians on the patterns of paternal diversity in modern Chinese populations.

Ancient genomes revealed the complex human interactions of the ancient western Tibetans

The western Tibetan Plateau is the crossroad between the Tibetan Plateau, Central Asia, and South Asia, and it is a potential human migration pathway connecting these regions. However, the population history of the western Tibetan Plateau remains largely unexplored due to the lack of ancient genomes covering a long-time interval from this area. Here, we reported genome-wide data of 65 individuals dated to 3,500-300 years before present (BP) in the Ngari prefecture. The ancient western Tibetan Plateau populations share the majority of their genetic components with the southern Tibetan Plateau populations and have maintained genetic continuity since 3,500 BP while maintaining interactions with populations within and outside the Tibetan Plateau. Within the Tibetan Plateau, the ancient western Tibetan Plateau populations were influenced by the additional expansion from the south to the southwest plateau before 1,800 BP. Outside the Tibetan Plateau, the western Tibetan Plateau populations interacted with both South and Central Asian populations at least 2,000 years ago, and the South Asian-related genetic influence, despite being very limited, was from the Indus Valley Civilization (IVC) migrants in Central Asia instead of the IVC populations from the Indus Valley. In light of the new genetic data, our study revealed the complex population interconnections across and within the Tibetan Plateau.

Evolutionary profiles and complex admixture landscape in East Asia: New insights from modern and ancient Y chromosome variation perspectives

Human Y-chromosomes are characterized by nonrecombination and uniparental inheritance, carrying traces of human history evolution and admixture. Large-scale population-specific genomic sources based on advanced sequencing technologies have revolutionized our understanding of human Y chromosome diversity and its anthropological and forensic applications. Here, we reviewed and meta-analyzed the Y chromosome genetic diversity of modern and ancient people from China and summarized the patterns of founding lineages of spatiotemporally different populations associated with their origin, expansion, and admixture. We emphasized the strong association between our identified founding lineages and language-related human dispersal events correlated with the Sino-Tibetan, Altaic, and southern Chinese multiple-language families related to the Hmong-Mien, Tai-Kadai, Austronesian, and Austro-Asiatic languages. We subsequently summarize the recent advances in translational applications in forensic and anthropological science, including paternal biogeographical ancestry inference (PBGAI), surname investigation, and paternal history reconstruction. Whole-Y sequencing or high-resolution panels with high coverage of terminal Y chromosome lineages are essential for capturing the genomic diversity of ethnolinguistically diverse East Asians. Generally, we emphasized the importance of including more ethnolinguistically diverse, underrepresented modern and spatiotemporally different ancient East Asians in human genetic research for a comprehensive understanding of the paternal genetic landscape of East Asians with a detailed time series and for the reconstruction of a reference database in the PBGAI, even including new technology innovations of Telomere-to-Telomere (T2T) for new genetic variation discovery.

The Genetic Origin of the Indo-Europeans

The Yamnaya archaeological complex appeared around 3300BCE across the steppes north of the Black and Caspian Seas, and by 3000BCE reached its maximal extent from Hungary in the west to Kazakhstan in the east. To localize the ancestral and geographical origins of the Yamnaya among the diverse Eneolithic people that preceded them, we studied ancient DNA data from 428 individuals of which 299 are reported for the first time, demonstrating three previously unknown Eneolithic genetic clines. First, a "Caucasus-Lower Volga" (CLV) Cline suffused with Caucasus hunter-gatherer (CHG) ancestry extended between a Caucasus Neolithic southern end in Neolithic Armenia, and a steppe northern end in Berezhnovka in the Lower Volga. Bidirectional gene flow across the CLV cline created admixed intermediate populations in both the north Caucasus, such as the Maikop people, and on the steppe, such as those at the site of Remontnoye north of the Manych depression. CLV people also helped form two major riverine clines by admixing with distinct groups of European hunter-gatherers. A "Volga Cline" was formed as Lower Volga people mixed with upriver populations that had more Eastern hunter-gatherer (EHG) ancestry, creating genetically hyper-variable populations as at Khvalynsk in the Middle Volga. A "Dnipro Cline" was formed as CLV people bearing both Caucasus Neolithic and Lower Volga ancestry moved west and acquired Ukraine Neolithic hunter-gatherer (UNHG) ancestry to establish the population of the Serednii Stih culture from which the direct ancestors of the Yamnaya themselves were formed around 4000BCE. This population grew rapidly after 3750-3350BCE, precipitating the expansion of people of the Yamnaya culture who totally displaced previous groups on the Volga and further east, while admixing with more sedentary groups in the west. CLV cline people with Lower Volga ancestry contributed four fifths of the ancestry of the Yamnaya, but also, entering Anatolia from the east, contributed at least a tenth of the ancestry of Bronze Age Central Anatolians, where the Hittite language, related to the Indo-European languages spread by the Yamnaya, was spoken. We thus propose that the final unity of the speakers of the "Proto-Indo-Anatolian" ancestral language of both Anatolian and Indo-European languages can be traced to CLV cline people sometime between 4400-4000 BCE.

Origin and dispersal history of Hepatitis B virus in Eastern Eurasia

Hepatitis B virus is a globally distributed pathogen and the history of HBV infection in humans predates 10000 years. However, long-term evolutionary history of HBV in Eastern Eurasia remains elusive. We present 34 ancient HBV genomes dating between approximately 5000 to 400 years ago sourced from 17 sites across Eastern Eurasia. Ten sequences have full coverage, and only two sequences have less than 50% coverage. Our results suggest a potential origin of genotypes B and D in Eastern Asia. We observed a higher level of HBV diversity within Eastern Eurasia compared to Western Eurasia between 5000 and 3000 years ago, characterized by the presence of five different genotypes (A, B, C, D, WENBA), underscoring the significance of human migrations and interactions in the spread of HBV. Our results suggest the possibility of a transition from non-recombinant subgenotypes (B1, B5) to recombinant subgenotypes (B2 - B4). This suggests a shift in epidemiological dynamics within Eastern Eurasia over time. Here, our study elucidates the regional origins of prevalent genotypes and shifts in viral subgenotypes over centuries.

Differentiated adaptative genetic architecture and language-related demographical history in South China inferred from 619 genomes from 56 populations

BACKGROUND

The underrepresentation of human genomic resources from Southern Chinese populations limited their health equality in the precision medicine era and complete understanding of their genetic formation, admixture, and adaptive features. Besides, linguistical and genetic evidence supported the controversial hypothesis of their origin processes. One hotspot case was from the Chinese Guangxi Pinghua Han people (GPH), whose language was significantly similar to Southern Chinese dialects but whose uniparental gene pool was phylogenetically associated with the indigenous Tai-Kadai (TK) people. Here, we analyzed genome-wide SNP data in 619 people from four language families and 56 geographically different populations, in which 261 people from 21 geographically distinct populations were first reported here.

RESULTS

We identified significant population stratification among ethnolinguistically diverse Guangxi populations, suggesting their differentiated genetic origin and admixture processes. GPH shared more alleles related to Zhuang than Southern Han Chinese but received more northern ancestry relative to Zhuang. Admixture models and estimates of genetic distances showed that GPH had a close genetic relationship with geographically close TK compared to Northern Han Chinese, supporting their admixture origin hypothesis. Further admixture time and demographic history reconstruction supported GPH was formed via admixture between Northern Han Chinese and Southern TK people. We identified robust signatures associated with lipid metabolisms, such as fatty acid desaturases (FADS) and medically relevant loci associated with Mendelian disorder (GJB2) and complex diseases. We also explored the shared and unique selection signatures of ethnically different but linguistically related Guangxi lineages and found some shared signals related to immune and malaria resistance.

CONCLUSIONS

Our genetic analysis illuminated the language-related fine-scale genetic structure and provided robust genetic evidence to support the admixture hypothesis that can explain the pattern of observed genetic diversity and formation of GPH. This work presented one comprehensive analysis focused on the population history and demographical adaptative process, which provided genetic evidence for personal health management and disease risk prediction models from Guangxi people. Further large-scale whole-genome sequencing projects would provide the entire landscape of southern Chinese genomic diversity and their contributions to human health and disease traits.

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.

Evolutionary origin of germline pathogenic variants in human DNA mismatch repair genes

BACKGROUND

Mismatch repair (MMR) system is evolutionarily conserved for genome stability maintenance. Germline pathogenic variants (PVs) in MMR genes that lead to MMR functional deficiency are associated with high cancer risk. Knowing the evolutionary origin of germline PVs in human MMR genes will facilitate understanding the biological base of MMR deficiency in cancer. However, systematic knowledge is lacking to address the issue. In this study, we performed a comprehensive analysis to know the evolutionary origin of human MMR PVs.

METHODS

We retrieved MMR gene variants from the ClinVar database. The genomes of 100 vertebrates were collected from the UCSC genome browser and ancient human sequencing data were obtained through comprehensive data mining. Cross-species conservation analysis was performed based on the phylogenetic relationship among 100 vertebrates. Rescaled ancient sequencing data were used to perform variant calling for archeological analysis.

RESULTS

Using the phylogenetic approach, we traced the 3369 MMR PVs identified in modern humans in 99 non-human vertebrate genomes but found no evidence for cross-species conservation as the source for human MMR PVs. Using the archeological approach, we searched the human MMR PVs in over 5000 ancient human genomes dated from 45,045 to 100 years before present and identified a group of MMR PVs shared between modern and ancient humans mostly within 10,000 years with similar quantitative patterns.

CONCLUSION

Our study reveals that MMR PVs in modern humans were arisen within the recent human evolutionary history.

Reconstructing complex admixture history using a hierarchical model

Various methods have been proposed to reconstruct admixture histories by analyzing the length of ancestral chromosomal tracts, such as estimating the admixture time and number of admixture events. However, available methods do not explicitly consider the complex admixture structure, which characterizes the joining and mixing patterns of different ancestral populations during the admixture process, and instead assume a simplified one-by-one sequential admixture model. In this study, we proposed a novel approach that considers the non-sequential admixture structure to reconstruct admixture histories. Specifically, we introduced a hierarchical admixture model that incorporated four ancestral populations and developed a new method, called HierarchyMix, which uses the length of ancestral tracts and the number of ancestry switches along genomes to reconstruct the four-way admixture history. By automatically selecting the optimal admixture model using the Bayesian information criterion principles, HierarchyMix effectively estimates the corresponding admixture parameters. Simulation studies confirmed the effectiveness and robustness of HierarchyMix. We also applied HierarchyMix to Uyghurs and Kazakhs, enabling us to reconstruct the admixture histories of Central Asians. Our results highlight the importance of considering complex admixture structures and demonstrate that HierarchyMix is a useful tool for analyzing complex admixture events.

Ancient mitochondrial genome depicts sheep maternal dispersal and migration in Eastern Asia

Sheep have been one of the most important groups of animals since ancient times. However, the knowledge of their migration routes and genetic relationships is still poorly understood. To investigate sheep maternal migration histories alongside Eurasian communications routes, in this study, we obtain mitochondrial genomes (mitogenomes) from 17 sheep remains in 6 Chinese sites and 1 Uzbekistan site dated 4429-3100 years before present (BP). By obtaining the mitogenomes from the sheep (4429-3556 BP) found in the Tongtian Cave site in Xinjiang, Altai region of northwest China, our results support the emergence of haplogroup C sheep in Xinjiang as early as 4429-3556 BP. The combined phylogenetic analyses with extant ancient and modern sheep mitogenomes suggest that the Uzbekistan-Altai region may have been a migration hub for early sheep in eastern Asia. At least two migration events have taken place for sheep crossing Eurasia to China, one passing by Uzbekistan and Northwest China to the middle and lower reaches of the Yellow River at approximately 4000 BP and another following the Altai region to middle Inner Mongolia from 4429 BP to 2500 BP. Overall, this study provides further evidence for early sheep utilization and migration patterns in Eastern Asia.

Whole mitochondrial genome analysis in highland Tibetans: further matrilineal genetic structure exploration

Introduction: The Qinghai-Tibet Plateau is one of the last terrestrial environments conquered by modern humans. Tibetans are among the few high-altitude settlers in the world, and understanding the genetic profile of Tibetans plays a pivotal role in studies of anthropology, genetics, and archaeology. Methods: In this study, we investigated the maternal genetic landscape of Tibetans based on the whole mitochondrial genome collected from 145 unrelated native Lhasa Tibetans. Molecular diversity indices, haplotype diversity (HD), Tajima's D and Fu's Fs were calculated and the Bayesian Skyline Plot was obtained to determining the genetic profile and population fluctuation of Lhasa Tibetans. To further explore the genetic structure of Lhasa Tibetans, we collected 107 East Asian reference populations to perform principal component analysis (PCA), multidimensional scaling (MDS), calculated Fst values and constructed phylogenetic tree. Results: The maternal genetic landscape of Tibetans showed obvious East Asian characteristics, M9a (28.28%), R (11.03%), F1 (12.41%), D4 (9.66%), N (6.21%), and M62 (4.14%) were the dominant haplogroups. The results of PCA, MDS, Fst and phylogenetic tree were consistent: Lhasa Tibetans clustered with other highland Tibeto-Burman speakers, there was obvious genetic homogeneity of Tibetans in Xizang, and genetic similarity between Tibetans and northern Han people and geographically adjacent populations was found. In addition, specific maternal lineages of Tibetans also be determined in this study. Discussion: In general, this study further shed light on long-time matrilineal continuity on the Tibetan Plateau and the genetic connection between Tibetans and millet famers in the Yellow River Basin, and further revealed that multiple waves of population interaction and admixture during different historical periods between lowland and highland populations shaped the maternal genetic profile of Tibetans.

Inferring the population history of Kyrgyz in Xinjiang, Northwest China from genome-wide array genotyping

OBJECTIVES

Xinjiang plays a vital role in the trans-Eurasian population migration, language diffusion, and culture and technology exchange. However, the underrepresentation of Xinjiang's genomes has hindered a more comprehensive understanding of Xinjiang's genetic structure and population history.

MATERIALS AND METHODS

We collected and genotyped 70 southern Xinjiang's Kyrgyz (SXJK) individuals and combined the data with modern and ancient Eurasians published. We used allele-frequency methods, including PCA, ADMIXTURE, f-statistics, qpWave/qpAdm, ALDER, Treemix, and haplotype-shared methods including shared-IBD segments, fineSTRUCTURE, and GLOBETROTTER to unveil the fine-scale population structure and reconstruct admixture history.

RESULTS

We identified genetic substructure within the SXJK population with subgroups showing different genetic affinities to West and East Eurasians. All SXJK subgroups were suggested to have close genetic relationships with surrounding Turkic-speaking groups that is, Uyghur, Kyrgyz from north Xinjiang and Tajikistan, and Chinese Kazakh, suggesting a shared ancestry among those populations. Outgroup-f3 and symmetrical f4 statistics showed a high genetic affinity of SXJK to present-day Tungusic, Mongolic-speaking populations and Ancient Northeast Asian (ANA) related groups. Allele sharing and haplotype sharing profiles revealed the east-west admixture pattern of SXJK. The qpAdm-based admixture models showed that SXJK derived ancestry from East Eurasian (ANA and East Asian, 42.7%-83.3%) and West Eurasian (Western Steppe herders and Central Asian, 16.7%-57.3%), the recent east-west admixture event could be traced to 1000 years ago based on ALDER and GLOBETROTTER analysis.

DISCUSSION

The high genetic affinity of SXJK to present-day Tungusic and Mongolic-speaking populations and short-shared IBD segments indicated their shared common ancestry. SXJK harbored a close genetic affinity to ANA-related populations, indicating the Northeast Asian origin of SXJK. The West and East Eurasian admixture models observed in SXJK further provided evidence of the dynamic admixture history in Xinjiang. The east-west admixture pattern and the identified ancestral makeup of SXJK suggested a genetic continuity from some Iron Age Xinjiang populations to present-day SXJK.

Earliest systematic coal exploitation for fuel extended to ~3600 B.P

Coal has long fueled human civilizations. The history of systematic coal fuel exploitation has been traced back to the late third millennium before present (post-2500 B.P.). Although sporadic combustion of coal for fuel was reported in some prehistoric archaeological sites, evidence for the systematic exploitation of coal for fuel before 2500 B.P. remains lacking. Here, we report comprehensive understanding for the earliest systematic exploitation of coal for fuel at the Jirentaigoukou site in Xinjiang, northwestern China, at ~3600 B.P. The main body of the site witnessed systematic exploitation of bituminous coals, illustrating a complete chaîne opératoire with selective mining, planned storage, and extensive combustion. Our results transform the knowledge of energy history by extending the upper limit of the systematic exploitation of coal for fuel by approximately a millennium, and provide a precedent of energy transition under intense conflict between social demand and environmental deterioration.

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.

The genetic echo of the Tarim mummies in modern Central Asians

The diversity of Central Asians has been shaped by multiple migrations and cultural diffusion. Although ancient DNA studies have revealed the demographic changes of the Central Asian since the Bronze Age, the contribution of the ancient populations to the modern Central Asian remains opaque. Herein, we performed high-coverage sequencing of 131 whole genomes of Indo-European-speaking Tajik and Turkic-speaking Kyrgyz populations to explore their genomic diversity and admixture history. By integrating the ancient DNA data, we revealed more details of the origins and admixture history of Central Asians. We found that the major ancestry of present-day Tajik populations can be traced back to the admixture of the Bronze Age Bactria–Margiana Archaeological Complex and Andronovo-related populations. Highland Tajik populations further received additional gene flow from the Tarim mummies, an isolated ancient North Eurasian–related population. The West Eurasian ancestry of Kyrgyz is mainly derived from Historical Era populations in Xinjiang of China. Furthermore, the recent admixture signals detected in both Tajik and Kyrgyz are ascribed to the expansions of Eastern Steppe nomadic pastoralists during the Historical Era.

Genomic insight into the population history and biological adaptations of high-altitude Tibetan highlanders in Nagqu

Tibetan, one of the largest indigenous populations living in the high-altitude region of the Tibetan Plateau (TP), has developed a suite of physiological adaptation strategies to cope with the extreme highland environment in TP. Here, we reported genome-wide SNP data from 48 Kham-speaking Nagqu Tibetans and analyzed it with published data from 1,067 individuals in 167 modern and ancient populations to characterize the detailed Tibetan subgroup history and population substructure. Overall, the patterns of allele sharing and haplotype sharing suggested (1) the relatively genetic homogeny between the studied Nagqu Tibetans and ancient Nepalese as well as present-day core Tibetans from Lhasa, Nagqu, and Shigatse; and (2) the close relationship between our studied Kham-speaking Nagqu Tibetans and Kham-speaking Chamdo Tibetans. The fitted qpAdm models showed that the studied Nagqu Tibetans could be fitted as having the main ancestry from late Neolithic upper Yellow River millet farmers and deeply diverged lineages from Southern East Asians (represented by Upper Paleolithic Guangxi_Longlin and Laos_Hoabinhian), and a non-neglectable western Steppe herder-related ancestry (∼3%). We further scanned the candidate genomic regions of natural selection for our newly generated Nagqu Tibetans and the published core Tibetans via FST, iHS, and XP-EHH tests. The genes overlapping with these regions were associated with essential human biological functions such as immune response, enzyme activity, signal transduction, skin development, and energy metabolism. Together, our results shed light on the admixture and evolutionary history of Nagqu Tibetan populations.