Mitogenomics of the Koryaks and Evens of the northern coast of the Sea of Okhotsk

Due to the geographical proximity of the northern coast of the Sea of Okhotsk and Kamchatka Peninsula to the Beringia, the indigenous populations of these territories are of great interest for elucidating the human settlement history of northern Asia and America. Meanwhile, there is a clear shortage of genetic studies of the indigenous populations of the northern coast of the Sea of Okhotsk. Here, in order to examine their fine-scale matrilineal genetic structure, ancestry and relationships with neighboring populations, we analyzed 203 complete mitogenomes (174 of which are new) from population samples of the Koryaks and Evens of the northern coast of the Sea of Okhotsk and the Chukchi of the extreme northeast Asia. The patterns observed underscore the reduced level of genetic diversity found in the Koryak, Even, and Chukchi populations, which, along with the high degree of interpopulation differentiation, may be the result of genetic drift. Our phylogeographic analysis reveals common Paleo-Asiatic ancestry for 51.1% of the Koryaks and 17.8% of the Evens. About third of the mitogenomes found in the Koryaks and Evens might be considered as ethno-specific, as these are virtually absent elsewhere in North, Central and East Asia. Coalescence ages of most of these lineages coincide well with the emergence and development of the Tokarev and Old Koryak archaeological cultures associated with the formation of the Koryaks, as well as with the period of separation and split of the North Tungusic groups migrated northwards from the Lake Baikal or the Amur River area.

Response to Wyckelsma et al.: Loss of α-actinin-3 during human evolution provides superior cold resilience and muscle heat generation

The common loss-of-function mutation R577X in the structural muscle protein ACTN3 emerged as a potential target of positive selection from early studies and has been the focus of insightful physiological work suggesting a significant impact on muscle metabolism. Adaptation to cold climates has been proposed as a key adaptive mechanism explaining its global allele frequency patterns. Here, we re-examine this hypothesis analyzing modern (n = 3,626) and ancient (n = 1,651) genomic data by using allele-frequency as well as haplotype homozygosity-based methods. The presented results are more consistent with genetic drift rather than selection in cold climates as the main driver of the ACTN3 R577X frequency distribution in human populations across the world. This Matters Arising paper is in response to Wyckelsma et al. (2021),1 published in The American Journal of Human Genetics. See also the response by Wyckelsma et al. (2022),2 published in this issue.

Complete mitogenome data for the Serbian population: the contribution to high-quality forensic databases

Mitochondrial genome (mtDNA) is a valuable resource in resolving various human forensic casework. The usage of variability of complete mtDNA genomes increases their discriminatory power to the maximum and enables ultimate resolution of distinct maternal lineages. However, their wider employment in forensic casework is nowadays limited by the lack of appropriate reference database. In order to fill in the gap in the reference data, which, considering Slavic-speaking populations, currently comprises only mitogenomes of East and West Slavs, we present mitogenome data for 226 Serbians, representatives of South Slavs from the Balkan Peninsula. We found 143 (sub)haplogroups among which West Eurasian ones were dominant. The percentage of unique haplotypes was 85%, and the random match probability was as low as 0.53%. We support previous findings on both high levels of genetic diversity in the Serbian population and patterns of genetic differentiation among this and ten studied European populations. However, our high-resolution data supported more pronounced genetic differentiation among Serbians and two Slavic populations (Russians and Poles) as well as expansion of the Serbian population after the Last Glacial Maximum and during the Migration period (fourth to ninth century A.D.), as inferred from the Bayesian skyline analysis. Phylogenetic analysis of haplotypes found in Serbians contributed towards the improvement of the worldwide mtDNA phylogeny, which is essential for the interpretation of the mtDNA casework.

Whole mitochondrial genome diversity in two Hungarian populations

Complete mitochondrial genomics is an effective tool for studying the demographic history of human populations, but there is still a deficit of mitogenomic data in European populations. In this paper, we present results of study of variability of 80 complete mitochondrial genomes in two Hungarian populations from eastern part of Hungary (Szeged and Debrecen areas). The genetic diversity of Hungarian mitogenomes is remarkably high, reaching 99.9% in a combined sample. According to the analysis of molecular variance (AMOVA), European populations showed a low, but statistically significant level of between-population differentiation (Fst = 0.61%, p = 0), and two Hungarian populations demonstrate lack of between-population differences. Phylogeographic analysis allowed us to identify 71 different mtDNA sub-clades in Hungarians, sixteen of which are novel. Analysis of ancestry-informative mtDNA sub-clades revealed a complex genetic structure associated with the genetic impact of populations from different parts of Eurasia, though the contribution from European populations is the most pronounced. At least 8% of ancestry-informative haplotypes found in Hungarians demonstrate similarity with East and West Slavic populations (sub-clades H1c23a, H2a1c1, J2b1a6, T2b25a1, U4a2e, K1c1j, and I1a1c), while the influence of Siberian populations is not so noticeable (sub-clades A12a, C4a1a, and probably U4b1a4).

Mitogenomic diversity and differentiation of the Buryats

In this paper we present a results of first comprehensive study of the complete mitogenomes in the Buryats with regard to their belonging to the main regional (eastern and western Buryats); tribal (Khori, Ekhirid, Bulagad, and Khongodor), and ethno-territorial (Aginsk, Alar, Balagansk, Barguzin, Ida, Khorinsk, Kuda, Selenga, Verkholensk, Olkhon, Tunka, and Shenehen Buryats) groups. The analysis of molecular variation performed using regional, tribal, and ethno-territorial divisions of the Buryats showed lack of genetic differentiation at all levels. Nonetheless, the complete mitogenome analysis revealed a very high level of genetic diversity in the Buryats which is the highest among Siberian populations and comparable to that in populations of eastern and western Asia. The AMOVA and MDS analyses results imply to a strong genetic similarity between the Buryats and eastern Asian populations of Chinese and Japanese, suggesting their origin on the basis of common maternal ancestry components. Several new Buryat-specific branches of haplogroup G (G2a2a, G2a1i, G2a5a) display signals of dispersals dating to 2.6-6.6 kya with a possible origin in eastern Asia, thus testifying Bronze Age and Neolithic arrival of ancestral eastern Asian component to the South Siberia region.

