Integrating Linguistics,Social Structure,and Geography to Model Genetic Diversity within India

India represents an intricate tapestry of population substructure shaped by geography, language, culture, and social stratification. Although geography closely correlates with genetic structure in other parts of the world, the strict endogamy imposed by the Indian caste system and the large number of spoken languages add further levels of complexity to understand Indian population structure. To date, no study has attempted to model and evaluate how these factors have interacted to shape the patterns of genetic diversity within India. We merged all publicly available data from the Indian subcontinent into a data set of 891 individuals from 90 well-defined groups. Bringing together geography, genetics, and demographic factors, we developed Correlation Optimization of Genetics and Geodemographics to build a model that explains the observed population genetic substructure. We show that shared language along with social structure have been the most powerful forces in creating paths of gene flow in the subcontinent. Furthermore, we discover the ethnic groups that best capture the diverse genetic substructure using a ridge leverage score statistic. Integrating data from India with a data set of additional 1,323 individuals from 50 Eurasian populations, we find that Indo-European and Dravidian speakers of India show shared genetic drift with Europeans, whereas the Tibeto-Burman speaking tribal groups have maximum shared genetic drift with East Asians.     

Cucurbits, Sanskrit, and the Indo-Aryas 1

Vedic and early post-Vedic Sanskrit works describe Indo-Aryan life in northern India during 2000–200 B.C. Wild and cultivated plants of the Cucurbitaceae played an important role in Indo-Aryan food, medicine, and culture; over 300 words describing cucurbits are found in the Sanskrit texts. These words, their etymological relationships, and the geographies of the cucurbit taxa were compared. Results indicate that the Indo-Aryas knew only a handful of cucurbits before entering the Indian subcontinent from the west. They learned of at least 11 more species from their neighbors in India, who were speakers of Dravidian or Munda languages.     

The Austroasiatic Munda population from India and Its enigmatic origin: a HLA diversity study

The Austroasiatic linguistic family disputes its origin between two geographically distant regions of Asia, India, and Southeast Asia, respectively. As genetic studies based on classical and gender-specific genetic markers provided contradictory results to this debate thus far, we investigated the HLA diversity (HLA-A, -B, and -DRB1 loci) of an Austroasiatic Munda population from Northeast India and its relationships with other populations from India and Southeast Asia. Because molecular methods currently used to test HLA markers often provide ambiguous results due to the high complexity of this polymorphism, we applied two different techniques (reverse PCR-SSO typing on microbeads arrays based on Luminex technology, and PCR-SSP typing) to type the samples. After validating the resulting frequency distributions through the original statistical method described in our companion article ( Nunes et al. 2011 ), we compared the HLA genetic profile of the sampled Munda to those of other Asiatic populations, among which Dravidian and Indo-European-speakers from India and populations from East and Southeast Asia speaking languages belonging to different linguistic families. We showed that the Munda from Northeast India exhibit a peculiar genetic profile with a reduced level of HLA diversity compared to surrounding Indian populations. They also exhibit less diversity than Southeast Asian populations except at locus DRB1. Several analyses using genetic distances indicate that the Munda are much more closely related to populations from the Indian subcontinent than to Southeast Asian populations speaking languages of the same Austroasiatic linguistic family. On the other hand, they do not share a closer relationship with Dravidians compared with Indo-Europeans, thus arguing against the idea that the Munda share a common and ancient Indian origin with Dravidians. Our results do not favor either a scenario where the Munda would be representative of an ancestral Austroasiatic population giving rise to an eastward Austroasiatic expansion to Southeast Asia. Rather, their peculiar genetic profile is better explained by a decrease in genetic diversity through genetic drift from an ancestral population having a genetic profile similar to present-day Austroasiatic populations from Southeast Asia (thus suggesting a possible southeastern origin), followed by intensive gene flow with neighboring Indian populations. This conclusion is in agreement with archaeological and linguistic information. The history of the Austroasiatic family represents a fascinating example where complex interactions among culturally distinct human populations occurred in the past.     

Genetic affinities of the central Indian tribal populations

Background

The central Indian state Madhya Pradesh is often called as ‘heart of India’ and has always been an important region functioning as a trinexus belt for three major language families (Indo-European, Dravidian and Austroasiatic). There are less detailed genetic studies on the populations inhabited in this region. Therefore, this study is an attempt for extensive characterization of genetic ancestries of three tribal populations, namely; Bharia, Bhil and Sahariya, inhabiting this region using haploid and diploid DNA markers.

