A prehistory of Indian Y chromosomes

Note: Performing your original search, dna india y mitochondrial, in Pub Med Central will retrieve 300 citations. Journal List > Proc Natl Acad Sci U S A > v.103(4); Jan 24, 2006

Anthropology http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=1347984

A prehistory of Indian Y chromosomes: Evaluating demic diffusion scenarios

Sanghamitra Sahoo,† Anamika Singh,† G. Himabindu,† Jheelam Banerjee,† T. Sitalaximi,† Sonali Gaikwad,† R. Trivedi,† Phillip Endicott,‡ Toomas Kivisild,§ Mait Metspalu,§ Richard Villems,§ and V. K. Kashyap†¶

†National DNA Analysis Centre, Central Forensic Science Laboratory, Kolkata 700014, India; ‡Department of Zoology, University of Oxford, Oxford OX1 3PS, United Kingdom; §Estonian Biocentre, 51010 Tartu, Estonia; and ¶National Institute of Biologicals, Noida 201307, India

To whom correspondence should be addressed. E-mail: vkk2k@hotmail.com .

Edited by Colin Renfrew, University of Cambridge, Cambridge, United Kingdom, and approved November 23, 2005 Received September 5, 2005.

Abstract

Understanding the genetic origins and demographic history of Indian populations is important both for questions concerning the early settlement of Eurasia and more recent events, including the appearance of Indo-Aryan languages and settled agriculture in the subcontinent. Although there is general agreement that Indian caste and tribal populations share a common late Pleistocene maternal ancestry in India, some studies of the Y-chromosome markers have suggested a recent, substantial incursion from Central or West Eurasia. To investigate the origin of paternal lineages of Indian populations, 936 Y chromosomes, representing 32 tribal and 45 caste groups from all four major linguistic groups of India, were analyzed for 38 single-nucleotide polymorphic markers. Phylogeography of the major Y-chromosomal haplogroups in India, genetic distance, and admixture analyses all indicate that the recent external contribution to Dravidian- and Hindi-speaking caste groups has been low. The sharing of some Y-chromosomal haplogroups between Indian and Central Asian populations is most parsimoniously explained by a deep, common ancestry between the two regions, with diffusion of some Indian-specific lineages northward. The Y-chromosomal data consistently suggest a largely South Asian origin for Indian caste communities and therefore argue against any major influx, from regions north and west of India, of people associated either with the development of agriculture or the spread of the Indo-Aryan language family. The dyadic Y-chromosome composition of Tibeto-Burman speakers of India, however, can be attributed to a recent demographic process, which appears to have absorbed and overlain populations who previously spoke Austro-Asiatic languages.

Keywords: agriculture, genetic origins, India, paternal lineages

Archaeological evidence advocates the settlement of India by modern humans, using Middle Palaeolithic tools, during the Late Pleistocene (1–5). The large number of deep-rooting, Indian-specific mtDNA lineages of macro haplogroups M and N, whose presence cannot be explained by a recent introduction from neighboring regions (6), is consistent with the archaeological data. These two lines of evidence suggest that the initial settlement, followed by local differentiation, has left a predominantly Late Pleistocene genetic signature in the maternal heritage of India (7–11). The initial settlement of South Asia, between 40,000 and 70,000 years ago, was most likely over the southern route from Africa because haplogroup M, which is the most frequent mtDNA component in India, is virtually absent in the Near East and Southwest Asia (6, 11–14).

Linguistically, the four main language families spoken in India have strong regional patterns, with the largest group, Indo-European (IE), prevalent in northern India. The second largest, the Dravidian (DR) family, covers the majority of the languages in the south. Most of the IE speakers belong to castes, whereas the majority of the tribal populations (>450) speak languages from the other three families (15). The existence of both IE-speaking tribal groups and DR castes indicates the complexity of historical interactions between Indian populations, and that there is no one-to-one correlation between language and mode of subsistence or social system. Even so, the trend of finding farming (castes) and IE languages grouped together has led some to suggest a major demic diffusion, associated with the spread of agriculture, from West Asia and/or Central Asia to India, (16). The situation in Northeast India is less intensely studied, but the proximity to related language families in East and Southeast Asia suggests possible origins for the Tibeto-Burman (TB) clade of the Sino-Tibetan family outside India. The Mundari group of Austro-Asiatic (AA) languages is currently found almost exclusively in East India. Khasi, a major tribe of Meghalaya, forms a notable exception, being surrounded by TB speakers in the northeast. Other members of the AA family are located in Southeast Asia, but the ultimate source of these languages is currently unresolved.

