Segmental duplications and evolutionary acquisition of UV damage response in the SPATA31 gene family of primates and humans.

BACKGROUND: Segmental duplications are an abundant source for novel gene functions and evolutionary adaptations. This mechanism of generating novelty was very active during the evolution of primates particularly in the human lineage. Here, we characterize the evolution and function of the SPATA31 gene family (former designation FAM75A), which was previously shown to be among the gene families with the strongest signal of positive selection in hominoids. The mouse homologue for this gene family is a single copy gene expressed during spermatogenesis. RESULTS: We show that in primates, the SPATA31 gene duplicated into SPATA31A and SPATA31C types and broadened the expression into many tissues. Each type became further segmentally duplicated in the line towards humans with the largest number of full-length copies found for SPATA31A in humans. Copy number estimates of SPATA31A based on digital PCR show an average of 7.5 with a range of 5-11 copies per diploid genome among human individuals. The primate SPATA31 genes also acquired new protein domains that suggest an involvement in UV response and DNA repair. We generated antibodies and show that the protein is re-localized from the nucleolus to the whole nucleus upon UV-irradiation suggesting a UV damage response. We used CRISPR/Cas mediated mutagenesis to knockout copies of the gene in human primary fibroblast cells. We find that cell lines with reduced functional copies as well as naturally occurring low copy number HFF cells show enhanced sensitivity towards UV-irradiation. CONCLUSION: The acquisition of new SPATA31 protein functions and its broadening of expression may be related to the evolution of the diurnal life style in primates that required a higher UV tolerance. The increased segmental duplications in hominoids as well as its fast evolution suggest the acquisition of further specific functions particularly in humans.

Geographic distribution and adaptive significance of genomic structural variants: an anthropological genetics perspective.

Anthropological geneticists have successfully used single-nucleotide and short tandem repeat variations across human genomes to reconstruct human history. These markers have also been used extensively to identify adaptive and phenotypic variation. The recent advent of high-throughput genomic technologies revealed an overlooked type of genomic variation: structural variants (SVs). In fact, some SVs may contribute to human adaptation in substantial and previously unexplored ways. SVs include deletions, insertions, duplications, inversions, and translocations of genomic segments that vary among individuals from the same species. SVs are much less numerous than single-nucleotide variants but account for at least seven times more variable base pairs than do single-nucleotide variants when two human genomes are compared. Moreover, recent studies have shown that SVs have higher mutation rates than single-nucleotide variants when the affected base pairs are considered, especially in certain parts of the genome. The null hypothesis for the evolution of SVs, as for single-nucleotide variants, is neutrality. Hence, drift is the primary force that shapes the current allelic distribution of most SVs. However, due to their size, a larger proportion of SVs appear to evolve under nonneutral forces (mostly purifying selection) than do single-nucleotide variants. In fact, as exemplified by several groundbreaking studies, SVs contribute to anthropologically relevant phenotypic variation and local adaptation among humans. In this review, we argue that with the advent of affordable genomic technologies, anthropological scrutiny of genomic structural variation emerges as a fertile area of inquiry to better understand human phenotypic variation. To motivate potential studies, we discuss scenarios through which structural variants (SVs) affect phenotypic variation among humans within an anthropological context. We further provide a methodological workflow in which we analyzed 1000 Genomes deletion variants and identified 16 exonic deletions that are specific to the African continent. We analyzed two of these deletion variants affecting the keratin-associated protein (KAP) cluster in a locus-specific manner. Our analysis revealed that these deletions may indeed affect phenotype and likely evolved under geography-specific positive selection. We outline all the major software and data sets for these analyses and provide the basic R and Perl codes we used for this example workflow analysis. Overall, we hope that this review will encourage and facilitate incorporation of genomic structural variation in anthropological research programs.

Genome-wide gene-environment study identifies glutamate receptor gene GRIN2A as a Parkinson's disease modifier gene via interaction with coffee.

