Origin and identity of the Brokpa of Dah-Hanu, Himalayas - an NRY-HG L1a2 (M357) legacy.

Background: The Brokpas are an isolated tribal population of the Dah-Hanu villages of the Leh district of India. They speak Dardic, a sub-branch of the Indo-European language family, and are putatively identified as "pure Aryan," a hegemonic impression perpetuated by foreign tourism.Aim: To determine if the above is true by looking for an appreciable frequency of NRY-HG-R1a1(M17) signatures which are common to Indo-European language speakers of mainland India and elsewhere.Subjects and methods: We studied 75 random Brokpa males from the Dah-Hanu region, on the northern bank of the Indus river.Results: Interestingly, the Brokpa males possessed a high proportion of NRY-HG-L1a2(M357) (62.7%) that are found sporadically in India and her neighbourhood. A global analysis of this clade (present study, 214 of 3327 men from 63 populations; from the literature 56 of 873) suggested that they originated from southern India.Conclusion: The Y chromosomal studies suggest the Brokpa to be pre-Vedic settlers of the Himalayas, 9000 ybp, with an isolated evolution. The mtDNA profile shows a predominance of mtDNA HG A4 that must have arrived from outside the Indian subcontinent.

The light skin allele of SLC24A5 in South Asians and Europeans shares identity by descent.

Skin pigmentation is one of the most variable phenotypic traits in humans. A non-synonymous substitution (rs1426654) in the third exon of SLC24A5 accounts for lighter skin in Europeans but not in East Asians. A previous genome-wide association study carried out in a heterogeneous sample of UK immigrants of South Asian descent suggested that this gene also contributes significantly to skin pigmentation variation among South Asians. In the present study, we have quantitatively assessed skin pigmentation for a largely homogeneous cohort of 1228 individuals from the Southern region of the Indian subcontinent. Our data confirm significant association of rs1426654 SNP with skin pigmentation, explaining about 27% of total phenotypic variation in the cohort studied. Our extensive survey of the polymorphism in 1573 individuals from 54 ethnic populations across the Indian subcontinent reveals wide presence of the derived-A allele, although the frequencies vary substantially among populations. We also show that the geospatial pattern of this allele is complex, but most importantly, reflects strong influence of language, geography and demographic history of the populations. Sequencing 11.74 kb of SLC24A5 in 95 individuals worldwide reveals that the rs1426654-A alleles in South Asian and West Eurasian populations are monophyletic and occur on the background of a common haplotype that is characterized by low genetic diversity. We date the coalescence of the light skin associated allele at 22-28 KYA. Both our sequence and genome-wide genotype data confirm that this gene has been a target for positive selection among Europeans. However, the latter also shows additional evidence of selection in populations of the Middle East, Central Asia, Pakistan and North India but not in South India.

Unusual selection on the KIR3DL1/S1 natural killer cell receptor in Africans.

Interactions of killer cell immunoglobulin-like receptors (KIRs) with major histocompatibility complex (MHC) class I ligands diversify natural killer cell responses to infection. By analyzing sequence variation in diverse human populations, we show that the KIR3DL1/S1 locus encodes two lineages of polymorphic inhibitory KIR3DL1 allotypes that recognize Bw4 epitopes of protein">HLA-A and HLA-B and one lineage of conserved activating KIR3DS1 allotypes, also implicated in Bw4 recognition. Balancing selection has maintained these three lineages for over 3 million years. Variation was selected at D1 and D2 domain residues that contact HLA class I and at two sites on D0, the domain that enhances the binding of KIR3D to HLA class I. HLA-B variants that gained Bw4 through interallelic microconversion are also products of selection. A worldwide comparison uncovers unusual KIR3DL1/S1 evolution in modern sub-Saharan Africans. Balancing selection is weak and confined to D0, KIR3DS1 is rare and KIR3DL1 allotypes with similar binding sites predominate. Natural killer cells express the dominant KIR3DL1 at a high frequency and with high surface density, providing strong responses to cells perturbed in Bw4 expression.

