Austro-Asiatic Tribes of Northeast India Provide Hitherto Missing Genetic Link between South and Southeast Asia

Abstract

Northeast India, the only region which currently forms a land bridge between the Indian subcontinent and Southeast Asia, has been proposed as an important corridor for the initial peopling of East Asia. Given that the Austro-Asiatic linguistic family is considered to be the oldest and spoken by certain tribes in India, Northeast India and entire Southeast Asia, we expect that populations of this family from Northeast India should provide the signatures of genetic link between Indian and Southeast Asian populations. In order to test this hypothesis, we analyzed mtDNA and Y-Chromosome SNP and STR data of the eight groups of the Austro-Asiatic Khasi from Northeast India and the neighboring Garo and compared with that of other relevant Asian populations. The results suggest that the Austro-Asiatic Khasi tribes of Northeast India represent a genetic continuity between the populations of South and Southeast Asia, thereby advocating that northeast India could have been a major corridor for the movement of populations from India to East/Southeast Asia.

About the Authors:

B. Mohan Reddy

*To whom correspondence should be addressed. E-mail: [email protected]

Affiliation: Biological Anthropology Unit, Indian Statistical Institute, Hyderabad, India

B. T. Langstieh

Affiliations Biological Anthropology Unit, Indian Statistical Institute, Hyderabad, India, Department of Anthropology, Northeast Hill University, Shillong, India

Vikrant Kumar

Affiliations Biological Anthropology Unit, Indian Statistical Institute, Hyderabad, India, Genome Institute of Singapore, Singapore, Singapore

T. Nagaraja

Affiliation: Biological Anthropology Unit, Indian Statistical Institute, Hyderabad, India

A. N. S. Reddy

Affiliation: Biological Anthropology Unit, Indian Statistical Institute, Hyderabad, India

Aruna Meka

Affiliation: Biological Anthropology Unit, Indian Statistical Institute, Hyderabad, India

A. G. Reddy

Affiliation: Centre for Cellular and Molecular Biology, Hyderabad, India

K. Thangaraj

Affiliation: Centre for Cellular and Molecular Biology, Hyderabad, India

Lalji Singh

Affiliation: Centre for Cellular and Molecular Biology, Hyderabad, India

Introduction

Two major routes have been proposed for the initial peopling of East Asia; one via Central Asia to Northeast Asia, which subsequently expanded towards Southeast Asia and beyond, and the other through India to Southeast Asia and further to different regions of East Asia [1]. It is pertinent in this context that the Indian subcontinent has been considered as a major corridor for the migration of human populations to East Asia [2]–[4]. Given its unique geographic position, Northeast India is the only region which currently forms a land bridge between the Indian subcontinent and Southeast Asia, hence hypothesized as an important passage for the initial peopling of East Asia. This region is inhabited by populations belonging to Indo-European, Tibeto-Burman and Austro-Asiatic linguistic families. Whereas Indo-European populations are also found in other parts of India, West Asia and Europe but absent in East Asia, Tibeto-Burman populations are otherwise found only in East Asia. However, Austro-Asiatic speakers, hypothesized as probably the earliest settlers in the Indian subcontinent ([5] and references their in), are also found in other parts of India as well as in East/Southeast Asia. Therefore, if Northeast India had served as an initial corridor, it is likely that the Austro-Asiatic tribes of this region should provide hitherto missing genetic link, which may reflect genetic continuity between Indian and East/Southeast Asian populations. Based on mitochondrial DNA (mtDNA) and Y-chromosome markers, Cordaux et al. [6] observed genetic discontinuity between the Indian and southeast Asian populations and inferred that Northeast India might have acted as a barrier rather than the facilitator of the movement of populations both into and out of India. However, this study included only a few Tibeto-Burman populations of Northeast India whose distribution is restricted only to this region in India [7]–[8], besides a few other populations from other parts of India, possibly with no genetic link with East Asians. It is therefore imperative that the framework of testing such a hypothesis should include adequate representation of these people from Northeast India. Further evidence is needed by way of determining the mtDNA and Y-chromosome haplogroups/lineages of the Austro-Asiatic tribes of the northeastern region and their comparison with appropriate set of South and Southeast Asian populations.

The Northeastern part of India is sandwiched by the marked presence of the young fold mountains of Eastern Himalayas on the northern side and the Indian Ocean on its southern side. The archaeological evidences, based on stone tools from the Garo hills of the Meghalaya region of Northeast India (Figure 1), suggest that this region might have been inhabited as early as in the Paleolithic period [9]–[12]. Within this Meghalaya region, one finds Khasi tribes whose language belongs to the Khasi-Khmic subfamily of the Austro-Asiatics [13], the other two branches of Austro-Asiatics being Mundari and Mon-Khmer, amidst the ethnic majority of Tibeto-Burman populations. Further, the Austro-Asiatic populations of Meghalaya and one of the contiguous Tibeto-Burman tribes namely, Garo, practice matrilineal pattern of descent and matrilocal pattern of residence while the other populations of this region are patrilineal and patrilocal. Such extreme cultural and linguistic diversity of this region may also imply high degree of genetic heterogeneity possibly due to passage of diverse populations through this region. Despite a possible major role played in the population dispersal by the Northeastern region, as transect between India and East Asia, the extent and nature of mtDNA and Y-chromosome diversity of this region is not adequately studied. Therefore, we present results based on the analyses of mtDNA and Y-Chromosome Single Nucleotide Polymorphisms (SNPs) and Short Tandem Repeats (STRs) data of the 8 subgroups of the Austro-Asiatic Khasi from Northeast India (Table 1 and Figure 1), probably for the first time, along with that of the other relevant populations to explore the missing genetic link between the Indian and Southeast Asian populations.

[Figure omitted. See PDF.]

Figure 1. Map of Meghlaya showing the core areas of geographic distribution of the Khasi and Garo populations.

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

[Figure omitted. See PDF.]

Table 1. Areas of sampling and the linguistic affiliations of the Meghalayan populations along with the number of samples typed for Y-chromosome and mtDNA.

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

Results

Distribution and diversity of Y-chromosome

