Next Article in Journal
Understanding the Mechanisms of Diet and Outcomes in Colon, Prostate, and Breast Cancer; Malignant Gliomas; and Cancer Patients on Immunotherapy
Next Article in Special Issue
Translation and Validation of the Caffeine Expectancy Questionnaire in Brazil (CaffEQ-BR)
Previous Article in Journal
Markers of Bone Health and Impact of Whey Protein Supplementation in Army Initial Entry Training Soldiers: A Double-Blind Placebo-Controlled Study
Previous Article in Special Issue
The Effect of Caffeine on the Risk and Progression of Parkinson’s Disease: A Meta-Analysis
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Nutrients
Volume 12
Issue 8
10.3390/nu12082228
nutrients-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Interactions of Habitual Coffee Consumption by Genetic Polymorphisms with the Risk of Prediabetes and Type 2 Diabetes Combined

by
Taiyue Jin
1,
Jiyoung Youn
1,
An Na Kim
1,
Moonil Kang
2,
Kyunga Kim
3,4,
Joohon Sung
5 and
Jung Eun Lee
1,6,*
1
Department of Food and Nutrition, College of Human Ecology, Seoul National University, Seoul 08826, Korea
2
Institute of Health and Environment, Graduate School of Public Health, Seoul National University, Seoul 08826, Korea
3
Statistics and Data Center, Research Institute for Future Medicine, Samsung Medical Center, Seoul 03181, Korea
4
Department of Digital Health, Samsung Advanced Institute for Health Sciences & Technology, Sungkyunkwan University, Seoul 06351, Korea
5
Department of Epidemiology, Graduate School of Public Health, Seoul National University, Seoul 08826, Korea
6
The Research Institute of Human Ecology, Seoul National University, Seoul 08826, Korea
*
Author to whom correspondence should be addressed.
Nutrients 2020, 12(8), 2228; https://doi.org/10.3390/nu12082228
Submission received: 2 July 2020 / Revised: 23 July 2020 / Accepted: 23 July 2020 / Published: 26 July 2020
(This article belongs to the Special Issue Coffee and Caffeine Consumption for Human Health)
Browse Figures
Versions Notes

Abstract

:
Habitual coffee consumption and its association with health outcomes may be modified by genetic variation. Adults aged 40 to 69 years who participated in the Korea Association Resource (KARE) study were included in this study. We conducted a genome-wide association study (GWAS) on coffee consumption in 7868 Korean adults, and examined whether the association between coffee consumption and the risk of prediabetes and type 2 diabetes combined was modified by the genetic variations in 4054 adults. In the GWAS for coffee consumption, a total of five single nucleotide polymorphisms (SNPs) located in 12q24.11-13 (rs2074356, rs11066015, rs12229654, rs11065828, and rs79105258) were selected and used to calculate weighted genetic risk scores. Individuals who had a larger number of minor alleles for these five SNPs had higher genetic risk scores. Multivariate logistic regression models were used to estimate the odds ratios (ORs) and 95% confidence intervals (95% CIs) to examine the association. During the 12 years of follow-up, a total of 2468 (60.9%) and 480 (11.8%) participants were diagnosed as prediabetes or type 2 diabetes, respectively. Compared with non-black-coffee consumers, the OR (95% CI) for ≥2 cups/day by black-coffee consumers was 0.61 (0.38–0.95; p for trend = 0.023). Similarly, sugared coffee showed an inverse association. We found a potential interaction by the genetic variations related to black-coffee consumption, suggesting a stronger association among individuals with higher genetic risk scores compared to those with lower scores; the ORs (95% CIs) were 0.36 (0.15–0.88) for individuals with 5 to 10 points and 0.87 (0.46–1.66) for those with 0 points. Our study suggests that habitual coffee consumption was related to genetic polymorphisms and modified the risk of prediabetes and type 2 diabetes combined in a sample of the Korean population. The mechanisms between coffee-related genetic variation and the risk of prediabetes and type 2 diabetes combined warrant further investigation.

1. Introduction

The incidence rate and prevalence of type 2 diabetes have steadily increased in Asian populations. The International Diabetes Federation (IDF) estimated that 163 million people (35.2% of global diabetic population) in the Western Pacific region had prevalent type 2 diabetes in 2019 [1], contributing the most to type 2 diabetes in the world. In Korea, the prevalence of type 2 diabetes has increased from 6.9% in 1998 to 10.8% in 2017 [2,3]. Additionally, type 2 diabetes contributed to 17.1% of the total deaths in Korea in 2018 [4].
Coffee consumption has been suggested to lower several chronic diseases, including type 2 diabetes [5], metabolic syndrome [6], coronary heart disease [7], liver disorders [8], and several types of cancers [9]. The bioactive compounds in coffee, such as caffeine and chlorogenic acids, have been investigated as potential compounds that lower the risk of type 2 diabetes. Caffeine has been shown to stimulate the metabolic rate [10,11], and its thermogenic effect has been hypothesized to decrease the risk of metabolic disease development. Antioxidants, including chlorogenic acids, commonly found in coffee, have also been highlighted as a preventing factor for type 2 diabetes by inhibiting the generation of free radicals and removing hyperglycemia-induced oxidative stress [12,13,14].
A heritability study on caffeine [15] and genome-wide association studies (GWASs) have suggested that coffee consumption behavior may be linked to genetic polymorphisms. The first genome-wide meta-analysis of coffee consumption was conducted in a European population and identified two independent loci, rs4410790 nears AHR and rs2470893 between CYP1A1 and CYP1A2, which are caffeine metabolism-related genes [16]. Additional European/Caucasian GWASs also discovered single nucleotide polymorphisms (SNPs) located in AHR, CYP1A1 and CYP1A2 as well as ABCG2, POR, BDNF, SLC6A4, GCKR and MLXIPL [17], near NRCAM or ULK3 [18]. A Japanese GWAS, the first GWAS of coffee consumption in Asia, identified 24 SNPs on chromosome 12, showing rs2074356 in HECTD4 as the strongest significant SNP [19]. Another Japanese GWAS found two loci located in CUX2 (rs7910258) and AHR (rs10251701) [20]. The few GWASs in Asia warrant further investigation of the coffee-related genetic polymorphisms in Asian populations because there has been an increase in coffee consumption and type 2 diabetes in Asia. Because coffee consumption has the potential to prevent type 2 diabetes [21], it is important to investigate whether coffee consumption is linked to a lower risk of type 2 diabetes and whether this association is modified by genetic variations common in Asian populations.
The objective of this study was to identify genetic polymorphisms associated with habitual coffee consumption and examine whether the association between coffee consumption and the risk of prediabetes and type 2 diabetes combined was modified by these genetic variants in Korean population.

2. Materials and Methods

2.1. Study Population

