The Phylogeography of Y-Chromosome Haplogroup H1a1a-M82 Reveals the Likely Indian Origin of the European Romani Populations

Abstract

Linguistic and genetic studies on Roma populations inhabited in Europe have unequivocally traced these populations to the Indian subcontinent. However, the exact parental population group and time of the out-of-India dispersal have remained disputed. In the absence of archaeological records and with only scanty historical documentation of the Roma, comparative linguistic studies were the first to identify their Indian origin. Recently, molecular studies on the basis of disease-causing mutations and haploid DNA markers (i.e. mtDNA and Y-chromosome) supported the linguistic view. The presence of Indian-specific Y-chromosome haplogroup H1a1a-M82 and mtDNA haplogroups M5a1, M18 and M35b among Roma has corroborated that their South Asian origins and later admixture with Near Eastern and European populations. However, previous studies have left unanswered questions about the exact parental population groups in South Asia. Here we present a detailed phylogeographical study of Y-chromosomal haplogroup H1a1a-M82 in a data set of more than 10,000 global samples to discern a more precise ancestral source of European Romani populations. The phylogeographical patterns and diversity estimates indicate an early origin of this haplogroup in the Indian subcontinent and its further expansion to other regions. Tellingly, the short tandem repeat (STR) based network of H1a1a-M82 lineages displayed the closest connection of Romani haplotypes with the traditional scheduled caste and scheduled tribe population groups of northwestern India.

About the Authors:

Niraj Rai

Contributed equally to this work with: Niraj Rai, Gyaneshwer Chaubey, Rakesh Tamang, George van Driem

Affiliation: CSIR-Centre for Cellular and Molecular Biology, Uppal Road, Hyderabad, India

Gyaneshwer Chaubey

Contributed equally to this work with: Niraj Rai, Gyaneshwer Chaubey, Rakesh Tamang, George van Driem

* E-mail: [email protected] (GC); [email protected] (KT)

Affiliation: Evolutionary Biology Group, Estonian Biocentre, Tartu, Estonia

Rakesh Tamang

Contributed equally to this work with: Niraj Rai, Gyaneshwer Chaubey, Rakesh Tamang, George van Driem

Affiliation: CSIR-Centre for Cellular and Molecular Biology, Uppal Road, Hyderabad, India

Ajai Kumar Pathak

Affiliation: Department of Evolutionary Biology, Institute of Molecular and Cell Biology, University of Tartu, Tartu, Estonia

Vipin Kumar Singh

Affiliation: CSIR-Centre for Cellular and Molecular Biology, Uppal Road, Hyderabad, India

Monika Karmin

Affiliations Evolutionary Biology Group, Estonian Biocentre, Tartu, Estonia, Department of Evolutionary Biology, Institute of Molecular and Cell Biology, University of Tartu, Tartu, Estonia

Manvendra Singh

Affiliation: CSIR-Centre for Cellular and Molecular Biology, Uppal Road, Hyderabad, India

Deepa Selvi Rani

Affiliation: CSIR-Centre for Cellular and Molecular Biology, Uppal Road, Hyderabad, India

Sharath Anugula

Affiliation: CSIR-Centre for Cellular and Molecular Biology, Uppal Road, Hyderabad, India

Brijesh Kumar Yadav

Affiliation: CSIR-Centre for Cellular and Molecular Biology, Uppal Road, Hyderabad, India

Ashish Singh

Affiliation: CSIR-Centre for Cellular and Molecular Biology, Uppal Road, Hyderabad, India

Ramkumar Srinivasagan

Affiliation: CSIR-Centre for Cellular and Molecular Biology, Uppal Road, Hyderabad, India

Anita Yadav

Affiliation: CSIR-Centre for Cellular and Molecular Biology, Uppal Road, Hyderabad, India

Manju Kashyap

Affiliation: CSIR-Centre for Cellular and Molecular Biology, Uppal Road, Hyderabad, India

Sapna Narvariya

Affiliation: CSIR-Centre for Cellular and Molecular Biology, Uppal Road, Hyderabad, India

Alla G. Reddy

Affiliation: CSIR-Centre for Cellular and Molecular Biology, Uppal Road, Hyderabad, India

