mtDNA analysis reveals a major late Paleolithic population expansion from southwestern to northeastern Europe.

mtDNA sequence variation was studied in 419 individuals from nine Eurasian populations, by high-resolution RFLP analysis, and it was followed by sequencing of the control region of a subset of these mtDNAs and a detailed survey of previously published data from numerous other European populations. This analysis revealed that a major Paleolithic population expansion from the "Atlantic zone" (southwestern Europe) occurred 10,000-15,000 years ago, after the Last Glacial Maximum. As an mtDNA marker for this expansion we identified haplogroup V, an autochthonous European haplogroup, which most likely originated in the northern Iberian peninsula or southwestern France at about the time of the Younger Dryas. Its sister haplogroup, H, which is distributed throughout the entire range of Caucasoid populations and which originated in the Near East approximately 25,000-30,000 years ago, also took part in this expansion, thus rendering it by far the most frequent (40%-60%) haplogroup in western Europe. Subsequent migrations after the Younger Dryas eventually carried those "Atlantic" mtDNAs into central and northern Europe. This scenario, already implied by archaeological records, is given overwhelming support from both the distribution of the autochthonous European Y chromosome type 15, as detected by the probes 49a/f, and the synthetic maps of nuclear data.     

Reproductive isolation and environmental adaptation shape the phylogeography of mountain pine beetle (Dendroctonus ponderosae)

Chromosomal rearrangement can be an important mechanism driving population differentiation and incipient speciation. In the mountain pine beetle (MPB, Dendroctonus ponderosae), deletions on the Y chromosome that are polymorphic among populations are associated with reproductive incompatibility. Here, we used RAD sequencing across the entire MPB range in western North America to reveal the extent of the phylogeographic differences between Y haplotypes compared to autosomal and X‐linked loci. Clustering and geneflow analyses revealed three distinct Y haplogroups geographically positioned within and on either side of the Great Basin Desert. Despite close geographic proximity between populations on the boundaries of each Y haplogroup, there was extremely low Y haplogroup mixing among populations, and gene flow on the autosomes was reduced across Y haplogroup boundaries. These results are consistent with a previous study suggesting that independent degradation of a recently evolved neo‐Y chromosome in previously isolated populations causes male sterility or inviability among Y haplotype lineages. Phylogeographic results supported historic contraction of MPB into three separate Pleistocene glacial refugia followed by postglacial range expansion and secondary contact. Distinct sets of SNPs were statistically associated with environmental data among the most genetically distinct sets of geographic populations. This finding suggests that the process of adaptation to local climatic conditions is influenced by population genetic structure, with evidence for largely independent evolution in the most genetically isolated Y haplogroup.     

Y-chromosome analysis in Egypt suggests a genetic regional continuity in Northeastern Africa

The geographic location of Egypt, at the interface between North Africa, the Middle East, and southern Europe, prompted us to investigate the genetic diversity of this population and its relationship with neighboring populations. To assess the extent to which the modern Egyptian population reflects this intermediate geographic position, ten Unique Event Polymorphisms (UEPs), mapping to the nonrecombining portion of the Y chromosome, have been typed in 164 Y chromosomes from three North African populations. The analysis of these binary markers, which define 11 Y-chromosome lineages, were used to determine the haplogroup frequencies in Egyptians, Moroccan Arabs, and Moroccan Berbers and thereby define the Y-chromosome background in these regions. Pairwise comparisons with a set of 15 different populations from neighboring European, North African, and Middle Eastern populations and geographic analysis showed the absence of any significant genetic barrier in the eastern part of the Mediterranean area, suggesting that genetic variation and gene flow in this area follow the "isolation-by-distance" model. These results are in sharp contrast with the observation of a strong north-south genetic barrier in the western Mediterranean basin, defined by the Gibraltar Strait. Thus, the Y-chromosome gene pool in the modern Egyptian population reflects a mixture of European, Middle Eastern, and African characteristics, highlighting the importance of ancient and recent migration waves, followed by gene flow, in the region.     

Differentiation of mitochondrial DNA and Y chromosomes in Russian populations

