Corded Ware and Bell Beaker related groups defined by patrilocality and female exogamy

tumulus-culture-eba-danube

Two new interesting papers concerning Corded Ware and Bell Beaker peoples appeared last week, supporting yet again what is already well-known since 2015 about West Uralic and North-West Indo-European speakers and their expansion.

Below are relevant excerpts (emphasis mine) and comments.

#UPDATE (27 OCT 2019): I have updated Y-DNA and mtDNA maps of Corded Ware, Bell Beaker, EBA, MBA, and LBA migrations. I have also updated PCA plots, which now include the newly reported samples and those from the Tollense valley, and I have tried some qpAdm models (see below).

I. Corded Ware and Battle Axe cultures

Open access The genomic ancestry of the Scandinavian Battle Axe Culture people and their relation to the broader Corded Ware horizon, by Malmström, Günther, et al. Philos. Trans. R. Soc. (2019).

I.1. Origins of Corded Ware peoples

The discovery of the Alexandria outlier represented a clear support for a long-lasting genomic difference between the two distinct cultural groups, Yamnaya and Corded Ware, already visible in an opposition Khvalynsk vs. late Sredni Stog ca. 4000 BC, i.e. well before the formation of both Late Eneolithic/Early Bronze Age groups.

However, the realization that it may not have been an Eneolithic individual, but rather a (Middle?) Bronze Age one, suggests that Sredni Stog was possibly not directly related to Corded Ware, and a potential direct connection with Yamnaya might have to be reevaluated, e.g. through the Carpathian Basin, as Anthony (2017) proposed.

pca-yamnaya-corded-ware-oblaczkowo
Principal component analysis of modern Europeans (grey) and projected ancient Europeans.

This new paper shows two early Corded Ware individuals from Obłaczkowo, Poland (ca. 2900-2600 BC) – hence close to the supposed original Proto-Corded Ware community – with an apparently (almost) full “Steppe-like” ancestry, clustering (almost) with Yamnaya individuals:

Similar to the BAC individuals, the newly sequenced individuals from the present-day Karlova in Estonia and Obłaczkowo in Poland appear to have strong genetic affinities to other individuals from BAC and CWC contexts across the Baltic Sea region. Some individuals from CWC contexts, including the two from Obłaczkowo, cluster closely with the potential source population of steppe-related ancestry, the Yamnaya herders. Notably, these individuals appear to be those with the earliest radiocarbon dates among all genetically investigated individuals from CWC contexts. Overall, for CWC-associated individuals, there is a clear trend of decreasing affinity to Yamnaya herders with time.

NOTE. Interestingly, this sample is almost certainly attributed to the skeleton E8-A, which had been supposedly already investigated by the Copenhagen group as the RISE1 sample:

We note that RISE1 is also described as the individual from Obłaczkowo feature E8-A. However, their genetic results differ from ours. They present this individual as a molecularly determined male that belongs to Y-chromosomal haplogroup (hg) R1b and to mtDNA hg K1b1a1 while our results show this individual to be female, carrying a mtDNA hg U3a’c profile

Since the typical Steppe_MLBA ancestry of Corded Ware groups does not show good fits for (Pre-)Yamnaya-derived ancestry, it is almost certain that these individuals will show no (or almost no) direct Yamnaya-related contribution, but rather a contribution of East European sub-Neolithic groups, more or less close to the steppe-forest region.

NOTE. They might show contributions from Pre-Yamnaya-influenced Sredni Stog, though, but if they show a contribution of Yamnaya, then they are probably outliers, related to Yamnaya vanguard groups (see image below). And for them to show it, then both sources, Yamnaya and Corded Ware, should be clearly distinguishable from each other and their relative contribution quantifiable in formal stats, something difficult (if not impossible) to ascertain today.

trypillian-yamnaya-influence-baltic
Trypillian routes of influence and Yamnaya culture influences in Central and Central-East Europe during the Late Eneolithic / Early Bronze Age. Images by Klochko (2009).

Their position in the published PCA – a plot apparently affected by projection bias – suggests a cluster in common with early Baltic samples, which are known to show contributions from East European sub-Neolithic populations (see qpAdm values for Baltic CWC samples).

NOTE. Results for previous samples labelled as Poland CWC are unreliable due to their low coverage.

The most interesting aspect about the ancestry shown by these early samples is their further support for an origin of the culture different than Sredni Stog, and for a rejection of the Alexandria outlier as ancestral to them, hence for a Volhynian-Podolian homeland of Proto-Corded Ware peoples, with an ancestry probably more closely related to the late Maykop Steppe- and Trypillian/GAC groups admixed with sub-Neolithic populations of the Eastern European Late Eneolithic.

NOTE. That is, unless there is a reason for the apparent increase in so-called “Steppe-ancestry” during the northward and westward migration of CWC peoples that represents another thing entirely…

#UPDATE (27 OCT 2019): Apparently, the PCA was actually not affected by projection bias:

Sample poz44 clusters ‘to the south’, with other early German ones, but also close to Yamnaya. Its poor coverage makes qpAdm results unreliable, but its common cluster close to central European and eastern CWC groups – despite belonging to the same Obłaczkowo site – supports that it is more representative of the Proto-CWC population than poz81.

Sample poz81 clusters with Yamnaya samples – or at least with the wider, Steppe-related cluster. Nevertheless, analyses with qpAdm – in combination with values obtained for other early Baltic samples – support that the ancestry of poz81 is more closely related to a core Corded Ware population admixed with sub-Neolithic peoples (similar to Samara LN).

NOTE. I have selected Czech CWC as a potential source closer to the Proto-CWC population, similar to models with Baltic samples. Since Czech CWC samples are later than these from Obłaczkowo, I have also checked the reverse model, with Poz81 and GAC Poland as a source for Czech CWC, and the fits are slightly worse. Anyway, ‘better’ or ‘worse’ p-values can’t determine the direction of migration

pca-corded-ware-poland-oblaczkowo-baltic-yamnaya
Detail of the PCA of Eurasian samples, including Corded Ware groups and related clusters, as well as outliers. Also marked is poz81.

I.2. CWC expansion under R1a bottlenecks

The two males in our dataset (ber1 and poz81) belonged to Y-chromosome R1a haplogroups, as do the majority of males (16/24) from the previously published CWC contexts, while a smaller fraction belonged to R1b [3/24] or I2a [3/24] lineages. The R1a haplogroup has not been found among Neolithic farmer populations nor in hunter–gatherer groups in central and western Europe, but it has been reported from eastern European hunter–gatherers and Eneolithic groups. Individuals from the Pontic–Caspian steppe, associated with the Yamnaya Culture, carry mostly R1b and not R1a haplotypes.

Sample poz81 is of basal hg. R1a-CTS4385*, an R1a-M417 subclade, supporting once again that most Corded Ware individuals from western and central European groups expanded under R1a-M417 (xZ645) lineages. The Battle Axe sample from Bergsgraven (ca. 2620-2470 BC) shows a basal hg. R1a-Y2395*, a R1a-Z283 subclade leading to the typically Fennoscandian R1a-Z284.

Both findings further support that typical lineages of West CWC groups, including R1a-M417 (xZ645) subclades, were fully replaced by incoming East Bell Beakers, and that the limited expansion of R1a-Z284 and I1 (the latter found in one newly reported Late Neolithic sample from Sweden) was the outcome of later regional bottlenecks within Scandinavia, after the creation of a maritime dominion by the Bell Beaker elites during the Dagger Period.

I.3. CWC and lactase persistence

(…) one of these individuals (kar1) carried at least one allele (-13910 C->T) associated with lactose tolerance, while the other two individuals (ber1 and poz81) carried at least one ancestral variant each, consistent with previous observations of low levels of lactose tolerance variants in the Neolithic and a slight increase among individuals from CWC contexts.

The fact that two early CWC individuals carry ancestral variants could be said to support the improbability of the individual from Alexandria representing a community ancestral to the Corded Ware community. On the other hand, the late CWC individual from Estonia carries one allele, but it still seems that only Bell Beakers and Steppe-related groups show the necessary two alleles during the Early Bronze Age, which is in line with a late Repin/early Yamnaya-related origin of the successful selection of the trait, consistent with the expansion of their specialized semi-nomadic cattle-breeding economy through the steppe biome during the Late Eneolithic.

rs4988235-lactase-persistence-history
Maps part of the public data used for the post by Iain Mathieson on Lactase Persistence. “By 2500 BP, the allele is present over a band stretching from Ireland to Central Asia at around 50 degrees latitude. This probably reflects the spread of Steppe ancestry populations in which the allele originated. However, the allele is still rare (say less than 1% frequency) over this entire range. It does not become common anywhere until some time in the past 2500 years – when it reaches its present-day high frequency in Britain and Central Europe”.

I.4. West Uralic spread from the East

The BAC groups fit as a sister group to the CWC-associated group from Estonia but not as a sister group to the CWC groups from Poland or Lithuania (|Z| > 3), indicating some differences in ancestry between these CWC groups and BAC. Supervised admixture modelling suggests that BAC may be the CWC-related group with the lowest YAM-related ancestry and with more ancestry from European Neolithic groups.

While the results of the paper are compatible with a migration from either the Eastern or the Western Baltic into Scandinavia, phylogeography and archaeology support that Battle Axe peoples emerged as a Baltic Corded Ware group close to the Vistula that expanded first to the north-east, and then to the west from Finland, continuing mostly unscathed during the whole Bronze Age mostly in eastern Fennoscandia with the development of Balto-Finnic- and Samic-speaking communities.

corded-ware-culture-ancestry-over-time
Correlation between f4(Chimp, LBK, YAM, X), where X is a CWC or BAC individual, and the date (BCE) of each individual. This statistic measures shared drift between CWC and Linear Pottery Culture (LBK) as opposed to YAM and should increase with the higher proportion of Neolithic farmer ancestry in CWC and BAC.

Radiocarbon dating showed that the three individuals from the Öllsjö megalithic tomb derived from later burials, where oll007 (2860–2500 cal BCE) overlaps with the time interval of the BAC, and oll009 and oll010 (1930–1650 cal BCE) fall within the Scandinavian Late Neolithic and Early Bronze Age

For more on how the Pitted Ware culture may have influenced Uralic-speaking Battle Axe peoples earlier than Indo-European-speaking Bell Beakers in Scandinavia, read more about Early Bronze Age Scandinavia and about the emergence of the Pre-Proto-Germanic community.

II. Bell Beakers through the Bronze Age

New paper (behind paywall) Kinship-based social inequality in Bronze Age Europe, by Mittnik et al. Science (2019).

II.1. Yamnaya vanguard settlers

In my last post, I showed how the ancestry of Corded Ware from Esperstedt is consistent with influence by incoming Yamnaya vanguard settlers or early Bell Beakers, stemming ultimately from the Carpathian Basin, something that could be inferred from the position of the Esperstedt outlier in the PCA, and by the knowledge of Yamnaya archaeological influences up to Saxony-Anhalt.

Yamnaya settlers are strongly suspected to have migrated in small so-called vanguard groups to the west and north of the Carpathians in the first half of the 3rd millennium BC, well before the eventual adoption of the Proto-Beaker package and their expansion ca. 2500 BC as East Bell Beakers.

Tauber Valley infiltration

As I mentioned in the books, one of the known – among the many more unknown – sites displaying Yamnaya-related traits and suggesting the expansion of Yamnaya settlers into Central Europe is Lauda-Königshofen, in the Tauber Valley.

From Diet and Mobility in the Corded Ware of Central Europe, by Sjögren, Price, & Kristiansen PLoS One (2017):

A series of CW cemeteries have been excavated in the Tauber valley. There are three large cemeteries known and some 30 smaller sites. The larger ones are Tauberbischofsheim-Dittingheim with 62 individuals, Tauberbischofsheim-Impfingen with 40 individuals, and Lauda-Königshofen with 91 individuals. The cemeteries are dispersed rather regularly along the Tauber valley, on both sides of the river, suggesting a quite densely settled landscape.

The Lauda-Königshofen graves consisted mostly of single inhumations in contracted position, usually oriented E-W or NE-SW. A total of 91 individuals were buried in 69 graves. At least 9 double graves and three graves with 3–4 individuals were present. In contrast to the common CW pattern, sexes were not distinguished by body position, only by grave goods. This trait is common in the Tauber valley and suggests a local burial tradition in this area. Stone axes were restricted to males, pottery to females, while other artifacts were common to both sexes. About a third of the graves were surrounded by ring ditches, suggesting palisade enclosures and possibly over-plowed barrows.

In particular, Frînculeasa, Preda, & Heyd (2015) used Lauda-Königshofen as representative of the mobility of horse-riding Yamnaya nomadic herders migrating into southern Germany, referring to the findings in Trautmann (2012) about the nomadic herders from the Tauber Valley, and their already known differences with other Corded Ware groups.

The likely influence of Yamnaya in the region has been reported at least since the 2000s, repeatedly mentioned by Jozef Bátora (2002, 2003, 2006), who compiled Yamnaya influences in a map that has been copied ever since, with little improvement over time. Heyd believes that there are potentially many Yamnaya remains along the Middle and Lower Danube and tributaries not yet found, though.

NOTE. Looking for this specific site, I realized that Bátora (and possibly many after him who, like me, copied his map) located Lauda-Königshofen in a more south-western position within Baden-Württemberg than its actual location. I have now corrected it in the maps of Chalcolithic migrations.