George van Driem

Contributed equally to this work with: Niraj Rai, Gyaneshwer Chaubey, Rakesh Tamang, George van Driem

Affiliation: Himalayan Languages Project, Institut für Sprachwissenschaft, Universität Bern, Bern, Switzerland

Peter A. Underhill

Affiliation: Department of Genetics, Stanford University School of Medicine, Stanford, California, United States of America

Richard Villems

Affiliations Evolutionary Biology Group, Estonian Biocentre, Tartu, Estonia, Department of Evolutionary Biology, Institute of Molecular and Cell Biology, University of Tartu, Tartu, Estonia, Estonian Academy of Sciences, Tallinn, Estonia

Toomas Kivisild

Affiliations Department of Evolutionary Biology, Institute of Molecular and Cell Biology, University of Tartu, Tartu, Estonia, Department of Biological Anthropology, University of Cambridge, Cambridge, United Kingdom

Lalji Singh

Affiliations CSIR-Centre for Cellular and Molecular Biology, Uppal Road, Hyderabad, India, Genome Foundation, Hyderabad, India, Banaras Hindu University, Varanasi, India

Kumarasamy Thangaraj

* E-mail: [email protected] (GC); [email protected] (KT)

Affiliation: CSIR-Centre for Cellular and Molecular Biology, Uppal Road, Hyderabad, India

Introduction

The Roma in England are traditionally known as gypsies because it was thought that they came from Egypt and were therefore ‘gypcians’. German Zigeuner, French tzigane and names in several other European languages derive from a designation for a Manichaean sect that practiced sorcery and soothsaying in the last centuries of the Byzantine Empire. On linguistic grounds, Grellmann [1] pointed out that the Roma must have originated in the Indian subcontinent. A detailed linguistic study by Pott [2] established that the various dialects spoken by the Roma derive specifically from North India. Roma populations are distributed widely within Europe including the Balkans and Scandinavia as well as throughout the Near East.

The name by which Roma designate themselves is Rroma (singular Rrom), whereby the double rr in Romani orthography represents a uvular ‘r’ [R] as opposed to an apical ‘r’ [r]. The autonym Rroma is held to be cognate with oma, a collective term for the ancient aboriginal populations of the Indian subcontinent. Many oma remained outcastes or tribals, whereas some were assimilated into the lower strata of the caste system by the Indo-European speaking Indians [3], [4].

The Roma route of migration has long been the object of linguistic study. Grierson [5] propagated the idea that the ‘Gipsy languages’ were of ‘Dardic origin’, but Turner [6] demonstrated that the Romani languages were not Dardic, but belonged to the same central Indo-Aryan subgroup as Hindi. The presence of Burushaski loans in Romani [7], the lack of Arabic loans and the presence of Dardic, Georgian, Ossetian, Armenian and mediaeval Greek loans [8] indicate that the Roma migrated to Europe by a northerly route, beginning around Gilgit in the northernmost Hindu Kush, thence along the southern Caspian littoral, the southern flank of the Caucasus, the southern shoreline of the Black Sea, across the Bosporus, and subsequently spreading across Europe since 13th century. A legacy of the migration is that some Roma refer to themselves as Sinti, an adjectival form derived from Sindh, the name of the Indus river. Morgenstierne [9] argued that the endangered omākī language spoken by several hundred oma in Gilgit and Yasin, belonging to the minstrel and blacksmith castes, represents an ethnolinguistic remnant of the early Roma migration through what today is northern Pakistan.

The Indian ancestry of Roma has been studied using mtDNA, Y-chromosomal and autosomal studies [10]–[16]. The South Asian-specific mtDNA haplogroups M5a1, M18 and M35b [13], [16], single Y-chromosomal haplogroup (hg) H1a1a-M82 [10], [17] and the pathogenic 1267delG mutation in CHRNE, i.e. cholinergic receptor, nicotinic, epsilon [18], root Roma ancestry in South Asia. It is also evident that the Roma exchanged a significant amount of genes with contemporary populations on their way to Europe [16], [17] (and references therein). The latter finding corroborates the admixture en route explicitly proposed by linguists in the 19th century [2], [19]. The genetic analysis of European Roma populations identified them as a suitable founder population to study Mendelian disorder and a valuable part of the European genetic landscape [12], [14], [15], [20], [21]. Recently, a study on primary congenital glaucoma reported a common mutation in the Roma and Jatt populations [21].

