New insights into the fine-scale history of western-eastern admixture of the northwestern Chinese population in the Hexi Corridor via genome-wide genetic legacy

Differentiated genomic footprints suggest isolation and long-distance migration of Hmong-Mien populations

BACKGROUND

The underrepresentation of Hmong-Mien (HM) people in Asian genomic studies has hindered our comprehensive understanding of the full landscape of their evolutionary history and complex trait architecture. South China is a multi-ethnic region and indigenously settled by ethnolinguistically diverse HM, Austroasiatic (AA), Tai-Kadai (TK), Austronesian (AN), and Sino-Tibetan (ST) people, which is regarded as East Asia's initial cradle of biodiversity. However, previous fragmented genetic studies have only presented a fraction of the landscape of genetic diversity in this region, especially the lack of haplotype-based genomic resources. The deep characterization of demographic history and natural-selection-relevant genetic architecture of HM people was necessary.

RESULTS

We reported one HM-specific genomic resource and comprehensively explored the fine-scale genetic structure and adaptative features inferred from the genome-wide SNP data of 440 HM individuals from 33 ethnolinguistic populations, including previously unreported She. We identified solid genetic differentiation between HM people and Han Chinese at 7.64‒15.86 years ago (kya) and split events between southern Chinese inland (Miao/Yao) and coastal (She) HM people in the middle Bronze Age period and the latter obtained more gene flow from Ancient Northern East Asians. Multiple admixture models further confirmed that extensive gene flow from surrounding ST, TK, and AN people entangled in forming the gene pool of Chinese coastal HM people. Genetic findings of isolated shared unique ancestral components based on the sharing alleles and haplotypes deconstructed that HM people from the Yungui Plateau carried the breadth of previously unknown genomic diversity. We identified a direct and recent genetic connection between Chinese inland and Southeast Asian HM people as they shared the most extended identity-by-descent fragments, supporting the long-distance migration hypothesis. Uniparental phylogenetic topology and network-based phylogenetic relationship reconstruction found ancient uniparental founding lineages in southwestern HM people. Finally, the population-specific biological adaptation study identified the shared and differentiated natural selection signatures among inland and coastal HM people associated with physical features and immune functions. The allele frequency spectrum of cancer susceptibility alleles and pharmacogenomic genes showed significant differences between HM and northern Chinese people.

CONCLUSIONS

Our extensive genetic evidence combined with the historical documents supported the view that ancient HM people originated from the Yungui regions associated with ancient "Three-Miao tribes" descended from the ancient Daxi-Qujialing-Shijiahe people. Then, some have recently migrated rapidly to Southeast Asia, and some have migrated eastward and mixed respectively with Southeast Asian indigenes, Liangzhu-related coastal ancient populations, and incoming southward ST people. Generally, complex population migration, admixture, and adaptation history contributed to the complicated patterns of population structure of geographically diverse HM people.

A proof-of-principle study: The potential application of MiniHap biomarkers in ancestry inference based on the QNome nanopore sequencing

Haplotyped SNPs convey forensic-related information, and microhaplotypes (MHs), as the most representative of this kind of marker, have proved the potential value for human forensics. In recent years, nanopore sequencing technology has developed rapidly, with its outstanding ability to sequence long continuous DNA fragments and obtain phase information, making the detection of longer haplotype marker possible. In this proof-of-principle study, we proposed a new type of forensic marker, MiniHap, based on five or more SNPs within a molecular distance less than 800 bp, and investigated the haplotype data of 56 selected MiniHaps in five Chinese populations using the QNome nanopore sequencing. The sequencing performance, allele (haplotype) frequencies, forensic parameters, effective number of alleles (Ae), and informativeness (In) were subsequently calculated. In addition, we performed principal component analysis (PCA), phylogenetic tree, and structure analysis to investigate the population genetic relationships and ancestry components among the five investigated populations and 26 worldwide populations. MiniHap-04 exhibited remarkable forensic efficacy, with 148 haplotypes reported and the Ae was 66.9268. In addition, the power of discrimination (PD) was 0.9934, the probability of exclusion (PE) was 0.9898, and the In value was 0.7893. Of the 56 loci, 85.71% had PD values above 0.85, 66.07% had PE values above 0.54, 67.86% had Ae values over 7.0%, and 55.36% were with In values above 0.2 across all samples, indicating that most of the MiniHaps are suitable for individual identification, paternity testing, mixture deconvolution, and ancestry inference. Moreover, the results of PCA, phylogenetic tree and structure analysis demonstrated that this MiniHap panel had the competency in continental population ancestry inference, but the differentiation within intracontinental/linguistically restricted subpopulations was not ideal. Such findings suggested that the QNome device for MiniHap detection was feasible and this novel marker has the potential in ancestry inference. Yet, the establishment of a more comprehensive database with sufficient reference population data remains necessary to screen more suitable MiniHaps.