Nei's [14] diversity statistic, h, based on the frequency of different haplogroups (Table 2), ranges from 77% in the Maram to 86.2% in the Pnar among the Khasi-Khmuic Austro-Asiatic groups, whereas it is 77.5% in the Tibeto-Burman Garo. For Y-STR haplotypes, while it ranges from 96.1% in Nongtrai to 99.9% in Khynriam in the Khasi-Khmic populations, it is 99.3% for the Garo. Out of the 26 potential haplogroups defined by the markers used in this study a total of 12 haplogroups were found in these populations (Figure 2). O-M95, with its frequency ranging from 17% in War-Khasi to 42% in War-Jaintia, was the most common haplogroup in all the Austro-Asiatic populations followed by the undifferentiated O-M122 (ranging from 11% in Nongtrai to 34% in Bhoi) where as in the Tibeto-Burman Garo the frequency of O-M134 and undifferentiated O-M122 haplogroups (23% and 17%, respectively) were the most common. H-M69 and its subclade H-M82 which is reported to be in high frequency in most of the Indo-European populations [15] are present with an average frequency of only 3% among them.

[Figure omitted. See PDF.]

Figure 2. Rooted maximum-parsimony tree of Y-chromosome haplogroups defined by binary markers along with their frequency in Nine Meghalayan Populations.

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

[Figure omitted. See PDF.]

Table 2. Genetic Diversity (in percentage) based on Y-Chromosome and mtDNA analysis of populations from Meghalaya

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

Population structure based on Y-chromosome

Based on the multidimensional scaling (MDS) of the Pairwise FST distances computed using haplogroup frequencies of Austro-Asiatic (Khasi from northeast India and others) and neighboring non-Austro-Asiatic populations, the two-dimensional MDS plot is furnished in Figure 3. A good fit between the two-dimensional MDS plot and the source data (pairwise value of FST) was obtained (stress value of 18%). Broadly speaking, most of the Austro-Asiatic populations, including all the three linguistic sub-families of Austro-Asiatics i.e Mundari, Khasi-Khmuic and Mon-Khmer tribes, irrespective of their geographic affiliations, are placed in the upper right quadrant; Nicobarese, Ho, Santhal, She and Zhuang are somewhat removed from the others. On the other hand, most of the Tibeto-Burman populations are differentiated from the Austro-Asiatic populations and the Indo-European populations (clustered in the lower right quadrant) on the 1st and 2nd dimension, respectively. The Khasi-Khmuic populations, which form a compact cluster near the centroid, do not cluster with the Tibeto-Burman populations of Northeast India, barring the Garo of Meghalaya which has contiguous geographic distribution and marital interaction with them. Overall, the populations of the same linguistic family seem to cluster together, with few exceptions such as the Austro-Asiatic Lodha, which is placed among the Indo-European populations.

[Figure omitted. See PDF.]

Figure 3. Plot on the first two dimensions derived from the multidimensional scaling of the pairwise FST distances of the populations based on Y-haplogroups.

Reference to comparative data is given at Table 5. SEA, Southeast Asian; AA, Austro-Asiatic; IE, Indo-European, TB, Tibeto-Burman; Khy, Khynriam; Lyn, Lyngngam; Wkhasi, War Khasi; Viet, Vietnamese;

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

The analysis of molecular variance (AMOVA) yielded a significant but low FST values for both Y-SNPs (0.02) and STRs (0.02), suggesting a probable recent differentiation of the Khasi-Khmuic populations (Table 3). For Y-SNPs, whereas among group differentiation between the Khasi and Southeast Asian Austro-Asiatic populations is low (0.03) and non-significant it is relatively high and significant between the Khasi and Mundari populations (0.08). On the other hand, the FCT value between Khasi-Khmuic and Indian Tibeto-Burman populations is very high and significant (0.30) while between Khasi-Khmuic and Southeast Asian Tibeto-Burman it was relatively low and non-significant (0.03). Although there is virtually no difference in the haplogroupic composition of the Tibeto-Burman Garo from Meghalaya and Southeast Asian Tibeto-Burman populations as suggested by the FCT (−0.01627), it is surprisingly high (0.17975) between the Garo and the other Indian Tibeto-Burman populations.

[Figure omitted. See PDF.]

Table 3. Analysis of Molecular Variance using Y-SNPs/STRs between groups of populations categorized on the basis of geography and languages

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

Profile of new mtDNA haplogroups

Based on Hypervariable segment (HVS) I and the known coding region SNPs most of the individuals could be assigned to specific haplogroups/lineages. However, there were still many individuals who could not be assigned to any existing lineages. Based on their HVS-I motif we could group these samples into 6 broad clades, and resequenced complete mtDNA of 1-2 samples from each of those clades to assign them to a known or new haplogroups (Fig 4). We also resequenced complete mtDNA for the samples falling in haplogroup B as none of the defining mutations for the subhaplogroups of B were found. The analysis of complete mtDNA suggests the presence of four new haplogroups which we have designated as M48, M49, M50 and B7. All the motifs in the coding region of the M48, except for 6336, which defines M30a [16] have not been reported and therefore we assign all these samples a new lineage. While the average frequency of M48 is 11% among the Austro-Asiatic Khasi groups, ranging from zero in War-Jaintia to as high as 26% in Lyngngam, it is present with a frequency of 4% among the Garo. Although haplogroups M49 and M50 are found with an average frequency of about 3% each in the Khasi populations, they could not be traced in the Garo as well as in some of the subgroups of Khasi. A subset of mutations at 150-9452-12950-13928C of our B-haplogroup samples has been reported in one of the samples (SD10313) of Han Chinese [17] which also falls in undifferentiated haplogroup B. We have proposed to name it as haplogroup B7 including the Han Chinese samples.

[Figure omitted. See PDF.]

Figure 4. Phylogenetic tree of new haplogroups based on full mtDNA along with the TMRCA and associated 95% Confidence Interval.

The samples names are on the tip of the haplogroups. The comparative data, for which sample names are written in capital letters, are from Kong et al. [17], Thangaraj et al. [3], Sun et al. [16] and Palanichamy et al. [19].

https://doi.org/10.1371/journal.pone.0001141.g004

