Article Open Access
Received: 14 March 2025 Accepted: 10 April 2025 Published: 18 April 2025
© 2025 The authors. This is an open access article under the Creative Commons Attribution 4.0 International License (https://creativecommons.org/licenses/by/4.0/).
The Okunev culture, primarily located in South Siberia, dates back to Afanasievo and Andronovo cultures [1]. The Okunev people had advanced animal husbandry, metallurgy, and unique burial practices and were known for their mysterious carvings and art. Scholars believe these cultural elements reflect a mysterious spiritual world and belief system [2,3]. Okunev cultural elements have been found in Central Asia, North China, and the Indus Valley Civilization, indicating its broad influence on Eurasian Bronze Age cultures [3]. However, the development and spread of the Okunev people remain unclear.
Ancient DNA research has offered insights into the origin and evolution of the Okunev culture. The ancient Okunev people can be seen as a blend of the indigenous Mongoloid population in South Siberia and a certain proportion of Indo-European- related ancestral components [4]. Their maternal lines include East Eurasian haplogroups such as A, A8, A8a, A8a1, A10, C5, C5c, D4, D4b2a2a, and West Eurasian ones like H, U5, H15b1, H1bv1, H6a, H6a1b, J [5,6]. Their paternal lines consist of Q-L330-YP761, R1b-Z2103, Q1b-L330-BZ180, Q-L940, Q-M25, N-Tat [5]. Scholars consider the presence of paternal line R1b-Z2103 as an indicator of the influence of the earlier Afanasievo culture, while N-Tat reflects the contemporaneous Uralic-speaking population’s impact [4,5]. Thus, the diverse branches of Q-M242 can be seen as the primary paternal type of the ancient Okunev people, with Q-L330-YP761 being the dominant one.
This paper will use modern-human sequences to explore the origin and evolution of the ancient Okunev people’s primary paternal lineage (Q-M242). The main aims are threefold. First, to track the flow of this lineage in later ancient populations. Second, to analyze the evolution and migration of the Okunev people’s descendants. Third, to identify the present-day descendants of the ancient Okunev people.
Samples and sequencing: Blood or saliva samples were collected from unrelated healthy male individuals from Chinese populations. All participants provided signed consent prior to participating in the study. Before signing the informed consent form, the contents of the form were clearly explained to the participants in either Chinese or their native language (like the language of minority). The study and sample collection procedures were reviewed and approved by the Medical Ethics Committee of Fudan University Approval NO.15137, adhering to the ethical principles of the 2013 Declaration of Helsinki by the World Medical Association.
Genomic DNA was extracted using the DP-318 Kit 128 (Tiangen Biotechnology, Beijing, China), following the manufacturer’s protocol. Firstly, routine Y-SNP tests were conducted to determine the general haplogroup of all male samples. DNA extracted from twenty-five Q-L330-YP761 male individuals’ biological samples was sent for next-generation Y-chromosome sequencing using the Illumina HiSeq2000 platform (Illumina, San Diego, CA, USA). Besides, the genome and information of twenty-five ancient individuals of haplogroup Q-L330-YP761 were extracted from previous publications [4,6,7,8,9,10,11]. Also, the variant situation of 10 haplogroup samples was adopted from publicly available resources, like previous literature [12] and the website of the human phylogenetic tree https://www.yfull.com/tree/Q-Y18330/ (accessed on 10 March 2025).
Data availability statement: Following the regulations of the Human Genetic Resources Administration of China (HGRAC), the raw sequence data reported in this paper are available at the request of the corresponding author.
Data analysis: In total, we summarized a final dataset of 60 samples, including 25 newly generated sequences (from the research group’s internal database, unpublished data, these individuals were preliminarily classified as haplogroup Q-L330-YP761), 25 sequences of ancient individuals sourced from published literature across distinct archaeological sites (detailed in Supplementary Table S1, all categorized as Q-L330-YP761), and 10 modern individuals from public available resources (all belonging to Q-L330-YP761 haplogroup). Read mapping and SNP calling from next-generation sequencing data were conducted using standard procedures (BWA and SAMtools) and the human reference genome sequence hg38 [13,14]. The final SNP list of 25 newly generated sequences was used to construct a phylogenetic tree and age estimation analyses. Bayesian evolutionary analyses were conducted using BEAST (v.2.4.3) [15]. To calculate divergence times in the phylogenetic tree, a point mutation rate of 0.74 × 10−9 per site per year was applied [16], which was very similar to the rate of 0.78 × 10−9 per site per year reported in analyses of the ancient genome of the 45,000-year-old Ust’Ishim male [17]. We refer to the protocol used in our previous publications for data processing and Bayesian evolutionary analyses [18,19,20]. Detailed sample information can be found in Supplementary Table S1.
In the next step, we placed the ancient available modern individuals onto the revised phylogenetic tree according to the variants of ancient genomes. Due to the low coverage of genomes of some ancient samples, the mutation situation of some SNPs of some samples is unclear. Therefore, some samples may be placed on a root location when they actually belong to a downstream sub-branch.
A schematic phylogenetic tree labeled ancient DNA is shown in Figure 1. The detailed tree can be found in Supplementary Table S1. We presented clear distribution patterns for sub-branches of haplogroup Q-L330-YP761 across Eurasia (Figure 2). The possible migration routes of sub-branches of haplogroup Q-L330-YP761 are shown in Figure 3.
The haplogroups provided by the original ancient DNA literature are general ones. Therefore, we re-analyzed the data and determined the most downstream sub-branch for ancient genomes (Table 1). Ancient individuals of the same haplogroup with Okunev samples were also summarized from public literature to trace the origin and descendants of Okunev lineages [4,5]. Haplogroup Q-L330-YP761 was dominant among remains from the core sites of Okunev culture. It was also the only paternal lineage found at the Tenisei Early Bronze Age site, which was proposed as the direct predecessor of Okunev culture. Therefore, we discussed this haplogroup in the next section.
Table 1. list of revised Y-chromosome haplogroup of ancient DNA.
ID | Location |
Y
Haplogroup |
Mt
Haplogroup |
Year |
---|---|---|---|---|
RISE664 | Russia | Q-YP761 | A8a1 | 4272 ybp |
RISE662 | Russia | Q-YP761 | H6a1b | 4063 ybp |
RISE718 | Russia | Q-YP761 | C5c | 4434 ybp |
DA73 ERS2374338 | Kyrgyzstan | Q-YP761 | G2a1 | 4000 ybp |
UKR109 | Ukraine | Q-YP761 | H5a1 | 2400 ybp |
DA96 ERS2374351 | Kyrgyzstan | Q-L332 | D4j5 | 2000 ybp |
DA161 ERS2374384 | Russia | Q-L332 | X2f | 2300 ybp |
I6231 | Mongolia | Q-L332 | F1b1b | 2300 ybp |
IS2 | Russia | Q-L332 | HV-CRS | |
DA221 ERS2374408 | Kazakhstan | Q-BZ425 | J1b1a1 | |
ERR5211315 | Mongolia | Q-YP4547 | 2180 ybp | |
ALN003.A0101 ERS6072385 | Kyrgyzstan | Q-YP4547 | U5a1g | 1700 ybp |
I6230 | Mongolia | Q-YP4547 | U2e1h | 2200 ybp |
UKR113 | Ukraine | Q-YP4547 | HV0a | 2600 ybp |
ERR5211347 | Mongolia | Q-Y145452 | 2200 ybp | |
CHN003 | Mongolia | Q-Y145452 | U4a1 | 2300 ybp |
CHN007 | Mongolia | Q-Y145452 | U2e1h | Mongolia |
DEL001 | Mongolia | Q-Y145452 | D4j11 | 2000 ybp |
SKT008 | Mongolia | Q-Y145452 | Hz | |
I7029 | Mongolia | Q-Y145452 | U4a1 | 2300 ybp |
CHN016 | Mongolia | Q-FTA37701 | A8a | |
ERR5211222 | Mongolia | Q-TY520082 | 2250 ybp | |
ERR5211316 | Mongolia | Q-FTA37701 | 2200 ybp | |
I6232 ERR5211317 | Mongolia | Q-TY520082 | 2250 ybp | |
I12970 | Mongolia | Q-TY520082 | A8a | 2300 ybp |
Haplogroup R1b-Z2103 was found in two samples (Kh7, RISE675) from the Uybat V site in Khakasia, South Siberia [4,5]. This haplogroup is one of the dominant paternal lineages of earlier Afanasievo people from the same region [5,21,22]. Therefore, it is clear that this paternal lineage, along with those west Eurasian mtDNA lineages, are results of admixture with local Afanasievo ancestry.
Haplogroup N-Tat was found in three samples (Kh12, Kh13, and Kh15) from the Itkol’II site and the Uybat-Charkov site in Khakasia, South Siberia. This haplogroup is one of the founder paternal lineages of Uralic-speaking populations and Turkic-speaking populations [21,22]. The duration of the Okunev culture (2500–1700 BC) was roughly concurrent with the Seima-Turbino phenomenon e (2200–1900 BC), which underwent significant expansion in the northern part of Eurasia during this period [23]. Therefore, this lineage among Okunev people is also likely a result of admixture with surrounding Uralic-speaking populations [24,25].
One individual of haplogroup Q1b-L330-BZ180 was found from Syda 5, Tumen site in Russia, South Siberia. The lineage was dominant among Early Bronze Age remains from the Cis-Baikal region and later in many sites across the Mongolian Plateau [8,11,26,27]. The appearance of this lineage among Okunev people may result from admixture. Four individuals related to Q-L940-Y2679 were found at the Verkhni Askiz site in Russia South Siberia [4]. They may belong to an extinct or a minor root lineage.
One individual of Q-M25-L712-YP832 was found in two samples (Kh7, RISE675) from the Uybat V site in Khakasia, South Siberia [4]. The earliest sample of this haplogroup, a 7500-year-old individual FRS002, was found in the Altai-Sayan region, South Siberia [28]. The haplogroup has also been found in large numbers among ancient Karasuk people [4], ancient Xiongnu [10], ancient Türk people, and ancient Avar people in Hungary [8]. This paternal lineage is likely flourishing among nomadic populations across the Eurasian Steppe.
In general, the paternal gene pool of ancient Okunev people was formed as a result of the rising of lineages (like Q-L330-YP761, Q-L940-Y2679, and Q-M25-L712-YP832) and admixture from surrounding people (like R1b-Z2103, N-Tat, and Q1b-L330-BZ180).
Previous studies have found many individuals of haplogroup Q-L330-YP761 from remains of Okunev culture and other archaeological cultures from the central part of the Eurasian Steppe at a later age (Table 1, Figure 1, and Table S1). The Okunev individuals’ variants suggested they are the direct predecessor of downstream sub-branch Q-L330-YP761-YP780, L332, which experienced continuous expansion about 4000 years ago. Samples of two root sub-linages, Y81743 and MF59110, were found among Han populations.
The revised phylogenetic tree suggested that sub-branch BZ425 is like a paternal lineage of ancient Xiongnu and Türk people. Samples of sub-lineage BZ513 were found among individuals from Kazakhstan, Hungary, and Chinese populations. The age of BZ513 (~2000 years ago) surprisingly coincides with the westward migration of ancient Xiongnu tribes.
The East-West dichotomous structure was also found among the topology of sub-lineage YP4564. Sub-lineage L329 was exclusively found in Europe while sub-lineage MF57829 was only found in China until now. This paternal is also likely the result of long-distance migration across the Eurasian Steppe.
Many Y145421 individuals were found among remains from the Chandman site of Uyuk-Sagly culture (900–300 BC) and several sites of early Xiongnu in the same region. According to the studies of archaeology and history, the northwest part of Mongolia and South Siberia became parts of the Xiongnu Empire after the rising of the Xiongnu people. Two modern individuals of this lineage were found in Tajik and Russia, suggesting a trace of westward migration.
Current aDNA research shows that Eurasian steppe populations during the Copper and Bronze Ages had highly diverse genetic lineages. In the late Paleolithic and Neolithic, NEA populations from the middle and lower Amur River diffused westward to the Baikal region and the central and eastern Mongolian Plateau [8,11]. Ancestral Uralic-speaking populations rapidly moved through the northern Eurasian taiga [29]. Between the Altai and Baikal, ANE populations differentiated rapidly, producing many genetic lineages and descendant groups [28]. These population movements and mixtures were very complex. Existing aDNA research has gradually revealed this demographic history [8,11,28,29]. This is reflected in the paternal genetic structure by the emergence of dozens of downstream lineages of C-L1373, Q-M242, and N-Tat and their extensive spread across northern Eurasia and the Eurasian steppe [18,21,22,30,31].
During the Copper and Bronze Age, the eastward expansion of Indo-European-related groups brought significant population movements and mixtures. West and East Eurasian lineages blended in Central Asia, shaping the region’s genetic makeup. Strong nomadic groups emerged in the eastern Eurasian steppe after the mid-Bronze Age (around 3500 BP), later establishing nomadic empires and shifting the population’s direction westward. This study’s detailed analysis of the ancient Okunev people’s paternal genetic structure supports the rise and spread of local lineages in the eastern Eurasian steppe, offering a fine-grained example of Eurasian steppe populations’ migration and genetic evolution.
The accumulation of aDNA studies has made the genetic structure of ancient archaeological cultures clearer over space and time. Meanwhile, the genetic makeup of modern populations has also become more evident. Therefore, integrating ancient and modern DNA research could be an effective way to explore the population history of Eurasian steppe peoples since ancient times.
In this paper, we explore the formation of the Okunev culture and the spread of its people. Ancient Okunev genetic lineages mainly flowed into later groups like the Xiongnu, Tiele, modern Chinese, and central-west Eurasians. Our high-resolution update of the Q-L330-YP761 phylogeography reveals multiple east-west differentiation patterns, strongly supporting the substantial westward migration of some genetic lineages from the eastern Eurasian steppe. However, ancient DNA data from northern China and northern South Asia are still scarce, and no Okunev-related genetic lineages have been identified in these regions so far. It remains unclear whether the influence of Okunev culture on these two regions originated from the spread of cultural factors or population movements. It remains unclear whether the influence of Okunev culture on these two regions is due to the spread of cultural factors or population mobility.
This research presents the largest dataset of ancient and modern sequences of the patrilineal haplogroup Q-L330-YP761 to date, the most revised phylogeny, and related evolutionary background information. The analysis in this paper shows that the patrilineal haplogroup Q-L330-YP761 originated in South Siberia, supporting the idea that the ancient Okunev people were direct descendants of the Upper Paleolithic Ancient Northern Eurasians (ANE). In subsequent stages, populations with diverse sub-branches of haplogroup Q-L330-YP761 as their main patrilineage were mainly active on the Mongolian Plateau, and they were likely part of the ancient Xiongnu, Turkic, and Uyghur peoples. From here, some of them spread to North China and eventually integrated into the Chinese population. Some other parts spread westward to Central Asia and Europe and participated in the formation of local populations in subsequent historical periods. The appearance of Okunev-like cultural elements in archaeological remains in the Mongolian Plateau, North China, India, Central Asia, and Eastern Europe may be related to the spread of the patrilineal haplogroup Q-L330-YP761. For the first time, we have provided a complete description of the origin and evolution of patrilineal haplogroup Q-L330-YP761 and its related ancient and modern populations over several thousand years.
In summary, this study comprehensively analyzed ancient and modern DNA data related to the paternal formation process of ancient Okunev people and their descendants. We provided a clear scenario for the rise of the Okunev people and the demographic history of the descendants’ lineages of the Okunev people in the late history period. This study’s data and analysis results will assist in further future research on the history of nomadic populations across the Eurasian Steppe.
The following supporting information can be found at: https://www.sciepublish.com/article/pii/504, Supplementary Table S1. The detailed patrilineal phylogenetic tree of haplogroup Q-L330-YP761.
We thank all donors for providing DNA samples and/or DNA sequences.
L.L.: Conceptualization, Methodology, Software, Formal analysis, Visualization, Writing—Original draft, Writing—Review & Editing. S.X.: Writing—Original draft, Validation, Investigation. H.L.: Data curation. L.W.: Data curation, Writing—Review & Editing, Supervision. H.Y.: Data curation, Writing—Review & Editing, Supervision, Project administration, Funding acquisition.
The study and sample collection procedures were reviewed and approved by the Medical Ethics Committee of Fudan University Approval NO.15137 (approved on 2 March 2015), adhering to the ethical principles of the 2013 Declaration of Helsinki by the World Medical Association.
Informed consent was obtained from all participants in the study.
Following the regulations of the Human Genetic Resources Administration of China (HGRAC), the raw data reported in this paper are available upon reasonable request from the corresponding author.
This work was supported by the Annual Open Project of the Key Research Institute for Humanities and Social Sciences at Inner Universities of Inner Mongolia Autonomous Region (Grant No. KFSM-SDSK0305).
The authors report no conflicts of interest.
Liu L, Xie S, Li H, Wei L, Yu H. Demographic History of Ancient Okunev People and Their Kin across Eurasia: A Patrilineal Perspective. Nature Anthropology 2025, 3, 10006. https://doi.org/10.70322/natanthropol.2025.10006
Liu L, Xie S, Li H, Wei L, Yu H. Demographic History of Ancient Okunev People and Their Kin across Eurasia: A Patrilineal Perspective. Nature Anthropology. 2025; 3(2):10006. https://doi.org/10.70322/natanthropol.2025.10006