Structure of the gene pool of eastern Ukrainians from Y-chromosome haplogroups

Y chromosomes from representative sample of Eastern Ukrainians (94 individuals) were analyzed for composition and frequencies of haplogroups, defined by 11 biallelic loci located in non-recombining part of the chromosome (SRY1532, YAP, 92R7, DYF155S2, 12f2, Tat, M9, M17, M25, M89, and M56). In the Ukrainian gene, pool six haplogroups were revealed: E, F (including G and I), J, N3, P, and R1a1. These haplogroups were earlier detected in a study of Y-chromosome diversity on the territory of Europe as a whole. The major haplogroup in the Ukrainian gene pool, haplogroup R1a1 (earlier designated HG3), accounted for about 44% of all Y chromosomes in the sample examined. This haplogroup is thought to mark the migration patterns of the early Indo-Europeans and is associated with the distribution of the Kurgan archaeological culture. The second major haplogroup is haplogroup F (21.3%), which is a combination of the lineages differing by the time of appearance. Haplogroup P found with the frequency of 9.6%, represents the genetic contribution of the population originating from the ancient autochthonous population of Europe. Haplogroups J and E (11.7 and 4.2%, respectively) mark the migration patterns of the Middle-Eastern agriculturists during the Neolithic. The presence of the N3 lineage (9.6%) is likely explained by a contribution of the assimilated Finno-Ugric tribes. The data on the composition and frequencies of Y-chromosome haplogroups in the sample studied substantially supplement the existing picture of the male lineage distribution in the Eastern Slav population.     

Phylogeographic history of mitochondrial haplogroup J in Scandinavia

Background

Mitochondrial DNA haplogroup J is the third most frequent haplogroup in modern-day Scandinavia, although it did not originate there. To infer the genetic history of haplogroup J in Scandinavia, we examined worldwide mitogenome sequences using a maximum-likelihood phylogenetic approach.

Methods

Haplogroup J mitogenome sequences were gathered from GenBank (n = 2245) and aligned against the ancestral Reconstructed Sapiens Reference Sequence. We also analyzed haplogroup J Viking Age sequences from the European Nucleotide Archive (n = 54). Genetic distances were estimated from these data and projected onto a maximum likelihood rooted phylogenetic tree to analyze clustering and branching dates.

Results

Haplogroup J originated approximately 42.6 kya (95% CI: 30.0–64.7), with several of its earliest branches being found within the Arabian Peninsula and Northern Africa. J1b was found most frequently in the Near East and Arabian Peninsula, while J1c occurred most frequently in Europe. Based on phylogenetic dating, subhaplogroup J1c has its early roots in the Mediterranean and Western Balkans. Otherwise, the majority of the branches found in Scandinavia are younger than those seen elsewhere, indicating that haplogroup J dispersed relatively recently into Northern Europe, most plausibly with Neolithic farmers.

Conclusions

Haplogroup J appeared when Scandinavia was transitioning to agriculture over 6 kya, with J1c being the most common lineage there today. Changes in the distribution of haplogroup J mtDNAs were likely driven by the expansion of farming from West Asia into Southern Europe, followed by a later expansion into Scandinavia, with other J subhaplogroups appearing among Scandinavian groups as early as the Viking Age.

     

Y chromosomal haplogroup J as a signature of the post-neolithic colonization of Europe

In order to attain a finer reconstruction of the peopling of southern and central-eastern Europe from the Levant, we determined the frequencies of eight lineages internal to the Y chromosomal haplogroup J, defined by biallelic markers, in 22 population samples obtained with a fine-grained sampling scheme. Our results partially resolve a major multifurcation of lineages within the haplogroup. Analyses of molecular variance show that the area covered by haplogroup J dispersal is characterized by a significant degree of molecular radiation for unique event polymorphisms within the haplogroup, with a higher incidence of the most derived sub-haplogroups on the northern Mediterranean coast, from Turkey westward; here, J diversity is not simply a subset of that present in the area in which this haplogroup first originated. Dating estimates, based on simple tandem repeat loci (STR) diversity within each lineage, confirmed the presence of a major population structuring at the time of spread of haplogroup J in Europe and a punctuation in the peopling of this continent in the post-Neolithic, compatible with the expansion of the Greek world. We also present here, for the first time, a novel method for comparative dating of lineages, free of assumptions of STR mutation rates.     