Our aim was to identify genes that influence the inverse association of coffee with the risk of developing Parkinson's disease (PD). We used genome-wide genotype data and lifetime caffeinated-coffee-consumption data on 1,458 persons with PD and 931 without PD from the NeuroGenetics Research Consortium (NGRC), and we performed a genome-wide association and interaction study (GWAIS), testing each SNP's main-effect plus its interaction with coffee, adjusting for sex, age, and two principal components. We then stratified subjects as heavy or light coffee-drinkers and performed genome-wide association study (GWAS) in each group. We replicated the most significant SNP. Finally, we imputed the NGRC dataset, increasing genomic coverage to examine the region of interest in detail. The primary analyses (GWAIS, GWAS, Replication) were performed using genotyped data. In GWAIS, the most significant signal came from rs4998386 and the neighboring SNPs in GRIN2A. GRIN2A encodes an NMDA-glutamate-receptor subunit and regulates excitatory neurotransmission in the brain. Achieving P(2df) = 10(-6), GRIN2A surpassed all known PD susceptibility genes in significance in the GWAIS. In stratified GWAS, the GRIN2A signal was present in heavy coffee-drinkers (OR = 0.43; P = 6×10(-7)) but not in light coffee-drinkers. The a priori Replication hypothesis that "Among heavy coffee-drinkers, rs4998386_T carriers have lower PD risk than rs4998386_CC carriers" was confirmed: OR(Replication) = 0.59, P(Replication) = 10(-3); OR(Pooled) = 0.51, P(Pooled) = 7×10(-8). Compared to light coffee-drinkers with rs4998386_CC genotype, heavy coffee-drinkers with rs4998386_CC genotype had 18% lower risk (P = 3×10(-3)), whereas heavy coffee-drinkers with rs4998386_TC genotype had 59% lower risk (P = 6×10(-13)). Imputation revealed a block of SNPs that achieved P(2df)<5×10(-8) in GWAIS, and OR = 0.41, P = 3×10(-8) in heavy coffee-drinkers. This study is proof of concept that inclusion of environmental factors can help identify genes that are missed in GWAS. Both adenosine antagonists (caffeine-like) and glutamate antagonists (GRIN2A-related) are being tested in clinical trials for treatment of PD. GRIN2A may be a useful pharmacogenetic marker for subdividing individuals in clinical trials to determine which medications might work best for which patients.

Uncovering the genetic architecture and evolutionary roots of androgenetic alopecia in African men.

Androgenetic alopecia is a highly heritable trait. However, much of our understanding about the genetics of male pattern baldness comes from individuals of European descent. Here, we examined a novel dataset comprising 2,136 men from Ghana, Nigeria, Senegal, and South Africa that were genotyped using a custom array. We first tested how genetic predictions of baldness generalize from Europe to Africa, finding that polygenic scores from European GWAS yielded AUC statistics that ranged from 0.513 to 0.546, indicating that genetic predictions of baldness in African populations performed notably worse than in European populations. Subsequently, we conducted the first African GWAS of androgenetic alopecia, focusing on self-reported baldness patterns at age 45. After correcting for present age, population structure, and study site, we identified 266 moderately significant associations, 51 of which were independent (p-value < 10-5, r2 < 0.2). Most baldness associations were autosomal, and the X chromosomes does not appear to have a large impact on baldness in African men. Finally, we examined the evolutionary causes of continental differences in genetic architecture. Although Neanderthal alleles have previously been associated with skin and hair phenotypes, we did not find evidence that European-ascertained baldness hits were enriched for signatures of ancient introgression. Most loci that are associated with androgenetic alopecia are evolving neutrally. However, multiple baldness-associated SNPs near the EDA2R and AR genes have large allele frequency differences between continents. Collectively, our findings illustrate how evolutionary history contributes to the limited portability of genetic predictions across ancestries.

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

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

Unraveling a fine-scale high genetic heterogeneity and recent continental connections of an Arabian Peninsula population.

Recent studies have showed the diverse genetic architecture of the highly consanguineous populations inhabiting the Arabian Peninsula. Consanguinity coupled with heterogeneity is complex and makes it difficult to understand the bases of population-specific genetic diseases in the region. Therefore, comprehensive genetic characterization of the populations at the finest scale is warranted. Here, we revisit the genetic structure of the Kuwait population by analyzing genome-wide single nucleotide polymorphisms data from 583 Kuwaiti individuals sorted into three subgroups. We envisage a diverse demographic genetic history among the three subgroups based on drift and allelic sharing with modern and ancient individuals. Furthermore, our comprehensive haplotype-based analyses disclose a high genetic heterogeneity among the Kuwaiti populations. We infer the major sources of ancestry within the newly defined groups; one with an obvious predominance of sub-Saharan/Western Africa mostly comprising Kuwait-B individuals, and other with West Eurasia including Kuwait-P and Kuwait-S individuals. Overall, our results recapitulate the historical population movements and reaffirm the genetic imprints of the legacy of continental trading in the region. Such deciphering of fine-scale population structure and their regional genetic heterogeneity would provide clues to the uncharted areas of disease-gene discovery and related associations in populations inhabiting the Arabian Peninsula.