Genome-wide signatures of male-mediated migration shaping the Indian gene pool.

Multiple questions relating to contributions of cultural and demographical factors in the process of human geographical dispersal remain largely unanswered. India, a land of early human settlement and the resulting diversity is a good place to look for some of the answers. In this study, we explored the genetic structure of India using a diverse panel of 78 males genotyped using the GenoChip. Their genome-wide single-nucleotide polymorphism (SNP) diversity was examined in the context of various covariates that influence Indian gene pool. Admixture analysis of genome-wide SNP data showed high proportion of the Southwest Asian component in all of the Indian samples. Hierarchical clustering based on admixture proportions revealed seven distinct clusters correlating to geographical and linguistic affiliations. Convex hull overlay of Y-chromosomal haplogroups on the genome-wide SNP principal component analysis brought out distinct non-overlapping polygons of F*-M89, H*-M69, L1-M27, O2a-M95 and O3a3c1-M117, suggesting a male-mediated migration and expansion of the Indian gene pool. Lack of similar correlation with mitochondrial DNA clades indicated a shared genetic ancestry of females. We suggest that ancient male-mediated migratory events and settlement in various regional niches led to the present day scenario and peopling of India.

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.

Herders of Indian and European cattle share their predominant allele for lactase persistence.

Milk consumption and lactose digestion after weaning are exclusively human traits made possible by the continued production of the enzyme lactase in adulthood. Multiple independent mutations in a 100-bp region--part of an enhancer--approximately 14-kb upstream of the LCT gene are associated with this trait in Europeans and pastoralists from Saudi Arabia and Africa. However, a single mutation of purported western Eurasian origin accounts for much of observed lactase persistence outside Africa. Given the high levels of present-day milk consumption in India, together with archaeological and genetic evidence for the independent domestication of cattle in the Indus valley roughly 7,000 years ago, we sought to determine whether lactase persistence has evolved independently in the subcontinent. Here, we present the results of the first comprehensive survey of the LCT enhancer region in south Asia. Having genotyped 2,284 DNA samples from across the Indian subcontinent, we find that the previously described west Eurasian -13910 C>T mutation accounts for nearly all the genetic variation we observed in the 400- to 700-bp LCT regulatory region that we sequenced. Geography is a significant predictor of -13910*T allele frequency, and consistent with other genomic loci, its distribution in India follows a general northwest to southeast declining pattern, although frequencies among certain neighboring populations vary substantially. We confirm that the mutation is identical by descent to the European allele and is associated with the same>1 Mb extended haplotype in both populations.

Recombination networks as genetic markers in a human variation study of the Old World.

We have analyzed human genetic diversity in 33 Old World populations including 23 populations obtained through Genographic Project studies. A set of 1,536 SNPs in five X chromosome regions were genotyped in 1,288 individuals (mostly males). We use a novel analysis employing subARG network construction with recombining chromosomal segments. Here, a subARG is constructed independently for each of five gene-free regions across the X chromosome, and the results are aggregated across them. For PCA, MDS and ancestry inference with STRUCTURE, the subARG is processed to obtain feature vectors of samples and pairwise distances between samples. The observed population structure, estimated from the five short X chromosomal segments, supports genome-wide frequency-based analyses: African populations show higher genetic diversity, and the general trend of shared variation is seen across the globe from Africa through Middle East, Europe, Central Asia, Southeast Asia, and East Asia in broad patterns. The recombinational analysis was also compared with established methods based on SNPs and haplotypes. For haplotypes, we also employed a fixed-length approach based on information-content optimization. Our recombinational analysis suggested a southern migration route out of Africa, and it also supports a single, rapid human expansion from Africa to East Asia through South Asia.

Leprosy and the adaptation of human toll-like receptor 1.