The genetic composition of the Russian population was investigated by analyzing both mitochondrial DNA (mtDNA) and Y-chromosome loci polymorphisms that allow for the different components of a population gene pool to be studied, depending on the mode of DNA marker inheritance. mtDNA sequence variation was examined by using hypervariable segment I (HVSI) sequencing and restriction analysis of the haplogroup-specific sites in 325 individuals representing 5 Russian populations from the European part of Russia. The Y-chromosome variation was investigated in 338 individuals from 8 Russian populations (including 5 populations analyzed for mtDNA variation) using 12 binary markers. For both uniparental systems most of the observed haplogroups fell into major West Eurasian haplogroups (97.9% and 99.7% for mtDNA and Y-chromosome haplogroups, respectively). Multidimensional scaling analysis based on pairwise F(ST) values between mtDNA HVSI sequences in Russians compared to other European populations revealed a considerable heterogeneity of Russian populations; populations from the southern and western parts of Russia are separated from eastern and northern populations. Meanwhile, the multidimensional scaling analysis based on Y-chromosome haplogroup F(ST) values demonstrates that the Russian gene pool is close to central-eastern European populations, with a much higher similarity to the Baltic and Finno-Ugric male pools from northern European Russia. This discrepancy in the depth of penetration of mtDNA and Y-chromosome lineages characteristic for the most southwestern Russian populations into the east and north of eastern Europe appears to indicate that Russian colonization of the northeastern territories might have been accomplished mainly by males rather than by females.     

From East to West: patterns of genetic diversity of populations living in four Eurasian regions

We have analyzed the distribution and patterns of the genetic diversity of eight Alu loci (ACE, ApoA1, PV92, TPA25, NBC27, NBC102, NBC148, and NBC182) in 1,049 individuals representing 16 populations of the Volga-Ural region (Bashkirs, Tatars, Komis, Maris, Mordvins, and Udmurts), Central Asia (Kazakhs, Uzbeks, and Uighurs), the North Caucasus (Karachays, Kumyks, Kuban Nogays, and Karanogays), and Central South Siberia (Yakuts, Kalmyks and Evenks). Geographic divide between Europe and Asia, e.g. the Ural Mountains and the Caspian Sea, can also be considered as a genetic boundary. The data indicates that the populations of the two boundary regions between Europe and Asia, the Volga-Ural region of Russia, and populations of the North Caucasus are more similar to European than to Asian populations. Finally, Siberian and Central Asian populations are genetically closely related to each other.     

High-resolution phylogenetic analysis of southeastern Europe traces major episodes of paternal gene flow among Slavic populations

The extent and nature of southeastern Europe (SEE) paternal genetic contribution to the European genetic landscape were explored based on a high-resolution Y chromosome analysis involving 681 males from seven populations in the region. Paternal lineages present in SEE were compared with previously published data from 81 western Eurasian populations and 5,017 Y chromosome samples. The finding that five major haplogroups (E3b1, I1b* (xM26), J2, R1a, and R1b) comprise more than 70% of SEE total genetic variation is consistent with the typical European Y chromosome gene pool. However, distribution of major Y chromosomal lineages and estimated expansion signals clarify the specific role of this region in structuring of European, and particularly Slavic, paternal genetic heritage. Contemporary Slavic paternal gene pool, mostly characterized by the predominance of R1a and I1b* (xM26) and scarcity of E3b1 lineages, is a result of two major prehistoric gene flows with opposite directions: the post-Last Glacial Maximum R1a expansion from east to west, the Younger Dryas-Holocene I1b* (xM26) diffusion out of SEE in addition to subsequent R1a and I1b* (xM26) putative gene flows between eastern Europe and SEE, and a rather weak extent of E3b1 diffusion toward regions nowadays occupied by Slavic-speaking populations.     

Similarities and distinctions in Y chromosome gene pool of Western Slavs

Analysis of Y chromosome Y‐STRs has proven to be a useful tool in the field of population genetics, especially in the case of closely related populations. We collected DNA samples from 169 males of Czech origin, 80 males of Slovakian origin, and 142 males dwelling Northern Poland. We performed Y‐STR analysis of 12 loci in the samples collected (PowerPlex Y system from Promega) and compared the Y chromosome haplotype frequencies between the populations investigated. Also, we used Y‐STR data available from the literature for comparison purposes. We observed significant differences between Y chromosome pools of Czechs and Slovaks compared to other Slavic and European populations. At the same time we were able to point to a specific group of Y‐STR haplotypes belonging to an R1a haplogroup that seems to be shared by Slavic populations dwelling in Central Europe. The observed Y chromosome diversity may be explained by taking into consideration archeological and historical data regarding early Slav migrations. Am J Phys Anthropol 142:540–548, 2010. © 2010 Wiley‐Liss, Inc.     

