Genetic diversity in India and the inference of Eurasian population expansion

The genetic admixture and assimilation of Ahom: a historic migrant from Thailand to India

The Northeastern region of India is considered a gateway for modern humans' dispersal throughout Asia. This region is a mixture of various ethnic and indigenous populations amalgamating multiple ancestries. One reason for such amalgamation is that, South Asia experienced multiple historic migrations from various parts of the world. A few examples explored genetically are Jews, Parsis and Siddis. Ahom is a dynasty that historically migrated to India during the 12th century. However, this putative migration has not been studied genetically at high resolution. Therefore, to validate this historical evidence, we genotyped autosomal data of the Modern Ahom population residing in seven sister states of India. Principal Component and Admixture analyses haave suggested a substantial admixture of the Ahom population with the local Tibeto-Burman populations. Moreover, the haplotype-based analysis has linked these Ahom individuals mainly with the Kusunda (a language isolated from Nepal) and Khasi (an Austroasiatic population of Meghalaya). Such unexpected presence of widespread population affinities suggests that Ahom mixed and assimilated a wide variety of Trans-Himalayan populations inhabiting this region after the migration. In summary, we observed a significant deviation of Ahom from their ancestral homeland (Thailand) and extensive admixture and assimilation with the local South Asian populations.

Characterizing substructure via mixture modeling in large-scale genetic summary statistics

Genetic summary data are broadly accessible and highly useful including for risk prediction, causal inference, fine mapping, and incorporation of external controls. However, collapsing individual-level data into groups masks intra- and inter-sample heterogeneity, leading to confounding, reduced power, and bias. Ultimately, unaccounted substructure limits summary data usability, especially for understudied or admixed populations. Here, we present Summix2, a comprehensive set of methods and software based on a computationally efficient mixture model to estimate and adjust for substructure in genetic summary data. In extensive simulations and application to public data, Summix2 characterizes finer-scale population structure, identifies ascertainment bias, and identifies potential regions of selection due to local substructure deviation. Summix2 increases the robust use of diverse publicly available summary data resulting in improved and more equitable research.

Access and barriers to genomic classifiers for breast cancer and prostate cancer in India

The incidence of cancer in general, including breast and prostate cancer specifically, is increasing in India. Breast and prostate cancers have genomic classifiers developed to guide therapy decisions. However, these genomic classifiers are often inaccessible in India due to high cost. These classifiers may also be less suitable to the Indian population, as data primarily from patients in wealthy Western countries were used in developing these genomic classifiers. In addition to the limitations in using these existing genomic classifiers, developing and validating new genomic classifiers for breast and prostate cancer in India is challenging due to the heterogeneity in the Indian population. However, there are steps that can be taken to address the various barriers that currently exist for accurate, accessible genomic classifiers for cancer in India.

Development of a decentralized cohort for studying post-acute sequelae of COVID-19 in India in the Data4life Study

BACKGROUND

Decentralized, digital health studies can provide real-world evidence of the lasting effects of COVID-19 on physical, socioeconomic, psychological, and social determinant factors of health in India. Existing research cohorts, however, are small and were not designed for longitudinal collection of comprehensive data from India's diverse population. Data4Life is a nationwide, digitally enabled, health research initiative to examine the post-acute sequelae of COVID-19 across individuals, communities, and regions. Data4Life seeks to build an ethnically and geographically diverse population of at least 100,000 participants in India.

METHODS

Here we discuss the feasibility of developing a completely decentralized COVID-19 cohort in India through qualitative analysis of data collection procedures, participant characteristics, participant perspectives on recruitment and reported study motivation.

RESULTS

As of June 13th, 2022, more than 6,000 participants from 17 Indian states completed baseline surveys. Friend and family referral were identified as the most common recruitment method (64.8%) across all demographic groups. Helping family and friends was the primary reason reported for joining the study (61.5%).

CONCLUSIONS

Preliminary findings support the use of digital technology for rapid enrollment and data collection to develop large health research cohorts in India. This demonstrates the potential for expansion of digitally enabled health research in India. These findings also outline the value of person-to-person recruitment strategies when conducting digital health research in modern-day India. Qualitative analysis reveals opportunities to increase diversity and retention in real time. It also informs strategies for improving participant experiences in the current Data4Life initiative and future studies.

Prevalence of Traditional Risk Factors in First-Degree Relatives of Patients With Established Cardiovascular Disease

INTRODUCTION

World Health Organization (WHO)/International Society of Hypertension (ISH) risk prediction charts are useful for predicting 10-year combined myocardial infarction and stroke risk (fatal and non-fatal). Hence the current study was conducted to assess the 10-year risk of cardiovascular disease among adults in Ahmedabad, India.

AIMS

The primary aim of the study was to assess the cardiovascular risk among first-degree relatives of patients attending the outpatient clinic. Also, to create awareness regarding assessment of cardiovascular risk among the studied group.

METHODS AND MATERIALS

A cross-sectional study was carried out among 372 first-degree relatives of patients at an out-patient cardiology clinic present in Vadaj, Ahmedabad. The WHO/ISH risk prediction chart for South-East Asia Region D (SEAR D) was used for calculating the 10-year cardiovascular risk.

RESULTS

A maximum (80.10%) of the study participants were in the low-risk (<10%) category followed by 8.33% for moderate-risk (10-20%), 7.25% for moderately high-risk (20-30%), 2.42% for high-risk (30-40%) and 1.88% for very high-risk (>40%).

CONCLUSION

WHO/ISH risk prediction charts provide a quick and effective way to assess and categorize the population in a low-resource setting which in turn helps in delivering focused intervention to the high-risk groups.

Identifying signatures of natural selection in Indian populations

In this study, we present the results of a genome-wide scan for signatures of positive selection using data from four tribal groups (Kokana, Warli, Bhil, and Pawara) and two caste groups (Deshastha Brahmin and Kunbi Maratha) from West of the Maharashtra State In India, as well as two samples of South Asian ancestry from the 1KG project (Gujarati Indian from Houston, Texas and Indian Telugu from UK). We used an outlier approach based on different statistics, including PBS, xpEHH, iHS, CLR, Tajima’s D, as well as two recently developed methods: Graph-aware Retrieval of Selective Sweeps (GRoSS) and Ascertained Sequentially Markovian Coalescent (ASMC). In order to minimize the risk of false positives, we selected regions that are outliers in all the samples included in the study using more than one method. We identified putative selection signals in 107 regions encompassing 434 genes. Many of the regions overlap with only one gene. The signals observed using microarray-based data are very consistent with our analyses using high-coverage sequencing data, as well as those identified with a novel coalescence-based method (ASMC). Importantly, at least 24 of these genomic regions have been identified in previous selection scans in South Asian populations or in other population groups. Our study highlights genomic regions that may have played a role in the adaptation of anatomically modern humans to novel environmental conditions after the out of Africa migration.

Extending precision medicine tools to populations at high risk of type 2 diabetes