In addition to these four new haplogroups we propose two new sub-haplogroups –M33b- within M33, and M31c within M31 (Fig 4). The samples falling in M31c has all the defining mutations of M31 but do not share any of the coding region motifs with either M31a which has been reported in the Andamanese of Andaman and Nicobar island [3] and other Tribal populations of India [18], or M31b found in Rajbanshis (SW1) of Northeast India [19]. Therefore, we propose a new haplogroup, i.e. M31c. While this haplogroup is absent in the Garo, it is found with an average frequency of ∼5% in the Austro-Asiatic Khasi populations with a maximum frequency of ∼17% among the Bhoi. The samples of M33b have mutations which define M33 and it also shares mutations at positions 1719-3221-16293-16324 with the Rajbanshi sample (SW23) which is now re-designated as M33b. The frequency of M33b, with the exception of Pnar (∼22%) is low and found only in Lyngngam, Khynriam and Garo (∼2, ∼3 and ∼3%, respectively). On the other hand, M33a which were found to be in extremely high frequency in the Garo (∼55%) and with an average frequency of∼5% in Khasi-Khumic populations has been also reported in the Brahmins of Uttar Pradesh, India [16] and in the two populations of South India [20]. It is interesting to note that all the samples of this study, except one Khynriam sample, forms a single sublineage defined by 16316 HVS-I motif which distinguishes it from other M33a lineages found in other parts of India.

Distribution and diversity of mtDNA haplotypes/lineages

In the 444 samples representing the 8 Khasi-Khmuic Austro-Asiatic tribes and a Tibeto-Burman Garo a total of 117 distinct HVSI haplotypes were found. Among these, 67 haplotypes are unique, each represented by single individual. Of the remaining, 37 are shared at least by two different tribes out of which only 10 are shared between Garo and Khasi subtribes. Based on the phylogenetic analysis of mtDNA control and coding region SNPs, 37 distinct haplogroups and subhaplogroups were observed among the studied populations (Fig 5). The samples that still remained unclassified in M and R are only ∼6%, and 0.5%, respectively. Among the Austro-Asiatic Khasi, ∼80% of the variation is accounted for by a set of 10 haplogroups–M*, M4a, M9a, M31c, M33a, M33b, M48, MD, MD4 and U2, whereas in the Garo a subset of only 3 haplogroups–M*, M33a and U2-accounted for ∼80% of the total sample. However, these 3 haplogroups account for only ∼18% of the sampled individuals from the neighbouring Austro-Asiatic Khasi populations.

[Figure omitted. See PDF.]

Figure 5. Tree Drawn from a Median-Joining Network of mtDNA Haplogroups Observed in Nine Meghalayan Populations.

The haplogroups box in yellow colour are either new or redefined haplogroups.

https://doi.org/10.1371/journal.pone.0001141.g005

The mtDNA haplogroup diversity (Table 2) among the Austro-Asiatic groups is low and ranges from 83.1% in War-Jaintia to 93.6% in Bhoi whereas in Garo the diversity is extremely low (66.9%). Similarly, the haplotype diversity (Table 2) for the Austro-Asiatic groups ranges from 86.8% in War-Jaintia to 96.1% in Khynriam where as in the Garo it is 68.1%.

Population relationships based on mtDNA haplogroups

The two dimensional plot of the multidimensional scaling of the genetic distance matrix of the 40 populations, including 8 Khasi subtribes and Garo of the present study and other relevant populations from the South and southeast Asia, is shown in Figure 6. The plot depicts the Tibeto-Burman Garo and Austro-Asiatic Nicobarese (a Mon-Khmer population) and Sakai as extreme outliers. As expected, the Mundari Austro-Asiatic populations, with predominantly South Asian mtDNA haplogroups, are placed as outliers aligning with the two Indian Indo-European populations on the extreme right corner of the plot. Although the Khasi-Khmuic Austro-Asiatic populations, except for Nongtrai, Lyngngam, form a constellation near to the left of centroid, it also has other populations such as Han, Lisu and Bai as part of this constellation. The Southeast Asian Tibeto-Burman populations is scattered along the 1st axis. Similarly, the Indian Tibeto-Burmans do not form its own cluster. Overall, the three different sub-families of Austro-Asiatic populations do not form a homogeneous cluster, unlike in the case of Y-chromosome.

[Figure omitted. See PDF.]

Figure 6. Plot on the first two dimensions derived from the multidimensional scaling of the pairwise FST distances of the populations based on mtDNA haplogroups.

Reference to comparative data is given at Table 6. SEA, Southeast Asian; AA, Austro-Asiatic; IE, Indo-European; TB, Tibeto-Burman; Wk, WarKhasi; Ch, Chong;

https://doi.org/10.1371/journal.pone.0001141.g006

Although the AMOVA suggests low FST value (0.05), hence low differentiation among the Khasi-Khmuic populations (Table 4), it is quite high between them and the Garo (0.12). The differentiation of Khasi-Khmuic tribes with Southeast Asian Austro-Asiatic populations is moderate (0.05) but is much higher with those of Mundari populations (0.12). Surprisingly, the Tibeto-Burman Garo of Meghalaya shows high degree of differentiation with the other Tibeto-Burman populations of India (0.17) as well as Southeast Asia (0.13).

[Figure omitted. See PDF.]

Table 4. Analysis of Molecular Variance using mtDNA haplogroups between groups of populations categorized on the basis of geography and languages

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

Time to Most Recent Common Ancestors (TMRCA)

The TMRCA was calculated based on mtDNA coding region (nucleotide position 577-16023) with the average sequence evolution rate as 1.26±0.08×10−8 base substitutions per nucleotide per year [21]. The TMRCA of the haplogroups based on the full mtDNA sequence (Fig 4) suggest a younger age of Khasi/Northeast Indian haplogroup M (41,000 YBP) compared to what has been obtained in the other studies [16] for Indian M haplogroup (54,000 YBP). This is because of the very low age contribution from the M48 haplogroup. Reanalyzing the data by removing M48 increases the age to ∼50,000 YBP which is close to what has been obtained in the other studies. The TMRCA of haplogroup M31 and M33 is ∼40,000 YBP and ∼50,000 YBP suggesting that M33 like M31 is an archaic lineage. The age of B7 suggests that this haplogroup has originated ∼28,000 YBP in East Asia where all the other sub-haplogroups of B have been hypothesized to have originated.

Discussion

Origin and spread of predominant Y-haplogroups of Northeast India

