Shared paternal ancestry of Han, Tai-Kadai-speaking, and Austronesian-speaking populations as revealed by the high resolution phylogeny of O1a-M119 and distribution of its sub-lineages within China

Development and validation of YARN: A novel SE-400 MPS kit for East Asian paternal lineage analysis

Y-chromosomal short tandem repeat polymorphisms (Y-STRs) and Y-chromosomal single nucleotide polymorphisms (Y-SNPs) are valuable genetic markers used in paternal lineage identification and population genetics. Currently, there is a lack of an effective panel that integrates Y-STRs and Y-SNPs for studying paternal lineages, particularly in East Asian populations. Hence, we developed a novel Y-chromosomal targeted panel called YARN (Y-chromosome Ancestry and Region Network) based on multiplex PCR and a single-end 400 massive parallel sequencing (MPS) strategy, consisting of 44 patrilineage Y-STRs and 260 evolutionary Y-SNPs. A total of 386 reactions were validated for the effectiveness and applicability of YARN according to SWGDAM validation guidelines, including sensitivity (with a minimum input gDNA of 0.125 ng), mixture identification (ranging from 1:1-1:10), PCR inhibitor testing (using substances such as 50 μM hematin, 100 μM hemoglobin, 100 μM humic acid, and 2.5 mM indigo dye), species specificity (successfully distinguishing humans from other animals), repeatability study (achieved 100% accuracy), and concordance study (with 99.91% accuracy for 1121 Y-STR alleles). Furthermore, we conducted a pilot study using YARN in a cohort of 484 Han Chinese males from Huaiji County, Zhaoqing City, Guangdong, China (GDZQHJ cohort). In this cohort, we identified 52 different Y-haplogroups and 73 different surnames. We found weak to moderate correlations between the Y-haplogroups, Chinese surnames, and geographical locations of the GDZQHJ cohort (with λ values ranging from 0.050 to 0.340). However, when we combined two different categories into a new independent variable, we observed stronger correlations (with λ values ranging from 0.617 to 0.754). Overall, the YARN panel, which combines Y-STR and Y-SNP genetic markers, meets forensic DNA quality assurance guidelines and holds potential for East Asian geographical origin inference and paternal lineage analysis.

Investigating demic versus cultural diffusion and sex bias in the spread of Austronesian languages in Vietnam

Austronesian (AN) is the second-largest language family in the world, particularly widespread in Island Southeast Asia (ISEA) and Oceania. In Mainland Southeast Asia (MSEA), groups speaking these languages are concentrated in the highlands of Vietnam. However, our knowledge of the spread of AN-speaking populations in MSEA remains limited; in particular, it is not clear if AN languages were spread by demic or cultural diffusion. In this study, we present and analyze new data consisting of complete mitogenomes from 369 individuals and 847 Y-chromosomal single nucleotide polymorphisms (SNPs) from 170 individuals from all five Vietnamese Austronesian groups (VN-AN) and five neighboring Vietnamese Austroasiatic groups (VN-AA). We found genetic signals consistent with matrilocality in some, but not all, of the VN-AN groups. Population affinity analyses indicated connections between the AN-speaking Giarai and certain Taiwanese AN groups (Rukai, Paiwan, and Bunun). However, overall, there were closer genetic affinities between VN-AN groups and neighboring VN-AA groups, suggesting language shifts. Our study provides insights into the genetic structure of AN-speaking communities in MSEA, characterized by some contact with Taiwan and language shift in neighboring groups, indicating that the expansion of AN speakers in MSEA was a combination of cultural and demic diffusion.

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.

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.

The formation of proto-austronesians: insights from a revised phylogeography of the paternal founder lineage

OBJECTIVES

Previous studies suggested that the Y-chromosome haplogroups O2-N6-B451-AM01756 and O1a-M119 are two founder lineages of proto-Austronesians at about five thousand years ago. The objective of this study was to investigate the formation of proto-Austronesians from the perspective of the paternal gene pool.

MATERIALS AND METHODS

In this study, we developed a highly evised phylogenetic tree with age estimates for haplogroup O2-N6 and early branches of O1a-M119 (M110, F1036, and F819). In addition, we also explored the geographical distribution of eight sub-branches of O2-N6 and O1a-M119, and spatial autocorrelation analysis was conducted for each sub-branch.

RESULTS

The paternal lineage combination of proto-Austronesians is a small subset of a diverse gene pool of populations from the mainland of East Asia. The distribution map and results of the spatial autocorrelation analysis suggested that the eastern coastal region of northern China is likely the source of lineage O2-N6 while the coastal region of southeastern China is likely the diffusion center of early branches of O1a-M119. We developed an evolutionary diagram for Austronesians and their ancestors in the past 18,000 years.

DISCUSSION

We proposed that the millet farming community in North China is the common ancestor group of the Austronesians and the Han people, while the diverse ancient people in the southeast coastal areas of East Asia form the common ancestor group of the Austronesians and the East Asian mainland population. The demographic history of multiple ancestral groups of the most recent common ancestor group itself in the more ancient period is helpful to understand the deep roots of the genetic components and cultural traditions of Austronesians.

