Total Pageviews

Showing posts with label Z93. Show all posts
Showing posts with label Z93. Show all posts

Saturday, August 18, 2012

R1a1a comparison STR111 Part II

First of all, I want to thank Humata, the blogger from http://vaedhya.blogspot.com/. He helped me analyzing R1a1a; I used the same tools that he used to analyze haplogroup Q in his recent post.

So, with Humata's help I analyzed all R1a1a individuals that have STR111 data. Again, I used adjusted distances to perform this analysis, which I described and used  here and here.

In order to analyze the data I had to increase IDs to 10 digits to prevent malfunction of the used software Fitch.

There are various versions to illustrate the data. Here, I am presenting two layouts.

Rectangular tree layout (see high resolution image):


Polar tree layout (see high resolution image):




Again, to better see "what is what" I annotate each ID with the proposed group used in the R1a1a and Subclades Project at FTDNA, and I used the same colors for the subclades as in this figure from FTDNA. Even with 111 STR values the main R1a1a SNPs (Z93, Z283, etc.) are overlapping.

So what does give more accurate results, unrooted network analysis or rooted tree analysis?

It is quiet obvious that even STR111 data are not sufficient enough to differentiate between the major R1a1a subclades. There are multiple overlapping haplogroups in the tree (presented below) and in the network (presented previously), i.e. Z283 and Z93 are overlapping when focusing only on the STR111 data. All previously presented trees are adding SNP information to the tree to correct this obvious overlapping. As an example: 
179005     Krikor Mirijanian, Arapkir, Turkey who is Z93+, L342+, L657-. 

Based on his STR111 values 179005 is closest to Z283+, Z280+ and Z283+, Z284+ individuals and not to other Z93+ individuals.


Update:
After Semargl was asked how he generated his R1a1a STR111 tree and why his tree shows clear clusters along the SNP branches, he responded that he is using not only STR but also SNP information to generate the tree.  Additionally, his current phylogenetic tree is a cladogram, that means that the cladogram tree does not have any information about the age or diversity of the R1a1a branches, e.g. the Ashkenazi-Jewish Z93+, L342+, L657- cluster and the MacDonald cluster takes a large part of the tree, even though these clusters are known to be very narrow (low diversity). Hopefully, he can generate a new tree that includes all this information.



From the Fitch software manual:

In Fitch you can also randomize the input order of the sequences with option "j", jumble. Often the input order of the sequences affects the outcome of the analysis. This can be assessed by randomizing the input order. The program also asks you to specify the number of times you want to randomize the input order of the sequences. It is advisable to do jumbling at least 10 times, because it almost certainly improves the results.

This is why I repeated the analysis with 10 runs as advised. Indeed, the R1a1a STR111 tree looks a little bit better now.


Rectangular tree layout (see high resolution image):


Polar tree layout (see high resolution image):





One of the new discovered SNPs is Z1282, downstream of L342. It was found in N77532 Sundardas Tulsyan, India. He is in the tree as # "N77532-2C*". Based on the presented tree and the previously presented network analysis #184336 , SAUD ABDUL AZIZ, Qatar would be a good candidate for L1282, too (1844336-2C* in the tree and in the network analysis).
 

Monday, August 13, 2012

R1a1a comparison STR111

Today, I want to present results I got by combining two methods that I used before. The first method is described here and here; the second method is based on SplitsTree that I also used here and here.

I used all R1a1a individuals with STR111 data from the R1a1a and Subclades Project at FTDNA, a total of 203 individuals (and lots of people with unknown SNP status). The two methods don't use any SNP information, so the clusters are enterily based on STR values and the mutation rate of each STR.

Here is the network as a high resolution jpg.



Update:
To better see "what is what" I annotate each ID with the proposed group used in the R1a1a and Subclades Project at FTDNA, and I used the same colors for the subclades as in this figure from FTDNA.


Here is the color-coded network as a high resolution jpg image.


 

Sunday, August 12, 2012

L342+ comparison STR67

Today, I want to present results I got by combining two methods that I used before. The first method is described here and here; the second method is based on SplitsTree that I also used here and here.

I used all L342+ individuals with at least STR67 data from the R1a1a and Subclades Project at FTDNA, so L657- and L657+ individuals (and lots of people with unknown SNP status).

