Protocol for a comprehensive pipeline to study ancient human genomes.

Ancient genomics has revolutionized our understanding of human evolution and migration history in recent years. Here, we present a protocol to prepare samples for ancient genomics research. We describe steps for releasing DNA from human remains, DNA library construction, hybridization capture, quantification, and sequencing. We then detail procedures for mapping sequence reads and population genetics analysis. This protocol also outlines challenges in extracting ancient DNA samples and authenticating ancient DNA to uncover the genetic history and diversity of ancient populations. For complete details on the use and execution of this protocol, please refer to Tao et al.1.

Tracing the fine-scale demographic history and recent admixture in Hmong-Mien speakers.

The linguistic, historical, and subsistent uniqueness of Hmong-Mien (HM) speakers offers a wonderful opportunity to investigate how these factors impact the genetic structure. The genetic differentiation among HM speakers and their population history are not well characterized. Here, we generate genome-wide data from 65 Yao ethnicity samples and analyze them with published data, particularly by leveraging haplotype-based methods. We determined that the fine-scale genetic substructure of HM speakers corresponds better with linguistic classification than with geography. Particularly, parallels between serial founder events and language differentiations can be observed in West Hmongic speakers. Multiple lines of evidence indicate that ~500-year-old GaoHuaHua individuals are most closely related to West Hmongic-speaking Bunu. The strong genetic bottleneck of some HM-speaking groups, especially Bunu, could potentially be associated with their long-term practice of swidden agriculture to some degree. The inferred admixture dates for most of the HM speakers overlap with the reign of the Ming dynasty (1368-1644 CE). Besides a common genetic origin for HM speakers, their genetic ancestry is shared primarily with neighboring Han Chinese and Tai-Kadai speakers in south China. In conclusion, our analyses reveal that recent isolation and admixture events have contributed to the genetic population history of present-day HM speakers.

Malaria resistance-related biological adaptation and complex evolutionary footprints inferred from one integrative Tai-Kadai-related genomic resource.

Pathogen‒host adaptative interactions and complex population demographical processes, including admixture, drift, and Darwen selection, have considerably shaped the Neolithic-to-Modern Western Eurasian population structure and genetic susceptibility to modern human diseases. However, the genetic footprints of evolutionary events in East Asia remain unknown due to the underrepresentation of genomic diversity and the design of large-scale population studies. We reported one aggregated database of genome-wide SNP variations from 796 Tai-Kadai (TK) genomes, including that of Bouyei first reported here, to explore the genetic history, population structure, and biological adaptative features of TK people from southern China and Southeast Asia. We found geography-related population substructure among TK people using the state-of-the-art population genetic structure reconstruction techniques based on the allele frequency spectrum and haplotype-resolved phased fragments. We found that the northern TK people from Guizhou harbored one TK-dominant ancestry maximized in the Bouyei people, and the southern TK people from Thailand were more influenced by Southeast Asians and indigenous people. We reconstructed fitted admixture models and demographic graphs, which showed that TK people received gene flow from ancient southern rice farmer-related lineages related to the Hmong-Mien and Austroasiatic people and from northern millet farmers associated with the Sino-Tibetan people. Biological adaptation focused on our identified unique TK lineages related to Bouyei, which showed many adaptive signatures conferring Malaria resistance and low-rate lipid metabolism. Further gene enrichment, the allele frequency distribution of derived alleles, and their correlation with the incidence of Malaria further confirmed that CR1 played an essential role in the resistance of Malaria in the ancient "Baiyue" tribes.

Ancient genome of the Chinese Emperor Wu of Northern Zhou.

Emperor Wu (, Wudi) of the Xianbei-led Northern Zhou dynasty, named Yuwen Yong (, 543-578 CE), was a highly influential emperor who reformed the system of regional troops, pacified the Turks, and unified the northern part of the country. His genetic profile and physical characteristics, including his appearance and potential diseases, have garnered significant interest from the academic community and the public. In this study, we have successfully generated a 0.343×-coverage genome of Wudi with 1,011,419 single-nucleotide polymorphisms (SNPs) on the 1240k panel. By analyzing pigmentation-relevant SNPs and conducting cranial CT-based facial reconstruction, we have determined that Wudi possessed a typical East or Northeast Asian appearance. Furthermore, pathogenic SNPs suggest Wudi faced an increased susceptibility to certain diseases, such as stroke. Wudi shared the closest genetic relationship with ancient Khitan and Heishui Mohe samples and modern Daur and Mongolian populations but also showed additional affinity with Yellow River (YR) farmers. We estimated that Wudi derived 61% of his ancestry from ancient Northeast Asians (ANAs) and nearly one-third from YR farmer-related groups. This can likely be attributed to continuous intermarriage between Xianbei royal families, and local Han aristocrats.1,2 Furthermore, our study has revealed genetic diversities among available ancient Xianbei individuals from different regions, suggesting that the formation of the Xianbei was a dynamic process influenced by admixture with surrounding populations.

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.

