Next Article in Journal
Tryptophan Intake and Metabolism in Older Adults with Mood Disorders
Next Article in Special Issue
Animal and Cellular Studies Demonstrate Some of the Beneficial Impacts of Herring Milt Hydrolysates on Obesity-Induced Glucose Intolerance and Inflammation
Previous Article in Journal
High-Sugar Diet Disrupts Hypothalamic but Not Cerebral Cortex Redox Homeostasis
Previous Article in Special Issue
Statistical and Machine-Learning Analyses in Nutritional Genomics Studies
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Genome-Wide Association Study Identifies the Association between the 12q24 Locus and Black Tea Consumption in Japanese Populations

1
Health Nutrition, Department of Applied Biological Chemistry, Graduate School of Agricultural and Life Sciences, The University of Tokyo, 1-1-1 Yayoi, Bunkyo-ku, Tokyo 113-8657, Japan
2
Department of Molecular Endocrinology, National Research Institute for Child Health and Development, 2-10-1 Okura, Setagaya-ku, Tokyo 157-8535, Japan
3
Research and Development Department, Genequest Inc., 5-29-11 Siba, Minato-ku, Tokyo 108-0014, Japan
4
Department of Genomic Data Analysis Service, Genome Analytics Japan Inc., 15-1-3205 Toyoshima-cho, Shinjuku-ku, Tokyo 162-0067, Japan
*
Authors to whom correspondence should be addressed.
Submission received: 15 September 2020 / Revised: 9 October 2020 / Accepted: 10 October 2020 / Published: 18 October 2020
(This article belongs to the Special Issue Genomics and Personalized Nutrition)

Abstract

:
Several genome-wide association studies (GWASs) have reported the association between genetic variants and the habitual consumption of foods and drinks; however, no association data are available regarding the consumption of black tea. The present study aimed to identify genetic variants associated with black tea consumption in 12,258 Japanese participants. Data on black tea consumption were collected by a self-administered questionnaire, and genotype data were obtained from a single nucleotide polymorphism array. In the discovery GWAS, two loci met suggestive significance (p < 1.0 × 10−6). Three genetic variants (rs2074356, rs144504271, and rs12231737) at 12q24 locus were also significantly associated with black tea consumption in the replication stage (p < 0.05) and during the meta-analysis (p < 5.0 × 10−8). The association of rs2074356 with black tea consumption was slightly attenuated by the additional adjustment for alcohol drinking frequency. In conclusion, genetic variants at the 12q24 locus were associated with black tea consumption in Japanese populations, and the association is at least partly mediated by alcohol drinking frequency.

Graphical Abstract

1. Introduction

Black tea, which is derived from the Camellia sinensis plant, is one of the most widely consumed beverages in the world owing to its several beneficial compounds, including polyphenols, amino acids, vitamins, and minerals. Nutritional studies have demonstrated that theaflavin, a polyphenol, has positive effects against oxidative stress, endoplasmic reticulum (ER) stress, hypertension, aging, and diabetes [1,2,3,4], and although it remains controversial, a meta-analysis of epidemiological studies indicated that consuming three cups per day of black tea could prevent the onset of ischemic stroke [5].
Recently, genome-wide association studies (GWASs) have been conducted to investigate the association between genetic variants and the habitual intake of foods and drinks [6,7]. Previous studies identified the association of genetic variants with the intake frequency of fish [8], alcohol consumption [9,10,11], and coffee consumption [12,13,14]. Furthermore, dietary trait-associated genetic variants were linked to health outcomes in Japanese populations [7]. Thus, it is crucial to accumulate data from GWASs to identify the relationship between genetic variants and dietary habits. However, no data are available regarding black tea consumption-related genetic variants in individuals of other ethnicities, including Asian, European, and American. Therefore, this study performed a GWAS to search for genetic variants associated with the habitual consumption of black tea in 12,258 Japanese participants.

2. Materials and Methods

2.1. Study Design

The participants were enrolled from customers of a direct-to-consumer genetic testing service in Japan, “HealthData Lab,” provided by Yahoo! Japan Corporation (Tokyo, Japan) and Genequest Inc. (Tokyo, Japan). The inclusion criterion was age ≥ 18 years. Eligible participants were asked to complete internet-based questionnaires covering socio-demographic factors, medical history, and lifestyle habits at the time of enrollment. Among the 12,621 participants, one participant who opted out was excluded from the present study.
This study was approved by the Ethics Committee of Genequest Inc. (2015-0907-1) and was conducted according to the principles expressed in the Declaration of Helsinki. Written informed consent was obtained from all the participants for the general use of the questionnaires and their genotype data. An additional agreement was obtained with the opportunity to opt-out after informing the participants of this study’s aim.

2.2. Black Tea Consumption