Although the Indian origin of the European Roma populations is linguistically and genetically well-established, accurate identification of their South Asian source has remained a matter of debate. Some linguistic studies argue that the proto-Romani founder population must have been in northwestern India [22]–[24], although their own origin myth suggests an origin on the Gangetic plain [25], [26]. The classical and mtDNA genetic markers suggested the closest affinity of the Roma with Rajput and Punjabi populations from Rajasthan and the Punjab respectively [16], [20], although these studies were compromised by low level phylogenetic resolution and a limited coverage of the Indian populations. These findings are now reconciled in light of our present genetic study thoroughly, comparing paternal lineages of Roma and a large number of South Asian populations, throwing light on the relationships and recent split of these groups.

In contrast with the maternal lineages, where three South Asian-specific lineages were reported [16], the paternal lineage of Roma carries only a single currently discernible South Asian-specific Y-chromosomal founder hg H1a1a-M82, ranging from 10 to 60% frequency in various Roma populations [10], [17], [27]–[31]. The paternal ancestor of hg H1a1a-M82, viz. hg H-M69, is an autochthonous component of the South Asian Y-chromosome gene pool, accounting, on average, for 22% of the studied paternal lineage pools [32]–[36]. Its relatively deep phylogenetic position within the overall Y-chromosome phylogeny as well as its extremely low frequency elsewhere, including the East and Southeast Asia, Middle East and Europe [32], [37]–[43], suggests that it arose in South Asia sometime prior to the Last Glacial Maximum, perhaps even earlier between 30 and 40 KYA ago [33].

There are examples where a single Y-chromosomal hg is associated with the expansion of a population from a homeland, e.g. hg C3 with Genghis Khan's army [44], spread of hg O1-M110 bearing Austronesians from Taiwan to the Admiralty Islands of Melanesia [45] and hg O2a-M95 with the incursion of Austroasiatic speakers into India [35]. The reduced diversity and expansion of H1a1a-M82 lineages in all Roma groups argues for a descent from a single paternal ancestor in the Indian subcontinent [10], [29]. The most popular model proposes an origin in northwestern India, but this model has not been tested using an adequately large collection of data from Y-chromosomes of the Indian subcontinent. Here, we address this issue by generating the largest dataset to date of H1a1a-M82 lineages including 214 ethnic populations from India that fill a major geographical gap (Table S1).

To detect possible source populations of European Roma paternal lineages within, we genotyped first M82-H1a1a hg) in a large number of Indian populations to produce a high-resolution dataset of hg H1a1a-M82 and thereafter, 17 Y-chromosome short tandem repeats (STRs) (Text S1 and Table S2). To reduce sampling bias, we have analysed more than 7000 samples from 214 ethnic populations covering the whole length and breadth of India (Table S1 and S3).

Results and Discussion

Hg H1a1a-M82 accounts for more than 12% of all male lineages in South Asia, with a highest incidence of 20% amongst South Indians and a lowest incidence of 0.2% among northeast Indian populations (Table S1 and S3). Hg H1a1a-M82 has a decreasing gradient from its peak frequency in southern India toward both the northwestern and the northeastern peripheries and a virtual absence further to the extreme east and west (Table S1, S3 and S4). Interestingly, substantial frequency changes occur despite the short geographical distance between the East Indian and their adjacent northeast Indian populations, suggesting different paternal population histories for these two neighboring regions. Consistent with the frequency distribution of H1a1a-M82 within India, the Y-STR variance and expansion ages basically reflect a similar pattern (Fig. 1 and Table 1). The overall TMRCA estimate of the H1a1a-M82 lineage in India using ρ-statistics [46], was ∼22KYA. The regional expansion times within India range between 16 and 24 KYA (Table 1). The higher expansion time with the associated mean pair-wise difference (MPD) and haplotype diversity (HD) (Fig. 1 and, Table 1 and 2) suggests that hg H1a1a-M82 originated in the same South Asian pre-LGM gene pool of Y-chromosomes.

