Y-chromosomal evidence for a limited Greek contribution to the Pathan population of Pakistan

24grammata.com/ αγγλικά/ ιατρική/ φυλή Καλάς

Sadaf Firasat1, Shagufta Khaliq1, Aisha Mohyuddin1, Myrto Papaioannou2, Chris Tyler-Smith3, Peter A Underhill4 and Qasim Ayub1

1Biomedical and Genetic Engineering Division, Dr. AQ Khan Research Laboratories, Islamabad, Pakistan
2Unit of Prenatal Diagnosis, Center for Thalassemia, Laiko General Hospital, Athens, Greece
3The Wellcome Trust Sanger Institute, Hinxton, Cambridge, UK
4Department of Genetics, Stanford University, Stanford, CA, USA

Correspondence: Dr Q Ayub, Biomedical and Genetic Engineering Division, Dr AQ Khan Research Laboratories, GPO Box 2891, Islamabad 44000, Pakistan. Tel: +92 51 926 1142; Fax: +92 51 926 1144; E-mail:[email protected]

Received 16 November 2005; Revised 9 August 2006; Accepted 1 September 2006; Published online 18 October 2006.

Abstract

Three Pakistani populations residing in northern Pakistan, the Burusho, Kalash and Pathan claim descent from Greek soldiers associated with Alexander’s invasion of southwest Asia. Earlier studies have excluded a substantial Greek genetic input into these populations, but left open the question of a smaller contribution. We have now typed 90 binary polymorphisms and 16 multiallelic, short-tandem-repeat (STR) loci mapping to the male-specific portion of the human Y chromosome in 952 males, including 77 Greeks in order to re-investigate this question. In pairwise comparisons between the Greeks and the three Pakistani populations using genetic distance measures sensitive to recent events, the lowest distances were observed between the Greeks and the Pathans. Clade E3b1 lineages, which were frequent in the Greeks but not in Pakistan, were nevertheless observed in two Pathan individuals, one of whom shared a 16 Y-STR haplotype with the Greeks. The worldwide distribution of a shortened (9 Y-STR) version of this haplotype, determined from database information, was concentrated in Macedonia and Greece, suggesting an origin there. Although based on only a few unrelated descendants, this provides strong evidence for a European origin for a small proportion of the Pathan Y chromosomes.
Keywords:

population genetics, Pakistan, Greek, Y chromosome polymorphism

Introduction

Pakistan lies within a region that was invaded by Alexander the Great in 327–323 BC,1 although archeological evidence in northern Pakistan suggests that the Greek influence predates this invasion.2 Many ethnically and linguistically distinct populations inhabit this region, three of which (Burusho, Kalash and the Pathan) claim to be descendents of Greek soldiers who invaded the Indian subcontinent.3, 4, 5 A preliminary study using a limited number of Y-chromosomal markers found no evidence for admixture between the Greeks and Burusho or Pathan, and provided ambiguous evidence of genetic admixture between the Greek and Kalash populations.6 A subsequent analysis of autosomal loci gave no indication that the Kalash or other populations were genetically related to the Greeks,7 and in another study the Kalash population were shown as a genetic isolate8 who might therefore have developed unusual genetic characteristics by drift.9

The genetic material comprising the great bulk of the patrilinearally inherited human Y-chromosome is effectively haploid and does not undergo inter-chromosomal recombination, making it useful for evolutionary studies and forensic investigations.10 In particular, it can provide very high resolution haplotypes with known phylogenetic relationships, and these can reveal low levels of admixture that would not be detected by other methods. In the present study, we have further investigated the origin and the genetic relationship of these three Pakistani populations with the extant Greek population by typing a large set of markers from the male-specific region of the Y chromosome in 77 Greeks and 875 Pakistani individuals, and applying analytical methods that should be more sensitive to low levels of admixture.

Materials and methods
DNA samples

DNA samples of 952 unrelated males were analyzed in this study. The DNA was extracted directly from peripheral blood mononuclear cells in case of the Greek samples (n=77) and from EBV-transformed lymphoblastoid cell lines for the Pakistani samples (n=875). The Pakistani samples included Burusho, Kalash and Pathan individuals. Informed consent was obtained from all participants in this study.
Genotyping

All samples were typed for 90 binary markers11 following a phylogenetic hierarchical approach to define evolutionary Y lineages or haplogroups.10 The markers and the method used for their genotyping are listed in Table 1. In addition, 16 multiallelic Y-STR (short-tandem-repeat) markers were analyzed in these populations by three multiplex PCR reactions12, 13 to identify Y haplotypes.