Several studies have argued that, in contrast to the relative uniformity of mtDNA, the Y chromosomes of Indian populations display relatively small genetic distances to those of West Eurasians (17), linking this finding to hypothetical migrations by Indo-Aryan speakers. Wells et al. (18) highlighted M17 (R1a) as a potential marker for one such event, as it demonstrates decreasing frequencies from Central Asia toward South India. Departing from the "one haplogroup equals one migration" scenario, Cordaux et al. (19) defined, heuristically, a package of haplogroups (J2, R1a, R2, and L) to be associated with the migration of IE people and the introduction of the caste system to India, again from Central Asia, because they had been observed at significantly lower proportions in South Indian tribal groups, with the high frequency of R1a among Chenchus of Andhra Pradesh (6) considered as an aberrant phenomenon (19). Conversely, haplogroups H, F*, and O2a, which were observed at significantly higher proportions among tribal groups of South India, led the same authors to single them out as having an indigenous Indian origin. Only O3e was envisaged as originating (recently) east of India (20), substantiating a linguistic correlation with the TB speakers of Southeast Asia.

The present study significantly increases the available sample size for India by typing 936 individuals from 77 populations, representing all four major linguistic groups (Fig. 1). The increased range of informative SNPs typed permits more detailed resolution of geographic patterns and the identification of some region-specific subsets of lineages. These Y chromosomes are analyzed in the context of available data from West Asia, East Asia, Southeast Asia, Central Asia, Europe, the Near East, and Ethiopia. Measures of genetic distance, admixture, and factor analysis drawn from the Y-chromosome data are used to investigate three themes central to population genetics in India: demographic links to West and Central Asia, the genetic relationship between castes and tribes, and geographic versus linguistic grouping for the current populations of the Indian subcontinent.

Fig. 1.

Map of India showing sample locations. Regional groupings of populations as used in the text are highlighted in different colors.

Results

A total of 18 haplogroups were detected in 936 Indian Y chromosomes (Fig. 3A, which is published as supporting information on the PNAS web site). Together, haplogroups R1, R2, L, O, H, J2, and C characterize >90% of the Y-chromosomal variation in all socio-linguistic groups of India (Tables 2 and 3, which are published as supporting information on the PNAS web site). Both IE- and DR-speaking populations show a high combined frequency of haplogroups C*, L1, H1, and R2. The total frequency of these four haplogroups outside of India is marginally low. In turn, haplogroups E, I, G, J*, and R1* have a combined frequency of 53% in the Near East among the Turks and 24% in Central Asia, but they are rare or absent in India (0.86% in all populations and almost solely because of R1*). Similarly, haplogroups C3, D, N, and O specific to Central Asian (36%) and Southeast Asian populations (subclades of haplogroup O; 85%) are virtually absent in India (Fig. 3A). Only haplogroups J2 and R1a have interregional frequency patterns west of India with J2 being most common in Afro-Asiatic-speaking (and IE-speaking) populations of the Near East and Middle East, whereas R1a occurs at the highest frequencies in populations of India, East Europe, and Central Asia. The O2a and O3e subclades of haplogroup O in India also have interregional distributions, overlapping with those of Southeast Asia and East Asia.

Principal component analysis (Fig. 3B) investigates the phylogeography of the Y haplogroups with respect to each other, illustrating the associations of haplogroups, irrespective of regional or cultural categories. The first two components account for 75% of the variation observed, and within India delineate R*, R2, F*, and H, within the sphere of L, K, P*, and R1a. Of all of the R lineages, only R1* is separated from this grouping, forming a cluster together with G, I, and J, consistent with their common and widespread distribution throughout (Western) Europe. The O lineages fall out with C* and D (the latter tending to derive from Sino-Tibetan speakers). Once the third and fourth factors are considered, the ambiguity of A, B, and E (typically African in origin) is resolved, and the positions of C3 and N, also non-Indian in their distribution, are delineated to Central Asia.