The questionnaires included a question regarding black tea consumption, i.e., “How many cups of black tea do you drink?” with the following seven answers: (i) hardly any, (ii) less than or equal to two cups per week, (iii) three to four cups per week, (iv) five to six cups per week, (v) one to two cups per day, (vi) three to four cups per day, and (vii) more than or equal to five cups per day. These answers were converted into continuous variables (cups per day): (i) 0 cups, (ii) 0.14 (=1/7) cups, (iii) 0.5 (=3.5/7) cups, (iv) 0.79 (=5.5/7) cups, (v) 1.5 cups, (vi) 3.5 cups, and (vii) 5 cups. This continuous value was then used to define as the amount of habitual black tea consumption.

2.3. Adjustment Variables

Data on alcohol consumption, alcohol drinking frequency, habitual coffee consumption, sweet preference, and body-mass index (BMI) were used as adjustment variables. These data were obtained from the questionnaire, as described previously [8,12,15]. Briefly, alcohol drinking frequency was defined according to the question, “How frequently do you drink alcohol?” with the following seven answers: (i) barely ever, (ii) less than once a month, (iii) one to three times per month, (iv) one to two times per week, (v) three to four times per week, (vi) five to six times per week, and (vii) every day. Alcohol consumption was calculated by adding the amount of alcohol in grams per day obtained from beer, red wine, white wine, highballs/cocktails, rice wine, and distilled spirits.
Coffee consumption was defined based on the response to the question, “How many cups of coffee (instant or regular) do you drink?” and “How many cups of coffee (can, PET bottles, or paper pack) do you drink?” with a choice of seven answers: (i) hardly any, (ii) less than or equal to two cups per week, (iii) three to four cups per week, (iv) five to six cups per week, (v) one to two cups per day, (vi) three to four cups per day, and (vii) more than or equal to five cups per day. The sum of these values was used to define as the amount of habitual coffee consumption [12].
Sweet taste preference was determined using the question, “Please tell us about your taste for sweet food,” with the following five-point scale: (i) 1, dislike a lot, (ii) 2, dislike, (iii) 3, neither like nor dislike, (iv) 4, like, and (v) 5, like a lot.

2.4. DNA Sampling, Genotype, Quality Control, and Genotype Imputation

The collection, stabilization, and transportation of saliva samples were performed using an Oragene DNA Collection Kit (DNA Genotek Inc., Ottawa, Ontario, Canada). Genotype analysis was performed using The Illumina HumanCore-12 + Custom BeadChip (Illumina, San Diego, CA, USA), which contains 302,073 markers, and the HumanCore-24 + Custom BeadChip, which contains 309,725 markers. We used a total of 296,675 markers in the present study.
After excluding 17 subjects who lived outside of Japan, 12,603 participants remained. They were divided into two groups: those from the Hondo region (except for Okinawa region) and those from the Okinawa region, and they were referred to as the discovery and replication cohorts, respectively. Then, we applied the quality control and association analysis procedure separately for each cohort.
In the quality control analysis, we filtered out single nucleotide polymorphism (SNP) markers with low call rates (<0.95), low Hardy–Weinberg equilibrium exact test p-values (<1 × 10−6), or low minor allele frequencies (<0.01). We also excluded subjects who had inconsistent sex information between the genotype and questionnaire, who had an estimated non-Japanese ancestry, or who had a low call rate (<0.95) [16,17]. Furthermore, we excluded the close relationship pairs determined using the identity-by-descent method (PI_HAT > 0.1875) according to previous studies [8,13,18,19]. Quality control analyses were carried out using PLINK [20,21] (version 1.90b3.42) and Eigensoft [16] (version 6.1.3).
The imputation of the genotype was performed using the 1000 Genome Phase 3 (version 5) reference panel [22]. Pre-phasing was carried out using EAGLE2 (version 2.4) [23], and genotype imputation was performed using Minimac3 (version 2.0.1) software [24]. We excluded variants with a low imputation quality (R2 < 0.8) and a low minor allele frequency (<1%) from further analysis. Finally, we used the dosage data for the 5,256,047 variants for the GWAS in the discovery phase.

2.5. Genome-Wide Association and Meta-Analysis

The association between genotype dosage and the amount of black tea consumption was examined using a linear regression model under the assumption of additive genetic effects. For the discovery GWAS, the association was tested with an adjustment for age, sex, and five principal components calculated using PLINK software. Variants that met suggestive significance (p < 1.0 × 10−6) in the discovery phase were further investigated in the replication stage to determine their associations with habitual black tea consumption. Thereafter, we combined the statistical data from both stages using a fixed-effects model and the inverse-variance weighting method with METAL software (version 2011-03-25) [25]. Variants achieving genome-wide significance (p < 5.0 × 10−8) in the meta-analysis were considered to be associated with habitual black tea consumption.

2.6. Confounding Factor Adjustment and Subgroup Analysis

Multivariate linear regression analysis was carried out to test the associations between SNPs and black tea consumption with an additional adjustment for alcohol frequency, alcohol consumption, coffee consumption, sweet preference, and BMI. The subgroup analysis was carried out according to sex and age, and P < 0.05 was considered statistically significant in the subgroup analysis. These analyses were conducted based on the discovery populations.