Methodology/Principal Findings

Mitochondrial DNA analysis showed high diversity, including some of the older sublineages of M haplogroup and prominent R lineages in all the three tribes. Y-chromosomal biallelic markers revealed high frequency of Austroasiatic-specific M95-O2a haplogroup in Bharia and Sahariya, M82-H1a in Bhil and M17-R1a in Bhil and Sahariya. The results obtained by haploid as well as diploid genetic markers revealed strong genetic affinity of Bharia (a Dravidian speaking tribe) with the Austroasiatic (Munda) group. The gene flow from Austroasiatic group is further confirmed by their Y-STRs haplotype sharing analysis, where we determined their founder haplotype from the North Munda speaking tribe, while, autosomal analysis was largely in concordant with the haploid DNA results.

Conclusions/Significance

Bhil exhibited largely Indo-European specific ancestry, while Sahariya and Bharia showed admixed genetic package of Indo-European and Austroasiatic populations. Hence, in a landscape like India, linguistic label doesn''t unequivocally follow the genetic footprints.     

Population genetic analysis among five Indian population groups using six microsatellite markers

Genetic variation at six tetranucleotide microsatellites (HUMTHO1, HUMVWA, F13A01, D3S1359, D12S66, and D12S67) has heen determined in five endogamous ethnic population groups of India belonging to two major linguistic families. The populations analyzed were Konkanastha Brahmins and Marathas (Maharashtra state) from the Indo-Aryan linguistic family and Nairs, Ezhavas, and Muslims (Kerala state) from the Dravidian family. All six loci show high gene diversity, ranging from 0.63 +/- 0.04 to 0.84 +/- 0.02. The average GST value observed was 1.7%, indicating that the differences between the populations account for less than 2% of the diversity, while the genetic variation is high within the five population groups studied (>98%). The phylogenetic tree fails to show any clear cluster. The absence of any cluster along with low average GST is suggestive of substantial genetic similarity among the studied populations, in spite of clear geographical, linguistic, and cultural barriers. This similarity indicates either a greater gene flow between these groups or, alternatively, may reflect a recent evolution for them, considering that the Indian caste system evolved only about 3000 years ago.     

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.     

Shared and unique components of human population structure and genome-wide signals of positive selection in South Asia

South Asia harbors one of the highest levels genetic diversity in Eurasia, which could be interpreted as a result of its long-term large effective population size and of admixture during its complex demographic history. In contrast to Pakistani populations, populations of Indian origin have been underrepresented in previous genomic scans of positive selection and population structure. Here we report data for more than 600,000 SNP markers genotyped in 142 samples from 30 ethnic groups in India. Combining our results with other available genome-wide data, we show that Indian populations are characterized by two major ancestry components, one of which is spread at comparable frequency and haplotype diversity in populations of South and West Asia and the Caucasus. The second component is more restricted to South Asia and accounts for more than 50% of the ancestry in Indian populations. Haplotype diversity associated with these South Asian ancestry components is significantly higher than that of the components dominating the West Eurasian ancestry palette. Modeling of the observed haplotype diversities suggests that both Indian ancestry components are older than the purported Indo-Aryan invasion 3,500 YBP. Consistent with the results of pairwise genetic distances among world regions, Indians share more ancestry signals with West than with East Eurasians. However, compared to Pakistani populations, a higher proportion of their genes show regionally specific signals of high haplotype homozygosity. Among such candidates of positive selection in India are MSTN and DOK5, both of which have potential implications in lipid metabolism and the etiology of type 2 diabetes.     

A microsatellite study to disentangle the ambiguity of linguistic,geographic, ethnic and genetic influences on tribes of India to get a better clarity of the antiquity and peopling of South Asia

An understanding of the genetic affinity and the past history of the tribal populations of India requires the untangling of the confounding influences of language, ethnicity, and geography on the extant diverse tribes. The present study examines the genetic relationship of linguistically (Dravidian, Austro‐Asiatic, and Tibeto‐Burman) and ethnically (Australian and East Asian) diverse tribal populations (46) inhabiting different regions of the Indian subcontinent. For the purpose, we have utilized the published data on allele frequency of 15 autosomal STR loci of our study on six Adi sub‐tribes of Arunachal Pradesh and compared the same with the reported allele frequency data, for nine common autosomal STR loci, of 40 other tribes. Phylogenetic and principal component analyses exhibit geography based clustering of Tibeto‐Burman speakers and separation of the Mundari and Mon‐Khmer speaking Austro‐Asiatic populations. The combined analyses of all 46 populations show clustering of the groups belonging to same ethnicity and inhabiting contiguous geographic regions, irrespective of their different languages. These results help us to reconstruct and understand three plausible scenarios of the antiquity of Indian tribal populations: the Dravidian and Austro‐Asiatic (Mundari) tribes were possibly derived from common early settlers; the Tibeto‐Burman tribes possibly belonged to a different ancestry and the Mon‐Khmer speaking Austro‐Asiatic populations share a common ancestry with some of the Tibeto‐Burman speakers. Am J Phys Anthropol, 2009. © 2009 Wiley‐Liss, Inc.     