Leprosy is an infectious disease caused by the obligate intracellular pathogen Mycobacterium leprae and remains endemic in many parts of the world. Despite several major studies on susceptibility to leprosy, few genomic loci have been replicated independently. We have conducted an association analysis of more than 1,500 individuals from different case-control and family studies, and observed consistent associations between genetic variants in both TLR1 and the HLA-DRB1/DQA1 regions with susceptibility to leprosy (TLR1 I602S, case-control P = 5.7 x 10(-8), OR = 0.31, 95% CI = 0.20-0.48, and HLA-DQA1 rs1071630, case-control P = 4.9 x 10(-14), OR = 0.43, 95% CI = 0.35-0.54). The effect sizes of these associations suggest that TLR1 and HLA-DRB1/DQA1 are major susceptibility genes in susceptibility to leprosy. Further population differentiation analysis shows that the TLR1 locus is extremely differentiated. The protective dysfunctional 602S allele is rare in Africa but expands to become the dominant allele among individuals of European descent. This supports the hypothesis that this locus may be under selection from mycobacteria or other pathogens that are recognized by TLR1 and its co-receptors. These observations provide insight into the long standing host-pathogen relationship between human and mycobacteria and highlight the key role of the TLR pathway in infectious diseases.

HLA-B57 and gender influence the occurrence of tuberculosis in HIV infected people of south India.

BACKGROUND: Substantial evidence exists for HLA and other host genetic factors being determinants of susceptibility or resistance to infectious diseases. However, very little information is available on the role of host genetic factors in HIV-TB coinfection. Hence, a longitudinal study was undertaken to investigate HLA associations in a cohort of HIV seropositive individuals with and without TB in Bangalore, South India. METHODS: A cohort of 238 HIV seropositive subjects were typed for HLA-A, B, and DR by PCR-SSP and followed up for 5 years or till manifestation of Tuberculosis. HLA data of 682 HIV Negative healthy renal donors was used as control. RESULTS: The ratio of males and females in HIV cohort was comparable (50.4% and 49.6%). But the incidence of TB was markedly lower in females (12.6%,) than males (25.6%). Further, HLA-B*57 frequency in HIV cohort was significantly higher among females without TB (21.6%, 19/88) than males (1.7%, 1/59); P = 0.0046; OR = 38. CD4 counts also were higher among females in this cohort. CONCLUSION: This study suggests that HIV positive women with HLA-B*57 have less occurrence of TB as compared to males.

Not all the infected develop the disease - A "Lotus and Cactus" model.

The immunogenetic dictum "not all the infected develop the disease" can best be explained by a "Lotus and Cactus" model. Lotuses grow in ponds and cacti in deserts: analogously, we can say that tubercle patient's lung (genetic makeup) functions as an ideal 'broth' for Mycobacterium tuberculosis (M.tb) germs to grow, but not the lungs of an endemic control. HLA association studies from Europe to Asia since 1983 till date, have shown a persistent HLA DR2 (15) association. Further, HLA DR2 and non-DR2 endemic controls showed disparate patterns of immune responses and gene expressions. The host and pathogen MHC diversities, Th1-Th2 paradigm and cytokine circuits all may play a crucial role in TB susceptibility. It is possible to decipher the protective immunity by controlling the known confounders - epidemiological, demographic, socio-biological and also host and pathogen diversities. This has become significant with our understanding on the 'out of Africa' migration and neolithic co-dispersal of M.tb with modern human. Divergence and expansion of various MHCs (eg HLA-DRB1*15, HLA-B*57) and non-MHC alleles in various continents might be responsible for the skewed transmission and distribution of the infectious diseases around the globe. The 'Lotus and Cactus' model proposed here exemplifies this. A holistic genetic epidemiology approach employing modern tools is the need of the hour to better understand infectious disease susceptibility.

Human leukocyte antigens in hypertrophic cardiomyopathy patients in South India.

