The Tungusic Ulchi population probably linked to haplogroup C2b1a

ulchi-marital

New paper (behind paywall) Demographic and Genetic Portraits of the Ulchi Population, by Balanovska et al. Russian Journal of Genetics (2018) 54(10):1245–1253.

Interesting excerpts (emphasis mine):

Marital structure. The intensity of interethnic marriages puts the existence of the Ulchi population at risk. The colorful ethnic composition of the Ulchi settlements is reflected in the marriage structure [see featured image]. We found that the proportion of single-ethnic marriages of the Ulchi is on average 51%. The greatest number of such marriages takes place in the village of Bulava. Marriages of Ulchi with Russians are in second place. Marriages with indigenous peoples of the Far East, Nanais, Nivkhs, Evenks, and others, are in third place. Thus, almost half of the Ulchi marriages are with representatives of other nationalities. Such a significant level of interethnic mixing makes it possible to talk about intense processes of assimilation of this indigenous people and puts to the forefront the problem of loss of the unique gene pool of the Ulchi.

Haplogroup C (its branch M48) was genotyped for its five subbranches with markers M86, B470, F13686, B93, and the marker at position 16645386 (GRCh37), which was found by our team for the first time. Variant B93 is rare in the Ulchi, and 14 samples (that is, more than a quarter of the entire gene pool of the Ulchi, Fig. 2) belong to M86 and its subvariants. Therefore, we genotyped STR markers of C-M86 carriers for the Ulchi and neighboring Amur populations and analyzed the relationships of detected haplotypes on the phylogenetic network (Fig. 3, STR haplotypes are available from authors upon request).

(…) On the network, different clusters are associated with different populations: most Mongols belong to F13686, all Evenks of the Amur River region with this haplogroup form a subcluster within F13686, and part of Upper Nanais is the basis of cluster B470.

ulchi-y-chromosome
Frequencies of haplogroups of Y chromosome in the Ulchi population. The nomenclature of haplogroups is given according to [9]. Markers that are not in bold type were not typed, but are ancestral for these nodes.

An estimate of the age of the entire haplogroup C-F12355 obtained from the data of genome-wide sequencing of seven specimens is 2400 ± 500 years (O.P. Balanovsky, unpublished data). That is, the common ancestor of all the studied representatives of various peoples with this haplogroup lived not so long ago, the first millennium BC. The formation time of cluster F13686 is somewhat later: 1990 ± 600 years.

(…) obvious traces of the interaction of the gene pool of the Ulchi with neighboring and remote peoples of the Far East and Central Asia in the time range of the last one to three thousand years were revealed. This shows that the results of work [4] on the similarity of the gene pool of the ancient (age of 7500 years) Neolithic genomes of the Amur River region to the Ulchi probably indicate not the uniqueness of the Ulchi, but the fact that this ancient gene pool was preserved in a vast circle of populations of the Far East interwoven with gene flows both with each other and, to a lesser extent, with populations of Central Asia.

The expansion of C2b1a2a-M86 (among many basal C2-M217 samples) is thus possibly associated with the spread of Tungusic, which puts C2b1a at the root of the Micro-Altaic expansion, with a formation date ca. 12700 BC, TMRCA 12500 BC (and not only Mongolian). This shows that Micro-Altaic is connected with a local population which shows a clear continuity since at least 3500 BC. This, however, tells us little about the origin of the language.

See also the recent ISBA presentation on the Houtaomuga site, Neolithic transition in Northeast Asia; and also Bronze Age population dynamics and rise of dairy pastoralism in Mongolia, Impact of colonization in north-eastern Siberia

That leaves the ancestral N lineages found among Far East Asians as Palaeo-Siberian in origin, and their late expansions to the west not particularly linked with any of the known Palaeo-Siberian ethnolinguistic groups, let alone a supposed “Uralo-Altaic” language…

Related

Dzudzuana, Sidelkino, and the Caucasus contribution to the Pontic-Caspian steppe

hunter-gatherer-pottery

It has been known for a long time that the Caucasus must have hosted many (at least partially) isolated populations, probably helped by geographical boundaries, setting it apart from open Eurasian areas.

David Reich writes in his book the following about India:

The genetic data told a clear story. Around a third of Indian groups experienced population bottlenecks as strong or stronger than the ones that occurred among Finns or Ashkenazi Jews. We later confirmed this finding in an even larger dataset that we collected working with Thangaraj: genetic data from more than 250 jati groups spread throughout India (…)

Rather than an invention of colonialism as Dirks suggested, long-term endogamy as embodied in India today in the institution of caste has been overwhelmingly important for millennia. (…)

The Han Chinese are truly a large population. They have been mixing freely for thousands of years. In contrast, there are few if any Indian groups that are demographically very large, and the degree of genetic differentiation among Indian jati groups living side by side in the same village is typically two to three times higher than the genetic differentiation between northern and southern Europeans. The truth is that India is composed of a large number of small populations.

There is little doubt now, based on findings spanning thousands of years, that the Mesolithic and Neolithic Caucasus hosted various very small populations, even if the ancestral components may be reduced to the few known to date (such as ANE, EHG, AME*, ENA, CHG, and other “deep” ancestral components).

NOTE. I will call the ancestral component of Dzudzuana/Anatolian hunter-gatherers Ancient Middle Easterner (AME), to give a clear idea of its likely extension during the Late Upper Palaeolithic, and to avoid using the more simplistic Dzudzuana, unless it is useful to mention these specific local samples.

dzudzuana-pca
Image modified from Lazaridis et al. (2018), including Caucasus, Don-Volga-Ural, and North Pontic Mesolithic-Neolithic populations. “Ancient West Eurasian population structure. (a) Geographical distribution of key ancient West Eurasian populations. (b) Temporal distribution of key ancient West Eurasian populations (approximate date in ky BP). (c) PCA of key ancient West Eurasians, including additional populations (shown with grey shells), in the space of outgroup f4-statistics (Methods).”

Genetic labs have a strong fixation with ancestry. I guess the use of complex statistical methods gives professionals and laymen alike the feeling of dealing with “Science”, as opposed to academic fields where you have to interpret data. I think language reveals a lot about the way people think, and the fact that ancestral components are called ‘lineages’ – while not wrong per se – is a clear symptom of the lack of interest in the true lineages: Y-DNA haplogroups.

Y-DNA bottlenecks

It has become quite clear that male-biased migrations are often the ones which can be confidently followed for actual population movements and ethnolinguistic identification, at least until the Iron Age. The frequently used Palaeolithic clusters offer a clear example of why ancestry does not represent what some people believe: They merely give a basic idea of sizeable population replacements by distant peoples.

Both concepts are important: sizeable and distant peoples. For example, during the Upper Palaeolithic in Europe there was a sizeable population replacement of the Aurignacian Goyet cluster by the Gravettian Vestonice cluster (probably from populations of far eastern Russia) coupled with the arrival of haplogroup I, although during the thousands of years that this material culture lasted, the previously expanded C1a2 lineages did not disappear, and there were probably different resurgence and admixture events.

Haplogroup I certainly expanded with the Gravettian culture to Iberia, where the Goyet ancestry did not change much – probably because of male-driven migrations -, to the extent that during the Magdalenian expansions haplogroup I expanded with an ancestry closer to Goyet, in what is called a ‘resurge’ of the Goyet cluster – even though there is a clear replacement of male lines.

The Villabruna (WHG) cluster is another good example. It probably spread with haplogroup R1b-L754, which – based on the extra ‘East Asian’ affinity of some samples and on modern samples from the Middle East – came probably from the east through a southern route, and not too long before the expansion of WHG likely from around the Black Sea, although this is still unclear. The finding of haplogroup I in samples of mostly WHG ancestry could confuse people that do not care about timing, sub-structured populations, and gene flow.

palaeolithic-expansions-reich
Image from David Reich’s Who We Are and How We Got Here. Having migrated out of Africa and the Near East, modern human pioneer populations spread throughout Eurasia (1). By at least thirty-nine thousand years ago, one group founded a lineage of European hunter-gatherers that persisted largely uninterrupted for more than twenty thousand years (2). Eventually, groups derived from an eastern branch of this founding population of European huntergatherers spread west (3), displaced previous groups, and were eventually themselves pushed out of northern Europe by the spread of glacial ice, shown at its maximum extent (top right). As the glaciers receded, western Europe was repeopled from the southwest (4) by a population that had managed to persist for tens of thousands of years and was related to an approximately thirty-five-thousand-year old individual from far western Europe. A later human migration, following the first strong warming period, had an even larger impact, with a spread from the southeast (5) that not only transformed the population of western Europe but also homogenized the populations of Europe and the Near East. At a single site—Goyet Caves in Belgium—ancient DNA from individuals spread over twenty thousand years reflects these transformations, with representatives from the Aurignacian, Gravettian, and Magdalenian periods.

NOTE. If you don’t understand why ‘clusters’ that span thousands of years don’t really matter for the many Palaeolithic population expansions that certainly happened among hunter-gatherers in Europe, just take a look at what happened with Bell Beakers expanding from Yamna into western Europe within 500 years.

If we don’t thread carefully when talking about population migrations, these terms are bound to confuse people. Just as the fixation on “steppe ancestry” – which marks the arrival in Chalcolithic Europe of peoples from the Pontic-Caspian region – has confused a lot of researchers to this day.

When I began to write about the Indo-European demic diffusion model, my concern was to find a single spot where a North-West Indo-European proto-language could have expanded from ca. 2000 BC (our most common guesstimate). Based on the 2015 papers, and in spite of their conclusions, I thought it had become clear that Corded Ware was not it, and it was rather Bell Beakers. I assumed that Uralic was spoken to the north (as was the traditional belief), and thus Corded Ware expanded from the forest zone, hence steppe ancestry would also be found there with other R1a lineages.

With the publication of Mathieson et al. (2017) and Olalde et al. (2017), I changed my mind, seeing how “steppe ancestry” did in fact appear quite late, hence it was likely to be the result of very specific population movements, probably directly from the Caucasus. Later, Mathieson published in a revision the sample from Alexandria of hg R1a-M417 (probably R1a-Z645, possibly Z93+), which further supported the idea that the migration of Corded Ware peoples started near the North Pontic forest-steppe (as I included in a the next revision).

The question remains the same I repeated recently, though: where do the extra Caucasus components (i.e. beyond EHG) of Eneolithic Ukraine/Corded Ware and Khvalynsk/Yamna come from?

Steppe ancestry: “EHG” + “CHG”?

About EHG ancestry

From Lazaridis et al. (2018):

Considering 2-way mixtures, we can model Karelia_HG as deriving 34 ± 2.8% of its ancestry from a Villabruna-related source, with the remainder mainly from ANE represented by the AfontovaGora3 (AG3) sample from Lake Baikal ~17kya.

AG3 was likely of haplogroup Q1a (as reported by YFull, see Genetiker), and probably the ANE ancestry found in Eastern Europe accompanied a Palaeolithic migration of Q1a2-M25 (formed ca. 22600 BC, TMRCA ca. 14300 BC).

NOTE. You can read more about the expansion of Q lineages during the Palaeolithic.

Combined with what we know about the Eneolithic Steppe and Caucasus populations – it is likely that ANE ancestry remained the most important component of some of the small ghost populations of the Caucasus until their emergence with the Lola culture.

pca-caucasus-dzudzuana
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. To understand the drawn potential Caucasus Mesolithic cluster, see above the PCA from Lazaridis et al. (2018).

The first sample we have now attributed to the EHG cluster is Sidelkino, from the Samara region (ca. 9300 BC), mtDNA U5a2. In Damgaard et al. (Science 2018), Yamnaya could be modelled as a CHG population related to Kotias Klde (54%) and the remaining from ANE population related to Sidelkino (>46%), with the following split events:

  1. 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).
  2. A split event, where the ANE component of Yamnaya splits from Sidelkino. This was inferred at about about 11 kya.
  3. 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.
  4. An ancestral split event between the CHG and ANE ancestral populations. This was inferred to occur around 40 kya.

Other samples classified as of the EHG cluster:

  • Popovo2 (ca. 6250 BC) of hg J1, mtDNA U4d – Po2 and Po4 from the same site (ca. 6550 BC) show continuity of mtDNA.
  • Karelia_HG, from Juzhnii Oleni Ostrov (ca. 6300 BC): I0211/UzOO40 (ca. 6300 BC) of hg J1(xJ1a), mtDNA U4a; and I0061/UzOO74 of hg R1a1(xR1a1a), mtDNA C1
  • UzOO77 and UzOO76 from Juzhnii Oleni Ostrov (ca. 5250 BC) of mtDNA R1b.
  • Samara_HG from Lebyanzhinka (ca. 5600 BC) of hg R1b1a, mtDNA U5a1d.

From the analysis of Lazaridis et al. (2018), we have some details about their admixture:

dzudzuana-admixture-sidelkino
Image modified from Lazaridis et al. (2018). Modeling present-day and ancient West-Eurasians. Mixture proportions computed with qpAdm (Supplementary Information section 4). The proportion of ‘Mbuti’ ancestry represents the total of ‘Deep’ ancestry from lineages that split prior to the split of Ust’Ishim, Tianyuan, and West Eurasians and can include both ‘Basal Eurasian’ and other (e.g., Sub-Saharan African) ancestry. (Left) ‘Conservative’ estimates. Each population 367 cannot be modeled with fewer admixture events than shown. (Right) ‘Speculative’ estimates. The highest number of sources (≤5) with admixture estimates within [0,1] are shown for each population. Some of the admixture proportions are not significantly different from 0 (Supplementary Information section 4).

About Anatolia_Neolithic ancestry

About the enigmatic Anatolia_Neolithic-related ancestry found in Pontic-Caspian steppe samples, this is what Wang et al. (2018) had to say:

We focused on model of mixture of proximal sources such as CHG and Anatolian Chalcolithic for all six groups of the Caucasus cluster (Eneolithic Caucasus, Maykop and Late Makyop, Maykop-Novosvobodnaya, Kura-Araxes, and Dolmen LBA), with admixture proportions on a genetic cline of 40-72% Anatolian Chalcolithic related and 28-60% CHG related (Supplementary Table 7). When we explored Romania_EN and Greece_Neolithic individuals as alternative southeast European sources (30-46% and 36-49%), the CHG proportions increased to 54-70% and 51-64%, respectively. We hypothesize that alternative models, replacing the Anatolian Chalcolithic individual with yet unsampled populations from eastern Anatolia, South Caucasus or northern Mesopotamia, would probably also provide a fit to the data from some of the tested Caucasus groups.

Also:

The first appearance of ‘Near Eastern farmer related ancestry’ in the steppe zone is evident in Steppe Maykop outliers. However, PCA results also suggest that Yamnaya and later groups of the West Eurasian steppe carry some farmer related ancestry as they are slightly shifted towards ‘European Neolithic groups’ in PC2 (Fig. 2D) compared to Eneolithic steppe. This is not the case for the preceding Eneolithic steppe individuals. The tilting cline is also confirmed by admixture f3-statistics, which provide statistically negative values for AG3 as one source and any Anatolian Neolithic related group as a second source

yamnaya-caucasus-dzudzuana
Modified image from Wang et al. (2018). In blue, Yamna-related populations. In red, Corded Ware-related populations, and two elevated Anatolia_Neolithic values in Yamna. Notice how only GAC-related admixture increases the Anatolian_N-related ancestry in the Yamna outlier from Ozero, and the late Yamna sample from Hungary, related to the homogeneous Yamna population. “Supplementary Table 14. P values of rank=3 and admixture proportions in modelling Steppe ancestry populations as a four-way admixture of distal sources EHG, CHG, Anatolian_Neolithic and WHG using 14 outgroups.Left populations: Steppe cluster, EHG, CHG, WHG, Anatolian_Neolithic. Right populations: Mbuti.DG, Ust_Ishim.DG, Kostenki14, MA1, Han.DG, Papuan.DG, Onge.DG, Villabruna, Vestonice16, ElMiron, Ethiopia_4500BP.SG, Karitiana.DG, Natufian, Iran_Ganj_Dareh_Neolithic.”

Detailed exploration via D-statistics in the form of D(EHG, steppe group; X, Mbuti) and D(Samara_Eneolithic, steppe group; X, Mbuti) show significantly negative D values for most of the steppe groups when X is a member of the Caucasus cluster or one of the Levant/Anatolia farmer-related groups (Supplementary Figs. 5 and 6). In addition, we used f- and D-statistics to explore the shared ancestry with Anatolian Neolithic as well as the reciprocal relationship between Anatolian- and Iranian farmer-related ancestry for all groups of our two main clusters and relevant adjacent regions (Supplementary Fig. 4). Here, we observe an increase in farmer-related ancestry (both Anatolian and Iranian) in our Steppe cluster, ranging from Eneolithic steppe to later groups. In Middle/Late Bronze Age groups especially to the north and east we observe a further increase of Anatolian farmer related ancestry consistent with previous studies of the Poltavka, Andronovo, Srubnaya and Sintashta groups and reflecting a different process not especially related to events in the Caucasus.

