Y-chromosome haplogroup N dispersals from south Siberia to Europe (2024)

Molecular phylogeography of Y-chromosome haplogroups provides intriguing clues into the history of human populations, especially in reconstructing of ancient migrations and more recent gene flow episodes (Underhill 2003). One of the most widespread Y-haplogroups in north Eurasia is haplogroup N defined by markers M231 and LLY22g (Jobling and Tyler-Smith 1995; Cinnioglu et al. 2004). This haplogroup is represented by subclades N^*, N1, N2 and N3, which are thought to be historically informative in north Eurasian populations (Rootsi et al. 2007). It has been generally accepted that the hg N ancestor is of east Asian ancestry, because haplotypes belonging to ancient paragroup N^* are still widely distributed in southeastern Asia and, rarely, even in south Siberia (Rootsi et al. 2007; Sengupta et al. 2006). However, the origin and migration ways of hgs N2 and N3 are not fully understood. According to the first studies, the most likely place of origin of these Y-chromosomes is the area of Mongolia and north China (Zerjal et al. 1997; Karafet et al. 2002). However, substantial frequency of hgs N2 and N3 in Eastern Europe, mainly in Finno-Ugric and Turkic-speaking populations, indicates that these haplogroups may have originated in Eastern Europe (Rosser et al. 2000; Underhill et al. 2001). Moreover, STR diversity within hg N3 has been detected on higher level in Eastern Europeans than in Siberians, which can be interpreted as supporting an origin of hg N3a in Eastern Europe (Villems et al. 2002). Recently, phylogeographic study of hg N structure in north Eurasian populations has demonstrated that south Siberia could be a place of transition of hg N (including N3a and N2) westward to Eastern Europe (Rootsi et al. 2007). However, there is no evidence still that the age of accumulated STR variation within hg N3a is higher in south Siberia.

In this study, we improve our understanding of hg N phylogeographic structure by examining STR variation in a large number of individuals of south Siberian and Eastern European origin. The results obtained give a clear evidence of that hg N3a is older in south Siberia than in Eastern Europe. In addition, the substructure of hg N2-A revealed in south Siberia suggests that this region was a place of several expansions of hg N haplotypes at different times during the Holocene.

Subjects and DNA typing

A total of 1,438 samples (whole blood and hair root samples) from unrelated males were collected in populations of south Siberia (south Altaians, Shors, Tuvinians, Tofalars, Sojots, Khakassians, Buryats, Kalmyks, Evenks, Yakuts, Evens, Koryaks), central and southwest Asia (Mongolians, Koreans, Tajiks, Iranians) and Eastern Europe (Russians) (Table1). All samples were collected with appropriate ethical approval and informed consent.

Hg N markers LLY22g (for the whole hg N) (Jobling and Tyler-Smith 1995), Tat (for N3) (Zerjal et al. 1997), M178 (for N3a) (Underhill et al. 2000; Kharkov et al. 2005), P43 (for N2) (Karafet et al. 2002) were assayed as described in the referred papers. The Y-SNP haplogroup nomenclature used here is according to the recommendations of the Y Chromosome Consortium (2002). A total of 234 samples belonging to haplogroups N2 and N3 were analyzed at 11 STR loci (DYS19, DYS385a/b, DYS389I, DYS389II, DYS390, DYS391, DYS392, DYS393, DYS437, DYS438, DYS439) using PowerPlex^® Y System (Promega Corporation, Madison, WI). Products of amplification were analyzed on ABI 3100 Genetic Analyzer (Applied Biosystems) with an appropriate filter set prepared on the basis of the relevant PowerPlex^® Matrix Standards provided by Promega Corporation. Electrophoresis results were analyzed using Genscan v. 3.7 and Genotyper v. 3.7 software (Applied Biosystems). Alleles were designated by repeat numbers.

Data analysis