3. Results

3.1. Research Flow and Characteristics of the Study Participants

Since the population of the Okinawa region forms a genetic cluster different from that of the other regions (Hondo region) in Japan [17], we investigated the effects of common variants on black tea consumption in 12,140 individuals for the discovery GWAS (Hondo region) and 118 individuals for the replication stage (Okinawa region), as illustrated in Figure 1. Thereafter, a meta-analysis from both the discovery and replication populations was carried out to identify black tea consumption-associated loci.
The characteristics of the Japanese participants from the discovery and replication stages are shown in Table 1. The mean age of the discovery and replication cohorts was 50.3 ± 13.2 years and 49.0 ± 12.3 years, respectively, and 46.8% and 45.8% of cohorts were female, respectively. The mean amount of black tea consumption per day was 0.20 ± 0.51 cups for the discovery population and 0.15 ± 0.41 cups for the replication population.

3.2. Discovery GWAS

The quantile-quantile plot indicated that the inflation factor (λ) was 1.005 and the 95% confidence interval was 1.003–1.007 in the discovery GWAS (Figure 2), indicating that the population structure was well-adjusted.
As illustrated in Figure 3, the discovery GWAS found that the 12q15 and 12q24 loci met suggestive significance (p < 1 × 10−6). Eleven SNPs at 12q24 locus and one SNP at 12q15 locus were detected as candidates for black tea consumption-associated variants (Table 2). The variants were positively associated with black tea consumption in the discovery stage.

3.3. Replication Stage and Meta-Analysis

The candidate variants found in the discovery GWAS were further examined in the replication stage. The eleven SNPs at 12q24 locus were associated with habitual black tea consumption (p < 0.05), whereas no significant association was observed in the case of rs1981764 (p = 0.81). A meta-analysis of both populations showed that three variants (rs2074356, rs144504271, and rs12231737) met genome-wide significance (p < 5 × 10−8), but other variants did not. With regard to the top SNP, rs2074356, the frequency of the T allele was 24.3%, and the effect size was estimated to be 0.042 (SE = 0.008) cups/day per allele in our meta-analysis. Significant heterogeneity in the genetic effects of rs2074356 was observed in the present study (I2 = 76.1 and p-value for heterogeneity = 0.04).

3.4. Adjustment for Potential Confounding Factors

We focused on the top SNP, rs2074356, for further analysis. The 12q24 locus is a strong long-range linkage disequilibrium (LD) region [26]. Previous GWASs in Asian populations showed the association of the 12q24 locus with dietary behaviors, preferences, and health outcomes [8,12,27,28]. Thus, the association of rs2074356 with black tea consumption might be confounded by the above-mentioned variables. Accordingly, we performed a multivariate linear regression analysis with further adjustment for dietary behaviors, preferences, and BMI. As shown in Table 3, the association between the 12q24 locus and black tea consumption was slightly attenuated by adding drinking frequency as an adjustment variable (p-value was changed from 1.2 × 10−7 to 3.8 × 10−5), although the association was still significant after the additional adjustment. In contrast, the addition of other variables hardly attenuated the association between the 12q24 locus and black tea consumption.

3.5. Subgroup Analysis According to Sex and Age

Subgroup analysis was conducted to examine whether the genetic effect of the variants at the 12q24 locus differed according to sex or age. In the subgroup analysis based on age, the median age of 51 years was used as the cutoff to categorize the individuals into the younger and older subgroups.
The results showed that the effect size of the rs2074356 variant was 0.028 (SE = 0.009) and 0.054 (SE = 0.013) cups per day in males and females, respectively. The genetic effect of the variant was slightly higher in females than in males, but the difference did not reach significance (p for interaction = 0.13; Figure 4a). No significant change was observed in the effect size between the younger subgroup (0.034, SE = 0.011) and the older subgroup (0.046, SE = 0.010) (p for interaction = 0.26; Figure 4b).

4. Discussion

