Paternal origin of Tungusic‐speaking populations: Insights from the updated phylogenetic tree of Y‐chromosome haplogroup C2a‐M86

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.

Genetic polymorphism of Y-chromosome in Kazakh populations from Southern Kazakhstan

BACKGROUND

The Kazakhs are one of the biggest Turkic-speaking ethnic groups, controlling vast swaths of land from the Altai to the Caspian Sea. In terms of area, Kazakhstan is ranked ninth in the world. Northern, Eastern, and Western Kazakhstan have already been studied in relation to genetic polymorphism 27 Y-STR. However, current information on the genetic polymorphism of the Y-chromosome of Southern Kazakhstan is limited only by 17 Y-STR and no geographical study of other regions has been studied at this variation.

RESULTS

The Kazakhstan Y-chromosome Haplotype Reference Database was expanded with 468 Kazakh males from the Zhambyl and Turkestan regions of South Kazakhstan by having their 27 Y-STR loci and 23 Y-SNP markers analyzed. Discrimination capacity (DC = 91.23%), haplotype match probability (HPM = 0.0029) and haplotype diversity (HD = 0.9992) are defined. Most of this Y-chromosome variability is attributed to haplogroups C2a1a1b1-F1756 (2.1%), C2a1a2-M48 (7.3%), C2a1a3-F1918 (33.3%) and C2b1a1a1a-M407 (6%). Median-joining network analysis was applied to understand the relationship between the haplotypes of the three regions. In three genetic layer can be described the position of the populations of the Southern region of Kazakhstan-the geographic Kazakh populations of Kazakhstan, the Kazakh tribal groups, and the people of bordering Asia.

CONCLUSION

The Kazakhstan Y-chromosome Haplotype Reference Database was formed for 27 Y-STR loci with a total sample of 1796 samples of Kazakhs from 16 regions of Kazakhstan. The variability of the Y-chromosome of the Kazakhs in a geographical context can be divided into four main clusters-south, north, east, west. At the same time, in the genetic space of tribal groups, the population of southern Kazakhs clusters with tribes from the same region, and genetic proximity is determined with the populations of the Hazaras of Afghanistan and the Mongols of China.

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.

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.

Paternal genetic structure of the Qiang ethnic group in China revealed by high-resolution Y-chromosome STRs and SNPs

The Qiang population mainly lived in Beichuan Qiang Autonomous County of Sichuan Province. It is one of the nomads in China, distributed along the Minjiang River. The Qiang population was assumed to have great affinity with the Han, the largest ethnic group in China, when it refers to the genetic origin. Whereas, it is deeply understudied, especially from the Y chromosome. Here in this study, we used validated high-resolution Y-chromosome single nucleotide polymorphisms (Y-SNPs) and short tandem repeats (Y-STRs) panels to study the Qiang ethnic group to unravel their paternal genetic, forensic and phylogenetic characteristics. A total of 422 male samples of the Qiang ethnic group were genotyped by 233 Y-SNPs and 29 Y-STRs. Haplogroup O-M175 (N = 312) was the most predominant haplogroup in the Qiang ethnic group, followed by D-M174 (N = 32) and C-M130 (N = 32), N-M231 (N = 27), and Q-M242 (N = 15). After further subdivision, O2a-M324 (N = 213) accounted for the majority of haplogroup O. Haplogroup C2b-Z1338 (N = 29), D1a-CTS11577 (N = 30). O2a2b1a1a1-F42 (N = 48), O2a1b1a1a1a-F11 (N = 35), and O2a2b1a1-M117 (N = 21) represented other large terminal haplogroups. The results unveiled that Qiang ethnic group was a population with a high percentage of haplogroup O2a2b1a1a1-F42 (48/422) and O2a1b1a1a1a-F11 (35/422), and O2a2b1a1-M117 (21/422), which has never been reported. Its haplogroup distribution pattern was different from any of the Han populations, implying that the Qiang ethnic group had its unique genetic pattern. Mismatch analysis indicated that the biggest mismatch number in haplogroup O2a2b1a1a1-F42 was 21, while that of haplogroup O2a1b1a1a1a-F11 was 20. The haplotype diversity of the Qiang ethnic group equaled 0.999788, with 392 haplotypes observed, of which 367 haplotypes were unique. The haplogroup diversity of the Qiang ethnic group reached 0.9767, and 53 terminal haplogroups were observed (The haplogroup diversity of the Qiang ethnic group was the highest among Qiang and all Han subgroups, indicating the larger genetic diversity of the Qiang ethnic group.). Haplogroup O2a2b1a1a1-F42 was the most predominant haplogroup, including 11.37 % of the Qiang individuals. Median-joining trees showed gene flow between the Qiang and Han individuals. Our results indicated that 1) the highest genetic diversity was observed in the Qiang ethnic group compared to any of the former studied Chinese population, suggesting that the Qiang might be an older paternal branch; 2) the haplogroup D-M174 individuals of Qiang, Tibetans and Japanese distributed in three different subclades, which was unable to identify through low-resolution Y-SNP panel; and 3) the Qiang had lower proportion of haplogroup D compared to Yi and Tibetan ethnic groups, showing that the Qiang had less genetic communication with them than with Han Chinese.

