Breuss, Martin | Heng, Julian Ik-Tsen | Poirier, Karine | Tian, Guoling | Jaglin, Xavier Hubert | Qu, Zhengdong | Braun, Andreas | Gstrein, Thomas | Ngo, Linh | Haas, Matilda | Bahi-Buisson, Nadia | Moutard, Marie-Laure | Passemard, Sandrine | Verloes, Alain | Gressens, Pierre | Xie, Yunli | Robson, Kathryn J.H. | Rani, Deepa Selvi | Thangaraj, Kumarasamy | Clausen, Tim | Chelly, Jamel | Cowan, Nicholas Justin | Keays, David Anthony
SUMMARY
The formation of the mammalian cortex requires the generation, migration, and differentiation of neurons. The vital role that the microtubule cytoskeleton plays in these cellular processes is reflected by the discovery that mutations in various tubulin isotypes cause different neurodevelopmental diseases, including lissencephaly (TUBA1A), polymicrogyria (TUBA1A, TUBB2B, TUBB3), and an ocular motility disorder (TUBB3). Here, we show that Tubb5 is expressed in neurogenic progenitors in the mouse and that its depletion in vivo perturbs the cell cycle of progenitors and alters the position of migrating neurons. We report the occurrence of three microcephalic patients with structural brain abnormalities harboring de novo mutations in TUBB5 (M299V, V353I, and E401K). These mutant proteins, which affect the chaperone-dependent assembly of tubulin heterodimers in different ways, disrupt neurogenic division and/or migration in vivo. Our results provide insight into the functional repertoire of the tubulin gene family, specifically implicating TUBB5 in embryonic neurogenesis and microcephaly.
doi:10.1016/j.celrep.2012.11.017
PMCID: PMC3595605
PMID: 23246003
Rai, Niraj | Chaubey, Gyaneshwer | Tamang, Rakesh | Pathak, Ajai Kumar | Singh, Vipin Kumar | Karmin, Monika | Singh, Manvendra | Rani, Deepa Selvi | Anugula, Sharath | Yadav, Brijesh Kumar | Singh, Ashish | Srinivasagan, Ramkumar | Yadav, Anita | Kashyap, Manju | Narvariya, Sapna | Reddy, Alla G. | van Driem, George | Underhill, Peter A. | Villems, Richard | Kivisild, Toomas | Singh, Lalji | Thangaraj, Kumarasamy | Yao, Yong-Gang
Linguistic and genetic studies on Roma populations inhabited in Europe have unequivocally traced these populations to the Indian subcontinent. However, the exact parental population group and time of the out-of-India dispersal have remained disputed. In the absence of archaeological records and with only scanty historical documentation of the Roma, comparative linguistic studies were the first to identify their Indian origin. Recently, molecular studies on the basis of disease-causing mutations and haploid DNA markers (i.e. mtDNA and Y-chromosome) supported the linguistic view. The presence of Indian-specific Y-chromosome haplogroup H1a1a-M82 and mtDNA haplogroups M5a1, M18 and M35b among Roma has corroborated that their South Asian origins and later admixture with Near Eastern and European populations. However, previous studies have left unanswered questions about the exact parental population groups in South Asia. Here we present a detailed phylogeographical study of Y-chromosomal haplogroup H1a1a-M82 in a data set of more than 10,000 global samples to discern a more precise ancestral source of European Romani populations. The phylogeographical patterns and diversity estimates indicate an early origin of this haplogroup in the Indian subcontinent and its further expansion to other regions. Tellingly, the short tandem repeat (STR) based network of H1a1a-M82 lineages displayed the closest connection of Romani haplotypes with the traditional scheduled caste and scheduled tribe population groups of northwestern India.
doi:10.1371/journal.pone.0048477
PMCID: PMC3509117
PMID: 23209554
Jha, Aditya Nath | Singh, Vipin Kumar | Kumari, Namrata | Singh, Ashish | Antony, Justin | van Tong, Hoang | Singh, Sakshi | Pati, Sudhanshu S. | Patra, Pradeep K. | Singh, Rajender | Toan, Nguyen L. | Song, Le H. | Assaf, Amal | Messias–Reason, Iara J. T. | Velavan, Thirumalaisamy P. | Singh, Lalji | Thangaraj, Kumarasamy | Ahmed, Niyaz
Background
Interleukin 4 (IL-4) is an anti-inflammatory cytokine, which regulates balance between TH1 and TH2 immune response, immunoglobulin class switching and humoral immunity. Polymorphisms in this gene have been reported to affect the risk of infectious and autoimmune diseases.
Methods
We have analyzed three regulatory IL-4 polymorphisms; -590C>T, -34C>T and 70 bp intron-3 VNTR, in 4216 individuals; including: (1) 430 ethnically matched case-control groups (173 severe malaria, 101 mild malaria and 156 asymptomatic); (2) 3452 individuals from 76 linguistically and geographically distinct endogamous populations of India, and (3) 334 individuals with different ancestry from outside India (84 Brazilian, 104 Syrian, and 146 Vietnamese).
Results
The -590T, -34T and intron-3 VNTR R2 alleles were found to be associated with reduced malaria risk (P<0.001 for -590C>T and -34C>T, and P = 0.003 for VNTR). These three alleles were in strong LD (r2>0.75) and the TTR2 (-590T, -34T and intron-3 VNTR R2) haplotype appeared to be a susceptibility factor for malaria (P = 0.009, OR = 0.552, 95% CI = 0.356 –0.854). Allele and genotype frequencies differ significantly between caste, nomadic, tribe and ancestral tribal populations (ATP). The distribution of protective haplotype TTR2 was found to be significant (χ23 = 182.95, p-value <0.001), which is highest in ATP (40.5%); intermediate in tribes (33%); and lowest in caste (17.8%) and nomadic (21.6%).
Conclusions
Our study suggests that the IL-4 polymorphisms regulate host susceptibility to malaria and disease progression. TTR2 haplotype, which gives protection against malaria, is high among ATPs. Since they inhabited in isolation and mainly practice hunter-gatherer lifestyles and exposed to various parasites, IL-4 TTR2 haplotype might be under positive selection.
doi:10.1371/journal.pone.0048136
PMCID: PMC3480467
PMID: 23110190
India is known for its vast human diversity, consisting of more than four and a half thousand anthropologically well-defined populations. Each population differs in terms of language, culture, physical features and, most importantly, genetic architecture. The size of populations varies from a few hundred to millions. Based on the social structure, Indians are classified into various caste, tribe and religious groups. These social classifications are very rigid and have remained undisturbed by emerging urbanisation and cultural changes. The variable social customs, strict endogamy marriage practices, long-term isolation and evolutionary forces have added immensely to the diversification of the Indian populations. These factors have also led to these populations acquiring a set of Indian-specific genetic variations responsible for various diseases in India. Interestingly, most of these variations are absent outside the Indian subcontinent. Thus, this review is focused on the peopling of India, the caste system, marriage practice and the resulting health and forensic implications.
doi:10.1186/2041-2223-3-20
PMCID: PMC3514343
PMID: 23020857
Admixture; caste; Indians; mtDNA; tribe; Y-chromosome
Background
Troponin I (TNNI3) is the inhibitory subunit of the thin filament regulatory complex Troponin, which confers calcium-sensitivity to striated muscle actomyosin ATPase activity. Mutations (2-7%) in this gene had been reported in hypertrophic cardiomyopathy patients (HCM). However, the frequencies of mutations and associated clinical presentation have not been established in cardiomyopathy patients of Indian origin, hence we have undertaken this study.
Methods
We have sequenced all the exons, including the exon-intron boundaries of TNNI3 gene in 101 hypertrophic cardiomyopathy patients (HCM), along with 160 healthy controls, inhabited in the same geographical region of southern India.
Results
Our study revealed a total of 16 mutations. Interestingly, we have observed Arginine to Glutamine (R to Q) mutation at 3 positions 98, 141 and 162, exclusively in HCM patients with family history of sudden cardiac death. The novel R98Q was observed in a severe hypertrophic obstructive cardiomyopathy patient (HOCM). The R141Q mutation was observed in two familial cases of severe asymmetric septal hypertrophy (ASH++). The R162Q mutation was observed in a ASH++ patient with mean septal thickness of 29 mm, and have also consists of allelic heterogeneity by means of having one more synonymous (E179E) mutation at g.4797: G → A: in the same exon 7, which replaces a very frequent codon (GAG: 85%) with a rare codon (GAA: 14%). Screening for R162Q mutation in all the available family members revealed its presence in 9 individuals, including 7 with allelic heterogeneity (R162Q and E179E) of which 4 were severely affected. We also found 2 novel SNPs, (g.2653; G → A and g.4003 C → T) exclusively in HCM, and in silico analysis of these SNPs have predicted to cause defect in recognition/binding sites for proteins responsible for proper splicing.
Conclusion
Our study has provided valuable information regarding the prevalence of TNNI3 mutations in Indian HCM patients and its risk assessment, these will help in genetic counseling and to adopt appropriate treatment strategies.
doi:10.1186/1471-2350-13-69
PMCID: PMC3495047
PMID: 22876777
TNNI3-Troponin I; Cardiomyopathy; SNPs; HCM; Indians; Mutations
Chaubey, Gyaneshwer | Metspalu, Mait | Choi, Ying | Mägi, Reedik | Romero, Irene Gallego | Soares, Pedro | van Oven, Mannis | Behar, Doron M. | Rootsi, Siiri | Hudjashov, Georgi | Mallick, Chandana Basu | Karmin, Monika | Nelis, Mari | Parik, Jüri | Reddy, Alla Goverdhana | Metspalu, Ene | van Driem, George | Xue, Yali | Tyler-Smith, Chris | Thangaraj, Kumarasamy | Singh, Lalji | Remm, Maido | Richards, Martin B. | Lahr, Marta Mirazon | Kayser, Manfred | Villems, Richard | Kivisild, Toomas
The geographic origin and time of dispersal of Austroasiatic (AA) speakers, presently settled in south and southeast Asia, remains disputed. Two rival hypotheses, both assuming a demic component to the language dispersal, have been proposed. The first of these places the origin of Austroasiatic speakers in southeast Asia with a later dispersal to south Asia during the Neolithic, whereas the second hypothesis advocates pre-Neolithic origins and dispersal of this language family from south Asia. To test the two alternative models, this study combines the analysis of uniparentally inherited markers with 610,000 common single nucleotide polymorphism loci from the nuclear genome. Indian AA speakers have high frequencies of Y chromosome haplogroup O2a; our results show that this haplogroup has significantly higher diversity and coalescent time (17–28 thousand years ago) in southeast Asia, strongly supporting the first of the two hypotheses. Nevertheless, the results of principal component and “structure-like” analyses on autosomal loci also show that the population history of AA speakers in India is more complex, being characterized by two ancestral components—one represented in the pattern of Y chromosomal and EDAR results and the other by mitochondrial DNA diversity and genomic structure. We propose that AA speakers in India today are derived from dispersal from southeast Asia, followed by extensive sex-specific admixture with local Indian populations.
doi:10.1093/molbev/msq288
PMCID: PMC3355372
PMID: 20978040
Austroasiatic; mtDNA; Y chromosome; autosomes; admixture
Sharma, Gunjan | Tamang, Rakesh | Chaudhary, Ruchira | Singh, Vipin Kumar | Shah, Anish M. | Anugula, Sharath | Rani, Deepa Selvi | Reddy, Alla G. | Eaaswarkhanth, Muthukrishnan | Chaubey, Gyaneshwer | Singh, Lalji | Thangaraj, Kumarasamy | Kivisild, Toomas
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.
doi:10.1371/journal.pone.0032546
PMCID: PMC3290590
PMID: 22393414
Mehrotra, Sanjana | Fakiola, Michaela | Oommen, Joyce | Jamieson, Sarra E | Mishra, Anshuman | Sudarshan, Medhavi | Tiwary, Puja | Rani, Deepa Selvi | Thangaraj, Kumarasamy | Rai, Madhukar | Sundar, Shyam | Blackwell, Jenefer M
Background
IL8RA and IL8RB, encoded by CXCR1 and CXCR2, are receptors for interleukin (IL)-8 and other CXC chemokines involved in chemotaxis and activation of polymorphonuclear neutrophils (PMN). Variants at CXCR1 and CXCR2 have been associated with susceptibility to cutaneous and mucocutaneous leishmaniasis in Brazil. Here we investigate the role of CXCR1/CXCR2 in visceral leishmaniasis (VL) in India.
Methods
Three single nucleotide polymorphisms (SNPs) (rs4674259, rs2234671, rs3138060) that tag linkage disequilibrium blocks across CXCR1/CXCR2 were genotyped in primary family-based (313 cases; 176 nuclear families; 836 individuals) and replication (941 cases; 992 controls) samples. Family- and population-based analyses were performed to look for association between CXCR1/CXCR2 variants and VL. Quantitative RT/PCR was used to compare CXCR1/CXCR2 expression in mRNA from paired splenic aspirates taken before and after treatment from 19 VL patients.
Results
Family-based analysis using FBAT showed association between VL and SNPs CXCR1_rs2234671 (Z-score = 2.935, P = 0.003) and CXCR1_rs3138060 (Z-score = 2.22, P = 0.026), but not with CXCR2_rs4674259. Logistic regression analysis of the case-control data under an additive model of inheritance showed association between VL and SNPs CXCR2_rs4674259 (OR = 1.15, 95%CI = 1.01-1.31, P = 0.027) and CXCR1_rs3138060 (OR = 1.25, 95%CI = 1.02-1.53, P = 0.028), but not with CXCR1_rs2234671. The 3-locus haplotype T_G_C across these SNPs was shown to be the risk haplotype in both family- (TRANSMIT; P = 0.014) and population- (OR = 1.16, P = 0.028) samples (combined P = 0.002). CXCR2, but not CXCR1, expression was down regulated in pre-treatment compared to post-treatment splenic aspirates (P = 0.021).
Conclusions
This well-powered primary and replication genetic study, together with functional analysis of gene expression, implicate CXCR2 in determining outcome of VL in India.
doi:10.1186/1471-2350-12-162
PMCID: PMC3260103
PMID: 22171941
Ayyasamy, Vanniarajan | Owens, Kjerstin M. | Desouki, Mohamed Mokhtar | Liang, Ping | Bakin, Andrei | Thangaraj, Kumarasamy | Buchsbaum, Donald J. | LoBuglio, Albert F. | Singh, Keshav K. | Lewin, Alfred
The Warburg Effect is characterized by an irreversible injury to mitochondrial oxidative phosphorylation (OXPHOS) and an increased rate of aerobic glycolysis. In this study, we utilized a breast epithelial cell line lacking mitochondrial DNA (rho0) that exhibits the Warburg Effect associated with breast cancer. We developed a MitoExpress array for rapid analysis of all known nuclear genes encoding the mitochondrial proteome. The gene-expression pattern was compared among a normal breast epithelial cell line, its rho0 derivative, breast cancer cell lines and primary breast tumors. Among several genes, our study revealed that over-expression of mitochondrial uncoupling protein UCP2 in rho0 breast epithelial cells reflects gene expression changes in breast cancer cell lines and in primary breast tumors. Furthermore, over-expression of UCP2 was also found in leukemia, ovarian, bladder, esophagus, testicular, colorectal, kidney, pancreatic, lung and prostate tumors. Ectopic expression of UCP2 in MCF7 breast cancer cells led to a decreased mitochondrial membrane potential and increased tumorigenic properties as measured by cell migration, in vitro invasion and anchorage independent growth. Consistent with in vitro studies, we demonstrate that UCP2 over-expression leads to development of tumors in vivo in an orthotopic model of breast cancer. Genipin, a plant derived small molecule, suppressed the UCP2 led tumorigenic properties, which were mediated by decreased reactive oxygen species and down-regulation of UCP2. However, UCP1, 3, 4 and 5 gene expression was unaffected. UCP2 transcription was controlled by SMAD4. Together, these studies suggest a tumor-promoting function of UCP2 in breast cancer. In summary, our studies demonstrate that i) the Warburg Effect is mediated by UCP2; ii) UCP2 is over-expressed in breast and many other cancers; iii) UCP2 promotes tumorigenic properties in vitro and in vivo and iv) genipin suppresses the tumor promoting function of UCP2.
doi:10.1371/journal.pone.0024792
PMCID: PMC3174207
PMID: 21935467
Underhill, Peter A | Myres, Natalie M | Rootsi, Siiri | Metspalu, Mait | Zhivotovsky, Lev A | King, Roy J | Lin, Alice A | Chow, Cheryl-Emiliane T | Semino, Ornella | Battaglia, Vincenza | Kutuev, Ildus | Järve, Mari | Chaubey, Gyaneshwer | Ayub, Qasim | Mohyuddin, Aisha | Mehdi, S Qasim | Sengupta, Sanghamitra | Rogaev, Evgeny I | Khusnutdinova, Elza K | Pshenichnov, Andrey | Balanovsky, Oleg | Balanovska, Elena | Jeran, Nina | Augustin, Dubravka Havas | Baldovic, Marian | Herrera, Rene J | Thangaraj, Kumarasamy | Singh, Vijay | Singh, Lalji | Majumder, Partha | Rudan, Pavao | Primorac, Dragan | Villems, Richard | Kivisild, Toomas
doi:10.1038/ejhg.2010.65
PMCID: PMC2987416
Mehrotra, Sanjana | Oommen, Joyce | Mishra, Anshuman | Sudharshan, Medhavi | Tiwary, Puja | Jamieson, Sarra E | Fakiola, Michaela | Rani, Deepa Selvi | Thangaraj, Kumarasamy | Rai, Madhukar | Sundar, Shyam | Blackwell, Jenefer M
Background
SLC11A1 has pleiotropic effects on macrophage function and remains a strong candidate for infectious disease susceptibility. 5' and/or 3' polymorphisms have been associated with tuberculosis, leprosy, and visceral leishmaniasis (VL). Most studies undertaken to date were under-powered, and none has been replicated within a population. Association with tuberculosis has replicated variably across populations. Here we investigate SLC11A1 and VL in India.
Methods
Nine polymorphisms (rs34448891, rs7573065, rs2276631, rs3731865, rs17221959, rs2279015, rs17235409, rs17235416, rs17229009) that tag linkage disequilibrium blocks across SLC11A1 were genotyped in primary family-based (313 cases; 176 families) and replication (941 cases; 992 controls) samples. Family- and population-based analyses were performed to look for association between SLC11A1 variants and VL. Quantitative RT/PCR was used to compare SLC11A1 expression in mRNA from paired splenic aspirates taken before and after treatment from 24 VL patients carrying different genotypes at the functional promoter GTn polymorphism (rs34448891).
Results
No associations were observed between VL and polymorphisms at SLC11A1 that were either robust to correction for multiple testing or replicated across primary and replication samples. No differences in expression of SLC11A1 were observed when comparing pre- and post-treatment samples, or between individuals carrying different genotypes at the GTn repeat.
Conclusions
This is the first well-powered study of SLC11A1 as a candidate for VL, which we conclude does not have a major role in regulating VL susceptibility in India.
doi:10.1186/1471-2350-12-71
PMCID: PMC3128845
PMID: 21599885
SLC11A1; visceral leishmaniasis; genetic susceptibility
Underhill, Peter A | Myres, Natalie M | Rootsi, Siiri | Metspalu, Mait | Zhivotovsky, Lev A | King, Roy J | Lin, Alice A | Chow, Cheryl-Emiliane T | Semino, Ornella | Battaglia, Vincenza | Kutuev, Ildus | Järve, Mari | Chaubey, Gyaneshwer | Ayub, Qasim | Mohyuddin, Aisha | Mehdi, S Qasim | Sengupta, Sanghamitra | Rogaev, Evgeny I | Khusnutdinova, Elza K | Pshenichnov, Andrey | Balanovsky, Oleg | Balanovska, Elena | Jeran, Nina | Augustin, Dubravka Havas | Baldovic, Marian | Herrera, Rene J | Thangaraj, Kumarasamy | Singh, Vijay | Singh, Lalji | Majumder, Partha | Rudan, Pavao | Primorac, Dragan | Villems, Richard | Kivisild, Toomas
Human Y-chromosome haplogroup structure is largely circumscribed by continental boundaries. One notable exception to this general pattern is the young haplogroup R1a that exhibits post-Glacial coalescent times and relates the paternal ancestry of more than 10% of men in a wide geographic area extending from South Asia to Central East Europe and South Siberia. Its origin and dispersal patterns are poorly understood as no marker has yet been described that would distinguish European R1a chromosomes from Asian. Here we present frequency and haplotype diversity estimates for more than 2000 R1a chromosomes assessed for several newly discovered SNP markers that introduce the onset of informative R1a subdivisions by geography. Marker M434 has a low frequency and a late origin in West Asia bearing witness to recent gene flow over the Arabian Sea. Conversely, marker M458 has a significant frequency in Europe, exceeding 30% in its core area in Eastern Europe and comprising up to 70% of all M17 chromosomes present there. The diversity and frequency profiles of M458 suggest its origin during the early Holocene and a subsequent expansion likely related to a number of prehistoric cultural developments in the region. Its primary frequency and diversity distribution correlates well with some of the major Central and East European river basins where settled farming was established before its spread further eastward. Importantly, the virtual absence of M458 chromosomes outside Europe speaks against substantial patrilineal gene flow from East Europe to Asia, including to India, at least since the mid-Holocene.
doi:10.1038/ejhg.2009.194
PMCID: PMC2987245
PMID: 19888303
Y chromosome; haplogroup R1a; human evolution; population genetics
Eaaswarkhanth, Muthukrishnan | Haque, Ikramul | Ravesh, Zeinab | Romero, Irene Gallego | Meganathan, Poorlin Ramakodi | Dubey, Bhawna | Khan, Faizan Ahmed | Chaubey, Gyaneshwer | Kivisild, Toomas | Tyler-Smith, Chris | Singh, Lalji | Thangaraj, Kumarasamy
Islam is the second most practiced religion in India, next to Hinduism. It is still unclear whether the spread of Islam in India has been only a cultural transformation or is associated with detectable levels of gene flow. To estimate the contribution of West Asian and Arabian admixture to Indian Muslims, we assessed genetic variation in mtDNA, Y-chromosomal and LCT/MCM6 markers in 472, 431 and 476 samples, respectively, representing six Muslim communities from different geographical regions of India. We found that most of the Indian Muslim populations received their major genetic input from geographically close non-Muslim populations. However, low levels of likely sub-Saharan African, Arabian and West Asian admixture were also observed among Indian Muslims in the form of L0a2a2 mtDNA and E1b1b1a and J*(xJ2) Y-chromosomal lineages. The distinction between Iranian and Arabian sources was difficult to make with mtDNA and the Y chromosome, as the estimates were highly correlated because of similar gene pool compositions in the sources. In contrast, the LCT/MCM6 locus, which shows a clear distinction between the two sources, enabled us to rule out significant gene flow from Arabia. Overall, our results support a model according to which the spread of Islam in India was predominantly cultural conversion associated with minor but still detectable levels of gene flow from outside, primarily from Iran and Central Asia, rather than directly from the Arabian Peninsula.
doi:10.1038/ejhg.2009.168
PMCID: PMC2859343
PMID: 19809480
Indian Muslims; mtDNA; Y chromosome; Middle East; sub-Saharan; gene flow
Thangaraj, Kumarasamy | Naidu, B. Prathap | Crivellaro, Federica | Tamang, Rakesh | Upadhyay, Shashank | Sharma, Varun Kumar | Reddy, Alla G. | Walimbe, S. R. | Chaubey, Gyaneshwer | Kivisild, Toomas | Singh, Lalji | Cordaux, Richard
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.
doi:10.1371/journal.pone.0015283
PMCID: PMC3004917
PMID: 21187967
Metspalu, Mait | Chaubey, Gyaneshwer | Yunusbayev, Bayazit | Romero, Irene Gallego | Karmin, Monika | Mallick, Chandana Basu | Metspalu, Ene | Shanmugalakshmi, Sadagopal | Balakrishnan, Karuppiah | Thangaraj, Kumarasamy | Singh, Lalji | Pitchappan, Ramasamy | Kivisild, Toomas | Villems, Richard
doi:10.1186/gb-2010-11-s1-o8
PMCID: PMC3026235
Eaaswarkhanth, Muthukrishnan | Haque, Ikramul | Ravesh, Zeinab | Romero, Irene Gallego | Meganathan, Poorlin Ramakodi | Dubey, Bhawna | Khan, Faizan Ahmed | Chaubey, Gyaneshwer | Kivisild, Toomas | Tyler-Smith, Chris | Singh, Lalji | Thangaraj, Kumarasamy
Islam is the second-most practiced religion in India, next to Hinduism. It is still unclear whether the spread of Islam in India has been only a cultural transformation or was associated with detectable levels of gene flow. To estimate the contribution of West Asian and Arabian admixture to Indian Muslims we have assessed genetic variation in mtDNA, Y-chromosomal and LCT/MCM6 markers in 472, 431 and 476 samples, respectively, representing six Muslim communities from different geographical regions of India. We found that most of the Indian Muslim populations received their major genetic input from geographically close non-Muslim populations. However, low levels of likely Arabian and West Asian admixture were also observed among Indian Muslims in the form of L0a2a2 mtDNA and E1b1b1a and J*(xJ2) Y-chromosomal lineages. The distinction between Iranian and Arabian sources was difficult to make with mtDNA and the Y chromosome as the estimates were highly correlated due to similar gene pool compositions in the sources. In contrast, the LCT/MCM6 locus, which shows a clear distinction between the two sources, enabled us to rule out significant gene flow from Arabia. Overall, our results support a model according to which the spread of Islam in India was predominantly cultural conversion associated with minor but still detectable levels of gene flow from outside, primarily from Iran and Central Asia, rather than directly from the Arabian Peninsula.
doi:10.1038/ejhg.2009.168
PMCID: PMC2859343
PMID: 19809480
Indian Muslims; mtDNA; Y chromosome; Middle East; sub-Saharan; gene flow
The present study was carried out to assess the role of androgen receptor CAG repeat polymorphism and X chromosome inactivation (XCI) pattern among Indian PCOS women and controls which has not been hitherto explored and also to test the hypothesis that shorter CAG alleles would be preferentially activated in PCOS. CAG repeat polymorphism and X chromosome methylation patterns were compared between PCOS and non-PCOS women. 250 PCOS women and 299 controls were included for this study. Androgen receptor CAG repeat sizes, XCI percentages, and clinical and biochemical parameters were measured. The mean CAG repeat number is similar between the cases (18.74±0.13) and controls (18.73±0.12). The obese PCOS women were significantly more frequent in the <18 and >20 CAG repeat category than the lean PCOS women, yielding a highly significant odds (p = 0.001). Among the women with non-random X-inactivation, alleles with <19 repeats were more frequently activated among cases than controls (p = 0.33). CAG repeat polymorphism by itself cannot be considered as a useful marker for discriminating PCOS. We observed a trend of preferential activation of the shorter allele among the PCOS cases with non random XCI pattern. In the obese PCOS women, this microsatellite variation may account for the hyperandrogenicity to a larger extent than the lean PCOS women.
doi:10.1371/journal.pone.0012401
PMCID: PMC2928732
PMID: 20865044
Nature
2009;461(7263):489-494.
India has been underrepresented in genome-wide surveys of human variation. We analyze 25 diverse groups to provide strong evidence for two ancient populations, genetically divergent, that are ancestral to most Indians today. One, the “Ancestral North Indians” (ANI), is genetically close to Middle Easterners, Central Asians, and Europeans, while the other, the “Ancestral South Indians” (ASI), is as distinct from ANI and East Asians as they are from each other. By introducing methods that can estimate ancestry without accurate ancestral populations, we show that ANI ancestry ranges from 39-71% in India, and is higher in traditionally upper caste and Indo-European speakers. Groups with only ASI ancestry may no longer exist in mainland India. However, the Andamanese are an ASI-related group without ANI ancestry, showing that the peopling of the islands must have occurred before ANI-ASI gene flow on the mainland. Allele frequency differences between groups in India are larger than in Europe, reflecting strong founder effects whose signatures have been maintained for thousands of years due to endogamy. We therefore predict that there will be an excess of recessive diseases in India, different in each group, which should be possible to screen and map genetically.
doi:10.1038/nature08365
PMCID: PMC2842210
PMID: 19779445
Aruna, Meka | Nagaraja, Theeya | Andal, Sadaranga | Tarakeswari, Surapaneni | Sirisha, Pisapati V. S. | Reddy, Alla G. | Thangaraj, Kumarasamy | Singh, Lalji | Reddy, B. Mohan | Aziz, Syed A.
Background
We attempt to ascertain if the 3 linked single nucleotide polymorphisms (SNPs) of the Progesterone Receptor (PR) gene (exon 1: G 1031 C; S344T, exon 4: G 1978 T; L660V and exon 5: C 2310 T; H770H) and the PROGINS insertion in the intron G, between exons 7 and 8, are associated with Recurrent Spontaneous Abortion (RSA) in the Indian population.
Methodology/Principal Findings
A total of 143 women with RSA and 150 controls were sequenced for all the 8 exons looking for the above 3 linked SNPs of the PR gene earlier implicated in the RSA, as well as for any new SNPs that may be possibly found in the Indian population. PROGINS insertion was screened by electrophoresis. We did not find any new mutations, not observed earlier, in our population. Further, we did not find significant role of the *2 allele (representing the mutant allele at the three SNP loci) or the T2 allele (PROGINS insertion) in the manifestation of RSA. We also did not find an LD pattern between each of the 3 SNPs and the PROGINS insertion.
Conclusions/Significance
The results suggest that the PR gene mutations may not play any exclusive role in the manifestation of RSA, and instead, given significantly higher frequency of the *2 allele among the normal women, we surmise if it does not really confer a protective role among the Indian populations, albeit further studies are required in the heterogeneous populations of this region before making any conclusive statement.
doi:10.1371/journal.pone.0008712
PMCID: PMC2806831
PMID: 20090851
Thangaraj, Kumarasamy | Nandan, Amrita | Sharma, Vishwas | Sharma, Varun Kumar | Eaaswarkhanth, Muthukrishnan | Patra, Pradeep Kumar | Singh, Sandhya | Rekha, Sashi | Dua, Monika | Verma, Narendra | Reddy, Alla G. | Singh, Lalji | Ahmed, Niyaz
Background
The phylogeny of the indigenous Indian-specific mitochondrial DNA (mtDNA) haplogroups have been determined and refined in previous reports. Similar to mtDNA superhaplogroups M and N, a profusion of reports are also available for superhaplogroup R. However, there is a dearth of information on South Asian subhaplogroups in particular, including R8. Therefore, we ought to access the genealogy and pre-historic expansion of haplogroup R8 which is considered one of the autochthonous lineages of South Asia.
Methodology/Principal Findings
Upon screening the mtDNA of 5,836 individuals belonging to 104 distinct ethnic populations of the Indian subcontinent, we found 54 individuals with the HVS-I motif that defines the R8 haplogroup. Complete mtDNA sequencing of these 54 individuals revealed two deep-rooted subclades: R8a and R8b. Furthermore, these subclades split into several fine subclades. An isofrequency contour map detected the highest frequency of R8 in the state of Orissa. Spearman's rank correlation analysis suggests significant correlation of R8 occurrence with geography.
Conclusions/Significance
The coalescent age of newly-characterized subclades of R8, R8a (15.4±7.2 Kya) and R8b (25.7±10.2 Kya) indicates that the initial maternal colonization of this haplogroup occurred during the middle and upper Paleolithic period, roughly around 40 to 45 Kya. These results signify that the southern part of Orissa currently inhabited by Munda speakers is likely the origin of these autochthonous maternal deep-rooted haplogroups. Our high-resolution study on the genesis of R8 haplogroup provides ample evidence of its deep-rooted ancestry among the Orissa (Austro-Asiatic) tribes.
doi:10.1371/journal.pone.0006545
PMCID: PMC2718812
PMID: 19662095
Dhandapany, Perundurai S | Sadayappan, Sakthivel | Xue, Yali | Powell, Gareth T | Rani, Deepa Selvi | Nallari, Prathiba | Rai, Taranjit Singh | Khullar, Madhu | Soares, Pedro | Bahl, Ajay | Tharkan, Jagan Mohan | Vaideeswar, Pradeep | Rathinavel, Andiappan | Narasimhan, Calambur | Ayapati, Dharma Rakshak | Ayub, Qasim | Mehdi, S Qasim | Oppenheimer, Stephen | Richards, Martin B | Price, Alkes L | Patterson, Nick | Reich, David | Singh, Lalji | Tyler-Smith, Chris | Thangaraj, Kumarasamy
Heart failure is a leading cause of mortality in South Asians. However, its genetic etiology remains largely unknown1. Cardiomyopathies due to sarcomeric mutations are a major monogenic cause for heart failure (MIM600958). Here, we describe a deletion of 25 bp in the gene encoding cardiac myosin binding protein C (MYBPC3) that is associated with heritable cardiomyopathies and an increased risk of heart failure in Indian populations (initial study OR = 5.3 (95% CI = 2.3–13), P = 2 × 10−6; replication study OR = 8.59 (3.19–25.05), P = 3 × 10−8; combined OR = 6.99 (3.68–13.57), P = 4 × 10−11) and that disrupts cardiomyocyte structure in vitro. Its prevalence was found to be high (~4%) in populations of Indian subcontinental ancestry. The finding of a common risk factor implicated in South Asian subjects with cardiomyopathy will help in identifying and counseling individuals predisposed to cardiac diseases in this region.
doi:10.1038/ng.309
PMCID: PMC2697598
PMID: 19151713
The caste system has persisted in Indian Hindu society for around 3,500 years. Like the Y chromosome, caste is defined at birth, and males cannot change their caste. In order to investigate the genetic consequences of this system, we have analysed male-lineage variation in a sample of 227 Indian men of known caste, 141 from the Jaunpur district of Uttar Pradesh and 86 from the rest of India. We typed 131 Y-chromosomal binary markers and 16 microsatellites. We find striking evidence for male substructure: in particular, Brahmins and Kshatriyas (but not other castes) from Jaunpur each show low diversity and the predominance of a single distinct cluster of haplotypes. These findings confirm the genetic isolation and drift within the Jaunpur upper castes, which are likely to result from founder effects and social factors. In the other castes, there may be either larger effective population sizes, or less strict isolation, or both.
doi:10.1007/s00439-006-0282-2
PMCID: PMC2590678
PMID: 17075717
Y chromosome; haplotype; human population substructure; Indian caste system
Thangaraj, Kumarasamy | Chaubey, Gyaneshwer | SelviRani, Deepa | Singh, Vijay Kumar | Ismail, Thanseem | Carvalho-Silva, Denise | Metspalu, Mait | Bhaskar, LVKS | Reddy, Alla G. | Chandra, Sharat | Pande, Veena | Naidu, Prathap | Adarsh, Niharika | Verma, Abhilasha | Jyothi, Inaganti Amara | Mallick, Chandana Basu | Shrivastava, Nidhi | Devasena, Ragala | Kumari, Babita | Singh, Amit Kumar | Diwedi, Shailendra Kumar Dhar | Singh, Shefali | Rao, Geeta | Gupta, Pranav | Sonvane, Vartika | Kumari, Kavita | Basha, Afsar | Bhargavi, K. R | Lalremruata, Albert | Gupta, Aravind Kumar | Kaur, Gurukamal | Reddy, K. K. | Rao, A. Papa | Villems, Richard | Tyler-Smith, Chris | Singh, Lalji
We have analyzed 7137 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/tRNALys 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.8% and 0.5%, respectively. We have also defined a novel deep-rooting haplogroup M43 and identified the rare haplogroup H14 in Indian populations carrying the 9bp-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.
doi:10.1159/000114160
PMCID: PMC2588665
PMID: 18223312
South Asia; 9bp indel; mtDNA; Haplogroup
Chaubey, Gyaneshwer | Karmin, Monika | Metspalu, Ene | Metspalu, Mait | Selvi-Rani, Deepa | Singh, Vijay Kumar | Parik, Jüri | Solnik, Anu | Naidu, B Prathap | Kumar, Ajay | Adarsh, Niharika | Mallick, Chandana Basu | Trivedi, Bhargav | Prakash, Swami | Reddy, Ramesh | Shukla, Parul | Bhagat, Sanjana | Verma, Swati | Vasnik, Samiksha | Khan, Imran | Barwa, Anshu | Sahoo, Dipti | Sharma, Archana | Rashid, Mamoon | Chandra, Vishal | Reddy, Alla G | Torroni, Antonio | Foley, Robert A | Thangaraj, Kumarasamy | Singh, Lalji | Kivisild, Toomas | Villems, Richard
Background
Human genetic diversity observed in Indian subcontinent is second only to that of Africa. This implies an early settlement and demographic growth soon after the first 'Out-of-Africa' dispersal of anatomically modern humans in Late Pleistocene. In contrast to this perspective, linguistic diversity in India has been thought to derive from more recent population movements and episodes of contact. With the exception of Dravidian, which origin and relatedness to other language phyla is obscure, all the language families in India can be linked to language families spoken in different regions of Eurasia. Mitochondrial DNA and Y chromosome evidence has supported largely local evolution of the genetic lineages of the majority of Dravidian and Indo-European speaking populations, but there is no consensus yet on the question of whether the Munda (Austro-Asiatic) speaking populations originated in India or derive from a relatively recent migration from further East.
Results
Here, we report the analysis of 35 novel complete mtDNA sequences from India which refine the structure of Indian-specific varieties of haplogroup R. Detailed analysis of haplogroup R7, coupled with a survey of ~12,000 mtDNAs from caste and tribal groups over the entire Indian subcontinent, reveals that one of its more recently derived branches (R7a1), is particularly frequent among Munda-speaking tribal groups. This branch is nested within diverse R7 lineages found among Dravidian and Indo-European speakers of India. We have inferred from this that a subset of Munda-speaking groups have acquired R7 relatively recently. Furthermore, we find that the distribution of R7a1 within the Munda-speakers is largely restricted to one of the sub-branches (Kherwari) of northern Munda languages. This evidence does not support the hypothesis that the Austro-Asiatic speakers are the primary source of the R7 variation. Statistical analyses suggest a significant correlation between genetic variation and geography, rather than between genes and languages.
Conclusion
Our high-resolution phylogeographic study, involving diverse linguistic groups in India, suggests that the high frequency of mtDNA haplogroup R7 among Munda speaking populations of India can be explained best by gene flow from linguistically different populations of Indian subcontinent. The conclusion is based on the observation that among Indo-Europeans, and particularly in Dravidians, the haplogroup is, despite its lower frequency, phylogenetically more divergent, while among the Munda speakers only one sub-clade of R7, i.e. R7a1, can be observed. It is noteworthy that though R7 is autochthonous to India, and arises from the root of hg R, its distribution and phylogeography in India is not uniform. This suggests the more ancient establishment of an autochthonous matrilineal genetic structure, and that isolation in the Pleistocene, lineage loss through drift, and endogamy of prehistoric and historic groups have greatly inhibited genetic homogenization and geographical uniformity.
doi:10.1186/1471-2148-8-227
PMCID: PMC2529308
PMID: 18680585
Background
The domestic goat is one of the important livestock species of India. In the present study we assess genetic diversity of Indian goats using 17 microsatellite markers. Breeds were sampled from their natural habitat, covering different agroclimatic zones.
Results
The mean number of alleles per locus (NA) ranged from 8.1 in Barbari to 9.7 in Jakhrana goats. The mean expected heterozygosity (He) ranged from 0.739 in Barbari to 0.783 in Jakhrana goats. Deviations from Hardy-Weinberg Equilibrium (HWE) were statistically significant (P < 0.05) for 5 loci breed combinations. The DA measure of genetic distance between pairs of breeds indicated that the lowest distance was between Marwari and Sirohi (0.135). The highest distance was between Pashmina and Black Bengal. An analysis of molecular variance indicated that 6.59% of variance exists among the Indian goat breeds. Both a phylogenetic tree and Principal Component Analysis showed the distribution of breeds in two major clusters with respect to their geographic distribution.
Conclusion
Our study concludes that Indian goat populations can be classified into distinct genetic groups or breeds based on the microsatellites as well as mtDNA information.
doi:10.1186/1471-2156-9-11
PMCID: PMC2268706
PMID: 18226239