Mitogenomic diversity in Russians and Poles

Complete mtDNA genome sequencing improves molecular resolution for distinguishing variation between individuals and populations, but there is still deficiency of mitogenomic population data. To overcome this limitation, we used Sanger-based protocol to generate complete mtDNA sequences of 376 Russian individuals from six populations of European part of Russia and 100 Polish individuals from northern Poland. Nearly complete resolution of mtDNA haplotypes was achieved - about 97% of haplotypes were unique both in Russians and Poles, and no haplotypes overlapped between them when indels were considered. While European populations showed a low, but statistically significant level of between-population differentiation (Fst=0.66%, p=0), Russians demonstrate lack of between-population differences (Fst=0.22%, p=0.15). Results of the Bayesian skyline analysis of Russian mitogenomes demonstrate not only post-Last Glacial Maximum expansion, but also rapid population growth starting from about 4.3kya (95% CI: 2.9-5.8kya), i.e. in the Bronze Age. This expansion strongly correlates with the Kurgan model established by archaeologists and confirmed by paleogeneticists.

Mitochondrial super-haplogroup U diversity in Serbians

BACKGROUND

Available mitochondrial (mtDNA) data demonstrate genetic differentiation among South Slavs inhabiting the Balkan Peninsula. However, their resolution is insufficient to elucidate the female-specific aspects of the genetic history of South Slavs, including the genetic impact of various migrations which were rather common within the Balkans, a region having a turbulent demographic history.

AIM

The aim was to thoroughly study complete mitogenomes of Serbians, a population linking westward and eastward South Slavs.

SUBJECTS AND METHODS

Forty-six predominantly Serbian super-haplogroup U complete mitogenomes were analysed phylogenetically against ∼4000 available complete mtDNAs of modern and ancient Western Eurasians.

RESULTS

Serbians share a number of U mtDNA lineages with Southern, Eastern-Central and North-Western Europeans. Putative Balkan-specific lineages (e.g. U1a1c2, U4c1b1, U5b3j, K1a4l and K1a13a1) and lineages shared among Serbians (South Slavs) and West and East Slavs were detected (e.g. U2e1b1, U2e2a1d, U4a2a, U4a2c, U4a2g1, U4d2b and U5b1a1).

CONCLUSION

The exceptional diversity of maternal lineages found in Serbians may be associated with the genetic impact of both autochthonous pre-Slavic Balkan populations whose mtDNA gene pool was affected by migrations of various populations over time (e.g. Bronze Age pastoralists) and Slavic and Germanic newcomers in the early Middle Ages.

Selective sweep on human amylase genes postdates the split with Neanderthals

Humans have more copies of amylase genes than other primates. It is still poorly understood, however, when the copy number expansion occurred and whether its spread was enhanced by selection. Here we assess amylase copy numbers in a global sample of 480 high coverage genomes and find that regions flanking the amylase locus show notable depression of genetic diversity both in African and non-African populations. Analysis of genetic variation in these regions supports the model of an early selective sweep in the human lineage after the split of humans from Neanderthals which led to the fixation of multiple copies of AMY1 in place of a single copy. We find evidence of multiple secondary losses of copy number with the highest frequency (52%) of a deletion of AMY2A and associated low copy number of AMY1 in Northeast Siberian populations whose diet has been low in starch content.

East Eurasian ancestry in the middle of Europe: genetic footprints of Steppe nomads in the genomes of Belarusian Lipka Tatars

Medieval era encounters of nomadic groups of the Eurasian Steppe and largely sedentary East Europeans had a variety of demographic and cultural consequences. Amongst these outcomes was the emergence of the Lipka Tatars—a Slavic-speaking Sunni-Muslim minority residing in modern Belarus, Lithuania and Poland, whose ancestors arrived in these territories via several migration waves, mainly from the Golden Horde. Our results show that Belarusian Lipka Tatars share a substantial part of their gene pool with Europeans as indicated by their Y-chromosomal, mitochondrial and autosomal DNA variation. Nevertheless, Belarusian Lipkas still retain a strong genetic signal of their nomadic ancestry, witnessed by the presence of common Y-chromosomal and mitochondrial DNA variants as well as autosomal segments identical by descent between Lipkas and East Eurasians from temperate and northern regions. Hence, we document Lipka Tatars as a unique example of former Medieval migrants into Central Europe, who became sedentary, changed language to Slavic, yet preserved their faith and retained, both uni- and bi-parentally, a clear genetic echo of a complex population interplay throughout the Eurasian Steppe Belt, extending from Central Europe to northern China.

POPULATION GENETICS. Genomic evidence for the Pleistocene and recent population history of Native Americans

How and when the Americas were populated remains contentious. Using ancient and modern genome-wide data, we found that the ancestors of all present-day Native Americans, including Athabascans and Amerindians, entered the Americas as a single migration wave from Siberia no earlier than 23 thousand years ago (ka) and after no more than an 8000-year isolation period in Beringia. After their arrival to the Americas, ancestral Native Americans diversified into two basal genetic branches around 13 ka, one that is now dispersed across North and South America and the other restricted to North America. Subsequent gene flow resulted in some Native Americans sharing ancestry with present-day East Asians (including Siberians) and, more distantly, Australo-Melanesians. Putative "Paleoamerican" relict populations, including the historical Mexican Pericúes and South American Fuego-Patagonians, are not directly related to modern Australo-Melanesians as suggested by the Paleoamerican Model.