Genetic Polymorphism and Population Genetic Structure Analysis of 21 Autosomal STR Loci for a Han-Chinese Population from Luzhou of Southwest China

The Han nationality is an ancient and populous people, and different places in China may have their distinct group relationships. Luzhou area, as a crossroads of several provinces in Southwest China, lacks autosomal short tandem repeat (STR) research and population genetics research. In this study, 21 autosomal STR loci were evaluated in 1959 Han-Chinese individuals from Luzhou area. There was no substantial linkage disequilibrium (LD) among the 21 autosomal STR markers, and all markers were in Hardy-Weinberg equilibrium (HWE). The total discrimination power (TDP) and cumulative probability of exclusion (CPE) of the 21 autosomal STR loci were calculated to be 1-9.8505 × 10^-16 and 1-1.9406 × 10^-9, respectively. There were 333 alleles for 21 STRs with allelic frequencies ranging from 0.00026 to 0.51302, and the number of alleles ranged from 7 in locus TPOX to 29 in locus Penta E. According to the results of population comparison and population differentiation, historical influences, geographical distribution, cultural integration, and economic development may have an impact on the Luzhou Han population and other Chinese populations. These 21 STR loci were found to enrich autosomal STR information in forensic databases and provide highly informative polymorphisms for our forensic practice in China, including personal identification and parentage testing.

Genome-wide allele and haplotype-sharing patterns suggested one unique Hmong-Mein-related lineage and biological adaptation history in Southwest China

BACKGROUND

Fine-scale genetic structure of ethnolinguistically diverse Chinese populations can fill the gap in the missing diversity and evolutionary landscape of East Asians, particularly for anthropologically informed Chinese minorities. Hmong-Mien (HM) people were one of the most significant indigenous populations in South China and Southeast Asia, which were suggested to be the descendants of the ancient Yangtze rice farmers based on linguistic and archeological evidence. However, their deep population history and biological adaptative features remained to be fully characterized.

OBJECTIVES

To explore the evolutionary and adaptive characteristics of the Miao people, we genotyped genome-wide SNP data in Guizhou HM-speaking populations and merged it with modern and ancient reference populations via a comprehensive population genetic analysis and evolutionary admixture modeling.

RESULTS

The overall genetic admixture landscape of Guizhou Miao showed genetic differentiation between them and other linguistically diverse Guizhou populations. Admixture models further confirmed that Miao people derived their primary ancestry from geographically close Guangxi Gaohuahua people. The estimated identity by descent and effective population size confirmed a plausible population bottleneck, contributing to their unique genetic diversity and population structure patterns. We finally identified several natural selection candidate genes associated with several biological pathways.

CONCLUSIONS

Guizhou Miao possessed a specific genetic structure and harbored a close genetic relationship with geographically close southern Chinese indigenous populations and Guangxi historical people. Miao people derived their major ancestry from geographically close Guangxi Gaohuahua people and experienced a plausible population bottleneck which contributed to the unique pattern of their genetic diversity and structure. Future ancient DNA from Shijiahe and Qujialing will provide new insights into the origin of the Miao people.

