Next year the Armenian Plateau hypothesis will collapse

It's been a great year for population genetics and paleogenomics, and also for this blog. I ran a lot of analyses in 2015 and managed to make a few discoveries that were subsequently confirmed, or at least, backed up by academia. For instance:

- first to show with ancient genomes that the Anglo-Saxons made a significant genetic impact on England. See here. Eventually confirmed here.

- first to show that the southern admixture in the Yamnaya pastoralists of the Early Bronze Age steppe was Georgian-related rather than Armenian-related. See here. Confirmed here.

- first to show that Anatolian Neolithic farmers were very similar to European Neolithic farmers, and lacked Ancient North Eurasian (ANE) ancestry. See here. Confirmed here.

- first to show using ancient DNA and formal statistics that South Asia experienced massive gene flow originating in Late Neolithic/Bronze Age Europe. See here. Backed up with my help here.

The fact that Caucasus hunter-gatherers (CHG) like Kotias are essentially an ideal fit for the southern ancestry in the Yamnaya is a big problem for the Armenian Plateau Indo-European homeland hypothesis. This TreeMix graph shows why.

Basically, it looks like the Kotias-related ancestry in the Yamnaya came from the North Caucasus, rather than any place closer to the Near East than Georgia. Unless, of course, the southern Caucasus was populated by unadmixed CHG right until the 4th Millennium BC, when the hypothetical Proto-Indo-Europeans from the Armenian Plateau set off on their journey to Northern Europe around the Caspian Sea. But let's be honest, that's extremely unlikely.

Indeed, I expect that next year we'll see the first Neolithic and Copper Age samples from Armenia and/or surrounds, and even though they will be in large part CHG, they'll be nowhere near unadmixed. This will essentially kill the Armenian Plateau hypothesis, and thus leave the Kurgan or steppe hypothesis as the only plausible choice.

In any case, 2016 will probably be the year when ancient DNA helps to settle the Indo-European homeland question once and for all. So get ready for more ancient DNA from the steppe, but also, among others, from Mesolithic and Neolithic Iran, Mycenaean Greece and the Maikop Culture of the North Caucasus. I'm also pretty sure that the Varna man with the golden codpiece will make an appearance in a paper about Neolithic and Copper Age Bulgaria. Bring it on!

At least three genetically distinct Indo-European migrations into South Asia

First came the Indo-Aryans, probably in a couple of waves. Historical linguistics and archeology tell us that they originated on the Trans-Urals steppe in the Sintashta-Andronovo horizon, and pushed south around 2,000 BC to establish themselves as the ruling elite over Central Asian agriculturalists, who were probably in large part of West Asian origin.

There are multiple lines of genetic evidence suggesting that this is indeed what happened, which I discussed in detail in several earlier blog posts, like here.

But arguably the easiest way to show it is with D-stats of the form D(Indo-Aryan,Southeast_Asian; X,Outgroup), where the Indo-Aryans are the Kalash, a population isolate from the Hindu Kush with a relatively low level of extra-West Eurasian admixture and speaking an archaic form of Indo-Aryan. The Southeast Asians are the Dai from southern China, one of the best proxies for the South and East Asian admixture in the Kalash, while X represents a wide variety of present-day and ancient populations in my dataset. The top five D-stats, each based on well over 500K SNPs, are listed below:

Kalash Dai Kotias Ju_hoan_North 0.0684 22.704
Kalash Dai Sintashta Ju_hoan_North 0.0632 25.036
Kalash Dai Georgian Ju_hoan_North 0.0625 30.991
Kalash Dai Afanasievo Ju_hoan_North 0.0612 24.496
Kalash Dai Yamnaya_Samara Ju_hoan_North 0.0611 27.97

Really cool results. Obviously, Kotias is the recently published Caucasus hunter-gatherer (CHG) genome. The Kalash appear to carry the highest level of Kotias-related ancestry among present-day populations, which they probably acquired from both the Central Asian agriculturists and Indo-Aryan invaders. At the same time, however, Georgians show the highest affinity to Kotias because they harbor less extra-West Eurasian admixture.