yamnaya-corded-ware-europe
Yamnaya influences in Central Europe suggestive of vanguard settlements, contemporary with Corded Ware groups. See full map.

Althäuser Hockergrab…Bell Beakers

Unfortunately, though, it is very difficult to attribute the reported R1b-L51 sample from the Tauber valley to a population preceding the arrival of East Bell Beakers in the region, so there is no uncontroversial smoking gun of Yamnaya vanguard settlers – yet. Reasons to doubt a Pre-Beaker origin are as follows:

1. This family of the Tauber valley shows a late radiocarbon date (ca. 2500 BC), i.e. from a time where East Bell Beakers are known to have been already expanding in all directions from the Middle and Upper Danube and its tributaries.

tauber-valley-althauser-hockergrab
Crouched burial from Althausen (Althäuser Hockergrab), dated ca. 2500 BC.

2. Archaeological information is scarce. Remains of these four individuals were discovered in 1939 and officially reported together with other findings in 1950, without any meaningful data that could distinguish between Bell Beakers and Corded Ware individuals.

This site is located in the Tauber valley, ca. 100 km to the northwest of the Lech valley. The site was discovered during the construction of a sports field in 1939 and was subsequently excavated by G. Müller and O. Paret. Four individuals in crouched position were found in the burial pit of a flat grave. The burial did not contain any grave goods, but due to the type of grave and positioning of the bodies (with heads pointing towards southwest) the site was attributed to the Corded Ware complex.

The classification of this burial as of CWC and not BBC seems to have been based entirely on the numerous CWC findings in the Tauber valley, rather than on its particular burial orientation following a regional custom (foreign to the described standard of both cultures), and on its grave type that was also found among Bell Beaker groups. Like many human remains recovered in dubious circumstances in the 20th century, these samples should have probably been labelled (at least in the genetic paper) more properly as Tauber_LN or Tauber_EBA.

yamnaya-bias-tauber-lech-valley
Changes in ancestry over time. (A) Median ages of individuals plotted against z scores of f4 (Mbuti, Test; Yamnaya_Samara, Anatolia_Neolithic) show increase of Anatolian farmer-related ancestry (indicated by more positive z-scores) and decrease of variation in ancestry over time. Grey shading indicates significant z scores, red line shonw near correlation (r = -0.35971; P = 0.003) and dotted lines the 95% confidence interval. (B) ancestry proportions on autosomes calculated with qpAdm. (C) Sex-bias z scores between autosomes and X chromosomes show significant male bias for steppe-related ancestry in the Tauber samples. Image modified from the paper: Surrounded with a blue circle in (A) are females with more Steppe-related ancestry, and in (C) surrounded by squares are the distinct sex biases found in the earliest BBC from the Tauber valley vs. later groups from the Lech valley.

3. In terms of ancestry, there seem to be no gross differences between the Lech Valley BBC individuals and previously reported South German Beakers, originally Yamnaya-like settlers admixing through exogamy with locals, including Corded Ware peoples, as the sex bias of the Lech Valley Beakers proves (see PCA plot below). In other words, northern and eastern Beakers admixed with regional (Epi-)Corded Ware females during their respective expansions, similar to how southern and western Beakers admixed with regional EEF-related females.

The two available Tauber Valley samples (“Tauber_CWC”) show the same pattern: a quite recent Steppe-related male bias and Anatolia_Neolithic-related female bias. Nevertheless, the male sample clusters ‘to the south’ in the PCA relative to all sampled Corded Ware individuals (see PCA plot below), and shows less Yamnaya-like ancestry than what is reported (or can be inferred) for Yamnaya from Hungary or early Bell Beakers of elevated Steppe-related ancestry.

yamnaya-ancestry-tauber-cwc-bbc-lech-eba-mba
Table S9. Three-way qpAdm admixture model for European MN/Chalcolithic group+Yamnaya_Samara. P-values greater than 0.05 (model is not rejected) marked in green.

The ancestry and position of the Althäuser male in the PCA is thus fully compatible with recently incoming East Bell Beakers admixing with local peoples (including Corded Ware) through exogamy, but not so much with a sample that would be expected from Yamanaya vanguard + Corded Ware-related ancestry (more like the Esperstedt outlier or the early France Beaker). Compared to the more ‘northern’ (fully Corded Ware-like) position ancestry of his female counterpart, there is little to support that both are part of the same native Tauber valley community after generations of ancestry levelling…

#UPDATE (27 OCT 2019): The PCA shows that the Althäuser male clusters, in fact, ‘to the north’ of the female one, almost on the same spot as a Bell Beaker sample from the Lech Valley.

Despite their reported damage and poor coverage, there seems to be a trend for qpAdm values to prefer a source population for the male (Alt_4) close to Germany Beakers, whereas the female sample (Alt_3) shows ‘better’ fits when a Corded Ware source is selected.

Also relevant is the Corded Ware ancestry of the male – closer to a Czech rather than German CWC source – compatible with an eastern origin, hence supporting a recent arrival via the Danube, in contrast to the local source of the CWC admixture of the female. The poorer coverage of the female sample makes these results questionable, though.

pca-bell-beaker-tauber-lech-valley-yamnaya-cwc
Detail of the PCA of Eurasian samples, including Bell Beaker groups and related clusters, as well as outliers. Also marked are the Tauber Valley male (M) and female (F).

4. The haplogroup inference is also unrevealing: whereas the paper reports that it is R1b-P310* (xU106, xP312), there is no data to support a xP312 call, so it may well be even within the P312 branch, like most sampled Bell Beaker males. Similarly, the paper also reports that HUGO_180Sk1 (ca. 2340 BC) shows a positive SNP for the U106 trunk, which would make it the earliest known U106 sample and originally from Central Europe, but there is no clear support for this SNP call, either. At least not in their downloadable BAM files, as far as I can tell. Even if both were true, they would merely confirm the path of expansion of Yamnaya / East Bell Beakers through the Danube, already visible in confirmed genomic data:

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

II.2. Proto-Celts and the Tumulus culture

The most interesting data from Mittnik et al. (2019) – overshadowed by the (at first sight) striking “CWC” label of the Althäuser male – is the finding that the most likely (Pre-)Proto-Celtic community of Southern Germany shows, as expected, major genetic continuity over time with Yamnaya/East Bell Beaker-derived patrilineal families, which suggests an almost full replacement of other Y-chromosome haplogroups in Southern German Bronze Age communities, too.

Sampled families form part of an evolving Bell Beaker-derived European BA cluster in common with other Indo-European-speaking cultures from Western, Southern, and Northern Europe, also including early Balto-Slavs, clearly distinct from the Corded Ware-related clusters surviving in the Eastern Baltic and the forest zone.

This Central European Bronze Age continuity is particularly visible in many generations of different patrilocal families practising female exogamy, showing patrilineal inheritance mainly under R1b-P312 (mostly U152+) lineages proper of Central European bottlenecks, all of them apparently following a similar sociopolitical system spanning roughly a thousand years, since the arrival of East Bell Beakers in the region (ca. 2500 BC) until – at least – the end of the Middle Bronze Age (ca. 1300 BC):

Here, we show a different kind of social inequality in prehistory, i.e., complex households that consisted of i) a higher-status core family, passing on wealth and status to descendants, ii) unrelated, wealthy and high-status non-local women and iii) local, low-status individuals. Based on comparisons of grave goods, several of the high-status non-local females could have come from areas inhabited by the Unetice culture, i.e., from a distance of at least 350 km. As the EBA evidence from most of Southern Germany is very similar to the Lech valley, we suggest that social structures comparable to our microregion existed in a much broader area. The EBA households in the Lech valley, however, seem similar to the later historically known oikos, the household sphere of classic Greece, as well as the Roman familia, both comprising the kin-related family and their slaves.

pca-lech-valley-bell-beaker-eba
Genetic structure of Late Neolithic and Bronze Age individuals from southern Germany. (A) Ancient individuals (covered at 20,000 or more SNPs) projected onto principal components defined by 1129 present day west Eurasians (shown in fig. S6); individuals in this study shown with outlines corresponding to their 87Sr/86Sr isotope value (black: consistent with local values, orange: uncertain/intermediate, red: inconsistent with local values). Selected published ancient European individuals are shown without outlines. Image modified from the paper. Surrounded by triangles in cyan, Corded Ware-like females; with a blue triangle, Yamnaya/Early BBC-like sample from the Tauber valley.

NOTE. For those unfamiliar with the usual clusters formed by the different populations in the PCA, you can check similar graphics: PCA with Bell Beaker communities, PCA with Yamnaya settlers from the Carpathians, a similar one from Wang et al. (2019) showing the Yamnaya-Hungary cline, or the chronological PCAs prepared by me for the books.

The gradual increase in local EEF-like ancestry among South Germany EBA and MBA communities over the previous BBC period offers a reasonable explanation as to how Italic and Celtic communities remained in loose contact (enough to share certain innovations) despite their physical separation by the Alps during the Early Bronze Age, and probably why sampled Bell Beakers from France were found to be the closest source of Celts arriving in Iberia during the Urnfield period.

Furthermore, continued contacts with Únětice-related peoples through exogamy also show how Celtic-speaking communities closer to the Danube might have influenced (and might have been influenced by) Germanic-speaking communities of the Nordic Late Neolithic and Bronze Age, helping explain their potentially long-lasting linguistic exchange.

Like other previous Neolithic or Chalcolithic groups that Yamnaya and Bell Beakers encountered in Europe, ancestry related to the Corded Ware culture became part of Bell Beaker groups during their expansion and later during the ancestry levelling in the European Early Bronze Age, which helps us distinguish the evolution of Indo-European-speaking communities in Europe, and suggests likely contacts between different cultural groups separated hundreds of km. from each other.

All in all, there is nothing to support that (epi-)Corded Ware groups might have survived in any way in Central or Western Europe: whether through their culture, their Y-chromosome haplogroups, or their ancestry, they followed the fate of other rapidly expanding groups before them, viz. Funnelbeaker, Baden, or Globular Amphorae cultural groups. This is very much unlike the West Uralic-speaking territory in the Eastern Baltic and the Russian forests, where Corded Ware-related cultures thrived during the Bronze Age.

lech-valley-yamnaya-ancestry-over-time
f4-statistics showing differences in ancestry in populations grouped by period. An increase in affinity to ancestry related to Anatolia Neolithic over time. Males and females grouped together shown as upward and downward pointing triangles, respectively.

Conclusion

It was about time that geneticists caught up with the relevance of Y-DNA bottlenecks when assessing migrations and cultural developments.

From Malmström et al. (2019):

The paternal lineages found in the BAC/CWC individuals remain enigmatic. The majority of individuals from CWC contexts that have been genetically investigated this far for the Y-chromosome belong to Y-haplogroup R1a, while the majority of sequenced individuals of the presumed source population of Yamnaya steppe herders belong to R1b. R1a has been found in Mesolithic and Neolithic Ukraine. This opens the possibility that the Yamnaya and CWC complexes may have been structured in terms of paternal lineages—possibly due to patrilineal inheritance systems in the societies — and that genetic studies have not yet targeted the direct sources of the expansions into central and northern Europe.

From Gibbons (2019), a commentary to Mittnik et al. (2019):

Some of the early farmers studied were part of the Neolithic Bell Beaker culture, named for the shape of their pots. Later generations of Bronze Age men who retained Bell Beaker DNA were high-ranking, buried with bronze and copper daggers, axes, and chisels. Those men carried a Y chromosome variant that is still common today in Europe. In contrast, low-ranking men without grave goods had different Y chromosomes, showing a different ancestry on their fathers’ side, and suggesting that men with Bell Beaker ancestry were richer and had more sons, whose genes persist to the present.

There was no sign of these women’s daughters in the burials, suggesting they, too, were sent away for marriage, in a pattern that persisted for 700 years. The only local women were girls from high-status families who died before ages 15 to 17, and poor, unrelated women without grave goods, probably servants, Mittnik says. Strontium levels from three men, in contrast, showed that although they had left the valley as teens, they returned as adults.

Also, from Scientific American:

(…) it has long been assumed that prior to the Athenian and Roman empires,—which arose nearly 2,500 and more than 2,000 years ago, respectively—human social structure was relatively straightforward: you had those who were in power and those who were not. A study published Thursday in Science suggests it was not that simple. As far back as 4,000 years ago, at the beginning of the Bronze Age and long before Julius Caesar presided over the Forum, human families of varying status levels had quite intimate relationships. Elites lived together with those of lower social classes and women who migrated in from outside communities. It appears early human societies operated in a complex, class-based system that propagated through generations.

It seems wrong (to me, at least) that the author and – as he believes – archaeologists and historians had “assumed” a different social system for the European Bronze Age, which means they hadn’t read about how Indo-European societies were structured. For example, long ago Benveniste (1969) already drew some coherent picture of these prehistoric peoples based on their reconstructed language alone: regarding their patrilocal and patrilineal family system; regarding their customs of female exogamy and marriage system; and regarding the status of foreigners and slaves as movable property in their society.

A long-lasting and pervasive social system of Bronze Age elites under Yamnaya lineages strikingly similar to this Southern German region can be easily assumed for the British Isles and Iberia, and it is likely to be also found in the Low Countries, Northern Germany, Denmark, Italy, France, Bohemia and Moravia, etc., but also (with some nuances) in Southern Scandinavia and Central-East Europe during the Bronze Age.

