Hui-Xin Yu1†
Cheliger Ao1†
Xiao-Peng Wang2†
Xian-Peng Zhang1Jin Sun3
Hui Li4,5Kai-Jun Liu6*
Lan-Hai Wei1,5*- 1School of Ethnology and Anthropology, Institute of Humanities and Human Sciences, Inner Mongolia Normal University, Hohhot, China
- 2School of Management, Dalian University of Technology, Dalian, China
- 3School of Literature and Media, Xingyi Normal University for Nationalities, Xingyi, China
- 4MOE Key Laboratory of Contemporary Anthropology, School of Life Sciences, Fudan University, Shanghai, China
- 5B&R International Joint Laboratory for Eurasian Anthropology, Fudan University, Shanghai, China
- 6School of International Tourism and culture, Guizhou Normal University, Guiyang, China
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.
Introduction
Bronze Age globalization led to large-scale archaeological and cultural changes in much of Eurasia (de Barros Damgaard et al., 2018; Jeong et al., 2018; Jeong et al., 2019). The widespread spread of Indo-Europeans and related populations has been accompanied by strong replacement and mixing of population in many regions (Allentoft et al., 2015), such as South Siberia (Yu et al., 2020), Central Asia (Ning et al., 2019; Zhang et al., 2021; Kumar et al., 2022), South Asia (Narasimhan et al., 2019), the Middle East (Lazaridis et al., 2016), and Europe (Haak et al., 2015; Damgaard et al., 2018; Mathieson et al., 2018; Olalde et al., 2018). However, the Bronze Age globalization does not necessarily lead to the large-scale replacement of the population. We proposed that the subsistence strategy may be the key element. In northern Eurasia, Middle East, and East Asia, the natural environment may be not conducive to the original subsistence strategy, a. ka. The typical nomadic lifestyle, of Indo-Europeans. If the original subsistence strategy of Indo-Europeans can be firmly established in a certain region after the initial diffusion, significant population replacement may occur in this region and the local population may revive again.
In the central part of East Asia, during the Chalcolithic Age and the Bronze Age, early states and early civilizations arose (Liu and Chen, 2017). However, studies of bioarchaeology and ancient DNA have so far not observed obvious Indo-European-related genetic components in the central part of East Asia (Zhao et al., 2014; Ning et al., 2020; Yang et al., 2020; Mao et al., 2021; Robbeets et al., 2021; Wang et al., 2021; Xue et al., 2022). The demographic context of the relevant historical processes is unclear.
Previous studies have provided clues to population exchanges between the eastern Eurasian steppe region and the central part of East Asia. The Bronze Age populations of the Mongolian Plateau (∼4000–2500 YBP) is mainly a mixture of northern East Asian, Ancient North Eurasians, Southern Siberian indigenous and Indo-European related populations (Jeong et al., 2018; Jeong et al., 2020; Robbeets et al., 2021; Wang et al., 2021). The genetic contribution of these populations to the population of ancient central China was limited. Ancient DNA studies show that the genetic component of the Bronze Age populations in northwest China are basically local types, and there are no obvious foreign branches, like what had observed in Shimao (Xue et al., 2022), Jinchankou (Ning et al., 2020), and Xiajiadian sites (Li et al., 2011). On the other hand, archaeological studies confirm that the Mongolian Plateau is the intermediary zone where the bronze culture of central Eurasia spread to the central part of East Asia (Новгородова, 1989; Liu, 2004; Cai et al., 2007; Houle and Board., 2014; Liu and Chen, 2017). The important question raised by this is, is the spread of bronze culture to the Central Plains purely cultural transmission? What ancient people’s activities led to the emergence and prosperity of the bronze culture in northwest China? The relevant ancient DNA evidence is limit up till now.
Studies on the paternal Y chromosome provide some clues to the above questions. There is a high diversity of paternal components in the Chinese populations, and dozens of major haplogroups have been found (Zhong et al., 2010; Yan et al., 2014; Lu et al., 2020). Most of them originated from the central region of East Asia (Zhong et al., 2011). However, the demographic history of two of them, Q1a1a-M120 and N1a2a-F1101, shows links to Siberia and the eastern Eurasian steppe (Jeong et al., 2018; Sun et al., 2019; Jeong et al., 2020). Haplogroup Q-M242 has dozens of downstream clades, and only one Q1a1a-M120 clade appears frequently in East Asian populations (Huang et al., 2018; Wei et al., 2018; Sun et al., 2019). Ancient DNA studies suggest that haplogroup Q1a1a-M120 is the main paternal lineage of Bronze Age populations in northwest China (Zhao et al., 2014; Sun et al., 2019) and this lineage also appeared on the Mongolia Plateau (Jeong et al., 2018; Jeong et al., 2020). Previously, we had published a research to explore the migration from southern Siberia to East Asia and the role Q1a1a-M120 in the rise of the Bronze Culture populations in northwest China (Sun et al., 2019).
Previously, ancient DNA studies showed that N1a2a-F1101 was likely the paternal lineage of the royal family of the third dynasty of ancient China, the Zhou Dynasty (1027–256 BC) (Ma et al., 2021; Wei et al., 2022). The Zhou dynasty was the last dynasty of the Bronze Age in central China. Interestingly, N1a2b-P43, the mostly close branch of N1a2a-F1101, is the founder paternal lineage of the Uralic populations in northern Eurasia (Rootsi et al., 2007; Karmin et al., 2015; Ilumäe et al., 2016). The initial divergence of these two clades should have been between modern northern China and Siberia. Haplogroup N1a2a-F1101 is one of the major paternal lineage of modern Chinese, especially the Han Chinese (Liu et al., 2022; Tao et al., 2023). However, it remains unclear that how haplogroup N1a2a-F1101, a lineage likely originated in the north boundary region of China, become an important paternal component among ancient and modern populations in central China.
In this study, we analyzed the sequences and distribution of haplogroup N1a2a-F1101 within China. Our first objective was to construct a detailed phylogenetic tree for this paternal lineage with age estimation. Second, we analyzed the demographic history of this lineage and its role during the formation of early state and ancient civilization in ancient China. Specifically, we discussed the special demographic history in East Asia during the Bronze age globalization.
Materials and methods
Samples and sequencing
Saliva samples were collected from unrelated healthy males in East Asian populations over the past few decades. All participants provided written informed consent prior to participating. The study and sample collection process were reviewed and approved by the Medical Ethics Committee of Fudan University and Inner Mongolian Normal University, and complied with the ethical principles of the 2013 Helsinki Declaration of the World Medical Association. Genomic DNA was extracted using the DP-318 Kit (Tiangen Biotechnology, Beijing, China) according to the manufacturer’s protocol. A series of SNaPshot panels were used to determine the downstream Y-SNP (Single Nucleotide Polymorphism) haplgroup (Cai et al., 2011; Hu et al., 2015). As indicated by previous study, the major paternal Y-chromosome haplogroup of East Asia populations (like C, D, N, O. Q, and R) and downstream sub-branches can be identified by these panels (Cai et al., 2011; Hu et al., 2015). After that, we can select samples belonging to a certain haplogroup for further research. Extracted DNA of 229 individuals of haplogroup N1a2a-F1101 were sent for next-generation sequencing of the Y-chromosome using the Illumina HiSeq2000 platform (San Diego, CA, USA). All of these 229 individuals came from Chinese populations and most of them are Han Chinese (Supplementary Table S1). All of samples included in this study are healthy unrelated male individuals (Supplementary Table S2).
Data availability statement
Following the regulations of the Human Genetic Resources Administration of China (HGRAC), the raw sequence data reported in this paper are available on request from the corresponding author. A list of variants of Y-SNP analyzed in this study was included in Supplementary Table S2.
Data analysis
In total, 229 sequences of N1a2a-F1101 were analyzed (Supplementary Table S1). Read mapping and SNP (single nucleotide polymorphism) calling from next-generation sequencing data were conducted using standard procedures (BWA and SAMtools) and the human reference genome sequence hg38 (Li and Durbin, 2009; 2010). Bayesian evolutionary analyses were conducted using BEAST (v.2.0.0) (Bouckaert et al., 2014). To calculate divergence times in the phylogenetic tree, a point mutation rate of 0.74 × 10-9 per site per year (Karmin et al., 2015), inferred from the ∼12,000-year-old Anzick-1 male infant genome (Rasmussen et al., 2014), was applied. We referred to the regulations established by the Y Chromosome Consortium for the assignment of SNP and revised haplogroup names (The Y Chromosome Consortium, 2002) (The Y Chromosome Consortium, 2002). New haplogroup names for sub-branches of N1a2a-F1101 are listed in Supplementary Table S2.
To analyze the distribution of N1a2a-F1101 sub-lineages, Chinese participants were drawn from our inhouse database, which contains various types of information for each sample, including the Y haplogroup, native province, and native city. All steps followed those described in our previous study (Sun et al., 2021). In total, 1,601 Q1a1a-M120 individuals were identified among 140,900 Chinese males (total frequency: 2.54%, Supplementary Table S3). Frequency data were used to generate distribution maps using ArcGIS (version 10.3; Environmental Systems Research Institute Inc. Redlands, CA, USA). Please refer to the Supplementary Text of our previous study (Sun et al., 2021) for details about the ArcGIS was also used for a spatial autocorrelation analysis. The HotSpots plot indicate the clustering region of hop spots and cold spots which are generally corresponding to the center of diffusion and its mirror image region shown by the general distribution pattern. The Moran’s I Index indicate the degree of correlation between distribution frequency and spatial geographical distribution. The Lisa Cluster show the cluster and outlier of general distribution pattern of high-low frequency.
Results
A revised phylogenetic tree of haplogroup N1a2a-F1101 were constructed with age estimation (Figure 1 and Supplementary Table S1). The haplogroups N1a2b-P43 and N1a2a-F1101 split at about 9300 years ago. There are similarities in the early history of the two haplogroups. They all experienced a very significant expansion after a bottleneck period of nearly 5,000 years and became the dominant paternal lineage of descendant populations. The main downstream branch of N1a2a-F1101 is N1a2a1-F1154, and the main differentiation node time is 4400 and 4000 years ago, and dozens of downstream branches are born. Among them, N1a2a1a1a1a1-F710 has undergone significant expansion after 3,350 years ago, giving birth to more than 70 downstream clades (Figure 1 and Supplementary Table S1). This topology suggests that the population expansion experienced by this paternal line around 3,000 years ago was the most significant of all paternal lineages in ancient East Asian populations at the same history period. Previously, ancient DNA studies suggested that this paternal line may be the paternal lineage of the Zhou Dynasty, the third dynasty of ancient China (Ma et al., 2021; Wei et al., 2022). The differentiation topology of this study supports the results of ancient DNA findings.
FIGURE 1. Schematic diagram of the revised phylogeny of haplogroup N1a2a-F1101. The red SNP labels indicated the definition marker of categories of sub-branches in this study. The green SNP labels indicated other important marker of the phylogenetic tree.
The proportion of early sub-branch N1a2a-F1101 clades in modern populations is very low and the distribution is discrete (Figure 2A, and Supplementary Table S2). Therefore, the possible diffusion centers shown by the HotSpots analysis in Figure 2A need to be treated with caution. Because sub-branch N1a2a1b-MF55705 is an early major branch of N1a2a-F1101 that emerged during the Chalcolithic Age (∼4400), we analyzed it separately. N1a2a1b-MF55705 is mainly distributed in eastern China (Figure 2B, and Supplementary Table S2). Results of spatial autocorrelation analysis suggests that the early spread center of this branch may have been the lower Yellow River region. Sub-branch N1a2a1a1-CTS2959 is the major downstream branch of N1a2a-F1101. High frequency of this lineage tend to found in northern China (Figure 2C, and Supplementary Table S2). HotSpots analysis in spatial autocorrelation analysis suggests that the early diffusion center of this clade may have been the central part of north China, consistent with the theory that this paternal line is the paternal lineage of the Zhou royal family.
FIGURE 2. Geographic distribution of sub-branches of N1a2a-F1101 and results of a spatial autocorrelation analysis of frequencies. (A): N1a2a*-F1101×(CTS2959,MF55705), (B): N1a2a1b-MF55705, (C): N1a2a1a1-CTS2959.
Discussion
Early history between 9.3 kya and 4.4 kya
As the only two downstream clades of N1a2-L666, the geographical distribution of N1a2a-F1101 and N1a2b-P43 is very different from each other. Ancient DNA studies have identified early branches of N1a2a-F1101 and N1a2b-P43 in sites in the Baikal region (de Barros Damgaard et al., 2018; Kilinc et al., 2021; Ma et al., 2021). The most recent branch of N1a2-L666 is N1a1-M46, the main paternal type of the Uralic population (Ilumäe et al., 2016). The first two early branches under N1a1-M46, N1a1b-Y149447 and N1a1a3-F4065, are mainly distributed in northeast China (https://www.yfull.com/tree/N/) (Hu et al., 2015). Therefore, we speculate that the initial spread of haplogroup N1a2-L666 may have been in the southwestern part of northeastern China (Figure 3). We proposed that this region is also the initial diffusion center of N1a1-M46, while the diffusion of N1a1-M46 (>12 kya) happened earlier than that of N1a2-L666 (<9.3 kya) (Hu et al., 2015). In the early Holocene (about 11.2kya-8kya), with climate change and the rise of early agricultural populations in northern China, a part of the descendants of the ancestor group, representing by sub-lineage N1a2b-P43, spread to the high latitude region of Siberia, eventually becoming part of the Ural-speaking populations. The other part, representing by sub-lineage N1a2a-F1101, remained in the local area and participated in the formation of the northern Chinese populations in the later historical period (Figure 3).
FIGURE 3. Possible diffusion pattern of N1a2a-F1101 in the past 9,000 years.
A bottleneck period of 5,000 years was observed early in the evolution of N1a2a-F1101 (Figure 1, Supplementary Table S1). Similar lengthy bottleneck periods were observed in downstream structures of N1a2b-P43, N1a1-M46, and Q1a1a-M120 (Ilumäe et al., 2016; Sun et al., 2019). This evolutionary pattern is very different from the expansion pattern of ancient agricultural populations in East Asia, which continued to expand since the beginning of Neolithic age (Yan et al., 2014). The differentiation of the downstream clades of Q-M242 and N-231 presents a similar structure, i.e., downstream clades with high frequency distribution both in East Asia and Siberia, respectively. Therefore, we speculate that in the bottleneck interval, ancient populations with Q1a1a-M120 and N1a2a-F1101 as the main paternal lineages are likely to exist in the form of prehistoric hunter-gatherer populations in the border between the eastern Eurasian steppe and the northern-northeastern China. The drought and harsh natural environment of this area had a great influence on the evolution of the two paternal lineages in later historical periods.
Expansion during the chalcolithic age and bronze age
During the Chalcolithic age (about 4.5 kya-4.0 kya) in East Asia, copper, cattle and wheat were introduced to the East Asian heartland (Liu and Chen, 2003; Liu, 2004; Liu and Chen, 2017). Archaeologists have suggested that the elements may have spread from northern boundary of China through the Eurasian steppe. However, the demographic context of this important cultural process is very ambiguous. Around 4,000 years ago, the Bronze culture arose in the agro-pastoral region of northwestern China and later spread across East Asia and Southeast Asia. The mixing of the bronze culture of agriculture and animal husbandry with the people of the middle and lower reaches of the Yellow River contributed to the establishment of three dynasties of the Bronze Age in ancient China, namely the Xia, Shang and Zhou dynasties (Liu and Chen, 2003; Liu, 2004; Liu and Chen, 2017).
As discussed above, ancient populations with Q1a1a-M120 and N1a2a-F1101 as the main paternal lineages may have played a mediating role in the spread of the Copper and Bronze cultures from the eastern Eurasian steppe to the central East Asian region, due to their area of activity in the junction zone. Due to the same reason, these two paternal lines experienced a very significant spread during the Bronze Age, becoming important patrilineal lineages that occupied an upper political position in the Bronze Age, and were frequently detected in the tombs of chiefs and nobles of the time (Zhao et al., 2014; Sun et al., 2019; Ma et al., 2021; Wei et al., 2022). An interesting thing is that the significant expansion of N1a2a-F1101 occurred after 3,300 years ago, significantly later than the major expansion period of Q1a1a-M120 (4.2 kya-3 kya, Figure 1). Nevertheless, several downstream clades of Q1a1a-M120, like F4759 and F4689, exhibit simultaneous expansion with N1a2a1a1a1a1-F710 (Sun et al., 2019). Ancient DNA data suggest that these two paternal lineages were concentrated in ancient populations in northwest China, and co-occurred in some tombs (Zhao et al., 2014; Ma et al., 2021; Wei et al., 2022). These ancient DNA studies also suggest that N1a2a-F1101 is likely the paternal lineage of the royal family of the Zhou Dynasty, while Q1a1a-M120 is the main paternal lineage of the Rong-Di populations (Means “Barbarians” in ancient Chinese). Both paternal lineages became the main paternal component of the Chinese group in later generations. In conclusion, we speculate that Q1a1a-M120 and N1a2a-F1101 together constitute the main paternal lineages of the populations that worked as farmers and pastoralists in northwest China during the Copper-Bronze Age. They played a key role in the emergence of bronze culture, early states, and early civilizations in central region of ancient China.
Bronze age globalization in East Asia
As, discussed in the Introduction section, Bronze age globalization has led to mass replacement and mixing of populations in multiple parts of Eurasia (Allentoft et al., 2015). In East Asia, however, the situation is quite different. Ancient DNA shows that during the Copper-Bronze Age, the populations in the central East Asian region did not experience large-scale replacement, and the genetic components from Indo-Europeans are nearly absent. Based on previous literature and the results of this paper, we suggest that the Gobi Desert on the border between China and Mongolia may have hindered the spread of the Bronze culture and Indo-European-related populations. The hunter-gatherer communities that originally operated in the north and south of the Gobi Desert relied on their familiarity with the environment and long-distance material exchange networks to spread relevant cultural elements as intermediaries. In later historical periods, they became the main founders of the bronze culture populations in northwest China. These demographic histories led to the spread of Bronze culture into central East Asia as a form of cultural diffusion, unlike what happened in other parts of Eurasia during the Bronze Age period of globalization.
In summary, we constructed a high-resolution phylogeny for Y-chromosome haplogroup N1a2a-F1101, one of main paternal lineages of modern Chinese. We explored the demographic of this paternal haplogroup in the past 9,000 years. We also discussed the activity of ancient populations with this lineage and their role during the appearance of Bronze Age culture, the formation of early state and early civilizations in central region of China. The newly-discovered sub-branches and variants will assist in exploring the formation process of gene pool of Chinese populations and their cultural traditions.
Data availability statement
Following the regulations of the Human Genetic Resources Administration of China (HGRAC), the raw sequence data reported in this paper are available on request from the corresponding author. A list of variants of Y-SNP analyzed in this study was included in Supplementary Table S2, which is sufficient for reader to repeat the analyses of this study.
Ethics statement
The studies involving human participants were reviewed and approved by the study and sample collection process were reviewed and approved by the Medical Ethics Committee of Fudan University and Inner Mongolian Normal University, and complied with the ethical principles of the 2013 Helsinki Declaration of the World Medical Association. The patients/participants provided their written informed consent to participate in this study.
Author contributions
Conceptualization, H-XY. Cheliger, X-PW, X-PZ, JS, HL, K-JL, and L-HW; literature search, H-XY. Cheliger, X-PW, X-PZ, and JS; writing—original draft preparation, H-XY. Cheliger, X-PW, K-JL, and L-HW; writing—review and editing, H-XY. Cheliger, X-PW, X-PZ, JS, HL, K-JL, and L-HW. All authors have read and agreed to the published version of the manuscript.
Funding
This study was also supported by Key projects of strategic international scientific and technological innovation cooperation of the Chinese Ministry of Science and Technology (2020YFE0201600 to HL) and Fundamental Research Funds for the Universities of Inner Mongolia Autonomous Region (2022JBBJ004 to L-HW).
Acknowledgments
We thank all donors for providing DNA samples.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The handling editor IY declared a past co-authorship with the author L-HW.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2023.1139722/full#supplementary-material
References
Allentoft, M. E., Sikora, M., Sjogren, K. G., Rasmussen, S., Rasmussen, M., Stenderup, J., et al. (2015). Population genomics of bronze age Eurasia. Nature 522 (7555), 167–172. doi:10.1038/nature14507
Bouckaert, R., Heled, J., Kuhnert, D., Vaughan, T., Wu, C. H., Xie, D., et al. (2014). Beast 2: A software platform for bayesian evolutionary analysis. PLoS Comput. Biol. 10(4), e1003537. doi:10.1371/journal.pcbi.1003537
Cai, D., Han, L., Xie, C.-Z., Li, S.-N., Zhou, H., and Zhu, H. (2007). Mitochondrial DNA analysis of bronze age horses recovered from chifeng region, inner Mongolia, China. Prog. Nat. Sci. 17 (5), 544–550. doi:10.1080/10020070708541034
Cai, X., Qin, Z., Wen, B., Xu, S., Wang, Y., Lu, Y., et al. (2011). Human migration through bottlenecks from Southeast Asia into East Asia during last glacial maximum revealed by Y chromosomes. PLoS One 6 (8), e24282. doi:10.1371/journal.pone.0024282
Damgaard, P. B., Marchi, N., Rasmussen, S., Peyrot, M., Renaud, G., Korneliussen, T., et al. (2018). 137 ancient human genomes from across the Eurasian steppes. Nature 557 (7705), 369–374. doi:10.1038/s41586-018-0094-2
de Barros Damgaard, P., Martiniano, R., Kamm, J., Moreno-Mayar, J. V., Kroonen, G., Peyrot, M., et al. (2018). The first horse herders and the impact of early Bronze Age steppe expansions into Asia. Science 360 (6396), eaar7711. doi:10.1126/science.aar7711
Haak, W., Lazaridis, I., Patterson, N., Rohland, N., Mallick, S., Llamas, B., et al. (2015). Massive migration from the steppe was a source for Indo-European languages in Europe. Nature 522 (7555), 207–211. doi:10.1038/nature14317
Houle, J.-L., and Board, O. H. E. (2014). “Bronze age Mongolia,” in The oxford handbook of topics in archaeology (Oxford, United Kingdom: Oxford University Press).
Hu, K., Yan, S., Liu, K., Ning, C., Wei, L.-H., Li, S.-L., et al. (2015). The dichotomy structure of Y chromosome Haplogroup N. arXiv, 1504.06463. doi:10.48550/arXiv.1504.06463
Huang, Y. Z., Pamjav, H., Flegontov, P., Stenzl, V., Wen, S. Q., Tong, X. Z., et al. (2018). Dispersals of the siberian Y-chromosome haplogroup Q in Eurasia. Mol. Genet. Genomics 293 (1), 107–117. doi:10.1007/s00438-017-1363-8
Новгородова, Э. А. (1989). Древняя ?онголия (?екоторые проблемы хронологии в этнокультурной истории). Москва: Наука.
Ilumäe, A. M., Reidla, M., Chukhryaeva, M., Jarve, M., Post, H., Karmin, M., et al. (2016). Human Y chromosome haplogroup N: A non-trivial time-resolved phylogeography that cuts across language families. Am. J. Hum. Genet. 99 (1), 163–173. doi:10.1016/j.ajhg.2016.05.025
Jeong, C., Balanovsky, O., Lukianova, E., Kahbatkyzy, N., Flegontov, P., Zaporozhchenko, V., et al. (2019). The genetic history of admixture across inner Eurasia. Nat. Ecol. Evol. 3 (6), 966–976. doi:10.1038/s41559-019-0878-2
Jeong, C., Wang, K., Wilkin, S., Taylor, W. T. T., Miller, B. K., Bemmann, J. H., et al. (2020). A dynamic 6,000-year genetic history of Eurasia's Eastern Steppe. Cell 183 (4), 890–904. doi:10.1016/j.cell.2020.10.015
Jeong, C., Wilkin, S., Amgalantugs, T., Bouwman, A. S., Taylor, W. T. T., Hagan, R. W., et al. (2018). Bronze Age population dynamics and the rise of dairy pastoralism on the eastern Eurasian steppe. Proc. Natl. Acad. Sci. U. S. A. 115 (48), E11248–E11255. doi:10.1073/pnas.1813608115
Karmin, M., Saag, L., Vicente, M., Wilson Sayres, M. A., Jarve, M., Talas, U. G., et al. (2015). A recent bottleneck of Y chromosome diversity coincides with a global change in culture. Genome Res. 25 (4), 459–466. doi:10.1101/gr.186684.114
Kilinc, G. M., Kashuba, N., Koptekin, D., Bergfeldt, N., Donertas, H. M., Rodriguez-Varela, R., et al. (2021). Human population dynamics and Yersinia pestis in ancient northeast Asia. Sci. Adv. 7 (2), eabc4587. doi:10.1126/sciadv.abc4587
Kumar, V., Wang, W., Zhang, J., Wang, Y., Ruan, Q., Yu, J., et al. (2022). Bronze and Iron Age population movements underlie Xinjiang population history. Science 376 (6588), 62–69. doi:10.1126/science.abk1534
Lazaridis, I., Nadel, D., Rollefson, G., Merrett, D. C., Rohland, N., Mallick, S., et al. (2016). Genomic insights into the origin of farming in the ancient Near East. Nature 536 (7617), 419–424. doi:10.1038/nature19310
Li, H., and Durbin, R. (2010). Fast and accurate long-read alignment with Burrows-Wheeler transform. Bioinformatics 26 (5), 589–595. doi:10.1093/bioinformatics/btp698
Li, H., and Durbin, R. (2009). Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 25 (14), 1754–1760. doi:10.1093/bioinformatics/btp324
Li, H., Zhao, X., Zhao, Y., Li, C., Si, D., Zhou, H., et al. (2011). Genetic characteristics and migration history of a bronze culture population in the West Liao-River valley revealed by ancient DNA. J. Hum. Genet. 56. 815, 822. doi:10.1038/jhg.2011.102
Liu, J., Jiang, L., Zhao, M., Du, W., Wen, Y., Li, S., et al. (2022). Development and validation of a custom panel including 256 Y-SNPs for Chinese Y-chromosomal haplogroups dissection. Forensic Sci. Int. Genet. 61, 102786. doi:10.1016/j.fsigen.2022.102786
Liu, L., and Chen, X.-C. (2017). The archaeology of China: From the late paleolithic to the early bronze age. Beijing: SDX Joint Publishing Company.
Liu, L. (2004). The Chinese neolithic: Trajectories to early states. Cambridge: Cambridge University Press.
Lu, Q., Cheng, H. Z., Li, L., Yao, H. B., Ru, K., Wen, S. Q., et al. (2020). Paternal heritage of the han Chinese in henan province (central China): High diversity and evidence of in situ neolithic expansions. Ann. Hum. Biol. 47 (3), 294–299. doi:10.1080/03014460.2020.1748226
Ma, P., Yang, X., Yan, S., Li, C., Gao, S., Han, B., et al. (2021). Ancient Y-DNA with reconstructed phylogeny provides insights into the demographic history of paternal haplogroup N1a2-F1360. J. Genet. Genomics 48 (12), 1130–1133. doi:10.1016/j.jgg.2021.07.018
Mao, X., Zhang, H., Qiao, S., Liu, Y., Chang, F., Xie, P., et al. (2021). The deep population history of northern East Asia from the late pleistocene to the Holocene. Cell 184 (12), 3256–3266.e13. doi:10.1016/j.cell.2021.04.040
Mathieson, I., Alpaslan-Roodenberg, S., Posth, C., Szecsenyi-Nagy, A., Rohland, N., Mallick, S., et al. (2018). The genomic history of southeastern Europe. Nature 555 (7695), 197–203. doi:10.1038/nature25778
Narasimhan, V. M., Patterson, N., Moorjani, P., Rohland, N., Bernardos, R., Mallick, S., et al. (2019). The formation of human populations in South and Central Asia. Science 365 (6457), eaat7487. doi:10.1126/science.aat7487
Ning, C., Li, T., Wang, K., Zhang, F., Li, T., Wu, X., et al. (2020). Ancient genomes from northern China suggest links between subsistence changes and human migration. Nat. Commun. 11 (1), 2700. doi:10.1038/s41467-020-16557-2
Ning, C., Wang, C. C., Gao, S., Yang, Y., Zhang, X., Wu, X., et al. (2019). Ancient genomes reveal yamnaya-related ancestry and a potential source of Indo-European speakers in iron age tianshan. Curr. Biol. 29(15), 2526–2532. doi:10.1016/j.cub.2019.06.044
Olalde, I., Brace, S., Allentoft, M. E., Armit, I., Kristiansen, K., Booth, T., et al. (2018). The Beaker phenomenon and the genomic transformation of northwest Europe. Nature 555 (7695), 190–196. doi:10.1038/nature25738
Rasmussen, M., Anzick, S. L., Waters, M. R., Skoglund, P., DeGiorgio, M., Stafford, T. W., et al. (2014). The genome of a Late Pleistocene human from a Clovis burial site in western Montana. Nature 506 (7487), 225–229. doi:10.1038/nature13025
Robbeets, M., Bouckaert, R., Conte, M., Savelyev, A., Li, T., An, D. I., et al. (2021). Triangulation supports agricultural spread of the Transeurasian languages. Nature 599 (7886), 616–621. doi:10.1038/s41586-021-04108-8
Rootsi, S., Zhivotovsky, L. A., Baldovic, M., Kayser, M., Kutuev, I. A., Khusainova, R., et al. (2007). A counter-clockwise northern route of the Y-chromosome haplogroup N from Southeast Asia towards Europe. Eur. J. Hum. Genet. 15 (2), 204–211. doi:10.1038/sj.ejhg.5201748
Sun, J., Li, Y. X., Ma, P. C., Yan, S., Cheng, H. Z., Fan, Z. Q., et al. (2021). 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. Am. J. Phys. Anthropol. 174 (4), 686–700. doi:10.1002/ajpa.24240
Sun, N., Ma, P. C., Yan, S., Wen, S. Q., Sun, C., Du, P. X., et al. (2019). Phylogeography of Y-chromosome haplogroup Q1a1a-M120, a paternal lineage connecting populations in Siberia and East Asia. Ann. Hum. Biol. 46 (3), 261–266. doi:10.1080/03014460.2019.1632930
Tao, R., Li, M., Chai, S., Xia, R., Qu, Y., Yuan, C., et al. (2023). Developmental validation of a 381 Y-chromosome SNP panel for haplogroup analysis in the Chinese populations. Forensic Sci. Int. Genet. 62, 102803. doi:10.1016/j.fsigen.2022.102803
The Y Chromosome Consortium (2002). A nomenclature system for the tree of human Y-chromosomal binary haplogroups. Genome Res. 12 (2), 339–348. doi:10.1101/gr.217602
Wang, C. C., Yeh, H. Y., Popov, A. N., Zhang, H. Q., Matsumura, H., Sirak, K., et al. (2021). Genomic insights into the formation of human populations in East Asia. Nature 591 (7850), 413–419. doi:10.1038/s41586-021-03336-2
Wei, L. H., Wang, L. X., Wen, S. Q., Yan, S., Canada, R., Gurianov, V., et al. (2018). Paternal origin of paleo-Indians in Siberia: Insights from Y-chromosome sequences. Eur. J. Hum. Genet. 26 (11), 1687–1696. doi:10.1038/s41431-018-0211-6
Wei, L. H., Wen, S. Q., Zhou, B. Y., Meng, H. L., and Li, H. (2022). "A study report on the ancient DNA of the remains of the xi-nan-cheng cemetery of western Zhou period in zhang-zi county, shanxi province of China," in A comprehensive study of the xi-nan-cheng cemetery of western Zhou period in zhang-zi county, shanxi province of China eds. B. H. Han (China: Shanxi Provincial Archaeological Research Institute. (Shanghai: Shanghai Ancient Books Publishing House). 135–145.
Xue, J., Wang, W., Shao, J., Dai, X., Sun, Z., Gardner, J. D., et al. (2022). Ancient mitogenomes reveal the origins and genetic structure of the neolithic Shimao population in northern China. Front. Genet. 13, 909267. doi:10.3389/fgene.2022.909267
Yan, S., Wang, C. C., Zheng, H. X., Wang, W., Qin, Z. D., Wei, L. H., et al. (2014). Y chromosomes of 40% Chinese descend from three Neolithic super-grandfathers. PLoS One 9 (8), e105691. doi:10.1371/journal.pone.0105691
Yang, M. A., Fan, X., Sun, B., Chen, C., Lang, J., Ko, Y. C., et al. (2020). Ancient DNA indicates human population shifts and admixture in northern and southern China. Science 369 (6501), 282–288. doi:10.1126/science.aba0909
Yu, H., Spyrou, M. A., Karapetian, M., Shnaider, S., Radzeviciute, R., Nagele, K., et al. (2020). Paleolithic to bronze age siberians reveal connections with first Americans and across Eurasia. Cell 181 (6), 1232–1245. doi:10.1016/j.cell.2020.04.037
Zhang, F., Ning, C., Scott, A., Fu, Q., Bjorn, R., Li, W., et al. (2021). The genomic origins of the bronze age tarim basin mummies. Nature 599 (7884), 256–261. doi:10.1038/s41586-021-04052-7
Zhao, Y.-B., Zhang, Y., Li, H. J., Cui, Y. Q., Zhu, H., and Zhou, H. (2014). Ancient DNA evidence reveals that the Y chromosome haplogroup Q1a1 admixed into the Han Chinese 3,000 years ago. Am. J. Hum. Biol. 26 (6), 813–821. doi:10.1002/ajhb.22604
Zhong, H., Shi, H., Qi, X. B., Duan, Z. Y., Tan, P. P., Jin, L., et al. (2011). Extended Y chromosome investigation suggests post-glacial migrations of modern humans into East Asia via the northern route. Mol. Biol. Evol. 28 (1), 717–727. doi:10.1093/molbev/msq247
Keywords: bronze age, Y-chromosome, N1a2a-F1101, Chinese, phylogeny
Citation: Yu H-X, Ao C, Wang X-P, Zhang X-P, Sun J, Li H, Liu K-J and Wei L-H (2023) The impacts of bronze age in the gene pool of Chinese: Insights from phylogeographics of Y-chromosomal haplogroup N1a2a-F1101. Front. Genet. 14:1139722. doi: 10.3389/fgene.2023.1139722
Received: 07 January 2023; Accepted: 23 February 2023;
Published: 10 March 2023.
Edited by:
Ivy Yeh, Nanyang Technological University, SingaporeReviewed by:
Longli Kang, Xizang Minzu University, ChinaAngélica G. Martínez-H, National Institute of Genomic Medicine (INMEGEN), Mexico
Copyright © 2023 Yu, Ao, Wang, Zhang, Sun, Li, Liu and Wei. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Kai-Jun Liu, jspzliukaijun@foxmail.com; Lan-Hai Wei, Ryan.lh.wei@gmail.com
†These authors have contributed equally to this work