The genetic legacy of the expansion of Turkic-speaking nomads across Eurasia

The Turkic peoples represent a diverse collection of ethnic groups defined by the Turkic languages. These groups have dispersed across a vast area, including Siberia, Northwest China, Central Asia, East Europe, the Caucasus, Anatolia, the Middle East, and Afghanistan. The origin and early dispersal history of the Turkic peoples is disputed, with candidates for their ancient homeland ranging from the Transcaspian steppe to Manchuria in Northeast Asia. Previous genetic studies have not identified a clear-cut unifying genetic signal for the Turkic peoples, which lends support for language replacement rather than demic diffusion as the model for the Turkic language’s expansion. We addressed the genetic origin of 373 individuals from 22 Turkic-speaking populations, representing their current geographic range, by analyzing genome-wide high-density genotype data. In agreement with the elite dominance model of language expansion most of the Turkic peoples studied genetically resemble their geographic neighbors. However, western Turkic peoples sampled across West Eurasia shared an excess of long chromosomal tracts that are identical by descent (IBD) with populations from present-day South Siberia and Mongolia (SSM), an area where historians center a series of early Turkic and non-Turkic steppe polities. While SSM matching IBD tracts (> 1cM) are also observed in non-Turkic populations, Turkic peoples demonstrate a higher percentage of such tracts (p-values ≤ 0.01) compared to their non-Turkic neighbors. Finally, we used the ALDER method and inferred admixture dates (~9th–17th centuries) that overlap with the Turkic migrations of the 5th–16th centuries. Thus, our results indicate historical admixture among Turkic peoples, and the recent shared ancestry with modern populations in SSM supports one of the hypothesized homelands for their nomadic Turkic and related Mongolic ancestors.

Centuries of nomadic migrations have ultimately resulted in the distribution of Turkic languages over a large area ranging from Siberia, across Central Asia to Eastern Europe and the Middle East. Despite the profound cultural impact left by these nomadic peoples, little is known about their prehistoric origins. Moreover, because contemporary Turkic speakers tend to genetically resemble their geographic neighbors, it is not clear whether their nomadic ancestors left an identifiable genetic trace. In this study, we show that Turkic-speaking peoples sampled across the Middle East, Caucasus, East Europe, and Central Asia share varying proportions of Asian ancestry that originate in a single area, southern Siberia and Mongolia. Mongolic- and Turkic-speaking populations from this area bear an unusually high number of long chromosomal tracts that are identical by descent with Turkic peoples from across west Eurasia. Admixture induced linkage disequilibrium decay across chromosomes in these populations indicates that admixture occurred during the 9th–17th centuries, in agreement with the historically recorded Turkic nomadic migrations and later Mongol expansion. Thus, our findings reveal genetic traces of recent large-scale nomadic migrations and map their source to a previously hypothesized area of Mongolia and southern Siberia.

A recent bottleneck of Y chromosome diversity coincides with a global change in culture

It is commonly thought that human genetic diversity in non-African populations was shaped primarily by an out-of-Africa dispersal 50–100 thousand yr ago (kya). Here, we present a study of 456 geographically diverse high-coverage Y chromosome sequences, including 299 newly reported samples. Applying ancient DNA calibration, we date the Y-chromosomal most recent common ancestor (MRCA) in Africa at 254 (95% CI 192–307) kya and detect a cluster of major non-African founder haplogroups in a narrow time interval at 47–52 kya, consistent with a rapid initial colonization model of Eurasia and Oceania after the out-of-Africa bottleneck. In contrast to demographic reconstructions based on mtDNA, we infer a second strong bottleneck in Y-chromosome lineages dating to the last 10 ky. We hypothesize that this bottleneck is caused by cultural changes affecting variance of reproductive success among males.

Mitochondrial DNA perspective of Serbian genetic diversity

Although south-Slavic populations have been studied to date from various aspects, the population of Serbia, occupying the central part of the Balkan Peninsula, is still genetically understudied at least at the level of mitochondrial DNA (mtDNA) variation. We analyzed polymorphisms of the first and the second mtDNA hypervariable segments (HVS-I and HVS-II) and informative coding-region markers in 139 Serbians to shed more light on their mtDNA variability, and used available data on other Slavic and neighboring non-Slavic populations to assess their interrelations in a broader European context. The contemporary Serbian mtDNA profile is consistent with the general European maternal landscape having a substantial proportion of shared haplotypes with eastern, central, and southern European populations. Serbian population was characterized as an important link between easternmost and westernmost south-Slavic populations due to the observed lack of genetic differentiation with all other south-Slavic populations and its geographical positioning within the Balkan Peninsula. An increased heterogeneity of south Slavs, most likely mirroring turbulent demographic events within the Balkan Peninsula over time (i.e., frequent admixture and differential introgression of various gene pools), and a marked geographical stratification of Slavs to south-, east-, and west-Slavic groups, were also found. A phylogeographic analyses of 20 completely sequenced Serbian mitochondrial genomes revealed not only the presence of mtDNA lineages predominantly found within the Slavic gene pool (U4a2a*, U4a2a1, U4a2c, U4a2g, HV10), supporting a common Slavic origin, but also lineages that may have originated within the southern Europe (H5*, H5e1, H5a1v) and the Balkan Peninsula in particular (H6a2b and L2a1k).

A novel multiplex assay amplifying 13 Y-STRs characterized by rapid and moderate mutation rate