Reconstructing the Indian origin and dispersal of the European Roma: a maternal genetic perspective

Previous genetic, anthropological and linguistic studies have shown that Roma (Gypsies) constitute a founder population dispersed throughout Europe whose origins might be traced to the Indian subcontinent. Linguistic and anthropological evidence point to Indo-Aryan ethnic groups from North-western India as the ancestral parental population of Roma. Recently, a strong genetic hint supporting this theory came from a study of a private mutation causing primary congenital glaucoma. In the present study, complete mitochondrial control sequences of Iberian Roma and previously published maternal lineages of other European Roma were analyzed in order to establish the genetic affinities among Roma groups, determine the degree of admixture with neighbouring populations, infer the migration routes followed since the first arrival to Europe, and survey the origin of Roma within the Indian subcontinent. Our results show that the maternal lineage composition in the Roma groups follows a pattern of different migration routes, with several founder effects, and low effective population sizes along their dispersal. Our data allowed the confirmation of a North/West migration route shared by Polish, Lithuanian and Iberian Roma. Additionally, eleven Roma founder lineages were identified and degrees of admixture with host populations were estimated. Finally, the comparison with an extensive database of Indian sequences allowed us to identify the Punjab state, in North-western India, as the putative ancestral homeland of the European Roma, in agreement with previous linguistic and anthropological studies.     

Microsatellite diversity reveals the interplay of language and geography in shaping genetic differentiation of diverse Proto-Australoid populations of west-central India

Microsatellite diversity was analyzed in four Proto-Australoid tribes, including Indo-European (Marathi)-speaking Katkari, Pawara, Mahadeo-Koli, and Dravidian (Gondi)-speaking groups of Maharashtra, west-central India, to understand their genetic structure and to identify the congruence between language and gene pool. Allele frequency data at 15 short tandem repeat (STR) loci in studied tribes was compared with data of 22 Indo-European- and Dravidian-speaking caste and tribal populations using heterozygosity, allele size variance, analysis of molecular variance (AMOVA), G(ST) estimate, PC plot, and Mantel correlation test. Our results demonstrate that "Gondi" tribes comprising the Madia-Gond, a hunter-gatherer population, and the agriculturist Dheria-Gond harbor lower diversity than "Marathi" tribal groups, which are culturally and genetically distinct. Katkari, a hunter-gatherer tribe, showed greater diversity and the presence of a large number of unique alleles, genetically distinct from all others except the Pawara, supporting their old cultural links. The agriculturist Pawara tribe represents a splinter subgroup of the Bhil tribe and has experienced gene flow. The Mahadeo-Koli, an agriculturally oriented tribe, displayed significant heterozygote deficiency, attributable to the practice of high endogamy. The Proto-Australoid tribal populations were genetically differentiated from castes of similar morphology, suggesting different evolutionary mechanisms operating upon the populations. The populations showed genetic and linguistic similarity, barring a few groups with varied migratory histories. The microsatellite variation clearly demonstrates the interplay of sociocultural factors including linguistic, geographical contiguity, and microevolutionary processes in shaping the genetic diversity of populations in contemporary India. This study supports the ethno-historical relationships of Indian populations.     

Genetic landscape of the people of India: a canvas for disease gene exploration

Analyses of frequency profiles of markers on disease or drug-response related genes in diverse populations are important for the dissection of common diseases. We report the results of analyses of data on 405 SNPs from 75 such genes and a 5.2 Mb chromosome, 22 genomic region in 1871 individuals from diverse 55 endogamous Indian populations. These include 32 large (> 10 million individuals) and 23 isolated populations, representing a large fraction of the people of India. We observe high levels of genetic divergence between groups of populations that cluster largely on the basis of ethnicity and language. Indian populations not only overlap with the diversity of HapMap populations, but also contain population groups that are genetically distinct. These data and results are useful for addressing stratification and study design issues in complex traits especially for heterogeneous populations. For complete authors’ list of the Indian Genome Variation Database Consortium please see Appendix. For correspondence. E-mail: skb@igib.res.in; Functional Genomics Unit, Institute of Genomics and Integrative Biology (CSIR), Mall Road, New Delhi 110 007, India.     