After the Indo-Aryans came the Iranians, in all likelihood also from the steppe. They were either an offshoot of Sintashta-Andronovo or the more westerly Srubnaya Culture. I'd say the D-stats below, of the form D(Eastern_Iranian,Southeast_Asian)(X,Outgroup), are inconclusive, because the differences are small, and the outcome possibly affected by the methodology and/or sampling bias.

Tajik_Shugnan Dai Sintashta Ju_hoan_North 0.0716 26.427
Tajik_Shugnan Dai Poltavka Ju_hoan_North 0.0695 25.234
Tajik_Shugnan Dai Afanasievo Ju_hoan_North 0.0691 24.703
Tajik_Shugnan Dai Srubnaya Ju_hoan_North 0.069 28.266
Tajik_Shugnan Dai Corded_Ware_Germany Ju_hoan_North 0.0684 27.328

But again, the top five results make a lot of sense in the context of historical linguistics and archeology. By the way, Tajik Shugnans are a population isolate in the Pamir Mountains, like the Kalash with low level extra-West Eurasian admixture, and thus likely to be among the best available reference groups for early Eastern Iranians.

Interestingly, based on that list the Shugnans look more European than the Kalash. In large part this might be a reflection of the sharp rise in the level of European-specific Western hunter-gatherer (WHG) admixture on the steppe during the Middle Bronze Age, probably caused by population movements originating at the western edge of the steppe and/or in East Central Europe.

As far as I can tell, the fact that the Shugnans and Kalash have around the same level of extra-West Eurasian admixture means that I can try to hone in on the differences between their steppe-derived ancestry with D-stats of the form D(Kalash,Tajik_Shugnan)(Kotias,X). The top result seems to confirm my hunch, because Loschbour is, of course, a Western hunter-gatherer.

Loschbour 0.0149 3.874
Basque_Spanish 0.0113 4.232
Anatolia_Neolithic 0.0112 4.257
Karelia_HG 0.0105 3.005
Poltavka 0.01 3.539
Corded_Ware_Germany 0.0099 3.734
Afanasievo 0.0094 3.213
Srubnaya 0.0094 3.538
Yamnaya_Kalmykia 0.0091 3.362
Albanian 0.0088 3.419
Altai_IA 0.0088 3.087
Sintashta 0.0088 3.146
Greek 0.0076 3.094

Full output available here

More recently, during historic times, large parts of northern South Asia were settled by the Balochi, a Western Iranian people from the South Caspian region, whose ancestors were probably Indo-Europeanized a couple millennia earlier by Proto-Iranians from the steppe moving west across the Iranian Plateau. D-stats comparing the Balochi to the Kalash and Shugnans, respectively, clearly reflect the Near Eastern origins of the Balochi.

BedouinB 0.0104 6.151
Anatolia_Neolithic 0.0094 5.495
Druze 0.0084 5.228
Cypriot 0.0082 4.839
Syrian 0.0079 4.714
Armenian 0.0063 3.935
Satsurblia 0.0059 2.472
Georgian 0.0055 3.443
Iranian 0.0055 3.345
Abkhasian 0.0053 3.279
Greek 0.0052 3.166

...

Okunevo -0.0081 -3.552
Karelia_HG -0.0104 -4.666

Full output available here

Satsurblia 0.007 2.078
BedouinB 0.0051 2.277
...

Basque_Spanish -0.0073 -3.156
Mezhovskaya -0.0085 -3.045
Altai_IA -0.0092 -3.677
Scythian_IA -0.0092 -3.108
Yamnaya_Samara -0.0095 -4.092
Karitiana -0.0098 -3.501
Karasuk -0.0099 -4.322
Andronovo -0.01 -4.09
Sintashta -0.01 -3.951
Corded_Ware_Germany -0.0102 -4.34
Srubnaya -0.0106 -4.605
Yamnaya_Kalmykia -0.011 -4.511
MA1 -0.0118 -3.691
Okunevo -0.0122 -3.783
Poltavka -0.0125 -5.043
Afanasievo -0.0136 -5.235
Loschbour -0.0148 -4.201
Karelia_HG -0.0208 -6.537

Full output available here

In this analysis I used ancient samples from the recently published Jones et al. and Mathieson et al. studies, available on request from the authors and at the Reich lab website here, respectively. The present-day samples are from the Human Origins dataset, also available at the Reich lab website.