Gene Pool Structure of Eastern Ukrainians as Inferred from the Y-Chromosome Haplogroups

Y chromosomes from representative sample of Eastern Ukrainians (94 individuals) were analyzed for composition and frequencies of haplogroups, defined by 11 biallelic loci located in non-recombining part of the chromosome (SRY1532, YAP, 92R7, DYF155S2, 12f2, Tat, M9, M17, M25,M89, andM56). In the Ukrainian gene, pool six haplogroups were revealed: E, F (including G and I), J, N3, P, and R1a1. These haplogroups were earlier detected in a study of Y-chromosome diversity on the territory of Europe as a whole. The major haplogroup in the Ukrainian gene pool, haplogroup R1a1 (earlier designated HG3), accounted for about 44% of all Y chromosomes in the sample examined. This haplogroup is thought to mark the migration patterns of the early Indo-Europeans and is associated with the distribution of the Kurgan archaeological culture. The second major haplogroup is haplogroup F (21.3%), which is a combination of the lineages differing by the time of appearance. Haplogroup P found with the frequency of 9.6%, represents the genetic contribution of the population originating from the ancient autochthonous population of Europe. Haplogroups J and E (11.7 and 4.2%, respectively) mark the migration patterns of the Middle-Eastern agriculturists during the Neolithic. The presence of the N3 lineage (9.6%) is likely explained by a contribution of the assimilated Finno–Ugric tribes. The data on the composition and frequencies of Y-chromosome haplogroups in the sample studied substantially supplement the existing picture of the male lineage distribution in the Eastern Slav population.     

Explaining the Imperfection of the Molecular Clock of Hominid Mitochondria

The molecular clock of mitochondrial DNA has been extensively used to date various genetic events. However, its substitution rate among humans appears to be higher than rates inferred from human-chimpanzee comparisons, limiting the potential of interspecies clock calibrations for intraspecific dating. It is not well understood how and why the substitution rate accelerates. We have analyzed a phylogenetic tree of 3057 publicly available human mitochondrial DNA coding region sequences for changes in the ratios of mutations belonging to different functional classes. The proportion of non-synonymous and RNA genes substitutions has reduced over hundreds of thousands of years. The highest mutation ratios corresponding to fast acceleration in the apparent substitution rate of the coding sequence have occurred after the end of the Last Ice Age. We recalibrate the molecular clock of human mtDNA as 7990 years per synonymous mutation over the mitochondrial genome. However, the distribution of substitutions at synonymous sites in human data significantly departs from a model assuming a single rate parameter and implies at least 3 different subclasses of sites. Neutral model with 3 synonymous substitution rates can explain most, if not all, of the apparent molecular clock difference between the intra- and interspecies levels. Our findings imply the sluggishness of purifying selection in removing the slightly deleterious mutations from the human as well as the Neandertal and chimpanzee populations. However, for humans, the weakness of purifying selection has been further exacerbated by the population expansions associated with the out-of Africa migration and the end of the Last Ice Age.     

A signal, from human mtDNA, of postglacial recolonization in Europe