The present study provides a comprehensive genetic analysis of the populations of Austro-Asiatic linguistic family inhabiting the Northeast Indian region, which has probably served as the corridor for the peopling of Southeast Asia. Two major haplogroups O-M95 and O-M122 and its subhaplogroups dominate the haplogroup composition of the Austro-Asiatic and Tibeto-Burman populations of Northeast India and East and Southeast Asia (Table 5). Indo-European groups of Northeast India lack these haplogroups. Kumar et al. [22] suggest that the haplogroup O-M95 had its origin probably in the ancestors of Mundari, one of the major subfamilies of the Indian Austro-Asiatics inhabiting Eastern and Central regions of the country, and correlated its spread with the movement of these populations to Southeast Asia and to other parts inhabited by the Austro-Asiatics via the Northeast Indian corridor. Our results are consistent with this observation as we not only find relatively high frequency O-M95 (Table 5) in all the subgroups of Khasi from Northeast India, but also find a decreasing gradient of O-M95 average frequency from Mundari (53%) to Khasi (∼31%) to Southeast Asian Austro-Asiatics (∼23%), suggesting diffusion of O-M95 from India to Southeast Asia. Further, O-M95 is either absent or, if present, only with negligible frequency in the other linguistic groups of India. It is found to be either absent or has a negligible presence in the other Tibeto-Burman populations of Northeast India and the moderate frequency of O-M95 in the Tibeto-Burman Garo may be because of the gene flow from the neighboring Khasi which is facilitated by the practice of matrilocality in these populations. This suggests that the haplogroup O-M95 is, by and large, restricted to Austro-Asiatic populations in India and supports the contention of Kumar et al. [22] that probably there was a concomitant spread of the ancestors of the present day AA people with proto-Austro-Asiatic language and haplogroup O-M95.

[Figure omitted. See PDF.]

Table 5. Y-haplogroup frequencies in different linguistic populations of South and Southeast Asia.

https://doi.org/10.1371/journal.pone.0001141.t005

Haplogroup O-M122 is found to be in high frequency in the Garo as well as in Khasi-Khmuic populations. However, further typing of O-M122 chromosomes suggests a high frequency of undifferentiated O-M122 among the Khasi-Khmic populations, whereas the frequency of O-M134 (one of the subhaplogroups of O-M122) is found to be much higher among the Garo. Incidentally, O-M134 is found to be in much higher frequency compared to the undifferentiated O-M122 in the other Indian Tibeto-Burman populations as well. Further, Tibeto-Burman populations of Southeast Asia also have relatively much higher frequency O-M134 compared to the Austro-Asiatics there whose undifferentiated O-M122 samples fall mostly into subhaplogroup O-M159 (Table 5). The presence of O-M134 in high frequency among the Tibeto-Burman populations, both from India and East/southeast Asia, strongly suggests possibility of its correlation with the migration and spread of Tibeto-Burman populations into India.

Distinct origin of Khasi-Khmuic tribes

Khasi-Khumic groups are surrounded by the predominant but ethnically similar Tibeto-Burman populations in the region and a comparison with them, particularly with the Garo which is also matrilocal like the Khasi and with whom they have contiguous distribution in Meghalaya, suggests a distinct genetic origin of the Khasi-Khmuic populations. For example, unlike the Khasi, the Tibeto-Burman Garo shows a very low frequency of O-M95 and undifferentiated O-M122 Y-haplogroups but a high frequency of differentiated O-M122 i.e. O-M134 and O-M133. Similarly, the other Tibeto-Burman groups from India show a low frequency of O-M95 and undifferentiated O-M122 but a very high frequency of O-M134 (Table 5). Further, while the Garo shows a high frequency of M33a (55%), it lacks M9a, MD and MD4 mtDNA haplogroups which are in high frequency among the Khasi. The other Tibeto-Burman groups from this region, although has 30% of unclassified M*, show a high presence F and its subhaplogroup, M8c, A and its subhaplogroups most of which have a negligible presence in the Austro-Asiatic Khasi populations. Moreover, M9a, MD and MD4 mtDNA haplogroups which account for ∼30% of Khasi-Khmuic samples are found only with a frequency of ∼10% among the Tibeto-Burman groups. Most importantly, Khasi-khumic group has ∼25% new haplogroups (M48, M49, M50, M31c and M33b) which has not been reported so far from any of the Northeast Indian groups except Garo which has M48 and M33b with a combined frequency of ∼7%. Thus, the composition of both the mtDNA and Y haplogroups in the Austro-Asiatic Khasi as a whole suggests their distinct origin and a separate migration vis-à-vis the Tibeto-Burman groups of this region.

mtDNA landscape of the Meghalaya populations

The Khasi-Khmuic and the Garo populations of Meghalaya essentially have three kinds of haplogroups: the commonly found South Asian haplogroups, East Asian haplogroups and new haplogroups (Table 6 and Fig 5). Among the new haplogroups, B7 is found only in Khasi-Khmuic populations and it has been otherwise reported as unclassified B in Han Chinese [17] suggesting its probable origin in East Asia, as is the case with other haplogroups of B. However, all the other new haplogroups viz. M48, M49 and M50 or the redefined subhaplogroups i.e. M31c and M33b have not yet been reported from East Asia and the neighboring East Asian populations and these haplogroups needs to be assessed, especially from Myanmar region, to trace their origin and movement. Further, these haplogroups have not been reported among the 23 Indo-European and 25 Dravidian populations of the Indian subcontinent [23]–[25] confirming the absence of this haplogroups in other parts of India.

[Figure omitted. See PDF.]

Table 6. mtDNA haplogroup frequencies in different linguistic populations of South and Southeast Asia.

https://doi.org/10.1371/journal.pone.0001141.t006

