search this blog

Wednesday, June 20, 2012

First direct evidence of genetic continuity in West and Central Poland from the Iron Age to the present


I've just been sent a fascinating thesis on the mtDNA of Iron Age and Medieval samples from Poland. It suggests direct genetic continuity between Iron Age samples belonging to the Przeworsk and Wielbark Cultures, of what is now West and Central Poland, and present-day Poles. Here's the English summary, and a map of the sites under study:

For many years the origin of the Slavs has been the subject-matter in archaeology, anthropology, history, linguistics and recently also modern human population genetics. By now there is no unambiguous answer to a question where, when and in what way the Slavs originated. For the purposes of this dissertation, the analysis of ancient human mitochondrial DNA was applied. The ancient DNA was isolated from 72 specimens which came from Iron-Age and medieval graveyards from the area of current Poland. Ancient mtDNA was extracted from two teeth from each individual and reproducible sequence results were obtained for 20 medieval and 23 Iron-Age specimens. On the basis of HVR I mtDNA mutation motifs and coding region SNPs each specimen was assigned to a mitochondrial haplogroup. The obtained results were used together with other ancient and modern populations to analyse shared haplotypes and population genetic distances illustrated by multidimentional scaling plots (MDS). The differences on genetic level and quite high genetic distances (FST) between medieval and Iron-Age populations as well as significant number of shared informative haplotypes with Belarus, Ukraine and Bulgaria may evidence genetic discontinuity between medieval and Iron Ages. From the other side, the highest number of shared informative haplotypes between Iron-Age and extant Polish population as well as the presence of subhaplogroup N1a1a2, can confirm that some genetic lines show continuity at least from Iron Age or even Neolithic in the areas of present day Poland. The results obtained in this work are considered to be the first ancient contribution in genetic history of the Slavs.


Below is an MDS from the thesis, based on data corrected for the effects of potential relatives in the Iron Age sample. I don't think it's a particularly useful way of judging the intra-European affinity of the two ancient Polish groups, mostly because the samples are small, and contemporary North, Central and East Europeans don't differ very much in terms of mtDNA. Nevertheless, we can see that both the Iron Age (Okres Rzymski) and Medieval (Sredniowiecze) samples fall within the range of modern European mtDNA diversity. On the other hand, the German Neolithic LBK sample (Neolit LBK Niemcy) clearly does not, because it's sitting at the far right of the plot, away from the main European cluster. This dichotomy between the genetic structure of the LBK farmers and modern Europeans has been demonstrated in previous studies, but the reasons for it are still a mystery.



Interestingly, modern Poles are closer to an Iron Age sample from Denmark (Okres Zelaza Dania) than to the Polish Iron Age set. However, as per the summary above, the author also compared the frequencies of the most informative haplotypes among the modern and ancient samples, and found that extant Poles are the closest group to the Polish Iron Age remains, followed by Balts, Swedes and Baltic Finns. Below is a table showing those results.




According to the author, these matches might hint at Baltic, Germanic and Finno-Ugric influences in the Polish Iron Age population. Perhaps, but in my opinion, they're simply in line with geography, and reflect the general North European character of maternal lineages shared by populations from around the Baltic, both today and during the Iron Age.

The results for the Medieval Polish sample are more intriguing, because they're somewhat out of whack with geography. Its best matching modern groups are Belorussians, Ukrainians and Bulgarians. This might suggest that, during the early middle ages, the territory of present day Poland experienced an influx of groups from what are now Belarus and Ukraine, who then melted into the gene pool of the natives of Polish Iron Age descent. However, conversely, it might mean that Belorussians, Ukrainians and Bulgarians descend in large part from fairly specific medieval groups from the area of modern Poland.




In any case, whether present day Polish territory saw some migrations from the immediate east during the Medieval period or not, this preliminary look at ancient Polish mtDNA suggests long-standing genetic continuity in the region. What it clearly doesn't show is a complete, or almost complete, population replacement in the areas between the Oder and Bug rivers during the migration period.

Indeed, the thesis results put into doubt past notions that the Przeworsk and Wielbark cultures were of Germanic origin.

The (mtDNA) haplogroup missing from both the Iron Age and medieval samples from the territory of modern Poland was haplogroup I. In contemporary Slavic populations, this haplogroup is found at levels ranging from 1.2% in Bulgarians to 4.8% in Slovaks. It was also recorded at high levels in ancient remains from Denmark. It showed a frequency of 12.5% in an Iron Age sample, and 13.8% in a medieval sample. Melchior et al. 2008 suggest that haplogroup I might have been more common in Denmark and Northern Europe during that period. Therefore, the lack of this haplogroup in ancient DNA from the territory of modern Poland, might mean that the Przeworsk and Wielbark cultures should not be identified with Germanic populations.