Genetic Affinity of the Bhil,Kol and Gond Mentioned in Epic Ramayana

Kol, Bhil and Gond are some of the ancient tribal populations known from the Ramayana, one of the Great epics of India. Though there have been studies about their affinity based on classical and haploid genetic markers, the molecular insights of their relationship with other tribal and caste populations of extant India is expected to give more clarity about the the question of continuity vs. discontinuity. In this study, we scanned >97,000 of single nucleotide polymorphisms among three major ancient tribes mentioned in Ramayana, namely Bhil, Kol and Gond. The results obtained were then compared at inter and intra population levels with neighboring and other world populations. Using various statistical methods, our analysis suggested that the genetic architecture of these tribes (Kol and Gond) was largely similar to their surrounding tribal and caste populations, while Bhil showed closer affinity with Dravidian and Austroasiatic (Munda) speaking tribes. The haplotype based analysis revealed a massive amount of genome sharing among Bhil, Kol, Gond and with other ethnic groups of South Asian descent. On the basis of genetic component sharing among different populations, we anticipate their primary founding over the indigenous Ancestral South Indian (ASI) component has prevailed in the genepool over the last several thousand years.     

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.     

Whole genome variant analysis in three ethnically diverse Indians

India represents an amazing confluence of geographically, linguistically and socially disparate ethnic populations (Indian Genome Variation Consortium, J Genet 87:3–20, 2008). Understanding the genetic diversity of Indian population remains a daunting task. In this paper we present detailed analysis of genomic variations (high-depth coverage (~?30×) using Illumina Hiseq 2000 platform) from three healthy Indian male individuals each belonging to three geographically delineated regions and linguistic phylum viz. high altitude region of Ladakh (Tibeto-Burman linguistic phylum), sub mountainous region of Kumaun (Indo-European linguistic phylum) and sea level region of Telangana (Dravidian linguistic phylum) for probing the extent of genetic diversity in our population. The sequencing analysis provided high quality data (~?95% of the total reads aligned to the human reference genome for each sample) and very good alignment quality (>?80% of the filtered mapped reads had a quality score of 60). A total of 4.3, 3.7 and 4.3 million single nucleotide variations were identified in the genome of high altitude, sub mountainous and sea level respectively by comparing with human reference genome. Approximately 17.3, 18.2, 17.4% of the variants were unique in the three genomes. The study identified many novel variations in the three diverse genomes (132,970 in Ladakh, 112,317 in Kumaun and 128,881 in Telangana individual) and is an important resource for creating a baseline and a comprehensive catalogue of human genomic variation across the Indian as well as the Asian continent.     

Deep common ancestry of indian and western-Eurasian mitochondrial DNA lineages

About a fifth of the human gene pool belongs largely either to Indo-European or Dravidic speaking people inhabiting the Indian peninsula. The 'Caucasoid share' in their gene pool is thought to be related predominantly to the Indo-European speakers. A commonly held hypothesis, albeit not the only one, suggests a massive Indo-Aryan invasion to India some 4,000 years ago [1]. Recent limited analysis of maternally inherited mitochondrial DNA (mtDNA) of Indian populations has been interpreted as supporting this concept [2] [3]. Here, this interpretation is questioned. We found an extensive deep late Pleistocene genetic link between contemporary Europeans and Indians, provided by the mtDNA haplogroup U, which encompasses roughly a fifth of mtDNA lineages of both populations. Our estimate for this split is close to the suggested time for the peopling of Asia and the first expansion of anatomically modern humans in Eurasia [4] [5] [6] [7] [8] and likely pre-dates their spread to Europe. Only a small fraction of the 'Caucasoid-specific' mtDNA lineages found in Indian populations can be ascribed to a relatively recent admixture.     

Haplotype analysis of the polymorphic 17 YSTR markers in Kerala nontribal populations