The phylogeny of M33 and M31 is quite intriguing. It is striking that these two lineages have evolved into ethnic specific branches, separated by a number of mutations, suggesting their deep antiquity. Haplogroup M31 has three subhpalogroups-M31a, M31b and M31c (Fig 4). While M31a is reported with a high frequency in the Greater Andmanese [3] and Lodha, Chenchu and Lambadi tribal groups of India [18], M31b and M31c are found in Northeast India, M31b among the Rajbanshis [19] and M31c with an average frequency of ∼5% in the Khasi-Khmuic populations (Fig 5). However, this haplogroup is absent in the Tibeto-Burmans of Northeast India or among the other Indian populations. Haplogroup M33b has a total frequency of ∼4% among the Khasi and otherwise it has been reported only from Rajbanshi from this region. On the other hand, haplogroup M33a, which has been reported in the Brahmins of Uttar Pradesh, India [16] and in the two populations of South India [20] bifurcates into two branches-one, without the motif 16316, identified in only 1 Khynriam sample and another defined by 16316 HVS-I motif and found in all the samples of M33a of Garo (∼55%) and Khasi-Khumic populations of Meghalaya (∼5%). However, none of the other Indian Austro-Asiatic populations (Mundari) shows M33 or its subhaplogroup (Kumar et al., unpublished results). Although we need to screen more populations in order to resolve the origin and distribution of M31 and M33, both these lineages are very old (>40,000 YBP) and probably originated in mainland India as their presence has been reported only from the Indian subcontinent. Overall, the new and the redefined haplogroups, excluding M33a, account for ∼25% of the total frequency in the Khasi-Khmuic populations and ∼4% in the Garo, while M33a accounts for ∼55% in Garo. The rest of the haplogroups are of either South Asian or East Asian types.

The East Asian haplogroups with high frequency in the Khasi-Khmuic populations (M9a, MD and MD4) account for ∼35% of the samples. Along with B7, the presence of very different set of East Asian haplogroups among them, as compared to the Tibeto-Burmans of India, may suggest a strong possibility of their admixture with or assimilation of certain East Asian populations, other than the Indian Tibeto-Burmans. This is also reflected in case of Y-chromosome with high frequency of undifferentiated O-M122. The South Asian haplogroups found in Khasi-Khmuic populations are M2 and M2a, M18, M25, M3, M4a, M5, R, R8, U2 and U2abc, and U3 which account for a total of ∼22% among them. However, the frequency of all these haplogroups is quite low except for M4a and U2, which account for ∼7 and ∼8%, respectively. The low frequency of South Asian haplogroups might be either due to admixture with the surrounding Indo-European populations or due to initial splintering of their gene pool from that of the Mundari tribes as both these linguistic groups show predominantly South Asian mtDNA haplogroups.

The Austro-Asiatic tribes of Northeast India: Genetic link/continuity between South and Southeast Asian populations

A rapid human migration through Southern route ∼60,000 YBP is suggested to have brought undifferentiated M and R into South and East Asia and subsequently differentiated into different subhaplogroups in different regions [17]. For example, the undifferentiated M evolved into different haplogroups such as M2, M3, M5 M6, M18 etc. in South Asia, whereas in East Asia, primarily haplogroup M9, MD, MG etc. evolved. In this backdrop, it is interesting that the Mundari Austro-Asiatic tribes of central and eastern India have mostly South Asian haplogroups (Kumar et al. unpublished results) with a high frequency of M2b, M40a, R6 and R7, which are considered to be the old lineages. This may support the migration of these tribes into India by Southern route and using the Western Indian corridor, bringing with them the undifferentiated M and R ∼50,000 YBP. The predominant lineages found among them might have evolved subsequently. However, these mtDNA haplogroups are with very low frequencies in the Khasi-Khmuic Austro-Asiatic populations from Northeast India who, nonetheless, have ∼25% of their haplogroups as new ones (Table 6). This would be feasible only if a section of them had separated soon after the ancestors of Austro-Asiatics had come to India, and migrated quite rapidly to Northeast India and further to Southeast Asia before the differentiation of M and R into any of the haplogroups found in South Asia. The undifferentiated M and R haplogroups carried by this group might have evolved into certain characteristic haplogroups such as M48, M49 and M50 in the Khasi but not in the Mundari groups around 40,000 YBP. A section of these initial migrants to Northeast India had probably moved further to Southeast Asia in rapid succession carrying with them the undifferentiated M (before the Khasi- specific haplogroups could evolve), which might have evolved into typical East Asian haplogroups. It is evident from the foregoing discussion that the Austro-Asiatic Khasi of Northeast India represents genetic continuity, linking the populations of South and Southeast Asia. Therefore, our findings reinforce the suggestion that Northeast India has acted as a corridor for initial movement of populations, not as a barrier as suggested in a recent study [6].

Materials and Methods

Collection of blood samples

About 5 ml of intravenous blood samples were collected in 5 ml Tarson tubes containing EDTA as an anticoagulant from a total of 444 healthy unrelated volunteers, both males and females, after obtaining the informed written consent. Prior approval for the study was obtained from the ‘Indian Statistical Institute Review Committee for Protection of Research Risk to Humans’. These samples represent 8 subgroups of the Austro-Asiatic Khasi, besides the neighbouring Tibeto-Burman Garo tribes from Meghalaya in the Northeastern part of India (Table 1 and Figure 1). While the samples for the different dialectical groups of Khasi tribe were drawn from almost all the areas of their distribution in Khasi hills, Garo samples were represented in bulk from the contiguous areas of the Lyngngam tribe, such as Rongjeng, Khonjoy, Shallang, etc. with a few samples drawn from all over the capital city of Meghalaya, Shillong.

mtDNA typing

The mtDNA genomes were amplified and sequenced by means of the procedures described in a recent study [17]. Sequences were edited and mutations scored relative to the revised Cambridge Reference Sequence [26] (rCRS). Initially, HVS-I (nucleotide positions [nps] 160001–16400) was sequenced, besides typing SNPs at 10398 and 10400. The published HVS-I sequences [16], [24], [27]–[30] were compared with our HVS-I sequences from Meghalaya in combination with the typed SNPs of this study to identify the relevant coding region SNPs, which are diagnostic of the main haplogroups and subhaplogroups within the mtDNA phylogeny. These coding regions were then selectively assayed by sequencing to obtain haplogroups and their derivatives (Fig 5). Representative samples from those that could not be assigned into known haplogroups were carefully chosen and complete mtDNA sequencing was done to assign/designate them to the new (sub)haplogroups (Fig 4).

Y-chromosome typing

The following 25 Y-SNPs which are known to detect variations in Asia were screened: RPS4Y, M210, M38, M8, M217, M93, M86, M89, M69, M82, M172, M9, M175, M122, M7, M164, M159, M121, M134, M133, M162, M95, M11, M45 and M173 [31]–[33]. The rooted maximum-parsimony trees of the haplogroups defined by these markers are presented in Figure 2. Many of the samples were typed with all the binary markers for internal check on the reliability of the typing and also to detect recurrent mutations. The nomenclature as suggested and followed by Y-chromosome consortium [32] was used. The following 6 Y-STRs loci were also typed: DYS19, DYS389I, DYS389b, DYS390, DYS391 and DYS393 (for data refer to Dataset S1) and the details of these loci are given at Butler et al [34].