By considering all haplogroup frequencies simultaneously, an indication of the relatedness between regions is obtained (Table 1). Here, for the sake of comparison only, the categories used by a previous study (19) are retained, but the tribal population is split into two because of the close association identified here between Hg O and tribal groups of the east and northeast of India (O2a represents 77% of AA speakers and 47% of TB speakers), which are combined to form the east and northeast tribes. In contrast to the earlier study (19), the caste populations of "north" and "south" India are not particularly more closely related to each other (average Fst value = 0.07) than they are to the tribal groups (average Fst value = 0.06). The multidimensional scaling plot of these values (Fig. 4, which is published as supporting information on the PNAS web site) demonstrates that the combined data set for the tribal peoples (derived from all regions of India, excluding those of the east and northeast) actually falls midway between those for northern and southern castes, whereas the tribal populations of the east and northeast are confirmed as a separate category. The position of the reduced tribal category, comprising groups from Southern, Northern, and Western India, is suggestive of geographical structuring north to south.

Table 1.

Genetic distances between populations estimated from Y-haplogroup frequencies

This geographical structure is displayed with greater precision by dividing the data set according to the regions of India presented in Fig. 1, except that the Punjab (caste only) is considered as a separate entity because of its isolation relative to the rest of the west (see Tables 2 and 3) and proximity to Central Asia. Considering individual haplogroup frequencies within each of these geographical regions, no consistent pattern (at the 95% level of certainty) was detected in the distribution of the Y haplogroups to distinguish either the castes from the tribes, or DRs from IEs (Fig. 5 B and C, which is published as supporting information on the PNAS web site). Therefore, it is appropriate to consider the distributions at the regional level, omitting Northeast India because of the dominance of haplogroup O there (Fig. 5A). The potential clines centered on North India (R1a), Northwest India (J2), South India (H), and East India (R2), identified in Fig. 5A, are illustrated by the distribution maps (Fig. 2 and Fig. 6, which is published as supporting information on the PNAS web site). These clines display distinct regional concentrations of J2, H, R1a, R2, O3, and O2a, confirming the primarily geographic nature of Y-chromosome frequency distribution in India.

Fig. 2.

Spatial frequency distribution maps of major Y-chromosome haplogroups in South Asia. For India, the data on tribal populations are shown in the inset maps and excluded from the main maps. The data for caste populations are averaged to the level of states (more ...)

Admixture analysis (21) evaluates the potential parental contributions to northwestern castes (Punjab) and southern castes (Table 4, which is published as supporting information on the PNAS web site) and Central Asia (Table 5, which is published as supporting information on the PNAS web site). For South Indian populations this analysis revealed reciprocally high local admixture contributions for both caste and tribal populations (0.91 contribution of tribes to castes, SE 0.1; 0.98 contribution of castes to tribes, SE 0.11) over the contributions from outside of India. It should be stressed that these values do not necessarily reflect actual admixture proportions between the tribes and castes, as the algorithm that is used to estimate the admixture proportions divides the whole genetic composition of a hybrid between given parental populations. Rather, these findings confirm the results obtained above from the Fst analyses, that Southern castes and tribals are very similar to each other in their Y-chromosomal haplogroup compositions, and that their gene pool is significantly related to the castes of Northwest India (Fig. 5A), among whom a South Indian tribal contribution of 0.48 (SE 0.12) was observed. In contrast, the potential contribution from Central Asia to the Indian Y-chromosomal pool is minor. In the case of Northwest India, there is nothing to choose between two opposing scenarios: (i) the flow of Y chromosomes from Central Asia, and (ii) the flow of Y chromosomes in the opposite direction, to Central Asia from Northwest India. Meanwhile, the West Asian contribution to the Indian Y-chromosomal pool was significantly smaller in all three admixture tests.