Mitochondrial haplogroup J associated with higher risk of obesity in the Qatari population.

Obesity, a major risk factor for metabolic disorders, is highly prevalent in Qatari population. Maternal transmission of obesity traits can be significant; for example, X haplogroup is known to be associated with lower BMI and body fat mass in Northern Europeans and T haplogroup which is a sister haplogroup of J is known to be associated with obesity in Caucasian subjects from Austria and Southern Italy. We aimed to delineate the mitochondrial haplogroups and variants associated with obesity in Qatari population. Mitochondrial genomes of 864 Qatari individuals were extracted from whole exome sequencing data with an average coverage of 77X. We distributed the participants into 2 sub-cohorts: obese (BMI ≥ 30) and non-obese (BMI < 30); the mean value of BMI from these two groups were 36.5 ± 5.7 and 26.5 ± 2.6, respectively. Mitochondrial haplogroup profiling followed by uni- and multivariant association tests adjusted for covariates were performed. Qatari individuals with mitochondrial haplogroup J had an increased (twofold) risk of obesity (odds ratio [OR] 1.925; 95% CI 1.234-3.002; P = 0.0038; the Bonferroni adjusted P value threshold is 0.0041), whereas the individuals with haplogroup X were at low risk of obesity (OR 0.387; 95% CI 0.175-0.857; P = 0.019). Further, a set of 38 mitochondrial variants were found to be associated (at P ≤ 0.05) with obesity in models adjusted for age, sex and haplogroup.

Salivary Microbiome Diversity in Kuwaiti Adolescents with Varied Body Mass Index-A Pilot Study.

The potential role of the salivary microbiome in human diseases has increasingly been explored. The salivary microbiome has been characterized in several global populations, except the Arabian Gulf region. Hence, in this pilot study, we profiled the salivary microbiome of Kuwaiti adolescents with varied body mass indexes (BMI). The analyses of core microbiome composition showed Firmicutes, Bacteroidota, Proteobacteria, Patescibacteria, Fusobacteriota, Actinobacteriota, and Campylobacterota as the common phylum found in the Kuwaiti adolescent population. We also illustrated a diverse microbial community among the sampled individuals grouped according to their BMI. Notably, the overweight group was found with a higher number of distinct taxa than other groups. As such, the core microbiome composition was found to be significantly different (p-value < 0.001) across different BMI groups. Overall, this pilot investigation outlined the microbial diversity and suggested that changes in salivary microbiome composition in people with obese or overweight BMI might reflect their susceptibility to oral diseases.

ACE2 and FURIN variants are potential predictors of SARS-CoV-2 outcome: A time to implement precision medicine against COVID-19.

The severity of the new COVID-19 pandemic caused by the SARS-CoV-2 virus is strikingly variable in different global populations. SARS-CoV-2 uses ACE2 as a cell receptor, TMPRSS2 protease, and FURIN peptidase to invade human cells. Here, we investigated 1,378 whole-exome sequences of individuals from the Middle Eastern populations (Kuwait, Qatar, and Iran) to explore natural variations in the ACE2, TMPRSS2, and FURIN genes. We identified two activating variants (K26R and N720D) in the ACE2 gene that are more common in Europeans than in the Middle Eastern, East Asian, and African populations. We postulate that K26R can activate ACE2 and facilitate binding to S-protein RBD while N720D enhances TMPRSS2 cutting and, ultimately, viral entry. We also detected deleterious variants in FURIN that are frequent in the Middle Eastern but not in the European populations. This study highlights specific genetic variations in the ACE2 and FURIN genes that may explain SARS-CoV-2 clinical disparity. We showed structural evidence of the functionality of these activating variants that increase the SARS-CoV-2 aggressiveness. Finally, our data illustrate a significant correlation between ACE2 variants identified in people from Middle Eastern origins that can be further explored to explain the variation in COVID-19 infection and mortality rates globally.

Delineation of Mitochondrial DNA Variants From Exome Sequencing Data and Association of Haplogroups With Obesity in Kuwait.

