On Maykop ancestry in Yamnaya

What Maykop ancestry in Yamnaya? There is none, or at least not enough worth discussing, except in one highly unusual female outlier from a burial in what is now eastern Ukraine. But apparently this is still up for debate? Well it shouldn't be.

To anyone with even a passing interest in the Yamnaya culture, it should be rather obvious that it formed during the tail end of the Eneolithic on the Pontic-Caspian steppe, as basically a direct offshoot of the earlier Repin culture, but perhaps also with significant influences from the earlier still Khvalynsk and Sredny Stog cultures. So why should its population history be much different from this?

It isn't, and this is fairly easy to demonstrate now despite the still rather poor sampling of Eneolithic remains from the Pontic-Caspian steppe.

Below is a series of qpAdm analyses in which I modeled several Yamnaya groups, as well as the closely related Afanasievo and Poltavka populations, exclusively and successfully as two- and three-way mixtures of a few Eneolithic singletons from various parts of the Pontic-Caspian steppe (obviously, I'd love to use homogeneous population sets instead, but, as per my point above, that's not possible yet). The models are sorted by their statistical fits, best to worst. Also note the large number and wide range of right pops or outgroups. I wanted to make sure that I wasn't missing anything.

Yamnaya_Samara
Dereivka_I_I4110 0.324±0.035
Progress_Eneolithic_PG2004 0.676±0.035
chisq 6.797
tail prob 0.976979
Full output

Afanasievo
Progress_Eneolithic_PG2004 0.638±0.038
Sredny_Stog_II_I6561 0.362±0.038
chisq 10.855
tail prob 0.818366
Full output

Yamnaya_Ukraine
Progress_Eneolithic_PG2001 0.655±0.073
Sredny_Stog_II_I6561 0.345±0.073
chisq 12.676
tail prob 0.696277
Full output

Poltavka
Dereivka_I_I4110 0.324±0.038
Progress_Eneolithic_PG2004 0.676±0.038
chisq 12.895
tail prob 0.680437
Full output

Yamnaya_Caucasus
Khvalynsk_Eneolithic_I0122 0.086±0.054
Sredny_Stog_II_I6561 0.221±0.070
Vonyuchka_Eneolithic_VJ1001 0.693±0.101
chisq 13.113
tail prob 0.593562
Full output

So, you might ask, is there any way to add Maykop to these models? Nope, it's pointless, because it doesn't improve the stats (for instance, see here, here and here). In other words, the situation is this: I already have awesome models, and I can't readily fit Maykop into my framework, so why do it? But if anyone out there wants to try, then by all means, and feel free to share the results with us in the comments.

Of course, the fact that most of these Yamnaya and Yamnaya-related populations are best modeled with somewhat different Eneolithic steppe singletons doesn't mean that they have radically different origins. In fact, they're all very closely related and they're basically like one Bronze Age steppe family. They just harbor somewhat different ratios of the same ancient ancestral components.

For the sake of being thorough, as per scientific literature, I pooled all of the above Afanasievo, Poltavka and Yamnaya samples into a Steppe_EMBA set and analyzed it with several genetically and geographically matching pairs of the Eneolithic singletons. This was one of the best fitting models, which I think is interesting, because the region roughly between the burial sites of these pairs of Eneolithic individuals was the home of the Repin culture.

Steppe_EMBA
North_Pontic_Eneolithic_I4110-I656 0.313±0.027
Progress_Eneolithic_PG2001-PG2004 0.687±0.027
chisq 15.378
tail prob 0.497157
Full output

Again, adding Maykop to this model makes no sense (see here, here and here). Clearly, I'd have to come up with a very different framework to successfully model Steppe_EMBA with a Maykop population. However, it's unlikely that such a model would make much sense in the context of various other types of genetic analyses and archeological data.

See also...

Yamnaya: home-grown

Big deal of 2018: Yamnaya not related to Maykop