Genetic structure and demographic history of Northern Han people in Liaoning Province inferred from genome-wide array data

In this study, we used typical and advanced population genetic analysis methods [principal component analysis (PCA), ADMIXTURE, FST, f3-statistics, f4-statistics, qpAdm/qpWave, qpGraph, ALDER (Admixture-induced Linkage Disequilibrium for Evolutionary Relationships) and TreeMix] to explore the genetic structure of 80 Han individuals from four different cities in Liaoning Province and reconstruct their demographic history based on the newly generated genome-wide data. We found that Liaoning Han people have genetic similarities with other northern Han people (Shandong, Henan, and Shanxi) and Liaoning Manchu people. Millet farmers in the Yellow River Basin (YRB) and the West Liao River Basin (WLRB) (57–98%) and hunter-gatherers in the Mongolian Plateau (MP) and the Amur River Basin (ARB) (40–43%) are the main ancestral sources of the Liaoning Han people. Our study further supports the “northern origin hypothesis”; YRB-related ancestry accounts for 83–98% of the genetic makeup of the Liaoning Han population. There are clear genetic influences of northern East Asian populations in the Liaoning Han people, ancient Northeast Asian-related ancestry is another dominant ancestral component, and large-scale population admixture has happened between Tungusic Manchu people and Han people. There are genetic differences among the Liaoning Han people, and we found that these differences are associated with different migration routes of Hans during the “Chuang Guandong” period in historical records.

Genetic admixture history and forensic characteristics of Tibeto-Burman-speaking Qiang people explored via the newly developed Y-STR panel and genome-wide SNP data

Fine-scale patterns of population genetic structure and diversity of ethnolinguistically diverse populations are important for biogeographical ancestry inference, kinship testing, and development and validation of new kits focused on forensic personal identification. Analyses focused on forensic markers and genome-wide single nucleotide polymorphism (SNP) data can provide new insights into the origin, admixture processes, and forensic characteristics of targeted populations. Qiang people had a large sample size among Tibeto-Burmanspeaking populations, which widely resided in the middle latitude of the Tibetan Plateau. However, their genetic structure and forensic features have remained uncharacterized because of the paucity of comprehensive genetic analyses. Here, we first developed and validated the forensic performance of the AGCU-Y30 Y-short tandem repeats (STR) panel, which contains slowly and moderately mutating Y-STRs, and then we conducted comprehensive population genetic analyses based on Y-STRs and genome-wide SNPs to explore the admixture history of Qiang people and their neighbors. The validated results of this panel showed that the new Y-STR kit was sensitive and robust enough for forensic applications. Haplotype diversity (HD) ranging from 0.9932 to 0.9996 and allelic frequencies ranging from 0.001946 to 0.8326 in 514 Qiang people demonstrated that all included markers were highly polymorphic in Tibeto-Burman people. Population genetic analyses based on Y-STRs [RST, FST, multidimensional scaling (MDS) analysis, neighboring-joining (NJ) tree, principal component analysis (PCA), and median-joining network (MJN)] revealed that the Qiang people harbored a paternally close relationship with lowland Tibetan-Yi corridor populations. Furthermore, we conducted a comprehensive population admixture analysis among modern and ancient Eurasian populations based on genome-wide shared SNPs. We found that the Qiang people were a genetically admixed population and showed closest relationship with Tibetan and Neolithic Yellow River farmers. Admixture modeling showed that Qiang people shared the primary ancestry related to Tibetan, supporting the hypothesis of common origin between Tibetan and Qiang people from North China.

Genomic Insights Into the Unique Demographic History and Genetic Structure of Five Hmong-Mien-Speaking Miao and Yao Populations in Southwest China