Statistical Analysis

Since the DYS389II allele length also contains DYS389I, for all statistical analyses a simple subtraction of DYS389I allele length from that of DYS389II was done to avoid the double-counting variation at DYS389I. The subtracted DYS389II allele is named as DYS389b. The Y-SNP and modified Y-STR data were then analyzed for haplogroup and haplotype diversity, respectively, along with their associated Standard Error by means of the software package ARLEQUIN 3.01 [35]. The frequencies of haplogroups constructed by binaray markers were used to compute pairwise FST genetic distance matrix. Based on the distance matrix, MDS analysis was performed using SPSS package. The genetic structure as reflected in the distribution of Y-SNPs and STRs was further explored through AMOVA, by grouping populations based on their geography and linguistic affiliations

Haplotype and Haplogroup diversity of mtDNA and AMOVA based on haplogroup frequencies were calculated using the ARLEQUIN 3.01 [35]. Phylogenetic relationships between the observed haplogroups were first drawn by hand and then confirmed by using the NETWORK program [36]. The frequencies of haplogroups were used as input vector to compute pairwise FST genetic distance matrix using Arlequin 3.01. Based on this distance matrix, MDS analysis was performed using SPSS package. The TMRCA of the clades and subclades and their associated SEs were calculated using their full mtDNA sequences [21].

Supporting Information

[Figure omitted. See PDF.]

Dataset S1.

Six Y-STR data of all the 9 populations

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

(0.05 MB DOC)

Acknowledgments

This study is essentially a part of the Indian Statistical Institute project undertaken by BMR in collaboration with Centre for Cellular & Molecular Biology, Hyderabad. BMR thanks Directors of both the Institutes for logistic support.

Author Contributions

Conceived and designed the experiments: BR BL. Performed the experiments: AR BL TN AR. Analyzed the data: VK AM. Contributed reagents/materials/analysis tools: KT BR LS. Wrote the paper: BR VK. Other: Collected blood samples: BL. Commented on the manuscript and participated in the discussions: KT LS.

Citation: Reddy BM, Langstieh BT, Kumar V, Nagaraja T, Reddy ANS, Meka A, et al. (2007) Austro-Asiatic Tribes of Northeast India Provide Hitherto Missing Genetic Link between South and Southeast Asia. PLoS ONE2(11): e1141. https://doi.org/10.1371/journal.pone.0001141

References

1. Cavalli-Sforza LL, Menozzi P, Piazza A (1994) The history and geography of human genes. Princeton, NJ: University Press. LL Cavalli-SforzaP. MenozziA. Piazza1994The history and geography of human genes.Princeton, NJUniversity Press

2. Jin L, Su B (2000) Natives or immigrants: modern human origin in East Asia. Nat Rev Genet 1: 126–133.L. JinB. Su2000Natives or immigrants: modern human origin in East Asia.Nat Rev Genet1126133

3. Thangaraj K, Chaubey G, Kivisild T, Reddy AG, Singh VK, et al. (2005) Reconstructing the origin of Andaman Islanders. Science 308: 996.K. ThangarajG. ChaubeyT. KivisildAG ReddyVK Singh2005Reconstructing the origin of Andaman Islanders.Science308996

4. Macaulay V, Hill C, Achilli A, Rengo C, Clarke D, et al. (2005) Single, rapid coastal settlement of Asia revealed by analysis of complete mitochondrial genomes. Science 308: 1034–1036.V. MacaulayC. HillA. AchilliC. RengoD. Clarke2005Single, rapid coastal settlement of Asia revealed by analysis of complete mitochondrial genomes.Science30810341036

5. Kumar V, Reddy BM (2003) Status of Austro-Asiatic groups in the peopling of India: An exploratory study based on the available prehistoric, Linguistic and Bioligcal evidences. J Biosci 28: 507–522.V. KumarBM Reddy2003Status of Austro-Asiatic groups in the peopling of India: An exploratory study based on the available prehistoric, Linguistic and Bioligcal evidences.J Biosci28507522

6. Cordaux R, Weiss G, Saha N, Stoneking M (2004) The northeast Indian passageway: a barrier or corridor for human migrations? Mol Biol Evol 21: 1525–1533.R. CordauxG. WeissN. SahaM. Stoneking2004The northeast Indian passageway: a barrier or corridor for human migrations?Mol Biol Evol2115251533

7. Su B, Xiao C, Deka R, Seielstad MT, Kangwanpong D, et al. (2000) Y chromosome haplotypes reveal prehistorical migrations to the Himalayas. Hum. Genet 107: 582–590.B. SuC. XiaoR. DekaMT SeielstadD. Kangwanpong2000Y chromosome haplotypes reveal prehistorical migrations to the Himalayas.Hum. Genet107582590

8. Krithika S, Trivedi R, Kashyap VK, Bharati P, Vasulu TS (2006) Antiquity, geographic contiguity and genetic affinity among Tibeto-Burman populations of India: a microsatellite study Ann Hum Biol 33: 26–42.S. KrithikaR. TrivediVK KashyapP. BharatiTS Vasulu2006Antiquity, geographic contiguity and genetic affinity among Tibeto-Burman populations of India: a microsatellite studyAnn Hum Biol332642

9. Sharma TC (1966) Researches on the prehistoric archaeology of Assam. Jr Assam Sc Soc IX: 1–11.TC Sharma1966Researches on the prehistoric archaeology of Assam.Jr Assam Sc SocIX111

10. Sharma TC (1980) Prehistoric archaeology in North Eastern India-A review in progress. In: Sharma TC, Majumdar DN, editors. Eastern Himalayas. New Delhi: Cosmo Publications. pp. 102–107.TC Sharma1980Prehistoric archaeology in North Eastern India-A review in progress.TC SharmaDN MajumdarEastern HimalayasNew DelhiCosmo Publications102107

11. Sharma TC (1984) North East India in prehistoric times. In: Karotemprel S, editor. Tribes of North East India. Shillong: Vendrame Missiological Institute, Centre of Indigenous Culture. pp. 1–23.TC Sharma1984North East India in prehistoric times.S. KarotemprelTribes of North East IndiaShillongVendrame Missiological Institute, Centre of Indigenous Culture123