Yamnaya isn't from Iran just like R1a isn't from India

Yamnaya: home-grown

I have some interesting news. It looks like Khvalynsk_Eneolithic I0434 can be used as essentially a perfect proxy for the Eneolithic steppe trio from Wang et al. 2018 when modeling the ancestry of the Yamnaya people of what is now the Samara region of Russia. Consider the qpAdm mixture models below, sorted by taildiff.

One of the best fitting models that also fairly closely matches archeological data, which suggest that Yamnaya was an amalgamation of the Khvalynsk, Repin and Sredny Stog cultures, is in bold. The worst fitting, and basically failed, models are listed below the dotted line. Note that almost all of these models feature reference populations from West and Central Asia.

Khvalynsk_I0434 + Iberia_ChL 0.681534184 > full output

Khvalynsk_I0434 + Globular_Amphora 0.525961242 > full output

Khvalynsk_I0434 + Iberia_Central_CA 0.515960444 > full output

Khvalynsk_I0434 + Sredny_Stog_I6561 0.485311962 > full output

Khvalynsk_I0434 + Varna 0.430411416 > full output

Khvalynsk_I0434 + Blatterhole_MN 0.328782809 > full output

Khvalynsk_I0434 + Baden_LCA 0.234307235 > full output

Khvalynsk_I0434 + Protoboleraz_LCA 0.231310724 > full output

Khvalynsk_I0434 + ALPc_MN 0.200002422 > full output

Khvalynsk_I0434 + Trypillia 0.193900977 > full output

Khvalynsk_I0434 + Balaton_Lasinja_CA 0.187031564 > full output

Khvalynsk_I0434 + Tiszapolgar_ECA 0.153940224 > full output

Khvalynsk_I0434 + Tisza_LN 0.145465993 > full output

Khvalynsk_I0434 + Balkans_ChL 0.111720163 > full output

...

Khvalynsk_I0434 + Armenia_EBA 0.0108890099 > full output

Khvalynsk_I0434 + Armenia_ChL 0.00882375703 > full output

Khvalynsk_I0434 + Levant_BA_North 0.0078751978 > full output

Khvalynsk_I0434 + Minoan_Lasithi 0.0675240088 > full output

Khvalynsk_I0434 + Peloponnese_N 0.046998906 > full output

Khvalynsk_I0434 + Hajji_Firuz_ChL 0.00269860335 > full output

Khvalynsk_I0434 + Shahr_I_Sokhta_BA1 0.00261908387 > full output

Khvalynsk_I0434 + Sarazm_Eneolithic 0.00120345503 > full output

Khvalynsk_I0434 + Seh_Gabi_ChL 0.00111898703 > full output

Khvalynsk_I0434 + Geoksiur_Eneolithic 0.000178295163 > full output

Khvalynsk_I0434 + Tepe_Hissar_ChL 0.000163698274 > full output

Khvalynsk_I0434 + Bustan_BA 0.000151088148 > full output

Why is this potentially important? Because unless Khvalynsk_Eneolithic I0434 was a recent migrant from the North Caucasus piedmont steppe, which is where the remains of the Eneolithic steppe trio were excavated, then Yamnaya's ethnogenesis might not have anything at all to do with Asia or even the Caucasus region. At least not within any reasonable time frame anyway. Here's a map showing the geographic locations of all of the populations relevant to the highlighted mixture model above.

I won't be fussed if it turns out that the majority of the ancestry of the Yamnaya, Corded Ware and other closely related ancient peoples was sourced from the Eneolithic populations of the North Caucasus piedmont steppe. But I think it's useful to make the point that there are still very few ancient samples available from the steppes between the Black and Caspian seas, so we don't yet have much of a clue how the groups living throughout this region during the Eneolithic and earlier fit into the grand scheme of things.