As microsatellites located on Y chromosome mutate with different rates, they may be exploited in evolutionary studies, genealogical testing of a variety of populations and even, as proven recently, aid individual identification. Currently available commercial Y-STR kits encompass mostly low to moderately mutating loci, making them a perfect choice for the first two applications. Some attempts have been made so far to utilize Y-STRs to provide a discriminatory tool for forensic purposes. Although all 13 rapidly mutating Y-STRs were already multiplexed, no single assay based on single-copy markers allowing at least a portion of close male relatives to be differentiated from one another is available. To fill in the blanks, we constructed and validated an assay comprised of single-copy Y-STR markers only with a mutation rate ranging from 8×10(-3) to 1×10(-2). Performance of the resulting combination of nine RM Y-STRs and four moderately mutating ones was tested on 361 father-son pairs and 1326 males from 9 populations revealing an overall mutation rate of 1.607×10(-1) for the assay as a whole. Application of the proposed 13 Y-STR set to differentiation of haplotypes present among homogenous population of Buryats resulted in a threefold increase of discrimination as compared with 10 Y-STRs from the PowerPlex(®) Y.

A Selective Sweep on a Deleterious Mutation in CPT1A in Arctic Populations

Arctic populations live in an environment characterized by extreme cold and the absence of plant foods for much of the year and are likely to have undergone genetic adaptations to these environmental conditions in the time they have been living there. Genome-wide selection scans based on genotype data from native Siberians have previously highlighted a 3 Mb chromosome 11 region containing 79 protein-coding genes as the strongest candidates for positive selection in Northeast Siberians. However, it was not possible to determine which of the genes might be driving the selection signal. Here, using whole-genome high-coverage sequence data, we identified the most likely causative variant as a nonsynonymous G>A transition (rs80356779; c.1436C>T [p.Pro479Leu] on the reverse strand) in CPT1A, a key regulator of mitochondrial long-chain fatty-acid oxidation. Remarkably, the derived allele is associated with hypoketotic hypoglycemia and high infant mortality yet occurs at high frequency in Canadian and Greenland Inuits and was also found at 68% frequency in our Northeast Siberian sample. We provide evidence of one of the strongest selective sweeps reported in humans; this sweep has driven this variant to high frequency in circum-Arctic populations within the last 6-23 ka despite associated deleterious consequences, possibly as a result of the selective advantage it originally provided to either a high-fat diet or a cold environment.

Western Eurasian ancestry in modern Siberians based on mitogenomic data

Background

Although the genetic heritage of aboriginal Siberians is mostly of eastern Asian ancestry, a substantial western Eurasian component is observed in the majority of northern Asian populations. Traces of at least two migrations into southern Siberia, one from eastern Europe and the other from western Asia/the Caucasus have been detected previously in mitochondrial gene pools of modern Siberians.

Results

We report here 166 new complete mitochondrial DNA (mtDNA) sequences that allow us to expand and re-analyze the available data sets of western Eurasian lineages found in northern Asian populations, define the phylogenetic status of Siberian-specific subclades and search for links between mtDNA haplotypes/subclades and events of human migrations. From a survey of 158 western Eurasian mtDNA genomes found in Siberia we estimate that nearly 40% of them most likely have western Asian and another 29% European ancestry. It is striking that 65 of northern Asian mitogenomes, i.e. ~41%, fall into 19 branches and subclades which can be considered as Siberian-specific being found so far only in Siberian populations. From the coalescence analysis it is evident that the sequence divergence of Siberian-specific subclades was relatively small, corresponding to only 0.6-9.5 kya (using the complete mtDNA rate) and 1–6 kya (coding region rate).

Conclusions

The phylogeographic analysis implies that the western Eurasian founders, giving rise to Siberian specific subclades, may trace their ancestry only to the early and mid-Holocene, though some of genetic lineages may trace their ancestry back to the end of Last Glacial Maximum (LGM). We have not found the modern northern Asians to have western Eurasian genetic components of sufficient antiquity to indicate traces of pre-LGM expansions.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-014-0217-9) contains supplementary material, which is available to authorized users.

A mitogenomic phylogeny and genetic history of sable (Martes zibellina)

We assessed phylogeny of sable (Martes zibellina, Linnaeus, 1758) by sequence analysis of nearly complete, new mitochondrial genomes in 36 specimens from different localities in northern Eurasia (Primorye, Khabarovsk and Krasnoyarsk regions, the Kamchatka Peninsula, the Kuril Islands and the Urals). Phylogenetic analysis of mtDNA sequences demonstrates that two clades, A and BC, radiated about 200-300 thousandyears ago (kya) according to results of Bayesian molecular clock and RelTime analyses of different mitogenome alignments (nearly complete mtDNA sequences, protein-coding region, and synonymous sites), while the age estimates of clades A, B and C fall within the Late Pleistocene (~50-140 kya). Bayesian skyline plots (BSPs) of sable population size change based on analysis of nearly complete mtDNAs show an expansion around 40 kya in the warm Karganian time, without a decline of population size around the Last Glacial Maximum (21 kya). The BSPs based on synonymous clock rate indicate that M. zibellina experienced demographic expansions later, approximately 22 kya. The A2a clade that colonized Kamchatka ~23-50 kya (depending on the mutation rate used) survived the last glaciation there as demonstrated by the BSP analysis. In addition, we have found evidence of positive selection acting at ND4 and cytochrome b genes, thereby suggesting adaptive evolution of the A2a clade in Kamchatka.

Complete mitochondrial DNA diversity in Iranians