The present study indicated an association of the 12q24 locus with black tea consumption in Japanese populations. Given that the 12q24 locus is specific to East Asian populations, this association cannot be found in studies of non-East Asian populations.
Regarding caffeine beverages, previous GWASs mainly investigated participants’ coffee consumption and identified several genetic variants in various populations [12,13,14]. Several studies also identified the variants related to caffeine intake, which was calculated by summing the consumption of caffeinated coffee, tea, soft drinks, and chocolate in individuals of European descent in populations from the United States and Costa Rica [29,30], whereas no GWAS data were available on the variants that are specific to black tea consumption. To the best of our knowledge, the present study is the first GWAS of habitual black tea consumption.
Furthermore, the 12q24 locus affects several dietary habits in East Asian populations [7,8,12,13,31]. Our analysis indicated that although the association between the 12q24 locus and black tea consumption was slightly attenuated by the further adjustment for drinking frequency, the association was still significant after the adjustment. Previous studies have indicated that the association of the 12q24 variant with fish intake frequency [8] (p-value was changed from 4.3 × 10−11 to 2.5 × 10−4) and sweet preference [15] (p-value from 1.2 × 10−69 to 6.5 × 10−10) was largely attenuated by adding drinking frequency as an adjustment variable. Moreover, the association of the 12q24 locus with black tea consumption was not evidently attenuated by additional adjustment for alcohol consumption. Thus, the association of the 12q24 locus with black tea consumption may be partly confounded by drinking frequency; regardless, the 12q24 locus certainly had a direct effect on black tea consumption.
In our meta-analysis, the rs2074356 SNP at the 12q24 locus had a beta value of 0.042, which indicated that Japanese participants’ black tea intake increased by 0.042 cups per day, 1.26 cups per month, and more than 15 cups per year for each copy of the minor allele. Considering that black tea contains several beneficial nutrients and other compounds, it is conceivable that the variants have health-promoting effects in East Asian populations. Particularly, caffeine is widely known to have health benefits related to type 2 diabetes mellitus and glucose tolerance [31]. Previous studies indicated that rs2074356 SNP was associated with coffee consumption in Japanese [13] and Korean populations [32], suggesting that the SNP could be associated with caffeine intake in East Asian populations. Theaflavin, a polyphenol in black tea, was also reported to suppress oxidative stress and hyperglycemia in a streptozotocin-induced diabetic model [4] and ER stress and hypertension in an angiotensin II-induced hypertension model [2]. Furthermore, a recent GWAS in Japanese populations indicated a strong association of the 12q24 locus with liver disorder-related parameters (alanine transaminase and aspartate transaminase) and blood pressure [7]. With regard to black tea-associated SNPs, HECTD4 is a ubiquitin ligase, and TRAFD1 relates to the toll-like receptor 4 signaling pathway [33], whereas there are limited data on the association between health outcomes and the three SNPs. Therefore, it is still important to study the complicated relationship between genetic background at 12q24 locus, dietary and drinking habits, and health outcomes.
Recently, GWASs have demonstrated genetic variants associated with dietary habits in several populations, and the results of these studies could be utilized for personalized nutritional recommendations [34]. In a Japanese population, a previous GWAS indicated that the 12q24 locus was associated with several dietary habits and health outcomes [7]. Thus, the genotype of the 12q24 locus may be a useful factor for personalized nutritional recommendations. In future studies, it is pertinent to accumulate GWAS data to develop novel personalized nutritional recommendations.
This study had the following two limitations. First, the 12q24 locus is specific to East Asian populations; however, our study investigated only Japanese populations. Second, we could not deny the possibility that the association between the 12q24 locus and black tea consumption may be confounded by unmeasured factors, including health outcomes and dietary factors. Therefore, further research with several populations and unmeasured factors are still required to determine the impacts of the genetic variants on black tea consumption and its related health outcomes.

5. Conclusions

We found that genetic variants at the 12q24 locus were associated with black tea consumption in Japanese populations. This finding not only provided fundamental information with regard to the association between genetic variants and dietary habits but also contributed to the development of novel strategies for personalized nutritional recommendations.

Author Contributions

K.F. wrote the manuscript and M.I., S.N., T.H., and H.J. revised the manuscript. K.S. oversaw the data management and prepared the genotype dataset for the GWAS. S.N. and K.K. performed the statistical analysis, GWAS, and meta-analysis. T.H. provided valuable suggestions for statistical analysis and interpretation of the obtained results. H.K. and H.J. supervised the work. H.K., H.J., M.I., and S.T. designed the project. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Grants-in-Aid (16K14920) from the Japan Society for the Promotion of Science and by the Honjo International Scholarship Foundation.

Conflicts of Interest

S.N. and K.K. are employees of Genequest Inc., and S.T. and K.S. are board members of Genequest Inc. T.H. is a board member of Genome Analytics Japan Inc. and an advisor of Genequest Inc. The other authors declare that they have no conflict of interest.

