Sex-Biased Population Admixture Mediated Subsistence Strategy Transition of Heishuiguo People in Han Dynasty Hexi Corridor

Abstract

The Hexi Corridor was an important arena for culture exchange and human migration between ancient China and Central and Western Asia. During the Han Dynasty (202 BCE–220 CE), subsistence strategy along the corridor shifted from pastoralism to a mixed pastoralist-agriculturalist economy. Yet the drivers of this transition remain poorly understood. In this study, we analyze the Y-chromosome and mtDNA of 31 Han Dynasty individuals from the Heishuiguo site, located in the center of the Hexi Corridor. A high-resolution analysis of 485 Y-SNPs and mitogenomes was performed, with the Heishuiguo population classified into Early Han and Late Han groups. It is revealed that (1) when dissecting genetic lineages, the Yellow River Basin origin haplogroups (i.e., Oα-M117, Oβ-F46, Oγ-IMS-JST002611, and O2-P164+, M134-) reached relatively high frequencies for the paternal gene pools, while haplogroups of north East Asian origin (e.g., D4 and D5) dominated on the maternal side; (2) in interpopulation comparison using PCA and Fst heatmap, the Heishuiguo population shifted from Southern-Northern Han cline to Northern-Northwestern Han/Hui cline with time, indicating genetic admixture between Yellow River immigrants and natives. By comparison, in maternal mtDNA views, the Heishuiguo population was closely clustered with certain Mongolic-speaking and Northwestern Han populations and exhibited genetic continuity through the Han Dynasty, which suggests that Heishuiguo females originated from local or neighboring regions. Therefore, a sex-biased admixture pattern is observed in the Heishuiguo population. Additionally, genetic contour maps also reveal the same male-dominated migration from the East to Hexi Corridor during the Han Dynasty. This is also consistent with historical records, especially excavated bamboo slips. Combining historical records, archeological findings, stable isotope analysis, and paleoenvironmental studies, our uniparental genetic investigation on the Heishuiguo population reveals how male-dominated migration accompanied with lifestyle adjustments brought by these eastern groups may be the main factor affecting the subsistence strategy transition along the Han Dynasty Hexi Corridor.

Introduction

Human history can be seen as a history of dealing with new challenges caused by changes of factors including resource distribution and social relationships. To overcome them, human beings have been using extrasomatic ways, including subsistence-, socio-, and ideo-technologies to create new niches for survival (Binford, 1962; Zhang 2021). Among the studies of adaptive changes, the evaluation of the factors leading to significant changes in subsistence strategy in human prehistory and history is a fascinating topic. Climate change, technological innovation, rapidly growing population, trans-continental cultural exchange, human migration, and geopolitics are clear potential candidates (Cohen, 1975; Bonsall et al., 2002; Gao et al., 2007; Pokharia et al., 2017; Cheung et al., 2019; Petraglia et al., 2020; Yang et al., 2020; Li et al., 2021). Climate change is well-studied and considered one of the most guiding factors of shifting subsistence strategies. A suitable climate will be conducive to the promotion of agricultural development, as rainy weather conditions facilitated the development of oasis agriculture in South-east Arabia around 5,100 cal BP (Preston, 2011), to take one example. On the other hand, an extraordinarily harsh climate can result in civilizational collapse. Humans are required to adopt different strategies, as well as adjust subsistence strategies, even migrate in search of better conditions, in order to cope with and adapt to such environmental changes (Polyak and Asmerom, 2001; Nunez et al., 2002; Haug et al., 2003; An et al., 2005; Preston et al., 2012; King et al., 2013; d’Alpoim Guedes et al., 2015; Jia et al., 2016; Pokharia et al., 2017).

The Hexi Corridor, located centrally on the eastern Silk Road, once played a crucial role in cultural exchange between east and west. The area was also a crossroads of agricultural and nomadic populations within China. Multiple lines of evidence showed that a significant shift in subsistence strategies along the Hexi Corridor occurred during the Han Dynasty (Ma et al., 2016; Yang et al., 2019a; Yang et al., 2020). Historical records from the Shiji (Records of the Grand Historian, 史记) to the Hanshu (Book of Han, 汉书) stated that the Hexi Corridor was occupied by nomadic pastoralists (Xiongnu, Yuezhi, and Wusun) before the Han Dynasty, a claim supported by excavated relics and faunal remains from Shajing culture (2,900–2,100 cal yr BP) and Shanma culture (2,700–2,100 cal yr BP) (Yang et al., 2019b) sites. The importance of domesticated pig declined precipitously during this same period (GPICRA, 2001), while sheep/goat, cattle, horse, and camel emerged as the major domesticated animals along the Hexi Corridor (GPICRA and SAMPU, 2011). A wealth of leather and woolen products had also been found at Shajing culture sites, suggesting how residents had initiated their own secondary products revolution with these domesticated pastoral animals (Sherratt, 1981). These were all indications suggesting that Hexi populations lived highly mobile nomadic lifestyles. During and after the Han Dynasty (202 BCE–220 CE), however, agriculture developed rapidly to become the subsistence strategy in this region. Numerous domesticated crop types had been found at archaeological sites, including barley, wheat, millet, highland barley, and pea (Sun and Liu, 2014), while advanced iron implements including iron plough, iron sickle, and iron spade had been excavated from numerous Han Dynasty sites (Yang, 2015). This strongly suggested that an advanced agricultural technology became widespread in the Hexi Corridor during this period of time. Remains of chicken, pigs, dogs, sheep/goat, cattle, and horse identified at the Heishuiguo Han Dynasty tombs should be placed in this context. Chicken would emerge as the most common domestic animal in this period, followed by pig (Li, 2021). This Heishuiguo population of domesticates resembled what we find with Central Plain farmers while being vastly different from that of nomads (the assemblage of domestic animals was sheep/goat, horse, and cattle) (Deng, 2015), who generally fed the camel and horse for long-distance migration, while pig and chicken may be more likely to appear in settled-peoples’ homes (Yang et al., 2019a). Finally, large numbers of painted murals depicting farming and animal-grazing had been found at Hexi Corridor sites, revealing this predominant Han-Jin Dynasties mixed economy (Zheng and Gao, 2019). We can be sure, therefore, that subsistence strategy in the Hexi Corridor shifted from nomadic style (pre-Han Dynasty) to a mixed style (during and post-Han Dynasty). The argument was also supported by stable isotope data (Li, 2021).

