The main angle of the recently released epic manuscript Haak et al. 2015 is that ancient DNA supports the steppe origin of at least some of Europe's Indo-European languages. That's certainly a move in the right direction, so that we can eventually do away with the Anatolian hypothesis, which was always a failed proposition.
But it's clear that the authors are holding back. They've obviously decided to be very cautious until they've looked at more ancient DNA, particularly from the Near East, Central Asia and India, before backing fully any one Proto-Indo-European (PIE) urheimat model.
That's understandable, considering how much opposition there is still to the steppe hypothesis, even though it does by and large have the support of historical linguists, which is what really counts. Nevertheless, my feeling is that Haak et al. are underselling their data, particularly the stuff from Eastern Europe.
I'm of the opinion that the steppe or Kurgan PIE model works just fine, and also not surprised by the ancient DNA evidence pointing to a massive expansion of people from the western steppe during the Late Neolithic/Early Bronze Age. So for me, the really big news in this paper is that the only two Eastern European forager samples belong to basal lineages of Y-chromosome haplogroups R1a and R1b. What this suggests, Id' say, is that ancient Eastern Europe was a key bifurcation region for R1.
Remarkably, it's possible to basically lay out the history and phylogeny of R1a in Europe using just three R1a samples from the paper. This can't be a coincidence.
- Mesolithic Hunter-Gatherer from Karelia: R1a (xM198)
- Late Neolithic Corded Ware pastoralist from Germany: R1a (M198, M417, xZ282)
- Late Bronze Age Urnfielder from Germany: R1a (M198, M417, Z282, Z280)
What we can see there is the progression from a basal R1a in pre-Neolithic Northeastern Europe to a derived R1a in late prehistoric Central Europe. The derived R1a is actually R1a1a1b1a2, which is by far the most common subclade of R1a in Europe today, and closely related to the Asian and Indo-Iranian-specific R1a1a1b2.
Interestingly, all seven of the Yamnaya males sampled by Haak et al., mostly from the Samara Valley, belong to R1b-M269, the most common subclade of R1b today. However, five belong to the West Asian-specific R1b-Z1203, but none to the West European-specific R1b-M412. Also, all nine Yamnaya samples show Near Eastern admixture, described in the paper as Armenian-like.
Does this perhaps mean that the Proto-Indo-Europeans (and thus Yamnaya) originated in the Near East, as per the Armenian Plateau hypothesis?
I doubt it. The aforementioned Eastern European R1b forager is also from the Samara Valley, and he clearly lacks Near Eastern admixture. So what are the chances that a Near Eastern population with a frequency of R1b-M269 of around 100% moved into an area of Eastern Europe where a more basal R1b was already present, and in fact in a population with no Near Eastern ancestry? Very slim, I'd say.
So how did the Yamnaya herders acquire their Near Eastern admixture? The answer is obvious if we look at their mtDNA haplogroups. These include H, T and W, all of which might have come to Eastern Europe from the Near East.
Of course this doesn't mean that the Eastern European steppe was overrun by Near Eastern Amazons. It's generally accepted that during the Neolithic the steppe was settled by farmers from the Near East, just like much of the rest of Europe, and I'd say that it was mostly the women from these groups who were incorporated into the later pastoralist societies of the steppe. The men, who probably belonged to Near Eastern haplogroups like G or T, might have been killed or marginalized in some way, so that their reproductive success was seriously hampered.
This is not a far fetched scenario. Typical hunter-gatherer Y-haplogroups like I2 and C6 have already been recorded alongside Near Eastern-specific mtDNA lineages at several Neolithic sites in Western and Central Europe. The social mechanisms for this might have been different there than on the steppe, but in any case, it seems that European hunter-gatherer males shacking up with farm girls of largely Near Eastern ancestry was not an unusual occurrence back in the day.
Now, if Eastern Europe was indeed a bifurcation hotspot for R1, then a large proportion, or even the majority of R1a and R1b in Eurasia today, might well be of Eastern European origin. If so, there should be some support for this in genome-wide DNA of present-day Asians, and indeed I think there is.
Below are a couple of principal component analyses (PCA). The first is from Haak et al. and the second from my own West Eurasia K8 analysis (see here). Unfortunately, I don't yet have access to the Yamnaya genomes, but I think it's petty easy to guesstimate where they will land on my plot when I run them in the K8. I marked this spot with an X.
Note that most of the Near Eastern and Caucasian populations are clearly shifted east towards ANE, and also up towards Europe. Moreover, I'd say many of these groups are specifically pushing up towards the Volga-Ural samples and thus the Yamnaya herders.
There's really no other way to explain this outcome. Quite simply, the vast majority of West Asians have relatively recent (post-Neolithic?) ancestry from the Ural or Kazakh steppe, which manifests itself as a west to east cline on PCA, running from the southern Levant to the north Caucasus. This result is easily reproduced on any decent PCA with West Eurasian populations, and can be seen on the Haak et al. plot.
I'm yet to find solid evidence that Indo-European speakers from the Near East, like Armenians, Kurds and Iranians, don't harbor fairly significant ancestry from this northeastern source.
For instance, unlike many people, I don't find unsupervised ADMIXTURE analyses very convincing when they show these groups to be entirely of Near Eastern ancestry. That's because when ADMIXTURE creates a modern Near Eastern/West Asian cluster, it usually lumps within it all of the ancient ancestral components that are today ubiquitous in the Near East. In other words, the steppe admixture which shows up amongst most West Asians on the PCA above is classified as native to the Near East, even though this is unlikely to be true.
See also...
High female mobility in Bronze Age Europe
Ust'-Ishim belongs to K-M526
I'll probably end up writing a whole series of posts on this paper. But for now, here's the abstract and a PCA.
We generated genome-wide data from 69 Europeans who lived between 8,000-3,000 years ago by enriching ancient DNA libraries for a target set of almost four hundred thousand polymorphisms. Enrichment of these positions decreases the sequencing required for genome-wide ancient DNA analysis by a median of around 250-fold, allowing us to study an order of magnitude more individuals than previous studies and to obtain new insights about the past. We show that the populations of western and far eastern Europe followed opposite trajectories between 8,000-5,000 years ago. At the beginning of the Neolithic period in Europe, ~8,000-7,000 years ago, closely related groups of early farmers appeared in Germany, Hungary, and Spain, different from indigenous hunter-gatherers, whereas Russia was inhabited by a distinctive population of hunter-gatherers with high affinity to a ~24,000 year old Siberian6. By ~6,000-5,000 years ago, a resurgence of hunter-gatherer ancestry had occurred throughout much of Europe, but in Russia, the Yamnaya steppe herders of this time were descended not only from the preceding eastern European hunter-gatherers, but from a population of Near Eastern ancestry. Western and Eastern Europe came into contact ~4,500 years ago, as the Late Neolithic Corded Ware people from Germany traced ~3/4 of their ancestry to the Yamnaya, documenting a massive migration into the heartland of Europe from its eastern periphery. This steppe ancestry persisted in all sampled central Europeans until at least ~3,000 years ago, and is ubiquitous in present-day Europeans. These results provide support for the theory of a steppe origin of at least some of the Indo-European languages of Europe.
Haak et al., Massive migration from the steppe is a source for Indo-European languages in Europe, bioRxiv, Posted February 10, 2015, doi: https://dx.doi.org/10.1101/013433
Here's the latest teaser for the new David Reich et al. paper on the ethnogenesis of present-day Europeans. It's part of an abstract for a seminar to be held by Professor Reich at Jesus College, Oxford, on February 9. Interestingly, it argues that migrations from the steppe resulted in a ~50% population turnover across northern Europe from the late Neolithic onwards, which is very much in agreement with recent discussions on the topic at Eurogenes (for instance, see here).
By ~6,000-5,000 years ago, a resurgence of hunter-gatherer ancestry had occurred throughout much of Europe, but in Russia, the Yamnaya steppe herders of this time were descended not only from the preceding eastern European hunter-gatherers, but also from a population of Near Eastern ancestry. Western and Eastern Europe came into contact ~4,500 years ago, as the Late Neolithic Corded Ware people from Germany traced ~3/4 of their ancestry to the Yamnaya, documenting a massive migration into the heartland of Europe from its eastern periphery. This steppe ancestry persisted in all sampled central Europeans until at least ~3,000 years ago, and comprises about half the ancestry of today’s northern Europeans. These results support the theory of a steppe origin of at least some of the Indo-European languages of Europe, and show the power of genome-wide ancient DNA studies to document human migrations.
Source: Ancient DNA documents three ancestral populations for present-day Europeans
Update 11/02/2015: Massive migration from the steppe is a source for Indo-European languages in Europe (Haak et al. 2015 preprint).
Haak et al., Massive migration from the steppe is a source for Indo-European languages in Europe, bioRxiv, Posted February 10, 2015, doi: https://dx.doi.org/10.1101/013433
I'm posting a new entry about the upcoming Corded Ware/Yamnaya paper because the last entry (see here) now has over 400 comments which aren't easy to load for many people.
One of the authors of this eagerly awaited paper, Nick Patterson of the Broad Institute, briefly joined our discussion. Nick's contribution is much appreciated. He wasn't able to reveal a great deal, because the manuscript is in submission, but he did make a couple of interesting points:
- the paper will feature Y-haplogroup results from the Yamnaya culture, represented by nine samples in all, including seven males
- the population with Near Eastern ancestry that mixed with the Eastern Hunter-Gatherers (EHG) on the Russian steppe to form the Yamnaya pastoralists by 5,000 YBP was also "rich" in ANE
- ancient DNA from the Caucasus, Iran and India is probably necessary to work out how the Indo-Europeans got to India, but the paper won't feature such data
It's nice to hear that Y-haplogroups aren't being ignored. My opinion is that they're at least as important as genome-wide data when tracking the movements across vast space and time of highly patriarchal and patrilineal groups like the ancient Indo-Europeans.
Indeed, we already know that the Slavic, Baltic and Norse-specific R1a1a1b1, defined by the Z282 mutation, is the sister clade of the Indo-Iranian-specific R1a1a1b2, defined by Z93. Thus, if the Yamnaya males were found to belong to these or upstream markers, this would suggest that they were the paternal ancestors of many Balts, Scandinavians, Slavs and Indo-Iranians, and correlate very nicely with the linguistic and archeological "steppe hypothesis" of Indo-European origins.
In fact, even if analyses based on high density genome-wide data suggest that Indians don't harbor any genome-wide European ancestry, we'd still have to accept the likelihood of gene flow - albeit perhaps very indirect gene flow - from the European steppe to India because many Indians belong to R1a1a1b2.
The second point made by Nick is perhaps surprising, but at least for me not totally unexpected. That's because we've already known for a while that the Yamnaya genomes can be successfully modeled as half Karelian EHG and half present-day Armenian (see here), and according to my own estimates Armenians carry an average of 15.5% ANE.
The fact that these Armenian-like, ANE-rich newcomers dampened the genome-wide affinity to ANE-proxy MA-1 on the Russian steppe might look like a contradiction, but not if we remember that the higher the Near Eastern ancestry the lower the genome-wide affinity to MA-1, and also consider that the steppe foragers probably carried a lot more ANE than the newcomers.
Actually, as far as I know, all of the Yamnaya samples in this study come from the Samara Valley, which is some distance north of the Caspian Sea near the southern Urals. So it makes senses that the pseudo Armenians who turned up there more than 5,000 years ago were not like the Neolithic farmers of Western and Central Europe, who lacked ANE.
I'd say that this as yet unidentified group (wild guess: immediate ancestors of the Repin culture people?) was the result of an admixture event, or perhaps a series of admixture events, with ANE-rich foragers somewhere on the steppe south of the Samara. If so, I won't be surprised if it turns out that R1a only appeared in the Samara Valley after their arrival.
In any case, it looks like even after this paper comes out, we'll still need a lot more ancient DNA from across Eurasia to help map out the early Indo-European dispersals with any confidence.
Update 11/02/2015: Massive migration from the steppe is a source for Indo-European languages in Europe (Haak et al. 2015 preprint) .
This is yet another teaser for the upcoming Corded Ware/Yamnaya paper from the Reich lab. Sadly, it doesn't mention Y-chromosome haplogroups, so perhaps the authors are going to tackle this issue later. However, check out what they say about the German and Spanish farmers being of the same stock, and the resurgence of hunter-gatherer ancestry in Western Europe after the early Neolithic. Fascinating stuff.
Ancient DNA points to the Eurasian steppe as a proximate source for Indo-European migrations into Europe
David Reich and Nick Patterson
Abstract: We generated genome-wide data from 65 Europeans who lived between 8,000-3,000 years ago by enriching ancient DNA libraries for a target set of about 390,000 single nucleotide polymorphisms. This strategy decreases the sequencing required to obtain genome-wide data from ancient DNA samples by around 1000-fold, allowing us to study an order of magnitude more individuals than previous studies and to obtain new insights about the past. We show that in western Europe, the farmers of both Germany and Spain >7,000 years ago were descended from a common ancestral stock. These farmers did not replace the earlier hunter-gatherers, but continued to mix with them, leading to a resurgence of hunter-gatherer ancestry in both Germany and Spain ~1,000-2,000 years later. In eastern Europe, the hunter-gatherers of Russia >7,000 years ago were distinct from those of the west, having an increased affinity to a ~24,000 year old individual from Siberia, but this affinity was reduced by ~5,000 years ago in the Yamnaya steppe pastoralists because of admixture with a population of Near Eastern ancestry. Western and Eastern Europe collided ~4,500 years ago with the appearance of the Corded Ware people in Central Europe, who derived at least two thirds of their ancestry from an eastern population closely related to the Yamnaya. The evidence for mass migration into Europe thousands of years after the arrival of agriculture, in combination with linguistic and archaeological data, makes a compelling case for the steppe as a proximate source for the spread of Indo-European languages into Europe.
Source: INA Kolloquium Ws 2014/15
Update 11/02/2015: Massive migration from the steppe is a source for Indo-European languages in Europe (Haak et al. 2015 preprint) .
Not long ago Lazaridis et al. proposed that most present-day Europeans were derived from three distinct ancestral populations: Ancient North Eurasians (ANE), Early European Farmers (EEF) and Western European Hunter-Gatherers (WHG).
However, this is essentially a stop-gap model, which will in all likelihood be replaced by a partly revised and more robust model once someone manages to sequence a genome or two from the Neolithic Near East. That's because EEF is clearly a hybrid component, largely made up of ancient Near Eastern ancestry and something very WHG-like, sometimes in very different proportions depending on the location and archeological context of the EEF genomes being analyzed.
So what will this new model look like, you might ask? Probably like this, where EEF is replaced by an Early Neolithic Farmer (ENF) component from the ancient Near East, or something very similar:
The diagram above is basically a Principal Component Analysis (PCA) based on output from my new West Eurasia K8 test (see here), in which the Near Eastern component is synonymous with ENF.
I'm quite certain that these results are very close to the truth. However, just in case the Near Eastern ancestry proportions are a little bit too high (and we won't know until we see those ancient genomes from the Near East), I've got another version that offers lower bound Near Eastern estimates.
It might be useful to keep in mind that I rotated the plots to fit geography. As a result, Component 1, which packs around 85% of the variance on both plots, appears smaller than Component 2, which only carries around 10% of the variance.
A spreadsheet with West Eurasia K8 results for a wide variety of populations is available here. Please note that there are two sheets, with the second sheet showing the lower bound Near Eastern ancestry proportions.
We'll probably learn of more ancient European meta-populations as many more genomes are sequenced from across Eurasia. Nevertheless, I doubt this will affect the model outlined above. That's because I'm expecting all such meta-populations to be mixtures of ANE, ENF and/or WHG, as well as, in some cases, extra-West Eurasian components.
However, I suspect that West Eurasia will have to be modeled in a different way from Europe, with, amongst other things, the so called Basal Eurasian component replacing ENF. But for this to happen we'll need at least one ancient genome that is in large-part of Basal Eurasian origin. In any case, that's a whole different subject.
See also...
4mix: four-way mixture modeling in R
I've seen quite a few comments on this blog suggesting that most of the Ancient North Eurasian (ANE) admixture found in Northern Europe today might come from Scandinavian hunter-gatherers like Motala12 and Ajvide58. It's probably obvious to most that this is not realistic, because the Scandinavian forager genomes sequenced to date show very high ratios of Western European Hunter-Gatherer (WHG) ancestry (>80%), so basically the math doesn't add up.
Nevertheless, I thought it might be useful to drive the point home using this Principal Component Analysis (PCA) based on my new West Eurasia K8 test. The datasheet is available here. You can view a spreadsheet of the results with extra samples here.
Please note that neither Motala12 nor Gokhem2, a late Neolithic farmer from south Sweden belonging to the Funnelbeaker culture, can pass for present-day Swedes. Moreover, mixing Gokhem2 with Motala12, in any proportions, will not produce a result even vaguely similar to present-day Swedes (ie. the outcome will fall somewhere along the dotted line).
I'd say one of the most obvious ways to get the right result would be to blend the Scandinavian forager and farmer with at least one other sample from somewhere below (ie. geographically speaking, east or southeast) of the Swedish cluster.
It might be possible to come up with a more precise plot location, and thus perhaps geographic origin, for this putative third source of Swedish ancestry by running some complex tests with the PCA datasheet. If anyone wants to have a go at that, and you actually manage to come up with a coherent outcome, then feel free to post your findings in the comments below.
I've decided not to bother, because as far as I can see, the options are infinite. What we really need are more genomes from the Swedish late Neolithic/early Bronze Age (LN/EBA), preferably belonging to one of the local spin-offs of the Corded Ware culture, which is thought to have originated in Eastern Europe, to provide more datapoints and help narrow down the options.
On a related note, I'm catching up on some reading this holiday season, and currently going through this book chapter which discusses the upheavals during the LN/EBA in south Scandinavia as seen through its archeology.
Rune Iversen, Beyond the Neolithic transition - the "de-Neolithisation" of south Scandinavia
See also...
Bell Beaker, Corded Ware, EHG and Yamnaya genomes in the fateful triangle
Here's a Principal Component Analysis (PCA) and an accompanying biplot based on output from an improved version of my ANE K7 ancestry test. Let's call it the West Eurasia K8. This one gives more accurate estimates of Western European Hunter-Gatherer (WHG) and Near Eastern admixture proportions, thanks to the use of new ancient samples.
When rotated accordingly (like here), the results are basically indistinguishable from those I get with genotype data (for instance, see here and here), which suggests that they're correct and based on ancestry proportions that are close to the truth. The Past3 data sheet used to create the PCA is available here. You can view a spreadsheet of the results with extra samples here.
Clearly, ANE is the main agent causing the west to east differentiation in dimension 2. Note that even a small rise in ANE, say, 4-5%, creates significant distance between samples on the PCA plot.
East and South Eurasian admixture has a similar effect, but must be more considerable to make an impact on a West Eurasian-specific PCA like this (and it does with the obvious Volga-Ural outliers, who come from Chuvashia and Tatarstan).
On the other hand, Near Eastern admixture without ANE creates almost the opposite effect. Note, for instance, that Neolithic genomes Stuttgart and NE1 show much higher levels of Near Eastern ancestry than most Europeans, and yet they're amongst the most western samples on the plot.
This suggests that the Near East, and in particular the Caucasus, experienced a significant rush of ANE admixture after early Neolithic farmers left the region for Europe. Alternatively, Caucasus populations may have carried even higher levels of ANE than they do today, before newcomers from the Near East mixed with them. But either way, a lot of ANE arrived in the Near East at some point.
It also suggests that, overall, the populations that moved west across northern Europe after the Neolithic, and shifted northern European genetic structure to the east, did not carry high ratios of Near Eastern ancestry. Instead, they harbored high ratios of ANE and WHG. What these ratios were exactly I haven't a clue, but ancient DNA should tell us that soon.
Below are the ancestry proportions for the five ancient genomes in this analysis, in chronological order. It's interesting to note (yet again) the rising and falling Near Eastern admixture, from the Mesolithic to Neolithic and then from the Neolithic to Bronze Age, respectively, as well as the steady rise of ANE from the Bronze Age to the Iron Age.
Loschbour (Mesolithic)
ANE 0
South_Eurasian 0
Near_Eastern 0
East_Eurasian 0
WHG 99.5
Oceanian 0.5
Pygmy 0
Sub-Saharan 0
Stuttgart (Neolithic)
ANE 0
South_Eurasian 0
Near_Eastern 72.19
East_Eurasian 0
WHG 27.8
Oceanian 0
Pygmy 0
Sub-Saharan 0
NE1 (Neolithic)
ANE 0
South_Eurasian 0
Near_Eastern 69.82
East_Eurasian 0
WHG 30.17
Oceanian 0
Pygmy 0
Sub-Saharan 0
BR2 (Bronze Age)
ANE 9.62
South_Eurasian 0.08
Near_Eastern 43.96
East_Eurasian 0
WHG 45.44
Oceanian 0.48
Pygmy 0.23
Sub-Saharan 0.19
Hinxton4 (Iron Age)
ANE 15.08
South_Eurasian 0.06
Near_Eastern 35.44
East_Eurasian 0.46
WHG 48.5
Oceanian 0
Pygmy 0
Sub-Saharan 0.46
See also...
The fateful triangle
Bell Beaker, Corded Ware, EHG and Yamnaya genomes in the fateful triangle
The map below is based on data from Warinner et al. 2014. It shows the consumption of milk, or lack of, among Late Neolithic/Bronze Age (LN/BA) individuals from across West Eurasia. Admittedly, the sampling is very sparse, but like I've said before on these blogs, the LN/BA was a time of profound changes in Europe, so every scrap of data from this period is very valuable.
Note the lack of milk consumption among the samples from north of the Alps, where today the vast majority of people consume milk as adults, and can do so because they carry the Lactase Persistence Allele (T-13910). This doesn't look like a coincidence, considering the mounting evidence of a major population turnover across much of Europe during the LN/BA, mostly as a result of migrations from the east.
Citation...
Warinner, C. et al. Direct evidence of milk consumption from ancient human dental calculus. Sci. Rep. 4, 7104; DOI:10.1038/srep07104 (2014).
See also...
Lactase persistence and ancient DNA
Ancient genomes from the Great Hungarian Plain
Update 20/05/2015: Large-scale recent expansion of European patrilineages
...
I wonder what the hardcore Y-DNA genetic genealogists will say about this effort? I know that many of those guys have been working with full Y-chromosome sequences for a while now. It's open access with lots of supplementary info.
Abstract: Many studies of human populations have used the male-specific region of the Y chromosome (MSY) as a marker, but MSY sequence variants have traditionally been subject to ascertainment bias. Also, dating of haplogroups has relied on Y-specific short tandem repeats (STRs), involving problems of mutation rate choice, and possible long-term mutation saturation. Next-generation sequencing can ascertain single nucleotide polymorphisms (SNPs) in an unbiased way, leading to phylogenies in which branch-lengths are proportional to time, and allowing the times-to-most-recent-common-ancestor (TMRCAs) of nodes to be estimated directly. Here we describe the sequencing of 3.7 Mb of MSY in each of 448 human males at a mean coverage of 51x, yielding 13,261 high-confidence SNPs, 65.9% of which are previously unreported. The resulting phylogeny covers the majority of the known clades, provides date estimates of nodes, and constitutes a robust evolutionary framework for analysing the history of other classes of mutation. Different clades within the tree show subtle but significant differences in branch lengths to the root. We also apply a set of 23 Y-STRs to the same samples, allowing SNP- and STR-based diversity and TMRCA estimates to be systematically compared. Ongoing purifying selection is suggested by our analysis of the phylogenetic distribution of non-synonymous variants in 15 MSY single-copy genes.
Here are a couple of interesting quotes. You can see the samples they're talking about on the tree below. As per the second paragraph, it seems there's a paper about to be published at Nature Communications on European Y-chromosome haplogroups based on some heavy resequencing data (see Batini et al. in the references list). Can't wait for that.
(viii) Rare deep-rooting hg Q lineages in NW Europe: Hg Q has been most widely investigated in terms of the peopling of the Americas from NE Asia (Karafet et al. 1999). Here, as well as an example of the common native American Q-M3 lineage, we included examples of rare European hg Q chromosomes. One of the English chromosomes belongs to the deepest-rooting lineage within Q (Q-M378) and may reflect the Jewish diaspora (Hammer et al. 2009); the other is distantly related, shares a deep node with the Mexican Q-M3 chromosome, and has an STR-haplotype closely related to those of scarce Scandinavian hg Q chromosomes (unpublished data).
(ix) Structure within the west Eurasian hg R: The TMRCA of hg R is 19 KYA, and within it both hgs R1a and R1b comprise young, star-like expansions discussed extensively elsewhere (Batini et al. submitted). The addition of Central Asian chromosomes here contributes a sequence to the deepest subclade of R1b-M269, while another, in a Bhutanese individual, forms an outgroup almost as old as the R1a/R1b split.
Citation...
Hallast et al., The Y-chromosome tree bursts into leaf: 13,000 high-confidence SNPs covering the majority of known clades, Molecular Biology & Evolution, published online December 2, 2014, doi: 10.1093/molbev/msu327