(…) Surprisingly, we found that a minimum of four streams of ancestry is needed to explain all eleven steppe ancestry groups tested, including previously published ones (Fig. 2; Supplementary Table 12). Importantly, our results show a subtle contribution of both Anatolian farmer-related ancestry and WHG-related ancestry (Fig.4; Supplementary Tables 13 and 14), which was likely contributed through Middle and Late Neolithic farming groups from adjacent regions in the West. The discovery of a quite old AME ancestry has rendered this probably unnecessary, because this admixture from an Anatolian-like ghost population could be driven even by small populations from the Caucasus.

yamna-caucasus-cwc-anatolia-neolithic
Image modified from Wang et al. (2018). Marked are: 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 1128 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 (see also Supplementary Tables 10, 14 and 20).”

NOTE. For a detailed account of the possibilities regarding this differential admixture in the North Pontic area in contrast to the Don-Volga-Ural region, you can read the posts Sredni Stog, Proto-Corded Ware, and their “steppe admixture”, and Corded Ware culture origins: The Final Frontier.

While it is not yet fully clear, the increased Anatolian_Neolithic-like ancestry in Ukraine_Eneolithic samples (see below) makes it unlikely that all such ancestry in Corded Ware groups comes from a GAC-related contribution. It is likely that at least part of it represents contributions from populations of the Caucasus, based on the mostly westward population movements in the steppe from ca. 4600 BC on, including the Suvorovo-Novodanilovka expansion, and especially the Kuban-Maykop expansion during the final Eneolithic into the North Pontic area.

NOTE. Since CHG-like groups from the Caucasus may have combinations of AME and ANE ancestry similar to Yamna (which may thus appear as ‘steppe ancestry’ in the North Pontic area), it is impossible to interpret with precision the following ADMIXTURE graphic:

ukraine-whg-ehg-steppe
Modified image from Mathieson et al. (2018). Supervised ADMIXTURE analysis, modelling each ancient individual (one per row) as a mixture of population clusters constrained to contain northwestern-Anatolian Neolithic (grey), Yamnaya from Samara (yellow), EHG (pink) and WHG (green) populations. Dates in parentheses indicate approximate range of individuals in each population.

North-Eastern Technocomplex

The East Asian contribution to samples from the WHG samples (like Loschbour or La Braña), as specified in Fu et al. (2016), does not seem to be related to Baikal_EN, and appears possibly (in the ADMIXTURE analysis) integrated into he Villabruna component. I guess this implies that the shared alleles with East Asians are quite early, and potentially due to the expansion of R1b-L754 from the East.

It would be interesting to know the specific material culture Sidelkino belonged to – i.e. if it was related to the expansion of the North-Eastern Technocomplex – , and its Y-DNA. The Post-Swiderian expansion into eastern Europe, probably associated with the expansion of R1b-P297 lineages (including R1b-M73, found later in Botai and in Baltic HG) is supposed to have begun during the 11th millennium BC, but migrations to the Urals and beyond are probably concentrated in the 9th millennium, so this sample is possibly slightly early for R1b.

NOTE. User Rozenfeld at Anthrogenica posted this, which I think is interesting (in case anyone wants to try a Y-SNP call):

there is something strange with Sidelkino EHG: first, its archaeological context is not described in the supplementary. Second, its sex is not listed in the supplementary tables. Third, after looking for info about this sample, I found that: “Сиделькино-3. Для снятия вопроса о половой принадлежности индивида была проведена генетическая экспертиза, выявившая принадлежность останков мужчине.”(translation: Sidelkino-3. To resolve the question about sex of the remains, the genetic analysis was conducted, which showed that remains belonged to male), source: http://static.iea.ras.ru/books/7487_Traditsii.pdf

So either they haven’t mentioned his Y-DNA in the paper for some reason, or there are more than one Sidelkino sample and the male one has not yet been published. The coverage of the Sidelkino sample from the paper is 2.9, more than enough to tell Y-DNA haplogroup.

zaliznyak-post-swiderian
The map of spreading of Post-Swiderian and Post-Krasnosillian sites in Mesolithic of Eastern Europe in the 8th millennia BC. From Zaliznyak (see here).

My speculative guess right now about specific population movements in far eastern Europe, based on the few data we have:

  • The expansion of the North-Eastern Technocomplex first around the 9th millennium BC, most likely expanded R1b-P279 ca. 11300 BC, judging by its TMRCA, with both R1b-M73 (TMRCA 5300) and R1b-M269 (TMRCA 4400 BC) info (with extra El Mirón ancestry) back, and thus Eurasiatic.
  • The expansion of haplogroup J1 to the north may have happened before or after the R1b-P279 expansion. Judging by the increase in AG3-related ancestry near Karelia compared to Baltic_HG, it is possible that it expanded just after R1b-P279 (hence possibly J1-Y6304? TMRCA 9700 BC). Its long-lasting presence in the Caucasus is supported by the Satsurblia (ca. 11300 BC) and the Dolmen BA (ca. 1300 BC) samples.
  • The expansion of R1a-M17 ca. 6600 BC is still likely to have happened from the east, based on the R1a-M17 samples found in Baikalic cultures slightly later (ca. 5300 BC). The presence of elevated Baikal_EN ancestry in Karelia HG and in Samara HG, and the finding of R1a-M417 samples in the Forest Zone after the Mesolithic suggests a connection with the expansion of Hunter-Gatherer pottery, from the Elshanka culture in the Samara region northward into the Forset Zone and westward into the North Pontic area.
  • The expansion of R1b-M73 ca. 5300 BC is likely to be associated with the emergence of a group east of the Urals (related to the later Botai culture, and potentially Pre-Yukaghir). Its presence in a Narva sample from Donkalnis (ca. 5200 BC) suggest either an early split and spread of both R1b-P297 lineages (M73 and M269) through Eastern Europe, or maybe a back-migration with hunter-gatherer pottery.
  • R1b-M269 spread successfully ca. 4400 BC (and R1b-L23 ca. 4100 BC, both based on TMRCA), and this successful expansion is probably to be associated with the Khvalynsk-Novodanilovka expansion. We already know that Samara_HG ca. 5600 was R1b1a, so it is likely that R1b-M269 appeared (or ‘resurged’) in the Volga-Ural region shortly after the expansion of R1a-M17, whose expansion through the region may be inferred by the additional AG3 and Baikal_EN ancestry. Interesting from Samara_HG compared to the previous Sidelkino sample is the introduction of more El Mirón-related ancestry, typical of WHG populations (and thus proper of Baltic groups).

NOTE. The TMRCA dates are obviously gross approximations, because a) the actual rate of mutation is unknown and b) TMRCA estimates are based on the convergence of lineages that survived. The potential finding of R1a-Z645 (possibly Z93+) in Ukraine Eneolithic (ca. 4000 BC), and the potential finding of R1b-L23 in Khvalynsk ca. 4250 BC complicates things further, in terms of dates and origins of any subclade.

The question thus remains as it was long ago: did R1b-M269 lineages expand (‘return’) from the east, near the Urals, or directly from the north? Were they already near Samara at the same time as the expansion of hunter-gatherer pottery, and were not much affected by it? Or did they ‘resurge’ from populations admixed with Caucasus-related ancestry after the expansion of R1a-M17 with this pottery (since there are different stepped expansions from the Samara region)? We could even ask, did R1a-M17 really expand from the east, i.e. are the dates on Baikalic subclades from Moussa et al. (2016) reliable? Or did R1a-M17 expand from some pockets in the Pontic-Caspian steppe, taking over the expansion of HG pottery at some point?

hunger-gatherer-pottery
Early Neolithic cultures in eastern and central Europe: 1–Yelshanian; 2–North Caspian; 3–Rakushechnyj Yar; 4–Surskian; 5–Dnieper-Donetsian; 6– Bug-Dniesterian; 7–Upper Volga; 8–Narvian; 9–Linear Pottery. White arrows: expansion of early farming; black arrows: spread of pottery-making traditions. From Dolukhanov et al. (2009).

Maglemose-related migrations

The most interesting aspect from the new paper (regarding Indo-Uralic migrations) is that Ancestral Middle Easterner ancestry will probably be a better proxy for the Anatolia_Neolithic component found in Ukraine Mesolithic to Eneolithic, and possibly also for some of the “more CHG-like” component found among Pontic-Caspian steppe populations, all likely derived from different admixture events with groups from the Caucasus.

NOTE. Even the supposed gene flow of Neolithic Iranian ancestry into the Caucasus can be put into question, since that means possibly a Dzudzuana-like population with greater “deep ancestry” proportion than the one found in CHG, which may still be found within the Caucasus.

If it was not clear already that following ‘steppe ancestry’ wherever it appears is a rather lame way of following Indo-European migrations, every single sample from the Caucasus and their admixture with Pontic-Caspian steppe populations will probably show that “steppe ancestry” is in fact formed by a variety of steppe-related ancestral components, impossible to follow coherently with a single population. Exactly what is happening already with the Siberian ancestry.

If the paper on the Dzudzuana samples has shown something, is that the expansion of an ANE-like population shook the entire Caucasus area up to the Zagros Mountains, creating this ANE – AME cline that are CHG and Iran_N, with further contributions of “deep ancestries” (probably from the south) complicating the picture further.

If this happens with few known samples, and we know of an ANE-like ghost population in the Caucasus (appearing later in the Lola culture), we can already guess that the often repeated “CHG component” found in Ukraine_Eneolithic and Khvalynsk will not be the same (except the part mediated by the Novodanilovka expansion).

This ANE-like expansion happened probably in the Late Upper Palaeolithic, and reached Northern Europe probably after the expansion of the Villabruna cluster (ca. 12000 BC), judging by the advance of AG3-like and ENA-like ancestry in later WHG samples.

The population movements during the Mesolithic and Early Neolithic in the North Pontic area are quite complicated: the extra AME ancestry is probably connected to the admixture with populations from the Caucasus, while the close similarity of Ukraine populations with Scandinavian ones (with an increase in Villabruna ancestry from Mesolithic to Neolithic samples), probably reveal population movements related to the expansion of Maglemose-related groups.

maglemose-mesolithic
Etno-cultural situation in Central and Eastern Europe in the Late Mesolithic — Early Neolithic (VI—V Mill. BC) (after Конча 2004: 201, карта 1; made after ideas by L. L. Zaliznyak). Legend: 1 — Maglemose circle in the VII Mill. BC (after Gr. Clark); 2—7 — Mesolithic cultures of the Post-Maglemose tradition, VI Mill. BC (after S. Kozłowsky, L. L. Zaliznyak): 2 — de Leyen-Wartena; 3 — Oldesloe — Godenaa; 4 — Chojnice — Peńki; 5 — Janisłavice; 6 — finds of Janisłavice artefacts outside of the main area; 7 — Donets culture; 8 — directions of the settling of Janisłavice people (after S. Kozłowsky and L. L. Zaliznyak); 9 — the south border of Mesolithic and Early Neolithic cultures of post-Swidrian and post-Arensburgian traditions; 10 — northern border of settlement of the Balkan-Danubian farmers; 11 — Bug- Dniester culture; 12 — Neolithic cultures emerged on the ethno-cultural basis of post-Maglemose: Э — Ertebölle-Ellerbeck, Н — Neman, Д — Dnieper-Donets, М — Mariupol (western variants). From Klein (2017).

These Maglemose-related groups were probably migrants from the north-west, originally from the Northern European Plains, who occupied the previous Swiderian territory, and then expanded into the North Pontic area. The overwhelming presence of I2a (likely all I2a2a1b1b) lineages in Ukraine Neolithic supports this migration.

The likely picture of Mesolithic-Neolithic migrations in the North Pontic area right now is then:

  1. Expansion of R1a-M459 from the east ca. 12000 BC – probably coupled with AG3 and also some Baikal_EN ancestry. First sample is I1819 from Vasilievka (ca. 8700 BC), another is from Dereivka ca. 6900 BC.
  2. Expansion of R1b-V88 from the Balkans in the west ca. 9700 BC, based on its TMRCA and also the Balkan hunter-gatherer population overwhemingly of this haplogroup from the 10th millennium until the Neolithic. First sample is I1734 from Vasilievka (ca. 7252 BC), which suggests that it replaced the male population there, based on their similar EHG-like adxmixture (and lack of sizeable WHG increase), and shared mtDNA U5b2, U5a2.
  3. Expansion of I2a-Y5606 probably ca. 6800 based on its TMRCA with Janislawice culture. Supporting this is the increase in WHG contribution to Neolithic samples, including the spread of U4 subclades compared to the previous period.
  4. Expansion of R1a-M17 starting probably ca. 6600 BC in the east (see above).

NOTE. The first sample of haplogroup I appears in the Mesolithic: I1763 (ca. 8100 BC) of haplogroup I2a1, probably related to an older Upper Palaeolithic expansion.

janislawice
Distribution of archeological cultures in the North Pontic Region during the Mesolithic (7th – 6th millennium BCE). Dotted, dashed and solid lines with corresponding arrows indicate alternative models of the spread of the Grebenyky culture groups. (After Bryuako IV., Samojlova TL., Eds, Drevnie kul’tury Severo-­‐Zapadnogo Prichernomor’ya, Odessa: SMIL, 2013.) Nikitin – Ivanova 2017.

Conclusion

It is becoming more and more clear with each new paper that – unless the number of very ancient samples increases – the use of Y-chromosome haplogroups remains one of the most important tools for academics; this is especially so in the steppes, in light of the diversity found in populations from the Caucasus. A clear example comes from the Yamna – Corded Ware similarities:

After the publication of the 2015 papers, it was likely that Yamna expanded with haplogroup R1b-L23, but it has only become crystal clear that Yamna expanded through the steppes into Bell Beakers, now that we have data about the strict genetic homogeneity of the whole Yamna population from west to east (including Afanasevo), in contrast with contemporary Corded Ware peoples which expanded from a different forest-steppe population.

The presence of haplogroups Q and R1a-M459 (xM17) in Khvalynsk along with a R1b1a sample, which some interpreted as being akin to modern ‘mixed’ populations in the past, is likely to point instead to a period of Khvalynsk-Novodanilovka expansion with R1b-M269, where different small populations from the steppe were being integrated into the common Khvalynsk stock, but where differences are seen in material culture surrounding their burials, as supported by the finding of R1b1 in the Kuban area already in the first half of the 5th millennium. The case would be similar to the early ‘mixed’ Icelandic population.

Only after the emergence of the Samara culture (in the second half of the 6th millennium BC), with a sample of haplogroup R1b1a, starts then the obvious connection with Early Proto-Indo-Europeans; and only after the appearance of late Sredni Stog and haplogroup R1a-M417 (ca. 4000 BC) is its connection with Uralic also clear. In previous population movements, I think more haplogroups were involved in migrations of small groups, and only some communities among them were eventually successful, expanding to be dominant, creating ever growing cultures during their expansions.

Indeed, if you think in terms of Uralic and Indo-European just as converging languages, and forget their potential genetic connection, then the genetic + linguistic picture becomes simplified, and the upper frontier of the 6th millennium BC with a division North Pontic (Mariupol) vs. Volga-Ural (Samara) is enough. However, tracing their movements backwards – with cultural expansions from west to east (with the expansion of farming), and earlier east to west (with hunter-gatherer pottery), and still earlier west to east (with the north-eastern technocomplex), offers an interesting way to prove their potential connection to macrofamilies, at least in terms of population movements.

corded-ware-uralic-qpgraph
Modified image from Tambets et al. (2018) Proportions of ancestral components in studied European and Siberian populations and the tested qpGraph model. a 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. The NeolL (Neolithic Levant) ancestry selected in this qpGraph is likely to correspond (at least in part) to a specific Dzudzuana-like component present in the CHG-like population that admixed in the North Pontic area.

I am quite convinced right now that it would be possible to connect the expansion of R1b-L754 subclades with a speculative Nostratic (given the R1b-V88 connection with Afroasiatic, and the obvious connection of R1b-L297 with Eurasiatic). Paradoxically, the connection of an Indo-Uralic community in the steppes (after the separation of Yukaghir) with any lineage expansion (R1a-M17, R1b-M269, or even Q, I or J1) seems somehow blurrier than one year ago, possibly just because there are too many open possibilities.

David Reich says about the admixture with Neanderthals, which he helped discover:

At the conclusion of the Neanderthal genome project, I am still amazed by the surprises we encountered. Having found the first evidence of interbreeding between Neanderthals and modern humans, I continue to have nightmares that the finding is some kind of mistake. But the data are sternly consistent: the evidence for Neanderthal interbreeding turns out to be everywhere. As we continue to do genetic work, we keep encountering more and more patterns that reflect the extraordinary impact this interbreeding has had on the genomes of people living today.

I think this is a shared feeling among many of us who have made proposals about anything, to fear that we have made a gross, evident mistake, and constantly look for flaws. However, it seems to me that geneticists are more preoccupied with being wrong in their developed statistical methods, in the theoretical models they are creating, and not so much about errors in the true ancient ethnolinguistic picture human population genetics is (at least in theory) concerned about. Their publications are, after all, constantly associating genetic finds with cultures and (whenever possible) languages, so this aspect of their research should not be taken lightly.

Seeing how David Anthony or Razib Khan (among many others) have changed their previously preferred migration models as new data was published, and they continue to be respected in their own fields, I guess we can be confident that professionals with integrity are going to accept whatever new picture appears. While I don’t think that genetic finds can change what we can reconstruct with comparative grammar, I am also ready to revise guesstimates and routes of expansion of certain dialects if R1a-Z645 is shown to have accompanied Late Proto-Indo-Europeans during their expansion with Yamna, and later integrated somehow with Corded Ware.

However, taking into account the obsession of some with an ancestral, uninterrupted R1a—Indo-European association, and the lack of actual political repercussion of Neanderthal admixture, I think the most common nightmare that all genetic researchers should be worried about is to keep inflating this “Yamnaya ancestry”-based hornet’s nest, which has been constantly stirred up for the past two years, by rejecting it – or, rather, specifying it into its true complex nature.

This succession of corrections and redefinitions, coupled with the distinct Y-DNA bottleneck of each steppe population, will eventually lead to a completely different ethnolinguistic picture of the Pontic-Caspian region during the Eneolithic, which is likely to eventually piss off not only reasonable academics stubbornly attached to the CWC-IE idea, but also a part of those interested in daydreaming about their patrilineal ancestors.

Sometimes it’s better to just rip off the band-aid once and for all…

Featured image from The oldest pottery in hunter-gatherer communitiesand models of Neolithisation of Eastern Europe (2015), by Andrey Mazurkevich and Ekaterina Dolbunova.

Related

Interesting is today’s post in Ancient DNA Era: Is Male-driven Genetic Replacement always meaning Language-shift?

Corded Ware—Uralic (II): Finno-Permic and the expansion of N-L392/Siberian ancestry

finno-ugric-samoyedic

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

I read from time to time that “we have not sampled Uralic speakers yet”, and “we are waiting to see when Uralic-speaking peoples are sampled”. Are we, though?

Proto-language homelands are based on linguistic data, such as guesstimates for dialectal evolution, loanwords and phonetic changes for language contacts, toponymy for ancient territories, etc. depending on the available information. The trace is then followed back, using available archaeological data, from the known historic speakers and territory to the appropriate potential prehistoric cultures. Only then can genetic analyses help us clarify the precise prehistoric population movements that better fit the models.

uralic-language-family
The traditional family tree of the Uralic branches. Kallio (2014)

The linguistic homeland

We thought – using linguistic guesstimates and fitting prehistoric cultures and their expansion – that Yamna was the Late Proto-Indo-European culture, so when Yamna was sampled, we had Late Proto-Indo-Europeans sampled. Simple deduction.

We thought that north-eastern Europe was a Uralic-speaking area during the Neolithic:

  • For those supporting a western continuity (and assuming CWC was Indo-European), the language was present at least since the Comb Ware culture, potentially since the Mesolithic.
  • For those supporting a late introduction into Finland, Uralic expanded the latest with Abashevo-related movements after its incorporation of Volosovo and related hunter-gatherers.

The expansion to the east must have happened through progressive infiltrations with Seima-Turbino / Andronovo-related expansions.

uralic-time-space
Some datings for the traditional proto-stages from Uralic to Finnic. Kallio (2014).

Finding the linguistic homeland going backwards can be described today as follows:

I. Proto-Fennic homeland

Based on the number of Baltic loanwords, not attested in the more eastern Uralic branches (and reaching only partially Mordvinic), the following can be said about western Finno-Permic languages (Junttila 2014):

The Volga-Kama Basin lies still too far east to be included in a list of possible contact locations. Instead, we could look for the contact area somewhere between Estonia in the west and the surroundings of Moscow in the east, a zone with evidence of Uralic settlement in the north and Baltic on the south side.

The only linguistically well-grounded version of the Stone Age continuation theory was presented by Mikko Korhonen in 1976. Its validity, however, became heavily threatened when Koivulehto 1983a-b proved the existence of a Late Proto-Indo-European or Pre-Baltic loanword layer in Saami, Finnic, and Mordvinic. Since this layer must precede the Baltic one and it was presumably acquired in the Baltic Sea region, Koivulehto posited it on the horizon of the Battle Axe period. This forces a later dating for the Baltic–Finnic contacts.

Today the Battle Axe culture is dated at 3200 to 3000 BC, a period far too remote to correspond linguistically with Proto-Baltic (Kallio 1998a).

Since the Baltic contacts began at a very initial phase of Proto-Finnic, the language must have been relatively uniform at that time. Hence, if we consider that the layer of Baltic loanwords may have spread over the Gulf of Finland at that time, we could also insist that the whole of the Proto-Finnic language did so.

migration-theory
Prehistoric Balts as the southern neighbours of Proto-Finnic speakers. 1 = The approximated area of Proto-Uralic. 2 = The approximated area of Finnic during the Iron Age. 3 = The area of ancient Baltic hydronyms. 4 = The area of Baltic languages in about 1200 AD. 5 = The problem: When did Uralic expand westwards and when did it meet Baltic? Junntila (2012).

II. Proto-Finno-Saamic homeland

The evidence of continued Palaeo-Germanic loanwords (from Pre- to Proto-Germanic stages) is certainly the most important data to locate the Finno-Saamic homeland, and from there backwards into the true Uralic homeland. Following Kallio (2017):

(…) the loanword evidence furthermore suggests that the ancestors of Finnic and Saamic had at least phonologically remained very close to Proto-Uralic as late as the Bronze Age (ca. 1700–500 BC). In particular, certain loanwords, whose Baltic and Germanic sources point to the first millennium BC, after all go back to the Finno-Saamic proto-stage, which is phonologically almost identical to the Uralic proto-stage (see especially the table in Sammallahti 1998: 198–202). This being the case, Dahl’s wave model could perhaps have some use in Uralic linguistics, too.

The presence of Pre-Germanic loanwords points rather to the centuries around the turn of the 2nd – 1st millennium BC or earlier. Proto-Germanic words must have been borrowed before the end of Germanic influence in the eastern Baltic at the beginning of the Iron Age, which sets a clear terminus ante quem ca. 800 BC.

The arrival of Bell Beaker peoples in Scandinavia ca. 2350 BC, heralding the formation of the Dagger Period, as well as the development of Pre-Germanic in common with Finnic-like populations point to the late 3rd / early 2nd millennium BC as the first time of close interaction through the Baltic region.

III. Proto-Uralic homeland

(…) the earliest Indo-European loanwords in the Uralic languages (…) show that Proto-Uralic cannot have been spoken much earlier than Proto-Indo-European dated about 3500 BC (Koivulehto 2001: 235, 257). As the same loanword evidence naturally also shows that the Uralic and Indo-European homelands were not located far from one another, the Uralic homeland can most likely be located in the Middle and Upper Volga region, right north of the Indo-European homeland*. From the beginning of the Subneolithic period about 5900 BC onwards, this region was an important innovation centre, from where several cultural waves spread to the Finnish Gulf area, such as the Sperrings Ware wave about 4900 BC, the Combed Ware wave about 3900 BC, and the Netted Ware wave about 1900 BC (Carpelan & Parpola 2001: 78–90).

The mainstream position is nowadays trying to hold together the traditional views of Corded Ware as Indo-European, and a Uralic Fennoscandia during the Bronze Age.

The following is an example of how this “Volosovo/Forest Zone hunter-gatherer theory” of Uralic origins looks like, as a ‘mixture’ of cultures and languages that benefits from the lack of genetic data for certain regions and periods (taken from Parpola 2018):

asbestos-ware
The extent of Typical Comb Ware (TCW), Asbestos- and Organic-tempered Wares (AOW) and Volosovo and Garino-Bor cultures; areas with deposits of native copper in Karelia and copperbearing sandstone in Volga-Kama-area are marked dark gray (after Zhuravlev 1977; Krajnov 1987; Nagovitsyn 1987; Chernykh 1992; Carpelan 1999; Zhul´nikov 1999). From Nordqvist et al. (2012).

The Corded Ware (or Battle Axe) culture intruded into the Eastern Baltic and coastal Finland already around 3100 BCE. The continuity hypothesis maintains that the early Proto-Finnic speakers of the coastal regions, who had come to Finland in the 4th millennium BCE with the Comb-Pitted Ware, coexisted with the Corded Ware newcomers, gradually adopting their pastoral culture and with it a number of NW-IE loanwords, but assimilating the immigrants linguistically.

The fusion of the Corded Ware and the local Comb-Pitted Ware culture resulted into the formation of the Kiukais culture (c. 2300–1500) of southwestern Finland, which around 2300 received some cultural impulses from Estonia, manifested in the appearance of the Western Textile Ceramic (which is different from the more easterly Textile Ceramic or Netted Ware, and which is first attested in Estonia c. 2700 BCE, cf. Kriiska & Tvauri 2007: 88), and supposed to have been accompanied by an influx of loanwords coming from Proto-Baltic. At the same time, the Kiukais culture is supposed to have spread the custom of burying chiefs in stone cairns to Estonia.

The coming of the Corded Ware people and their assimilation created a cultural and supposedly also a linguistic split in Finland, which the continuity hypothesis has interpreted to mean dividing Proto-Saami-Finnic unity into its two branches. Baltic Finnic, or simply Finnic, would have emerged in the coastal regions of Finland and in the northern East Baltic, while preforms of Saami would have been spoken in the inland parts of Finland.

The Nordic Bronze Age culture, correlated above with early Proto-Germanic, exerted a strong influence upon coastal Finland and Estonia 1600–700 BCE. Due to this, the Kiukais culture was transformed into the culture of Paimio ceramics (c. 1600–700 BCE), later continued by Morby ceramics (c. 700 BCE – 200 CE). The assumption is that clear cultural continuity was accompanied by linguistic continuity. Having assimilated the language of the Germanic traders and relatively few settlers of the Bronze Age, the language of coastal Finland is assumed to have reached the stage of Proto-Finnish at the beginning of the Christian era. In Estonia, the Paimio ceramics have a close counterpart in the contemporaneous Asva ceramics.

Eastern homelands?

I will not comment on Siberian or Central Asian homeland proposals, because they are obviously not mainstream, still less today when we know that Uralic was certainly in contact with Proto-Indo-European, and then with Pre- and Proto-Indo-Iranian, as supported even by the Copenhagen group in Damgaard et al. (2018).

This is what Kallio (2017) has to say about the agendas behind such proposals:

Interestingly, the only Uralicists who generally reject the Central Russian homeland are the Russian ones who prefer the Siberian homeland instead. Some Russians even advocate that the Central Russian homeland is only due to Finnish nationalism or, as one of them put it a bit more tactfully, “the political and ideological situation in Finland in the first decades of the 20th century” (Napolskikh 1995: 4).

Still, some Finns (and especially those who also belong to the “school who wants it large and wants it early”) simultaneously advocate that exactly the same Central Russian homeland is due to Finnlandisierung (Wiik 2001: 466).

Hence, for those of you willing to learn about fringe theories not related to North-Eastern Europe, you also have then the large and early version of the Uralic homeland, with Wiik’s Palaeolithic continuity of Uralic peoples spread over all of eastern and central Europe (hence EHG and R1a included):

atlantic-finnic-theory
Palaeolithic boat peoples and Finno-Ugric. Source

These fringe Finnish theories look a lot like the Corded Ware expansion… Better not go the Russian or Finnish nationalist ways? Agreed then, let’s discuss only rational proposals based on current data.

The archaeological homeland

For a detailed account of the Corded Ware expansion with Battle Axe, Fatyanovo-Balanovo, and Abashevo groups into the area, you can read my recent post on the origin of R1a-Z645.

1. Textile ceramics

During the 2nd millennium BC, textile impressions appear in pottery as a feature across a wide region, from the Baltic area through the Volga to the Urals, in communities that evolve from late Corded Ware groups without much external influence.

While it has been held that this style represents a north-west expansion from the Volga region (with the “Netted Ware” expansion), there are actually at least two original textile styles, one (earlier) in the Gulf of Finland, common in the Kiukainen pottery, which evolves into the Textile ware culture proper, and another which seems to have an origin in the Middle Volga region to the south-east.

The Netted ware culture is the one that apparently expands into inner Finland – a region not densely occupied by Corded Ware groups until then. There are, however, no clear boundaries between groups of both styles; textile impressions can be easily copied without much interaction or population movement; and the oldest textile ornamentation appeared on the Gulf of Finland. Hence the tradition of naming all as groups of Textile ceramics.

textile-ware-cultures
Maximum distribution of Textile ceramics during the Bronze Age (ca. 2000-800 BC). Asbestos-tempered ware lies to the north (and is also continued in western Fennoscandia).

The fact that different adjacent groups from the Gulf of Finland and Forest Zone share similar patterns making it very difficult to differentiate between ‘Netted Ware’ or ‘Textile Ware’ groups points to:

  • close cultural connections that are maintained through the Gulf of Finland and the Forest Zone after the evolution of late Corded Ware groups; and
  • no gross population movements in the original Battle Axe / Fatyanovo regions, except for the expansion of Netted Ware to inner Finland, Karelia, and the east, where the scattered Battle Axe finds and worsening climatic conditions suggest most CWC settlements disappeared at the end of the 3rd millennium BC and recovered only later.

NOTE. This lack of population movement – or at least significant replacement by external, non-CWC groups – is confirmed in genetic investigation by continuity of CWC-related lineages (see below).

The technology present in Textile ceramics is in clear contrast to local traditions of sub-Neolithic Lovozero and Pasvik cultures of asbestos-tempered pottery to the north and east, which point to a different tradition of knowledge and learning network – showing partial continuity with previous asbestos ware, since these territories host the main sources of asbestos. We have to assume that these cultures of northern and eastern Fennoscandia represent Palaeo-European (eventually also Palaeo-Siberian) groups clearly differentiated from the south.

The Chirkovo culture (ca. 1800-700 BC) forms on the middle Volga – at roughly the same time as Netted Ware formed to the west – from the fusion of Abashevo and Balanovo elites on Volosovo territory, and is also related (like Abashevo) to materials of the Seima-Turbino phenomenon.

Bronze Age ethnolinguistic groups

In the Gulf of Finland, Kiukainen evolves into the Paimio ceramics (in Finland) — Asva Ware (in Estonia) culture, which lasts from ca. 1600 to ca. 700 BC, probably representing an evolving Finno-Saamic community, while the Netted Ware from inner Finland (the Sarsa and Tomitsa groups) and the groups from the Forest Zone possibly represent a Volga-Finnic community.

NOTE. Nevertheless, the boundaries between Textile ceramic groups are far from clear, and inner Finland Netted Ware groups seem to follow a history different from Netted Ware groups from the Middle and Upper Volga, hence they could possibly be identified as an evolving Pre-Saamic community.

Based on language contacts, with Early Baltic – Early Finnic contacts starting during the Iron Age (ca. 500 BC onwards), this is a potential picture of the situation at the end of this period, when Germanic influence on the coast starts to fade, and Lusatian culture influence is stronger:

aikio-finnic-saamic
The linguistic situation in Lapland and the northern Baltic Sea Area in the Early Iron Age prior to the expansion of Saami languages; the locations of the language groups are schematic. The black line indicates the distribution of Saami languages in the 19th century, and the gray line their approximate maximal distribution before the expansion of Finnic. Aikio (2012)

The whole Finno-Permic community remains thus in close contact, allowing for the complicated picture that Kallio mentions as potentially showing Dahl’s wave model for Uralic languages.

Genetic data shows a uniform picture of these communities, with exclusively CWC-derived ancestry and haplogroups. So in Mittnik et al. (2018) all Baltic samples show R1a-Z645 subclades, while the recent session on Estonian populations in ISBA 8 (see programme in PDF) clearly states that:

[Of the 24 Bronze Age samples from stone-cist graves] all 18 Bronze Age males belong to R1a.

Regarding non-Uralic substrates found in Saami, supposedly absorbed during the expansion to the north (and thus representing languages spoken in northern Fennoscandia during the Bronze Age) this is what Aikio (2012) has to say:

The Saami substrate in the Finnish dialects thus reveals that also Lakeland Saami languages had a large number of vocabulary items of obscure origin. Most likely many of these words were substrate in Lakeland Saami, too, and ultimately derive from languages spoken in the region before Saami. In some cases the loan origin of these words is obvious due to their secondary Proto-Saami vowel combinations such as *ā–ë in *kāvë ‘bend; small bay’ and *šāpšë ‘whitefish’. This substrate can be called ‘Palaeo-Lakelandic’, in contrast to the ‘Palaeo-Laplandic’ substrate that is prominent in the lexicon of Lapland Saami. As the Lakeland Saami languages became extinct and only fragments of their lexicon can be reconstructed via elements preserved in Finnish place-names and dialectal vocabulary, we are not in a position to actually study the features of this Palaeo-Lakelandic substrate. Its existence, however, appears evident from the material above.

If we wanted to speculate further, based on the data we have now, it is very likely that two opposing groups will be found in the region:

A) The central Finnish group, in this hypothesis the Palaeo-Lakelandic group, made up of the descendants of the Mesolithic pioneers of the Komsa and Suomusjärvi cultures, and thus mainly Baltic HG / Scandinavian HG ancestry and haplogroups I / R1b(xM269) (see more on Scandinavian HG).

siberian-ancestry-map
Frequency map of the so-called ‘Siberian’ component. From Tambets et al. (2018).

B) Lapland and Kola were probably also inhabited by similar Mesolithic populations, until it was eventually assimilated by expanding Siberian groups (of Siberian ancestry and N1c-L392 lineages) from the east – entering the region likely through the Kola peninsula – , forming the Palaeo-Laplandic group, which was in turn later replaced by expanding Proto-Saamic groups.

Siberian ancestry appears first in Fennoscandia at Bolshoy Oleni Ostrov ca. 1520 BC, with haplogroup N1c-L392 (2 samples, BOO002 and BOO004), and with Siberian ancestry. This is their likely movement in north-eastern Europe, from Lamnidis et al (2018):

The large Siberian component in the Bolshoy individuals from the Kola Peninsula provides the earliest direct genetic evidence for an eastern migration into this region. Such contact is well documented in archaeology, with the introduction of asbestos-mixed Lovozero ceramics during the second millenium BC, and the spread of even-based arrowheads in Lapland from 1,900 BCE. Additionally, the nearest counterparts of Vardøy ceramics, appearing in the area around 1,600-1,300 BCE, can be found on the Taymyr peninsula, much further to the east. Finally, the Imiyakhtakhskaya culture from Yakutia spread to the Kola Peninsula during the same period.

saamic-lovozero-pca
PCA plot of 113 Modern Eurasian populations, with individuals from this study projected on the principal components. Uralic speakers are highlighted in light purple. Image modified from Lamnidis et al. (2018)

Obviously, these groups of asbestos-tempered ware are not connected to the Uralic expansion. From the same paper:

The fact that the Siberian genetic component is consistently shared among Uralic-speaking populations, with the exceptions of Hungarians and the non-Uralic speaking Russians, would make it tempting to equate this component with the spread of Uralic languages in the area. However, such a model may be overly simplistic. First, the presence of the Siberian component on the Kola Peninsula at ca. 4000 yBP predates most linguistic estimates of the spread of Uralic languages to the area. Second, as shown in our analyses, the admixture patterns found in historic and modern Uralic speakers are complex and in fact inconsistent with a single admixture event. Therefore, even if the Siberian genetic component partly spread alongside Uralic languages, it likely presented only an addition to populations carrying this component from earlier.

2. The Early Iron Age

The Ananino culture appears in the Vyatka-Kama area, famed for its metallurgy, with traditions similar to the North Pontic area, by this time developing Pre-Sauromatian traditions. It expanded to the north in the first half of the first millennium BC, remaining in contact with the steppes, as shown by the ‘Scythian’ nature of its material culture.

NOTE. The Ananino culture can be later followed through its zoomorphic styles into Iron Age Pjanoborskoi and Gljadenovskoi cultures, later to Ural-Siberian Middle Age cultures – Itkuska, Ust’-Poluiska, Kulaiska cultures –, which in turn can be related as prototypes of medieval Permian styles.

ananino-culture-homeland
Territory of (early and maximum) Ananino material culture. Vasilyev (2002).

At the same time as the Ananino culture begins to expand ca. 1000 BC, the Netted Ware tradition from the middle Oka expanded eastwards into the Oka-Vyatka interfluve of the middle Volga region, until then occupied by the Chirkovo culture. Eventually the Akozino or Akhmylovo group (ca. 800-300 BC) emerged from the area, showing a strong cultural influence from the Ananino culture, by that time already expanding into the Cis-Urals region.

The Akozino culture remains nevertheless linked to the western Forest Zone traditions, with long-ranging influences from as far as the Lusatian culture in Poland (in metallurgical techniques), which at this point is also closely related with cultures from Scandinavia (read more on genetics of the Tollense Valley).

malar-celts-ananino
Mälar celts and molds for casting (a) and the main distribution area (в) of Mälar-type celts of the Mälar type in the Volga-Kama region (according to Kuzminykh 1983: figure 92) and Scandinavia (according to Baudou 1960: Karte 10); Ananino celts and molds for casting (б) and the main distribution area (г) of the distribution of the celts of the Ananino type in the Volga-Kama area (according to Kuzminykh 1983: figure 9); dagger of Ananino type (д).Map from (Yushkova 2010)

Different materials from Akozino reach Fennoscandia late, at the end of the Bronze Age and beginning of the Early Iron Age, precisely when the influence of the Nordic Bronze Age culture on the Gulf of Finland was declining.

This is a period when Textile ceramic cultures in north-eastern Europe evolve into well-armed chiefdom-based groups, with each chiefdom including thousands or tens of thousands, with the main settlements being hill forts, and those in Fennoscandia starting ca. 1000-400 BC.

Mälar-type celts and Ananino-type celts appear simultaneously in Fennoscandia and the Forest Zone, with higher concentrations in south-eastern Sweden (Mälaren) and the Volga-Kama region, supporting the existence of a revived international trade network.

akozino-malar-axes-fennoscandia
Distribution of the Akozino-Mälar axes according to Sergej V. Kuz’minykh (1996: 8, Abb. 2).

The Paimio—Asva Ware culture evolves (ca. 700-200 BC) into the Morby (in Finland) — Ilmandu syle (in Estonia, Latvia, and Mälaren) culture. The old Paimio—Asva tradition continues side by side with the new one, showing a clear technical continuity with it, but with ornamentation compared to the Early Iron Age cultures of the Upper Volga area. This new south-eastern influence is seen especially in:

  • Akozino-Mälar axes (ca. 800-500 BC): introduced into the Baltic area in so great numbers – especially south-western Finland, the Åland islands, and the Mälaren area of eastern Sweden – that it is believed to be accompanied by a movement of warrior-traders of the Akozino-Akhmylovo culture, following the waterways that Vikings used more than a thousand years later. Rather than imports, they represent a copy made with local iron sources.
  • Tarand graves (ca. 500 BC – AD 400): these ‘mortuary houses’ appear in the coastal areas of northern and western Estonia and the islands, at the same time as similar graves in south-western Finland, eastern Sweden, northern Latvia and Courland. Similar burials are found in Akozino-Akhmylovo, with grave goods also from the upper and middle Volga region, while grave goods show continuity with Textile ware.

The use of asbestos increases in mainland Finnish wares with Kjelmøy Ware (ca. 700 BC – AD 300), which replaced the Lovozero Ware; and in the east in inner Finland and Karelia with the Luukonsaari and Sirnihta wares (ca. 700-500 BC – AD 200), where they replaced the previous Sarsa-Tomitsa ceramics.

The Gorodets culture appears during the Scythian period in the forest-steppe zone north and west of the Volga, shows fortified settlements, and there are documented incursions of Gorodets iron makers into the Samara valley, evidenced by deposits of their typical pottery and a bloom or iron in the region.

Iron Age ethnolinguistic groups

According to (Koryakova and Epimakhov 2007):

It is commonly accepted by archaeology, ethnography, and linguistics that the ancestors of the Permian peoples (the Udmurts, Komi-Permians, and Komi-Zyryans) left the sites of Ananyino cultural intercommunity.

NOTE. For more information on the Late Metal Ages and Early Medieval situation of Finno-Ugric languages, see e.g. South-eastern contact area of Finnic languages in the light of onomastics (Rahkonen 2013).

finno-saamic-mordvin
Yakhr-, -khra, yedr-, -dra and yer-/yar, -er(o), -or(o) names of lakes in Central and North Russia and the possible boundary of the proto-language words *jäkra/ä and *järka/ä. Rahkonen (2011)

Certain innovations shared between Proto-Fennic (identified with the Gulf of Finland) and Proto-Mordvinic (from the Gorodets culture) point to their close contact before the Proto-Fennic expansion, and thus to the identification of Gorodets as Proto-Mordvinic, hence Akozino as Volgaic (Parpola 2018):

  • the noun paradigms and the form and function of individual cases,
  • the geminate *mm (foreign to Proto-Uralic before the development of Fennic under Germanic influence) and other non-Uralic consonant clusters.
  • the change of numeral *luka ‘ten’ with *kümmen.
  • The presence of loanwords of non-Uralic origin, related to farming and trees, potentially Palaeo-European in nature (hence possibly from Siberian influence in north-eastern Europe).
ananino-textile-ware-cultures
Map of archaeological cultures in north-eastern Europe ca. 8th-3rd centuries BC. [The Mid-Volga Akozino group not depicted] Shaded area represents the Ananino cultural-historical society. Purple area show likely zones of predominant Siberian ancestry and N1c-L392 lineages. Blue areas likely zones of predominant CWC ancestry and R1a-Z645 lineages. Fading purple arrows represent likely stepped movements of haplogroup N1c-L392 for centuries (Siberian → Ananino → Akozino → Fennoscandia), found eventually in tarand graves. Blue arrows represent eventual expansions of Fennic and (partially displaced) Saamic. Modified image from Vasilyev (2002).

The introduction of a strongly hierarchical chiefdom system can quickly change the pre-existing social order and lead to a major genetic shift within generations, without a radical change in languages, as shown in Sintashta-Potapovka compared to the preceding Poltavka society (read more about Sintashta).

Fortified settlements in the region represented in part visiting warrior-traders settled through matrimonial relationships with local chiefs, eager to get access to coveted goods and become members of a distribution network that could guarantee them even military assistance. Such a system is also seen synchronously in other cultures of the region, like the Nordic Bronze Age and Lusatian cultures (Parpola 2013).

The most likely situation is that N1c subclades were incorporated from the Circum-Artic region during the Anonino (Permic) expansion to the north, later emerged during the formation of the Akozino group (Volgaic, under Anonino influence), and these subclades in turn infiltrated among the warrior traders that spread all over Fennoscandia and the eastern Baltic (mainly among Fennic, Saamic, Germanic, and Balto-Slavic peoples), during the age of hill forts, creating alliances partially based on exogamy strategies (Parpola 2013).

Over the course of these events, no language change is necessary in any of the cultures involved, since the centre of gravity is on the expanding culture incorporating new lineages:

  • first on the Middle Volga, when Ananino expands to the north, incorporatinig N1c lineages from the Circum-Artic region.
  • then with the expansion of the Akozino-Akhmylovo culture into Ananino territory, admixing with part of its population;
  • then on the Baltic region, when materials are imported from Akozino into Fennoscandia and the eastern Baltic (and vice versa), with local cultures being infiltrated by foreign (Akozino) warrior-traders and their materials;
  • and later with the different population movements that led eventually to a greater or lesser relevance of N1c in modern Finno-Permic populations.

To argue that this infiltration and later expansion of lineages changed the language in one culture in one of these events seems unlikely. To use this argument of “opposite movement of ethnic and language change” for different successive events, and only on selected regions and cultures (and not those where the greatest genetic and cultural impact is seen, like e.g. Sweden for Akozino materials) is illogical.

NOTE. Notice how I write here about “infiltration” and “lineages”, not “migration” or “populations”. To understand that, see below the next section on autosomal studies to compare Bronze Age, Iron Age, Medieval and Modern Estonians, and see how little the population of Estonia (homeland of Proto-Fennic and partially of Proto-Finno-Saamic) has changed since the Corded Ware migrations, suggesting genetic continuity and thus mostly close inter-regional and intra-regional contacts in the Forest Zone, hence a very limited impact of the absorbed N1c lineages (originally at some point incorporated from the Circum-Artic region). You can also check on the most recent assessment of R1a vs. N1c in modern Uralic populations.

Iron Age and later populations

From the session on Estonian samples on ISBA 8, by Tambets et al.:

[Of the 13 samples from the Iron Age tarand-graves] We found that the Iron Age individuals do in fact carry chrY hg N3 (…) Furthermore, based on their autosomal data, all of the studied individuals appear closer to hunter-gatherers and modern Estonians than Estonian CWC individuals do.

EDIT (16 OCT) A recent abstract with Saag as main author (Tambets second) cites 3 out of 5 sampled Iron Age individuals as having haplogroup N3.

EDIT (28 OCT): Notice also the appearance of N1a1a1a1a1a1a1-L1025 in Lithuania (ca. 300 AD), from Damgaard (Nature 2018); the N1c sample of the Krivichi Pskov Long Barrows culture (ca. 8th-10th c. AD), and N1a1a1a1a1a1a7-Y4341 among late Vikings from Sigtuna (ca. 10th-12th c. AD) in Krzewinska (2018).

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

Looking at the plot, the genetic inflow marking the change from the Bronze Age to the Iron Age looks like an obvious expansion of nearby peoples with CWC-related ancestry, i.e. likely from the south-east, near the Middle Volga, where influence of steppe peoples is greater (hence likely Akozino) into a Proto-Fennic population already admixed (since the arrival of Corded Ware groups) with Comb Ware-like populations.

All of these groups were probably R1a-Z645 (likely R1a-Z283) since the expansion of Corded Ware peoples, with an introduction of some N1c lineages precisely during this Iron Age period. This infiltration of N1c-L392 with Akozino is obviously not directly related to Siberian cultures, given what we know about the autosomal description of Estonian samples.

Rather, N1c-L392 lineages were likely part of the incoming (Volgaic) Akozino warrior-traders, who settled among developing chiefdoms based on hill fort settlements of cultures all over the Baltic area, and began to appear thus in some of the new tarand graves associated with the Iron Age in north-eastern Europe.f

A good way to look at this is to realize that no new cluster appears compared to the data we already have from Baltic LN and BA samples from Mittnik et al. (2018), so the Estonian BA and IA clusters must be located (in a proper PCA) in the cline from Pit-Comb Ware culture through Baltic BA to Corded Ware groups:

baltic-samples
PCA and ADMIXTURE analysis reflecting three time periods in Northern European prehistory. a Principal components analysis of 1012 present-day West Eurasians (grey points, modern Baltic populations in dark grey) with 294 projected published ancient and 38 ancient North European samples introduced in this study (marked with a red outline). Population labels of modern West Eurasians are given in Supplementary Fig. 7 and a zoomed-in version of the European Late Neolithic and Bronze Age samples is provided in Supplementary Fig. 8. b Ancestral components in ancient individuals estimated by ADMIXTURE (k = 11)

This genetic continuity from Corded Ware (the most likely Proto-Uralic homeland) to the Proto-Fennic and Proto-Saamic communities in the Gulf of Finland correlates very well with the known conservatism of Finno-Saamic phonology, quite similar to Finno-Ugric, and both to Proto-Uralic (Kallio 2017): The most isolated region after the expansion of Corded Ware peoples, the Gulf of Finland, shielded against migrations for almost 1,500 years, is then the most conservative – until the arrival of Akozino influence.

NOTE. This has its parallel in the phonetic conservatism of Celtic or Italic compared to Finno-Ugric-influenced Germanic, Balto-Slavic, or Indo-Iranian.

Only later would certain regions (like Finland or Lappland) suffer Y-DNA bottlenecks and further admixture events associated with population displacements and expansions, such as the spread of Fennic peoples from their Estonian homeland (evidenced by the earlier separation of South Estonian) to the north and east:

diversification-finnic
The Finnic family tree. Kallio (2014).

The initial Proto-Fennic expansion was probably coupled with the expansion of Proto-Saami to the north, with the Kjelmøy Ware absorbing the Siberian population of Lovozero Ware, and potentially in inner Finland and Karelia with the Luukonsaari and Sirnihta wares (Carpelan and Parpola 2017).

This Proto-Saami population expansion from the mainland to the north, admixing with Lovozero-related peoples, is clearly reflected in the late Iron Age Saamic samples from Levänluhta (ca. 400-800 AD), as a shift (of 2 out of 3 samples) to Siberian-like ancestry from their original CWC_Baltic-like situation (see PCA from Lamnidis et al. 2018 above).

Also, Volgaic and Permic populations from inner Finland and the Forest Zone to the Cis-Urals and Circum-Artic regions probably incorporate Siberian ancestry and N1c-L392 lineages during these and later population movements, while the westernmost populations – Estonian, Mordvinic – remain less admixed (see PCA from Tambets et al. 2018 below).

We also have data of N1c-L392 in Nordic territory in the Middle Ages, proving its likely strong presence in the Mälaren area since the Iron Age, with the arrival of Akozino warrior traders. Similarly, it is found among Balto-Slavic groups along the eastern Baltic area. Obviously, no language change is seen in Nordic Bronze Age and Lusatian territory, and none is expected in Estonian or Finnish territory, either.

Therefore, no “N1c-L392 + Siberian ancestry” can be seen expanding Finno-Ugric dialects, but rather different infiltrations and population movements with limited effects on ancestry and Y-DNA composition, depending on the specific period and region.

estonians-hungarians-mordvinian
Selection of the PCA, with the group of Estonians, Mordovians, and Hungarians selected. See Tambets et al. (2018) for more information.

An issue never resolved

Because N1c-L392 subclades & Siberian ancestry, which appear in different proportions and with different origins among some modern Uralic peoples, do not appear in cultures supposed to host Uralic-speaking populations until the Iron Age, people keep looking into any direction to find the ‘true’ homeland of those ‘Uralic N1c peoples’? Kind of a full circular reasoning, anyone? The same is valid for R1a & steppe ancestry being followed for ‘Indo-Europeans’, or R1b-P312 & Neolithic farmer ancestry being traced for ‘Basques’, because of their distribution in modern populations.

I understand the caution of many pointing to the need to wait and see how samples after 2000 BC are like, in every single period, from the middle and upper Volga, Kama, southern Finland, and the Forest Zone between Fennoscandia and the steppe. It’s like waiting to see how people from Western Yamna and the Carpathian Basin after 3000 BC look like, to fill in what is lacking between East Yamna and Bell Beakers, and then between them and every single Late PIE dialect.

But the answer for Yamna-Bell Beaker-Poltavka peoples during the Late PIE expansion is always going to be “R1b-L23, but with R1a-Z645 nearby” (we already have a pretty good idea about that); and the answer for the Forest Zone and northern Cis- and Trans-Urals area – during the time when Uralic languages are known to have already been spoken there – is always going to be “R1a-Z645, but with haplogroup N nearby”, as is already clear from the data on the eastern Baltic region.

So, without a previously proposed model as to where those amateurs expressing concern about ‘not having enough data’ expect to find those ‘Uralic peoples’, all this waiting for the right data looks more like a waiting for N1c and Siberian ancestry to pop up somewhere in the historic Uralic-speaking area, to be able to say “There! A Uralic-speaking male!”. Not a very reasonable framework to deal with prehistoric peoples and their languages, I should think.

But, for those who want to do that, let me break the news to you already:

ananino-culture-balto-slavic
First N1c – Finno-Ugric person arrives in Estonia to teach Finno-Saamic to Balto-Slavic peoples.

And here it is, an appropriate fantasy description of the ethnolinguistic groups from the region. You are welcome:

  • During the Bronze Age, late Corded Ware groups evolve as the western Textile ware Fennic Balto-Slavic group in the Gulf of Finland; the Netted Ware Saamic Balto-Slavic group of inner Finland; the south Netted Ware / Akozino Volgaic Balto-Slavic groups of the Middle Volga; and the Anonino Permic Balto-Slavic group in the north-eastern Forest Zone; all developing still in close contact with each other, allowing for common traits to permeate dialects.
  • These Balto-Slavic groups would then incorporate west of the Urals during and after the Iron Age (ca. 800-500 BC first, and also later during their expansion to the north) limited ancestry and lineages from eastern European hunter-gatherer groups of Palaeo-European Fennic and Palaeo-Siberian Volgaic and Permic languages from the Circum-Artic region, but they adopted nevertheless the language of the newcomers in every single infiltration of N1c lineages and/or admixture with Siberian ancestry. Oh and don’t forget the Saamic peoples from central Sweden, of course, the famous N1c-L392 ‘Rurikid’ lineages expanding Saamic to the north and replacing Proto-Germanic…

The current model for those obsessed with modern Y-DNA is, therefore, that expanding Neolithic, Bronze Age and Iron Age cultures from north-eastern Europe adopted the languages of certain lineages originally from sub-Neolithic (Scandinavian and Siberian) hunter-gatherer populations of the Circum-Artic region; lineages that these cultures incorporated unevenly during their expansions. Hmmmm… Sounds like an inverse Western movie, where expanding Americans end up speaking Apache, and the eastern coast speaks Spanish until Italian migrants arrive and make everyone speak English… or something. A logic, no-nonsense approach to ethnolinguistic identification.

I kid you not, this is the kind of models we are going to see very soon. In 2018 and 2019, with ancient DNA able to confirm or reject archaeological hypotheses based on linguistic data, people will keep instead creating new pet theories to support preconceived ideas based on the Y-DNA prevalent among modern populations. That is, information available in the 2000s.

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

[Next post on the subject: Corded Ware—Uralic (III): Seima-Turbino and the Ugric and Samoyedic expansion]

Related

Munda admixture happened probably during the ANI-ASI mixture

language-tree-munda

Preprint The genetic legacy of continental scale admixture in Indian Austroasiatic speakers, by Tätte et al. bioRxiv (2018).

Interesting excerpts:

Studies analysing mtDNA and Y chromosome markers have revealed a sex-specific admixture pattern of admixture of Southeast and South Asian ancestry components for Munda speakers. While close to 100% of mtDNA lineages present in Mundas match those in other Indian populations, around 65% of their paternal genetic heritage is more closely related to Southeast Asian than South Asian variation. Such a contrasting distribution of maternal and paternal lineages among the Munda speakers is a classic example of ‘father tongue hypothesis’. However, the temporality of this expansion is contentious. Based on Y-STR data the coalescent time of Indian O2a-M95 haplogroup was estimated to be >10 KYA. Recently, the reconstructed phylogeny of 8.8 Mb region of Y chromosome data showed that Indian O2a-M95 lineages coalesce within a clade nested within East/Southeast Asian within the last ~5-7 KYA. This date estimate sets the upper boundary for the main episode of gene flow of Y chromosomes from Southeast Asia to India.

munda-pca
Supplementary Figure S4. First two components of principal component analysis (PCA). Individuals and population medians (circles) are marked with abbreviations from population names. Different colours represent populations from different geographic areas and/or linguistic groups as shown on the legend on the right. For the full names of populations see Supplementary Table S1. PCA was performed using software EIGENSOFT 6.1.42 on the whole filtered dataset (1072 individuals), previously LD pruned as described in the title of Supplementary Figure S1. The first two principal components describe 5.13% and 2.57% of total variance.

Admixture proportions suggest a novel scenario

Regardless of which West Asian population we used, we found that Munda speakers can be described on average as a mixture of ~19% Southeast Asian, 15% West Asian and 66% Onge (South Asian) components. Alternatively, the West and South Asian components of Munda could be modelled using a single South Asian population (Paniya), accounting on average to 77% of the Munda genome. When rescaling the West and South Asian (Onge) components to 1 to explore the Munda genetic composition prior to the introduction of the Southeast Asian component, we note that the West Asian component is lower (~19%) in Munda compared to Paniya (27%) (Supplementary Table S4: *Average_Lao=0). Consistently with qpGraph analyses in Narasimhan et al. (2018), this may point to an initial admixture of a Southeast Asian substrate with a South Asian substrate free of any West Asian component, followed by the encounter of the resulting admixed population with a Paniya-like population. Such a scenario would imply an inverse relationship between the Southeast and West Asian relative proportions in Munda or, in other words, the increase of Southeast Asian component should cause a greater reduction of the West Asian compared to the reduction in the South Asian component in Munda.

admixture-munda-india
The distribution of genetic components (K=13) based on the global ADMIXTURE analysis (Supplementary Figure S1, S2, S3) for a subset of populations on a map of South and Southeast Asia. The circular legend in the bottom left corner shows the ancestral components corresponding to the colours on pie charts. The sector sizes correspond to population median.

Dating the admixture event

In this study, we have replicated a result previously reported in Chaubey et al. (2011)7 that the Mundas lack one ancestral component (k2) that is characteristic to Indian Indo-European and Dravidian speaking populations. If this component came to India through one of the Indo-Aryan migrations then it would be fair to presume that the Munda admixture happened before this component reached India or at least before it spread all over the country. However, the admixture time computed here, falls in the exact same timeframe as the ANI-ASI mixture has been estimated to have happened in India through which the k2 component probably spread. Therefore, we propose that if the Munda admixture happened at the same time, it is possible for it to have happened in the eastern part of the country, east of Bangladesh, and later when populations from East Asia moved to the area, the Mundas migrated towards central India. Such a scenario, which may be further clarified by ancient DNA analyses, seems to be further supported by the fact that Mundas harbor a smaller fraction of West Asian ancestry compared to contemporary Paniya (Supplementary Table S4) and cannot therefore be seen as a simple admixture product of Southern Indian populations with incoming Southeast Asian ancestries.

namazga-expansion-south-asia
Image from Damgaard et al. (2018). 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

Linguistics and genome-wide data

(…) by and large, the linguistic classification justifies itself but Kharia and Juang do not fit in this simplification perfectly.

Once again, with the current level of detail in genetic studies, there is often no clear dialectal division possible for certain groups without fine-scale population studies, and the help from linguistics and archaeology.

Featured image from open access paper by Chaubey et al. (2011).

Related

Haplogroup R1a and CWC ancestry predominate in Fennic, Ugric, and Samoyedic groups

uralic-languages

Open access Genes reveal traces of common recent demographic history for most of the Uralic-speaking populations, by Tambets et al. Genome Biology (2018).

Interesting excerpts (emphasis mine):

Methods

A total of 286 samples of Uralic-speaking individuals, of those 121 genotyped in this study, were analysed in the context of 1514 Eurasian samples (including 14 samples published for the first time) based on whole genome single nucleotide polymorphisms (SNPs) (Additional file 1: Table S1). All these samples, together with the larger sample set of Uralic speakers, were characterized for mtDNA and chrY markers.

The question as which material cultures may have co-spread together with proto-Uralic and Uralic languages depends on the time estimates of the splits in the Uralic language tree. Deeper age estimates (6,000 BP) of the Uralic language tree suggest a connection between the spread of FU languages from the Volga River basin towards the Baltic Sea either with the expansion of the Neolithic culture of Combed Ware, e.g. [6, 7, 17, 26] or with the Neolithic Volosovo culture [7]. Younger age estimates support a link between the westward dispersion of Proto-Finno-Saamic and eastward dispersion of Proto-Samoyedic with a BA Sejma-Turbino (ST) cultural complex [14, 18, 27, 28] that mediated the diffusion of specific metal tools and weapons from the Altai Mountains over the Urals to Northern Europe or with the Netted Ware culture [23], which succeeded Volosovo culture in the west. It has been suggested that Proto-Uralic may have even served as the lingua franca of the merchants involved in the ST phenomenon [18]. All these scenarios imply that material culture of the Baltic Sea area in Europe was influenced by cultures spreading westward from the periphery of Europe and/or Siberia. Whether these dispersals involved the spread of both languages and people remains so far largely unknown.

The population structure of Uralic speakers

To contextualize the autosomal genetic diversity of Uralic speakers among other Eurasian populations (Additional file 1: Table S1), we first ran the principal component (PC) analysis (Fig. 2a, Additional file 3: Figure S1). The first two PCs (Fig. 2a, Additional file 3: Figure S1A) sketch the geography of the Eurasian populations along the East-West and North-South axes, respectively. The Uralic speakers, along with other populations speaking Slavic and Turkic languages, are scattered along the first PC axis in agreement with their geographic distribution (Figs. 1 and 2a) suggesting that geography is the main predictor of genetic affinity among the groups in the given area. Secondly, in support of this, we find that FST-distances between populations (Additional file 3: Figure S2) decay in correlation with geographical distance (Pearson’s r = 0.77, p < 0.0001). On the UPGMA tree based on these FST-distances (Fig. 2b), the Uralic speakers cluster into several different groups close to their geographic neighbours.

uralic-pca
Principal component analysis (PCA) and genetic distances of Uralic-speaking populations. a PCA (PC1 vs PC2) of the Uralic-speaking populations.

We next used ADMIXTURE [48], which presents the individuals as composed of inferred genetic components in proportions that maximize Hardy-Weinberg and linkage equilibrium in the overall sample (see the ‘Methods’ section for choice of presented K). Overall, and specifically at lower values of K, the genetic makeup of Uralic speakers resembles that of their geographic neighbours. The Saami and (a subset of) the Mansi serve as exceptions to that pattern being more similar to geographically more distant populations (Fig. 3a, Additional file 3: S3). However, starting from K = 9, ADMIXTURE identifies a genetic component (k9, magenta in Fig. 3a, Additional file 3: S3), which is predominantly, although not exclusively, found in Uralic speakers. This component is also well visible on K = 10, which has the best cross-validation index among all tests (Additional file 3: S3B). The spatial distribution of this component (Fig. 3b) shows a frequency peak among Ob-Ugric and Samoyed speakers as well as among neighbouring Kets (Fig. 3a). The proportion of k9 decreases rapidly from West Siberia towards east, south and west, constituting on average 40% of the genetic ancestry of FU speakers in Volga-Ural region (VUR) and 20% in their Turkic-speaking neighbours (Bashkirs, Tatars, Chuvashes; Fig. 3a). The proportion of this component among the Saami in Northern Scandinavia is again similar to that of the VUR FU speakers, which is exceptional in the geographic context. It is also notable that North Russians, sampled from near the White Sea, differ from other Russians by sporting higher proportions of k9 (10–15%), which is similar to the values we observe in their Finnic-speaking neighbours. Notably, Estonians and Hungarians, who are geographically the westernmost Uralic speakers, virtually lack the k9 cluster membership.

siberian-ancestry
Population structure of Uralic-speaking populations inferred from ADMIXTURE analysis on autosomal SNPs in Eurasian context. a Individual ancestry estimates for populations of interest for selected number of assumed ancestral populations (K3, K6, K9, K11). Ancestry components discussed in a main text (k2, k3, k5, k6, k9, k11) are indicated and have the same colours throughout. The names of the Uralic-speaking populations are indicated with blue (Finno-Ugric) or orange (Samoyedic). The full bar plot is presented in Additional file 3: Figure S3. b Frequency map of component k9

We also tested the different demographic histories of female and male lineages by comparing outgroup f3 results for autosomal and X chromosome (chrX) data for pairs of populations (Estonians, Udmurts or Khanty vs others) with high versus low probability to share their patrilineal ancestry in chrY hg N (see the ‘Methods’ section, Additional file 3: Figure S13). We found a minor but significant excess of autosomal affinity relative to chrX for pairs of populations that showed a higher than 10% chance of two randomly sampled males across the two groups sharing their chrY ancestry in hg N3-M178, compared to pairs of populations where such probability is lower than 5% (Additional file 3: Figure S13).

In sum, these results suggest that most of the Uralic speakers may indeed share some level of genetic continuity via k9, which, however, also extends to the geographically close Turkic speakers.

uralic-modern-europe

Identity-by-descent

We found that it is the admixture with the Siberians that makes the Western Uralic speakers different from the tested European populations (Additional file 3: Figure S4A-F, H, J, L). Differentiating between Estonians and Finns, the Siberians share more derived alleles with Finns, while the geographic neighbours of Estonians (and Finns) share more alleles with Estonians (Additional file 3: Figure S4M). Importantly, Estonians do not share more derived alleles with other Finnic, Saami, VUR FU or Ob-Ugric-speaking populations than Latvians (Additional file 3: Figure S4O). The difference between Estonians and Latvians is instead manifested through significantly higher levels of shared drift between Estonians and Siberians on the one hand and Latvians and their immediate geographic neighbours on the other hand. None of the Uralic speakers, including linguistically close Khanty and Mansi, show significantly closer affinities to the Hungarians than any non-FU population from NE Europe (Additional file 3: Figure S4R).

ibd-uralic-genetics
Share of ~ 1–2 cM identity-by-descent (IBD) segments within and between regional groups of Uralic speakers. For each Uralic-speaking population representing lines in this matrix, we performed permutation test to estimate if it shows higher IBD segment sharing with other population (listed in columns) as compared to their geographic control group. Empty rectangles indicate no excess IBD sharing, rectangles filled in blue indicate comparisons when statistically significant excess IBD sharing was detected between one Uralic-speaking population with another Uralic-speaking population (listed in columns), rectangles filled in green mark the comparisons when a Uralic-speaking population shows excess IBD sharing with a non-Uralic-speaking population. For each tested Uralic speaker (matrix rows) populations in the control group that were used to generate permuted samples are indicated using small circles. For example, the rectangle filled in blue for Vepsians and Komis (A) implies that the Uralic-speaking Vepsians share more IBD segments with the Uralic-speaking Komis than the geographic control group for Vepsians, i.e. populations indicated with small circles (Central and North Russians, Swedes, Latvians and Lithuanians). The rectangle filled in green for Vepsians and Dolgans shows that the Uralic-speaking Vepsians share more IBD segments with the non-Uralic-speaking Dolgans than the geographic control group

Time of Siberian admixture

The time depth of the Globetrotter (Fig. 5b) inferred admixture events is relatively recent—500–1900 AD (see also complementary ALDER results, in Additional file 13: Table S12 and Additional file 3: Figure S7)—and agrees broadly with the results reported in Busby et al. [55]. A more detailed examination of the ALDER dates, however, reveals an interesting pattern. The admixture events detected in the Baltic Sea region and VUR Uralic speakers are the oldest (800–900 AD or older) followed by those in VUR Turkic speakers (∼1200–1300 AD), while the admixture dates for most of the Siberian populations (>1500 AD) are the most recent (Additional file 3: Figure S7). The West Eurasian influx into West Siberia seen in modern genomes was thus very recent, while the East Eurasian influx into NE Europe seems to have taken place within the first millennium AD (Fig. 5b, Additional file 3: Figure S7).

Affinities of the Uralic speakers with ancient Eurasians

We next calculated outgroup f3-statistics [48] to estimate the extent of shared genetic drift between modern and ancient Eurasians (Additional file 14: Table S13, Additional file 3: Figures S8-S9). Consistent with previous reports [45, 50], we find that the NE European populations including the Uralic speakers share more drift with any European Mesolithic hunter-gatherer group than Central or Western Europeans (Additional file 3: Figure S9A-C). Contrasting the genetic contribution of western hunter-gatherers (WHG) and eastern hunter-gatherers (EHG), we find that VUR Uralic speakers and the Saami share more drift with EHG. Conversely, WHG shares more drift with the Finnic and West European populations (Additional file 3: Figure S9A). Interestingly, we see a similar pattern of excess of shared drift between VUR and EHG if we substitute WHG with the aDNA sample from the Yamnaya culture (Additional file 3: Figure S9D). As reported before [2, 45], the genetic contribution of European early farmers decreases along an axis from Southern Europe towards the Ural Mountains (Fig. 6, Additional file 3: Figure S9E-F).

yamna-cwc-qpgraph-admixture-uralic
Proportions of ancestral components in studied European and Siberian populations and the tested qpGraph model. a 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. b Admixture proportions (%) of ancestral components. We calculated the admixture proportions summing up the relative shares of a set of intermediate populations to explain the full spectrum of admixture components in the test population. We further did the same for the intermediate node CWC’ and present the proportions of the mixing three components in the stacked column bar of CWC’. Colour codes for ancestral components are as follows: dark green—Western hunter gatherer (WHG’); light green—Eastern hunter gatherer (EHG’); grey—European early farmer (LBK’); dark blue—carriers of Corded Ware culture (CWC’); and dark grey—Siberian. CWC’ consists of three sub-components: blue—Caucasian hunter-gatherer in Yamnaya (CHGinY’); light blue—Eastern hunter-gatherer in Yamnaya (EHGinY’); and light grey—Neolithic Levant (NeolL’)

We then used the qpGraph software [48] to test alternative demographic scenarios by trying to fit the genetic diversity observed in a range of the extant Finno-Ugric populations through a model involving the four basic European ancestral components: WHG, EHG, early farmers (LBK), steppe people of Yamnaya/Corded Ware culture (CWC) and a Siberian component (Fig. 6, Additional file 3: Figure S10). We chose the modern Nganasans to serve as a proxy for the latter component because we see least evidence for Western Eurasian admixture (Additional file 3: Figure S3) among them. We also tested the Khantys for that proxy but the model did not fit (yielding f2-statistics, Z-score > 3). The only Uralic-speaking population that did not fit into the tested model with five ancestral components were Hungarians. The qpGraph estimates of the contributions from the Siberian component show that it is the main ancestry component in the West Siberian Uralic speakers and constitutes up to one third of the genomes of modern VUR and the Saami (Fig. 6). It drops, however, to less than 10% in most of NE Europe, to 5% in Estonians and close to zero in Latvians and Lithuanians.

Discussion

uralic-groups-haplogroup-r1a
Additional file 6: Table S5. Y chromosome haplogroup frequencies in Eurasia. Modified by me: in bold haplogroup N1c and R1a from Uralic-speaking populations, with those in red showing where R1a is the major haplogroup. Observe that all Uralic subgroups – Finno-Permic, Ugric, and Samoyedic – have some populations with a majority of R1a lineages.

One of the notable observations that stands out in the fineSTRUCTURE analysis is that neither Hungarians nor Estonians or Mordovians form genetic clusters with other Uralic speakers but instead do so with a broad spectrum of geographically adjacent samples. Despite the documented history of the migration of Magyars [63] and their linguistic affinity to Khantys and Mansis, who today live east of the Ural Mountains, there is nothing in the present-day gene pool of the sampled Hungarians that we could tie specifically to other Uralic speakers.

Perhaps even more surprisingly, we found that Estonians, who show close affinities in IBD analysis to neighbouring Finnic speakers and Saami, do not share an excess of IBD segments with the VUR or Siberian Uralic speakers. This is eIn this context, it is important to remind that the limited (5%, Fig. 6) East Eurasian impact in the autosomal gene pool of modern Estonians contrasts with the fact that more than 30% of Estonian (but not Hungarian) men carry chrY N3 that has an East Eurasian origin and is very frequent among NE European Uralic speakers [36]. However, the spread of chrY hg N3 is not language group specific as it shows similar frequencies in Baltic-speaking Latvians and Lithuanians, and in North Russians, who in all our analyses are very similar to Finnic-speakers. The latter, however, are believed to have either significantly admixed with their Uralic-speaking neighbours or have undergone a language shift from Uralic to Indo-European [38].ven more striking considering that the immediate neighbours—Finns, Vepsians and Karelians—do.

With some exceptions such as Estonians, Hungarians and Mordovians, both IBD sharing and Globetrotter results suggest that there are detectable inter-regional haplotype sharing ties between Uralic speakers from West Siberia and VUR, and between NE European Uralic speakers and VUR. In other words, there is a fragmented pattern of haplotype sharing between populations but no unifying signal of sharing that unite all the studied Uralic speakers.

Comments

The paper is obviously trying to find a “N1c/Siberian ancestry = Uralic” link, but it shows (as previous papers using ancient DNA) that this identification is impossible, because it is not possible to identify “N1c=Siberian ancestry”, “N1c=Uralic”, or “Siberian ancestry = Uralic”. In fact, the arrival of N subclades and Siberian ancestry are late, both events (probably multiple stepped events) are unrelated to each other, and represent east-west demic diffusion waves (as well as founder effects) that probably coincide in part with the Scythian and Turkic (or associated) expansions, i.e. too late for any model of Proto-Uralic or Proto-Finno-Ugric expansion.

On the other hand, it shows interesting data regarding ancestry of populations that show increased Siberian influence, such as those easternmost groups admixed with Yeniseian-like populations (Samoyedic), those showing strong founder effects (Finnic), or those isolated in the Circum-Artic region with neighbouring Siberian peoples in Kola (Saami). All in all, Hungarians, Estonians and Mordovians seem to show the original situation better than the other groups, which is also reflected in part in Y-DNA, conserved as a majority of R1a lineages precisely in these groups. Just another reminder that CWC-related ancestry is found in every single Uralic group, and that it represents the main ancestral component in all non-Samoyedic groups.

estonians-hungarians-mordvinian
Selection of the PCA, with the group of Estonians, Mordovians, and Hungarians selected.

The qpGraph shows the ancestor of Yamna (likely Khvalynsk) and Corded Ware stemming as different populations from a common (likely Neolithic) node – whose difference is based on the proportion of Anatolian-related ancestry – , that is, probably before the Indo-Hittite expansion; and ends with CWC groups forming the base for all Uralic peoples. Below is a detail of the qpGraph on the left, and my old guess (2017) on the right, for comparison:

yamna-corded-ware-qpgraph

#EDIT (22 sep 2018): I enjoyed re-reading it, and found this particular paragraph funny:

Despite the documented history of the migration of Magyars [63] and their linguistic affinity to Khantys and Mansis, who today live east of the Ural Mountains, there is nothing in the present-day gene pool of the sampled Hungarians that we could tie specifically to other Uralic speakers.

They are so obsessed with finding a link to Siberian ancestry and N1c, and so convinced of Kristiansen’s idea of CWC=Indo-European, that they forgot to examine their own data from a critical point of view, and see the clear link between all Uralic peoples with Corded Ware ancestry and R1a-Z645 subclades… Here is a reminder about Hungarians and R1a-Z282, and about the expansion of R1a-Z645 with Uralic peoples.

Related

Neolithic and Bronze Age Anatolia, Urals, Fennoscandia, Italy, and Hungary (ISBA 8, 20th Sep)

jena-isba8

I will post information on ISBA 8 sesions today as I see them on Twitter (see programme in PDF, and sessions from yesterday).

Official abstracts are listed first (emphasis mine), then reports and images and/or link to tweets. Here is the list for quick access:

Russian colonization in Yakutia

Exploring the genomic impact of colonization in north-eastern Siberia, by Seguin-Orlando et al.

Yakutia is the coldest region in the northern hemisphere, with winter record temperatures below minus 70°C. The ability of Yakut people to adapt both culturally and biologically to extremely cold temperatures has been key to their subsistence. They are believed to descend from an ancestral population, which left its original homeland in the Lake Baykal area following the Mongol expansion between the 13th and 15th centuries AD. They originally developed a semi-nomadic lifestyle, based on horse and cattle breeding, providing transportation, primary clothing material, meat, and milk. The early colonization by Russians in the first half of the 17th century AD, and their further expansion, have massively impacted indigenous populations. It led not only to massive epidemiological outbreaks, but also to an important dietary shift increasingly relying on carbohydrate-rich resources, and a profound lifestyle transition with the gradual conversion from Shamanism to Christianity and the establishment of new marriage customs. Leveraging an exceptional archaeological collection of more than a hundred of bodies excavated by MAFSO (Mission Archéologique Française en Sibérie Orientale) over the last 15 years and naturally kept frozen by the extreme cold temperatures of Yakutia, we have started to characterize the (epi)genome of indigenous individuals who lived from the 16th to the 20th century AD. Current data include the genome sequence of approximately 50 individuals that lived prior to and after Russian contact, at a coverage from 2 to 40 fold. Combined with data from archaeology and physical anthropology, as well as microbial DNA preserved in the specimens, our unique dataset is aimed at assessing the biological consequences of the social and biological changes undergone by the Yakut people following their neolithisation by Russian colons.

NOTE: For another interesting study on Yakutian tribes, see Relationships between clans and genetic kin explain cultural similarities over vast distances.

Ancient DNA from a Medieval trading centre in Northern Finland

Using ancient DNA to identify the ancestry of individuals from a Medieval trading centre in Northern Finland, by Simoes et al.

Analyzing genomic information from archaeological human remains has proved to be a powerful approach to understand human history. For the archaeological site of Ii Hamina, ancient DNA can be used to infer the ancestries of individuals buried there. Situated approximately 30 km from Oulu, in Northern Finland, Ii Hamina was an important trade place since Medieval times. The historical context indicates that the site could have been a melting pot for different cultures and people of diversified genetic backgrounds. Archaeological and osteological evidence from different individuals suggest a rich diversity. For example, stable isotope analyses indicate that freshwater and marine fish was the dominant protein source for this population. However, one individual proved to be an outlier, with a diet containing relatively more terrestrial meat or vegetables. The variety of artefacts that was found associated with several human remains also points to potential differences in religious beliefs or social status. In this study, we aimed to investigate if such variation could be attributed to different genetic ancestries. Ten of the individuals buried in Ii Hamina’s churchyard, dating to between the 15th and 17th century AD, were screened for presence of authentic ancient DNA. We retrieved genome-wide data for six of the individuals and performed downstream analysis. Data authenticity was confirmed by DNA damage patterns and low estimates of mitochondrial contamination. The relatively recent age of these human remains allows for a direct comparison to modern populations. A combination of population genetics methods was undertaken to characterize their genetic structure, and identify potential familiar relationships. We found a high diversity of mitochondrial lineages at the site. In spite of the putatively distant origin of some of the artifacts, most individuals shared a higher affinity to the present-day Finnish or Late Settlement Finnish populations. Interestingly, different methods consistently suggested that the individual with outlier isotopic values had a different genetic origin, being more closely related to reindeer herding Saami. Here we show how data from different sources, such as stable isotopes, can be intersected with ancient DNA in order to get a more comprehensive understanding of the human past.

A closer look at the bottom left corner of the poster (the left columns are probably the new samples):

finland-medieval-admixture

Plant resources processed in HG pottery from the Upper Volga

Multiple criteria for the detection of plant resources processed in hunter-gatherer pottery vessels from the Upper Volga, Russia, by Bondetti et al.

In Northern Eurasia, the Neolithic is marked by the adoption of pottery by hunter-gatherer communities. The degree to which this is related to wider social and lifestyle changes is subject to ongoing debate and the focus of a new research programme. The use and function of early pottery by pre-agricultural societies during the 7th-5th millennia BC is of central interest to this debate. Organic residue analysis provides important information about pottery use. This approach relies on the identification and isotopic characteristics of lipid biomarkers, absorbed into the pores of the ceramic or charred deposits adhering to pottery vessel surfaces, using a combined methodology, namely GC-MS, GC-c-IRMS and EA-IRMS. However, while animal products (e.g., marine, freshwater, ruminant, porcine) have the benefit of being lipid-rich and well-characterised at the molecular and isotopic level, the identification of plant resources still suffers from a lack of specific criteria for identification. In huntergatherer contexts this problem is exacerbated by the wide range of wild, foraged plant resources that may have been potentially exploited. Here we evaluate approaches for the characterisation of terrestrial plant food in pottery through the study of pottery assemblages from Zamostje 2 and Sakhtysh 2a, two hunter-gatherer settlements located in the Upper Volga region of Russia.

GC-MS analysis of the lipids, extracted from the ceramics and charred residues by acidified methanol, suggests that pottery use was primarily oriented towards terrestrial and aquatic animal products. However, while many of the Early Neolithic vessels contain lipids distinctive of freshwater resources, triterpenoids are also present in high abundance suggesting mixing with plant products. When considering the isotopic criteria, we suggest that plants were a major commodity processed in pottery at this time. This is supported by the microscopic identification of Viburnum (Viburnum Opulus L.) berries in the charred deposits on several vessels from Zamostje.

The study of Upper Volga pottery demonstrated the importance of using a multidisciplinary approach to determine the presence of plant resources in vessels. Furthermore, this informs the selection of samples, often subject to freshwater reservoir effects, for 14C dating.

Studies on hunter-gatherer pottery – appearing in eastern Europe before Middle Eastern Neolithic pottery – may be important to understand the arrival of R1a-M17 lineages to the region before ca. 7000 BC. Or not, right now it is not very clear what happened with R1b-P297 and R1a-M17, and with WHG—EHG—ANE ancestry

Bronze Age population dynamics and the rise of dairy pastoralism on the eastern Eurasian steppe

Bronze Age population dynamics and the rise of dairy pastoralism on the eastern Eurasian steppe, by Warinner et al.

Recent paleogenomic studies have shown that migrations of Western steppe herders (WSH), beginning in the Eneolithic (ca. 3300-2700 BCE), profoundly transformed the genes and cultures of Europe and Central Asia. Compared to Europe, the eastern extent of this WSH expansion is not well defined. Here we present genomic and proteomic data from 22 directly dated Bronze Age khirigsuur burials from Khövsgöl, Mongolia (ca. 1380-975 BCE). Only one individual showed evidence of WSH ancestry, despite the presence of WSH populations in the nearby Altai-Sayan region for more than a millennium. At the same time, LCMS/ MS analysis of dental calculus provides direct protein evidence of milk consumption from Western domesticated livestock in 7 of 9 individuals. Our results show that dairy pastoralism was adopted by Bronze Age Mongolians despite minimal genetic exchange with Western steppe herders.

Detail of the images:

mongol-bronze-age-pca

mongol-bronze-age-f4-ancestry

Early Medieval Alemannic graveyard shows diverse cultural and genetic makeup

alemannic-niederstotzingen

Open access Ancient genome-wide analyses infer kinship structure in an Early Medieval Alemannic graveyard, by O’Sullivan et al., Science (2018) 4(9):eaao1262

Interesting excerpts:

Introduction

The Alemanni were a confederation of Germanic tribes that inhabited the eastern Upper Rhine basin and surrounding region (Fig. 1) (1). Roman ethnographers mentioned the Alemanni, but historical records from the 3rd to the 6th century CE contain no regular description of these tribes (2). The upheaval that occurred during the European Migration Period (Völkerwanderung) partly explains the interchangeability of nomenclature with the contemporaneous Suebi people of the same region and periods of geographic discontinuity in the historical record (3). This diverse nomenclature reflects centuries of interactions between Romans and other Germanic groups such as the Franks, Burgundians, Thuringians, Saxons, and Bavarians. With the defeat of the Alemanni by Clovis I of the Franks in 497 CE, Alamannia became a subsumed Duchy of the Merovingian Kingdom. This event solidified the naming of the inhabitants of this region as Alemanni (3). From the 5th to the 8th century CE, integration between the Franks and the Alemanni was reflected by changed burial practices, with households (familia) buried in richly furnished graves (Adelsgrablege) (4). The splendor of these Adelsgräber served to demonstrate the kinship structure, wealth, and status of the familia and also the power of the Franks (Personenverbandstaaten, a system of power based on personal relations rather than fixed territory). Because inclusion in familia during the Merovingian period was not necessarily based on inheritance or provenance, debate continues on the symbolism of these burial rites (5).

The 7th century CE Alemannic burial site at Niederstotzingen in southern Germany, used circa 580 to 630 CE, represents the best-preserved example of such an Alemannic Adelsgrablege. (…)

alemannic-haplogroup

Strontium and oxygen isotope data from the enamel showed that most individuals are local rather than migrants (Table 1, table S2, and fig. S2), except for individuals 10 and 3B. (…)

Analysis of uniparental markers

mtDNA haplogroups were successfully assigned to all 13 individuals (Table 1). Notably, there are three groups of individuals that share, among the assigned positions, identical haplotypes: individuals 4, 9, and 12B in haplogroup X2b4; individuals 1 and 3A in haplogroup K1a; and individuals 2 and 5 in haplogroup K1a1b2a1a.

Most individuals belong to the R1b haplogroup (individuals 1, 3A, 3C, 6, 9, 12A, 12B, and 12C), which has the highest frequency (>70%) in modern western European populations (20). Five individuals (1, 3A, 9, 12B, and 12C) share the same marker (Z319) defining haplogroup R1b1a2a1a1c2b2b1a1 [=ISOGG R1b1a1a2a1a1c2b2b1a1a] (…) individuals 1, 3A, and 6 have R1b lineage and marker Z347 (R1b1a2a1a1c2b2b) [=ISOGG R1b1a1a2a1a1c2b2b], which belongs to the same male ancestral lineage as marker Z319 [i.e. all R1b-U106]. Individual 3B instead carries NRY haplogroup G2a2b1, which is rare in modern north, west, and east European populations (<5%), only reaching common abundance in the Caucasus (>70%), southern Europe, and the Near East (10 to 15%)

Genome-wide capture

alemannic-pca
PCA plot of Niederstotzingen individuals, modern west Eurasians, and selected ancient Europeans. Genome-wide ancient data were projected against modern west Eurasian populations. Colors on PCA indicate more general Eurasian geographic boundaries than countries: dark green, Caucasus; bright green, eastern Europe; yellow, Sardinia and Canary Islands; bright blue, Jewish diaspora; bright purple, western and central Europe; red, southern Europe; dark brown, west Asia; light purple, Spain; dark purple, Russia; pale green, Middle East; orange, North Africa. The transparent circles serve to highlight the genetic overlap between regions of interest.

Genomically, the individuals buried at Niederstotzingen can be split into two groups: Niederstotzingen North (1, 3A, 6, 9, 12B, and 12C), who have genomic signals that most resemble modern northern and eastern European populations, and Niederstotzingen South (3B and 3C), who most resemble modern-day Mediterraneans, albeit with recent common ancestry to other Europeans. Niederstotzingen North is composed of those buried with identifiable artifacts: Lombards (individual 6), Franks (individual 9), and Byzantines (individuals 3A and 12B), all of whom have strontium and oxygen isotope signals that support local provenance (fig. S2) (8). Just two individuals, 3B (Niederstotzingen South) and 10 (no sufficient autosomal data, with R1 Y-haplogroup), have nonlocal strontium isotope signals. The δ18O values suggest that individuals 10 and 3B may have originated from a higher-altitude region, possibly the Swiss-German Alpine foothills (8). Combined with the genome affinity of individual 3B to southern Europeans, these data provide direct evidence for incoming mobility at the site and for contact that went beyond exchange of grave goods (4). Familia had holdings across the Merovingian Kingdom and traveled long distances to maintain them; these holdings could have extended from northern Italy to the North Sea. Nobles displayed and accrued power by recruiting outside individuals into the household as part of their traveling retinue. Extravagant burial rites of these familia are symbolic evidence of the Frankish power systems based on people Personenverbandstaaten imposed from the 5th until the 8th century CE (4). The assignment of grave goods and the burial pattern do not follow any apparent pattern with respect to genetic origin or provenance, suggesting that relatedness and fellowship were held in equal regard at this burial.

Kinship

Both kinship estimates show first-degree relatedness for pairs 1/3A, 1/6, 1/9, 3A/9, and 9/12B and second-degree relatedness for 1/12B, 3A/6, 3A/12B, and 6/9. Except for 12C, all of the Niederstotzingen North individuals are detectably and closely related. The Niederstotzingen South individuals are not detectably related to each other or any other members of the cohort. (…)

We demonstrated that five of the individuals (1, 3A, 6, 9, and 12B) were kin to at least second degree (Fig. 3 and tables S15 and S16); four of these were buried with distinguishable grave goods (discussed above and in fig. S1). These data show that at Niederstotzingen, at least in death, diverse cultural affiliations could be appropriated even within the same family across just two generations. This finding is somewhat similar to the burial of the Frankish King Childeric in the 5th century CE with a combination of Frankish and Byzantine grave goods that symbolized both his provenance and military service to the Romans (4). The burial of three unrelated individuals (3B, 3C, and 12C) in multiple graves beside the rest of the cohort would imply that this Alemannic group buried their dead based on a combination of familial ties and fellowship. One explanation could be that they were adopted as children from another region to be trained as warriors, which was a common practice at the time; these children were raised with equal regard in the familia (2, 4).

alemannic-family
Reconstruction of first- and second-degree relatedness among all related individuals. Bold black lines and blue lines indicate first- and second-degree relatedness, respectively. Dark blue squares are identified males with age-at-death estimates years old (y.o.), mtDNA haplotypes, and NRY haplogroups. Red circles represent unidentified females that passed maternal haplotypes to their offspring. The light square represents one male infant that shares its maternal haplotype with individuals 12B and 9. N.D., not determined.

Conclusion

The 7th century CE burial in Niederstotzingen represents the best-preserved example of an Alemannic Adelsgrablege. The observation that burial of the remains was close to a Roman crossroads, orientated in a considered way, and associated with rich grave goods points to a noble gravesite of an Alemannic familia with external cultural influences. The high percentage of males in the burial site suggests that this site was intended for a ranked warrior group, meaning that the individuals are not representative of the population existing in 7th century CE Alemannia. The kinship estimates show that kinship structure was organized around the familia, which is defined by close association of related and unrelated individuals united for a common purpose. The apparent kinship structure is consistent with the hypothesized Personenverbandstaaten, which was a system by which Merovingian nobles enforced rule in the Duchies of Alemannia, Thuringia, Burgundy, and elsewhere. Beyond the origin of the grave goods, we show isotopic and genetic evidence for contact with communities external to the region and evidence for shared ancestry between northern and southern Europeans. This finding invites debate on the Alemannic power system that may have been highly influenced by mobility and personal relations.

Texts and images distributed under the terms of the Creative Commons Attribution-NonCommercial license.

Related

On the origin of haplogroup R1b-L51 in late Repin / early Yamna settlers

steppe-eneolithic-migrations

A recent comment on the hypothetical Central European origin of PIE helped me remember that, when news appeared that R1b-L51 had been found in Khvalynsk ca. 4250-4000 BC, I began to think about alternative scenarios for the expansion of this haplogroup, with one of them including Central Europe.

Because, if YFull‘s (and Iain McDonald‘s) estimation of the split of R1b-L23 in L51 and Z2103 (ca. 4100 BC, TMRCA ca. 3700 BC) was wrong, by as much as the R1a-Z645 estimates proved wrong, and both subclades were older than expected, then maybe R1b-L51 was not part of the Yamna expansion, but rather part of an earlier expansion with Suvorovo-Novodanilovka into central Europe.

That is, R1b-L51 and R1b-Z2103 would have expanded wih Khvalynsk-Novodanilovka migrants, and they would have either disappeared among local populations, or settled and expanded with successful lineages in certain regions. I think this may give rise to two potential models.

A hidden group in the European east-central steppes?

Here is what Heyd (2011), for example, has to say about the effect of the Khvalynsk-Novodanilovka expansion in the 4th millennium BC, with the first Kurgan wave that shuttered the social, economic, and cultural foundations of south-eastern Europe (before the expansion of west Yamna migrants in the region):

indo-european-anatolian-uralic-migrations
Proto-Anatolian migrations with Khvalynsk-Novodanilovka expansion, including ADMIXTURE data from Wang et al. (2018).

As the Boleraz and Baden tumuli cases in Serbia and Hungary demonstrate, there are earlier, 4th millennium cal. B.C. round tumuli in the Carpathian basin. There are also earlier north-Pontic steppe populations who infiltrated similar environments west of the Black Sea prior to the rise of the Yamnaya culture. This situation can be traced back to the 2nd half of the 5th millennium cal. B.C. to a group of distinct burials, zoomorphic maceheads, long flint blades, triangular flint points, etc., summarized under the term Suvurovo-Novodanilovka (Govedarica 2004; Rassamakin 2004; Anthony 2007; Heyd forthcoming 2011). They also erected round personalized tumuli, though smaller in size and height, above inhumations of single individuals. Suvorovo and Casimcea are the key examples in the lower Danube region of Romania. In northeast Bulgaria, the primary grave of Polska Kosovo (ochre-stained supine extended body position: information communicated by S. Alexandrov) can also be seen as such, as should the Targovishte-“Gonova mogila” primary grave 1 in the Thracian plain with a burial arranged in a supine position with flexed legs, southeast-northwest orientated, and strewed with ochre (Kanchev 1991 , p. 56- 57; Ivanova Gaydarska 2007). In addition to the many copper and shell beads, the 17.4cm long obsidian blade is exceptional, which links this grave to the Csongrád-“Kettoshalom” grave in the south Hungarian plain (Ecsedy 1979). It also yielded an obsidian blade ( 13.2cm long) and copper, shell and limestone beads.

suvorovo-novodanilovka-expansion-europe
The Southeast European distribution of graves of the Suvorovo-Novodanilovka group and such unequipped ones mentioned in the text which can be attributed by burial custom and stratigraphic position in the barrow, plus zoomorphic and abstract animal head sceptres as well as specific maceheads with knobs as from Decea Maresului (mid-5th millennium until around 4000 BC). Heyd (2016).

However, no traces of a tumulus have been recorded above the Kettoshalom tomb. Conventionally, it is dated to the Bodrogkeresztur-period in east Hungary, shortly after 4000 cal. B.C., which would correspond very well with the suggested Cernavodă I (or its less known cultural equivalent in the Thracian plain) attribution for the “Gonova mogila” grave, a cultural background to which the Csongrád grave should have also belonged. Bodrogkeresztur and Cernavodă I periods are not the only examples of 4th millennium cal. B.C. tumuli and burials displaying this steppe connection. Indeed we can find this early steppe impact throughout the 4th millennium cal. B.C. These include adscriptions to the Horodiștea II (Corlateni-Dealul Stadole, grave I: Burtanescu l 998, p. 37; Holbocai, grave 34: Coma 1998, p. 16); to Gordinești-Cernavodă 11 (Liești-Movila Arbănașu, grave 22: Brudiu 2000); to Gorodsk-Usatovo (Corlăteni Dealul Cetăţii, grave I: Comșa 1998, p. 17- 18, in Romania; Durankulak, grave 982: Vajsov 2002, in Bulgaria); and to Cernavodă III(Golyama Detelina, tum. 4: Leshtakov, Borisov 1995), and early (end of 4th millennium cal. B.C.) Ezero in Ovchartsi, primary grave (Kalchev 1994, p. 134-138) and Golyama Detelina, tum. 2 (Kanchev 1991) in Bulgaria. Also the Boleráz and Baden tumuli of Banjevac-Tolisavac and Mokrin in the south Carpathian basin account for this, since one should perhaps take into account primary grave 12 of the Sárrédtudavari-Orhalom tumulus in the Hungarian Alfold: a left-sided crouched juvenile ( 15- 17 y) individual in an oval, NW-SE orientated grave pit 14C dated to 3350-3100 cal. B.C. at 2 sigma (Dani, Ncpper 2006). Neither the burial custom (no ochre strewing or depositing a lump of ochre has been recorded), nor date account for its ascription to the Yamnaya!

All of these tumuli and burials demonstrate, though, that there is already a constant but perhaps low-level 4th millennium cal. B.C. steppe interaction, linking the regions of the north of the Black Sea with those of the west, and reaching deep into the Carpathian basin. This has to be acknowledged. even if these populations remain small, bounded to their steppe habitat with an economy adapted to this special environment, and are not always visible in the record. Indirect hints may help in seeing them, such as the frequent occurrence of horse bones, regarded as deriving from domesticated horses, in Hungarian Baden settlements (Bokonyi 1978; Benecke 1998), and in those of the south German Cham Culture (Matuschik 1999, p. 80-82) and the east German Bernburg Culture (Becker 1999; Benecke 1999). These occur, however, always in low numbers, perhaps not enough to maintain and regenerate a herd. Does this point us towards otherwise archaeologically hidden horsebreeders in the Carpathian basin, before the Yamnaya? In any case, I hope to make one case clear: these are by no means Yamnaya burials in the strict definition! Attribution to the Yamnaya in its strict definition applies.

pit-graves-central-europe
Distribution of Pit-Grave burials west of the Black Sea likely dating to the 2nd half of the 4th millennium BC (triangles: side-crouched burials; filled circles: supine extended burials; open circles: suspected). In Alin Frînculeasa, Bianca Preda, Volker Heyd, Pit-Graves, Yamnaya and Kurgans along the Lower Danube.

Also, about the expansion of Yamna settlers along the steppes:

However, it should have been made clear by the distribution map of the Western Yamnaya that they were confining themselves solely to their own, well-known, steppe habitat and therefore not occupying, or pushing away and expelling, the locally settled farming societies. Also, living solely in the steppes requires another lifestyle, and quite different economic and social bases, most likely very different to the established farming societies. Although surely regarded as incoming strangers, they may therefore not have been seen as direct competitors. This argument can be further enforced when remembering that the lowlands and the steppes in the southeast of Europe had already been populated throughout the 4th millennium cal. B.C., as demonstrated above, by societies with a similar north-Pontic steppe origin and tradition, albeit in lower numbers. It is only for these groups that the Yamnaya may have become a threat, but their common origin and perhaps a similar economic/ social background with comparable lifestyles would surely have assisted to allow rapid assimilation. More important, though, is that farming societies in this region may therefore have been accustomed to dealing and interacting with different people and ethnic strangers for a long time. (…)

When assessing farming and steppe societies’ interaction from a general point of view, attitudes can diverge in three main directions:

  1. the violent one; with raids, fights, struggles, warfare, suppression and finally the superiority and exploitation of the one over the other;
  2. the peaceful one; with a continuous exchange of gifts, goods, work, information and genes in a balanced reciprocal system, leading eventually to the merging of the two societies and creation of a new identity;
  3. the neutral one; with the two societies ignoring each other for a long time.

What we see from trying to understand the record of the Yamnaya, based on their tumuli and burials, and the local and neighbouring contemporary societies, based on their settlements, hoards, and graves, is likely a mixture of all three scenarios, with the balance perhaps more towards exchange in a highly dynamic system with alterations over time. However, violence and raids cannot be ruled out; they would be difficult to see in the archaeological record; or only indirectly, such as the building of hill forts, particularly the defence-like chain of Vucedol hillforts along the south shore of the Danube on the Serbian/Croatian border zone (Tasic 1995a), and the retreat of people into them (Falkenstein 1998, p. 261-262), with other interpretations also possible. And finally, we are dealing here with very different local and neighbouring societies, as well as with more distant contemporary ones, looking, in reality, rather like a chequer board of societies and archaeological cultures (see Parzinger 1993 for the overview). These display different regional backgrounds and traditions leading to different social and settlement organizations, different economic bases and material cultures in the wide areas between Prut and Maritza rivers, and Black Sea and Tisza river. They surely found their individual way of responding to the incoming and settling Yamnaya people.

yamna-tumuli-west-carpathians
Yamnaya tumuli signalling the expansion of West Yamna from ca. 3100 BC (especially after ca. 2950 BC). Heyd (2011).

The best data we have about this potential non-Yamna origin of R1b-L51 – and thus in favour of its admixture in the Carpathian basin – lies in:

  1. The majority of R1a-Z2103 subclades found to date among Yamna samples.
  2. The presence of R1b-Z2103 in the Catacomb culture – in the Northern Caucasus and in Ukraine.
  3. The limited presence of (ancient and modern) R1b-L51 in eastern Europe and India, whose isolated finds are commonly (and simplistically) attributed to ‘late migrations’.
  4. The presence of R1b-L51 (xZ2103) in cultures related to the ‘Yamna package’, but supposedly not to Yamna settlers. So for example I7043, of haplogroup R1b-L151(xU106,xP312), ca. 2500-2200 BC from Szigetszentmiklós-Üdülősor, probably from the Bell Beaker (Csepel group), but maybe from the early Nagýrev culture.
  5. The expansion of its subclades apparently only from a single region, around the Carpathian basin, in contrast to R1b-Z2103.
  6. The already ‘diluted’ steppe admixture found in the earliest samples with respect to Yamna, which points to the appearance after the Yamna admixture with the local population.
  7. Ukrainian archaeologists (in contrast to their Russian colleagues) point to the relevance of North Pontic cultures like Kvitjana and Lower Mikhailovka in the development of Early Yamna in the west, and some eastern European researchers also believe in this similarity.
  8. If R1b-Z2103 and R1b-L51 had expanded with Suvorovo-Novodanilovka migrants to the west, and had admixed later as Hungary_LCA-LBA-like peoples with Yamna migrants during the long-term contacts with other ‘kurganized cultures’ ca. 2900-2500 BC in the Great Hungarian Plains, it could explain some peculiar linguistic traits of North-West Indo-European, and also why R1b-Z2103 appears in cultures associated with this earlier ‘steppe influence’ (i.e. not directly related to Yamna) such as Vučedol (with a R1b-Z2103 sample, see below). That could also explain the presence of R1b-L151(xP312, xU106) in similar Balkan cultures, possibly not directly related to Yamna.
PCA-r1b-l51
Image modified from Wang et al. (2018). PCA of ancient and modern samples. Red circle in dashed line around Varna, Greece Neolithic, and (approximate position of) Smyadovo outliers, part of Khvalynsk-Novodanilovka settlers.

A hidden group among north or west Pontic Eneolithic steppe cultures?

The expansion of Khvalynsk as Novodanilovka into the North Pontic area happened through the south across the steppe, near the coast, with the forest-steppe region working as a clear natural border for this culture of likely horse-riding chieftains, whose economy was probably based on some rudimentary form of mobile pastoralism.

Although archaeologists are divided as to the origin of each individual Middle Eneolithic group near the Black Sea after the end of the Khvalynsk-Novodanilovka period, it seems more or less clear that steppe cultures like Cernavodă, Lower Mikhailovka, or Kvitjana are closer (or “more archaic”) in their steppe features, which connects them to Volga–Ural and Northern Caucasus cultures, like Northern Caucasus, Repin or Khvalynsk.

On the other hand, forest-steppe cultures like Dereivka (including Alexandria) show innovative traits and contacts with para- or sub-Neolithic cultures to the north, like Comb-Pit Ware groups, apart from corded decoration influenced by Trypillian groups to the west, especially in their later (‘Proto-Corded Ware‘) stage after ca. 3500 BC.

If Ukrainian researchers like Rassamakin are right, Early Yamna expanded not only from Repin settlers, but also from local steppe cultures adopting Repin traits to develop an Early Yamna culture, similar to how eastern (Volga–Ural groups) seem to have synchronously adopted Early Yamna without massive affluence of Repin settlements.

Furthermore, local traits develop in southern groups, like anthropomorphic stelae (shared with Kemi-Oba, direct heir of Lower Mikhailovka), and rich burials featuring wagons. These traits are seen in west Yamna settlers.

north-pontic-kvityana-dereivka-repin
Modified from Rassamakin (1999), adding red color to Repin expansion. The system of the latest Eneolithic Pointic cultures and the sites of the Zhivotilovo-Volchanskoe type: 1) Volchanskoe; 2) Zhivotilovka; 3) Vishnevatoe; 4) Koisug.

Problems of this model include:

  1. On the North Pontic area – in contrast to the Volga–Ural region – , there was a clear “colonization” wave of Repin settlers, also supported by Ukrainian researchers, based on the number of new settlements and burials, and on the progressive retreat of Dereivka, Kvitjana, as well as (more recent) Maykop- and Trypillia-related groups from the North Pontic area ca. 3350/3300 BC. It seems unlikely that these expansionist, semi-nomadic, cattle-breeding, patrilineally-related steppe clans that were driving all native populations out of their territories suddenly decided, at some point during their spread into the North Pontic area ca. 3300-3100 BC, to join forces with some foreign male lineages from the area, and then continue their expansion to the west…
  2. Similar to the fate of R1b-P297 subclades in the Baltic after the expansion of Corded Ware migrants, previous haplogropus of the North Pontic region – such as R1a, R1b-V88, and I2 subclades basically disappeared from the ancient DNA record after the expansion of Khvalynsk-Novodanilovka, and then after the expansion of Yamna, as is clear from Yamna, Afanasevo, and Bell Beaker samples obtained to date. This, in combination with what we know about Y-chromosome bottlenecks in post-Neolithic expansions, leaves little space to think that a big enough territorial group with a majority of “native” haplogroups could survive later expansions (be it R1b-L51 or R1a-Z645).
  3. Supporting an expansion of the same male (and partly female) population, the Yamna admixture from east to west is quite homogeneous, with the only difference found in (non-significant) EEF-like proportion which becomes elevated in distant areas [apart from significant ‘southern’ contribution to certain outlier samples]. Based on the also homogeneous Y-DNA picture, the heterogeneity must come, in general, from the female exogamy practiced by expanding groups.
  4. There is a short period, spanning some centuries (approximately 3300-2700 BC), in which the North Pontic area – especially the forest-steppe territories to the west of the Dnieper, i.e. the Upper Dniester, Boh, and Prut-Siret areas – are a chaos of incoming and emigrating, expanding and shrinking groups of different cultures, such as late Trypillian groups, Maykop-related traits, TRB, GAC, (Proto-)Corded Ware, and Early Yamna settlements. No natural geographic frontier can be delimited between these groups, which probably interacted in different ways. Nevertheless, based on their cultural traits, admixture, and especially on their Y-DNA, it seems that they never incorporated foreign male lineages, beyond those they probably had during their initial expansion trends.
  5. The further expansionist waves of Early Yamna seen ca. 3100 BC, from the Danube Delta to the west, give an overall image of continuously expanding patrilineal clans of R1b-M269 subclades since the Khvalynsk-Novodanilovka migration, in different periodic steps, mostly from eastern Pontic-Caspian nuclei, usually overriding all encountered cultures and (especially male) populations, rather than showing long-term collaboration and interaction. Such interaction is seen only in exceptional cases, e.g. the long-term admixture between Abashevo and Poltavka, as seen in Proto-Indo-Iranian peoples and their language.
PCA-Ukraine-r1b-l51
Image modified from Wang et al. (2018). PCA of ancient and modern samples. Arrows depicting Khvalynsk -> Yamna drift (blue), and hypothetic approximate Ukraine Eneolithic -> Yamna drift accompanying R1b-L51 (red).

Consequences

We are living right now an exemplary ego-, (ethno-)nationalism-, and/or supremacy-deflating moment, for some individuals of eastern and northern European descent who believed that R1a or ‘steppe ancestry proportions’ meant something special. The same can be said about those who had interiorized some social or ethnolinguistic meaning for the origin of R1b in western Europe, N1c in north-eastern Europe, as well as Greeks, Iranians, Armenians, or Mediterranean peoples in general of ‘Near Eastern’ ancestry or haplogroups, or peoples of Near Eastern origin and/or language.

These people had linked their haplogroups or ancestry with some fantasy continuity of ‘their’ ancestral populations to ‘their’ territories or languages (or both), and all are being proven wrong.

Apart from teaching such people a lesson about what simplistic views are useful for – whether it is based on ABO or RH group, white skin, blond hair, blue eyes, lactase persistence, or on the own ancestry or Y-DNA haplogroup -, it teaches the rest of us what can happen in the near future among western Europeans. Because, until recently, most western Europeans were comfortably settled thinking that our ancestors were some remnant population from an older, Palaeolithic or Mesolithic population, who acquired Indo-European languages by way of cultural diffusion in different periods, including only minor migrations.

Judging by what we can see now among some individuals of Northern and Eastern European descent, the only thing that can worsen the air of superiority among western Europeans is when they realize (within a few years, when all these stupid battles to control the narrative fade) that not only are they the cultural ‘heirs’ of the Graeco-Roman tradition that began with the Roman Empire, but that most of them are the direct patrilineal descendants of Khvalynsk, Yamna, Bell Beaker, and European Bronze Age peoples, and thus direct descendants of Middle PIE, Late PIE, and NWIE speakers.

steppe-chalcolithic-migrations
Steppe-related migrations ca. 3100-2600 BC with tentative linguistic identification.

The finding of R1b-L51 and R1b-Z2103 among expanding Suvorovo-Novodanilovka chieftains, with pockets of R1b-L51 remaining in steppe-like societies of the Balkans and the Carpathian Basin, would have beautifully complemented what we know about the East Yamna admixture with R1a-Z93 subclades (Uralic speakers) ca. 2600-2100 BC to form Proto-Indo-Iranian, and about the regional admixtures seen in the Balkans, e.g. in Proto-Greeks, with the prevalent J subclades of the region.

It would have meant an end to any modern culture or nation identifying themselves with the ‘true’ Late PIE and Yamna heirs, because these would be exclusively associated with the expansion of R1b-Z2103 subclades with late Repin, and later as the full-fledged Late PIE with Yamna settlers to south-east and central Europe, and to the southern Urals. The language would have had then obviously undergone different language changes in all these territories through long-lasting admixture with other populations. In that sense, it would have ended with the ideas of supremacy in western Europe before they even begin.

The most likely future

However limited the evidence, it seems that R1b-L51 expanded with Yamna, though, based on the estimates for the haplogroups involved, and on marginal hints at the variability of L23 subclades within Yamna and neighbouring populations. If R1b-L51 expanded with West Repin / Early Yamna settlers, this is why they have not yet been found among Yamna samples:

steppe-eneolithic-migrations
Simplified map of Repin expansions from ca. 3500/3400 BC.
  • The subclade division of Yamna settlers needs not be 50:50 for L51:Z2103, either in time or in space. I think this is the simplistic view underlying many thoughts on this matter. Many different expanding patrilineal clans of L23 subclades may have been more or less successful in different areas, and non-Z2103 may have been on the minority, or more isolated relative to Z2103-clans among expanding peoples on the steppe, especially on the east. In fact, we usually talk in terms of “Z2103 vs. L51” as if
    1. these two were the only L23 subclades; and
    2. both had split and succeeded (expanding) synchronously;

    that is, as if there had not been multiple subclades of both haplogroups, and as if there had not been different expansion waves for hundreds of years stemming from different evolving nuclei, involving each time only limited (successful) clans. Many different subclades of haplogroups L23 (xZ2103, xL51), Z2103, and L51 must have been unsuccessful during the ca. 1,500 years of late Khvalynsk and late Repin-Early Yamna expansions in which they must have participated (for approximately 60-75 generations, based on a mean 20-25 years).

  • If we want to imagine a pocket of ‘hidden’ L51 for some region of the North Pontic or Carpathian region, the same can be imagined – and much more likely – for any unsampled territory of expanding late Repin/Early Yamna settlers from the Lower Don – Lower Volga region (probably already a mixed society of L51 and Z2103 subclades since their beginning, as the early Repin culture, ca. 3800 BC), with L51 clans being probably successful to the west.
  • The Repin culture expanded only in small, mobile settlements from the Lower Don – Lower Volga to the north, east, and south, starting ca. 3500/3400 BC, in the waves that eventually gave a rather early distant offshoot in the Altai region, i.e. Afanasevo. Starting ca. 3300 BC in the archaeological record, the majority of R1b-Z2103 subclades found to date in Afanasevo also supports either
    • a mixed Repin society, with Z2103-clans predominating among eastern settlers; or
    • a Repin society marked by haplogroup L51, and thus a cultural diffusion of late Repin/Early Yamna traits among neighbouring (Khvalynsk, Samara, etc.) groups of essentially the same (early Khvalynsk-Novodanilovka) genetic stock in the Volga–Ural region.

    Both options could justify a majority of Z2103 in the Lower Volga–Ural region, with the latter being supported by the scattered archaeological remains of late Repin in the region before the synchronous emergence of Early Yamna findings in the whole Pontic-Caspian steppe.

  • Most Z2103 from Yamna samples to date are from around 3100 BC (in average) onward, and from the right bank of the Lower Don to the east, particularly from the Lower Volga–Ural area (especially the Samara region), which – based on the center of expansion of late Repin settlers – may be depicting an artificially high Z2103-distribution of the whole Yamna community.
repin-expansion-khvalynsk-cultures
Repin expansion into the Volga–Ural region from ca. 3500/3400 BC. Map made by me based on maps and data from Morgunova (2014, 2016). Lopatino is marked with number 64.
  • Yamna sample I0443, R1b-L23 (Y410+, L51-), ca. 3300-2700 BCE from Lopatino II, points to an intermediate subclade between L23 and L51, near one of the supposed late Repin sites (based on kurgan burials with late Repin cultural traits) in the Samara region.
  • Other Balkan cultures potentially unrelated to the Yamna expansion also show Z2103 (and not only L51) subclades, like I3499 (ca. 2884-2666 calBC), of the Vučedol culture, from Beli Manastir-Popova zemlja, which points to the infiltration of Yamna peoples in other cultures. In any case, the appearance of R1b-L23 subclades in the region happens only after the Yamna expansion ca. 3100 BC, probably through intrusions into different neighbouring regions, if these Balkan cultures are not directly derived from Yamna settlements (which is probably the case of the Csepel Bell Beaker or early Nagýrev sample, see above).
  • The diversity of haplogroups found in or around the Carpathian Basin in Late Chalcolithic / Early Bronze Age samples, including L151(xP312, xU106), P312, U106, Z2103, makes it the most likely sink of Yamna settlers, who spread thus with expanding family clans of different R1b-L23 subclades.
  • Even though some Yamna vanguard groups are known to have expanded up to Saxony-Anhalt before ca. 2700 BC, haplogroup Z2103 seems to be restricted to more eastern regions, which suggests that R1b-L51 was already successful among expanding West Yamna clans in Hungary, which gave rise only later to expanding East Bell Beakers (overwhelmingly of L151 subclades). The source of R1b-L51 and L151 expansion over Z2103 must lie therefore in the West Yamna period, and not in the Bell Beaker expansion.
indo-european-uralic-migrations-yamna-gac
Yamna migrants ca. 3300-2600. Most likely site of admixture with GAC circled in red.
  • The R1b-Z2103 found in Poltavka, Catacomb, and to the south point to a late migration displacing the western R1b-L51, only after the late Repin expansion. This is also seen in the steppe ancestry and R1b-Z2103 south of the Caucasus, in Hajji Firuz, which points to this route as a potential source of the supposed “Earliest Proto-Indo-Iranian” (the mariannu term) of the Near East. A similar replacement event happened some centuries later with expanding R1a-Z93 subclades from the east wiping out haplogroup R1b-Z2103 from the Pontic-Caspian steppe.
  • Many ancient samples from Khvalynsk, Northern Caucasus, Yamna, or later ones are reported simply as R1b-M269 or L23, without a clear subclade, so the simplistic ‘Yamna–Z2103’ picture is not real: if one takes into account that Z2103 might have been successful quite early in the eastern region, it is more likely to obtain a successful Y-SNP call of a Z2103 subclade in the Volga–Ural region than a xZ2103 one.
  • There are some modern samples of R1b-L51 in eastern Europe and Asia, whose common simplistic attribution to “late expansions” is usually not substantiated; and also ancient R1b-L51 samples might be confirmed soon for Asia.
  • ‘Western’ features described by archaeologists for West Yamna settlers, associated with Kemi Oba and southern Yamna groups in the North Pontic area – like rich burials with anthropomorphic stelae and wagons – are actually absent in burials from settlers beyond Bulgaria, which does not support their affiliation with these local steppe groups of the Black Sea. Also, a mix with local traditions is seen accross all Early Yamna groups of the Pontic-Caspian steppe, and still genetics and common cultural traits point to their homogeneization under the same patrilineal clans expanding continuously for centuries. The maintenance of local traditions (as evidenced by East Bell Beakers in Iberia related to Iberian Proto-Beakers) is often not a useful argument in genetics, especially when the female population is not replaced.
yamna-settlers-hungary
Yamna settlers in the Great Pannonian Plain, showing only kurgans of Hungary ca. 2950-2500 BC. Yamna Hungary was one of the biggest West Yamna provinces. From Hórvath et al. (2013).

Conclusion

This is what we know, using linguistics, archaeology, and genetics:

  • Middle Proto-Indo-European expanded with Khvalynsk-Novodanilovka after ca. 4800 BC, with the first Suvorovo settlements dated ca. 4600 BC.
  • Archaic Late Proto-Indo-European expanded with late Repin (or Volga–Ural settlers related to Khvalynsk, influenced by the Repin expansion) into Afanasevo ca. 3500/3400 BC.
  • Late Proto-Indo-European expanded with Early Yamna settlers to the west into central Europe and the Balkans ca. 3100 BC; and also to the east (as Pre-Proto-Indo-Iranian) into the southern Urals ca. 2600 BC.
  • North-West Indo-European expanded with Yamna Hungary -> East Bell Beakers, from ca. 2500 BC.
  • Proto-Indo-Iranian expanded with Sintashta, Potapovka, and later Andronovo and Srubna from ca. 2100 BC.

It seems that the subclades from Khvalynsk ca. 4250-4000 BC were wrongly reported – like those of Narasimhan et al. (2018). However, even if they are real and YFull estimates have to be revised, and even if the split had happened before the expansion of Suvorovo-Novodanilovka, the most likely origin of R1b-L51 among Bell Beakers will still be the expansion of late Repin / Early Yamna settlers, and that is what ancient DNA samples will most likely show, whatever the social or political consequences.

The only relevance of the finding of R1b-L51 in one place or another – especially if it is found to be a remnant of a Middle PIE expansion coupled with centuries of admixture and interaction in the Carpathian Basin – is the potential influence of an archaic PIE (or non-IE) layer on the development of North-West Indo-European in Yamna Hungary -> East Bell Beaker. That is, more or less like the Uralic influence related to the appearance of R1a-Z93 among Proto-Indo-Iranians, of R1a-Z284 among Pre-Germanic peoples, and of R1a-Z282 among Balto-Slavic peoples.

I think there is little that ancient DNA samples from West Yamna could add to what we know in general terms of archaeology or linguistics at this point regarding Late PIE migrations, beyond many interesting details. I am sure that those who have not attributed some random 6,000-year-old paternal ancestor any magical (ethnic or nationalist) meaning are just having fun, enjoying more and more the precise data we have now on European prehistoric populations.

As for those who believe in magical consequences of genetic studies, I don’t think there is anything for them to this quest beyond the artificially created grand-daddy issues. And, funnily enough, those who played (and play) the ‘neutrality’ card to feel superior in front of others – the “I only care about the truth”-type of lie, while secretly longing for grandpa’s ethnolinguistic continuity – are suffering the hardest fall.

Related