N1c-L392 associated with expanding Turkic lineages in Siberia


Second in popularity for the expansion of haplogroup N1a-L392 (ca. 4400 BC) is, apparently, the association with Turkic, and by extension with Micro-Altaic, after the Uralic link preferred in Europe; at least among certain eastern researchers.

New paper in a recently created journal, by the same main author of the group proposing that Scythians of hg. N1c were Turkic speakers: On the origins of the Sakhas’ paternal lineages: Reconciliation of population genetic / ancient DNA data, archaeological findings and historical narratives, by Tikhonov, Gurkan, Demirdov, and Beyoglu, Siberian Research (2019).

Interesting excerpts:

According to the views of a number of authoritative researchers, the Yakut ethnos was formed in the territory of Yakutia as a result of the mixing of people from the south and the autochthonous population [34].

These three major Sakha paternal lineages may have also arrived in Yakutia at different times and/ or from different places and/or with a difference in several generations instead, or perhaps Y-chromosomal STR mutations may have taken place in situ in Yakutia. Nevertheless, the immediate common ancestor(s) from the Asian Steppe of these three most prevalent Sakha Y-chromosomal STR haplotypes possibly lived during the prominence of the Turkic Khaganates, hence the near-perfect matches observed across a wide range of Eurasian geography, including as far as from Cyprus in the West to Liaoning, China in the East, then Middle Lena in the North and Afghanistan in the South (Table 3 and Figure 5). There may also be haplotypes closely-related to ‘the dominant Elley line’ among Karakalpaks, Uzbeks and Tajiks, however, limitations in the loci coverage for the available dataset (only eight Y-chromosomal STR loci) precludes further conclusions on this matter [25].

17-loci median-joining network analysis of the original/dominant Elley, Unknown and Omogoy Y-chromosomal STR haplotypes with the YHRD matches from outside Yakutia populations.

According to the results presented here, very similar Y-STR haplotypes to that of the original Elley line were found in the west: Afghanistan and northern Cyprus, and in the east: Liaoning Province, China and Ulaanbaator, Northern Mongolia. In the case of the dominant Omogoy line, very closely matching haplotypes differing by a single mutational step were found in the city of Chifen of the Jirin Province, China. The widest range of similar haplotypes was found for the Yakut haplotype Unknown: In Mongolia, China and South Korea. For instance, haplotypes differing by a single step mutation were found in Northern Mongolia (Khalk, Darhad, Uryankhai populations), Ulaanbaator (Khalk) and in the province of Jirin, China (Han population).

14-loci median-joining network analysis for the original/dominant Elley (Ell), Unknown Clan
(Vil), Omogoy (Omo), Eurasian (Eur) and Xiongnu (Xuo) Y-chromosomal STR haplotypes and that for a representative ancient DNA sample (Ch0 or DSQ04) from the Upper Xiajiadian Culture
recovered from the Inner Mongolia Autonomous Region, China.

Notably, Tat-C-bearing Y-chromosomes were also observed in ancient DNA samples from the 2700-3000 years-old Upper Xiajiadian culture in Inner Mongolia, as well as those from the Serteya II site at the Upper Dvina region in Russia and the ‘Devichyi gory’ culture of long barrow burials at the Nevel’sky district of Pskovsky region in Russia. A 14-loci Y-chromosomal STR median-joining network of the most prevalent Sakha haplotypes and a Tat-C-bearing haplotype from one of the ancient DNA samples recovered from the Upper Xiajiadian culture in Inner Mongolia (DSQ04) revealed that the contemporary Sakha haplotype ‘Xuo’ (Table 2, Haplotype ID “Xuo”) classified as that of ‘the Xiongnu clan’ in our current study, was the closest to the ancient Xiongnu haplotype (Figure 6). TMRCA estimate for this 14-loci Y-chromosomal STR network was 4357 ± 1038 years or 2341 ± 1038 BCE, which correlated well with the Upper Xiajiadian culture that was dated to the Late Bronze Age (700-1000 BCE).

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

NOTE. Also interesting from the paper seems to be the proportion of E1b1b among admixed Russian populations, in a proportion similar to R1a or I2a(xI2a1).

It is tempting to associate the prevalent presence of N1c-L392 in ancient Siberian populations with the expansion of Altaic, by simplistically linking the findings (in chronological order) near Lake Baikal (Damgaard et al. 2018), Upper Xiajiadian (Cui et al. 2013), among Khövsgöl (Jeong et al. 2018), in Huns (Damgaard et al. 2018), and in Mongolic-speaking Avars (Csáky et al. 2019).

However, its finding among Palaeo-Laplandic peoples in the Kola peninsula ca. 1500 BC (Lamnidis et al. 2018) and among Palaeo-Siberian populations near the Yana River (Sikora et al. 2018) ca. AD 1200 should be enough to accept the hypothesis of ancestral waves of expansion of the haplogroup over northern Eurasia, with acculturation and further expansions in the different regions since the Iron Age (see more on its potential expansion waves).

Also, a simple look at the TMRCA and modern distribution was enough to hypothesize long ago the lack of connection of N1c-L392 with Altaic or Uralic peoples. From Ilumäe et al. (2016):

Previous research has shown that Y chromosomes of the Turkic-speaking Yakuts (Sakha) belong overwhelmingly to hg N3 (formerly N1c1). We found that nearly all of the more than 150 genotyped Yakut N3 Y chromosomes belong to the N3a2-M2118 clade, just as in the Turkic-speaking Dolgans and the linguistically distant Tungusic-speaking Evenks and Evens living in Yakutia (Table S2). Hence, the N3a2 patrilineage is a prime example of a male population of broad central Siberian ancestry that is not intrinsic to any linguistically defined group of people. Moreover, the deepest branch of hg N3a2 is represented by a Lebanese and a Chinese sample. This finding agrees with the sequence data from Hallast et al., where one Turkish Y chromosome was also assigned to the same sub-clade. Interestingly, N3a2 was also found in one Bhutan individual who represents a separate sub-lineage in the clade. These findings show that although N3a2 reflects a recent strong founder effect primarily in central Siberia (Yakutia, Sakha), the sub-clade has a much wider distribution area with incidental occurrences in the Near East and South Asia.

Frequency-Distribution Maps of Individual Sub-clades of hg N3a2, by Ilumäe et al. (2016).

The most striking aspect of the phylogeography of hg N is the spread of the N3a3’6-CTS6967 lineages. Considering the three geographically most distant populations in our study—Chukchi, Buryats, and Lithuanians—it is remarkable to find that about half of the Y chromosome pool of each consists of hg N3 and that they share the same sub-clade N3a3’6. The fractionation of N3a3’6 into the four sub-clades that cover such an extraordinarily wide area occurred in the mid-Holocene, about 5.0 kya (95% CI = 4.4–5.7 kya). It is hard to pinpoint the precise region where the split of these lineages occurred. It could have happened somewhere in the middle of their geographic spread around the Urals or further east in West Siberia, where current regional diversity of hg N sub-lineages is the highest (Figure 1B). Yet, it is evident that the spread of the newly arisen sub-clades of N3a3’6 in opposing directions happened very quickly. Today, it unites the East Baltic, East Fennoscandia, Buryatia, Mongolia, and Chukotka-Kamchatka (Beringian) Eurasian regions, which are separated from each other by approximately 5,000–6,700 km by air. N3a3’6 has high frequencies in the patrilineal pools of populations belonging to the Altaic, Uralic, several Indo-European, and Chukotko-Kamchatkan language families. There is no generally agreed, time-resolved linguistic tree that unites these linguistic phyla. Yet, their split is almost certainly at least several millennia older than the rather recent expansion signal of the N3a3’6 sub-clade, suggesting that its spread had little to do with linguistic affinities of men carrying the N3a3’6 lineages.