Median-joining (MJ) networks of hg N STR-haplotypes were constructed using the Network 4.1.1.2 program (http://www.fluxus-engineering.com). Networks were calculated by the median-joining method after processing the data with the reduced median network method in order to reduce reticulations within the network. For network construction, STR loci were weighted according to the average of their variability in the corresponding haplogroups. All of these approaches are described elsewhere (Bandelt et al. 2000).

The age of microsatellite variation within Y-chromosome haplogroups was estimated as the average squared difference in the number of repeats between all current chromosomes and the founder haplotype, averaged over microsatellite loci and divided by w=6.9×10⁻⁴ per 25years (Zhivotovsky et al. 2004). We followed the recommendation of Sengupta et al. (2006) and used the median haplotypes (formed by the median values of the repeat scores at each microsatellite locus within each haplogroup) as the founding ones.

Published hg N STR-data for Eastern Europeans (Komis, Tatars, Maris, Udmurts, Bashkirs, Karelians, Vepsas, Estonians, Ukrainians and Slovaks) and Asians (Altaians, Tuvinians, Khakassians, Evenks, Eskimo and Chinese) (Rootsi et al. 2007) were also included in the analyses.

SNP analysis of 1,438 Y-chromosomes from 17 ethnic groups representing populations of north Asia and Eastern Europe demonstrates that hgs N3a and N2 are relatively frequent in some populations (Table1). Haplogroup N3a, present at its highest frequency in Yakuts (80%), Tofalars (27%), Koryaks (25%) and Buryats (18%), is relatively frequent in all other Siberian populations as well as in the Russians (11%) studied. Haplogroup N2 chromosomes are present in the majority of the Siberian samples, reaching the highest frequencies in Tofalars (43%), Khakassians (34%) and Evenks (24%). However, this haplogroup is absent in northeastern Siberian populations (such as Evens, Koryaks and Yakuts). In addition, hg N2 is very rare in Russians (only 0.2%).

In order to obtain a better understanding of phylogenetic relationships between Y-chromosomes found in individuals carrying N2- and N3-lineages, we have determined their haplotypes at 11 microsatellite loci (Table2). In general, 234 samples were genotyped. The age of accumulated STR variation within hg N3a estimated using the method of Zhivotovsky et al. (2004) indicates that this haplogroup is much more diverse in Eastern Europe than in Siberia (giving ages of 9.0ky and 7.9ky ago, respectively) (Table3), which is in agreement with earlier observations made by Villems et al. (2002). However, median network analysis has shown that, in fact, there are two subclusters of N3a, both present in Siberia and Eastern Europe—N3a1 and N3a2 (Fig.1). This subdivision of hg N3a has been suggested recently (Rootsi et al. 2007); however, phylogenetic and statistical analysis of hg N3a subclusters had not been performed.

The presumed founder haplotypes for both subclusters (they are shown in bold in Table2) are one-step neighbors to each other because there is only one difference in the numbers of repeats at the DYS391 locus between the two haplotypes. Nevertheless, this subdivision appears to be phylogenetically stable, because only four haplotypes (two in Russians, one in Yakuts and one in Khakassians) out of 40 N3a1-haplotypes demonstrated back-mutations from DYS392-11 to DYS392-10 allele (Table2). Meanwhile, no back mutations from DYS392-10 to DYS392-11 allele were revealed in subcluster N3a2. Our analysis has shown that subcluster N3a1 is much older than N3a2, giving the age of 9.1ky and 5.0ky, respectively. Importantly, the subcluster N3a1 appears to be older in south Siberia than in Europe (10ky and 8.2ky ago, respectively), allowing the suggestion that the first expansion of N3a1 occurred in the south Siberian region, and then this subcluster spread to the West, into the Urals and Eastern Europe. The ancestral haplotype for N3a1 is present in Eastern Europeans (Russians and Estonians) and in south Siberians from the east Sayan regions (such as Tuvinians, Tofalars and Sojots), giving also rise to regional-specific branches of N3a1 in Yakuts and Khakassians/Tuvinians (Fig.1). In general, this subcluster encompasses 58% of N3a individuals studied here (Table1). It is noteworthy that the overwhelming majority of Russians (83%) as well as other Slavs and Baltic-Finnic individuals taken from the study of Rootsi et al. (2007) belong to subcluster N3a1.

Subcluster N3a2 is characterized by different tree topology. It has an almost star-like branching pattern in south Siberia, where it is present mostly in the Baikal region among Buryats. But it shows a very complex picture of the phylogenetic tree branching into its Eastern European part represented predominantly by Finno-Ugric and Turkic-speaking populations of the Volga-Ural region. High STR variation of subcluster N3a2 observed in the Volga-Ural populations, in comparison with south Siberians (6.6ky vs. 3.7ky ago), seems to be a consequence of serial bottlenecks and/or the small sample sizes of Eastern European N3a2 chromosomes studied. In addition, the Eastern European part of the N3a2 network displayed clear non-starlike features, which greatly reduce the precision of genetic age estimates (Sailard et al. 2000). To elucidate this question, much more detailed information on STR variation should be obtained from Finno-Ugric and Turkic-speaking populations of the Volga-Ural region. Meanwhile, traces of recent expansions of subcluster N3a2 are obvious in different populations of south Siberia—mainly, in Buryats and south Altaians.

It has been recently demonstrated that N2 haplotypes form two separate subgroups, N2-A and N2-E, characterized by an Asian and European pattern of distribution, respectively (Rootsi et al. 2007). Analyzing STR variation of N2 haplotypes in a large sample of Siberians, we have found that all of them belong to cluster N2-A. According to our data, the age of cluster N2-A is 4.8±1.6ky, which is somewhat less than estimates made by Rootsi et al. (2007). Median network analysis of N2-A haplotypes showed an existence of two subclusters determined by different STR composition in their founder haplotypes (with differences at loci DYS19, DYS391 and DYS439). Subcluster N2-A1, with founder haplotype present in different populations (such as Khakassians, Shors, Mongols, Altaians, Evenks and Tofalars), appears to be an ancestral subcluster, later giving rise to N2-A2 and N2-E. This follows from its central position in the median joining tree constructed with taking into account the location of N^* lineages (Fig.2). The age of STR variation within hg N2-A1 is a little higher than that of N2-A2 (2.4 and 1.8kys, respectively). Interestingly, N2-A2 expansion is restricted mainly to the Tuva region.

Overall, the data obtained indicate that south Siberia was a place of several expansions of hgs N2 and N3 haplotypes at different times during the Holocene. Future examinations of Y-chromosome variation using a combined STR-SNP approach in additional population samples may enable a better definition of the differences between north Eurasian population groups.

Authors and Affiliations

1. Genetics Laboratory, Institute of Biological Problems of the North, Russian Academy of Sciences, Portovaya str. 18, 685000, Magadan, Russia

   Miroslava Derenko,Boris Malyarchuk&Galina Denisova

2. Forensic Medicine Institute, the Ludwik Rydygier Medical College, The Nicolaus Copernicus University in Torun, 85-094, Bydgoszcz, Poland

   Marcin Wozniak&Tomasz Grzybowski

3. Institute of General and Experimental Biology, Russian Academy of Sciences, 670047, Ulan-Ude, Russia

   Irina Dambueva

4. Vavilov Institute of General Genetics, Russian Academy of Sciences, Gubkin Str. 3, 119991, Moscow, Russia

   Ilia Zakharov

Corresponding author

Correspondence to Miroslava Derenko.

Reprints and permissions