References

  1. Leung, L.K.; Su, Y.; Chen, R.; Zhang, Z.; Huang, Y.; Chen, Z.-Y. Theaflavins in Black Tea and Catechins in Green Tea Are Equally Effective Antioxidants. J. Nutr. 2001, 131, 2248–2251. [Google Scholar] [CrossRef]
  2. Cheang, W.S.; Ngai, C.Y.; Tam, Y.Y.; Tian, X.Y.; Wong, W.T.; Zhang, Y.; Lau, C.W.; Chen, Z.Y.; Bian, Z.X.; Huang, Y.; et al. Black tea protects against hypertension-associated endothelial dysfunction through alleviation of endoplasmic reticulum stress. Sci. Rep. 2015, 5, 10340. [Google Scholar] [CrossRef]
  3. Cameron, A.R.; Anton, S.; Melville, L.; Houston, N.P.; Dayal, S.; McDougall, G.J.; Stewart, D.; Rena, G. Black tea polyphenols mimic insulin/insulin-like growth factor-1 signalling to the longevity factor FOXO1a. Aging Cell 2008, 7, 69–77. [Google Scholar] [CrossRef] [PubMed]
  4. Gothandam, K.; Ganesan, V.S.; Ayyasamy, T.; Ramalingam, S. Antioxidant potential of theaflavin ameliorates the activities of key enzymes of glucose metabolism in high fat diet and streptozotocin-induced diabetic rats. Redox Rep. Commun. Free Radic. Res. 2019, 24, 41–50. [Google Scholar] [CrossRef] [Green Version]
  5. Arab, L.; Liu, W.; Elashoff, D. Green and black tea consumption and risk of stroke: A meta-analysis. Stroke 2009, 40, 1786–1792. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Robino, A.; Concas, M.P.; Catamo, E.; Gasparini, P. A Brief Review of Genetic Approaches to the Study of Food Preferences: Current Knowledge and Future Directions. Nutrients 2019, 11, 1735. [Google Scholar] [CrossRef] [Green Version]
  7. Matoba, N.; Akiyama, M.; Ishigaki, K.; Kanai, M.; Takahashi, A.; Momozawa, Y.; Ikegawa, S.; Ikeda, M.; Iwata, N.; Hirata, M.; et al. GWAS of 165,084 Japanese individuals identified nine loci associated with dietary habits. Nat. Hum. Behav. 2020, 4, 308–316. [Google Scholar] [CrossRef] [PubMed]
  8. Igarashi, M.; Nogawa, S.; Kawafune, K.; Hachiya, T.; Takahashi, S.; Jia, H.; Saito, K.; Kato, H. Identification of the 12q24 locus associated with fish intake frequency by genome-wide meta-analysis in Japanese populations. Genes Nutr. 2019, 14, 21. [Google Scholar] [CrossRef] [PubMed]
  9. Jorgenson, E.; Thai, K.K.; Hoffmann, T.J.; Sakoda, L.C.; Kvale, M.N.; Banda, Y.; Schaefer, C.; Risch, N.; Mertens, J.; Weisner, C.; et al. Genetic contributors to variation in alcohol consumption vary by race/ethnicity in a large multi-ethnic genome-wide association study. Mol. Psychiatry 2017, 22, 1359–1367. [Google Scholar] [CrossRef] [Green Version]
  10. Quillen, E.E.; Chen, X.-D.; Almasy, L.; Yang, F.; He, H.; Li, X.; Wang, X.-Y.; Liu, T.-Q.; Hao, W.; Deng, H.-W.; et al. ALDH2 is associated to alcohol dependence and is the major genetic determinant of “daily maximum drinks” in a GWAS study of an isolated rural Chinese sample. Am. J. Med. Genet. B Neuropsychiatr. Genet. 2014, 165B, 103–110. [Google Scholar] [CrossRef] [Green Version]
  11. Yang, X.; Lu, X.; Wang, L.; Chen, S.; Li, J.; Cao, J.; Chen, J.; Hao, Y.; Li, Y.; Zhao, L.; et al. Common variants at 12q24 are associated with drinking behavior in Han Chinese. Am. J. Clin. Nutr. 2013, 97, 545–551. [Google Scholar] [CrossRef] [PubMed]
  12. Jia, H.; Nogawa, S.; Kawafune, K.; Hachiya, T.; Takahashi, S.; Igarashi, M.; Saito, K.; Kato, H. GWAS of habitual coffee consumption reveals a sex difference in the genetic effect of the 12q24 locus in the Japanese population. BMC Genet. 2019, 20, 61. [Google Scholar] [CrossRef] [PubMed]
  13. Nakagawa-Senda, H.; Hachiya, T.; Shimizu, A.; Hosono, S.; Oze, I.; Watanabe, M.; Matsuo, K.; Ito, H.; Hara, M.; Nishida, Y.; et al. A genome-wide association study in the Japanese population identifies the 12q24 locus for habitual coffee consumption: The J-MICC Study. Sci. Rep. 2018, 8, 1493. [Google Scholar] [CrossRef]
  14. Pirastu, N.; Kooyman, M.; Robino, A.; van der Spek, A.; Navarini, L.; Amin, N.; Karssen, L.C.; Van Duijn, C.M.; Gasparini, P. Non-additive genome-wide association scan reveals a new gene associated with habitual coffee consumption. Sci. Rep. 2016, 6, 31590. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Kawafune, K.; Hachiya, T.; Nogawa, S.; Takahashi, S.; Jia, H.; Saito, K.; Kato, H. Strong association between the 12q24 locus and sweet taste preference in the Japanese population revealed by genome-wide meta-analysis. J. Hum. Genet. 2020. [Google Scholar] [CrossRef] [PubMed]
  16. Price, A.L.; Patterson, N.J.; Plenge, R.M.; Weinblatt, M.E.; Shadick, N.A.; Reich, D. Principal components analysis corrects for stratification in genome-wide association studies. Nat. Genet. 2006, 38, 904–909. [Google Scholar] [CrossRef]
  17. Yamaguchi-Kabata, Y.; Nakazono, K.; Takahashi, A.; Saito, S.; Hosono, N.; Kubo, M.; Nakamura, Y.; Kamatani, N. Japanese population structure, based on SNP genotypes from 7003 individuals compared to other ethnic groups: Effects on population-based association studies. Am. J. Hum. Genet. 2008, 83, 445–456. [Google Scholar] [CrossRef] [Green Version]
  18. Hachiya, T.; Komaki, S.; Hasegawa, Y.; Ohmomo, H.; Tanno, K.; Hozawa, A.; Tamiya, G.; Yamamoto, M.; Ogasawara, K.; Nakamura, M.; et al. Genome-wide meta-analysis in Japanese populations identifies novel variants at the TMC6–TMC8 and SIX3–SIX2 loci associated with HbA1c. Sci. Rep. 2017, 7, 16147. [Google Scholar] [CrossRef]
  19. Shimanoe, C.; Hachiya, T.; Hara, M.; Nishida, Y.; Tanaka, K.; Sutoh, Y.; Shimizu, A.; Hishida, A.; Kawai, S.; Okada, R.; et al. A genome-wide association study of coping behaviors suggests FBXO45 is associated with emotional expression. Genes Brain Behav. 2019, 18, e12481. [Google Scholar] [CrossRef]
  20. Chang, C.C.; Chow, C.C.; Tellier, L.C.; Vattikuti, S.; Purcell, S.M.; Lee, J.J. Second-generation PLINK: Rising to the challenge of larger and richer datasets. GigaScience 2015, 4, 7. [Google Scholar] [CrossRef]
  21. Purcell, S.; Neale, B.; Todd-Brown, K.; Thomas, L.; Ferreira, M.A.R.; Bender, D.; Maller, J.; Sklar, P.; de Bakker, P.I.W.; Daly, M.J.; et al. PLINK: A Tool Set for Whole-Genome Association and Population-Based Linkage Analyses. Am. J. Hum. Genet. 2007, 81, 559–575. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. The Genomes Project Consortium; Auton, A.; Abecasis, G.R.; Altshuler, D.M.; Durbin, R.M.; Abecasis, G.R.; Bentley, D.R.; Chakravarti, A.; Clark, A.G.; Donnelly, P.; et al. A global reference for human genetic variation. Nature 2015, 526, 68. [Google Scholar] [CrossRef] [Green Version]
  23. Loh, P.-R.; Danecek, P.; Palamara, P.F.; Fuchsberger, C.; Reshef, Y.A.; Finucane, H.K.; Schoenherr, S.; Forer, L.; McCarthy, S.; Abecasis, G.R.; et al. Reference-based phasing using the Haplotype Reference Consortium panel. Nat. Genet. 2016, 48, 1443–1448. [Google Scholar] [CrossRef] [Green Version]
  24. Das, S.; Forer, L.; Schönherr, S.; Sidore, C.; Locke, A.E.; Kwong, A.; Vrieze, S.I.; Chew, E.Y.; Levy, S.; McGue, M.; et al. Next-generation genotype imputation service and methods. Nat. Genet. 2016, 48, 1284–1287. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Willer, C.J.; Li, Y.; Abecasis, G.R. METAL: Fast and efficient meta-analysis of genomewide association scans. Bioinformatics (Oxford) 2010, 26, 2190–2191. [Google Scholar] [CrossRef]
  26. Kim, Y.J.; Go, M.J.; Hu, C.; Hong, C.B.; Kim, Y.K.; Lee, J.Y.; Hwang, J.-Y.; Oh, J.H.; Kim, D.-J.; Kim, N.H.; et al. Large-scale genome-wide association studies in East Asians identify new genetic loci influencing metabolic traits. Nat. Genet. 2011, 43, 990–995. [Google Scholar] [CrossRef]
  27. Zhang, H.; Mo, X.B.; Xu, T.; Bu, X.Q.; Lei, S.F.; Zhang, Y.H. Novel Genes Affecting Blood Pressure Detected Via Gene-Based Association Analysis. G3 (Bethesda) 2015, 5, 1035–1042. [Google Scholar] [CrossRef] [Green Version]
  28. Wen, W.; Zheng, W.; Okada, Y.; Takeuchi, F.; Tabara, Y.; Hwang, J.Y.; Dorajoo, R.; Li, H.; Tsai, F.J.; Yang, X.; et al. Meta-analysis of genome-wide association studies in East Asian-ancestry populations identifies four new loci for body mass index. Hum. Mol. Genet. 2014, 23, 5492–5504. [Google Scholar] [CrossRef] [Green Version]
  29. Cornelis, M.C.; Monda, K.L.; Yu, K.; Paynter, N.; Azzato, E.M.; Bennett, S.N.; Berndt, S.I.; Boerwinkle, E.; Chanock, S.; Chatterjee, N.; et al. Genome-wide meta-analysis identifies regions on 7p21 (AHR) and 15q24 (CYP1A2) as determinants of habitual caffeine consumption. PLoS Genet. 2011, 7, e1002033. [Google Scholar] [CrossRef]
  30. Josse, A.R.; Da Costa, L.A.; Campos, H.; El-Sohemy, A. Associations between polymorphisms in the AHR and CYP1A1-CYP1A2 gene regions and habitual caffeine consumption. Am. J. Clin. Nutr. 2012, 96, 665–671. [Google Scholar] [CrossRef] [Green Version]
  31. Higdon, J.V.; Frei, B. Coffee and Health: A Review of Recent Human Research. Crit. Rev. Food Sci. Nutr. 2006, 46, 101–123. [Google Scholar] [CrossRef]
  32. Jin, T.; Youn, J.; Kim, A.N.; Kang, M.; Kim, K.; Sung, J.; Lee, J.E. Interactions of Habitual Coffee Consumption by Genetic Polymorphisms with the Risk of Prediabetes and Type 2 Diabetes Combined. Nutrients 2020, 12, 2228. [Google Scholar] [CrossRef]
  33. Fang, C.M.; Roy, S.; Nielsen, E.; Paul, M.; Maul, R.; Paun, A.; Koentgen, F.; Raval, F.M.; Szomolanyi-Tsuda, E.; Pitha, P.M. Unique contribution of IRF-5-Ikaros axis to the B-cell IgG2a response. Genes Immun. 2012, 13, 421–430. [Google Scholar] [CrossRef] [Green Version]
  34. Drabsch, T.; Holzapfel, C. A Scientific Perspective of Personalised Gene-Based Dietary Recommendations for Weight Management. Nutrients 2019, 11, 617. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Research flow of the present study. GWAS: Genome-wide association study.