[Figure omitted. See PDF.]

Figure 1. The most parsimonious route of prehistoric expansion of Y-chromosomal haplogroup H1a1a-M82 and the recent out-of -India migration of European Roma ancestors.

https://doi.org/10.1371/journal.pone.0048477.g001

[Figure omitted. See PDF.]

Table 1. Measure of Genetic Diversity and expansion ages of haplogroup H1a1a-M82 among different population groups.

https://doi.org/10.1371/journal.pone.0048477.t001

[Figure omitted. See PDF.]

Table 2. The average number of pairwise differences are shown within (PiX along diagonal) and between (PiXY above diagonal) populations; pairwise (δμ)2 genetic distance values are depicted below the diagonal.

https://doi.org/10.1371/journal.pone.0048477.t002

Within India, language group-wise expansion times and variance were observed to be highest among Dravidian speakers (25±3.5 KYA, 0.59), followed by Austroasiatic speakers (22.4±3.3 KYA, 0.52), and lowest among Indo-European speakers (19.5±2.2 KYA, 0.47). Just as with the expansion time of hg O2a-M95 among Austroasiatic speakers [35], the H1a1a-M82 expansion time was established to be substantially higher among South Munda (22±3.5 KYA) than in the North Munda (16.3±3 KYA) group. The South Munda groups live in close proximity to the Dravidian groups, and it is likely that the already diverse Dravidian H1a1a-M82 might have assimilated into the Munda group after the Austroasiatic incursion into the Indian subcontinent.

The geographical distribution of hg H1a1a-M82 is largely restricted to South Asia, and its significant occurrence among European Roma populations strongly links the Roma to the Indian subcontinent. The high frequency of H1a1a-M82 among all the Roma groups and their reduced genetic diversity relative to South Asian populations can most likely be attributed to their recent migration from India (Table 1 and 2). In the network analysis (Fig. 2), the Roma Y-STR haplotypes cluster predominantly close to the northwestern Indian haplotypes. In addition, the northwestern Indian haplotypes, while diverse, generally radiate from the core of the network while the Roma haplotypes being distributed further away. These patterns point to northwestern India as the source of the Roma H1a1a-M82 chromosomes. The average age estimate of Roma founders considering their distance from Northwest Indian founders is 1405±688 YBP (Table 3), which is largely in agreement with the time frames suggested previously.

[Figure omitted. See PDF.]

Figure 2. Phylogenetic network relating Y-STR haplotypes within haplogroup H1a1a -M82.

The network was constructed using a median joining with MP (maximum parsimony) algorithm as implemented in the Network 4.6 program. The size of the circles is proportional to the number of samples. The data used for comparison have been taken from the literature (See Text S1).

https://doi.org/10.1371/journal.pone.0048477.g002

[Figure omitted. See PDF.]

Table 3. Y-chromosomal haplogroup H1a1a-M82 founder analysis for Roma.

https://doi.org/10.1371/journal.pone.0048477.t003

In order to ascertain the closest population group among northwestern Indians, we redrew the network of Roma haplotypes exclusively within the northwestern Indian variation (Fig. 3). It is highly revealing that the closest or matching haplotypes with the Roma haplotypes were found in scheduled caste and scheduled tribe populations, while the middle and upper caste haplotypes were more distant to the Roma haplotypes (Fig. 3). Scheduled castes and Scheduled tribes are the endogamous groups in India that are given a special status by the Government of India to uplift their social status (for more details, refer [47]). Historically, the assimilation of so-called tribals into the caste system generally did little to ameliorate the socio-economic barriers or enhance the marriageability of former outcastes to members of the middle or high castes. However their language and means of subsistence were often affected, e.g. assimilation to an Indo-Aryan language and the shift from foraging, hunting and fishing to a more sedentary existence. Not surprisingly, the genetic differences between scheduled tribes and scheduled castes are not found to be substantial [47]. On the basis of our findings, it is therefore most parsimonious to conclude that the genealogically closest patrilineal ancestors of the Roma were among the ancestors of the present scheduled tribes and scheduled caste populations of northwestern India. The genetic data analysed here for the first time provide strong population genetic support for the linguistic based identification of the ancestral Roma with the presumed aboriginal oma of northwestern India and the Gangetic plain.