Therefore, only the modern genetic pool of some border North-West Indo-European-speaking communities of Europe need further information to describe a precise chain of events before their eventual expansion in more recent times:

  1. the relative geographical isolation causing the visible regional founder effects in Scandinavia, proper of the maritime dominion of the Nordic Late Neolithic (related thus to the Island Biogeography Theory); and
  2. the situation of the (Pre-)Proto-Balto-Slavic community close to the Western Baltic which, I imagine, will be shown to be related to a resurge of local lineages, possibly due to a shift of power structures similar to the case described for Babia Góra.

NOTE. Rumour has it that R1b-L23 lineages have already been found among Mycenaeans, while they haven’t been found among sampled early West European Corded Ware groups, so the westward expansion of Indo-European-speaking Yamnaya-derived peoples mainly with R1b-L23 lineages through the Danube Basin merely lacks official confirmation.

Related

Bell Beakers and Mycenaeans from Yamnaya; Corded Ware from the forest steppe

eba-yamnaya-ancestry-hungary

I have recently written about the spread of Pre-Yamnaya or Yamnaya ancestry and Corded Ware-related ancestry throughout Eurasia, using exclusively analyses published by professional geneticists, and filling in the gaps and contradictory data with the most reasonable interpretations. I did so consciously, to avoid any suspicion that I was interspersing my own data or cherry picking results.

Now I’m finished recapitulating the known public data, and the only way forward is the assessment of these populations using the available datasets and free tools.

Understanding the complexities of qpAdm is fairly difficult without a proper genetic and statistical background, which I won’t pretend to have, so its tweaking to get strictly correct results would require an unending game of trial and error. I have sadly little time for this, even taking my tendency to procrastination into account… so I have used a simple model akin to those published before – in particular, the outgroup selection by Ning, Wang et al. (2019), who seem to be part of the only group interested in distinguishing Yamnaya-related from Corded Ware-related ancestry, probably the most relevant question discussed today in population genomics regarding the Proto-Indo-European and Proto-Uralic homelands.

eneolithic-steppe-best-fits
Supplementary Table 13. P values of rank=2 and admixture proportions in modelling Steppe ancestry populations as a three-way admixture of Eneolithic steppe Anatolian_Neolithic and WHG using 14 outgroups.
Left populations: Test, Eneolithic_steppe, Anatolian_Neolithic, WHG.
Right populations: Mbuti.DG, Ust_Ishim.DG, Kostenki14, MA1, Han.DG, Papuan.DG, Onge.DG, Villabruna, Vestonice16, ElMiron, Ethiopia_4500BP.SG, Karitiana.DG, Natufian, Iran_Ganj_Dareh_Neolithic.

I have used for all analyses below a merged dataset including the curated one of the Reich Lab, the latest on Central and South Asia by Narasimhan, Patterson et al. (2019), on Iberia by Olalde et al. (2019), and on the East Baltic by Saag et al. (2019), as well as datasets including samples from Wang et al. (2019) and Lamnidis et al. (2018). I used (and intend to use) the same merged dataset in all cases, despite its huge size, to avoid adding one more uncontrolled variable to the analyses, so that all results obtained can be compared.

I try to prepare in advance a bunch of relevant files with left pops and right pops for each model:

  1. It seems a priori more reasonable to use geographically and chronologically closer proxy populations (say, Trypillia or GAC for Steppe-related peoples) than hypothetic combinations of ancestral ones (viz. Anatolian farmer, WHG, and EHG).
  2. This also means using subgroups closer to the most likely source population, such as (Don-Volga interfluve) Yamnaya_Kalmykia rather than (Middle Volga) Yamnaya_Samara for the western expansion of late Repin/early Yamnaya, or the early Germany_Corded_Ware.SG or Czech_Corded Ware for the group closest to the Proto-Corded Ware population (see below), likely neighbouring the Upper Vistula region.
  3. I usually test two source populations for different targets, which seems like a much more efficient way of using computer resources, whenever I know what I want to test, since I need my PC back for its normal use; whenever I don’t know exactly what to test, I use three-way admixture models and look for subsets to try and improve the results.

I have probably left out some more complex models by individualizing the most relevant groups, but for the time being this would have to do. Also, no other formal stats have been used in any case, which is an evident shortcoming, ruling out an interpretation drawn directly and only from the results below.

Full qpAdm results for each batch of samples are presented in a Google Spreadsheet, with each tab (bottom of the page) showing a different combination of sources, usually in order of formally ‘best’ (first to the left) to ‘worst’ (last to the right) fits, although the order is difficult to select in highly heterogeneous target groups, as will be readily visible.

maykop-trypillia-intrusion-steppes
Disintegration, migration, and imports of the Azov–Black Sea region. First migration event (solid arrows): Gordineşti–Maikop expansion (groups: I – Bursuchensk; II – Zhyvotylivka; III – Vovchans’k; IV – Crimean; V – Lower Don; VI – pre-Kuban). Second migration event (hollow arrows): Repin expansion. After Rassamakin (1999), Demchenko (2016).

Corded Ware origins

The latest publications on the Yampil barrow complex have not improved much our understanding of the complexity of Corded Ware origins from an archaeological point of view, involving multiple cultural (hence likely population) influences. This bit is from Ivanova et al., Baltic-Pontic Studies (2015) 20:1, and most hypotheses of the paper remain unanswered (except maybe for the relevance of the Złota group):

In the light of the above outline therefore one should argue that the ‘architecture of barrows’ associated in the ‘Yampil landscape’ of the Middle Dniester Area with the Eneolithic (specifically, mainly with the TC), precedes the development of a similar phenomenon that can be observed from 2900/2800 BC in the Upper Dniester Area and drainage basin of the Upper Vistula, associated with the CWC [Goslar et al. 2015; Włodarczak 2006; 2007; 2008; Jarosz, Włodarczak 2007]. The most consuming research question therefore is whether ritual customs making use of Eneolithic (Tripolye) ‘barrow architecture’ could have penetrated northwards along the Dniester route, where GAC communities functioned. One could also ask what role the rituals played among the autochthons [Kośko 2000; Włodarczak 2008; 2014: 335; Ivanova, Toshchev 2015b].

This issue has already been discussed with a resulting tentative systemic taxonomy in the studies of Włodarczak, arguing for the Złota culture (ZC) in the Vistula region as an illustration of one of the (Małopolska) reception centres of civilization inspirations from the oldest Pontic ‘barrow culture’ circle associated with the Eneolithic and Early Bronze Age [Włodarczak 2008]. Notably, it is in the ZC that one can notice a set of cultural traits (catacomb grave construction, burial details, forms and decoration of vessels) analogous to those shared by the north-western Black Sea Coast groups of the forest-steppe Eneolithic (chiefly Zhyvotilovka-Volchansk) and the Late Tripolye circle (chiefly Usatovo-Gordinești-Horodiștea-Kasperovtsy).

gac-trypillia-usatovo-corded-ware
Globular Amphorae culture „exodus” to the Danube Delta: a – Globular Amphorae culture; b – GAC (1), Gorodsk (2), Vykhvatintsy (3) and Usatovo (4) groups of Trypillia culture; c – Coţofeni culture; d – northern border of the late phase of Baden culture;red arrows – direction of Globular Amphora culture expansion; blue arrow – direction of „reflux” of Globular Amphora culture (apud Włodarczak, 2008, with changes).

Taking into account that I6561 might be wrongly dated, we cannot include the Corded Ware-like sample of the end-5th millennium BC in the analysis of Corded Ware origins. That uncertainty in the chronology of the appearance of “Steppe ancestry” in Proto-Corded Ware peoples complicates the selection of any potential source population from the CHG cline.

Nevertheless, the lack of hg. R1a-M417 and sizeable Pre-Yamnaya-related ancestry in the sampled Pontic forest-steppe Eneolithic populations (represented exclusively by two samples from Dereivka ca. 3600-3400 BC) would leave open the interesting possibility that a similar ancestry got to the forest-steppe region between modern Poland and Ukraine during the known complex population movements of the Late Eneolithic.

It is known that Corded Ware-derived groups and Steppe Maykop show bad fits for Pre-Yamnaya/Yamnaya ancestry, and also that Steppe Maykop is a potential source of “Steppe-related ancestry” within the Eneolithic CHG mating network of the Pontic-Caspian steppes and forest-steppes. Testing Corded Ware for recent Trypillia and Maykop influences, proper of Late Trypillia and Late Maykop groups in the North Pontic area (such as Zhyvotylivka–Vovchans’k and Gordineşti) side by side with potential Pre-Yamnaya and Yamnaya sources makes thus sense:

Now, the main obvious difference between Khvalynsk-Yamnaya and Corded Ware is the long-lasting, pervasive Y-chromosome bottlenecks under R1b lineages in the former, compared to the haplogroup variability and late bottleneck under R1a-M417 in the latter, which speaks in favour – on top of everything else – of a different community of sub-Neolithic hunter-gatherers including hg. R1a-M417 hijacking the expansion of Steppe_Maykop-related ancestry around the Volhynian-Podolian Upland.

Akin to how Yamnaya patrilineal descendants hijacked regional EEF (±CWC) ancestry components mainly through exogamy, dragging them into the different expanding Bell Beaker groups (see below), but kept their Indo-European languages, these hunter-gatherers that admixed with peoples of “Steppe ancestry” were the most likely vector of expansion of Uralic languages in Eastern Europe.

corded-ware-from-trypillia-maykop
PCA of ancient Eurasian samples. Marked likely Proto-Corded Ware samples and potential origin of its PCA cluster based on qpAdm results. See full PCA and more related files.

Baltic Corded Ware

One of the most interesting aspects of the results above is the surprising heterogeneity of the different regional groups, which is also reflected in the Y-DNA variability of early Corded Ware samples.

Seeing how Baltic CWC groups, especially the early Latvia_LN sample, show particularly bad fits with the models above, it seems necessary to test how this population might have come to be. My first impression in 2017 was that they could represent early Corded Ware groups admixed with Yamnaya settlers through their interactions along the Dnieper-Dniester corridor.

However, I recently predicted that the most likely admixture leading to their ancestry and PCA cluster would involve a Corded Ware-like group and a group related to sub-Neolithic cultures of eastern Europe, whose best proxy to date are EHG-like Khvalynsk samples (i.e. excluding the outlier with Pre-Yamnaya ancestry, I0434):

corded-ware-pca-sub-neolithic-europe
Detail of the PCA of the Corded Ware expansion. See full PCA and more related files.

Late Corded Ware + Yamnaya vanguard

Relevant are also the mixtures of Corded Ware from Esperstedt, and particularly those of the sample I0104, which I have repeated many times in this blog I suspected to be influenced by vanguard Yamnaya settlers:

The infeasible models of CWC + Yamnaya_Kalmykia ± Hungary_Baden (see below for Bell Beakers) and the potential cluster formed with other samples from the Baltic suggest that it could represent a more complex set of mixtures with sub-Neolithic populations. On the other hand, its location in Germany, late date (ca. 2500 BC or later), and position in the PCA, together with the good fits obtained for Germany_Beaker as a source, suggest that the increase in Steppe-related ancestry + EEF makes it impossible for the model (as I set it) to directly include Yamnaya_Kalmykia, despite this excess Steppe-related ancestry actually coming from Yamnaya vanguard groups.

I think it is very likely that the future publication of EEF-admixed Yamnaya_Hungary samples (or maybe even Yamnaya vanguard samples) will improve the fits of this model.

These results confirm at least the need to distrust the common interpretation of mixtures including late Corded Ware samples from Esperstedt (giving rise to the “up to 75% Yamnaya ancestry of CWC” in the 2015 papers) as representative of the Corded Ware culture as a whole, and to keep always in mind that an admixture of European BA groups including Corded Ware Esperstedt as a source also includes East BBC-like ancestry, unless proven otherwise.

yamnaya-vanguard-corded-ware-chalcolithic-early
Yamnaya vanguard groups in Corded Ware territory before the expansion of Bell Beakers (ca. 2500 BC). See full map.

Bell Beaker expansion

A hotly (re)debated topic in the past 6 months or so, and for all the wrong reasons, is the origin of the Bell Beaker folk. Archaeology, linguistics, and different Y-chromosome bottlenecks clearly indicate that Bell Beakers were at the origin of the North-West Indo-European expansion in Europe, while the survival of Corded Ware-related groups in north-eastern Europe is clearly related to the expansion of Uralic languages.

NOTE. For the interesting case of Proto-Indo-Iranians expanding with Corded Ware-like ancestry, see more on the formation of Sintashta-Potapovka-Filatovka from East Uralic-speaking Abashevo and Pre-Proto-Indo-Iranian-speaking Poltavka herders. See also more on R1a in Indo-Iranians and on the social complexity of Sintashta.

Nevertheless, every single discarded theory out there seems to keep coming back to life from time to time, and a new wave of interest in “Bell Beaker from the Single Grave culture” somehow got revived in the process, too, because this obsession – unlike the “Bell Beakers from Iberia Chalcolithic” – is apparently acceptable in certain circles, for some reason.

We know that Iberian Beakers, British Beakers, or Sicilian EBA – representing the most likely closest source population of speakers of Proto-Galaico-Lusitanian, Pre-Celtic Indo-European, and Proto-Elymian, respectively – have already been successfully tested for a direct origin among Western European Beakers in Olalde et al. (2018), Olalde et al. (2019), and Fernandes et al. (2019).

