Most of the extant mtDNA boundaries in South and Southwest Asia were likely shaped during the initial settlement of Eurasia by anatomically modern humans
© Metspalu et al; licensee BioMed Central Ltd. 2004
Received: 07 May 2004
Accepted: 31 August 2004
Published: 31 August 2004
Recent advances in the understanding of the maternal and paternal heritage of south and southwest Asian populations have highlighted their role in the colonization of Eurasia by anatomically modern humans. Further understanding requires a deeper insight into the topology of the branches of the Indian mtDNA phylogenetic tree, which should be contextualized within the phylogeography of the neighboring regional mtDNA variation. Accordingly, we have analyzed mtDNA control and coding region variation in 796 Indian (including both tribal and caste populations from different parts of India) and 436 Iranian mtDNAs. The results were integrated and analyzed together with published data from South, Southeast Asia and West Eurasia.
Four new Indian-specific haplogroup M sub-clades were defined. These, in combination with two previously described haplogroups, encompass approximately one third of the haplogroup M mtDNAs in India. Their phylogeography and spread among different linguistic phyla and social strata was investigated in detail. Furthermore, the analysis of the Iranian mtDNA pool revealed patterns of limited reciprocal gene flow between Iran and the Indian sub-continent and allowed the identification of different assemblies of shared mtDNA sub-clades.
Since the initial peopling of South and West Asia by anatomically modern humans, when this region may well have provided the initial settlers who colonized much of the rest of Eurasia, the gene flow in and out of India of the maternally transmitted mtDNA has been surprisingly limited. Specifically, our analysis of the mtDNA haplogroups, which are shared between Indian and Iranian populations and exhibit coalescence ages corresponding to around the early Upper Paleolithic, indicates that they are present in India largely as Indian-specific sub-lineages. In contrast, other ancient Indian-specific variants of M and R are very rare outside the sub-continent.
Two mtDNA macro-haplogroups (M and N) that arose from the African haplogroup L3 encompass virtually all mtDNAs outside Africa [1–4]. The phylogenetic node N (including R) has spread its branches all over Eurasia, in contrast to haplogroup M, which is found in Eastern Eurasia but is virtually absent in Europe. The numerous branches of N are, however, generally segregated to either the eastern (e.g. A, B , Y , R9  or western (e.g. N1 [8, 9], N2 (comprising of W and its sister-clade identified by ), TJ, HV, U ) Eurasian-specific pools.
The majority of Indian mtDNAs belong to macro-haplogroup M [8, 12–21]. While the topology of the M sub-haplogroups that are common in mainland East Asia (M7, M8 (including C, Z), M9 (including E), D, G [7, 22, 23]) and in Africa (M1 ) is established in detail, the internal haplogroup structure of M in India has remained largely undefined. We have previously demonstrated that transitions at nps. 477G, 1780, 8502 and 16319 designate Indian-specific haplogroup M2, the most frequent M clade in India . Another Indian-specific M clade supported by HVS-I variation as well as coding region markers, is M6 . Haplogroups M3, M4 and M5 have been discriminated preliminarily by their characteristic HVS-I mutations , but since their defining positions, 16126, 16311, and 16129, respectively, are phylogenetically unstable [25, 26], it is unlikely that the proposed haplogroups are monophyletic. Most numerous sub-groups of macro-haplogroup N in India are the Indian-specific variants of the phylogenetic node R including haplogroups R5, R6, U2(a, b, c) [8, 13, 27].
The overwhelming majority of the Iranian mtDNAs have been shown to lie in the West Eurasian domain of the global human mtDNA pool [27, 28]. Here we focus on the analysis of mtDNA lineages that are shared between Indians and Iranians and bear signals of pre-Holocene expansion in the region.
India congregates four linguistic domains (Indo-European, Dravidic, Austro-Asiatic and Tibeto-Burman) that occupy non-random spheres of the geographic distribution of its populations. The majority of the recent studies based on mtDNA variation have, in contrast to some , provided evidence that linguistic groups of India do not represent genetically homogeneous units and are not, therefore, traceable to different immigration waves from distinct sources [8, 13, 19]. The complexity that arises in defining populations and groups of populations in India based on genetic and cultural criteria has been recently demonstrated in South Indian tribal and caste populations. The combined data from mtDNA, Y-chromosome and autosomal genes indicated that the tribes and castes derive largely from the same genetic heritage of Late Pleistocene southern and southwestern Asians, and have received limited gene flow from external sources since the Holocene . Similar results were obtained by Cordaux et al. , who demonstrated that caste and tribal groups exhibit similar levels of molecular variance. However, genetic distances indicated that the Tibeto-Burman speaking tribal populations (from eastern India) were more closely related to East Asians than to other Indians . This is consistent with an earlier suggestion placing the origin of these tribal groups east of India – in Tibet and Myanmar .
Characteristics of the Indian and Iranian population samples whose mtDNA variation has been determined in the course of this study.
State / Region
Mixed caste people
Mixed caste people from Cochin
Mixed caste people
Mixed caste people
Uttar Pradesh Brahmin
ca. 68 Million
Results and Discussion
Geographic distribution of macro-clades M and N in India
Geographic, linguistic and socio-cultural distribution of major Indian-specific mtDNA haplogroups
HAPLOGROUP FREQUENCY (95% CR FOR PROPORTION)
R2; U7; W
Socio-cultural affiliation (Indian data only)
Brahmins & Kshatriyas
Language groups of India
Indian States and Bangladesh (excluding tribal populations) k
With the exception of the diverse set of largely Indian-specific R lineages, the most frequent mtDNA haplogroup in India that derives from the phylogenetic node N is haplogroup W . The frequency peak of haplogroup W is 5% in the northwestern states – Gujarat, Punjab and Kashmir. Elsewhere in India its frequency is very low (from 0 to 0.9%) (Table 2) forming a significant spatial cline (Figure 4).
MtDNA haplogroups in Iran
Over 90% of the mtDNAs found in Iran belong to haplogroups HV, TJ, U, N1, N2 and X, commonly found in West Eurasia (Table 2). In contrast to Europe, where H is predominant among the mtDNA haplogroups, in Iran the frequency of haplogroup U (29%) is higher than that of haplogroup H (17%) (Table 9, see Additional file 4). This difference accounts, at least partly for the presence in Iran of U sub-groups, such as U7 (9.4%), that are virtually absent in Europe.
Compared to India, haplogroup M frequency in Iran is marginally low (5.3%) and there are no distinguished Iranian-specific sub-clades of haplogroup M. All Iranian haplogroup M lineages can be seen as derived from other regional variants of the haplogroup: eleven show affiliation to haplogroup M lineages found in India, twelve in East and Central Asia (D, G, and M8) and one in northeast Africa (M1).
Indian-specific (R5 and Indian-specific M and U2 variants) and East Asian-specific (A, B and East Asian-specific M subgroups) mtDNAs, both, make up less than 4% of the Iranian mtDNA pool. We used Turkey (88.8 ± 4.0%) as the third parental population for evaluating the relative proportions of admixture from India (2.2 ± 1.7%) and China (9.1 ± 4.1%) into Iran. Therefore we can conclude that historic gene flow from India to Iran has been very limited.
The package of the most ancient mtDNA haplogroups in India
Indian-specific sub-clades of mtDNA haplogroups M and R.
DIAGNOSTIC CODING REGION MARKERS
ANCESTRAL HVS-I MOTIF
COALESCENCE b (years)
Proportion (95% CR)
70,600 ± 21,000
70,100 ± 20,700
48,300 ± 20,100
46,700 ± 22,800
54,800 ± 25,000
57,600 ± 22,300
17,300 ± 7,400
17,300 ± 7,600
3.0f (2.4 – 3.8)
19,200 ± 9,000
19,100 ± 9,000
33,000 ± 13,900
30,000 ± 13,600
19,100 ± 7,600
15,700 ± 8,100 d
6,000 ± 2,100
12,100 ± 4,500
9,400 ± 3,200
17,100 ± 4,700
19,400 ± 7,200
22,300 ± 8,600
1.0g (0.7 – 1.5)
40,400 ± 14,300
52,500 ± 21,700
66,100 ± 22,000
69,800 ± 24,800
30,000 ± 11,000
30,300 ± 11,600
Together with the Indian-specific sub-clades of haplogroup U2 [13, 27], haplogroups M2 and R5 can be discriminated as a package of Indian-specific mtDNA clades harboring extremely deep coalescence times (around 50,000 – 70,000 ybp). Together they constitute nearly 15% of the Indian mtDNAs. Importantly, these haplogroups are virtually absent elsewhere in Eurasia [13, 15], this study]. Because most of Indian varieties of haplogroup M are still unclassified (M*), this package is likely to be extended when more mtDNA coding region information will become available for the M* lineages in India.
The geographic distribution of the M2, R5 and U2i package seems to be rather uniform in the context of the Indian-specific maternal lineages (SAA p > 0.05). When excluding the mtDNAs that are likely to have arrived more recently from West or East Eurasia, the share of the package among the caste populations in northern and southern India is roughly similar (CR: 16.1–24.5% and 19.5–26.0 %, respectively). However, in accordance with expectations from the individual haplogroup distributions, the tribal groups speaking Austro-Asiatic and Tibeto-Burman languages are characterized by considerably lower values (CR: 3.9–14.1% and 3.8–15.1%, respectively).
The quest for finding the origin of haplogroup M and a plausible scenario for the peopling of Eurasia
Based on the high frequency and diversity of haplogroup M in India and elsewhere in Asia, some authors have suggested (versus ) that M may have arisen in Southwest Asia [16, 17, 31]. Finding M1 or a lineage ancestral to M1 in India, could help to explain the presence of M1 in Africa as a result of a back migration from India. Yet, to date this has not been achieved , this study). Therefore, one cannot rule out the still most parsimonious scenario that haplogroup M arose in East Africa . Furthermore, the lack of L3 lineages other than M and N (indeed, L3M and L3N) in India is more consistent with the African launch of haplogroup M. On the other hand, one also observes that: i) M1 is the only variant of haplogroup M found in Africa; ii) M1 has a fairly restricted phylogeography in Africa, barely penetrating into sub-Saharan populations, being found predominantly in association with the Afro-Asiatic linguistic phylum – a finding that appears to be inconsistent with the distribution of sub-clades of haplogroups L3 and L2 that have similar time depths. That, plus the presence of M1 without accompanying L lineages in the Caucasus  and [our unpublished data], leaves the question about the origin of haplogroup M still open.
The improved structure of autochthonous Indian mtDNA clades
Nearly a third of Indian mtDNAs belonging to haplogroup M could be assigned to its existing boughs and limbs with the current knowledge of the mtDNA coding region polymorphisms (Table 9, see Additional file 4). It is likely that the unclassified Indian M* and R* mtDNAs are also to a large extent autochthonous because neither the East nor West Eurasian mtDNA pools include such lineages at notable frequencies.
Coalescence estimates for these Indian-specific mtDNA haplogroups (M3a, M4a, M6, M25 and R6) fall largely between 20,000 and 30,000 ybp. These estimates overlap with those of many West Eurasian-specific (e.g. H, HV, preHV, U3, U4, K, X [9, 34]) and East Eurasian-specific (A, F2, D4, M7c1, M7a1, M8a [7, 22]) mtDNA clades, suggesting a rather synchronic worldwide demographic expansion event in the late Pleistocene, during an interglacial period preceding the LGM.
Several Indian-specific mtDNA clades demonstrate a similar spread-pattern in southern India. We found haplogroups M4a, M6a and M18 in southeastern Tamil Nadu and Andhra Pradesh while they were absent from neighboring Karnataka and Kerala (Figure 1 panel M6a and Figure 2 panels M4a and M18). One possible explanation is that admixture has been facilitated along the coastlines of the Arabian Sea and the Bay of Bengal. On the other hand, because the absolute frequencies of these haplogroups are rather low, it cannot be ruled out that an increase of sample sizes would disrupt the observed spread-pattern.
Were the Austro-Asiatic speaking tribal people the earliest inhabitants of India?
By calculating nucleotide diversities and expansion times (using the method from ) for different linguistic groups of India, some previous studies on mtDNA variation have distinguished the Austro-Asiatic speaking tribal groups as the carriers of the genetic legacy of the earliest settlers of the subcontinent [17, 38]. However, because the linguistic groups of India do not cluster into distinct branches of the Indian mtDNA tree [13, 15, 19], this study], calculating the beginning of expansion for those groupings is problematic and likely controversial as well.
Recently, Basu et al. (2003) supported the conclusions of [17, 38] by reporting that the frequency of the ancient haplogroup M2 among the Austro-Asiatic tribal populations is as high as 19%, and that they lack the slightly younger haplogroup M4 (the likely paraphyletic mother-clade of M4a). The authors have regarded the HVS-I transition at np 16319 as sufficient in defining haplogroup M2. This assumption, however, might lead to an overestimation of M2 frequency and age. Indeed, the 16319 transition has arisen several times on the background of other Indian haplogroup M lineages (Table 7, see Additional file 2), more specifically, in mtDNAs lacking the coding region markers that define M2 .
Although two out of the four M2 sequences reported by Basu et al. (2003) among the Lodha, Mundas and Santals (HVS-I sequences originally published by ) do harbor the characteristic M2a HVS-I motif (16223-16270-16274-16319-16352), without information from the coding region it is not clear whether the other two sequences (HVS-I motifs: 16092-16179-16223-16289-16294-16319 and 16147G-16172-16223-16319) represent novel M2 sub-clades (because these sequences cannot be affiliated with M2a or M2b) or derive from two independent branches of haplogroup M where 16319 transition has arisen recurrently. HVS-I motif 16147G-16172-16223, for example, is commonly associated with haplogroup N1a. Since sequence data on the five M2s among the Austro-Asiatic speaking tribe Ho, reported by Basu et al. (2003), have not been made available in the publication, we cannot rely on their haplogroup classification. Thus, we are left with one Munda and one Santal mtDNA belonging to haplogroup M2. They make up just 5% of the Austro-Asiatic tribal sample of 37 subjects (excluding the ten Ho). Interestingly, we found no instances of haplogroup M2 among the 56 Lodhas analyzed in this study. Consequently, when excluding the recurrences of the 16319 transition on the background of other sequence motifs, the frequency of M2 among the Austro-Asiatic speaking tribal groups from West Bengal in the combined dataset (Table 7, see Additional file 2) is significantly reduced to about 2%. The corrected value is comparable to the M2 frequency (>3%) in tribal populations speaking Indo-European languages of Punjab and Uttar Pradesh, but is significantly lower than its frequency (>14%) among the Dravidic-speaking tribal groups of Andhra Pradesh (Table 8, see Additional file 3).
Language families present today in India, such as Indo-European, Dravidic and Austro-Asiatic, are all much younger than the majority of indigenous mtDNA lineages found among their present-day speakers at high frequencies (see Additional file 9). It would make it highly speculative to infer, from the extant mtDNA pools of their speakers, whether one of the listed above linguistically defined group in India should be considered more "autochthonous" than any other in respect of its presence in the subcontinent.
Additionally, we note that some recent linguistic and archaeological evidence place the spread of the Austro-Asiatic languages in the Neolithic, in conjunction with the dispersal of rice cultivation from the Yangtze River basin . If this were the case, it would imply that the arrival of this linguistic phylum in India was not associated with female gene flow.
Gene flow from West Eurasia
In comparison to an overall frequency of 17% (CR 15.1–19.3%) among the caste populations, only 7% (CR 6.2 – 9.0%) of the mtDNAs from the tribal groups show affiliation to the West Eurasian haplogroups. The observed difference could be caused by differences in the extent of gene flow from the west to different social layers of the Indian society , and/or a more pronounced genetic drift among the tribal groups.
Coalescence estimates and diversity values for mtDNA haplogroups U7, W and R2 in India, Central Asia, Near and Middle East.
NEAR AND MIDDLE EAST
41,400 ± 15,800
41,200 ± 14,800
34,400 ± 15,500
37,900 ± 11,100
32,000 ± 8,700
27,400 ± 8,300
40,400 ± 14,300
36,100 ± 11,100
25,200 ± 13,300
Similarly to HV, TJ, N1, the spread of U7, W and R2 in India is geographically uneven – the three haplogroups are much more frequent in the northwestern states (SAA p < 0.05, Figure 4). Together, they constitute nearly 14% of the maternal gene pool in Indian Punjab, Gujarat and Rajasthan. Their frequency is also high in Iran (13%) and Pakistan (10%) but declines in Central Asia (5%).
Although haplogroup W is not highly frequent in European populations, it is nevertheless quite common , reaching the highest frequencies among the central-northern Finns (9% ). Yet, it is virtually absent from the Finno-Ugric speaking populations of the Volga basin . In Central Asia the frequency of haplogroup W stays below 2%.
Many European populations are lacking haplogroup U7 , but its frequency climbs over 4% in the Near East and up to 5% in Pakistan, reaching nearly 10% level in Iranians (Table 9, see Additional file 4). In India, haplogroup U7 frequency peaks at over 12% in Gujarat, in the westernmost state of India, while for the whole of India its frequency stays around 2%. Expansion times and haplotype diversities for the Indian and Near and Middle Eastern U7 mtDNAs are strikingly similar (Table 4). The possible homeland of this haplogroup spans likely in Gujarat and Iran because from there its frequency declines steeply both to the east and to the west. If the origin were in Iran rather than in India, then its equally high frequency as well as diversity in Gujarat favors a scenario whereby U7 has been introduced to the coastal western India either very early, or by multiple founders. Notably, the overlap of the Indian and Iranian lineages is largely restricted to the ancestral nodes while the coalescence age estimates for the nearly exclusively Indian (16207-16309-16318T) and West Eurasian (16126-16309-16318T) founder HVS-I motifs both yield time depths of about 20,000 to 30,000 ybp.
Haplogroup R2 appears at low frequencies in Near and Middle East and India and is virtually absent elsewhere. The spread of haplogroup R2 in Europe is restricted to a few populations in the Volga basin where it is represented by nodal haplotypes and by a region-specific subclade characterized by the HVS-I motif 16037–16172 . The coalescence estimate of this sub-clade is 11,400 ± 9,000 ybp. However, its wide error range prevents us from drawing any firm conclusions.
To summarize, the West and South Asian phylogeography of haplogroups W, U7 and R2 can be viewed as a genetic continuum that spans from the Near East into India, extending north into Central Asia. The coalescence times of these haplogroups suggest that this continuum took shape somewhere between 30,000 to 50,000 ybp (Table 4), thus falling within the climatically favorable interglacial period. We notice that the extant U7 and W frequencies along the proposed continuum are not uniform. U7 is more predominant in Iran, Pakistan, northwestern India and the Arabian peninsula, while W is more frequent in the western Near-East, Anatolia and the Caucasus. The coalescence ages of the Indian- and Iranian-specific U7 clades suggest that the time-window of this continuum was closed by ca. 20,000 ybp. The inferred extreme aridity of eastern Iran and western India during the last glacial maximum, which is well documented in paleovegetation reconstructions  may explain the observed segregation.
It has been suggested that the Jews settled in southwest India on the coast of the Arabian Sea sometime during the early Middle Ages. However, the mtDNA pool of the extant Cochin Jews is overwhelmingly Indian-specific (Table 10, see Additional file 5). We found exact or close matches to the fourteen HVS-I haplotypes observed among the Cochin Jews in other Indian populations. It is not clear whether the Near Eastern mtDNA lineages have been lost or the initial Jewish settlers did not include women.
Gene flow from East Eurasia
The East Eurasian-specific mtDNA haplogroups are less common in India and more sharply geographically segregated than the haplogroups of western Eurasian ancestry (Table 2; Figure 11, panel C). Indian caste populations harbor only about 4% of such mtDNAs, compared to 17% of the West Eurasian ones (Table 2). Elevated frequencies of haplogroups common in eastern Eurasia are observed in Bangladesh (17%) and Indian Kashmir (21%) and may be explained by admixture with the adjacent populations of Tibet and Myanmar (and possibly further east: from China and perhaps Thailand). On the other hand, the high frequencies of East Eurasian-specific mtDNAs found in the southern Indian state of Tamil Nadu (21%) are unexpected when considering their relatively low frequencies (~1%) in West Bengal and Andhra Pradesh. We notice, however, that the haplogroup assignments used here for the Tamil Nadu sample (A4, B4, F1a and M7) (Table 7, see Additional file 2) are based on HVS-I sequences alone . As shown and discussed elsewhere , such type of assignment is prone to mistakes.
Two varieties of haplogroup M, D4c and G2a, were recently identified as largely specific to Central Asia . In spite of geographical proximity we did not find these haplogroups in northern or northwestern India. Haplogroup G2a did, though in marginally low frequency, turn up in Iran (CR: 0.1 – 1.6%) and in southern Indian states Andhra Pradesh (CR: 0.3 – 2.9%) and Sri Lanka (CR: 0.2 – 4.1%).
Tibeto-Burman speaking tribal populations of eastern and northern India exhibit the highest frequencies of East Eurasian-specific mtDNA haplogroups. As inferred from the published HVS-I sequences , their share sums up to approximately two thirds of mtDNAs among the tribal groups in Assam, Nagaland, Arunachal Pradesh and Tripura, (Table 8, see Additional file 3; Figure 11, panel D). MtDNA haplogroups native to East Eurasia are also highly frequent in the northern states of India, reaching a peak of nearly 50% among the Kanet of Himachal Pradesh. Papiha and colleagues have previously demonstrated through the typing of immunoglobin allotypes that the Tibetan admixture among the regional Kanet groups decreases as the distance from the Tibetan border increases . Thus, mtDNA data are consistent with an ancestral origin of the Tibeto-Burman speaking tribal populations outside (east of) India in the neighboring Tibet and Myanmar [30, 45].
Haplotype sharing between populations
The majority (70%) of the 1136 mtDNA haplotypes found among continental Indians (including Pakistan and Bangladesh) are singletons and 41% of those that occur more than once are restricted to a single population. Only a few haplotypes are shared among five or more populations.
Proportion of mtDNA haplotype sharing between population groups of South Asia
1256 ± 215
617 ± 109
.399 – .453
1426 ± 215
643 ± 114
.396 – .447
132 (112 – 154)
.411 – .467
.324 – .372
107 (90 – 128)
.277 – .329
.345 – .394
96 (80 – 116)
.453 – .511
.331 – .379
107 (90 – 128)
.428 – 490
.273 – .317
90 (74 – 109)
.328 – .379
.393 – .446
106 (89 – 127)
.552 – .746
.014 – .025
9 (5 – 17)
.174 – .277
.091 – .115
24 (16 – 35)
.295 – .346
.425 – .487
87 (71 – 106)
An alternative method that assesses the degree of haplotype sharing between populations is to investigate the combined frequency of the shared haplotypes in two population groups. Thus, amongst the northern and the southern population groups the combined frequency of the haplotypes present also in the other group is significantly lower than that which we observed in the case of random groups. This is not surprising because West Eurasian-specific mtDNA haplogroups are rather frequent in northwest India. Because the Indo-European and the Dravidic speakers of India are largely concentrated to the northern and southern parts of the subcontinent, respectively, the differences arising from geographic division of the Indian populations also correspond to these linguistic groupings (Table 5).
Three Indian-specific haplogroups, M2, U2i and R5, which encompass about 15% of the Indian mtDNA pool, exhibit equally deep coalescence ages of about 50,000 – 70,000 years. Thus, their spread can be associated with the initial peopling of South Asia.
Haplogroups U7, W and R2 harbor a number of similar traits. Their overlapping geographic distributions and coalescence times suggest some degree of genetic continuum in the area spanning from the Near and Middle East through northwest India and reaching north into Central Asia. The coalescence estimates for these haplogroups are equally deep (around 30,000 – 50,000 years) in these different regions. That may be a result of either relatively more recent albeit large in scale migrations that brought along most of the diversity or may indeed reflect the region-specific expansions of these haplogroups. The former explanation could be ruled out since it is impossible to envisage a substantial movement of mtDNAs from South Asia that would not include haplogroup M. The same is true for the opposite – the share of U7, W and R2 within the West Eurasian-specific mtDNA haplogroups is two-fold higher in India than it is in Iran. Moreover, the South- and West Asian-specific sub-branches of haplogroup U7 predate the last glacial maximum. Therefore, deep autochthonous history of these haplogroups in the region remains to be the most parsimonious explanation.
Through the use of mtDNA coding region markers, we were able to classify altogether a quarter of the Indian M and R mtDNAs into a number of Indian-specific mtDNA haplogroups, four of which we newly identified. Several of these are characterized by clear patterns in their geographic distribution and/or different frequencies among different socio-cultural groups of India. Additional efforts should be undertaken to identify new coding region markers in order to further improve the classification of Indian mtDNAs. Indeed, much of the information is still hidden by the poor resolution of the Indian mtDNA tree and further piling of HVS-I datasets would add little to deepen our understanding of the demographic history of South and Southwest Asia.
We found that haplogroup M frequency drops abruptly from about 60% in India to about 5% in Iran, marking the western border of the haplogroup M distribution. A similarly sharp border cuts the distribution of Indian-specific mtDNA haplogroups to the east and to the north of the subcontinent. We therefore propose that the initial mtDNA pool established upon the peopling of South Asia has not been replaced but has rather been reshaped in situ by major demographic episodes in the past and garnished by relatively minor events of gene flow both from the West and the East during more recent chapters of the demographic history in the region.
MtDNA sequence variation in a total of 796 Indian samples most of which are held in a collection at Newcastle University  was analyzed. The samples cover a wide geographical range that spans from Himachal Pradesh in the north, Sri Lanka in the south, West Bengal in the east and Gujarat in the west (Table 1).
Tribal populations constitute 15% of the total sample size. The Lodha (n = 56) live mostly in the western part of Midnapore district of West Bengal where they are also known as Kheria or Kharia. Their total population size was ~59000 according to the 1981 census. Their language belongs to the Mundari branch of the Austro-Asiatic language family . The Kanet (n = 37) make up two thirds of the ~50000 inhabitants of the Kinnaur district in Himachal Pradesh . Their language belongs to the Himalayan group of the Tibeto-Burman language family . Five Bhoksa and twenty-six Tharu individuals were included in the present study in addition to those from the same populations that we have previously reported . Most of the Tharu live in the Terai areas (a belt of marshy land at the foot of the Himalayas) of Nepal (n = 720000). Approximately 96000 reside in Uttar Pradesh and Uttaranchal, Indian states adjacent to Nepal. Around 32000 Bhoksas live in the lowland areas of Uttaranchal and the neighboring Bijnor district of Uttar Pradesh. Both these tribal groups speak languages belonging to the Indo-European phylum .
The non-tribal Indian samples analyzed contained 105 West Bengalis of different caste rank, 58 Konkanastha Brahmins from Bombay, 53 Gujaratis, 50 Moors and 82 Sinhalese from Sri Lanka, 109 Punjabis of different caste rank from the Punjab, 25 Brahmins from Uttar Pradesh, 35 Rajputs from Rajasthan, 55 Parsees from Maharashtra and 100 subjects from Cochin, Kerala (including 45 Jews who have moved to Israel) (Table 1).
The Iranian sample of 436 individuals was collected in different locations mainly from southwestern and northwestern Iran (Table 1).
The new Indian mtDNA sequence data was combined with that previously published on Indian populations [3, 13, 15, 17, 18, 29, 35, 49–51] to produce a pooled dataset (n = 2572) in which the tribal populations constitute slightly over 50% (Table 6, see Additional file 1). By including also the data on Iranian and published data on Pakistani [13, 27, 29], Bangladeshi , Chinese [7, 22, 52–54] and Thai [53, 55, 56] populations (n = 145, 29, 919 and 552, respectively), a dataset of over 4600 samples spanning from West to East Asia (Table 6 and 7, see Additional files 1 and 2) was obtained. In many cases the published data was reanalyzed. Analogy with the newly obtained HVS-I motifs, that were classified into haplogroups using diagnostic markers of the mtDNA coding region, served as basis for haplogroup assignments of the published HVS-I sequences.
An even wider range of data was used for the production of the isofrequency maps: the published data from Armenia (n = 192)  and Kyrgyz (n = 95) , and the unpublished data from Yemen (n = 118), Oman (n = 79), Saudi Arabia (n = 205), Ethiopia (n = 270), Uzbekistan (n = 114) and Tajikistan (n = 41). Further details about the mtDNA variation of these populations will be published elsewhere.
MtDNA molecular analyses
DNA was extracted using standard phenol-chloroform methods . The hypervariable segment (HVS)-I (between nucleotide positions (np) 16024–16400) of the control region was sequenced in all the 796 Indian and 436 Iranian samples. Preparation of sequencing templates was carried out following Kaessmann et al. . Purified products were sequenced with the DYEnamic™ ET terminator cycle sequencing kit (Amersham Pharmacia Biotech) and analyzed on ABI 377 DNA or Megabace2000 Sequencers. Sequences were aligned and analyzed with the Wisconsin Package (GCG). In addition, informative mtDNA coding region positions [3, 11, 27] were assayed (Table 7, see Additional file 2) in selected individuals from different HVS-I haplotypes to determine haplogroup affiliations.
Median networks [60, 61] were constructed using Network 3.111 and Network 2.10B programs  with default settings (r = 2; ε = 0). MtDNA coding region markers were given five times the weight of the HVS-I positions. Coalescence of mtDNA haplogroups and sub-haplogroups was calculated using ρ (the averaged distance to a specified founder haplotype) and a mutation rate of one transition per 20,180 years between nps 16090–16365 . Standard errors for coalescence estimates and efficiency of a sample for coalescence time calculation (ρ/nδ2) were calculated following . Coalescence times were further calculated with the exclusion of tribal populations sequences, yet preserving the cluster topology (coalescence times marked with an asterisk). A more intense genetic drift (particularly founder effects) could introduce a bias into the coalescence time calculation (for example, see below the coalescence time of haplogroup M18 with and without the Lodha sample).
The software kindly provided by Vincent Macaulay was adopted in order to calculate the 95% credible regions (CR) from the posterior distribution of the proportion of a haplogroup/sub-haplogroup in the population.
Haplotype diversity was estimated as
where n is the number of sequences, k the number of distinct haplotypes, and ni number of sequences with one distinct haplotype.
Haplogroup isofrequency maps were generated using the Kriging method in Surfer 7 program of Golden Software. Haplogroup frequencies were averaged over populations from the same state in India, provinces in China and Thailand and geographic divisions in Iran. The data points for Kriging are shown as black dots, while the sample size applicable to the data point is given adjacent to the dot. In relatively small and isolated groups (e.g. tribal groups) random genetic drift might seriously affect the haplogroup frequencies, which may become uninformative when a whole region (e.g. state) is considered (e.g. M18 among the Lodha, see below). Therefore, the tribal data were excluded from the haplogroup isofrequency maps calculation. When illustrating the spread of mtDNA haplogroups native to West Eurasia, East Eurasia and India (Figure 11, panel D) we present these data as pie diagrams. The respective sample size and origin are indicated adjacent to the diagrams.
Spatial autocorrelation analysis was done using the PASSAGE software packet [65, 66]. The correlograms of Moran's I autocorrelation coefficient were calculated using binary weight matrix with five distance classes. Data on tribal populations was not used (see previous paragraph for explanation).
Haplogroup frequencies based admixture proportions were calculated using the ADMIX 2.0 software .
Haplotypes were defined as the HVS-I motif combined with the haplogroup label for haplotype analysis. In this way it was possible to discriminate between mtDNAs with identical HVS-I motifs but otherwise known to belong to different haplogroups (e.g. sequences with the CRS motif in HVS-I belonging to haplogroups H or R). In order to assess the results of haplotype sharing analysis between population groups (e.g. tribal and caste populations) we divided the Indian populations ten times randomly into two sets and analyzed the level of haplotype sharing between these sets.
We thank Lluis Quintana-Murci for providing data prior to publication, Vincent Macaulay for the algorithm for calculating credible regions for haplogroup frequencies, and Jaan Lind and Hille Hilpus for technical assistance. This work was supported by Estonian Science Foundation grants 514 (to RV), 5574 (to TK), 5807 (to EM), the Italian Ministry of the University: Progetti Ricerca Interesse Nazionale 2002 and 2003 (to AT), and European Commission grants ICA1CT20070006 and QLG2-CT-2002-90455 (to RV).
- Torroni A, Miller JA, Moore LG, Zamudio S, Zhuang J, Droma T, Wallace DC: Mitochondrial DNA analysis in Tibet: implications for the origin of the Tibetan population and its adaptation to high altitude. Am J Phys Anthropol. 1994, 93 (2): 189-199.View ArticlePubMedGoogle Scholar
- Chen YS, Torroni A, Excoffier L, Santachiara-Benerecetti AS, Wallace DC: Analysis of mtDNA variation in African populations reveals the most ancient of all human continent-specific haplogroups. Am J Hum Genet. 1995, 57 (1): 133-149.PubMed CentralPubMedGoogle Scholar
- Quintana-Murci L, Semino O, Bandelt H-J, Passarino G, McElreavey K, Santachiara-Benerecetti AS: Genetic evidence of an early exit of Homo sapiens sapiens from Africa through eastern Africa. Nat Genet. 1999, 23 (4): 437-441. 10.1038/70550.View ArticlePubMedGoogle Scholar
- Alves-Silva J, da Silva Santos M, Guimarães PE, Ferreira AC, Bandelt H-J, Pena SD, Prado VF: The ancestry of Brazilian mtDNA lineages. Am J Hum Genet. 2000, 67 (2): 444-461. 10.1086/303004.PubMed CentralView ArticlePubMedGoogle Scholar
- Torroni A, Schurr TG, Cabell MF, Brown MD, Neel JV, Larsen M, Smith DG, Vullo CM, Wallace DC: Asian affinities and continental radiation of the four founding Native American mtDNAs. Am J Hum Genet. 1993, 53 (3): 563-590.PubMed CentralPubMedGoogle Scholar
- Schurr TG, Sukernik RI, Starikovskaya YB, Wallace DC: Mitochondrial DNA variation in Koryaks and Itel'men: population replacement in the Okhotsk Sea-Bering Sea region during the Neolithic. Am J Phys Anthropol. 1999, 108 (1): 1-39. 10.1002/(SICI)1096-8644(199901)108:1<1::AID-AJPA1>3.0.CO;2-1.View ArticlePubMedGoogle Scholar
- Kivisild T, Tolk H-V, Parik J, Wang Y, Papiha SS, Bandelt H-J, Villems R: The emerging limbs and twigs of the East Asian mtDNA tree. Mol Biol Evol. 2002, 19 (10): 1737-1751.View ArticlePubMedGoogle Scholar
- Kivisild T, Kaldma K, Metspalu M, Parik J, Papiha SS, Villems R: The place of the Indian mitochondrial DNA variants in the global network of maternal lineages and the peopling of the Old World. In: Genomic Diversity. Edited by: Papiha SS, Deka R, Chakraborty R. 1999, Kluwer Academic/Plenum Publishers, 135-152.View ArticleGoogle Scholar
- Richards M, Macaulay V, Hickey E, Vega E, Sykes B, Guida V, Rengo C, Sellitto D, Cruciani F, Kivisild T, et al: Tracing European founder lineages in the Near Eastern mtDNA pool. Am J Hum Genet. 2000, 67 (5): 1251-1276.PubMed CentralView ArticlePubMedGoogle Scholar
- Derbeneva OA, Starikovskaia EB, Volod'ko NV, Wallace DC, Sukernik RI: [Mitochondrial DNA variation in Kets and Nganasans and the early peoples of Northern Eurasia]. Genetika. 2002, 38 (11): 1554-1560.PubMedGoogle Scholar
- Macaulay VA, Richards MB, Hickey E, Vega E, Cruciani F, Guida V, Scozzari R, Bonné-Tamir B, Sykes B, Torroni A: The emerging tree of West Eurasian mtDNAs: a synthesis of control-region sequences and RFLPs. Am J Hum Genet. 1999, 64 (1): 232-249. 10.1086/302204.PubMed CentralView ArticlePubMedGoogle Scholar
- Ballinger SW, Schurr TG, Torroni A, Gan YY, Hodge JA, Hassan K, Chen KH, Wallace DC: Southeast Asian mitochondrial DNA analysis reveals genetic continuity of ancient mongoloid migrations. Genetics. 1992, 130 (1): 139-152.PubMed CentralPubMedGoogle Scholar
- Kivisild T, Bamshad MJ, Kaldma K, Metspalu M, Metspalu E, Reidla M, Laos S, Parik J, Watkins WS, Dixon ME, et al: Deep common ancestry of Indian and western-Eurasian mitochondrial DNA lineages. Curr Biol. 1999, 9 (22): 1331-1334. 10.1016/S0960-9822(00)80057-3.View ArticlePubMedGoogle Scholar
- Kivisild T, Papiha SS, Rootsi S, Parik J, Kaldma K, Reidla M, Laos S, Metspalu M, Pielberg G, Adojaan M, et al: An Indian Ancestry: a key for understanding human diversity in Europe and beyond. In: Archaeogenetics: DNA and the population prehistory of Europe. Edited by: Renfrew C, Boyle K. 2000, Cambridge: McDonald Institute for Archaeological Research University of Cambridge, 267-279.Google Scholar
- Kivisild T, Rootsi S, Metspalu M, Mastana S, Kaldma K, Parik J, Metspalu E, Adojaan M, Tolk H-V, Stepanov V, et al: The genetic heritage of the earliest settlers persists both in Indian tribal and caste populations. Am J Hum Genet. 2003, 72: 313-332. 10.1086/346068.PubMed CentralView ArticlePubMedGoogle Scholar
- Roychoudhury S, Roy S, Dey B, Chakraborty M, Roy M, Roy B, Ramesh A, Prabhakaran N, Rani MVU, Vishwanathan HMM, et al: Fundamental genomic unity of ethnic India is revealed by analysis of mitochondrial DNA. Current Science. 2000, 79 (9): 1182-1192.Google Scholar
- Roychoudhury S, Roy S, Basu A, Banerjee R, Vishwanathan H, Usha Rani MV, Sil SK, Mitra M, Majumder PP: Genomic structures and population histories of linguistically distinct tribal groups of India. Hum Genet. 2001, 109 (3): 339-350. 10.1007/s004390100577.View ArticlePubMedGoogle Scholar
- Bamshad MJ, Watkins WS, Dixon ME, Jorde LB, Rao BB, Naidu JM, Prasad BV, Rasanayagam A, Hammer MF: Female gene flow stratifies Hindu castes. Nature. 1998, 395 (6703): 651-652. 10.1038/27103.View ArticlePubMedGoogle Scholar
- Bamshad M, Kivisild T, Watkins WS, Dixon ME, Ricker CE, Rao BB, Naidu JM, Prasad BV, Reddy PG, Rasanayagam A, et al: Genetic evidence on the origins of Indian caste populations. Genome Res. 2001, 11 (6): 994-1004. 10.1101/gr.GR-1733RR.PubMed CentralView ArticlePubMedGoogle Scholar
- Passarino G, Semino O, Bernini LF, Santachiara-Benerecetti AS: Pre-Caucasoid and Caucasoid genetic features of the Indian population, revealed by mtDNA polymorphisms. Am J Hum Genet. 1996, 59 (4): 927-934.PubMed CentralPubMedGoogle Scholar
- Basu A, Mukherjee N, Roy S, Sengupta S, Banerjee S, Chakraborty M, Dey B, Roy M, Roy B, Bhattacharyya N, et al: Ethnic India: a genomic view, with special reference to peopling and structure. Genome Res. 2003, 13 (10): 2277-2290. 10.1101/gr.1413403.PubMed CentralView ArticlePubMedGoogle Scholar
- Yao YG, Kong QP, Bandelt HJ, Kivisild T, Zhang YP: Phylogeographic differentiation of mitochondrial DNA in Han Chinese. Am J Hum Genet. 2002, 70 (3): 635-651. 10.1086/338999.PubMed CentralView ArticlePubMedGoogle Scholar
- Kong QP, Yao YG, Sun C, Bandelt HJ, Zhu CL, Zhang YP: Phylogeny of East Asian mitochondrial DNA lineages inferred from complete sequences. Am J Hum Genet. 2003, 73 (3): 671-676. 10.1086/377718.PubMed CentralView ArticlePubMedGoogle Scholar
- Passarino G, Semino O, Quintana-Murci L, Excoffier L, Hammer M, Santachiara-Benerecetti AS: Different genetic components in the Ethiopian population, identified by mtDNA and Y-chromosome polymorphisms. Am J Hum Genet. 1998, 62 (2): 420-434. 10.1086/301702.PubMed CentralView ArticlePubMedGoogle Scholar
- Hasegawa M, Di Rienzo A, Kocher TD, Wilson AC: Toward a more accurate time scale for the human mitochondrial DNA tree. J Mol Evol. 1993, 37 (4): 347-354. 10.1007/BF00178865.View ArticlePubMedGoogle Scholar
- Wakeley J: Substitution rate variation among sites in hypervariable region 1 of human mitochondrial DNA. J Mol Evol. 1993, 37 (6): 613-623. 10.1007/BF00182747.View ArticlePubMedGoogle Scholar
- Quintana-Murci L, Chaix R, Wells S, Behar D, Sayar H, Scozzari R, Rengo C, Al-Zahery N, Semino O, Santachiara-Benerecetti A, et al: Where West meets East: The complex mtDNA landscape of the Southwest and Central Asian corridor. Am J Hum Genet. 2004, 74: 827-845. 10.1086/383236.PubMed CentralView ArticlePubMedGoogle Scholar
- Kivisild T, Rootsi S, Metspalu M, Metspalu E, Parik J, Kaldma K, Usanga E, Mastana S, Papiha S, Villems R: The Genetics of Language and Farming Spread in India. In: Examining the farming/language dispersal hypothesis. Edited by: Bellwood P, Renfrew C. 2003, Cambridge: The McDonald Institute for Archaeological Research, 215-222.Google Scholar
- Cordaux R, Saha N, Bentley G, Aunger R, Sirajuddin S, Stoneking M: Mitochondrial DNA analysis reveals diverse histories of tribal populations from India. Eur J Hum Genet. 2003, 3: 253-264. 10.1038/sj.ejhg.5200949.View ArticleGoogle Scholar
- Guha B: The racial affinities of the people of India. In: Census of India, 1931, Part III Ethnographical. 1935, Simla: Government of India PressGoogle Scholar
- Richards M, Rengo C, Cruciani F, Gratrix F, Wilson J, Scozzari R, Macaulay V, Torroni A: Extensive female-mediated gene flow from sub-Saharan Africa into near eastern Arab populations. Am J Hum Genet. 2003, 72 (4): 1058-1064. 10.1086/374384.PubMed CentralView ArticlePubMedGoogle Scholar
- Tambets K, Kivisild T, Metspalu E, Parik J, Kaldma K, Laos S, Tolk H-V, Gölge M, Demirtas H, Geberhiwot T, et al: The topology of the maternal lineages of the Anatolian and Trans-Caucasus populations and the peopling of the Europe: some preliminary considerations. In: Archaeogenetics: DNA and the population prehistory of Europe. Edited by: Renfrew C, Boyle K. 2000, Cambridge: Cambridge University Press, 219-235.Google Scholar
- Underhill PA, Passarino G, Lin AA, Shen P, Mirazon Lahr M, Foley R, Oefner PJ, Cavalli-Sforza LL: The phylogeography of Y chromosome binary haplotypes and the origins of modern human populations. Ann Hum Genet. 2001, 65 (1): 43-62. 10.1046/j.1469-1809.2001.6510043.x.View ArticlePubMedGoogle Scholar
- Reidla M, Kivisild T, Metspalu E, Kaldma K, Tambets K, Tolk H, Parik J, Loogvali E, Derenko M, Malyarchuk B, et al: Origin and Diffusion of mtDNA Haplogroup X. Am J Hum Genet. 2003, 73 (6): 1178-1190. 10.1086/379380.PubMed CentralView ArticlePubMedGoogle Scholar
- Mountain JL, Hebert JM, Bhattacharyya S, Underhill PA, Ottolenghi C, Gadgil M, Cavalli-Sforza LL: Demographic history of India and mtDNA-sequence diversity. Am J Hum Genet. 1995, 56 (4): 979-992.PubMed CentralPubMedGoogle Scholar
- Herrnstadt C, Elson JL, Fahy E, Preston G, Turnbull DM, Anderson C, Ghosh SS, Olefsky JM, Beal MF, Davis RE, et al: Reduced-median-network analysis of complete mitochondrial DNA coding-region sequences for the major African, Asian, and European haplogroups. Am J Hum Genet. 2002, 70 (5): 1152-1171. 10.1086/339933.PubMed CentralView ArticlePubMedGoogle Scholar
- Slatkin M, Hudson RR: Pairwise comparisons of mitochondrial DNA sequences in stable and exponentially growing populations. Genetics. 1991, 129 (2): 555-562.PubMed CentralPubMedGoogle Scholar
- Majumder PP: Ethnic populations of India as seen from an evolutionary perspective. J Biosci. 2001, 26 (4 Suppl): 533-545.View ArticlePubMedGoogle Scholar
- Higham C: Languages and Farming Dispersals: Austroasiatic Languages and Rice Cultivation. In: Examining the farming/language dispersal hypothesis. Edited by: Bellwood P, Renfrew C. 2003, Cambridge: The McDonald Institute for Archaeological ResearchGoogle Scholar
- Finnilä S, Lehtonen MS, Majamaa K: Phylogenetic network for European mtDNA. Am J Hum Genet. 2001, 68 (6): 1475-1484. 10.1086/320591.PubMed CentralView ArticlePubMedGoogle Scholar
- Bermisheva M, Tambets K, Villems R, Khusnutdinova E: Diversity of mitochondrial DNA haplotypes in ethnic populations of the Volga-Ural region of Russia. Mol Biol (Mosk). 2002, 36 (6): 990-1001.View ArticleGoogle Scholar
- Global land environments since the last interglacial:. [http://members.cox.net/quaternary/.]
- Comas D, Plaza S, Wells RS, Yuldaseva N, Lao O, Calafell F, Bertranpetit J: Admixture, migrations, and dispersals in Central Asia: evidence from maternal DNA lineages. Eur J Hum Genet. 2004, 12 (6): 495-504. 10.1038/sj.ejhg.5201160.View ArticlePubMedGoogle Scholar
- Papiha SS, Chahal SM, Mastana SS: Variability of genetic markers in Himachal Pradesh, India: variation among the subpopulations. Hum Biol. 1996, 68 (5): 629-654.PubMedGoogle Scholar
- Passarino G, Semino O, Pepe G, Shrestha SL, Modiano G, Santachiara Benerecetti AS: MtDNA polymorphisms among Tharus of eastern Terai (Nepal). Gene Geogr. 1992, 6 (3): 139-147.PubMedGoogle Scholar
- Papiha SS: Genetic variation in India. Hum Biol. 1996, 68 (5): 607-628.PubMedGoogle Scholar
- Singh KS, ed: The Scheduled Tribes. 1997, Oxford: Oxford University PressGoogle Scholar
- Papiha SS, Schanfield MS, Chakraborty R: Immunoglobulin allotypes and estimation of genetic admixture among populations of Kinnaur District, Himachal Pradesh, India. Hum Biol. 1996, 68 (5): 777-794.PubMedGoogle Scholar
- Thangaraj K, Ramana GV, Singh L: Y-chromosome and mitochondrial DNA polymorphisms in Indian populations. Electrophoresis. 1999, 20 (8): 1743-1747. 10.1002/(SICI)1522-2683(19990101)20:8<1743::AID-ELPS1743>3.0.CO;2-V.View ArticlePubMedGoogle Scholar
- Prasad B, Ricker C, Watkins W, Dixon M, Rao B, Naidu J, Jorde L, Bamshad M: Mitochondrial DNA variation in Nicobarese Islanders. Hum Biol. 2001, 73 (5): 715-725.View ArticlePubMedGoogle Scholar
- Thangaraj K, Singh L, Reddy A, Rao V, Sehgal S, Underhill P, Pierson M, Frame I, Hagelberg E: Genetic affinities of the andaman islanders, a vanishing human population. Current Biology. 2003, 13 (2): 86-93. 10.1016/S0960-9822(02)01336-2.View ArticlePubMedGoogle Scholar
- Melton T, Clifford S, Martinson J, Batzer M, Stoneking M: Genetic evidence for the proto-Austronesian homeland in Asia: mtDNA and nuclear DNA variation in Taiwanese aboriginal tribes. Am J Hum Genet. 1998, 63 (6): 1807-1823. 10.1086/302131.PubMed CentralView ArticlePubMedGoogle Scholar
- Yao Y, Nie L, Harpending H, Fu Y, Yuan Z, Zhang Y: Genetic relationship of Chinese ethnic populations revealed by mtDNA sequence diversity. Am J Phys Anthropol. 2002, 118 (1): 63-76. 10.1002/ajpa.10052.View ArticlePubMedGoogle Scholar
- Horai S, Murayama K, Hayasaka K, Matsubayashi S, Hattori Y, Fucharoen G, Harihara S, Park KS, Omoto K, Pan IH: mtDNA polymorphism in East Asian Populations, with special reference to the peopling of Japan. Am J Hum Genet. 1996, 59 (3): 579-590.PubMed CentralPubMedGoogle Scholar
- Oota H, Settheetham-Ishida W, Tiwawech D, Ishida T, Stoneking M: Human mtDNA and Y-chromosome variation is correlated with matrilocal versus patrilocal residence. Nat Genet. 2001, 29 (1): 20-21. 10.1038/ng711.View ArticlePubMedGoogle Scholar
- Fucharoen G, Fucharoen S, Horai S: Mitochondrial DNA polymorphisms in Thailand. J Hum Genet. 2001, 46 (3): 115-125. 10.1007/s100380170098.View ArticlePubMedGoogle Scholar
- Comas D, Calafell F, Mateu E, Perez-Lezaun A, Bosch E, Martinez-Arias R, Clarimon J, Facchini F, Fiori G, Luiselli D, et al: Trading genes along the silk road: mtDNA sequences and the origin of Central Asian populations. Am J Hum Genet. 1998, 63 (6): 1824-1838. 10.1086/302133.PubMed CentralView ArticlePubMedGoogle Scholar
- Sambrook J, Fritsch EF, Maniatis T: Molecular cloning: a laboratory manual.,. 1989, Cold Spring Harbor, NY: Cold Spring Harbor Laboratory Press, 2Google Scholar
- Kaessmann H, Heissig F, von Haeseler A, Paabo S: DNA sequence variation in a non-coding region of low recombination on the human X chromosome. Nat Genet. 1999, 22 (1): 78-81. 10.1038/8785.View ArticlePubMedGoogle Scholar
- Bandelt H-J, Forster P, Sykes BC, Richards MB: Mitochondrial portraits of human populations using median networks. Genetics. 1995, 141 (2): 743-753.PubMed CentralPubMedGoogle Scholar
- Bandelt H-J, Forster P, Röhl A: Median-joining networks for inferring intraspecific phylogenies. Mol Biol Evol. 1999, 16 (1): 37-48.View ArticlePubMedGoogle Scholar
- Free Phylogenetic Network Analysis Shareware Software: NETWORK,. [http://www.fluxus-engineering.com.]
- Forster P, Harding R, Torroni A, Bandelt H-J: Origin and evolution of Native American mtDNA variation: a reappraisal. Am J Hum Genet. 1996, 59 (4): 935-945.PubMed CentralPubMedGoogle Scholar
- Saillard J, Forster P, Lynnerup N, Bandelt H-J, Nørby S: mtDNA variation among Greenland Eskimos: the edge of the Beringian expansion. Am J Hum Genet. 2000, 67 (3): 718-726. 10.1086/303038.PubMed CentralView ArticlePubMedGoogle Scholar
- Rosenberg MS: PASSAGE. Pattern Analysis, Spatial Statistics, and Geographic Exegesis. Version 1.1. In: Department of Biology, Arizona State University, Tempe, AZ. 2001Google Scholar
- Rosenberg M: Passage software packet:. [http://lsweb.la.asu.edu/rosenberg/Passage/.]
- Dupanloup I, Bertorelle G: Inferring admixture proportions from molecular data: extension to any number of parental populations. Mol Biol Evol. 2001, 18 (4): 672-675.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open-access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.