Hypertrophic cardiomyopathy is characterized by massive ventricular hypertrophy, reduced diastolic function, and excessive ventricular contraction. The human leukocyte antigens HLA-A, HLA-B, and HLA-DR were studied in 14 hypertrophic cardiomyopathy patients with left ventricular obstruction from South India. They were compared with 81 normal age- and sex-matched individuals from the same ethnic background. The human leucocyte antigens were identified using the standard serological assay with a longer incubation for DR antigens. The odds ratio, frequency, chi-squared value, p-value, etiological fraction, preventive fraction, and haplotype frequency estimates were calculated. The HLA-B51 and HLA-DR2 levels were significantly increased in hypertrophic cardiomyopathy patients compared to controls, whereas HLA-A19, HLA-B7, and HLA-DR4 were decreased when compared to the controls. It was noticed that haplotype B51-DR2-DQ3 was significantly associated with hypertrophic cardiomyopathy patients from South India. Hypertrophic cardiomyopathy may be associated with genes in the human leukocyte antigen region, and immunogenetic factors linked to human leukocyte antigens appear to play a major role in the pathogenesis.

Geographic population structure analysis of worldwide human populations infers their biogeographical origins.

The search for a method that utilizes biological information to predict humans' place of origin has occupied scientists for millennia. Over the past four decades, scientists have employed genetic data in an effort to achieve this goal but with limited success. While biogeographical algorithms using next-generation sequencing data have achieved an accuracy of 700 km in Europe, they were inaccurate elsewhere. Here we describe the Geographic Population Structure (GPS) algorithm and demonstrate its accuracy with three data sets using 40,000-130,000 SNPs. GPS placed 83% of worldwide individuals in their country of origin. Applied to over 200 Sardinians villagers, GPS placed a quarter of them in their villages and most of the rest within 50 km of their villages. GPS's accuracy and power to infer the biogeography of worldwide individuals down to their country or, in some cases, village, of origin, underscores the promise of admixture-based methods for biogeography and has ramifications for genetic ancestry testing.

Population differentiation of southern Indian male lineages correlates with agricultural expansions predating the caste system.

Previous studies that pooled Indian populations from a wide variety of geographical locations, have obtained contradictory conclusions about the processes of the establishment of the Varna caste system and its genetic impact on the origins and demographic histories of Indian populations. To further investigate these questions we took advantage that both Y chromosome and caste designation are paternally inherited, and genotyped 1,680 Y chromosomes representing 12 tribal and 19 non-tribal (caste) endogamous populations from the predominantly Dravidian-speaking Tamil Nadu state in the southernmost part of India. Tribes and castes were both characterized by an overwhelming proportion of putatively Indian autochthonous Y-chromosomal haplogroups (H-M69, F-M89, R1a1-M17, L1-M27, R2-M124, and C5-M356; 81% combined) with a shared genetic heritage dating back to the late Pleistocene (10-30 Kya), suggesting that more recent Holocene migrations from western Eurasia contributed <20% of the male lineages. We found strong evidence for genetic structure, associated primarily with the current mode of subsistence. Coalescence analysis suggested that the social stratification was established 4-6 Kya and there was little admixture during the last 3 Kya, implying a minimal genetic impact of the Varna (caste) system from the historically-documented Brahmin migrations into the area. In contrast, the overall Y-chromosomal patterns, the time depth of population diversifications and the period of differentiation were best explained by the emergence of agricultural technology in South Asia. These results highlight the utility of detailed local genetic studies within India, without prior assumptions about the importance of Varna rank status for population grouping, to obtain new insights into the relative influences of past demographic events for the population structure of the whole of modern India.

Meiotic recombination generates rich diversity in NK cell receptor genes, alleles, and haplotypes.

Natural killer (NK) cells contribute to the essential functions of innate immunity and reproduction. Various genes encode NK cell receptors that recognize the major histocompatibility complex (MHC) Class I molecules expressed by other cells. For primate NK cells, the killer-cell immunoglobulin-like receptors (KIR) are a variable and rapidly evolving family of MHC Class I receptors. Studied here is KIR3DL1/S1, which encodes receptors for highly polymorphic human HLA-A and -B and comprises three ancient allelic lineages that have been preserved by balancing selection throughout human evolution. While the 3DS1 lineage of activating receptors has been conserved, the two 3DL1 lineages of inhibitory receptors were diversified through inter-lineage recombination with each other and with 3DS1. Prominent targets for recombination were D0-domain polymorphisms, which modulate enhancer function, and dimorphism at position 283 in the D2 domain, which influences inhibitory function. In African populations, unequal crossing over between the 3DL1 and 3DL2 genes produced a deleted KIR haplotype in which the telomeric "half" was reduced to a single fusion gene with functional properties distinct from its 3DL1 and 3DL2 parents. Conversely, in Eurasian populations, duplication of the KIR3DL1/S1 locus by unequal crossing over has enabled individuals to carry and express alleles of all three KIR3DL1/S1 lineages. These results demonstrate how meiotic recombination combines with an ancient, preserved diversity to create new KIR phenotypes upon which natural selection acts. A consequence of such recombination is to blur the distinction between alleles and loci in the rapidly evolving human KIR gene family.

Distinct diversity of KIR genes in three southern Indian populations: comparison with world populations revealed a link between KIR gene content and pre-historic human migrations.

By interacting with polymorphic HLA class I molecules, the killer cell immunoglobulin-like receptors (KIR) influence the innate and adaptive immune response to infection. The KIR family varies in gene content and sequence polymorphism, thereby, distinguishing individuals and populations. To investigate KIR diversity in the earliest settlers of India, we have characterized the KIR gene content in three Dravidian-speaking populations (Mollukurumba, Kanikar, and Paravar) from the state of Tamil Nadu, southern India. The activating KIR genes and putative group-B KIR haplotypes were frequent in Paravar and Kanikar, a scenario analogous to those seen previously in other populations of Indian origin, indicating that predominance of group-B KIR haplotypes is the characteristic feature of Indian populations. In contrast, the KIR gene profile of Mollukurumba was more related to Caucasian type. It is not clear whether a local-specific selection or a recent admixture from Iran is responsible for such discrete profile in Mollukurumba. Each southern Indian population had distinct KIR genotype profile. Comparative analyses with world populations revealed that group-B KIR haplotypes were frequent in the natives of India, Australia, and America, the populations associated with those involved in extensive prehistoric human migrations. Whether or not natural selection has acted to enrich group-B KIR haplotypes in these migratory descendants is an issue that requires objective testing.

Variation in MICA and MICB genes and enhanced susceptibility to paucibacillary leprosy in South India.

In a study of mainly paucibacillary leprosy-affected sib-pair families from South India, in addition to the expected associations with the HLA-DRB1 locus, we have identified significant association with a functional variant of the MICA gene as well as a microsatellite in the flanking region of the MICB gene. The associations with MICA and MICB cannot be accounted for by linkage disequilibrium with the HLA class II locus indicating a role in genetic susceptibility to leprosy that is independent of HLA-DRB1. Previous studies have shown that MICA and MICB are expressed on the surface of cells in response to infection, where they are recognized by the NKG2D receptor on gammadelta T cells, CD8+ alphabeta T cells and natural killer cells, all of which contribute to defense against mycobacteria. The MICA*5A5.1 allele, associated here with leprosy susceptibility, encodes a protein lacking a cytoplasmic tail providing a possible mechanism for defective immune surveillance against mycobacteria.

A region of chromosome 20 is linked to leprosy susceptibility in a South Indian population.

A major susceptibility locus for leprosy has recently been mapped on chromosome 10 (10p13) by genome-wide linkage analysis. Microsatellite markers from this genome screen that showed suggestive evidence of linkage to leprosy were evaluated in an additional 140 families with affected sib pairs. A second region of linkage has thus been identified on chromosome 20 (20p12). The peak of linkage lies at marker D20S115, which has a significant single-point maximum logarithm of odds score of 3.48 (P=.00003). Transmission disequilibrium testing of the microsatellite markers in 20p12 showed that the marker D20S835 is associated with protection against leprosy (P=.021), which suggests that a locus controlling susceptibility lies close to this marker.