search this blog

Friday, July 22, 2016

The Basal-rich K7


Update 25/01/2018: The Basal-rich K7 is now available to personal genomics customers for $6 USD a pop (see here).

...

I've got a new test. Currently I'm only using it to explore ancient genomes, but at some point I'll make another version available to personal genomics customers, one way or another. However, that might take a little bit of work and time to mitigate the effects of the calculator effect and so on.

Below is a spreadsheet featuring a wide range of ancient and present-day samples from recent papers. A table with the Fst genetic distances between the seven ancestral populations is available here.


Please note that the Basal-rich component is unlikely to be a perfect representation of the hypothetical Basal Eurasian population. At the same time, it's likely that the two hunter-gatherer components, Ancient North Eurasian (or AG3-related) and Villabruna-related, contain some Basal Eurasian admixture.

Here's a Principal Component Analysis (PCA) of the West Eurasian populations based on their K7 ancestry proportions. It captures all of the main features of West Eurasian genetic diversity, including the two parallel clines made up of Europeans and Near Easterners, and the intermediate position of South Central Asians between the ancient samples from Neolithic Iran and Bronze Age Europe.


An extra large version of the same PCA, with the samples labeled individually, can be downloaded here.

Also, using the K7 ancestry proportions, I modeled the ancient ancestry of a few present-day populations from the Near east, Northern Europe and South Central Asia with the nMonte R script. Bronze Age steppe admixture in groups from the latter two regions is usually inferred at 40-50% with tools based on formal stats, such as qpAdm and TreeMix, so I wanted to check if I could reproduce such results.

Iranian_Mazandarani
Iran_Chalcolithic 42.05
Iran_IA:F38 31.35
Iran_Neolithic 14.95
Andronovo_Kytmanovo 7.6
Yamnaya-Catacomb_Ulan 1.95
Han 1.35
Andamanese_Onge 0.75
Papuan 0

distance%=0.285 / distance=0.00285

Iranian_Zoroastrian
Iran_Chalcolithic 40.95
Iran_IA:F38 38.25
Yamnaya-Catacomb_Ulan 11.25
Iran_Neolithic 6.45
Andronovo_Kytmanovo 2.8
Andamanese_Onge 0.15
Han 0.15
Papuan 0

distance%=0.2979 / distance=0.002979

Kurdish
Iran_IA:F38 96.3
Andronovo_Kytmanovo 2.35
Han 0.6
Andamanese_Onge 0.55
Papuan 0.2
Iran_Chalcolithic 0
Iran_Neolithic 0
Yamnaya-Catacomb_Ulan 0

distance%=1.3845 / distance=0.013845

Latvian
Yamnaya_Peshany 44
Loschbour 30.65
Sweden_MN:Gokhem4 25.35
Barcin_Neolithic 0
Ulchi 0

distance%=1.2569 / distance=0.012569

Polish
Sweden_MN:Gokhem4 39.4
Yamnaya_Peshany 39
Loschbour 21.6
Barcin_Neolithic 0
Ulchi 0

distance%=0.5208 / distance=0.005208

Swedish
Sweden_MN:Gokhem4 42.15
Yamnaya_Peshany 36.55
Loschbour 21.3
Barcin_Neolithic 0
Ulchi 0

distance%=0.939 / distance=0.00939

Kalash
Iran_Neolithic 54.05
Yamnaya-Catacomb_Ulan 25.5
Andronovo_Kytmanovo 9.35
Han 7.6
Andamanese_Onge 1.95
Papuan 1.55

distance%=0.4592 / distance=0.004592

Pashtun_Afghanistan
Iran_Neolithic 54.8
Andronovo_Kytmanovo 31
Han 7.5
Yamnaya-Catacomb_Ulan 3.5
Andamanese_Onge 1.75
Papuan 1.45

distance%=0.4921 / distance=0.004921

Pathan
Iran_Neolithic 55.25
Yamnaya-Catacomb_Ulan 19.1
Andronovo_Kytmanovo 12.3
Han 8.45
Andamanese_Onge 3.4
Papuan 1.5

distance%=0.4848 / distance=0.004848

Tajik_Pomiri
Iran_Neolithic 41.7
Andronovo_Kytmanovo 30.65
Yamnaya-Catacomb_Ulan 19.15
Han 6.95
Andamanese_Onge 1.05
Papuan 0.5

distance%=0.5318 / distance=0.005318

Admittedly, these estimates look very conservative, but certainly not out of the ballpark. I suspect that I'll be able to improve the models and statistical fits as new Bronze Age steppe samples become available. Indeed, I'll be updating the spreadsheet above regularly.

Thursday, July 14, 2016

Early Neolithic genomes from the eastern Fertile Crescent (Broushaki et al. 2016)


Open access at Science:

