Global population genetics and diversity in the TAS2R bitter taste receptor family

Bitter taste receptors (TAS2Rs) are noted for their role in perception, and mounting evidence suggests that they mediate responses to compounds entering airways, gut, and other tissues. The importance of these roles suggests that TAS2Rs have been under pressure from natural selection. To determine the extent of variation in TAS2Rs on a global scale and its implications for human evolution and behavior, we analyzed patterns of diversity in the complete 25 gene repertoire of human TAS2Rs in ∼2,500 subjects representing worldwide populations. Across the TAS2R family as a whole, we observed 721 single nucleotide polymorphisms (SNPs) including 494 nonsynonymous SNPs along with 40 indels and gained and lost start and stop codons. In addition, computational predictions identified 169 variants particularly likely to affect receptor function, making them candidate sources of phenotypic variation. Diversity levels ranged widely among loci, with the number of segregating sites ranging from 17 to 41 with a mean of 32 among genes and per nucleotide heterozygosity (π) ranging from 0.02% to 0.36% with a mean of 0.12%. F_ST ranged from 0.01 to 0.26 with a mean of 0.13, pointing to modest differentiation among populations. Comparisons of observed π and F_ST values with their genome wide distributions revealed that most fell between the 5th and 95th percentiles and were thus consistent with expectations. Further, tests for natural selection using Tajima’s D statistic revealed only two loci departing from expectations given D’s genome wide distribution. These patterns are consistent with an overall relaxation of selective pressure on TAS2Rs in the course of recent human evolution.

Bitter taste receptors: Genes, evolution and health

Bitter taste perception plays vital roles in animal behavior and fitness. By signaling the presence of toxins in foods, particularly noxious defense compounds found in plants, it enables animals to avoid exposure. In vertebrates, bitter perception is initiated by TAS2Rs, a family of G protein-coupled receptors expressed on the surface of taste buds. There, oriented toward the interior of the mouth, they monitor the contents of foods, drinks and other substances as they are ingested. When bitter compounds are encountered, TAS2Rs respond by triggering neural pathways leading to sensation. The importance of this role placed TAS2Rs under selective pressures in the course of their evolution, leaving signatures in patterns of gene gain and loss, sequence polymorphism, and population structure consistent with vertebrates' diverse feeding ecologies. The protective value of bitter taste is reduced in modern humans because contemporary food supplies are safe and abundant. However, this is not always the case. Some crops, particularly in the developing world, retain surprisingly high toxicity and bitterness remains an important measure of safety. Bitter perception also shapes health through its influence on preference driven behaviors such as diet choice, alcohol intake and tobacco use. Further, allelic variation in TAS2Rs is extensive, leading to individual differences in taste sensitivity that drive these behaviors, shaping susceptibility to disease. Thus, bitter taste perception occupies a critical intersection between ancient evolutionary processes and modern human health.

Worldwide diversity, association potential, and natural selection in the superimposed taste genes, CD36 and GNAT3

CD36 and GNAT3 mediate taste responses, with CD36 acting as a lipid detector and GNAT3 acting as the α subunit of gustducin, a G protein governing sweet, savory, and bitter transduction. Strikingly, the genes encoding CD36 and GNAT3 are genomically superimposed, with CD36 completely encompassing GNAT3. To characterize genetic variation across the CD36-GNAT3 region, its implications for phenotypic diversity, and its recent evolution, we analyzed from ~2,500 worldwide subjects sequenced by the 1000 Genomes Project (1000GP). CD36-GNAT3 harbored extensive diversity including 8,688 single-nucleotide polymorphisms (SNPs), 414 indels, and other complex variants. Sliding window analyses revealed that nucleotide diversity and population differentiation across CD36-GNAT3 were consistent with genome-wide trends in the 1000GP (π = 0.10%, P = 0.64; FST = 9.0%, P = 0.57). In addition, functional predictions using SIFT and PolyPhen-2 identified 60 variants likely to alter protein function, and they were in weak linkage disequilibrium (r2 < 0.17), suggesting their effects are largely independent. However, the frequencies of predicted functional variants were low (P¯ = 0.0013), indicating their contributions to phenotypic variance on population scales are limited. Tests using Tajima's D statistic revealed that pressures from natural selection have been relaxed across most of CD36-GNAT3 during its recent history (0.39 < P < 0.67). However, CD36 exons showed signs of local adaptation consistent with prior reports (P < 0.035). Thus, CD36 and GNAT3 harbor numerous variants predicted to affect taste sensitivity, but most are rare and phenotypic variance on a population level is likely mediated by a small number of sites.

