Saturday, July 13, 2013

A Hidden Gem in Central Asia: Previously Unknown Y-DNA R1b Haplotype [Original Work]

1. Introduction

Central Asian Y-DNA diversity has been an area of constant intrigue in the genetics community. Wells et
al.'s The Eurasian Heartland: A continental perspective on Y-chromosome diversity paved the way, with several others following in their regard. Members of the same team (including Dr. Wells) produced another paper - A Genetic Landscape Reshaped by Recent Events: Y-Chromosomal Insights into Central Asia - on the same topic in the following year, this time headed by Dr. Tatania Zerjal. I noted a greater emphasis on East-Central Asian populations as well as a mentioning of Y-STR analysis in the study itself. However, none of this data was supplied, with only Y-SNP information included (shown sporadically in this entry). The age of this paper is apparent through the nomenclature used (see Method section).

Several months ago, I made a request to obtain the Y-STR data from this study to one of the co-authors, Dr. Tyler-Smith, who kindly replied with the results of all sampled populations (Data Sink > Zerjal et al. Raw Data).

In this blog entry, the Y-STR data is showcased with a special emphasis on the Y-DNA R1b-M269 which was discovered.


2. Method
Y-SNP Phylogeny in original paper (Zerjal et al.) [1]


The maximum number of compatible Y-STR's were utilised for processing in Urasin's YPredictor for easier haplogroup identification (14 of a possible 16, DYS434 and 435 were excluded). All data was run through YPredictor. Only samples with ≥70% probability were included in the final results (Data Sink > Processed Data). As discussed below, relevant findings are compared with the basic Y-SNP haplogroups shown in the original study (on right).

One point which needs to be addressed immediately is the high frequency of "_DE-M1" and "P-M45". It appears that the STR selection has led to a phantom result, rendering many of the samples useless. For instance, the original study shows the Kazakhs belong overwhelmingly to C3c-M48, [1] although the probable results shown here are mostly "_DE-M1".  The exclusion of DYS434 and 435 from my level of processing likely contributed to this; if one assigns equal weight to the statistical strength of a prediction, removal of two STR's from a panel numbering 16, accuracy is reduced by 12.5%. Additionally, some conversion error seems to have applied with DYS437 (i.e. a value <12 is unusual). Therefore, "_DE-M1" and "P-M45" results were dismissed on account of the mismatch between predicted and likely confirmed haplogroups probably due to a compatibility issue between the study's STR panel and YPredictor..


3. Results

As the majority of samples were removed owing to the caveat described above, this entry will take a qualitative rather than quantitative approach to analysis on the general picture formed. Much of the remaining results are congruent with findings in other papers. Populations around the Caucasus are signified by plenty of R1b-M269, J2a-M410 and G2a-P15. Tajiks and the Kyrgyz were predominantly R1a1a-M17. Mongolians and other East-Central Asian ethnic groups yielded the most O3-M122 and "NO-M14" (likely to be Y-DNA N or O suffering from the STR restrictions described in the Method section).


Y-SNP distribution in Central Asia (Zerjal et al.) [1]


3i. The R1b Signal 

R1b-M269 was found across Central Asia and not only in the Caucasus (Armenians, Azeris, Georgians, Ossetians). It was mostly detected among the Turkmen (trk1, trk2, trk4, trk6, trk7, trk22, T29, T32) with a single sample among the Uzbek (uz-s110). [1]

Analysis of the haplotypes (including DYS434 and DYS435) revealed the nine Central Asian R1b samples belonged to a secure haplotype (Data Sink > R1b Results). trk6 diverged greatest, albeit with two 1-step mutations on DYS393 and DYS434. The rest match this haplotype exactly or have single 1-step mutations. [1] When this Central Asian R1b haplotype is compared with the other Caucasian samples, a mixed picture emerges, with the poorest being an Armenian (arm47) at 8/16, whereas the best are another Armenian (arm12) and Azeri (az48), both at 15/16. [1]
One interesting point is the Kurds sampled in this study (some of whom also belong to R1b-M269) are actually the displaced population positioned on the Iranian-Turkmenistani border. All of whom match the Central Asian R1b haplotype with a similar value (12-13/16). This definitively rules out the Kurds as a source for the haplotype, particularly as better matches can be found further to the west. It should be noted the Kurds themselves formed their own R1b haplotype (defined here by DYS389II=27, DYS391=10). [1]

