The 9-bp deletion between the mitochondrial lysine tRNA and COII genes in tribal populations of India

A 9-base-pair (bp) deletion located between the lysine tRNA (MTTK) and COII (MTCOX*2) genes in the human mitochondrial genome is a valuable marker for tracing population relationships. Previous research has shown that the 9-bp deletion is associated with two major clusters of control region sequences; one occurs in sub-Saharan Africa, while the other is associated with Asian populations and populations of Asian origin. We surveyed 898 individuals from 16 tribal populations in India and found 6 individuals with the 9-bp deletion. Sequences of the first hypervariable segment (HV1) of the mtDNA control region from these 9-bp deletion-bearing mtDNAs were compared to those previously reported from Asian and African populations. Phylogenetic analysis indicates three distinct clusters of tribal Indian 9-bp deletion mtDNA types. One cluster, found in northeast India, includes southeast Asian and Indonesian mtDNA types. The remaining two clusters appear to have unique origins in southern India. These data provide further evidence of past migrations from Asia into the northeast corner of the Indian subcontinent.     

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.     

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.     

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.      

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.     

Peopling of South Asia: investigating the caste-tribe continuum in India

In recent years, mtDNA and Y chromosome studies involving human populations from South Asia and the rest of the world have revealed new insights about the peopling of the world by anatomically modern humans during the late Pleistocene, some 40,000-60,000 years ago, over the southern coastal route from Africa. Molecular studies and archaeological record are both largely consistent with autochthonous differentiation of the genetic structure of the caste and tribal populations in South Asia. High level of endogamy created by numerous social boundaries within and between castes and tribes, along with the influence of several evolutionary forces such as genetic drift, fragmentation and long-term isolation, has kept the Indian populations diverse and distant from each other as well as from other continental populations. This review attempts to summarize recent genetic studies on Indian caste and tribal populations with the focus on the information embedded in the socially defined structure of Indian populations.     

Mitochondrial DNA variation in Nicobarese Islanders.

The aboriginal populations living in the Nicobar Islands are hypothesized to be descendants of people who were part of early human dispersals into Southeast Asia. However, analyses of ethnographic histories, languages, morphometric data, and protein polymorphisms have not yet resolved which worldwide populations are most closely related to the Nicobarese. Thus, to explore the origins and affinities of the Nicobar Islanders, we analyzed mitochondrial DNA (mtDNA) hypervariable region 1 sequence data from 33 Nicobarese Islanders and compared their mtDNA haplotypes to those of neighboring East Asians, mainland and island Southeast Asians, Indians, Australian aborigines, Pacific Islanders, and Africans. Unique Nicobarese mtDNA haplotypes, including five Nicobarese mtDNA haplotypes linked to the COII/tRNA(Lys) 9-bp deletion, are most closely related to mtDNA haplotypes from mainland Southeast Asian Mon-Kmer-speaking populations (e.g., Cambodians). Thus, the dispersal of southern Chinese into mainland Southeast Asia may have included a westward expansion and colonization of the islands of the Andaman Sea.     

Genetic Structure of Tibeto-Burman Populations of Bangladesh: Evaluating the Gene Flow along the Sides of Bay-of-Bengal

Human settlement and migrations along sides of Bay-of-Bengal have played a vital role in shaping the genetic landscape of Bangladesh, Eastern India and Southeast Asia. Bangladesh and Northeast India form the vital land bridge between the South and Southeast Asia. To reconstruct the population history of this region and to see whether this diverse region geographically acted as a corridor or barrier for human interaction between South Asia and Southeast Asia, we, for the first time analyzed high resolution uniparental (mtDNA and Y chromosome) and biparental autosomal genetic markers among aboriginal Bangladesh tribes currently speaking Tibeto-Burman language. All the three studied populations; Chakma, Marma and Tripura from Bangladesh showed strikingly high homogeneity among themselves and strong affinities to Northeast Indian Tibeto-Burman groups. However, they show substantially higher molecular diversity than Northeast Indian populations. Unlike Austroasiatic (Munda) speakers of India, we observed equal role of both males and females in shaping the Tibeto-Burman expansion in Southern Asia. Moreover, it is noteworthy that in admixture proportion, TB populations of Bangladesh carry substantially higher mainland Indian ancestry component than Northeast Indian Tibeto-Burmans. Largely similar expansion ages of two major paternal haplogroups (O2a and O3a3c), suggested that they arose before the differentiation of any language group and approximately at the same time. Contrary to the scenario proposed for colonization of Northeast India as male founder effect that occurred within the past 4,000 years, we suggest a significantly deep colonization of this region. Overall, our extensive analysis revealed that the population history of South Asian Tibeto-Burman speakers is more complex than it was suggested before.     

The influence of natural barriers in shaping the genetic structure of Maharashtra populations

Background

The geographical position of Maharashtra state makes it rather essential to study the dispersal of modern humans in South Asia. Several hypotheses have been proposed to explain the cultural, linguistic and geographical affinity of the populations living in Maharashtra state with other South Asian populations. The genetic origin of populations living in this state is poorly understood and hitherto been described at low molecular resolution level.

Methodology/Principal Findings

To address this issue, we have analyzed the mitochondrial DNA (mtDNA) of 185 individuals and NRY (non-recombining region of Y chromosome) of 98 individuals belonging to two major tribal populations of Maharashtra, and compared their molecular variations with that of 54 South Asian contemporary populations of adjacent states. Inter and intra population comparisons reveal that the maternal gene pool of Maharashtra state populations is composed of mainly South Asian haplogroups with traces of east and west Eurasian haplogroups, while the paternal haplogroups comprise the South Asian as well as signature of near eastern specific haplogroup J2a.

