Multiple origins of the mtDNA 9-bp deletion in populations of South India

The origins and genetic affinities of the more than 500 tribal populations living in South Asia are widely disputed. This may reflect differential contributions that continental populations have made to tribal groups in South Asia. We assayed for the presence of the intergenic COII/tRNALys 9-bp deletion in human mtDNA in 646 individuals from 12 caste and 14 tribal populations of South India and compared them to individuals from Africa, Europe, and Asia. The 9-bp deletion is observed in four South Indian tribal populations, the Irula, Yanadi, Siddi, and Maria Gond, and in the Nicobarese. Length polymorphisms of the 9-bp motif are present in the Santal, Khonda Dora, and Jalari, all of whom live in a circumscribed region on the eastern Indian coast. Phylogenetic analyses of mtDNA control region sequence from individuals with the 9-bp deletion indicate that it has arisen independently in some Indian tribal populations. Other 9-bp deletion haplotypes are likely to be of Asian and African origin, implying multiple origins of the 9-bp deletion in South India. These results demonstrate varying genetic affinities of different South Indian tribes to continental populations and underscore the complex histories of the tribal populations living in South Asia. Am J Phys Anthropol 109:147–158, 1999. © 1999 Wiley-Liss, Inc.     

Different population histories of the Mundari- and Mon-Khmer-speaking Austro-Asiatic tribes inferred from the mtDNA 9-bp deletion/insertion polymorphism in Indian populations

Length variation in the human mtDNA intergenic region between the cytochrome oxidase II (COII) and tRNA lysine (tRNA^lys) genes has been widely studied in world populations. Specifically, Austronesian populations of the Pacific and Austro-Asiatic populations of southeast Asia most frequently carry the 9-bp deletion in that region implying their shared common ancestry in haplogroup B. Furthermore, multiple independent origins of the 9-bp deletion at the background of other mtDNA haplogroups has been shown in populations of Africa, Europe, Australia, and India. We have analyzed 3293 Indian individuals belonging to 58 populations, representing different caste, tribal, and religious groups, for the length variation in the 9-bp motif. The 9-bp deletion (one copy) and insertion (three copies) alleles were observed in 2.51% (2.15% deletion and 0.36% insertion) of the individuals. The maximum frequency of the deletion (45.8%) was observed in the Nicobarese in association with the haplogroup B5a D-loop motif that is common throughout southeast Asia. The low polymorphism in the D-loop sequence of the Nicobarese B5a samples suggests their recent origin and a founder effect, probably involving migration from southeast Asia. Interestingly, none of the 302 (except one Munda sample, which has 9-bp insertion) from Mundari-speaking Austro-Asiatic populations from the Indian mainland showed the length polymorphism of the 9-bp motif, pointing either to their independent origin from the Mon-Khmeric-speaking Nicobarese or to an extensive admixture with neighboring Indo-European-speaking populations. Consistent with previous reports, the Indo-European and Dravidic populations of India showed low frequency of the 9-bp deletion/insertion. More than 18 independent origins of the deletion or insertion mutation could be inferred in the phylogenetic analysis of the D-loop sequences.     

Polynesian genetic affinities with Southeast Asian populations as identified by mtDNA analysis.

Polynesian genetic affinities to populations of Asia were studied using mtDNA markers. A total of 1,037 individuals from 12 populations were screened for a 9-bp deletion in the intergenic region between the COII and tRNA(Lys) genes that approaches fixation in Polynesians. Sequence-specific oligonucleotide probes that identify specific mtDNA control region nucleotide substitutions were used to describe variation in individuals with the 9-bp deletion. The 9-bp deletion was not observed in northern Indians, Bangladeshis, or Pakistanis but was seen at low to moderate frequencies in the nine other Southeast Asian populations. Three substitutions in the control region at positions 16217, 16247, and 16261 have previously been observed at high frequency in Polynesian mtDNAs; this "Polynesian motif" was observed in 20% of east Indonesians with the 9-bp deletion but was observed in only one additional individual. mtDNA types related to the Polynesian motif are highest in frequency in the corridor from Taiwan south through the Philippines and east Indonesia, and the highest diversity for these types is in Taiwan. These results are consistent with linguistic evidence of a Taiwanese origin for the proto-Polynesian expansion, which spread throughout Oceania by way of Indonesia.     