Update 24/12/2018: I decided to repeat the analysis, but this time with Caucasus Hunter-Gatherers (CHG) as one of the outgroups (or right pops). The reason I initially didn't include CHG in the outgroups was because I didn't want to discriminate, perhaps unfairly, against West and Central Asians with high levels of CHG-related ancestry, and in favor of Europeans with no or minimal CHG-related input. But in my opinion, the new results clearly make more sense, with Sredny Stog and Varna at the top of the list.

Khvalynsk_I0434 + Sredny_Stog_I6561 0.410719649 > full output

Khvalynsk_I0434 + Varna 0.394089365 > full output

Khvalynsk_I0434 + Iberia_ChL 0.16554258 > full output

Khvalynsk_I0434 + Globular_Amphora 0.128348823 > full output

Khvalynsk_I0434 + Iberia_Central_CA 0.126100242 > full output

Khvalynsk_I0434 + Trypillia 0.135306664 > full output

Khvalynsk_I0434 + Baden_LCA 0.0853031796 > full output

Khvalynsk_I0434 + Protoboleraz_LCA 0.0766892008 > full output

Khvalynsk_I0434 + Tisza_LN 0.0661622403 > full output

Khvalynsk_I0434 + Tiszapolgar_ECA 0.0626469042 > full output

Khvalynsk_I0434 + Balaton_Lasinja_CA 0.0536293042 > full output

Khvalynsk_I0434 + ALPc_MN 0.0505788809 > full output

...

Khvalynsk_I0434 + Minoan_Lasithi 0.0439451605 > full output

Khvalynsk_I0434 + Balkans_ChL 0.0436885241 > full output

Khvalynsk_I0434 + Blatterhole_MN 0.0329758292 > full output

Khvalynsk_I0434 + Peloponnese_N 0.0181930605 > full output

Khvalynsk_I0434 + Armenia_EBA 0.014715999 > full output

Khvalynsk_I0434 + Armenia_ChL 0.0060437014 > full output

Khvalynsk_I0434 + Levant_BA_North 0.00514574731 > full output

Khvalynsk_I0434 + Shahr_I_Sokhta_BA1 0.00350059625 > full output

Khvalynsk_I0434 + Hajji_Firuz_ChL 0.00228771991 > full output

Khvalynsk_I0434 + Seh_Gabi_ChL 0.00117061206 > full output

Khvalynsk_I0434 + Sarazm_Eneolithic 0.001118931 > full output

Khvalynsk_I0434 + Bustan_BA 0.00021203609 > full output

Khvalynsk_I0434 + Tepe_Hissar_ChL 0.000200643323 > full output

Khvalynsk_I0434 + Geoksiur_Eneolithic 0.000175941977 > full output

Update 17/02/2019: I basically managed to confirm my analysis with samples from the Wang et al. Caucasus paper. See here.

See also...

Big deal of 2018: Yamnaya not related to Maykop

"The Homeland: In the footprints of the early Indo-Europeans" time map

Late PIE ground zero now obvious; location of PIE homeland still uncertain, but...

The Poltavka outlier

Anyone who still thinks that Y-chromosome haplogroup R1a originated in South Asia should burn this map into their brains. It'll come in useful over the next few years as we learn from ancient DNA about the conquest of the Indian subcontinent, and indeed much of Asia, by pastoralists from the western Russian and Ukrainian steppes.

X marks the spot of the burial site of Poltavka sample I0432 from the Mathieson et al. 2015 dataset. This individual belongs to Y-chromosome haplogroup R1a-Z93(Z94+), which today accounts for well over 90% of the R1a lineages in Asia and peaks in frequency at over 60% in the northern parts of South Asia.

Moreover, the dating of his burial site, 2925-2536 calBCE, suggests that he lived not long after the Z93 and Z94 mutations came into existence. That's because Z93 doesn't appear to be much older than 5,000 years based on full Y-chromosome sequence data (see here and here, including the comments).

So I0432 could well turn out to be a crucial piece in the puzzle of the peopling of South Asia.