[Figure omitted. See PDF.]

Figure 3. Phylogenetic network relating Y-STR haplotypes within haplogroup H1a1a -M82- of Roma with northwest Indian populations.

The network was constructed using a median joining with MP (maximum parsimony) algorithm as implemented in the Network 4.6 program. The size of the circles is proportional to the number of samples. SC = Scheduled Caste and ST = Scheduled Tribe.

https://doi.org/10.1371/journal.pone.0048477.g003

This suggestion also enables a resolution of the controversy on the Punjab and the Gangetic plain. The oma are indigenous to the Punjab as well as to the entire Gangetic plain as far east as the lower course of the Brahmaputra. Yet not all oma groups were historically ancestral to the European Roma, and in fact most descendants of the aboriginal oma reside in the Indian subcontinent today [4]. It is held that a large number of aboriginal oma were recruited in the Punjab to repel the Ghaznavid invasions of the kingdom of Jayapāla between 1001 and 1026. These oma were rewarded by nominal promotion to Katriya ‘warrior’ or Rājpūt caste status, but such nominal promotions in the social context of the Indian subcontinent generally did not lead to a genuine enhancement of eligibility in marriage or social status because of the enduring nature of local memory [4]. With the fall of the Hindu polities in what today is Pakistan, the westward migration of the oma was set into motion. Arguably, the humble social position of the oma could frequently have led to circumstances, which could have promoted their geographical mobility. Yet the timing of the Ghaznavid invasion presents the probable temporal and spatial historical correlation.

In order to conduct an intra-group comparison of Y-chromosome variation, we estimated the genetic distances (RST values) from 15 loci Y-STR haplotypes and made a multidimensional scaling plot (Fig. S1). The first dimension largely separated all the Indian groups from one another, while the second dimension had a maximum impact in separating east Indian from north-central Indian and northwestern Indian from the Roma groups as well as Croatian Roma from Portuguese and Serbian Roma. We also noted other significant genetic differences between Croatian Roma and Portuguese and Serbian Roma. Portuguese and Serbian Roma share a closer affinity than Croatian Roma (Fig. S1 and Table S5). South and south-central Indian groups occupy a pivotal position in the plot, with other groups scattered around them in accordance with their geographical affinity. The modal haplotypes of different population groups are given in Table S6 and their distance to it in Table S7. The Analysis of Molecular Variance (AMOVA) based on Y-STR data shows that the Roma populations are closest to northwestern Indian populations (Table 4). The average mutational distance from the Roma modal haplotype shows a consistent pattern of northwestern Indian populations as representing the closest Indian groups to European Roma (Table S7).

[Figure omitted. See PDF.]

Table 4. Analysis of Molecular Variance (AMOVA) using Y-STRs between groups of populations categorized on the basis of geography.

https://doi.org/10.1371/journal.pone.0048477.t004

In conclusion, the analysis of Y-chromosome hg H1a1a-M82 variation in 214 ethnic groups from India shows that northwest Indian populations are the closest to the hg H1a1a-M82 variants observed in the present-day European Roma populations. Although only entire genome analyses can provide the richest narrative of history and migration, the Y-chromosome hg H1a provides an exceptional record of Indian and European Roma-specific paternal heritage, including their exodus from northwestern India and subsequent recent expansion. Our phylogeographical study of hg H1a1a-M82 enables now to better understand the temporal and spatial parameters of this migration. This first genetic evidence of this nature allows us to develop a more detailed picture of the paternal genetic history of European Roma, revealing that the ancestors of present scheduled tribes and scheduled caste populations of northern India, traditionally referred to collectively as the oma, are the likely ancestral populations of modern European Roma. Our findings corroborate the hypothesized cognacy of the terms Rroma and oma and resolve the controversy about the Gangetic plain and the Punjab in favour of the northwestern portion of the diffuse widespread range of the oma ancestral population of northern India.

Materials and Methods