The Scythian

It's time to have a look at the Scythian steppe warrior from the Mathieson et al. dataset. This is the first Scythian individual to be genotyped.

He comes from the eastern end of the Pontic-Caspian steppe, is dated to 380-200 calBCE, and belongs to Y-chromosome haplogroup R1a, which is the dominant Y-haplogroup in Scythian and related remains tested to date.

His genome-wide data puts him closest to Northeast and Northwest Europeans from among present-day populations, rather than West and South Asians, who should, in theory, carry significant Scythian ancestry. We can probably put this down to the complex ancestry of West and South Asians.

Moreover, he can be modeled as a mixture of the Middle Bronze Age Potapovka people of the Pontic-Caspian steppe and present-day Nganasans of Siberia. This gels rather nicely with archaeological evidence, which suggests that Scythians were the descendants of Bronze Age Eastern European migrants to South Siberia, who expanded west across the Eurasian steppe during the Iron Age and eventually ended up back in Europe.

Identical-by-State (IBS) similarity

Lithuanian 0.645247
Estonian 0.645233
Latvian 0.645024
Russian_Kostroma 0.644946
Irish 0.644902
Orcadian 0.644792
Norwegian 0.644754
Belorussian 0.644727
Swedish 0.644667
Polish 0.644664
Austrian 0.644639
Danish 0.644587
English_Cornwall 0.644556
Belgian 0.644552
Scottish_Argyll 0.644548

Full output available here

Outgroup f3 shared drift statistics

Estonian 0.313726
Latvian 0.313664
Lithuanian 0.313574
Russian_Orel 0.313346
Finnish_Southwest 0.312997
Orcadian 0.312768
Norwegian 0.312768
Belorussian 0.312676
Russian_Kostroma 0.312669
Swedish 0.312608
Karelian 0.312567
Polish 0.31243
Irish 0.312281
Polish_Estonian 0.312156
Finnish 0.312102

Full output available here

qpAdm mixture model

Scythian_IA
Potapovka 0.913
Nganasan 0.087
chisq 5.815 tail prob 0.213365

Full output available here

Citation...

Mathieson et al., Genome-wide patterns of selection in 230 ancient Eurasians, Nature, Published online 23 November, 2015doi:10.1038/nature16152. Genotype dataset available here.

See also...

Cimmerians, Scythians and Sarmatians came from...

Genetic origins and legacy of the Scythians and Sarmatians

Mixed marriages on the early Eneolithic steppe

It looks like the Sredny Stog culture was the early vector for the spread of both Anatolian Neolithic and Caucasus hunter-gatherer (CHG) admixture onto the steppe, from the west and east, respectively:

These data testify the assumption about the existence of mixed Tripolye-Sredniy Stog marriages, because Tripolye population represented the Mediterranean anthropological type according to the not numerous Tripolye burials (Потехина 1999, c.154). It is interesting, that the massive Protoeuropoid type was typical for the oldest and the most eastern monuments of Sredniy Stog, while mesomorphic Mediterranean type was typical for the Igren cemetery, which was one of the youngest monuments related to the second and third periods of the Sredniy Stog culture and synchronous to the Tripolye B I and B I-II.

...

Appearance of pottery with pearls at the settlements of the third period of Sredniy Stog culture and glossy ceramics without ornamentation in the eastern variant sites, as well as the group of vessels with the steppe traces at the Svobodnoe settlement, allow me to assume the existence of mixed marriages between the Sredniy Stog and Northern Caucasus population.

Source: Early Eneolithic in the Pontic Steppes, book by Nadezhda Sergeenva Kotova, available at Academia.edu here.

See also...

The beast among Y-haplogroups

Late PIE ground zero now obvious; location of PIE homeland still uncertain, but...

Who's your (proto) daddy Western Europeans?

The Khvalynsk men

This is where the three Samara Eneolithic or Khvalynsk samples from the recent Mathieson et al. paper plot on my Principal Component Analysis (PCA) of ancient West Eurasia. They're labeled as Steppe_CA (steppe Copper Age). I've also marked them with their Y-chromosome haplogroups.