Forager-farmer transition at the crossroads of East and Southeast Asia 4900 years ago.

The southward expansion of East Asian farmers profoundly influenced the social evolution of Southeast Asia by introducing cereal agriculture. However, the timing and routes of cereal expansion in key regions are unclear due to limited empirical evidence. Here we report macrofossil, microfossil, multiple isotopic (C/N/Sr/O) and paleoproteomic data directly from radiocarbon-dated human samples, which were unearthed from a site in Xingyi in central Yunnan and which date between 7000 and 3300 a BP. Dietary isotopes reveal the earliest arrival of millet ca. 4900 a BP, and greater reliance on plant and animal agriculture was indicated between 3800 and 3300 a BP. The dietary differences between hunter-gatherer and agricultural groups are also evident in the metabolic and immune system proteins analysed from their skeletal remains. The results of paleoproteomic analysis indicate that humans had divergent biological adaptations, with and without farming. The combined application of isotopes, archaeobotanical data and proteomics provides a new approach to documenting dietary and health changes across major subsistence transitions.

Unraveling the paternal genetic structure and forensic traits of the Hui population in Liaoning Province, China using Y-chromosome analysis.

The Hui people are the second-largest ethnic minority in China, and they are distributed throughout the country. A previous study explored the paternal genetic structure of the Hui population in nine different regions of China, but it overlooked the Liaoning province. In this study, we examined the paternal genetic makeup and forensic traits of the Hui population in Liaoning province by analyzing 157 Y-chromosome single nucleotide polymorphisms (Y-SNPs) and 26 short tandem repeats (Y-STRs). We successfully genotyped 282 unrelated male individuals from the Hui population of Liaoning province using the SNaPshot® single base extension assay and Goldeneye™ Y26 system kit (PEOPLESPOT R&D, Beijing, China). The results revealed high haplotypic diversity (0.9998) and identified 46 terminal haplogroups for the Hui population. Additional analyses, such as heat maps, principal component analysis (PCA), genetic distance (FST), Multidimensional scaling (MDS) analysis, and median-joining network (MJ) analysis, showed that the Hui population could be classified into three groups: Northwest Hui populations (NWH), including Liaoning, Xinjiang, Qinghai, Gansu, Ningxia, Shaanxi, and Henan; Hui populations from Sichuan and Shandong (SSH); and Yunnan Hui populations (YNH). Pairwise genetic distance (Rst) comparisons with other Chinese populations revealed that the Hui population displayed genetic affinity with the Han population. The comprehensive understanding of the Hui population in Liaoning province, explored by Y-SNPs and Y-STRs, can be utilized to interpret their genetic structure and enhance the accuracy of forensic databases.

Genomic formation of Tibeto-Burman speaking populations in Guizhou, Southwest China.

Sino-Tibetan is the most prominent language family in East Asia. Previous genetic studies mainly focused on the Tibetan and Han Chinese populations. However, due to the sparse sampling, the genetic structure and admixture history of Tibeto-Burman-speaking populations in the low-altitude region of Southwest China still need to be clarified. We collected DNA from 157 individuals from four Tibeto-Burman-speaking groups from the Guizhou province in Southwest China. We genotyped the samples at about 700,000 genome-wide single nucleotide polymorphisms. Our results indicate that the genetic variation of the four Tibeto-Burman-speaking groups in Guizhou is at the intermediate position in the modern Tibetan-Tai-Kadai/Austronesian genetic cline. This suggests that the formation of Tibetan-Burman groups involved a large-scale gene flow from lowland southern Chinese. The southern ancestry could be further modelled as deriving from Vietnam's Late Neolithic-related inland Southeast Asia agricultural populations and Taiwan's Iron Age-related coastal rice-farming populations. Compared to the Tibeto-Burman speakers in the Tibetan-Yi Corridor reported previously, the Tibeto-Burman groups in the Guizhou region received additional gene flow from the southeast coastal area of China. We show a difference between the genetic profiles of the Tibeto-Burman speakers of the Tibetan-Yi Corridor and the Guizhou province. Vast mountain ranges and rivers in Southwest China may have decelerated the westward expansion of the southeast coastal East Asians. Our results demonstrate the complex genetic profile in the Guizhou region in Southwest China and support the multiple waves of human migration in the southern area of East Asia.

Ancient genomes reveal millet farming-related demic diffusion from the Yellow River into southwest China.