Southern China was the original center of multiple ancestral populations related to modern Hmong-Mien, Tai-Kadai, Austroasiatic, and Austronesian people. More recent genetic surveys have focused on the fine-scale genetic structure and admixture history of southern Chinese populations, but the genetic formation and diversification of Hmong-Mien speakers are far from clear due to the sparse genetic sampling. Here, we reported nearly 700,000 single-nucleotide polymorphisms (SNPs) data from 130 Guizhou Miao and Yao individuals. We used principal component analysis, ADMIXTURE, f-statistics, qpAdm, phylogenetic tree, fineSTRUCTURE, and ALDER to explore the fine-scale population genetic structure and admixture pattern of Hmong-Mien people. The sharing allele patterns showed that our studied populations had a strong genetic affinity with ancient and modern groups from southern and southeastern East Asia. We identified one unique ancestry component maximized in Yao people, which widely existed in other Hmong-Mien-speaking populations in southern China and Southeast Asia and ancient samples of Guangxi. Guizhou Hmong-Mien speakers harbored the dominant proportions of ancestry related to southern indigenous East Asians and minor proportions of northern ancestry related to Yellow River farmers, suggesting the possibility of genetic admixture between Hmong-Mien people and recent southward Sino-Tibetan-related populations. Furthermore, we found a genetic substructure among geographically different Miao and Yao people in Leishan and Songtao. The Yao and Miao people in Leishan harbored more southern East Asian ancestry, but Miao in Songtao received more northern East Asian genetic influence. We observed high mtDNA but low Y-chromosome diversity in studied Hmong-Mien groups, supporting the role of sex-specific residence in influencing human genetic variation. Our data provide valuable clues for further exploring population dynamics in southern China.

Forensic Features and Genetic Structure Analyses of the Beijing Han Nationality Disclosed by a Self-Developed Panel Containing a Series of Ancestry Informative Deletion/Insertion Polymorphism Loci

The utilization of the ancestry informative markers to disclose the ancestral composition of a certain population and explore the genetic affinities between diverse populations is beneficial to inferring the biogeographic ancestry of unknown individuals and assisting in case detection, as well as avoiding the impacts of population stratification during genome-wide association analysis studies. In the present study, we applied an in-house ancestry informative deletion/insertion polymorphic multiplex amplification system to investigate the ancestral compositions of the Beijing Han population and analyze the genetic relationships between the Beijing Han population and 31 global reference populations. The results demonstrated that 32 loci of this self-developed panel containing 39 loci significantly contributed to the inference of genetic information for the Beijing Han population. The results of multiple population genetics statistical analyses indicated that the ancestral component and genetic architecture of the Beijing Han population were analogous to the reference East Asian populations, and that the Beijing Han population was genetically close to the reference East Asian populations.

Fine-Scale Population Admixture Landscape of Tai–Kadai-Speaking Maonan in Southwest China Inferred From Genome-Wide SNP Data

Guizhou Province harbors extensive ethnolinguistic and cultural diversity with Sino-Tibetan-, Hmong–Mien-, and Tai–Kadai-speaking populations. However, previous genetic analyses mainly focused on the genetic admixture history of the former two linguistic groups. The admixture history of Tai–Kadai-speaking populations in Guizhou needed to be characterized further. Thus, we genotyped genome-wide SNP data from 41 Tai–Kadai-speaking Maonan people and made a comprehensive population genetic analysis to explore their genetic origin and admixture history based on the pattern of the sharing alleles and haplotypes. We found a genetic affinity among geographically different Tai–Kadai-speaking populations, especially for Guizhou Maonan people and reference Maonan from Guangxi. Furthermore, formal tests based on the f₃/f₄-statistics further identified an adjacent connection between Maonan and geographically adjacent Hmong–Mien and Sino-Tibetan people, which was consistent with their historically documented shared material culture (Zhang et al., iScience, 2020, 23, 101032). Fitted qpAdm-based two-way admixture models with ancestral sources from northern and southern East Asians demonstrated that Maonan people were an admixed population with primary ancestry related to Guangxi historical people and a minor proportion of ancestry from Northeast Asians, consistent with their linguistically supported southern China origin. Here, we presented the landscape of genetic structure and diversity of Maonan people and a simple demographic model for their evolutionary process. Further whole-genome-sequence–based projects can be presented with more detailed information about the population history and adaptative history of the Guizhou Maonan people.