In summary, the data reveals that the Turkmen are particularly abundant in R1b-M269 and all belong to the same haplotype as one of the Uzbek samples. This haplotype matched some Caucasians very well, but others not so well. The Kurds living in Turkmenistan belonged to their own haplotype.


3ii. Is This Actually R1b-M269?

Attention must first be shown to the original paper again; any potential R1b-M269 here will be present as P(xR1a)-92R7 (shown in the paper as "Haplogroup 1"). [1] Evidently, this makes up approximately half of the Turkmen lines and a quarter of Uzbek ones. Other haplogroups (such as other forms of R1b, R2a-M124, various Q subclades) presumably make up the rest of "Haplogroup 1" shown.

The next step is to verify whether or not this Central Asian R1b haplotype matches other R1b haplotypes online. As Y-DNA R1b-M269 is fortunately well-represented in the world of genetic genealogy, searching for the haplotype's matches on ySearch is a reasonable enterprise. DYS437 had to be excluded here due to a conversion issue, leaving the haplotype at 15 STR's. A genetic distance (GD) of 3 was allowed on these 15 markers. Results are shown on the right.

ySearch results for Central Asian R1b haplotype
With some confidence, the search has demonstrated that the Central Asian R1b haplotype does indeed belong to R1b-M269, as all the seven matches shown (one of whom is Armenian) belong to it.

Expanding the line of inquiry one further step came through comparing this haplotype with Iranian haplotypes [2] which were readily available. Due to differences in STR panels (an overlap of only 11) this proved to be inconclusive, aside from the observation that DYS389i+ii was completely different between the Central Asian modal (10-26) and the Iranian values. At this point I suspect that, much like DYS437, there is a conversion issue with DYS389 also.

Finally, a comparison was made with the R1b found in Afghanistan last year [3]. Interestingly, if DYS389i+ii and DYS437 are excluded, the two Uzbeks (samples 35 and 181) match the Central Asian R1b haplotype almost exactly based on the remaining 11 STR's. The one Tajik (sample 32) is less likely to be related due to two 1-step mutations on different STR's.


4. Conclusion

The inferences made from the data hang by a metaphorical thread due to the persistent STR issue; different labs have used different panels in the past decade, making it excruciatingly difficult to use materials from older papers. Fortunately, the presence of a specific strain of R1b-M269 in Central Asian (in Turkmen and Uzbeks) has successfully been demonstrated after select exclusions and no modifications to the data.

However, some larger questions remain. If STR limitations were not an issue, how would the Iranians from Haber et al. have compared? Would the Tajik from the other Haber et al. paper have belonged to the same haplotype in the end?

The origin of this Central Asian R1b haplotype will, I anticipate, also be a point discussed heavily among interested parties. At this point in time, I must stress that none of the evidence thus far points to anything in particular without ruling other theories out, although it leaves the door for interpretation wide open.

Having given this cautionary statement, the main thrust of this entry should be emphasised; R1b-M269 in Central Asia is a confirmed reality and here to stay. I will defer any subsequent analyses to the experts on Y-DNA R1b which grace several genetic genealogy boards for their take on the flavour of this haplotype.


5. Acknowledgement

I publicly extend my gratitude to Dr. Tyler-Smith for being so kind in sending me the raw STR's from this important paper for my research, as well as co-authoring the other two excellent studies I have cited here and in the past.


6. References

1. Zerjal T, Wells RS, Yuldasheva N, Ruzibakiev R, Tyler-Smith C. A genetic landscape reshaped by recent events: Y-chromosomal insights into central Asia. Am J Hum Genet. 2002 Sep;71(3):466-82. Epub 2002 Jul 17.

2. Haber M, Platt DE, Badro DA, Xue Y, El-Sibai M, Bonab MA. Influences of history, geography, and religion on genetic structure: the Maronites in Lebanon. Eur J Hum Genet. 2011 Mar;19(3):334-40. doi: 10.1038/ejhg.2010.177. Epub 2010 Dec 1.

3. Haber M, Platt DE, Ashrafian Bonab M, Youhanna SC, Soria-Hernanz DF, Martínez-Cruz B. Afghanistan's ethnic groups share a Y-chromosomal heritage structured by historical events. PLoS One. 2012;7(3):e34288. doi: 10.1371/journal.pone.0034288. Epub 2012 Mar 28.