Participants were recruited from the Korea Association Resource (KARE) study, which is part of the Korean Genome and Epidemiology Study (KoGES), a community-based cohort study. A total of 5012 and 5018 participants aged 40–69 years were enrolled from Ansan and Ansung, respectively, between 2001 and 2002. Socio-demographic status, anthropometric indices, dietary lifestyle, physical activity, and disease history were examined at baseline and follow-up phase every 2 years. In this study, we included data from baseline to the 6th follow-up (2001–2014). The details of the KoGES project have been described elsewhere [22].
Among the 10,030 participants of this study, DNA samples from a total of 8840 participants were genotyped at baseline. We excluded participants who did not have genetic information (n = 1190); those who were diagnosed with a type 2 diabetes by physicians or had been treated with oral hypoglycemic medication or insulin therapy at baseline (n = 619); those who had a history of cardiovascular disease including myocardial infarction and stroke (n = 150) or cancer (n = 20) at baseline; and those who did not provide daily coffee consumption frequency (n = 129) or amount (n = 54). In summary, 7868 participants (3685 men and 4183 women) were included in the GWAS (Table S1).
When we examined the association between habitual coffee consumption and the incident risk of prediabetes and type 2 diabetes combined, we further excluded participants who had prevalent type 2 diabetes (n = 666) or prediabetes (n = 2099) determined by the fasting plasma glucose (FPG) test, a 2-h oral glucose tolerance test (OGTT) and a hemoglobin A1c (HbA1c) test at baseline. We also excluded those who attended neither fifth nor sixth follow-up examination (n = 1049). As a result, a total of 4054 participants (1904 men and 2150 women) were included in the association analysis. This study was approved by the Institutional Review Board of Seoul National University (IRB No. E1911/003-011).

2.2. Dietary Assessment

A 103-item semi-quantitative food frequency questionnaire (FFQ) for the KoGES was used to assess the habitual coffee consumption at baseline. The validity of the FFQ was evaluated by 3-day diet records [23,24]. Participants were asked to answer the frequency of coffee consumption over the last year (almost none, 1 time/month, 2–3 times/month, 1–2 times/week, 3–4 times/week, 5–6 times/week, 1 time/day, 2 times/day or 3 times/day) as well as the amount of coffee (0.5, 1 or 2 cups per time). The amount of sugar added in coffee was also investigated (1, 2 or 3 teaspoons per time). We multiplied the frequency and the amount of the daily consumption of coffee to obtain the number of cups consumed per day. Coffee consumption was categorized into non-coffee consumers, <1 cup/day, 1 to <2 cups/day and ≥2 cups/day. Among coffee consumers, participants who consumed coffee without sugar were defined as black-coffee consumers and those who consumed coffee with sugar as sugared-coffee consumers.
We calculated the amount of caffeine consumption of each participant by using the caffeine database published by the Korea Ministry of Food and Drug Safety (MFDS) [25] or the United States Department of Agriculture (USDA) [26] or by contacting the product companies. The foods considered to contain caffeine were as follows: beverages (coffee, tea and soft drinks) and foods made with cocoa (chocolate, chocolate candies, chocolate ice cream, chocolate milk, chocolate pies, and chocolate cakes).

2.3. Ascertainment of Cases

Participants underwent a 2-h 75 g OGTT and HbA1c test at baseline and each follow-up. The concentrations of FPG and 2-h plasma glucose were measured by the hexokinase method (ADVIA 1650; Bayer, Berkley, MI, USA). Incident cases of prediabetes and type 2 diabetes were defined according to the American Diabetes Association (ADA) criteria [27]. Incident type 2 diabetes was ascertained if a participant had one of the following during the follow-up examination: had been diagnosed with a type 2 diabetes by physicians or had been treated with oral hypoglycemic medication or insulin therapy; had an FPG ≥ 126 mg/dL (≥7.0 mmol/L); had a 2-h plasma glucose ≥ 200 mg/dL (≥11.1 mmol/L); or had a HbA1c ≥ 6.5%. Prediabetes incidence was defined as an FPG ≥ 100 mg/dL (≥5.6 mmol/L) or 2-h plasma glucose ≥ 140 mg/dL (≥7.8 mmol/L) or HbA1c ≥ 6.0%.

2.4. Genotyping and Quality Control

DNA samples were genotyped using the Affymetrix Genome-Wide Human SNP Array 5.0 (Affymetrix, Santa Clara, CA, USA). The Bayesian Robust Linear Modeling using Mahalanobis Distance (BRLMM) Genotyping Algorithm was used for the genotype calling of 500,568 SNPs [28]. After quality control filtering, a total of 352,228 SNPs remained for analysis. Details about quality control criteria have been described elsewhere [29].

2.5. GWAS on Coffee Consumption