Mitochondrial HVS-I sequences from 10,365 subjects belonging to 56 populations/geographical regions of western Eurasia and northern Africa were first surveyed for the presence of the T-->C transition at nucleotide position 16298, a mutation which has previously been shown to characterize haplogroup V mtDNAs. All mtDNAs with this mutation were then screened for a number of diagnostic RFLP sites, revealing two major subsets of mtDNAs. One is haplogroup V proper, and the other has been termed "pre*V," since it predates V phylogenetically. The rather uncommon pre*V tends to be scattered throughout Europe (and northwestern Africa), whereas V attains two peaks of frequency: one situated in southwestern Europe and one in the Saami of northern Scandinavia. Geographical distributions and ages support the scenario that pre*V originated in Europe before the Last Glacial Maximum (LGM), whereas the more recently derived haplogroup V arose in a southwestern European refugium soon after the LGM. The arrival of V in eastern/central Europe, however, occurred much later, possibly with (post-)Neolithic contacts. The distribution of haplogroup V mtDNAs in modern European populations would thus, at least in part, reflect the pattern of postglacial human recolonization from that refugium, affecting even the Saami. Overall, the present study shows that the dissection of mtDNA variation into small and well-defined evolutionary units is an essential step in the identification of spatial frequency patterns. Mass screening of a few markers identified using complete mtDNA sequences promises to be an efficient strategy for inferring features of human prehistory.     

百合目叶绿体rbcL基因核苷酸替代速率检测

分析了百合目主要类群叶绿体中编码核酮糖1,5二磷酸羧化氧化酶大亚基rbcL基因的42条序列,使用RRTree相对速率检测方法,详细研究rbcL基因在百合目7科间同义替代速率和非同义替代速率的变化.相对速率检测显示:百合目内秋水仙科(Colchicaceae)的同义替代速率和非同义替代速率均最快,金梅草科(Campynemat-aceae)同义替代速率最慢,百合科(Liliaceae)的非同义替代速率最慢,但在百合目各科间,无论同义替代速率还是非同义替代速率差异均不显著.     

Estimate of within population incremental selection through branch imbalance in lineage trees

Incremental selection within a population, defined as limited fitness changes following mutation, is an important aspect of many evolutionary processes. Strongly advantageous or deleterious mutations are detected using the synonymous to non-synonymous mutations ratio. However, there are currently no precise methods to estimate incremental selection. We here provide for the first time such a detailed method and show its precision in multiple cases of micro-evolution. The proposed method is a novel mixed lineage tree/sequence based method to detect within population selection as defined by the effect of mutations on the average number of offspring. Specifically, we propose to measure the log of the ratio between the number of leaves in lineage trees branches following synonymous and non-synonymous mutations. The method requires a high enough number of sequences, and a large enough number of independent mutations. It assumes that all mutations are independent events. It does not require of a baseline model and is practically not affected by sampling biases. We show the method''s wide applicability by testing it on multiple cases of micro-evolution. We show that it can detect genes and inter-genic regions using the selection rate and detect selection pressures in viral proteins and in the immune response to pathogens.     

Population history of the Red Sea--genetic exchanges between the Arabian Peninsula and East Africa signaled in the mitochondrial DNA HV1 haplogroup

Archaeological studies have revealed cultural connections between the two sides of the Red Sea dating to prehistory. The issue has still not been properly addressed, however, by archaeogenetics. We focus our attention here on the mitochondrial haplogroup HV1 that is present in both the Arabian Peninsula and East Africa. The internal variation of 38 complete mitochondrial DNA sequences (20 of them presented here for the first time) affiliated into this haplogroup testify to its emergence during the late glacial maximum, most probably in the Near East, with subsequent dispersion via population expansions when climatic conditions improved. Detailed phylogeography of HV1 sequences shows that more recent demographic upheavals likely contributed to their spread from West Arabia to East Africa, a finding concordant with archaeological records suggesting intensive maritime trade in the Red Sea from the sixth millennium BC onwards. Closer genetic exchanges are apparent between the Horn of Africa and Yemen, while Egyptian HV1 haplotypes seem to be more similar to the Near Eastern ones.     

藏鸡主要组织相容性复合体B-LBⅡ基因第二外显子序列遗传变异分析