Frequency-Distribution Maps of Individual Subclade N3a3 / N1a1a1a1a1a-CTS2929/VL29.

It was thus clear long ago that N1c-L392 lineages must have expanded explosively in the 5th millennium through Northern Eurasia, probably from a region to the north of Lake Baikal, and that this expansion – and succeeding ones through Northern Eurasia – may not be associated to any known language group until well into the common era.


“Steppe ancestry” step by step: Khvalynsk, Sredni Stog, Repin, Yamna, Corded Ware


Wang et al. (2018) is obviously a game changer in many aspects. I have already written about the upcoming Yamna Hungary samples, about the new Steppe_Eneolithic and Caucasus Eneolithic keystones, and about the upcoming Greece Neolithic samples with steppe ancestry.

An interesting aspect of the paper, hidden among so many relevant details, is a clearer picture of how the so-called Yamnaya or steppe ancestry evolved from Samara hunter-gatherers to Yamna nomadic pastoralists, and how this ancestry appeared among Proto-Corded Ware populations.

Image modified from Wang et al. (2018). Marked are in orange: equivalent Steppe_Maykop ADMIXTURE; in red, approximate limit of Anatolia_Neolithic ancestry found in Yamna populations; in blue, Corded Ware-related groups. “Modelling results for the Steppe and Caucasus cluster. Admixture proportions based on (temporally and geographically) distal and proximal models, showing additional Anatolian farmer-related ancestry in Steppe groups as well as additional gene flow from the south in some of the Steppe groups as well as the Caucasus groups.”

Please note: arrows of “ancestry movement” in the following PCAs do not necessarily represent physical population movements, or even ethnolinguistic change. To avoid misinterpretations, I have depicted arrows with Y-DNA haplogroup migrations to represent the most likely true ethnolinguistic movements. Admixture graphics shown are from Wang et al. (2018), and also (the K12) from Mathieson et al. (2018).

1. Samara to Early Khvalynsk

The so-called steppe ancestry was born during the Khvalynsk expansion through the steppes, probably through exogamy of expanding elite clans (eventually all R1b-M269 lineages) originally of Samara_HG ancestry. The nearest group to the ANE-like ghost population with which Samara hunter-gatherers admixed is represented by the Steppe_Eneolithic / Steppe_Maykop cluster (from the Northern Caucasus Piedmont).

Steppe_Eneolithic samples, of R1b1 lineages, are probably expanded Khvalynsk peoples, showing thus a proximate ancestry of an Early Eneolithic ghost population of the Northern Caucasus. Steppe_Maykop samples represent a later replacement of this Steppe_Eneolithic population – and/or a similar population with further contribution of ANE-like ancestry – in the area some 1,000 years later.


This is what Steppe_Maykop looks like, different from Steppe_Eneolithic:


NOTE. This admixture shows how different Steppe_Maykop is from Steppe_Eneolithic, but in the different supervised ADMIXTURE graphics below Maykop_Eneolithic is roughly equivalent to Eneolithic_Steppe (see orange arrow in ADMIXTURE graphic above). This is useful for a simplified analysis, but actual differences between Khvalynsk, Sredni Stog, Afanasevo, Yamna and Corded Ware are probably underestimated in the analyses below, and will become clearer in the future when more ancestral hunter-gatherer populations are added to the analysis.

2. Early Khvalynsk expansion

We have direct data of Khvalynsk-Novodanilovka-like populations thanks to Khvalynsk and Steppe_Eneolithic samples (although I’ve used the latter above to represent the ghost Caucasus population with which Samara_HG admixed).

We also have indirect data. First, there is the PCA with outliers:


Second, we have data from north Pontic Ukraine_Eneolithic samples (see next section).

Third, there is the continuity of late Repin / Afanasevo with Steppe_Eneolithic (see below).

3. Proto-Corded Ware expansion

It is unclear if R1a-M459 subclades were continuously in the steppe and resurged after the Khvalynsk expansion, or (the most likely option) they came from the forested region of the Upper Dnieper area, possibly from previous expansions there with hunter-gatherer pottery.

Supporting the latter is the millennia-long continuity of R1b-V88 and I2a2 subclades in the north Pontic Mesolithic, Neolithic, and Early Eneolithic Sredni Stog culture, until ca. 4500 BC (and even later, during the second half).

Only at the end of the Early Eneolithic with the disappearance of Novodanilovka (and beginning of the steppe ‘hiatus’ of Rassamakin) is R1a to be found in Ukraine again (after disappearing from the record some 2,000 years earlier), related to complex population movements in the north Pontic area.

NOTE. In the PCA, a tentative position of Novodanilovka closer to Anatolia_Neolithic / Dzudzuana ancestry is selected, based on the apparent cline formed by Ukraine_Eneolithic samples, and on the position and ancestry of Sredni Stog, Yamna, and Corded Ware later. A good alternative would be to place Novodanilovka still closer to the Balkan outliers (i.e. Suvorovo), and a source closer to EHG as the ancestry driven by the migration of R1a-M417.


The first sample with steppe ancestry appears only after 4250 BC in the forest-steppe, centuries after the samples with steppe ancestry from the Northern Caucasus and the Balkans, which points to exogamy of expanding R1a-M417 lineages with the remnants of the Novodanilovka population.


4. Repin / Early Yamna expansion

We don’t have direct data on early Repin settlers. But we do have a very close representative: Afanasevo, a population we know comes directly from the Repin/late Khvalynsk expansion ca. 3500/3300 BC (just before the emergence of Early Yamna), and which shows fully Steppe_Eneolithic-like ancestry.


Compared to this eastern Repin expansion that gave Afanasevo, the late Repin expansion to the west ca. 3300 BC that gave rise to the Yamna culture was one of colonization, evidenced by the admixture with north Pontic (Sredni Stog-like) populations, no doubt through exogamy:


This admixture is also found (in lesser proportion) in east Yamna groups, which supports the high mobility and exogamy practices among western and eastern Yamna clans, not only with locals:


5. Corded Ware

Corded Ware represents a quite homogeneous expansion of a late Sredni Stog population, compatible with the traditional location of Proto-Corded Ware peoples in the steppe-forest/forest zone of the Dnieper-Dniester region.


We don’t have a comparison with Ukraine_Eneolithic or Corded Ware samples in Wang et al. (2018), but we do have proximate sources for Abashevo, when compared to the Poltavka population (with which it admixed in the Volga-Ural steppes): Sintashta, Potapovka, Srubna (with further Abashevo contribution), and Andronovo:


The two CWC outliers from the Baltic show what I thought was an admixture with Yamna. However, given the previous mixture of Eneolithic_Steppe in north Pontic steppe-forest populations, this elevated “steppe ancestry” found in Baltic_LN (similar to west Yamna) seems rather an admixture of Baltic sub-Neolithic peoples with a north Pontic Eneolithic_Steppe-like population. Late Repin settlers also admixed with a similar population during its colonization of the north Pontic area, hence the Baltic_LN – west Yamna similarities.

NOTE. A direct admixture with west Yamna populations through exogamy by the ancestors of this Baltic population cannot be ruled out yet (without direct access to more samples), though, because of the contacts of Corded Ware with west Yamna settlers in the forest-steppe regions.


A similar case is found in the Yamna outlier from Mednikarovo south of the Danube. It would be absurd to think that Yamna from the Balkans comes from Corded Ware (or vice versa), just because the former is closer in the PCA to the latter than other Yamna samples. The same error is also found e.g. in the Corded Ware → Bell Beaker theory, because of their proximity in the PCA and their shared “steppe ancestry”. All those theories have been proven already wrong.