We identified SNPs related to habitual coffee consumption in the GWAS. Participants were grouped into non-coffee consumers and coffee consumers, and this group information was used as a binary outcome variable in a logistic regression model for each SNP. We adjusted for age (years, continuous), sex, and alcohol consumption (g/day, continuous) in the GWAS. We also adjusted for baseline BMI in the secondary GWAS. In addition, we replicated a GWAS using a continuous variable of coffee consumption. A GWAS on continuous caffeine intake was also conducted. A box-cox transformed coffee or caffeine intake and a linear regression model were used for continuous analysis. A Manhattan plot and quantile-quantile (Q-Q) plot were generated, and the inflation factor (λ) was calculated. The GWAS was performed using PLINK version 1.07, and SNPs with a p-value < 1 × 10−5 were considered suggestive significant. A regional association plot of the significant SNPs was generated using the LocusZoom program (http://locuszoom.org). The linkage disequilibrium (LD) between the significant SNPs was examined using the HaploView program [30].

2.6. Statistical Analysis

To calculate the genetic risk scores, each coffee-related SNP identified from the GWAS was assigned 0, 1, or 2 according to the number of minor alleles and then weighted by its relative effective size (β coefficient obtained from the GWAS). The genetic risk scores were calculated as following the equation: 5 × (β1 × SNP1 + β2 × SNP2 + β3 × SNP3 + β4 × SNP4 + β5 × SNP5)/(β1 + β2 + β3 + β4 + β5) [31]. The genetic risk scores ranged from 0 to 10 points, and the median value was 5 points among those with at least one minor allele. Participants were categorized into 3 groups: 0 point, 0.1 to <5 points, and 5 to 10 points. We examined the association between coffee consumption and the risk of prediabetes and type 2 diabetes combined in black coffee and sugared coffee. Non-coffee consumers and sugared-coffee consumers combined were regarded as reference group for black-coffee consumers. For sugared-coffee consumers, non-coffee consumers and black-coffee consumers combined were regarded as a reference group. Odds ratios (ORs) and 95% confidence intervals (CIs) for the association between coffee consumption and the risk of prediabetes and type 2 diabetes combined were calculated using multivariate logistic regression models among men and women combined or separately. The median cups of coffee consumed per day were assigned to each coffee group and used to test the linear trends. We also examined whether the association varied by genetic risk scores. The interaction analysis was performed by comparing the models with or without interaction term using a likelihood ratio test. All of the analyses were adjusted for age (years, continuous), sex (for men and women combined), body mass index (BMI, <23, 23 to <25, 25 to <30 and ≥30 kg/m2), smoking status (never smokers, ≤10 and >10 pack-years for black-coffee consumers; never smokers, ≤10, 10.1 to ≤20, 20.1 to ≤30 and >30 pack-years for sugared-coffee consumers), alcohol consumption (non-drinkers, ≤5, 5.1 to ≤10 and >10 g/day for black-coffee consumers; non-drinkers, ≤5, 5.1 to ≤10, 10.1 to ≤20 and >20 g/day for sugared-coffee consumers), family history of type 2 diabetes (yes or no), and total energy intake (kcal/day, continuous). Additionally, we also adjusted for the amount of sugar added in coffee when analyzing the association between black coffee consumption and the risk of prediabetes and type 2 diabetes combined (0, ≤5, 5.1 to ≤10, 10.1 to ≤15 and >15 g/day). SAS version 9.4 (SAS Institute, Cary, NC, USA) was used for all analyses, and p-value < 0.05 in two-sided tests was defined as a significant difference.

3. Results

3.1. Baseline Characteristics

Of the 4054 participants in the association analysis, a total of 480 (11.8%) and 2468 (60.9%) participants developed type 2 diabetes and prediabetes during the 12-year follow-up period, respectively. Table 1 shows the baseline characteristics of the study population according to habitual coffee consumption. Compared with non-coffee consumers, participants who consumed either black coffee or sugared coffee were more likely to be younger, smoke, drink alcohol, and have a higher BMI.

3.2. Genetic Polymorphisms Associated with Coffee Consumption

In the GWAS of habitual coffee consumption, a total of 18 SNPs located in 12q24 achieved a suggestive significance (p < 1 × 10−5) (Table S2). When we additionally adjusted for baseline BMI, the same 18 significant SNPs were identified as well (Table S3). SNPs identified based on habitual coffee consumption were the same when we additionally conducted a GWAS for continuous caffeine intake.
Figure 1 and Figure 2 show the Manhattan plot and regional association plot, respectively. Figure 3 shows the Q-Q plot, and the genomic inflation factor (λ) of the GWAS was 1.0103, suggesting that the study population structure was well-adjusted [32]. Among the SNPs achieved suggestive significance at p < 1 × 10−5, we selected SNPs as follows, which were used later for calculating the genetic risk scores. First, pairwise correlations were examined and high correlation was declared when r2 > 0.8 (Figure 4). Then, we excluded all imputed SNPs having a high correlation with any genotyped SNP, and selected the SNPs among highly-correlated genotyped SNPs. As a result, three genotyped SNPs (rs2074356, rs11066015, and rs12229654) and two imputed SNPs (rs11065828 and rs79105258) were considered as coffee-related SNPs and used for calculating the genetic risk scores. Table 2 presents the information on these five SNPs related to coffee consumption. The SNPs were located in genes HECTD4, ACAD10, MYL2, and CUX2, and the minor allele frequency ranged from 0.143 to 0.172 in our population.

3.3. Coffee Consumption and the Risk of Prediabetes and Type 2 Diabetes Combined

We examined the association between black coffee consumption and the risk of prediabetes and type 2 diabetes combined (Table 3). Compared with non-black-coffee consumers, participants who consumed ≥2 cups/day of black coffee had a 39% lower risk of prediabetes and type 2 diabetes combined among men and women combined (95% CI = 0.38–0.95; p for trend = 0.023). When we separated men and women, compared with non-black-coffee consumers, participants who consumed ≥2 cups/day of black coffee had a 54% lower risk of prediabetes and type 2 diabetes combined among men (95% CI = 0.23–0.94; p for trend = 0.026). Although the association was not statistically significant among women, an inverse trend was observed (OR = 0.74; 95% CI = 0.41–1.34).
Table 4 presents the association between sugared coffee consumption and the risk of prediabetes and type 2 diabetes combined. The ORs (95% CIs) of prediabetes and type 2 diabetes combined comparing sugared-coffee consumers with non-sugared-coffee consumers were 0.73 (0.60–0.89; p for trend = 0.005) for men and women combined, 0.71 (0.52–0.97; p for trend = 0.015) for men, and 0.75 (0.57–0.99; p for trend = 0.080) for women.

3.4. Associations Modified by Genetic Risk Scores

We analyzed whether the association between habitual coffee consumption and the risk of prediabetes and type 2 diabetes combined varied by genetic polymorphisms. For black-coffee consumers, we found that an inverse association was more pronounced among individuals with high genetic risk scores, but the interaction was not statistically significant (p for interaction = 0.261) (Table 5). The ORs (95% CIs) for ≥2 cups/day of black coffee vs. non-black-coffee consumers were 0.87 (0.46–1.66) for participants with 0 points, 0.49 (0.17–1.44) for those with 0.1 to <5 points, and 0.36 (0.15–0.88) for those with 5 to 10 points for their genetic risk scores. For sugared coffee consumption, the associations with the risk of prediabetes and type 2 diabetes combined were similar across genetic risk scores (p for interaction = 0.608) (Table 6). When we conducted an interaction analysis for each coffee-related SNP for black coffee consumption, the association between black coffee consumption and the risk of prediabetes and type 2 diabetes combined was inverse for the minor allele of each SNP (Table 7).

4. Discussion

Our first GWAS of coffee consumption in the Korean population identified five SNPs (rs2074356 in HECTD4, rs11066015 in ACAD10, rs12229654 in MYL2, rs11065828 and rs79105258 in CUX2) related to habitual coffee consumption. Compared with non-coffee consumers, the risk of prediabetes and type 2 diabetes being combined was inversely associated with habitual coffee consumption, either black coffee or sugared coffee. Individuals with black coffee consumption had a lower risk of prediabetes and type 2 diabetes combined compared with non-black-coffee consumers among those with multiple minor alleles for these five SNPs.
The significant SNPs discovered to be related to habitual coffee consumption in our GWAS were all introns. Although the introns were noncoding regions of genes, they may affect the transcription rate and translation efficiency, further regulating gene expression [33]. Recent GWASs have identified several loci on the AHR, CYP1A1, and CYP1A2 genes associated with coffee consumption [16,17]. However, most of the GWASs on coffee consumption were conducted in European populations, only two GWASs, to our knowledge, reported SNPs related to coffee consumption in Asian populations. A GWAS in the Japan Multi-Institutional Collaborative Cohort (J-MICC) study, the first GWAS on coffee consumption in Asia, found that rs2074356 located in 12q24 was most strongly associated with habitual coffee consumption (p = 2.2 × 10−6) [19]. Similarly, in our GWAS, a total of 18 SNPs were associated with coffee consumption at p < 1 × 10−5, all of which were found to be related to habitual coffee consumption in the J-MICC study, and the strongest significant variant was rs2074356 as well. Another Japanese coffee GWAS identified two independent loci (rs79105258 in 12q24 and rs10252701 in 7p21) that were associated with coffee consumption [20]. rs79105258 was also selected as a coffee-related variant in our study. In addition to the association with habitual coffee consumption, these SNPs were associated with type 2 diabetes [34], blood glucose levels [35], blood pressure levels [36], and obesity [37]. rs2074356 in HECTD4 was associated with prevalent type 2 diabetes and blood glucose level in a Korean population [34,35]. Three SNPs (rs12229654, rs11066015 and rs2074356) were also identified to be linked with both systolic and diastolic blood pressure in a Japanese GWAS [36].
Previous epidemiologic studies have shown that coffee consumption was inversely associated with the risk of type 2 diabetes. A meta-analysis of 28 cohort and nested case-control studies reported that participants who consumed 5 cups/day of coffee had a 30% lower risk of type 2 diabetes compared with almost non-consumers, and the associations were similar between men and women [5]. Although we found a stronger inverse association among men than among women, further investigation is needed to explore a larger amount of coffee consumption, e.g., 3 or more cups/day, and whether the inverse association holds for women.
The lower risk of type 2 diabetes linked to coffee consumption could be linked to several biological mechanisms. As a main polyphenolic compound in coffee, chlorogenic acids have been shown as inhibitors of hepatic glucose-6-phosphatase, the rate-limiting enzyme of glucose hydrolysis [38]. Reduced hepatic glucose-6-phosphatase may affect the glucose output and thus decrease the blood glucose concentration. In addition, chlorogenic acids act as antioxidants to lower oxidative stress shown in both in vitro and in vivo studies [39]. Additionally, caffeine and magnesium, both of which are commonly found in coffee, have been suggested to have roles in type 2 diabetes prevention by improving insulin resistance. Previous studies suggested that caffeine could improve insulin resistance by stimulating insulin secretion from pancreatic β cells [40]. In addition to insulin secretion, caffeine increases thermogenesis, lipolysis, and β-oxidation [41]. Magnesium supplementation has reduced the development of type 2 diabetes and improved glucose disposal in experimental studies [42,43], and cohort studies have reported a significant inverse association between magnesium intake and type 2 diabetes risk [44].
In this study, we observed that both black coffee and sugared coffee decreased the risk of prediabetes and type 2 diabetes combined for individuals who consumed more than 2 cups of coffee per day. But participants who consumed 1 to <2 cups/day of sugared coffee had a 45% higher risk of prediabetes and type 2 diabetes combined among men alone, but not among women. Although sugar-sweetened beverage intake including carbonated beverages has been positively associated with the risk of type 2 diabetes [45,46], the results of coffee with sugar remained equivocal. In a French cohort study, compared to non-coffee consumers, the incident type 2 diabetes decreased by 40% and 31% among participants consuming more than 1.1 cups of coffee per lunch with or without sugar, respectively [47]. In a small clinical trial, where eight lean, young and healthy adults drank six types of beverages or water 1 h before a potato-based meal, postprandial hyperglycemia, an early abnormality of type 2 diabetes [48], was significantly reduced when they drank sweetened coffee before their meals [49]. Further studies are needed to examine whether the benefit of coffee remains even after adding a small amount of sugar.
To our knowledge, this study is the first GWAS of habitual coffee consumption in a Korean population. The strengths of this study include good ascertainment of prediabetes and type 2 diabetes, adjustment for potential confounding factors, and a 12-year follow-up. The incidence of prediabetes and type 2 diabetes were identified based on the circulating levels of FPG, 2-h plasma glucose, and HbA1c, which could minimize the misclassification. Adjustment for potential confounding factors, including smoking status and alcohol consumption, may enable us to remove the effect of the confounding factors. However, we cannot rule out the possibility that residual confounding factors remained. There are several limitations to our study. First, the rate of revisits to the clinic for the blood draw decreased to 60% in the last sixth follow-up. Therefore, our study may not be representative of the full cohort of KARE study. However, the internal validity may not be impaired as we obtained relatively accurate information on the incidence of prediabetes and type 2 diabetes during the 12-year follow-up. Second, we could not distinguish caffeinated and decaffeinated coffee or boiled and filtered coffee. However, previous studies have shown that the associations between coffee and type 2 diabetes were similar by the amount of caffeine [5] or preparation methods [50]. Third, we were not able to examine high coffee consumption (e.g., 3 or more cups/day) because only a few participants consumed more than 3 cups/day. Fourth, we did not consider medicinal caffeine intake when we performed a GWAS. However, beverages may mainly contribute to daily caffeine intake in Korea.

5. Conclusions

In our study, we conducted a GWAS and discovered five SNPs (rs2074356 in HECTD4, rs11066015 in ACAD10, rs12229654 in MYL2, rs11065828 and rs79105258 in CUX2) associated with habitual coffee consumption in a Korean population. We observed that moderate black coffee and sugared coffee consumption reduced the risk of prediabetes and type 2 diabetes combined. We found that an inverse association was stronger among black-coffee consumers with minor alleles of five SNPs related to coffee consumption compared to those with major alleles. Further Asian GWASs and epidemiological studies are needed to elucidate the effects of coffee-related genetic variation and high coffee consumption on chronic disease risk.

Supplementary Materials

The following are available online at https://www.mdpi.com/2072-6643/12/8/2228/s1. Table S1: Baseline characteristics of the participants included in the GWAS on coffee consumption; Table S2: The significant SNPs discovered from the GWAS on coffee consumption; Table S3: The significant SNPs discovered from the GWAS on coffee consumption after additionally adjusted for BMI.

Author Contributions

Conceptualization, T.J. and J.E.L.; data curation, T.J., J.Y., K.K., J.S. and J.E.L.; formal analysis, T.J., J.Y. and J.E.L.; investigation, T.J., J.Y., A.N.K., M.K., K.K., J.S. and J.E.L.; methodology, T.J., J.Y., M.K., K.K., J.S. and J.E.L.; project administration, J.S. and J.E.L.; supervision, J.E.L.; writing—original draft, T.J. and J.E.L.; writing—review & editing, T.J., J.Y., A.N.K., M.K., K.K., J.S. and J.E.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

This study was conducted with bioresources from National Biobank of Korea, the Center for Disease Control and Prevention, Republic of Korea (KBN-2018-044).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Saeedi, P.; Petersohn, I.; Salpea, P.; Malanda, B.; Karuranga, S.; Unwin, N.; Colagiuri, S.; Guariguata, L.; Motala, A.A.; Ogurtsova, K.; et al. Global and regional diabetes prevalence estimates for 2019 and projections for 2030 and 2045: Results from the International Diabetes Federation Diabetes Atlas, 9th edition. Diabetes Res. Clin. Pract. 2019, 157, 107843. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Ko, B.; Lim, J.; Kim, Y.Z.; Park, H.S. Trends in type 2 diabetes prevalence according to income levels in Korea (1998–2012). Diabetes Res. Clin. Pract. 2016, 115, 137–139. [Google Scholar] [CrossRef] [PubMed]
  3. Shin, J.-Y. Trends in the prevalence and management of diabetes in Korea: 2007–2017. Epidemiol. Health 2019, 41, e2019029. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Statistics Korea. Causes of Death Statistics in 2018; Statistics Korea: Daejeon, Korea, 2019.
  5. Ding, M.; Bhupathiraju, S.N.; Chen, M.; van Dam, R.M.; Hu, F.B. Caffeinated and decaffeinated coffee consumption and risk of type 2 diabetes: A systematic review and a dose-response meta-analysis. Diabetes Care 2014, 37, 569–586. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Shang, F.; Li, X.; Jiang, X. Coffee consumption and risk of the metabolic syndrome: A meta-analysis. Diabetes Metab. 2016, 42, 80–87. [Google Scholar] [CrossRef]
  7. Wu, J.-N.; Ho, S.C.; Zhou, C.; Ling, W.-H.; Chen, W.-Q.; Wang, C.-L.; Chen, Y.-M. Coffee consumption and risk of coronary heart diseases: A meta-analysis of 21 prospective cohort studies. Int. J. Cardiol. 2009, 137, 216–225. [Google Scholar] [CrossRef]
  8. Saab, S.; Mallam, D.; Cox Ii, G.A.; Tong, M.J. Impact of coffee on liver diseases: A systematic review. Liver Int. 2014, 34, 495–504. [Google Scholar] [CrossRef]
  9. Arab, L. Epidemiologic Evidence on Coffee and Cancer. Nutr. Cancer 2010, 62, 271–283. [Google Scholar] [CrossRef] [Green Version]
  10. Acheson, K.J.; Zahorska-Markiewicz, B.; Pittet, P.; Anantharaman, K.; Jéquier, E. Caffeine and coffee: Their influence on metabolic rate and substrate utilization in normal weight and obese individuals. Am. J. Clin. Nutr. 1980, 33, 989–997. [Google Scholar] [CrossRef]
  11. Astrup, A.; Toubro, S.; Cannon, S.; Hein, P.; Breum, L.; Madsen, J. Caffeine: A double-blind, placebo-controlled study of its thermogenic, metabolic, and cardiovascular effects in healthy volunteers. Am. J. Clin. Nutr. 1990, 51, 759–767. [Google Scholar] [CrossRef]
  12. Liang, N.; Kitts, D.D. Role of Chlorogenic Acids in Controlling Oxidative and Inflammatory Stress Conditions. Nutrients 2015, 8, 16. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Ludwig, I.A.; Clifford, M.N.; Lean, M.E.J.; Ashihara, H.; Crozier, A. Coffee: Biochemistry and potential impact on health. Food Funct. 2014, 5, 1695–1717. [Google Scholar] [CrossRef] [PubMed]
  14. Ceriello, A.; Testa, R.; Genovese, S. Clinical implications of oxidative stress and potential role of natural antioxidants in diabetic vascular complications. Nutr. Metab. Cardiovasc. Dis. 2016, 26, 285–292. [Google Scholar] [CrossRef] [PubMed]
  15. Yang, A.; Palmer, A.A.; de Wit, H. Genetics of caffeine consumption and responses to caffeine. Psychopharmacology 2010, 211, 245–257. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. 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] [PubMed]
  17. Cornelis, M.C.; Byrne, E.M.; Esko, T.; Nalls, M.A.; Ganna, A.; Paynter, N.; Monda, K.L.; Amin, N.; Fischer, K.; Renstrom, F.; et al. Genome-wide meta-analysis identifies six novel loci associated with habitual coffee consumption. Mol. Psychiatry 2015, 20, 647–656. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Amin, N.; Byrne, E.; Johnson, J.; Chenevix-Trench, G.; Walter, S.; Nolte, I.M.; kConFab, I.; Vink, J.M.; Rawal, R.; Mangino, M.; et al. Genome-wide association analysis of coffee drinking suggests association with CYP1A1/CYP1A2 and NRCAM. Mol. Psychiatry 2012, 17, 1116–1129. [Google Scholar] [CrossRef]
  19. 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]
  20. 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]
  21. Carlström, M.; Larsson, S.C. Coffee consumption and reduced risk of developing type 2 diabetes: A systematic review with meta-analysis. Nutr. Rev. 2018, 76, 395–417. [Google Scholar] [CrossRef]
  22. Kim, Y.; Han, B.-G.; KoGES Group. Cohort Profile: The Korean Genome and Epidemiology Study (KoGES) Consortium. Int. J. Epidemiol. 2017, 46, e20. [Google Scholar] [CrossRef] [PubMed]
  23. Ahn, Y.; Kwon, E.; Shim, J.E.; Park, M.K.; Joo, Y.; Kimm, K.; Park, C.; Kim, D.H. Validation and reproducibility of food frequency questionnaire for Korean genome epidemiologic study. Eur. J. Clin. Nutr. 2007, 61, 1435. [Google Scholar] [CrossRef] [PubMed]
  24. Kim, J.; Kim, Y.; Ahn, Y.O.; Paik, H.Y.; Ahn, Y.; Tokudome, Y.; Hamajima, N.; Inoue, M.; Tajima, K. Development of a food frequency questionnaire in Koreans. Asia Pac. J. Clin. Nutr. 2003, 12, 243–250. [Google Scholar] [PubMed]
  25. Ministry of Food and Drug Safety (KR). Food Composition Database. Available online: http://www.foodsafetykorea.go.kr/fcdb/ (accessed on 5 March 2020).
  26. United States Department of Agriculture (USDA), A.R.S. Abridged List Ordered by Nutrient Content in Household Measure, Nutrients: Caffeine(mg). Available online: https://www.nal.usda.gov/sites/www.nal.usda.gov/files/caffeine.pdf (accessed on 17 February 2020).
  27. American Diabetes Association. Diagnosis and Classification of Diabetes Mellitus. Diabetes Care 2014, 37, S81. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Rabbee, N.; Speed, T.P. A genotype calling algorithm for affymetrix SNP arrays. Bioinformatics 2005, 22, 7–12. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Cho, Y.S.; Go, M.J.; Kim, Y.J.; Heo, J.Y.; Oh, J.H.; Ban, H.J.; Yoon, D.; Lee, M.H.; Kim, D.J.; Park, M.; et al. A large-scale genome-wide association study of Asian populations uncovers genetic factors influencing eight quantitative traits. Nat. Genet. 2009, 41, 527–534. [Google Scholar] [CrossRef] [PubMed]
  30. Barrett, J.C.; Fry, B.; Maller, J.; Daly, M.J. Haploview: Analysis and visualization of LD and haplotype maps. Bioinformatics 2005, 21, 263–265. [Google Scholar] [CrossRef] [Green Version]
  31. Wang, T.; Huang, T.; Heianza, Y.; Sun, D.; Zheng, Y.; Ma, W.; Jensen, M.K.; Kang, J.H.; Wiggs, J.L.; Pasquale, L.R.; et al. Genetic Susceptibility, Change in Physical Activity, and Long-term Weight Gain. Diabetes 2017, 66, 2704–2712. [Google Scholar] [CrossRef] [Green Version]
  32. Zeng, P.; Zhao, Y.; Qian, C.; Zhang, L.; Zhang, R.; Gou, J.; Liu, J.; Liu, L.; Chen, F. Statistical analysis for genome-wide association study. J. Biomed. Res. 2015, 29, 285–297. [Google Scholar] [CrossRef] [Green Version]
  33. Shaul, O. How introns enhance gene expression. Int. J. Biochem. Cell Biol. 2017, 91, 145–155. [Google Scholar] [CrossRef]
  34. Go, M.J.; Hwang, J.Y.; Kim, Y.J.; Hee Oh, J.; Kim, Y.J.; Heon Kwak, S.; Soo Park, K.; Lee, J.; Kim, B.J.; Han, B.G.; et al. New susceptibility loci in MYL2, C12orf51 and OAS1 associated with 1-h plasma glucose as predisposing risk factors for type 2 diabetes in the Korean population. J. Hum. Genet. 2013, 58, 362–365. [Google Scholar] [CrossRef] [PubMed]
  35. Go, M.J.; Hwang, J.Y.; Park, T.J.; Kim, Y.J.; Oh, J.H.; Kim, Y.J.; Han, B.G.; Kim, B.J. Genome-wide association study identifies two novel Loci with sex-specific effects for type 2 diabetes mellitus and glycemic traits in a korean population. Diabetes Metab. J. 2014, 38, 375–387. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Yamada, Y.; Sakuma, J.; Takeuchi, I.; Yasukochi, Y.; Kato, K.; Oguri, M.; Fujimaki, T.; Horibe, H.; Muramatsu, M.; Sawabe, M.; et al. Identification of polymorphisms in 12q24.1, ACAD10, and BRAP as novel genetic determinants of blood pressure in Japanese by exome-wide association studies. Oncotarget 2017, 8, 43068–43079. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. 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]
  38. Arion, W.J.; Canfield, W.K.; Ramos, F.C.; Schindler, P.W.; Burger, H.J.; Hemmerle, H.; Schubert, G.; Below, P.; Herling, A.W. Chlorogenic acid and hydroxynitrobenzaldehyde: New inhibitors of hepatic glucose 6-phosphatase. Arch. Biochem. Biophys. 1997, 339, 315–322. [Google Scholar] [CrossRef]
  39. Sato, Y.; Itagaki, S.; Kurokawa, T.; Ogura, J.; Kobayashi, M.; Hirano, T.; Sugawara, M.; Iseki, K. In vitro and in vivo antioxidant properties of chlorogenic acid and caffeic acid. Int. J. Pharm. 2011, 403, 136–138. [Google Scholar] [CrossRef]
  40. Keijzers, G.B.; De Galan, B.E.; Tack, C.J.; Smits, P. Caffeine Can Decrease Insulin Sensitivity in Humans. Diabetes Care 2002, 25, 364. [Google Scholar] [CrossRef] [Green Version]
  41. Greenberg, J.A.; Boozer, C.N.; Geliebter, A. Coffee, diabetes, and weight control. Am. J. Clin. Nutr. 2006, 84, 682–693. [Google Scholar] [CrossRef] [Green Version]
  42. Schulze, M.B.; Schulz, M.; Heidemann, C.; Schienkiewitz, A.; Hoffmann, K.; Boeing, H. Fiber and Magnesium Intake and Incidence of Type 2 Diabetes: A Prospective Study and Meta-analysis. Arch. Intern. Med. 2007, 167, 956–965. [Google Scholar] [CrossRef] [Green Version]
  43. Balon, T.W.; Gu, J.L.; Tokuyama, Y.; Jasman, A.P.; Nadler, J.L. Magnesium supplementation reduces development of diabetes in a rat model of spontaneous NIDDM. Am. J. Physiol.-Endocrinol. Metab. 1995, 269, E745–E752. [Google Scholar] [CrossRef]
  44. Larsson, S.C.; Wolk, A. Magnesium intake and risk of type 2 diabetes: A meta-analysis. J. Intern. Med. 2007, 262, 208–214. [Google Scholar] [CrossRef] [PubMed]
  45. Drouin-Chartier, J.P.; Zheng, Y.; Li, Y.; Malik, V.; Pan, A.; Bhupathiraju, S.N.; Tobias, D.K.; Manson, J.E.; Willett, W.C.; Hu, F.B. Changes in Consumption of Sugary Beverages and Artificially Sweetened Beverages and Subsequent Risk of Type 2 Diabetes: Results from Three Large Prospective U.S. Cohorts of Women and Men. Diabetes Care 2019, 42, 2181–2189. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Schulze, M.B.; Manson, J.E.; Ludwig, D.S.; Colditz, G.A.; Stampfer, M.J.; Willett, W.C.; Hu, F.B. Sugar-Sweetened Beverages, Weight Gain, and Incidence of Type 2 Diabetes in Young and Middle-Aged Women. JAMA 2004, 292, 927–934. [Google Scholar] [CrossRef] [PubMed]
  47. Sartorelli, D.S.; Fagherazzi, G.; Balkau, B.; Touillaud, M.S.; Boutron-Ruault, M.-C.; de Lauzon-Guillain, B.; Clavel-Chapelon, F. Differential effects of coffee on the risk of type 2 diabetes according to meal consumption in a French cohort of women: The E3N/EPIC cohort study. Am. J. Clin. Nutr. 2010, 91, 1002–1012. [Google Scholar] [CrossRef] [Green Version]
  48. American Diabetes Association. Postprandial Blood Glucose. Diabetes Care 2001, 24, 775. [Google Scholar] [CrossRef] [Green Version]
  49. Louie, J.C.; Atkinson, F.; Petocz, P.; Brand-Miller, J.C. Delayed effects of coffee, tea and sucrose on postprandial glycemia in lean, young, healthy adults. Asia Pac. J. Clin. Nutr. 2008, 17, 657–662. [Google Scholar]
  50. Tuomilehto, J.; Hu, G.; Bidel, S.; Lindström, J.; Jousilahti, P. Coffee Consumption and Risk of Type 2 Diabetes Mellitus Among Middle-aged Finnish Men and Women. JAMA 2004, 291, 1213–1219. [Google Scholar] [CrossRef] [Green Version]
