Yamnaya replaced Europeans, but admixed heavily as they spread to Asia


Recent papers The formation of human populations in South and Central Asia, by Narasimhan, Patterson et al. Science (2019) and An Ancient Harappan Genome Lacks Ancestry from Steppe Pastoralists or Iranian Farmers, by Shinde et al. Cell (2019).

NOTE. For direct access to Narasimhan, Patterson et al. (2019), visit this link courtesy of the first author and the Reich Lab.

I am currently not on holidays anymore, and the information in the paper is huge, with many complex issues raised by the new samples and analyses rather than solved, so I will stick to the Indo-European question, especially to some details that have changed since the publication of the preprint. For a summary of its previous findings, see the book series A Song of Sheep and Horses, in particular the sections from A Clash of Chiefs where I discuss languages and regions related to Central and South Asia.

I have updated the maps of the Preshistory Atlas, and included the most recently reported mtDNA and Y-DNA subclades. I will try to update the Eurasian PCA and related graphics, too.

NOTE. Many subclades from this paper have been reported by Kolgeh (download), Pribislav and Principe at Anthrogenica on this thread. I have checked some out for comparison, but even if it contradicted their analyses mine would be the wrong ones. I will upload my spreadsheets and link to them from this page whenever I find the time.

Ancestry clines (1) before and (2) after the advent of farming. Colour modified from the original to emphasize the CHG cline: notice the apparent relevance of forest-steppe groups in the formation of this CHG mating network from which Pre-Yamnaya peoples emerged.


I think the Narasimhan, Patterson et al. (2019) paper is well-balanced, and unexpectedly centered – as it should – on the spread of Yamnaya-related ancestry (now Western_Steppe_EMBA) as the marker of Proto-Indo-European migrations, which stretched ca. 3000 BC “from Hungary in the west to the Altai mountains in the east”, spreading later Indo-European dialects after admixing with local groups, from the Atlantic to South Asia.

I. Afanasievo

I.1. East or West PIE?

I expected Afanasievo to show (1) R1b-L23(xZ2103, xL51) and (2) R1b-L51 lineages, apart from (3) the known R1b-Z2103 ones, pointing thus to an ancestral PIE community before the typical Yamnaya bottlenecks, and with R1b-L51 supporting a connection with North-West Indo-European. The presence of some samples of hg. Q pointed in this direction, too.

However, Afanasievo samples show overwhelmingly R1b-Z2103 subclades (all except for those with low coverage), all apparently under R1b-Z2108 (formed ca. 3500 BC, TMRCA ca. 3500 BC), like most samples from East Yamnaya.

This necessarily shifts the split and spread of R1b-L23 lineages to Khvalynsk/early Repin-related expansions, in line with what TMRCA suggested, and what advances by Anthony (2019) and Khokhlov (2018) on future samples from the Reich Lab suggest.

Given the almost indistinguishable ancestry between Afanasievo and Early Yamnaya, there seems to be as of yet little potential information to support in population genomics that Pre-Tocharians were more closely related to North-West Indo-Europeans than to Graeco-Aryans, as it is proposed in linguistics based on the few shared traits between them, and the lack of innovations proper of the Graeco-Aryan community.

NOTE. A new issue of Wekʷos contains an abstract from a relevant paper by Blažek on vocabulary for ‘word’, including the common NWIE *wrdʰo-/wordʰo-, but also a new (for me, at least) Northern Indo-European one: *rēki-/*rēkoi̯-, shared by Slavic and Tocharian.

The fact that bottlenecks happened around the time of the late Repin expansion suggests that we might be able to see different clans based on the predominant lineages developing around the Don-Volga area in the 4th millennium BC. The finding of Pre-R1b-L51 in Lopatino (see below), and of a Catacomb sample of hg. R1b-Z2103(Z2105-) in the North Caucasus steppe near Novoaleksandrovskij also support a star-like phylogeny of R1b-L23 stemming from the Don-Volga area.

NOTE. Interestingly, a dismissal of a common trunk between Tocharian and North-West Indo-European would mean that shared similarities between such disparate groups could be traced back to a Common Late PIE trunk, and not to a shared (western) Repin community. For an example of such a ‘pure’ East-West dialectal division, see the diagram of Adams & Mallory (2007) at the end of the post. It would thus mean a fatal blow to Kortlandt’s Indo-Slavonic group among other hypothetical groupings (remade versions of the ancient Centum-Satem division), as well as to certain assumptions about laryngeal survival or tritectalism that usually accompany them. Still, I don’t think this is the case, so the question will remain a linguistic one, and maybe some similarities will be found with enough number of samples that differentiate Northern Indo-Europeans from the East Yamna/Catacomb-Poltavka-Balkan_EBA group.

Y-chromosome haplogroups of Afanasievo samples and neighbouring groups. See full maps.

I.2. Expansion or resurgence of hg. Q1b?

Haplogroup Q1b-Y6802(xY6798) seems to be the main lineage that expanded with Afanasievo, or resurged in their territory. It’s difficult to tell, because the three available samples are family, and belong to a later period.

NOTE. I have finally put some order to the chaos of Q1a vs. Q1b subclades in my spreadsheet and in the maps. The change of ISOGG 2016 to 2017 has caused that many samples reported as of Q1 subclades from papers prepared during the 2017-2018 period, and which did not provide specific SNP calls, were impossible to define with certainty. By checking some of them I could determine the specific standard used.

In favour of the presence of this haplogroup in the Pre-Yamnaya community are:

  • The statement by Anthony (2019) that Q1a [hence maybe Q1b in the new ISOGG nomenclature] represented a significant minority among an R1b-rich community.
  • The sample found in a Sintastha WSHG outlier (see below), of hg. Q1b-Y6798, and the sample from Lola, of hg. Q1b-L717, are thus from other lineage(s) separated thousands of years from the Afanasievo subclade, but might be related to the Khvalynsk expansion, like R1b-V1636 and R1b-M269 are.

These are the data that suggest multiple resurgence events in Afanasievo, rather than expanding Q1b lineages with late Repin:

  • Overwhelming presence of R1b in early Yamnaya and Afanasievo samples; one Q1(xQ1b) sample reported in Khvalynsk.
  • The three Q1b samples appear only later, although wide CI for radiocarbon dates, different sites, and indistinguishable ancestry may preclude a proper interpretation of the only available family.
    • Nevertheless, ancestry seems unimportant in the case of Afanasievo, since the same ancestry is found up to the Iron Age in a community of varied haplogroups.
  • Another sample of hg. Q1b-Y6802(xY6798) is found in Aigyrzhal_BA (ca. 2120 BC), with Central_Steppe_EMBA (WSHG-related) ancestry; however, this clade formed and expanded ca. 14000 BC.
  • The whole Altai – Baikal area seems to be a Q1b-L54 hotspot, although admittedly many subclades separated very early from each other, so they might be found throughout North Eurasia during the Neolithic.
  • One Afanasievo sample is reported as of hg. C in Shin (2017), and the same haplogroup is reported by Hollard (2014) for the only available sample of early Chemurchek to date, from Kulala ula, North Altai (ca. 2400 BC).
Y-chromosome haplogroups of late Afanasievo – early Chemurchek samples and neighbouring groups. See full maps.

I.3. Agricultural substrate

Evidence of continuous contacts of Central_Steppe_MLBA populations with BMAC from ca. 2100 BC on – visible in the appearance of Steppe ancestry among BMAC samples and BMAC ancestry among Steppe pastoralists – supports the close interaction between Indo-Iranian pastoralists and BMAC agriculturalists as the origin of the Asian agricultural substrate found in Proto-Indo-Iranian, hence likely related to the language of the Oxus Civilization.

Similar to the European agricultural substrate adopted by West Yamnaya settlers (both NWIE and Palaeo-Balkan speakers), Tocharian shows a few substrate terms in common with Indo-Iranian, which can be explained by contacts in different dialectal stages through phonetic reconstruction alone.

The recent Hermes et al. (2019) supports the early integration of pastoralism and millet cultivation in Central Asia (ca. 2700 BC or earlier), with the spread of agriculture to the north – through the Inner Asian Mountain Corridor – being thus unrelated to the Indo-Iranian expansions, which might support independent loans.

However, compared to the huge number of parallel shared loans between NWIE and Palaeo-Balkan languages in the European substratum, Indo-Iranians seem to have been the first borrowers of vocabulary from Asian agriculturalists, while Proto-Tocharian shows just one certain related word, with phonetic similarities that warrant an adoption from late Indo-Iranian dialects.

Y-chromosome haplogroups of Sintashta, Central Asia, and neighbouring groups in the Early Bronze Age. See full maps.

The finding of hg. (pre-)R1b-PH155 in a BMAC sample from Dzharkutan (to the west of Xinjiang) together with hg. R1b in a sample from Central Mongolia previously reported by Shin (2017) support the widespread presence of this lineage to the east and west of Xinjiang, which means it might have become incorporated to Indo-Iranian migrants into the Xiaohe horizon, to the Afanasievo-Chemurchek-derived groups, or the later from the former. In other words, the Island Biogeography Theory with its explanation of founder effects might be, after all, applicable to the whole Xinjiang area, not only during the Chemurchek – Tianshan-Beilu – Xiaohe interaction.

Of course, there is no need for too complicated models of haplogroup resurgence events in Central and South Asia, seeing how the total amount of hg. R1a-L657 (today prevalent among Indo-Aryan speakers from South Asia) among ancient Western/Central_Steppe_MLBA-related samples amounts to a total of 0, and that many different lineages survived in the region. Similar cases of haplogroup resurgence and Y-DNA bottleneck events are also found in the Central and Eastern Mediterranean, and in North-Eastern Europe. From the paper:

[It] could reflect stronger ecological or cultural barriers to the spread of people in South Asia than in Europe, allowing the previously established groups more time to adapt and mix with incoming groups. A second difference is the smaller proportion of Steppe pastoralist– related ancestry in South Asia compared with Europe, its later arrival by ~500 to 1000 years, and a lower (albeit still significant) male sex bias in the admixture (…).

Y-chromosome haplogroups of samples from the Srubna-Andronovo and Andronovo-related horizon, Xiaohe, late BMAC, and neighbouring groups. See full maps.

II. R1b-Beakers replaced R1a-CWC peoples

II.1. R1a-M417-rich Corded Ware

Newly reported Corded Ware samples from Radovesice show hg. R1a-M417, at least some of them xZ645, ‘archaic’ lineages shared with the early Bergrheinfeld sample (ca. 2650 BC) and with the coeval Esperstedt family, hence supporting that it eventually became the typical Western Corded Ware lineage(s), probably dominating over the so-called A-horizon and the Single Grave culture in particular. On the other hand, R1a-Z645 was typical of bottlenecks among expanding Eastern Corded Ware groups.

Interestingly, it is supported once again that known bottlenecks under hg. R1a-M417 happened during the Corded Ware expansion, evidenced also by the remarkable high variability of male lineages among early Corded Ware samples. Similarly, these Corded Ware samples from Bohemia form part of the typical ‘Central European’ cluster in the PCA, which excludes once again not only the ‘official’ Espersted outlier I1540, but also the known outlier with Yamnaya ancestry.

NOTE. The fact that Esperstedt is closely related geographically and in terms of ancestry to later Únětice samples further complicates the assumption that Únětice is a mixture of Bell Beakers and Corded Ware, being rather an admixture of incoming Bell Beakers with post-Yamnaya vanguard settlers who admixed with Corded Ware (see more on the expansion of Yamnaya ancestry). In other words, Únětice is rather an admixture of Yamnaya+EEF with Yamnaya+(CWC+EEF).

Y-chromosome haplogroups of samples from Catacomb, Poltavka, Balkan EBA, and Bell Beaker, as well as neighbouring groups. See full maps.

On Ukraine_Eneolithic I6561

If the bottlenecks are as straightforward as they appear, with a star-like phylogeny of R1a-M417 starting with the Pre-Corded Ware expansion, then what is happening with the Alexandria sample, so precisely radiocarbon dated to ca. 4045-3974 BC? The reported hg. R1a-M417 was fully compatible, while R1a-Z645 could be compatible with its date, but the few positive SNPs I got in my analysis point indeed to a potential subclade of R1a-Z94, and I trust more experienced hobbyists in this ‘art’ of ascertaining the SNPs of ancient samples, and they report hg. R1a-Z93 (Z95+, Y26+, Y2-).

Seeing how Y-DNA bottlenecks worked in Yamnaya-Afanasievo and in Corded Ware and related groups, and if this sample really is so deep within R1a-Z93 in a region that should be more strongly affected by the known Neolithic Y-chromosome bottlenecks and forest-steppe ecotone, someone from the lab responsible for this sample should check its date once again, before more people keep chasing their tails with an individual that (based on its derived SNPs’ TMRCA) might actually be dated to the Bronze Age, where it could make much more sense in terms of ancestry and position in the PCA.

EDIT (14 SEP 2019): … and with the fact that he is the first individual to show the genetic adaptation for lactase persistence (I3910-T), which is only found later among Bell Beakers, and much later in Sintashta and related Steppe_MLBA peoples (see comments below).

This is also evidenced by the other Ukraine_Eneolithic (likely a late Yamnaya) sample of hg. R1b-Z2103 from Dereivka (ca. 2800 BC) and who – despite being in a similar territory 1,000 years later – shows a wholly diluted Yamnaya ancestry under typically European HG ancestry, even more so than other late Sredni Stog samples from Dereivka of ca. 3600-3400 BC, suggesting a decrease in Steppe ancestry rather than an increase – which is supposedly what should be expected based on the ancestry from Alexandria…

Like the reported Chalcolithic individual of Hajji Firuz who showed an apparently incompatible subclade and Yamnaya ancestry at least some 1,000 years before it should, and turned out to be from the Iron Age (see below), this may be another case of wrong radiocarbon dating.

NOTE. It would be interesting, if this turns out to be another Hajji Firuz-like error, to check how well different ancestry models worked in whose hands exactly, and if anyone actually pointed out that this sample was derived, and not ancestral, to many different samples that were used in combination with it. It would also be a great control to check if those still supporting a Sredni Stog origin for PIE would shift their preference even more to the north or west, depending on where the first “true” R1a-M417 samples popped up. Such a finding now could be thus a great tool to discover whether haplogroup-based bias plays a role in ancestry magic as related to the Indo-European question, i.e. if it really is about “pure statistics”, or there is something else to it…

II.1. R1b-L51-rich Bell Beakers

The overwhelming majority of R1b-L51 lineages in Radovesice during the Bell Beaker period, just after the sampled Corded Ware individuals from the same site, further strengthen the hypothesis of an almost full replacement of R1a-M417 lineages from Central Europe up to southern Scandinavia after the arrival of Bell Beakers.

Yet another R1b-L151* sample has popped up in Central Europe, in the individual classified as Bilina_BA (ca. 2200-800 BC), which clusters with Bell Beakers from Bohemia, with the outlier from Turlojiškė, and with Early Slavs, suggesting once again that a group of central-east European Beakers represented the Pre-Proto-Balto-Slavic community before their spread and admixture events to the east.

The available ancient distribution of R1b-L51*, R1b-L52* or R1b-L151* is getting thus closer to the most likely origin of R1b-L51 in the expansion of East Bell Beakers, who trace their paternal ancestors to Yamnaya settlers from the Carpathian Basin:

NOTE. Some of these are from other sources, and some are samples I have checked in a hurry, so I may have missed some derived SNPs. If you send me a corrected SNP call to dismiss one of these, or more ‘archaic’ samples, I’ll correct the map accordingly. See also maps of modern distributionof R1b-M269 subclades.

Distribution of ‘archaic’ R1b-L51 subclades in ancient samples, overlaid over a map of Yamnaya and Bell Beaker migrations. In blue, Yamnaya Pre-L51 from Lopatino (not shown) and R1b-L52* from BBC Augsburg. In violet, R1b-L51 (xP312,xU106) from BBC Prague and Poland. In maroon, hg. R1b-L151* from BBC Hungary, BA Bohemia, and (not shown) a potential sample from BBC at Mondelange, which is certainly xU106, maybe xP312. Interestingly, the earliest sample of hg. R1b-U106 (a lineage more proper of northern Europe) has been found in a Bell Beaker from Radovesice (ca. 2350 BC), between two of these ‘archaic’ R1b-L51 samples; and a sample possibly of hg. R1b-ZZ11+ (ancestral to DF27 and U152) was found in a Bell Beaker from Quedlinburg, Germany (ca. 2290 BC), to the north-west of Bohemia. The oldest R1b-U152 are logically from Central Europe, too.

III. Proto-Indo-Iranian

Before the emergence of Proto-Indo-Iranian, it seems that Pre-Proto-Indo-Iranian-speaking Poltavka groups were subjected to pressure from Central_Steppe_EMBA-related peoples coming from the (south-?)east, such as those found sampled from Mereke_BA. Their ‘kurgan’ culture was dated correctly to approximately the same date as Poltavka materials, but their ancestry and hg. N2(pre-N2a) – also found in a previous sample from Botai – point to their intrusive nature, and thus to difficulties in the Pre-Proto-Indo-Iranian community to keep control over the previous East Yamnaya territory in the Don-Volga-Ural steppes.

We know that the region does not show genetic continuity with a previous period (or was not under this ‘eastern’ pressure) because of an Eastern Yamnaya sample from the same site (ca. 3100 BC) showing typical Yamnaya ancestry. Before Yamnaya, it is likely that Pre-Yamnaya ancestry formed through admixture of EHG-like Khvalynsk with a North Caspian steppe population similar to the Steppe_Eneolithic samples from the North Caucasus Piedmont (see Anthony 2019), so we can also rule out some intermittent presence of a Botai/Kelteminar-like population in the region during the Khvalynsk period.

It is very likely, then, that this competition for the same territory – coupled with the known harsher climate of the late 3rd millennium BC – led Poltavka herders to their known joint venture with Abashevo chiefs in the formation of the Sintashta-Potapovka-Filatovka community of fortified settlements. Supporting these intense contacts of Poltavka herders with Central Asian populations, late ‘outliers’ from the Volga-Ural region show admixture with typical Central_Steppe_MLBA populations: one in Potapovka (ca. 2220 BC), of hg. R1b-Z2103; and four in the Sintashta_MLBA_o1 cluster (ca. 2050-1650 BC), with two samples of hg. R1b-L23 (one R1b-Z2109), one Q1b-L56(xL53), one Q1b-Y6798.

Outlier analysis reveals ancient contacts between sites. We plot the average of principal component 1 (x axis) and principal component 2 (y axis) for the West Eurasian and All Eurasian PCA plots (…). In the Middle to Late Bronze Age Steppe, we observe, in addition to the Western_Steppe_MLBA and Central_Steppe_MLBA clusters (indistinguishable in this projection), outliers admixed with other ancestries. The BMAC-related admixture in Kazakhstan documents northward gene flow onto the Steppe and confirms the Inner Asian Mountain Corridor as a conduit for movement of people.

Similar to how the Sintashta_MLBA_o2 cluster shows an admixture with central steppe populations and hg. R1a-Z645, the WSHG ancestry in those outliers from the o1 cluster of typically (or potentially) Yamnaya lineages show that Poltavka-like herders survived well after centuries of Abashevo-Poltavka coexistence and admixture events, supporting the formation of a Proto-Indo-Iranian community from the local language as pronounced by the incomers, who dominated as elites over the fortified settlements.

The Proto-Indo-Iranian community likely formed thus in situ in the Don-Volga-Ural region, from the admixture of locals of Yamnaya ancestry with incomers of Corded Ware ancestry – represented by the ca. 67% Yamnaya-like ancestry and ca. 33% ancestry from the European cline. Their community formed thus ca. 1,000 years later than the expansion of Late PIE ca. 3500 BC, and expanded (some 500 years after that) a full-fledged Proto-Indo-Iranian language with the Srubna-Andronovo horizon, further admixing with ca. 9% of Central_Steppe_EMBA (WSHG-related) ancestry in their migration through Central Asia, as reported in the paper.

IV. Armenian

The sample from Hajji Firuz, of hg. R1b-Z2103 (xPF331), has been – as expected – re-dated to the Iron Age (ca. 1193-1019 BC), hence it may offer – together with the samples from the Levant and their Aegean-like ancestry rapidly diluted among local populations – yet another proof of how the Late Bronze Age upheaval in Europe was the cause of the Armenian migration to the Armenoid homeland, where they thrived under the strong influence from Hurro-Urartian.

Y-chromosome haplogroups of the Middle East and neighbouring groups during the Late Bronze Age / Iron Age. See full maps.

Indus Valley Civilization and Dravidian

A surprise came from the analysis reported by Shinde et al. (2019) of an Iran_N-related IVC ancestry which may have split earlier than 10000 BC from a source common to Iran hunter-gatherers of the Belt Cave.

For the controversial Elamo-Dravidian hypothesis of the Muscovite school, this difference in ancestry between both groups (IVC and Iran Neolithic) seems to be a death blow, if population genomics was even needed for that. Nevertheless, I guess that a full rejection of a recent connection will come down to more recent and subtle population movements in the area.

EDIT (12 SEP): Apparently, Iosif Lazaridis is not so sure about this deep splitting of ‘lineages’ as shown in the paper, so we may be talking about different contributions of AME+ANE/ENA, which means the Elamo-Dravidian game is afoot; at least in genomics:

I shared the idea that the Indus Valley Civilization was linked to the Proto-Dravidian community, so I’m inclined to support this statement by Narasimhan, Patterson, et al. (2019), even if based only on modern samples and a few ancient ones:

The strong correlation between ASI ancestry and present-day Dravidian languages suggests that the ASI, which we have shown formed as groups with ancestry typical of the Indus Periphery Cline moved south and east after the decline of the IVC to mix with groups with more AASI ancestry, most likely spoke an early Dravidian language.

Natural neighbour interpolation of qpAdm results – Maximum A Posteriori Estimate from the Hierarchical Model (estimates used in the Narasimhan, Patterson et al. 2019 figures) for Central_Steppe_MLBA-related (left), Indus_Periphery_West-related (center) and Andamanese_Hunter-Gatherer-related ancestry (right) among sampled modern Indian populations. In blue, peoples of IE language; in red, Dravidian; in pink, Tibeto-Burman; in black, unclassified. See full image.

I am wary of this sort of simplistic correlation with modern speakers, because we have seen what happened with the wrong assumptions about modern Balto-Slavic and Finno-Ugric speakers and their genetic profile (see e.g. here or here). In fact, I just can’t differentiate as well as those with deep knowledge in South Asian history the social stratification of the different tribal groups – with their endogamous rules under the varna and jati systems – in the ancestry maps of modern India. The pattern of ancestry and language distribution combined with the findings of ancient populations seem in principle straightforward, though.


The message to take home from Shinde et al. (2019) is that genomic data is fully at odds with the Anatolian homeland hypothesis – including the latest model by Heggarty (2014)* – whose relevance is still overvalued today, probably due in part to the shift of OIT proponents to more reasonable Out-of-Iran models, apparently more fashionable as a vector of Indo-Aryan languages than Eurasian steppe pastoralists?
*The authors listed this model erroneously as Heggarty (2019).

The paper seems to play with the occasional reference to Corded Ware as a vector of expansion of Indo-European languages, even after accepting the role of Yamnaya as the most evident population expanding Late PIE to western Europe – and the different ancestry that spread with Indo-Iranian to South Asia 1,000 years later. However, the most cringe-worthy aspect is the sole citation of the debunked, pseudoscientific glottochronological method used by Ringe, Warnow, and Taylor (2002) to support the so-called “steppe homeland”, a paper and dialectal scheme which keeps being referenced in papers of the Reich Lab, probably as a consequence of its use in Anthony (2007).

On the other hand, these are the equivalent simplistic comments in Narasimhan, Patterson et al. (2019):

The Steppe ancestry in South Asia has the same profile as that in Bronze Age Eastern Europe, tracking a movement of people that affected both regions and that likely spread the unique features shared between Indo-Iranian and Balto-Slavic languages. (…), which despite their vast geographic separation share the “satem” innovation and “ruki” sound laws.

Indo-European dialectal relationships, from Mallory and Adams (2006).

The only academic closely related to linguistics from the list of authors, as far as I know, is James P. Mallory, who has supported a North-West Indo-European dialect (including Balto-Slavic) for a long time – recently associating its expansion with Bell Beakers – opposed thus to a Graeco-Aryan group which shared certain innovations, “Satemization” not being one of them. Not that anyone needs to be a linguist to dismiss any similarities between Balto-Slavic and Indo-Iranian beyond this phonetic trend, mind you.

Even Anthony (2019) supports now R1b-rich Pre-Yamnaya and Yamnaya communities from the Don-Volga region expanding Middle and Late Proto-Indo-European dialects.

So how does the underlying Corded Ware ancestry of eastern Europe (where Pre-Balto-Slavs eventually spread to from Bell Beaker-derived groups) and of the highly admixed (“cosmopolitan”, according to the authors) Sintashta-Potapovka-Filatovka in the east relate to the similar-but-different phonetic trends of two unrelated IE dialects?

If only there was a language substrate that could (as Shinde et al. put it) “elegantly” explain this similar phonetic evolution, solving at the same time the question of the expansion of Uralic languages and their strong linguistic contacts with steppe peoples. Say, Eneolithic populations of mainly hunter-fisher-gatherers from the North Pontic forest-steppes with a stronger connection to metalworking


Migrations in the Levant region during the Chalcolithic, also marked by distinct Y-DNA


Open access Ancient DNA from Chalcolithic Israel reveals the role of population mixture in cultural transformation, by Harney et al. Nature Communications (2018).

Interesting excerpts (emphasis mine, reference numbers deleted for clarity):


The material culture of the Late Chalcolithic period in the southern Levant contrasts qualitatively with that of earlier and later periods in the same region. The Late Chalcolithic in the Levant is characterized by increases in the density of settlements, introduction of sanctuaries, utilization of ossuaries in secondary burials, and expansion of public ritual practices as well as an efflorescence of symbolic motifs sculpted and painted on artifacts made of pottery, basalt, copper, and ivory. The period’s impressive metal artifacts, which reflect the first known use of the “lost wax” technique for casting of copper, attest to the extraordinary technical skill of the people of this period.

The distinctive cultural characteristics of the Late Chalcolithic period in the Levant (often related to the Ghassulian culture, although this term is not in practice applied to the Galilee region where this study is based) have few stylistic links to the earlier or later material cultures of the region, which has led to extensive debate about the origins of the people who made this material culture. One hypothesis is that the Chalcolithic culture in the region was spread in part by immigrants from the north (i.e., northern Mesopotamia), based on similarities in artistic designs. Others have suggested that the local populations of the Levant were entirely responsible for developing this culture, and that any similarities to material cultures to the north are due to borrowing of ideas and not to movements of people.

Previous genome-wide ancient DNA studies from the Near East have revealed that at the time when agriculture developed, populations from Anatolia, Iran, and the Levant were approximately as genetically differentiated from each other as present-day Europeans and East Asians are today. By the Bronze Age, however, expansion of different Near Eastern agriculturalist populations — Anatolian, Iranian, and Levantine — in all directions and admixture with each other substantially homogenized populations across the region, thereby contributing to the relatively low genetic differentiation that prevails today. Showed that the Levant Bronze Age population from the site of ‘Ain Ghazal, Jordan (2490–2300 BCE) could be fit statistically as a mixture of around 56% ancestry from a group related to Levantine Pre-Pottery Neolithic agriculturalists (represented by ancient DNA from Motza, Israel and ‘Ain Ghazal, Jordan; 8300–6700 BCE) and 44% related to populations of the Iranian Chalcolithic (Seh Gabi, Iran; 4680–3662 calBCE). Suggested that the Canaanite Levant Bronze Age population from the site of Sidon, Lebanon (~1700 BCE) could be modeled as a mixture of the same two groups albeit in different proportions (48% Levant Neolithic-related and 52% Iran Chalcolithic-related). However, the Neolithic and Bronze Age sites analyzed so far in the Levant are separated in time by more than three thousand years, making the study of samples that fill in this gap, such as those from Peqi’in, of critical importance.

This procedure produced genome-wide data from 22 ancient individuals from Peqi’in Cave (4500–3900 calBCE) (…)


We find that the individuals buried in Peqi’in Cave represent a relatively genetically homogenous population. This homogeneity is evident not only in the genome-wide analyses but also in the fact that most of the male individuals (nine out of ten) belong to the Y-chromosome haplogroup T, a lineage thought to have diversified in the Near East. This finding contrasts with both earlier (Neolithic and Epipaleolithic) Levantine populations, which were dominated by haplogroup E, and later Bronze Age individuals, all of whom belonged to haplogroup J.

Detailed sample background data for each of the 22 samples from which we successfully obtained ancient DNA. Additionally, background information for all samples from Peqi’in that were screened is included in Supplementary Data 1. *Indicates that Y-chromosome haplogroup call should be interpreted with caution, due to low coverage data.

Our finding that the Levant_ChL population can be well-modeled as a three-way admixture between Levant_N (57%), Anatolia_N (26%), and Iran_ChL (17%), while the Levant_BA_South can be modeled as a mixture of Levant_N (58%) and Iran_ChL (42%), but has little if any additional Anatolia_N-related ancestry, can only be explained by multiple episodes of population movement. The presence of Iran_ChL-related ancestry in both populations – but not in the earlier Levant_N – suggests a history of spread into the Levant of peoples related to Iranian agriculturalists, which must have occurred at least by the time of the Chalcolithic. The Anatolian_N component present in the Levant_ChL but not in the Levant_BA_South sample suggests that there was also a separate spread of Anatolian-related people into the region. The Levant_BA_South population may thus represent a remnant of a population that formed after an initial spread of Iran_ChL-related ancestry into the Levant that was not affected by the spread of an Anatolia_N-related population, or perhaps a reintroduction of a population without Anatolia_N-related ancestry to the region. We additionally find that the Levant_ChL population does not serve as a likely source of the Levantine-related ancestry in present-day East African populations.

These genetic results have striking correlates to material culture changes in the archaeological record. The archaeological finds at Peqi’in Cave share distinctive characteristics with other Chalcolithic sites, both to the north and south, including secondary burial in ossuaries with iconographic and geometric designs. It has been suggested that some Late Chalcolithic burial customs, artifacts and motifs may have had their origin in earlier Neolithic traditions in Anatolia and northern Mesopotamia. Some of the artistic expressions have been related to finds and ideas and to later religious concepts such as the gods Inanna and Dumuzi from these more northern regions. The knowledge and resources required to produce metallurgical artifacts in the Levant have also been hypothesized to come from the north.

Our finding of genetic discontinuity between the Chalcolithic and Early Bronze Age periods also resonates with aspects of the archeological record marked by dramatic changes in settlement patterns, large-scale abandonment of sites, many fewer items with symbolic meaning, and shifts in burial practices, including the disappearance of secondary burial in ossuaries. This supports the view that profound cultural upheaval, leading to the extinction of populations, was associated with the collapse of the Chalcolithic culture in this region.

Genetic structure of analyzed individuals. a Principal component analysis of 984 present-day West Eurasians (shown in gray) with 306 ancient samples projected onto the first two principal component axes and labeled by culture. b ADMIXTURE analysis of 984 and 306 ancient samples with K = 11
ancestral components. Only ancient samples are shown


I think the most interesting aspect of this paper is – as usual – the expansion of peoples associated with a single Y-DNA haplogroup. Given that the expansion of Semitic languages in the Middle East – like that of Anatolian languages from the north – must have happened after ca. 3100 BC, coinciding with the collapse of the Uruk period, these Chalcolithic north Levant peoples are probably not related to the posterior Semitic expansion in the region. This can be said to be supported by their lack of relationship with posterior Levantine migrations into Africa. The replacement of haplogroup E before the arrival of haplogroup J suggests still more clearly that Natufians and their main haplogroup were not related to the Afroasiatic expansions.

Distribution of Semitic languages. From Wikipedia.

On the other hand, while their ancestry points to neighbouring regional origins, their haplogroup T1a1a (probably T1a1a1b2) may be closely related to that of other Semitic peoples to the south, as found in east Africa and Arabia. This may be due either to a northern migration of these Chalcolithic Levantine peoples from southern regions in the 5th millennium BC, or maybe to a posterior migration of Semitic peoples from the Levant to the south, coupled with the expansion of this haplogroup, but associated with a distinct population. As we know, ancestry can change within certain generations of intense admixture, while Y-DNA haplogroups are not commonly admixed in prehistoric population expansions.

Without more data from ancient DNA, it is difficult to say. Haplogroup T1a1 is found in Morocco (ca. 3780-3650 calBC), which could point to a recent expansion of a Berbero-Semitic branch; but also in a sample from Balkans Neolithic ca. 5800-5400 calBCE, which could suggest an Anatolian origin of the specific subclades encountered here. In any case, a potential origin of Proto-Semitic anywhere near this wide Near Eastern region ca. 4500-3500 BC cannot be discarded, knowing that their ancestors came probably from Africa.

Distribution of haplogroup T of Y-chromosome. From Wikipedia.

Interesting from this paper is also that we are yet to find a single prehistoric population expansion not associated with a reduction of variability and expansion of Y-DNA haplogroups. It seems that the supposedly mixed Yamna community remains the only (hypothetical) example in history where expanding patrilineal clans will not share Y-DNA haplogroup…


Mitogenomes show continuity of Neolithic populations in Southern India

New paper (behind paywall) Neolithic phylogenetic continuity inferred from complete mitochondrial DNA sequences in a tribal population of Southern India, by Sylvester et al. Genetica (2018).

This paper used a complete mtDNA genome study of 113 unrelated individuals from the Melakudiya tribal population, a Dravidian speaking tribe from the Kodagu district of Karnataka, Southern India.

Some interesting excerpts (emphasis mine):

Autosomal genetic evidence indicates that most of the ethnolinguistic groups in India have descended from a mixture of two divergent ancestral populations: Ancestral North Indians (ANI) related to People of West Eurasia, the Caucasus, Central Asia and the Middle East, and Ancestral South Indians (ASI) distantly related to indigenous Andaman Islanders (Reich et al. 2009). It is presumed that proto-Dravidian language, most likely originated in Elam province of South Western Iran, and later spread eastwards with the movement of people to the Indus Valley and later the subcontinent India (McAlpin et al. 1975; Cavalli-Sforza et al. 1988; Renfrew 1996; Derenko et al. 2013). West Eurasian haplogroups are found across India and harbor many deep-branching lineages of Indian mtDNA pool, and most of the mtDNA lineages of Western Eurasian ancestry must have a recent entry date less than 10 Kya (Kivisild et al. 1999a). The frequency of these lineages is specifically found among the higher caste groups of India (Bamshad et al. 1998, 2001; Basu et al. 2003) and many caste groups are direct descendants of Indo-Aryan immigrants (Cordaux et al. 2004). These waves of various invasions and subsequent migrations resulted in major demographic expansions in the region, which added new languages and cultures to the already colonized populations of India. Although previous genetic studies of the maternal gene pools of Indians had revealed a genetic connection between Iranian populations and the Arabian Peninsula, likely the result of both ancient and recent gene flow (Metspalu et al. 2004; Terreros et al. 2011).


Haplogroup HV14

mtDNA haplogroup HV14 has prominence in North/Western Europe, West Eurasia, Iran, and South Caucasus to Central Asia (Malyarchuk et al. 2008; Schonberg et al. 2011; Derenko et al. 2013; De Fanti et al. 2015). Although Palanichamy identified haplogroup HV14a1 in three Indian samples (Palanichamy et al. 2015), it is restricted to limited unknown distribution. In the present study, by the addition of considerable sequences from the Melakudiya population, a unique novel subclade designated as HV14a1b was found with a high frequency (43%) allowed us to reveal the earliest diverging sequences in the HV14 tree prior to the emergence of HV14a1b in Melakudiya. (…) The coalescence age for haplogroup HV14 in this study is dated ~ 16.1 ± 4.2 kya and the founder age of haplogroup HV14 in Melakudiya tribe, which is represented by a novel clade HV14a1b is ~ 8.5 ± 5.6 kya

Maximum Parsimonious tree of complete mitogenomes constructed using 38 sequences from Melakudiya tribe and 11 previously published sequences belonging to haplogroup HV14 [Supplementary file Table S2] Suffixes @ indicate back mutation, a plus sign (+) an insertion. Control region mutations are underlined, and synonymous transitions are shown in normal font and non-synonymous mutations are shown in bold font. Coalescence ages (Kya) for complete coding region are shown in normal font and synonymous transitions are shown in Italics

Haplogroup U7a3a1a2

The coalescence age of haplogroup U7a3a1a2 dates to ~ 13.3 ± 4.0 kya. (…)

Although, haplogroup U7 has its origin from the Near East and is widespread from Europe to India, the phylogeny of Melakudiya tribe with subclade U7a3a1a2 clusters with populations of India (caste and tribe) and neighboring populations (Irwin et al. 2010; Ranaweera et al. 2014; Sahakyan et al. 2017), hint about the in-situ origin of the subclade in India from Indo-Aryan immigrants.

I am not a native English speaker, but this paper looks like it needs a revision by one.

Also – without comparison with ancient DNA – it is not enough to show coalescence age to prove an origin of haplogroup expansion in the Neolithic instead of later bottlenecks. However, since we are talking about mtDNA, it is likely that their analysis is mostly right.

Finally, one thing is to prove that the origin of the Indus Valley Civilization lies (in part) in peoples from the Iranian plateau, and to show with ASI ancestry that they are probably the origin of Proto-Dravidian expansion, and another completely different thing is to prove an Elamo-Dravidian connection.

Since that group is not really accepted in linguistics, it is like talking about proving – through that Iran Neolithic ancestry – a Sumero-Dravidian, or a Hurro-Dravidian connection…


Expansion of domesticated goat echoes expansion of early farmers


New paper (behind paywall) Ancient goat genomes reveal mosaic domestication in the Fertile Crescent, by Daly et al. Science (2018) 361(6397):85-88.

Interesting excerpts (emphasis mine):

Thus, our data favor a process of Near Eastern animal domestication that is dispersed in space and time, rather than radiating from a central core (3, 11). This resonates with archaeozoological evidence for disparate early management strategies from early Anatolian, Iranian, and Levantine Neolithic sites (12, 13). Interestingly, our finding of divergent goat genomes within the Neolithic echoes genetic investigation of early farmers. Northwestern Anatolian and Iranian human Neolithic genomes are also divergent (14–16), which suggests the sharing of techniques rather than large-scale migrations of populations across Southwest Asia in the period of early domestication. Several crop plants also show evidence of parallel domestication processes in the region (17).

PCA affinity (Fig. 2), supported by qpGraph and outgroup f3 analyses, suggests that modern European goats derive from a source close to the western Neolithic; Far Eastern goats derive from early eastern Neolithic domesticates; and African goats have a contribution from the Levant, but in this case with considerable admixture from the other sources (figs. S11, S16, and S17 and tables S26 and 27). The latter may be in part a result of admixture that is discernible in the same analyses extended to ancient genomes within the Fertile Crescent after the Neolithic (figs. S18 and S19 and tables S20, S27, and S31) when the spread of metallurgy and other developments likely resulted in an expansion of inter-regional trade networks and livestock movement.

Maximumlikelihood phylogeny and geographical distributions of ancient mtDNA haplogroups. (A) A phylogeny placing ancient whole mtDNA sequences in the context of known haplogroups. Symbols denoting individuals are colored by clade membership; shape indicates archaeological period (see key). Unlabeled nodes are modern bezoar and outgroup sequence (Nubian ibex) added for reference.We define haplogroup T as the sister branch to the West Caucasian tur (9). (B and C) Geographical distributions of haplogroups show early highly structured diversity in the Neolithic period (B) followed by collapse of structure in succeeding periods (C).We delineate the tiled maps at 7250 to 6950 BP, a period >bracketing both our earliest Chalcolithic sequence (24, Mianroud) and latest Neolithic (6, Aşağı Pınar). Numbered archaeological sites also include Direkli Cave (8), Abu Ghosh (9), ‘Ain Ghazal (10), and Hovk-1 Cave (11) (table S1) (9).

Our results imply a domestication process carried out by humans in dispersed, divergent, but communicating communities across the Fertile Crescent who selected animals in early millennia, including for pigmentation, the most visible of domestic traits.


Kura-Araxes implicated in the transformation of regional trade in the Near East


Craft production at Köhne Shahar, a Kura-Araxes settlement in Iranian Azerbaijan, by Alizadeh et al. J Anthropol Arch (2018) 51:127-143.

Interesting excerpts (emphasis mine):


Kura-Araxes communities first emerged throughout the southern Caucasus in the mid-4th millennium BC (Sagona, 1984; Rothman, 2005; Kohl, 2009) or possibly earlier in Nakhchivan (Marro et al., 2014; Palumbi and Chataigner, 2014: 250; Marro et al., 2015; Palumbi and Chataigner, 2015). By the late 4th-early 3rd millennium BC, their characteristic material culture, particularly hand-made black burnished pottery, spread throughout much of Southwest Asia after 2900 BCE (Fig. 1). The widespread dissemination of this material culture, along with the small size of most sites, the ephemeral nature of their architectural remains in these smaller sites, and their presence in both fertile lowlands and seasonally-inhospitable highlands, have been used to portray Kura-Araxes communities as small, egalitarian communities of mobile pastoralists or sedentary agriculturalists; economically undifferentiated and socially non-hierarchical (Smith, 2005: 258; Frangipane and Palumbi, 2007; Kohl, 2007: 113; 2009: 250). Limited evidence for craft production and trade among Kura-Araxes communities has further strengthened the argument that Kura-Araxes economies were dominated by domestic and subsistence-related activities (Palumbi, 2008: 53). With some rare exceptions (Marro et al., 2010; Stöllner, 2014; Simonyan and Rothman, 2015), Kura-Araxes settlements lack any evidence of craft production, mining, or resource extraction.

Distribution of Kura-Araxes material culture in the Near East (modified from Wikimedia).

Kura-Araxes communities, however, are also implicated in the evolution and transformation of regional trade in the Near East. Cause and effect of the spread of Kura-Araxes material culture beyond the Caucasus “homeland” and the establishment of diaspora is hotly debated. Among proponents of emigration, the strongest arguments for movement out of the Caucasus include the presence of strong pull factors, notably productive activities like meat and wool production, viticulture, and metals and metallurgy (Rothman, 2003). Kura-Araxes populations primarily inhabited mountains and intermontane valleys of the highland zone surrounding Mesopotamia. Kura-Araxes communities had access to metals, precious and semi-precious stones, stones for tool making, wood, and animal products; resources that were abundant in the mountain zone, yet critical to the evolution of Mesopotamian societies. The frequent appearance of simple bronze and copper objects at temporary camps of Kura-Araxes herders suggests that mobile agropastoralists engaged in metallurgy and trade in metals, especially with societies of the Upper Euphrates (Frangipane et al., 2001; Hauptmann et al., 2002; Rothman, 2003; Connor and Sagona, 2007; Frangipane, 2014). Wool and textiles products from sheep herded by mountainous communities may have been major exports of the mountain zone to Mesopotamia (Anthony, 2007: 284; Nosch et al., 2013; Breniquet and Michel, 2014).

It is argued that by the second half of the 4th millennium BC (Surenhagen, 1986; Algaze, 1989, 2004, 2007), Uruk polities of southern Mesopotamia established colonies across northern Mesopotamia, southern Anatolia, and western Iran to better control regional trade. Although the nature of these colonies and local developments is still debated (Stein, 2002, 2014), co-occurrence of the sudden abandonment of these colonies and regional expansion of Kura-Araxes communities by the end of the 4th millennium BC has led some scholars to argue that Kura-Araxes communities were emergent competitors of Mesopotamia whose economic activities possibly contributed to the decline and eventual collapse of the Uruk system (Algaze, 2001: 76; Kohl, 2007: 97–98; Lamberg-Karlovsky, 2008: 10).

Major Kura-Araxes sites in the Caucasus region and location of Köhne Shahar (modified map from wikimedia.org).


The abundant evidence of craft specialization at Köhne Shahar clearly shows that Kura-Araxes communities were not all homogenous and undifferentiated. Excavations and a geophysical survey at Köhne Shahar demonstrate that multi-craft production activities were practiced at a community-level inside the citadel at the site, and that a large portion of the population may have engaged in this specialized, extrahousehold craft economy. The possible involvement of a political apparatus with a specialized craft economy at Köhne Shahar may have necessitated control over various aspects of production such as labor, commodities, resource procurement, exchange, and grain storage. As Adam Smith (Smith, 2015: 106) argues, all of these point to complex labor coordination at Köhne Shahar.

Although excavations exposed a limited area, the scale of craft production at Köhne Shahar and the scarcity of finished products may suggest that consumers of finished goods were not necessarily residents of Köhne Shahar, but instead occupied other areas on the landscape. Communication between these nodes of production and consumption necessitated a network of exchange and interaction. The miniature sumptuary container at Köhne Shahar points to possible interaction with regions of Central Asia and the Persian Gulf, while the bitumen used to mend vessels points to interaction with northern Mesopotamia or the Zagros mountains in western Iran. It is possible that long-distance interaction brought Köhne Shahar chiefs into contact with other complex societies in the region, connecting them to a larger inter-regional exchange and trade network.

Archaeological and geophysical evidence for community-level production documents Köhne Shahar’s emergence as a regional economic center. The extent of Köhne Shahar’s regional engagements and ambitions, however, have yet to be fully understood. Köhne Shahar’s economic focus on production may have enabled its producers to contribute to regional transformations. When trade became a significant part of the economy of early complex societies in the Near East in the second half of the 4th millennium BC (Surenhagen, 1986; Algaze, 1989, 2004), Kura-Araxes communities like Köhne Shahar may have emerged as a primary center of specialized craft production in the late 4th/early 3rd millennium BC. Alternatively, Köhne Shahar’s economic success may have been due to its ability to satisfy regional demand (highlands of NW Iran, eastern Anatolia, or northern Mesopotamia) by filling a supply vacuum created following the collapse of Uruk colonies. Political and entrepreneurial ambitions of Köhne Shahar chiefs may have also provided the impetus for the selection of the site’s naturally defensible area and the construction of a large and defensive fortification wall; two barriers intended to safeguard the production machinery of the citadel from the onset of the site’s occupation in the late 4th millennium BC (Alizadeh et al., 2015).

I don’t have much to add to what I recently wrote about potential intrusive steppe admixture in the Caucasus.


On the Maykop – Upper Mesopotamia cultural province, distinct from the steppe


New paper (behind paywall) The Production of Thin‐Walled Jointless Gold Beads from the Maykop Culture Megalithic Tomb of the Early Bronze Age at Tsarskaya in the North Caucasus: Results of Analytical and Experimental Research, by Trifonov et al. Archaeometry (2018)

Interesting excerpts (emphasis mine):

In 1898, two megalithic tombs containing graves of a local social elite dated to the Early Bronze Age were discovered by N. I. Veselovsky near the village of Tsarskaya (modern Novosvobodnaya, Republic of Adygeya) (Fig. 1 (a)) (Baye 1900, 43–59; IAC 1901, 33–8; Sagona 2018, 281–97).

Radiocarbon dates place both tombs within the Novosvobodnaya phase of the Maykop culture, between c. 3200 and 2900 BC (Trifonov et al. 2017). Along with the human remains (one adult individual was interred in each dolmen), the tombs yielded rich funerary offerings, including artefacts made of gold, silver and semi-precious stones. (…) This paper presents results of a technical analysis of just one type of artefact, from kurgan 2 at Tsarskaya: thin-walled jointless beads made from gold.

(a) A map of the Caucasus and part of Western Asia, showing the locations of sites mentioned in the text: 1, Tsarskaya (modern Novosvobodnaya); 2, Maykop; 3, Staromyshastovskaya; 4, Andryukovskaya; 5, Psebaiskaya; 6, Inozemtsevo; 7, Kudakhurt; 8, Soyuq Bulaq; 9, Sé Girdan; 10, Tepe Gawra. (b) The string of thin-walled jointless gold beads, silver and carnelian beads from the dolmen in kurgan 2 at Tsarskaya, Western Caucasus (1898).

Ever since M. I. Rostovtzeff noted a stylistic similarity between Maykop art and Sumerian art (Rostovtzeff 1920) and M. V. Andreeva described this phenomenon within a broad cultural and chronological context (Andreeva 1977), new archaeological studies have only extended this picture of a vast cultural province that appeared between the Caucasus and the northern fringe of Western Asia (Trifonov 1987). The discovery of the Leyla-Tepe culture (Narimanov 1987) and Maykop-type kurgans in Azerbaijan (Lyonnet et al. 2008) and adjacent Iran (Muscarella 1969, 1971, 2003; Trifonov 2000) has confirmed the spatial and temporal unity of this phenomenon as a precondition for free circulation of cultural patterns and technical innovations across vast areas of the Caucasus and Western Asia. Jewellery made of gemstones and precious metals, primarily gold, was probably one such innovation.

Attempts to demarcate the historical region where the Maykop culture emerged and developed have emphasized the role of Upper Mesopotamia in the development of the Sumerian civilization and the definition of a northern centre of urbanization, independent from the centres of the south (Rothman 2002; Oats et al. 2007). The turn of the fourth millennium BC saw the development of various cultural traditions in south-east Anatolia, north-east Syria and north-west Iran; on the northern fringe, these traditions manifested themselves through the Maykop culture. Perhaps it is no coincidence that the first high-status burials containing gold and gemstone jewellery (including carnelian, turquoise and lapis lazuli) appear in these northern, rather than southern, centres in the first quarter of 4000 BC (e.g., Tepe Gawra, graves 109, 110) (Piasnall 2002). With regard to funeral rites and stylistic characteristics of jewellery pieces, these graves have many parallels with early Maykop burials (Munchaev 1975, 329; Trifonov 1987, 20).

It still remains unclear if the goldsmiths of Upper Mesopotamia mastered the technique of making thin-walled jointless beads. The gold beads from Tepe Gawra are described as spherical or ball-shaped, but their maximum diameter (5–8mm) always exceeds the length of the bore (3–4mm) (Tobler 1950, 89, 199, pl. LV, a). On the whole, these measurements are consistent with the proportions and sizes of some Maykop beads.(…)

It is quite possible that a distinctive technique of making thin-walled jointless beads from gold was a regional technological development of Maykop culture goldsmiths, within a wider tradition of Near East metalwork, as a type of production regulated by ritual beliefs (Gell 1992; Benzel 2013).

These deep-rooted Near East traditions of ritualization of the production and use of jewellery pieces made of gold, silver and gemstones in the Maykop culture, on the one hand, maintained familiar canons of ritual behaviour and, on the other, made perception of sophisticated symbolism of gemstones more difficult for neighbouring cultures with different living standards, levels of social development and value systems to understand. The jewellery traditions of the Maykop culture had no successors in the Caucasus or the adjacent steppes. In the third millennium BC, the goldsmiths of Europe and Asia had to reinvent the technique of making thin-walled jointless gold beads from scratch (Born et al. 2009).

Also interesting is Holocene environmental history and populating of mountainous Dagestan (Eastern Caucasus, Russia), by Ryabogina et al., Quaternary International (2018).

The combination of Holocene environment changes and the settlement of the territory of Dagestan.

Related excerpts, about the climate of an adjacent region of the Caucasus before, during, and after the Maykop culture:

The 7th millennium BC featured a warm and arid climate, so that time corresponds to the steppe landscapes in the final stage of the Mesolithic. It is likely that the formation of a producing economy in the mountainous zone of Dagestan gradually emerged against this background. In the Neolithic period, the area remained almost treeless, as it was still warm and quite dry. However, archaeological data indicates that long-term settlements with well-developed farming spread in the mountainous zone around 6200-5500 BC.

The beginning of increasing humidity and the appearance of deciduous forests corresponds to the early Chalcolithic period of the Eastern Caucasus. It is the most poorly studied period in the history of this region. Covering a time span of 2000 years, this period was the least saturated by archaeological sites. At the start of this period, only the stands of herdsman in the mountain zone are known, dating to the second half of the 6th millennium BC (Gadgiev, 1991). It is still not clear whether the mountains were not settled in such a favorable climatic stage. The uncertainty may be due to the fact that people have chosen other ecological niches, or it could be we simply do not have data due to the insufficient archaeological survey of the territory. It is surprising that the turn to drier climate and the reduction of deciduous forests in the inner mountainous part of Dagestan, the large, long-term settlements like Ginchi emerge with pronounced specialization in agriculture (Fig. 7 panel (2)) (Gadgiev, 1991).

After the dry climate, simultaneously with cooling, the subsequent spread of pine forests coincides with the beginning of expansion of Kura-Araxes culture from the territory of Georgia through Chechnya to the mountainous Dagestan. Debates on the impact of past climate on Kura-Araxes societies in Transcaucasus have a long history (for the comprehensive review see, for example, Connor and Kvavadze, 2014 and references therein). In general, it is clear that after 3000 BC, forest cover in most areas of the Kura-Araxes region in the Transcaucasia reached its maximum extent in the Holocene (Connor and Kvavadze, 2014). However, at the same time lakes in Central Anatolia began to dry out and Caspian Sea levels fell (Roberts et al. 2011; Leroy et al. 2013), and arid conditions were identified in mountainous Dagestan in the 4th millennium. Clearly the regional moisture balance shifted in the Eastern Caucasus only in the late 4th to early 3rd millennium BC (this study). The only available radiocarbon dating of Dagestan confirms that the agricultural settlements of the Early Bronze Age appear not in the middle of the 4th millennium BC, but in the early 3rd millennium BC; that is not earlier than the stage of increasing moistening and the appearance of pine forests.

See also:

Steppe and Caucasus Eneolithic: the new keystones of the EHG-CHG-ANE ancestry in steppe groups


Some interesting excerpts from Wang et al. (2018):

An interesting observation is that steppe zone individuals directly north of the Caucasus (Eneolithic Samara and Eneolithic steppe) had initially not received any gene flow from Anatolian farmers. Instead, the ancestry profile in Eneolithic steppe individuals shows an even mixture of EHG and CHG ancestry, which argues for an effective cultural and genetic border between the contemporaneous Eneolithic populations in the North Caucasus, notably Steppe and Caucasus. Due to the temporal limitations of our dataset, we currently cannot determine whether this ancestry is stemming from an existing natural genetic gradient running from EHG far to the north to CHG/Iran in the south or whether this is the result of farmers with Iranian farmer/ CHG-related ancestry reaching the steppe zone independent of and prior to a stream of Anatolian farmer-like ancestry, where they mixed with local hunter-gatherers that carried only EHG ancestry.

Image modified from Wang et al. (2018). Samples projected in PCA of 84 modern-day West Eurasian populations (open symbols). Previously known clusters have been marked and referenced. An Eastern European (blue) and a Caucasus (brown) ‘clouds’ have been drawn in dotted circles, leaving Pontic-Caspian steppe and derived groups between them.See the original file here.

Concerning the influences from the south, our oldest dates from the immediate Maykop predecessors Darkveti-Meshoko (Eneolithic Caucasus) indicate that the Caucasus genetic profile was present north of the range ~6500 BP, 4500 calBCE. This is in accordance with the Neolithization of the Caucasus, which had started in the flood plains of the great rivers in the South Caucasus in the 6th millennium BCE from where it spread to the West and Northwest Caucasus during the 5th millennium BCE9, 49. It remains unclear whether the local CHG ancestry profile (represented by Late Upper Palaeolithic/Mesolithic individuals from Kotias Klde and Satsurblia in today’s Georgia) was also present in the North Caucasus region before the Neolithic. However, if we take the Caucasus hunter-gatherer individuals from Georgia as a local baseline and the oldest Eneolithic Caucasus individuals from our transect as a proxy for the local Late Neolithic ancestry, we notice a substantial increase in Anatolian farmer-related ancestry. This in all likelihood is linked to the process of Neolithization, which also brought this type of ancestry to Europe. As a consequence, it is possible that Neolithic groups could have reached the northern flanks of the Caucasus earlier50 (Supplementary Information 1) and in contact with local hunter gatherers facilitated the exploration of the steppe environment for pastoralist economies. Hence, additional sampling from older individuals is needed to fill this temporal and spatial gap.

The newest paper of the Reich/Jena group has brought samples (probably) much nearer to the actual CHG and ANE contribution seen in Eneolithic steppe peoples than the previously available Kotias Klde, Satsurblia, Afontova Gora 3, or Mal’ta.

It is impossible to say without direct access to the samples, but it is very likely that we will soon be able to break down different gross contributions from groups similar to these Steppe/Caucasus Neolithic ancestral groups into the diverse Eneolithic cultures of the Pontic-Caspian steppe, and thus trace more precisely each of these cultures to their genetic (and thus ethnolinguistic) heirs.

Admixture Graph modelling of the population history of the Caucasus region. We started with a skeleton tree without admixture including Mbuti, Loschbour and MA1. We grafted onto this EHG, CHG, Globular_Amphora, Eneolithic_steppe, Maykop, and Yamnaya_Caucasus, adding them consecutively to all possible edges in the tree and retaining only graph solutions that provided no differences of |Z|>3 between fitted and estimated statistics. The worst match is |Z|=2.824 for this graph. We note that the maximum discrepancy is f4(MA1, Maykop; EHG, Eneolithic_steppe) = -3.369 if we do not add the 4% EHG ancestry to Maykop. Drifts along edges are multiplied by 1000 and dashed lines represent admixture.”

Some more representative samples from Eneolithic steppe, steppe-forest and forest zone cultures of Eastern Europe will probably help with the fine-scale structure of different Chalcolithic groups, especially the homeland of early Corded Ware groups.

These new samples seem another good reason (like the Botai and R1b-M73) to rethink the role of (what I assumed were) different westward Mesolithic Eurasian waves of expansion influencing the formation of an Indo-Uralic and Indo-European community in Eastern Europe, and return to the simpler idea of local contributions from North Caucasus and steppe peoples absorbed by expanding EHG-like groups.


No large-scale steppe migration into Anatolia; early Yamna migrations and MLBA brought LPIE dialects in Asia


Another, simultaneous paper with the Eurasian samples from Nature, The first horse herders and the impact of early Bronze Age steppe expansions into Asia, by de Barros Damgaard et al., Science (2018).

A lot of interesting data, I will try to analyse its main implications, if only superficially, in sections.

Anatolian samples

Anatolia_EBA from Ovaören, and Anatolia_MLBA (this including Assyrian and Old Hittite samples), all from Kalehöyük, show almost no change in Y-DNA lineages (three samples J2a, one G2a), and therefore an origin of these people in common with CHG and Iranian Neolithic populations is likely. No EHG ancestry is found. And PCA cluster is just somehow closer to Europe, but not to EHG populations.

NOTE. Hittite is attested only in the late first half of the 2nd millennium, although the authors cite (in the linguistic supplement) potential evidence from the palatial archives of the ancient city of Ebla in Syria to argue that Indo-European languages may have been already spoken in the region in the late 3rd millennium BCE.

Regarding the Assyrian samples (one J2a) from Ovaören:

Layer V of GT-137 was the richest in terms of architectural finds and dates to the Early Bronze Age II. In this layer, 2 different structures and a well were uncovered. The well was filled with stones, pottery, and human skeletons (Figs. S2 and S3). In total, skeletons belonging to 22 individuals, including adults, young adults, and children, must belong to the disturbed Early Bronze Age II graves adjacent to the well (103). Pottery and stones found below the skeletons demonstrate that the water well was consciously filled and closed. The fill consists of dumped stones, sherds and skeletons, and the closing stones demonstrate that the water well was consciously filled and cancelled.

Regarding the site most likely associated with the emergence of Old Hittite (two samples J2a1, one G2a2b1), this is what we know:

The Middle Bronze Age at Kaman-Kalehöyük represented by stratum IIIc yields material remains (seals and ceramics) contemporary with the international trade system managed by expatriate Assyrian merchants evidenced at the nearby site of Kültepe/Kanesh. It is therefore also referred to as belonging to the “Assyrian Colony Period” (98). The stratum has revealed three burned architectural units, and it has been suggested that the seemingly site-wide conflagration might be connected to a destruction event linked with the emergence of the Old Hittite state (99). (…) Omura (100) suggests that the rooms could belong to a public building, and that it might even be a small trade center based on the types of artifacts recovered. Omura (100) has concluded that the evidence from the first complex indicates a battle between 2 groups took place at the site. It is possible that a group died inside the buildings, mostly perishing in the fire, while another group died in the courtyard.

NOTE. For more on the Old Hittite period, you can read this for example.

Regarding PCA:

The PCA (Fig. 2B) indicates that all the Anatolian genome sequences from the Early Bronze Age ( -2200 BCE) and Late Bronze Age (-1600 BCE) cluster with a previously sequenced Copper Age ( -3900- 3700 BCE) individual from Northwestern Anatolia and lie between Anatolian Neolithic (Anatolia_ N) samples and CHG samples but not between Anatolia_N and EHG samples.

(…) we are not able to reject a two-population qpAdm model in which these groups derive -60% of their ancestry from Anatolian farmers and -40% from CHG-related ancestry (p-value = 0.5). This signal is not driven by Neolithic Iranian ancestry.

Principal Component Analysis estimated with ancient and modern Eurasians.

NOTE. Anatolian Iron Age samples, from the Hellenistic period, which was obviously greatly influenced by different, later Indo-European migrations, does show a change in PCA.

Regarding CHG ancestry:

Ancient DNA findings suggest extensive population contact between the Caucasus and the steppe during the Copper Age (-5000-3000 BCE) (1, 2, 42). Particularly, the first identified presence of Caucasian genomic ancestry in steppe populations is through the Khvalynsk burials (2, 47) and that of steppe ancestry in the Caucasus is through Armenian Copper Age individuals (42). These admixture processes likely gave rise to the ancestry that later became typical of the Yamnaya pastoralists (7), whose IE language may have evolved under the influence of a Caucasian language, possibly ‘from the Maykop culture (50, 55). This scenario is consistent with both the “Copper Age steppe” (4) and the “Caucasian” models for the origin of the Proto-Anatolian language (56).

The CHG specific ancestry and the absence of EHG-related ancestry in Bronze Age Anatolia would be in accordance with intense cultural interactions between populations in the Caucasus and Anatolia observed during the late 5th millennium BCE that seem to come to an end in the first half of the 4th millennium BCE with the village-based egalitarian Kura-Araxes society (59, 60), thus preceding the emergence and dispersal of Proto-Anatolian.

Our results indicate that the early spread of IE languages into Anatolia was not associated with any large-scale steppe-related migration, as previously suggested (61). Additionally, and in agreement with the later historical record of the region (62), we find no correlation between genetic ancestry and exclusive ethnic or political identities among the populations of Bronze Age Central Anatolia, as has previously been hypothesized ( 63).

The Anatolian question

There is no steppe ancestry or R1b-M269 lineages near early historic Hittites. Yet.

Nevertheless, we already know about potentially similar cases:

So there seems to be thus no theoretical problem in accepting:

  • That neither steppe ancestry nor R1b-M269 subclades, already diminished in Bulgaria in the mid-5th millennium, did reach Anatolia, but only those Common Anatolian-speaking Aegean groups over whose ancestors Proto-Anatolians (marked by incoming EHG ancestry) would have previously dominated in the Balkans.
  • That steppe ancestry and R1b-M269 subclades did in fact arrive in the Aegean, but EHG was further diluted among the CHG-related population by the time of the historic Anatolian-speaking peoples in central Anatolia. Or, the most likely option, that their trace have not been yet found. Probably the western Luwian peoples, near Troy, were genetically closer to Common Anatolians.

Both of these scenarios are interesting, in that they show potential links between Pre-Greek peoples of Hellas (related to Anatolians) and the Pelasgian substrate of early Greek dialects, since they show a similar recent CHG-related wave from the East.

What we can assert right now is that Proto-Anatolian must have separated quite early for this kind of data to show up. This should mean an end to the Late PIE origin of Anatolian, if there was some lost soul from the mid-20th century still rooting for this.

As I said in my review of Lazaridis’ latest preprint, we will have to wait for the appropriate potential routes of expansion of Proto-Anatolian to be investigated. As he answered, the lack of EHG poses a problem for steppe expansion into Anatolia, but there is still no better alternative model proposed.

Model-based clustering analysis of present-day and ancient individuals assuming K = 6 ancestral components. The main ancestry components at K = 6 correlate well with CHG (turquoise), a major component of Iran_N, Namazga_CA and South Asian dines; EHG (pale blue), a component of the steppe dine and present in South Asia; East Asia (yellow ochre), the other component of the steppe d ine also in Tibeto-Burman South Asian populations; South Indian (pink), a core component of South Asian populations; Anatolian_N (purple), an important component of Anatolian Bronze Age and Steppe_MLBA; Onge (dark pink) forms its own component.

This is what the authors have to say:

Our findings are thus consistent with historical models of cultural hybridity and “Middle Ground” in a multi-cultural and multi-lingual but genetically homogeneous Bronze Age Anatolia (68, 69). Current linguistic estimations converge on dating the Proto-Anatolian split from residual PIE to the late 5th or early 4th millennia BCE (58, 70) and place the breakup of Anatolian IE inside Turkey prior to the mid-3rd millennium (53, 71,72).

We cannot at this point reject a scenario in which the introduction of the Anatolian IE languages into Anatolia was coupled with the CHG-derived admixture prior to 3700 BCE, but note that this is contrary to the standard view that PIE arose in the steppe north of the Caucasus (4) and that CHG ancestry is also associated with several non-IE-speaking groups, historical and current. Indeed, our data are also consistent with the first speakers of Anatolian IE coming to the region by way of commercial contacts and small-scale movement during the Bronze Age. Among comparative linguists, a Balkan route for the introduction of Anatolian IE is generally considered more likely than a passage through the Caucasus, due, for example, to greater Anatolian IE presence and language diversity in the west (73). Further discussion of these options is given in the archaeological and linguistic supplementary discussions (48, 49).

If you are asking yourselves why the Danish school (of Allentoft, Kristiansen, and Kroonen, co-authors of this paper) was not so fast to explain the findings the same way the proposed their infamous Indo-European – steppe ancestry association (i.e. ancestry = language, ergo CHG = PIE in this case), and resorted to mainstream anthropological models instead to explain the incongruence, I can think of two main reasons:

The possibility of having an early PIE around the Caucasus, potentially closely related not only to Uralic to the north, but also to Caucasian languages, Sumerian, Afroasiatic, Elamo-Dravidian, etc. could be a good reason for those excited with these few samples to begin dealing with macro-language proposals, such as Eurasiatic and Nostratic. If demonstrated to be true, a Northern Iranian origin of Middle PIE would also help relieve a little bit the pressure that some are feeling about the potentially male-driven Indo-European continuity (even if not “autochthonous”) associated with the expansion of R1b-L23 subclades.

On the other hand, I am a firm supporter of solid anthropological models of migration, and of “late and small” language expansions, usually accompanied by demic diffusion, which has been demonstrated to be linked with haplogroup expansion and reduction in variability.

Therefore, for the moment, even if it is weak – as weak as it always was (but still stronger than Gimbutas’ Maykop route) – the Balkan route seems like the best fit for all the data combined.

In fact, we already have steppe ancestry moving into the Lower Danube and Bulgaria in the mid-5th millennium. Let’s not forget that.

Yamna expansion to the East

Interesting data from an early East Yamna offshoot at Karagash, ca. 3018-2887 BC, of R1b-Z2106 lineage, which shows some ancestry, lineage, and cultural continuity in Sholpan, ca. 2620-2468 BC, in Kazakhstan.

This sample might be part of another descendant group from the migration waves that reached Afanasevo, and can thus be related to other early Asian R1b-L23 samples found in Narasimhan et al. (2018).

On the formation of Yamna and its CHG contribution, from the supplementary material:

  1. An admixture event, where Yamnaya is formed from a CHG population related to KK1 [=Kotias, dated ca. 7800 BC] and an ANE population related to Sidelkino and Botai. We inferred 54% of the Yamnaya ancestry to come from CHG and the remaining 46% to come from ANE.
  2. A split event, where the CHG component of Yamnaya splits from KK1. The model inferred this time at 27 kya (though we note the larger models in Sections S2.12.4 and S2.12.5 inferred a more recent split time [see below graphic]).
  3. A split event, where the ANE component of Yamnaya splits from Sidelkino. This was inferred at about about 11 kya.
  4. A split event, where the ANE component of Yamnaya splits from Botai. We inferred this to occur 17 kya. Note that this is above the Sidelkino split time, so our model infers Yamnaya to be more closely related to the EHG Sidelkino, as expected.
  5. An ancestral split event between the CHG and ANE ancestral populations. This was inferred to occur around 40 kya.
A 10-leaf model based on combining the models in Fig. S16 and Fig. S19 and re-estimating the model parameters.

On the expansion of domestication

CHG is not found in Botai, no gene flow from Yamna is found in its samples, and they are more related to East Asians, while Yamna is related to West Eurasians:

The lack of evidence of admixture between Botai horse herders and western steppe pastoralists is consistent with these latter migrating through the central steppe but not settling until they reached the Altai to the east (4). More significantly, this lack of admixture suggests that horses were domesticated by hunter-gatherers not previously familiar with farming, as were the cases for dogs (38) and reindeer (39). Domestication of the horse thus may best parallel that of the reindeer, a food animal that can be milked and ridden, which has been proposed to be domesticated by hunters via the “prey path” (40); indeed anthropologists note similarities in cosmological beliefs between hunters and reindeer herders (41). In contrast, most animal domestications were achieved by settled agriculturalists (5).

NOTE. I am not sure, but they seem to hint that there were separate events of horse domestication and horse-riding technique by the Botai and Yamna populations due to their lack of genetic contribution from the latter to the former. I guess they did not take into account farming spreading to the steppe without genetic contribution beyond the Dnieper… In fact, the superiority in horse-riding shown by the expanding Yamna peoples – as they state – should also serve to suggest from where the original technique expanded.

Indo-Iranian migrations

On the expansion of Yamna, and the different expansion of Steppe MLBA (with Indo-Iranian speakers) into Asia, further supporting Narasimhan et al. (2018), they have this to say:

However, direct influence of Yamnaya or related cultures of that period is not visible in the archaeological record, except perhaps for a single burial mound in Sarazm in present-day Tajikistan of contested age (44, 45). Additionally, linguistic reconstruction of proto-culture coupled with the archaeological chronology evidences a Late (-2300-1200 BCE) rather than Early Bronze Age (-3000-2500 BCE) arrival of the Indo-Iranian languages into South Asia (16, 45, 46). Thus, debate persists as to how and when Western Eurasian genetic signatures and IE languages reached South Asia.

Samples from the Namazga region (current Turkmenistan) from the Iron Age show an obvious influence from steppe MLBA (ca. 2300-1200 BC), and not steppe EBA (i.e. Yamna), population, in contrast with samples from the Chalcolithic (ca. 3300 BC), which don’t show this influence. This helps distinguish prior contacts with Iran Neolithic from the actual steppe population that expanded Indo-Iranian into Asia.

Very interesting therefore the Namazga CA sample (ca. 855 BC), of R1a-Z93 subclade, showing the sign of immigrant Indo-Aryans in the region. For more on this we will need an evaluation in common with the corrected data from Narasimhan et al. (2018), and all, including de Barros (Nature 2018), in combination with statistical methods to ascertain differences between early Indo-Aryans and Iranians.

A summary of the four qpAdm models fitted for South Asian populations. For each modern South Asian population. we fit different models with qpAdm to explain their ancestry composition using ancient groups and present the f irst model that we could not reject in the following priority order: 1. Namazga_CA + Onge, 2. Namazga_CA + Onge + Late Bronze Age Steppe, 3. Namazga_CA + Onge + Xiongnu_lA (East Asian proxy). and 4. Turkmenistan_lA + Xiongnu_lA. Xiongnu_lA were used here to represent East Asian ancestry. We observe that while South Asian Dravidian speakers can be modeled as a mixture of Onge and Namazga_CA. an additional source related to Late Bronze Age steppe groups is required for IE speakers. In Tibeto-Burman and Austro-Asiatic speakers. an East Asian rather than a Steppe_MLBA source is required.

Siberian peoples and N1c lineages

We have already seen how the paper on Eurasian steppe samples tries to assign Uralic to Neolithic peoples east of the Urals. The association with Okunevo is unlikely, since most are of haplogroup Q1a2, but they seem to suggest (combining both papers) that they accompanied N lineages from Siberian hunter-gatherers (present e.g. in Botai or Shamanka II, during the Early Neolithic), and formed part of (or suffered from) different demic diffusion waves:

These serial changes in the Baikal populations are reflected in Y-chromosome lineages (Fig. SA; figs. S24 to S27, and tables S13 and SI4). MAI carries the R haplogroup, whereas the majority of Baikal_EN males belong to N lineages, which were widely distributed across Northern Eurasia (29), and the Baikal_LNBA males all carry Q haplogroups, as do most of the Okunevo_EMBA as well as some present-day Central Asians and Siberians.

NOTE. Also interesting to see no R1a in Baikal hunter-gatherers after ca. 3500 BC, and a prevalence of N lineages as supported in a previous paper on the Kitoi culture, which some had questioned in the past.

In fact, the only N1c1 sample comes from Ust’Ida Late Neolithic, 180km to the north of Lake Baikal, apparently before the expansion of Q1a2a lineages during the EBA period. While this sample may be related to those expanded later in Finno-Ugric territory (although it may only be related to those expanded much later with Yakuts), other samples are not clearly from those found widely distributed among North-East Europeans only after the Iron Age, or – as in the case of Shamanka II (N1c2), they are clearly not of the same haplogroup.

Geographical location of ancient samples belonging to major clade N of the Y-chromosome.


It is great to see the paper and the supplementary material deal with Y-DNA haplogroups and their relevance for migrations with such detail. Especially because this paper comes from the same Copenhagen-based research group that originally associated ancestry with language, creating thus today’s mess based on steppe ancestry.

Regarding Y-DNA data, once again almost 100% of samples from late Khvalynsk/Yamna and derived cultures (like Afanasevo and Bell Beaker) are R1b-L23, no single R1a-M417 lineage found, and few expected by now, if any, within Late Proto-Indo-European territory.

While they claim to take Y-DNA into account to assess migrations – as they do for example with Asian cultures – , their previous model of a Yamna “R1a-R1b community” remains oddly unchanged, and they even insist on it in the supplementary materials, as they do in their parallel Nature paper.

They have also expressly mitigated the use of ancestral components to assess populations, citing the ancestral and modern association of CHG ancestry with different ethnolinguistic groups in the Middle East, to dismiss any rushed conclusions on the origin of Anatolian, and consequently of Middle PIE. And they did so evidently because it did not fit the anthropological data that is mainstream today (supporting a Balkan route), which is the right thing to do.

However, they have apparently not stopped to reconsider the links of CWC and steppe ancestry to ancestral and modern Uralic peoples – although they expressly mention the strong connection with modern Karelians in the supplementary material.

Also, after Narasimhan et al. (2018), there is a clear genetic continuity with East Yamna (in ancestry as in R1b-L23 subclades), so their interpretations about Indo-Iranian in this paper and especially de Barros (Nature 2018) – regarding the Abashevo -> Sintashta/Srunba/Andronovo connection – come, again, too late.