Polymorphism of hypervariable segment I of mitochondrial DNA in three ethnic groups of the Volga-Ural region

A population genetic study of the polymorphism in the first hypervariable segment (HVSI) of mitochondrial DNA control region was carried out for three ethnic populations of the Volga-Ural region, Bashkirs, Russians, and Komi-Permyaks. This analysis showed that most of the mtDNA HVSI haplotypes detected in the populations of Bashkirs, Russians and Komi-Permyaks contained the combinations of nucleotide substitutions detected earlier in Asian, European, and Finno-Ugric populations. These findings are consistent with historical, anthropological, and ethnographical data suggesting the presence of European and Mongoloid components of different geographical descent in the gene pool of the contemporary population of the Volga-Ural region. The data on the genetic structure and the phylogenetic relationships between populations of the Volga-Ural region based on modern molecular genetic methods of mitochondrial genome investigation would be a substantial addition to the already existing information for some other regions of Europe and Asia. These data would provide more complete examination of the development of interethnic diversity of mitochondrial gene pools of contemporary ethnic populations with the purpose of reconstructing the genetic demographic processes that accompanied peopling of the Middle Ural and Volga region.     

Human-Mediated Admixture and Selection Shape the Diversity on the Modern Swine (Sus scrofa) Y Chromosomes

Throughout its distribution across Eurasia, domestic pig (Sus scrofa) populations have acquired differences through natural and artificial selection, and have often interbred. We resequenced 80 Eurasian pigs from nine different Asian and European breeds; we identify 42,288 reliable SNPs on the Y chromosome in a panel of 103 males, among which 96.1% are newly detected. Based on these new data, we elucidate the evolutionary history of pigs through the lens of the Y chromosome. We identify two highly divergent haplogroups: one present only in Asia and one fixed in Europe but present in some Asian populations. Analyzing the European haplotypes present in Asian populations, we find evidence of three independent waves of introgression from Europe to Asia in last 200 years, agreeing well with the literature and historical records. The diverse European lineages were brought in China by humans and left significant imprints not only on the autosomes but also on the Y chromosome of geographically and genetically distinct Chinese pig breeds. We also find a general excess of European ancestry on Y chromosomes relative to autosomes in Chinese pigs, an observation that cannot be explained solely by sex-biased migration and genetic drift. The European Y haplotype is associated with leaner meat production, and we hypothesize that the European Y chromosome increased in frequency in Chinese populations due to artificial selection. We find evidence of Y chromosomal gene flow between Sumatran wild boar and Chinese pigs. Our results demonstrate how human-mediated admixture and selection shaped the distribution of modern swine Y chromosomes.     

Tracing past human male movements in northern/eastern Africa and western Eurasia: new clues from Y-chromosomal haplogroups E-M78 and J-M12

Detailed population data were obtained on the distribution of novel biallelic markers that finely dissect the human Y-chromosome haplogroup E-M78. Among 6,501 Y chromosomes sampled in 81 human populations worldwide, we found 517 E-M78 chromosomes and assigned them to 10 subhaplogroups. Eleven microsatellite loci were used to further evaluate subhaplogroup internal diversification. The geographic and quantitative analyses of haplogroup and microsatellite diversity is strongly suggestive of a northeastern African origin of E-M78, with a corridor for bidirectional migrations between northeastern and eastern Africa (at least 2 episodes between 23.9-17.3 ky and 18.0-5.9 ky ago), trans-Mediterranean migrations directly from northern Africa to Europe (mainly in the last 13.0 ky), and flow from northeastern Africa to western Asia between 20.0 and 6.8 ky ago. A single clade within E-M78 (E-V13) highlights a range expansion in the Bronze Age of southeastern Europe, which is also detected by haplogroup J-M12. Phylogeography pattern of molecular radiation and coalescence estimates for both haplogroups are similar and reveal that the genetic landscape of this region is, to a large extent, the consequence of a recent population growth in situ rather than the result of a mere flow of western Asian migrants in the early Neolithic. Our results not only provide a refinement of previous evolutionary hypotheses but also well-defined time frames for past human movements both in northern/eastern Africa and western Eurasia.     

Mitochondrial DNA variation in Russian populations of Stavropol krai,Orel and Saratov oblasts