Here is the network as a jpg. (Update: New results of N101746 are included in jpg link.)
I started to annotate the different clusters and I will finish it soon, still I want to share the current status.

Update:
N101746 (Central&Southwest Asian; India) is in a cluster with:
M6851 (Arabic II; Saudi-Arabia),
M6699 (Arabic II; Unknown),
M6698 (Arabic II; Unknown),
M6853 (Arabic II; Unknown),
197670 (Central&Southwest Asian; India),
M6132 (Central&Southwest Asian; UAE), and
160543 (Central&Southwest Asian; Iraq).
The "Arabic II" individuals are very close to each other, while the "Central&Southwest Asian" individuals in this cluster show more diversity.

This described cluster above is close to another cluster that is very narrow. The follwoing individuals belong to this cluster:
157103 (Arabic; Saudi-Arabia),
160271 (Arabic; Qatar),
162855 (Arabic; Saudi-Arabia),
178907 (Arabic Saudi-Arabia),
178905 (Arabic; Saudi-Arabia),
157621 (Arabic; Saudi-Arabia),
157621 (Arabic; Saudi-Arabia),
157619 (Arabic; Saudi-Arabia),
M6895 (Arabic; Saudi-Arabia),
178906 (Arabic; Saudi-Arabia),
M7066 (Arabic; Unknown),
M6183 (Arabic; Unknown),
M6679 (Arabic; Unknown),
M6458 (Arabic; Kuwait),
M7013 (Arabic; Kuwait),
M6982 (Arabic; Unknown), and
M6285 (Arabic; Qatar).

Zoom in:


Wednesday, May 16, 2012

Deep look into R1a1a

I tried to analyze the currently available STR data for Kurds having the Y-haplogroup R1a1a.
I will focus on these STRs:
DYS393    DYS390    DYS19    DYS391    DYS385a    DYS385b    DYS426    DYS388    DYS439    DYS389i    DYS392    DYS389ii    DYS458    DYS459a    DYS459b    DYS455    DYS454    DYS447    DYS437    DYS448    DYS449    DYS464a    DYS464b    DYS464c    DYS464d    DYS460    Y-GATA-H4    YCAIIa    YCAIIb    DYS456

The ranking is sorted by similarity: (number of equal STRs/number of comparable STRs**) in percentage (%)

**comparable STRs= STRs that got tested by both sides that were compared

Let's start with H1483, he is being tested as Z93+ L342+ (L657-)





In this case the STR results fit with the known SNP results. H1483 is sharing the highest similarity with other people that are Z93+ L342+.


The next individual that I want to compare is 214352, he is being tested as Z283+ M417+ (Z93- Z280- Pk5- P98- M64.2- M56- M458- M434- M334- M157.1- L260- L176.1- L175- Z284-)






Here, the STR results do not fit with the known SNP results. 214352 is sharing a relatively high similarity not only with other people that are Z283+ but also with some that are Z283- and/or are not even R1a. He belongs to the root of Z283+. We already know that his genetic results are special and intriguing.


Besides the two previously described individuals we now have a 3rd individual with R1a1a. He is a Zaza Alevi from Dersim. He is AXYH2. He did not get any SNPs tested, so I have to focus on his STR results.

Let's compare him to others:






The picture is not clear but if I would have to predict his SNPs, I would guess Z93+.



Same exercise as before but this time with 37 STRs:
DYS393    DYS390    DYS19    DYS391    DYS385a    DYS385b    DYS426    DYS388    DYS439    DYS389i    DYS392    DYS389ii    DYS458    DYS459a    DYS459b    DYS455    DYS454    DYS447    DYS437    DYS448    DYS449    DYS464a    DYS464b    DYS464c    DYS464d    DYS460    Y-GATA-H4    YCAIIa    YCAIIb    DYS456    DYS 607    DYS 576    DYS 570    CDY a    CDY b    DYS 442    DYS 438

Let's start with H1483 again, he is being tested as Z93+ L342+ (L657-). Here, I could not compare 37 STRs but only 32 STRs.




Next is 214352, he is being tested as Z283+ M417+ (Z93- Z280- Pk5- P98- M64.2- M56- M458- M434- M334- M157.1- L260- L176.1- L175- Z284-).








Latest R1a1a individual is AXYH2. He did not get any SNPs tested, so I have to focus on his STR results.