The origin of the Kerala non tribal population has been a matter of contention for centuries. While some claim that Negritos were the first inhabitants, some historians suggest a Dravidian origin for all Keralites. The aim of our study has been to provide sufficient scientific evidence based on Y chromosome short tandem repeat (Y STR) analysis for tracing the paternal lineage and also to create a database of the Y STR haplotype of the male population for future forensic analysis. Whole blood samples (n = 168) were collected from unrelated healthy men of the Kerala non-tribal population over a period of 2 years from October 2009. Genomic DNA was extracted by salting out method. All samples were genotyped for the 17 Y STR loci by the AmpFLSTR Y-filer PCR Amplification Kit. The haplotype and allele frequencies were determined by direct counting and analyzed using Arlequin 3.1 software, and molecular variance was calculated with the Y chromosome haplotype reference database online analysis tool, . Haplotype diversity was calculated using HaPYDive (). The majority of haplotypes were unique (149/168). The variant allele 17.1 was observed in DYS 385 loci in three samples. Fifteen samples (8.93%) showed the presence of alleles that are not within the established marker range denoted as outside marker range (OMR). The allele frequency of Kerala non tribal population ranged from 0.00003 to 0.5809. The most polymorphic single locus marker was DYS 458. The haplotype diversity value for Kerala non tribal population was 0.9978. The pairwise difference value ranged from 0.0531 to 0.0854 on comparison of the haplotypes of the Kerala non tribals with other Indian populations. The multi dimensional scaling plot depicted the proximity of Kerala non tribal population with Vasterbotten population (Swedish) and Paiwan, Patyal population of Taiwan, Thailand, and Zhuang population of China. The results of the study indicate towards a European paternal lineage in the non tribal Kerala population.     

Anadromy and heterogenous population of a tropical shad Tenualosa ilisha in Malaysia,as revealed by otolith microchemistry and molecular evidence

Tenualosa ilisha was found recently in the Perak River in western Peninsular Malaysia. Molecular phylogenetic and haplotype network analyses suggest that T. ilisha has two genetically distinct populations/groups: (i) Peninsular Malaysia (Malaysia population), and (ii) Peninsular Malaysia, Thailand, India and Bangladesh (Indian Ocean population). The results also suggest that the T ilisha population in Peninsular Malaysia is genetically heterogeneous with a typical anadromous migration pattern.     

Phylogeographic Analysis of Mitochondrial DNA in the Nogays: A Strong Mixture of Maternal Lineages from Eastern and Western Eurasia

Analysis of mtDNA markers in a population of the Nogays (n = 206), the people inhabiting the North Caucasus and speaking a Turkic language of the Altaic linguistic family, has revealed a high level of genetic diversity (H = 0.99). The identified haplotypes include all major West Eurasian haplogroups, with the prevalence of H and U clusters (22 and 21%, respectively), but the percentage of lineages specific for East Eurasian populations is the highest (40%). Some other mtDNA variants in the Nogay population belong to the M1 haplogroups typical of northeastern Africa and U2 characteristic of Indian populations. Thus, components of different origin have contributed to the gene pool of Nogays. An erratum to this article is available at .     

Gametophytic self-incompatibility in Lycium parishii (Solanaceae): allelic diversity, genealogical structure, and patterns of molecular evolution at the S-RNase locus

We characterized allelic diversity at the locus controlling self-incompatibility (SI) for a population of Lycium parishii (Solanaceae) from Organ Pipe National Monument, Arizona. Twenty-four partial sequences of S-RNase alleles were recovered from 25 individuals. Estimates of allelic diversity range from 23 to 27 alleles and, consistent with expectations for SI, individuals are heterozygous. We compare S-RNase diversity, patterns of molecular evolution, and the genealogical structure of alleles from L. parishii to a previously studied population of its congener L. andersonii. Gametophytic SI is well characterized for Solanaceae and although balancing selection is hypothesized to be responsible for high levels of allelic divergence, the pattern of selection varies depending on the portion of the gene considered. Site-specific models investigating patterns of selection for L. parishii and L. andersonii indicate that positive selection occurs in those regions of the S-RNase gene hypothesized as important to the recognition response, whereas positive selection was not detected for any position within regions previously characterized as conserved. A 10-species genealogy including S-RNases from a pair of congeners from each of five genera in Solanaceae reveals extensive transgeneric evolution of L. parishii S-RNases. Further, within Lycium, the Dn/Ds ratios for pairs of closely related alleles for intraspecific versus interspecific comparisons were not significantly different, suggesting that the S-RNase diversity recovered in these two species was present prior to the speciation event separating them. Despite this, two S-RNases from L. parishii are identical to two previously reported alleles for L. andersonii, suggesting gene flow between these species.