Bitter Fruit: Inverse Associations Between PTC and Antidesma bunius Perception

Ability to perceive the bitter compound phenylthiocarbamide (PTC) is inherited via a dominant "taster" allele of the TAS2R38 gene, whereas inability is inherited via a recessive "non-taster" allele. This raises a question: Is the non-taster allele functionless, or does it mediate perception of compounds other than PTC? New evidence supports speculation that it is indeed functional. Associations between TAS2R38 mutations and bitter sensitivity to the tropical berry Antidesma bunius are the inverse of those PTC, suggesting that the non-taster allele enables perception to compounds in the fruit.

Copy Number Variation in TAS2R Bitter Taste Receptor Genes: Structure, Origin, and Population Genetics

Bitter taste receptor genes (TAS2Rs) harbor extensive diversity, which is broadly distributed across human populations and strongly associated with taste response phenotypes. The majority of TAS2R variation is composed of single-nucleotide polymorphisms. However, 2 closely positioned loci at 12p13, TAS2R43 and -45, harbor high-frequency deletion (Δ) alleles in which genomic segments are absent, resulting in copy number variation (CNV). To resolve their chromosomal structure and organization, we generated maps using long-range contig alignments and local sequencing across the TAS2R43-45 region. These revealed that the deletion alleles (43Δ and 45Δ) are 37.8 and 32.2kb in length, respectively and span the complete coding region of each gene (~1kb) along with extensive up- and downstream flanking sequence, producing separate CNVs at the 2 loci. Comparisons with a chimpanzee genome, which contained intact homologs of TAS2R43, -45, and nearby TAS2Rs, indicated that the deletions evolved recently, through unequal recombination in a cluster of closely related loci. Population genetic analyses in 946 subjects from 52 worldwide populations revealed that copy number ranged from 0 to 2 at both TAS2R43 and TAS2R45, with 43Δ and 45Δ occurring at high global frequencies (0.33 and 0.18). Estimated recombination rates between the loci were low (ρ = 2.7×10(-4); r = 6.6×10(-9)) and linkage disequilibrium was high (D' = 1.0), consistent with their adjacent genomic positioning and recent origin. Geographic variation pointed to an African origin for the deletions. However, no signatures of natural selection were found in population structure or integrated haplotype scores spanning the region, suggesting that patterns of diversity at TAS2R43 and -45 are primarily due to genetic drift.

Receptor Polymorphism and Genomic Structure Interact to Shape Bitter Taste Perception

The ability to taste bitterness evolved to safeguard most animals, including humans, against potentially toxic substances, thereby leading to food rejection. Nonetheless, bitter perception is subject to individual variations due to the presence of genetic functional polymorphisms in bitter taste receptor (TAS2R) genes, such as the long-known association between genetic polymorphisms in TAS2R38 and bitter taste perception of phenylthiocarbamide. Yet, due to overlaps in specificities across receptors, such associations with a single TAS2R locus are uncommon. Therefore, to investigate more complex associations, we examined taste responses to six structurally diverse compounds (absinthin, amarogentin, cascarillin, grosheimin, quassin, and quinine) in a sample of the Caucasian population. By sequencing all bitter receptor loci, inferring long-range haplotypes, mapping their effects on phenotype variation, and characterizing functionally causal allelic variants, we deciphered at the molecular level how a subjects’ genotype for the whole-family of TAS2R genes shapes variation in bitter taste perception. Within each haplotype block implicated in phenotypic variation, we provided evidence for at least one locus harboring functional polymorphic alleles, e.g. one locus for sensitivity to amarogentin, one of the most bitter natural compounds known, and two loci for sensitivity to grosheimin, one of the bitter compounds of artichoke. Our analyses revealed also, besides simple associations, complex associations of bitterness sensitivity across TAS2R loci. Indeed, even if several putative loci harbored both high- and low-sensitivity alleles, phenotypic variation depended on linkage between these alleles. When sensitive alleles for bitter compounds were maintained in the same linkage phase, genetically driven perceptual differences were obvious, e.g. for grosheimin. On the contrary, when sensitive alleles were in opposite phase, only weak genotype-phenotype associations were seen, e.g. for absinthin, the bitter principle of the beverage absinth. These findings illustrate the extent to which genetic influences on taste are complex, yet arise from both receptor activation patterns and linkage structure among receptor genes.

Human bitter taste is believed to protect us from the ingestion of poisonous substances, thereby shaping food rejections. Bitter perception differs, however, across individuals, due to genetic variations in the ~25 bitter taste receptor (TAS2R) genes. A famous example known since the 1930s is the inherited bitter taste sensitivity to phenylthiocarbamide, which is associated with genetic polymorphisms in a single TAS2R gene. Yet, such simple receptor-substance associations do not reflect the full complexity of bitter perception, since individual bitter substances frequently activate several TAS2Rs. Here, we provide an in-depth analysis of the genetic variability influencing human bitter taste. While each study subject carried a different set of genetic polymorphisms, we found that most variations reside in just six blocks, each harboring only one to five haplotypes. Thus, besides simple associations between taste and TAS2R gene polymorphisms, we revealed complex associations dependent on linkage between several high- and low-sensitivity alleles. Indeed, subjects carried either sensitive or insensitive alleles for receptors sensitive to grosheimin, a bitter compound in artichoke, or at least one sensitive allele for receptors specific for absinthin, the bitter principle of absinth. In short, simple associations and complex genomic linkage determine sensitivity to selected dietary bitter compounds.

Association of a bitter taste receptor mutation with Balkan Endemic Nephropathy (BEN)

Background

Balkan Endemic Nephropathy (BEN) is late-onset kidney disease thought to arise from chronic exposure to aristolochic acid, a phytotoxin that contaminates wheat supplies in rural areas of Eastern Europe. It has recently been demonstrated that humans are capable of perceiving aristolochic acid at concentrations below 40 nM as the result of high-affinity interactions with the TAS2R43 bitter taste receptor. Further, TAS2R43 harbors high-frequency loss-of-function mutations resulting in 50-fold variability in perception. This suggests that genetic variation in TAS2R43 might affect susceptibility to BEN, with individuals carrying functional forms of the receptor being protected by an ability to detect tainted foods.

Methods

To determine whether genetic variation in TAS2R43 predicts BEN susceptibility, we examined genotype-phenotype associations in a case–control study. A cohort of 88 affected and 99 control subjects from western Bulgaria were genotyped with respect to two key missense variants and a polymorphic whole-gene deletion of TAS2R43 (W35S, H212R, and wt/Δ), which are known to affect taste sensitivity to aristolochic acid. Tests for association between haplotypes and BEN status were then performed.

Results

Three major TAS2R43 haplotypes observed in previous studies (TAS2R43-W35/H212, -S35/R212 and –Δ) were present at high frequencies (0.17, 0.36, and 0.47 respectively) in our sample, and a significant association between genotype and BEN status was present (P = 0.020; odds ratio 1.18). However, contrary to expectation, BEN was positively associated with TAS2R43-W35/H212, a highly responsive allele previously shown to confer elevated bitter sensitivity to aristolochic acid, which should drive aversion but might also affect absorption, altering toxin activation.

Conclusions

Our findings are at strong odds with the prediction that carriers of functional alleles of TAS2R43 are protected from BEN by an ability to detect and avoid aristolochic acid exposure. Evidence for a positive association between high-sensitivity alleles and BEN status suggests instead that possession of toxin-responsive receptor variants may paradoxically increase vulnerability, possibly by shifting attractive responses associated with low-intensity bitter sensations. The broad-spectrum tuning of the ~25-member TAS2R family as a whole toward xenobiotics points to a potentially far-reaching relevance of bitter responses to exposure-related disease in both individuals and populations.

Blindly using Wald's test can miss rare disease-causal variants in case-control association studies

There are four tests--the likelihood ratio (LR) test, Wald's test, the score test and the exact test--commonly employed in genetic association studies. On comparison of the four tests, we found that Wald's test, popular in genome-wide screens due to its low computational demands, exhibited a paradoxical behaviour in that the test statistic decreased as the effect size of the variant increased, resulting in a loss of power. The LR test always achieved the most significant P-values, followed by the exact test. We further examined the results in a real data set composed of high- and low-cholesterol subjects from the Dallas Heart Study (DHS). We also compared the single-variant LR test with two multi-variant analysis approaches--the burden test and the C-alpha test--in analysing the sequencing data by simulation. Our results call for caution in using Wald's test in genome-wide case-control association studies and suggest that the LR test is a better alternative in spite of its computational demands.

Genomic, genetic and functional dissection of bitter taste responses to artificial sweeteners

Bitter taste perception is initiated by TAS2R receptors, which respond to agonists by triggering depolarization of taste bud cells. Mutations in TAS2Rs are known to affect taste phenotypes by altering receptor function. Evidence that TAS2Rs overlap in ligand specificity suggests that they may also contribute joint effects. To explore this aspect of gustation, we examined bitter perception of saccharin and acesulfame K, widely used artificial sweeteners with aversive aftertastes. Both substances are agonists of TAS2R31 and -43, which belong to a five-member subfamily (TAS2R30-46) responsive to a diverse constellation of compounds. We analyzed sequence variation and linkage structure in the ∼140 kb genomic region encoding TAS2R30-46, taste responses to the two sweeteners in subjects, and functional characteristics of receptor alleles. Whole-gene sequences from TAS2R30-46 in 60 Caucasian subjects revealed extensive diversity including 34 missense mutations, two nonsense mutations and high-frequency copy-number variants. Thirty markers, including non-synonymous variants in all five genes, were associated (P< 0.001) with responses to saccharin and acesulfame K. However, linkage disequilibrium (LD) in the region was high (D', r(2) > 0.95). Haplotype analyses revealed that most associations were spurious, arising from LD with variants in TAS2R31. In vitro assays confirmed the functional importance of four TAS2R31 mutations, which had independent effects on receptor response. The existence of high LD spanning functionally distinct TAS2R loci predicts that bitter taste responses to many compounds will be strongly correlated even when they are mediated by different genes. Integrative approaches combining phenotypic, genetic and functional analysis will be essential in dissecting these complex relationships.

Human population genetic structure and diversity inferred from polymorphic L1(LINE-1) and Alu insertions

BACKGROUND/AIMS

The L1 retrotransposable element family is the most successful self-replicating genomic parasite of the human genome. L1 elements drive replication of Alu elements, and both have had far-reaching impacts on the human genome. We use L1 and Alu insertion polymorphisms to analyze human population structure.

METHODS

We genotyped 75 recent, polymorphic L1 insertions in 317 individuals from 21 populations in sub-Saharan Africa, East Asia, Europe and the Indian subcontinent. This is the first sample of L1 loci large enough to support detailed population genetic inference. We analyzed these data in parallel with a set of 100 polymorphic Alu insertion loci previously genotyped in the same individuals.

RESULTS AND CONCLUSION

The data sets yield congruent results that support the recent African origin model of human ancestry. A genetic clustering algorithm detects clusters of individuals corresponding to continental regions. The number of loci sampled is critical: with fewer than 50 typical loci, structure cannot be reliably discerned in these populations. The inclusion of geographically intermediate populations (from India) reduces the distinctness of clustering. Our results indicate that human genetic variation is neither perfectly correlated with geographic distance (purely clinal) nor independent of distance (purely clustered), but a combination of both: stepped clinal.

Sequencing and haplotype analysis of the Activator of CREM in the Testis (ACT) gene in populations of fertile and infertile males

cAMP-responsive element modulator (CREM) is a key transcription factor in the differentiation of round spermatids into mature spermatozoa. During spermiogenesis, CREM is regulated in part by activator of CREM in the testis (ACT), which activates CREM in a phosphorylation-independent fashion. We hypothesized that the ACT gene, which is expressed exclusively in the testis, could be involved in male factor infertility in patients with idiopathic-impaired spermatogenesis. To test this hypothesis, we sequenced the coding regions and flanking intronic regions of the ACT gene in 96 azoo- or oligospermic patients and 69 fertile controls. A total of 12 single-nucleotide polymorphisms (SNPs) was identified, and four of them leading to amino acid substitutions. An association study was performed based on calculated haplotype frequencies, and statistically significant differences were found between the patient and control populations for some haplotypes. To help establish the evolutionary relationships between the haplotypes, the coding regions of both the chimpanzee and the gorilla ACT gene were sequenced and evaluated. To test whether the different haplotypes conferred a functional change to the ACT protein, a yeast two-hybrid assay was designed to test the interaction between the two most divergent ACT haplotypes and their known binding partners, CREM and KIF17b. We identified one ACT haplotype that had a 45% reduction in its interaction with CREM. Our results suggest that different haplotypes within the ACT gene may contribute to male factor subfertility.

Identification of polymorphisms and balancing selection in the male infertility candidate gene, ornithine decarboxylase antizyme 3

Background

The antizyme family is a group of small proteins that play a role in cell growth and division by regulating the biosynthesis of polyamines (putrescine, spermidine, spermine). Antizymes regulate polyamine levels primarily through binding ornithine decarboxylase (ODC), an enzyme key to polyamine production, and targeting ODC for destruction by the 26S proteosome. Ornithine decarboxylase antizyme 3 (OAZ3) is a testis-specific antizyme paralog and the only antizyme expressed in the mid to late stages of spermatogenesis.

Methods

To see if mutations in the OAZ3 gene are responsible for some cases of male infertility, we sequenced and evaluated the genomic DNA of 192 infertile men, 48 men of known paternity, and 34 African aborigines from the Mbuti tribe in the Democratic Republic of the Congo. The coding sequence of OAZ3 was further screened for polymorphisms by SSCP analysis in the infertile group and an additional 250 general population controls. Identified polymorphisms in the OAZ3 gene were further subjected to a haplotype analysis using PHASE 2.02 and Arlequin 2.0 software programs.

Results

A total of 23 polymorphisms were identified in the promoter, exons or intronic regions of OAZ3. The majority of these fell within a region of less than two kilobases. Two of the polymorphisms, -239 A/G in the promoter and 4280 C/T, a missense polymorphism in exon 5, may show evidence of association with male infertility. Haplotype analysis identified 15 different haplotypes, which can be separated into two divergent clusters.

Conclusion

Mutations in the OAZ3 gene are not a common cause of male infertility. However, the presence of the two divergent haplotypes at high frequencies in all three of our subsamples (infertile, control, African) suggests that they have been maintained in the genome by balancing selection, which was supported by a test of Tajima's D statistic. Evidence for natural selection in this region implies that these haplotypes may be associated with a trait other than infertility. This trait may be related to another function of OAZ3 or a region in tight linkage disequilibrium to the gene.