Maternal footprints of Southeast Asians in North India

We have analyzed 7,137 samples from 125 different caste, tribal and religious groups of India and 99 samples from three populations of Nepal for the length variation in the COII/tRNA(Lys) region of mtDNA. Samples showing length variation were subjected to detailed phylogenetic analysis based on HVS-I and informative coding region sequence variation. The overall frequencies of the 9-bp deletion and insertion variants in South Asia were 1.9 and 0.6%, respectively. We have also defined a novel deep-rooting haplogroup M43 and identified the rare haplogroup H14 in Indian populations carrying the 9-bp deletion by complete mtDNA sequencing. Moreover, we redefined haplogroup M6 and dissected it into two well-defined subclades. The presence of haplogroups F1 and B5a in Uttar Pradesh suggests minor maternal contribution from Southeast Asia to Northern India. The occurrence of haplogroup F1 in the Nepalese sample implies that Nepal might have served as a bridge for the flow of eastern lineages to India. The presence of R6 in the Nepalese, on the other hand, suggests that the gene flow between India and Nepal has been reciprocal.     

mtDNA control-region sequence variation suggests multiple independent origins of an "Asian-specific" 9-bp deletion in sub-Saharan Africans.

The intergenic COII/tRNA(Lys) 9-bp deletion in human mtDNA, which is found at varying frequencies in Asia, Southeast Asia, Polynesia, and the New World, was also found in 81 of 919 sub-Saharan Africans. Using mtDNA control-region sequence data from a subset of 41 individuals with the deletion, we identified 22 unique mtDNA types associated with the deletion in Africa. A comparison of the unique mtDNA types from sub-Saharan Africans and Asians with the 9-bp deletion revealed that sub-Saharan Africans and Asians have sequence profiles that differ in the locations and frequencies of variant sites. Both phylogenetic and mismatch-distribution analysis suggest that 9-bp deletion arose independently in sub-Saharan Africa and Asia and that the deletion has arisen more than once in Africa. Within Africa, the deletion was not found among Khoisan peoples and was rare to absent in western and southwestern African populations, but it did occur in Pygmy and Negroid populations from central Africa and in Malawi and southern African Bantu-speakers. The distribution of the 9-bp deletion in Africa suggests that the deletion could have arisen in central Africa and was then introduced to southern Africa via the recent "Bantu expansion."     

Mitochondrial DNA polymorphisms in nine aboriginal groups of Taiwan: implications for the population history of aboriginal Taiwanese

Mitochondrial DNA (mtDNA) polymorphisms in the D-loop region and the intergenic COII/tRNA(Lys) 9-bp deletion were examined in 180 individuals from all nine aboriginal Taiwanese groups: Atayal, Saisiat, Bunun, Tsou, Rukai, Paiwan, Ami, Puyuma, and Yami. A comparison of 563-bp sequences showed that there were 61 different sequence types, of which 42 types were specific to respective aboriginal groups. D-loop sequence variation and phylogenetic analysis enabled the 180 aboriginal lineages to be classified into eight monophyletic clusters (designated C1-C8). Phylogeographic analysis revealed that two (C2 and C4) of the eight clusters were new characteristic clusters of aboriginal Taiwanese and accounted for 8.3% and 13.9% of the aboriginal lineages, respectively. From the estimated coalescent times for the two unique clusters, the mtDNA lineages leading to such clusters were inferred to have been introduced into Taiwan approximately 11,000-26,000 years ago, suggesting ancient immigrations of the two mtDNA lineages. Genetic distances, based on net nucleotide diversities between populations, revealed three distinct clusters that were comprised of northern mountain (Atayal and Saisiat), southern mountain (Rukai and Paiwan), and middle mountain/east coast (Bunun, Tsou, Ami, Puyuma, and Yami) groups, respectively. Furthermore, phylogenetic analysis of 16 human populations (including six other Asian populations and one African population) confirmed that the three clusters for aboriginal Taiwanese had remained largely intact. Each of the clusters (north, south, and middle-east coast) was characterized by a high frequency of a particular lineage (C4, C2, and 9-bp deletion, respectively). This may result from random genetic drift among the aboriginal groups after a single introduction of all the mtDNA lineages into Taiwan, but another plausible explanation is that at least three genetically distinct ancestral populations have contributed to the maternal gene pool of aboriginal Taiwanese.     