12. Hussain Z (1991) Who are the pre-historic dwellers of the Meghalaya plateau? In: Singh JP, Sengupta G, editors. Archaeology of North East India. New Delhi: Vikas Publishing House Pvt. Ltd. pp. 74–85.Z. Hussain1991Who are the pre-historic dwellers of the Meghalaya plateau?JP SinghG. SenguptaArchaeology of North East IndiaNew DelhiVikas Publishing House Pvt. Ltd7485

13. Diffloth G (2005) The contribution of linguistic palaeontology to the homeland of Austro-Asiatic. In: Sagart L, Blench R, Sanchez-Mazas A, editors. The Peopling of East Asia: Putting Together Archaeology, Linguistics and Genetics. London: Routledge Curzon. pp. 77–81.G. Diffloth2005The contribution of linguistic palaeontology to the homeland of Austro-Asiatic.L. SagartR. BlenchA. Sanchez-MazasThe Peopling of East Asia: Putting Together Archaeology, Linguistics and GeneticsLondonRoutledge Curzon7781

14. Nei M (1987) Molecular evolutionary genetics. New York: Columbia University Press. M. Nei1987Molecular evolutionary genetics.New YorkColumbia University Press

15. 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–21.S. SenguptaLA ZhivotovskyR. KingSQ MehdiCA Edmonds2006Polarity 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 Genet7820221

16. Sun C, Kong QP, Palanichamy MG, Agrawal S, Bandelt HJ, et al. (2006) The dazzling array of basal branches in the mtDNA macrohaplogroup M from India as inferred from complete genomes. Mol Biol Evol 23: 683–690.C. SunQP KongMG PalanichamyS. AgrawalHJ Bandelt2006The dazzling array of basal branches in the mtDNA macrohaplogroup M from India as inferred from complete genomes.Mol Biol Evol23683690

17. Kong QP, Yao YG, Sun C, Bandelt HJ, Zhu CL, et al. (2003) Phylogeny of East Asian mitochondrial DNA lineages inferred from complete sequences. Am J Hum Genet 73: 671–676.QP KongYG YaoC. SunHJ BandeltCL Zhu2003Phylogeny of East Asian mitochondrial DNA lineages inferred from complete sequences.Am J Hum Genet73671676

18. Endicott P, Metspalu M, Stringer C, Macaulay V, Cooper A, et al. (2006) Multiplexed SNP typing of ancient DNA clarifies the origin of andaman mtDNA haplogroups amongst south Asian tribal populations. PLoS ONE 1: e81.P. EndicottM. MetspaluC. StringerV. MacaulayA. Cooper2006Multiplexed SNP typing of ancient DNA clarifies the origin of andaman mtDNA haplogroups amongst south Asian tribal populations.PLoS ONE1e81

19. Palanichamy MG, Agrawal S, Yao YG, Kong QP, Sun C, et al. (2006) Comment on “Reconstructing the origin of Andaman islanders”. Science 311: 470.MG PalanichamyS. AgrawalYG YaoQP KongC. Sun2006Comment on “Reconstructing the origin of Andaman islanders”.Science311470

20. Thanseem I, Thangaraj K, Chaubey G, Singh VK, Bhaskar LV, et al. (2006) Genetic affinities among the lower castes and tribal groups of India: inference from Y chromosome and mitochondrial DNA. BMC Genet 7: 42.I. ThanseemK. ThangarajG. ChaubeyVK SinghLV Bhaskar2006Genetic affinities among the lower castes and tribal groups of India: inference from Y chromosome and mitochondrial DNA.BMC Genet742

21. Mishmar D, Ruiz-Pesini E, Golik P, Macaulay V, Clark AG, et al. (2003) Natural selection shaped regional mtDNA variation in humans. Proc Natl Acad Sci U S A 7: 171–6.D. MishmarE. Ruiz-PesiniP. GolikV. MacaulayAG Clark2003Natural selection shaped regional mtDNA variation in humans.Proc Natl Acad Sci U S A71716

22. Kumar V, Reddy AN, 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.V. KumarAN ReddyJP BabuTN RaoBT Langstieh2007Y-chromosome evidence suggests a common paternal heritage of Austro-Asiatic populations.BMC Evol Biol.747

23. Roychoudhury S, Roy S, Basu A, Banerjee R, Vishwanathan H, et al. (2001) Genomic structures and population histories of linguistically distinct tribal groups of India. Hum Genet 3: 339–350.S. RoychoudhuryS. RoyA. BasuR. BanerjeeH. Vishwanathan2001Genomic structures and population histories of linguistically distinct tribal groups of India.Hum Genet3339350

24. Kivisild T, Rootsi S, Metspalu M, Mastana S, Kaldma K, et al. (2003) The genetic heritage of the earliest settlers persists both in Indian. tribal and caste populations. Am J Hum Genet 72: 313–332.T. KivisildS. RootsiM. MetspaluS. MastanaK. Kaldma2003The genetic heritage of the earliest settlers persists both in Indian. tribal and caste populations.Am J Hum Genet72313332

25. Metspalu M, Kivisild T, Metspalu E, Parik J, Hudjashov G, et al. (2004) Most of the extant mtDNA boundaries in south and southwest Asia were likely shaped during the initial settlement of Eurasia by anatomically modern humans. BMC Genet 5: e26.M. MetspaluT. KivisildE. MetspaluJ. ParikG. Hudjashov2004Most of the extant mtDNA boundaries in south and southwest Asia were likely shaped during the initial settlement of Eurasia by anatomically modern humans.BMC Genet5e26

26. Andrews RM, Kubacka I, Chinnery PF, Lightowlers RN, Turnbull DM, et al. (1999) Reanalysis and revision of the Cambridge reference sequence for human mitochondrial DNA. Nat Genet 23: 147.RM AndrewsI. KubackaPF ChinneryRN LightowlersDM Turnbull1999Reanalysis and revision of the Cambridge reference sequence for human mitochondrial DNA.Nat Genet23147