Samples were collected with the informed written consent from 3498 unrelated healthy individuals belonging to 57 populations from all the four linguistic groups of India (Table S1 and S3). Y-chromosome marker M82 (defining hg H1a) was genotyped in all the samples. From one to five M82 derived samples of each population were randomly selected for Y-STR genotyping based upon the frequency of H1a1a-M82 in the respective population. Our main aim was to cover highest geographic area (Table S1 and S3). In total, 204, M82 derived samples were genotyped by using the AmpFℓSTR® Y-filer™ PCR amplification Kit (Applied Biosystems). The detailed material and method section is given in Text S1.

Supporting Information

[Figure omitted. See PDF.]

Text S1.

Detailed materials and methods .

https://doi.org/10.1371/journal.pone.0048477.s001

(DOC)

Figure S1.

Multidimensional scaling plot of the Rst distances of the population groups based on haplogroup H1a1a-M82 Y-STR data.

https://doi.org/10.1371/journal.pone.0048477.s002

(TIF)

Table S1.

Details of the Indian and Nepali samples included in the present study and haplogroup H1a1a-M82 frequencies.

https://doi.org/10.1371/journal.pone.0048477.s003

(DOC)

Table S2.

The 15 loci Y-STR profile of haplogroup H1a1a-M82 belonging to Indian, Afghani and European Roma population, used in the present analysis.

https://doi.org/10.1371/journal.pone.0048477.s004

(DOC)

Table S3.

Regionwise haplogroup frequency in India and Nepal.

https://doi.org/10.1371/journal.pone.0048477.s005

(DOC)

Table S4.

The frequency of haplogroup H1a1a-M82 among different world populations.

https://doi.org/10.1371/journal.pone.0048477.s006

(DOC)

Table S5.

Mean pairwise Fst between different studied groups for haplogroup H1a1a-M82.

https://doi.org/10.1371/journal.pone.0048477.s007

(DOC)

Table S6.

Modal H1a1a-M82 Y-STR haplotype of different population groups.

https://doi.org/10.1371/journal.pone.0048477.s008

(DOC)

Table S7.

Average mutational distances from Roma Modal haplotype.

https://doi.org/10.1371/journal.pone.0048477.s009

(DOC)

Acknowledgments

We thank all the donors for participating in this study. We are grateful to the anonymous reviewers for their constructive suggestions in improving the quality of our manuscript.

Author Contributions

Conceived and designed the experiments: GC KT. Performed the experiments: NR GC RT AKP VKS M. Karmin MS DSR SA BKY AS RS AY M. Kashyap SN AGR. Analyzed the data: NR GC RT AKP. Contributed reagents/materials/analysis tools: AGR RV LS KT. Wrote the paper: GC KT NR RT GVD. Provided input on manuscript writing: TK PAU RV M. Karmin.

Citation: Rai N, Chaubey G, Tamang R, Pathak AK, Singh VK, Karmin M, et al. (2012) The Phylogeography of Y-Chromosome Haplogroup H1a1a-M82 Reveals the Likely Indian Origin of the European Romani Populations. PLoS ONE 7(11): e48477. https://doi.org/10.1371/journal.pone.0048477

References

1. Grellman H (1783) Historischer Versuch über die Zigeuner betreffend die Lebensart und Verfassung, Sitten und Schicksale dieses Volks seit seiner Erscheinung in Europa, und dessen Ursprung: Goöttingen, bey Johann Christian Dieterich. 183p

2. Pott A (1844) Die Zigeuner in Europa und Asien, Erster Theil. Einleitung und Grammatik. Halle: Edouard Heynemann. 540p.

3. Trail G (1828) Statistical sketch of Kamaon. Asiatick Researches 16: 137–234.

4. van Driem G (2001) Languages of the Himalayas: An Ethnolinguistic Handbook of the Greater Himalayan Region, containing an Introduction to the Symbiotic Theory of Language, Leiden: Brill. 1401p.

5. Grierson G (1922) Linguistic Survey of India. Gypsy languages, Calcutta: Superintendent of Government Printing. 213p.

6. Turner R (1926) The position of Romani in Indo-Aryan. Journal of the Gypsy Lore Society 3: 145–94.

7. Berger H (1959) Die Burušaski-Lehnwörter in der Zigeunersprache. Indo-Iranian Journal 3: 17–43.