Due to its pivotal geographical location and proximity to transcontinental migratory routes, Iran has played a key role in subsequent migrations, both prehistoric and historic, between Africa, Asia and Europe. To shed light on the genetic structure of the Iranian population as well as on the expansion patterns and population movements which affected this region, the complete mitochondrial genomes of 352 Iranians were obtained. All Iranian populations studied here exhibit similarly high diversity values comparable to the other groups from the Caucasus, Anatolia and Europe. The results of AMOVA and MDS analyses did not associate any regional and/or linguistic group of populations in the Anatolia/Caucasus and Iran region pointing to close genetic positions of Persians and Qashqais to each other and to Armenians, and Azeris from Iran to Georgians. By reconstructing the complete mtDNA phylogeny of haplogroups R2, N3, U1, U3, U5a1g, U7, H13, HV2, HV12, M5a and C5c we have found a previously unexplored genetic connection between the studied Iranian populations and the Arabian Peninsula, India, Near East and Europe, likely the result of both ancient and recent gene flow. Our results for Persians and Qashqais point to a continuous increase of the population sizes from ∼24 kya to the present, although the phase between 14-24 kya is thought to be hyperarid according to the Gulf Oasis model. Since this would have affected hunter-gatherer ranges and mobility patterns and forced them to increasingly rely on coastal resources, this transition can explain the human expansion across the Persian Gulf region.

Phylogeny and genetic history of the Siberian salamander (Salamandrella keyserlingii, Dybowski, 1870) inferred from complete mitochondrial genomes

We assessed phylogeny of the Siberian salamander (Salamandrella keyserlingii, Dybowski, 1870), the most northern ectothermic, terrestrial vertebrate in Eurasia, by sequence analysis of complete mitochondrial genomes in 26 specimens from different localities (China, Khabarovsk region, Sakhalin, Yakutia, Magadan region, Chukotka, Kamchatka, Ural, European part of Russia). In addition, a complete mitochondrial genome of the Schrenck salamander, Salamandrella schrenckii, was determined for the first time. Bayesian phylogenetic analysis of the entire mtDNA genomes of S. keyserlingii demonstrates that two haplotype clades, AB and C, radiated about 1.4 million years ago (Mya). Bayesian skyline plots of population size change through time show an expansion around 250 thousand years ago (kya) and then a decline around the Last Glacial Maximum (25 kya) with subsequent restoration of population size. Climatic changes during the Quaternary period have dramatically affected the population genetic structure of the Siberian salamanders. In addition, complete mtDNA sequence analysis allowed us to recognize that the vast area of Northern Eurasia was colonized only by the Siberian salamander clade C1b during the last 150 kya. Meanwhile, we were unable to find evidence of molecular adaptation in this clade by analyzing the whole mitochondrial genomes of the Siberian salamanders.

The history of Slavs inferred from complete mitochondrial genome sequences

To shed more light on the processes leading to crystallization of a Slavic identity, we investigated variability of complete mitochondrial genomes belonging to haplogroups H5 and H6 (63 mtDNA genomes) from the populations of Eastern and Western Slavs, including new samples of Poles, Ukrainians and Czechs presented here. Molecular dating implies formation of H5 approximately 11.5–16 thousand years ago (kya) in the areas of southern Europe. Within ancient haplogroup H6, dated at around 15–28 kya, there is a subhaplogroup H6c, which probably survived the last glaciation in Europe and has undergone expansion only 3–4 kya, together with the ancestors of some European groups, including the Slavs, because H6c has been detected in Czechs, Poles and Slovaks. Detailed analysis of complete mtDNAs allowed us to identify a number of lineages that seem specific for Central and Eastern Europe (H5a1f, H5a2, H5a1r, H5a1s, H5b4, H5e1a, H5u1, some subbranches of H5a1a and H6a1a9). Some of them could possibly be traced back to at least ∼4 kya, which indicates that some of the ancestors of today's Slavs (Poles, Czechs, Slovaks, Ukrainians and Russians) inhabited areas of Central and Eastern Europe much earlier than it was estimated on the basis of archaeological and historical data. We also sequenced entire mitochondrial genomes of several non-European lineages (A, C, D, G, L) found in contemporary populations of Poland and Ukraine. The analysis of these haplogroups confirms the presence of Siberian (C5c1, A8a1) and Ashkenazi-specific (L2a1l2a) mtDNA lineages in Slavic populations. Moreover, we were able to pinpoint some lineages which could possibly reflect the relatively recent contacts of Slavs with nomadic Altaic peoples (C4a1a, G2a, D5a2a1a1).

Y-chromosome variation in Tajiks and Iranians

AIM

The purpose of this study was to characterize Y-chromosome diversity in Tajiks from Tajikistan and in Persians and Kurds from Iran.

METHOD

Y-chromosome haplotypes were identified in 40 Tajiks, 77 Persians and 25 Kurds, using 12 short tandem repeats (STR) and 18 binary markers.

RESULTS

High genetic diversity was observed in the populations studied. Six of 12 haplogroups were common in Persians, Kurds and Tajiks, but only three haplogroups (G-M201, J-12f2 and L-M20) were the most frequent in all populations, comprising together ~60% of the Y-chromosomes in the pooled data set. Analysis of genetic distances between Y-STR haplotypes revealed that the Kurds showed a great distance to the Iranian-speaking populations of Iran, Afghanistan and Tajikistan. The presence of Indian-specific haplogroups L-M20, H1-M52 and R2a-M124 in both Tajik samples from Afghanistan and Tajikistan demonstrates an apparent genetic affinity between Tajiks from these two regions.

CONCLUSIONS

Despite the marked similarities between Y-chromosome gene pools of Iranian-speaking populations, there are differences between them, defined by many factors, including geographic and linguistic relationships.

On the Y-chromosome haplogroup C3c classification

As there are ambiguities in classification of the Y-chromosome haplogroup C3c, relatively frequent in populations of Northern Asia, we analyzed all three haplogroup-defining markers M48, M77 and M86 in C3-M217-individuals from Siberia, Eastern Asia and Eastern Europe. We have found that haplogroup C3c is characterized by the derived state at M48, whereas mutations at both M77 and M86 define subhaplogroup C3c1. The branch defined by M48 alone would belong to subhaplogroup C3c*, characteristic for some populations of Central and Eastern Siberia, such as Koryaks, Evens, Evenks and Yukaghirs. Subhaplogroup C3c* individuals could be considered as remnants of the Neolithic population of Siberia, based on the age of C3c*-short tandem repeat variation amounting to 4.5 ± 2.4 thousand years.