The study of southwest China is vital for understanding the dispersal and development of farming because of the coexistence of millet and rice in this region since the Neolithic period.1,2 However, the process of the Neolithic transition in southwest China is largely unknown, mainly due to the lack of ancient DNA from the Neolithic period. Here, we report genome-wide data from 11 human samples from the Gaoshan and Haimenkou sites with mixed farming of millet and rice dating to between 4,500 and 3,000 years before present in southwest China. The two ancient groups derived approximately 90% of their ancestry from the Neolithic Yellow River farmers, suggesting a demic diffusion of millet farming to southwest China. We inferred their remaining ancestry to be derived from a Hòabìnhian-related hunter-gatherer lineage. We did not detect rice farmer-related ancestry in the two ancient groups, which indicates that they likely adopted rice farming without genetic assimilation. We, however, observed rice farmer-related ancestry in the formation of some present-day Tibeto-Burman populations. Our results suggested the occurrence of both demic and cultural diffusion in the development of Neolithic mixed farming in some parts of southwest China.

Distinguished biological adaptation architecture aggravated population differentiation of Tibeto-Burman-speaking people.

Tibeto-Burman (TB) people have endeavored to adapt to the hypoxic, cold, and high-UV high-altitude environments in the Tibetan Plateau and complex disease exposures in lowland rainforests since the late Paleolithic period. However, the full landscape of genetic history and biological adaptation of geographically diverse TB-speaking people, as well as their interaction mechanism, remain unknown. Here, we generate a whole-genome meta-database of 500 individuals from 39 TB-speaking populations and present a comprehensive landscape of genetic diversity, admixture history, and differentiated adaptative features of geographically different TB-speaking people. We identify genetic differentiation related to geography and language among TB-speaking people, consistent with their differentiated admixture process with incoming or indigenous ancestral source populations. A robust genetic connection between the Tibetan-Yi corridor and the ancient Yellow River people supports their Northern China origin hypothesis. We finally report substructure-related differentiated biological adaptative signatures between highland Tibetans and Loloish speakers. Adaptative signatures associated with the physical pigmentation (EDAR and SLC24A5) and metabolism (ALDH9A1) are identified in Loloish people, which differed from the high-altitude adaptative genetic architecture in Tibetan. TB-related genomic resources provide new insights into the genetic basis of biological adaptation and better reference for the anthropologically informed sampling design in biomedical and genomic cohort research.

Extensive genetic admixture between Tai-Kadai-speaking people and their neighbours in the northeastern region of the Yungui Plateau inferred from genome-wide variations.

BACKGROUND: Yungui Plateau in Southwest China is characterized by multi-language and multi-ethnic communities and is one of the regions with the wealthiest ethnolinguistic, cultural and genetic diversity in East Asia. There are numerous Tai-Kadai (TK)-speaking populations, but their detailed evolutionary history and biological adaptations are still unclear. RESULTS: Here, we genotyped genome-wide SNP data of 77 unrelated TK-speaking Zhuang and Dong individuals from the Yungui Plateau and explored their detailed admixture history and adaptive features using clustering patterns, allele frequency differentiation and sharing haplotype patterns. TK-speaking Zhuang and Dong people in Guizhou are closely related to geographically close TK and Hmong-Mien (HM)-speaking populations. Besides, we identified that Guizhou TK-speaking people have a close genetic relationship with Austronesian (AN)-speaking Atayal and Paiwan people, which is supported by the common origin of the ancient Baiyue tribe. We additionally found subtle genetic differences among the newly studied TK people and previously reported Dais via the fine-scale genetic substructure analysis based on the shared haplotype chunks. Finally, we identified specific selection candidate signatures associated with several essential human immune systems and neurological disorders, which could provide evolutionary evidence for the allele frequency distribution pattern of genetic risk loci. CONCLUSIONS: Our comprehensive genetic characterization of TK people suggested the strong genetic affinity within TK groups and extensive gene flow with geographically close HM and Han people. We also provided genetic evidence that supported the common origin hypothesis of TK and AN people. The best-fitted admixture models further suggested that ancestral sources from northern millet farmers and southern inland and coastal people contributed to the formation of the gene pool of the Zhuang and Dong people.

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.

The genetic structure and admixture of Manchus and Koreans in northeast China.

BACKGROUND: The fine-scale genetic profiles and population history of Manchus and Koreans remain unclear. AIM: To infer a fine-scale genetic structure and admixture of Manchu and Korean populations. SUBJECTS AND METHODS: We collected and genotyped 16 Manchus from Liaoning and 18 Koreans from Jilin province with about 700K genome-wide SNPs. We analysed the data using principal component analysis (PCA), ADMIXTURE, Fst, TreeMix, f-statistics, qpWave, and qpAdm. RESULTS: Manchus and Koreans showed a genetic affinity with northern East Asians. Chinese Koreans showed a long-term genetic continuity with Bronze Age populations from the West Liao River and had a strong affinity with Koreans in South Korea and Japan. Manchus had a different genetic profile compared with other Tungusic populations since the Manchus received additional genetic influence from the southern Chinese but didn't have West Eurasian-related admixture. CONCLUSIONS: The genetic formation of Manchus involving southern Chinese was consistent with the extensive interactions between Manchus and populations from central and southern China. The large-scale genetic continuity between ancient West Liao River farmers and Koreans highlighted the role farming expansion played in the peopling of the Korean Peninsula.