Whole mitochondrial genome analysis of the Daur ethnic minority from Hulunbuir in the Inner Mongolia Autonomous Region of China

BACKGROUND

Mitochondrial DNA (mtDNA) variations are often associated with bioenergetics, disease, and speciation and can be used to track the history of women. Although advances in massively parallel sequencing (MPS) technology have greatly promoted our understanding of the population's history (especially genome-wide data and whole Y chromosome sequencing), the whole mtDNA sequence of many important groups has not been fully studied. In this study, we employed whole mitogenomes of 209 healthy and unrelated individuals from the Daur group, a Mongolic-speaking representative population of the indigenous groups in the Heilongjiang River basin (also known as the Amur River basin).

RESULTS

The dataset presented 127 distinct mtDNA haplotypes, resulting in a haplotype diversity of 0.9933. Most of haplotypes were assigned to eastern Eurasian-specific lineages, such as D4 (19.62%), B4 (9.09%), D5 (7.66%) and M7 (4.78%). Population comparisons showed that the Daurians do have certain connections with the ancient populations in the Heilongjiang River basin but the matrilineal genetic composition of the Daur group was also greatly influenced by other non-Mongolic groups from neighboring areas.

CONCLUSIONS

Collectively, the whole mtDNA data generated in the present study will augment the existing mtDNA database. Our study provides genetic links between the Daur population and the aborigine peoples from Siberia and the Amur-Ussuri Region. But on the whole, compared with other Mongolic-speaking groups, the modern Daur population is closer to the East Asian ancestry group.

Genetic insights into the paternal admixture history of Chinese Mongolians via high-resolution customized Y-SNP SNaPshot panels

The Mongolian people, one of the Mongolic-speaking populations, are native to the Mongolian Plateau in North China and southern Siberia. Many ancient DNA studies recently reported extensive population transformations during the Paleolithic to historic periods in this region, while little is known about the paternal genetic legacy of modern geographically different Mongolians. Here, we genotyped 215 Y-chromosomal single nucleotide polymorphisms (Y-SNPs) and 37 Y-chromosomal short tandem repeats (Y-STRs) among 679 Mongolian individuals from Hohhot, Hulunbuir, and Ordos in North China using the AGCU Y37 kit and our developed eight Y-SNP SNaPshot panels (including two panels first reported herein). The C-M130 Y-SNP SNaPshot panel defines 28 subhaplogroups, and the N/O/Q complementary Y-SNP SNaPshot panel defines 30 subhaplogroups of N1b-F2930, N1a1a1a1a3-B197, Q-M242, and O2a2b1a1a1a4a-CTS4658, which improved the resolution our developed Y-SNP SNaPshot panel set and could be applied for dissecting the finer-scale paternal lineages of Mongolic speakers. We found a strong association between Mongolian-prevailing haplogroups and some observed microvariants among the newly generated Y-STR haplotype data, suggesting the possibility of haplogroup prediction based on the distribution of Y-STR haplotypes. We identified three main ancestral sources of the observed Mongolian-dominant haplogroups, including the local lineage of C2*-M217 and incoming lineages from other regions of southern East Asia (O2*-M122, O1b*-P31, and N1*-CTS3750) and western Eurasia (R1*-M173). We also observed DE-M145, D1*-M174, C1*-F3393, G*-M201, I-M170, J*-M304, L-M20, O1a*-M119, and Q*-M242 at relatively low frequencies (< 5.00%), suggesting a complex admixture history between Mongolians and other incoming Eurasians from surrounding regions. Genetic clustering analyses indicated that the studied Mongolians showed close genetic affinities with other Altaic-speaking populations and Sinitic-speaking Hui people. The Y-SNP haplotype/haplogroup-based genetic legacy not only revealed that the stratification among geographically/linguistically/ethnically different Chinese populations was highly consistent with the geographical division and language classification, but also demonstrated that patrilineal genetic materials could provide fine-scale genetic structures among geographically different Mongolian people, suggesting that our developed high-resolution Y-SNP SNaPshot panels have the potential for forensic pedigree searches and biogeographical ancestry inference.

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.

Medieval Super-Grandfather founder of Western Kazakh Clans from Haplogroup C2a1a2-M48

Western Kazakhstan is populated by three clans totaling 2 million people. Since the clans are patrilineal, the Y-chromosome is the most informative genetic system for tracing their origin. We genotyped 40 Y-SNP and 17 Y-STR markers in 330 Western Kazakhs. High phylogenetic resolution within haplogroup C2a1a2-M48 was achieved by using additional SNPs. Three lines of evidence indicate that the Alimuly and Baiuly clans (but not the Zhetiru clan) have a common founder placed 700 ± 200 years back by the STR data and 500 ± 200 years back by the sequencing data. This supports traditional genealogy claims about the descent of these clans from Emir Alau, who lived 650 years ago and whose lineage might be carried by two-thirds of Western Kazakhs. There is accumulation of specific haplogroups in the subclans representing other lineages, confirming that the clan structure corresponds with the paternal genetic structure of the steppe population.