North-West Indo-Europeans of Iberian Beaker descent and haplogroup R1b-P312

iron-age-early-mediterranean

The recent data on ancient DNA from Iberia published by Olalde et al. (2019) was interesting for many different reasons, but I still have the impression that the authors – and consequently many readers – focused on not-so-relevant information about more recent population movements, or even highlighted the least interesting details related to historical events.

I have already written about the relevance of its findings for the Indo-European question in an initial assessment, then in a more detailed post about its consequences, then about the arrival of Celtic languages with hg. R1b-M167, and later in combination with the latest hydrotoponymic research.

This post is thus a summary of its findings with the help of natural neighbour interpolation maps of the reported Germany_Beaker and France_Beaker ancestry for individual samples. Even though maps are not necessary, visualizing geographically the available data facilitates a direct comprehension of the most relevant information. What I considered key points of the paper are highlighted in bold, and enumerated.

NOTE. To get “more natural” maps, extrapolation for the whole Iberian Peninsula is obtained by interpolation through the use of external data from the British Isles, Central Europe, and Africa. This is obviously not ideal, but – lacking data from the corners of the Iberian Peninsula – this method gives a homogeneous look to all maps. Only data in direct line between labelled samples in each map is truly interpolated for the Iberian Peninsula, while the rest would work e.g. for a wider (and more simplistic) map of European Bronze Age ancestry components.

Chalcolithic

iberia-chalcolithic
Iberian Chalcolithic groups and expansion of the Proto-Beaker package. See full map.

The Proto-Beaker package may or may not have expanded into Central Europe with typical Iberia_Chalcolithic ancestry. A priori, it seems a rather cultural diffusion of traits stemming from west Iberia roughly ca. 2800 BC.

iberia-y-dna-map-chalcolithic
Map of Y-DNA haplogroups among Iberia Chalcolithic samples. See full map.

The situation during the Chalcolithic is only relevant for the Indo-European question insofar as it shows a homogeneous Iberia_Chalcolithic-like ancestry with typical Y-chromosome (and mtDNA) haplogroups of the Iberian Neolithic dominating over the whole Peninsula until about 2500 BC. This might represent an original Basque-Iberian community.

iberia-mtdna-map-chalcolithic
Map of mtDNA haplogroups among Iberia Chalcolithic samples. See full map.

Bell Beaker period

iberia-bell-beaker-period
Iberian Bell Beaker groups and potential routes of expansion. See full map.

The expansion of the Bell Beaker folk brought about a cultural and genetic change in all Europe, to the point where it has been rightfully considered by Mallory (2013) – the last one among many others before him – the vector of expansion of North-West Indo-European languages. Olalde et al. (2019) proved two main points in this regard, which were already hinted in Olalde et al. (2018):

(1) East Bell Beakers brought hg. R1b-L23 and Yamnaya ancestry to Iberia, ergo the Bell Beaker phenomenon was not a (mere) local development in Iberia, but involved the expansion of peoples tracing their ancestry to the Yamnaya culture who eventually replaced a great part of the local population.

iberia-ancestry-bell-beaker-germany_beaker
Natural neighbor interpolation of Germany_Beaker ancestry in Iberia during the Bell Beaker period (ca. 2600-2250 BC). See full map.

(2) Classical Bell Beakers have their closest source population in Germany Beakers, and they reject an origin close to Rhine Beakers (i.e. Beakers from the British Isles, the Netherlands, or northern France), ergo the Single Grave culture was not the origin of the Bell Beaker culture, either (see here).

iberia-y-dna-map-bell-beaker-period
Map of Y-DNA haplogroups among Iberian Bell Beaker samples. See full map.
iberia-mtdna-map-bell-beaker-period
Map of mtDNA haplogroups among Iberian Bell Beaker samples. See full map.

Early Bronze Age

iberia-early-bronze-age
Iberian Early Bronze Age groups and likely population and culture expansions. See full map.

Interestingly, the European Early Bronze Age in Iberia is still a period of adjustments before reaching the final equilibrium. Unlike the situation in the British Isles, where Bell Beakers brought about a swift population replacement, Iberia shows – like the Nordic Late Neolithic period – centuries of genomic balancing between Indo-European- and non-Indo-European-speaking peoples, as could be suggested by hydrotoponymic research alone.

(3) Palaeo-Indo-European-speaking Old Europeans occupied first the whole Iberian Peninsula, before the potential expansion of one or more non-Indo-European-speaking groups, which confirms the known relative chronology of hydrotoponymic layers of Iberia.

iberia-ancestry-early-bronze-age-germany_beaker
Natural neighbor interpolation of Germany_Beaker ancestry in Iberia during the Early Bronze Age period (ca. 2250-1750 BC). See full map.

This balancing is seen in terms of Germany_Beaker vs. Iberia_Chalcolithic ancestry, but also in terms of Y-chromosome haplogroups, with the most interesting late developments happening in southern Iberia, around the territory where El Argar eventually emerged in radical opposition to the Bell Beaker culture.

iberia-y-dna-map-early-bronze-age
Map of Y-DNA haplogroups among Iberia Early Bronze Age samples. See full map.

(4) Bell Beakers and descendants expanded under male-driven migrations, proper of the Indo-European patrilineal tradition, seen in Yamnaya and even earlier in Khvalynsk:

We obtained lower proportions of ancestry related to Germany_Beaker on the X-chromosome than on the autosomes (Table S14), although the Z-score for the differences between the estimates is 2.64, likely due to the large standard error associated to the mixture proportions in the X-chromosome.

germany-beaker-x-chromosome

iberia-mtdna-map-early-bronze-age
Map of mtDNA haplogroups among Iberia Early Bronze Age samples. See full map.

Regarding the PCA, Iberia Bronze Age samples occupy an intermediate cluster between Iberia Chalcolithic and Bell Beakers of steppe ancestry, with Yamnaya-rich samples from the north (Asturias, Burgos) representing the likely source Old European population whose languages survived well into the Roman Iron Age:

iberia-pca-bronze-age
PCA of ancient European samples. Marked and labelled are Bronze Age groups and relevant samples. See full image.

Middle Bronze Age

iberia-middle-bronze-age
Iberian Middle Bronze Age groups and likely population and culture expansions. See full map.

During the Middle Bronze Age, the equilibrium reached earlier is reversed, with a (likely non-Indo-European-speaking) Argaric sphere of influence expanding to the west and north featuring Iberia Chalcolithic and lesser amount of Germany_Beaker ancestry, present now in the whole Peninsula, although in varying degrees.

iberia-ancestry-middle-bronze-age-germany_beaker
Natural neighbor interpolation of Germany_Beaker ancestry in Iberia during the Middle Bronze Age period (ca. 1750-1250 BC). See full map.

All Iberian groups were probably already under a bottleneck of R1b-DF27 lineages, although it is likely that specific subclades differed among regions:

iberia-y-dna-map-middle-bronze-age
Map of Y-DNA haplogroups among Iberia Middle Bronze Age samples. See full map.
iberia-mtdna-map-middle-bronze-age
Map of mtDNA haplogroups among Iberia Middle Bronze Age samples. See full map.

Late Bronze Age

iberia-late-bronze-age
Iberian Late Bronze Age groups and likely population and culture expansions. See full map.

The Late Bronze Age represents the arrival of the Urnfield culture, which probably expanded with Celtic-speaking peoples. A Late Bronze Age transect before their genetic impact still shows a prevalent Germany_Beaker-like Steppe ancestry, probably peaking in north/west Iberia:

iberia-ancestry-late-bronze-age-germany_beaker
Natural neighbor interpolation of Germany_Beaker ancestry in Iberia during the Late Bronze Age period (ca. 1250-750 BC). See full map.

(5) Galaico-Lusitanians were descendants of Iberian Beakers of Germany_Beaker ancestry and hg. R1b-M269. Autosomal data of samples I7688 and I7687, of the Final Bronze (end of the reported 1200-700 BC period for the samples), from Gruta do Medronhal (Arrifana, Coimbra, Portugal) confirms this.

In the 1940s, human bones, metallic artifacts (n=37) and non-human bones were discovered in the natural cave of Medronhal (Arrifana, Coimbra). All these findings are currently housed in the Department of Life Sciences of the University of Coimbra and are analyzed by a multidisciplinary team. The artifacts suggest a date at the beginning of the 1st millennium BC, which is confirmed by radiocarbon date of a human fibula: 890–780 cal BCE (2650±40 BP, Beta–223996). This natural cave has several rooms and corridors with two entrances. No information is available about the context of the human remains. Nowadays these remains are housed mixed and correspond to a minimum number of 11 individuals, 5 adults and 6 non-adults.

In particular, sample I7687 shows hg. R1b-M269, with no available quality SNPs, positive or negative, under it (see full report). They represent thus another strong support of the North-West Indo-European expansion with Bell Beakers.

iberia-y-dna-map-late-bronze-age
Map of Y-DNA haplogroups among Iberian Late Bronze Age samples. See full map.
iberia-mtdna-map-late-bronze-age
Map of mtDNA haplogroups among Iberian Late Bronze Age samples. See full map.

NOTE. To understand how the region around Coimbra was (Proto-)Lusitanian – and not just Old European in general – until the expansion of the Turduli Oppidani, see any recent paper on Bronze Age expansion of warrior stelae, hydrotoponymy, anthroponymy, or theonymy (see e.g. about Spear-vocabulary).

Iron Age

iberia-iron-age-early
Iberian Pre-Roman Iron Age groups and likely population and culture expansions. See full map.

In a complex period of multiple population movements and language replacements, the temporal transect in Olalde et al. (2019) offers nevertheless relevant clues for the Pre-Roman Iron Age:

(6) The expansion of Celtic languages was associated with the spread of France_Beaker-like ancestry, most likely already with the LBA Urnfield culture, since a Tartessian and a Pre-Iberian samples (both dated ca. 700-500 BC) already show this admixture, in regions which some centuries earlier did not show it. Similarly, a BA sample from Álava ca. 910–840 BC doesn’t show it, and later Celtiberian samples from the same area (ca. 4th c. BC and later) show it, depicting a likely north-east to west/south-west routes of expansion of Celts.

iberia-ancestry-iron-age-france_beaker
Natural neighbor interpolation of France_Beaker ancestry in Iberia during the Pre-Roman Iron Age period (ca. 750-250 BC). See full map.

(7) The distribution of Germany_Beaker ancestry peaked, by the Iron Age, among Old Europeans from west Iberia, including Galaico-Lusitanians and probably also Astures and Cantabri, in line with what was expected before genetic research:

iberia-ancestry-iron-age-germany_beaker
Natural neighbor interpolation of Germany_Beaker ancestry in Iberia during the Pre-Roman Iron Age period (ca. 750-250 BC). See full map.

A probably more precise picture of the Final Bronze – Early Iron Age transition is obtained by including the Final Bronze samples I2469 from El Sotillo, Álava (ca. 910-875 BC) as Celtic ancestry buffer to the west, and the sample I3315 from Menorca (ca. 904-861 BC), lacking more recent ones from intermediate regions:

iberia-ancestry-ia-germany_beaker
Natural neighbor interpolation of Germany_Beaker ancestry in Iberia during the Final Bronze Age – Early Iron Age transition. See full map.
iberia-ancestry-ia-france_beaker
Natural neighbor interpolation of France_Beaker ancestry in Iberia during the Final Bronze Age – Early Iron Age transition. See full map.

In terms of Y-DNA and mtDNA haplogroups, the situation is difficult to evaluate without more samples and more reported subclades:

iberia-y-dna-map-iron-age
Map of Y-DNA haplogroups among Iberian Iron Age samples. See full map.
iberia-mtdna-map-iron-age
Map of mtDNA haplogroups among Iberian Iron Age samples. See full map.

In the PCA, Proto-Lusitanian samples occupy an intermediate cluster between Iberian Bronze Age and Bronze Age North (see above), including the Final Bronze sample from Álava, while Celtic-speaking peoples (including Pre-Iberians and Iberians of Celtic descent from north-east Iberia) show a similar position – albeit evidently unrelated – due to their more recent admixture between Iberian Bronze Age and Urnfield/Hallstatt from Central Europe:

iberia-pca-iron-age
PCA of ancient European samples. Marked and labelled are Iron Age groups and relevant samples. See full image.

(8) Iberian-speaking peoples in north-east Iberia represent a recent expansion of the language from the south, possibly accompanied by an increase in Iberia_Chalcolithic/Germany_Beaker admixture from east/south-east Iberia.