Abstract: We sequenced Early Neolithic genomes from the Zagros region of Iran (eastern Fertile Crescent), where some of the earliest evidence for farming is found, and identify a previously uncharacterized population that is neither ancestral to the first European farmers nor has contributed significantly to the ancestry of modern Europeans. These people are estimated to have separated from Early Neolithic farmers in Anatolia some 46-77,000 years ago and show affinities to modern day Pakistani and Afghan populations, but particularly to Iranian Zoroastrians. We conclude that multiple, genetically differentiated hunter-gatherer populations adopted farming in SW-Asia, that components of pre-Neolithic population structure were preserved as farming spread into neighboring regions, and that the Zagros region was the cradle of eastward expansion.


Broushaki et al., Early Neolithic genomes from the eastern Fertile Crescent, Science 14 Jul 2016, DOI: 10.1126/science.aaf7943

See also...

Economic overhaul + population shift in Late Neolithic Iran

Modeling Steppe_EMBA

Yamnaya =/= Eastern Hunter-Gatherers + Iran Chalcolithic

Monday, July 11, 2016

Genome-wide variants of Eurasian facial shape differentiation


Very interesting new preprint at bioRxiv:

Abstract: It is a long standing question as which genes define the characteristic facial features among different ethnic groups. In this study, we use Uyghurs, an ancient admixed population to query the genetic bases why Europeans and Han Chinese look different. Facial trait variations were analyzed based on high dense 3D facial images; numerous biometric spaces were examined for divergent facial features between European and Han Chinese, ranging from inner-landmarks to dense shape geometrics. A series of genome-wide association analyses were conducted on a discovery panel of Uyghurs. Six significant loci were identified and four of which, rs1868752, rs118078182, rs60159418 at or near UBASH3B, COL23A1, PCDH7 and rs17868256 were replicated in two independent cohorts of Uyghurs or Southern Han Chinese. We further developed a quantitative model to predict 3D faces based on 277 top GWAS SNPs. In hypothetic forensic scenarios, this model was found to significantly enhance the rate of suspect verification, suggesting a practical potential of related research.

Lu Qiao et al., Detecting Genome-wide Variants of Eurasian Facial Shape Differentiation: DNA based Face Prediction Tested in Forensic Scenario, bioRxiv, posted July 11, 2016, doi: http://dx.doi.org/10.1101/062950

Layers of Ancient North Eurasian-related ancestry in East Asia


Back in May I hypothesized that present-day East Asians were prehistoric hybrids of partly Ancient North Eurasian (ANE) origin. I got the idea from a series of TreeMix runs (see here).

This was essentially confirmed recently in the Lazaridis et al. 2016 preprint. Refer to page 147 in the paper's supplementary information PDF here.

However, based on more recent TreeMix runs featuring data from Lazaridis et al., I'd say the situation is more complex than just some minor ANE-related admixture in East Asians. I suspect now that all East Asians, including even the Onge, an ancient isolate population from the Andaman Islands, harbor significant ANE-related ancestry that may have arrived in East Asia in separate waves.

Here's what I'm talking about. Note that all of the samples on the East Asian node - Upper Paleolithic west Siberian forager Ust-Ishim, Han Chinese and Onge - are influenced by a massive migration edge from the base of the AG3-MA1 or ANE branch. However, as per the second graph, only the ancestors of more northerly East Asians, like those of the Han, appear to have been recipients of the latest ANE-related admixture into East Asia.



Indeed, when I add the Natufians from the Epipaleolithic Levant to the analysis, Ust-Ishim and the East Asians join AG3-MA1 on the same branch, but now receive a 36% migration edge from a point basal to all Eurasians. This is not admixture from the hypothesized Basal Eurasian clade, but probably from another basal clade, specific to East Asians, which I'd say occasionally shows up as pseudo Sub-Saharan admixture in East Asians.


But obviously, we'll need a solid selection of ancient genomes from across space and time in East Asia to confirm these results. Rumor has it that they're on their way.

See also...

East and West Eurasians seperated at least 45,000 years ago, but...

Saturday, July 9, 2016

Modeling Steppe_EMBA


Lazaridis et al. showed that their Steppe_EMBA grouping, which included Afanasievo, Poltavka and Yamnaya, as well as two Potapovka samples, one Russia_EBA sample and one Srubnaya_outlier sample, were best modeled in the following two ways using qpAdm:

Steppe_EMBA
Eatern Hunter-Gatherer (EHG) 0.568
Iran Chalcolithic (Iran_ChL) 0.432

Steppe_EMBA
Caucasus Hunter-Gatherer (CHG) 0.181
Eastern Hunter-Gatherer (EHG) 0.527
Iran Chalcolithic (Iran_ChL) 0.292

I'm not a huge fan of either of these models, but especially the first one, even though I understand that they're both statistically very sound. For one, the uniparental markers don't match, and two, TreeMix seems to disagree (see here).


So let's try something a little different and see what happens when I model Steppe_EMBA as EHG, CHG, and Anatolia Chalcolithic.

Outgroups
Anatolia_Neolithic
Andamanese_Onge
Chukchi
Han
Israel_Natufian
Karitiana
Kostenki14
Levant_Neolithic
MA1
Mbuti.DG
Papuan
WHG

Steppe_EMBA
Anatolia Chalcolithic (Anatolia_ChL) 0.128
Caucasus Hunter-Gatherer (CHG) 0.375
Eastern Hunter-Gatherer (EHG) 0.497

As far as I can tell, it's a very decent fit, especially considering that I'm using 12 outgroups and three reference populations. To me, at least, the standard errors look surprisingly low for such a complex model: 0.033, 0.046 and 0.020, respectively.

Now, I'm not arguing here that Chalcolithic Anatolia is the answer. What I'm saying is that multiple lines of evidence do not support Chalcolithic Iran as a real source of admixture for Steppe_EMBA, and I'm offering what I see as a plausible alternative among the currently available samples.

I know that this is a work in progress for the Broad MIT/Harvard team, and we'll have to wait for more ancient samples and another paper or two before a consensus is reached on the topic.

But here's my prediction: Steppe_EMBA only has 10-15% admixture from the post-Mesolithic Near East not including the North Caucasus, and basically all of this comes via female mediated gene flow from farming communities in the Caucasus and perhaps present-day Ukraine.

See also...

Another look at the genetic structure of Yamnaya

Friday, July 8, 2016

Khazar shmazar #2


Open access at Genome Biology and Evolution:

Abstract: In a recent interdisciplinary study, Das and co-authors have attempted to trace the homeland of Ashkenazi Jews and of their historical language, Yiddish (Das et al. 2016. Localizing Ashkenazic Jews to Primeval Villages in the Ancient Iranian Lands of Ashkenaz. Genome Biology and Evolution). Das and co-authors applied the geographic population structure (GPS) method to autosomal genotyping data and inferred geographic coordinates of populations supposedly ancestral to Ashkenazi Jews, placing them in Eastern Turkey. They argued that this unexpected genetic result goes against the widely accepted notion of Ashkenazi origin in the Levant, and speculated that Yiddish was originally a Slavic language strongly influenced by Iranian and Turkic languages, and later remodeled completely under Germanic influence. In our view, there are major conceptual problems with both the genetic and linguistic parts of the work. We argue that GPS is a provenancing tool suited to inferring the geographic region where a modern and recently unadmixed genome is most likely to arise, but is hardly suitable for admixed populations and for tracing ancestry up to 1000 years before present, as its authors have previously claimed. Moreover, all methods of historical linguistics concur that Yiddish is a Germanic language, with no reliable evidence for Slavic, Iranian, or Turkic substrata.

Flegontov et al., Pitfalls of the geographic population structure (GPS) approach applied to human genetic history: A case study of Ashkenazi Jews, Genome Biol Evol (2016) doi: 10.1093/gbe/evw162

See also...

Khazar shmazar

Irano-Turko-Slavic roots of Ashkenazi Jews?

Monday, July 4, 2016

Economic overhaul + population shift in Late Neolithic Iran


Courtesy of Arbuckle et al. at the Journal of Archaeological Science. Emphasis is mine:

Abstract: In this paper we address the timing of and mechanisms for the appearance of domestic cattle in the Eastern Fertile Crescent (EFC) region of SW Asia through the analysis of new and previously published species abundance and biometric data from 86 archaeofaunal assemblages. We find that Bos exploitation was a minor component of animal economies in the EFC in the late Pleistocene and early Holocene but increased dramatically in the sixth millennium BC. Moreover, biometric data indicate that small sized Bos, likely representing domesticates, appear suddenly in the region without any transitional forms in the early to mid sixth millennium BC. This suggests that domestic cattle were imported into the EFC, possibly associated with the spread of the Halaf archaeological culture, several millennia after they first appear in the neighboring northern Levant.

These findings more or less correlate with the results in the new Lazaridis et al. preprint:

During subsequent millennia, the early farmer populations of the Near East expanded in all directions and mixed, as we can only model populations of the Chalcolithic and subsequent Bronze Age as having ancestry from two or more sources. The Chalcolithic people of western Iran can be modelled as a mixture of the Neolithic people of western Iran, the Levant, and Caucasus Hunter Gatherers (CHG), consistent with their position in the PCA (Fig. 1b).

In other words, the small cows weren't just imported into the Eastern Fertile Crescent; they came with people who also made a major genetic impact on the region.

Here's my own PCA featuring the relevant Lazaridis et al. samples. Key: Caucasus_HG = Caucasus Hunter-Gatherer; Iran_ChL = Iran Chalcolithic; Iran_HG = Iran Hunter-Gatherer; Iran_N = Iran Neolithic; Levant_N = Levant Neolithic.


Thus, it would seem that after the early Neolithic farmers from Iran migrated to South Asia, they were largely replaced in their own homeland by Halaf pastoralists and/or related groups. Moreover, their descendents in South Asia, and especially South Central Asia, were then largely replaced by pastoralists from the Bronze Age Eurasian steppe (for instance, see here).

Obviously, this doesn't square too well with the idea of a Proto-Indo-European homeland in the Zagros Mountains of western Iran, does it?

See also...

Yamnaya =/= Eastern Hunter-Gatherers + Iran Chalcolithic

Zarathushtra and his steppe posse