Discussion

Leaving aside, for the moment, TB and AA speakers, the distributions of Y haplogroups between India and West and Central Asia display a clear patterning. J is the predominant Y-chromosome haplogroup in populations living west of India. The frequency and subgroup variation of J in West Asia, in the context of the complete absence of J1 and most J2 subgroups within the Indian sample, is consistent with an influx of a subset of J2 lineages to India from the Near East, followed by their subsequent diffusion from India's northwest toward the south and east. In contrast, within India, the complete absence of the derived C3 lineages, which represent >95% of haplogroup C variation in Central Asia (22), suggests that Indian C lineages cannot be ascribed to a recent admixture from the north.

Similarly, the proposition that a high frequency of R1a in India is caused by admixture with populations of Central Asian origin is difficult to substantiate, as the proposed source region does not meet the expectation of containing high frequencies of the other components of haplogroup R, with no examples of R* and generally low incidence of R2, which, unlike J2, does not show evidence of a recent diffusion throughout India from the northwest. Second, it is notable that the results from the admix2 program gave relatively high reciprocal admixture (0.3–0.35) proportions for Northwest Indian and Central Asian populations, despite the incompatibility of the respective haplogroup frequency pools; our Northwest Indian sample totally lacks haplogroups C3, DE, J*, I, G, N, and O, which cover almost half of the Central Asian Y chromosomes, whereas the Central Asian sample is poor in haplogroups C*, F*, H, L, and R2 (with a combined frequency of 10%). Hence, the admixture proportions are driven solely by the shared high frequency of R1a. In other words, if the source of R1a variation in India comes from Central Asia, as claimed by Wells et al. (18) and Cordaux et al. (19), then, under a recent gene flow scenario, one would expect to find the other Central Asian-derived NRY haplogroups (C3, DE, J*, I, G, N, O) in Northwest India at similarly elevated frequencies, but that is not the case.

Alternatively, although the simple admixture scenario does not hold, one could nevertheless argue that the other haplogroups were lost during a hypothetical bottleneck (lineage sorting among the early Indo-Aryans arriving to India). But in line with this scenario, one should expect to observe dramatically lower genetic variation among Indian R1a lineages. In fact, the opposite is true: the STR haplotype diversity on the background of R1a in Central Asia (and also in Eastern Europe) has already been shown to be lower than that in India (6). Rather, the high incidence of R1* and R1a throughout Central Asian and East European populations (without R2 and R* in most cases) is more parsimoniously explained by gene flow in the opposite direction, possibly with an early founder effect in South or West Asia. Note that the admixture method reports positive admixture proportions in cases where just one haplogroup is shared between populations (possibly because of shared deep common ancestry), even if other haplogroup frequencies strongly argue against a recent simple admixture scenario.

Even though more than one explanation could exist for genetic differentiation between castes and tribes in India, the Indo-Aryan migration scenario advocated in ref. 19 rested on the suggestion that all Indian caste groups are similar to each other while being significantly different from the tribes. Using a much more representative data set, numerically, geographically, and definitively, it was not possible to confirm any of the purported differentiations between the caste and tribal pools. Although differences could be found to occur within particular regions, between particular caste and tribal groups, consistent and statistically significant variations at the subcontinental scale were not detected. Although it is arguable that assimilation of tribal populations into the caste system could skew distributions in any particular region, it cannot explain the persistence and prevalence of those lineages put forward as being typical of incoming IEs (J2, R1a, R2, and L) among many of those populations who are still designated as tribals [see also the credibility gap in the groupings of Corduax et al. (19) illustrated in the factor analysis, Fig. 2]. Rather, taken together with the evidence from Fst values, the elements discussed so far (i.e., admixture, factor analysis, and frequency distributions) are more parsimoniously explained by a predominantly pre-IE, pre-Neolithic presence in India, for the majority of those Y lineages considered here (R1a, R2, L1), which occur together with strictly Indian-specific haplogroups and paragroups (C*, F*, H) among both caste and tribal groups.

