On the trail of the Proto-Uralic speakers (work in progress)

Historical linguists have long posited that Fennoscandia was a busy contact zone between early Germanic and Uralic languages. The first ancient DNA samples from what is now Finland have corroborated their inferences, by showing that during the Iron Age the western part of the country was inhabited by a genetically heterogeneous population closely related to both the Uralic-speaking Saami and Germanic-speaking southern Scandinavians.
The samples were sequenced and analyzed by two different teams of researches, and their findings published recently in Lamnidis et al. and Sikora et al. (see here and here, respectively).
This is how most of these ancients, whose remains were excavated from the Levanluhta burial site dated to 300–800 CE, behave in a Principal Component Analysis (PCA) based on my Global25 data. Levanluhta_IA are the Saami-related samples, while Levanluhta_IA_o is an Scandinavian-like outlier. Baltic_IA is an early Iron Age individual from what is now Lithuania from the recent Damgaard et al. paper (see here). Note the accuracy of the Global25 data in pinpointing their genetic affinities and also the trajectory of the Levanluhta_IA cluster, which seems to be “pulling” towards Levanluhta_IA_o.

The Saami and Levanluhta_IA are clear outliers from the main Northern European cluster. There are two reasons for this: excess East Asian/Siberian-related ancestry and Saami-specific genetic drift. However, this eastern admixture and genetic drift are shared in varying degrees by other North European populations, especially those that also speak Uralic languages, and this is why they appear to be “pulling” towards the Saami/Levanluhta_IA clusters in my PCA. Thus, what this suggests is that the expansion of Uralic languages across Northeastern Europe was intimately linked with the spread of Siberian-related ancestry into the region.
This idea has been around for a long time and is now becoming even more widely accepted (see here). However, Lamnidis et al. also featured samples from a likely pre-Uralic (1523±87 calBCE) burial site at Bolshoy Oleni Ostrov in the Kola Peninsula, present-day northern Russia, and, perhaps surprisingly, found that they showed even more Siberian-related ancestry than Levanluhta_IA. So what’s going on?
I’m confident that this discrepancy can be explained by multiple waves of migrations from the east into Northeastern Europe, possibly before, during and after the time of the people buried at Bolshoy Oleni Ostrov, by pre-Uralic, para-Uralic and/or Proto-Uralic-speaking populations.
Consider the following qpAdm output, in which Levanluhta_IA is just barely modeled successfully as a two-way mixture between Levanluhta_IA_o and Bolshoy_Oleni_Ostrov. The statistical fit improves significantly with the addition of Glazkovo_EBA as a third mixture source. This is an ancient population from near Lake Baikal dated to 4597-3726 BC from the aforementioned Damgaard et al. paper.

Levanluhta_IA
Bolshoy_Oleni_Ostrov 0.468±0.036
Levanluhta_IA_o 0.532±0.036
chisq 19.129
tail prob 0.0854706
Full output
Levanluhta_IA
Bolshoy_Oleni_Ostrov 0.241±0.092
Glazkovo_EBA 0.162±0.059
Levanluhta_IA_o 0.597±0.046
chisq 7.756
tail prob 0.734966
Full output

For the sake of being complete, I also tested whether Levanluhta_IA_o could be substituted by other similar ancient samples from the neighborhood, including those associated with the Battle-Axe and Corded Ware cultures. There’s not much to report; qpAdm returned poor statistical fits and/or implausible ancestry proportions (for the full output from my runs, see here). Baltic_IA did produce a statistically sound model, but with excess Glazkovo_EBA-related ancestry. I also had to drop Bolshoy_Oleni_Ostrov from the analysis to make things work, which suggests to me that the result shouldn’t be taken too literally.

Levanluhta_IA
Baltic_IA 0.677±0.034
Glazkovo_EBA 0.323±0.034
chisq 8.547
tail prob 0.741095
Full output

So as far as I can see, the western ancestry in Levanluhta_IA is likely to be mostly of Germanic origin, and thus Indo-European, meaning that it’s logical to look east, perhaps far to the east, for the source of its Uralic ancestry. This might seem like a complicated and uncertain task, considering that Levanluhta_IA could well be at least a thousand years younger than the first entry of Uralic speakers into Fennoscandia. However, take a look what happens when I substitute Glazkovo_EBA with a variety of Uralic-speaking populations from around the Ural Mountains, which is where the Proto-Uralic homeland is generally considered to have been located.

Levanluhta_IA
Bolshoy_Oleni_Ostrov 0.210±0.091
Khanty 0.283±0.090
Levanluhta_IA_o 0.507±0.035
chisq 7.007
tail prob 0.798532
Full output
Levanluhta_IA
Bolshoy_Oleni_Ostrov 0.193±0.098
Levanluhta_IA_o 0.495±0.035
Mansi 0.312±0.100
chisq 7.884
tail prob 0.7237
Full output
Levanluhta_IA
Bolshoy_Oleni_Ostrov 0.300±0.065
Levanluhta_IA_o 0.337±0.072
Mari 0.363±0.121
chisq 8.393
tail prob 0.677705
Full output
Levanluhta_IA
Bolshoy_Oleni_Ostrov 0.238±0.084
Levanluhta_IA_o 0.553±0.036
Nenets 0.209±0.067
chisq 7.210
tail prob 0.78181
Full output
Levanluhta_IA
Bolshoy_Oleni_Ostrov 0.302±0.069
Levanluhta_IA_o 0.324±0.081
Udmurt 0.373±0.135
chisq 9.195
tail prob 0.60393
Full output

All of these models look great, and easily rival the best model with Glazkovo_EBA. Moreover, they make good sense in terms of linguistics. The only problem is that they’re anachronistic, because the Uralic-speaking reference populations are younger than Levanluhta_IA. So I can’t be certain that they reflect reality without corroboration from ancient DNA. It might turn out, for instance, that an Glazkovo_EBA-like population was already present somewhere deep in Europe before or during the time of Bolshoy_Oleni_Ostrov, while no such population existed around the Ural Mountains until the time of Levanluhta_IA.
By the way, it might be important to note that the present-day Finnish samples in my dataset can’t be modeled as a mixture between Levanluhta_IA and Levanluhta_IA_o. But they can be modeled as a mixture between Baltic_IA and Levanluhta_IA. I don’t know which part of Finland they’re from exactly; probably all over the place, so it’d be useful to test regional Finnish populations to see how they behave in such models. Of course, Finns aren’t Saamic speakers, they’re Finnic speakers, and they’re probably the result of a more recent Uralic expansion into Fennoscandia than the one that gave rise to the Saami.

Finnish
Baltic_IA 0.671±0.076
Levanluhta_IA 0.329±0.076
chisq 14.114
tail prob 0.293508
Full output

Damgaard et al. didn’t report the Y-haplogroup for Baltic_IA, but the word round the campfire is that this individual belonged to N1c, which is today the most common Y-haplogroup among Uralic speakers. Obviously, we need a lot more ancient DNA to sort all of this out, but things are already looking pretty much as expected. Stay tuned for new posts in this series following the publication of more ancient DNA relevant to this fascinating topic.
See also…
The Uralic cline in the Global25
Late PIE ground zero now obvious; location of PIE homeland still uncertain, but…
Genetic and linguistic structure across space and time in Northern Europe

Source