BACKGROUND/OBJECTIVES: Whole-exome sequencing is a valuable tool to determine genetic variations that are associated with rare and common health conditions. A limited number of studies demonstrated that mitochondrial DNA can be captured using whole-exome sequencing. Previous studies have suggested that mitochondrial DNA variants and haplogroup lineages are associated with obesity. Therefore, we investigated the role of mitochondrial variants and haplogroups contributing to the risk of obesity in Arabs in Kuwait using exome sequencing data. SUBJECTS/METHODS: Indirect mitochondrial genomes were extracted from exome sequencing data from 288 unrelated native Arab individuals from Kuwait. The cohort was divided into obese [body mass index (BMI) ≥ 30 kg/m2] and non-obese (BMI < 30 kg/m2) groups. Mitochondrial variants were identified, and haplogroups were classified and compared with other sequencing technologies. Statistical analysis was performed to determine associations and identify mitochondrial variants and haplogroups affecting obesity. RESULTS: Haplogroup R showed a protective effect on obesity [odds ratio (OR) = 0.311; P = 0.006], whereas haplogroup L individuals were at high risk of obesity (OR = 2.285; P = 0.046). Significant differences in mitochondrial variants between the obese and non-obese groups were mainly haplogroup-defining mutations and were involved in processes in energy generation. The majority of mitochondrial variants and haplogroups extracted from exome were in agreement with technical replica from Sanger and whole-genome sequencing. CONCLUSIONS: This is the first to utilize whole-exome data to extract entire mitochondrial haplogroups to study its association with obesity in an Arab population.

Could the D614G substitution in the SARS-CoV-2 spike (S) protein be associated with higher COVID-19 mortality?

The increasing number of deaths due to the COVID-19 pandemic has raised serious global concerns. Increased testing capacity and ample intensive care availability could explain lower mortality in some countries compared to others. Nevertheless, it is also plausible that the SARS-CoV-2 mutations giving rise to different phylogenetic clades are responsible for the apparent death rate disparities around the world. Current research literature linking the genetic make-up of SARS-CoV-2 with fatalities is lacking. Here, we suggest that this disparity in fatality rates may be attributed to SARS-CoV-2 evolving mutations and urge the international community to begin addressing the phylogenetic clade classification of SARS-CoV-2 in relation to clinical outcomes.

Genome-Wide Selection Scan in an Arabian Peninsula Population Identifies a TNKS Haplotype Linked to Metabolic Traits and Hypertension.

Despite the extreme and varying environmental conditions prevalent in the Arabian Peninsula, it has experienced several waves of human migrations following the out-of-Africa diaspora. Eventually, the inhabitants of the peninsula region adapted to the hot and dry environment. The adaptation and natural selection that shaped the extant human populations of the Arabian Peninsula region have been scarcely studied. In an attempt to explore natural selection in the region, we analyzed 662,750 variants in 583 Kuwaiti individuals. We searched for regions in the genome that display signatures of positive selection in the Kuwaiti population using an integrative approach in a conservative manner. We highlight a haplotype overlapping TNKS that showed strong signals of positive selection based on the results of the multiple selection tests conducted (integrated Haplotype Score, Cross Population Extended Haplotype Homozygosity, Population Branch Statistics, and log-likelihood ratio scores). Notably, the TNKS haplotype under selection potentially conferred a fitness advantage to the Kuwaiti ancestors for surviving in the harsh environment while posing a major health risk to present-day Kuwaitis.

Diverse genetic origin of Indian Muslims: evidence from autosomal STR loci.

The origin and relationships of Indian Muslims is still dubious and are not yet genetically well studied. In the light of historically attested movements into Indian subcontinent during the demic expansion of Islam, the present study aims to substantiate whether it had been accompanied by any gene flow or only a cultural transformation phenomenon. An array of 13 autosomal STR markers that are common in the worldwide data sets was used to explore the genetic diversity of Indian Muslims. The austere endogamy being practiced for several generations was confirmed by the genetic demarcation of each of the six Indian Muslim communities in the phylogenetic assessments for the markers examined. The analyses were further refined by comparison with geographically closest neighboring Hindu religious groups (including several caste and tribal populations) and the populations from Middle East, East Asia and Europe. We found that some of the Muslim populations displayed high level of regional genetic affinity rather than religious affinity. Interestingly, in Dawoodi Bohras (TN and GUJ) and Iranian Shia significant genetic contribution from West Asia, especially Iran (49, 47 and 46%, respectively) was observed. This divulges the existence of Middle Eastern genetic signatures in some of the contemporary Indian Muslim populations. Our study reveals that the spread of Islamic faith in the Indian subcontinent was predominantly cultural transformation associated with minor gene flow from West Asia.