mtDNA polymorphism in East Asian Populations,with special reference to the peopling of Japan.

Nucleotide sequences of the major noncoding (D-loop) region of human mtDNA from five East Asian populations including mainland Japanese, Ainu, Ryukyuans, Koreans, and Chinese were analyzed. On the basis of a comparison of 482-bp sequences in 293 East Asians, 207 different sequence types were observed. Of these, 189 were unique to their respective populations, whereas 18 were shared between two or three populations. Among the shared types, eight were found in common between the mainland Japanese and Koreans, which is the largest number in the comparison. The intergenic COII/tRNA(Lys) 9-bp deletion was observed in every East Asian population with varying frequencies. The D-loop sequence variation suggests that the deletion event occurred only once in the ancestry of East Asians. Phylogenetic analysis revealed that East Asian lineages were classified into at least 18 monophyletic clusters, though lineages from the five populations were completely intermingled in the phylogenetic tree. However, we assigned 14 of the 18 clusters for their specificity on the basis of the population from which the maximum number of individuals in each cluster was derived. Of note is the finding that 50% of the mainland Japanese had continental specificity in which Chinese or Koreans were dominant, while < 20% of either Ryukyuans or Ainu possessed continental specificity. Phylogenetic analysis of the entire human population revealed the closest genetic affinity between the mainland Japanese and Koreans. Thus, the results of this study are compatible with the hybridization model on the origin of modern Japanese. It is suggested that approximately 65% of the gene pool in mainland Japanese was derived from the continental gene flow after the Yayoi Age.     

Genetic evidence for the proto-Austronesian homeland in Asia: mtDNA and nuclear DNA variation in Taiwanese aboriginal tribes.

Previous studies of mtDNA variation in indigenous Taiwanese populations have suggested that they held an ancestral position in the spread of mtDNAs throughout Southeast Asia and Oceania (Melton et al. 1995; Sykes et al. 1995), but the question of an absolute proto-Austronesian homeland remains. To search for Asian roots for indigenous Taiwanese populations, 28 mtDNAs representative of variation in four tribal groups (Ami, Atayal, Bunun, and Paiwan) were sequenced and were compared with each other and with mtDNAs from 25 other populations from Asia and Oceania. In addition, eight polymorphic Alu insertion loci were analyzed, to determine if the pattern of mtDNA variation is concordant with nuclear DNA variation. Tribal groups shared considerable mtDNA sequence identity (P>.90), where gene flow is believed to have been low, arguing for a common source or sources for the tribes. mtDNAs with a 9-bp deletion have considerable mainland-Asian diversity and have spread to Southeast Asia and Oceania through a Taiwanese bottleneck. Only four Taiwanese mtDNA haplotypes without the 9-bp deletion were shared with any other populations, but these shared types were widely dispersed geographically throughout mainland Asia. Phylogenetic and principal-component analyses of Alu loci were concordant with conclusions from the mtDNA analyses; overall, the results suggest that the Taiwanese have temporally deep roots, probably in central or south China, and have been isolated from other Asian populations in recent history.     

A Russian family of Slavic origin carrying mitochondrial DNA with a 9-bp deletion in region V and a long C-stretch in D-loop

A 9-bp deletion first described in the mitochondrial DNA (mtDNA) for East Asian, Polynesian or Indian American populations of the B haplogroup is now discovered in Slavs. The Russian family carrying that deletion belongs to a new branch of the T haplogroup as deduced from D-loop sequence and haplogroup-specific restriction fragment length polymorphism analysis. One family member had a Kearns-Sayre syndrome with a 5.5 kb mtDNA deletion. This family also presented a long C-stretch in the D-loop. Whether or not the formation of the 5.5 kb deletion might be related to the 9-bp deletion or to the long C-stretch in the D-loop is discussed.     