In this Perspective, Shivani Misra and Jose C Florez discuss the application of precision medicine tools in under-represented populations.

An Alu insertion map of the Indian population: identification and analysis in 1021 genomes of the IndiGen project

Actively retrotransposing primate-specific Alu repeats display insertion-deletion (InDel) polymorphism through their insertion at new loci. In the global datasets, Indian populations remain under-represented and so do their Alu InDels. Here, we report the genomic landscape of Alu InDels from the recently released 1021 Indian Genomes (IndiGen) (available at https://clingen.igib.res.in/indigen). We identified 9239 polymorphic Alu insertions that include private (3831), rare (3974) and common (1434) insertions with an average of 770 insertions per individual. We achieved an 89% PCR validation of the predicted genotypes in 94 samples tested. About 60% of identified InDels are unique to IndiGen when compared to other global datasets; 23% of sites were shared with both SGDP and HGSVC; among these, 58% (1289 sites) were common polymorphisms in IndiGen. The insertions not only show a bias for genic regions, with a preference for introns but also for the associated genes showing enrichment for processes like cell morphogenesis and neurogenesis (P-value &lt; 0.05). Approximately, 60% of InDels mapped to genes present in the OMIM database. Finally, we show that 558 InDels can serve as ancestry informative markers to segregate global populations. This study provides a valuable resource for baseline Alu InDels that would be useful in population genomics.

Contrasting maternal and paternal genetic histories among five ethnic groups from Khyber Pakhtunkhwa, Pakistan

Northwest Pakistan has served as a point of entry to South Asia for different populations since ancient times. However, relatively little is known about the population genetic history of the people residing within this region. To better understand human dispersal in the region within the broader history of the subcontinent, we analyzed mtDNA diversity in 659 and Y-chromosome diversity in 678 individuals, respectively, from five ethnic groups (Gujars, Jadoons, Syeds, Tanolis and Yousafzais), from Swabi and Buner Districts, Khyber Pakhtunkhwa Province, Pakistan. The mtDNAs of all individuals were subject to control region sequencing and SNP genotyping, while Y-chromosomes were analyzed using 54 SNPs and 19 STR loci. The majority of the mtDNAs belonged to West Eurasian haplogroups, with the rest belonging to either South or East Asian lineages. Four of the five Pakistani populations (Gujars, Jadoons, Syeds, Yousafzais) possessed strong maternal genetic affinities with other Pakistani and Central Asian populations, whereas one (Tanolis) did not. Four haplogroups (R1a, R1b, O3, L) among the 11 Y-chromosome lineages observed among these five ethnic groups contributed substantially to their paternal genetic makeup. Gujars, Syeds and Yousafzais showed strong paternal genetic affinities with other Pakistani and Central Asian populations, whereas Jadoons and Tanolis had close affinities with Turkmen populations from Central Asia and ethnic groups from northeast India. We evaluate these genetic data in the context of historical and archeological evidence to test different hypotheses concerning their origins and biological relationships.

COVID-19: Impact on linguistic and genetic isolates of India

The rapid expansion of coronavirus SARS-CoV-2 has impacted various ethnic groups all over the world. The burden of infectious diseases including COVID-19 are generally reported to be higher for the Indigenous people. The historical knowledge have also suggested that the indigenous populations suffer more than the general populations in the pandemic. Recently, it has been reported that the indigenous groups of Brazil have been massively affected by COVID-19. Series of studies have shown that many of the indigenous communities reached at the verge of extinction due to this pandemic. Importantly, South Asia also has several indigenous and smaller communities, that are living in isolation. Till date, despite the two consecutive waves in India, there is no report on the impact of COVID-19 for indigenous tribes. Since smaller populations experiencing drift may have greater risk of such pandemic, we have analysed Runs of Homozygosity (ROH) among South Asian populations and identified several populations with longer homozygous segments. The longer runs of homozygosity at certain genomic regions may increases the susceptibility for COVID-19. Thus, we suggest extreme careful management of this pandemic among isolated populations of South Asia.

Lung cancer in the older population:Interactive effects of angiotensin converting enzyme gene polymorphism (rs 4340 ID) and tobacco addiction in risk assessment

BACKGROUND

rs4340ID polymorphism of angiotensin-converting enzyme (ACE) correlates with serum ACE levels in many known cancers. This study analyzed ACE rs4340 ID polymorphism in lung cancer (LC) in older patients of North India and correlated it with addiction status.

METHODS

The study enrolled all subjects aged 60 years and above with 154 LC and 205 healthy controls. Genotyping was done by polymerase chain reaction (PCR) and validated by sequencing of 10% of the sample. Statistical analysis was done by SPSS Statistics 21.

RESULTS

Genotype II was observed to have a significant 2.21-fold increased risk of LC as compared to the DD genotype and 3.43-folds enhanced risk with interaction of I allele with tobacco consumption habits as compared to D allele in LC was seen.

CONCLUSION

The risk of LC was higher with II genotype as compared to DD genotype. Interactive effect showed that I allele with tobacco habits may increase the risk of LC.

Disease characteristics and clinical outcomes of adults and children with anti-MDA-5 antibody-associated myositis: a prospective observational bicentric study

To study the demographic, clinical and serologic characteristics of anti-MDA5-positive DM from two geographically and ethnically disparate inception cohorts in India. To identify the clinical and serological parameters at inception that could predict mortality among these individuals. Individuals with anti-MDA5 antibody-positive DM diagnosed between 2017 and 2020 from two centres in India were prospectively followed up. The clinical and serological characteristics at baseline and the treatment outcome at follow-up were assessed for this study. Anti-MDA5 antibody was positive in 25 (7.5%) out of the 330 individuals with myositis. These 25 (21 adults, 4 juvenile) patients were followed up for a median duration of 14 months. Among adults, a majority had cutaneous manifestations 21 (84%) followed by, arthritis 17 (80%), and interstitial lung disease 12 (ILD, 57.1%). Four (19%) had rapidly progressive ILD (RP-ILD). Eight (38%) presented as clinically amyopathic DM. Among cutaneous manifestations, majority (62%) had classic features (gottron's papules/sign, heliotrope rash) while 8 (38%) had cutaneous ulceration and 2 each had periorbital edema and tendon rupture. Eight (38%) were positive for anti-Ro-52 antibody. Out of 21 adults, 8 (38%) succumbed to the diseases. RP-ILD (n = 4; 19%), ulcerative gottron's (n = 5) and anti-Ro-52 (n = 8) were significantly associated with mortality (p < 0.05). Upon binary logistic regression, positive anti-Ro-52 antibody predicted mortality [HR 17.3 (95%CI 1.4-210, p = 0.025)]. All juvenile anti-MDA5-positive DMs had classic cutaneous features with 2 of them having ulcerative gottron's. None of the juvenile patients had ILD and everyone survived till the last follow-up. Indian adults with anti-MDA5 DM have high mortality. Rarer atypical features like tendon rupture or periorbital edema could assist in diagnosis. Ulcerative gottron's, positive anti-Ro 52 antibodies, and RP-ILD are valuable clinical-serological markers that portend poor prognosis.

A data set of variants derived from 1455 clinical and research exomes is efficient in variant prioritization for early-onset monogenic disorders in Indians

Given the genomic uniqueness, a local data set is most desired for Indians, who are underrepresented in existing public databases. We hypothesize patients with rare monogenic disorders and their family members can provide a reliable source of common variants in the population. Exome sequencing (ES) data from families with rare Mendelian disorders was aggregated from five centers in India. The dataset was refined by excluding related individuals and removing the disease-causing variants (refined cohort). The efficiency of these data sets was assessed in a new set of 50 exomes against gnomAD and GenomeAsia. Our original cohort comprised 1455 individuals from 1203 families. The refined cohort had 836 unrelated individuals that retained 1,251,064 variants with 181,125 population-specific and 489,618 common variants. The allele frequencies from our cohort helped to define 97,609 rare variants in gnomAD and 44,520 rare variants in GenomeAsia as common variants in our population. Our variant dataset provided an additional 1.7% and 0.1% efficiency for prioritizing heterozygous and homozygous variants respectively for rare monogenic disorders. We observed additional 19 genes/human knockouts. We list carrier frequency for 142 recessive disorders. This is a large and useful resource of exonic variants for Indians. Despite limitations, datasets from patients are efficient tools for variant prioritization in a resource-limited setting.

Most frequent South Asian haplotypes of ACE2 share identity by descent with East Eurasian populations

It was shown that the human Angiotensin-converting enzyme 2 (ACE2) is the receptor of recent coronavirus SARS-CoV-2, and variation in this gene may affect the susceptibility of a population. Therefore, we have analysed the sequence data of ACE2 among 393 samples worldwide, focusing on South Asia. Genetically, South Asians are more related to West Eurasian populations rather than to East Eurasians. In the present analyses of ACE2, we observed that the majority of South Asian haplotypes are closer to East Eurasians rather than to West Eurasians. The phylogenetic analysis suggested that the South Asian haplotypes shared with East Eurasians involved two unique event polymorphisms (rs4646120 and rs2285666). In contrast with the European/American populations, both of the SNPs have largely similar frequencies for East Eurasians and South Asians, Therefore, it is likely that among the South Asians, host susceptibility to the novel coronavirus SARS-CoV-2 will be more similar to that of East Eurasians rather than to that of Europeans.

The genetic legacy of continental scale admixture in Indian Austroasiatic speakers

Surrounded by speakers of Indo-European, Dravidian and Tibeto-Burman languages, around 11 million Munda (a branch of Austroasiatic language family) speakers live in the densely populated and genetically diverse South Asia. Their genetic makeup holds components characteristic of South Asians as well as Southeast Asians. The admixture time between these components has been previously estimated on the basis of archaeology, linguistics and uniparental markers. Using genome-wide genotype data of 102 Munda speakers and contextual data from South and Southeast Asia, we retrieved admixture dates between 2000–3800 years ago for different populations of Munda. The best modern proxies for the source populations for the admixture with proportions 0.29/0.71 are Lao people from Laos and Dravidian speakers from Kerala in India. The South Asian population(s), with whom the incoming Southeast Asians intermixed, had a smaller proportion of West Eurasian genetic component than contemporary proxies. Somewhat surprisingly Malaysian Peninsular tribes rather than the geographically closer Austroasiatic languages speakers like Vietnamese and Cambodians show highest sharing of IBD segments with the Munda. In addition, we affirmed that the grouping of the Munda speakers into North and South Munda based on linguistics is in concordance with genome-wide data.

General obesity and Cardiovascular diseases among Gaur Brahmin population of NCR/Delhi

AIMS

Cardiovascular diseases are one of the leading causes of mortality and morbidity among human beings. The presence of endemic Cardiovascular diseases and their risk factors differ from population to population. The Cardiovascular diseases associated risk factors are sub-categorised into two forms, one is traditional and the other is non-traditional risk factors. The present study shows the prevalence of both risk factors and its association with Cardiovascular diseases, especially with reference to general obesity.

MATERIALS AND METHODS

The present study includes a total of 506 Gaur Brahmins residing in Delhi and National Capital Region India. Household survey was conducted and data were collected by using pre-tested interview schedule. Somatometric measurements were taken following the international standard techniques. Approx 5 ml blood was collected from each individual unrelated up to the first cousion. The serum was used to analyse the lipid profiles and fasting glucose level. All necessary statistical analyses were performed using the Statistical Package for Social Sciences (SPSS) and MS Excel. The ethical clearance was obtained from the Ethical Committee of the Department of Anthropology, University of Delhi, Delhi.

RESULTS AND CONCLUSION

The mean value of Somatometric variables such as Body Mass Index, Waist circumference and Waist-hip ratio and physiological variables DBP and SBP (diastolic blood pressure and systolic blood pressure) were found to be higher than their respective ranges in the studied population. General obesity, though found to be less common in this population as compared to abdominal obesity, but it is found to be contributing to dyslipidemia.

Activation of Epithelial-Mesenchymal Transition and Altered β-Catenin Signaling in a Novel Indian Colorectal Carcinoma Cell Line

Colorectal cancer is the third major cause of cancer-related mortality worldwide. The upward trend in incidence and mortality rates, poor sensitivity to conventional therapies and a dearth of early diagnostic parameters pose a huge challenge in the management of colorectal cancer in India. Due to the high level of genetic diversity present in the Indian population, unraveling the genetic contributions toward pathogenesis is key for understanding the etiology of colorectal cancer and in reversing this trend. We have established a novel cell line, MBC02, from an Indian colorectal cancer patient and have carried out extensive molecular characterization to unravel the pathological alterations in this cell line. In-depth molecular analysis of MBC02 revealed suppression of E-cadherin expression, concomitant with overexpression of EMT related molecules, which manifested in the form of highly migratory and invasive cells. Loss of membrane-tethered E-cadherin released β-catenin from the adherens junction resulting in its cytoplasmic and nuclear accumulation and consequently, upregulation of c-Myc. MBC02 also showed dramatic transcriptional upregulation of β-catenin. Remarkably, we observed significantly elevated proteasome activity that perhaps co-evolved to compensate for the unnaturally high mRNA level of β-catenin to regulate the increased protein load. In addition, there was substantial misregulation of other clinically relevant signaling pathways that have clinical relevance in the pathogenesis of colorectal cancer. Our findings pave the way toward understanding the molecular differences that could define pathogenesis in cancers originating in the Indian population.

The Genetic Ancestry of Modern Indus Valley Populations from Northwest India

The Indus Valley has been the backdrop for several historic and prehistoric population movements between South Asia and West Eurasia. However, the genetic structure of present-day populations from Northwest India is poorly characterized. Here we report new genome-wide genotype data for 45 modern individuals from four Northwest Indian populations, including the Ror, whose long-term occupation of the region can be traced back to the early Vedic scriptures. Our results suggest that although the genetic architecture of most Northwest Indian populations fits well on the broader North-South Indian genetic cline, culturally distinct groups such as the Ror stand out by being genetically more akin to populations living west of India; such populations include prehistorical and early historical ancient individuals from the Swat Valley near the Indus Valley. We argue that this affinity is more likely a result of genetic continuity since the Bronze Age migrations from the Steppe Belt than a result of recent admixture. The observed patterns of genetic relationships both with modern and ancient West Eurasians suggest that the Ror can be used as a proxy for a population descended from the Ancestral North Indian (ANI) population. Collectively, our results show that the Indus Valley populations are characterized by considerable genetic heterogeneity that has persisted over thousands of years.

The genetic makings of South Asia

South Asia is home for more than a billion people culturally structured into innumerable groups practicing different levels of endogamy. Linguistically South Asia is broadly characterized by four major language families which has served as access way for disentangling the genetic makings of South Asia. In this review we shall give brief account on the recent developments in the field. Advances are made in two fronts simultaneously. Whole genome characterisation of many extant South Asians paint the picture of the genetic diversity and its implications to health-care. On the other hand ancient DNA studies, which are finally reaching South Asia, provide new incites to the demographic history of the subcontinent. Before the spread of agriculture, South Asia was likely inhabited by hunter-gatherer groups deriving much of their ancestry from a population that split from the rest of humanity soon after expanding from Africa. Early Iranian agriculturalists mixing with these local hunter-gatherers probably formed the population that flourished during the blossoming of the Indus Valley Civilisation. Further admixture with the still persisting HG groups and population(s) from the Eurasian Steppe, formed the two ancestral populations (ANI and ASI), the north-south mixing pattern of whom is known today as the 'Indian Cline'. Studies on natural selection in South Asia have so far revealed strong signals of sweeps that are shared with West Eurasians. Future studies will have to fully unlock the aDNA promise for South Asia.

Carriers of mitochondrial DNA macrohaplogroup L3 basal lineages migrated back to Africa from Asia around 70,000 years ago

Background

The main unequivocal conclusion after three decades of phylogeographic mtDNA studies is the African origin of all extant modern humans. In addition, a southern coastal route has been argued for to explain the Eurasian colonization of these African pioneers. Based on the age of macrohaplogroup L3, from which all maternal Eurasian and the majority of African lineages originated, the out-of-Africa event has been dated around 60-70 kya. On the opposite side, we have proposed a northern route through Central Asia across the Levant for that expansion and, consistent with the fossil record, we have dated it around 125 kya. To help bridge differences between the molecular and fossil record ages, in this article we assess the possibility that mtDNA macrohaplogroup L3 matured in Eurasia and returned to Africa as basal L3 lineages around 70 kya.

Results

The coalescence ages of all Eurasian (M,N) and African (L3 ) lineages, both around 71 kya, are not significantly different. The oldest M and N Eurasian clades are found in southeastern Asia instead near of Africa as expected by the southern route hypothesis. The split of the Y-chromosome composite DE haplogroup is very similar to the age of mtDNA L3. An Eurasian origin and back migration to Africa has been proposed for the African Y-chromosome haplogroup E. Inside Africa, frequency distributions of maternal L3 and paternal E lineages are positively correlated. This correlation is not fully explained by geographic or ethnic affinities. This correlation rather seems to be the result of a joint and global replacement of the old autochthonous male and female African lineages by the new Eurasian incomers.

Conclusions

These results are congruent with a model proposing an out-of-Africa migration into Asia, following a northern route, of early anatomically modern humans carrying pre-L3 mtDNA lineages around 125 kya, subsequent diversification of pre-L3 into the basal lineages of L3, a return to Africa of Eurasian fully modern humans around 70 kya carrying the basal L3 lineages and the subsequent diversification of Eurasian-remaining L3 lineages into the M and N lineages in the outside-of-Africa context, and a second Eurasian global expansion by 60 kya, most probably, out of southeast Asia. Climatic conditions and the presence of Neanderthals and other hominins might have played significant roles in these human movements. Moreover, recent studies based on ancient DNA and whole-genome sequencing are also compatible with this hypothesis.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1211-4) contains supplementary material, which is available to authorized users.

The Geographic Origins of Ethnic Groups in the Indian Subcontinent: Exploring Ancient Footprints with Y-DNA Haplogroups

Several studies have evaluated the movements of large populations to the Indian subcontinent; however, the ancient geographic origins of smaller ethnic communities are not clear. Although historians have attempted to identify the origins of some ethnic groups, the evidence is typically anecdotal and based upon what others have written before. In this study, recent developments in DNA science were assessed to provide a contemporary perspective by analyzing the Y chromosome haplogroups of some key ethnic groups and tracing their ancient geographical origins from genetic markers on the Y-DNA haplogroup tree. A total of 2,504 Y-DNA haplotypes, representing 50 different ethnic groups in the Indian subcontinent, were analyzed. The results identified 14 different haplogroups with 14 geographic origins for these people. Moreover, every ethnic group had representation in more than one haplogroup, indicating multiple geographic origins for these communities. The results also showed that despite their varied languages and cultural differences, most ethnic groups shared some common ancestors because of admixture in the past. These findings provide new insights into the ancient geographic origins of ethnic groups in the Indian subcontinent. With about 2,000 other ethnic groups and tribes in the region, it is expected that more scientific discoveries will follow, providing insights into how, from where, and when the ancestors of these people arrived in the subcontinent to create so many different communities.

SAGE: a comprehensive resource of genetic variants integrating South Asian whole genomes and exomes

South Asia is home to \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$\sim $\end{document}20% of the world population and characterized by distinct ethnic, linguistic, cultural and genetic lineages. Only limited representative samples from the region have found its place in large population-scale international genome projects. The recent availability of genome scale data from multiple populations and datasets from South Asian countries in public domain motivated us to integrate the data into a comprehensive resource. In the present study, we have integrated a total of six datasets encompassing 1213 human exomes and genomes to create a compendium of 154 814 557 genetic variants and adding a total of 69 059 255 novel variants. The variants were systematically annotated using public resources and along with the allele frequencies are available as a browsable-online resource South Asian genomes and exomes. As a proof of principle application of the data and resource for genetic epidemiology, we have analyzed the pathogenic genetic variants causing retinitis pigmentosa. Our analysis reveals the genetic landscape of the disease and suggests subset of genetic variants to be highly prevalent in South Asia.

Reference genome assessment from a population scale perspective: an accurate profile of variability and noise

MOTIVATION

Current plant and animal genomic studies are often based on newly assembled genomes that have not been properly consolidated. In this scenario, misassembled regions can easily lead to false-positive findings. Despite quality control scores are included within genotyping protocols, they are usually employed to evaluate individual sample quality rather than reference sequence reliability. We propose a statistical model that combines quality control scores across samples in order to detect incongruent patterns at every genomic region. Our model is inherently robust since common artifact signals are expected to be shared between independent samples over misassembled regions of the genome.

RESULTS

The reliability of our protocol has been extensively tested through different experiments and organisms with accurate results, improving state-of-the-art methods. Our analysis demonstrates synergistic relations between quality control scores and allelic variability estimators, that improve the detection of misassembled regions, and is able to find strong artifact signals even within the human reference assembly. Furthermore, we demonstrated how our model can be trained to properly rank the confidence of a set of candidate variants obtained from new independent samples.

AVAILABILITY AND IMPLEMENTATION

This tool is freely available at http://gitlab.com/carbonell/ces.

CONTACT

jcarbonell.cipf@gmail.com or joaquin.dopazo@juntadeandalucia.es.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Y-STR Haplogroup Diversity in the Jat Population Reveals Several Different Ancient Origins

The Jats represent a large ethnic community that has inhabited the northwest region of India and Pakistan for several thousand years. It is estimated the community has a population of over 123 million people. Many historians and academics have asserted that the Jats are descendants of Aryans, Scythians, or other ancient people that arrived and lived in northern India at one time. Essentially, the specific origin of these people has remained a matter of contention for a long time. This study demonstrated that the origins of Jats can be clarified by identifying their Y-chromosome haplogroups and tracing their genetic markers on the Y-DNA haplogroup tree. A sample of 302 Y-chromosome haplotypes of Jats in India and Pakistan was analyzed. The results showed that the sample population had several different lines of ancestry and emerged from at least nine different geographical regions of the world. It also became evident that the Jats did not have a unique set of genes, but shared an underlying genetic unity with several other ethnic communities in the Indian subcontinent. A startling new assessment of the genetic ancient origins of these people was revealed with DNA science.

Carriers of mitochondrial DNA macrohaplogroup R colonized Eurasia and Australasia from a southeast Asia core area

BACKGROUND

The colonization of Eurasia and Australasia by African modern humans has been explained, nearly unanimously, as the result of a quick southern coastal dispersal route through the Arabian Peninsula, the Indian subcontinent, and the Indochinese Peninsula, to reach Australia around 50 kya. The phylogeny and phylogeography of the major mitochondrial DNA Eurasian haplogroups M and N have played the main role in giving molecular genetics support to that scenario. However, using the same molecular tools, a northern route across central Asia has been invoked as an alternative that is more conciliatory with the fossil record of East Asia. Here, we assess as the Eurasian macrohaplogroup R fits in the northern path.

RESULTS

Haplogroup U, with a founder age around 50 kya, is one of the oldest clades of macrohaplogroup R in western Asia. The main branches of U expanded in successive waves across West, Central and South Asia before the Last Glacial Maximum. All these dispersions had rather overlapping ranges. Some of them, as those of U6 and U3, reached North Africa. At the other end of Asia, in Wallacea, another branch of macrohaplogroup R, haplogroup P, also independently expanded in the area around 52 kya, in this case as isolated bursts geographically well structured, with autochthonous branches in Australia, New Guinea, and the Philippines.

CONCLUSIONS

Coeval independently dispersals around 50 kya of the West Asia haplogroup U and the Wallacea haplogroup P, points to a halfway core area in southeast Asia as the most probable centre of expansion of macrohaplogroup R, what fits in the phylogeographic pattern of its ancestor, macrohaplogroup N, for which a northern route and a southeast Asian origin has been already proposed.

Extremely low-coverage whole genome sequencing in South Asians captures population genomics information

BACKGROUND

The cost of Whole Genome Sequencing (WGS) has decreased tremendously in recent years due to advances in next-generation sequencing technologies. Nevertheless, the cost of carrying out large-scale cohort studies using WGS is still daunting. Past simulation studies with coverage at ~2x have shown promise for using low coverage WGS in studies focused on variant discovery, association study replications, and population genomics characterization. However, the performance of low coverage WGS in populations with a complex history and no reference panel remains to be determined.

RESULTS

South Indian populations are known to have a complex population structure and are an example of a major population group that lacks adequate reference panels. To test the performance of extremely low-coverage WGS (EXL-WGS) in populations with a complex history and to provide a reference resource for South Indian populations, we performed EXL-WGS on 185 South Indian individuals from eight populations to ~1.6x coverage. Using two variant discovery pipelines, SNPTools and GATK, we generated a consensus call set that has ~90% sensitivity for identifying common variants (minor allele frequency ≥ 10%). Imputation further improves the sensitivity of our call set. In addition, we obtained high-coverage for the whole mitochondrial genome to infer the maternal lineage evolutionary history of the Indian samples.

CONCLUSIONS

Overall, we demonstrate that EXL-WGS with imputation can be a valuable study design for variant discovery with a dramatically lower cost than standard WGS, even in populations with a complex history and without available reference data. In addition, the South Indian EXL-WGS data generated in this study will provide a valuable resource for future Indian genomic studies.

A genetic chronology for the Indian Subcontinent points to heavily sex-biased dispersals

BACKGROUND

India is a patchwork of tribal and non-tribal populations that speak many different languages from various language families. Indo-European, spoken across northern and central India, and also in Pakistan and Bangladesh, has been frequently connected to the so-called "Indo-Aryan invasions" from Central Asia ~3.5 ka and the establishment of the caste system, but the extent of immigration at this time remains extremely controversial. South India, on the other hand, is dominated by Dravidian languages. India displays a high level of endogamy due to its strict social boundaries, and high genetic drift as a result of long-term isolation which, together with a very complex history, makes the genetic study of Indian populations challenging.

RESULTS

We have combined a detailed, high-resolution mitogenome analysis with summaries of autosomal data and Y-chromosome lineages to establish a settlement chronology for the Indian Subcontinent. Maternal lineages document the earliest settlement ~55-65 ka (thousand years ago), and major population shifts in the later Pleistocene that explain previous dating discrepancies and neutrality violation. Whilst current genome-wide analyses conflate all dispersals from Southwest and Central Asia, we were able to tease out from the mitogenome data distinct dispersal episodes dating from between the Last Glacial Maximum to the Bronze Age. Moreover, we found an extremely marked sex bias by comparing the different genetic systems.

CONCLUSIONS

Maternal lineages primarily reflect earlier, pre-Holocene processes, and paternal lineages predominantly episodes within the last 10 ka. In particular, genetic influx from Central Asia in the Bronze Age was strongly male-driven, consistent with the patriarchal, patrilocal and patrilineal social structure attributed to the inferred pastoralist early Indo-European society. This was part of a much wider process of Indo-European expansion, with an ultimate source in the Pontic-Caspian region, which carried closely related Y-chromosome lineages, a smaller fraction of autosomal genome-wide variation and an even smaller fraction of mitogenomes across a vast swathe of Eurasia between 5 and 3.5 ka.

Chronic pancreatitis in Eastern India: Experience from a tertiary care center

There is a wide variation in the clinical presentation of chronic pancreatitis (CP) in the different parts of India. Data regarding the clinical profile of CP from eastern India are scarce. We describe the clinical and demographic profiles of patients with CP in eastern India. Consecutive patients were evaluated for the clinical presentation, etiology and complication of CP. One hundred and thirty-nine patients with CP (mean age 39.57±14.88 years; M/F 3.48:1) were included. Idiopathic CP (50.35%) was the most common etiology followed by alcohol (33.81%); 68.34% had calcific CP and 31.65% had noncalcific CP. The median duration of symptoms was 24 (1-240) months. Pain was the most common symptom, being present in 93.52% of the patients. Diabetes, steatorrhea and pseudocyst were present in 45.32%, 14.38% and 7.19% of the cases, respectively. Moderate to severe anemia was revealed in 16.53% of the patients. Benign biliary stricture was diagnosed in 19.42% of the cases (symptomatic in 6.47%). The common radiological findings were the following: pancreatic calculi (68.34%), dilated pancreatic duct (PD) (58.99%), parenchymal atrophy (25.89%) and PD stricture (23.74%). In our center, idiopathic CP followed by alcoholic CP was the most frequent form of CP. Tropical CP was distinctly uncommon.

Characterising private and shared signatures of positive selection in 37 Asian populations

The Asian Diversity Project (ADP) assembled 37 cosmopolitan and ethnic minority populations in Asia that have been densely genotyped across over half a million markers to study patterns of genetic diversity and positive natural selection. We performed population structure analyses of the ADP populations and divided these populations into four major groups based on their genographic information. By applying a highly sensitive algorithm haploPS to locate genomic signatures of positive selection, 140 distinct genomic regions exhibiting evidence of positive selection in at least one population were identified. We examined the extent of signal sharing for regions that were selected in multiple populations and observed that populations clustered in a similar fashion to that of how the ancestry clades were phylogenetically defined. In particular, populations predominantly located in South Asia underwent considerably different adaptation as compared with populations from the other geographical regions. Signatures of positive selection present in multiple geographical regions were predicted to be older and have emerged prior to the separation of the populations in the different regions. In contrast, selection signals present in a single population group tended to be of lower frequencies and thus can be attributed to recent evolutionary events.

Molecular Characterization and Meta-Analysis of Gut Microbial Communities Illustrate Enrichment of Prevotella and Megasphaera in Indian Subjects

The gut microbiome has varied impact on the wellbeing of humans. It is influenced by different factors such as age, dietary habits, socio-economic status, geographic location, and genetic makeup of individuals. For devising microbiome-based therapies, it is crucial to identify population specific features of the gut microbiome. Indian population is one of the most ethnically, culturally, and geographically diverse, but the gut microbiome features remain largely unknown. The present study describes gut microbial communities of healthy Indian subjects and compares it with the microbiota from other populations. Based on large differences in alpha diversity indices, abundance of 11 bacterial phyla and individual specific OTUs, we report inter-individual variations in gut microbial communities of these subjects. While the gut microbiome of Indians is different from that of Americans, it shared high similarity to individuals from the Indian subcontinent i.e., Bangladeshi. Distinctive feature of Indian gut microbiota is the predominance of genus Prevotella and Megasphaera. Further, when compared with other non-human primates, it appears that Indians share more OTUs with omnivorous mammals. Our metagenomic imputation indicates higher potential for glycan biosynthesis and xenobiotic metabolism in these subjects. Our study indicates urgent need of identification of population specific microbiome biomarkers of Indian subpopulations to have more holistic view of the Indian gut microbiome and its health implications.

Frailty and geography: should these two factors be added to the ABCDE contemporary guide to diabetes therapy?

On the road towards personalized treatments for type 2 diabetes, we suggest here that two parameters could be added to the ABCDE algorithm, 'F' for frailty and 'G' for geography. Indeed, the progressive ageing of population is causing a simultaneous increase of frailty worldwide. The identification of the optimal therapeutic approach is often difficult in frail subjects because of the complexity of 'frailty syndrome'. Nevertheless, given the relevance of diabetes in the development and progression of frailty, a safe and effective cure of diabetes is extremely important to guarantee a good medical outcome. There are few data about diabetes treatment in this delicate category of patients, and the choice of the appropriate therapy mostly remains a challenge. Moreover, type 2 diabetes affects more than 382 million people of different countries, races and ethnicities. To face the lack of solid evidence-based medicine for the treatment of diabetes in different ethnic groups, it is extremely important to increase knowledge about the different pathophysiology of diabetes according to ethnicity. In this way, a tailored approach to treatment of various ethnic groups living in the same or different regions can eventually be developed. Copyright © 2016 John Wiley & Sons, Ltd.

Genetic factors affecting drug disposition in Asian cancer patients

INTRODUCTION

In the era of genomic medicine, it is increasingly recognized that ethnogeographic differences in drug pharmacology exist between Asian and other populations. This is particularly pertinent to oncology, where drugs forming the backbone of chemotherapy often have narrow therapeutic windows and are frequently dosed close to maximally tolerable levels.

AREAS COVERED

At the population level, ancestry is important because historical-biogeographical confluences have shaped population genetics and pharmacoethnicity in the Asian race through allelic differentiation and interethnic differences in inheritance patterns of linkage disequilibrium. At the individual level, cis- and trans-acting germline polymorphisms and somatic mutations in genes encoding drug-metabolizing enzymes and transporters act in a multifactorial manner to determine drug disposition phenotype and clinical response in Asian cancer patients. A growing body of evidence also finds that complex genetic interactions and regulation, including a multiplicity of gene control mechanisms, are increasingly implicated in genotype-phenotype correlates than has hitherto been appreciated--potentially serving as the mechanistic links between hits in non-coding regions of genome-wide association studies and drug toxicity. Together, these genetic factors contribute to the clinical heterogeneity of drug disposition in Asian cancer patients.

EXPERT OPINION

This topic has broad relevance for the optimization and individualization of anticancer strategies in Asians.

Implications of Cardiovascular Disease Risk Assessment Using the WHO/ISH Risk Prediction Charts in Rural India

Cardiovascular disease (CVD) risk in India is currently assessed using the World Health Organization/International Society for Hypertension (WHO/ISH) risk prediction charts since no population-specific models exist. The WHO/ISH risk prediction charts have two versions-one with total cholesterol as a predictor (the high information (HI) model) and the other without (the low information (LI) model). However, information on the WHO/ISH risk prediction charts including guidance on which version to use and when, as well as relative performance of the LI and HI models, is limited. This article aims to, firstly, quantify the relative performance of the LI and HI WHO/ISH risk prediction (for WHO-South East Asian Region D) using data from rural India. Secondly, we propose a pre-screening (simplified) point-of-care (POC) test to identify patients who are likely to benefit from a total cholesterol (TC) test, and subsequently when the LI model is preferential to HI model. Analysis was performed using cross-sectional data from rural Andhra Pradesh collected in 2005 with recorded blood cholesterol measurements (N = 1066). CVD risk was computed using both LI and HI models, and high risk individuals who needed treatment(THR) were subsequently identified based on clinical guidelines. Model development for the POC assessment of a TC test was performed through three machine learning techniques: Support Vector Machine (SVM), Regularised Logistic Regression (RLR), and Random Forests (RF) along with a feature selection process. Disagreement in CVD risk predicted by LI and HI WHO/ISH models was 14.5% (n = 155; p<0.01) overall and comprised 36 clinically relevant THR patients (31% of patients identified as THR by using either model). Using two patient-specific parameters (age, systolic blood pressure), our POC assessment can pre-determine the benefit of TC testing and choose the appropriate risk model (out-of-sample AUCs:RF-0.85,SVM-0.84,RLR:0.82 and maximum sensitivity-98%). The identification of patients benefitting from a TC test for CVD risk stratification can aid planning for resource-allocation and save costs for large-scale screening programmes.

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

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

Distance-dependent patterns of molecular divergences in Tuatara mitogenomes

Population genetic models predict that populations that are geographically close to each other are expected to be genetically more similar to each other compared to those that are widely separate. However the patterns of relationships between geographic distance and molecular divergences at neutral and constrained regions of the genome are unclear. We attempted to clarify this relationship by sequencing complete mitochondrial genomes of the relic species Tuatara (Sphenodon punctatus) from ten offshore islands of New Zealand. We observed a positive relationship that showed a proportional increase in the neutral diversity at synonymous sites (dS), with increasing geographical distance. In contrast we showed that diversity at evolutionarily constrained sites (dC) was elevated in the case of comparisons involving closely located populations. Conversely diversity was reduced in the case of comparisons between distantly located populations. These patterns were confirmed by a significant negative relationship between the ratio of dC/dS and geographic distance. The observed high dC/dS could be explained by the abundance of deleterious mutations in comparisons involving closely located populations, due to the recent population divergence times. Since distantly related populations were separated over long periods of time, deleterious mutations might have been removed by purifying selection.

Sampling strategies for frequency spectrum-based population genomic inference

Background

The allele frequency spectrum (AFS) consists of counts of the number of single nucleotide polymorphism (SNP) loci with derived variants present at each given frequency in a sample. Multiple approaches have recently been developed for parameter estimation and calculation of model likelihoods based on the joint AFS from two or more populations. We conducted a simulation study of one of these approaches, implemented in the Python module δaδi, to compare parameter estimation and model selection accuracy given different sample sizes under one- and two-population models.

Results

Our simulations included a variety of demographic models and two parameterizations that differed in the timing of events (divergence or size change). Using a number of SNPs reasonably obtained through next-generation sequencing approaches (10,000 - 50,000), accurate parameter estimates and model selection were possible for models with more ancient demographic events, even given relatively small numbers of sampled individuals. However, for recent events, larger numbers of individuals were required to achieve accuracy and precision in parameter estimates similar to that seen for models with older divergence or population size changes. We quantify i) the uncertainty in model selection, using tools from information theory, and ii) the accuracy and precision of parameter estimates, using the root mean squared error, as a function of the timing of demographic events, sample sizes used in the analysis, and complexity of the simulated models.