Data analysis

Population pairwise PhiST values were estimated by using the Arlequin package14 based on STR variation within haplogroups weighted as described previously.6 Population pairwise rho genetic distances were calculated according to Helgason et al.15 Median-joining networks were constructed by Network 4.1.1.216 using the following five-fold range weighting scheme. The weights assigned were specific for the haplogroup and took into account the Y-STR variation across the haplogroup in the Pakistani and Greek populations: variance 0.00–0.14, weight 5; variance 0.15–0.29, weight 4; variance 0.30–0.44, weight 3; variance 0.45–0.59, weight 2; variance >0.59, weight 1. Network was also used to estimate the time to the most recent common ancestor (TMRCA). The rho genetic distance matrix was used to construct phylogenetic tree by neighbor-joining method17, 18 using the MEGA software package.19

Results

The Y-chromosomal lineages and their frequencies in the Greeks, Burusho, Kalash, Pathan and the rest of the Pakistani population are shown in Figure 1. The combination of biallelic markers identified 12 Y-chromosomal haplogroups or lineages in the Greeks, 17 in the Burusho and 15 in the Pathan populations. Only eight Y lineages were found in the Kalash population. Principal component analysis of Y haplogroup frequencies incorporating published data from European20 and West Asian21, 22 populations (Figure 2) revealed that the Pakistani populations cluster together, separately from the Europeans, consistent with the previous conclusion that none of the Pakistani populations had a large male contribution from Greece, and demonstrating that this conclusion was not an artifact of the low phylogenetic resolution used before.

A rooted maximum-parsimony tree of Y lineages found in the Greek, Burusho, Kalash, Pathan and Pakistan. The lineages were defined by binary markers whose designations and population frequencies are given below each branch. The YCC lineage names10 are shown below the frequencies. Branch lengths are arbitrary.
Full figure and legend (37K)

A plot of the first two principal coordinates based upon the analysis of Y haplogroup frequencies in (a) Pakistani and Greek populations using 33 biallelic markers. (b) Pakistani, Greek and selected populations from Europe, Africa, Middle East and West Asia.
Full figure and legend (24K)

The genetic distances between the populations were then calculated using measures that are more sensitive to recent events (Table 2). Pakistani–Greek population pairwise PhiST values based on the variation of STRs within haplogroups6 ranged from 0.131 to 0.213, with the lowest value between the Pathan and the Greeks. Pairwise rho genetic distances (the number of steps between a haplotype in one population and the closest haplotype in the second population, averaged over all comparisons)15 ranged from 4.3 to 8.1, with the lowest value again between the Pathan and the Greeks. Phylogenetic analysis using the matrix of rho genetic distances between populations with tree validation carried out by bootstrap resampling (10 000 replicates) also demonstrated that of the three Pakistani populations, the Pathans were closest to the Greek (Figure 3). These results therefore suggest that there might have been a low degree of recent Pathan–Greek admixture. To investigate this possibility further, we have examined individual lineages.
Figure 3.
Figure 3 – Unfortunately we are unable to provide accessible alternative text for this. If you require assistance to access this image, please contact [email protected] or the author

Neighbor-joining tree showing the relationship between the Greek and three Pakistani ethnic groups. The tree is based on rho genetic distances. Bootstrap values from 10 000 replicates are shown.
Full figure and legend (4K)

Clade E lineages were more frequent in the Greeks (21%) as compared to Pakistan (4%). The majority of haplogroup E chromosomes belonged to clade E3b and all Greek and Pakistani samples were resolved into the branches E3b1 (M78) and E3b3 (M123). Among the three Pakistani populations claiming Greek descent, this clade was observed only in the Pathans. The Pathan samples belonged to clade E3b1 that constituted 17% of the Greek samples.

A median-joining network of clade E Y chromosomes was constructed in order to examine the genetic relationship between these Greek and Pathan samples. A duplication of 10 and 13 repeat units was observed in the clade-E-derived Y chromosomes for the trinucleotide repeat DYS425 and this locus was, therefore, excluded from the network. The most striking feature of this network was the sharing of haplotypes between the Pathan and Greek samples (Figure 4). One Pathan individual shared the same Y-STR haplotype with three Greek individuals, and the other Pathan sample was separated from this cluster by a single mutation at the DYS436 locus. This demonstrates a very close relationship between the Pathan and Greek E lineages, but how surprising is this?

Median-joining network of clade E lineages in Pakistan (open circles) and Greece (hatched circles). Circles represent haplotypes and have an area proportional to frequency. The Pathan individuals are shown in black.
Full figure and legend (12K)