I'm sure more ancient DNA studies are on the way looking at the origins of Slavs and Poles. Indeed, if the Y-chromosomes of Przeworsk and Wielbark remains are successfully tested, I won't be surprised if they look fairly typical of modern Poles, with a decent representation of R1a1a-M458, which is the most common Y-chromosome haplogroup in Poland today.

Anna Juras, Etnogeneza Słowian w świetle badań kopalnego DNA, Praca doktorska wykonana w Zakładzie Biologii Ewolucyjnej Człowieka Instytutu Antropologii UAM w Poznaniu pod kierunkiem Prof. dr hab. Janusza Piontka


Friday, June 15, 2012

Ancient mtDNA from the Dnieper-Donets cultural complex


A new paper at the Journal of Human Genetics reports on the mtDNA gene pool of the Dnieper-Donets (DD) cultural complex of Neolithic Ukraine. The authors were able to confirm the presence of the following haplogroups in the ancient remains from the Mariupol-type sites along the Dnieper: two H, two C, one C4a2, one U5a1a and one U3. So, three out of the seven samples belonged to haplogroup C, which is a Siberian-specific marker.

We've already had a sneak peek at these results thanks to a thesis abstract published last year by the Grand Valley State University (see here). The presence of mtDNA C among the DD remains suggests that the North Pontic steppe was formerly inhabited by genetically heterogeneous groups, and part of their ancestry came from Siberia.

The fact that Siberian-specific mtDNA lineages are very rare in Ukraine today, means that something must have happened there since the Bronze Age that basically wiped them out.


Alexey G Nikitin et al., Mitochondrial haplogroup C in ancient mitochondrial DNA from Ukraine extends the presence of East Eurasian genetic lineages in Neolithic Central and Eastern Europe, Journal of Human Genetics, advance online publication, 7 June 2012; doi:10.1038/jhg.2012.69


Saturday, April 21, 2012

So who's the most (indigenous) European of us all?


Basically, the first map below reveals the answer. It shows the spread of a European specific cluster from a global-wide ADMIXTURE analysis at K=8 (eight ancestral populations assumed), which I call "North European". Thus, genetically, the most European populations are found around the Baltic Sea, and in particular in the East Baltic region. In my genome collection, samples from Lithuania clearly and consistently score the highest percentages in ADMIXTURE clusters specific to Europe. However, I suspect that if I had Latvians with no known foreign ancestry going back more than four generations, they'd come out the "most European". Hopefully we can test that in the near future.


Below are the fifteen Eurogenes sample sets that scored the highest levels of membership in the North European cluster. The list only includes groups with five or more individuals present in the analysis, so some populations, like Estonians or Danes, weren't included, even though they easily made the cut. The spreadsheet with all the results from this run can be seen here. A table of Fst (genetic) distances between the eight clusters is available here.

Lithuanians 77%
Finns 74%
Belorussians 70%
Swedes 69%
Norwegians 68%
Kargopol Russians 68%
Russians 68%
Poles 68%
Erzya 66%
Ukrainians 66%
Moksha 66%
Orcadians 63%
HapMap Utah Americans (CEU) 63%
Irish 63%
British 62%

So why did I pick the results from K=8, and not some other K, like 2, 10, or 25? Well, it's not possible to evaluate who is more European without a European-specific cluster (ie. modal in Europeans, with a low frequency outside of Europe). Provided that a decent number and range of global and West Eurasian samples are used in the analysis, such clusters begin appearing at around K=5 or K=6, and start breaking up into local clusters from about K=9. I found that runs below K=8 produced European clusters that spilled too generously outside of the borders of Europe. On the other hand, runs above K=8 produced European clusters that weren't representative of enough European groups (ie. too localized). But the European cluster from K=8 was pretty much perfect, and I think that's obvious from the map. In fact, I can hardly believe how well it fits the modern geographic concept of Europe - north of the Mediterranean and west of the Urals. Amazing stuff.

There are two other clusters that show up across Europe in non-trivial amounts - Mediterranean and Caucasus (see maps below). These can also be thought of as native European clusters, since they've been on the continent for thousands of years. However, their peak frequencies are found in West Asia, so they're not particularly useful signals of European-specific ancestry.