Updating Phylogeny of Mitochondrial DNA Macrohaplogroup M in India: Dispersal of Modern Human in South Asian Corridor

To construct maternal phylogeny and prehistoric dispersals of modern human being in the Indian sub continent, a diverse subset of 641 complete mitochondrial DNA (mtDNA) genomes belonging to macrohaplogroup M was chosen from a total collection of 2,783 control-region sequences, sampled from 26 selected tribal populations of India. On the basis of complete mtDNA sequencing, we identified 12 new haplogroups - M53 to M64; redefined/ascertained and characterized haplogroups M2, M3, M4, M5, M6, M8′C′Z, M9, M10, M11, M12-G, D, M18, M30, M33, M35, M37, M38, M39, M40, M41, M43, M45 and M49, which were previously described by control and/or coding-region polymorphisms. Our results indicate that the mtDNA lineages reported in the present study (except East Asian lineages M8′C′Z, M9, M10, M11, M12-G, D ) are restricted to Indian region.The deep rooted lineages of macrohaplogroup ‘M’ suggest in-situ origin of these haplogroups in India. Most of these deep rooting lineages are represented by multiple ethnic/linguist groups of India. Hierarchical analysis of molecular variation (AMOVA) shows substantial subdivisions among the tribes of India (Fst = 0.16164). The current Indian mtDNA gene pool was shaped by the initial settlers and was galvanized by minor events of gene flow from the east and west to the restricted zones. Northeast Indian mtDNA pool harbors region specific lineages, other Indian lineages and East Asian lineages. We also suggest the establishment of an East Asian gene in North East India through admixture rather than replacement.     

The genetic heritage of the earliest settlers persists both in Indian tribal and caste populations

Two tribal groups from southern India--the Chenchus and Koyas--were analyzed for variation in mitochondrial DNA (mtDNA), the Y chromosome, and one autosomal locus and were compared with six caste groups from different parts of India, as well as with western and central Asians. In mtDNA phylogenetic analyses, the Chenchus and Koyas coalesce at Indian-specific branches of haplogroups M and N that cover populations of different social rank from all over the subcontinent. Coalescence times suggest early late Pleistocene settlement of southern Asia and suggest that there has not been total replacement of these settlers by later migrations. H, L, and R2 are the major Indian Y-chromosomal haplogroups that occur both in castes and in tribal populations and are rarely found outside the subcontinent. Haplogroup R1a, previously associated with the putative Indo-Aryan invasion, was found at its highest frequency in Punjab but also at a relatively high frequency (26%) in the Chenchu tribe. This finding, together with the higher R1a-associated short tandem repeat diversity in India and Iran compared with Europe and central Asia, suggests that southern and western Asia might be the source of this haplogroup. Haplotype frequencies of the MX1 locus of chromosome 21 distinguish Koyas and Chenchus, along with Indian caste groups, from European and eastern Asian populations. Taken together, these results show that Indian tribal and caste populations derive largely from the same genetic heritage of Pleistocene southern and western Asians and have received limited gene flow from external regions since the Holocene. The phylogeography of the primal mtDNA and Y-chromosome founders suggests that these southern Asian Pleistocene coastal settlers from Africa would have provided the inocula for the subsequent differentiation of the distinctive eastern and western Eurasian gene pools.     

Length variations in the COII-tRNA(Lys) intergenic region of mitochondrial DNA in Indonesian populations.

The prevalence of a 9-base-pair (bp) deletion between the mitochondrial cytochrome oxidase II (MTCOX*2) and lysine tRNA (MTTK) genes (region V) has been used to estimate the genetic relationships among Asian and Pacific populations. Many East Asian and Pacific Island populations have been examined previously, but the mitochondrial DNA (mtDNA) diversity of the intervening Indonesian archipelago has not previously been systematically examined. The 17,500 islands of Indonesia currently contain nearly 213 million people and extensive cultural, linguistic, and, presumably, genetic diversity. This study of 1091 individuals representing 15 ethnic groups is the most extensive mtDNA survey to date of the Indonesian archipelago. Six distinct length polymorphisms in region V were observed within these 15 populations. The 9-bp deletion was found in every population examined at frequencies comparable to those of previously examined East Asian populations and substantially lower than those in most Pacific Island populations. Despite the inclusion of Austronesian-speaking populations and a Papuan-speaking population, there was no statistically significant heterogeneity in the frequency of the 9-bp deletion among the 15 populations (p = 0.09). These data indicate that substantial gene flow occurred among the populations at some time in the past. Our observations of no significant correlations between genetic and geographic distances (r = -0.04, p = 0.53) coupled with the extensive cultural and linguistic differences currently within the archipelago suggest that little gene flow among neighboring populations has occurred recently.     