Individual 10433, belonging to Y-chromosome haplogroup R1a, is almost a pure Eastern European Hunter-Gatherer, which is perhaps surprising, considering he was buried with copper artifacts. On the other hand, sample 10434, the one belonging to haplogroup Q1a, and positioned further east than the other two, appears to have been whacked over the head a few times and simply thrown into a ditch.

The PCA also has most of the other samples featured in Mathieson et al., including Neolithic Anatolians (labeled Anatolia_N), as well as extra samples from Allentoft et al. and Jones et al.

See also...

The Khvalynsk men #2

CHG and the Indo-European question

The recent Jones et al. palaeogenomics paper focusing on Caucasus hunter-gatherers (CHG) has this to say about the Indo-European and Indo-Aryan expansions:

CHG left their imprint on modern populations from the Caucasus and also central and south Asia possibly marking the arrival of Indo-Aryan languages.

...

It has been proposed that modern Indians are a mixture of two ancestral components, an Ancestral North Indian component related to modern West Eurasians and an Ancestral South Indian component related more distantly to the Onge [25]; here Kotias proves the majority best surrogate for the former [28,29] (Supplementary Table 10). It is estimated that this admixture in the ancestors of Indian populations occurred relatively recently, 1,900–4,200 years BP, and is possibly linked with migrations introducing Indo-European languages and Vedic religion to the region (28).

...

Finally, we found that CHG ancestry was also carried east to become a major contributor to the Ancestral North Indian component found in the Indian subcontinent. Exactly when the eastwards movement occurred is unknown, but it likely included migration around the same time as their contribution to the western European gene pool and may be linked with the spread of Indo-European languages. However, earlier movements associated with other developments such as that of cereal farming and herding are also plausible.

To their credit, in that last quote the authors leave open the possibility that CHG arrived in South Asia in multiple waves and with a variety of groups, including Neolithic farmers. Nevertheless, I'd say their comments are still confusing and perhaps also incredibly naive, because essentially they appear to be hoping that in CHG they've identified the Proto-Indo-European and Proto-Indo-Aryan genetic component.

Indeed, a lot of people actually believe that the overwhelming part of the West Eurasian admixture in South Asia should be attributed to the Indo-Aryans. But that's just stupid.

After all, many Dravidian groups that in all likelihood never spoke Indo-Aryan languages carry significant ratios of West Eurasian ancestry. Some of this influence can be explained by admixture with Indo-Aryans, but uniparental markers suggest that much of it was brought from West Asia by the Proto-Dravidians (see here).

Below is the aforementioned Supplementary Table 10. Note that two of the Indian populations that are best modeled with D-stats as mixtures of Kotias (one of the two CHG genomes) and Onge are Dravidian speakers (Mala and Vishwabrahmin or Viskwakarma, a Malayali community). Another three are Indo-Aryans (GujaratiC, GujaratiD and Lodhi), but with high levels of Ancestral South Indian (ASI) admixture, which suggests their ancestors might have been language shifters.

On the other hand, the three populations that are best modeled as Afanasievo (a pastoralist group from the Early Bronze Age steppe) and Onge are all Indo-Aryans (GujaratiA, GujaratiB and Tiwari).

But like I say, South Asia is a complex melting pot of Indo-Aryans, Dravidians, and several other linguistic groups, so a more comprehensive analysis than a comparison of a few D-stats is needed to unravel the origins of its people in a meaningful way.

By the way, Jones et al. also argue that CHG is basically an offshoot of the so called Basal Eurasian clade, which was first described in Lazaridis et al. 2014. I'm highly skeptical of this claim, and I might check it out after I get my hands on the CHG genomes.

Citation...

Jones, E. R. et al. Upper palaeolithic genomes reveal deep roots of modern eurasians. Nat. Commun. 6:8912 doi: 10.1038/ncomms9912 (2015).

See also...

The "fourth strand" of European ancestry came from the Caucasus

The "fourth strand" of European ancestry came from the Caucasus

From a news feature about a forthcoming palaeogenomics paper:

"The question of where the Yamnaya come from has been something of a mystery up to now," said one of the lead senior authors Dr Andrea Manica, from Cambridge's Department of Zoology.