Genomic Insights Into the Population History and Biological Adaptation of Southwestern Chinese Hmong-Mien People

Hmong-Mien (HM) -speaking populations, widely distributed in South China, the north of Thailand, Laos, and Vietnam, have experienced different settlement environments, dietary habits, and pathogenic exposure. However, their specific biological adaptation remained largely uncharacterized, which is important in the population evolutionary genetics and Trans-Omics for regional Precision Medicine. Besides, the origin and genetic diversity of HM people and their phylogenetic relationship with surrounding modern and ancient populations are also unknown. Here, we reported genome-wide SNPs in 52 representative Miao people and combined them with 144 HM people from 13 geographically representative populations to characterize the full genetic admixture and adaptive landscape of HM speakers. We found that obvious genetic substructures existed in geographically different HM populations; one localized in the HM clines, and others possessed affinity with Han Chinese. We also identified one new ancestral lineage specifically existed in HM people, which spatially distributed from Sichuan and Guizhou in the north to Thailand in the south. The sharing patterns of the newly identified homogenous ancestry component combined the estimated admixture times via the decay of linkage disequilibrium and haplotype sharing in GLOBETROTTER suggested that the modern HM-speaking populations originated from Southwest China and migrated southward in the historic period, which is consistent with the reconstructed phenomena of linguistic and archeological documents. Additionally, we identified specific adaptive signatures associated with several important human nervous system biological functions. Our pilot work emphasized the importance of anthropologically informed sampling and deeply genetic structure reconstruction via whole-genome sequencing in the next step in the deep Chinese Population Genomic Diversity Project (CPGDP), especially in the regions with rich ethnolinguistic diversity.

The genomic history of southwestern Chinese populations demonstrated massive population migration and admixture among proto-Hmong-Mien speakers and incoming migrants

Southwest China was the crossroad for the initial settler people of East Asia, which shows the highest diversity in languages and genetics. This region played a significant role in the formation of the genetic makeup of the proto-Hmong-Mien-speaking people and in the north-to-south human expansion during the Neolithic-to-historic transformation. Their genetic history covering migration events and the admixture processes still needs to be further explored. Therefore, in the current study, we have generated genome-wide data from three genomic aspects covering autosomal, mitochondrial and Y-chromosomal regions in 260 Hmong-Mien, Tibeto-Burman, and Sinitic people from 29 different southwestern Chinese groups, and further analyzed them with 2676 published modern and ancient Eurasian genomes. Here, we have noticed a new southwestern East Asian genetic cline composed of the Hmong-Mien-specific ancestry enriched in modern Hmong and Pathen. This newly identified southern inland East Asian lineage contributed to a great extent of the gene pool in the modern southern East Asians. We also have observed genetic substructure among Hmong-Mien-speaking populations. The southern Hmong-Mien-speaking people showed more genetic affinity with modern Tai-Kadai/Austroasiatic people, while the northern Hmong-Mien speakers expressed a closer genetic connection with the Neolithic-to-modern northern East Asians. Moreover, southwestern Sinitic populations had a strong genomic affinity with the adjacent Hmong-Mien-speaking populations and the lowlander Tibeto-Burman-speaking populations, which suggested the large-scale genetic admixture occurred between them. Allele-sharing-based qpAdm/qpGraph results further confirmed that all included southwestern Chinese populations could be modeled as a mixed result of the major ancestry component from the northern millet farmers in the Yellow River basin and the minor ancestry component from the southern rice farmers in the Yangtze River basin. Usually, this newly identified Hmong-Mien-associated southern East Asian ancestry could improve our understanding of the full-scale genetic landscape of the evolutionary and admixture history of southwestern East Asians. Further ancient genomic studies from southeastern China are required to shed deeper light on our established phylogeny context.