This success in ascertaining a closer Beaker source is probably due to the physical isolation of the specific groups (related to Germany_Beaker, Netherlands_Beaker, and NE_Mediterranean_Beaker samples, respectively) after their migration into regions dominated by peoples without Steppe-related ancestry. Furthermore, Celtic-speaking populations expanding with Urnfield south of the Pyrenees also show a good fit with a source close to France_Beaker.

So I decided to test sampled Bell Beaker populations, to see if it could shed light to the most likely source population of individual Beaker groups and the direction of migration within Central Europe, i.e. roughly eastwards or westwards. As it was to be expected for closely related populations (see the relevant discussion here), an attempt to offer a simplistic analysis of direction based on formal stats does not make any sense, because most of the alternative hypotheses cannot be rejected:

Not only because of the similar values obtained, but because it is absurd to take p-values as a measure of anything, especially when most of these conflicting groups with slightly ‘better’ or ‘worse’ p-values represent multiple different mixtures of the type (Yamnaya + EEF) + (Corded Ware + EEF ± Yamnaya), impossible to distinguish without selecting proper, direct ancestral populations…

A further example of how explosive the Bell Beaker expansion was into different territories, and of their extensive local admixture, is shown by the unsuccessful attempt by Olalde et al. (2018) to obtain an origin of the EEF source for all Beaker groups (excluding Iberian Beakers):

bell-beaker-local-population-iberia
Investigating the genetic makeup of Beaker-complex-associated individuals. Testing different populations as a source for the Neolithic ancestry component in Beaker-complex-associated individuals. The table shows P values (* indicates values > 0.05) for the fit of the model: ‘Steppe_EBA + Neolithic/Copper Age’ source population.
burials-yamnaya-hungary
Map of attested Yamnaya pit-grave burials in the Hungarian plains; superimposed in shades of blue are common areas covered by floods before the extensive controls imposed in the 19th century; in orange, cumulative thickness of sand, unfavourable loamy sand layer. Marked are settlements/findings of Boleráz (ca. 3500 BC on), Baden (until ca. 2800 BC), Kostolac (precise dates unknown), and Yamna kurgans (from ca. 3100/3000 BC on).

Now, there is a simpler way to understand what kind of Steppe-related ancestry is proper of Bell Beakers. I tested two simple models for some Beaker groups: Yamnaya + Hungary Baden vs. Corded Ware + GAC Poland. After all, the Bell Beaker folk should prefer a source more closely related to either Yamnaya Hungary or Central European Corded Ware:

Interestingly, models including Yamnaya + Baden show good fits for the most important groups related to North-West Indo-Europeans, including Bell Beakers from Germany, the Netherlands, Italy, and Poland, representing the most likely closest source populations of speakers of Pre-Proto-Celtic, Pre-Proto-Germanic, Proto-Italo-Venetic, and Pre-Proto-Balto-Slavic, respectively.

The admixed Yamnaya samples from Hungary that will hopefully be published soon by the Jena Lab will most likely further improve these fits, especially in combination with intermediate Chalcolithic populations of the Middle and Upper Danube and its tributaries, to a point where there will be an absolute chronological and geographical genomic trail from the fully Yamnaya-like Yamnaya settlers from Hungary to all North-West Indo-European-speaking groups of the Early Bronze Age.

The only difference between groups will be the gradual admixture events of their source Beaker group with local populations on their expansion paths, including peoples of mainly EEF, CWC+EEF, or CWC+EEF+Yamnaya related ancestry. There is ample evidence beyond ancestry models to support this, in particular continued Y-DNA bottlenecks under typical Yamnaya paternal lineages, mainly represented by R1b-L51 subclades.

east-bell-beaker-group-expansion
Distribution of the Bell Beaker East Group, with its regional provinces, as of c. 2400 cal BC (after Heyd et al. 2004, modified). See full maps.

European Early Bronze Age

European EBA groups that might show conflicting results due to multiple admixture events with Corded Ware-related populations are the Únětice culture and the Nordic Late Neolithic.

The results for Únětice groups seem to be in line with what is expected of a Central European EBA population derived from Bell Beakers admixed with surrounding poulations of East Bell Beaker and/or late (Epi-)Corded Ware descent.

Potential models of mixture for Nordic Late Neolithic samples – despite the bad fits due to the lack of direct ancestral CWC and BBC groups from Denmark – seem to be impossible to justify as derived exclusively from Single Grave or (even less) from Battle Axe peoples, supporting immigration waves of Bell Beakers from the south and further admixture events with local groups through maritime domination.

PCA of ancient European samples. Marked are Bronze Age clusters. See full PCAs.

Balkans Bronze Age

The potential origin of the typical Corded Ware Steppe-related ancestry in the social upheaval and population movements of the Dnieper-Dniester forest-steppe corridor during the 4th millennium BC raises the question: how much do Balkan Bronze Age groups owe their ancestry to a population different than the spread of Pre-Yamnaya-like Suvorovo-Novodanilovka chieftains? Furthermore, which Bronze Age groups seem to be more likely derived exclusively from Pre-Yamnaya groups, and which are more likely to be derived from a mixture of Yamnaya and Pre-Yamnaya? Do the formal stats obtained correspond to the expected results for each group?

Since the expansion of hg. I2a-L699 (TMRCA ca. 5500 BC) need not be associated with Yamnaya, some of these values – together with the assessment of each individual archaeological culture – may question their origin in a Yamnaya-related expansion rather than in a Khvalynsk-related one.

NOTE. These are the last ones I was able to test yesterday, and I have not thought these models through, so feel free to propose other source and target groups. In particular, complex movements through the North Pontic area during the Late Eneolithic would suggest that there might have been different Steppe-ancestry-related vs. EEF-related interactions in the north-west and west Pontic area before and during the expansion of Yamnaya.

Mycenaeans

One of the key Indo-European populations that should be derived from Yamnaya to confirm the Steppe hypothesis, together with North-West Indo-Europeans, are Proto-Greeks, who will in turn improve our understanding of the preceding Palaeo-Balkan community. Unfortunately, we only have Mycenaean samples from the Aegean, with slight contributions of Steppe-related ancestry.

Still, analyses with potential source populations for this Steppe ancestry show that the Yamnaya outlier from Bulgaria is a good fit:

The comparison of all results makes it quite evident the why of the good fits from (Srubnaya-related) Bulgaria_MLBA I2163 or of Sintashta_MLBA relative to the only a priori reasonable Yamnaya and Catacomb sources: it is not about some hypothetical shared ancestor in Graeco-Aryan-speaking East Yamnaya– or even Catacomb-Poltavka-related groups, because all available Yamnaya-related peoples are almost indistinguishable from each other (at least with the sampling available today). These results reflect a sizeable contribution of similar EEF-related populations from around the Carpathians in both Steppe-related groups: Corded Ware and Yamnaya settlers from the Balkans.

mycenaeans-minyan-ware-greece-minoan
Cultural groups in and around the Balkans during the Early Bronze Age. See full maps.

qpAdm magic

In hobby ancestry magic, as in magic in general, it is not about getting dubious results out of thin air: misdirection is the key. A magician needs to draw the audience attention to ‘remarkable’ ancestry percentages coupled with ‘great’ (?) p-values that purportedly “prove” what the audience expects to see, distracting everyone from the true interesting aspects, like statistical design, the data used (and its shortcomings), other opposing models, a comparison of values, a proper interpretation…you name it.

I reckon – based on the examples above – that the following problems lie at the core of bad uses of qpAdm:

  1. In the formal aspect, the poor understanding of what p-values and other formal stats obtained actually mean, and – more importantly – what they don’t mean. The simplistic trend to accept results of a few analyses at face value is necessarily wrong, in so far as there is often no proper reasoning of what is being assessed and how, and there is never a previous opinion about what could be expected if the alternative hypotheses were true.
  2. In the interpretation aspect, the poor judgement of accompanying any results with simplistic, superficial, irrelevant, and often plainly wrong archaeological or linguistic data selected a posteriori; the inclusion of some racial or sociopolitical overtones in the mixture to set a propitious mood in the target audience; and a sort of ritualistic theatrics with the main theme of ‘winning’, that is best completed with ad hominems.

If you get rid of all this, the most reasonable interpretation of the output of a model proposed and tested should be similar to Nick Patterson’s words in his explanation of qpWave and qpAdm use:

Here we see that, at least in this analysis there are reasonable models with CordedWareNeolithic is a mix of either WHG or LBKNeolithic and YamnayaEBA. (…) The point of this note is not to give a serious phylogenetic analysis but the results here certainly support a major Steppe contribution to the Corded Ware population, which is entirely concordant with the archaeology [?].

Very far, as you can see, from the childish “Eureka! I proved the source!”-kind of thinking common among hobbyists.

The Mycenaean case is an illustrative example: if the Yamnaya outlier from Bulgaria were not available, and if one were not careful when designing and assessing those mixture models, the interpretation would range from erroneous (viz. a Graeco-Aryan substrate, as I initially thought) to impossible (say, inventing migration waves of Sintashta or Srubnaya peoples into Crete). The models presented above show that a contribution of Yamnaya to Mycenaeans couldn’t be rejected, and this alone should have been enough to accept Yamnaya as the most likely source population of “Steppe ancestry” in Proto-Greeks, pending intermediate samples from the Balkans. In other words, one could actually find that ‘the best’ p-values for source populations of Mycenaeans is a combination of modern Poles + Turks, despite the impracticality of such a model…

I haven’t been able to reproduce results which supposedly showed that Corded Ware is more likely to be derived from (Pre-)Yamnaya than other source population, or that Corded Ware is better suited as the ancestral population of Bell Beakers. The analyses above show values in line with what has been published in recent scientific papers, and what should be expected based on linguistics and archaeology. So I’ll go out on a limb here and say that it’s only through a careful selection of outgroups and samples tested, and of as few compared models as possible, that you could eventually get this kind of results and interpretation, if at all.

Whether that kind of special care for outgroups and samples is about (a) an acceptable fine-tuning of the analyses, (b) a simplistic selection dragged from the first papers published and applied indiscriminately to all models, or (c) cherry picking analyses until results fit the expected outcome, is a question that will become mostly irrelevant when future publications continue to support an origin of the expansion of ancient Indo-European languages in Khvalynsk- and Yamnaya-related migrations.

Feel free to suggest (reasonable) modifications to correct some of these models in the comments. Also, be sure to check out other values such as proportions, SD or SNPs of the different results that I might have not taken into account when assessing ‘good’ or ‘bad’ fits.

Related

European hydrotoponymy (V): Etruscans and Rhaetians after Italic peoples

italy-mediterranean-bronze-age

There is overwhelming evidence that the oldest hydrotoponymic layer in Italy (and especially Etruria) is of Old European nature, which means that non-Indo-European-speaking (or, at least, non-Old-European-speaking) Etruscans came later to the Apennine Peninsula.

Furthermore, there is direct and indirect linguistic, archaeological, and palaeogenomic data supporting that the intrusive Tursānoi came from the Aegean during the Late Bronze Age, possibly through the Adriatic, and that their languages spread to Etruria and probably also to the eastern Alps.

Hydrotoponymic layer

The following are translated excerpts (emphasis mine) from Lenguas, genes y culturas en la Prehistoria de Europa y Asia suroccidental, by Villar et al. Universidad de Salamanca (2007):

villar-vascos
Lenguas, genes y culturas en la Prehistoria de Europa y Asia suroccidental (2007). Buy the ebook online (or the printed version, if available).

‘(Indo-)Mediterranean’ substrate?

The name Indo-Mediterranean substrate was spread in Italy by the work of V. Pisani. Other Italian scholars continued this idea, such as W. Belardi, L. Heilmann, D. Silvestri, etc. In their hands, the nuclear area of ​​the Indo-Mediterranean substratum was established as follows: “il mondo culturale indomediterráneo trova i suoi più importanti centri di gravitazione (e, soltanto secondariamente, di espansione) nel Mediterràneo Orientale (Creta, Cipro, Asia Minore), nella ‘regione dei due fiumi’ (area di espansione subarea) e nella valle dell’Indo (civiltà de Harappa e Mohenjo Daro)”. From there they could have spread to other areas, such as the western Mediterranean. Even at one point there was talk of “a Mediterranean oasis in the Baltic”, whose main basis was the existence of numerous lexical elements, real or supposedly pre-Indo-European in the Baltic languages.

One of the paradoxes of the theory of the Mediterranean substrate is that the lexical or toponymic components that are attributed to it can rarely be explained etymologically from the surviving languages ​​of said supposed substrate; sometimes they are not even very compatible with what we know of the non-Indo-European languages ​​of the corresponding area. For example, neither Basque nor Iberian have an ancestral and autochthonous phoneme /p/, while that phoneme is frequent in substrate words (cf. among the few mentioned above *pal- and *lap-). In fact, for these three languages ​​other alternative origins have been imagined, so that they would not be representatives of the local substrate: Basque (North Africa, the Caucasus), Iberian (North Africa), Etruscan (Asia Minor). Thus, under such hypotheses the non-Indo-European languages ​​attested in Italy and the Iberian Peninsula would not be autochthonous, but as immigrant as the Indo-European languages.

akwa-hydronyms
Akʷa hydronyms. The majority of old serial elements are found in Italy, with 9, where they don’t appear as second element. Different to the southern areas, they are found in especially frequent compounds in the acha-Namen in Germany, and hyper-represented (as usual) in Lithuania, which shows strictly 8 ancient names.

Italy and Iberia

Let’s review data on Italy:

I. Serial tponyms and hydronyms of Italy:

  1. ub-: Caecubus, Egubium, Litubium, Marrubium, Olobia, Rutuba, Tardoba, Tardubius, Verubius, etc.
  2. uc-: Aluca, Arucia, Arugus, Ausucum, Ausugum, Motuca, Uccia.
  3. ur-: Orinos, Stura, Stura, Astura, Tibur, Caburrum, Calorem.
  4. urc-: Coturga, Orgus, Urcia, Urcinia, Urgo.
  5. bai-: Baebiani.
  6. tuc-: Tucianus (pagus).
  7. murc-: Murcia, Murgantia, Murgantia.
  8. *war: Varduli, Barduli.

ub-hydronyms

II. Non-serial toponyms and hydronyms of Italy: Aesis, Aisis, Ana, Ania, Anios, Arsia, Astura, Ausa, Ausonia, Ausculum, Bardinisca vallis, Barduli, Basentius, Basta, Boron, Cabienses (Cabia), Caburrum, Cales, Cales, Casta Ballenis, Ceresium, Cerili, Corsica, Cortona, Curicum, Ispelum, Ispila, Isporos, Istonium, Istria, lacus, Latis, Latium, Laurentum, Laurentes, Luca, Lucania, Lucera, Maleventum, mare, Marrucini, Minio, Minius, Oscela, Osci, Ossa, Ostia, Paestum, Pisaurum, Pisaurus, Sabini, Sagis, Savo, Sila, Silarus, Silis, Soletum, etc.

italy-iberia-hydronymy-toponymy

Not few of the coincident place names between the southern Iberian and Italic material are rigorous cognates. We understand by such the names that not only coincide in the root or in the serial element, but in the whole root set plus suffixes, or – if it is a compound – in the two sets of roots plus suffixes. In addition to the ones that we are going to present below, there are others that we did not mention because the Iberian correlate was not found within the southern group, but in other geographical areas, as is the case, for example, with the Italian Mantua and the Spanish Mantua (Carpetania).

As can be seen, the parallels between the southern Iberian toponymic area and the Italic one are so wide and strict that the mere calculation of probabilities makes any attempt to attribute them to the mere chance of random homophony irrational. And the improbability of chance increases as coincidences are added in new places in Europe. What will not prevent, for sure, that some would resort to it as an explanation, in particular those who are reluctant to abandon the conception of the prehistory of the European continent that underlies their usual approaches, which suffer an irreparable strike when they are confronted with these data.

The second aspect, the compatibility of this material with Indo-European etymology, offers another significant correlation: the “southern” series that are also found in the Ibero-Pyrenean region and in Italy (and the rest of western Europe) are compatible with Indo-European etymologies; (…)

I will spare the reader of all proposed Indo-European etymologies, most of which are fairly evident. Those interested should buy one of the books, or both.

or-hydronyms

Etruria

(…) in the whole of Italy there is a considerable collection of toponyms and hydronyms of “Southern Iberian” type, whose joint inventory we have contributed to above. From them we find in Etruria Ause, Veturris / Bituriza, Castola, Hasta, Cortona, Luca, Minio, Osa / Ossa, Pissai, Pistoria. The Hispanic and Italian correlates of those names are:

iberian-etruscan-indo-european

However, the inventory of ancient names and hydronyms of Etruria compatible without discussion with well-known Indo-European etymologies is much wider: Albina, Alma, Alsium, Arnine, Arnos, Arnus, Aventia, Marta, Pallia, Umbro, Vetulonium, Volsinii. Furthermore, the majority of Etrurian hydronyms have non-Latin Indo-European etymology: Albina, Alma, Arnine, Arnos, Arnus, Auser, Aventia, Marta, Minio, Osa, Ossa, Pallia, Umbro. And very few of the others (Clusinus, Cremera, Lingeus, Trasumenus, Vesidia) could claim an Etruscan etymology, if only one could do so.

In summary, the territory occupied by Etruscans presents a hydro-toponymic situation very similar to that of the rest of Italy and Western Europe: it exhibits a very deep toponymic stratum of Indo-European character to which most hydronyms attested in antiquity belong. As we know the history of Etruria from the end of the 1st millennium BC, and we know that no other Indo-European peoples mediated between the Etruscans and the Romanization of the territory, we must conclude that this ancient toponymy was there before the Etruscans arrived or emerged in that place. And, when the Etruscans settled there, they did not have the opportunity to put names of their language to the rivers in general, because they had already received them from a previous people and the Etruscans limited themselves to learning them, adapting them to their language, and transmitting them in turn to the Romans. When the latter Romanized Etruria, they limited themselves to incorporating those names and adapting them to Latin.

maro-maranto

Etruscans

The ‘foreign’ Tyrsenians

Here is a recapitulation of the main reasons why Etruscans were recently intrusive to Italy, as they appeared in The Origin of the Etruscans, by Beekes (2003):

NOTE. You can read another version of the text in PDF, as the main paper from Biblioteca Orientalis LIX(3-4) 2002.

  1. The tradition as given by Herodotus and Dionysius of Halikarnassos.
  2. The story that the Etruscans were Pelasgians.
  3. The use of the term ‘Tyrsēnoi’ for both Etruscans and a people in north-western Asia Minor. Above we argued that the eastern Tyrsēnoi are the remnant of a population. This means that the Tyrsēnoi/Etruscans came from this area.
  4. The Lemnos inscription.
  5. To the testimony of Lemnos must now be added that Herodotus says that the people of Plakiê and Skylakê spoke the same language as the Etruscans.
  6. etruscan-homeland

  7. The kumdanlı inscription. (…) lake Egridir (of which the old name is unknown, unless it was just Limnai). This is just over the border of classical Lydia. The inscription dates from the second century ad and is given by Ramsay (i883); the same inscription is cited by Sundwall (i9i3, 22i). It mentions three people as Tyrsēnoi(67, 68, i02). Though very late, the inscription is of great interest, as it is the only time that we have inscriptional evidence for Tyrsēnoi in Asia Minor. (And nobody will argue that these were Etruscans from Italy.) (…)
  8. The suffix -ānos. The suffix -ānos in the name Tyrsēnoi (with ē from ā) points to the north-west of Asia Minor. It has long since been recognized that this suffix for ethnic names is at home in north-west Asia Minor; some think that it is of non-Greek origin; cf. Αβυδηνός , Ολυμπιηνός, Περγαμηνός, Σαρδηνός; (see Chantraine i933, 206; Schwyzer 490 (6); De Simone i993, 88ff.). This proves that the name Tyrsēnoi originated in the north-west of Asia Minor. (…)
  9. Loanwords. As to the language, Steinbauer (i999, 367) observes that Etruscan shows most connections (loanwords) with Lydian (…)
  10. Tarchon. The definite proof of the oriental origin of the Etruscans is that a ‘hero’ of great significance is Tarchon (Briquel i99i). He is clearly the Stormgod Tarhun(t)-, the highest god of the Luwians and Hittites.
  11. Nanas. This identification is strongly confirmed by the story that the Etruscans were Pelasgians who came from Greece under Nanas (Nanos), mentioned by Hellanikos. This name was long ago recognized as an Anatolian ‘Lallname’.
  12. The triumphus complex. In his study of the Roman triumphus Versnel has shown that (i970, 293): ‘the Etruscans brought the New Year festival with them from Asia Minor, together with the god who formed the centre of it, a god whom the Greeks called Dionysos, the Etruscans Tinia (or by an Italic name Voltumna), a figure of the ‘dying and rising’ type, who was invoked by the cry *thriambe and who on New Year’s Day was represented by the king.’ And on p. 300: ‘The Etruscans brought the New Year festival with them from Asia Minor and gave Rome two ceremonies: the ludi Romani as the festival of the New Year, the triumph as the festival of the victory. … Only along this way is it possible to explain the data: i. the Dionysiac call to epiphany triumpe, introduced via Etruria; 2. the identification of the Roman victorious general and of the magistrate leading the games with the god Iuppiter; 3. the typological and historic relation between the ludi Romani and the triumph.’
  13. The double axe. On a smaller issue Versnel concludes (p. 299): ‘When this bipennis [‘double axe’], property of ‘Zeus Bakchos’, carried as symbol of sacred power by Lydian kings, is encountered again as the symbol of the royal authority of the Etruscan kings, particularly of the supreme king of the federation of cities, this may be considered an important indication of the Asia Minor origin of the entire underlying ideology, and of the ceremony of investiture in which the bipennis played a part.’ These conclusions are of primary importance, as they concern a deeprooted complex of religious views that cannot have been taken over from elsewhere.
  14. The Kabeiroi. One might also recall the Latin word camillus, which means a young boy of noble birth who assists with ritual actions. (…) Probably more evidence can be found in the field of religion, such as the much discussed hepatoscopy. It seems quite probable to me that the lituus, the crosier used by the Roman priests, is Anatolian (see e.g. Wainwright i959, 2i0; cf. Haas i99i, Abb. 75, the Stormgod standing on an animal with his lituus over his shoulder).
  15. The Etruscan way of life. There was in antiquity much criticism on Etruscan customs, concerning cruelty, sexual behaviour, and the behaviour of women. (…) Dionysius concluded from the fact that they were so strange that they had always lived in Italy, whereas it is of course much more natural to explain it by assuming that they were strangers.
  16. No withdrawal area. We have seen above that Tuscany is not a ‘withdrawal area’, where an ancient people may hold out when the country is invaded. On the contrary, it is a desirable area which the Indo-European peoples, had they come later, would certainly have occupied. (But it went the other way: the Etruscans came long after the Indo-Europeans and settled there/conquered the country.)
  17. sea-peoples-expansion-tyrsenians
    The Sea Peoples in the Eastern Mediterranean c. 1200 BC. Map by Ian Mladjov.
  18. Archaeology. Many scholars would like to see archaeological evidence, but I think that it is quite possible that we shall never find any.
  19. The 1200 crisis. In 1200 the whole Mediterranean was in commotion; the Mycenaean and Hittite worlds, between which the TyrseOEnoi lived, disappeared. So the movement of the Etruscans fits very well in the general picture. That this was the setting of the migration of the Etruscans has been assumed by many earlier scholars.
  20. The ten saecula. As to the time, it has been argued that the Etruscans thought that their world would last ten saecula (Briquel i999, 58; Pfiffig i975, i59ff.). The way of counting provides several problems, however (…) If we accept it, we arrive at 968 bc. Now we do not know from when one started counting. This might have been a decisive victory over the Umbrians, or a kind of unification of the Etruscans, or the founding of an important city. It could well be that this was some 200 years after the arrival of the Etruscans, which would take us to 1168 bc. (…)
  21. The famine. Herodotus states that the reason for the departure of the Tyrsēnoi was a long famine. This has been identified as the famine about i200. (…)
  22. The sea-peoples. (…) The phenomenon as a whole stands, it seems; the problem is the details: which peoples took part in which movements? In our case, as the Lukka are mentioned (which were very probably the Lycians), the Tyrsēnoi may have been involved as well. So the question is whether the T(w)r(w)š, mentioned by Merneptah, were the Tyrsēnoi. We have no confirmation, but it seems quite possible.
  23. The journey. We know from the abundant finds of ceramics in the i3th century that the Mycenaeans knew the sea-route to Italy. (…)
  24. The Umbrians. Pliny (3, ii2) states that the Etruscans conquered 300 cities from the Umbrians (Trecenta eorum oppida Tusci debellasse reperiuntur.). This clearly refers to the ‘Landnahme’. This statement is confirmed by the river Umbro (mod. Ombrone), which flows in its full length in Etruscan territory. The river will have given its name to the people, or vice versa. Anyhow, the river will have flowed in Umbrian territory; so the Etruscans must have pushed the Umbrians out.
  25. The name Sergestus, of a prominent friend of Aeneas, seems identical with Lydian Srkastu- and Phrygian Surkastos (…) it is excluded that (Virgil) got it from Lydia or Phrygia, or Asa Minor in general. So he must have got it at home, from a source that was acqainted with Etruscan traditions. This means that the name was known to the Etruscans (or those who studied their traditions). Above I proposed that it lives on in Etr. Sekst-alu-.

You can read the full text (and its appendices) for further evidences adduced by Beekes, who considers the matter mostly settled.

Local Italic peoples

Another main reason for the intrusion of Tyrsenians among local groups is the ancient connection between Italic languages, which most likely formed an ancient Apennine dialect continuum:

  • the core Italic group with Latino-Faliscan and Palaeo-Sabellic – probably also including an Ausonian-Siculian branch – separated ca. 1500-1000 BC;
  • NOTE. Sicel is believed to have arrived in Sicily with Ausonian-Siculian speakers either around the 13th c. or in the middle of the 11th c. BC (or in both waves), from their ancient settlements in the mainland, driving prior inhabitants (Elymians) to the east of the island, which sets another clear terminus ante quem for the expansion of Italic languages in southern Italy.

  • and the possibly more distantly related North Picene and Venetic, connecting all roughly to an early to mid-2nd millennium BC language.

This continuum was probably broken (with language replacement and displacement events) with the 12th c. BC turmoil and the emergence of new social hierarchies. The adoption of older place and river names, as well as the lack of long-lasting influence on neighbouring languages, suggests that the predominance of the Etruscan language in its proto-historic territory was probably gradual and quite recent.

NOTE. For more on guesstimates, relative chronological expansions and potential archaeological identifications, see e.g. “Ausgliederung und Aufgliederung der italischen Sprachen”, by Helmut Rix In: Languages in Prehistoric Europe (2003). Or, basically, any recent (linguistic) text on the distribution and attribution of ancient Apennine languages to the Ital(o-Venet)ic group.

Italic-venetic-etruscan-languages-map
Languages of pre-Roman Italy and nearby islands. Italo-Venetic languages surrounded with shadowed red border. I1, South Picene; I2, Umbrian; I3, Sabine; I4, Faliscan; I5, Latin; I6, Volscian and Hernican; I7, Central Italic (Marsian, Aequian, Paeligni, Marrucinian, Vestinian); I8, Oscan, Sidicini, Pre-Samnite; I9, Sicel; IE1, Venetic; IE2, North Picene; IE3, Ligurian; IE4, Elymian; IE5, Messapian; C1, Lepontic; C2, Gaulish; G1-G2-G3, Greek dialects (G1: Ionic, G2: Aeolic, G3: Doric); P1, Punic; N1, Rhaetian; N2, Etruscan; N3, Nuragic. Image modified from Davius Sanctex.

Archaeology

The main criticism against this ethnolinguistic model of foreign Tyrsenians comes, surprisingly, from the lack of archaeological data to support this arrival. Or, rather, fitting anthropological interpretations of a culture of Asia Minor with similar hierarchical societies (?). From Review of R. S. P. Beekes, The Origin of the Etruscans, by Mahoney, Etruscan Studies (2008):

A crucial part of Beekes’ argument, however, is that there is a significant cultural break in Etruria around 1200, at the beginning of the Early Bronze Age or Proto-Villanovan period (p. 34, citing Briquel and Torelli). The introduction of cremation can be dated to around this period, and there is also evidence for a new hierarchical social organization (convenient summary in Barker and Rasmussen, p. 53-60). Beekes simply says that there is a change, and changes of this sort can come about when new people move in to an existing society, so therefore this change is consistent with his theory. That is correct as far as it goes, but what is missing is any consideration of how and why people coming in from Asia Minor would cause the particular changes that take place in Etruria. Can we argue that the society of the pre-migration Tyrsenians was hierarchical in the same way as those of the various Indo-European-speaking peoples in the region? Beekes simply says “what we still would like to have is material objects, or art traditions etc., from Etruria agreeing with their homeland” (p. 34). What we would really like to have is evidence for the organization of society in this alleged homeland.

Weird as this criticism is, here it is yet another example of the social change brought about under Eastern Mediterranean influences during the Final Bronze Age, from a recent paper (behind paywall) Mobile elites at Frattesina: flows of people in a Late Bronze Age ‘port of trade’ in northern Italy, by Cavazzuti et al. Antiquity (2019):

Introduction

The collapse of the Terramare system c. 1150 BC was followed by a sudden and substantial depopulation of the central part of the Po Plain (Cardarelli 2009). At the beginning of the Final Bronze Age, the southern part of the Po Valley was almost abandoned. In contrast, in the northern part of the Po Valley, some villages survived (…) Concurrently, a new territorial system arose, pivoting around the socio-economic pole of Frattesina (Calzavara Capuis et al. 1984; Bietti Sestieri et al. 2015; Cupitò et al. 2015). Therefore, within the area of the wider Terramare ‘culture’, local responses to the crisis led to different outcomes, some of which were relatively successful and others catastrophic. Economic factors—both in terms of internal carrying capacity and degree of openness to external relations—probably played a key role in determining different responses to the tensions.

The communities of the Terramare, especially in the southern area, were probably not flexible enough to adapt their political structure and modes of production to the needs of a rapidly changing world. Moreover, the domino effect from the overall geo-political instability of the twelfth century BC, in a highly interconnected system such as the Mediterranean, was undoubtedly another factor (Cardarelli 2009). The lack of evidence in the southern Terramare area for connections with the Aegean and the Levant suggests a more ‘closed’ system located on the edge of the ‘globalised’ world of the Late Bronze Age. In contrast, there is well-documented evidence from the largest terramare on the northern side of the Po River for possible incipient institutionalised, well-connected elites—particularly at Fondo Paviani, which has yielded locally produced pottery in Levantine and Late Helladic IIIC Aegean-Mycenaean styles (Bettelli et al. 2015).

The display of austere equality that dominated the Middle and Late Bronze Age ‘urnfields’ (Salzani 2005; Cardarelli 2014) strongly limited funerary expressions of social differentiation. Internal inequalities nonetheless existed between different co-resident extended families and lineages comprising tens of individuals at most (e.g. at Casinalbo; Cardarelli et al. 2014: 722–28), and, above all, between large centres, such as the terramara at Fondo Paviani and dependent satellite settlements (Balista et al. 2005; Cupitò et al. 2015). It seems reasonable therefore to hypothesise that groups based at nodal sites in the system attracted more prestige goods from exotic places, along with individuals from distant areas, while small villages attracted people mainly from within a local radius (Cavazzuti et al. 2019a). Within this dynamic cultural context, the Final Bronze Age funerary evidence from Frattesina documents a more elaborate display of power and wealth concentrated in the hands of elites. At Le Narde (Frattesina’s cemetery), this privileged segment of society, probably with its own entourage, is clearly represented by a small number of burials with several indicators of prestige.

bow-fibula-italy-aegean
Distribution of the violin-bow fibulae with two temple knots in the different source categories. Map by Sabine Pabst (2018).

Results

(…) the individual in burial Narde1-168 may have achieved the status of a ‘warrior-chief’, as symbolised by the presence of an Allerona-type sword (Bianco Peroni 1970). This was ritually broken and deposited in pieces inside the grave, along with a bronze pin, a pair of tweezers and other ornaments (Figure 8). (…) yielded a strontium isotope ratio (0.70983) that is incompatible with the local 0–20km baseline (Table 3), but fits within the 20–50km range. By contrast, the value obtained from the femoral cortical bone (0.70924) is consistent with the local range of Frattesina. This means that this individual moved to the site after early childhood—possibly during youth or early adulthood—and he probably spent the last years of his life there, at the apex of the community.

Marshall Sahlins (1981), in his famous article ‘The stranger-king: or Dumézil among the Fijians’, compares the dynamics of power in the Fiji Islands to the Indo-European tradition, arguing that human societies tend to locate power as originating from the outside (Sahlins 1981, 2008; see also Ling & Rowlands 2015). Sahlins focuses on origin myths across ancient polities in the Indo-European language area, which systematically feature a dichotomy between what the Romans called gravitas and celeritas. The former refers to the conservative, peaceful and productive character of an established native community, while celeritas represents the disruptive, transformative violence personified in the stranger king, who “erupts upon a pastoral scene of peaceful husbandry and political equality (or at least limited authority)” (Sahlins 1981: 112).

grave-goods-frattesina-warrior-chief-allerona-sword
The grave goods and cremated bones of burial Narde1-168 (after Salzani 1989). Urn height is 0.26m, sword length is 0.46m.

The individual buried in grave Narde1-168 at Frattesina was probably neither a true ‘king’, nor a true ‘stranger’. Despite its uniqueness, his grave resembles those of the rest of the community and is included within a large collective—or at least not evidently exclusive—burial mound. ‘Warrior-chief’ perhaps would be a more appropriate definition for this individual. Moreover, his place of origin was not so distant as to define him as a ‘stranger’. Nonetheless, Sahlins’s archetype of the ‘stranger-king’ evokes the power of alterity; burial Narde1-168 perfectly embodies celeritas, which breaks with the gravitas of the former Terramare tradition and guided whatever survived the collapse towards a new social model. Since the discovery of Frattesina and its cemeteries, Italian scholars have debated the mechanisms underlying the origin and economic success of the settlement, and the degree of foreign (i.e. Cypriot and Levantine) involvement in this process as suggested by archaeological finds (Cupitò et al. 2015). The new isotopic data presented here demonstrate that even though some individuals may have come from the Levant—where the available 87Sr/86Sr baseline ranges from 0.7079–0.7086 (Sheridan & Gregoricka 2015; Gregoricka & Guise Sheridan 2016)—or were from other exotic places, they nonetheless represent a minority of the population and, in any case, not the upper elite. The latter appear quite mobile, although probably from within the broader hinterland radius.

Adriatic or western route?

One of the interesting questions, and probably non-trivial for the correct interpretation of ancestry in future ancient DNA samples, is from where exactly did Tyrsenians come from, and more importantly where exactly did the arrive, and when. I have the impression that a Tyrrhenian Sea route is more commonly depicted (as in my maps) due to the historical predominance of Etruscans in the west, but that archaeologists usually consider the Adriatic – and thus a spread from the Po River Valley and/or Pannonia – a more likely route for Tyrsenian speakers, and probably rightly so.

NOTE. The tentative (and highly speculative) classification of fragmentary Rhaetian as more archaic than Etruscan relative to Lemnian may give further support to this route.

Failing a precise time transect from a population geographically close to the origin of their expansion in central or northern Italy, we are bound to see the same misinterpretations of the data we have seen in the case of Sea Peoples of hg. R1b behind Philistines. Nevertheless, here are some interesting predictions of population movements by Pabst (2013) based particularly on the Stätzling-/Allerona-sword from Narde in Veneto, which have been confirmed for the moment with isotope analyses, showing that some peoples of Frattesina had previously lived in the eastern Mediterranean, and that local elites had a much closer origin:

staetzling-swords
Distribution of the various blade profiles of the Stätzling (l) and Casale (H) type of leaf blade sword: 1 White symbols: blade with rapier-like ribs. – 2 black symbols: flat rhombic blade profile.- 3 Large gray symbols: a blade with a narrow midrib and longitudinal grooves.- Small gray symbols: lenticular or indefinite blade profile. (Map S. Pabst).

An Ingot fragment from the hoard of Hočko Pohorje in Styria, Slovenia indicates that possibly also Pannonia was involved in the 12th century BC (or during stage Ha A1) in the East and Central Mediterranean copper trade. According to the chemical composition or the high iron content, it is particularly close to individual finds from Sardinia, Italy and Anatolia.

The people behind the Stätzling swords could have been the potentates of this supraregional trade in the Adriatic and Ionian seas. This is also to be expected from the presence of late Mycenaean populations on the upper Adriatic. This is indicated – in addition to individual Mycenaean ceramics imports – especially in the Aegean Stätzling sword from the fly cave of Škocjan in the hinterland of Trieste, in this exchange network of the 12th century BC. However, not only people from the late Mycenaean cultural area were involved in the process. For native elites are suspected behind the mostly locally manufactured Stätzling swords in Pannonia and Italy, according to the analysis of the grave find 227 of Narde; perhaps local organizers of the trade, as allies of the Mycenaean chiefs.

Palaeogenomics

Palaeogenomics might help shed light upon the complex matter of the Tyrsenian emergence in Europe. Even though Rhaetian is a fragmentary language, it seems that it is related to Etruscan, and neither are remnant languages from the Bronze Age, but rather intrusive languages to Italy and Central Europe.

It is more than likely, then, that ancient DNA will show an increase in Aegean ancestry during the Late/Final Bronze Age in central and/or northern Italy, even if this change is found rapidly diluted within generations, as happened with the Aegean ancestry among Philistines, who – in spite of this dilution – also left their prolonged linguistic mark in the Levant.

This is the summary I made of an online report from oral communication A 12,000-year Genetic History of Rome and the Italian Peninsula, by Hannah Moots, the 6th February 2019, with 134 samples from Lazio and surrounding areas:

bronze-iron-age-romans-etruscans-osco-umbrians-map
Bronze Age – Iron Age evolution of Italy Top Left: Early Bronze Age cultures. Bottom left: PCA of groups from the Bronze Age; marked in red are previous Italy Bell Beakers. Top Right: Early Iron Age cultures. Bottom right: PCA of groups from the Iron Age – Middle Ages; marked in red are the approximate location of described ancient Italian clusters, one including Etruscans, Osco-Umbrians, Picentes, etc., and the wider cluster of Romans (dates unknown). See full maps and PCAs.

While Bronze Age samples of west-central Italy show a clear homogenisation of the genetic pool, with a shift in the PCA towards central Europe (away from the previous CHG/Iran Neolithic influence), and thus close to the modern Sardinian cluster, the few investigated Iron Age samples from the Republican period (ca. 700–20 BC) show a widespread genetic cluster encompassing the modern Italian ones, overlapping North Italian (ca. 60%) or South Italian/Sicilian (ca. 40%) clusters. The arrival or increase of EHG-, Levant Neolithic-, or CHG/IN-related ancestry in samples from this period suggest influence from previous population movements during the LBA from the north or through the Mediterranean, respectively. The Imperial Period shows influence from CHG/IN-related ancestry, but only sporadically Levant Neolithic.

NOTE. For more on the referred northern and southern Italian clusters, see Population structure of modern-day Italians reveals patterns of ancient and archaic ancestries in Southern Europe, by Raveane et al. bioRxiv (2018).

italian-north-south-clusters
Principal component analysis projecting 63 ancient individuals onto the components inferred from modern individuals. A) Principal component analysis projecting 63 ancient individuals onto the components inferred from 3,282 modern individuals assigned, through a CP/fS analysis, to European West Asian and Caucasian clusters.

The alternative view

Kristiansen is among those who offer an alternative view in the archaeological question, supporting the opposite direction of population movements: of Terramare migrants in Greece, a theory which is not to be lightly dismissed, in the complex setting of population movements across the Mediterranean during the Final Bronze Age.

As a weak linguistic support for such a movement, one can find the hypothesis of Eteo-Cretans as Osco-Umbrian speakers, based on de Ligt’s speculative interpretation of the Praisos inscription (Talanta 2008-2009).

It seems that, even if these views are also correct, the overwhelming evidence is for a foreign origin of Tyrsenians:

  • lack of Tyrsenian hydrotoponymic layers in Italy or central Europe;
  • guesstimates and “split” distribution of Italo-Venetic languages;
  • foreign culture and influences of Etruscans;
  • recent predominance and influence of Etruscan language and culture;
  • East Mediterranean peoples in the LBA Po Valley (isotope analyses);
  • genetics of Sea Peoples from the Aegean.

Related

Villabruna cluster in Late Epigravettian Sicily supports South Italian corridor for R1b-V88

epipalaeolithic-whg-expansion

New preprint Late Upper Palaeolithic hunter-gatherers in the Central Mediterranean: new archaeological and genetic data from the Late Epigravettian burial Oriente C (Favignana, Sicily), by Catalano et al. bioRxiv (2019).

Interesting excerpts (emphasis mine):

Grotta d’Oriente is a small coastal cave located on the island of Favignana, the largest (~20 km2) of a group of small islands forming the Egadi Archipelago, ~5 km from the NW coast of Sicily.

The Oriente C funeral pit opens in the lower portion of layer 7, specifically sublayer 7D. Two radiocarbon dates on charcoal from the sublayers 7D (12149±65 uncal. BP) and 7E, 12132±80 uncal. BP are consistent with the associated Late Epigravettian lithic assemblages (Lo Vetro and Martini, 2012; Martini et al., 2012b) and refer the burial to a period between about 14200-13800 cal. BP, when Favignana was connected to the main island (Agnesi et al., 1993; Antonioli et al., 2002; Mannino et al. 2014).

sicily-grotta-oriente
A-B) Geographic location of Grotta d’Oriente.

The anatomical features of Oriente C are close to those of Late Upper Palaeolithic populations of the Mediterranean and show strong affinity with other Palaeolithic individuals of Sicily. As suggested by Henke (1989) and Fabbri (1995) the hunter-gatherer populations were morphologically rather uniform.

Genetic analysis

We confirmed the originally reported mitochondrial haplogroup assignment of U2’3’4’7’8’9. This haplogroup is present in both pre- and post-LGM populations, but is rare by the Mesolithic, when U5 dominates (Posth et al.2016).

Lipson et al. (2018) (their supplementary Figure S5.1) and Villalba-Mouco et al. (2019) (their Figure 2A) showed that European Late Palaeolithic and Mesolithic hunter-gatherers fall along two main axes of genetic variation. Multidimensional scaling (MDS) of f3-statistics shows that these axes form a “V” shape (Fig. 3). (…)

Focusing further on Oriente C, we find that it shares most drift with individuals from Northern Italy, Switzerland and Luxembourg, and less with individuals from Iberia, Scandinavia, and East and Southeast Europe (Fig. 4A-B). Shared drift decreases significantly with distance (Fig. 4C) and with time (Fig. 4D) although in a linear model of drift with distance and time as a covariate, only distance (p=1.3×10-6) and not time (p=0.11) is significant. Consistent with the overall E-W cline in hunter-gatherer ancestry, genetic distance to Oriente C increases more rapidly with longitude than latitude, although this may also be affected by geographic features. For example, Oriente C shares significantly more drift with the 8,000 year-old 1,400 km distant individual from Loschbour in Luxembourg (Lazaridis et al.,2014), than with the 9,000 year old individual from Vela Spila in Croatia (Mathieson et al.,2018) only 700 km away as shown by the D-statistic (Patterson et al.,2012) D (Mbuti, Oriente C, Vela Spila, Villabruna); Z=3.42. Oriente C’s heterozygosity was slightly lower than Villabruna (14% lower at 1240k transversion sites), but this difference is not significant (bootstrap P=0.12).

oriente-c-villabruna-f3-statistics
Multidimensional scaling of outgroup f3-statistics for Late 531 Upper Palaeolithic and Mesolithic hunter-gatherers.

Discussion and Conclusion

The robust record of radiocarbon dates proves that they reached Sicily not before 15-14 ka cal. BP, several millennia after the LGM peak. In our opinion, in fact, the hypothesis about an early colonization of Sicily by Aurignacians (Laplace, 1964; Chilardi et al., 1996) must be rejected, on the basis of a recent reinterpretation of the techno-typological features of the lithic industries from Riparo di Fontana Nuova (Martini et al., 2007; Lo Vetro and Martini, 2012; on this topic see also Di Maida et al., 2019).

These analyses have implications for understanding the origin and diffusion of the hunter-gatherers that inhabited Europe during the Late Upper Palaeolithic and Mesolithic. Our findings indicate that Oriente C shows a strong genetic relationship with Western European Late Upper Palaeolithic and Mesolithic hunter-gatherers, suggesting that the “Western hunter-gatherers” was a homogeneous population widely distributed in the Central Mediterranean, presumably as a consequence of continuous gene flow among different groups, or a range expansion following the LGM.

shared-drift-whg-villabruna-oriente-c
The same statistic as in A plotted with geographic position

The South Italian corridor

Once again, a hypothesis based on phylogeography – apart from scarce archaeological and palaeolinguistic data (“Semitic”-like topo-hydronymy and substrates in Europe) – seems to be confirmed step by step. Since the finding of the Villabruna individual of hg. R1b-L754 (likely R1b-V88, like south-eastern European lineages expanded with WHG ancestry), it was quite likely to find out that southern Europe was the origin of the expansion of R1b-V88 into Africa.

The most likely explanation for the presence of “archaic” R1b-V88 subclades among modern Sardinians was, therefore, that they represented a remnant from a Late Upper Palaeolithic/Early Mesolithic population that had not been replaced in subsequent migrations, and thus that the migration of these lineages into Northern Africa and the Green Sahara happened during a period when Italy was connected by a shallower Mediterranean (and more land connections) to Northern Africa.

late-epigravettian
Likely Late Epigravettian/Mesolithic expansion of R1b-V88 into Northern Africa. See full map.

Nevertheless, the arguments for a quite recent expansion of R1b-V88 through the Mediterranean and into Africa keep being repeated, probably based on ancestry from the few ancient (and many modern) populations that have been investigated to date, a simplistic approach prone to important errors that overarch whole migration models.

For example, in the recent paper by Marcus et al. (2019) the presence of these lineages among ancient Sardinians (from the late 4th millennium BC on) is interpreted as an expansion of R1b-V88 with the Cardial Neolithic based on their ancestry, disregarding the millennia-long gap between these samples and the presence of this haplogroup in Palaeolithic/Mesolithic Northern Iberia and Northern Italy, and the comparatively much earlier splits in the phylogenetic tree and dispersal among African populations.

Afroasiatic and Nostratic

I was asked recently if I really believed that we could reconstruct Proto-Nostratic and connect it with any ancestral population. My answer is simple: until the Chalcolithic – when the whole picture of Indo-Europeans, Uralians, Egyptians or Semites becomes quite clear – we have just very few (linguistic, archaeological, genetic) dots which we would like to connect, and we do so the best we can. The earlier the population and proto-language, the more difficult this task becomes.

NOTE. 1) I tentatively connected hg. R with Nostratic in a previous text – when it appeared that R1a expanded from around Lake Baikal, hence Eurasiatic; R1b from the south with AME-WHG ancestry, hence Afroasiatic; and R2 with Dravidian.

2) After that, I though it was more likely to be connected to AME ancestry and the Middle East, because of the apparent expansion of WHG from south-eastern Europe, and the potential association of Afroasiatic and (Elamo-?)Dravidian to Middle Eastern populations.

3) However, after finding more and more R1b samples expanding through northern Eurasia, spreading through the (then wider) steppe regions; and R1a essentially surviving among other groups in eastern Europe for thousands of years without being associated to significant migrations (like, say, hg. C after the Palaeolithic), it didn’t seem like this division was accurate, hence my most recent version.

But, in essence, it’s all about connecting the dots, and we have very few of them…

eurasiatic-phylum-ultraconserved-words
Phylogenetic tree from Pagel et al. (2013), partially in agreement with Kortlandt’s view on Eurasiatic. “Consensus phylogenetic tree of Eurasiatic superfamily (A) superimposed on Eurasia and (B) rooted tree with estimated dates of origin of families and of superfamily. (A) Unrooted consensus tree with branch lengths (solid lines) shown to scale and illustrating the correspondence between the tree and the contemporary north-south and east-west geographical positions of these language families. Abbreviations: P (proto) followed by initials of language family: PD, proto-Dravidian; PK, proto-Kartvelian; PU, proto-Uralic; PIE, proto–Indo-European; PA, proto-Altaic; PCK, proto–Chukchi-Kamchatkan; PIY, proto–Inuit-Yupik. The dotted line to PIY extends the inferred branch length into the area in which Inuit-Yupik languages are currently spoken: it is not a measure of divergence. The cross-hatched line to PK indicates that branch has been shortened (compare with B). The branch to proto-Dravidian ends in an area that Dravidian populations are thought to have occupied before the arrival of Indo-Europeans (see main text). (B) Consensus tree rooted using proto-Dravidian as the outgroup. The age at the root is 14.45 ± 1.75 kya (95% CI = 11.72–18.38 kya) or a slightly older 15.61 ± 2.29 kya (95% CI = 11.72–20.40 kya) if the tree is rooted with proto-Kartvelian. The age assumes midpoint rooting along the branch leading to proto-Dravidian (rooting closer to PD would produce an older root, and vice versa), and takes into account uncertainty around proto–Indo-European date of 8,700 ± 544 (SD) y following ref. 35 and the PCK date of 692 ± 67 (SD) y ago.”

In linguistics, I trust traditional linguists who tend to trust other more experimental linguists (like Hyllested or Kortlandt) who consider that – in their experience – an Indo-Uralic and a Eurasiatic phylum can be reconstructed. Similarly, linguists like Kortlandt are apparently (partially) supportive of attempts like that of Allan Bomhard with Nostratic – although almost everyone is critic of the Muscovite school‘s attachment to the Brugmannian reconstruction, stuck in pre-laryngeal Proto-Indo-Anatolian and similar archaisms.

I mostly use Nostratic as a way to give a simplistic ethnolinguistic label to the genetically related prehistoric peoples whose languages we will probably never know. I think it’s becoming clear that the strongest connection right now with the expansion of potential Eurasiatic dialects is offered by ANE-related populations (hence Y-chromosome bottlenecks under hg. R, Q, probably also N), however complicated the reconstruction of that hypothetic community (and its dialectalization) may be.

Therefore, the multiple expansions of lineages more or less closely associated to ANE-related peoples – like R1b-V88 in the case of Afrasian, or R2 in the case of Dravidians – are the easiest to link to the traditionally described Nostratic dialects and their highly hypothetic relationship.

green-sahara-neolithic
Reconstruction of North African vegetation during past green Sahara periods. Estimated and reconstructed MAP for the Holocene GSP (6–10 kyr BP) projected onto a cross-section along the eastern Sahara (left panel) and map view of reconstructed MAP, vegetation and physiographic elements [7,8,11,45] (right panel). Image from Larrasoaña et al. (2013).

What should be clear to anyone is that the attempt of many modern Afroasiatic speakers to connect their language to their own (or their own community’s main) haplogroups, frequently E and/or J, is flawed for many reasons; it was simplistic in the 2000s, but it is absurd after the advent of ancient DNA investigation and more recent investigation on SNP mutation rates. R1b-V88 should have been on the table of discussions about the expansion of Afroasiatic communities through the Green Sahara long ago, whether one supports a Nostratic phylum or not.

The fact that the role of R1b bottlenecks and expansions in the spread of Afroasiatic is usually not even discussed despite their likely connection with the most recent population expansions through the Green Sahara fitting a reasonable time frame for Proto-Afroasiatic reconstruction, a reasonable geographical homeland, and a compatible dialectal division – unlike many other proposed (E or J) subclades – reveals (once again) a lot about the reasons behind amateur interest in genetics.

Just like seeing the fixation in (and immobility of) recent writings about the role of I1, I2, or (more recently) R1a in the Proto-Indo-European expansion, R1b with Vasconic, or N1c with Proto-Uralic.

NOTE. That evident interest notwithstanding, it is undeniable that we have a much better understanding of the expansions of R1b subclades than other haplogroups, probably due in great part to the easier recovery of ancient DNA from Eurasia (and Europe in particular), for many different – sociopolitical, geographical, technological – reasons. It is quite possible that a more thorough temporal transect of ancient DNA from the Middle East and Africa might radically change our understanding of population movements, especially those related to the Afroasiatic expansion. I am referring in this post to interpretations based on the data we currently have, despite that potential R1b-based bias.

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

Global demographic history inferred from mitogenomes

Open access Global demographic history of human populations inferred from whole mitochondrial genomes, by Miller, Manica, and Amos, Royal Society Open Science (2018).

Relevant excerpts (emphasis mine):

Material

The Phase 3 sequence data from 20 populations, comprising five populations for each of the four main geographical regions of Europe, East Asia, South Asia and Africa, were downloaded from the 1000 Genomes Project website (www.1000genomes.org/data, [8]), including whole mitochondrial genome data for 1999 individuals. We decided not to analyse populations from the Americas due to the region’s complex history of admixture [13,14].

The European populations were as follows: Finnish sampled in Finland (FIN); European Caucasians resident in Utah, USA (CEU); British in England and Scotland (GBR); an Iberian population from Spain (IBS) and Toscani from Italy (TSI). Representing East Asia were the Han Chinese in Beijing (CHB); Southern Han Chinese (CHS); Dai Chinese from Xishuangbanna, China (CDX); Kinh population from Ho Chi Minh City, Vietnam (KHV) and Japanese from Tokyo (JPT). The South Asian populations were Punjabi Indians from Lahore, Pakistan (PJL); Gujarati Indians in Houston, USA (GIH) as well as Indian Telugu sampled in the UK (ITU); Bengali from Bangladesh (BEB) and Sri Lankan Tamil from the UK (STU). (…)

Method

We analysed our mtDNA data with the extended Bayesian skyline plot (EBSP) method, a Bayesian, non-parametric technique for inferring past population size fluctuations from genetic data. Building on the previous Bayesian skyline plot (BSP) approach, EBSP uses a piecewise-linear model and Markov chain Monte Carlo (MCMC) methods to reconstruct a populations’ demographic history [17] and is implemented in the software package BEAST v. 2.3.2 [11]. Alignments for each of the 20 populations were loaded separately into the Bayesian Evolutionary Analysis Utility tool (BEAUti v. 2.3.2) in NEXUS format.

1000-genomes-similarity-fst
Relationship between profile similarity and genetic distance, measured as Fst. Comparisons between regions, circles, are colour-coded: black ¼ AFR-EA; yellow ¼ AFR-EUR; blue ¼ AFR-SA; orange ¼ EUR-EA; green ¼ EA-SA; red ¼ EUR-SA. Comparisons within regions, squares, are coded: peach ¼ EUR; pink ¼ EA; dark blue ¼ EA; light blue ¼ AFR. Profile similarity is calculated as inferred size difference summed over 20 evenly spaced intervals (see Material and methods).

Regional demographic histories

Europe:

The five European profiles are presented in figure 2. The four southerly populations all show profiles with a stable size up to approximately 14 ka followed by a sudden, rapid increase that becomes progressively less steep towards the present. There is also a north-south trend, with confidence intervals becoming broader towards the north, particularly for the oldest time-points. The Finnish population profile appears rather different, but this is to be expected both because it is so far north and because previous studies have identified Finns as a strong genetic outlier in Europe [19–22].

europe-mtdna
Inferred demographic histories of five European populations. Dotted line is the median estimate of Ne and the thin grey lines show the boundary of the 95% CPD interval. The x-axis represents time from the present in years and all plots are on the same scale. Map shows origins of sampled populations.

South Asia:

The five profiles for South Asia are shown in figure 3. All populations reveal a period of rapid growth approximately 45–40 ka which then slows. Near the present the two southerly populations, GIH and STU both show evidence of a decline. However, this may be due to these samples being drawn from populations no longer living on the subcontinent, with the downward trend capturing a bottleneck associated with moving to Europe/America, perhaps accentuated by the tendency for immigrant populations to group by region, religion and race [23].

asia-mtdna
Inferred South Asian population demographic histories. Dotted line is the median Ne estimate and the thin grey lines show the boundary of the 95% CPD intervals. The x-axis represents time from the present in thousands of years and all plots are on the same scale. The map shows location of sampled populations.

Related

Long-term matrilineal continuity in a nonisolated region of Tuscany

iron_age_europe_mediterranean

New paper (behind paywall) The female ancestor’s tale: Long‐term matrilineal continuity in a nonisolated region of Tuscany, by Leonardi et al. Am J Phys Anthr (2018).

EDIT (10 SEP 2018): The main author has shared an open access link to read the PDF.

Interesting excerpts:

Here we analyze North-western Tuscany, a region that was a corridor of exchanges between Central Italy and the Western Mediterranean coast.

We newly obtained mitochondrial HVRI sequences from 28 individuals, and after gathering published data, we collected genetic information for 119 individuals from the region. Those span five periods during the last 5,000 years: Prehistory, Etruscan age, Roman age, Renaissance, and Present-day. We used serial coalescent simulations in an approximate Bayesian computation framework to test for continuity between the mentioned groups.

In all cases, a simple model of a long-term genealogical continuity proved to fit the data better, and sometimes much better, than the alternative hypothesis of discontinuity.

The low number of samples analyzed requires some caution in the interpretation. Because we did not test for gene flow, it is at this stage impossible to reject it, but our results suggest at least significant levels of genealogical continuity. Moreover, as it has not been possible to obtain more precise information on the age of the Eneolithic samples, they were grouped together considering the average archaeological period of interest, which may cause a bias in the analyses. (…)

tuscany
Geographic location of the samples considered in this work

(…) clearly, our samples show high levels of continuity when considering the whole Tuscan region as a genetic reservoir during the Iron Age.

The posterior distributions of the parameters confirm a high degree of genetic isolation in the sampled population, with very small values for the female effective population sizes across time. Such values, in particular the Neolithic ones, are in accord with the estimates obtained in similar studies, both in Tuscany (Ghirotto et al., 2013) and in France (Rivollat et al., 2017).

tuscany-genetic-diversity

Taken at their face value, our results do not show any major shift in the composition of the maternal ancestry of the population, across 50 centuries. This does not mean that no demographic process of relevance has affected the population, and indeed the higher diversity accumulating in time is the likely consequence of immigrating people, enriching the mitochondrial gene pool.

(…) the population of the current Lucca province appears to have retained very ancient mitochondrial features, despite occupying a geographical corridor between the Ligurian and the Tyrrhenian coast, and despite not showing the persistence of unique cultural traits through the centuries.

tuscany-genetic-diversity-hap

Another possibility is that that the different populations passing through the area (Etruscans, Romans, and Lombards) had a consistent social and/or sex bias. An example of similar patterns has been observed several times. Between the Late Neolithic and the Early Bronze Age, female exogamy in patrilocal society has been observed in Southern Germany (Knipper et al., 2017); during the Bronze Age the migrations toward Europe from the steppes appears to have consisted prevalently of males (Goldberg, Günther, Rosenberg, & Jakobsson, 2017); and in more recent periods in the Canary Islands, the female ancestry maintains a significant amount of autochthonous lineages, while the male ancestry was strongly influenced by the European colonization (Fregel et al., 2009, b).

It is well known that military invasions may not have a significant genetic impact upon the invaded population (Schiffels et al., 2016; Sokal, Oden, Walker, Di Giovanni, & Thomson, 1996;Weale,Weiss, Jager, Bradman, & Thomas, 2002), especially at the mitochondrial level, because of the limited size of a sustainable army, and of the fact that armies are generally composed mostly or only of males. Even if a substantial share of invaders decided to remain and settle the region, this form of gene flow would affect mostly or only the paternal lineages, rather than the maternal ones. We can also hypothesize the immigration of a number of people (e.g., Romans, Lombards) that may have acted as ruler of the region, remaining socially (and so genetically) separated by the local population, and leaving few (if any) traces in the gene pools of the local population.

tuscany-mtdna-haplogroups
Supporting Information, Table S1 New ancient samples genotyped

We expect to see that certain migrations since the Iron Age – like the Celtic and Roman ones – were somehow different from previous ones, where, at least since the Neolithic, male-dominated expansions were the rule.

If, however, male-biased expansions are also seen during the Iron Age – probably driven by particular subclades then – , this would certainly justify the continuity of admixture in certain regions in spite of these population expansions, and thus the importance of Y-DNA to track more recent language changes.

One of the most interesting details of the upcoming paper of Italic peoples will be the Y-DNA (and admixture) of Etruscans compared to other neighbouring peoples, given the known conflicting theories regarding their recent vs. older origin in the East before the historical record.

Related

Y-chromosome mixture in the modern Corsican population shows different migration layers

mesolithic-europe

Open access Prehistoric migrations through the Mediterranean basin shaped Corsican Y-chromosome diversity, by Di Cristofaro et al. PLOS One (2018).

Interesting excerpts:

This study included 321 samples from men throughout Corsica; samples from Provence and Tuscany were added to the cohort. All samples were typed for 92 Y-SNPs, and Y-STRs were also analyzed.

Haplogroup R represented approximately half of the lineages in both Corsican and Tuscan samples (respectively 51.8% and 45.3%) whereas it reached 90% in Provence. Sub-clade R1b1a1a2a1a2b-U152 predominated in North Corsica whereas R1b1a1a2a1a1-U106 was present in South Corsica. Both SNPs display clinal distributions of frequency variation in Europe, the U152 branch being most frequent in Switzerland, Italy, France and Western Poland. Calibrated branch lengths from whole Y chromosome sequencing [44,45] and ancient DNA studies [46] both indicated that R1a and R1b diversification began relatively recently, about 5 Kya, consistent with Bronze Age and Copper Age demographic expansion. TMRCA estimations are concordant with such expansion in Corsica.

corsica-haplogroups
Spatial frequency maps for haplogroups with frequencies above 3%, their Y-STR based phylogenetic networks in Corsican populations (Blue: North, Green: West, Orange: South, Black: Center and Purple: East) and their TMRCA (in years, +/- SE).

Haplogroup G reached 21.7% in Corsica and 13.3% in Tuscany. Sub-clade G2a2a1a2-L91 accounted for 11.3% of all haplogroups in Corsica yet was not present in Provence or in Tuscany. Thirty-four out of the 37 G2a2a1a2-L91 displayed a unique Y-STR profile, illustrated by the star-like profile of STR networks (Fig 1). G2a2a1a2-L91 and G2a2a-PF3147(xL91xM286) show their highest frequency in present day Sardinia and southern Corsica compared to low levels from Caucasus to Southern Europe, encompassing the Near and Middle East [21,47–50]. Ancient DNA results from Early and Middle Neolithic samples reported the presence of haplogroup G2a-P15 [51–53], consistent with gene flow from the Mediterranean region during the Neolithic transition. Td expansion time estimated by STR for P15-affiliated chromosomes was estimated to be 15,082+/-2217 years ago [49]. Ötzi, the 5,300-year-old Alpine mummy, was derived for the L91 SNP [21]. A genetic relationship between G haplogroups from Corsica and Sardinia is further supported by DYS19 duplication, reported in North Sardinia [14], and observed in the southern part of the Corsica in 9 out of 37 G2a2a1a2-L91 chromosomes and in 4 out of 5 G2a2a-PF3147(xL91xM286) chromosomes, 3 of which displayed an identical STR profile (S4 Table).

This lineage has a reported coalescent age estimated by whole sequencing in Sardinian samples of about 9,000 years ago. This could reflect common ancestors coming from the Caucasus and moving westward during the Neolithic period [48], whereas their continental counterparts would have been replaced by rapidly expanding populations associated with the Bronze Age [46,54,55]. Estimated TMRCA for L91 lineage in Corsica is 4529 +/- 853 years. G-L497 showed high frequencies in Corsica compared to Provence and Tuscany, and this haplogroup was common in Europe, but rare in Greece, Anatolia and the Middle East. Fifteen out of the 17 Corsican G2a2b2a1a1b-L497 displayed a unique Y-STR profile (S4 Table) with an estimated TMRCA of 6867 +/- 1294 years. Haplogroup G2a2b1-M406, associated with Impressed Ware Neolithic markers, along with J2a1-DYS445 = 6 and J2a1b1-M92 [22,49], had very low levels in Corsica. Conversely, G2a2b2a-P303was highly represented and seemed to be independent of the G2a2b1-M406 marker. The 7 G2a2b2a-P303(xL497xM527) Corsican chromosomes displayed a unique Y-STR profile (S4 Table).

pca-corsica
First and second axes of the PCA based on 12 Y-chromosome haplogroup frequencies in 83 west Mediterranean populations.

Haplogroup J, mainly represented by J2a1b-M67(xM92), displayed intermediate frequencies in Corsica compared to Tuscany and Provence. J2a1b-M67(xM92) derived STR network analysis displayed a quite homogeneous profile across the island with an estimated TMRCA of 2381 +/- 449 years (Fig 1) and individuals displaying M67 were peripheral compared to Northwestern Italians (S2 Fig). The haplogroup J2a1-Page55(xM67xM530), characteristic of non-Greek Anatolia [22], was found in the north-west of Corsica. Haplogroup J2a1-DYS445 = 6 was found in the north-west with DYS391 = 10 repeats, and in the far south with DYS391 = 9 repeats, the former was associated with Anatolian Greek samples, whereas the second was found in central Anatolia [22]. The 7 J2b2a-M241 displayed a unique Y-STR profile (S4 Table), they were only detected in the Cap Corse region, this sub-haplogroup shows frequency peaks in both the southern Balkans and northern-central Italy [56] and is associated with expansion from the Near East to the Balkans during Neolithic period [57].

Haplogroup E, mainly represented by E1b1b1a1b1a-V13, displayed intermediate frequencies in Corsica compared to Tuscany and Provence. E1b1b1a1b1a-V13 was thought to have initiated a pan-Mediterranean expansion 7,000 years ago starting from the Balkans [52] and its dispersal to the northern shore of the Mediterranean basin is consistent with the Greek Anatolian expansion to the western Mediterranean [22], characteristic of the region surrounding Alaria, and consistent with the TMRCA estimated in Corsica for this haplogroup. A few E1b1a-V38 chromosomes are also observed in the same regions as V13.

Related: