Jianxue Xiong1†
Panxin Du1†Guoke Chen2
Yichen Tao1Boyan Zhou3Yishi Yang2Hui Wang4,5Yao Yu4Xin Chang4Edward Allen4Chang Sun1Juanjuan Zhou1Yetao Zou1Yiran Xu4Hailiang Meng1Jingze Tan1*
Hui Li1*
Shaoqing Wen4,5*- 1Ministry of Education Key Laboratory of Contemporary Anthropology, Department of Anthropology and Human Genetics, School of Life Sciences, Fudan University, Shanghai, China
- 2Institute of Cultural Relics and Archaeology in Gansu Province, Lanzhou, China
- 3Division of Biostatistics, Department of Population Health, School of Medicine, New York University, New York, NY, United States
- 4Institute of Archaeological Science, Fudan University, Shanghai, China
- 5Center for the Belt and Road Archaeology and Ancient Civilizations (BRAAC), Fudan University, Shanghai, China
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 14C 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. The geographical location of the Heishuiguo site and the Hexi Corridor.
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.
TABLE 1. Ancient individuals sampled in this study.
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. The phylogenetic tree of 485 Y-SNPs in this study.
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. The phylogenetic relationship of Y-chromosome haplogroups in this study and their haplogroup-based frequencies in the sampled populations (Early Han, the early Han group; Late Han, the late Han group). Marker names are shown along the branches, and haplogroup names are shown to the right, based on ISOGG Y-DNA Haplogroup Tree 2019. Asterisks distinguish potentially paraphyletic undefined subgroups from recognized haplogroups. The markers in red are key Y-SNPs. a-f: These markers are not designed in Y-SNP panel but very common in Y phylogenetic trees.
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. The phylogenetic relationship of mtDNA haplogroups surveyed in this study and their haplogroup-based frequencies in the sampled populations. Haplogroup names are shown to the right, according to the PhyloTree mtDNA tree Build 17.
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. Principal component plots of Heishuiguo population and reference populations for Y-chromosome haplogroups (A) and mtDNA haplogroups (B).
FIGURE 6. The genetic distance (Fst) heatmap plot of Heishuiguo population and reference populations for Y-chromosome.
FIGURE 7. The genetic distance (Fst) heatmap plot of Heishuiguo population and reference populations for mtDNA.
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. Genetic contour maps of Heishuiguo population and reference populations for Y-chromosome with early Han group (A), Y-chromosome with late Han group (B), mtDNA with early Han group (C), mtDNA with late Han group (D).
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. Immigrants from 20 counties in middle and lower Yellow River to Hexi Corridor according to the contents of unearthed bamboo slips.
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.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Ethics Statement
Approval for the use of ancient human individuals was curated by co-authors and obtained with permission from the respective provincial archaeology institutes or universities that managed the samples. The permission and oversight were also provided by the Ethics Committee of Fudan University of Life Sciences to study their ancient genomes.
Author Contributions
JT, HL, and SW designed this study. HL and SW supervised the study. GC, YiY, and HW provided materials and resources. HM and JT collected the samples. YiY performed archaeological data analysis. JX, PD, JZ, YZ, and YX performed genetic laboratory work. YT, BZ, CS, XC, and YaY performed genetic data analysis. JX and PD integrated the genetic data. JX, PD, and SW wrote the paper, with contributions from EA.
Funding
This work was supported by the National Key R&D Program of China (2020YFE0201600 and 2020YFC1521607), B&R Joint Laboratory of Eurasian Anthropology (18490750300), National Natural Science Foundation of China (31771325 and 32070576), Major Research Program of National Natural Science Foundation of China (91731303), Major Project of National Social Science Foundation of China (20&ZD212), Shanghai Municipal Science and Technology Major Project (2017SHZDZX01), the 111 Project (B13016), and European Research Council (ERC) grant to Dan Xu (ERC-2019-ADG-883700-TRAM).
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.
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.2022.827277/full#supplementary-material
Supplementary Figure S1 | Ancient DNA damage patterns. Nucleotide misincorporation patterns cause by C deamination in ancient DNA sequences.
Supplementary Figure S2 | Principal component plot of Heishuiguo population and published ancient populations for mtDNA haplogroups.
Supplementary Figure S3 | The genetic distance (Fst) heatmap plot of Heishuiguo population and published ancient populations for mtDNA.
References
An, C.-B., Tang, L., Barton, L., and Chen, F.-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. doi:10.1016/j.yqres.2005.02.004
Bonsall, C., Macklin, M. G., Anderson, D. E., and Payton, R. W. (2002). Climate Change and the Adoption of Agriculture in north-west Europe. Eur. J. Archaeol. 5 (1), 9–23. doi:10.1179/eja.2002.5.1.9
Buikstra, J. E. (1994). Standards for Data Collection from Human Skeletal Remains. Fayetteville, Arkansas: Arkansas archeological survey. Research series no. 44.
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
Calafell, F., and Larmuseau, M. H. D. (2017). The Y Chromosome as the Most Popular Marker in Genetic Genealogy Benefits Interdisciplinary Research. Hum. Genet. 136, 559–573. doi:10.1007/s00439-016-1740-0
Chen, G. K., Yang, Y., and Liu, B. (2019). Ganzhou, Zhangye: The Excavation Report of Han Dynasty Cemetery in Heishuiguo Site (Volume II). Lanzhou: Gansu Education Publishing House.
Cheung, C., Zhang, H., Hepburn, J. C., Yang, D. Y., and Richards, M. P. (2019). Stable Isotope and Dental Caries Data Reveal Abrupt Changes in Subsistence Economy in Ancient China in Response to Global Climate Change. PLOS ONE 14 (7), e0218943. doi:10.1371/journal.pone.0218943
Cohen, M. N. (1975). Archaeological Evidence for Population Pressure in Pre-agricultural Societies. Am. Antiq. 40 (4), 471–475. doi:10.2307/279335
Dabney, J., Knapp, M., Glocke, I., Gansauge, M.-T., Weihmann, A., Nickel, B., 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. doi:10.1073/PNAS.1314445110
d’Alpoim Guedes, J., Jin, G., and Bocinsky, R. 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 one 10 (6), e0130430. doi:10.1371/journal.pone.0130430
Derenko, M., Malyarchuk, B., Grzybowski, T., Denisova, G., Rogalla, U., Perkova, M., et al. (2010). Origin and post-glacial Dispersal of Mitochondrial DNA Haplogroups C and D in Northern Asia. PloS one 5 (12), e15214. doi:10.1371/journal.pone.0015214
Gansu Bamboo Slips Museum (GBSM) (2013). Han Dynasty Wooden Slips from Jianshuijinguan (Third). Shanghai: Zhongxi Book Company, 52.
Gao, H., Zhu, C., and Xu, W. (2007). Environmental Change and Cultural Response Around 4200 Cal. Yr BP in the Yishu River Basin, Shandong. J. Geogr. Sci. 17, 285–292. doi:10.1007/s11442-007-0285-5
Ge, J. X., Wu, S. D., and Chao, S. J. (1997). Zhongguo Yimin Shi (The Migration History of China). Fuzhou: Fuzhou: Fujian People’s Publishing House.
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.
GPICRA (Gansu Provincial Institute of Cultural Relics and Archaeology) (2001). Excavation Report of Xigangchaiwangang Tomb in Yongchang. Lanzhou: Gansu People’s Publishing House.
Haug, G. H., Günther, D., Peterson, L. C., Sigman, D. M., Hughen, K. A., and Aeschlimann, B. (2003). Climate and the Collapse of Maya Civilization. Science 299, 1731–1735. doi:10.1126/science.1080444
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. arXiv:1504.06463.
Ilumäe, A.-M., Reidla, M., Chukhryaeva, M., Järve, 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
Jia, X., Sun, Y., Wang, L., Sun, W., Zhao, Z., Lee, H. 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 Holocene 26 (5), 781–789. doi:10.1177/0959683615618262
Jin, H.-J., Tyler-Smith, C., and Kim, W. (2009). The Peopling of Korea Revealed by Analyses of Mitochondrial DNA and Y-Chromosomal Markers. PloS one 4 (1), e4210. doi:10.1371/journal.pone.0004210
Jónsson, H., Ginolhac, A., Schubert, M., Johnson, P. L. F., and Orlando, L. (2013). mapDamage2.0: Fast Approximate Bayesian Estimates of Ancient DNA Damage Parameters. Bioinformatics (Oxford, England) 29 (13), 1682–1684. doi:10.1093/bioinformatics/btt193
Karafet, T. M., Hallmark, B., Cox, M. P., Sudoyo, H., Downey, S., Lansing, J. S., et al. (2010). Major East-West Division Underlies Y Chromosome Stratification across Indonesia. Mol. Biol. Evol. 27 (8), 1833–1844. doi:10.1093/molbev/msq063
King, C. L., Bentley, R. A., Tayles, N., Viðarsdóttir, U. S., Nowell, G., and Macpherson, C. 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. doi:10.1016/J.JAS.2012.11.013
Klales, A. R., Ousley, S. D., and Vollner, J. 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. doi:10.1002/ajpa.22102
Knapp, M., Clarke, A. C., Horsburgh, K. A., and Matisoo-Smith, E. A. (2012). Setting the Stage - Building and Working in an Ancient DNA Laboratory. Ann. Anat. - Anatomischer Anzeiger 194 (1), 3–6. doi:10.1016/j.aanat.2011.03.008
Ko, A. M.-S., Chen, C.-Y., Fu, Q., Delfin, F., Li, M., Chiu, H.-L., et al. (2014). Early Austronesians: into and Out of Taiwan. Am. J. Hum. Genet. 94 (3), 426–436. doi:10.1016/j.ajhg.2014.02.003
Kong, Q.-P., Yao, Y.-G., Sun, C., Bandelt, H.-J., Zhu, C.-L., and Zhang, Y.-P. (2003). Phylogeny of East Asian Mitochondrial DNA Lineages Inferred from Complete Sequences. Am. J. Hum. Genet. 73 (3), 671–676. doi:10.1086/377718
Kutanan, W., Kampuansai, J., Srikummool, M., Kangwanpong, D., Ghirotto, S., Brunelli, A., 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. doi:10.1007/s00439-016-1742-y
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., Liu, Z., James, N., Li, X., Hu, Z., Shi, H., 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. doi:10.3389/feart.2021.661684
Li, T. Y. (2018). A Study about the Garrisons’ Native Place at Zhangye Area in Bamboo Slips of Han Dynasty. Master’s thesis. Changchun: Jilin University.
Li, X. (2021). Human Diets and its Influencing Factors during Han and Jin Periods in the Hexi Corridor and its Adjacent Areas. Dissertation. Lanzhou: Lanzhou University.
Li, Y.-C., Ye, W.-J., Jiang, C.-G., Zeng, Z., Tian, J.-Y., Yang, L.-Q., et al. (2019). River Valleys Shaped the Maternal Genetic Landscape of Han Chinese. Mol. Biol. Evol. 36 (8), 1643–1652. doi:10.1093/molbev/msz072
Liu, Y. J. (2012). Study of Immigrants of Hexi in Han Dynasty. Master’s thesis. Lanzhou: Northwest Normal University.
Luo, X. Q., Du, P. X., Wang, L. X., Zhou, B. Y., Li, Y. C., Zheng, H. 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. doi:10.13110/humanbiology.91.4.05
Ma, M., Dong, G., Jia, X., Wang, H., Cui, Y., and Chen, F. (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. doi:10.1016/j.quascirev.2016.05.041
Ma, P., Yang, X., Yan, S., Li, C., Gao, S., Han, B., et al. (2021). Ancient Y-DNA with Reconstructed Phylogeny Provide 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
Meyer, M., and Kircher, M. (2010). Illumina Sequencing Library Preparation for Highly Multiplexed Target Capture and Sequencing. Cold Spring Harb Protoc. 2010 (6), pdb.prot5448–prot5448. doi:10.1101/pdb.prot5448
Núñez, L., Grosjean, M., and Cartajena, I. (2002). Human occupations and climate change in the Puna de Atacama, Chile. Science 298, 821–824. doi:10.1126/science.1076449
Oota, H., Settheetham-Ishida, W., Tiwawech, D., Ishida, T., and Stoneking, M. (2001). Human mtDNA and Y-Chromosome Variation Is Correlated with Matrilocal versus Patrilocal Residence. Nat. Genet. 29 (1), 20–21. doi:10.1038/ng711
Peltzer, A., Jäger, G., Herbig, A., Seitz, A., Kniep, C., Krause, J., et al. (2016). EAGER: Efficient Ancient Genome Reconstruction. Genome Biol. 17, 60. doi:10.1186/s13059-016-0918-z
Petraglia, M. D., Groucutt, H. S., Guagnin, M., Breeze, P. S., and Boivin, N. (2020). Human Responses to Climate and Ecosystem Change in Ancient Arabia. Proc. Natl. Acad. Sci. USA 117, 8263–8270. doi:10.1073/pnas.1920211117
Pokharia, A. K., Agnihotri, R., Sharma, S., Bajpai, S., Nath, J., Kumaran, R. 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 one 12 (10), e0185684. doi:10.1371/journal.pone.0185684
Polyak, V. J., and Asmerom, Y. (2001). Late Holocene Climate and Cultural Changes in the Southwestern United States. Science 294, 148–151. doi:10.1126/science.1062771
Preston, G. 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.
Preston, G. W., Parker, A. G., Walkington, H., Leng, M. J., and Hodson, M. J. (2012). From Nomadic Herder-hunters to Sedentary Farmers: the Relationship between Climate Change and Ancient Subsistence Strategies in South-Eastern Arabia. J. Arid Environments 86, 122–130. doi:10.1016/j.jaridenv.2011.11.030
Qi, C. (1983). The Rise and Fall of Ancient Hexi. The Archaeology Team of Silk Road. Lanzhou: Gansu people's Publishing House.
Ralf, A., Montiel González, D., Zhong, K., and Kayser, M. (2018). Yleaf: Software for Human Y-Chromosomal Haplogroup Inference from Next-Generation Sequencing Data. Mol. Biol. Evol. 35 (5), 1291–1294. doi:10.1093/molbev/msy032
Renaud, G., Slon, V., Duggan, A. T., and Kelso, J. (2015). Schmutzi: Estimation of Contamination and Endogenous Mitochondrial Consensus Calling for Ancient DNA. Genome Biol. 16, 224. doi:10.1186/s13059-015-0776-0
Sawyer, S., Renaud, G., Viola, B., Hublin, J. J., Gansauge, M. T., Shunkov, M. V., et al. (2015). Nuclear and Mitochondrial DNA Sequences from Two Denisovan Individuals. Proc. Nat. Acad. Sci. USA 112 (51), 15696–15700. doi:10.1073/pnas.1519905112
Shao-Qing, W., Ruo-Yu, B., Bo-Yan, Z., Pan-Xin, D., Chang, S., Liang, C., et al. (2019). China National DNA Martyry: a beacon of hope for the Martyrs' Coming home. J. Hum. Genet. 64 (10), 1045–1047. doi:10.1038/s10038-019-0649-6
Sherratt, A. (1981). “Plough and Pastoralism: Aspects of the Secondary Products Revolution,” in Pattern of the Past: Studies in Honour of David Clarke. Editors I. Hodder, G. Isaac, and N. Hammond (Cambridge: Cambridge University Press), 261–305.
Shi, Z., Chen, T., Storozum, M. J., and Liu, F. (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. doi:10.1016/j.quaint.2018.12.016
Si, M. Q. (2002). Shiji. The Biographies of Dawan. Changsha: Yuelu Publishing House, 697–698. doi:10.1007/s00106-002-0694-9
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
Sun, X.-f., Wen, S.-q., Lu, C.-q., Zhou, B.-y., Curnoe, D., Lu, H.-y., et al. (2021). Ancient DNA and Multimethod Dating Confirm the Late Arrival of Anatomically Modern Humans in Southern China. Proc. Natl. Acad. Sci. USA 118 (8), e2019158118. doi:10.1073/pnas.2019158118
Sun, Z., and Liu, S. (2014). Study of Hexi Diet from Unearthed Slips on the Han Dynasty. Gansu Soc. Sci. 000 (006), 90–94.
Thorvaldsdottir, H., Robinson, J. T., and Mesirov, J. P. (2013). Integrative Genomics Viewer (IGV): High-Performance Genomics Data Visualization and Exploration. Brief. Bioinform. 14, 178–192. doi:10.1093/bib/bbs017
Trejaut, J. A., Poloni, E. S., Yen, J.-C., Lai, Y.-H., Loo, J.-H., Lee, C.-L., et al. (2014). Taiwan Y-Chromosomal DNA Variation and its Relationship with Island Southeast Asia. BMC Genet. 15 (1), 77–23. doi:10.1186/1471-2156-15-77
Volodko, N. V., Starikovskaya, E. B., Mazunin, I. O., Eltsov, N. P., Naidenko, P. V., Wallace, D. 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. doi:10.1016/j.ajhg.2008.03.019
Weissensteiner, H., Pacher, D., Kloss-Brandstätter, A., Forer, L., Specht, G., Bandelt, H.-J., et al. (2016). HaploGrep 2: Mitochondrial Haplogroup Classification in the Era of High-Throughput Sequencing. Nucleic Acids Res. 44 (W1), W58–W63. doi:10.1093/nar/gkw233
Wen, B., Li, H., Gao, S., Mao, X., Gao, Y., Li, F., 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. doi:10.1093/molbev/msi055
Wen, B., Li, H., Lu, D., Song, X., Zhang, F., He, Y., et al. (2004a). Genetic Evidence Supports Demic Diffusion of Han Culture. Nature 431, 302–305. doi:10.1038/nature02878
Wen, B., Xie, X., Gao, S., Li, H., Shi, H., Song, X., 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. doi:10.1086/386292
Wen, S.-Q., Tong, X.-Z., and Li, H. (2016). Y-chromosome-based Genetic Pattern in East Asia Affected by Neolithic Transition. Quat. Int. 426, 50–55. doi:10.1038/s10038-020-0775-110.1016/j.quaint.2016.03.027
Wu, Q., Cheng, H.-Z., Sun, N., Ma, P.-C., Sun, J., Yao, H.-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. doi:10.1038/s10038-020-0775-1
Xu, K., and Hu, S. (2015). Population Data of Mitochondrial DNA HVS-I and HVS-II Sequences for 208 Henan Han Chinese. Leg. Med. 17 (4), 287–294. doi:10.1016/j.legalmed.2015.02.003
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
Yan, S., Wang, C. C., Wang, C.-C., Li, H., Li, S.-L., and Jin, L. (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. doi:10.1038/ejhg.2011.64
Yang, J. (2015). Study on Food Culture in Hexi Areas of Han to Jin Dynasties from Archaeological Perspective (Lanzhou: Northwest Normal University). [Dissertation].
Yang, L., Shi, Z., Zhang, S., and Lee, H. F. (2020). Climate Change, Geopolitics, and Human Settlements in the Hexi Corridor over the Last 5,000 Years. Acta Geologica Sinica ‐ English Edition 94, 612–623. doi:10.1111/1755-6724.14529
Yang, Y., Ren, L., Dong, G., Cui, Y., Liu, R., Chen, G., et al. (2019a). Economic Change in the Prehistoric Hexi Corridor (4800-2200bp), North-West China. Archaeometry 61, 957–976. doi:10.1111/arcm.12464
Yang, Y., Zhang, S., Oldknow, C., Qiu, M., Chen, T., Li, H., 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. doi:10.1007/S11430-018-9375-4
Yao, Y.-G., Kong, Q.-P., Bandelt, H.-J., Kivisild, T., and Zhang, Y.-P. (2002). Phylogeographic Differentiation of Mitochondrial DNA in Han Chinese. Am. J. Hum. Genet. 70 (3), 635–651. doi:10.1086/338999
Zhang, H. L. (2017). A Study on Han and Jin Tombs in the Luoyang Area (Zhengzhou: Zhengzhou University). [Dissertation].
Zhang, M. (2021). Late Pleistocene and Early Holocene Microblade-Based Industries in Northeastern Asia: A Macroecological Approach to Foraging Societies. Oxford: British Archaeological Reports.
Zhang, X., Liao, S., Qi, X., Liu, J., Kampuansai, J., Zhang, H., 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. doi:10.1038/srep15486
Zheng, B., and Gao, G. (2019). The Murals Unearthed from Wei, Jin and Tang Tombs in Gansu Province. Lanzhou: Gansu University Press.
Zhong, H., Shi, H., Qi, X.-B., Xiao, C.-J., Jin, L., Ma, R. 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. doi:10.1038/jhg.2010.40
Keywords: Hexi Corridor, subsistence strategy, sex-biased admixture, Y chromosome, mitogenome, ancient DNA
Citation: Xiong J, Du P, Chen G, Tao Y, Zhou B, Yang Y, Wang H, Yu Y, Chang X, Allen E, Sun C, Zhou J, Zou Y, Xu Y, Meng H, Tan J, Li H and Wen S (2022) Sex-Biased Population Admixture Mediated Subsistence Strategy Transition of Heishuiguo People in Han Dynasty Hexi Corridor. Front. Genet. 13:827277. doi: 10.3389/fgene.2022.827277
Received: 01 December 2021; Accepted: 10 January 2022;
Published: 10 March 2022.
Edited by:
Marc Via, University of Barcelona, SpainReviewed by:
Levon Yepiskoposyan, Armenian National Academy of Sciences, ArmeniaJiang Huang, Guizhou Medical University, China
Copyright © 2022 Xiong, Du, Chen, Tao, Zhou, Yang, Wang, Yu, Chang, Allen, Sun, Zhou, Zou, Xu, Meng, Tan, Li and Wen. 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: Jingze Tan, jztan@fudan.edu.cn; Hui Li, lihui.fudan@gmail.com; Shaoqing Wen, wenshaoqing@fudan.edu.cn
†These authors have contributed equally to this work