8. Hancock I (1995) A Handbook of Vlax Romani. Ohio: Slavic Publishers. 178p.

9. Morgenstierne G (1973) Irano-Dardica. University of Michigan: Dr. Ludwig Reichert Verlag. 388p.

10. Gresham D, Morar B, Underhill PA, Passarino G, Lin AA, et al. (2001) Origins and divergence of the Roma (Gypsies). Am J Hum Genet 69: 1314–1331.

11. Zhivotovsky LA, Underhill PA, Cinnioglu C, Kayser M, Morar B, et al. (2004) The effective mutation rate at Y-chromosome short tandem repeats, with application to human population-divergence time. Am J Hum Genet 74: 50–61.

12. Kalaydjieva L, Morar B, Chaix R, Tang H (2005) A newly discovered founder population: the Roma/Gypsies. Bioessays 27: 1084–1094.

13. Malyarchuk BA, Perkova MA, Derenko MV, Vanecek T, Lazur J, et al. (2008) Mitochondrial DNA variability in Slovaks, with application to the Roma origin. Ann Hum Genet 72: 228–240.

14. Mihaylova V, Hantke J, Sinigerska I, Cherninkova S, Raicheva M, et al. (2007) Highly variable neural involvement in sphingomyelinase-deficient Niemann-Pick disease caused by an ancestral Gypsy mutation. Brain 130: 1050–1061.

15. Bouwer S, Angelicheva D, Chandler D, Seeman P, Tournev I, et al. (2007) Carrier rates of the ancestral Indian W24X mutation in GJB2 in the general Gypsy population and individual subisolates. Genet Test 11: 455–458.

16. Mendizabal I, Valente C, Gusmão A, Alves C, Gomes V, et al. (2011) Reconstructing the Indian origin and dispersal of the European Roma: a maternal genetic perspective. PLoS ONE 6: e15988.

17. Klarić IM, Salihović MP, Lauc LB, Zhivotovsky LA, Rootsi S, et al. (2009) Dissecting the molecular architecture and origin of Bayash Romani patrilineages: genetic influences from South-Asia and the Balkans. Am J Phys Anthropol 138: 333–342.

18. Abicht A, Stucka R, Karcagi V, Herczegfalvi A, Horváth R, et al. (1999) A common mutation (epsilon1267delG) in congenital myasthenic patients of Gypsy ethnic origin. Neurology 53: 1564–1569.

19. Crawfurd J (1865) “On the Origin of the Gypsies”. Transactions of the Ethnological Society of London: Royal Anthropological Institute of Great Britain and Ireland. 25–36p

20. Mastana SS, Papiha SS (1992) Origin of the Romany gypsies–genetic evidence. Z Morphol Anthropol 79: 43–51.

21. Ali M, McKibbin M, Booth A, Parry DA, Jain P, et al. (2009) Null mutations in LTBP2 cause primary congenital glaucoma. Am J Hum Genet 84: 664–671.

22. Hancock IF (2002) We Are the Romani People. United Kingdom: University Of Hertfordshire Press. 288p.

23. Fraser A (1992) The Gypsies. Oxford: Blackwell Publishers Ltd. 359p.

24. Iovita RP, Schurr TG (2004) Reconstructing the Origins and Migrations of Diasporic Populations: The Case of the European Gypsies. American Anthropologist 106: 267–281.