Microsatellite diversity delineates genetic relationships of Shia and Sunni Muslim populations of Uttar Pradesh, India.

In this study we characterize the genetic diversity and relationships between the Shia and Sunni Muslim populations of North India and geographically targeted neighboring and global populations. We examined a number of parameters of population genetic and forensic interest based on the allele frequencies from 15 autosomal STR loci (D8S1179, D21S11, D7S820, CSF1PO, D19S433, VWA, TPOX, D18S51, D3S1358, THO1, D13S317, D16S539, D2S1338, D5S818, and FGA). All the studied loci were consistent with Hardy-Weinberg equilibrium, except loci D18S51 and FGA for both Muslim populations, even after applying the Bonferroni correction. The combined power of exclusion and combined power of discrimination values for all 15 STR loci were 0.9999 and >0.99999, respectively, in both Muslim populations. Gene diversity values ranged from 0.6784 (TPOX) to 0.9027 (FGA) for Shia Muslims and from 0.7152 (CSF1PO) to 0.9120 (D18S51) for Sunni Muslims. The observed heterozygosity (H(o)) ranged from 0.5833 (D18S51) to 0.8595 (VWA) in Shia Muslims and from 0.6818 (CSF1PO) to 0.8333 (D21S11) in Sunni Muslims and was lower than the expected heterozygosity (H(e)) for 11 out of the 15 STRs typed. We analyzed the genetic affinities of the Shia and Sunni Muslim populations with their geographically closest neighboring North Indian, Middle Eastern, East Asian, and European populations using distance-based methods, including neighbor-joining trees and multidimensional scaling. In addition, we estimated the genetic contribution of the putative parental populations included in the analysis to the Shia and Sunni Muslim gene pool using admixture analysis. Although we observed a certain degree of genetic contribution from Iran to both Muslim populations, the results of the phylogenetic analyses based on autosomal STRs suggest genetic relatedness with some of the geographically closest neighboring Hindu religious populations.

Mitochondrial DNA D-loop sequencing reveals obesity variants in an Arab population.

Background: The association of mitochondrial DNA (mtDNA) variations with obesity has been investigated in diverse populations across the world. However, such obesity-associated mtDNA examinations are rarely conducted in Arab populations. Materials and methods: We re-sequenced mtDNA displacement loop (D-loop) region of 395 Arab individuals of Kuwait. We categorized the individuals based on their BMI scores as obese (n=232; BMI ≥30 kg/m2), overweight (n=110; BMI ≥25 kg/m2 and <30 kg/m2), and lean (n=53; BMI <25 kg/m2). We performed all the statistical tests by combining obese and overweight individuals in one group. Association analyses were conducted applying Fisher's exact test and logistic regression model. Results: We identified that the mtDNA variations m.73A>G, and m.523delAC were positively correlated with obesity, while m.310T>C, and m.16318A>T were negatively associated. All these variants, except m.16318A>T, remain statistically significant after adjusting for age and gender. We found that the variant m.73A>G increases the likelihood of being obese by 6-fold, whereas haplogroup H decreases the probability of being obese in Arab individuals of Kuwait. Haplotype analysis revealed that a haplotype, A263G-C309CT-T310C, defining the H2a clade of H haplogroup, reduces the probability of being obese. Conclusion: Our study reports, for the first time, the obesity-related mtDNA variants in Arabs of Kuwait. Based on the mtDNA D-loop region variations, we detected particular variants and haplogroup that are related with increased and decreased probability of being obese in the Kuwait Arab population.

Assessment of coding region variants in Kuwaiti population: implications for medical genetics and population genomics.

Consanguineous populations of the Arabian Peninsula have been underrepresented in global efforts that catalogue human exome variability. We sequenced 291 whole exomes of unrelated, healthy native Arab individuals from Kuwait to a median coverage of 45X and characterised 170,508 single-nucleotide variants (SNVs), of which 21.7% were 'personal'. Up to 12% of the SNVs were novel and 36% were population-specific. Half of the SNVs were rare and 54% were missense variants. The study complemented the Greater Middle East Variome by way of reporting many additional Arabian exome variants. The study corroborated Kuwaiti population genetic substructures previously derived using genome-wide genotype data and illustrated the genetic relatedness among Kuwaiti population subgroups, Middle Eastern, European and Ashkenazi Jewish populations. The study mapped 112 rare and frequent functional variants relating to pharmacogenomics and disorders (recessive and common) to the phenotypic characteristics of Arab population. Comparative allele frequency data and carrier distributions of known Arab mutations for 23 disorders seen among Arabs, of putative OMIM-listed causal mutations for 12 disorders observed among Arabs but not yet characterized for genetic basis in Arabs, and of 17 additional putative mutations for disorders characterized for genetic basis in Arab populations are presented for testing in future Arab studies.

