Association between taste receptor (TAS) genes and the perception of wine characteristics

Several studies have suggested a possible relationship between polymorphic variants of the taste receptors genes and the acceptance, liking and intake of food and beverages. In the last decade investigators have attempted to link the individual ability to taste 6-n-propylthiouracil (PROP) and the sensations, such as astringency and bitterness, elicited by wine or its components, but with contradictory results. We have used the genotype instead of the phenotype (responsiveness to PROP or other tastants), to test the possible relation between genetic variability and the perception of wine characteristic in 528 subjects from Italy and the Czech Republic. We observed several interesting associations, among which the association between several TAS2R38 gene single nucleotide polymorphisms (P = 0.002) and the TAS2R16-rs6466849 polymorphism with wine sourness P = 0.0003). These associations were consistent in both populations, even though the country of origin was an important factor in the two models, thus indicating therefore that genetics alongside cultural factors also play a significant role in the individual liking of wine.

For example, lower alcohol consumption was observed in individuals carrying the PAV haplotype of the TAS2R38 gene, whereas individuals with the ancestral T allele of rs846664-TAS2R16 gene showed an increased risk of alcohol dependence [36][37][38] . In addition Lanier and colleagues reported that sweet and bitter tastes of alcoholic beverages may mediate alcohol intake 38 . Genetic variability has also been suggested to be associated with alcohol abuse 39,40 .
To further our understanding on wine acceptance or preference we have used, for the first time, a Mendelian randomization approach. We tested the possible relation between 29 taste receptor (TAS) single nucleotide polymorphisms (SNP) and the perception of wine characteristics. This approach is also supported by a recent report in which the authors identified several TAS genes, such as TAS2R4, TAS2R5, TAS2R39 and TAS2R7, to be activated by polyphenol compounds that are frequently present in red wine 41 .

Results
Data filtering and quality controls. All SNPs were in Hardy Weinberg equilibrium (HWE) with the exception of TAS2R7 -rs2588350) and TAS2R1 -rs2234233 that were therefore excluded from further analyses. TAS1R2 -rs4920566 and CA6 -rs2274333 were not genotyped successfully in the Italian population, therefore all the results regarding these two SNPs refer to the Czech population only. All samples that had a genotyping call rate lower than 75% (N = 71) were not further analysed. For the remaining samples (N = 528), the median genotyping call was 96.77%. The final number of participants who were eligible for association studies between genotype and phenotype (fully characterized by the questionnaire, stimuli perception and genotyping), included 272 Czechs (average age 43.4 ± 11.7; male 58%, female 42%) and 235 Italians (average age 37.7 ± 15.2; male 47%, female 53%). The association analyses were carried out using 27 SNPs belonging to 20 genes in a total number of 507 individuals. The study-wide Meff was 25.196, resulting in a threshold for significance, after the correction for multiple comparisons, of p = 0.05/25.196 = 0.002.