Complete mitochondrial DNA analysis of eastern Eurasian haplogroups rarely found in populations of northern Asia and eastern Europe

With the aim of uncovering all of the most basal variation in the northern Asian mitochondrial DNA (mtDNA) haplogroups, we have analyzed mtDNA control region and coding region sequence variation in 98 Altaian Kazakhs from southern Siberia and 149 Barghuts from Inner Mongolia, China. Both populations exhibit the prevalence of eastern Eurasian lineages accounting for 91.9% in Barghuts and 60.2% in Altaian Kazakhs. The strong affinity of Altaian Kazakhs and populations of northern and central Asia has been revealed, reflecting both influences of central Asian inhabitants and essential genetic interaction with the Altai region indigenous populations. Statistical analyses data demonstrate a close positioning of all Mongolic-speaking populations (Mongolians, Buryats, Khamnigans, Kalmyks as well as Barghuts studied here) and Turkic-speaking Sojots, thus suggesting their origin from a common maternal ancestral gene pool. In order to achieve a thorough coverage of DNA lineages revealed in the northern Asian matrilineal gene pool, we have completely sequenced the mtDNA of 55 samples representing haplogroups R11b, B4, B5, F2, M9, M10, M11, M13, N9a and R9c1, which were pinpointed from a massive collection (over 5000 individuals) of northern and eastern Asian, as well as European control region mtDNA sequences. Applying the newly updated mtDNA tree to the previously reported northern Asian and eastern Asian mtDNA data sets has resolved the status of the poorly classified mtDNA types and allowed us to obtain the coalescence age estimates of the nodes of interest using different calibrated rates. Our findings confirm our previous conclusion that northern Asian maternal gene pool consists of predominantly post-LGM components of eastern Asian ancestry, though some genetic lineages may have a pre-LGM/LGM origin.

Ancient links between Siberians and Native Americans revealed by subtyping the Y chromosome haplogroup Q1a

To investigate the structure of Y chromosome haplogroups R-M207 and Q-M242 in human populations of North Asia, we have performed high-resolution genotyping using both single nucleotide polymorphisms and short tandem repeat (STR)-based approaches of 121 M207- and M242-derived samples from 885 males of 16 ethnic groups of Siberia and East Asia. As a result, the following Y chromosome haplogroups were revealed: R1b1b1-M73 (2.0%), R1b1b2-M269 (0.7%), R2-M124 (1.1%), Q1a*-MEH2 (0.5%), Q1a2-M25 (0.1%), Q1a3*-M346 (9.2%) and Q1a3a-M3 (0.2%). Despite the low coalescence age of haplogroup Q1a3*-M346, which is estimated in South Siberia as about 4.5±1.5 thousand years ago (Ka), divergence time between these Q1a3*-M346 haplotypes and Amerindian-specific haplogroup Q1a3a-M3 is equal to 13.8±3.9 Ka, pointing to a relatively recent entry date to America. In addition, unique cluster of haplotypes belonging to Q1a*-MEH2 was found in Koryaks inhabiting the Sea of Okhotsk coast (at a frequency of 10.3%). Although the level of STR diversity associated with Q1a*-MEH2 is very low, this lineage appears to be closest to the extinct Palaeo-Eskimo individuals belonging to the Saqqaq culture arisen in the New World Arctic about 5.5 Ka. This finding suggests that Q1a*-MEH2 likely traces a population migration originating in Northeast Siberia across the Bering Strait.

Origin and post-glacial dispersal of mitochondrial DNA haplogroups C and D in northern Asia

More than a half of the northern Asian pool of human mitochondrial DNA (mtDNA) is fragmented into a number of subclades of haplogroups C and D, two of the most frequent haplogroups throughout northern, eastern, central Asia and America. While there has been considerable recent progress in studying mitochondrial variation in eastern Asia and America at the complete genome resolution, little comparable data is available for regions such as southern Siberia--the area where most of northern Asian haplogroups, including C and D, likely diversified. This gap in our knowledge causes a serious barrier for progress in understanding the demographic pre-history of northern Eurasia in general. Here we describe the phylogeography of haplogroups C and D in the populations of northern and eastern Asia. We have analyzed 770 samples from haplogroups C and D (174 and 596, respectively) at high resolution, including 182 novel complete mtDNA sequences representing haplogroups C and D (83 and 99, respectively). The present-day variation of haplogroups C and D suggests that these mtDNA clades expanded before the Last Glacial Maximum (LGM), with their oldest lineages being present in the eastern Asia. Unlike in eastern Asia, most of the northern Asian variants of haplogroups C and D began the expansion after the LGM, thus pointing to post-glacial re-colonization of northern Asia. Our results show that both haplogroups were involved in migrations, from eastern Asia and southern Siberia to eastern and northeastern Europe, likely during the middle Holocene.

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.

Mitochondrial haplogroup N1a phylogeography, with implication to the origin of European farmers

BACKGROUND

Tracing the genetic origin of central European farmer N1a lineages can provide a unique opportunity to assess the patterns of the farming technology spread into central Europe in the human prehistory. Here, we have chosen twelve N1a samples from modern populations which are most similar with the farmer N1a types and performed the complete mitochondrial DNA genome sequencing analysis. To assess the genetic and phylogeographic relationship, we performed a detailed survey of modern published N1a types from Eurasian and African populations.

RESULTS