Conclusions

Here, we illustrate the utility of the genome-wide AFS for estimating demographic history and provide recommendations to guide sampling in population genomics studies that seek to draw inference from the AFS. Our results indicate that larger samples of individuals (and thus larger AFS) provide greater power for model selection and parameter estimation for more recent demographic events.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-014-0254-4) contains supplementary material, which is available to authorized users.

Characterizing the genetic differences between two distinct migrant groups from Indo-European and Dravidian speaking populations in India

Background

India is home to many ethnically and linguistically diverse populations. It is hypothesized that history of invasions by people from Persia and Central Asia, who are referred as Aryans in Hindu Holy Scriptures, had a defining role in shaping the Indian population canvas. A shift in spoken languages from Dravidian languages to Indo-European languages around 1500 B.C. is central to the Aryan Invasion Theory. Here we investigate the genetic differences between two sub-populations of India consisting of: (1) The Indo-European language speaking Gujarati Indians with genome-wide data from the International HapMap Project; and (2) the Dravidian language speaking Tamil Indians with genome-wide data from the Singapore Genome Variation Project.

Results

We implemented three population genetics measures to identify genomic regions that are significantly differentiated between the two Indian populations originating from the north and south of India. These measures singled out genomic regions with: (i) SNPs exhibiting significant variation in allele frequencies in the two Indian populations; and (ii) differential signals of positive natural selection as quantified by the integrated haplotype score (iHS) and cross-population extended haplotype homozygosity (XP-EHH). One of the regions that emerged spans the SLC24A5 gene that has been functionally shown to affect skin pigmentation, with a higher degree of genetic sharing between Gujarati Indians and Europeans.