Evolutionary history of the COII/tRNALys intergenic 9 base pair deletion in human mitochondrial DNAs from the Pacific

Length changes in human mitochondrial DNA (mtDNA) are potentially useful markers for inferring the evolutionary history of populations. One such length change is a nine base pair (9-bp) deletion that is located in the intergenic region between the COII gene and the Lysine tRNA gene (COII/tRNALys intergenic region). This deletion has been used as a genetic marker to trace descent from peoples of East Asian origin. A geographic cline of the deletion frequency across modern Pacific Islander populations suggests that the deletion may be useful for tracing prehistoric Polynesian origins and affinities. Mitochondrial DNA sequence variation within two variable segments of the control region (CR) permits a number of inferences regarding the evolutionary history of the 9-bp deletion that cannot be determined from frequency data alone. We obtained CR sequences from 74 mtDNAs with the 9-bp deletion from Indonesia, coastal Papua New Guinea (PNG), and American Samoa. Phylogenetic and pairwise distribution analysis of these CR sequences pooled with previously published CR sequences reveals that the deletion arose independently in Africa and Asia and suggests possible multiple origins of the deletion in Asia. A clinal increase of the frequency of the 9-bp deletion across the three Pacific populations is associated with a decrease in CR sequence diversity, consistent with founder events. Furthermore, analysis of pairwise difference distributions indicates an expansion time of proto-Polynesians that began 5,500 yr ago from Southeast Asia. These results are consistent with the express train model of Polynesian origins.      

C3 polymorphism in some Indian populations.

The distribution of C'3 phenotypes was studied in one tribal and three urban populations from India. The C'3F gene was found low in frequency compared to European and West Asian populations. Quantitatively also, the concentration of the C3 component in the Indian region was found significantly low to the European and West Asian populations reported previously.     

Structure and diversity of the mitochondrial gene pool of the aboriginal population of Tuva and Buriatia from restriction polymorphism data

The populations of Tuvinians (N = 36) and Buryats (N = 105) were characterized by using the data on mitochondrial DNA (mtDNA) polymorphism. The gene pools of both ethnic groups possessed the mtDNA types belonging to the four main haplogroups, A, B, C, and D, found only in the indigenous populations of Asia and America. The total frequencies of the A, B, C, and D haplogroups in Tuvinians and Buryats were 72.3% and 52.4%, respectively. These values, along with the frequency for Altai populations (57.2%), were highest in the Asian populations studied, indicating that the populations Southern and Eastern Siberia can be considered as ancestral relatives to the ethnic groups of the New World. Analysis of the mtDNA region V polymorphism showed the presence of 9-bp deletion and 4-bp insertion in both populations with frequencies respectively of 13.9 and 5.56% in Tuvinians and 4.8 and 1.9% in Buryats. The frequency of the +AvaII/8249 variant was 11.1% in Tuvinians and 3.81% in Buryats. Analysis of the association between the region V deletion-insertion polymorphism and certain restriction haplogroups pointed to repeated and independent emergence of the 4-bp insertion in Siberia.     

Austro-Asiatic tribes of Northeast India provide hitherto missing genetic link between South and Southeast Asia

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.     

Most of the extant mtDNA boundaries in South and Southwest Asia were likely shaped during the initial settlement of Eurasia by anatomically modern humans

Background

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.

Results

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.

Conclusions

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.     

Peopling of Sahul: mtDNA variation in aboriginal Australian and Papua New Guinean populations.