So what do these three clusters show exactly? They represent certain allele frequencies in modern populations, and in fact, these can change fairly rapidly due to admixture, selection, and genetic drift. So claiming that such clusters represent pure ancient populations is unlikely to be true in most cases, if ever. However, I don't think there's anything wrong in saying that, when robust enough, they can be thought of as signals of ancestry from relatively distinct ancestral groups.

Indeed, anyone who's read up on the prehistory of Europe, knows that there are three general Neolithic archeological waves to consider when trying to untangle the story of the peopling of Europe. These are Mediterranean Neolithic, Anatolian Neolithic and Forest Neolithic (for example, see here).

Mediterranean Neolithic refers to a series of migrations from West Asia via the Mediterranean and its coasts. The areas most profoundly affected by these movements include the islands of Sardinia and Corsica, and the Southwest European mainland. Anatolian Neolithic describes migrations into Europe from modern day Turkey, mostly into the Balkans, but also as far as Germany and France. At the moment, Forest Neolithic of Northeastern Europe is something of a mystery. However, the general opinion is that it was largely the result of native Mesolithic hunter-gatherers adopting agriculture.

Obviously, it's very difficult to dismiss the correlations between these three broad archeological groups and the European and two European/West Asian clusters produced in my K=8 ADMIXTURE analysis. Is it a coincidence that the Mediterranean cluster today peaks in Sardinia, which has been largely shielded from foreign admixture since the Neolithic, and today forms a very distinct Southern European isolate? Why does the North European cluster show the highest peaks in classic Forest Neolithic territory? And why does the Caucasus cluster radiate in Europe from the southeast, which is where Anatolian farmers had the greatest impact? These can't all be coincidences, and I'm willing to bet that none of them are. I'm convinced that the three clusters from my K=8 run are strong signals from the Neolithic, and the North European cluster also from the Mesolithic.

Eventually, these issues will be settled with ancient DNA data, in a much more comprehensive way than ever possible using modern genomes. We've already seen some preliminary results, mostly from Mesolithic, Neolithic and Bronze Age sites around Europe, so perhaps it's useful to ask whether my ADMIXTURE analysis and commentary here mirror these early findings? I think they do. For instance, here's an interesting conclusion regarding the East Baltic area from a study on ancient Scandinavian mtDNA by Malmström et al.

Through analysis of DNA extracted from ancient Scandinavian human remains, we show that people of the Pitted Ware culture were not the direct ancestors of modern Scandinavians (including the Saami people of northern Scandinavia) but are more closely related to contemporary populations of the eastern Baltic region. Our findings support hypotheses arising from archaeological analyses that propose a Neolithic or post-Neolithic population replacement in Scandinavia [7]. Furthermore, our data are consistent with the view that the eastern Baltic represents a genetic refugia for some of the European hunter-gatherer populations.

I suppose there will be people wondering why I didn't take Sub-Saharan African, East Asian, and South Asian admixtures into account in my analysis. The reason is that I wasn't looking at which group was most West Eurasian, or Caucasoid. Based on everything I've seen to date, in my own work as well as elsewhere, the most West Eurasian group would probably be the French Basques from the HGDP. However, the differences between them, and certain groups from Northeastern Europe, like Northern Poles and Lithuanians, really wouldn't be that great anyway. I might do a write up about that at some point.


Credits...

- Maps by Eurogenes project member FR7

- Additional stats by Eurogenes project member DESEUK1


References...

Helena Malmström et al., Ancient DNA Reveals Lack of Continuity between Neolithic Hunter-Gatherers and Contemporary Scandinavians, Current Biology, 24 September 2009, doi:10.1016/j.cub.2009.09.017

Noreen von Cramon-Taubadel and Ron Pinhasi, Craniometric data support a mosaic model of demic and cultural Neolithic diffusion to outlying regions of Europe, Proc. R. Soc. B published online 23 February 2011, doi: 10.1098/rspb.2010.2678


Thursday, March 29, 2012

What's the point of peer review (in reference to Marc Haber et al. 2012)?


I was just reading a new study in PLoS One about Y-DNA in Afghanistan, and stopped when I saw this idiotic paragraph...

R1a1a-M17 diversity declines toward the Pontic-Caspian steppe where the mid-Holocene R1a1a7-M458 sublineage is dominant [46]. R1a1a7-M458 was absent in Afghanistan, suggesting that R1a1a-M17 does not support, as previously thought [47], expansions from the Pontic Steppe [3], bringing the Indo-European languages to Central Asia and India.