Associations between PROP sensitivity and wine descriptors.
We observed a statistically significant inverse correlation between PROP non-taster status and the level of bitterness perceived in the wine. This association was consistent in the Italian population, in the Czech population and in the pooled analysis, but it was statistically significant at the 0.05 level only in the Italian and in pooled groups, with coefficients −0.26 (95% CI −0.46; −0.06, P = 0.01) and −0.20 (95% CI −0.34; −0.06, P = 0.005) respectively. The Czech population showed the same trend but did not reach statistical significance: coefficient −0.17 (95% CI −0.37; 0.03), P = 0.09. There was no association between PROP tasting phenotype and the perceived astringency and sourness of the wine in any of the strata analysed. The results are shown in Table 1.
Associations between PROP sensitivity and TAS2R38 SNPs. We tested, as a proof of principle, the associations between the selected SNPs and PROP phenotype. As expected, we observed strong statistical associations for TAS2R38-rs10246939, TAS2R38-rs172686 and TAS2R38-rs713598, as shown in Table 2. We also observed an association between TAS2R5-rs2227264 and PROP phenotype, however this SNP and the variants in the TAS2R38 gene are in linkage disequilibrium (D' = 0.67 In 1000 G for Caucasian) and reflect the same signal Supplementary Table S1.
Associations between SNPs and wine descriptors. We observed several promising associations between the selected polymorphic variants and the wine descriptors. The most consistent was the one between the TAS2R38 SNPs and wine bitterness. All the variant alleles were inversely associated with bitterness perception. We recorded the strongest effect when comparing the variant homozygous T/T with the common C/C of the rs1726866 SNP: coefficient −1.16 (95% CI −2.13; −0.18), P = 0.02 (Czechs), coefficient −0.94 (95% CI −1.94; −0.05), P = 0.06 (Italians) and coefficient −1.10 (95% CI −1.80; −0.40), P = 0.002 (pooled data). The results are shown in Table 3 and in Supplementary Table S2. In addition, the variant allele of the TAS2R16-rs6466849 was associated with wine sourness with the strongest association observed in the pooled data: coefficient 0.94 (95% CI 0.43; 1.46), P = 0.0003. This association was consistent across all three strata but did not reach statistical significance among the Italians (Table 3 and in Supplementary Table S3). For TAS1R2-rs4920566 we observed a statistically significant association with wine astringency: coefficient −1.24 (95% CI −1.97; −0.50), P = 0.001. However this SNP was successfully genotyped only in the Czech population (Table 3 and in Supplementary Table S4). Additionally, using an unconditional logistic regression, adjusting for age, gender, nationality, we observed that individuals with at least one A allele of the TAS2R16 − rs6466849 had a decreased, although not significant, tendency of drinking wine OR 0.65 (95% CI; 0.40-1.04, P = 0.071). By stratifying for gender this association was statistically significant only in females subjects (p = 0.005) suggesting a possible interaction between gender and wine drinking.
Possible functional effects. We used several bioinformatic tools to predict possible functional relevance for TAS2R16-rs6466849. RegulomeDB showed a score of 6 suggesting minimal binding evidence. HaploReg suggested that this polymorphism alters the sequence recognized by the transcription factors SMC3_disc2, AIRE_1 and Myb_1. No significant association between rs6466849 and expression of any gene in any of the available tissues was reported in the GTEx project.

Discussion
In this study we found a statistically significant association between PROP and wine bitterness, confirming previous observations by Pickering and colleagues 8 . We did not, however, identify any association between PROP sensitivity and wine astringency and sourness as reported in that same study 8 . This divergence in the results may be explained by the different methods used to rank the wine descriptors or by the different number of volunteers involved in the two studies with ours being, by far, the larger of the two. The novelty of the present study is the use of a Mendelian randomization approach to directly assess the effect of genetic variability on the perception of wine tastes. We observed two potentially relevant signals, namely the associations between the TAS2R38 SNPs and TAS2R16-rs6466849 with wine bitterness and sourness respectively. The association between TAS2R38 polymorphisms and wine bitterness is consistent in the three strata analyzed (Czechs alone, Italians alone and the two population together) indicating that TAS2R38 "non tasters" have a tendency to judge the wine less bitter. Among the Italians, this association did not reach statistical significance suggesting a stronger effect of cultural and lifestyle habits in that country, which certainly play a relevant role in wine perception. Our findings are consistent with two studies 42,43 that reported that the bitterness of ethanol solutions differed by TAS2R38 polymorphisms. There are overwhelming evidences supporting a strong functional effect of the TAS2R38 SNPs we studied, which have been correlated with multiple human traits. Our report is, however the first one to highlight the possible   Table 3. Statistically significant associations between TAS genes polymorphisms and wine descriptors.
involvement of the polymorphic variants in wine sensing. The other novel finding of this study is represented by the association between the TAS2R16-rs6466849 SNP and wine sourness. The above gene is known to encode a receptor able to bind to a variety of natural and synthetic bitter compounds 14 . Its polymorphic variants have been associated with a variety of human traits and phenotypes, such as nicotine dependence, ageing and alcohol consumption 7,20,37 . TAS2R16-rs6466849 is situated approximately 1000 base pairs from the 3' end of the gene and bioinformatic tools used to assess its potential functional relevance did not help in elucidating its role in gene function. Considering also that wine is a complex blend of various molecules it is difficult to understand the exact mechanism underlying the association. However, given the strong statistical effect, it might be worth performing further studies to identify which molecules the receptor binds to, and how the polymorphic variant modifies the binding. Given that in a report by Hayes and colleagues 44 the genetic variability of the TAS2R16 gene was associated with alcohol intake, we tested, as an explanatory analysis, the possible relations between wine consumption, sourness and TAS2R16 genotypes. We observed that carriers of the G allele of the TAS2R16-rs6466849 showed a borderline tendency of drinking less wine. The polymorphic variant analyzed by Hayes and colleagues (rs846672) is in Linkage Disequilibrium (LD) with rs6466849 therefore the findings of the two studies likely represent the same signal and support each other. Assuming that sourness in red wine is generally perceived as a negative sensation, the different sourness experienced by carriers of the different genotypes of rs6466849 might influence alcohol liking and subsequently drinking behavior. However since we could not observe a direct association between sourness sensing and alcohol intake we lack statistical evidence to support the association between rs6466849 and alcohol liking. This hypothesis, therefore, remains highly speculative. A limitation of the study is that we performed several tests and this mitigates the robustness of the outcome. In summary using a Mendelian randomization approach we have shown a suggestive direct relation between genetic variability of taste receptors genes and wine perception. In this study we have used two populations from Europe that are similar in allelic frequencies for taste receptors, but have different dietary habits (Mediterranean Vs. central European) and we observed that one locus (TAS2R38-SNP variants) was associated with wine bitterness and another one (TAS2R16-rs6466849) with wine sourness. These associations were consistent in both populations, even if the country of origin was an important factor in the two models, thus indicating that genetics alongside cultural factors also play a significant role in individual liking of wine.

Material and Methods
Study population. A sample of 599 individuals of legal drinking age composed by students, university staff and blood donors, was recruited in two different European countries: 299 at the Institute of Experimental Medicine, Czech Academy of Sciences, Prague, Czech Republic and 300 at the Biology Department, Pisa University, Italy. Volunteers were enrolled to obtain two groups with similar mean age (40 ± 15 years), and an equal male: female ratio. Participants were asked to complete a brief questionnaire including height, weight, age, gender, smoking (number of cigarette/day, cigar, pipe), drinking (wine, beer, spirits, and glass/day or by week) and ethnicity (i.e. origin of grandparents). All participants signed a written informed consent. The study was approved by ethical review boards of the institutions responsible for subject recruitment in each of the recruitment center (University of Pisa and the Czech Academy of Sciences). All procedures performed were in accordance with the 1964 Helsinki declaration and its later amendments or comparable ethical standards.
Determination of the PROP phenotype. The individual taste determination of the PROP phenotype was determined using the staircase method. Perception response was determined using filter paper disks (2.3 cm diameter) impregnated with 100 µl solutions containing seven increasing concentrations of PROP (0.05 mM, 0.1 mM, 0.25 mM, 0.5 mM, 3 mM, 5 mM, 10 mM). The disks were dried at room temperature in the dark. Paper disks were stored in sealed plastic bags, in the dark, at room temperature up to two weeks before use; all paper disks were prepared in Pisa and half of them were immediately sent to Prague in order to avoid any bias in the preparation process. Volunteers were asked not to smoke, eat and drink coffee at least two hours before the sessions that were held in a dedicated silent room, that was adequately wide in order for volunteers not to disturb each other. Each volunteer had his/her own desk, and was assisted by qualified personnel throughout the session. After a brief description of the taste anatomy and physiology, volunteers (about 10-15 per session) were invited to place the paper disk on their tongues, to moisten it with saliva and move the paper in different parts of the mouth using the tongue mimicking chewing for at least 10 seconds. Firstly, they started with an empty paper disk in order to be trained to assess the "empty paper" taste. Then they were invited to use the same procedure to taste paper disk with increasing concentration (as described above) and to stop as soon as they had an unambiguous bitter perception of a taste. In case of substantial uncertainty (no perception or questionably detectable), volunteers were invited to taste the next concentration, and if perceived, it represented the threshold level.

Wine selection and scaling methodology. A panel of experts from the CREA-Research Centre for
Viticulture at Conegliano, chose Piave Raboso DOC red wine, 2008 vintage, produced from the native grape variety Raboso Piave (V. vinifera L.). It is cultivated in the Piave Denomination of Controlled Origin area (Veneto region, Italy), for its defined taste personality and mildness in aroma, suitable to decrease the confounding effects of olfaction. These characteristics make it easier for non-experts to evaluate the three chosen indicators: astringency, bitterness and sourness of wine. After an explanation of the three tastes to be assessed in the wine, the volunteers were provided with some examples: strong coffee and dark chocolate without sugar for bitterness; vinegar for sourness and raw artichoke and immature apple for astringency 45 . The participants were then invited to taste the red wine keeping it in the mouth for at least 7-10 seconds and moving it forward and backward. Next, they were asked to indicate the degree of bitterness, sourness and astringency using a visual analogic scale (VAS) from 0 (not detectable) to 10 (strong taste impression). Selection of tagging SNPs. Our aim was to consider genes likely to be involved in the human bitter (TAS2Rs) and sweet (TAS1Rs) taste receptor system, including functional polymorphisms belonging to genes involved in transduction pathway such as guanine nucleotide binding protein, alpha transducing 3 (GNAT3) and carbonic anhydrase VI (CA6). Therefore, according to the algorithm described by Carlson and coworkers 46 , we followed a hybrid tagging-functional method, conceived to select the maximum informative set of tag SNPs in a candidate gene/candidate region for an association study.  Table S5 shows the genes and SNPs included in the study, their position in the genome, in the gene, the amino acidic change specified and the gene function.
DNA extraction and genotyping. DNA was extracted from buccal swab following Invisorb ® Spin Swab Kit protocol and from blood samples with standard proteinase K digestion followed by phenol/chloroform extraction and ethanol precipitation. Biological material was collected during one research visit. The Invisorb ® Spin Swab Kits simple procedure comprises the following steps: lysis of cells with standard proteinase K, binding the genomic DNA to the membrane of a Spin Filter, washing the membrane with elimination of ethanol and finally elution of genomic DNA. DNA samples were transferred on 384 well PCR plates with duplicates (12%) as randomized quality controls on the plates, in order to guarantee and monitor the effective success of the PCR. All genotyping was carried out using KASPar (Kbioscience, Heddesdon, UK) assay, the PCR plates were read on a ViiA Real-Time PCR system using the 7 RUO Software(Applied Biosystems).

Statistical analysis.
Hardy Weinberg equilibrium (HWE) was tested in the two populations separately using a Pearson chi-square test. The associations between PROP phenotype, genetic polymorphisms and wine descriptors were calculated using a general linear model computing coefficients and 95% confidence intervals (95% CI), chi-square, ordinal, logistic and multinomial regression. Genetic data were coded using a codominant inheritance model with the most common genotype as the reference category. We performed the analysis considering the two centers separately and pooled them together. All the analyses were adjusted for age and gender. In the pooled analysis we also adjusted for country of origin. We additionally performed an unconditional logistic regression to assess the possible effect of TAS2R16 genotypes on alcohol intake. Considering the large number of tests performed in this study, we calculated the number of effective independent variables (Meff) for each gene 49 . We obtained a study-wide Meff = 25.196, by adding up each individual gene Meff.
Bioinformatic analysis. We used several bioinformatic tools to assess the possible functional relevance for the SNP TAS2R16-rs6466849, showing the most significant associations. RegulomeDB (http://regulome.stanford.edu/) 50 and HaploReg v2B 51 were used to identify the regulatory potential of the region nearby each SNP. The GTEx portal web site 52 was used to identify potential associations between the SNP and expression levels of nearby genes (eQTL) in all the available tissues.