We examined genetic affinities of Aboriginal Australian and New Guinean populations by using nucleotide variation in the two hypervariable segments of the mtDNA control region (CR). A total of 318 individuals from highland Papua New Guinea (PNG), coastal PNG, and Aboriginal Australian populations were typed with a panel of 29 sequence-specific oligonucleotide (SSO) probes. The SSO-probe panel included five new probes that were used to type an additional 1,037 individuals from several Asian populations. The SSO-type data guided the selection of 78 individuals from Australia and east Indonesia for CR sequencing. A gene tree of these CR sequences, combined with published sequences from worldwide populations, contains two previously identified highland PNG clusters that do not include any Aboriginal Australians; the highland PNG clusters have coalescent time estimates of approximately 80,000 and 122,000 years ago, suggesting ancient isolation and genetic drift. SSO-type data indicate that 84% of the sample of PNG highlander mtDNA belong to these two clusters. In contrast, the Aboriginal Australian sequences are intermingled throughout the tree and cluster with sequences from multiple populations. Phylogenetic and multidimensional-scaling analyses of CR sequences and SSO types split PNG highland and Aboriginal Australian populations and link Aboriginal Australian populations with populations from the subcontinent of India. These mtDNA results do not support a close relationship between Aboriginal Australian and PNG populations but instead suggest multiple migrations in the peopling of Sahul.     

Genomic structures and population histories of linguistically distinct tribal groups of India

There are various conflicting hypotheses regarding the origins of the tribal groups of India, who belong to three major language groups--Austro-Asiatic, Dravidian and Tibeto-Burman. To test some of the major hypotheses we designed a genetic study in which we sampled tribal populations belonging to all the three language groups. We used a set of autosomal DNA markers, mtDNA restriction-site polymorphisms (RSPs) and mtDNA hypervariable segment-1 (HVS-1) sequence polymorphisms in this study. Using the unlinked autosomal markers we found that there is a fair correspondence between linguistic and genomic affinities among the Indian tribal groups. We reconstructed mtDNA RSP haplotypes and found that there is extensive haplotype sharing among all tribal populations. However, there is very little sharing of mtDNA HVS-1 sequences across populations, and none across language groups. Haplogroup M is ubiquitous, and the subcluster U2i of haplogroup U occurs in a high frequency. Our analyses of haplogroup and HVS-1 sequence data provides evidence in support of the hypothesis that the Austro-Asiatic speakers are the most ancient inhabitants of India. Our data also support the earlier finding that some of the western Eurasian haplogroups found in India may have been present in India prior to the entry of Aryan speakers. However, we do not find compelling evidence to support the theory that haplogroup M was brought into India on an "out of Africa" wave of migration through a southern exit route from Ethiopia. On the contrary, our data raise the possibility that this haplogroup arose in India and was later carried to East Africa from India.     

The northeast Indian passageway: a barrier or corridor for human migrations?

The northeast Indian passageway connecting the Indian subcontinent to East/Southeast Asia is thought to have been a major corridor for human migrations. Because it is also an important linguistic contact zone, it is predicted that northeast India has witnessed extensive population interactions, thus, leading to high genetic diversity within groups and heterogeneity among groups. To test this prediction, we analyzed 14 biallelic and five short tandem-repeat Y-chromosome markers and hypervariable region 1 mtDNA sequence variation in 192 northeast Indians. We find that both northeast Indian Y chromosomes and mtDNAs consistently show strikingly high homogeneity among groups and strong affinities to East Asian groups. We detect virtually no Y-chromosome and mtDNA admixture between northeast and other Indian groups. Northeast Indian groups are also characterized by a greatly reduced Y-chromosome diversity, which contrasts with extensive mtDNA diversity. This is best explained by a male founder effect during the colonization of northeast India that is estimated to have occurred within the past 4,000 years. Thus, contrary to the prediction, these results provide strong evidence for a genetic discontinuity between northeast Indian groups and other Indian groups. We, therefore, conclude that the northeast Indian passage way acted as a geographic barrier rather than as a corridor for human migrations between the Indian subcontinent and East/Southeast Asia, at least within the past millennia and possibly for several tens of thousand years, as suggested by the overall distinctiveness of the Indian and East Asian Y chromosome and mtDNA gene pools.