Worldwide data for the 16-element haplotype are not available, but a subset of nine of the STRs are included in by the Y-STR Haplotype Reference Database (YHRD)23 and were used to search this. The haplotype DYS19=13; 389I=13; 389II=30; 390=24; 391=10; 392=11; 393=13; 438=10; 439=12 was found in 53 individuals in a worldwide population sample of 7897 haplotypes and was highly specific for the Balkans . The contour map shows a major concentration around Macedonia and Greece, with a low scattering in other European countries, Tunisia, West Africa and the Pathans. This gives a strong indication of a European, possibly Greek, origin of these Pathan Y chromosomes.

Contour map showing the 9 Y-STR haplotype frequency distribution in Eurasia and northern Africa. This haplotype was shared between three Greeks and a Pathan individual belonging to clade E3b1.
Full figure and legend (53K)
Discussion
An extensive analysis of Y diversity within Greeks and three Pakistani populations – the Burusho, Kalash and Pathan – who claim descent from Greek soldiers allowed us to compare Y lineages within these populations and re-evaluate their suggested Greek origins. This study as a whole seems to exclude a large Greek contribution to any Pakistani population, confirming previous observations.7 However, it provides strong evidence in support of the Greek origins for a small proportion of Pathans, as demonstrated by the clade E network  and the low pairwise genetic distances between these two populations.

The Pathans were the only population among the three that claim Greek ancestry in which clade E was present. This branch is observed in Europe, Middle East, North and East Africa with a suggested origin in East Africa.24 Sub-clade E3b is common in Europe and probably originated in Africa.25 Compelling evidence in support of the genetic relationship between the Pathan and Greek E3b1 Y chromosomes was provided by the median-joining network (Figure 4). One Pathan shared a Y-STR haplotype, that included a duplication of 10 and 13 repeat units for the DYS425 locus, with three Greek individuals and the other was separated from this cluster by a single mutation, which enabled us to estimate the TMRCA (meanplusminusSD) using the Network software as between 2000plusminus400 and 5000plusminus1200 YBP depending upon the observed26 or inferred mutation rates,27 respectively. This coincides with the period of Alexander’s invasion during 327–323 BC. This haplotype was not observed in any other E3b1-derived Pakistani Y chromosome but was highly specific for the Balkans – the highest frequency being in Macedonia.

It is worth emphasizing here that the chance of picking up rare events largely amplified by drift affecting a limited portion of the population cannot be discounted, and Cruciani et al28 also recommend caution when using microsatellite alleles as surrogates of unique event polymorphisms. The genetic data alone do not tell us when the Balkan chromosomes arrived in Pakistan: it is necessary to turn to the historical record for this. There has been no known Greek admixture within the last few generations, but in addition to Alexander’s armies, the possibility of admixture between the Greek slaves who were brought to this region by Xerxes around one hundred and fifty years before Alexander’s arrival, and the local population, cannot be discounted. At that time Afghanistan and present day Pakistan were part of the Persian Empire.1 Nevertheless, Alexander’s army of 25 000–30 000 mercenary foot soldiers from Persia and West Asia and 5000–7000 Macedonian cavalry29 perhaps provides a more likely explanation because of their elite status and substantial political impact on the region.

References