Comprehensive insights into the genetic background of Chinese populations using Y chromosome markers

China is located in East Asia. With a high genetic and cultural diversity, human migration in China has always been a hot topic of genetics research. To explore the origins and migration routes of Chinese males, 3333 Chinese individuals (Han, Hui, Mongolia, Yi and Kyrgyz) with 27 Y-STRs and 143 Y-SNPs from published literature were analysed. Our data showed that there are five dominant haplogroups (O2-M122, O1-F265, C-M130, N-M231, R-M207) in China. Combining analysis of haplogroup frequencies, geographical positions and time with the most recent common ancestor (TMRCA), we found that haplogroups C-M130, N-M231 and R1-M173 and O1a-M175 probably migrated into China via the northern route. Interestingly, we found that haplogroup C*-M130 in China may originate in South Asia, whereas the major subbranches C2a-L1373 and C2b-F1067 migrated from northern China. The results of BATWING showed that the common ancestry of Y haplogroup in China can be traced back to 17 000 years ago, which was concurrent with global temperature increases after the Last Glacial Maximum.

Unveiling 2,000 years of differentiation among Tungusic-speaking populations: a revised phylogeny of the paternal founder lineage C2a-M48-SK1061

Previous studies demonstrated Y chromosome haplogroup C2a-M48-SK1061 is the only founding paternal lineage of all Tungusic-speaking populations. To infer the differentiation history of these populations, we studied more sequences and constructed downstream structure of haplogroup C2a-M48-SK1061 with better resolution. In this study, we generated 100 new sequences and co-analyzed 140 sequences of C2a-M48-SK1061 to reconstruct a highly revised phylogenetic tree with age estimates. We also performed the analysis of the geographical distribution and spatial autocorrelation of sub-branches. Dozens of new sub-branches were discovered, many sub-branches were nearly unique for Ewenki, Evens, Oroqen, Xibe, Manchu, Daur, and Mongolian. The topology of these unique sub-branches is the key evidence for understanding the complex evolutionary relationship between different Tungusic-speaking populations. The revised phylogeny provided a clear pattern for the differentiation history of haplogroup C2a-M48-SK1061 in the past 2,000 years. This study showed that the divergence pattern of founder lineage is essential to understanding the differentiation history of populations.

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.

Multiple founding paternal lineages inferred from the newly-developed 639-plex Y-SNP panel suggested the complex admixture and migration history of Chinese people

BACKGROUND

Non-recombining regions of the Y-chromosome recorded the evolutionary traces of male human populations and are inherited haplotype-dependently and male-specifically. Recent whole Y-chromosome sequencing studies have identified previously unrecognized population divergence, expansion and admixture processes, which promotes a better understanding and application of the observed patterns of Y-chromosome genetic diversity.

RESULTS

Here, we developed one highest-resolution Y-chromosome single nucleotide polymorphism (Y-SNP) panel targeted for uniparental genealogy reconstruction and paternal biogeographical ancestry inference, which included 639 phylogenetically informative SNPs. We genotyped these loci in 1033 Chinese male individuals from 33 ethnolinguistically diverse populations and identified 256 terminal Y-chromosomal lineages with frequency ranging from 0.0010 (singleton) to 0.0687. We identified six dominant common founding lineages associated with different ethnolinguistic backgrounds, which included O2a2b1a1a1a1a1a1a1-M6539, O2a1b1a1a1a1a1a1-F17, O2a2b1a1a1a1a1b1a1b-MF15397, O2a2b2a1b1-A16609, O1b1a1a1a1b2a1a1-F2517, and O2a2b1a1a1a1a1a1-F155. The AMOVA and nucleotide diversity estimates revealed considerable differences and high genetic diversity among ethnolinguistically different populations. We constructed one representative phylogenetic tree among 33 studied populations based on the haplogroup frequency spectrum and sequence variations. Clustering patterns in principal component analysis and multidimensional scaling results showed a genetic differentiation between Tai-Kadai-speaking Li, Mongolic-speaking Mongolian, and other Sinitic-speaking Han Chinese populations. Phylogenetic topology inferred from the BEAST and Network relationships reconstructed from the popART further showed the founding lineages from culturally/linguistically diverse populations, such as C2a/C2b was dominant in Mongolian people and O1a/O1b was dominant in island Li people. We also identified many lineages shared by more than two ethnolinguistically different populations with a high proportion, suggesting their extensive admixture and migration history.

CONCLUSIONS

Our findings indicated that our developed high-resolution Y-SNP panel included major dominant Y-lineages of Chinese populations from different ethnic groups and geographical regions, which can be used as the primary and powerful tool for forensic practice. We should emphasize the necessity and importance of whole sequencing of more ethnolinguistically different populations, which can help identify more unrecognized population-specific variations for the promotion of Y-chromosome-based forensic applications.

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.

Ancient Components and Recent Expansion in the Eurasian Heartland: Insights into the Revised Phylogeny of Y-Chromosomes from Central Asia