Mitochondrial DNA (mtDNA) polymorphism was examined in three Russian populations from the European part of Russia (Stavropol krai, Orel oblast, and Saratov oblast). This analysis showed that mitochondrial gene pool of Russians was represented by the mtDNA types belonging to haplogroups H, V, HV*, J, T, U, K, I, W, and X. A mongoloid admixture (1.5%) was revealed in the form of mtDNA types of macrohaplogroup M. Comparative analysis of the mtDNA haplogroup frequency distribution patterns in six Russian populations from the European part of Russia indicated the absence of substantial genetic differences between them. However, in Russian populations from the southern and central regions the frequency of haplogroup V (average frequency 8%) was higher than in the populations from more northern regions. Based on the data on mtDNA HVS1 sequence variation, it was shown that the diversity of haplogroup V in Russians (h = 0.72) corresponded to the highest h values observed in Europe. The reasons for genetic differentiation of the Russian population (historical, ecological, and adaptive) are discussed.     

Mitochondrial DNA Variation in Russian Populations of Stavropol Krai, Orel and Saratov Oblasts

Mitochondrial DNA (mtDNA) polymorphism was examined in three Russian populations from the European part of Russia (Stavropol krai, Orel oblast, and Saratov oblast). This analysis showed that mitochondrial gene pool of Russians was represented by the mtDNA types belonging to haplogroups H, V, HV*, J, T, U, K, I, W, and X. A mongoloid admixture (1.5%) was revealed in the form of mtDNA types of macrohaplogroup M. Comparative analysis of the mtDNA haplogroup frequency distribution patterns in six Russian populations from the European part of Russia indicated the absence of substantial genetic differences between them. However, in Russian populations from the southern and central regions the frequency of haplogroup V (average frequency 8%) was higher than in the populations from more northern regions. Based on the data on mtDNA HVS1 sequence variation, it was shown that the diversity of haplogroup V in Russians (h= 0.742) corresponded to the highest h values observed in Europe. The reasons for genetic differentiation of the Russian population (historical, ecological, and adaptive) are discussed.     

A Genome-Wide Analysis of Populations from European Russia Reveals a New Pole of Genetic Diversity in Northern Europe

Several studies examined the fine-scale structure of human genetic variation in Europe. However, the European sets analyzed represent mainly northern, western, central, and southern Europe. Here, we report an analysis of approximately 166,000 single nucleotide polymorphisms in populations from eastern (northeastern) Europe: four Russian populations from European Russia, and three populations from the northernmost Finno-Ugric ethnicities (Veps and two contrast groups of Komi people). These were compared with several reference European samples, including Finns, Estonians, Latvians, Poles, Czechs, Germans, and Italians. The results obtained demonstrated genetic heterogeneity of populations living in the region studied. Russians from the central part of European Russia (Tver, Murom, and Kursk) exhibited similarities with populations from central–eastern Europe, and were distant from Russian sample from the northern Russia (Mezen district, Archangelsk region). Komi samples, especially Izhemski Komi, were significantly different from all other populations studied. These can be considered as a second pole of genetic diversity in northern Europe (in addition to the pole, occupied by Finns), as they had a distinct ancestry component. Russians from Mezen and the Finnic-speaking Veps were positioned between the two poles, but differed from each other in the proportions of Komi and Finnic ancestries. In general, our data provides a more complete genetic map of Europe accounting for the diversity in its most eastern (northeastern) populations.     

High levels of Y-chromosome differentiation among native Siberian populations and the genetic signature of a boreal hunter-gatherer way of life

We examined genetic variation on the nonrecombining portion of the Y chromosome (NRY) to investigate the paternal population structure of indigenous Siberian groups and to reconstruct the historical events leading to the peopling of Siberia. A set of 62 biallelic markers on the NRY were genotyped in 1432 males representing 18 Siberian populations, as well as nine populations from Central and East Asia and one from European Russia. A subset of these markers defines the 18 major NRY haplogroups (A-R) recently described by the Y Chromosome Consortium (YCC 2002). While only four of these 18 major NRY haplogroups accounted for -95% of Siberian Y-chromosome variation, native Siberian populations differed greatly in their haplogroup composition and exhibited the highest phiST value for any region of the world. When we divided our Siberian sample into four geographic regions versus five major linguistic groupings, analyses of molecular variance (AMOVA) indicated higher phiST and phiCT values for linguistic groups than for geographic groups. Mantel tests also supported the existence of NRY genetic patterns that were correlated with language, indicating that language affiliation might be a better predictor of the genetic affinity among Siberians than their present geographic position. The combined results, including those from a nested cladistic analysis, underscored the important role of directed dispersals, range expansions, and long-distance colonizations bound by common ethnic and linguistic affiliation in shaping the genetic landscape of Siberia. The Siberian pattern of reduced haplogroup diversity within populations combined with high levels of differentiation among populations may be a general feature characteristic of indigenous groups that have small effective population sizes and that have been isolated for long periods of time.     

