New add-on tool to MDLP World-22 calculator

Alex_Axe from Molgen forum announced the release of a new tool Oracle_AdMix4.

Basically, you can think of this tool as an useful add-on that extends the Oracle two-population results into four-population format.

Currently the tool support only DOS/Windows platform, and its usage is quite straightforward. The current distribution of  Oracle_AdMix4 includes the required file (data.txt) with the frequencies of 22 components in each of reference populations. The second file - input.txt - is the one you need to modify in order to use the tool on your own World-22 calculator results:

30     {Quantity of output results}
1,5    {Threshold of components to ignore noise} 
0,5    {Threshold of method to ignore noise}
2      {Power parameter for weighted method or zero for least-squares method}   
0      Pygmy
15     West-Asian
0      North-European-Mesolithic
0      Indo-Tibetan
0      Mesoamerican
0      Arctic-Amerind
0      South-America_Amerind
5      Indian
0      North-Siberean
25     Atlantic_Mediterranean_Neolithic
5      Samoedic
0      Indo-Iranian
0      East-Siberean
40     North-East-European
0      South-African
0      North-Amerind
0      Sub-Saharian
0      East-South-Asian
10     Near_East
0      Melanesian
0      Paleo-Siberian
0      Austronesian

All you need to do is to change the default values of components with your own results of World-22 calculator.

The output should appear in the following format:

1 Mordovian+Russian_South+Swedish+Swedish @ 20,749639
2 Estonian+Mordovian+Swedish+Swedish @ 23,267935
3 Finnish-North+Mordovian+Polish+Swedish @ 23,656645
4 Finnish-North+Mordovian+Russian_South+Swedish @ 23,83736
5 Estonian+Finnish-South+Mordovian+Swedish @ 24,192354
6 Mordovian+Polish+Swedish+Swedish @ 24,846391
7 Estonian+German-North+Mordovian+Swedish @ 25,176036
8 Lithuanian+Mordovian+Swedish+Swedish @ 25,511204
9 Finnish-North+Lithuanian+Mordovian+Swedish @ 25,523722
10 Finnish-South+Lithuanian_V+Mordovian+Swedish @ 25,539837

A quick update to SupportMix's Chromosome Painting

A thoughtful reader of our blog has noticed that some of chromosome_paintings (Chr 5 set 9, Chr 7 sets 4,5; Chr 9 set 7, Chr 11 set 4) were missing in the original tar.gz bundle  distributed via Google Data Drive. I've had to re-upload the archive (with missing files), the new location of  the archive is here.

Additional experiment.

In addition to that quick fix, i decided to test the accuracy of SupportMix's chromosome paintings by juxtaposing them over the MDLP-World22 chromosome graphs. Due to time limitations, i used only first 7 chromosomes of my own SNP data. At first, i ran the MDLP-World22 modification of DIYDodecad v2.1 in byseg mode on "windows" of 500 contiguous SNPs along a chromosome, slided  by increments of 50. After that i cut out chromosome paintings of each chromosome from SupportMix's graphic output and aligned them to the scale of corresponding DIYDodecad chromosome graphs:

 After the preliminary evaluation of results, i have mentioned an approximate correlation between the byseg-output MDLP World22-DIYcalculator and SupportMix for two major "components" in my genome (North-East-European and Atlantic-Mediterranean). Moreover, "Near-Eastern segments" (assigned by SupportMix"  partially overlaps with the peaks of "Near-East segments" in DIYDodecad output. However, the situation with the minor components is much less uncertain. The lack of correlation for the minor components could be explained by different factors:

1) DIYDodecad operates on the unphased raw data of genotypes

2) DIYDodecad program doesn't take into consideration genetic distance/recombination

3) last, but not least: small segments may appear more noisier than they are, because there may not be any informative SNPs in a particular region to distinguish between some of the minor ancestral components (Dienekes Pontikos' observation)

UPDATE: At first i thought that it would be a great idea to calculate the index of correlation between 'byseg' output and SupporMix' Tprobs output file. But it seems that the results are not directly comparable - the assignment of segments in 'byseg' is measured in frequencies, while the assignment of segments in SupportMix is expressed by probability of assignment. If someone has a solution to this problem, please let me know.

The component maps of MDLP-World22 calculator

I would like to express my gratitude to the fellow members of ABF - Loxias and Wojewoda- for creating amazing maps of components and plots showing interrelationships between 22 components.

Since the readers of my blog might be interested in visual inspecting of components , i have decided to upload them on-line. 

The first batch of maps  includes "the component portraits", created by Loxias:

West-Asian

Atlantic-Mediterranean

East-Siberian

Amerind

Indian

Indo-Iranian

Indo-Tibetan

Paleo-Siberian

Pygmy

Sub-Saharan

Near-East

North-East-European

North-Siberian

North-Mesolithic-European

 The second batch includes the PCA scatter plot of components, created by Wojewoda. PCA scatter plots one PCA component vs. another PCA component of data obtained from ancestry coefficients  collected for each of MDLP World22 components: the spots are connected over time to produce a trajectory for each component.