OK, first of all, the authors made a glaring error in claiming that R1a1a7-M458 dominates the Pontic-Caspian Steppe. It most certainly does not, and this information is available in the paper they referenced (Underhill et al. 2009). Secondly, who cares about R1a1a STR diversity? It's not relevant, because it tells us nothing about the origins of R1a1a. In fact, there's no way anyone can accurately estimate the ages/expansion times of Y-chromosome haplogroups. Scientists have attempted such feats on many occasions in recent years, and often came up with ridiculous results. So perhaps it's now time to admit there's a problem and move on?

I'd say there's no reason why R1a1a7-M458 should be present in Central Asia and India. Simply by looking at its modern distribution and frequencies, without even attempting any complex calculations, it seems to have expanded around East Central Europe well after the early Indo-European dispersals (see here). The most sensible claim anyone can make about R1a1a7-M458 is that it's a Slavic, or even West Slavic marker, that probably originated in the people who would become Slavs in or near modern-day Poland.

The European-specific R1a1a SNP that scientists should be looking for in Central Asia, in order to track the movements of the early Indo-Europeans, is Z280. This marker has a much wider distribution in Central and Eastern Europe than M458, and perhaps that was also the case during the relevant time frames - the Chalcolithic and early Bronze Age? Indeed, it has already been found in native Central Asian samples, both in private and academic tests, and the latter results will hopefully be published soon.

To widen the net, they should also test for Z283, which is upstream to both M458 and Z280. Its present distribution hints that it might have been a common marker within the late Neolithic Corded Ware cultural horizon of the North European Plain, which is usually thought of as an early Indo-European culture. If that was the case, then in theory, based on archeological data, it might have traveled with representatives of the Eastern Corded Ware, the Abashevo culture, past the Urals and as far as East Central Asia.

Haber M, Platt DE, Ashrafian Bonab M, Youhanna SC, Soria-Hernanz DF, et al., (2012) Afghanistan's Ethnic Groups Share a Y-Chromosomal Heritage Structured by Historical Events. PLoS ONE 7(3):e34288. doi:10.1371/journal.pone.0034288



See also...

First official attempt to divide R1a1 into multiple subclades since the discovery of R-M458


Wednesday, March 7, 2012

Southwest Eurasians + Northwest Eurasians + Mesolithic survivors = modern Europeans


Update 23/13/2013: Ancient human genomes suggest (more than) three ancestral populations for present-day Europeans

...

For a long time, it was generally accepted that Europeans were direct descendants of Palaeolithic settlers of the continent, with some Middle Eastern ancestry in the Mediterranean regions, courtesy of Neolithic farmers. However, in the last few years, largely thanks to ancient DNA, it dawned on most people that such a scenario was unrealistic. It now seems that Europe was populated after the Ice Age in a big way, by multiple waves of migrants from almost all directions, but especially from the southeast.

Getting to grips with the finer details of the peopling of Europe is going to be a difficult and painstaking process, and will require ancient DNA technology that probably isn’t even available at the moment. However, the mystery about the basic origins and genetic structure of Europeans was solved for me this week, after I completed a series of ADMIXTURE runs focusing on West Eurasia (see
K=10, K=11, K=12, and K=13). The map below, produced by one of my project members, surmises very nicely the most pertinent information from those runs (thanks FR7!). It shows the relative spread of three key genetic clusters, from the K=13, in a wide range of populations from Europe, North Africa, and West, Central and South Asia (i.e. the data represents the nature of West Eurasian alleles in the sampled groups, with only three clusters considered). The yellow cluster is best described as Mediterranean or Southwest Eurasian, while the cyan and magenta, which are sister clades, as Northwest Eurasian.





Thus, it appears as if modern Europeans are made up of two major Neolithic groups, which are related, but at some point became distinct enough to leave persistent signals of that split. They spread into different parts of Western Asia before moving into Europe. The Southwest Eurasians, possibly from the southern Levant, dominated the Mediterranean Basin, including North Africa, Southern Europe, and the Arabian Peninsula. I’m pretty sure that Otzi the Iceman is the best known representative of the ancient Southwest Eurasians (see here).

The Northwest Eurasians might have originated in the northern Levant, but that’s a pure guess. In fact, judging by the map above, their influence isn’t particularly strong in that part of the world today, and only becomes noticeable several hundred kilometers to the north and east, in the North Caucasus and Iran respectively. However, the northern Levant is actually dominated by a fourth West Eurasian cluster, tagged by me as "Caucasus" in the K=13 run, and not shown on the map above. Various calculations show that this can also be assigned to the Northwest Eurasian group, except that it seems to have split from the other Northwest Eurasian components at an early stage (see comments section here).

After their initial spread, it appears as if the Northwest Eurasians inhaled varying amounts of native Mesolithic groups in their newly acquired territories west, north and east of the Levant. This is being strongly suggested by the aforementioned ancient DNA results, at least as far as Europe is concerned. They also mixed heavily with Southwest Eurasians in Europe and nearby. That’s why, for instance, you’ll never find an Irishman who clusters closer genetically to an Indian than to other Europeans. However, even a basic analysis of their DNA, like my own ADMIXTURE runs, shows that a large subset of their genomes comes from the same, relatively recent, “Northwest Eurasian” source.


We can follow the same logic when talking about the differentiation between modern descendants of Southwest Eurasians. For instance, those in Iberia have significant admixture from Northwest Eurasians, while those in North Africa carry appreciable amounts of West and East African influence.

I’m convinced that the scenario of the peopling of Europe outlined above, by two basic stocks of migrants from Neolithic West Asia, is the only plausible one, because the signals from the data are too strong to argue against it. I’m sure you’ll be seeing the same story told by scientists over the next few years in peer reviewed papers. They’ll probably come up with different monikers for the Southwest and Northwest Eurasians, but the general concepts will be the same.

However, that was the easy part. The hard part is linking the myriad of movements of these Southwest and Northwest Eurasians with archaeological and linguistic groups. Perhaps the earliest Southwest Eurasians into Europe were Afro-Asiatic speakers? To be honest, I have no idea, because that’s not an area I’ve studied closely. But I would say that it’s almost certain that the proto-Indo-Europeans were of Northwest Eurasian stock. It’s an obvious conclusion, due to the trivial to non-existent amounts of Southwest Eurasian influence in regions associated with the early Indo-Europeans, like Eastern Europe and Central Asia.

Perhaps the simplest and most diplomatic thing to do for the time being, would be to associate the entire Northwest Eurasian group with an early (Neolithic) spread of Indo-European languages from somewhere on the border between West Asia and Europe? I know that would work for a lot of people, specifically those who’d like to see an Indo-European urheimat in Asia, as opposed to Europe. But it wouldn’t work for me, especially not after taking a closer look at that map above.

As already mentioned, the Northwest Eurasians can be reliably split into two clusters, marked on the map in cyan and magenta. I call the cyan cluster North Atlantic, because it peaks among the Irish and other Atlantic fringe groups, and the magenta Baltic, because it shows the highest frequencies among Lithuanians and nearby populations. The story suggested by the map is pretty awesome, with the Baltic cluster seemingly exploding from somewhere in the middle of the Northwest Eurasian range, and pushing its close relatives to the peripheries of that range. Thus, under such a dramatic model, the North Atlantic is essentially the remnant of the pre-Baltic Northwest Eurasians, and appears to have found refuge in Western and Northwestern Europe, in the valleys of the Caucasus Mountains, and in South Asia.

Indeed, there seems to be a correlation between the highest relative frequencies of the North Atlantic and regions that are still home to non-Indo-European speakers, or were known to have been home to such groups in historic times. For instance, France has the Basques, while the British Isles had the Picts, who are hypothesized to be of non-Indo-European stock. Note also the native, non-Indo European speakers in the Caucasus, like the Chechens, who show extreme relative frequencies of the North Atlantic component. Moreover, at the south-eastern end of the Northwest Eurasian range, in India, there are still many groups of Dravidian speakers.

Below are two maps that isolate the relative frequencies of the North Atlantic (cyan) and Baltic (magenta) components, versus each other and the Southwest Eurasian cluster, to better show the hole in the distribution of the North Atlantic. To be sure, this North Atlantic can be broken down further, but only with more a comprehensive sampling strategy, especially of Northern and Western Europe.






That’s my take on what the data is showing, and other explanations are possible. But I don’t really know what they might be? I should also mention that the potentially proto-Indo-European Baltic cluster shows a remarkable correlation with the spread of Y-chromosome haplogroup R1a, and ancient DNA rich in this haplogroup from supposed early Indo-Europeans. For more info on that, see the links below:


Best of 2008: Corded Ware DNA from Germany


Ancient Siberians carrying R1a1 had light eyes


Ancient Siberians carrying R1a1 had light eyes - take 2

Bronze Age Tarim Basin "Caucasoids" carried R1a1 (and European mtDNA lineages too)

