Bronze and Iron Age population movements underlie Xinjiang population history

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.

Abstract Figure

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.