Genetic diversity of two haploid markers in the Udegey population from southeastern Siberia

The Udegeys are a small ethnic group who live along the tributaries of the Amur River Basin of southeastern Siberia in Russia. They are thought to speak a language belonging to a subdivision of the Tungusic‐Manchu branch of the Altaic family. To understand the genetic features and genetic history of the Udegeys, we analyzed two haploid markers, mitochondrial DNA (mtDNA), and Y‐chromosomal variation, in 51 individuals (including 21 males) from the Udegey population. In general, the Udegeys' mtDNA profiles revealed similarities to Siberians and other northeastern Asian populations, although a moderate European contribution was also detected. Interestingly, pairwise values of F_ST and the MDS plots based on the mtDNA variation showed that the Orok and Nivkh inhabiting the very same region of the Udegey were significantly different from the Udegey, implying that they may have been isolated and undergone substantial genetic drift. The Udegeys were characterized by a high frequency (66.7%) of Y chromosome haplogroup C, indicating a close genetic relationship with Mongolians and Siberians. On the paternal side, however, very little admixture was observed between the Udegeys and Europeans. Thus, the combined haploid genetic markers of both mtDNA and the Y chromosome imply that the Udegeys are overall closest to Siberians and northeast Asians of the Altaic linguistic family, with a minor maternal contribution from the European part of the continent. Am J Phys Anthropol, 2010. © 2009 Wiley‐Liss, Inc.     

Genome-wide association study of bronchial asthma in the Volga-Ural region of Russia

Bronchial asthma is a chronic inflammatory respiratory disease that is caused by the complex interaction of environmental influences and genetic susceptibility. The first genome-wide association study of bronchial asthma discovered a significant association between SNPs within 17q12-21 genomic region and childhood bronchial asthma in individuals of European descent. Association with this genomic region was then replicated in a number of independent samples of European and Asian descent. Here we report results of the first genome-wide association study of bronchial asthma in the Volga-Ural region of Russia. The present study includes 358 unrelated patients with physician-diagnosed bronchial asthma and 369 disease-free control subjects of different ethnic origin (Russians, Tatars and Bashkirs). Genotyping of DNA samples was carried out using the Illumina Human610 quad array as a part of GABRIEL project (contract from the EC No LSHB-CT-2006-018996). After QC filtering procedures, a final set of 550915 SNPs genotyped in 330 cases and 348 controls was tested for association with bronchial asthma. Five markers on chromosome 17q12-21 showed statistically significant association with bronchial asthma (p < or = 4.79 x 10(-7)). SNP rs7216389 with the strongest evidence for association (p = 1.01 x 10(-7)) is located within the first intron of the GSDMB gene. Evidence for association was stronger with childhood-onset asthma (p = 1.97 x 10(-6) for SNP rs7216389) compared to late-onset asthma (p = 1.8 x 10(-4) for SNP rs7216389). Our replication study using three SNPs within GSDMB gene confirmed association with only childhood-onset asthma. In summary, these results suggest an important role for genetic variants within 17q12-q21 region in the development of bronchial asthma in the Volga-Ural region of Russia.     

Distribution of Y chromosomes among native North Americans: a study of Athapaskan population history

In this study, 231 Y chromosomes from 12 populations were typed for four diagnostic single nucleotide polymorphisms (SNPs) to determine haplogroup membership and 43 Y chromosomes from three of these populations were typed for eight short tandem repeats (STRs) to determine haplotypes. These data were combined with previously published data, amounting to 724 Y chromosomes from 26 populations in North America, and analyzed to investigate the geographic distribution of Y chromosomes among native North Americans and to test the Southern Athapaskan migration hypothesis. The results suggest that European admixture has significantly altered the distribution of Y chromosomes in North America and because of this caution should be taken when inferring prehistoric population events in North America using Y chromosome data alone. However, consistent with studies of other genetic systems, we are still able to identify close relationships among Y chromosomes in Athapaskans from the Subarctic and the Southwest, suggesting that a small number of proto-Apachean migrants from the Subarctic founded the Southwest Athapaskan populations.     