根据鸡主要组织相容性复合体B-LBⅡ基因序列设计特异性引物,在藏鸡基因组中扩增了一个包括其第二外显子和第二内含子在内长度为374 bp的片段,并通过克隆和PCR直接测序获得了该片段的核苷酸序列。发现了15个B-LBⅡ新等位基因。对18个B-LBⅡ等位基因核苷酸序列和其所编码的MHCB-LBⅡ分子β1结构域的氨基酸序列分析显示,第二外显子核苷酸序列遗传多态性异常丰富,存在着62个多态变异位点(共包括80个突变),其中41个为简约性多态位点;衡量该序列遗传多样性的π值为0.0718;反映其群体内遗传变异度的平均遗传距离为0.056±0.008,低于在5个外来品种所估算的平均遗传距离。该编码区核苷酸相对异义替换率(15.61±2.69%)显著高于其同义替换率(3.25±0.94%),进一步分析表明,基因重组和平衡选择机制可能是引起B-LBⅡ基因序列变异的主要因素。在β1结构域氨基酸序列中,存在11个同义替换和27个异义替换;在24个肽结合位点中有12个变异位点;与其他6个中国地方鸡品种和一个外来品种比较发现,有11个异义氨基酸替换仅出现在藏鸡群体中,并被认为与藏鸡的免疫特异性有关,可为鸡的抗病力研究提供分子依据。     

Excavating Y-chromosome haplotype strata in Anatolia

Analysis of 89 biallelic polymorphisms in 523 Turkish Y chromosomes revealed 52 distinct haplotypes with considerable haplogroup substructure, as exemplified by their respective levels of accumulated diversity at ten short tandem repeat (STR) loci. The major components (haplogroups E3b, G, J, I, L, N, K2, and R1; 94.1%) are shared with European and neighboring Near Eastern populations and contrast with only a minor share of haplogroups related to Central Asian (C, Q and O; 3.4%), Indian (H, R2; 1.5%) and African (A, E3*, E3a; 1%) affinity. The expansion times for 20 haplogroup assemblages was estimated from associated STR diversity. This comprehensive characterization of Y-chromosome heritage addresses many multifaceted aspects of Anatolian prehistory, including: (1) the most frequent haplogroup, J, splits into two sub-clades, one of which (J2) shows decreasing variances with increasing latitude, compatible with a northward expansion; (2) haplogroups G1 and L show affinities with south Caucasus populations in their geographic distribution as well as STR motifs; (3) frequency of haplogroup I, which originated in Europe, declines with increasing longitude, indicating gene flow arriving from Europe; (4) conversely, haplogroup G2 radiates towards Europe; (5) haplogroup E3b3 displays a latitudinal correlation with decreasing frequency northward; (6) haplogroup R1b3 emanates from Turkey towards Southeast Europe and Caucasia and; (7) high resolution SNP analysis provides evidence of a detectable yet weak signal (<9%) of recent paternal gene flow from Central Asia. The variety of Turkish haplotypes is witness to Turkey being both an important source and recipient of gene flow.     

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.     

Y-chromosome lineages in Cabo Verde Islands witness the diverse geographic origin of its first male settlers

The Y-chromosome haplogroup composition of the population of the Cabo Verde Archipelago was profiled by using 32 single-nucleotide polymorphism markers and compared with potential source populations from Iberia, west Africa, and the Middle East. According to the traditional view, the major proportion of the founding population of Cabo Verde was of west African ancestry with the addition of a minor fraction of male colonizers from Europe. Unexpectedly, more than half of the paternal lineages (53.5%) of Cabo Verdeans clustered in haplogroups I, J, K, and R1, which are characteristic of populations of Europe and the Middle East, while being absent in the probable west African source population of Guiné-Bissau. Moreover, a high frequency of J* lineages in Cabo Verdeans relates them more closely to populations of the Middle East and probably provides the first genetic evidence of the legacy of the Jews. In addition, the considerable proportion (20.5%) of E3b(xM81) lineages indicates a possible gene flow from the Middle East or northeast Africa, which, at least partly, could be ascribed to the Sephardic Jews. In contrast to the predominance of west African mitochondrial DNA haplotypes in their maternal gene pool, the major west African Y-chromosome lineage E3a was observed only at a frequency of 15.9%. Overall, these results indicate that gene flow from multiple sources and various sex-specific patterns have been important in the formation of the genomic diversity in the Cabo Verde islands.An erratum to this article can be found at     