The Genetic Structure and East-West Population Admixture in Northwest China Inferred From Genome-Wide Array Genotyping

Northwest China is a contacting region for East and West Eurasia and an important center for investigating the migration and admixture history of human populations. However, the comprehensive genetic structure and admixture history of the Altaic speaking populations and Hui group in Northwest China were still not fully characterized due to insufficient sampling and the lack of genome-wide data. Thus, We genotyped genome-wide SNPs for 140 individuals from five Chinese Mongolic, Turkic speaking groups including Dongxiang, Bonan, Yugur, and Salar, as well as the Hui group. Analysis based on allele-sharing and haplotype-sharing were used to elucidate the population history of Northwest Chinese populations, including PCA, ADMIXTURE, pairwise Fst genetic distance, f-statistics, qpWave/qpAdm and ALDER, fineSTRUCTURE and GLOBETROTTER. We observed Dongxiang, Bonan, Yugur, Salar, and Hui people were admixed populations deriving ancestry from both East and West Eurasians, with the proportions of West Eurasian related contributions ranging from 9 to 15%. The genetic admixture was probably driven by male-biased migration- showing a higher frequency of West Eurasian related Y chromosomal lineages than that of mtDNA detected in Northwest China. ALDER-based admixture and haplotype-based GLOBETROTTER showed this observed West Eurasian admixture signal was introduced into East Eurasia approximately 700 ∼1,000 years ago. Generally, our findings provided supporting evidence that the flourish transcontinental communication between East and West Eurasia played a vital role in the genetic formation of northwest Chinese populations.

ACE and ACTN3 Gene Polymorphisms and Genetic Traits of Rowing Athletes in the Northern Han Chinese Population

This investigation aimed to explore the effects of ACE I/D and ACTN3 R577X gene polymorphisms on specific quantitative variables, including height, weight, arm span, biacromial breadth, forced vital capacity (FVC), FVC/weight, maximal oxygen uptake (VO2max), prone bench pull (PBP), loaded barbell squat (LBS), and 3,000-m run, in 243 Chinese rowing athletes. The ACE and ACTN3 genotypes were obtained for each athlete via polymerase chain reaction on saliva samples, and the genotype frequency was analyzed. The ACE genotype frequency of rowing athletes were 45.8% II, 42.2% ID, and 12% DD for males and 33.6% II, 48% ID, and 18.4% DD for females. There were significant differences in weight in male athletes, PBP in female athletes, and ACE genotypes. A linear regression analysis using PBP and LBS as different dependent variables and ACE genotypes as independent variables based on the ACE I allele additive genetic effect showed a statistical significance in female athletes (p < 0.05). There was a significant difference in the distribution of the three genotypes among male athletes (36.7% XX, 38.5% RX, and 24.8% RR, χ2 = 5.191, df = 2, p = 0.022 < 0.05). There were no significant differences in the distribution of the three genotypes among female athletes (23.8% XX, 47.8% RX, 28.4% RR, χ 2 = 0.24, df = 2, p = 0.619 > 0.05). The ACTN3 gene polymorphism of male rowing athletes was dominated by the ACTN3 577X allele. There were significant differences in the χ 2 test between groups of male athletes. The ACTN3 R577 allele was dominant in female athletes. There were significant differences between PBP and FVC/body weight and ACTN3 genotypes in male athletes by ANOVA, respectively (p < 0.05). A linear regression analysis using FVC and FVC/body weight as dependent variables and ACTN3 genotypes as independent variables based on the ACTN3 577X allele recessive genetic effect showed statistical significance in male athletes (p < 0.05). These results suggested that ACE and ACTN3 gene polymorphisms may be used as biomarkers of genetic traits in Chinese rowing athletes.

Genomic Insight Into the Population Admixture History of Tungusic-Speaking Manchu People in Northeast China