Tuesday, March 26, 2013

Y-DNA Haplogroup N in India: Wayward Uralics or Lab Error? [Original Work]

Introduction
Y-DNA Haplogroup N Eurasian Distribution


Per ISOGG's 2013 SNP tree and as has been the case for years, Y-DNA Haplogroup N is defined by the M231 mutation (G->A at rs9341278) on the Y-Chromosome. With a predominantly North Eurasian distribution, it peaks in Europe among the Finnish people and various ethnic groups residing in Russia's far north through the N1c-Tat subclade. N1c-Tat specifically is frequently associated with Uralic-speaking populations in the literature.

Haplogroup N also appears to have an association with Central Asia as shown in the N Y-DNA Haplogroup Project (FTDNA) results, with several samples coming in from Kazakhstan, Uzbekistan and Mongolia. It has also been observed in Turkey (KurdishDNA blog entry) as well as appearing in 1.6% of Iran's Azeri population (Grugni et al. entry).

The finding of Haplogroup N in India through Sharma et al.'s The Indian origin of paternal haplogroup R1a1* substantiates the autochthonous origin of Brahmins and the caste system [1] is a curious one. Unfortunately, the paper did not include any Y-STR material to help understand the basis of N's presence in India.


Significance of Potential Haplogroup N in India

Linguistics provides us with a plausible scenario regarding how Haplogroup N may have arrived in the Indian Subcontinent. Contacts between early Finno-Ugric and Indo-Iranian groups took place around the Ural mountains, specifically between the forest and steppe zones. Evidence of transmission in horsekeeping techniques, economy, deities and common words are firmly established from Andronovo archaeological horizon on the steppes into the "Andronovoid" societies living in the nearby forests. [2]

The presence of Haplogroup N in India, if present in relevant populations and displaying MRCA values or STR clusters consistent with a Neolithic origin further north, would satisfy the likelihood of Haplogroup N representing an accompanying genetic signal from the steppe zone roughly four thousand years ago, as well as serving as a genetic remnant of the interactions that undoubtedly took place between Indo-Iranian and Finno-Ugric tribes.


Current Findings