Conclusions

Our finding points to a gene-flow from Europe to north India that provides an explanation for the lighter skin tones present in North Indians in comparison to South Indians.

Population and genomic lessons from genetic analysis of two Indian populations

Indian demographic history includes special features such as founder effects, interpopulation segregation, complex social structure with a caste system and elevated frequency of consanguineous marriages. It also presents a higher frequency for some rare mendelian disorders and in the last two decades increased prevalence of some complex disorders. Despite the fact that India represents about one-sixth of the human population, deep genetic studies from this terrain have been scarce. In this study, we analyzed high-density genotyping and whole-exome sequencing data of a North and a South Indian population. Indian populations show higher differentiation levels than those reported between populations of other continents. In this work, we have analyzed its consequences, by specifically assessing the transferability of genetic markers from or to Indian populations. We show that there is limited genetic marker portability from available genetic resources such as HapMap or the 1,000 Genomes Project to Indian populations, which also present an excess of private rare variants. Conversely, tagSNPs show a high level of portability between the two Indian populations, in contrast to the common belief that North and South Indian populations are genetically very different. By estimating kinship from mates and consanguinity in our data from trios, we also describe different patterns of assortative mating and inbreeding in the two populations, in agreement with distinct mating preferences and social structures. In addition, this analysis has allowed us to describe genomic regions under recent adaptive selection, indicating differential adaptive histories for North and South Indian populations. Our findings highlight the importance of considering demography for design and analysis of genetic studies, as well as the need for extending human genetic variation catalogs to new populations and particularly to those with particular demographic histories.