Manchu is the third-largest ethnic minority in China and has the largest population size among the Tungusic-speaking groups. However, the genetic origin and admixture history of the Manchu people are far from clear due to the sparse sampling and a limited number of markers genotyped. Here, we provided the first batch of genome-wide data of genotyping approximate 700,000 single-nucleotide polymorphisms (SNPs) in 93 Manchu individuals collected from northeast China. We merged the newly generated data with data of publicly available modern and ancient East Asians to comprehensively characterize the genetic diversity and fine-scale population structure, as well as explore the genetic origin and admixture history of northern Chinese Manchus. We applied both descriptive methods of ADMIXTURE, fineSTRUCTURE, F ST , TreeMix, identity by decedent (IBD), principal component analysis (PCA), and qualitative f-statistics (f 3, f 4, qpAdm, and qpWave). We found that Liaoning Manchus have a close genetic relationship and significant admixture signal with northern Han Chinese, which is in line with the cluster patterns in the haplotype-based results. Additionally, the qpAdm-based admixture models showed that modern Manchu people were formed as major ancestry related to Yellow River farmers and minor ancestry linked to ancient populations from Amur River Bain, or others. In summary, the northeastern Chinese Manchu people in Liaoning were an exception to the coherent genetic structure of Tungusic-speaking populations, probably due to the large-scale population migrations and genetic admixtures in the past few hundred years.

Fine-Scale Genetic Structure and Natural Selection Signatures of Southwestern Hans Inferred From Patterns of Genome-Wide Allele, Haplotype, and Haplogroup Lineages

The evolutionary and admixture history of Han Chinese have been widely discussed via traditional autosomal and uniparental genetic markers [e.g., short tandem repeats, low-density single nucleotide polymorphisms). However, their fine-scale genetic landscapes (admixture scenarios and natural selection signatures) based on the high-density allele/haplotype sharing patterns have not been deeply characterized. Here, we collected and generated genome-wide data of 50 Han Chinese individuals from four populations in Guizhou Province, one of the most ethnolinguistically diverse regions, and merged it with over 3,000 publicly available modern and ancient Eurasians to describe the genetic origin and population admixture history of Guizhou Hans and their neighbors. PCA and ADMIXTURE results showed that the studied four populations were homogeneous and grouped closely to central East Asians. Genetic homogeneity within Guizhou populations was further confirmed via the observed strong genetic affinity with inland Hmong-Mien people through the observed genetic clade in Fst and outgroup f₃/f₄-statistics. qpGraph-based phylogenies and f₄-based demographic models illuminated that Guizhou Hans were well fitted via the admixture of ancient Yellow River Millet farmers related to Lajia people and southern Yangtze River farmers related to Hanben people. Further ChromoPainter-based chromosome painting profiles and GLOBETROTTER-based admixture signatures confirmed the two best source matches for southwestern Hans, respectively, from northern Shaanxi Hans and southern indigenes with variable mixture proportions in the historical period. Further three-way admixture models revealed larger genetic contributions from coastal southern East Asians into Guizhou Hans compared with the proposed inland ancient source from mainland Southeast Asia. We also identified candidate loci (e.g., MTUS2, NOTCH4, EDAR, ADH1B, and ABCG2) with strong natural selection signatures in Guizhou Hans via iHS, nSL, and ihh, which were associated with the susceptibility of the multiple complex diseases, morphology formation, alcohol and lipid metabolism. Generally, we provided a case and ideal strategy to reconstruct the detailed demographic evolutionary history of Guizhou Hans, which provided new insights into the fine-scale genomic formation of one ethnolinguistically specific targeted population from the comprehensive perspectives of the shared unlinked alleles, linked haplotypes, and paternal and maternal lineages.

Peopling History of the Tibetan Plateau and Multiple Waves of Admixture of Tibetans Inferred From Both Ancient and Modern Genome-Wide Data