Figure 1. Research flow of the present study. GWAS: Genome-wide association study.
Nutrients 12 03182 g001
Figure 2. Quantile-quantile plot for the discovery GWAS. The x-axis and y-axis show the theoretical –log10 p-value and the observed –log10 p-value, respectively. The red line represents y = x. CI: confidence interval, GWAS: genome-wide association study.
Figure 2. Quantile-quantile plot for the discovery GWAS. The x-axis and y-axis show the theoretical –log10 p-value and the observed –log10 p-value, respectively. The red line represents y = x. CI: confidence interval, GWAS: genome-wide association study.
Nutrients 12 03182 g002
Figure 3. Manhattan plot for the discovery GWAS of habitual black tea consumption in the Japanese population. Single nucleotide polymorphism (SNPs) are ordered by chromosome and position in the x-axis. The y-axis shows the negative logarithm of the association of each SNP with black tea consumption. The red line represents p < 1.0 × 10−6. GWAS: genome-wide association study.
Figure 3. Manhattan plot for the discovery GWAS of habitual black tea consumption in the Japanese population. Single nucleotide polymorphism (SNPs) are ordered by chromosome and position in the x-axis. The y-axis shows the negative logarithm of the association of each SNP with black tea consumption. The red line represents p < 1.0 × 10−6. GWAS: genome-wide association study.
Nutrients 12 03182 g003
Figure 4. Effects of rs2074356 on habitual black tea consumption stratified by sex (a) and age (b). The black tea consumption is marked for each allele, CC, TC, and TT, for rs2074356. Error bars and β indicate the standard deviation and the estimated regression coefficient, respectively.
Figure 4. Effects of rs2074356 on habitual black tea consumption stratified by sex (a) and age (b). The black tea consumption is marked for each allele, CC, TC, and TT, for rs2074356. Error bars and β indicate the standard deviation and the estimated regression coefficient, respectively.
Nutrients 12 03182 g004aNutrients 12 03182 g004b
Table 1. Participant characteristics.
Table 1. Participant characteristics.
ParametersDiscoveryReplication
N12,140118
Female (%)46.845.8
Age, years (mean ± SD)50.3 ± 13.249.0 ± 12.3
Black tea consumption, cups/day (mean ± SD)0.20 ± 0.510.15 ± 0.41
Drinking frequency, times/week (mean ± SD)2.21 ± 2.672.63 ± 2.76
Alcohol consumption, g/day (means ± SD)7.10 ± 11.949.87 ± 15.02
Coffee consumption, cups/day (mean ± SD)1.70 ± 1.501.64 ± 1.28
Sweet taste preference, (mean ± SD)3.74 ± 0.903.72 ± 0.89
BMI, kg/m2 (mean ± SD)23.1 ± 3.723.9 ± 4.0
BMI: body mass index, SD: standard deviation.
Table 2. Variants associated with habitual black tea consumption.
Table 2. Variants associated with habitual black tea consumption.
SNPChrPositionGeneEANEAPopulationEAFBetaSE (Beta)pAssociation
rs207435612112645401HECTD4TCDiscovery0.2440.0400.0081.2 × 10−7
Replication0.1740.1840.0700.01
Meta-analysis0.2430.0420.0082.4 × 10−8
rs14450427112112627350HECTD4AGDiscovery0.2640.0410.0081.2 × 10−7
Replication0.1900.1830.0720.01
Meta-analysis0.2630.0420.0083.1 × 10−8
rs1223173712112574616TRAFD1TCDiscovery0.2660.0410.0081.4 × 10−7
Replication0.1920.1870.0720.01
Meta-analysis0.2650.0430.0084.3 × 10−8
rs11687308712112511913NAA25CGDiscovery0.2630.0410.0082.2 × 10−7
Replication0.1890.1860.0740.01
Meta-analysis0.2620.0430.0085.8 × 10−8
rs1106613212112468206NAA25TCDiscovery0.2610.0400.0083.2 × 10−7
Replication0.1880.1880.0740.01
Meta-analysis0.2600.0420.0089.6 × 10−8
rs7806906612112337924MAPKAPK5
TMEM116
AGDiscovery0.2670.0390.0083.5 × 10−7
Replication0.1890.1770.0720.01
Meta-analysis0.2660.0400.0081.1 × 10−7
rs464677612112230019ALDH2CGDiscovery0.2640.0370.0074.4 × 10−7
Replication0.1870.1710.0700.02
Meta-analysis0.2630.0390.0071.4 × 10−7
rs67112112241766ALDH2AGDiscovery0.2640.0370.0074.5 × 10−7
Replication0.1860.1710.0700.02
Meta-analysis0.2630.0390.0071.5 × 10−7
rs1106600112112119171BRAPCTDiscovery0.2620.0380.0086.2 × 10−7
Replication0.1850.1750.0720.02
Meta-analysis0.2610.0390.0081.6 × 10−7
rs1106601512112168009ACAD10AGDiscovery0.2630.0370.0075.7 × 10−7
Replication0.1870.1710.0700.02
Meta-analysis0.2620.0380.0071.8 × 10−7
rs378288612112110489BRAPCTDiscovery0.2640.0370.0087.9 × 10−7
Replication0.1860.1750.0720.02
Meta-analysis0.2630.0390.0082.6 × 10−7
rs19817641268217659DYRK2– IFNGGADiscovery0.1230.0490.0109.9 × 10−7
Replication0.119−0.0210.0870.81
Meta-analysis0.1270.0480.0101.4 × 10−6
Chromosomal positions are based on the human genome assembly version GRCh37/hg19. SNP: single nucleotide polymorphism, Chr: chromosome, EA: effect allele, NEA: non-effect allele, EAF: effect allele frequency, SE: standard error.
Table 3. Additional adjustment for confounding factors.
Table 3. Additional adjustment for confounding factors.
Adjustment VariablesBetaSE (Beta)p Association
Age, Sex, Population structure (5 PCs)0.0400.0081.2 × 10−7
Age, Sex, Population structure (5 PCs), Drinking frequency0.0330.0083.8 × 10−5
Age, Sex, Population structure (5 PCs), Alcohol consumption0.0370.0083.1 × 10−6
Age, Sex, Population structure (5 PCs), Coffee consumption0.0450.0089.9 × 10−9
Age, Sex, Population structure (5 PCs), Sweet preference0.0390.0083.0 × 10−7
Age, Sex, Population structure (5 PCs), BMI0.0400.0089.9 × 10−8
SE: standard error, BMI: body mass index, PC: principal component.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Furukawa, K.; Igarashi, M.; Jia, H.; Nogawa, S.; Kawafune, K.; Hachiya, T.; Takahashi, S.; Saito, K.; Kato, H. A Genome-Wide Association Study Identifies the Association between the 12q24 Locus and Black Tea Consumption in Japanese Populations. Nutrients 2020, 12, 3182. https://0-doi-org.brum.beds.ac.uk/10.3390/nu12103182

AMA Style

Furukawa K, Igarashi M, Jia H, Nogawa S, Kawafune K, Hachiya T, Takahashi S, Saito K, Kato H. A Genome-Wide Association Study Identifies the Association between the 12q24 Locus and Black Tea Consumption in Japanese Populations. Nutrients. 2020; 12(10):3182. https://0-doi-org.brum.beds.ac.uk/10.3390/nu12103182

Chicago/Turabian Style

Furukawa, Kyohei, Maki Igarashi, Huijuan Jia, Shun Nogawa, Kaoru Kawafune, Tsuyoshi Hachiya, Shoko Takahashi, Kenji Saito, and Hisanori Kato. 2020. "A Genome-Wide Association Study Identifies the Association between the 12q24 Locus and Black Tea Consumption in Japanese Populations" Nutrients 12, no. 10: 3182. https://0-doi-org.brum.beds.ac.uk/10.3390/nu12103182

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop