Graphing the truth

I haven't used TreeMix since qpGraph became freely available for Linux. Among other things, the latter offers greater control, reproducibility and transparency.

However, I'd say that in its current form qpGraph is not the most objective way to analyze data. That's because if you're really good with it, and you want a graph to work, then often you can make it work by tweaking whatever it is that needs to be tweaked.

It's not possible to do a lot of tweaking with TreeMix. Indeed, once the user picks the samples for the TreeMix run, the rest of the process can be totally unsupervised, and thus free from human interference. Obviously, that's not a guarantee of accuracy, but it can be useful.

I feel I need to run more unsupervised analyses, especially when exploring new data. So to that end, I've dusted off TreeMix and will be using it regularly again.

There's been some talk lately online about migrations from Central Asia giving rise to the Eneolithic populations of the North Caucasus Piedmont steppe. In my opinion, that sounds like nonsense. But let's see what TreeMix has to say on the matter. In the graphs below look for the samples labeled Progress_En and Vonyuchka_En, respectively.

As far as I can tell, both of these graphs essentially corroborate the results from my recent Principal Component Analyses (PCA) with many of the same ancients (see here). In other words, Progress_En and Vonyuchka_En can be described as mixtures of populations closely related to the hunter-gatherers of the Caucasus on one hand, and those of Eastern Europe on the other. How does Central Asia fit into this, you might ask? It doesn't, unless you really want it to.

See also...

Did South Caspian hunter-fishers really migrate to Eastern Europe?

A final note for the year

I feel like I've spent a good part of 2019 banging my head against a thicker than average brick wall.

Much of this feeling is tied to the controversy over the ethnogenesis of the Yamnaya people, and my often futile attempts to explain that their origin cannot be sought in what is now Iran, or, indeed, anywhere outside of Eastern Europe.

This post is my final attempt to lay out the facts in regards to this topic. Next year I'll have better things to do than to argue the bleeding obvious.

Below are two graphs from a Principal Component Analysis (PCA) based on relatively high quality ancient human genotype data from the Caucasus and surrounds. They include two typical Yamnaya individuals from burial sites north of the Caspian Sea. I made the graphs with the Vahaduo Custom PCA tool here. The relevant datasheet can be downloaded here.

Here's what I'm seeing:

- the Yamnaya individuals sit on genetic clines made up of hunter-gatherers native to the Caucasus and various parts of Eastern Europe, including a trio from the southernmost part of the Pontic-Caspian steppe (labeled Steppe_Eneolithic), with whom they form a distinct cluster

- the samples from the Caucasus and the Iranian Plateau form very different clusters, so there's no support here for the ancient Caucasus/Iranian grouping that is often haphazardly invoked in scientific literature

- there's no indication that the Yamnaya and/or Steppe_Eneolithic groups experienced recent gene flow, or, for that matter, any gene flow whatsoever, from what is now Iran.

Of course, analyses based on formal statistics suggest that the Yamnaya population harbors minor western ancestry that is missing in Steppe_Eneolithic. In fact, I was first to argue this point (see here). So let's add a couple of ancient farmers from Western Europe to my PCA to see how they affect the graphs. The relevant datasheet is available here.

Yep, the Yamnaya pair appears to be peeling away very slightly, but deliberately, from the Steppe_Eneolithic individuals towards the part of the plot occupied by the farmers.

Admittedly, I'm no Sherlock Holmes, but even with my fairly average sleuthing abilities, I'm pretty sure I know how the Yamnaya people came to be. They formed largely on the base of a population very much like Steppe_Eneolithic somewhere deep in Eastern Europe, well to the north of the Caucasus, and nowhere near the Iranian Plateau.

See also...

A note on Steppe Maykop

A note on Steppe Maykop

I'm reading a new book titled Dispersals and Diversification: Linguistic and Archaeological Perspectives on the Early Stages of Indo-European (see here). One of the chapters is authored by archeologist David Anthony, in which he makes the following claims:

A previously unknown genetic population actually was identified in Wang et al. (2019), but it was a peculiar relict-seeming group related to Paleo Siberians and American Indians (Kennewick) that had survived isolated somewhere in the Caspian steppes or perhaps in the North Caucasus Mountains. The Maykop people did admix with this previously isolated Siberian/Kennewick population in graves labeled "Steppe Maykop" in Wang et al. (2019).

But this just makes it clearer that a cultural choice motivated the Maykop people to exclude marriages with Yamnaya and pre-Yamnaya people specifically, even while exchanges of material goods, ideas, technologies continued. Neither the Maykop nor the North Caucasus/Siberian/Kennewick population can be the source of most of the CHG [Caucasus hunter-gatherer] ancestry in Yamnaya. In order to narrow down when and where CHG ancestry entered the steppes, we must widen our geographic frame beyond the Caucasus.

Unfortunately, this is way off the mark. Especially unsound is his inference that the CHG-related ancestry in the Yamnaya population may have come from beyond the Caucasus.

In fact, the chances that the Steppe Maykop people were derived from a relict Siberian/Kennewick-related group that survived into the Maykop era in the Caspian steppes or the North Caucasus are exactly zero.

The real story was surely more complicated. In my opinion, it initially involved the migration during the Eneolithic or earlier of a people rich in CHG ancestry from the southernmost steppes into the Volga Delta and surrounds, and then the back-migration during the Early Bronze Age (EBA) of their descendants with around 50% admixture from Central Asian foragers. If so, these foragers were very similar to indigenous West Siberians and also relatively closely related to Native Americans.

I don't know why such an exotic people migrated into the North Caucasus steppes to form the bulk of the Steppe Maykop population, but I'm certain they did, and one interesting possibility is that they were recruited by Maykop chiefs to create a buffer zone against hostile Yamnaya-related groups trying to push into the Caucasus, possibly from the lower Don region.

Of course, the same ancient northward migration of the CHG-rich population that may have eventually given rise to the Steppe Maykop people might also explain the deep origins of the Yamnaya people.

The key sample in all of this is VJ1001 from the Wang et al. paper. This female comes from an Eneolithic (4332-4238 calBCE) kurgan burial in the North Caucasus steppes. But despite her early date, she's genetically very similar to most Yamnaya individuals. And she's also a perfect proxy for half of the ancestry of three out of the six Steppe Maykop individuals. Here's a mixture model that I put together using the Broad MIT/Harvard software qpAdm:

RUS_Steppe_Maykop (3/6)
RUS_Eneolithic_steppe_VJ1001 0.452±0.023
RUS_Tyumen_HG 0.548±0.023
chisq 7.494
tail prob 0.874914
Full output

Indeed, these Steppe Maykop samples don't harbor any Maykop ancestry. They're simply a two-way mixture between a population closely resembling VJ1001 and another one similar to hunter-gatherers from Tyumen, West Siberia.

Importantly, a couple of Steppe Maykop-related populations were inadvertently discovered by Narasimhan et al. northeast of the Caspian Sea in what is now Kazakhstan. One of these groups is labeled Kumsay_EBA, after the location of its cemetery. It's roughly contemporaneous with Steppe Maykop and basically identical to the aforementioned Steppe Maykop trio.

KAZ_Kumsay_EBA
RUS_Eneolithic_steppe_VJ1001 0.440±0.022
RUS_Tyumen_HG 0.560±0.022
chisq 10.573
tail prob 0.646513
Full output

I suppose it's possible that Kumsay_EBA represents the migration of Steppe Maykop people into the Kazakh steppes. But even if this is true, then there had to have been an earlier migration of a group from the Kazakh steppes or West Siberia that mixed with the VJ1001-related natives of the North Caucasus steppes to give rise to Steppe Maykop.

I'm assuming that the Yamnaya-like VJ1001 and her people were the indigenous population of the North Caucasus steppes because there are no indications that they or their ancestors migrated there within any reasonable time frame from anywhere else, and certainly not from as far afield as, say, what is now Iran.

The other three Steppe Maykop individuals, who are genetic outliers in varying degrees from the main Steppe Makyop cluster, show variable levels of Maykop ancestry, with an average of about 50%. But they too harbor significant VJ1001-related ancestry. So despite the fact that there was some irregular mixing between the Maykop and Steppe Maykop peoples, this is not what created the typical Steppe Maykop genetic profile.

RUS_Steppe_Maykop_o
RUS_Eneolithic_steppe_VJ1001 0.234±0.074
RUS_Maykop_Novosvobodnaya 0.461±0.046
RUS_Tyumen_HG 0.305±0.033
chisq 7.378
tail prob 0.831667
Full output

And, of course, it should be obvious by now that the ancestry of the vast majority of Yamnaya individuals is better modeled without any input whatsoever from the Maykop or Steppe Maykop samples.

In fact, early indications are that the Yamnaya people flooded into Steppe Maykop territory from the north and completely replaced its population (see here). Despite this, in Dispersals and Diversification archeologist Kristian Kristiansen makes the following claim: "steppe Maykop expanded north, leading to the formation of the Yamnaya Culture and Proto-Indo-European". Not a chance in hell Professor.

See also...

A final note for the year

The PIE homeland controversy: August 2019 status report

Some myths die hard

An exceptional burial indeed, but not that of an Indo-European

Avalon vs Valhalla revisited

Pictured below is a new version of my Celtic vs Germanic genetic map. It's based on the same Principal Component Analysis (PCA) as the original (which can be seen here), but more focused on Northwestern Europe and produced with a different program.

To see the interactive online version, navigate to Vahaduo Custom PCA and copy paste the text from here into the empty space under the PCA DATA tab. Then press the PLOT PCA button under the PCA PLOT tab. For more guidance, refer to the screen caps here and here.

To include a wider range of populations in the key, just edit the data accordingly. For instance, to break up the ancient grouping into more specific populations, delete the Ancient: prefix in all of the relevant rows. This is what you should see:

Conversely, you can leave the ancient sample set intact and instead reorder the present-day linguistic groupings into, say, geographic groupings. To achieve this just delete all of the linguistic prefixes, such as Celtic:, Germanic:, and so on. You should end up with a datasheet like this and plot like this.

Of course, you can design your own plot by using any combination of the ancient and present-day individuals and populations that I've already run in this PCA. Their coordinates are listed here. Indeed, if you're in the possession of your own Celtic vs Germanic PCA coordinates, you can add yourself to the plot. And if you're not, see here.

It's also possible to re-process PCA data via the SOURCE tab. But I don't recommend doing this with the Celtic vs Germanic data, which are derived from a fine scale analysis and don't pack much variation. On the other hand, Global25 data are ideal for such re-processing. I made the plots below from subsets of Global25 coordinates available in a zip file here. To see how, refer to the screen caps here and here.

See also...

Modeling your ancestry has never been easier

Getting the most out of the Global25

Modeling genetic ancestry with Davidski: step by step

Big deal of 2019: ancient DNA confirms the link between Y-haplogroup N and Uralic expansions

The academic consensus is that Indo-European languages first spread into the Baltic region from the Eastern European steppes along with the Corded Ware culture (CWC) and its people during the Late Neolithic, well before the expansion of Uralic speakers into Fennoscandia and surrounds, probably from somewhere around the Ural Mountains.

On the other hand, the views that the Uralic language family is native to Northern Europe and/or closely associated with the CWC are fringe theories usually espoused by people not familiar with the topic or, unfortunately it has to be said, mentally unstable trolls.

The likely close relationship between the CWC expansion and the early spread of Indo-European languages was discussed in several papers in recent years (for instance, see here). This year, we saw the first ancient DNA paper focusing on the transition from the Bronze Age to the Iron Age in the East Baltic, including the likely first arrival of Uralic speech in what is now Estonia.

Published in Current Biology courtesy of Saag et al., the paper showed that the genetic structure of present-day East Baltic populations largely formed in the Iron Age (see here). It was during this time, the authors revealed, that the region experienced a sudden influx of Y-chromosome haplogroup N, which is today common in many Uralic speaking populations and often referred to as a Proto-Uralic marker. Little wonder then that Saag et al. linked this genetic shift in the East Baltic to the westward migrations of early Uralic speakers.

The table below, based on data from the Saag et al. paper, surely doesn't leave much to the imagination about what happened.

Unfortunately, I have to say that the genome-wide analysis in the paper was less informative than it could have been. The authors focused their attention on rather broad genetic components, and, as a result, missed an interesting fine scale distinction between their Bronze Age and Iron Age samples. The spatial maps below, based on my Global25 data for most of the ancients from Saag et al., show what I mean. The hotter the color the higher the genetic similarity between them and present-day West Eurasian populations.

Note that the Bronze Age (Baltic_EST_BA) samples are most similar to the Baltic-speaking, and thus also Indo-European-speaking, Latvians and Lithuanians, rather than the Uralic-speaking Estonians, even though they're from burial sites in Estonia. On the other hand, the Iron Age (Baltic_EST_IA) samples show strong similarity to a wider range of populations, including Estonians and many other Uralic-speaking groups.

See also...

It was always going to be this way

Fresh off the sledge

More on the association between Uralic expansions and Y-haplogroup N

Viking Age Iceland

I finally managed to get some of the Icelandic ancients from Ebenesersdóttir et al. 2018 into the Global25 datasheets (see here). Better late than never. Look for the"ISL_Viking_Age" prefix. Below is a screen cap of a Principal Component Analysis (PCA) with the new samples. It was done with an online Global25 PCA runner freely available here.

The individuals classified as unadmixed Gaels and Norse by Ebenesersdóttir et al. generally also look like it based on their Global25 coordinates.

The mixture models below, using all of the populations from the Global25 "modern pop averages scaled" datasheet, were run with an online tool freely available here. Note that the ADD DIST COL option is set to 1X. This is a useful feature for modeling the fine scale ancestry of samples that are derived from very similar populations.

See also...

They came, they saw, and they mixed

Commoner or elite?

Who were the people of the Nordic Bronze Age?

Etruscans, Latins, Romans and others

I've just added coordinates for more than 100 ancient genomes from the recently published Antonio et al. ancient Rome paper to the Global25 datasheets. Look for the population and individual codes listed here. Same links as always:

Global25 datasheet ancient scaled

Global25 pop averages ancient scaled

Global25 datasheet ancient

Global25 pop averages ancient

Thus far I've only managed to check a handful of the coordinates, so please let me know if you spot any issues. Below is a Principal Component Analysis (PCA) featuring the Etruscan and Italic speakers. I ran the PCA with an online tool specifically designed for Global25 coordinates freely available here.

Can we say anything useful about the origins of the Etruscan and early Italic populations thanks to these new genomes? Also, to reiterate my question from the last blog post, what are the genetic differences exactly between the Etruscans, early Latins, Romans and present-day Italians? Feel free to let me know in the comments below.

Update 13/11/2019: Here's another, similar PCA. This one, however, is based on genotype data, and it also highlights many more of the samples from the Antonio et al. paper. Considering these results, I'm tempted to say that the present-day Italian gene pool largely formed in the Iron Age, and that it was only augmented by population movements during later periods. The relevant datasheet is available here.

Update 13/11/2019: It seems to me that the two Latini-associated outliers show significant ancestry from the Levant, which possibly means that they're in part of Phoenician origin. These qpAdm models speak for themselves:

ITA_Ardea_Latini_IA_o
ITA_Proto-Villanovan 0.547±0.081
Levant_ISR_Ashkelon_IA2 0.453±0.081
chisq 7.573
tail prob 0.87027
Full output

ITA_Prenestini_tribe_IA_o
ITA_Proto-Villanovan 0.679±0.068
Levant_ISR_Ashkelon_IA2 0.321±0.068
chisq 7.222
tail prob 0.89033
Full output

The Proto-Villanovan singleton is also a key part of the models. Dating to the Bronze Age/Iron Age transition, she appears to be of western Balkan origin. Moreover, her steppe ancestry is probably derived directly from the Yamnaya horizon.

ITA_Proto-Villanovan
HRV_Vucedol 0.677±0.031
Yamnaya_RUS_Samara 0.323±0.031
chisq 10.397
tail prob 0.661174
Full output

The cluster made up of four early Italic speakers can be modeled with minor Proto-Villanovan-related ancestry, but, perhaps crucially, it doesn't need to be. Indeed, judging by the qpAdm output below, it's possible that almost all of its steppe ancestry came from the Bell Beaker complex, and, thus, the Corded Ware culture complex before that.

ITA_Italic_IA
Bell_Beaker_Mittelelbe-Saale 0.480±0.055
ITA_Grotta_Continenza_CA 0.411±0.042
ITA_Proto-Villanovan 0.109±0.084
chisq 10.294
tail prob 0.590205
Full output

Two out of the three available Etruscans look very similar to the Italic speakers in the above PCA plots, and yet they show a lot more Proto-Villanovan-related ancestry in my qpAdm run. The statistical fit is also relatively poor, perhaps suggesting that something important is missing.

ITA_Etruscan
Bell_Beaker_Mittelelbe-Saale 0.186±0.081
ITA_Grotta_Continenza_CA 0.283±0.064
ITA_Proto-Villanovan 0.531±0.126
chisq 17.175
tail prob 0.143143
Full output

Interestingly, the Etruscan outlier with significant North African admixture (proxied in my run by MAR_LN) doesn't need to be modeled with any Bell Beaker ancestry.

ITA_Etruscan_o
ITA_Proto-Villanovan 0.675±0.057
MAR_LN 0.325±0.057
chisq 14.864
tail prob 0.315912
Full output

Update 17/11/2019: The spatial maps below show how three groups of ancient Romans (from the Imperial, Late Antiquity and Medieval periods) compare to present-day West Eurasian populations in terms of their Global25 coordinates. The hotter the color, the higher the similarity. More here.

See also...

Getting the most out of the Global25

Modeling your ancestry has never been easier

An exceedingly simple, yet feature-packed, online tool ideal for modeling ancestry with Global25 coordinates is freely available HERE. It works offline too, after downloading the web page onto your computer. Just copy paste the coordinates of your choice under the "source" and "target" tabs, and then mess around with the buttons to see what happens. The screen caps below show me doing just that.

Another free, easy to use online tool that works with Global25 coordinates is the Principal Component Analysis (PCA) runner HERE. Below is a screen cap of me checking out one of the many PCA that it offers.

See also...

Getting the most out of the Global25

The Balkan connection

The hot topic at the moment is social inequality in Bronze Age Europe, thanks to a new paper by Mittnik et al. at Science. The full article is sitting behind an exceedingly robust paywall here.

However, the genotype dataset from the paper is freely available at the Max Planck Society's Edmond data repository here. Below is my Principal Component Analysis (PCA) of ancient West Eurasian genetic variation featuring 41 of the highest quality ancients from the new dataset. Almost all of them are from the Lech Valley in the Bavarian Alps, covering the period from the Bell Beaker culture (BBC) to the Middle Bronze Age (MBA). Two of the samples are from a mass Corded Ware culture (CWC) burial in the more northerly Tauber Valley.

I've also highlighted other ancients on the plot associated with the BBC and CWC from present-day Netherlands and Germany, respectively. The relevant PCA datasheet can be downloaded here.

Social stratification in ancient Europe is a fascinating topic, and it's an issue that I've started looking at myself (see here). However, I can't see any correlation between the inferred social standing of the individuals from the Lech and Tauber valleys and their positions in my PCA.

Nevertheless, the PCA is interesting in that it highlights considerable genetic heterogeneity within the Lech Valley BBC population. Indeed, how is this heterogeneity even possible, if, as per Mittnik et al., ancient DNA "has shown that the spread of the BBC throughout continental Europe did not involve large-scale migrations"?

Below is another version of my PCA, but this time focusing on three males: Lech Valley Beakers UNTA58_68Sk1 and WEHR_1192SkA, as well as ALT_4 from the aforementioned mass CWC grave in the Tauber Valley. Note that UNTA58_68Sk1 and WEHR_1192SkA represent genetically the most southern and northern, respectively, Lech Valley BBC samples that had enough data to be run in my analysis. I chose to focus on males because they carry the Y-chromosome, which can be informative about male-mediated ancient population expansions.

The PCA outcomes for these individuals are generally in line with their results in other types of genetic analyses, including those based on formal statistics. For instance, compared to the other two, ALT_4 harbors excess early steppe herder ancestry, UNTA58_68Sk1 excess early European farmer ancestry, and WEHR_1192SkA excess European hunter-gatherer ancestry. Moreover...

- UNTA58_68Sk1 shows a non-local isotopic signature and belongs to Y-haplogroup G2a, a marker essentially missing from BBC populations north of the Alps, and is best modeled as a two-way mixture between Bronze Age populations from the Balkans and the Pontic-Caspian steppe (see here), which probably means that he was a migrant to the Lech Valley from south of the Alps

- importantly, UNTA58_68Sk1 is not an isolated case, at least in the sense that several other BBC individuals from Bavaria, Bohemia, Hungary and Poland show varying ratios of Balkan-related ancestry, although almost all of these people are women

- WEHR_1192SkA is very similar to Bell Beakers from the northern Netherlands with whom he shares the R1b-P312 Y-haplogroup, suggesting that he was part of a population that moved into the Lech Valley from potentially as far away as the North Sea coast

- although ALT_4 probably shares the R1b-L51 Y-haplogroup with WEHR_1192SkA and many other BBC and Bronze Age individuals from the Bavarian Alps and surrounds, this can't be used as evidence of significant local genetic continuity after the CWC period, especially considering the comparatively eastern genome-wide structure of ALT_4.

Of course, archeological data suggest that the BBC was influenced in some important ways by the Copper and Bronze Age cultures of the Balkans and Carpathian Basin. So much so, in fact, that Marija Gimbutas, author of The Civilization of the Goddess, believed that the BBC originated in the Balkans from a synthesis of the local Vucedol culture and the intrusive Yamnaya culture from the Pontic-Caspian steppe.

Considering the ancient DNA evidence, however, the main demographic center of the early BBC could not have been south of the Alps.

Rather, it appears that early BBC and even CWC groups from north of the Alps moved into the Balkans and Carpathian Basin, where they may have established contacts with the local elites. If so, this might explain the significant southern cultural influences on the BBC, but limited accompanying genetic impact. This scenario also has support from archeological data (for instance, see here).

See also...

Is Yamnaya overrated?

The Boscombe Bowmen

Single Grave > Bell Beakers

Y-haplogroup R1a and mental health

I've updated my map of pre-Corded Ware culture R1a samples with a couple of new entries from Central and South Asia (the original is still here). However, before any of you get overly excited, please note that these samples aren't older than the Corded Ware culture. The reason I added them to my map is to counter the ongoing absurd claims online that South Asian R1a isn't derived from European R1a.

Just in case the map can't be viewed in all of its glory in some devices, here's what the fine print says:

The oldest example of R1a in ancient DNA from Central Asia is dated to 2132-1940 calBCE (ID I3770, Narasimhan 2019). Moreover, this sequence is closely related to much older R1a samples from Central, Eastern and Northern Europe, and phylogenetically nested within their diversity. Thus, it must surely represent a population expansion from Europe to Central Asia. Indeed, it's also associated with the Bronze Age Andronovo archeological culture, which is usually seen as an offshoot of the Corded Ware culture (CWC) of Late Neolithic Europe. The vast majority of present-day R1a lineages in Central Asia are closely related to that of I3770, and so must also ultimately derive from Europe.

The oldest instance of R1a in ancient DNA from South Asia is dated to just 1044-922 calBCE (ID I12457, Narasimhan 2019). This sequence, as well as the vast majority of present-day South Asian R1a lineages, are closely related to much older R1a samples from Central, Eastern and Northern Europe, and phylogenetically nested within their diversity. Thus, they must surely represent a population expansion from Europe to South Asia via Central Asia, in all likelihood during the Bronze Age. Even if R1a existed in South Asia before the Bronze Age, which is extremely unlikely, because it's found in samples from indigenous European hunter-gatherers, the vast majority of present-day R1a lineages in South Asia must be ultimately from Europe.

The idea that most, if not all, South Asian R1a is derived from European R1a seriously scares a lot of people. This is obvious in many online discussions on the topic. I suspect they're so frightened by it because, in their minds, it has the potential to encourage discrimination and even racism, perhaps by re-defining the colonization of much of the world by European nations in the recent past as the natural order of things?

In any case, clearly we're dealing with some sort of mass phobia here. I've got advice for those of you suffering from this problem: if you're honestly worried that the geographic provenance and expansion history of some Y-haplogroup is going to negatively impact on your life in any meaningful way, then it's time to find yourself a quality mental health professional. All the best with that.

See also...

The mystery of the Sintashta people

The Poltavka outlier

Yamnaya isn't from Iran just like R1a isn't from India