NOTE. A similar fallacy is found in potential Sintashta→Mycenaean connections, where we should distinguish statistically that result from an East/West Yamna + Balkans_BA admixture. In fact, genetic links of Mycenaeans with west Yamna settlers prove this (there are some related analyses in Anthrogenica, but the site is down at this moment). To try to relate these two populations (separated more than 1,000 years before Sintashta) is like comparing ancient populations to modern ones, without the intermediate samples to trace the real anthropological trail of what is found…Pure numbers and wishful thinking.


Yamna and Corded Ware show a similar “steppe ancestry” due to convergence. I have said so many times (see e.g. here). This was clear long ago, just by looking at the Y-chromosome bottlenecks that differentiate them – and Tomenable noticed this difference in ADMIXTURE from the supplementary materials in Mathieson et al. (2017), well before Wang et al. (2018).

This different stock stems from (1) completely different ancestral populations + (2) different, long-lasting Y-chromosome bottlenecks. Their similarities come from the two neighbouring cultures admixing with similar populations.

If all this does not mean anything, and each lab was going to support some pre-selected archaeological theories from the 1960s or the 1980s, coupled with outdated linguistic models no matter what – Anthony’s model + Ringe’s glottochronological tree of the early 2000s in the Reich Lab; and worse, Kristiansen’s CWC-IE + Germano-Slavonic models of the 1940s in the Copenhagen group – , I have to repeat my question again:

What’s (so much published) ancient DNA useful for, exactly?

See also


Iron Age bottleneck of the Proto-Fennic population in Estonia


Demographic data and figures derived from Estonian Iron Age graves, by Raili Allmäe, Papers on Anthropology (2018) 27(2).

Interesting excerpts (emphasis mine):


Inhumation and cremation burials were both common in Iron Age Estonia; however, the pattern which burials were prevalent has regional as well temporal peculiarities. In Estonia, cremation burials appear in the Late Bronze Age (1100–500 BC), for example, in stone-cist graves and ship graves, although inhumation is still characteristic of the period [28, 18]. Cremation burials have occasionally been found beneath the Late Bronze Age cists and the Early Iron Age (500 BC–450 AD) tarand graves [30, 28]. In south-eastern Estonia, including Setumaa, the tradition to bury cremated human remains in pit graves also appeared in the Bronze Age and lasted during the Pre-Roman period (500 BC–50 AD) and the Roman Iron Age (50–450 AD), and even up to the medieval times [30, 23, 33, 9]. During the Early Iron Age, cremations appear in cairn graves and have occasionally been found in many Pre-Roman early tarand graves where they appear with inhumations [27, 28, 19, 20, 21, 22, 24]. In Roman Iron Age tarand graves, cremation as well inhumation were practiced [28, 36, 37]; however, cremation was the prevailing burial practice during the Roman Iron Age, for example, in tarand graves in south-eastern Estonia [30, 28]. Roman Iron Age (50 AD–450 AD) burial sites have not been found in continental west Estonia [28, 38]). At the beginning of the Middle Iron Age (450–800 AD), burial sites, for example stone graves without a formal structure, like Maidla I, Lihula and Ehmja ‘Varetemägi’, appear in Läänemaa, west Estonia; in these graves cremations as well inhumations have been found [39, 48]. Like underground cremation burial, the stone grave without a formal structure was the most common grave type during the Late Iron Age (800– 1200 AD) in west Estonia [39, 35, 48]. In south-eastern and eastern Estonia, sand barrows with cremation burials appeared at the beginning of the Middle Iron Age. Cremation barrows are attributed to the Culture of Long Barrows and are most numerous in the villages Laossina and Rõsna in northern Setomaa, on the western shore of Lake Peipsi [8, 48]. Apparently during the Iron Age, the practiced burial customs varied in Estonia.

Typical prehistoric Estonian graves. Top: Cist-graves common during the Bronze Age, by Terker (GNU FDL 1.2). Bottom: Tarand graves of the Iron Age, by Marika Mägi (2017)


Three Iron Age cremation graves from south-eastern Estonia and four graves including cremations as well inhumations from western Estonia were analysed by osteological and palaeodemographic methods in order to estimate the age and sex composition of burial sites, and to propose some possible demographic figures and models for living communities.

The crude birth/death rate estimated on the basis of juvenility indices varied between 55.1‰ and 60.0‰ (58.5‰ on average) at Rõsna village in south-eastern Estonia in the Middle Iron Age. The birth/death rates based on juvenility indices for south eastern graves varied to a greater degree. The estimated crude birth/death rate was somewhat lower (38.9‰) at Maidla in the Late Iron Age and extremely high (92.1‰) at Maidla in the Middle Iron Age, which indicates an unsustainable community. High crude birth/death rates are also characteristic of Poanse tarand graves from the Pre-Roman Iron Age – 92.3‰ for the 1st grave and 69.6‰for the 2nd grave. Expectedly, newborn life expectancies are extremely low in both communities – 10.8 years at Poanse I and 14.4 years at Poanse II. Most likely, both Maidla I and Poanse I were unsustainable communities.

Locations of the investigated Estonian Iron Age graves. Map by R. Allmäe

According to the main model where the given period of grave usage is 150 years, the burial grounds were most likely exploited by communities of 3–14 people. In most cases, this corresponds to one family or household. In comparison with other graves, Maidla II stone grave in western Estonia and Rõsna-Saare I barrow cemetery in south-eastern Estonia could have been used by a somewhat larger community, which may mean an extended family, a larger household or usage by two nuclear families.

More papers on the same subject by the author – who participated in the recent Mittnik et al. (2018) paper – include Observations On Estonian Iron Age Cremations (2013), and The demography of Iron Age graves in Estonia (2014).

Fast life history in Iron Age Estonia

While the demographic situation in the Gulf of Finland during the Iron Age is not well known – and demography is always difficult to estimate based on burials, especially when cremation is prevalent – , there is a clear genetic bottleneck in Finns, which has been estimated precisely to this period, coincident with the expansion of Proto-Fennic.

PCA of Estonian samples from the Bronze Age, Iron Age and Medieval times. Tambets et al. (2018, upcoming).

The infiltration of N1c lineages in Estonia – the homeland of Proto-Fennic – happened during the Iron Age – as of yet found in 3 out of 5 sampled Tarand graves – , while the previous period was dominated by 100% R1a and Corded Ware + Baltic HG ancestry. With the Iron Age, a slight shift towards Corded Ware ancestry can be seen, which probably signals the arrival of warrior-traders from the Alanino culture, close to the steppe. They became integrated through alliances and intermarriages in a culture of chiefdoms dominated by hill forts. Their origin in the Mid-Volga area is witnessed by their material culture, such as Tarand-like graves (read here a full account of events).

This new political structure, reminiscent of the chiefdom system in Sintashta (with a similar fast life history causing a bottleneck of R1a-Z645 lineages), coupled with the expansion of Fennic (and displaced Saamic) peoples to the north, probably caused the spread of N1c-L392 among Finnic peoples. The linguistic influence of these early Iron Age trading movements from the Middle Volga region can be seen in similarities between Fennic and Mordvinic, which proves that the Fenno-Saamic community was by then not only separated linguistically, but also physically (unlike the period of long-term Palaeo-Germanic influence, where loanwords could diffuse from one to the other).

NOTE. Either this, or the alternative version: an increase in Corded Ware ancestry in Estonia during the Iron Age marks the arrival of the first Fennic speakers ca. 800 BC or later, splitting from Mordvinic? A ‘Mordvin-Fennic’ group in the Volga, of mainly Corded Ware ancestry…?? Which comes in turn from a ‘Volga-Saamic’ population of Siberian ancestry in the Artic region??? And, of course, Palaeo-Germanic widely distributed in North-Eastern Europe with R1a during the Bronze Age! Whichever model you find more logical.


Corded Ware—Uralic (I): Differences and similarities with Yamna


This is the first of four posts on the Corded Ware—Uralic identification:

I was reading The Bronze Age Landscape in the Russian Steppes: The Samara Valley Project (2016), and I was really surprised to find the following excerpt by David W. Anthony:

The Samara Valley links the central steppes with the western steppes and is a north-south ecotone between the pastoral steppes to the south and the forest-steppe zone to the north [see figure below]. The economic contrast between pastoral steppe subsistence, with its associated social organizations, and forest-zone hunting and fishing economies probably explains the shifting but persistent linguistic border between forest-zone Uralic languages to the north (today largely displaced by Russian) and a sequence of steppe languages to the south, recently Turkic, before that Iranian, and before that probably an eastern dialect of Proto-Indo-European (Anthony 2007). The Samara Valley represents several kinds of borders, linguistic, cultural, and ecological, and it is centrally located in the Eurasian steppes, making it a critical place to examine the development of Eurasian steppe pastoralism.

Language map of the middle Volga-Ural region. After “Geographical Distribution of the Uralic Languages” by Finno-Ugrian Society, Helsinki, 1993.

Khokhlov (translated by Anthony) further insists on the racial and ethnic divide between both populations, Abashevo to the north, and Poltavka to the south, during the formation of the Abashevo – Sintashta-Potapovka community that gave rise to Proto-Indo-Iranians:

Among all cranial series in the Volga-Ural region, the Potapovka population represents the clearest example of race mixing and probably ethnic mixing as well. The cultural advancements seen in this period might perhaps have been the result of the mixing of heterogeneous groups. Such a craniometric observation is to some extent consistent with the view of some archaeologists that the Sintashta monuments represent a combination of various cultures (principally Abashevo and Poltavka, but with other influences) and therefore do not correspond to the basic concept of an archaeological culture (Kuzmina 2003:76). Under this option, the Potapovka-Sintashta burial rite may be considered, first, a combination of traits to guarantee the afterlife of a selected part of a heterogeneous population. Second, it reflected a kind of social “caste” rather than a single population. In our view, the decisive element in shaping the ethnic structure of the Potapovka-Sintashta monuments was their extensive mobility over a fairly large geographic area. They obtained knowledge of various cultures from the populations with whom they interacted.

Late Middle Bronze Age cultures with the Proto-Indo-Iranian Sintashta-Potapovka-Filatovka group (shaded). After Anthony (2007 Figure 15.5), from Anthony (2016).

Interesting is also this excerpt about the predominant population in the Abashevo – Sintashta-Potapovka admixture (which supports what Chetan said recently, although this does not seemed backed by Y-DNA haplogroups found in the richest burials), coupled with the sign of incoming “Uraloid” peoples from the east, found in both Sintashta and eastern Abashevo:

The socially dominant anthropological component was Europeoid, possibly the descendants of Yamnaya. The association of craniofacial types with archaeological cultures in this period is difficult, primarily because of the small amount of published anthropological material of the cultures of steppe and forest belt (Balanbash, Vol’sko-Lbishche) and the eastern and southern steppes (Botai-Tersek). The crania associated with late MBA western Abashevo groups in the Don-Volga forest zone were different from eastern Abashevo in the Urals, where the expression of the Old Uraloid craniological complex was increased. Old Uraloid is found also on a single skull of Vol’sko-Lbishche culture (Tamar Utkul VII, Kurgan 4). Potentially related variants, including Mongoloid features, could be found among the Seima-Turbino tribes of the forest-steppe zone, who mixed with Sintashta and Abashevo. In the Sintashta Bulanova cemetery from the western Urals, some individuals were buried with implements of Seima-Turbino type (Khalyapin 2001; Khokhlov 2009; Khokhlov and Kitov 2009). Previously, similarities were noted between some individual skulls from Potapovka I and burials of the much older Botai culture in northern Kazakhstan (Khokhlov 2000a). Botai-Tersek is, in fact, a growing contender for the source of some “eastern” cranial features.

Facial reconstructions based on skulls from (a) Khvalynsk II Grave 24, a young adult male; (b) Poludin Grave 6, Yamnaya culture, a mature male (both by A. I. Nechvaloda); and (c) Luzanovsky cemetery, Srubnaya culture (by L. T. Yablonsky). In Khokhlov (2016).

The wave of peoples associated with “eastern” features can be seen in genetics in the Sintashta outliers from Narasimhan et al. (2018), and it probably will be eventually seen in Abashevo, too. These may be related to the Seima-Turbino international network – but most likely it is directly connected to Sintashta through the starting Andronovo and Seima-Turbino horizons, by admixing of prospective groups and small-scale back-migrations.

Corded Ware – Yamna similarities?

So, if peoples of north-eastern Europe have been assumed for a long time to be Uralic speakers, what is happening with the Corded Ware = IE obsession? Is it Gimbutas’ ghost possessing old archaeologists? Probably not.

It is about certain cultural similarities evident at first sight, which have been traditionally interpreted as a sign of cultural diffusion or migration. Not dissimilar to the many Bell Beaker models available, where each archaeologist is pushing certain differences, mixing what seemed reasonable, what still might seem reasonable, and what certainly isn’t anymore after the latest ancient DNA data.

“European dialect” expansion of Proto-Indo-European according to Gimbutas (1963)

The initial models of Gimbutas, Kristiansen, or Anthony – which are known to many today – were enunciated in the infancy of archaeological studies in the regions, during and just after the fall of the USSR, and before many radiocarbon dates that we have today were published (with radiocarbon dating being still today in need of refinement), so it is only logical that gross mistakes were made.

We have similar gross mistakes related to the origins of Bell Beakers, and studying them was certainly easier than studying eastern data.

  • Gimbutas believed – based mainly on Kurgan-like burials – that Bell Beaker formed from a combination of Yamna settlers with the Vučedol culture, so she was not that far from the truth.
  • The expansion of Corded Ware from peoples of the North Pontic forest-steppe area, proposed by Gimbutas and later supported also by Kristiansen (1989) as the main Indo-European expansion – , is probably also right about the approximate origins of the culture. Only its ‘Indo-European’ nature is in question, given the differences with Khvalynsk and Yamna evolution.
  • Anthony only claimed that Yamna migrants settled in the Balkans and along the Danube into the Hungarian steppes. He never said that Corded Ware was a Yamna offshoot until after the first genetic papers of 2015 (read about his newest proposal). He initially claimed that only certain neighbouring Corded Ware groups “adopted” Indo-European (through cultural diffusion) because of ‘patron-client’ relationships, and was never preoccupied with the fate of Corded Ware and related cultures in the east European forest zone and Finland.

So none of them was really that far from the true picture; we might say a lot people are more way off the real picture today than the picture these three researchers helped create in the 1990s and 2000s. Genetics is just putting the last nail in the coffin of Corded Ware as a Yamna offshoot, instead of – as we believed in the 2000s – to Vučedol and Bell Beaker.

So let’s revise some of these traditional links between Corded Ware and Yamna with today’s data:


Even more than genetics – at least until we have an adequate regional and temporary sampling – , archaeological findings lead what we have to know about both cultures.

It is essential to remember that Corded Ware, starting ca. 3000/2900 BC in east-central Europe, has been proposed to be derived from Early Yamna, which appeared suddenly in the Pontic-Caspian steppes ca. 3300 BC (probably from the late Repin expansion), and expanded to the west ca. 3000.

Early Yamna is in turn identified as the expanding Late Proto-Indo-European community, which has been confirmed with the recent data on Afanasevo, Bell Beaker, and Sintashta-Potapovka and derived cultures.

The question at hand, therefore, is if Corded Ware can be considered an offshoot of the Late PIE community, and thus whether the CWC ethnolinguistic community – proven in genetics to be quite homogeneous – spoke a Late PIE dialect, or if – alternatively – it is derived from other neighbouring cultures of the North Pontic region.

NOTE. The interpretation of an Indo-Slavonic group represented by a previous branching off of the group is untenable with today’s data, since Indo-Slavonic – for those who support it – would itself be a branch of Graeco-Aryan, and Palaeo-Balkan languages expanded most likely with West Yamna (i.e. R1b-L23, mainly R1b-Z2103) to the south.

The convoluted alternative explanation would be that Corded Ware represents an earlier, Middle PIE branch (somehow carrying R1a??) which influences expanding Late PIE dialects; this has been recently supported by Kortlandt, although this simplistic picture also fails to explain the Uralic problem.

Kurgans: The Yamna tradition was inherited from late Repin, in turn inherited from Khvalynsk-Novodanilovka proto-Kurgans. As for the CWC tradition, it is unclear if the tumuli were built as a tradition inherited from North and West Pontic cultures (in turn inherited or copied from Khvalynsk-Novodanilovka), such as late Trypillia, late Kvityana, late Dereivka, late Sredni Stog; or if they were built because of the spread of the ‘Transformation of Europe’, set in motion by the Early Yamna expansion ca. 3300-3000 BC (as found in east-central European cultures like Coţofeni, Lizevile, Șoimuș, or the Adriatic Vučedol). My guess is that it inherits an older tradition than Yamna, with an origin in east-central Europe, because of the mound-building distribution in the North Pontic area before the Yamna expansion, but we may never really know.

Distribution of Pit-Grave burials west of the Black Sea likely dating to the 2nd half of the IVth millennium BC (triangles: side-crouched burials; filled circles: supine extended burials; open circles: suspected). Frînculeasa, Preda, and Heyd (2015)

Burial rite: Yamna features (with regional differences) single burials with body on its back, flexed upright knees, poor grave goods, common orientation east-west (heads to the west) inherited from Repin, in turn inherited from Khvalynsk-Novodanilovka. CWC tradition – partially connected to Złota and surrounding east-central European territories (in turn from the Khvalynsk-Novodanilovka expansion) – features single graves, body in fetal position, strict gender differentiation – men on the right, women on the left -, looking to the south, graves with standardized assemblages (objects representing affirmation of battle, hunting, and feasting). The burial rites clearly represent different ideologies.

Left: Pit-Grave burial types expanded with Khvalynsk-Novodanilovka. Right: Pit-Grave burial types associated with the Yamna expansion and influence. Frînculeasa, Preda, and Heyd (2015)

Corded decoration: Corded ware decoration appears in the Balkans during the 5th millennium, and represents a simple technique whereby a cord is twisted, or wrapped around a stick, and then pressed directly onto the fresh surface of a vessel leaving a characteristic decoration. It appears in many groups of the 5th and 4th millennium BC, but it was Globular Amphorae the culture which popularized the drinking vessels and their corded ornamentation. It appears thus in some regional groups of Yamna, but it becomes the standard pottery only in Corded Ware (especially with the A-horizon), which shows continuity with GAC pottery.

Origins of the first Corded Ware horizon (5th millennium BC) after the Khvalynsk-Novodanilovka expansion. Corded Ware (circles) and horse-head scepters (rectangles) and other steppe elements (triangles). Image from Bulatović (2014).

Economy: Yamna expands from Repin (and Repin from Khvalynsk-Novodanilovka) as a nomadic or semi-nomadic purely pastoralist society (with occasional gathering of wild seeds), which naturally thrives in the grasslands of the Pontic-Caspian, lower Danube and Hungarian steppes. Corded Ware shows agropastoralism (as late Eneolithic forest-steppe and steppe groups of eastern Europe, such as late Trypillian, TRB, and GAC groups), inhabits territories north of the loess line, with heavy reliance of hunter-gathering depending on the specific region.

Cattle herding: Interestingly, both west Yamna and Corded Ware show more reliance on cattle herding than other pastoralist groups, which – contrasted with the previous Eneolithic herding traditions of the Pontic-Caspian steppe, where sheep-goats predominate – make them look alike. However, the cattle-herding economy of Yamna is essential for its development from late Repin and its expansion through the steppes (over western territories practising more hunter-gathering and sheep-goat herding economy), and it does not reach equally the Volga-Ural region, whose groups keep some of the old subsistence economy (read more about the late Repin expansion). Corded Ware, on the other hand, inherits its economic strategy from east European groups like TRB, GAC, and especially late Trypillian communities, showing a predominance of cattle herding within an agropastoral community in the forest-steppe and forest zones of Volhynia, Podolia, and surrounding forest-steppe and forest regions.

Scheme of interlinked socio-economic-ideological innovations forming the Yamnaya. Frînculeasa, Preda, and Heyd (2015)

Horse riding: Horse riding and horse transport is proven in Yamna (and succeeding Bell Beaker and Sintashta), assumed for late Repin (essential for cattle herding in the seas of grasslands that are the steppes, without nearby water sources), quite likely during the Khvalynsk expansion (read more here), and potentially also for Samara, where the predominant horse symbolism of early Khvalynsk starts. Corded Ware – like the north Pontic forest-steppe and forest areas during the Eneolithic – , on the other hand, does not show a strong reliance on horse riding. The high mobility and short-term settlements characteristic of Corded Ware, that are often associated with horse riding by association with Yamna, may or may not be correct, but there is no need for horses to explain their herding economy or their mobility, and the north-eastern European areas – the one which survived after Bell Beaker expansion – did certainly not rely on horses as an essential part of their economy.

NOTE: I cannot think of more supposed similarities right now. If you have more ideas, please share in the comments and I will add them here.

Genetic similarities

EHG: This is the clearest link between both communities. We thought it was related to the expansion of ANE-related ancestry to the west into WHG territory, but now it seems that it will be rather WHG expanding into ANE territory from the Pontic-Caspian region to the east (read more on recent Caucasus Neolithic, on , and on Caucasus HG).

NOTE. Given how much each paper changes what we know about the Palaeolithic, the origin and expansion of the (always developing) known ancestral components and specific subclades (see below) is not clear at all.

CHG: This is the key link between both cultures, which will delimit their interaction in terms of time and space. CHG is intermediate between EHG and Iran N (ca. 8000 BC). The ancestry is thus linked to the Caucasus south of the steppe before the emergence of North Pontic (western) and Don-Volga-Ural (eastern) communities during the Mesolithic. The real question is: when we have more samples from the steppe and the Caucasus during the Neolithic, how many CHG groups are we going to find? Will the new specific ancestral components (say CHG1, CHG2, CHG3, etc.) found in Yamna (from Khvalynsk, in the east) and Corded Ware (probably from the North Pontic forest-steppe) be the same? My guess is, most likely not, unless they are mediated by the Khvalynsk-Novodanilovka expansion (read more on CHG in the Caucasus).

Formation of Yamna and CHG contribution, in Damgaard et al. (Science 2018). A 10-leaf model based on combining the models in Fig. S16 and Fig. S19 and re-estimating the model parameters.

WHG/EEF: This is the obvious major difference – known today – in the formation of both communities in the steppe, and shows the different contacts that both groups had at least since the Eneolithic, i.e. since the expansion of Repin with its renewed Y-DNA bottleneck, and probably since before the early Khvalynsk expansion (read more on Yamna-Corded Ware differences contrasting with Yamna-Afanasevo, Yamna-Bell Beaker, and Yamna-Sintashta similarities).

NOTE 1. Some similarities between groups can be seen depending on the sampled region; e.g. Baltic groups show more similarities with southern Pontic-Caspian steppe populations, probably due to exogamy.

Tested qpGraph model in Tambets et al. (2018). The qpGraph model fitting the data for the tested populations. “Colour codes for the terminal nodes: pink—modern populations (‘Population X’ refers to test population) and yellow—ancient populations (aDNA samples and their pools). Nodes coloured other than pink or yellow are hypothetical intermediate populations. We putatively named nodes which we used as admixture sources using the main recipient among known populations. The colours of intermediate nodes on the qpGraph model match those on the admixture proportions panel.”

NOTE 2. We have this information on the differences in “steppe ancestry” between Yamna and Corded Ware, compared to previous studies, because now we have more samples of neighbouring, roughly contemporaneous Eneolithic groups, to analyse the real admixture processes. This kind of fine scale studies is what is going to show more and more differences between Khvalynsk-Yamna and Sredni Stog-Corded Ware as more data pours in. The evolution of both communities in archaeology and in PCA (see below) is probably witness to those differences yet to be published.

R1: Even though some people try very hard to think in terms of “R1” vs. (Caucasus) J or G or any other upper clade, this is plainly wrong. It is possible, given what we know now, that Q1a2-M242 expanded ANE ancestry to the west ca. 13000 BC, while R1b-P279 expanded WHG ancestry to the east with the expansion of post-Swiderian cultures, creating EHG as a WHG:ANE cline. The role of R1a-M459 is unknown, but it might be related to any of these migrations, or others (plural) along northern Eurasia (read more on the expansion of R1b-P279, on Palaeolithic Q1a2, and on R1a-M417).

NOTE. I am inclined to believe in a speculative Mesolithic-Early Neolithic community involving Eurasiatic movements accross North Eurasia, and Indo-Uralic movements in its western part, with the last intense early Uralic-PIE contacts represented by the forming west (Mariupol culture) and east (Don-Volga-Ural cultures, including Samara) communities developing side by side. Before their known Eneolithic expansions, no large-scale Y-DNA bottleneck is going to be seen in the Pontic-Caspian steppe, with different (especially R1a and R1b subclades) mixed among them, as shown in North Pontic Neolithic, Samara HG, and Khvalynsk samples.

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. Marked and labelled are the Balkan samples referenced in this text An EHG and a Caucasus ‘clouds’ have been drawn, leaving Pontic-Caspian steppe and derived groups between them. See the original file here.

Corded Ware and ‘steppe ancestry’

If we take a look at the evolution of Corded Ware cultures, the expansion of Bell Beakers – dominated over most previous European cultures from west to east Europe – influenced the development of the whole European Bronze Age, up to Mierzanowice and Trzciniec in the east.

The only relevant unscathed CWC-derived groups, after the expansion of Sintashta-Potapovka as the Srubna-Andronovo horizon in the Eurasian steppes, were those of the north-eastern European forest zone: between Belarus to the west, Finland to the north, the Urals to the east, and the forest-steppe region to the south. That is, precisely the region supposed to represent Uralic speakers during the Bronze Age.

This inconsistency of steppe ancestry and its relation with Uralic (and Balto-Slavic) peoples was observed shortly after the publication of the first famous 2015 papers by Paul Heggarty, of the Max-Planck Institute for Evolutionary Anthropology (read more):

Haak et al. (2015) make much of the high Yamnaya ancestry scores for (only some!) Indo-European languages. What they do not mention is that those same results also include speakers of other languages among those with the highest of all scores for Yamnaya ancestry. Only these are languages of the Uralic family, not Indo-European at all; and their Yamnaya-ancestry signals are far higher than in many branches of Indo-European in (southern) Europe. Estonian ranks very high, while speakers of the very closely related Finnish are curiously not shown, and nor are the Saami. Hungarian is relevant less directly since this language arrived only c. 900 AD, but also high.


These data imply that Uralic-speakers too would have been part of the Yamnaya > Corded Ware movement, which was thus not exclusively Indo-European in any case. And as well as the genetics, the geography, chronology and language contact evidence also all fit with a Yamnaya > Corded Ware movement including Uralic as well as Balto-Slavic.

Both papers fail to address properly the question of the Uralic languages. And this despite — or because? — the only Uralic speakers they report rank so high among modern populations with Yamnaya ancestry. Their linguistic ancestors also have a good claim to have been involved in the Corded Ware and Yamnaya cultures, and of course the other members of the Uralic family are scattered across European Russia up to the Urals.

NOTE. Although the author was trying to support the Anatolian hypothesis – proper of glottochronological studies often published from the Max Planck Institute – , the question remains equally valid: “if Proto-Indo-European expands with Corded Ware and steppe ancestry, what is happening with Uralic peoples?”

For my part, I claimed in my draft that ancestral components were not the only relevant data to take into account, and that Y-DNA haplogroups R1a and R1b (appearing separately in CWC and Yamna-Bell Beaker-Afanasevo), together with their calculated timeframes of formation – and therefore likely expansion – did not fit with the archaeological and linguistic description of the spread of Proto-Indo-European and its dialects.

In fact, it seemed that only one haplogroup (R1b-M269) was constantly and consistenly associated with the proposed routes of Late PIE dialectal expansions – like Anthony’s second (Afanasevo) and third (Lower Danube, Balkan) waves. What genetics shows fits seamlessly with Mallory’s association of the North-West Indo-European expansion with Bell Beakers (read here how archaeologists were right).

Map of the much beloved steppe (or “Yamnaya”) ancestry in modern populations, by Balanovsky. Modified from Klejn (2017).

More precise inconsistencies were observed after the publication of Olalde et al. (2017) and Mathieson et al. (2017), by Volker Heyd in Kossinna’s smile (2017). Letting aside the many details enumerated (you can read a summary in my latest draft), this interesting excerpt is from the conclusion:

NOTE. An open access ealier draft version of the paper is offered for download by the author.

Simple solutions to complex problems are never the best choice, even when favoured by politicians and the media. Kossinna also offered a simple solution to a complex prehistoric problem, and failed therein. Prehistoric archaeology has been aware of this for a century, and has responded by becoming more differentiated and nuanced, working anthropologically, scientifically and across disciplines (cf. Müller 2013; Kristiansen 2014), and rejecting monocausal explanations. The two aDNA papers in Nature, powerful and promising as they are for our future understanding, also offer rather straightforward messages, heavily pulled by culture-history and the equation of people with culture. This admittedly is due partly to the restrictions of the medium that conveys them (and despite the often relevant additional detail given as supplementary information, which is unfortunately not always given full consideration).

While I have no doubt that both papers are essentially right, they do not reflect the complexity of the past. It is here that archaeology and archaeologists contributing to aDNA studies find their role; rather than simply handing over samples and advising on chronology, and instead of letting the geneticists determine the agenda and set the messages, we should teach them about complexity in past human actions and interactions. If accepted, this could be the beginning of a marriage made in heaven, with the blessing smile of Gustaf Kossinna, and no doubt Vere Gordon Childe, were they still alive, in a reconciliation of twentieth- and twenty-first-century approaches. For us as archaeologists, it could also be the starting point for the next level of a new archaeology.

Main distribution of Yamnaya kurgans in the Pontic-Caspian steppe of modern day Russia, Ukraine, and Kazakhstan, and its western branch in modern south-east European countries of Romania, Bulgaria, Serbia, and Hungary, with numbers of excavated kurgans and graves given. Picture: Volker Heyd (2018).

The question was made painfully clear with the publication of Olalde et al. (2018) & Mathieson et al. (2018), where the real route of Yamna expansion into Europe was now clearly set through the steppes into the Carpathian basin, later expanded as Bell Beakers.

This has been further confirmed in more recent papers, such as Narasimhan et al. (2018), Damgaard et al. (2018), or Wang et al. (2018), among others.

However, the discussion is still dominated by political agendas based on prevalent Y-DNA haplogroups in modern countries and ethnic groups.


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…


Polygyny as a potential reason for Y-DNA bottlenecks among agropastoralists


Open access Greater wealth inequality, less polygyny: rethinking the polygyny threshold model by Ross et al. Journal of the Royal Society Interface (2018).

Interesting excerpts, from the discussion (emphasis mine):

We use cross-cultural data and a new mutual mate choice model to propose a resolution to the polygyny paradox. Following Oh et al. [17], we extend the standard polygyny threshold model to a mutual mate choice model that accounts for both female supply to, and male demand for, polygynous matchings, in the light of the importance of, and inequality in, rival and non-rival forms of wealth. The empirical results presented in figures 5 and 6 demonstrate two phenomena that are jointly sufficient to generate a transition to more frequent monogamy among populations with a co-occurring transition to a more unequal, highly stratified, class-based social structure. In such populations, fewer men can cross the wealth threshold required to obtain a second wife, and those who do may be fabulously wealthy, but—because of diminishing marginal fitness returns to increasing number of marriages—do not acquire wives in full proportion to their capacity to support them with rival wealth. Together, these effects reduce the population-level fraction of wives in polygynous marriages.

Our model demonstrates that a low population-level frequency of polygyny will be an equilibrium outcome among fitness maximizing males and females in a society characterized by a large class of wealth-poor peasants and a small class of exceptionally wealthy elite. Our mutual mate choice model thus provides an empirically plausible resolution to the polygyny paradox and the transition to monogamy which co-occurred with the rise of highly unequal agricultural populations.

(a) Mean frequency of married women who are married polygynously by production system (+2 s.e.) using the Standard Cross-Cultural Sample [30]. Rates of polygyny are measured with variable ]872, per cent of wives with co-wives. (b) Rates of monogamy and polygyny by production system are measured with variable ]861, the standard polygamy code. Data on subsistence come from variable ]858, categorized subsistence. In general, agricultural populations show reduced rates of polygyny and increased rates of monogamy relative to other subsistence systems. See electronic supplementary material for more information. (c) Gini of wealth by production system in our sample.

The reasons for this decrease in marginal fitness returns are explained as either a) a potential missing of important rival forms of wealth in the statistical model, or b) one or more of the following reasons:

  • [A] male’s time and attention are rival inputs to his own fitness (…) A single rich man will have to defend his 10 wives from nine unmarried men on average.”As the wealth ratio grows even more skewed, this situation could become increasingly difficult to manage (e.g. requiring the use of eunochs to defend harems [74]).
  • A related possibility is that a growing number of unmarried men could socially censure wealthy polygynous males, imposing costs on them that reduce male demand for and/or female supply to polygynous marriage [23,24]. (…)
  • A third possibility is that sexually transmitted infection (STI) burden [22,75] could diminish returns to polygyny, if polygyny enhances infection rates [76,77]. (…)
  • Finally, impediments to cooperation or even outright conflict among co-wives can be greater as the number of wives increases. Interference competition among co-wives could impose significant fitness costs in settings where effective child rearing benefits from cooperation [79,80].(…)
between the Gini coefficient on completed rival wealth and per cent completed female polygyny.

I have previously argued against some reasons traditionally given to explain the replacement of native male populations after migrations (i.e. polygyny, slavery, targeted male extermination, etc.), because I believe that a gradual successful expansion of patrilineal clans over some generations based on wealth alone is enough to explain the obvious Y-DNA bottlenecks that happened in many different prehistoric and historic cultures (especially among steppe pastoralists, including Indo-Europeans).

I realize that I haven’t really used any study to support my opinion, though, and data from modern and ancient pastoralists from different regions seem to contradict it, so maybe ancient DNA can show that Indo-Europeans had often children with more than one woman at the same time. I don’t remember seeing that kind of information in supplementary materials to date. From memory I can think of maybe two or three examples of agnate siblings published, but I doubt the archaeological age estimation (based on simple observation of skeletal remains) combined with radiocarbon age (usually given with broad CI) could be enough to prove a similar age of conception. Maybe a case of many siblings clearly of the same age and from many different mothers in the same burial could be a strong proof of this…

I recently read that theoretical models are actually trusted by no one except for the researchers who propose them, and experimental data are trusted by everyone except for the researchers who worked with them. I cannot agree more. However, we lack information about this question (as far as I know), so we may have to rely on indirect estimations, like the kind of models presented in the paper (or the one proposed for Post-Neolithic Y-chromosome bottlenecks).

The Late Proto-Indo-European word for bride comes from a root meaning ‘drive, lead’, hence literally ‘deportation’, so the bride was transferred from her father’s family to her husband’s house. Marriage was certainly an asymmetrical contract for its members, and the reconstructible word for ‘dowry’ further supports the weaker position of the wife in it. Also, ancient marriage could differ from a family agreement, because marriage by elopement, bride kidnapping or hostage was probably common (more or less socially regulated) for people belonging the same culture. Apart from this, I don’t know about reconstructed linguistic data pointing to polygyny, and I doubt archaeological data alone – without genetics – can help.


Male-biased expansions and migrations also observed in Northwestern Amazonia

Open access preprint Cultural Innovations influence patterns of genetic diversity in Northwestern Amazonia, by Arias et al., bioRxiv (2018).

Abstract (emphasis mine):

Human populations often exhibit contrasting patterns of genetic diversity in the mtDNA and the non-recombining portion of the Y-chromosome (NRY), which reflect sex-specific cultural behaviors and population histories. Here, we sequenced 2.3 Mb of the NRY from 284 individuals representing more than 30 Native-American groups from Northwestern Amazonia (NWA) and compared these data to previously generated mtDNA genomes from the same groups, to investigate the impact of cultural practices on genetic diversity and gain new insights about NWA population history. Relevant cultural practices in NWA include postmarital residential rules and linguistic-exogamy, a marital practice in which men are required to marry women speaking a different language. We identified 2,969 SNPs in the NRY sequences; only 925 SNPs were previously described. The NRY and mtDNA data showed that males and females experienced different demographic histories: the female effective population size has been larger than that of males through time, and both markers show an increase in lineage diversification beginning ~5,000 years ago, with a male-specific expansion occurring ~3,500 years ago. These dates are too recent to be associated with agriculture, therefore we propose that they reflect technological innovations and the expansion of regional trade networks documented in the archaeological evidence. Furthermore, our study provides evidence of the impact of postmarital residence rules and linguistic exogamy on genetic diversity patterns. Finally, we highlight the importance of analyzing high-resolution mtDNA and NRY sequences to reconstruct demographic history, since this can differ considerably between males and females.

MDS plots for mtDNA and NRY. Stress values (within parentheses) are indicated in percentages.

Looking more precisely at the different groups (even with the resampling approach), there are no significant differences between matrilocal and patrilocal groups. At best, as the study proposes, “this is just one of the factors at play in structuring the observed genetic variation”.

Interesting excerpts:

(…) we found evidence that the patterns of genetic differentiation depend on the geographical scale of the study. The magnitude of between-population differentiation in the NRY compared to the mtDNA is smaller when looking at the continental scale than in NWA (Figure 6). This is in agreement with the findings of Wilkins and Marlowe (2006), who showed that the excess of between-population differentiation for the NRY in comparison to the mtDNA decreases when comparing more geographically distant populations. Heyer et al. (2012) and Wilkins and Marlowe (2006) have proposed that at a local scale the patterns of genetic diversity reflect cultural practices over a relatively small number of generations, whereas at a larger geographic scale the genetic diversity reflects old migration and/or old common ancestry patterns(Heyer et al. 2012; Wilkins and Marlowe 2006).

BSPs for the mtDNA and NRY sequences from NWA. The dotted lines indicate the 95% HPD intervals. Ne was corrected for generation time according to (Fenner 2005), using 26 years for mtDNA and 31 years for NRY.

The BSP plots and the diversity statistics indicate that overall the Ne of males has been smaller than that of females. One tentative explanation for this difference is that it reflects larger differences in reproductive success among males than among females. Some support for this explanation comes from the shape of the phylogenies (Supplementary Figures 1 and 6), since differences in reproductive success and the cultural transmission of fertility lead to imbalance phylogenies (Blum et al. 2006; Heyer et al. 2015). We estimated a common index of tree imbalance (Colless index) and calculated whether the mtDNA and NRY trees were more unbalanced than 1000 simulated trees generated under a Yule process (Bortolussi et al. 2006) (i.e. a simple pure birth process that assumes that the birth rate of new lineages is the same along the tree). We found that the NRY tree is more unbalanced than predicted by the Yule model (p-value=0.001), whereas the mtDNA tree is not significantly different from trees generated by the Yule model (p-value=0.628). It has been suggested that highly mobile hunter-gatherer societies, such as those typical of most of human prehistory, were polygynous bands (Dupanloup et al. 2003); similarly, nomadic horticulturalist Amazonian societies exhibit strong differences in reproductive success due to the common practice of polygyny, especially among community chiefs, whose offspring also enjoy a high fertility (Neel 1970; 1980; Neel and Weiss 1975).

Furthermore, a more recent expansion can be observed in the BSP based on the NRY, but not in the mtDNA BSP (Figure 5), indicating an expansion specifically in the paternal line. The reasons behind this recent male-biased population expansion, which starts ~3.5 kya, are as yet unclear. However, similar male-biased expansions have been observed in other studies using high-resolution NRY sequences (Batini et al. 2017; Karmin et al. 2015).


Post-Neolithic Y-chromosome bottleneck explained by cultural hitchhiking and competition between patrilineal clans

Open access study Cultural hitchhiking and competition between patrilineal kin groups explain the post-Neolithic Y-chromosome bottleneck, by Zeng, Aw, and Feldman, Nature Communications (2018).

Abstract (emphasis mine):

In human populations, changes in genetic variation are driven not only by genetic processes, but can also arise from cultural or social changes. An abrupt population bottleneck specific to human males has been inferred across several Old World (Africa, Europe, Asia) populations 5000–7000 BP. Here, bringing together anthropological theory, recent population genomic studies and mathematical models, we propose a sociocultural hypothesis, involving the formation of patrilineal kin groups and intergroup competition among these groups. Our analysis shows that this sociocultural hypothesis can explain the inference of a population bottleneck. We also show that our hypothesis is consistent with current findings from the archaeogenetics of Old World Eurasia, and is important for conceptions of cultural and social evolution in prehistory.

Relevant excerpts:

Tree of Y-chromosome genotypes from samples found among cultures with hunter-gatherer subsistence, and agropastoralist subsistence. The blue background represents hunter-gatherer subsistence while the green background represents agropastoralist subsistence. Letters in red circles match individuals from sites with their archaeological context. Note that R1b-P321 is synonymous with R1b-S116. Adapted from Figs. 3, 4, 5 and 6 of Kivisild67, with addition of information from Olalde et al.64. The vertical axis represents time; the position of branch points represent the ages of branch-defining mutations, with nomenclature and age from yfull (https://www.yfull.com/tree/)

Our hypothesis explains the bottleneck as a consequence of intergroup competition between patrilineal kin groups, which caused cultural hitchhiking between Y-chromosomes and cultural groups and reduction in Y-chromosomal diversity. Competition between demes can dramatically reduce genetic diversity within a population1, especially if the population is structured such that variation is greater between demes than within demes. Culturally transmitted kinship ideals and norms can cause homophilous sorting and limit interdemic gene flow, creating homogeneous demes that differ strongly from one another. Patrilineal corporate kin groups, with coresiding male group members descending from a common male ancestor, would produce such an effect on Y-chromosomes only, as patrilineal corporate kin groups generally coexist with female exogamy40, which would homogenize the mitochondrial gene pools of different groups41,42.

With intergroup competition between patrilineal corporate kin groups, two mechanisms would operate to reduce Y-chromosomal diversity. First, patrilineal corporate kin groups produce high levels of Y-chromosomal homogeneity within each social group due to common descent, as well as high levels of between-group variation. Second, the presence of such groups results in violent intergroup competition preferentially taking place between members of male descent groups, instead of between unrelated individuals. Casualties from intergroup competition then tend to cluster among related males, and group extinction is effectively the extinction of lineages.

There is evidence that other analogous situations involving gene-culture hitchhiking in culturally-defined social groups may have affected genetic diversity. Central Asian pastoralists, who are organized into patriclans, have high levels of intergroup competition and demonstrate ethnolinguistic and population-genetic turnover down into the historical period59. They also have a markedly lower diversity in Y-chromosomal lineages than nearby agriculturalists42,60. In fact, Central Asians are the only population whose male effective population size has not recovered from the post-Neolithic bottleneck; it remains disproportionately reduced, compared to female estimates using mtDNA4. Central Asians are also the only population to have star-shaped expansions of Y-chromosomes within the historical period, which may be due to competitive processes that led to the disproportionate political success of certain patrilineal clans60.

The simulation offers an interesting graphic. I had been thinking for some time about developing an interactive image with waves of expansion showing how only few haplogroups expand and thus their variability is reduced in successive migration waves, because a lot of people seemed not to be willing to accept this:

Schematic of the steps in the simulation, according to the order described in the algorithm. a (i) Patrilineal (PT) starting conditions, where cultural groups strictly determine haplogroup type. a (ii) The non-patrilineal (NPT) condition where they are perfectly uncorrelated. b The killing step, with a more (PT) and less (NPT) patrilineal starting condition. The number of deaths in each group is inversely related to group size. The blue cultural group goes extinct in both cases. This causes the haplogroup represented by the diamonds to go extinct in PT, but no haplogroup extinction occurs in NPT. c The mutation step, where a small number of individuals in the largest haplogroup change their haplogroup. d The regeneration step, where (i) is a replica of (b) PT (iii), and (d) (ii) shows how the original number of individuals before the killing step is restored by proportionally increasing the number of individuals in all cells. e Group fission step. Where an empty row occurs, the largest cultural group splits, and half the individuals form a new cultural group in the empty row. The step in which we remove cultural groups that are too small—between (c, d) (see Methods)—is not shown

You only have to imagine this process happening in many successive waves of expansion (external as well as internal to each culture) since the first Neolithic expansions in the steppe in the late-6th millennium BC, even before the formation of the Khvalynsk-Sredni Stog cultural-historical community, to understand what happened in the next thousands of years with evolving patrilineal clans and their distinct cultures.

The whole paper is an interesting read. It’s great to see sociology and genetics finally catch up and interact to develop more complex anthropological hypotheses.

The fact that this paper appears in mid-2018 and geneticists are beginning to discuss this only now when their statistical methods fail to explain the obvious (see David Reich’s recent interview) seems anachronistic, though, because all this was quite clear already in 2015 – at least for those who were looking for mainstream Yamna – Bell Beaker connections, instead of inventing new migration pathways to justify the results of certain statistical analyses

Anyway, better late than never.

Also, they use YFull estimates, which vindicates my use of them in the Indo-European demic diffusion model (2017). On the other hand, their use of these estimates right now in 2018 for R1a-M417 and R1b-M269 – when we know of a R1a-Z93 case much older than YFull’s estimated 5,000 YBP for this subclade, and possibly for R1b-L23, too, is the biggest pitfall in their temporal assessment, although the bottlenecks seen in Chalcolithic expansions seem to have indeed began during the Mesolithic-Neolithic transition in the steppe.

So, say goodbye (if you haven’t already) to dat fantasy ‘steppe people’ of mixed R1a/R1b descent cooperating with the same mixed steppe language, all represented by the Yamnaya™ ancestral component, and say hello to distinct, competing ethnolinguistic steppe groups during the Neolithic.