The geographic origin and expansion of farmer lineages related N1a subclades have been deduced from combined analysis of 19 complete sequences with 166 N1a haplotypes. The phylogeographic analysis revealed that the central European farmer lineages have originated from different sources: from eastern Europe, local central Europe, and from the Near East via southern Europe.

CONCLUSIONS

The results obtained emphasize that the arrival of central European farmer lineages did not occur via a single demic diffusion event from the Near East at the onset of the Neolithic spread of agriculture into Europe. Indeed these results indicate that the Neolithic transition process was more complex in central Europe and possibly the farmer N1a lineages were a result of a 'leapfrog' colonization process.

Mitogenomic diversity in Tatars from the Volga-Ural region of Russia

To investigate diversity of mitochondrial gene pool of Tatars inhabiting the territory of the middle Volga River basin, 197 individuals from two populations representing Kazan Tatars and Mishars were subjected for analysis of mitochondrial DNA (mtDNA) control region variation. In addition, 73 mitochondrial genomes of individuals from Mishar population were sequenced completely. It was found that mitochondrial gene pool of the Volga Tatars consists of two parts, but western Eurasian component prevails considerably (84% on average) over eastern Asian one (16%). Eastern Asian mtDNAs detected in Tatars belonged to a heterogeneous set of haplogroups (A, C, D, G, M7, M10, N9a, Y, and Z), although only haplogroups A and D were revealed simultaneously in both populations. Complete mtDNA variation study revealed that the age of western Eurasian haplogroups (such as U4, HV0a, and H) is less than 18,000 years, thus suggesting re-expansion of eastern Europeans soon after the Last Glacial Maximum.

The peopling of Europe from the mitochondrial haplogroup U5 perspective

It is generally accepted that the most ancient European mitochondrial haplogroup, U5, has evolved essentially in Europe. To resolve the phylogeny of this haplogroup, we completely sequenced 113 mitochondrial genomes (79 U5a and 34 U5b) of central and eastern Europeans (Czechs, Slovaks, Poles, Russians and Belorussians), and reconstructed a detailed phylogenetic tree, that incorporates previously published data. Molecular dating suggests that the coalescence time estimate for the U5 is ∼25–30 thousand years (ky), and ∼16–20 and ∼20–24 ky for its subhaplogroups U5a and U5b, respectively. Phylogeographic analysis reveals that expansions of U5 subclusters started earlier in central and southern Europe, than in eastern Europe. In addition, during the Last Glacial Maximum central Europe (probably, the Carpathian Basin) apparently represented the area of intermingling between human flows from refugial zones in the Balkans, the Mediterranean coastline and the Pyrenees. Age estimations amounting for many U5 subclusters in eastern Europeans to ∼15 ky ago and less are consistent with the view that during the Ice Age eastern Europe was an inhospitable place for modern humans.

Phylogeography and molecular adaptation of Siberian salamander Salamandrella keyserlingii based on mitochondrial DNA variation

We assessed the phylogeographic pattern of Siberian salamander (Salamandrella keyserlingii, Dybowski, 1870), which appear to be the most northern ectothermic, terrestrial vertebrate in Northern Eurasia, by sequence analysis of a 611-bp fragment of the mitochondrial cytochrome b gene in 159 specimens from different localities (Khabarovsk region, Sakhalin, Yakutia, Magadan region, Chukotka, Kamchatka and others). The data revealed that cytochrome b lineages of S. keyserlingii are divided into haplogroups A, B and C. Haplogroup A and B sequences are widespread in the Far East region, whereas haplogroup C consisting of several phylogenetic clusters (C1, C2, C3) is present in the all range of S. keyserlingii. Among them, cluster C3 appears to be specific for Sakhalin; most likely, it has arisen in situ in this island, with the entry time of the founder mtDNA estimated at about 0.4 MY. Analysis of cytochrome b gene variation by using different neutrality tests (including those based on K(A)/K(S)-ratio) has shown that differences between haplogroups were statistically insignificant, thus suggesting selective neutrality. However, analysis of amino acid changes allowed us to detect a signature of molecular adaptation, which might have led to appearance of adaptive cytochrome b variants in haplogroup C, originating most likely at the end of Eopleistocene (about 0.64 MY based on the haplogroup C divergence level). It seems probable that this adaptive mechanism could promote subsequent populating of new regions.

Mitochondrial DNA phylogeny in Eastern and Western Slavs

To resolve the phylogeny of certain mitochondrial DNA (mtDNA) haplogroups in eastern Europe and estimate their evolutionary age, a total of 73 samples representing mitochondrial haplogroups U4, HV*, and R1 were selected for complete mitochondrial genome sequencing from a collection of about 2,000 control region sequences sampled in eastern (Russians, Belorussians, and Ukrainians) and western (Poles, Czechs, and Slovaks) Slavs. On the basis of whole-genome resolution, we fully characterized a number of haplogroups (HV3, HV4, U4a1, U4a2, U4a3, U4b, U4c, U4d, and R1a) that were previously described only partially. Our findings demonstrate that haplogroups HV3, HV4, and U4a1 could be traced back to the pre-Neolithic times ( approximately 12,000-19,000 years before present [YBP]) in eastern Europe. In addition, an ancient connection between the Caucasus/Europe and India has been revealed by analysis of haplogroup R1 diversity, with a split between the Indian and Caucasus/European R1a lineages occurring about 16,500 years ago. Meanwhile, some mtDNA subgroups detected in Slavs (such as U4a2a, U4a2*, HV3a, and R1a1) are definitely younger being dated between 6,400 and 8,200 YBP. However, robust age estimations appear to be problematic due to the high ratios of nonsynonymous to synonymous substitutions found in young mtDNA subclusters.

Phylogeographic analysis of mitochondrial DNA in northern Asian populations