Located in an arid climatic transitional belt, the Hexi Corridor is extremely sensitive to changes in its environment. Previous studies have primarily attempted to explain shifting subsistence strategies through the mirror of environmental change (Zhou et al., 2016; Shi et al., 2018; Yang et al., 2019a; Yang et al., 2020). But we believe that the influence of population migration should not be overlooked, especially when considering the Hexi Corridor’s unique geographical position. Uniparentally inherited markers (Y chromosome or mitochondria) have been widely used in human population migration studies. The non-recombining portion of the Y chromosome (NRY) is strictly inherited paternally, while mtDNA is inherited maternally (Calafell and Larmuseau, 2017). MtDNA and Y chromosome therefore provide a matrilineal and patrilineal demographic history, respectively, revealing pictures of sex-specific processes in the past. Notably, genetic history revealed by Y chromosomes need not be identical to that by mtDNA. By using both mtDNA and Y chromosome inherited markers, sex-biased migration has been frequently found in studies of human populations and therefore might truly reflect the influence of social behaviors (Oota et al., 2001). So as to discuss the relationship between the change of subsistence strategies in the Hexi Corridor and population migrations during the Han Dynasty, here we analyze the Y-chromosome (including 485 Y-SNP markers) and mitochondrial genomes from 31 samples from the Heishuiguo site covering the whole Han Dynasty.

Materials and Methods

Materials

The Heishuiguo site is located in Ganzhou county, Zhangye city, Gansu province, China (Figure 1). The site was excavated by the Institute of Cultural Relics and Archaeology of Gansu Province in 2018. Tomb distribution at Heishuiguo reveals a large cemetery consisting of family burial grounds and scattered burial groups (Chen et al., 2019). Based on ¹⁴C dating, tomb morphology and grave good assemblages (Li, 2021), burials at Heishuiguo have been divided into four phases: Phase 1, during the middle Western Han Dynasty (118–49 BCE); Phase 2, during the late Western Han Dynasty (48 BCE–6CE); Phase 3, from the Wangmang Xin Dynasty through the early Eastern Han Dynasty (7–67 CE); and Phase 4, from the middle to the late Eastern Han Dynasty (67–191 CE) (Zhang, 2017; Chen et al., 2019).

FIGURE 1

This study classified 31 individuals into two groups, the Early Han (Phases 1–2) and Late Han Dynasty (Phases 3–4), respectively. Heishuiguo individual sex was determined by pelvic (Klales et al., 2012) and skull morphology (Buikstra, 1994). Details are provided in Table 1 and Supplementary Table S1.

Methods

DNA Extraction

We extracted DNA from 31 samples using a dedicated aDNA facility at Fudan University, Shanghai, following the established precautions for working with ancient human DNA (Knapp et al., 2012; Sun et al., 2021). Negative extraction control samples (no sample powder used) were included, to monitor against contamination, as well as library negative controls (extract supplemented by water) in every batch of samples, which were processed and carried through the entire wet laboratory processing.

Prior to sampling, all samples were irradiated with UV-light for 30 min on each side and wiped with a 5% bleach solution. Teeth and other osseous materials were sandblasted to remove the outer surface before being ground to fine power with a mixer mill (Retsch, Germany). A dense section of temporal bones was cut around the cochlea by first removing the outer part, then grinding the clean inner part into fine powder. A total of 100 mg of bone powder was then used to extract DNA. Pre-lysis methods included the addition of a 1 ml extraction buffer, containing 0.5 M EDTA, 0.25 mg/ml Proteinase K (Merck, Germany), pH 8.0, after which samples were rotated for 1 hour at a temperature of 37°C. After centrifugation, the supernatant was discarded and 2.5 ml extraction buffer added, followed by overnight rotation at 37°C. We mixed 20 μl magnetic beads (Enlighten Biotech, China) with 12.5 ml binding buffer containing 5 M GuHCl, 40% Isopropanol, 25 mM sodium acetate, and 0.05% Tween-20 (Merck, Germany), (PH 5.2). We then transferred the supernatant (∼2.5 ml) to a binding buffer/bead mixture prior to the robotic extraction (Enlighten Biotech, China) procedure. At last, the DNA underwent elution through a 50 μl TET buffer (QIAGEN, Germany).

Library Preparation

We prepared double-stranded libraries in accordance with Meyer’s protocols (Meyer and Kircher, 2010), but with minor alterations. Libraries were amplified with indexing primers in two parallel polymerase chain reactions (PCR) using Q5 High-Fidelity DNA Polymerase (New England Biolabs, USA). Indexed products from the same library were pooled and purified using Agencourt AMPure XP beads (Beckman Coulter, Germany) and then eluted in 20 μl TET buffer. We quantitated the clean-up libraries using a Qubit 2.0 Fluorometer (Thermo Fisher, USA). We then sequenced the libraries on an Illumina HiSeq X10 instrument at the Annoroad Company (Beijing, China) using the 150-bp paired-end sequencing design.

Mitochondrial Capture and Sequencing

Target enrichment of the mitogenome was performed on each amplified library using a MyGenostics Human Mitochondria Capture Kit (MyGenostics Company, Beijing, China) as described by Sun et al. (2021), with hybridization and wash temperatures lowered to 60°C and 55°C to facilitate the enrichment of our short library molecules, in line with Dabney et al. (2013). A final post-enrichment amplification was performed for 15 cycles. The post-enrichment amplified product was then quantified via qPCR. Sequencing was performed using a Novaseq 6,000 platform at Mingma Technologies Company (Shanghai, China). Next, 150 bp paired-end reads were generated according to the manufacturer’s instructions.

Multiplex Polymerase Chain Reactions Targeted Amplification and Sequencing for Y Chromosome

Ancient DNA fragment lengths are generally skewed toward short fragments, the vast majority of which are typically shorter than 100 bp (Sawyer et al., 2015). Furthermore, in order to preliminary screen ancient samples quickly and inexpensively, we opted for multiplex PCR targeting enrichment with short amplicons based on the NGS (Next Generation Sequencing) platform. After amplification enrichment, a large number of samples could be detected and analyzed in parallel using the NGS platform. In view of the characteristics of highly degraded DNA involved in this study, we designed a more sensitive short amplifier primer system (Wen et al., 2019) and conducted tests on samples from the Heishuiguo site. The system comprises 485 Y-SNPs (Figure 2; Supplementary Tables S2, S3), covering common lineages in East Asia. Details can be found in the Supplementary Material. PCR amplification, sequencing, and data analysis can also be found in the Supplementary Materials.

FIGURE 2

Sequence Data Processing