European admixture among ancient East Asians (aka. two-rooted canines carried by early Indo-Europeans to China)


Sunday, February 26, 2012

Genetic substructures within the HapMap CEU sample (and Eurogenes' Northwest Europeans)


In this experiment I attempt to characterize more precisely the origins of some of the individuals from the HapMap CEU cohort. These samples are described by the HapMap project as Utah Americans of Western and Northern European descent. But this doesn't seem to be exactly true for at least two of them, who actually come out very Central European in all my tests. Moreover, it's obvious that some of the samples fit nicely into very specific areas of Western and Northern Europe. For instance, at this level of resolution, a few could pass as Irish, and others for Danes or even Swedes. Below is a quick and dirty ADMIXTURE analysis designed specifically for this experiment.



Key: Red = Sub-Saharan African, Yellow = Southern European, Green = North-Central European, Aqua = North Atlantic, Blue = Baltic, Pink = East Asian. See spreadsheet for details.


Based on the K=6 results it's fair to say that at least six of the CEU samples might pass for unmixed Scandinavians, most likely Danes or southern Swedes (NA12003, NA12057, NA12248, NA12249, NA12776 and NA12875). At least five could be confused for Irish or western British samples (NA10850, NA12005, NA12006, NA12386 and NA12812). The two Central European-like Utahns stick out from the CEU set due to their unusually high Baltic scores (NA11917 and NA12286). From the little I know about the CEU samples, I'd say that these two were of eastern or southeastern German origin. But they might have fairly recent ancestry from further east than that. My own MDS analysis (first image below) and a PCA plot from Lao et al. 2008 (second image, slightly edited by me to remove article text) confirm that such Scandinavian-like, German-like and Irish-like individuals do exist in the CEU set.



So the CEU set is not a homogeneous one, and carries clear substructures that can be picked up via fairly basic means. However, this doesn't make the CEU samples less valuable, but more so, due to the lack of public access to continental Northwestern European samples. Secondly, the test reveals some interesting information about the genetic substructures within Northwestern Europe. Here are some of my observations:

- Scandinavians often show very high levels of the North-Central European component, and moderately high levels of the North Atlantic component. Many also carry clear amounts of the Baltic component, but, as a rule, lower levels of the Southern European component.

- Germans mainly differ from the Scandinavians in that they carry the Southern European component at appreciable amounts. They show variable amounts of the Baltic component, with those from eastern Germany carrying the highest levels.

- Irish project members, especially those from western Ireland, show very high levels of the North Atlantic component, but low levels of the Southern European component.

- Western British samples, like those from Cornwall or western Scotland, are generally very similar to the Irish, mainly in that they carry the North Atlantic component at high levels. However, they often show somewhat higher levels of the Southern European component.

I'm eventually going to test these classifications of the CEU samples with ChromoPainter, which is by far the most accurate tool for such things at the moment. Unfortunately, it's also a lot of hard work and computationally intensive, so it might take a few weeks. I do have the allele frequencies from the above ADMIXTURE run, and it is possible to make a stand alone test from them. However, I'm not certain that's a good idea at present, due to the small number of samples involved. It might be worth doing when the right samples swell in number, so I can run a more robust analysis. In particular, I need more people from Ireland, Scotland and Scandinavia.

Reference...

Oscar Lao et al, Correlation between Genetic and Geographic Structure in Europe, Current Biology, Volume 18, Issue 16, 1241-1248, 26 August 2008, doi:10.1016/j.cub.2008.07.049


Sunday, January 22, 2012

Eurogenes' North Euro clusters - phase 2, final results


This is a continuation of my ChromoPainter analysis of Europeans from north of the Pyrenees, Alps and Balkans (see here). To obtain the most accurate results possible on my laptop, I increased the burn-ins and iterations in fineSTRUCTURE to 500K each (5 hour run in all, which is all I'm willing to put this machine through). The end product looks very similar to my initial analysis, in which I explored the data at 200K burn-ins and iterations. What I think this shows is that the results are robust, and I doubt they'd change much even after a couple of days of running fineSTRUCTURE.

Indeed, as mentioned in my previous blog entry, this appears to be the most detailed and accurate cluster analysis of this part of Europe produced anywhere to date. There are 21 clusters in all, with at least 20 looking like strong signals of genetic substructures across North, West, Central and East Europe (see spreadsheet for individual classifications). They include:

pop0 - West Finnish1: This is a pair of reference individuals, most likely from Western Finland, judging by their PCA and ADMIXTURE results. They are either from the same community, or have a very similar mix of very specific ancestries.

pop1 - Erzya + Moksha: This includes all of the Erzya and Moksha in the project, plus a Russian with recent Erzya ancestry. It's closely related to ethnic Russian clusters that stretch from Northwest Russia to near the Volga, and also to the Estonian cluster.

pop2 - South/Central Finnish: This is the largest Finnish cluster, and that's probably more than just the result of sampling bias. I would say that the greater part of the Finnish population would belong to this type of cluster, which occupies regions of highest population density within the country.

pop3 - Fenno-Scandian: This cluster includes a Northern Swede, a Swede with probable recent Finnish ancestry, and Finns with probable recent Swedish influence. I have a feeling that Finland Swedes and Aland Islanders would also be placed here more often than not.

pop4 - Northwest Russian/Southeast Finnish: Although this cluster includes only two individuals, it's definitely much more than just the result of two relatively closely related samples being in the same run. I'd hazard a guess that Northwest Russians with, say, significant Ingrian ancestry, would land here, and so would Finns with recent Russian ancestry.

pop5 - West Finnish2: Based on PCA and ADMIXTURE results, most of these Finns likely come from Western Finland, probably from places like Southern Ostrobothnia. They possibly also have some Swedish influence.

pop6 - West German: This cluster is based on individuals from Western and Northwestern Germany. It also includes a Dutchman, Austrian and people of mixed origin, like a Dane with French and German ancestry, and Americans with British, German, Scandinavian and/or Polish ancestry. In other words, this is where Northwestern Europe meets Central Europe.

pop7 - Vologda Russian: Most of the Vologda Russians from the HGDP land here, so this appears to be a local cluster. Judging from its phylogeny, it looks like a mix of North Slavic, Baltic and Finnic influences.

pop8 - East Finnish: All the project and reference Finns with substantial ancestry from new settlement areas of Eastern Finland appear in this cluster. No wonder then, that this is the cluster with the highest chunk count in this analysis.

pop9 - Estonian: This is a mixed cluster, including individuals from Estonia, and, as far as I know, Russians with substantial ancestry from near Estonia. As mentioned above, it's closely related to the Erzya + Moksha, Northwest Russian and Vologda clusters. However, it's clearly much more western than any of these clusters (for instance, see the PCA below), which suggests Germanic influence in its makeup.

pop10 - Cornish: Almost all of my Cornish samples from the 1000 Genomes Project feature in this very local cluster, which shows the highest chunk count among the Western European samples. The overall results suggest a lack of outbreeding in recent times.

pop11 - French/Belgian: Interestingly, this cluster includes the bulk of the French samples, a French Canadian, and two Belgians. On the other hand, the most northerly French are placed in the more cosmopolitan Northwest European cluster (see below).

pop12 - Lithuanian: All of the more or less pure Lithuanians fall in this cluster. Those that don't are a reference sample from Behar et al. 2009, who always appears very Belorussian like in other analyses, and here sits in the East Slavic cluster, and a project member with recent German ancestry (LIT3). The Western European influence carried by the latter pushes him into the Polish/West Ukrainian cluster, despite not having any documented Polish or Ukrainian ancestry.

pop13 - Northwest Russian: This cluster appears to be made up of Russians who have more Finnic, and/or perhaps Eastern Baltic, ancestry than the individuals in the East Slavic cluster. In other words, it's more northerly, less westerly, and more closely related to the Finnic-speaking Erzya, Moksha and Estonians.

pop14 - Irish + West British: Most Irish individuals fall in this cluster, as well as British samples from Western Scotland and Wales. It's tempting to correlate this cluster with Celtic genetic ancestry in the Isles.

pop15 - South/West Scandinavian: This is basically a Norwegian and Southern Swedish cluster. It also features Swedes from other parts of the country who most likely have some German, Walloon and/or French influence.

pop16 - East German: This cluster includes individuals with significant or even overwhelming Germanic ancestry, but also with very clear Western Slavic input. One of the individuals here is of mixed Polish, German and Swedish ancestry, which pretty much sums up the character of this cluster in a modern context. The presence of two Hungarians from Behar et al. 2009. isn't surprising, because Hungary was settled by both Germanic and Western Slavic groups from the early Middle Ages until modern times.

pop17 - Northwest European: I had reasonable hopes of breaking up this large cluster into a couple of units at least. However, that did not happen, and I don't think it will unless I obtain more samples from the relevant areas of Europe, like Holland and specific parts of the UK. I think the main reason this cluster failed to budge was because of its cosmopolitan nature. In other words, the samples here include some of the most outbred in the analysis, and this, coupled with the fact that they carry very similar ancestral components, means that fineSTRUCTURE doesn't have anything to latch onto to create divisions.

pop18 - East Scandinavian: This could also be called a Swedish cluster. It's almost entirely made up of Swedes, usually from Eastern or Southeastern Sweden, and/or occasionally with recent Finnish influence.

pop19 - Polish/West Ukrainian: The vast majority of the Poles fall in this cluster, and about half of the Ukrainians from Yunusbayev et al. 2011. Most of these Ukrainians appear to be from the Lviv district in the west, and some might even have fairly recent Polish and/or German ancestry. In fact, I would say the latter is a good bet for UkrLv240Y, who shows large Western European segments on several chromosomes.

pop20 - East Slavic: All of the Belorussians cluster here, and so do Russians from near Belorussia and Ukraine, and almost half of the Ukrainians from Yunusbayev et al. 2011 (those who show more easterly genetic characteristics). An individual of mixed Polish and Lithuanian ancestry also makes an appearance here, suggesting that one of the main factors differentiating this cluster from the Polish/West Ukrainian group is a higher level of Baltic admixture in the former.

pop21 - East Central European:
This cluster is based on most of the Hungarians in my dataset, but it also includes a number of Western and Southern Slavs, often with significant German ancestry. Not surprisingly, this cluster shows very high affinity with both the East German and Polish/West Ukrainian clusters.

Let's now move on to some graphics. Below, in order of appearance, are the following: raw data coancestry matrix, showing the placement of individual samples; aggregate coancestry matrix, showing the populations (or clusters) described above; pairwise coincidence matrix, which is useful for spotting very recent ancestral ties; a PCA plot of the 21 clusters. More detailed ChromoPainter/fineSTRUCTURE PCAs of Western Europe can be found at this link.





Finally, those of you who wish to run your own experiments with the ChromoPainter datasheets from this analysis can download them here. Please note, the sheets don't reveal any raw or traits/disease data.

Friday, January 13, 2012

Eurogenes' North Euro clusters - phase 1, exploring the data


I have some preliminary results from a new intra-North Euro cluster analysis, using a cutting edge tool called ChromoPainter. More than 400 samples and 270K SNPs were tested, in linkage mode, and then the output processed in fineSTRUCTURE at 200K burn-ins and iterations. Like I say, the results should be treated as preliminary, but they already look better than any other cluster analysis I've ever seen dealing with Europe north of the Alps, Pyrenees and Balkans. The algorithm identified 21 clusters, with most located in Eastern and Northeastern Europe (see spreadsheet for details). Below are two plots showing how the clusters relate to each other via a tree diagram and heat maps – the first shows an aggregate view, and the second the individual samples.





It's interesting that the Baltic Finns seem to create clusters at a drop of a hat, but they also share the highest number of chunks, and the longest chunks, than any other group. Indeed, all of the Finnish clusters are closely related, and many of the individuals, especially from East Finland, even look like distant relatives on the heat map (note the ultra-hot, blue squares). On the other hand, the large Northwestern European cluster, featuring samples from across the UK, as well as from several nearby countries, is holding firm, and might be tough to break up in this analysis.

I have some theories about the reasons for the obvious genetic homogeneity and diversity in Western Europe, and these include the effects of the Black Death. It decimated many populations in the western half of the continent, thus encouraging migrations into emptied areas, and eventually leading to more open, mobile societies. It's an interesting subject, and I might write much more on it in the future. Meantime, here's a PCA plot from the ChromoPainter chunk counts data. Note the large distances spanned by groups from Northern and Eastern Europe, and the tight bundle of samples from the west, mostly from the UK, Ireland, France and the Low Countries. Interestingly, and perhaps counter-intuitively, it's the closely related Finns who take up most of the space on the plot.



The first component picked up by this PCA appears to be an Atlantic one. It peaks in the Cornish samples, but shows similar levels in all the British, Irish, French, Dutch and Belgians (post-Black Death mobility?). If we are to assume that I identified the component correctly, then it appears as if the East Finns, Vologda Russians, Erzya from the Middle Volga, and Lithuanians are the least “Atlantic” samples in this analysis. These groups, especially the East Finns, also happen to act like relative genetic isolates in many of my experiments (such as ADMIXTURE and MDS analyses). Thus, it seems they've been sheltered from significant gene flow from outside in recent times, including from the west, like German emigration to East Central Europe and Scandinavian influence in Western and Southwestern Finland.