Effect of Si on the Oxidation Behaviors of Ti3 Al1-x Six C2 at 1000 °C.

In this work, Ti3 Al1-x Six C2 (x=0, 0.2, 0.4, and 0.6) with Al/Si solid solution structure are synthesized, and the effects of Si on their oxidation behaviors at 1000 °C are evaluated. The addition of Si not only contributes to the formation of Ti5 Si3 impurity but also affects the composition of the oxide scale. Particularly, the incorporation of Si in the TiO2 lattice is demonstrated, which alters the formation energy of the (110) plane in TiO2 , thus leading to the preferential growth of Si-doped TiO2 to dendritic congeries. Moreover, the Si addition is believed to affect mass transportation during the oxidation process, which accelerates the formation of a continuous Al2 O3 layer in the oxide scale. With an optimized Si content, the oxidation of Ti3 Al1-x Six C2 is restrained. However, with excess Si content, the continuity of the resulting Al2 O3 layer is destroyed, thus the oxidation rate rises again.

Evaluation of facial hair-associated SNPs: a pilot study on male Pakistani Punjabi population.

Variation in facial hair is one of the most conspicuous features of facial appearance, particularly in South Asia and Middle East countries. A genome-wide association study in Latin Americans has identified multiple genetic variants at distinct loci being associated with facial hair traits including eyebrow thickness, beard thickness, and monobrow. In this pilot study, we have evaluated 16 SNPs associated with facial hair traits in 58 male individuals from the Punjabi population of Pakistan. In our sample, rs365060 in EDAR and rs12597422 in FTO showed significant association with monobrow, rs6684877 in MACF1 showed significant association with eyebrow thickness, and two SNPs in LOC105379031 (rs9654415 and rs7702331) showed significant association with beard thickness. Our results also suggest that genetic association may vary between ethnic groups and geographic regions. Although more data are needed to validate our results, our findings are of value in forensic molecular photofitting research in Pakistan.

Cultural and demic co-diffusion of Tubo Empire on Tibetan Plateau.

A high point of Tibetan Plateau (TP) civilization, the expansive Tubo Empire (618-842 AD) wielded great influence across ancient western China. However, whether the Tubo expansion was cultural or demic remains unclear due to sparse ancient DNA sampling. Here, we reported ten ancient genomes at 0.017- to 0.867-fold coverages from the Dulan site with typical Tubo archaeological culture dating to 1308-1130 BP. Nine individuals from three different grave types have close relationship with previously reported ancient highlanders from the southwestern Himalayas and modern core-Tibetan populations. A Dulan-related Tubo ancestry contributed overwhelmingly (95%-100%) to the formation of modern Tibetans. A genetic outlier with dominant Eurasian steppe-related ancestry suggesting a potential population movement into the Tubo-controlled regions from Central Asia. Together with archeological evidence from burial styles and customs, our study suggested the impact of the Tubo empire on the northeast edge of the TP involved both cultural and demic diffusion.

Genetic diversity and forensic application of Y-filer STRs in four major ethnic groups of Pakistan.

17 Y-chromosomal STRs which are part of the Yfiler Amplification Kit were investigated in 493 unrelated Pakistani individuals belonging to the Punjabi, Sindhi, Baloch, and Pathan ethnic groups. We have assessed the forensic parameters and population genetic structure for each group. Among the 493 unrelated individuals from four ethnic groups (128 Baloch, 122 Pathan, 108 Punjabi, and 135 Sindhi), 82 haplotypes were observed with haplotype diversity (HD) of 0.9906 in Baloch, 102 haplotypes with HD value of 0.9957 in Pathans, 80 haplotypes with HD value of 0.9924 in Punjabi, and 105 haplotypes with HD value of 0.9945 in the Sindhi population. The overall gene diversity for Baloch, Pathan, Punjabi, and Sindhi populations was 0.6367, 0.6479, 0.6657, and 0.6112, respectively. The results had shown us that Pakistani populations do not have a unique set of genes but share the genetic affinity with regional (Central Asia and Northern India) populations. The observed low gene diversity (heterozygosity) values may be because of endogamy trends and this observation is equally supported by the results of forensic parameters which are mostly static across 4 combinations (minimal STRs, extended 11 Y-STRs, Powerplex 12 Y System, and Yfiler 17 Y-STRs) of STRs in these four populations.