MTHFR gene polymorphism in acute lymphoblastic leukemia among North Indian children: a case-control study and meta-analysis updated from 2011

Studies on the association of methylenetetrahydrofolate reductase (MTHFR) genotype in childhood acute lymphoblastic leukemia (ALL) have yielded conflicting results. The present study examines this association in north Indian children with ALL and includes an updated meta-analysis. MTHFR (677 and 1298) genotype of children with ALL and healthy adult controls were done by the PCR-restriction fragment length polymorphism (PCR-RFLP) method and were compared using various models of inheritance. A total of 150 patients and 300 controls were included. The 677T allele was found protective (odds ratio (OR) 0.21, 95% confidence interval (CI) 0.04-0.94), whereas 1298C allele led to an increase in risk (OR 4.44, 95% CI 2.19-8.99) of childhood ALL. Meta-analysis included 31 and 27 studies examining the association of 677 and 1298 genotypes, respectively. The 677 C -> T polymorphism was protective (OR 0.90, 95% CI 0.82-0.99). Protection was more pronounced in folate-sufficient populations as compared with those not covered by folate fortification guidelines. The 1298A->C polymorphism was associated with a marginal increase in risk (OR 1.19, 95% CI 1.01-1.40).

Forensic and population genetic analyses of eighteen non-CODIS miniSTR loci in the Korean population

We analyzed the variation of eighteen miniSTR loci in 411 randomly chosen individuals from Korea to increase the probability that a degraded sample can be typed, as well as to provide an expanded and reliable population database. Six multiplex PCR systems were developed (multiplex I: D1S1677, D2S441 and D4S2364; multiplex II: D10S1248, D14S1434 and D22S1045; multiplex III: D12S391, D16S3253 and D20S161; multiplex IV: D3S4529, D8S1115 and D18S853; multiplex V: D6S1017, D11S4463 and D17S1301; multiplex VI: D5S2500, D9S1122 and D21S1437). Allele frequencies and forensic parameters were calculated to evaluate the suitability and robustness of these non-CODIS miniSTR systems. No significant deviation from Hardy-Weinberg equilibrium expectations were observed, except for D4S2364, D5S2500 and D20S161 loci. A multidimensional scaling plot based on allele frequencies of the six miniSTR loci (D1S1677, D2S441, D4S2364, D10S1248, D14S1434 and D22S1045) showed that Koreans appeared to have most genetic affinity with Chinese and Japanese than to other Eurasian populations compared here. The combined probability of match calculated from the 18 miniSTR loci was 2.902 × 10(-17), indicating a high degree of polymorphism. Thus, the 18 miniSTR loci can be suitable for recovering useful information for analyzing degraded forensic casework samples and for adding supplementary genetic information for a variety of analyses involving closely related individuals where there is a need for additional genetic information.

NAT2 sequence polymorphisms and acetylation profiles in Indians

BACKGROUND

NAT2, a broad-spectrum drug-metabolizing gene, is of high pharmacogenetic interest. Based on seven different mutations in the NAT2 gene, an individual can either be categorized as a slow or fast acetylator.

MATERIALS & METHODS

In order to characterize acetylation profiles of Indians, where data are poorly available, we sequenced the 873 bp NAT2 coding region in 250 Indians, covering the whole of India including three tribes.

RESULTS

Altogether, 35 NAT2 alleles forming two acetylator phenotypes (distributed almost in equal proportion in India) were found; while the alleles determining slow acetylators were highly differentiated, the fast acetylator alleles were less in number but highly frequent.

CONCLUSION

Interestingly, distribution of two different acetylation phenotypes correlated well with historical dietary pattern in India. The neighbor-joining phylogenetic tree based on NAT2 gene polymorphisms in worldwide humans revealed genetic affinities among populations with similar acetylation phenotypes, which also placed Indians and Africans together in a single cluster.

HGDP and HapMap analysis by Ancestry Mapper reveals local and global population relationships

Knowledge of human origins, migrations, and expansions is greatly enhanced by the availability of large datasets of genetic information from different populations and by the development of bioinformatic tools used to analyze the data. We present Ancestry Mapper, which we believe improves on existing methods, for the assignment of genetic ancestry to an individual and to study the relationships between local and global populations. The principle function of the method, named Ancestry Mapper, is to give each individual analyzed a genetic identifier, made up of just 51 genetic coordinates, that corresponds to its relationship to the HGDP reference population. As a consequence, the Ancestry Mapper Id (AMid) has intrinsic biological meaning and provides a tool to measure similarity between world populations. We applied Ancestry Mapper to a dataset comprised of the HGDP and HapMap data. The results show distinctions at the continental level, while simultaneously giving details at the population level. We clustered AMids of HGDP/HapMap and observe a recapitulation of human migrations: for a small number of clusters, individuals are grouped according to continental origins; for a larger number of clusters, regional and population distinctions are evident. Calculating distances between AMids allows us to infer ancestry. The number of coordinates is expandable, increasing the power of Ancestry Mapper. An R package called Ancestry Mapper is available to apply this method to any high density genomic data set.

Genetic diversity and admixture patterns in Indian populations

India is a diverse land whose population holds the history of waves of human dispersal. Recent studies suggest two major ancestral contributions to most of the Indian sub-populations. However, present day Indians are thought to contain huge genetic diversity derived consequent to multiple cultural, linguistic and geographical variations. Genome-wide survey of individuals from current North (N-In) and South (S-In) India along with populations from HapMap Phase III and Indian sub-populations from HUGO Pan-Asian SNP Consortium is performed. Multivariate analysis (MDS and PCA) was carried out after merging data from the current study and other consortia. Indian sub-populations clustered separately from populations of major global geographical regions in MDS and PCA in a loose agglomeration except for two Indian subpopulations clustering near far eastern populations. F(st) values indicated diversity among Indian sub-populations which was substantiated by STRUCTURE analysis suggesting the possibility of additional admixture events.

Demographic inference using spectral methods on SNP data, with an analysis of the human out-of-Africa expansion

We present an implementation of a recently introduced method for estimating the allele-frequency spectrum under the diffusion approximation. For single-nucleotide polymorphism (SNP) frequency data from multiple populations, the method computes numerical solutions to the allele-frequency spectrum (AFS) under a complex model that includes population splitting events, migration, population expansion, and admixture. The solution to the diffusion partial differential equation (PDE) that mimics the evolutionary process is found by means of truncated polynomial expansions. In the absence of gene flow, our computation of frequency spectra yields exact results. The results are compared to those that use a finite-difference method and to forward diffusion simulations. In general, all the methods yield comparable results, although the polynomial-based approach is the most accurate in the weak-migration limit. Also, the economical use of memory attained by the polynomial expansions makes the study of models with four populations possible for the first time. The method was applied to a four-population model of the human expansion out of Africa and the peopling of the Americas, using the Environmental Genome Project (EGP) SNP database. Although our confidence intervals largely overlapped previous analyses of these data, some were significantly different. In particular, estimates of migration among African, European, and Asian populations were considerably lower than those in a previous study and the estimated time of migration out of Africa was earlier. The estimated time of founding of a human population outside of Africa was 52,000 years (95% confidence interval: 36,000–80,800 years).

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.