Archeologically attested human occupation on the Tibetan Plateau (TP) can be traced back to 160 thousand years ago (kya) via the archaic Xiahe people and 30∼40 kya via the Nwya Devu anatomically modern human. However, the history of the Tibetan populations and their migration inferred from the ancient and modern DNA remains unclear. Here, we performed the first ancient and modern genomic meta-analysis among 3,017 Paleolithic to present-day Eastern Eurasian genomes (2,444 modern individuals from 183 populations and 573 ancient individuals). We identified a close genetic connection between the ancient-modern highland Tibetans and lowland island/coastal Neolithic Northern East Asians (NEA). This observed genetic affinity reflected the primary ancestry of high-altitude Tibeto-Burman speakers originated from the Neolithic farming populations in the Yellow River Basin. The identified pattern was consistent with the proposed common north-China origin hypothesis of the Sino-Tibetan languages and dispersal patterns of the northern millet farmers. We also observed the genetic differentiation between the highlanders and lowland NEAs. The former harbored more deeply diverged Hoabinhian/Onge-related ancestry and the latter possessed more Neolithic southern East Asian (SEA) or Siberian-related ancestry. Our reconstructed qpAdm and qpGraph models suggested the co-existence of Paleolithic and Neolithic ancestries in the Neolithic to modern East Asian highlanders. Additionally, we found that Tibetans from Ü-Tsang/Ando/Kham regions showed a strong population stratification consistent with their cultural background and geographic terrain. Ü-Tsang Tibetans possessed a stronger Chokhopani-affinity, Ando Tibetans had more Western Eurasian related ancestry and Kham Tibetans harbored greater Neolithic southern EA ancestry. Generally, ancient and modern genomes documented multiple waves of human migrations in the TP's past. The first layer of local hunter-gatherers mixed with incoming millet farmers and arose the Chokhopani-associated Proto-Tibetan-Burman highlanders, which further respectively mixed with additional genetic contributors from the western Eurasian Steppe, Yellow River and Yangtze River and finally gave rise to the modern Ando, Ü-Tsang and Kham Tibetans.

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

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

Significant East Asian Affinity of the Sichuan Hui Genomic Structure Suggests the Predominance of the Cultural Diffusion Model in the Genetic Formation Process

The ancestral origin and genomic history of Chinese Hui people remain to be explored due to the paucity of genome-wide data. Some evidence argues that an eastward migration of Central Asians gave rise to modern Hui people, which is referred to as the demic diffusion hypothesis; other evidence favors the cultural diffusion hypothesis, which posits that East Asians adopted Muslim culture to form the modern culturally distinct populations. However, the extent to which the observed genetic structure of the Huis was mediated by the movement of people or the assimilation of Muslim culture also remains highly contentious. Analyses of over 700 K SNPs in 109 western Chinese individuals (49 Sichuan Huis and 60 geographically close Nanchong Hans) together with the available ancient and modern Eurasian sequences allowed us to fully explore the genomic makeup and origin of Hui and neighboring Han populations. The results from PCA, ADMIXTURE, and allele-sharing-based f-statistics revealed a strong genomic affinity between Sichuan Huis and Neolithic-to-modern Northern East Asians, which suggested a massive gene influx from East Asians into the Sichuan Hui people. Three-way admixture models in the qpWave/qpAdm analyses further revealed a small stream of gene influx from western Eurasians into the Sichuan Hui people, which was further directly confirmed via the admixture event from the temporally distinct Western sources to Sichuan Hui people in the qpGraph-based phylogenetic model, suggesting the key role of the cultural diffusion model in the genetic formation of the Sichuan Huis. ALDER-based admixture date estimation showed that this observed western Eurasian admixture signal was introduced into the Sichuan Huis during the historic periods, which was concordant with the extensive western-eastern communication along the Silk Road and historically documented Huis' migration history. In summary, although significant cultural differentiation exists between Hui people and their neighbors, our genomic analysis showed their strong genetic affinity with modern and ancient Northern East Asians. Our results support the hypothesis that the Sichuan Huis arose from a mixture of minor western Eurasian ancestry and predominant East Asian ancestry.