Wolpert S: A new history of India. New York: Oxford University Press, 2000.
Jettmar K: Beyond the Gorges of the Indus: Archaeology before excavation. Karachi, Pakistan: Oxford University Press, 2001.
Biddulph J: Tribes of the Hindoo Koosh. Karachi, Pakistan: Indus Publications, 1977.
Lines M: The Kalasha people of North-western Pakistan. Peshawar, Pakistan: Emjay Books International, 1999.
Caroe O: The Pathans. Karachi, Pakistan: Oxford University Press, 1958.
Qamar R, Ayub Q, Mohyuddin A et al: Y-chromosomal DNA variation in Pakistan. Am J Hum Genet 2002; 70: 1007–1024. | Article |
Mansoor A, Mazhar K, Khaliq S et al: Investigation of the Greek ancestry of populations from northern Pakistan. Hum Genet 2004; 114: 484–490. | Article | PubMed | ChemPort |
Rosenberg NA, Pritchard JK, Weber JL et al: Genetic structure of human populations. Science 2002; 298: 2381–2385. | Article | PubMed | ISI | ChemPort |
Rosenberg NA, Mahajan S, Ramachandran S, Zhao C, Pritchard JK, Feldman MW: Clines, clusters, and the effect of study design on the inference of human population structure. PLoS Genet 2005; 1: e70. | Article | PubMed | ChemPort |
Jobling MA, Tyler-Smith C: The human Y chromosome: an evolutionary marker comes of age. Nat Rev Genet 2003; 4: 598–612. | Article | PubMed | ISI | ChemPort |
Underhill PA, Shen P, Lin AA et al: Y chromosome binary haplotypes and origins of modern human populations. Ann Hum Genet 2001; 26: 358–361.
Thomas MG, Bradman N, Flinn HM: High throughput analysis of 10 microsatellite and 11 diallelic polymorphisms on the human Y-chromosome. Hum Genet 1999; 105: 577–581. | Article | PubMed | ChemPort |
Ayub Q, Mohyuddin A, Qamar R et al: Identification and characterisation of novel human Y-chromosomal microsatellites from sequence database information. Nucleic Acids Res 2000; 28: e8. | Article | PubMed | ChemPort |
Schneider S, Kueffer J-M, Roessli D, Excoffier L: Arlequin ver 1.1: a software for population genetic data analysis. Switzerland: Genetics and Biometry Laboratory, University of Geneva, 1997.
Helgason A, Sigurdardottir S, Nicholson J et al: Estimating Scandinavian and Gaelic ancestry in the male settlers of Iceland. Am J Hum Genet 2000; 67: 697–717. | Article | PubMed | ISI | ChemPort |
Bandelt HJ, Forster P, Rohl A: Median-joining networks for inferring intraspecific phylogenies. Mol Biol Evol 1999; 16: 37–48. | PubMed | ISI | ChemPort |
Nei M, Tajima F, Tateno Y: Accuracy of estimated phylogenetic trees from molecular data. II. Gene frequency data. J Mol Evol 1983; 19: 153–170. | Article | PubMed | ChemPort |
Nei M: Analysis of gene diversity in subdivided populations. Proc Natl Acad Sci USA 1973; 70: 3321–3323. | Article | PubMed | ChemPort |
Nei M, Kumar S: Molecular evolution and phylogenetics. New York: Oxford University Press, 2000.
Francalacci P, Morelli L, Underhill PA et al: Peopling of three Mediterranean islands (Corsica, Sardinia, and Sicily) inferred by Y-chromosome biallelic variability. Am J Phys Anthropol 2003; 121: 270–279. | Article | PubMed | ChemPort |
Flores C, Maca-Meyer N, Larruga JM et al: Isolates in a corridor of migrations: a high-resolution analysis of Y-chromosome variation in Jordan. J Hum Genet 2005; 50: 435–441. | Article | PubMed |
Cinnioglu C, King R, Kivisild T et al: Excavating Y-chromosome haplotype strata in Anatolia. Hum Genet 2004; 114: 127–148. | Article | PubMed |
Roewer L, Krawczak M, Willuweit S et al: Online reference database of European Y-chromosomal short tandem repeat (STR) haplotypes. Forensic Sci Int 2001; 118: 106–113. | Article | PubMed | ISI | ChemPort |
Semino O, Magri C, Benuzzi G et al: Origin, diffusion, and differentiation of Y-chromosome haplogroup E and J: inferences on the neolithization of Europe and later migratory events in the Mediterranean area. Am J Hum Genet 2004; 74: 1023–1034. | Article | PubMed | ISI | ChemPort |
Cruciani F, La Fratta R, Santolamazza PG et al: Phylogeographic analysis of haplogroup E3b (E-M215) Y chromosomes reveals multiple migratory events within and out of Africa. Am J Hum Genet 2004; 74: 1014–1022. | Article | PubMed | ISI | ChemPort |
Kayser M, Roewer L, Hedman M et al: Characteristics and frequency of germline mutations at microsatellite loci from the human Y chromosome, as revealed by direct observation in father/son pairs. Am J Hum Genet 2000; 66: 1580–1588. | Article | PubMed | ISI | ChemPort |
Zhivotovsky LA, Underhill PA, Cinnioglu C et al: The effective mutation rate at Y chromosome short tandem repeats, with application to human population-divergence time. Am J Hum Genet 2004; 74: 50–61. | Article | PubMed | ISI | ChemPort |
Cruciani F, La Fratta R, Torroni A et al: Molecular dissection of the Y chromosome haplogroup E-M78 (E3b1a): a posteriori evaluation of a microsatellite-network-based approach through six new biallelic markers. Hum Mutation 2006; 27: 831–832. | Article |
Engels DW: Alexander the Great and the logistics of the Macedonian Army. Berkeley, CA: University of California Press, 1981.

www.nature.com