Genetic variants associated with warfarin dosage in Kuwaiti population.

Assessing the distinct prevalence or absence of genetic variants associated with differential response to the anticoagulant medication of warfarin in different population groups is actively pursued by pharmacogenomics community. Populations from Arabian Peninsula are underrepresented in such studies. By way of examining exome- and genome-wide genotype data from 1395 Arab individuals in Kuwait, we report distinct occurrence of warfarin response-related variants rs12460590_A/CYP2A7, rs2108622_T/CYP4F2, rs2884737_C/VKORC1 and distinct absence of rs11150606_C/PRSS53 in Kuwaiti population. The presented results in conjunction with similar literature reports on Qatari population enhance the worldwide understanding on population-specific distributions of genetic variants associated with warfarin drug dosage.

Atopic Dermatitis Susceptibility Variants in Filaggrin Hitchhike Hornerin Selective Sweep.

Human skin has evolved rapidly, leaving evolutionary signatures in the genome. The filaggrin (FLG) gene is widely studied for its skin-barrier function in humans. The extensive genetic variation in this gene, especially common loss-of-function (LoF) mutations, has been established as primary risk factors for atopic dermatitis. To investigate the evolution of this gene, we analyzed 2,504 human genomes and genotyped the copy number variation of filaggrin repeats within FLG in 126 individuals from diverse ancestral backgrounds. We were unable to replicate a recent study claiming that LoF of FLG is adaptive in northern latitudes with lower ultraviolet light exposure. Instead, we present multiple lines of evidence suggesting that FLG genetic variation, including LoF variants, have little or no effect on fitness in modern humans. Haplotype-level scrutinization of the locus revealed signatures of a recent selective sweep in Asia, which increased the allele frequency of a haplotype group (Huxian haplogroup) in Asian populations. Functionally, we found that the Huxian haplogroup carries dozens of functional variants in FLG and hornerin (HRNR) genes, including those that are associated with atopic dermatitis susceptibility, HRNR expression levels and microbiome diversity on the skin. Our results suggest that the target of the adaptive sweep is HRNR gene function, and the functional FLG variants that involve susceptibility to atopic dermatitis, seem to hitchhike the selective sweep on HRNR. Our study presents a novel case of a locus that harbors clinically relevant common genetic variation with complex evolutionary trajectories.

Genetic affinities of the central Indian tribal populations.

BACKGROUND: The central Indian state Madhya Pradesh is often called as 'heart of India' and has always been an important region functioning as a trinexus belt for three major language families (Indo-European, Dravidian and Austroasiatic). There are less detailed genetic studies on the populations inhabited in this region. Therefore, this study is an attempt for extensive characterization of genetic ancestries of three tribal populations, namely; Bharia, Bhil and Sahariya, inhabiting this region using haploid and diploid DNA markers. METHODOLOGY/PRINCIPAL FINDINGS: Mitochondrial DNA analysis showed high diversity, including some of the older sublineages of M haplogroup and prominent R lineages in all the three tribes. Y-chromosomal biallelic markers revealed high frequency of Austroasiatic-specific M95-O2a haplogroup in Bharia and Sahariya, M82-H1a in Bhil and M17-R1a in Bhil and Sahariya. The results obtained by haploid as well as diploid genetic markers revealed strong genetic affinity of Bharia (a Dravidian speaking tribe) with the Austroasiatic (Munda) group. The gene flow from Austroasiatic group is further confirmed by their Y-STRs haplotype sharing analysis, where we determined their founder haplotype from the North Munda speaking tribe, while, autosomal analysis was largely in concordant with the haploid DNA results. CONCLUSIONS/SIGNIFICANCE: Bhil exhibited largely Indo-European specific ancestry, while Sahariya and Bharia showed admixed genetic package of Indo-European and Austroasiatic populations. Hence, in a landscape like India, linguistic label doesn't unequivocally follow the genetic footprints.