To elucidate the human colonization process of northern Asia and human dispersals to the Americas, a diverse subset of 71 mitochondrial DNA (mtDNA) lineages was chosen for complete genome sequencing from the collection of 1,432 control-region sequences sampled from 18 autochthonous populations of northern, central, eastern, and southwestern Asia. On the basis of complete mtDNA sequencing, we have revised the classification of haplogroups A, D2, G1, M7, and I; identified six new subhaplogroups (I4, N1e, G1c, M7d, M7e, and J1b2a); and fully characterized haplogroups N1a and G1b, which were previously described only by the first hypervariable segment (HVS1) sequencing and coding-region restriction-fragment-length polymorphism analysis. Our findings indicate that the southern Siberian mtDNA pool harbors several lineages associated with the Late Upper Paleolithic and/or early Neolithic dispersals from both eastern Asia and southwestern Asia/southern Caucasus. Moreover, the phylogeography of the D2 lineages suggests that southern Siberia is likely to be a geographical source for the last postglacial maximum spread of this subhaplogroup to northern Siberia and that the expansion of the D2b branch occurred in Beringia ~7,000 years ago. In general, a detailed analysis of mtDNA gene pools of northern Asians provides the additional evidence to rule out the existence of a northern Asian route for the initial human colonization of Asia.

Y-chromosome haplogroup N dispersals from south Siberia to Europe

In order to reconstruct the history of Y-chromosome haplogroup (hg) N dispersals in north Eurasia, we have analyzed the diversity of microsatellite (STR) loci within two major hg N clades, N2 and N3, in a total sample of 1,438 males from 17 ethnic groups, mainly of Siberian and Eastern European origin. Based on STR variance analysis we observed that hg N3a is more diverse in Eastern Europe than in south Siberia. However, analysis of median networks showed that there are two STR subclusters of hg N3a, N3a1 and N3a2, that are characterized by different genetic histories. Age calculation of STR variation within subcluster N3a1 indicated that its first expansion occurred in south Siberia [approximately 10,000 years (ky)] and then this subcluster spread into Eastern Europe where its age is around 8 ky ago. Meanwhile, younger subcluster N3a2 originated in south Siberia (probably in the Baikal region) approximately 4 ky ago. Median network and variance analyses of STR haplotypes suggest that south Siberian N3a2 haplotypes spread further into Volga-Ural region undergoing serial bottlenecks. In addition, median network analysis of STR data demonstrates that haplogroup N2-A is represented by two subclusters, showing recent expansion times. The data obtained allow us to suggest Siberian origin of haplogroups N3 and N2 that are currently widespread in some populations of Eastern Europe.

Allelic and haplotypic frequencies at 11 Y-STR loci in Buryats from South-East Siberia

We have obtained Y-STR haplotypes in 12 loci (DYS19, DYS385, DYS389I, DYS389II, DYS390, DYS391, DYS392, DYS393, DYS437, DYS438 and DYS439) from 215 Buryat males. We have found that one haplotype (15-11,18-13-28-23-10-11-14-14-10-12) comprises more than 30% of Y chromosomes in this population while another haplotype (14-11,13-14-30-23-10-14-14-14-10-10) comprises additional 14% of chromosomes. The population under study seems to be very homogenous as far as Y chromosome is regarded and the most frequent haplotype seems to be the modal haplotype for Buryats.

Contrasting patterns of Y-chromosome variation in South Siberian populations from Baikal and Altai-Sayan regions

In order to investigate the genetic history of autochthonous South Siberian populations and to estimate the contribution of distinct patrilineages to their gene pools, we have analyzed 17 Y-chromosomal binary markers (YAP, RPS4Y(711), SRY-8299, M89, M201, M52, M170, 12f2, M9, M20, 92R7, SRY-1532, DYS199, M173, M17, Tat, and LLY22 g) in a total sample of 1,358 males from 14 ethnic groups of Siberia (Altaians-Kizhi, Teleuts, Shors, Tuvinians, Todjins, Tofalars, Sojots, Khakassians, Buryats, Evenks), Central/Eastern Asia (Mongolians and Koreans) and Eastern Europe (Kalmyks and Russians). Based on both, the distribution pattern of Y-chromosomal haplogroups and results on AMOVA analysis we observed the statistically significant genetic differentiation between the populations of Baikal and Altai-Sayan regions. We suggest that these regional differences can be best explained by different contribution of Central/Eastern Asian and Eastern European paternal lineages into gene pools of modern South Siberians. The population of the Baikal region demonstrates the prevalence of Central/Eastern Asian lineages, whereas in the populations of Altai and Sayan regions the highest paternal contribution resulted from Eastern European descent is revealed. Yet, our data on Y-chromosome STRs variation demonstrate the clear differences between the South Siberian and Eastern European R1a1-lineages with the evolutionary ages compatible with divergence time between these two regional groups.

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.

Origin and diffusion of mtDNA haplogroup X

A maximum parsimony tree of 21 complete mitochondrial DNA (mtDNA) sequences belonging to haplogroup X and the survey of the haplogroup-associated polymorphisms in 13,589 mtDNAs from Eurasia and Africa revealed that haplogroup X is subdivided into two major branches, here defined as "X1" and "X2." The first is restricted to the populations of North and East Africa and the Near East, whereas X2 encompasses all X mtDNAs from Europe, western and Central Asia, Siberia, and the great majority of the Near East, as well as some North African samples. Subhaplogroup X1 diversity indicates an early coalescence time, whereas X2 has apparently undergone a more recent population expansion in Eurasia, most likely around or after the last glacial maximum. It is notable that X2 includes the two complete Native American X sequences that constitute the distinctive X2a clade, a clade that lacks close relatives in the entire Old World, including Siberia. The position of X2a in the phylogenetic tree suggests an early split from the other X2 clades, likely at the very beginning of their expansion and spread from the Near East.