27. Macaulay V, Richards M, Hickey E, Vega E, Cruciani F, et al. (1999) The emerging tree of West Eurasian mtDNAs: a synthesis of control-region sequences and RFLPs. Am J Hum Genet 64: 232–249.V. MacaulayM. RichardsE. HickeyE. VegaF. Cruciani1999The emerging tree of West Eurasian mtDNAs: a synthesis of control-region sequences and RFLPs.Am J Hum Genet64232249

28. Qunitana-Murci L, Semino O, Bandelt HJ, Passarino G, McElreavey K, et al. (1999) Genetic evidence of an early exit of Homo sapies sapiens from Africa through eastern Africa. Nat Genet 23: 437–441.L. Qunitana-MurciO. SeminoHJ BandeltG. PassarinoK. McElreavey1999Genetic evidence of an early exit of Homo sapies sapiens from Africa through eastern Africa.Nat Genet23437441

29. Richards M, Macaulay V, Hickey E, Vega E, Sykes B, et al. (2000) Tracing European founder lineages in the Near Eastern mtDNA pool. Am J Hum Genet 67: 1251–1276.M. RichardsV. MacaulayE. HickeyE. VegaB. Sykes2000Tracing European founder lineages in the Near Eastern mtDNA pool.Am J Hum Genet6712511276

30. Palanichamy MG, Sun C, Agrawal S, Bandelt HJ, Kong QP, et al. (2004) Phylogeny of mitochondrial DNA macrohaplogroup N in India, based on complete sequencing: implications for the peopling of South Asia. Am J Hum Genet 75: 966–978.MG PalanichamyC. SunS. AgrawalHJ BandeltQP Kong2004Phylogeny of mitochondrial DNA macrohaplogroup N in India, based on complete sequencing: implications for the peopling of South Asia.Am J Hum Genet75966978

31. Underhill PA, Passarino G, Lin AA, Shen P, Mirazon Lahr M, et al. (2001) The phylogeography of Y chromosome binary haplotypes and the origins of modern human populations. Ann Hum Genet 65: 43–62.PA UnderhillG. PassarinoAA LinP. ShenM. Mirazon Lahr2001The phylogeography of Y chromosome binary haplotypes and the origins of modern human populations.Ann Hum Genet654362

32. YCC (The Y Chromosome Consortium) (2002) A nomenclature system for the tree of human Y-chromosomal binary haplogroups. Genome Res 12: 339–348.YCC (The Y Chromosome Consortium)2002A nomenclature system for the tree of human Y-chromosomal binary haplogroups.Genome Res12339348

33. Jobling MA, Tyler-Smith C (2003) The human Y chromosome: an evolutionary marker comes of age. Nat Rev Genet 8: 598–612.MA JoblingC. Tyler-Smith2003The human Y chromosome: an evolutionary marker comes of age.Nat Rev Genet8598612

34. Butler JM, Schoske R, Vallone PM, Kline MC, Redd AJ, et al. (2002) A novel multiplex for simultaneous amplification of 20 Y chromosome STR markers. Forensic Sci Int 129: 10–24.JM ButlerR. SchoskePM ValloneMC KlineAJ Redd2002A novel multiplex for simultaneous amplification of 20 Y chromosome STR markers.Forensic Sci Int1291024

35. Excoffier L, Laval G, Schneider S (2005) Arlequin ver. 3.0: An integrated software package for population genetics data analysis. Evolutionary Bioinformatics Online 1: 47–50.L. ExcoffierG. LavalS. Schneider2005Arlequin ver. 3.0: An integrated software package for population genetics data analysis.Evolutionary Bioinformatics Online14750

36. Bandelt H, Forster P, Rohl A (1999) Median joining networks for inferring intraspecific phylogenies. Mol Biol Evol 16: 37–48.H. BandeltP. ForsterA. Rohl1999Median joining networks for inferring intraspecific phylogenies.Mol Biol Evol163748

37. Thangaraj K, Singh L, Reddy A, Rao V, Sehgal S, et al. (2003) Genetic affinities of the andaman islanders, a vanishing human population. Curr Biol 13: 86–93.K. ThangarajL. SinghA. ReddyV. RaoS. Sehgal2003Genetic affinities of the andaman islanders, a vanishing human population.Curr Biol138693

38. Karafet T, Xu L, Du R, Wang W, Feng S, et al. (2001) Paternal population history of East Asia: sources, patterns, and microevolutionary processes. Am J Hum Genet 69: 615–628.T. KarafetL. XuR. DuW. WangS. Feng2001Paternal population history of East Asia: sources, patterns, and microevolutionary processes.Am J Hum Genet69615628

39. Fucharoen G, Fucharoen S, Horai S (2001) Mitochondrial DNA polymorphisms in Thailand. J Hum Genet 3: 115–125.G. FucharoenS. FucharoenS. Horai2001Mitochondrial DNA polymorphisms in Thailand.J Hum Genet3115125

40. Oota H, Settheetham-Ishida W, Tiwawech D, Ishida T, Stoneking M (2001) Human mtDNA and Y-chromosome variation is correlated with matrilocal versus patrilocal residence. Nat Genet 29: 20–1.H. OotaW. Settheetham-IshidaD. TiwawechT. IshidaM. Stoneking2001Human mtDNA and Y-chromosome variation is correlated with matrilocal versus patrilocal residence.Nat Genet29201

41. Prasad BV, Ricker CE, Watkins WS, Dixon ME, Rao BB, et al. (2001) Mitochondrial DNA variation in Nicobarese Islanders. Hum Biol 5: 715–725.BV PrasadCE RickerWS WatkinsME DixonBB Rao2001Mitochondrial DNA variation in Nicobarese Islanders.Hum Biol5715725

42. Yao Y, Nie L, Harpending H, Fu Y, Yuan Z, et al. (2002) Genetic relationship of Chinese ethnic populations revealed by mtDNA sequence diversity. Am J Phys Anthropol 118: 63–76.Y. YaoL. NieH. HarpendingY. FuZ. Yuan2002Genetic relationship of Chinese ethnic populations revealed by mtDNA sequence diversity.Am J Phys Anthropol1186376

43. Yao YG, Kong QP, Bandelt HJ, Kivisild T, Zhang YP (2002) Phylogeographic differentiation of mitochondrial DNA in Han Chinese. Am J Hum Genet 70: 635–651.YG YaoQP KongHJ BandeltT. KivisildYP Zhang2002Phylogeographic differentiation of mitochondrial DNA in Han Chinese.Am J Hum Genet70635651