Diversity of mitochondrial DNA haplotypes in ethnic populations of the Volga-Ural region of Russia

The mtDNA polymorphism was analyzed in eight ethnic groups (N = 979) of the Volga-Ural region. Most mtDNA variants belonged to haplogroups H, U, T, J, W, I, R, and N1 characteristic of West Eurasian populations. The most frequent were haplogroups H (12-42%) and U (18-44%). East Eurasian mtDNA types (A, B, Y, F, M, N9) were also observed. Genetic diversity was higher in Turkic than in Finno-Ugric populations. The frequency of mtDNA types characteristic of Siberian and Central Asian populations substantially increased in the ethnic groups living closer to the Urals, a boundary between Europe and Asia. Geographic distances, rather than linguistic barriers, were assumed to play the major role in distribution of mtDNA types in the Volga-Ural region. Thus, as concerns the maternal lineage, the Finno-Ugric populations of the region proved to be more similar to their Turkic neighbors rather than to linguistically related Balto-Finnish ethnic groups.     

Phylogeography and Pleistocene refugia of the Little Owl Athene noctua inferred from mtDNA sequence data

Pleistocene glaciations greatly affected the distribution of genetic diversity in animal populations. The Little Owl is widely distributed in temperate regions and could have survived the last glaciations in southern refugia. To describe the phylogeographical structure of European populations, we sequenced the mitochondrial cytochrome c oxidase I (COI) and control region (CR1) in 326 individuals sampled from 22 locations. Phylogenetic analyses of COI identified two deeply divergent clades: a western haplogroup distributed in western and northwestern Europe, and an eastern haplogroup distributed in southeastern Europe. Faster evolving CR1 sequences supported the divergence between these two main clades, and identified three subgroups within the eastern clade: Balkan, southern Italian and Sardinian. Divergence times estimated from COI with fossil calibrations indicate that the western and eastern haplogroups split 2.01–1.71 Mya. Slightly different times for splits were found using the standard 2% rate and 7.3% mtDNA neutral substitution rate. CR1 sequences dated the origin of endemic Sardinian haplotypes at 1.04–0.26 Mya and the split between southern Italian and Balkan haplogroups at 0.72–0.21 Mya, coincident with the onset of two Pleistocene glaciations. Admixture of mtDNA haplotypes was detected in northern Italy and in central Europe. These findings support a model of southern Mediterranean and Balkan refugia, with postglacial expansion and secondary contacts for Little Owl populations. Central and northern Europe was predominantly recolonized by Little Owls from Iberia, whereas expansion out of the Balkans was more limited. Northward expansion of the Italian haplogroup was probably prevented by the Alps, and the Sardinian haplotypes remained confined to the island. Results showed a clear genetic pattern differentiating putative subspecies. Genetic distances between haplogroups were comparable with those recorded between different avian species.     

Genetic Evidence of an East Asian Origin and Paleolithic Northward Migration of Y-chromosome Haplogroup N

The Y-chromosome haplogroup N-M231 (Hg N) is distributed widely in eastern and central Asia, Siberia, as well as in eastern and northern Europe. Previous studies suggested a counterclockwise prehistoric migration of Hg N from eastern Asia to eastern and northern Europe. However, the root of this Y chromosome lineage and its detailed dispersal pattern across eastern Asia are still unclear. We analyzed haplogroup profiles and phylogeographic patterns of 1,570 Hg N individuals from 20,826 males in 359 populations across Eurasia. We first genotyped 6,371 males from 169 populations in China and Cambodia, and generated data of 360 Hg N individuals, and then combined published data on 1,210 Hg N individuals from Japanese, Southeast Asian, Siberian, European and Central Asian populations. The results showed that the sub-haplogroups of Hg N have a distinct geographical distribution. The highest Y-STR diversity of the ancestral Hg N sub-haplogroups was observed in the southern part of mainland East Asia, and further phylogeographic analyses supports an origin of Hg N in southern China. Combined with previous data, we propose that the early northward dispersal of Hg N started from southern China about 21 thousand years ago (kya), expanding into northern China 12–18 kya, and reaching further north to Siberia about 12–14 kya before a population expansion and westward migration into Central Asia and eastern/northern Europe around 8.0–10.0 kya. This northward migration of Hg N likewise coincides with retreating ice sheets after the Last Glacial Maximum (22–18 kya) in mainland East Asia.