For shotgun and mtDNA capture data, we clipped sequencing adapters and merged using sequences by ClipAndMerge v1.7.8 (Peltzer et al., 2016). Following this, we mapped merged reads to the human reference genome (hs37d5; GRCh37 with decoy sequences) using BWA v0.7.17 (Li and Durbin, 2010). PCR duplicates were removed using Dedup v0.12.3 (Peltzer et al., 2016). Using trimBam implemented in BamUtil v1.0.14 (https://github.com/statgen/bamUtil), we clipped four bases from both ends of each read to avoid an excess of remaining C->T and G->A transitions at the ends of the sequences.

Ancient DNA Authentication

The authenticity of the ancient genome sequence was mainly determined by the combination of two observations of the same specimen. Firstly, deep sequencing of the mitochondrial genome would show that the vast majority of the DNA fragments had come from a single individual. Second, the patterns of DNA degradation (Supplementary Figure S1), in particular nucleotide misincorporations resulting from deamination of cytosine residues at the ends of DNA fragments, would indicate that the mtDNA is ancient. We first checked DNA damage pattern and estimated the 5′ C>T and 3′ G>A misincorporation rate using mapDamage v 2.0.61 (Jónsson et al., 2013). We then made use of a Schmutzi program to test mitochondrial contamination rates for all individuals (Renaud et al., 2015).

Uniparental Haplogroup Assignment

For mtDNA, we employed a log2fasta program implemented in Schmutzi (Renaud et al., 2015) in order to call the mtDNA consensus sequences. Mutations that appeared when checked against rCRS were also re-checked in BAM (Binary Alignment Map) files through visual inspection with IGV software (Helga et al., 2013). Then, we used Haplogrep 2 (Weissensteiner et al., 2016) to assign the haplogroups.

Y chromosome haplogroups were examined by aligning a set of positions in the ISOGG (International Society of Genetic Genealogy, http://isogg.org/) and Y-full (https://www.yfull.com/tree/) databases, and analysis performed in the case of a base and mapping quality exceeding 30. Haplogroup determination was performed with the script Yleaf.py in Yleaf software (Ralf et al., 2018), which provides outputs for allele counts of ancestral and derived SNPs along a path of branches of the Y-chromosome tree. Finally, we re-checked the SNPs by visual inspection with IGV software (Helga et al., 2013).

Statistical Analyses

Principal component analysis (PCA) was performed using FactoMineR package of R 3.6.3. Reference populations are listed in Supplementary Tables S4–S6. The pairwise genetic Fst was calculated using Arlequin 3.5 software. Heatmaps were also used to illustrate the clusters based on Fst by stats package of R 3.6.3. At last, to visualize the origins, we used PC1 values from PCA plot to generate contour maps with Surfer 12.0 software (Golden Software, https://www.goldensoftware.com/).

Results

Y-Chromosome and mtDNA Haplogroup Profile

Y chromosome haplogroups of 30 Heishuiguo males were determined according to the ISOGG’s Y-DNA Haplogroup Tree 2019 (Figure 3; Supplementary Table S1). Overall, the Heishuiguo population was comprised of the haplogroups Oγ-IMS-JST002611(26.7%), Oα-M117 (13.3%), C2-M217(13.3%), N-F1206 (13.3%), O1b1a2-Page59 (10%), O1a-M119+, P203- (6.7%), O2-P164+, M134- (6.7%), Oβ-F46 (3.3%), O1b1a1a-M95 (3.3%), and Q-M120 (3.3%). Haplogroups Oγ-IMS-JST002611 (26.7%), Oα-M117 (13.3%), and Oβ-F46(3.3%) represent three major founder paternal lineages (Yan et al., 2014; Wen et al., 2016) and encompass more than 40% of the present-day Han Chinese (estimated 16% for Oα, 11% for Oβ, and 14% for Oγ) (Yan et al., 2011). These three lineages have been considered to be derived from Neolithic farmers in Yellow River Basin (Yan et al., 2014). The three haplogroups are also the most frequent among the Heishuiguo population, representing 43.3% of the total, similar to modern Han populations. Haplogroup O2-P164+, M134-, may have expanded similarly with the above three haplogroups. Haplogroups C2-M217, N-F1206, and Q-M120 are common in Chinese, Altaic, Uralic, and Northern Eurasian Indo-European populations. The vast majority of haplogroup C in China belongs to C2-M217 (Zhong et al., 2010), constituting ∼10% of Han Chinese, as well as a great portion of Altaic-speaking populations, i.e., Mongol, Manchu, and Kazakh people. In the Heishuiguo population, the haplogroup C2-M217 can be further classified into C2*-M217 and C-F1144. The latter is a southern clade (Yan et al., 2014) and one of the six major paternal lineages of the Han Chinese (Wen et al., 2016; Wu et al., 2020). N-F1206 is called the northern clade of N haplogroup and widely distributed across Northern Eurasia (Hu et al., 2015). N-F1206 can be sub-divided into N-TAT and N-F710. The highest frequency of N-TAT is found in Altaic, Uralic, and Indo-European speaking populations (Ilumäe et al., 2016), such as Vilyuy Yakuts (91.525%), Evenks (50.877%), Buryats (41.441%), Udmurts (66.667%), Finns (53.846%), and Latvians (43.023%), while the N-F710 is considered as having migrated from Northeast Asia southward into the Yellow River region at no later than 2.7 kya (Ma et al., 2021). Haplogroup Q-M120 originated in South Siberia and expanded across northwestern China between 5–3 kya (Sun et al., 2019). The lineage was absorbed into ancient Huaxia (Han Chinese) populations before 2 kya and would eventually become one of the six founder lineages in modern Han populations (Wen et al., 2016; Sun et al., 2019). Finally, the O1b1a2-Page59, O1a-M119+, P203-, and O1b1a1a-M95 haplogroups, making up 20% of the Heishuiguo haplogroup proportions, are of minor southern East Asia origin (Karafet et al., 2010; Cai et al., 2011; Trejaut et al., 2014; Zhang et al., 2015; Luo et al., 2020). The paternal diversity of the Late Han Heishuiguo group exceeds that of the Early Han group, and features the addition of haplogroups N-F710, C-F1144, Oβ-F46, and Q-M120.

FIGURE 3

MtDNA haplogroups of 27 Heishuiguo samples were determined based on PhyloTree mtDNA tree Build 17 (Figure 4; Supplementary Table S1). Different from the Y chromosome, the mtDNA gene pool is more heterogeneous. In the Heishuiguo population, these haplogroups consisted of D4 (25.93%), D5(18.52%), B5 (11.11%), R11(11.11%), B4 (3.7%), C4 (3.7%), F1(3.7%), G1(3.7%), G3 (3.7%), M11 (3.7%), M33 (3.7%), M9 (3.7%), and N9 (3.7%). Among these, haplogroups D4, D5, C4, G1, G3, M11, and M9 have origins in north East Asia. Haplogroup D4 occurs with the greatest frequency (25.93%) in the Heishuiguo population, and also retains a very high frequency in northern Asian (average 16.7%), central Asian (average 15.3%), and eastern Asian populations (average 22.5%) (Derenko et al., 2010). In this study, D4 could be further divided into sub-clades D4, D4a3b*, D4a6, D4b2b*, and D4j*. As for D5, this haplogroup spread across northern East Asia with moderate to low frequency and reaches its highest levels in Tubalar (25%) (Volodko et al., 2008), Beijing Han (15.38%) (Jin et al., 2009), Shannan Tibetan (15%) (Xu and Hu, 2015), Henan Han (11.4%) (Xu and Hu, 2015), Orochens (11.4%) (Kong et al., 2003), Koreans (10.4%) (Kong et al., 2003) from China, and Shandong Han (10%) (Yao et al., 2002). Here, D5 includes the sub-clades D5a2*, D5a2a, D5a2a1+@16172*, D5b1b*, and D5c*. These two dominant haplogroups, D4 and D5, with 44.44% of the Heishuiguo total, have a Northern Asian origin and distribution (Derenko et al., 2010). The remaining C4, G1, G3, M11, and M9 haplogroups are more common in northern East Asia than southern Eastern Asia. The B5, B4, and F1 haplogroups (occurring with 18.52% frequency at Heishuiguo) are relatively common in southern East Asia and may indicate gene influx from these regions. The haplogroup B5 retains high frequency in southern East Asia, especially in Tai-Kadai populations such as the Seak (69.23%) (Kutanan et al., 2017) and Kalueng (40%) (Kutanan et al., 2017) in Northeastern Thailand, Austro-Asiatic population like the Bru (45.83%) (Kutanan et al., 2017) from Northeastern Thailand, and Hmong-Mien population such as Hunan Yao (29.17%) (Wen et al., 2005) and Guizhou Han (10.19%) (Li et al., 2019). B5 is comprised of the haplogroups B5a*, B5a2a1a, and B5b2a2*. Likewise, the haplogroups F1 and B4 were common in southern East Asia and have southern Asian origins (Wen et al., 2005; Ko et al., 2014; Kutanan et al., 2017). Finally, the haplogroups R11, N9, and M33 occur sporadically in East Asian populations. We noted a similar degree of genetic diversity among both the Early Han and Late Han groups. In summary, the dominant mtDNA haplogroups at Heishuiguo, such as D4 and D5, are more frequent in northern East Asian populations and exhibit a northern Asian origin, while haplogroups B5, B4, and F1 might reflect gene flows from southern China.

FIGURE 4

Population Comparisons

In order to investigate the genetic relationships between Heishuiguo population and reference populations, we conducted a principal component analysis (PCA) based on haplogroup frequencies (Figure 5, Supplementary Figure S2; Supplementary Tables S4–S6). A genetic distance (Fst) heatmap was also visualized, to further explore population relationships (Figures 6, 7; Supplementary Figure S3). In this study, the Heishuiguo population was classified into three groups: Early Han, Late Han, and Overall group merging the former two groups. The PCA map (Figure 5) shows a division between the Early Han and Late Han groups according to PC1. Our Y-chromosome principal component plot (Figure 5A) reveals an overall Heishuiguo population is projected on to the Southern Han and Northern Han cline. More specifically, the Early Group appears to cluster around Southern Han populations, especially Fujian Han and Guangdong Han. According to previous studies, these two latter populations were descended from Northern China immigrations beginning in the Han Dynasty (202 BCE–220 CE) (Ge et al., 1997; Wen et al., 2004a). Our Early Heishuiguo population may therefore also reflect migration, in this case from the Central Plains to Hexi Corridor. The Late Group clusters closely with the Northern Han population in PCA plot and additionally clusters with Northern Han and Northwestern Hui in our heatmap (Figure 6), suggesting the genetic admixture with indigenous people over time. On the maternal side, the three Heishuiguo populations are tightly congregated (Figure 5B), mirroring genetic continuity from the Early Han to Late Han. Meanwhile, the Heishuiguo population maps out close to Mongolian-speaking (e.g. Buryat) and Northwestern Han (e.g., Ningxia Han) populations, which indicates an indigenous maternal origin.

FIGURE 5

FIGURE 6

FIGURE 7

In summary, an interpopulation comparison reveals that the Heishuiguo population shows close affinity with Han Chinese in terms of paternal structure and Mongolic and Northwestern Han populations in terms of maternal structure. This is a clear indication of sex-biased population admixture.

The Origins of the Heishuiguo Population

The above-mentioned PCA plot shows the eastern and western population cline according to PC1. We can make further use of PC1 values to generate genetic contour maps (Figure 8) and further visualize the possible origins of Heishuiguo population. On the paternal side, the contour map with early Heishuiguo group (Figure 8A) shows PC1 values gradually increasing from east-to-west across northwest China before dropping abruptly in the Hexi Corridor, a clue to the non-local and likely eastern origins of the population. Such a pattern does not repeat for the Late Han Heishuiguo group (Figure 8B). On the other hand, the maternal side paints a different picture (Figure 8C and Figure 8D), where we do not observe significant fluctuations between the Early Han and Late Han groups. This result points to the local roots of Heishuiguo females and their genetic continuity covering the whole Han Dynasty, coinciding with our PCA plot.

FIGURE 8

Discussion

Multiple lines of evidence (i.e., historical records, archeological finds, and stable isotope analysis) have demonstrated how subsistence strategies along the Hexi Corridor shifted from a nomadic economy to mixed economy (i.e., pastoralism and farming) from the Han Dynasty (Sun and Liu, 2014; Li 2021). Prior research has focused on exploring the relationship between climate change and subsistence strategy transition in this region during the Han, suggesting a cold and arid climate had dried up considerable stretches of river and impoverished fertile lands prior to the Han Dynasty, with resulting prevalent nomadism. A warmer and wetter climate from the Han onwards promoted the thriving of agriculture (Yang et al., 2020). Up till now, it has seemed farfetched to link this transition to climate factor simply.

On basis of Y chromosome and mtDNA profiles of Heishuiguo population, we have located a Heishuiguo patrilineal population consisting of Yellow River Basin origin haplogroups (i.e., Oα-M117, Oβ-F46, Oγ-IMS-JST002611, and O2-P164+, M134-) at an overall rate of over 50%, along with southern East Asian origin haplogroups (O1a and O1b) at ∼20% and northern East Asian origin haplogroups (C2-M217, N-F1206, and Q-M120) at ∼30%. The matrilineal Heishuiguo population consisted of northern East Asian haplogroups (e.g., haplogroups D4, D5, and C4), which accounted for ∼62.95%, alongside southern East Asian haplogroups (such as B5, B4, and F1) at 18.52%. By way of interpopulation comparisons (PCA and Fst heatmap) with reference populations (e.g., Southern Han, Northern Han, Hui, Mongolic, and Tibetan), we have been able to show closer paternal genetic affinity with Northern Han and Hui populations among the groups at Heishuiguo. Via PCA, we observed genetic structure changes from Southern-Northern Han cline to Northern-Northwestern Han/Hui cline with time (Figure 5A), indicating genetic admixture between Yellow River immigrants and natives. Historical records and archaeological finds add further credence to our results. According to historical documents (Qi, 1983; Si, 2002; Ban, 2008; Bamboo Slips Museum, 2013), the Han Dynasty government migrated population on a large scale from about 20 counties in Yellow River Basin to four counties (i.e. Zhangye 张掖郡, Jiuquan 酒泉郡, Wuwei 武威郡, and Dunhuang 敦煌郡) in Hexi to strengthen administration and control over this region. Bamboo Slips (简牍) provide more detail and even give us the possible source of these male immigrants in 21 counties from the Middle and Lower Yellow River (Hongnong 弘农郡, He’nei 河内郡, Langya琅邪郡, Changyi 昌邑国, Pinggan State平干国, Dahe 大河郡, Chenliu 陈留郡, Runan 汝南郡, Julu 巨鹿郡, Yingchuan 颍川郡, Shangdang 上党郡, Henan 河南郡, Jiyin 济阴郡, Nanyang 南阳郡, Hedong 河东郡, Zhao State赵国, Dong 东郡, Liang State梁国, Zhangye张掖郡, Huaiyang 淮阳郡, and Wei 魏郡) (Liu, 2012; Li, 2018, Figure 9). Our view from the maternal mtDNA (Figure 5B), however, shows the Heishuiguo population closely clustered with certain Mongolic and Northwestern Han populations and exhibiting a genetic continuity covering the whole Han Dynasty, suggesting a possible local origin for Heishuiguo females. This is in accordance with historical records (Liu 2012), where major migration events were often male-dominated migratory and frequently involved migration for garrison building, political migration, and migration of minority groups. Young males were usually the ones building garrisons, and most couldn’t bring their families outside of small local garrisons. Political migration involved political prisoners, ordinary crimes, and victims of natural calamities. Migration of minority groups targeted rebels from border areas such as Di and Qiang peoples in the Upper Yellow River Basin. Among them, the military migration was in the majority. This is likely why a sex-biased admixture pattern can be clearly observed in the Heishuiguo population. Such sex-biased admixture patterns have also been observed in population expansion of Han and Tibeto-Burman-speaking (Wen et al., 2004a; Wen et al., 2004b) people. We plotted genetic contour maps to visualize the possible origins of Heishuiguo population. As shown in Figure 8, we can easily observe that the primary eastern East Asian origin of male ancestry and native origins of female ancestry.

FIGURE 9

In this study, by means of a uniparental genetic analysis, we have observed a male-dominated admixture event occurring during this period, one additionally supported by historical records and archaeological findings. That shifting subsistence strategy along the Hexi Corridor residents that kept pace with human migration cannot be coincidental. Mass migration of individuals and transplantation of subsistence lifestyles would have impacted the former subsistence strategy in Hexi Corridor during the Han Dynasty. This study provides new insights and possibilities into how population admixture serves as a key factor in changes of subsistence strategy.

References

• 1

  AnC.-B.TangL.BartonL.ChenF.-H. (2005). Climate Change and Cultural Response Around 4000 Cal Yr B.P. In the Western Part of Chinese Loess Plateau. Quat. Res.63 (3), 347–352. 10.1016/j.yqres.2005.02.004

  • CrossRef
  • Google Scholar

• 2

  BanG. (2008). Hanshu. The Biography of Zhaochongguo. Taiyuan: Sanjin publishing House, 143.

  • Google Scholar

• 3

  BinfordL. R. (1962). Archaeology as Anthropology. Am. Antiq.28 (2), 217–225. 10.2307/278380

  • CrossRef
  • Google Scholar

• 4

  BonsallC.MacklinM. G.AndersonD. E.PaytonR. W. (2002). Climate Change and the Adoption of Agriculture in north-west Europe. Eur. J. Archaeol.5 (1), 9–23. 10.1179/eja.2002.5.1.9

  • CrossRef
  • Google Scholar

• 5

  BuikstraJ. E. (1994). Standards for Data Collection from Human Skeletal Remains. Fayetteville, Arkansas: Arkansas archeological survey. Research series no. 44.

  • Google Scholar

• 6

  CaiX.QinZ.WenB.XuS.WangY.LuY.et al (2011). Human Migration through Bottlenecks from Southeast Asia into East Asia during Last Glacial Maximum Revealed by Y Chromosomes. PloS one6 (8), e24282. 10.1371/JOURNAL.PONE.0024282

  • CrossRef
  • Google Scholar

• 7

  CalafellF.LarmuseauM. H. D. (2017). The Y Chromosome as the Most Popular Marker in Genetic Genealogy Benefits Interdisciplinary Research. Hum. Genet.136, 559–573. 10.1007/s00439-016-1740-0

  • CrossRef
  • Google Scholar

• 8

  ChenG. K.YangY.LiuB. (2019). Ganzhou, Zhangye: The Excavation Report of Han Dynasty Cemetery in Heishuiguo Site (Volume II). Lanzhou: Gansu Education Publishing House.

  • Google Scholar

• 9

  CheungC.ZhangH.HepburnJ. C.YangD. Y.RichardsM. P. (2019). Stable Isotope and Dental Caries Data Reveal Abrupt Changes in Subsistence Economy in Ancient China in Response to Global Climate Change. PLOS ONE14 (7), e0218943. 10.1371/journal.pone.0218943

  • CrossRef
  • Google Scholar

• 10

  CohenM. N. (1975). Archaeological Evidence for Population Pressure in Pre-agricultural Societies. Am. Antiq.40 (4), 471–475. 10.2307/279335

  • CrossRef
  • Google Scholar

• 11

  DabneyJ.KnappM.GlockeI.GansaugeM.-T.WeihmannA.NickelB.et al (2013). Complete Mitochondrial Genome Sequence of a Middle Pleistocene Cave bear Reconstructed from Ultrashort DNA Fragments. Proc. Natl. Acad. Sci.110 (39), 15758–15763. 10.1073/PNAS.1314445110

  • CrossRef
  • Google Scholar

• 12

  d’Alpoim GuedesJ.JinG.BocinskyR. K. (2015). The Impact of Climate on the Spread of Rice to North-Eastern China: A New Look at the Data from Shandong Province. PloS one10 (6), e0130430. 10.1371/journal.pone.0130430

  • CrossRef
  • Google Scholar

• 13

  DengH. (2015). Animal Remains in Han Danasty Tomb. Cult. Relics Southern China3, 58–69.

  • Google Scholar

• 14

  DerenkoM.MalyarchukB.GrzybowskiT.DenisovaG.RogallaU.PerkovaM.et al (2010). Origin and post-glacial Dispersal of Mitochondrial DNA Haplogroups C and D in Northern Asia. PloS one5 (12), e15214. 10.1371/journal.pone.0015214

  • CrossRef
  • Google Scholar

• 15

  Gansu Bamboo Slips Museum (GBSM) (2013). Han Dynasty Wooden Slips from Jianshuijinguan (Third). Shanghai: Zhongxi Book Company, 52.

  • Google Scholar

• 16

  GaoH.ZhuC.XuW. (2007). Environmental Change and Cultural Response Around 4200 Cal. Yr BP in the Yishu River Basin, Shandong. J. Geogr. Sci.17, 285–292. 10.1007/s11442-007-0285-5

  • CrossRef
  • Google Scholar

• 17

  GeJ. X.WuS. D.ChaoS. J. (1997). Zhongguo Yimin Shi (The Migration History of China). Fuzhou: Fuzhou: Fujian People’s Publishing House.

  • Google Scholar

• 18

  GPICRA and SAMPU (Gansu Provincial Institute of Cultural Relics and Archaeology and School of Archaeology and Museology of Peking University) (2011). The Report of Prehistoric Archaeology Survey in the Hexi Corridor. Beijing: Cultural Relics Publishing House.

  • Google Scholar

• 19

  GPICRA (Gansu Provincial Institute of Cultural Relics and Archaeology) (2001). Excavation Report of Xigangchaiwangang Tomb in Yongchang. Lanzhou: Gansu People’s Publishing House.

  • Google Scholar

• 20

  HaugG. H.GüntherD.PetersonL. C.SigmanD. M.HughenK. A.AeschlimannB. (2003). Climate and the Collapse of Maya Civilization. Science299, 1731–1735. 10.1126/science.1080444

  • CrossRef
  • Google Scholar

• 21

  HuK.YanS.LiuK.NingC.WeiL. H.LiS. L.et al (2015). The Dichotomy Structure of Y Chromosome Haplogroup N. arXiv. arXiv:1504.06463.

  • Google Scholar

• 22

  IlumäeA.-M.ReidlaM.ChukhryaevaM.JärveM.PostH.KarminM.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. 10.1016/j.ajhg.2016.05.025

  • CrossRef
  • Google Scholar

• 23

  JiaX.SunY.WangL.SunW.ZhaoZ.LeeH. F.et al (2016). The Transition of Human Subsistence Strategies in Relation to Climate Change during the Bronze Age in the West Liao River Basin, Northeast China. The Holocene26 (5), 781–789. 10.1177/0959683615618262

  • CrossRef
  • Google Scholar

• 24

  JinH.-J.Tyler-SmithC.KimW. (2009). The Peopling of Korea Revealed by Analyses of Mitochondrial DNA and Y-Chromosomal Markers. PloS one4 (1), e4210. 10.1371/journal.pone.0004210

  • CrossRef
  • Google Scholar

• 25

  JónssonH.GinolhacA.SchubertM.JohnsonP. L. F.OrlandoL. (2013). mapDamage2.0: Fast Approximate Bayesian Estimates of Ancient DNA Damage Parameters. Bioinformatics (Oxford, England)29 (13), 1682–1684. 10.1093/bioinformatics/btt193

  • CrossRef
  • Google Scholar

• 26

  KarafetT. M.HallmarkB.CoxM. P.SudoyoH.DowneyS.LansingJ. S.et al (2010). Major East-West Division Underlies Y Chromosome Stratification across Indonesia. Mol. Biol. Evol.27 (8), 1833–1844. 10.1093/molbev/msq063

  • CrossRef
  • Google Scholar

• 27

  KingC. L.BentleyR. A.TaylesN.ViðarsdóttirU. S.NowellG.MacphersonC. G. (2013). Moving Peoples, Changing Diets: Isotopic Differences Highlight Migration and Subsistence Changes in the Upper Mun River Valley, Thailand. J. Archaeological Sci.40, 1681–1688. 10.1016/J.JAS.2012.11.013

  • CrossRef
  • Google Scholar

• 28

  KlalesA. R.OusleyS. D.VollnerJ. M. (2012). A Revised Method of Sexing the Human Innominate Using Phenice's Nonmetric Traits and Statistical Methods. Am. J. Phys. Anthropol.149 (1), 104–114. 10.1002/ajpa.22102

  • CrossRef
  • Google Scholar

• 29

  KnappM.ClarkeA. C.HorsburghK. A.Matisoo-SmithE. A. (2012). Setting the Stage - Building and Working in an Ancient DNA Laboratory. Ann. Anat. - Anatomischer Anzeiger194 (1), 3–6. 10.1016/j.aanat.2011.03.008

  • CrossRef
  • Google Scholar

• 30

  KoA. M.-S.ChenC.-Y.FuQ.DelfinF.LiM.ChiuH.-L.et al (2014). Early Austronesians: into and Out of Taiwan. Am. J. Hum. Genet.94 (3), 426–436. 10.1016/j.ajhg.2014.02.003

  • CrossRef
  • Google Scholar

• 31

  KongQ.-P.YaoY.-G.SunC.BandeltH.-J.ZhuC.-L.ZhangY.-P. (2003). Phylogeny of East Asian Mitochondrial DNA Lineages Inferred from Complete Sequences. Am. J. Hum. Genet.73 (3), 671–676. 10.1086/377718

  • CrossRef
  • Google Scholar

• 32

  KutananW.KampuansaiJ.SrikummoolM.KangwanpongD.GhirottoS.BrunelliA.et al (2017). Complete Mitochondrial Genomes of Thai and Lao Populations Indicate an Ancient Origin of Austroasiatic Groups and Demic Diffusion in the Spread of Tai-Kadai Languages. Hum. Genet.136 (1), 85–98. 10.1007/s00439-016-1742-y

  • CrossRef
  • Google Scholar

• 33

  LiH.DurbinR. (2010). Fast and Accurate Long-Read Alignment with Burrows-Wheeler Transform. Bioinformatics26 (5), 589–595. 10.1093/bioinformatics/btp698

  • CrossRef
  • Google Scholar

• 34

  LiH.LiuZ.JamesN.LiX.HuZ.ShiH.et al (2021). Agricultural Transformations and Their Influential Factors Revealed by Archaeobotanical Evidence in Holocene Jiangsu Province, Eastern China. Front. Earth Sci.9, 661–684. 10.3389/feart.2021.661684

  • CrossRef
  • Google Scholar

• 35

  LiT. Y. (2018). A Study about the Garrisons’ Native Place at Zhangye Area in Bamboo Slips of Han Dynasty. Master’s thesis. Changchun: Jilin University.

  • Google Scholar

• 36

  LiX. (2021). Human Diets and its Influencing Factors during Han and Jin Periods in the Hexi Corridor and its Adjacent Areas. Dissertation. Lanzhou: Lanzhou University.

  • Google Scholar

• 37

  LiY.-C.YeW.-J.JiangC.-G.ZengZ.TianJ.-Y.YangL.-Q.et al (2019). River Valleys Shaped the Maternal Genetic Landscape of Han Chinese. Mol. Biol. Evol.36 (8), 1643–1652. 10.1093/molbev/msz072

  • CrossRef
  • Google Scholar

• 38

  LiuY. J. (2012). Study of Immigrants of Hexi in Han Dynasty. Master’s thesis. Lanzhou: Northwest Normal University.

  • Google Scholar

• 39

  LuoX. Q.DuP. X.WangL. X.ZhouB. Y.LiY. C.ZhengH. X.et al (2019). Uniparental Genetic Analyses Reveal the Major Origin of Fujian Tanka from Ancient Indigenous Daic Populations. Hum. Biol.91 (4), 257–277. 10.13110/humanbiology.91.4.05

  • CrossRef
  • Google Scholar

• 40

  MaM.DongG.JiaX.WangH.CuiY.ChenF. (2016). Dietary Shift after 3600 Cal Yr BP and its Influencing Factors in Northwestern China: Evidence from Stable Isotopes. Quat. Sci. Rev.145, 57–70. 10.1016/j.quascirev.2016.05.041

  • CrossRef
  • Google Scholar

• 41

  MaP.YangX.YanS.LiC.GaoS.HanB.et al (2021). Ancient Y-DNA with Reconstructed Phylogeny Provide Insights into the Demographic History of Paternal Haplogroup N1a2-F1360. J. Genet. Genomics48 (12), 1130–1133. 10.1016/j.jgg.2021.07.018

  • CrossRef
  • Google Scholar

• 42

  MeyerM.KircherM. (2010). Illumina Sequencing Library Preparation for Highly Multiplexed Target Capture and Sequencing. Cold Spring Harb Protoc.2010 (6), pdb.prot5448–prot5448. 10.1101/pdb.prot5448

  • CrossRef
  • Google Scholar

• 43

  NúñezL.GrosjeanM.CartajenaI. (2002). Human occupations and climate change in the Puna de Atacama, Chile. Science298, 821–824. 10.1126/science.1076449

  • CrossRef
  • Google Scholar

• 44

  OotaH.Settheetham-IshidaW.TiwawechD.IshidaT.StonekingM. (2001). Human mtDNA and Y-Chromosome Variation Is Correlated with Matrilocal versus Patrilocal Residence. Nat. Genet.29 (1), 20–21. 10.1038/ng711

  • CrossRef
  • Google Scholar

• 45

  PeltzerA.JägerG.HerbigA.SeitzA.KniepC.KrauseJ.et al (2016). EAGER: Efficient Ancient Genome Reconstruction. Genome Biol.17, 60. 10.1186/s13059-016-0918-z

  • CrossRef
  • Google Scholar

• 46

  PetragliaM. D.GroucuttH. S.GuagninM.BreezeP. S.BoivinN. (2020). Human Responses to Climate and Ecosystem Change in Ancient Arabia. Proc. Natl. Acad. Sci. USA117, 8263–8270. 10.1073/pnas.1920211117

  • CrossRef
  • Google Scholar

• 47

  PokhariaA. K.AgnihotriR.SharmaS.BajpaiS.NathJ.KumaranR. N.et al (2017). Altered Cropping Pattern and Cultural Continuation with Declined prosperity Following Abrupt and Extreme Arid Event at ∼4,200 Yrs BP: Evidence from an Indus Archaeological Site Khirsara, Gujarat, Western India. PloS one12 (10), e0185684. 10.1371/journal.pone.0185684

  • CrossRef
  • Google Scholar

• 48

  PolyakV. J.AsmeromY. (2001). Late Holocene Climate and Cultural Changes in the Southwestern United States. Science294, 148–151. 10.1126/science.1062771

  • CrossRef
  • Google Scholar

• 49

  PrestonG. W. (2011). From Nomadic Herder-hunters to Sedentary Farmers: The Relationship between Climate, Environment and Human Societies in the united arab emirates from the Neolithic to the Iron Age. Dissertation. Oxford: Oxford Brookes University.

  • Google Scholar

• 50

  PrestonG. W.ParkerA. G.WalkingtonH.LengM. J.HodsonM. J. (2012). From Nomadic Herder-hunters to Sedentary Farmers: the Relationship between Climate Change and Ancient Subsistence Strategies in South-Eastern Arabia. J. Arid Environments86, 122–130. 10.1016/j.jaridenv.2011.11.030

  • CrossRef
  • Google Scholar

• 51

  QiC. (1983). The Rise and Fall of Ancient Hexi. The Archaeology Team of Silk Road. Lanzhou: Gansu people's Publishing House.

  • Google Scholar

• 52

  RalfA.Montiel GonzálezD.ZhongK.KayserM. (2018). Yleaf: Software for Human Y-Chromosomal Haplogroup Inference from Next-Generation Sequencing Data. Mol. Biol. Evol.35 (5), 1291–1294. 10.1093/molbev/msy032

  • CrossRef
  • Google Scholar

• 53

  RenaudG.SlonV.DugganA. T.KelsoJ. (2015). Schmutzi: Estimation of Contamination and Endogenous Mitochondrial Consensus Calling for Ancient DNA. Genome Biol.16, 224. 10.1186/s13059-015-0776-0

  • CrossRef
  • Google Scholar

• 54

  SawyerS.RenaudG.ViolaB.HublinJ. J.GansaugeM. T.ShunkovM. V.et al (2015). Nuclear and Mitochondrial DNA Sequences from Two Denisovan Individuals. Proc. Nat. Acad. Sci. USA112 (51), 15696–15700. 10.1073/pnas.1519905112

  • CrossRef
  • Google Scholar

• 55

  Shao-QingW.Ruo-YuB.Bo-YanZ.Pan-XinD.ChangS.LiangC.et al (2019). China National DNA Martyry: a beacon of hope for the Martyrs' Coming home. J. Hum. Genet.64 (10), 1045–1047. 10.1038/s10038-019-0649-6

  • CrossRef
  • Google Scholar

• 56

  SherrattA. (1981). “Plough and Pastoralism: Aspects of the Secondary Products Revolution,” in Pattern of the Past: Studies in Honour of David Clarke. Editors HodderI.IsaacG.HammondN. (Cambridge: Cambridge University Press), 261–305.

  • Google Scholar

• 57

  ShiZ.ChenT.StorozumM. J.LiuF. (2019). Environmental and Social Factors Influencing the Spatiotemporal Variation of Archaeological Sites during the Historical Period in the Heihe River basin, Northwest China. Quat. Int.507, 34–42. 10.1016/j.quaint.2018.12.016

  • CrossRef
  • Google Scholar

• 58

  SiM. Q. (2002). Shiji. The Biographies of Dawan. Changsha: Yuelu Publishing House, 697–698. 10.1007/s00106-002-0694-9

  • CrossRef
  • Google Scholar

• 59

  SunN.MaP.-C.YanS.WenS.-Q.SunC.DuP.-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. 10.1080/03014460.2019.1632930

  • CrossRef
  • Google Scholar

• 60

  SunX.-f.WenS.-q.LuC.-q.ZhouB.-y.CurnoeD.LuH.-y.et al (2021). Ancient DNA and Multimethod Dating Confirm the Late Arrival of Anatomically Modern Humans in Southern China. Proc. Natl. Acad. Sci. USA118 (8), e2019158118. 10.1073/pnas.2019158118

  • CrossRef
  • Google Scholar

• 61

  SunZ.LiuS. (2014). Study of Hexi Diet from Unearthed Slips on the Han Dynasty. Gansu Soc. Sci.000 (006), 90–94.

  • Google Scholar

• 62

  ThorvaldsdottirH.RobinsonJ. T.MesirovJ. P. (2013). Integrative Genomics Viewer (IGV): High-Performance Genomics Data Visualization and Exploration. Brief. Bioinform.14, 178–192. 10.1093/bib/bbs017

  • CrossRef
  • Google Scholar

• 63

  TrejautJ. A.PoloniE. S.YenJ.-C.LaiY.-H.LooJ.-H.LeeC.-L.et al (2014). Taiwan Y-Chromosomal DNA Variation and its Relationship with Island Southeast Asia. BMC Genet.15 (1), 77–23. 10.1186/1471-2156-15-77

  • CrossRef
  • Google Scholar

• 64

  VolodkoN. V.StarikovskayaE. B.MazuninI. O.EltsovN. P.NaidenkoP. V.WallaceD. C.et al (2008). Mitochondrial Genome Diversity in Arctic Siberians, with Particular Reference to the Evolutionary History of Beringia and Pleistocenic Peopling of the Americas. Am. J. Hum. Genet.82 (5), 1084–1100. 10.1016/j.ajhg.2008.03.019

  • CrossRef
  • Google Scholar

• 65

  WeissensteinerH.PacherD.Kloss-BrandstätterA.ForerL.SpechtG.BandeltH.-J.et al (2016). HaploGrep 2: Mitochondrial Haplogroup Classification in the Era of High-Throughput Sequencing. Nucleic Acids Res.44 (W1), W58–W63. 10.1093/nar/gkw233

  • CrossRef
  • Google Scholar

• 66

  WenB.LiH.GaoS.MaoX.GaoY.LiF.et al (2005). Genetic Structure of Hmong-Mien Speaking Populations in East Asia as Revealed by mtDNA Lineages. Mol. Biol. Evol.22 (3), 725–734. 10.1093/molbev/msi055

  • CrossRef
  • Google Scholar

• 67

  WenB.LiH.LuD.SongX.ZhangF.HeY.et al (2004a). Genetic Evidence Supports Demic Diffusion of Han Culture. Nature431, 302–305. 10.1038/nature02878

  • CrossRef
  • Google Scholar

• 68

  WenB.XieX.GaoS.LiH.ShiH.SongX.et al (2004b). Analyses of Genetic Structure of Tibeto-Burman Populations Reveals Sex-Biased Admixture in Southern Tibeto-Burmans. Am. J. Hum. Genet.74 (5), 856–865. 10.1086/386292

  • CrossRef
  • Google Scholar

• 69

  WenS.-Q.TongX.-Z.LiH. (2016). Y-chromosome-based Genetic Pattern in East Asia Affected by Neolithic Transition. Quat. Int.426, 50–55. 10.1038/s10038-020-0775-110.1016/j.quaint.2016.03.027

  • CrossRef
  • Google Scholar

• 70

  WuQ.ChengH.-Z.SunN.MaP.-C.SunJ.YaoH.-B.et al (2020). Phylogenetic Analysis of the Y-Chromosome Haplogroup C2b-F1067, a Dominant Paternal Lineage in Eastern Eurasia. J. Hum. Genet.65, 823–829. 10.1038/s10038-020-0775-1

  • CrossRef
  • Google Scholar

• 71

  XuK.HuS. (2015). Population Data of Mitochondrial DNA HVS-I and HVS-II Sequences for 208 Henan Han Chinese. Leg. Med.17 (4), 287–294. 10.1016/j.legalmed.2015.02.003

  • CrossRef
  • Google Scholar

• 72

  YanS.WangC. C.ZhengH. X.WangW.QinZ. D.WeiL. H.et al (2014). Y Chromosomes of 40% Chinese Descend from Three Neolithic Super-grandfathers. PloS one9 (8), e105691. 10.1371/journal.pone.0105691

  • CrossRef
  • Google Scholar

• 73

  YanS.WangC. C.WangC.-C.LiH.LiS.-L.JinL. (2011). An Updated Tree of Y-Chromosome Haplogroup O and Revised Phylogenetic Positions of Mutations P164 and PK4. Eur. J. Hum. Genet.19 (9), 1013–1015. 10.1038/ejhg.2011.64

  • CrossRef
  • Google Scholar

• 74

  YangJ. (2015). Study on Food Culture in Hexi Areas of Han to Jin Dynasties from Archaeological Perspective (Lanzhou: Northwest Normal University). [Dissertation].

  • Google Scholar

• 75

  YangL.ShiZ.ZhangS.LeeH. F. (2020). Climate Change, Geopolitics, and Human Settlements in the Hexi Corridor over the Last 5,000 Years. Acta Geologica Sinica ‐ English Edition94, 612–623. 10.1111/1755-6724.14529

  • CrossRef
  • Google Scholar

• 76

  YangY.RenL.DongG.CuiY.LiuR.ChenG.et al (2019a). Economic Change in the Prehistoric Hexi Corridor (4800-2200bp), North-West China. Archaeometry61, 957–976. 10.1111/arcm.12464

  • CrossRef
  • Google Scholar

• 77

  YangY.ZhangS.OldknowC.QiuM.ChenT.LiH.et al (2019b). Refined Chronology of Prehistoric Cultures and its Implication for Re-evaluating Human-Environment Relations in the Hexi Corridor, Northwest China. Sci. China Earth Sci.62, 1578–1590. 10.1007/S11430-018-9375-4

  • CrossRef
  • Google Scholar

• 78

  YaoY.-G.KongQ.-P.BandeltH.-J.KivisildT.ZhangY.-P. (2002). Phylogeographic Differentiation of Mitochondrial DNA in Han Chinese. Am. J. Hum. Genet.70 (3), 635–651. 10.1086/338999

  • CrossRef
  • Google Scholar

• 79

  ZhangH. L. (2017). A Study on Han and Jin Tombs in the Luoyang Area (Zhengzhou: Zhengzhou University). [Dissertation].

  • Google Scholar

• 80

  ZhangM. (2021). Late Pleistocene and Early Holocene Microblade-Based Industries in Northeastern Asia: A Macroecological Approach to Foraging Societies. Oxford: British Archaeological Reports.

  • Google Scholar

• 81

  ZhangX.LiaoS.QiX.LiuJ.KampuansaiJ.ZhangH.et al (2015). Y-chromosome Diversity Suggests Southern Origin and Paleolithic Backwave Migration of Austro-Asiatic Speakers from Eastern Asia to the Indian Subcontinent. Sci. Rep.5 (1), 1–8. 10.1038/srep15486

  • CrossRef
  • Google Scholar

• 82

  ZhengB.GaoG. (2019). The Murals Unearthed from Wei, Jin and Tang Tombs in Gansu Province. Lanzhou: Gansu University Press.

  • Google Scholar

• 83

  ZhongH.ShiH.QiX.-B.XiaoC.-J.JinL.MaR. Z.et al (2010). Global Distribution of Y-Chromosome Haplogroup C Reveals the Prehistoric Migration Routes of African Exodus and Early Settlement in East Asia. J. Hum. Genet.55 (7), 428–435. 10.1038/jhg.2010.40

  • CrossRef
  • Google Scholar

• 84

  ZhouX.LiX.DodsonJ.ZhaoK. (2016). Rapid Agricultural Transformation in the Prehistoric Hexi Corridor, China. Quat. Int.426, 33–41. 10.1016/j.quaint.2016.04.021

  • CrossRef
  • Google Scholar