Interestingly, this individual was flagged as an outlier in the Poltavka sample set by Mathieson et al., hence his other moniker: the Poltavka outlier. However, this wasn't because of any ancestry from South or even Central Asia. In fact, it was because he was too western.

Principal Component Analyses (PCA) featuring a wide range of present-day and ancient samples from Europe and Asia, like the one below, show that Poltavka outlier clusters further west than most Corded Ware individuals from Germany. Right click and open in a new tab to view full size.

In the past, using qpAdm, I modeled Poltavka outlier as 63.7% Yamnaya Samara and 36.3% German Middle Neolithic. This is probably not very far from the truth, but qpAdm offers a supervised mixture test in which the results are heavily reliant on the choice of outgroups, so I thought I'd revisit the issue with TreeMix, which allows an unsupervised analysis.

In a dataset including seven relatively high coverage Copper Age (CA), Early Bronze Age and Middle Neolithic (MN) European genomes, TreeMix picked out Poltavka outlier as the most likely sample to be admixed, showing a mixture edge of 33% from the base of the branch leading to the Iberian MN individual to that of Poltavka outlier.

This outcome is very similar to my qpAdm model, but it suggests an even more western source of admixture in Poltavka outlier. Could this admixture actually be from Iberia? I wouldn't discount this possibility, considering the presence of Bell Beaker communities, possibly of Atlantic or even Iberian origin, as far east as present-day Poland. Indeed, according to Cassidy et al. 2015, German Beakers show high affinity to MN and CA Iberians (see page 51 in the supp info here).

I double checked my TreeMix result with D-stats, and yep, when placed in a clade with Poltavka or Samara Yamnaya, Poltavka outlier shows the strongest signal of admixture from the Iberia MN individual.

At the same time, however, the signal from the Early Neolithic (EN) Iberian fails to reach significance (Z=<3), which suggests that, in fact, TreeMix and D-stats might be seeing the Iberia MN sample as the most attractive mixture source due to her high level of Western European hunter-gatherer (WHG) ancestry, which Poltavka outlier also has plenty of, rather than anything specific to Iberia.

In any case, it's clear enough that Poltavka outlier was the result of mixture between Yamnaya-related western steppe pastoralists and the descendants of Middle Neolithic Europeans with a high ratio of WHG ancestry. Where this admixture actually took place and which archaeological cultures were involved will have to be resolved with further sampling of ancient remains from Central and Eastern Europe.

However, it's already impossible to place the origin of Poltavka outlier anywhere in Asia, which suggests that both Z93 and Z94 are also from well inside the generally accepted borders of Europe.

This obviously has implications for the origins of the Indo-Iranians, because the widespread presence of these mutations in Asia gels very nicely with the idea, and indeed academic consensus, that Indo-Iranian languages expanded rapidly from the Eurasian steppe into Asia during the Bronze Age.

Considering that Poltavka outlier came from a Kurgan burial, and was therefore an individual of some social standing, he might be the direct ancestor of many millions of present-day Asians. If so, this won't be very difficult to prove in the near future as ancient DNA research revs up a few notches.

On a related note, apparently there's a paper on the way with ancient DNA results from Rakhigarhi, a Harappan site in Haryana, northern India (see here). As far as I know, the results will include Y-chromosome haplogroups of three males, but I don't think we'll see any decent genome-wide data at this stage. However, hopefully I'm wrong and the paper will come out with full ancient genomes.

Feel free to post your predictions in the comments. I'm tentatively expecting a couple of instances of J2 and maybe an L or H. Razib made basically the same prediction recently so I'm not being original. What I do know is that we won't see any R1a-Z93. The only way that might happen is if, say, someone coughed or sneezed on the Harappan remains.

Data source and reference...

Mathieson et al., Genome-wide patterns of selection in 230 ancient Eurasians, Nature, 528, 499–503 (24 December 2015), doi:10.1038/nature16152

See also...

The beast among Y-haplogroups