In the past two decades, studies of Y chromosomal single nucleotide polymorphisms (Y-SNPs) and short tandem repeats (Y-STRs) have shed light on the demographic history of Central Asia, the heartland of Eurasia. However, complex patterns of migration and admixture have complicated population genetic studies in Central Asia. Here, we sequenced and analyzed the Y-chromosomes of 187 male individuals from Kazakh, Kyrgyz, Uzbek, Karakalpak, Hazara, Karluk, Tajik, Uyghur, Dungan, and Turkmen populations. High diversity and admixture from peripheral areas of Eurasia were observed among the paternal gene pool of these populations. This general pattern can be largely attributed to the activities of ancient people in four periods, including the Neolithic farmers, Indo-Europeans, Turks, and Mongols. Most importantly, we detected the consistent expansion of many minor lineages over the past thousand years, which may correspond directly to the formation of modern populations in these regions. The newly discovered sub-lineages and variants provide a basis for further studies of the contributions of minor lineages to the formation of modern populations in Central Asia.

Genetic substructure of Guizhou Tai-Kadai-speaking people inferred from genome-wide single nucleotide polymorphisms data

The genome-wide characteristics and admixture history of the Tai-Kadai-speaking populations are essential for understanding the population genetic diversity in southern China. We genotyped about 700,000 single nucleotide polymorphisms (SNPs) of 239 individuals from six Tai-Kadai-speaking populations residing in the mountainous Guizhou Province of southwestern China. We merged the genome-wide data with available populations and ancients in East and Southeast Asia to infer Tai-Kadai-speaking populations’ admixture history and genetic structure. We observed a genetic substructure within the studied six populations in the PCA, ADMIXTURE, ChromoPainter, GLOBETROTTER, f-statistics, and qpWave analysis. The Dong, Zhuang, and Bouyei people had a strong genetic affinity with other Tai-Kadai-speaking and Austronesian groups in the surrounding area. However, Gelao showed an affinity to Sino-Tibetan groups, and Mulao people were genetically close to Hmong-Mien populations. qpAdm further illuminated that Gelao and Dong_Tongren composited more Han-related ancestry than Dong, Zhuang, Bouyei, and Mulao people. Meanwhile, we observed high frequencies of Y-chromosome haplogroup O in studied Tai-Kadai-speaking groups except for Gelao people with a high haplogroup N frequency. From the maternal side, haplogroup M7 was frequent in studied populations except for Tongren Dong, who had a high frequency of haplogroup B5. Our newly reported data are helpful for further exploring population dynamics in southern China.

Whole mitochondrial genome analysis of Tai-Kadai-speaking populations in Southwest China

As a single matrilineal gene, human mitochondrial DNA plays a very important role in the study of population genetics. The whole mitogenome sequences of 287 individuals of the Tai-Kadai-speaking population in Guizhou were obtained. It was discovered that there were 82, 104, and 94 haplotypes in 83 Bouyei individuals, 107 Dong individuals, and 97 Sui individuals, respectively; and the haplotype diversity in Bouyei, Dong, and Sui groups was 1.000 ± 0.02, 0.9993 ± 0.0015, and 0.999 ± 0.002, respectively. The result of neutrality tests of the Tai-Kadai-speaking population in Guizhou showed significant negative values, and the analysis of mismatch distribution showed an obvious unimodal distribution. The results implied that Guizhou Tai-Kadai-speaking populations had high genetic diversities and may have experienced recent population expansion. In addition, the primary haplogroups of studied populations were M*, F, B, D, and R*, implying that they may origin from Southern China. The matrilineal genetic structure of the Tai-Kadai-speaking populations in Guizhou was analyzed by merging the mitogenome data of 79 worldwide populations as reference data. The results showed that there were close relationships between studied populations and other Tai-Kadai as well as some Austronesian populations in East and Southeast Asia. Overall, the mitogenome data generated in this study will provide important data for the study of genetic structure of Tai-Kadai speaking populations.

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.

Episodes of diversification and isolation in Island Southeast Asian and Near Oceanian male lineages

Island Southeast Asia (ISEA) and Oceania host one of the world’s richest assemblages of human phenotypic, linguistic, and cultural diversity. Despite this, the region’s male genetic lineages are globally among the last to remain unresolved. We compiled ∼9.7 Mb of Y chromosome (chrY) sequence from a diverse sample of over 380 men from this region, including 152 first reported here. The granularity of this data set allows us to fully resolve and date the regional chrY phylogeny. This new high-resolution tree confirms two main population bursts: multiple rapid diversifications following the region’s initial settlement ∼50 kya, and extensive expansions <6 kya. Notably, ∼40–25 kya the deep rooting local lineages of C-M130, M-P256, and S-B254 show almost no further branching events in ISEA, New Guinea, and Australia, matching a similar pause in diversification seen in maternal mitochondrial DNA lineages. The main local lineages start diversifying ∼25 kya, at the time of the last glacial maximum. This improved chrY topology highlights localized events with important historical implications, including pre-Holocene contact between Mainland and ISEA, potential interactions between Australia and the Papuan world, and a sustained period of diversification following the flooding of the ancient Sunda and Sahul continents as the insular landscape observed today formed. The high-resolution phylogeny of the chrY presented here thus enables a detailed exploration of past isolation, interaction, and change in one of the world’s least understood regions.

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.