Historical isolation versus recent long-distance connections between Europe and Africa in bifid toadflaxes (Linaria sect. Versicolores)

Background

Due to its complex, dynamic and well-known paleogeography, the Mediterranean region provides an ideal framework to study the colonization history of plant lineages. The genus Linaria has its diversity centre in the Mediterranean region, both in Europe and Africa. The last land connection between both continental plates occurred during the Messinian Salinity Crisis, in the late Miocene (5.96 to 5.33 Ma).

Methodology/Principal Findings

We analyzed the colonization history of Linaria sect. Versicolores (bifid toadflaxes), which includes c. 22 species distributed across the Mediterranean, including Europe and Africa. Two cpDNA regions (rpl32-trnL^UAG and trnK-matK) were sequenced from 66 samples of Linaria. We conducted phylogenetic, dating, biogeographic and phylogeographic analyses to reconstruct colonization patterns in space and time. Four major clades were found: two of them exclusively contain Iberian samples, while the other two include northern African samples together with some European samples. The bifid toadflaxes have been split in African and European clades since the late Miocene, and most lineage and speciation differentiation occurred during the Pliocene and Quaternary. We have strongly inferred four events of post-Messinian colonization following long-distance dispersal from northern Africa to the Iberian Peninsula, Sicily and Greece.

Conclusions/Significance

The current distribution of Linaria sect. Versicolores lineages is explained by both ancient isolation between African and European populations and recent events of long-distance dispersal over sea barriers. This result provides new evidence for the biogeographic complexity of the Mediterranean region.     

Discovery of Western European r1b1a2 y chromosome variants in 1000 genomes project data: an online community approach

The authors have used an online community approach, and tools that were readily available via the Internet, to discover genealogically and therefore phylogenetically relevant Y-chromosome polymorphisms within core haplogroup R1b1a2-L11/S127 (rs9786076). Presented here is the analysis of 135 unrelated L11 derived samples from the 1000 Genomes Project. We were able to discover new variants and build a much more complex phylogenetic relationship for L11 sub-clades. Many of the variants were further validated using PCR amplification and Sanger sequencing. The identification of these new variants will help further the understanding of population history including patrilineal migrations in Western and Central Europe where R1b1a2 is the most frequent haplogroup. The fine-grained phylogenetic tree we present here will also help to refine historical genetic dating studies. Our findings demonstrate the power of citizen science for analysis of whole genome sequence data.     

African female heritage in Iberia: a reassessment of mtDNA lineage distribution in present times

The Iberian peninsula is a peripheral region of Europe in close proximity to Africa. Its inhabitants have an overall mtDNA genetic landscape typical of European background, although with signs of some African influence, whose features we deemed to disclose by analyzing available mtDNA HVRI distributions and new data. We analyzed 1,045 sequences. The most relevant results are the following: (1) North African sequences (haplogroup U6) present an overall frequency of 2.39%, and sub-Saharan sequences reach 3.83%, values that are, in both cases, much higher than those generally observed in Europe; and (2) there is a substantial geographic heterogeneity in the distribution of these lineages (haplogroup L being the most frequent in the south, whereas haplogroup U6 is generally more common in the north). The analysis of the observed diversity within each haplogroup strongly suggests that both were recently introduced (in historical times). Although for haplogroup U6 the documented event that is demographically compatible is the Islamic period (beginning of the 8th century to the end of the 15th century), for haplogroup L the most probable origin is the modern slave trade (mid 15th century to the end of the 18th century). However, the observed geographic structuring for one of the haplogroups does not fit the expected distribution provided by simplistic historical considerations. In fact, although for haplogroup L the north-south increasing frequency is corroborated by historical data, the opposite trend, observed for haplogroup U6, is more difficult to reconcile with the magnitude and time span of the Islamic political and cultural influence, which lasted longer and was more intense in the south. To clarify this conundrum, we need not only a substantial increase in the amount of mtDNA data (particularly for North Africa) but also new historical data and interpretations.     