Nutrients 12 02228 g001 550
Figure 1. Manhattan plot of the genome-wide association study (GWAS) on coffee consumption. Each color represents a different chromosome. The strongest significant SNP was rs2074356 in chromosome 12 (p-value = 6.62 × 10−8).
Figure 1. Manhattan plot of the genome-wide association study (GWAS) on coffee consumption. Each color represents a different chromosome. The strongest significant SNP was rs2074356 in chromosome 12 (p-value = 6.62 × 10−8).
Nutrients 12 02228 g001
Nutrients 12 02228 g002 550
Figure 2. Regional association plot of the 18 significant single nucleotide polymorphisms (SNPs) discovered from the GWAS on coffee consumption. The strongest significant SNP, rs2074356, was shown in purple, and the gray dots represent chromosomal positions of other SNPs near the rs2074356.
Figure 2. Regional association plot of the 18 significant single nucleotide polymorphisms (SNPs) discovered from the GWAS on coffee consumption. The strongest significant SNP, rs2074356, was shown in purple, and the gray dots represent chromosomal positions of other SNPs near the rs2074356.
Nutrients 12 02228 g002
Nutrients 12 02228 g003 550
Figure 3. Quantile-quantile (Q-Q) plot of the GWAS on coffee consumption.
Figure 3. Quantile-quantile (Q-Q) plot of the GWAS on coffee consumption.
Nutrients 12 02228 g003
Nutrients 12 02228 g004 550
Figure 4. Linkage disequilibrium (LD) plot of the 18 significant SNPs discovered from the GWAS on coffee consumption. Values shown in the red boxes represent the LD (r2) between the SNPs.
Figure 4. Linkage disequilibrium (LD) plot of the 18 significant SNPs discovered from the GWAS on coffee consumption. Values shown in the red boxes represent the LD (r2) between the SNPs.
Nutrients 12 02228 g004
Table
Table 1. Baseline characteristics of study population according to habitual coffee consumption.
Table 1. Baseline characteristics of study population according to habitual coffee consumption.
Non-Coffee Consumers (n = 864)Black-Coffee Consumers
<1 Cup/Day (n = 100)1 to <2 Cups/Day (n = 70)≥2 Cups/Day (n = 105)
Age, mean ± SD (years)53.3 ± 8.749.7 ± 7.547.4 ± 6.946.9 ± 5.7
Sex, n (%)
Men299 (34.6)32 (32.0)21 (30.0)47 (44.8)
Women565 (65.4)68 (68.0)49 (70.0)58 (55.2)
BMI, mean ± SD (kg/m2)23.8 ± 3.024.3 ± 2.824.6 ± 3.124.7 ± 3.1
Smoking status, n (%)
Never smokers638 (73.8)76 (76.0)49 (70.0)57 (54.3)
Past smokers103 (11.9)9 (9.0)11 (15.7)14 (13.3)
Current smokers123 (14.2)15 (15.0)10 (14.3)34 (32.4)
Alcohol consumption, n (%)
Never drinkers588 (68.1)55 (55.0)44 (62.9)46 (43.8)
≤5 g/day104 (12.0)15 (15.0)12 (17.1)23 (21.9)
5 to ≤10 g/day36 (4.2)8 (8.0)4 (5.7)8 (7.6)
10 to ≤20 g/day46 (5.3)10 (10.0)6 (8.6)11 (10.5)
>20 g/day90 (10.4)12 (12.0)4 (5.7)17 (16.2)
Family history of type 2 diabetes, n (%)
Yes65 (7.5)11 (11.0)5 (7.1)14 (13.3)
No799 (92.5)89 (89.0)65 (92.9)91 (86.7)
Sugar added in coffee, mean ± SD (g/day)0.1 ± 0.90 ± 00 ± 00 ± 0
Total energy intake, mean ± SD (kcal/day)1882.1 ± 762.51981.1 ± 747.61863.9 ± 718.02143.7 ± 1121.4
Non-Coffee Consumers (n = 864)Sugared-coffee Consumers
<1 Cup/Day (n = 892)1 to <2 Cups/Day (n = 986)≥2 Cups/Day (n = 1037)
Age, mean ± SD (years)53.3 ± 8.751.4 ± 8.550.0 ± 8.449.2 ± 8.0
Sex, n (%)
Men299 (34.6)442 (49.6)425 (43.1)638 (61.5)
Women565 (65.4)450 (50.5)561 (56.9)399 (38.5)
BMI, mean ± SD (kg/m2)23.8 ± 3.024.4 ± 2.924.3 ± 2.924.3 ± 3.0
Smoking status, n (%)
Never smokers638 (73.8)555 (62.2)626 (63.5)458 (44.2)
Past smokers103 (11.9)147 (16.5)129 (13.1)188 (18.1)
Current smokers123 (14.2)190 (21.3)231 (23.4)391 (37.7)
Alcohol consumption, n (%)
Never drinkers588 (68.1)435 (48.8)488 (49.5)459 (44.3)
≤5 g/day104 (12.0)166 (18.6)196 (19.9)204 (19.7)
5 to ≤10 g/day36 (4.2)77 (8.6)59 (6.0)85 (8.2)
10 to ≤20 g/day46 (5.3)90 (10.1)88 (8.9)91 (8.8)
>20 g/day90 (10.4)124 (13.9)155 (15.7)198 (19.1)
Family history of type 2 diabetes, n (%)
Yes65 (7.5)104 (11.7)91 (9.2)96 (9.3)
No799 (92.5)788 (88.3)895 (90.8)941 (90.7)
Sugar added in coffee, mean ± SD (g/day)0.1 ± 0.91.6 ± 1.54.9 ± 1.211.9 ± 4.4
Total energy intake, mean ± SD (kcal/day)1882.1 ± 762.51881.2 ± 638.61969.5 ± 596.22115.8 ± 749.4
Abbreviations: SD, standard deviation; BMI, body mass index.
Table
Table 2. SNPs related to coffee consumption (p-value < 1 × 10−5).
Table 2. SNPs related to coffee consumption (p-value < 1 × 10−5).
ChrSNPPositionGeneAlleles 1MAF 2MAF 3Beta 4p-Value 5
12rs2074356112,645,401HECTD4A/G0.1490.1290.31856.62 × 10−8
12rs11066015112,168,009ACAD10A/G0.1720.1760.24697.79 × 10−6
12rs12229654111,414,461MYL2G/T0.1430.1590.28671.49 × 10−6
12rs11065828111,629,389CUX2A/C0.1720.2140.26541.26 × 10−6
12rs79105258111,718,231CUX2A/C0.1560.2160.29123.89 × 10−7
Abbreviations: Chr, chromosome; SNP, single nucleotide polymorphism; MAF, minor allele frequency. 1 Alleles are presented as minor allele/major allele; 2 minor allele frequency in this study population; 3 minor allele frequency in the 1000Genomes, East Asian; 4 the beta (β) coefficient was obtained from the GWAS; 5 the p-value was calculated using a Wald test from logistic regression model adjusted for age (years; continuous), sex, and alcohol consumption (g/day; continuous).
Table
Table 3. Multivariate-adjusted ORs and 95% CIs for the risk of prediabetes and type 2 diabetes combined according to black coffee consumption.
Table 3. Multivariate-adjusted ORs and 95% CIs for the risk of prediabetes and type 2 diabetes combined according to black coffee consumption.
Non-Black-Coffee Consumers 1Black-Coffee Consumersp for Trend
<1 Cup/Day1 to <2 Cups/Day≥2 Cups/Day
Men and women combined 0.023
Case/total2759/377973/10047/7069/105
ORs (95% CIs) 2Reference0.96 (0.59, 1.55)0.75 (0.44, 1.28)0.61 (0.38, 0.95)
Men 0.026
Case/total1394/180426/3215/2132/47
ORs (95% CIs) 2Reference1.18 (0.45, 3.15)0.69 (0.25, 1.93)0.46 (0.23, 0.94)
Women 0.271
Case/total1365/197547/6832/4937/58
ORs (95% CIs) 2Reference0.91 (0.52, 1.59)0.78 (0.41, 1.47)0.74 (0.41, 1.34)
Abbreviations: OR, odds ratio; 95% CI, 95% confidence interval. 1 For black-coffee consumers, non-coffee consumers and sugared-coffee consumers combined were regarded as reference group; 2 the ORs (95% CIs) were adjusted for age (years, continuous), sex (for men and women combined), body mass index (BMI, <23, 23 to <25, 25 to <30 and ≥30 kg/m2), smoking status (never smokers, ≤10 and >10 pack-years), alcohol consumption (non-drinkers, ≤5, 5.1 to ≤10 and >10 g/day), family history of type 2 diabetes (yes or no), total energy intake (kcal/day, continuous), and the amount of sugar added in coffee (0, ≤5, 5.1 to ≤10, 10.1 to ≤15 and >15 g/day).
Table
Table 4. Multivariate-adjusted ORs and 95% CIs for the risk of prediabetes and type 2 diabetes combined according to sugared coffee consumption.
Table 4. Multivariate-adjusted ORs and 95% CIs for the risk of prediabetes and type 2 diabetes combined according to sugared coffee consumption.
Non-Sugared-Coffee Consumers 1Sugared-Coffee Consumersp for Trend
<1 Cup/Day1 to <2 Cups/Day≥2 Cups/Day
Men and women combined 0.005
Case/total834/1139644/892749/986721/1037
ORs (95% CIs) 2Reference0.84 (0.68, 1.03)1.11 (0.90, 1.35)0.73 (0.60, 0.89)
Men 0.015
Case/total310/399335/442357/425465/638
ORs (95% CIs) 2Reference0.83 (0.59, 1.15)1.45 (1.01, 2.08)0.71 (0.52, 0.97)
Women 0.080
Case/total 524/740309/450392/561256/399
ORs (95% CIs) 2Reference0.84 (0.64, 1.09)0.96 (0.75, 1.23)0.75 (0.57, 0.99)
Abbreviations: OR, odds ratio; 95% CI, 95% confidence interval. 1 For sugared-coffee consumers, non-coffee consumers and black-coffee consumers combined were regarded as reference group; 2 the ORs (95% CIs) were adjusted for age (years, continuous), sex (for men and women combined), body mass index (BMI, <23, 23 to <25, 25 to <30 and ≥30 kg/m2), smoking status (never smokers, ≤10, 10.1 to ≤20, 20.1 to ≤30 and >30 pack-years), alcohol consumption (non-drinkers, ≤5, 5.1 to ≤10, 10.1 to ≤20 and >20 g/day), family history of type 2 diabetes (yes or no), and total energy intake (kcal/day, continuous).
Table
Table 5. Multivariate-adjusted ORs and 95% CIs for the risk of prediabetes and type 2 diabetes combined by genetic risk scores according to black coffee consumption.
Table 5. Multivariate-adjusted ORs and 95% CIs for the risk of prediabetes and type 2 diabetes combined by genetic risk scores according to black coffee consumption.
Genetic Risk Scores 1Non-Black-Coffee Consumers 2Black-Coffee Consumersp for Interaction
<1 Cup/Day1 to <2 Cups/Day≥2 Cups/Day
0 point 0.261
Case/total1639/220447/6125/3946/61
ORs (95% CIs) 3Reference1.03 (0.54, 1.97)0.67 (0.33, 1.35)0.87 (0.46, 1.66)
0.1 to <5 points
Case/total524/72112/188/119/17
ORs (95% CIs) 3Reference1.00 (0.33, 3.06)0.90 (0.21, 3.94)0.49 (0.17, 1.44)
5 to 10 points
Case/total596/85414/2114/2014/27
ORs (95% CIs) 3Reference0.78 (0.28, 2.15)0.89 (0.31, 2.58)0.36 (0.15, 0.88)
Abbreviations: OR, odds ratio; 95% CI, 95% confidence interval. 1 Genetic risk scores were calculated by 5 SNPs related to coffee consumption weighted by relative effect size (β coefficient); 2 for black-coffee consumers, non-coffee consumers and sugared-coffee consumers combined were regarded as reference group; 3 the ORs (95% CIs) were adjusted for age (years, continuous), sex (for men and women combined), body mass index (BMI, <23, 23 to <25, 25 to <30 and ≥30 kg/m2), smoking status (never smokers, ≤10 and >10 pack-years), alcohol consumption (non-drinkers, ≤5, 5.1 to ≤10 and >10 g/day), family history of type 2 diabetes (yes or no), total energy intake (kcal/day, continuous), and the amount of sugar added in coffee (0, ≤5, 5.1 to ≤10, 10.1 to ≤15 and >15 g/day).
Table
Table 6. Multivariate-adjusted ORs and 95% CIs for the risk of prediabetes and type 2 diabetes combined by genetic risk scores according to sugared coffee consumption.
Table 6. Multivariate-adjusted ORs and 95% CIs for the risk of prediabetes and type 2 diabetes combined by genetic risk scores according to sugared coffee consumption.
Genetic Risk Scores 1Non-Sugared-Coffee Consumers 2Sugared-Coffee Consumersp for Interaction
<1 Cup/Day1 to <2 Cups/Day≥2 Cups/Day
0 point 0.608
Case/total521/698397/544442/571397/552
ORs (95% CIs) 3Reference0.79 (0.61, 1.04)1.11 (0.85, 1.46)0.74 (0.56, 0.97)
0.1 to <5 points
Case/total139/196121/162151/196142/213
ORs (95% CIs) 3Reference1.02 (0.62, 1.67)1.28 (0.79, 2.07)0.71 (0.45, 1.12)
5 to 10 points
Case/total174/245126/186156/219182/272
ORs (95% CIs) 3Reference0.79 (0.51, 1.22)1.03 (0.68, 1.57)0.74 (0.49, 1.12)
Abbreviations: OR, odds ratio; 95% CI, 95% confidence interval. 1 Genetic risk scores were calculated by 5 SNPs related to coffee consumption weighted by relative effect size (β coefficient); 2 for sugared-coffee consumers, non-coffee consumers and black-coffee consumers combined were regarded as reference group; 3 the ORs (95% CIs) were adjusted for age (years, continuous), sex (for men and women combined), body mass index (BMI, <23, 23 to <25, 25 to <30 and ≥30 kg/m2), smoking status (never smokers, ≤10, 10.1 to ≤20, 20.1 to ≤30 and >30 pack-years), alcohol consumption (non-drinkers, ≤5, 5.1 to ≤10, 10.1 to ≤20 and >20 g/day), family history of type 2 diabetes (yes or no), and total energy intake (kcal/day, continuous).
Table
Table 7. Multivariate-adjusted ORs and 95% CIs for the risk of prediabetes and type 2 diabetes combined by 5 coffee-related SNPs according to black coffee consumption.
Table 7. Multivariate-adjusted ORs and 95% CIs for the risk of prediabetes and type 2 diabetes combined by 5 coffee-related SNPs according to black coffee consumption.
Non-Black-Coffee Consumers 1Black-Coffee Consumersp for Interaction
<1 Cup/Day1 to <2 Cups/Day≥2 Cups/Day
rs2074356 0.171
GG
Case/total1989/269454/7133/4852/73
ORs (95% CIs) 2Reference1.07 (0.59, 1.94)0.83 (0.43, 1.61)0.77 (0.44, 1.35)
GA+AA
Case/total770/108519/2914/2217/32
ORs (95% CIs) 2Reference0.72 (0.31, 1.71)0.65 (0.25, 1.70)0.37 (0.16, 0.84)
rs11066015 0.143
GG
Case/total1894/255754/6929/4351/70
ORs (95% CIs) 2Reference1.24 (0.67, 2.30)0.81 (0.41, 1.62)0.84 (0.47, 1.51)
GA+AA
Case/total865/122219/3118/2718/35
ORs (95% CIs) 2Reference0.56 (0.25, 1.27)0.68 (0.28, 1.65)0.35 (0.16, 0.75)
rs12229654 0.366
TT
Case/total2049/276555/7331/4752/72
ORs (95% CIs) 2Reference1.02 (0.57, 1.83)0.69 (0.36, 1.32)0.75 (0.42, 1.33)
TG+GG
Case/total710/101418/2716/2317/33
ORs (95% CIs) 2Reference0.82 (0.33, 2.01)0.94 (0.35, 2.52)0.42 (0.19, 0.93)
rs11065828 0.460
CC
Case/total1909/257551/7026/4252/71
ORs (95% CIs) 2Reference0.88 (0.49, 1.56)0.60 (0.31, 1.18)0.81 (0.45, 1.45)
CA+AA
Case/total850/120422/3021/2817/34
ORs (95% CIs) 2Reference1.17 (0.48, 2.87)1.19 (0.46, 3.06)0.36 (0.17, 0.80)
rs79105258 0.395
CC
Case/total1986/267252/7031/4851/70
ORs (95% CIs) 2Reference0.98 (0.55, 1.77)0.66 (0.35, 1.25)0.79 (0.44, 1.42)
CA+AA
Case/total773/110721/3016/2218/35
ORs (95% CIs) 2Reference0.91 (0.38, 2.17)1.04 (0.37, 2.92)0.40 (0.18, 0.87)
Abbreviations: OR, odds ratio; 95% CI, 95% confidence interval. 1 For black-coffee consumers, non-coffee consumers and sugared-coffee consumers combined were regarded as reference group; 2 the ORs (95% CIs) were adjusted for age (years, continuous), sex (for men and women combined), body mass index (BMI, <23, 23 to <25, 25 to <30 and ≥30 kg/m2), smoking status (never smokers, ≤10 and >10 pack-years), alcohol consumption (non-drinkers, ≤5, 5.1 to ≤10 and >10 g/day), family history of type 2 diabetes (yes or no), total energy intake (kcal/day, continuous), and the amount of sugar added in coffee (0, ≤5, 5.1 to ≤10, 10.1 to ≤15 and >15 g/day.

Share and Cite

MDPI and ACS Style

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. https://doi.org/10.3390/nu12082228

AMA Style

Jin T, Youn J, Kim AN, Kang M, Kim K, Sung J, Lee JE. Interactions of Habitual Coffee Consumption by Genetic Polymorphisms with the Risk of Prediabetes and Type 2 Diabetes Combined. Nutrients. 2020; 12(8):2228. https://doi.org/10.3390/nu12082228

Chicago/Turabian Style

Jin, Taiyue, Jiyoung Youn, An Na Kim, Moonil Kang, Kyunga Kim, Joohon Sung, and Jung Eun Lee. 2020. "Interactions of Habitual Coffee Consumption by Genetic Polymorphisms with the Risk of Prediabetes and Type 2 Diabetes Combined" Nutrients 12, no. 8: 2228. https://doi.org/10.3390/nu12082228

APA Style

Jin, T., Youn, J., Kim, A. N., Kang, M., Kim, K., Sung, J., & Lee, J. E. (2020). Interactions of Habitual Coffee Consumption by Genetic Polymorphisms with the Risk of Prediabetes and Type 2 Diabetes Combined. Nutrients, 12(8), 2228. https://doi.org/10.3390/nu12082228

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