"We can now answer that as we've found that their genetic make-up is a mix of Eastern European hunter-gatherers and a population from this pocket of Caucasus hunter-gatherers who weathered much of the last Ice Age in apparent isolation. This Caucasus pocket is the fourth major strand of ancient European ancestry, one that we were unaware of until now," he said

Read more at: 'Fourth strand' of European ancestry originated with hunter-gatherers isolated by Ice Age

Update: the paper is now out and open access at Nature Communications.

The two ancient Georgian genomes belong to Y-chromosome haplogroup J. So it looks like I was right when I said that this type of ancestry mostly entered the European steppe from the Caucasus via female mediated gene flow during the Bronze Age (see here and here).

However, one of the Eastern European Hunter-Gatherers (EHG) from Mathieson et al. 2015 also belonged to haplogroup J. This suggests that there was intermittent gene flow, including some paternal gene flow, between the Caucasus and the steppe well before the Bronze Age.

Nevertheless, it's now even more difficult to accept that Y-haplogroup R1 and the Proto-Indo-Europeans might have originated south of the steppe. Clearly, R1 appears to be a steppe marker from way back, and I seriously doubt that Indo-European languages were introduced into highly patriarchal steppe societies by female migrants from the Caucasus.

Image credit: Nature Communications, dx.doi.org/10.1038/ncomms9912

Mitochondrial DNA from Maykop + Wolfgang Haak on Near Eastern-related ancestry in Yamnaya

From page 166 of a report posted recently at Academia.edu:

Majkop verfügen sowohl über eine «paläolithische» Haplogruppe (U8) als auch über «neolithische» Haplogruppen: V (Недолужко u. a. 2014), T2, N1. Bei einem Objekt aus einem Grab bei der Staniza Novosvobodnaja fanden wir auch die Haplogruppe М52. Die gewonnenen Daten sprechen für eine (auf dem Niveau der mitochondrialen DNA) mögliche genetische Gemeinschaft der archäologischen Kulturen von Majkop und Novosvobodnaja.

The presence of Indian-specific mtDNA haplogroup M52 is surprising. Maykop territory was located just south of the steppe, but M52 isn't found in any of the Bronze and Iron Age samples from the steppe tested to date.

Here's the comment from Haak, from an abstract titled The role of the Caucasus in the formation of the Eurasia's genetic makeup: Insights and questions from ancient DNA research.

Recent genetic research on autosomal and uniparentally-inherited markers has shown a remarkable genetic uniformity of Caucasian populations despite the region’s notable linguistic and cultural diversity. When compared to neighbouring regions, the smooth genetic transition from the Near/Middle East to the Caucasus is in stark contrast to the marked differences to populations from the East European steppes. Flanked by the Black and the Caspian Seas, it remains unclear to what extent the Caucasus served as a corridor and whether and if so when ancient migrations had affected and shaped the region’s genetic profile. Ancient DNA research on Mesolithic, Neolithic and Bronze Age individuals from Western Eurasia have recently thrown fresh light on the Caucasus as region, which appears to have played a critical role in the formation of the genetic ancestry of the Yamnaya people, Bronze Age pastoralist of the east European steppes. The Yamnaya carry strong signals of eastern hunter-gatherer (EHG) ancestry and ancient Near Eastern ancestry that is different from the one that giving rise to early European farmers. While modern-day Armenians are the best proxy for the putative source population of the EHG dilution in the steppes, it is highly likely that prehistoric cultural groups from the Caucasus will provide a much better temporal and contextual fit.

Actually, I'd say western Georgians are the best proxy for the putative source population of the EHG dilution in the steppes. See here...

Yamnaya's exotic ancestry: The Kartvelian connection

See also...

Steppe Maykop: a buffer zone?

Genetic borders are usually linguistic borders too

On the genetic prehistory of the Greater Caucasus (Wang et al. 2018 preprint)

Plague germs may have facilitated Bronze Age expansions from the steppe

Update 07/12/2018: Europe's ancient proto-cities may have been ravaged by the plague

...

Open access at Cell:

Summary: The bacteria Yersinia pestis is the etiological agent of plague and has caused human pandemics with millions of deaths in historic times. How and when it originated remains contentious. Here, we report the oldest direct evidence of Yersinia pestis identified by ancient DNA in human teeth from Asia and Europe dating from 2,800 to 5,000 years ago. By sequencing the genomes, we find that these ancient plague strains are basal to all known Yersinia pestis. We find the origins of the Yersinia pestis lineage to be at least two times older than previous estimates. We also identify a temporal sequence of genetic changes that lead to increased virulence and the emergence of the bubonic plague. Our results show that plague infection was endemic in the human populations of Eurasia at least 3,000 years before any historical recordings of pandemics.

Rasmussen et al., Early Divergent Strains of Yersinia pestis in Eurasia 5,000 Years Ago, Cell, Volume 163, Issue 3, p571–582, 22 October 2015, DOI: https://dx.doi.org/10.1016/j.cell.2015.10.009

Also, some juicy quotes at ScienceDaily:

Study co-author Dr Marta Mirazón-Lahr, from Cambridge's Leverhulme Centre for Human Evolutionary Studies (LCHES), points out that a study earlier this year from Willerslev's Copenhagen group showed the Bronze Age to be a highly active migratory period, which could have led to the spread of pneumonic plague.

"The Bronze Age was a period of major metal weapon production, and it is thought increased warfare, which is compatible with emerging evidence of large population movements at the time. If pneumonic plague was carried as part of these migrations, it would have had devastating effects on small groups they encountered," she said.

"Well-documented cases have shown the pneumonic plague's chain of infection can go from a single hunter or herder to ravaging an entire community in two to three days."

University of Cambridge. "Plague in humans 'twice as old' but didn't begin as flea-borne, ancient DNA reveals." ScienceDaily, 22 October 2015. www.sciencedaily.com/releases/2015/10/151022124532.htm.

Basques are not simply a fusion of Iberian hunter-gatherers and early farmers

I thought I'd revisit the issue of Basque origins with my new Principal Component Analysis (PCA) of West Eurasian genetic variation. The useful thing about this PCA is that it gets around two problems that routinely affect PCA featuring ancient samples: projection bias, otherwise known as shrinkage, and exaggerated outcomes for individuals with high counts of homozygous genotypes.

A couple of recent papers argued that Basques were the direct descendants of local hunter-gatherers and early Neolithic farmers who arrived in Iberia from the eastern Mediterranean. This is probably correct for the most part, but it doesn't tell the whole story.

On the PCA above, Basques are quite distinct from Early Neolithic, Middle Neolithic and Copper Age Iberians (marked Iberia_EN, Iberia_MN and Iberia_CA, respectively), because they are significantly more eastern. In fact, they cluster with the only Bronze Age Iberian on the plot (Iberia_BA), which is the same individual that I found to harbor steppe-related ancestry (see here).

Thus, the story told by the PCA is that Basques are the progeny of Bronze Age Iberians, who, unlike their Copper Age predecessors, experienced a pulse of steppe-related admixture from the east.

Formal statistics back this up. For instance, here's a quote from the recently revised Mathieson et al. preprint:

However, the statistic f4(Basque, Iberia_Chalcolithic; Yamnaya_Samara,Chimp)=0.00168 is significantly positive (Z=8.1), as is the statistic f4(Spanish, Iberia_Chalcolithic; Yamnaya_Samara, Chimp)= 0.00092 (Z=4.6). This indicates that steppe ancestry occurs in present-day southwestern European populations, and that even the Basques cannot be considered as mixtures of early farmers and hunter-gatherers without it (4).

The key question now is who brought the steppe-related ancestry to Basque country. Were they Indo-Europeans or speakers of Proto-Basque? Also, did they actually come from the steppe, or somewhere nearby, like the Carpathian Basin?

The reason I mention the Carpathian Basin is because, as per the PCA, Basques more or less cluster between Copper Age Iberians and some of the Bronze Age Hungarians (marked Hungary_BA). But this is just one possibility, and I'm not sure at this stage how plausible it looks with, say, formal statistics.

In this analysis I used samples from the Allentoft et al., Gunther et al., Haak et al. and Lazaridis et al. datasets, all of which are publicly available. The latter two are found at the Reich Lab site here. If you're confused by some of the acronyms in the PCA key, see here.