(9) Modern Basques represent a recent isolation + Y-DNA bottlenecks after the Roman Iron Age population movements, probably from Aquitanians migrating south of the Pyrenees, admixing with local peoples, and later becoming isolated during the Early Middle Ages and thereafter:

[Modern Basques] overlap genetically with Iron Age populations showing substantial levels of Steppe ancestry.

Assuming that France_Beaker ancestry is associated with the Urnfield culture (spreading with Celtic-speaking peoples), Vasconic speakers were possibly represented by some population – most likely from France – whose ancestry is close to Rhine Beakers (see here).

Alternatively, a Vasconic language could have survived in some France/Iberia_Chalcolithic-like population that got isolated north of the Pyrenees close to the Atlantic Façade during the Bronze Age, and who later admixed with Celtic-speaking peoples south of the Pyrenees, such as the Vascones, to the point where their true ancestry got diluted.

In any case, the clear Celtic Steppe-like admixture of modern Basques supports for the time being their recent arrival to Aquitaine before the proto-historical period, which is in line with hydrotoponymic research.

Conclusion

The most interesting aspects to discuss after the publication of Olalde et al. (2019) would have been thus the nature of controversial Palaeohispanic peoples for which there is not much linguistic data, such as:

  • the Astures and the Cantabri, usually considered Pre-Celtic Indo-European (see here);
  • the Vaccaei, usually considered Celtic;
  • the Vettones, traditionally viewed as sharing the same language as Lusitanians due to their apparent shared hydrotoponymic, anthroponymic, and/or theonymic layers, but today mostly viewed as having undergone Celticization and helped the westward expansion of Celtic languages (and archaeologically clearly divided from Old European hostile neighbours to the west by their characteristic verracos);
  • the Pellendones or the Carpetani, who were once considered Pre-Celtic Indo-Europeans, too;
  • the nature of Tartessian as Indo-European, or maybe even as “Celtic”, as defended by Koch;
  • or the potential remote connection of Basque and Iberian languages in a common trunk featuring Iberian/France_Chalcolithic ancestry (also including Palaeo-Sardo).
pre-roman-palaeohispanic-languages-peoples-iberia-300bc
Pre-Roman Palaeohispanic peoples ca. 300 BC. See full map. Image modified from the version at Wikipedia, a good example of how to disseminate the wrong ideas about Palaeohispanic languages.

Despite these interesting questions still open for discussion, the paper remarked something already known for a long time: that modern Basques had steppe ancestry and Y-DNA proper of the Yamnaya 5,000 years ago, and that Bell Beakers had brought this steppe ancestry and R1b-P312 lineages to Iberia. This common Basque-centric interpretation of Iberian prehistory is the consequence of a 19th-century tradition of obsessively imagining Vasconic-speaking peoples in their medieval territories extrapolated to Cro-Magnons and Atapuerca (no, really), inhabiting undisturbed for millennia a large territory encompassing the whole Iberia and France, “reduced” or “broken” only with the arrival of Celts just before the Roman conquests. A recursive idea of “linguistic autochthony” and “genetic purity” of the peoples of Iberia that has never had any scientific basis.

Similarly, this paper offered the Nth proof already in population genomics that traditional nativist claims for the origin of the Bell Beaker folk in Western Europe were wrong, both southern (nativist Iberian origin) and northern European (nativist Lower Rhine origin). Both options could be easily rejected with phylogeography since 2015, they were then rejected in Olalde et al. and Mathieson et al (2017), then again with the update of many samples in Olalde et al. (2018) and Mathieson et al (2018), and it has most clearly been rejected recently with data from Wang et al. (2018) and its Yamnaya Hungary samples. Findings from Olalde et al. (2019) are just another nail to coffins that should have been well buried by now.

Even David Anthony didn’t have any doubt in his latest model (2017) about the Carpathian Basin origin of North-West Indo-Europeans (see here), and his latest update to the Proto-Indo-European homeland question (2019) shows that he is convinced now about R1b bottlenecks and proper Pre-Yamnaya ancestry stemming from a time well before the Bell Beaker expansion. This won’t be the last setback to supporters of zombie theories: like the hypotheses of an Anatolian, Armenian, or OIT origin of the PIE homeland, other mythical ideas are so entrenched in nationalist and/or nativist tradition that many supporters will no doubt prefer them to die hard, under the most numerous and shameful rejections of endlessly remade reactionary models.

Related

Volga Basin R1b-rich Proto-Indo-Europeans of (Pre-)Yamnaya ancestry

yamnaya-expansion

New paper (behind paywall) by David Anthony, Archaeology, Genetics, and Language in the Steppes: A Comment on Bomhard, complementing in a favourable way Bomhard’s Caucasian substrate hypothesis in the current issue of the JIES.

NOTE. I have tried to access this issue for some days, but it’s just not indexed in my university library online service (ProQuest) yet. This particular paper is on Academia.edu, though, as are Bomhard’s papers on this issue in his site.

Interesting excerpts (emphasis mine):

Along the banks of the lower Volga many excavated hunting-fishing camp sites are dated 6200-4500 BC. They could be the source of CHG ancestry in the steppes. At about 6200 BC, when these camps were first established at Kair Shak III and Varfolomievka (42 and 28 on Figure 2), they hunted primarily saiga antelope around Dzhangar, south of the lower Volga, and almost exclusively onagers in the drier desert-steppes at Kair-Shak, north of the lower Volga. Farther north at the lower/middle Volga ecotone, at sites such as Varfolomievka and Oroshaemoe hunter-fishers who made pottery similar to that at Kair-Shak hunted onagers and saiga antelope in the desert-steppe, horses in the steppe, and aurochs in the riverine forests. Finally, in the Volga steppes north of Saratov and near Samara, hunter-fishers who made a different kind of pottery (Samara type) and hunted wild horses and red deer definitely were EHG. A Samara hunter-gatherer of this era buried at Lebyazhinka IV, dated 5600-5500 BC, was one of the first named examples of the EHG genetic type (Haak et al. 2015). This individual, like others from the same region, had no or very little CHG ancestry. The CHG mating network had not yet reached Samara by 5500 BC.

morgunova-eneolithic-pontic-caspian
Eneolithic settlements (1–5, 7, 10–16, 20, 22–43, 48, 50), burial grounds (6, 8–9, 17–19, 21, 47, 49) and kurgans (44–46) of the steppe Ural-Volga region: 1 Ivanovka; 2 Turganik; 3 Kuzminki; 4 Mullino; 5 Davlekanovo; 6 Sjezheye (burial ground); 7 Vilovatoe; 8 Ivanovka; 9 Krivoluchye; 10–13 LebjazhinkaI-III-IV-V; 14 Gundorovka; 15–16 Bol. Rakovka I-II; 17–18 Khvalunsk I-II; 19 Lipoviy Ovrag; 20 Alekseevka; 21 Khlopkovskiy; 22 Kuznetsovo I; 23 Ozinki II; 24 Altata; 25 Monakhov I; 26 Oroshaemoe; 27 Rezvoe; 28 Varpholomeevka; 29 Vetelki; 30 Pshenichnoe; 31 Kumuska; 32 Inyasovo; 33 Shapkino VI; 34 Russkoe Truevo I; 35 Tsaritsa I-II; 36 Kamenka I; 37 Kurpezhe-Molla; 38 Istay; 39 Isekiy; 40 Koshalak; 41 Kara-Khuduk; 42 Kair-Shak VI; 43 Kombakte; 44 Berezhnovka I-II; 45 Rovnoe; 46 Politotdelskoe; 47 burial near s. Pushkino; 48 Elshanka; 49 Novoorsk; 50 Khutor Repin. Modified from Morgunova (2014).

But before 4500 BC, CHG ancestry appeared among the EHG hunter-fishers in the middle Volga steppes from Samara to Saratov, at the same time that domesticated cattle and sheep-goats appeared. The Reich lab now has whole-genome aDNA data from more than 30 individuals from three Eneolithic cemeteries in the Volga steppes between the cities of Saratov and Samara (Khlopkov Bugor, Khvalynsk, and Ekaterinovka), all dated around the middle of the fifth millennium BC. Many dates from human bone are older, even before 5000 BC, but they are affected by strong reservoir effects, derived from a diet rich in fish, making them appear too old (Shishlina et al 2009), so the dates I use here accord with published and unpublished dates from a few dated animal bones (not fish-eaters) in graves.

Only three individuals from Khvalynsk are published, and they were first published in a report that did not mention the site in the text (Mathieson et al. 2015), so they went largely unnoticed. Nevertheless, they are crucial for understanding the evolution of the Yamnaya mating network in the steppes. They were mentioned briefly in Damgaard et al (2018) but were not graphed. They were re-analyzed and their admixture components were illustrated in a bar graph in Wang et al (2018: figure 2c), but they are not the principal focus of any published study. All of the authors who examined them agreed that these three Khvalynsk individuals, dated about 4500 BC, showed EHG ancestry admixed substantially with CHG, and not a trace of Anatolian Farmer ancestry, so the CHG was a Hotu-Cave or Kotias-Cave type of un-admixed CHG. The proportion of CHG in the Wang et al. (2018) bar graphs is about 20-30% in two individuals, substantially less CHG than in Yamnaya; but the third Khvalynsk individual had more than 50% CHG, like Yamnaya. The ca. 30 additional unpublished individuals from three middle Volga Eneolithic cemeteries, including Khvalynsk, preliminarily show the same admixed EHG/CHG ancestry in varying proportions. Most of the males belonged to Y-chromosome haplogroup R1b1a, like almost all Yamnaya males, but Khvalynsk also had some minority Y-chromosome haplogroups (R1a, Q1a, J, I2a2) that do not appear or appear only rarely (I2a2) in Yamnaya graves.

eneolithic-steppes
Pontic-Caspian steppe and neighbouring groups in the Neolithic. See full map.

Wang et al. (2018) discovered that this middle Volga mating network extended down to the North Caucasian steppes, where at cemeteries such as Progress-2 and Vonyuchka, dated 4300 BC, the same Khvalynsk-type ancestry appeared, an admixture of CHG and EHG with no Anatolian Farmer ancestry, with steppe-derived Y-chromosome haplogroup R1b. These three individuals in the North Caucasus steppes had higher proportions of CHG, overlapping Yamnaya. Without any doubt, a CHG population that was not admixed with Anatolian Farmers mated with EHG populations in the Volga steppes and in the North Caucasus steppes before 4500 BC. We can refer to this admixture as pre-Yamnaya, because it makes the best currently known genetic ancestor for EHG/CHG R1b Yamnaya genomes. The Progress-2 individuals from North Caucasus steppe graves lived not far from the pre-Maikop farmers of the Belaya valley, but they did not exchange mates, according to their DNA.

The hunter-fisher camps that first appeared on the lower Volga around 6200 BC could represent the migration northward of un-admixed CHG hunter-fishers from the steppe parts of the southeastern Caucasus, a speculation that awaits confirmation from aDNA. After 5000 BC domesticated animals appeared in these same sites in the lower Volga, and in new ones, and in grave sacrifices at Khvalynsk and Ekaterinovka. CHG genes and domesticated animals flowed north up the Volga, and EHG genes flowed south into the North Caucasus steppes, and the two components became admixed. After approximately 4500 BC the Khvalynsk archaeological culture united the lower and middle Volga archaeological sites into one variable archaeological culture that kept domesticated sheep, goats, and cattle (and possibly horses). In my estimation, Khvalynsk might represent the oldest phase of PIE.

eneolithic-early-steppes
Pontic-Caspian steppe and neighbouring groups in the Early Eneolithic. See full map.

Anatolian Farmer ancestry and Yamnaya origins

The Eneolithic Volga-North Caucasus mating network (Khvalynsk/Progress-2 type) exhibited EHG/CHG admixtures and Y-chromosome haplogroups similar to Yamnaya, but without Yamnaya’s additional Anatolian Farmer ancestry. (…)

Like the Mesolithic and Neolithic populations here, the Eneolithic populations of Dnieper-Donets II type seem to have limited their mating network to the rich, strategic region they occupied, centered on the Rapids. The absence of CHG shows that they did not mate frequently if at all with the people of the Volga steppes, a surprising but undeniable discovery. Archaeologists have seen connections in ornament types and in some details of funeral ritual between Dnieper-Donets cemeteries of the Mariupol-Nikol’skoe type and cemeteries in the middle Volga steppes such as Khvalynsk and S’yez’zhe (Vasiliev 1981:122-123). Also their cranio-facial types were judged to be similar (Bogdanov and Khokhlov 2012:212). So it it surprising that their aDNA does not indicate any genetic admixture with Khvalynsk or Progress-2. Also, neither they nor the Volga steppe Eneolithic populations showed any Anatolian Farmer ancestry. (…)

All three of the steppe-admixed exceptions were from the Varna region (Mathieson et al. 2018). One of them was the famous “golden man’ at Varna (Krause et al. 2016), Grave 43, whose steppe ancestry was the most doubtful of the three. If he had steppe ancestry, it was sufficiently distant (five+ generations before him) that he was not a statistically significant outlier, but he was displaced in the steppe direction, away from the central values of the majority of typical Anatolian Farmers at Varna and elsewhere. The other two, at Varna (grave 158, a 5-7-year-old girl) and Smyadovo (grave 29, a male 20-25 years old), were statistically significant outliers who had recent steppe ancestry (consistent with grandparents or great-grandparents) of the EHG/CHG Khvalynsk/Progress-2 type, not of the Dnieper Rapids EHG/WHG type.

(…) I believe that the Suvorovo-Cernavoda I movement into the lower Danube valley and the Balkans about 4300 BC separated early PIE-speakers (pre-Anatolian) from the steppe population that stayed behind in the steppes and that later developed into late PIE and Yamnaya.

This archaeological transition marked the breakdown of the mating barrier between steppe and Anatolian Farmer mating networks. After this 4300-4200 BC event, Anatolian Farmer ancestry began to pop up in the steppes. The currently oldest sample with Anatolian Farmer ancestry in the steppes in an individual at Aleksandriya, a Sredni Stog cemetery on the Donets in eastern Ukraine. Sredni Stog has often been discussed as a possible Yamnaya ancestor in Ukraine (Anthony 2007: 239- 254). The single published grave is dated about 4000 BC (4045– 3974 calBC/ 5215±20 BP/ PSUAMS-2832) and shows 20% Anatolian Farmer ancestry and 80% Khvalynsk-type steppe ancestry (CHG&EHG). His Y-chromosome haplogroup was R1a-Z93, similar to the later Sintashta culture and to South Asian Indo-Aryans, and he is the earliest known sample to show the genetic adaptation to lactase persistence (I3910-T). Another pre-Yamnaya grave with Anatolian Farmer ancestry was analyzed from the Dnieper valley at Dereivka, dated 3600-3400 BC (grave 73, 3634–3377 calBC/ 4725±25 BP/ UCIAMS-186349). She also had 20% Anatolian Farmer ancestry, but she showed less CHG than Aleksandriya and more Dereivka-1 ancestry, not surprising for a Dnieper valley sample, but also showing that the old fifth-millennium-type EHG/WHG Dnieper ancestry survived into the fourth millennium BC in the Dnieper valley (Mathieson et al. 2018).

late-eneolithic-repin
Pontic-Caspian steppe and neighbouring groups in the Late Eneolithic. See full map.

Probably, late PIE (Yamnaya) evolved in the same part of the steppes—the Volga-Caucasus steppes between the lower Don, the lower and middle Volga, and the North Caucasus piedmont—where early PIE evolved, and where appropriate EHG/CHG admixtures and Y-chromosome haplogroups were seen already in the Eneolithic (without Anatolian Farmer). There have always been archaeologists who argued for an origin of Yamnaya in the Volga steppes, including Gimbutas (1963), Merpert (1974), and recently Morgunova (2014), who argued that this was where Repin-type ceramics, an important early Yamnaya pottery type, first appeared in dated contexts before Yamnaya, about 3600 BC. The genetic evidence is consistent with Yamnaya EHG/CHG origins in the Volga-Caucasus steppes. Also, if contact with the Maikop culture was a fundamental cause of the innovations in transport and metallurgy that defined the Yamnaya culture, then the lower Don-North Caucasus-lower Volga steppes, closest to the North Caucasus, would be where the earliest phase is expected.

I would still guess that the Darkveti-Meshoko culture and its descendant Maikop culture established the linguistic ancestor of the Northwest Caucasian languages in approximately the region where they remained. I also accept the general consensus that the appearance of the hierarchical Maikop culture about 3600 BC had profound effects on pre-Yamnaya and early Yamnaya steppe cultures. Yamnaya metallurgy borrowed from the Maikop culture two-sided molds, tanged daggers, cast shaft hole axes with a single blade, and arsenical copper. Wheeled vehicles might have entered the steppes through Maikop, revolutionizing steppe economies and making Yamnaya pastoral nomadism possible after 3300 BC.

For those who still hoped that Proto-Indo-Europeans of Yamnaya/Afanasievo ancestry from the Don-Volga region were associated with the expansion of hg. R1a-M417, in a sort of mythical “R1-rich” Indo-European society, it seems this is going to be yet another prediction based on ancestry magic that goes wrong.

Proto-Indo-Europeans were, however, associated with other subclades beyond R1b-M269, probably (as I wrote recently) R1b-V1636, I2a-L699, Q1a-M25, and R1a-YP1272, but also interestingly some J subclade, so let’s see what surprises the new study on Khvalynsk and Yamnaya settlers from the Carpathian Basin brings…

On the bright side, it is indirectly confirmed that late Sredni Stog formed part of the neighbouring Corded Ware-like populations of ca. 20-30%+ Anatolian farmer ancestry that gave Yamnaya its share (ca. 6-10%), relative to the comparatively unmixed Khvalynsk and late Repin population (as shown by Afanasevo).

In this steppe mating network that opened up after the Khvalynsk expansion, the increasing admixture of Anatolian farmer-related ancestry in Yamnaya from east (ca. 2-10%) to west (ca. 6-15%) points to an exogamy of late Repin males in their western/south-western regions with populations around the Don River basin and beyond (and endogamy within the Yamnaya community), in an evolution relevant for language expansions and language contacts during the Late Eneolithic.

NOTE. “Mating network” is my new preferred term for “ancestry”. Also great to see scholars finally talk about “Pre-Yamnaya” ancestry, which – combined with the distinction of Yamnaya from Corded Ware ancestry – will no doubt help differentiate fine-scale population movements of steppe- and forest-steppe-related populations.

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.

The whole issue of the JIES is centered on Caucasian influences on Early PIE as an Indo-Uralic dialect, and this language contact/substrate is useful to locate the most likely candidates for the Northeast and Northwest Caucasian and the Proto-Indo-European homelands.

On the other hand, it would also be interesting to read a discussion of how this Volga homeland of Middle PIE and Don-Volga-Ural homeland of Late PIE would be reconciled with the known continuous contacts of Uralic with Middle and Late PIE (see here) to locate the most likely Proto-Uralic homeland.

Especially because Corded Ware fully replaced all sub-Neolithic groups to the north and east of Khvalynsk/Yamnaya, like Volosovo, so no other population neighbouring Middle and Late Proto-Indo-Europeans survived into the Bronze Age…

EDIT: For those new to this blog, this information on unpublished samples from the Volga River basin is yet another confirmation of Khokhlov’s report on the R1b-L23 samples from Yekaterinovka, and its confirmation by a co-author of The unique elite Khvalynsk male from a Yekaterinovskiy Cape burial, apart from more support to the newest data placing Yekaterinovka culturally and probably chronologically between Samara and Khvalynsk.

Related

“Dinaric I2a” and the expansion of Common Slavs from East-Central Europe

late-iron-age-eastern-europe

A recently published abstract for an upcoming chapter about Early Slavs shows the generalized view among modern researchers that Common Slavs did not spread explosively from the east, an idea proper of 19th-century Romantic views about ancestral tribes of pure peoples showing continuity since time immemorial.

Migrations and language shifts as components of the Slavic spread, by Lindstedt and Salmela, In: Language contact and the early Slavs, Eds. Tomáš Klír, Vít Boček, Universitätsverlag Winter (2019):

The rapid spread of the Proto-Slavic language in the second half of the first millennium CE was long explained by the migration of its speakers out of their small primary habitat in all directions. Starting from the 1980s, alternative theories have been proposed that present language shift as the main scenario of the Slavic spread, emphasizing the presumed role of Slavic as the lingua franca of the Avar Khaganate. Both the migration and the language shift scenarios in their extreme forms suffer from factual and chronological inaccuracy. On the basis of some key facts about human population genetics (the relatively recent common ancestry of the East European populations), palaeoclimatology (the Late Antique Little Ice Age from 536 to around 660 CE), and historical epidemiology (the Justinianic Plague), we propose a scenario that includes a primary rapid demographic spread of the Slavs followed by population mixing and language shifts to and from Slavic in different regions of Europe. There was no single reason for the Slavic spread that would apply to all of the area that became Slavic-speaking. The northern West Slavic area, the East Slavic area, and the Avar sphere and South-Eastern Europe exhibit different kinds of spread: mainly migration to a sparsely populated area in the northwest, migration and language shift in the east, and a more complicated scenario in the southeast. The remarkable homogeneity of Slavic up to the jer shift was not attributable to a lingua-franca function in a great area, as is often surmised. It was a founder effect: Proto-Slavic was originally a small Baltic dialect with little internal variation, and it took time for the individual Slavic languages to develop in different directions.

While I would need to read the whole chapter, in principle it seems easier to agree with this summary than with Curta’s (sort of diffuse) Danubian origin of Common Slavic, based on the likely origin of the Balto-Slavic expansion with the Trzciniec and/or Lusatian culture, close to the Baltic.

A multi-ethnic Chernyakhov culture

In a sneak peek to the expected Järve et al. (2019) paper in review, there are three Chernyakhov samples (ca. calAD 350-550) with different ancestry probably corresponding to the different regions where they stem from (see image below), which supports the idea that Iron Age eastern Europe was a true melting pot where the eventual language of the different cultures depended on many different factors:

chernyakhov-samples-region
Map of the samples from Järve et al. (2019).

From the paper:

The Chernyakhiv culture was likely an ethnically heterogeneous mix based on Goths (Germanic tribes) but also including Sarmatians, Alans, Slavs, late Scythians and Dacians – the entire ancient population of the northern coast of the Black Sea.

Contacts with neighbouring regions were active, and the Chernyakhiv culture is associated with a number of historical events that took place in Europe at that time. In particular, during the Scythian or Gothic wars of the 230s and 270s, barbarians living in the territory of the Chernyakhiv culture (Goths, Ferules, Carps, Bastarns, etc.) carried out regular raids across the Danube Limes of the Roman Empire. However, from the end of the 3rd century the relations of the barbarians with the Roman Empire gained a certain stability. From the reign of Constantine I the Goths, who were part of the Chernyakhiv culture, became federates (military allies) of the Empire.

The Goths also interacted with the inhabitants of the East European forest zone. The Roman historian Jordanes described the military campaigns of the Gothic king Ermanaric against northern peoples (the ancestors of Vends, Slavs, etc., and the inhabitants of the northern Volga region).

NOTE. As it has become traditional in writings about eastern Europe, ‘Slavs’ are assumed – for no particular reason – to be part of the ‘northern peoples of the forest’ since who knows when exactly, and thus appear mentioned in this very text simultaneously as part of Chernyakhov, but also part of peoples to the north of Chernyakhov warring against them…

admixture-chernyakhov
Proportions of Eastern Hunter-Gatherer (EHG, blue), Natufian (red) and Altaian (green) ancestries in Scythian/Sarmatian groups and groups pre- and postdating them inferred using the a) qpAdm and b) ChromoPainter/NNLS method. c–e Correlation of qpAdm and CP/NNLS proportions for the three putative sources evaluated. Steppe populations predating the Scythians: Yamnaya_Ukraine [26], Yamnaya_Kalmykia [15], Ukr_BA (this study). Scythians and Sarmatians: Nomad_IA [15], Scythian_East and Sarmatian_SU [3], Hungarian Scythian, Sarmatian, Central Saka, Tian Shan Saka and Tagar [1], Scy_Ukr, ScySar_SU and Scy_Kaz (this study). Population postdating the Scythians: Chern (this study). See also Table S3.

Genetic variation

(…) the Chernyakhiv samples overlapped with modern Europeans, representing the most ‘western’ range of variation among the groups of this study.

After the end of the Scythian period in the western Eurasian Steppe, the Chernyakhiv culture samples have higher Near Eastern affinity compared to the Scythians preceding them, agreeing with the Gothic component in the multi-ethnic mix of the Chernyakhiv culture.

The higher proportion Near Eastern and (according to CP/NNLS) lower proportion of eastern ancestry in the Chernyakhiv culture samples were mirrored by f4 analyses where Chern showed lower affinity to Han (Z score –3.097) and EHG (Z score –3.643) than Ukrainian Scythian and Bronze Age samples, respectively, as well as higher Near Eastern (Levant_N and Anatolia_N) affinity than Ukrainian Scythians (Z scores 4.696 and 3.933, respectively). It is plausible to assume that this excess Near Eastern ancestry in Chern is related to European populations whose Near Eastern proportion has exceeded that in the steppe populations since the Neolithic expansion of early farmers. While the Chernyakhiv culture was likely ethnically heterogeneous, the three samples in our Chern group appear to represent its Gothic component.

chernyakhov-goths-uralic-clines
PCA obtained by projecting the ancient samples of this study together with published Scythian/Sarmatian and related samples onto a plot based on 537,802 autosomal SNPs in 1,422 modern Eurasians. To improve readability, the modern populations have been plotted as population medians (after outlier removal). Image modified from the paper, including Sredni Stog, Corded Ware/Uralic (with Srubna outliers) and Chernyakhov clusters.Notice the two new Late Yamna and Catacomb samples from Ukraine clustering with other published samples, despite being from the same region as Sredni Stog individuals.

Early Slavs of hg. I2-L621

A post in Anthrogenica shows some subclades of the varied haplogroups that are expected from medieval Poland:

KO_55, Kowalewko (100-300 AD), I1a3a1a1-Y6626
KO_45, Kowalewko (100-300 AD), I2a2a1b2a-L801
KO_22, Kowalewko (100-300 AD), G2a2b-L30
KO_57, Kowalewko (100-300 AD), G2a2b-L30

ME_7, Markowice (1000-1200 AD), I1a2a2a5-Y5384
NA_13, Niemcza, (900-1000 AD), I2a1b2-L621
NA_18, Niemcza, (900-1000 AD), J2a1a-L26

Just because of these samples among Early Slavs, and looking again more carefully at the modern distribution of I2a-L621 subclades, I think now I was wrong in assuming that I2a-L621 in early Hungarian Conquerors would mean they would appear around the Urals as a lineage integrated in Eastern Corded Ware groups. It seems rather a haplogroup with an origin in Central Europe. Whether it was part of a Baltic community that expanded south, or was incorporated during the expansions to the south is unclear. Like hg. E-V13, it doesn’t seem to have been incorporated precisely along the Danube, but closer to the north-east Carpathians.

Especially interesting is the finding of I2a-L621 among Early Slavs from Silesia, a zone of close interaction among early West Slavs. From Curta (2019):

On Common Slavs

In Poland, settlement discontinuity was postulated, to make room for the new, Prague culture introduced gradually from the southeast (from neighboring Ukraine). However, there is increasing evidence of 6th-century settlements in Lower Silesia (western Poland and the lands along the Middle Oder) that have nothing to do with the Prague culture. Nor is it clear how and when did the Prague culture spread over the entire territory of Poland.

On Great Moravia

Svatopluk’s remarkably strong position was immediately recognized by Pope John VIII, who ordered the immediate release of Methodius from his monastic prison in order to place him in 873 under Svatopluk’s protection. One year later (874), Louis the German himself was forced to recognize Svatopluk’s independence through the peace of Forchheim. By that time, the power of Svatopluk had extended into the upper Vistula Basin, over Bohemia, the lands between the Saale and the Elbe rivers, as well as the northern and northeastern parts of the Carpathian Basin.* The Czech prince Bořivoj, a member of the Přemyslid family which would unify and rule Bohemia in the following century, is believed to have been baptized in 874 by Methodius in Moravia together with his wife Ludmila (St. Wenceslas’s grandmother).

*Brather, Archäologie, p. 71. The expansion into the region of the Upper Vistula (Little Poland) results from one of St. Methodius’ prophecies, for which see the Life of Methodius 11, p. 72; Poleski, “Contacts between the Great Moravian empire and the tribes”; Poleski, “Contacts between the tribes in the basins.” Despite an early recognition of the Moravian influences on the material culture in 9th-century southern Poland and Silesia (e.g., Dostál, “Das Vordringen”), the question of Svatopluk’s expansion has triggered in the 1990s a fierce debate among Polish archaeologists. See Wachowski, “Problem”; Abłamowicz, “Górny Śląsk”; Wachowski, “Północny zasięg ekspansji”; Szydłowski, “Czy ślad”; Jaworski, “Elemente.”

On Piast Poland

Mieszko agreed to marry Oda, the daughter of the margrave of the North March, for his first wife had died in 977. The marriage signaled a change in the relations with the Empire, for Mieszko sent troops to help Otto II against the Slavic rebels of 983. He also attacked Bohemia and incorporated Silesia and Lesser Poland into the Piast realm, which prompted Bohemians to ally themselves with the Slavic rebels against whom Emperor Otto was now fighting. By 980, therefore, Mieszko was part of a broader configuration of power, and his political stature was recognized in Scandinavia as well. His daughter, Swietoslawa married first Erik Segersäll of Sweden (ca. 970–ca. 995) and then Sweyn Forkbeard of Denmark (986–1014).26 In the early 990s, together with his wife and children, Mieszko offered his state (called “civitas Schinesghe,” the state of Gniezno) to the pope as a fief, as attested by a unique document known as Dagome iudex and preserved in a late 11th-century summary. The document describes the inner boundaries of the state and peripheral provinces, as if Gniezno were a civitas (city) in Italy, with its surrounding territory. Regional centers, however, did indeed come into being shortly before AD 1000 in Lesser Poland (Cracow, Sandomierz), Pomerania (Gdańsk), and Silesia (Wrocław). Such regional centers came to be distinguished from other strongholds by virtue of the presence within their walls of some of the earliest churches built in stone. Mieszko got his own, probably missionary bishop.

In light of this recent find, which complements the Early Slav of the High Middle Ages from Sunghir (ca. AD 1100-1200), probably from the Vladimir-Suzdalian Rus’, we can assume now less speculatively that I2a-CTS10228 most likely expanded with Common Slavs, because alternative explanations for its emergence in the Carpathian Basin, among Early West Slavs, and among Early East Slavs within this short period of time requires too many unacceptable assumptions.

dinaric-i2a-distribution
Modern distribution of “Dinaric” I2a. Modified from Balanovsky et al. (2008)

Hungarian Conquerors

Knowing that R1a-Z280 was an Eastern Corded Ware lineage, found from Baltic Finns to Finno-Ugric populations of the Trans-Urals, we can probably assign expanding Magyars to at least R1a-Z280, R1a-Z93, and N1c-L392 (xB197) lineages.

From Curta (2019):

Earlier Latin sources, especially those of the first half of the 10th century, refer to Magyars as Huns or Avars. They most likely called themselves Magyars, a word indicating that the language they spoke was not Turkic, but Finno-Ugrian, related to a number of languages spoken in Western Siberia and the southern Ural region. The modern word—Hungarian—derives from the Slavic word for those people, U(n)gri, which is another indication of Ugric roots. This has encouraged the search for the origin of the Hungarian people in the lands to the east from the Ural Mountains, in western Siberia, where the Hungarian language is believed to have emerged between 1000 and 500 BC.

In looking for the Magyar primordial homeland, they draw comparisons with the assemblages found in Hungary that have been dated to the 10th century and attributed to the Magyars. Some of those comparisons had extraordinary results. For example, the excavation of the burial mound cemetery recently discovered near Lake Uelgi, in the Cheliabinsk region of Russia, has produced rosette-shaped harness mounts and silver objects ornamented with palmette and floral designs arranged in reticulated patterns, which are very similar to those of Hungary. But Uelgi is not dated to prehistory, and many finds from that site coincided in time with those found in burial assemblages in Hungary. In other words, although there can be no doubt about the relations between Uelgi and the sites in Hungary attributed to the first generations of Magyars, those relations indicate a migration directly from the Trans-Ural lands, and not gradually, with several other stops in the forest-steppe and steppe zones of Eastern Europe. In the lands west of the Ural Mountains, the Magyars are now associated with the Kushnarenkovo (6th to 8th century) and Karaiakupovo (8th to 10th century) cultures, and with such burial sites as Sterlitamak (near Ufa, Bashkortostan) and Bol’shie Tigany (near Chistopol, Tatarstan).14 However, the same problem with chronology makes it difficult to draw the model of a migration from the lands along the Middle Volga. Many parallels for the so typically Magyar sabretache plates found in Hungary are from that region. They have traditionally been dated to the 9th century, but more recent studies point to the coincidence in time between specimens found in Eastern Europe and those from Hungary.

Adding J2a and I1a samples to the Early Slavic stock, based on medieval samples from Poland – with G2a and E-V13 lineages probably shared with Goths from Wielbark/Chernyakhov, or becoming acculturated in the Carpathian Basin – one is left to wonder which of these lineages actually took part in Common Slavic migrations/acculturation events, whenever and wherever those actually happened.

I have tentatively re-assigned lineages of Hungarian conquerors according to their likely origins in a simplistic way – similar to how the paper classifies them – , now (I think) less speculatively, assuming that Early Slavs likely formed eventually part of them:

hungarian-conquerors-y-dna-slavs
Image modified from the paper, with drawn red square around lineages of likely East Slavic origin, and blue squares around R1a-Z93, R1a-Z283, N1a-Z1936, and N1a-M2004 samples, of likely Ugric origin Y-Hg-s determined from 46 males grouped according to sample age, cemetery and Hg. Hg designations are given according to ISOGG Tree 2019. Grey shading designate distinguished individuals with rich grave goods, color shadings denote geographic origin of Hg-s according to Fig. 1. For samples K3/1 and K3/3 the innermost Hg defining marker U106* was not covered, but had been determined previously.

NOTE. The ancestral origin of lineages is meaningless for an ethnolinguistic identification. The only reasonable assumption is that all the individuals sampled formed part of the Magyar polity, shared Magyar culture, and likely spoke Hungarian, unless there is a clear reason to deny this: which I guess should include at least a clearly ‘foreign’ ancestry (showing a distant cluster compared to the group formed by all other samples), ‘foreign’ isotopic data (showing that he was born and/or raised outside of the Carpathian Basin), and particularly ‘foreign’ cultural assemblage of the burial, if one really wants to risk assuming that the individual didn’t speak Hungarian as his mother tongue.

“Dinaric” or Slavic I2a?

I don’t like the use of “Dinaric I2a”, because it is reminiscent of the use of “Iberian R1b-DF27”, or “Germanic R1b-U106”, when ancient DNA has shown that this terminology is most often wrong, and turns out to be misleading. As misleading as “Slavic R1a”. Recently, a Spanish reader wrote me emails wondering how could I possibly say that R1b-DF27 came from Central Europe, because modern distribution maps (see below) made it evident that the haplogroup expanded from Iberia…

DF27-iberia-france-m167
Contour maps of the derived allele frequencies of the SNPs analyzed in Solé-Morata et al. (2017).

The obvious answer is, these maps show modern distributions, not ancient ones. In the case of R1b-DF27, different Iberian lineages are not even related to the same expansion. At least R1b-M167/SRY2627 lineages seem to have expanded from Central Europe into Iberia much more recently than other DF27 subclades associated with Bell Beakers. What’s more, if R1b-M167/SRY2627 appear densest in north-east Spain it is not because of the impact of Celts or Iberians before the arrival of Romans, but because of the impact of medieval expansions during the Reconquista from northern kingdoms expanding south in the Middle Ages:

iberian-medieval-kingdoms-expansion-population-genomics
Genetic differentiation and the footprints of historical migrations in the Iberian Peninsula. Image modified from Bycroft et al. (2018).

Similarly, the term “Dinaric I2a”, based on the higher density in the Western Balkans, is misleading because it is probably the result of later bottlenecks. Just like the density of different R1a subclades among Modern Slavs is most likely the result of acculturation of different groups, especially to the east and north-east, where language shift is known to have happened in historical times, with the cradle of Russians in particular being a Finno-Volgaic hotspot, later expanding with hg. R1a-Z280 and N1c-L392 lineages.

Now, one may think that maybe Slavs expanded with ALL of these different lineages. Since we are talking about late Iron Age / medieval expansions, there might be confederations of different peoples expanding with a single lingua franca… But no, not really. Not likely in linguistics, not likely in archaeology, and apparently not in population genomics, either.

How many ancient peoples from the Iron Age and Early Middle Ages expanded with so many different lineages? We see bottlenecks in expansions even in recent times: say, in Visigoths under E-V13 (probably recently incorporated during their migrations); in Moors (mostly Berbers) with E-M81 and J; in medieval Iberians under different DF27 bottlenecks during the Reconquista (including huge bottlenecks among Basques); similarly, huge bottlenecks are found in Finnic expansions under N1c…How likely is it that Proto-Slavs (and Common Slavs) expanded with all those attested lineages to date among Early Slavs (E-V13, I2a-L621, R1a-M458, I1, J2a) AND also with other R1a subclades prevalent today, but almost absent in sampled Early Slavs?

To sum up, I am not so sure anymore about the possibility of simplistically assigning R1a-M458 to expanding Common Slavs. R1a-M458 may well have been the prevalent R1a subclade in Central Europe among early Balto-Slavic – and possibly also neighbouring Northern Indo-European-speaking – peoples (let’s see what subclades Tollense and Unetice samples bring), but it is more and more likely that most of the density we see in modern R1a-M458 distribution maps is actually the effect of medieval bottlenecks of West Slavs, similar to the case of Iberia.

r1a-m458-underhill-2015
Modern distribution of R1a-M458, after Underhill et al. (2015).

Related

Złota a GAC-CWC transitional group…but not the origin of Corded Ware peoples

koszyce-gac-zlota-cwc

Open access Unraveling ancestry, kinship, and violence in a Late Neolithic mass grave, by Schroeder et al. PNAS (2019).

Interesting excerpts of the paper and supplementary materials, about the Złota group variant of Globular Amphora (emphasis mine):

A special case is the so-called Złota group, which emerged around 2,900 BCE in the northern part of the Małopolska Upland and existed until 2,600-2,500 BCE. Originally defined as a separate archaeological “culture” (15), this group is mainly defined by the rather local introduction of a distinct form of burial in the area mentioned. Distinct Złota settlements have not yet been identified. Nonetheless, because of the character of its burial practices and material culture, which both retain many elements of the GAC and yet point forward to the Corded Ware tradition, and because of its geographical location, the Złota group has attracted significant archaeological attention (15, 16).

The Złota group buried their dead in a new, distinct type of funerary structure; so-called niche graves (also called catacomb graves). These structures featured an entrance shaft or pit and, below that, a more or less extensive niche, sometimes connected to the entrance area by a narrow corridor. Local limestone was used to seal off the entrance shaft and to pave the floor of the niche, on which the dead were usually placed along with grave goods. This specific and relatively sophisticated form of burial probably reflects contacts between the northern Małopolska Upland and the steppe and forest-steppe communities further to the east, who also buried their dead in a form of catacomb graves. Individual cases of the use of ochre and of deformation of skulls in Złota burials provide further indications of such a connection (15). At the same time, the Złota niche grave practice also retains central elements of the GAC funerary tradition, such as the frequent practice of multiple burials in one grave, often entailing redeposition and violation of the anatomical order of corpses, and thus differs from the catacomb grave customs found on the steppes which are strongly dominated by single graves. Nonetheless, at Złota group cemeteries single burial graves appear, and even in multiple burial graves the identity of each individual is increasingly emphasized, e.g. by careful deposition of the body and through the personal nature of grave goods (16).

globular-amphorae-corded-ware-zlota-amphorae
Correspondence analysis of amphorae from the Złota-graveyards reveals that there is no typological break between Globular Amphorae and Corded Ware Amphorae, including ‘Strichbündelamphorae’ (after Furholt 2008)

Just like its burial practices, the material culture and grave goods of the Złota group combine elements of the GAC, such as amber ornaments and central parts of the ceramic inventory, with elements also found in the Corded Ware tradition, such as copper ornaments, stone shaft-hole axes, bone and shell ornaments, and other stylistic features of the ceramic inventory. In particular, Złota group ceramic styles have been seen as a clear transitional phenomenon between classical GAC styles and the subsequent Corded Ware ceramics, probably playing a key role in the development of the typical cord decoration patterns that came to define the latter (17).

As briefly summarized above, the Złota group displays a distinct funerary tradition and combination of material culture traits, which give the clear impression of a cultural “transitional situation”. While the group also appears to have had long-distance contacts directed elsewhere (e.g. to Baden communities to the south), it is the combination of Globular Amphora traits, on the one hand, and traits found among late Yamnaya or Catacomb Grave groups to the east as well as the closely related Corded Ware groups that emerged around 2,800 BCE, on the other hand, that is such a striking feature of the Złota group and which makes it interesting when attempting to understand cultural and demographic dynamics in Central and Eastern Europe during the early 3rd millennium BCE.

catacomb-grave-ksiaznice
Catacomb grave no. 2a/06 from Książnice, Złota culture (acc. to Wilk 2013). Image from Włodarczak (2017)

Książnice (site 2, grave 3ZC), Świętokrzyskie province. This burial, a so-called niche grave of the Złota type (with a vertical entrance shaft and perpendicularly situated niche), was excavated in 2006 and contained the remains of 8 individuals, osteologically identified as three adult females and five children, positioned on limestone pavement in the niche part of the grave. Radiocarbon dating of the human remains indicates that the grave dates to 2900-2630 BCE, 95.4% probability (Dataset S1). The grave had an oval entrance shaft with a diameter of 60 cm and depth of 130 cm; the depth of the niche reached to 170 cm (both measured from the modern surface), and it also contained a few animal bones, a few flint artefacts and four ceramic vessels typical of the Złota group. Książnice is located in the western part of the Małopolska Upland, which only has a few Złota group sites but a stronger presence of other, contemporary groups (including variants of the Baden culture).

Wilczyce (site 90, grave 10), Świętokrzyskie province. A rescue excavation in 2001 uncovered a niche grave of the Złota type, which had a round entrance shaft measuring 90 cm in diameter. The grave was some 60-65 cm deep below the modern surface and the bottom of the niche was paved with thin limestone plates, on which remains of three individuals had been placed; two adults, one female and one male, and one child. Four ceramic vessels of Złota group type were deposited in the niche along with the bodies. Wilczyce is located in the Sandomierz Upland, an area with substantial presence of both the Globular Amphora culture and Złota group, as well as the Corded Ware culture from 2800 BCE.

zlota-gac-cwc
Genetic affinities of the Koszyce individuals and other GAC groups (here including Złota) analyzed in this study. (A) Principal component analysis of previously published and newly sequenced ancient individuals. Ancient genomes were projected onto modern reference populations, shown in gray. (B) Ancestry proportions based on supervised ADMIXTURE analysis (K = 3), specifying Western hunter-gatherers, Anatolian Neolithic farmers, and early Bronze Age steppe populations as ancestral source populations. LP, Late Paleolithic; M, Mesolithic; EN, Early Neolithic; MN, Middle Neolithic; LN, Late Neolithic; EBA, Early Bronze Age; PWC, Pitted Ware culture; TRB, Trichterbecherkultur/Funnelbeaker culture; LBK, Linearbandkeramik/Linear Pottery culture; GAC, Globular Amphora culture; Złota, Złota culture. Image modified to outline in red GAC and Złota groups.

To further investigate the ancestry of the Globular Amphora individuals, we performed a supervised ADMIXTURE (6) analysis, specifying typical western European hunter-gatherers (Loschbour), early Neolithic Anatolian farmers (Barcın), and early Bronze Age steppe populations (Yamnaya) as ancestral source populations (Fig. 2B). The results indicate that the Globular Amphora/Złota group individuals harbor ca. 30% western hunter-gatherer and 70% Neolithic farmer ancestry, but lack steppe ancestry. To formally test different admixture models and estimate mixture proportions, we then used qpAdm (7) and find that the Polish Globular Amphora/Złota group individuals can be modeled as a mix of western European hunter-gatherer (17%) and Anatolian Neolithic farmer (83%) ancestry (SI Appendix, Table S2), mirroring the results of previous studies.

zlota-steppe-ancestry-cwc
Table S2. qpADM results. The ancestry of most Globular Amphora/Złota group individuals
can be modelled as a two-way mixture of Mesolithic western hunter-gatherers (WHG), and early Anatolian Neolithic farmers (Barcın). The five individuals from Książnice (Złota group) show evidence for additional gene flow, most likely from an eastern source.

The lack of a direct genetic connection of Corded Ware peoples with the Złota group despite their common “steppe-like traits” – shared with Yamna – reveals, once more, how the few “Yamna-like” traits of Corded Ware do not support a direct connection with Indo-Europeans, and are the result of the expansion of the so-called steppe package all over Europe, and particularly among cultures closely related to the Khvalynsk expansion, and later under the influence of expanding Yamna peoples.

The results from Książnice may support that early Corded Ware peoples were in close contact with GAC peoples in Lesser Poland during the complex period of GAC-Trypillia-CWC interactions, and especially close to the Złota group at the beginning of the 3rd millennium BC. Nevertheless, patrilineal clans of Złota apparently correspond to Globular Amphorae populations, with the only male sample available yet being within haplogroup I2a-L801, prevalent in GAC.

NOTE. The ADMIXTURE of Złota samples in common with GAC samples (and in contrast with the shared Sredni Stog – Corded Ware “steppe ancestry”) makes the possibility of R1a-M417 popping up in the Złota group from now on highly unlikely. If it happened, that would complicate further the available picture of unusually diverse patrilineal clans found among Uralic speakers expanding with early Corded Ware groups, in contrast with the strict patrilineal and patrilocal culture of Indo-Europeans as found in Repin, Yamna and Bell Beakers.

Once again the traditional links between groups hypothesized by archaeologists – like Gimbutas and Kristiansen in this case – are wrong, as is the still fashionable trend in descriptive archaeology, of supporting 1) wide cultural relationships in spite of clear-cut inter-cultural differences (and intra-cultural uniformity kept over long distances by genetically-related groups), 2) peaceful interactions among groups based on few common traits, and 3) regional population continuities despite cultural change. These generalized ideas made some propose a steppe language shared between Pontic-Caspian groups, most of which have been proven to be radically different in culture and genetics.

gimbutas-kurgan-indo-european
The background shading indicates the tree migratory waves proposed by Marija Gimbutas, and personally checked by her in 1995. Image from Tassi et al. (2017).

Furthermore, paternal lines show once again marked bottlenecks in expanding Neolithic cultures, supporting their relevance to follow the ethnolinguistic identity of different cultural groups. The steppe- or EHG-related ancestry (if it is in fact from early Corded Ware peoples) in Książnice was thus probably, as in the case of Trypillia, in the form of exogamy with females of neighbouring groups:

The presence of unrelated females and related males in the grave is interesting because it suggests that the community at Koszyce was organized along patrilineal lines of descent, adding to the mounting evidence that this was the dominant form of social organization among Late Neolithic communities in Central Europe. Usually, patrilineal forms of social organization go hand in hand with female exogamy (i.e., the practice of women marrying outside their social group). Indeed, several studies (11, 12) have shown that patrilocal residence patterns and female exogamy prevailed in several parts of Central Europe during the Late Neolithic. (…) the high diversity of mtDNA lineages, combined with the presence of only a single Y chromosome lineage, is certainly consistent with a patrilocal residence system.

funnelbeaker-trypillia-corded-ware
Map of territorial ranges of Funnel Beaker Culture (and its settlement concentrations in Lesser Poland), local Tripolyan groups and Corded Ware Culture settlements (■) at the turn of the 4th/3rd millennia BC.

Since ancient and modern Uralians show predominantly Corded Ware ancestry, and Proto-Uralic must have been in close contact with Proto-Indo-European for a very long time – given the different layers of influence that can be distinguished between them -, it follows as logical consequence that the North Pontic forest-steppes (immediately to the west of the PIE homeland in the Don-Volga-Ural steppes) is the most likely candidate for the expansion of Proto-Uralic, accompanying the spread of Sredni Stog ancestry and a bottleneck under R1a-M417 lineages.

The early TMRCAs in the 4th millennium BC for R1a-M417 and R1a-Z645 support this interpretation, like the R1a-M417 sample found in Sredni Stog. On the other hand, the resurgence of typical GAC-like ancestry in late Corded Ware groups, with GAC lineages showing late TMRCAs in the 3rd millennium BC, proves the disintegration of Corded Ware all over Europe (except in Textile Ceramics- and Abashevo-related groups) as the culture lost its cohesion and different local patrilineal clans used the opportunity to seize power – similar to how eventually I2a-L621 infiltrated eastern (Finno-Ugrian) groups.

Related

Magyar tribes brought R1a-Z645, I2a-L621, and N1a-L392(xB197) lineages to the Carpathian Basin

hungarian-conquerors-turks

The Nightmare Week of “N1c=Uralic” proponents (see here) continues, now with preprint Y-chromosome haplogroups from Hun, Avar and conquering Hungarian period nomadic people of the Carpathian Basin, by Neparaczki et al. bioRxiv (2019).

Abstract:

Hun, Avar and conquering Hungarian nomadic groups arrived into the Carpathian Basin from the Eurasian Steppes and significantly influenced its political and ethnical landscape. In order to shed light on the genetic affinity of above groups we have determined Y chromosomal haplogroups and autosomal loci, from 49 individuals, supposed to represent military leaders. Haplogroups from the Hun-age are consistent with Xiongnu ancestry of European Huns. Most of the Avar-age individuals carry east Eurasian Y haplogroups typical for modern north-eastern Siberian and Buryat populations and their autosomal loci indicate mostly unmixed Asian characteristics. In contrast the conquering Hungarians seem to be a recently assembled population incorporating pure European, Asian and admixed components. Their heterogeneous paternal and maternal lineages indicate similar phylogeographic origin of males and females, derived from Central-Inner Asian and European Pontic Steppe sources. Composition of conquering Hungarian paternal lineages is very similar to that of Baskhirs, supporting historical sources that report identity of the two groups.

Interesting excerpts (emphasis mine):

All N-Hg-s identified in the Avars and Conquerors belonged to N1a1a-M178. We have tested 7 subclades of M178; N1a1a2-B187, N1a1a1a2-B211, N1a1a1a1a3-B197, N1a1a1a1a4-M2118, N1a1a1a1a1a-VL29, N1a1a1a1a2-Z1936 and the N1a1a1a1a2a1c1-L1034 subbranch of Z1936. The European subclades VL29 and Z1936 could be excluded in most cases, while the rest of the subclades are prevalent in Siberia 23 from where this Hg dispersed in a counter-clockwise migratory route to Europe (…). All the 5 other Avar samples belonged to N1a1a1a1a3-B197, which is most prevalent in Chukchi, Buryats, Eskimos, Koryaks and appears among Tuvans and Mongols with lower frequency.

haplogroup-n-pca
First two components of PCA from Hg N1a subbranch distribution in 51 populations including Avars and Conquerors. Colors indicate geographic regions. Three letter codes are given in Supplementary Table S5.

By contrast two Conquerors belonged to N1a1a1a1a4-M2118, the Y lineage of nearly all Yakut males, being also frequent in Evenks, Evens and occurring with lower frequency among Khantys, Mansis and Kazakhs.

Three Conqueror samples belonged to Hg N1a1a1a1a2-Z1936 , the Finno-Permic N1a branch, being most frequent among northeastern European Saami, Finns, Karelians, as well as Komis, Volga Tatars and Bashkirs of the Volga-Ural region.Nevertheless this Hg is also present with lower frequency among Karanogays, Siberian Nenets, Khantys, Mansis, Dolgans, Nganasans, and Siberian Tatars.

The west Eurasian R1a1a1b1a2b-CTS1211 subclade of R1a is most frequent in Eastern Europe especially among Slavic people. This Hg was detected just in the Conqueror group (K2/18, K2/41 and K1/10). Though CTS1211 was not covered in K2/36 but it may also belong to this sub-branch of Z283.

Hg I2a1a2b-L621 was present in 5 Conqueror samples, and a 6th sample form Magyarhomorog (MH/9) most likely also belongs here, as MH/9 is a likely kin of MH/16 (see below). This Hg of European origin is most prominent in the Balkans and Eastern Europe, especially among Slavic speaking groups. It might have been a major lineage of the Cucuteni-Trypillian culture and it was present in the Baden culture of the Chalcolithic Carpathian Basin.

hungarian-conquerors-y-dna
Image modified from the paper, with drawn red square around lineages of likely Ugric origin, and squares around R1a-Z93, R1a-Z283, N1a-Z1936, and N1a-M2004 samples. Y-Hg-s determined from 46 males grouped according to sample age, cemetery and Hg. Hg designations are given according to ISOGG Tree 2019. Grey shading designate distinguished individuals with rich grave goods, color shadings denote geographic origin of Hg-s according to Fig. 1. For samples K3/1 and K3/3 the innermost Hg defining marker U106* was not covered, but had been determined previously.

We identified potential relatives within Conqueror cemeteries but not between them. The uniform paternal lineages of the small Karos3 (19 graves) and Magyarhomorog (17 graves) cemeteries approve patrilinear organization of these communities. The identical I2a1a2b Hg-s of Magyarhomorog individuals appears to be frequent among high-ranking Conquerors, as the most distinguished graves in the Karos2 and 3 cemeteries also belong to this lineage. The Karos2 and Karos3 leaders were brothers with identical mitogenomes 11 and Y-chromosomal STR profiles (Fóthi unpublished). The Sárrétudvari commoner cemetery seems distinct from the others, containing other sorts of European Hg-s. Available Y-chromosomal and mtDNA data from this cemetery suggest that common people of the 10th century rather represented resident population than newcomers. The great diversity of Y Hg-s, mtDNA Hg-s, phenotypes and predicted biogeographic classifications of the Conquerors indicate that they were relatively recently associated from very diverse populations.

Surprising about the Hungarian conquerors – although in line with the historical accounts – is the varied patrilineal origin of clans, including Q1a, G2a2b, I1, E1b1b, R1b, J1, or J2 – some of which (depending on specific lineages) may have appeared earlier in the Carpathian Basin or south-eastern Europe.

However, out of the 27 conqueror elite samples, 17 are of haplogroups most likely related to Ugric populations beyond the Urals: R1a-Z645, I2-L621, and two specific N1a-L392 lineages (see below). In fact, there are three high-ranking conqueror elites of hg. I2-L621 (one of them termed a “leader”, brother to an unpublished leader of Karos3, and all of them possibly family), one of hg. R1a-Z280, one of hg. R1a-Z93 (which should be added to the Árpáds), and one of hg. N1a-Z1936, which gives a good idea of the ruling class among the elite Ugric settlers.

NOTE. The Q1a sample is also likely to be found in the mixed population of the West Siberian forest-steppes, since it was found in Mesolithic-Neolithic samples from eastern Europe to Lake Baikal, and in Bronze Age Siberian groups, although admittedly it may have formed part of an Avar Transtisza group, or even earlier Hunnic or Scythian groups along the steppes. Without precise subclades it’s impossible to know.

arrival-of-hungarians-arpad
The seven chieftains of the Hungarians, detail of Arrival of the Hungarians, from Árpád Feszty’s and his assistants’ vast (1800 m2) cyclorama, painted to celebrate the 1000th anniversary of the Magyar conquest of Hungary, now displayed at the Ópusztaszer National Heritage Park in Hungary. Image from Wikipedia.

I2a-L621

I2a-L621 (xS17250) or I2a1b2 in the old nomenclature, is found in 6 early conquerors (including one leader), on a par with R1a and N samples. This haplogroup is found widely distributed in ancient samples, due to its early split (formed ca. 9200 BC, TMRCA ca. 4500 BC) and expansion, probably with Neolithic populations. I can’t seem to find samples of this early haplogroup from the Carpathian Basin, as mentioned in the text, although it wouldn’t be strange, because it appears also in Neolithic Iberia, and in modern populations from western Europe.

Nevertheless, I2a-L621 samples seem to be concentrated mainly in Mesolithic-Neolithic cultures of Fennoscandia, and appeared also in Sikora et al. (2017) in a sample of the High Middle Ages from Sunghir (ca. AD 1100-1200), probably from the Vladimir-Suzdalian Rus’, in a region where clearly tribes of Volga Finns were being assimilated at the time. The reported SNP call by Genetiker is A16681 (see Yfull), deep within I2a-CTS10228. It is possibly also behind a modern Saami from Chalmny Varre (ca. AD 1800) of hg. I2a in Lamnidis et al. (2018).

Lacking precise subclades from Hungarian conquerors this is pure speculation, but modern samples may also point to I2a-CTS10228 (formed ca. 3100 BC, TMRCA ca. 1800 BC) as a Finno-Ugric lineage in common with R1a, which must have expanded to the Urals and beyond with eastern Corded Ware groups or (more likely) succeeding cultures. This is in line with the association of certain I2a lineages with modern Uralic peoples or populations from their historical regions in eastern Europe, and linked thus to the most likely homeland of Uralians in the eastern European forests:

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, and also of I lineages. Data from Tambets et al. (2018).

R1a-Z645

Regarding the important question of the ethnic makeup of Ugric populations stemming from the Urals, the most interesting (and expected) data is the presence of R1a-Z645 lineages among high-ranking conquerors, in particular four R1a-Z280 subclades proper of Finno-Ugrians.

This proves that, in line with the old split and expansion of R1a-CTS1211 (formed ca. 2600 BC, TMRCA ca. 2400 BC), and its finding in Bronze Age Fennoscandian samples, only some late R1a-Z280 (xZ92) lineages (see Z280 on YFull) may show a clear identification with early acculturated Uralic speakers, with the main early acculturated Balto-Slavic R1a haplogroup remaining R1a-M458.

I recently hypothesized this late connection of Slavs with very specific R1a-Z280 (xZ92) lineages based on analyses of modern populations (like Slovenians), because the connection of ancient Finno-Ugrians with modern Z92 samples was already evident:

(…) subclades of hg. R1a1a1b1a2-Z280 (xR1a1a1b1a2a-Z92) seem to have also been involved in early Slavic expansions, like R1a1a1b1a2b3a-CTS3402 (formed ca. 2200 BC, TMRCA ca. 2200 BC), found among modern West, South, and East Slavic populations and in Fennoscandia, prevalent e.g. among modern Slovenians which points to a northern origin of its expansion (Maisano Delser et al. 2018).

This finding also supports the expected shared R1a-Z280 lineages among ancient Finno-Ugric populations, as predicted from the study of modern Permic and Ugric peoples in Dudás et al. (2019).

r1a-z282-z280-z2125-distribution
Modified image, from Underhill et al. (2015). Spatial frequency distributions of Z282 (green) and Z93 (blue) affiliated haplogroups. Notice the distribution of R1a-Z280 (xZ92), i.e. R1a-M558, compared to the ancient Finno-Ugric distribution.

Furthermore, while we don’t have precise R1a-Z93 lineages to compare with the new Hunnic sample reported, we already know that some archaic R1a-Z2124 subclades stem from the forest-steppe areas of the Cis- and Trans-Urals, and the two newly reported R1a-Z93 Hungarian conqueror elites, like those of the Árpád dynasty, probably belong to them.

There is an obvious lack of continuity in specific paternal lineages among the Hunnic, the Avar, and the Conqueror periods, which makes any simplistic identification of all R1a-Z93 lineages as stemming from Avars, Huns, or the Iron Age Pontic-Caspian steppes clearly flawed. Comparing R1a-Z93 in Hungarian Conquerors with Huns is like comparing them with samples of the Srubna or earlier periods… Similarly, comparing the Hunnic R1b-U106 or the early Avar I1 to later Hungarian samples is not warranted without precise subclades, because they most likely correspond to different Germanic populations: Goths among Huns, then Longobards, then likely peoples descended from Franks and Irish Monks (the latter with R1b-P312).

N1a-L392

Second behind R1a subclades are, as expected, N1a-L392 (N1c in the old nomenclature).

Avars are dominated by a specific N1a-L392 subclade, N1a-B197, as we recently discovered in Csáky et al. (2019).

Hungarian conquerors show three N1a-Z1936 subclades, which is known to stem from the northern Ural region, including the Arctic (likely Palaeo-Laplandic peoples) and cross-stamped cultures of the northern Eurasian forests.

haplogroup_n3a4
Frequency-Distribution Maps of Individual Subclade N3a4 / N1a1a1a1a2-Z1936, probably with the Samic (first) and Fennic (later) expansions into Paleo-Lakelandic and Palaeo-Laplandic territories.

On the other hand, the two N1a-M2118 lineages are more clearly associated with Palaeo-Siberian populations east of the Urals, but became incorporated into the Ugric stock in the Trans-Urals region probably in the same way as N1a-Z1936, by infiltration from (and acculturation of) hunter-gatherers of forest and taiga cultures.

NOTE. You can read more about the infiltration of N1a lineages in the recent post Corded Ware—Uralic (IV): Hg R1a and N in Finno-Ugric and Samoyedic expansions, and in the specific sections for each Uralic group in A Clash of Chiefs.

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

Conclusion

The picture offered by the paper on Hungarian Conquerors, while in line with historical accounts of multi-ethnic tribes incorporating regional lineages, shows nevertheless patrilineal clans clearly associated with Uralic peoples, in a distribution which could have been easily inferred from ancient Trans-Uralian forest-steppe cultures and modern samples (even regarding I2a-L621).

In spite of this, there is a great deal of discussion in the paper about specific N1a subclades in Hungarian conquerors, while the presence of R1a-Z280 (among early Magyar elites!) is interpreted, as always, as recently acculturated Slavs. This is sadly coupled with the simplistic identification of I2a-L621 as of local origin around the Carpathians.

The introduction of the paper to the history of Hungarians is also weird, for example giving credibility to the mythic accounts of the Árpád dynasty’s origin in Attila, which is in line, I guess, with what the authors intended to support all along, i.e. the association of Magyars with Turks from the Eurasian steppes, which they are apparently willing to achieve by relating them to haplogroup R1a-Z93

The conclusion is thus written to appease modern nation-building myths more than anything else, like many other papers before it:

It is generally accepted that the Hungarian language was brought to the Carpathian Basin by the Conquerors. Uralic speaking populations are characterized by a high frequency of Y-Hg N, which have often been interpreted as a genetic signal of shared ancestry. Indeed, recently a distinct shared ancestry component of likely Siberian origin was identified at the genomic level in these populations, modern Hungarians being a puzzling exception36. The Conqueror elite had a significant proportion of N Hgs, 7% of them carrying N1a1a1a1a4-M2118 and 10% N1a1a1a1a2-Z1936, both of which are present in Ugric speaking Khantys and Mansis. At the same time none of the examined Conquerors belonged to the L1034 subclade of Z1936, while all of the Khanty Z1936 lineages reported in 37 proved to be L1034 which has not been tested in the 23 study. Population genetic data rather position the Conqueror elite among Turkic groups, Bashkirs and Volga Tatars, in agreement with contemporary historical accounts which denominated the Conquerors as “Turks”. This does not exclude the possibility that the Hungarian language could also have been present in the obviously very heterogeneous, probably multiethnic Conqueror tribal alliance.

So, back to square one, and new circular reasoning: If ancient populations from north-eastern Europe believed to represent ancient Finno-Ugrians are of R1a-Z645 lineages, it’s because they were not Finno-Ugric speakers. If ancient and modern populations known to be of Finno-Ugric language show clear connections with R1a-Z645, it’s because they are “multi-ethnic”.

The only stable basis for discussion in genetic papers, apparently, is the own making of geneticists, with their traditional 2000s “R1a=Indo-European” and “N1c=Uralic”, coupled with national beliefs. It does not matter how many predictions based on that have been proven wrong, or how many predictions based on the Corded Ware = Uralic expansion have been proven right.

Related

Ancient Sardinia hints at Mesolithic spread of R1b-V88, and Western EEF-related expansion of Vasconic

nuragic-sardinia-neolithic

New preprint Population history from the Neolithic to present on the Mediterranean island of Sardinia: An ancient DNA perspective, by Marcus et al. bioRxiv (2019)

Interesting excerpts (emphasis mine, edited for clarity):

On the high frequency of R1b-V88

Our genome-wide data allowed us to assign Y haplogroups for 25 ancient Sardinian individuals. More than half of them consist of R1b-V88 (n=10) or I2-M223 (n=7).

Francalacci et al. (2013) identi fied three major Sardinia-specifi c founder clades based on present-day variation within the haplogroups I2-M26, G2-L91 and R1b-V88, and here we found each of those broader haplogroups in at least one ancient Sardinian individual. Two major present-day Sardinian haplogroups, R1b-M269 and E-M215, are absent.

Compared to other Neolithic and present-day European populations, the number of identi fied R1b-V88 carriers is relatively high.

(…)ancient Sardinian mtDNA haplotypes belong almost exclusively to macro-haplogroups HV (n = 16), JT (n = 17) and U (n = 9), a composition broadly similar to other European Neolithic populations.

r1b-v88-europe
Geographic and temporal distribution of R1b-V88 Y-haplotypes in ancient European samples. We plot the geographic position of all ancient samples inferred to carry R1b-V88 equivalent markers. Dates are given as years BCE (means of calibrated 2s radio-carbon dates). Multiple V88 individuals with similar geographic positions are vertically stacked. We additionally color-code the status of the R1b-V88 subclade R1b-V2197, which is found in most present-day African R1b-V88 carriers.

On the origin of a Vasconic-like Paleosardo with the Western EEF

(…) the Neolithic (and also later) ancient Sardinian individuals sit between early Neolithic Iberian and later Copper Age Iberian populations, roughly on an axis that differentiates WHG and EEF populations and embedded in a cluster that additionally includes Neolithic British individuals. This result is also evident in terms of absolute genetic differentiation, with low pairwise FST ~ 0.005 +- 0.002 between Neolithic Sardinian individuals and Neolithic western mainland European populations. Pairwise outgroup-f3 analysis shows a very similar pattern, with the highest values of f3 (i.e. most shared drift) being with Neolithic and Copper Age Iberia, gradually dropping off for temporally and geographically distant populations.

In explicit admixture models (using qpAdm, see Methods) the southern French Neolithic individuals (France-N) are the most consistent with being a single source for Neolithic Sardinia (p ~ 0:074 to reject the model of one population being the direct source of the other); followed by other populations associated with the western Mediterranean Neolithic Cardial Ware expansion.

sardinians-ancient-eef
Principal Components Analysis based on the Human Origins dataset. A: Projection of ancient individuals’ genotypes onto principal component axes de fined by modern Western Eurasians (gray labels).

Pervasive Western Hunter-Gatherer ancestry in Iberian/French/Sardinian population

Similar to western European Neolithic and central European Late Neolithic populations, ancient Sardinian individuals are shifted towards WHG individuals in the top two PCs relative to early Neolithic Anatolians Admixture analysis using qpAdm infers that ancient Sardinian individuals harbour HG ancestry (~ 17%) that is higher than early Neolithic mainland populations (including Iberia, ~ 8%), but lower than Copper Age Iberians (~ 25%) and about the same as Southern French Middle-Neolithic individuals (~ 21%).

sardinia-modern-ancient-nuragic-pca
Principal Components Analysis based on the Human Origins dataset. B: Zoom into the region most relevant for Sardinian individuals.

Continuity from Sardinia Neolithic through the Nuragic

We found several lines of evidence supporting genetic continuity from the Sardinian Neolithic into the Bronze Age and Nuragic times. Importantly, we observed low genetic differentiation between ancient Sardinian individuals from various time periods.

A qpAdm analysis, which is based on simultaneously testing f-statistics with a number of outgroups and adjusts for correlations, cannot reject a model of Neolithic Sardinian individuals being a direct predecessor of Nuragic Sardinian individuals (…) Our qpAdm analysis further shows that the WHG ancestry proportion, in a model of admixture with Neolithic Anatolia, remains stable at ~17% throughout three ancient time-periods.

sardinians-modern-ancient-pca-admixture
Present-day genetic structure in Sardinia reanalyzed with aDNA. A: Scatter plot of the rst two principal components trained on 1577 present-day individuals with grand-parental ancestry from Sardinia. Each individual is labeled with a location if at least 3 of the 4 grandparents were born in the same geographical location (\small” three letter abbreviations); otherwise with \x” or if grand-parental ancestry is missing with \?”. We calculated median PC values for each Sardinian province (large abbreviations). We also projected each ancient Sardinian individual on to the top two PCs (gray points). B/C: We plot f-statistics that test for admixture of modern Sardinian individuals (grouped into provinces) when using Nuragic Sardinian individuals as one source population. Uncertainty ranges depict one standard error (calculated from block bootstrap). Karitiana are used in the f-statistic calculation as a proxy for ANE/Steppe ancestry (Patterson et al., 2012).

Steppe influx in Modern Sardinians

While contemporary Sardinian individuals show the highest affinity towards EEF-associated populations among all of the modern populations, they also display membership with other clusters (Fig. 5). In contrast to ancient Sardinian individuals, present-day Sardinian individuals carry a modest “Steppe-like” ancestry component (but generally less than continental present-day European populations), and an appreciable broadly “eastern Mediterranean” ancestry component (also inferred at a high fraction in other present-day Mediterranean populations, such as Sicily and Greece).

Related

ASoSaH Reread (II): Y-DNA haplogroups among Uralians (apart from R1a-M417)

corded-ware-yamna-ancestry

This is mainly a reread of from Book Two: A Game of Clans of the series A Song of Sheep and Horses: chapters iii.5. Early Indo-Europeans and Uralians, iv.3. Early Uralians, v.6. Late Uralians and vi.3. Disintegrating Uralians.

“Sredni Stog”

While the true source of R1a-M417 – the main haplogroup eventually associated with Corded Ware, and thus Uralic speakers – is still not known with precision, due to the lack of R1a-M198 in ancient samples, we already know that the Pontic-Caspian steppes were probably not it.

We have many samples from the north Pontic area since the Mesolithic compared to the Volga-Ural territory, and there is a clear prevalence of I2a-M223 lineages in the forest-steppe area, mixed with R1b-V88 (possibly a back-migration from south-eastern Europe).

R1a-M459 (xR1a-M198) lineages appear from the Mesolithic to the Chalcolithic scattered from the Baltic to the Caucasus, from the Dniester to Samara, in a situation similar to haplogroups Q1a-M25 and R1b-L754, which supports the idea that R1a, Q1a, and R1b expanded with ANE ancestry, possibly in different waves since the Epipalaeolithic, and formed the known ANE:EHG:WHG cline.

y-dna-khvalynsk
Y-DNA samples from Khvalynsk and neighbouring cultures. See full version.

The first confirmed R1a-M417 sample comes from Alexandria, roughly coinciding with the so-called steppe hiatus. Its emergence in the area of the previous “early Sredni Stog” groups (see the mess of the traditional interpretation of the north Pontic groups as “Sredni Stog”) and its later expansion with Corded Ware supports Kristiansen’s interpretation that Corded Ware emerged from the Dnieper-Dniester corridor, although samples from the area up to ca. 4000 BC, including the few Middle Eneolithic samples available, show continuity of hg. I2a-M223 and typical Ukraine Neolithic ancestry.

NOTE. The further subclade R1a-Z93 (Y26) reported for the sample from Alexandria seems too early, given the confidence interval for its formation (ca. 3500-2500 BC); even R1a-Z645 could be too early. Like the attribution of the R1b-L754 from Khvalynsk to R1b-V1636 (after being previously classifed as of Pre-V88 and M73 subclade), it seems reasonable to take these SNP calls with a pinch of salt: especially because Yleaf (designed to look for the furthest subclade possible) does not confirm for them any subclade beyond R1a-M417 and R1b-L754, respectively.

The sudden appearance of “steppe ancestry” in the region, with the high variability shown by Ukraine_Eneolithic samples, suggests that this is due to recent admixture of incoming foreign peoples (of Ukraine Neolithic / Comb Ware ancestry) with Novodanilovka settlers.

The most likely origin of this population, taking into account the most common population movements in the area since the Neolithic, is the infiltration of (mainly) hunter-gatherers from the forest areas. That would confirm the traditional interpretation of the origin of Uralic speakers in the forest zone, although the nature of Pontic-Caspian settlers as hunter-gatherers rather than herders make this identification today fully unnecessary (see here).

EDIT (3 FEB 2019): As for the most common guesstimates for Proto-Uralic, roughly coinciding with the expansion of this late Sredni Stog community (ca. 4000 BC), you can read the recent post by J. Pystynen in Freelance Reconstruction, Probing the roots of Samoyedic.

eneolithic-ukraine-corded-ware
Late Sredni Stog admixture shows variability proper of recent admixture of forest-steppe peoples with steppe-like population. See full version here.

NOTE. Although my initial simplistic interpretation (of early 2017) of Comb Ware peoples – traditionally identified as Uralic speakers – potentially showing steppe ancestry was probably wrong, it seems that peoples from the forest zone – related to Comb Ware or neighbouring groups like Lublyn-Volhynia – reached forest-steppe areas to the south and eventually expanded steppe ancestry into east-central Europe through the Volhynian Upland to the Polish Upland, during the late Trypillian disintegration (see a full account of the complex interactions of the Final Eneolithic).

The most interesting aspect of ascertaining the origin of R1a-M417, given its prevalence among Uralic speakers, is to precisely locate the origin of contacts between Late Proto-Indo-European and Proto-Uralic. Traditionally considered as the consequence of contacts between Middle and Upper Volga regions, the most recent archaeological research and data from ancient DNA samples has made it clear that it is Corded Ware the most likely vector of expansion of Uralic languages, hence these contacts of Indo-Europeans of the Volga-Ural region with Uralians have to be looked for in neighbours of the north Pontic area.

sredni-stog-repin-contacts
Sredni Stog – Repin contacts representing Uralic – Late Indo-European contacts were probably concentrated around the Don River.

My bet – rather obvious today – is that the Don River area is the source of the earliest borrowings of Late Uralic from Late Indo-European (i.e. post-Indo-Anatolian). The borrowing of the Late PIE word for ‘horse’ is particularly interesting in this regard. Later contacts (after the loss of the initial laryngeal) may be attributed to the traditionally depicted Corded Ware – Yamna contact zone in the Dnieper-Dniester area.

NOTE. While the finding of R1a-M417 populations neighbouring R1b-L23 in the Don-Volga interfluve would be great to confirm these contacts, I don’t know if the current pace of more and more published samples will continue. The information we have right now, in my opinion, suffices to support close contacts of neighbouring Indo-Europeans and Uralians in the Pontic-Caspian area during the Late Eneolithic.

Classical Corded Ware

After some complex movements of TRB, late Trypillia and GAC peoples, Corded Ware apparently emerged in central-east Europe, under the influence of different cultures and from a population that probably (at least partially) stemmed from the north Pontic forest-steppe area.

Single Grave and central Corded Ware groups – showing some of the earliest available dates (emerging likely ca. 3000/2900 BC) – are as varied in their haplogroups as it is expected from a sink (which does not in the least resemble the Volga-Ural population):

Interesting is the presence of R1b-L754 in Obłaczkowo, potentially of R1b-V88 subclade, as previously found in two Central European individuals from Blätterhole MN (ca. 3650 and 3200 BC), and in the Iron Gates and north Pontic areas.

Haplogroups I2a and G have also been reported in early samples, all potentially related to the supposed Corded Ware central-east European homeland, likely in southern Poland, a region naturally connected to the north Pontic forest-steppe area and to the expansion of Neolithic groups.

corded-ware-haplogroups
Y-DNA samples from early Corded Ware groups and neighbouring cultures. See full version.

The true bottlenecks under haplogroup R1a-Z645 seem to have happened only during the migration of Corded Ware to the east: to the north into the Battle Axe culture, mainly under R1a-Z282, and to the south into Middle Dnieper – Fatyanovo-Balanovo – Abashevo, probably eventually under R1a-Z93.

This separation is in line with their reported TMRCA, and supports the split of Finno-Permic from an eastern Uralic group (Ugric and Samoyedic), although still in contact through the Russian forest zone to allow for the spread of Indo-Iranian loans.

This bottleneck also supports in archaeology the expansion of a sort of unifying “Corded Ware A-horizon” spreading with people (disputed by Furholt), the disintegrating Uralians, and thus a source of further loanwords shared by all surviving Uralic languages.

Confirming this ‘concentrated’ Uralic expansion to the east is the presence of R1a-M417 (xR1a-Z645) lineages among early and late Single Grave groups in the west – which essentially disappeared after the Bell Beaker expansion – , as well as the presence of these subclades in modern Central and Western Europeans. Central European groups became thus integrated in post-Bell Beaker European EBA cultures, and their Uralic dialect likely disappeared without a trace.

NOTE. The fate of R1b-L51 lineages – linked to North-West Indo-Europeans undergoing a bottleneck in the Yamna Hungary -> Bell Beaker migration to the west – is thus similar to haplogroup R1a-Z645 – linked to the expansion of Late Uralians to the east – , hence proving the traditional interpretation of the language expansions as male-driven migrations. These are two of the most interesting genetic data we have to date to confirm previous language expansions and dialectal classifications.

It will be also interesting to see if known GAC and Corded Ware I2a-Y6098 subclades formed eventually part of the ancient Uralic groups in the east, apart from lineages which will no doubt appear among asbestos ware groups and probably hunter-gatherers from north-eastern Europe (see the recent study by Tambets et al. 2018).

Corded Ware ancestry marked the expansion of Uralians

Sadly, some brilliant minds decided in 2015 that the so-called “Yamnaya ancestry” (now more appropriately called “steppe ancestry”) should be associated to ‘Indo-Europeans’. This is causing the development of various new pet theories on the go, as more and more data contradicts this interpretation.

There is a clear long-lasting cultural, populational, and natural barrier between Yamna and Corded Ware: they are derived from different ancestral populations, which show clearly different ancestry and ancestry evolution (although they did converge to some extent), as well as different Y-DNA bottlenecks; they show different cultures, including those of preceding and succeeding groups, and evolved in different ecological niches. The only true steppe pastoralists who managed to dominate over grasslands extending from the Upper Danube to the Altai were Yamna peoples and their cultural successors.

corded-ware-yamna-pca
Corded Ware admixture proper of expanding late Sredni Stog-like populations from the forest-steppe. See full version here.

NOTE. You can also read two recent posts by FrankN in the blog aDNA era, with detailed information on the Pontic-Caspian cultures and the formation of “steppe ancestry” during the Palaeolithic, Mesolithic and Neolithic: How did CHG get into Steppe_EMBA? Part 1: LGM to Early Holocene and How did CHG get into Steppe_EMBA? Part 2: The Pottery Neolithic. Unlike your typical amateur blogger on genetics using few statistical comparisons coupled with ‘archaeolinguoracial mumbo jumbo’ to reach unscientific conclusions, these are obviously carefully redacted texts which deserve to be read.

I will not enter into the discussion of “steppe ancestry” and the mythical “Siberian ancestry” for this post, though. I will just repost the opinion of Volker Heyd – an archaeologist specialized in Yamna Hungary and Bell Beakers who is working with actual geneticists – on the early conclusions based on “steppe ancestry”:

[A]rchaeologist Volker Heyd at the University of Bristol, UK, disagreed, not with the conclusion that people moved west from the steppe, but with how their genetic signatures were conflated with complex cultural expressions. Corded Ware and Yamnaya burials are more different than they are similar, and there is evidence of cultural exchange, at least, between the Russian steppe and regions west that predate Yamnaya culture, he says. None of these facts negates the conclusions of the genetics papers, but they underscore the insufficiency of the articles in addressing the questions that archaeologists are interested in, he argued. “While I have no doubt they are basically right, it is the complexity of the past that is not reflected,” Heyd wrote, before issuing a call to arms. “Instead of letting geneticists determine the agenda and set the message, we should teach them about complexity in past human actions.

Related

Resurge of local populations in the final Corded Ware culture period from Poland

poland-kujawy

Open access A genomic Neolithic time transect of hunter-farmer admixture in central Poland, by Fernandes et al. Scientific Reports (2018).

Interesting excerpts (emphasis mine, stylistic changes):

Most mtDNA lineages found are characteristic of the early Neolithic farmers in south-eastern and central Europe of the Starčevo-Kőrös-Criş and LBK cultures. Haplogroups N1a, T2, J, K, and V, which are found in the Neolithic BKG, TRB, GAC and Early Bronze Age samples, are part of the mitochondrial ‘Neolithic package’ (which also includes haplogroups HV, V, and W) that was introduced to Europe with farmers migrating from Anatolia at the onset of the Neolithic17,31.

A noteworthy proportion of Mesolithic haplogroup U5 is also found among the individuals of the current study. The proportion of haplogroup U5 already present in the earliest of the analysed Neolithic groups from the examined area differs from the expected pattern of diversity of mtDNA lineages based on a previous archaeological view and on the aDNA findings from the neighbouring regions which were settled by post-Linear farmers similar to BKG at that time. A large proportion of Mesolithic haplogroups in late-Danubian farmers in Kuyavia was also shown in previous studies concerning BKG samples based on mtDNA only, although these frequencies were derived on the basis of very small sample sizes.

y-dna-poland

A significant genetic influence of HG populations persisted in this region at least until the Eneolithic/Early Bronze Age period, when steppe migrants arrived to central Europe. The presence of two outliers from the middle and late phases of the BKG in Kuyavia associated with typical Neolithic burial contexts provides evidence that hunter-farmer contacts were not restricted to the final period of this culture and were marked by various episodes of interaction between two societies with distinct cultural and subsistence differences.

The identification of both mitochondrial and Y-chromosome haplogroup lineages of Mesolithic provenance (U5 and I, respectively) in the BKG support the theory that both male and female hunter-gatherers became part of these Neolithic agricultural societies, as has been reported for similar cases from the Carpathian Basin, and the Balkans. The identification of an individual with WHG affinity, dated to ca. 4300 BCE, in a Middle Neolithic context within a BKG settlement, provides direct evidence for the regional existence of HG enclaves that persisted and coexisted at least for over 1000 years, from the arrival of the LBK farmers ca. 5400 BCE until ca. 4300 BCE, in proximity with Neolithic settlements, but without admixing with their inhabitants.

poland-pca
Principal component analysis with modern populations greyed out on the background (top), ADMIXTURE results with K = 10 with samples from this study amplified (bottom).

The analysis of two Late Neolithic cultures, the GAC and CWC, shows that steppe ancestry was present only among the CWC individuals analysed, and that the single GAC individual had more WHG ancestry than previous local Neolithic individuals. (…) The CWC’s affinity to WHG, however, contrasts with results from published CWC individuals that identified steppe ancestry related to Yamnaya as the major contributor to the CWC genomes, while here we report also substantial contributions from WHG that could relate to the late persistence of pockets of WHG populations, as supported by the admixture results of N42 and the finding of the 4300-year-old N22 HG individual. These results agree with archaeological theories that suggest that the CWC interaction with incoming steppe cultures was complex and that it varied by region.

Some comments

About the analyzed CWC samples, it is remarkable that, even though they are somehow related to each other, they do not form a tight cluster. Also, their Y-DNA (I2a), and this:

When compared to previously published CWC data, our CWC group (not individuals) is genetically significantly closer to WHG than to steppe individuals (Z = −4.898), a result which is in contrast with those for CWC from Germany (Z = 2.336), Estonia (Z = 0.555), and Latvia (Z = 1.553).

ancestry-proportions-poland
Ancestry proportions based on qpAdm. Visual representation of the main results presented in Supplementary Table S5. Populations from this study marked with an asterisk. Values and populations in brackets show the nested model results marked in green in Supplementary Table S5.

Włodarczak (2017) talks about the CWC period in Poland after ca. 2600 BC as a time of emergence of an allochthnous population, marked by the rare graves of this area, showing infiltrations initially mainly from Lesser Poland, and later (after 2500 BC) from the western Baltic zone.

Since forest sub-Neolithic populations would have probably given more EHG to the typical CWC population, these samples support the resurge of ‘local’ pockets of GAC- or TRB-like groups with more WHG (and also Levant_Neolithic) ancestry.

The known presence of I2a2a1b lineages in GAC groups in Poland also supports this interpretation, and the subsistence of such pockets of pre-steppe-like populations is also seen with the same or similar lineages appearing in comparable ‘resurge’ events in Central Europe, e.g. in samples from the Únětice and Tumulus culture.

About the Bronze Age sample, we have at last official confirmation of haplogroup R1a1a (sadly no subclade*) at the very beginning of the Trzciniec period – in a region between western (Iwno) and eastern (Strzyżów) groups related to Mierzanowice – , which has to be put in relation with the samples from the final Trzciniec period in the Baltic published in Mittnik et al. (2018).

EDIT (8 OCT 2018): More specific subclades have been published, including a R1a-Z280 lineage for the Bronze Age sample (see spreadsheet).

This confirms the early resurge of R1a-Z645 (probably R1a-Z282) lineages at the core of the developing East European Bronze Age, a province of the European Bronze Age that emerged from evolving Bell Beaker groups in Poland.

bell-beakers-poland-kujawy
Arrival of Bell Beakers in Poland after ca. 2400 BC, and their origin in other BBC centres (Czebreszuk and Szmyt 2011).

I don’t have any hope that the Balto-Slavic evolution through BBC Poland → Mierzanowice/Iwno → Trzciniec → Lusatian cultures is going to be confirmed any time soon, until we have a complete trail of samples to follow all the way to historic Slavs of the Prague culture. However, I do think that the current data on central-east Europe – and the recent data we are receiving from north-east Europe and the Iranian steppes, at odds with the Indo-Slavonic alternative – supports this model.

I guess that, in the end, similar to how the Yamna vs. Corded Ware question is being solved, the real route of expansion of Proto-Balto-Slavic (supposedly spoken ca. 1500-1000 BC) is probably going to be decided by the expansion of either R1a-M458 (from the west) or R1a-Z280 lineages (from the east), because the limited precision of genetic data and analyses available today are going to show ‘modern Slavic’-like populations from the whole eastern half of Europe for the past 4,000 years…

Related