The distribution of R2, with its concentration in Eastern and Southern India, is not consistent with a recent demographic movement from the northwest. Instead, its prevalence among castes in these regions might represent a recent population expansion, perhaps associated with the transition to agriculture, which may have occurred independently in South Asia (23). A pre-Neolithic chronology for the origins of Indian Y chromosomes is also supported by the lack of a clear delineation between DR and IE speakers. Again, although appeals to language change are plausible for explaining the appearance of supposedly tribe-specific Y lineages among incoming IE speakers, it is much harder to conceive of a systematic movement of external Y-chromosome types in the opposite direction, via the uptake of DR languages. The near absence of L lineages within the IE speakers from Bihar (0%), Orissa (0%), and West Bengal (1.5%) further suggests that the current distribution of Y haplogroups in India is associated primarily with geographic rather than linguistic or cultural determinants.

In contrast, the situation with the TB- and AA-speaking populations is rather intriguing and warrants further discussion. The AA groups have a very clear association with O2a Y-chromosome haplogroup, both in India and Southeast Asia (24, 25), whereas the close association between TB groups and the O3e lineage may indicate a second case where a Y haplogroup is linked to a cultural entity. The present-day distribution of haplogroup O argues for a Southeast Asian homeland for the AA speakers of India (Mundari group), in distinct contrast to the suggestions, based on mtDNA, that the Mundari speakers represent the earliest settlers of India (9, 26). Yet, the contemporary distribution within India of Y-chromosomal haplogroup O2a, on one hand, and AA speakers on the other, cautions against simplistic interpretations of either linguistic or genetic correlations. AA languages, besides being concentrated in East India, also appear as outliers in Madhya Pradesh (Central India) and Maharashtra (West India), whereas O2a is present, sporadically, within other linguistic groups in both South and East India.

Among TB speakers the share of mtDNAs typical of East Asia increases to nearly two-thirds (64%), inferred from ref. 27. This scenario would be consistent with a more recent migration event or the continued movement of women into India through the maintenance of social links. The near total absence of AA-speaking groups between East India and Southeast Asia has been interpreted as representing a recent (mid-Holocene) influx of TB populations, bearing O3e Y chromosomes, into this region (20). Cordaux et al. (20), when considering different scenarios for the prehistory of this area, favored the view that it was previously an unoccupied territory that had acted as a barrier to human migrations, possibly since the late Pleistocene. However, the presence of the AA-speaking Khasi in Meghalaya provides an alternative explanation, namely that there were previous inhabitants in this region who had been predominantly AA speakers. This explanation is favored by the presence of both O3e and O2a Y haplogroups within the TB populations reported here. The parsimonious explanation for this is that AA speakers were formerly distributed from Southeast Asia to India and intermixed with TB speakers as they migrated to the area. This scenario is supported by the widespread presence of East Asian mtDNA lineages among TB groups. So, paradoxically, it is in Northeastern India that there is evidence, from the Y chromosome, for both large-scale immigration (TB speakers) and language change (former AA speakers). One of the reasons this is still detectable is the relatively shallow time depth proposed for this "event," a chronology that still covers the period proposed for the appearance of the caste system, the IE language family, and agriculture into India through the northwest (20).

Conclusions

It is not necessary, based on the current evidence, to look beyond South Asia for the origins of the paternal heritage of the majority of Indians at the time of the onset of settled agriculture. The perennial concept of people, language, and agriculture arriving to India together through the northwest corridor does not hold up to close scrutiny. Recent claims for a linkage of haplogroups J2, L, R1a, and R2 with a contemporaneous origin for the majority of the Indian castes' paternal lineages from outside the subcontinent are rejected, although our findings do support a local origin of haplogroups F* and H. Of the others, only J2 indicates an unambiguous recent external contribution, from West Asia rather than Central Asia. The current distributions of haplogroup frequencies are, with the exception of the O lineages, predominantly driven by geographical, rather than cultural determinants. Ironically, it is in the northeast of India, among the TB groups that there is clear-cut evidence for large-scale demic diffusion traceable by genes, culture, and language, but apparently not by agriculture.