Conclusions/Significance

Our analysis suggests that Indian populations, including Maharashtra state, are largely derived from Paleolithic ancient settlers; however, a more recent (∼10 Ky older) detectable paternal gene flow from west Asia is well reflected in the present study. These findings reveal movement of populations to Maharashtra through the western coast rather than mainland where Western Ghats-Vindhya Mountains and Narmada-Tapti rivers might have acted as a natural barrier. Comparing the Maharastrian populations with other South Asian populations reveals that they have a closer affinity with the South Indian than with the Central Indian populations.     

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.     

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.     

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."     

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.     

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.     

Land, language, and loci: mtDNA in Native Americans and the genetic history of Peru

Despite a long history of complex societies and despite extensive present-day linguistic and ethnic diversity, relatively few populations in Peru have been sampled for population genetic investigations. In order to address questions about the relationships between South American populations and about the extent of correlation between genetic distance, language, and geography in the region, mitochondrial DNA (mtDNA) hypervariable region I sequences and mtDNA haplogroup markers were examined in 33 individuals from the state of Ancash, Peru. These sequences were compared to those from 19 American Indian populations using diversity estimates, AMOVA tests, mismatch distributions, a multidimensional scaling plot, and regressions. The results show correlations between genetics, linguistics, and geographical affinities, with stronger correlations between genetics and language. Additionally, the results suggest a pattern of differential gene flow and drift in western vs. eastern South America, supporting previous mtDNA and Y chromosome investigations.     

Independent origins of Indian caste and tribal paternal lineages

The origins of the nearly one billion people inhabiting the Indian subcontinent and following the customs of the Hindu caste system are controversial: are they largely derived from Indian local populations (i.e. tribal groups) or from recent immigrants to India? Archaeological and linguistic evidence support the latter hypothesis, whereas recent genetic data seem to favor the former hypothesis. Here, we analyze the most extensive dataset of Indian caste and tribal Y chromosomes to date. We find that caste and tribal groups differ significantly in their haplogroup frequency distributions; caste groups are homogeneous for Y chromosome variation and more closely related to each other and to central Asian groups than to Indian tribal or any other Eurasian groups. We conclude that paternal lineages of Indian caste groups are primarily descended from Indo-European speakers who migrated from central Asia approximately 3,500 years ago. Conversely, paternal lineages of tribal groups are predominantly derived from the original Indian gene pool. We also provide evidence for bidirectional male gene flow between caste and tribal groups. In comparison, caste and tribal groups are homogeneous with respect to mitochondrial DNA variation, which may reflect the sociocultural characteristics of the Indian caste society.     

Mitochondrial diversity and phylogeographic structure of Chinese domestic goats

China has numerous native domestic goat breeds, but so far there has been no extensive study on genetic diversity, population demographic history, and origin of Chinese goats. Here, we examined the genetic diversity and phylogeographic structure of Chinese domestic goats by determining a 481-bp fragment of the first hypervariable region of mitochondrial DNA (mtDNA) control region from 368 individuals representing 18 indigenous breeds. Phylogenetic analyses revealed that there were four mtDNA lineages (A-D) identified in Chinese goats, in which lineage A was predominant, lineage B was moderate, and lineages C and D were at low frequency. These results further support the multiple maternal origins of domestic goats. The pattern of genetic variation in goat mtDNA sequences indicated that the two larger lineages A and B had undergone population expansion events. In a combined analysis of previously reported sequences and our sequences belonging to lineage B, we detected two subclades, in which one was unique to eastern Asia and another was shared between eastern and southern Asia. A larger genetic variation in eastern Asia than southern Asia and the pattern of phylogeographic variation in lineage B suggest that at least one subclade of lineage B originated from eastern Asia. There was no significant geographical structuring in Chinese goat populations, which suggested that there existed strong gene flow among goat populations caused by extensive transportation of goats in history.     

Evolutionary history of the mtDNA 9-bp deletion in Chinese populations and its relevance to the peopling of east and southeast Asia

In total, 1218 Chinese from twelve ethnic groups and nine Han geographic groups were screened for the mtDNA 9-bp deletion motif. The frequency of the 9-bp deletion in all samples was 14.7% but ranged from 0% to 32% in the various ethnic groups. Three individuals had a triplication of the 9-bp segment. Phylogenetic and demographic analyses of the mtDNA hypervariable segment 1 (HVS 1) sequences suggest that the 9-bp deletion occurred more than once in China. The majority of the Chinese deletion haplotypes (about 90%) have a common origin as a mutational event following an initial expansion of modern humans in eastern Asia. Other deletion haplotypes and the three haplotypes with a 9-bp triplication may have arisen independently in the Chinese, presumably by replication error. HVS1 haplotype analysis suggests two possible migration routes of the 9-bp deletion in east and southeast Asia. Both migrations originated in China with one route leading to the Pacific Islands via Taiwan, the other to southeast Asia and possibly the Nicobar Islands. Along both routes of peopling, a decrease in HVSI diversity of the mtDNA haplotypes is observed. The "Polynesian motif (16217T/C, 16247A/G, and 16261C/T)" and the 16140T/C, 16266C/A, or C/G polymorphisms appear specific to each migration route.