25. Chaubey G (2010) The demographic history of India: A perspective based on genetic evidence. Universitatis Tartuensis, Estonia. Ph. D Thesis.185p. (http://hdl.handle.net/10062/15240)

26. Moreau R (2002) The Rom: walking in the paths of the Gypsies. Firefly Books 301p.

27. Pericic M, Lauc LB, Klaric IM, Rootsi S, Janicijevic B, et al. (2005) High-Resolution Phylogenetic Analysis of Southeastern Europe Traces Major Episodes of Paternal Gene Flow Among Slavic Populations. Mol Biol Evol 22: 1964–1975.

28. Gusmão A, Gusmão L, Gomes V, Alves C, Calafell F, et al. (2008) A perspective on the history of the Iberian Gypsies provided by phylogeographic analysis of Y-chromosome lineages. Ann Hum Genet 72: 215–227.

29. Pamjav H, Zalán A, Béres J, Nagy M, Chang YM (2011) Genetic structure of the paternal lineage of the Roma people. Am J Phys Anthropol 145: 21–29.

30. Zalán A, Béres J, Pamjav H (2011) Paternal genetic history of the Vlax Roma. Forensic Sci Int: Genet 5: 109–113.

31. Regueiro M, Stanojevic A, Chennakrishnaiah S, Rivera L, Varljen T, et al. (2011) Divergent patrilineal signals in three Roma populations. Am J Phys Anthropol 144: 80–91.

32. Sengupta S, Zhivotovsky LA, King R, Mehdi SQ, Edmonds CA, et al. (2006) 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 78: 202–221.

33. Trivedi R, Sahoo S, Singh A, Bindu G, Banerjee J, et al. (2008) Genetic Imprints of Pleistocene origin of Indian populations: a comprehensive phylogeographic sketch of Indian Y-chromosomes. Int J Hum Genet 8: 97–118.

34. Kumar V, Reddy ANS, Babu JP, Rao TN, Langstieh BT, et al. (2007) Y-chromosome evidence suggests a common paternal heritage of Austro-Asiatic populations. BMC Evol Biol 7: 47.

35. Chaubey G, Metspalu M, Choi Y, Mägi R, Romero IG, et al. (2011) Population Genetic Structure in Indian Austroasiatic speakers: the role of landscape barriers and sex-specific admixture. Mol Biol Evol 28: 1013–1024.

36. Debnath M, Palanichamy MG, Mitra B, Jin JQ, Chaudhuri TK, et al. (2011) Y-chromosome haplogroup diversity in the sub-Himalayan Terai and Duars populations of East India. J Hum Genet 56: 765–771.

37. Underhill PA, Shen P, Lin AA, Jin L, Passarino G, et al. (2000) Y-chromosome sequence variation and the history of human populations. Nat Genet 26: 358–361.

38. Gayden T, Cadenas AM, Regueiro M, Singh NB, Zhivotovsky LA, et al. (2007) The Himalayas as a directional barrier to gene flow. Am J Hum Genet 80: 884–894.

39. Cadenas AM, Zhivotovsky LA, Cavalli-Sforza LL, Underhill PA, Herrera RJ (2008) Y-chromosome diversity characterizes the Gulf of Oman. Eur J Hum Genet 16: 374–386.

40. Battaglia V, Fornarino S, Al-Zahery N, Olivieri A, Pala M, et al. (2009) Y-chromosomal evidence of the cultural diffusion of agriculture in southeast Europe. Eur J Hum Genet 17: 820–830.

41. Zhong H, Shi H, Qi XB, Duan ZY, Tan PP, et al. (2011) Extended Y-chromosome investigation suggests postglacial migrations of modern humans into East Asia via the northern route. Mol Biol Evol 28: 717–727.

42. He JD, Peng MS, Quang HH, Dang KP, Trieu AV, et al. (2012) Patrilineal perspective on the Austronesian diffusion in mainland southeast Asia. PLoS ONE 7: e36437.

43. Grugni V, Battaglia V, Hooshiar Kashani B, Parolo S, Al-Zahery N, et al. (2012) Ancient migratory events in the middle East: new clues from the y-chromosome variation of modern iranians. PLoS ONE 7: e41252.

44. Zerjal T, Xue Y, Bertorelle G, Wells RS, Bao W, et al. (2003) The genetic legacy of the Mongols. Am J Hum Genet 72: 717–721.

45. Kayser M, Choi Y, van Oven M, Mona S, Brauer S, et al. (2008) The impact of the Austronesian expansion: evidence from mtDNA and Y-chromosome diversity in the Admiralty Islands of Melanesia. Mol Biol Evol 25: 1362–1374.

46. Forster P, Harding R, Torroni A, Bandelt HJ (1996) Origin and evolution of Native American mtDNA variation: a reappraisal. Am J Hum Genet 59: 935–945.

47. Chaubey G, Metspalu M, Kivisild T, Villems R (2007) Peopling of South Asia: investigating the caste-tribe continuum in India. Bioessays 29: 91–100.