The role of mtDNA background in disease expression: a new primary LHON mutation associated with Western Eurasian haplogroup J

Leber's hereditary optic neuropathy (LHON) is a maternally transmitted form of blindness caused by mitochondrial DNA (mtDNA) mutations. Approximately 90% of LHON cases are caused by 3460A, 11778A, or 14484C mtDNA mutations. These are designated "primary" mutations because they impart a high risk for LHON expression. Although the 11778A and 14484C mutations unequivocally predispose carriers to LHON, they are preferentially associated with mtDNA haplogroup J, one of nine Western Eurasian mtDNA lineages, suggesting a synergistic and deleterious interaction between these LHON mutations and haplogroup J polymorphism(s). We report here the characterization of a new primary LHON mutation in the mtDNA ND4L gene at nucleotide pair 10663. The homoplasmic 10663C mutation has been found in three independent LHON patients who lack a known primary mutation and all of which belong to haplogroup J. This mutation has not been found in a large number of haplotype-matched or non-haplogroup-J control mtDNAs. Phylogenetic analysis with primarily complete mtDNA sequence data demonstrates that the 10663C mutation has arisen at least three independent times in haplogroup J, indicating that it is not a rare lineage-specific polymorphism. Analysis of complex I function in patient lymphoblasts and transmitochondrial cybrids has revealed a partial complex I defect similar in magnitude to the 14484C mutation. Thus, the 10663C mutation appears to be a new primary LHON mutation that is pathogenic when co-occurring with haplogroup J. These results strongly support a role for haplogroup J in the expression of certain LHON mutations.     

Clustering of Caucasian Leber hereditary optic neuropathy patients containing the 11778 or 14484 mutations on an mtDNA lineage.

Leber hereditary optic neuropathy (LHON) is a type of blindness caused by mtDNA mutations. Three LHON mtDNA mutations at nucleotide positions 3460, 11778, and 14484 are specific for LHON and account for 90% of worldwide cases and are thus designated as "primary" LHON mutations. Fifteen other "secondary" LHON mtDNA mutations have been identified, but their pathogenicity is unclear. mtDNA haplotype and phylogenetic analysis of the primary LHON mutations in North American Caucasian patients and controls has shown that, unlike the 3460 and 11778 mutations, which are distributed throughout the European-derived (Caucasian) mtDNA phylogeny, patients containing the 14484 mutation tended to be associated with European mtDNA haplotype J. To investigate this apparent clustering, we performed chi2-based statistical analyses to compare the distribution of LHON patients on the Caucasian phylogenetic tree. Our results indicate that, unlike the 3460 and 11778 mutations, the 14484 mutation was not distributed on the phylogeny in proportion to the frequencies of the major Caucasian mtDNA haplogroups found in North America. The 14484 mutation was next shown to occur on the haplogroup J background more frequently that expected, consistent with the observation that approximately 75% of worldwide 14484-positive LHON patients occur in association with haplogroup J. The 11778 mutation also exhibited a moderate clustering on haplogroup J. These observations were supported by statistical analysis using all available mutation frequencies reported in the literature. This paper thus illustrates the potential importance of genetic background in certain mtDNA-based diseases, speculates on a pathogenic role for a subset of LHON secondary mutations and their interaction with primary mutations, and provides support for a polygenic model for LHON expression in some cases.     

Positive selection of synonymous mutations in vesicular stomatitis virus

Prevailing evolutionary forces are typically deduced from the pattern of differences in synonymous and non-synonymous mutations, under the assumption of neutrality in the absence of amino acid change. We determined the complete sequence of ten vesicular stomatitis virus populations evolving under positive selection. A significant number of the mutations occurred independently in two or more strains, a process known as parallel evolution, and a substantial fraction of the parallel mutations were silent. Parallel evolution was also identified in non-coding regions. These results indicate that silent mutations can significantly contribute to adaptation in RNA viruses, and relative frequencies of synonymous and non-synonymous substitutions may not be useful to resolve their evolutionary history.     

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.