In 2009, Sharma et al. published a paper highlighting the Y-Chromosome haplogroup differences between various upper caste (Brahmin) and tribal populations across India. The paper went on to deduce that Haplogroup R1a1a in India was autochthonous in origin based on their findings [1] (now disputable and improbable based on Underhill et al.'s landmark study on Y-DNA R1a1a and recent findings by the R1a Subclades FTDNA Project, although this topic is beyond the scope of this entry).

It was this very paper by Sharma et al. which revealed the presence of Y-DNA N in India. Haplogroup N1-LLY22g was found in Brahmins from Gujarat, Madhya Pradesh and Mahastra (3.13%, 2.38% and 3.33% respectively), as well as tribal populations from Uttar Pradesh (1.56%). Their results were extended to include greater caste differentiation (Brahmins vs. Scheduled Castes vs. Tribals); here, Brahmins were found to have five times greater the frequency of N1-LLY22g than tribal groups (0.5% vs. 0.1% respectively).  [1]

Although the frequencies were arguably insignificant, the inference stood - Y-DNA Haplogroup N showed an association with the upper caste practitioners of Hinduism in India, paving the way for the scenario described in the above chapter to be considered.

However, the strength of this conclusion is weakened greatly by cross-sectional data from numerous studies concerning the Indian Subcontinent produced in the past decade:


  • Sengupta et al. (2006) revealed that, out of 1090 samples, with the majority coming from the Indian Subcontinent, the only populations revealing any Haplogroup N (N-M231) and associated downstream subclades were either East Asian (Chinese ethnicities, Cambodian) or Siberian (Yakut). [3] No groups from India belonged to Haplogroup N-M231.
  • Furthermore, Sahoo et al. (2006) also sampled individuals from across the Indian Subcontinent (n=1074) and failed to find a single instance of N-M231. [4]
  • In a recent study on various populations in Tamil Nadu (South India), Haplogroup N was completely absent in the 1680 samples tested. [5]
  • Y-DNA N1c-Tat was absent in the 607 tribal populations tested from East and Northeast India. [6]
  • Returning to the north of the country, 560 men from various upper castes and Muslim groups were tested by Zhao et al. and N1c-Tat was absent from all. [7] 
  • Focused specifically on Brahmins from Saraswat (Jammu-Kashmir), Yadav et al. found none of the approximately 109 haplotypes to belong to any derivative of Haplogroup N. [8]

Finally, the N Y-DNA Haplogroup Project at FTDNA currently does not show any samples whatsoever from the Indian Subcontinent.


Possible Explanation

Despite over 4,000 samples over five studies representing various groups from across India, not a single trace of Haplogroup N has been detected. What explains this glaring discrepancy with Sharma et al.'s findings? Differences in sampling strategy between the other studies with Sharma et al. cannot account for this; there is enough regional overlap to rule this out.

As was the case with Sengupta et al. where several Hazara haplogroup classifications were allegedly due to a laboratory error, it is probable the Haplogroup N seen here follows the same suit. By reasonable deduction, if one study reveals a trend that several others covering thousands of samples cannot verify, there must be something intrinsically erroneous in the former.


Conclusion

Until I can physically view the purported Haplogroup N haplotypes reported in Sharma et al., it is the conclusion of this entry that they are most likely the result of a laboratory error given the complete absence of any flavour of N-M231 in India through other recent studies. If any Haplogroup N is found, it must be contrasted against Sharma et al. and should be investigated on a separate line of inquiry. As ever, details of any future cases of Haplogroup N in India should be taken into consideration. If of a Mughal background, the paternal origins are readily explained by Medieval Central Asian ancestry. If from the furthest northeast of the Indian Subcontinent, the possibility of Nepali ancestry should be sought. [9] Although prehistoric indirect influence from Finno-Ugric interactions in the second millennium BC onwards shouldn't be dismissed outright, other more recent explanations exist.


References


1. Sharma S, Rai E, Sharma P, Jena M, Singh S, Darvishi K. The Indian origin of paternal haplogroup R1a1* substantiates the autochthonous origin of Brahmins and the caste system. J Hum Genet. 2009 Jan;54(1):47-55. doi: 10.1038/jhg.2008.2. Epub 2009 Jan 9.

2. Kuz'mina EE. The Origin of the Indo-Iranians. Koninklijke Brill NV, Leiden, The Netherlands. 2007.

3. Sengupta S, Zhivotovsky LA, King R, Mehdi SQ, Edmonds CA, Chow CE. Polarity and temporality of high-resolution y-chromosome distributions in India identify both indigenous and exogenous expansions and reveal minor genetic influence of Central Asian pastoralists. Am J Hum Genet. 2006 Feb;78(2):202-21. Epub 2005 Dec 16.

4. Sahoo S, Singh A, Himabindu G, Banerjee J, Sitalaximi T, Gaikwad S. A prehistory of Indian Y chromosomes: evaluating demic diffusion scenarios. Proc Natl Acad Sci U S A. 2006 Jan 24;103(4):843-8. Epub 2006 Jan 13.

5. Arunkumar G, Soria-Hernanz DF, Kavitha VJ, Arun VS, Syama A, Ashokan KS. Population differentiation of southern Indian male lineages correlates with agricultural expansions predating the caste system. PLoS One. 2012;7(11):e50269. doi: 10.1371/journal.pone.0050269. Epub 2012 Nov 28.

6. Borkar M, Ahmad F, Khan F, Agrawal S. Paleolithic spread of Y-chromosomal lineage of tribes in eastern and northeastern India. Ann Hum Biol. 2011 Nov;38(6):736-46. doi: 10.3109/03014460.2011.617389. Epub 2011 Oct 6.

7. Zhao Z, Khan F, Borkar M, Herrera R, Agrawal S. Presence of three different paternal lineages among North Indians: a study of 560 Y chromosomes. Ann Hum Biol. 2009 Jan-Feb;36(1):46-59. doi: 10.1080/03014460802558522.

8. Yadav B, Raina A, Dogra TD. Genetic polymorphisms for 17 Y-chromosomal STR haplotypes in Jammu and Kashmir Saraswat Brahmin population. Leg Med (Tokyo). 2010 Sep;12(5):249-55. doi: 10.1016/j.legalmed.2010.05.003.

9. Gayden T, Chennakrishnaiah S, La Salvia J, Jimenez S, Regueiro M, Maloney T. Y-STR diversity in the Himalayas. Int J Legal Med. 2011 May;125(3):367-75. doi: 10.1007/s00414-010-0485-x. Epub 2010 Jul 21.