- Open Access
Association between hypertension and coffee drinking based on CYP1A2 rs762551 single nucleotide polymorphism in Taiwanese
Nutrition & Metabolism volume 18, Article number: 78 (2021)
Hypertension increases the likelihood of cardiovascular diseases (CVDs). Cytochrome P450 1A2 (CYP1A2) single nucleotide polymorphism (SNP) is related to caffeine metabolism and the risk of CVD among coffee drinkers. CYP1A2 rs762551 influenced the risk of stroke among hypertensive patients. We examined the relationship between hypertension and coffee drinking based on CYP1A2 rs762551 SNP in Taiwanese adults.
We used data contained in the Taiwan Biobank database (2011–2018) and included 19,133 participants having complete information on hypertension, rs762551 polymorphism, coffee intake, etc. The risk of hypertension was determined using multiple logistic regression.
Coffee intake was significantly associated with a lower risk of hypertension. The odds ratio (OR), 95% confidence interval (CI), and p-value were 0.877, 0.807–0.954, and 0.0032, respectively. CYP1A2 rs762551 was not significantly associated with the risk of hypertension, but it had a significant interactive association with coffee drinking (p value = 0.0303). After stratification by rs762551 genotypes, the inverse coffee drinking-hypertension association was retained, but significant results were observed only in those with the AC + CC genotype (OR 0.678, 95% CI 0.722–900, p value = 0.0001). According to the combination of coffee drinking and rs762551 genotypes (reference group: no coffee drinking and rs762551 AA), the coffee drinking-AC + CC group had a lower risk of hypertension (OR 0.888, 95% CI 0.789–0.999, p value = 0.0483).
Coffee drinking, particularly among individuals with the CYP1A2 rs762551 AC + CC genotype was associated with lower odds of hypertension.
Hypertension is a significant public health burden. Its prevalence among adults has been increasing. In 2000, about 972 million adults in the world were believed to be hypertensive . However, by 2010, hypertension among adults was estimated at 1.4 billion . It is projected at about 1.56 billion by 2025 . The disease is multiplying rapidly in Asia, and the prevalence rate in Taiwan is about 25% in men and 18% in women . Hypertension increases the likelihood of CVDs like heart failure, stroke, and coronary heart disease [1, 4,5,6], which account for a considerable number of global deaths [7, 8]. The burden of deaths attributable to CVDs was approximately 12.3 and 17.3 million in 1990 and 2013, respectively , and 17.92 million in 2015 . Some cardiovascular outcomes are believed to be the consequence of the effect of coffee on blood pressure (BP) .
Coffee consumption has been increasing in Taiwan and the world [10, 11]. Dietary and lifestyle factors play substantial roles in modifying the risk of hypertension . A previous study indicated that regular consumption of three or more cups of coffee a day was protective against hypertension . In the SUN project, regular consumption of coffee was associated with a decreased risk of hypertension (about 26%) only in women . A meta-analysis suggested that long-term consumption of coffee may lower the incidence of hypertension . Daily intake of green or black tea was associated with a modest increase in hypertension risk . However, the risk became non-significant when caffeine was adjusted for in the model. Marjan and colleagues found that regular intake of tea reduced systolic blood pressure (SBP) and diastolic blood pressure (DBP) by 4.81 mmHg and 1.98 mmHg, respectively .
So far, findings from epidemiological research on coffee drinking and hypertension have been controversial . Some short-term intervention studies revealed a positive association . However, one prospective study found that coffee intake was beneficial for hypertension only among non-smokers  while several prospective studies found that coffee intake could elevate the risk of hypertension [20,21,22,23].
Genetic factors contribute to variability in the effect of coffee on cardiovascular health . Specific genetic polymorphisms related to caffeine metabolism pose a risk of CVD among coffee drinkers. For instance, polymorphism at CYP1A2 rs762551 was associated with impaired caffeine metabolism, owing to low enzyme actions . This variant was also associated with stroke risk in a Chinese population with hypertension, cerebral infarction, and coronary heart disease . About 95% of caffeine metabolism is attributed to CYP1A2 . CYP1A2 modulated the effect of caffeinated coffee on myocardial infarction  and its polymorphism was associated with hypertension . Moreover, caffeine intake modified the association between CYP1A2 variants and blood pressure (BP) among European adults [30, 31]. Research on hypertension based on CYP1A2 polymorphism and coffee intake in Taiwanese is seemingly rare. Considering the controversial influence of coffee on hypertension and cardiovascular health, and the importance of genetic insights in the management of diseases in this era, we undertook this study to evaluate the relationship between hypertension and coffee drinking based on CYP1A2 rs762551 single nucleotide polymorphism in Taiwanese.
We conducted this study using data contained in the Taiwan Biobank database (2011–2018). The biobank contains genetic data as well as information like medical history, sex, age, diet, and other lifestyle factors. Enrolment into the Taiwan Biobank project is restricted to Taiwanese adults between 30 and 70 years without a personal history of cancer.
Study population and variable definitions
Variables included hypertension (yes/no), coffee drinking (yes/no), rs762551 genotypes (AA/AC + CC), tea drinking (yes/no), sex (men/women), age, body mass index or simply BMI (underweight/normal /overweight/obese), cigarette smoking (never/former/current), alcohol drinking (never/former/current), and exercise (no/yes). Hypertensive cases were those who confirmed through questionnaires to have ever been clinically diagnosed with the disease. Coffee drinking was defined according to self-reports on weekly consumption frequency. The coffee-drinking group included individuals who self-reported a habit of drinking coffee at least thrice per week while the non-drinking group included those who habitually drank coffee less than three times per week. Tea drinking was defined based on daily consumption frequency. The tea-drinking category included those who self-reported a habitual daily tea intake of at least one time while the non-drinking category included those who had a habit of not drinking tea every day. Sex and age were self-reported while BMI was derived from weight (kg) and height (m) using the standard formula: weight divided by (height)2. Underweight was defined as BMI < 18.5 kg/m2, normal BMI as 18.5 ≤ BMI < 24 kg/m2, overweight as 24 ≤ BMI < 27 kg/m2, and obesity as BMI ≥ 27 kg/m2. Participants were considered non-smokers if they had no history of cigarette smoking for continuously 6 months or more. Former smokers were those who smoked cigarettes continuously for at least 6 months but were not smoking during data collection. Current smokers were those who had a history of continuous cigarette smoking for at least 6 months and were still smoking during data collection. Participants were considered non-drinkers if they had no history of continuous drinking for 6 months or drank less than 150 cc of alcohol per week. Former drinkers were those who abstained from alcohol for more than 6 months while current drinkers were those who drank at least 150 cc of alcohol per week continuously for 6 months. Exercise was defined as regular exercise taken 3 times a week (with each session lasting over 30 min). This study included 3411 hypertensive and 15,722 non-hypertensive individuals. All participants signed an informed consent letter. We obtained ethical approval from the Institutional Review Board of Chung Shan Medical University (CS1-20,009).
CYP1A2 rs762551 was validated using quality control measures including, Hardy–Weinberg equilibrium p > 0.001, a minor allele frequency ≥ 0.05, and call rate ≥ 0.05. Quality control was performed using PLINK version 1.9. Chi-square test was used to differentiate the categorical variables while Student’s T test was used to differentiate the continuous variables between cases and non-cases. Categorical variables were presented as n (percentage) and continuous variables as mean ± standard error. Multivariate logistic regression was used to evaluate the association of coffee drinking and rs762551 with hypertension and also the interaction between coffee drinking and rs762551. Statistical analyses were performed using SAS version 9.4 (SAS Institute Inc., Cary, North Carolina, USA), and p values less than 0.05 were considered to be statistically significant.
Table 1 shows the demographic data of 19,133 participants stratified as coffee drinkers (n = 7766) and non-drinkers (n = 11,367). Of the 19,133 participants, 3411 were hypertensive while 15,722 were non-hypertensive. Coffee drinkers and non-drinkers significantly differed according to hypertension (p value < 0.0001) but not according to rs762551 genotypes (p value = 0.8659). Coffee drinkers and non-drinkers also differed in terms of tea drinking (p value < 0.0001), sex (p value = 0.0154), age, BMI, cigarette smoking, and alcohol drinking (p value < 0.0001) but not exercise (p value = 0.2581).
Table 2 shows the association of hypertension with coffee and other variables. Coffee intake was significantly associated with a lower (OR 0.877, 95% CI 0.807–0.954, p value = 0.0032) while the rs762551 variant was associated with a higher but not significant risk of hypertension (OR 1.021, 95% CI 0.941–1.107, p value = 0.6176).
The rs762551 variant had a significant interaction with coffee (p value = 0.0303) but not with tea (p value = 0.4183) as seen in Table 3. After stratification by rs762551 genotypes, coffee drinking was significantly associated with a lower risk of hypertension only in those with the AC + CC genotype (OR 0.806, 95% CI 0.722–0.900, p value = 0.0001). Further stratification of coffee drinking revealed a significant dose–response association (p value < 0.0001) between coffee drinking and hypertension in those carrying the AC + CC genotype (Table 4).
According to the combination of coffee drinking and the rs762551 genotypes (Table 5) with the reference group being no coffee drinking and the rs762551 AA genotype, the risk of hypertension was lower in the coffee drinking and AC + CC group (OR 0.888, CI 0.789–0.999, p value = 0.0483).
We used the 2011–2018 Taiwan Biobank cohort in the current study. Descriptive data showed that coffee drinkers and non-drinkers significantly differed according to hypertension but not according to the rs762551 genotypes. Logistic regression analysis also revealed a non-significant association between AC + CC genotype (reference: AA) and coffee drinking (data not shown). After multivariate logistic regression analysis, coffee drinking had a significant inverse association with hypertension. CYP1A2 rs762551 was not significantly associated with the risk of hypertension but had a significant interaction with coffee drinking. After stratifying participants by rs762551 genotypes, the inverse coffee-hypertension association was retained only among those with the AC + CC genotype. Combining coffee and rs762551 genotypes revealed a lower risk of hypertension among coffee drinkers with the AC + CC genotype. The inverse association between coffee drinking and hypertension in the AC + CC group was in a dose-dependent manner.
The non-significant association between coffee drinking and rs762551 genotype observed in the current study is in concordance with previous reports among Costa Ricans  and Italians . Contrarily, coffee drinking was significantly lower among Japanese with the AC + CC genotype . In a meta-analysis, coffee drinking among Asians was not significantly associated with the rs762551 genotypes . On the other hand, the aforementioned meta-analysis revealed that coffee drinking was significantly higher among Caucasians carrying the rs762551 AA genotype compared to the AC + CC genotype .
The relationship between coffee and hypertension has been controversial. For instance, a meta-analysis of randomized clinical trials demonstrated a positive but weak association between coffee intake and hypertension . In a follow-up study, coffee intake was associated with higher blood pressure but played a very minute role in the onset of hypertension . In another follow-up study, there was no significant linear association between coffee drinking and hypertension . While a lower risk of hypertension was observed in men who abstained from coffee compared to those who drank moderate amounts, heavy coffee consumption among women was associated with a lower risk of hypertension . In conformity to our observation, results from a Mendelian randomization analysis revealed an inverse association between coffee drinking and hypertension . The seemingly protective effect of coffee against hypertension could in part elucidate why coffee drinking is associated with a lower risk of stroke [37, 38]. Coffee drinking is also associated with other cardiovascular health benefits, including protection against cerebral infarction , coronary heart disease , and CVD-related/all-cause mortality . These beneficial effects of coffee are attributed to its richness in polyphenols which possess antihypertensive, antioxidant, and anti-inflammatory properties [16, 41].
The observed differences in the relationship between coffee drinking and hypertension are in part due to genetics . Coffee consumption was associated with an increase in CYP1A2 enzyme activity . Moreover, the association of hypertension with coffee drinking was also confirmed to vary according to CYP1A2 genotypes . AA homozygotes compared to AC individuals have greater enzyme activities and are believed to metabolize caffeine relatively fast [25, 44]. Both moderate and heavy coffee consumption among Italians carrying the C allele appeared to confer higher susceptibility to hypertension . Swedish individuals with higher caffeine intake carrying the A allele were more likely to have a lower BP and lower risk of hypertension than those carrying the C allele . In the present study, coffee drinking among Taiwanese adults aged 30–70 years carrying the AC + CC genotypes was associated with a lower risk of hypertension. The discrepancy in findings could be due to differences in sample sizes, age, and ethnicity. In the study of Palatino and colleagues (2009), only 553 individuals were included in the study. Moreover, the age range of the participants was 18–45 years . We had a relatively large sample size (n = 19,133) and the age-range was 30–70 years. It was proposed that the CYP1A2 enzyme affects BP by differentially inhibiting adenosine receptors (42).
Despite the large sample size, our study is limited in that the TWB questionnaires contained no information on caffeine intake and coffee type. Therefore, we could not determine the interaction between CYP1A2 polymorphism and caffeine or decaffeinated coffee on hypertension risk. Moreover, we defined hypertension based on the Taiwan Biobank questionnaires (as described in the Methods section). This might have resulted in information bias. However, it should be noted that in our initial analyses, we merged the Taiwan Biobank database with the National Health Insurance Research Database (NHIRD) using encrypted personal identification numbers with the resources at the Health and Welfare Data Science Center (HWDC). The NHIRD contained diagnostic codes based on the International Classification of Diseases, Ninth Revision, Clinical Modification (ICD-9-CM). Hypertension was identified based on either two outpatient visits or one-time hospitalization between 1998 and 2015 using the diagnostic codes ICD-9-CM 401-405, A260, and A269. Even though we defined hypertension based on self-reports, all cases were those whose diagnosis was confirmed by a medical doctor. Noteworthy, the results observed are similar to those observed when hypertension was defined using ICD codes.
Generally, coffee drinking was associated with a lower risk of hypertension, while the AC + CC genotype was associated with a higher but non-significant risk. Notwithstanding, coffee drinking and rs762551 had a significant interactive relationship. Grouped analyses revealed an inverse association between hypertension and coffee drinking, particularly among individuals with the CYP1A2 rs762551 AC + CC genotype. Our results indicate the possibility of a protective effect of coffee on hypertension and confirm the joint role of coffee and CYP1A2 rs762551 polymorphism on hypertension.
Body mass index
Cytochrome P450 1A2
Single nucleotide polymorphism
National Health Insurance Research Database
International Classification of Diseases, Ninth Revision, Clinical Modification
Health and Welfare Data Science Center
Kearney PM, Whelton M, Reynolds K, Muntner P, Whelton PK, He J. Global burden of hypertension: analysis of worldwide data. Lancet. 2005;365(9455):217–23.
Mills KT, Bundy JD, Kelly TN, Reed JE, Kearney PM, Reynolds K, et al. Global disparities of hypertension prevalence and control: a systematic analysis of population-based studies from 90 countries. Circulation. 2016;134(6):441–50.
Chiang C-E, Wang T-D, Ueng K-C, Lin T-H, Yeh H-I, Chen C-Y, et al. 2015 Guidelines of the Taiwan Society of Cardiology and the Taiwan Hypertension Society for the Management of Hypertension. J Chin Med Assoc. 2015;78(1):1–47.
Geleijnse JM. Habitual coffee consumption and blood pressure: an epidemiological perspective. Vascular Health Risk Manag. 2008;4(5):963.
Stamler J, Stamler R, Neaton JD. Blood pressure, systolic and diastolic, and cardiovascular risks: US population data. Arch Intern Med. 1993;153(5):598–615.
Lawes CM, Vander Hoorn S, Rodgers A. Global burden of blood-pressure-related disease, 2001. Lancet. 2008;371(9623):1513–8.
Roth GA, Huffman MD, Moran AE, Feigin V, Mensah GA, Naghavi M, et al. Global and regional patterns in cardiovascular mortality from 1990 to 2013. Circulation. 2015;132(17):1667–78.
Roth GA, Johnson C, Abajobir A, Abd-Allah F, Abera SF, Abyu G, et al. Global, regional, and national burden of cardiovascular diseases for 10 causes, 1990 to 2015. J Am Coll Cardiol. 2017;70(1):1–25.
De Caterina R, El-Sohemy A. Moving towards specific nutrigenetic recommendation algorithms: caffeine, genetic variation and cardiovascular risk. Lifestyle Genomics. 2016;9(2–4):106–15.
Loftfield E, Cornelis MC, Caporaso N, Yu K, Sinha R, Freedman N. Association of coffee drinking with mortality by genetic variation in caffeine metabolism: findings from the UK Biobank. JAMA Intern Med. 2018;178(8):1086–97.
Wann J-W, Kao C-Y, Yang Y-C. Consumer preferences of locally grown specialty crop: The Case of Taiwan Coffee. Sustainability. 2018;10(7):2396.
Slama M, Susic D, Frohlich ED. Prevention of hypertension. Curr Opin Cardiol. 2002;17(5):531–6.
D’Elia L, La Fata E, Galletti F, Scalfi L, Strazzullo P. Coffee consumption and risk of hypertension: a dose–response meta-analysis of prospective studies. Eur J Nutr. 2019;58(1):271–80.
Navarro AM, Martinez-Gonzalez MA, Gea A, Ramallal R, Ruiz-Canela M, Toledo E. Coffee consumption and risk of hypertension in the SUN Project. Clin Nutr. 2019;38(1):389–97.
Grosso G, Micek A, Godos J, Pajak A, Sciacca S, Bes-Rastrollo M, et al. Long-term coffee consumption is associated with decreased incidence of new-onset hypertension: a dose–response meta-analysis. Nutrients. 2017;9(8):890.
Larsson SC. Coffee, tea, and cocoa and risk of stroke. Stroke. 2014;45(1):309–14.
Mahdavi-Roshan M, Salari A, Ghorbani Z, Ashouri A. The effects of regular consumption of green or black tea beverage on blood pressure in those with elevated blood pressure or hypertension: a systematic review and meta-analysis. Complement Therap Med. 2020;51:102430.
Steffen M, Kuhle C, Hensrud D, Erwin PJ, Murad MH. The effect of coffee consumption on blood pressure and the development of hypertension: a systematic review and meta-analysis. J Hypertens. 2012;30(12):2245–54.
Miranda AM, Goulart AC, Benseñor IM, Lotufo PA, Marchioni DM. Coffee consumption and risk of hypertension: a prospective analysis in the cohort study. Clin Nutr. 2020;40:542.
Hamer M. Coffee and health: explaining conflicting results in hypertension. J Hum Hypertens. 2006;20(12):909–12.
Noordzij M, Uiterwaal CS, Arends LR, Kok FJ, Grobbee DE, Geleijnse JM. Blood pressure response to chronic intake of coffee and caffeine: a meta-analysis of randomized controlled trials. LWW. 2005;23:921.
Grosso G, Stepaniak U, Polak M, Micek A, Topor-Madry R, Stefler D, et al. Coffee consumption and risk of hypertension in the Polish arm of the HAPIEE cohort study. Eur J Clin Nutr. 2016;70(1):109–15.
Zhang Z, Hu G, Caballero B, Appel L, Chen L. Habitual coffee consumption and risk of hypertension: a systematic review and meta-analysis of prospective observational studies. Am J Clin Nutr. 2011;93(6):1212–9.
Yang A, Palmer AA, de Wit H. Genetics of caffeine consumption and responses to caffeine. Psychopharmacology. 2010;211(3):245–57.
Sachse C, Brockmöller J, Bauer S, Roots I. Functional significance of a C→ A polymorphism in intron 1 of the cytochrome P450 CYP1A2 gene tested with caffeine. Br J Clin Pharmacol. 1999;47(4):445–9.
Mao Y, Yang L, Chen Q, Li G, Sun Y, Wu J, et al. The influence of CYP1A1 and CYP1A2 polymorphisms on stroke risk in the Chinese population. Lipids Health Dis. 2020;19(1):1–9.
Rasmussen BB, Brix TH, Kyvik KO, Brøsen K. The interindividual differences in the 3-demthylation of caffeine alias CYP1A2 is determined by both genetic and environmental factors. Pharmacogenet Genomics. 2002;12(6):473–8.
Cornelis MC, El-Sohemy A, Kabagambe EK, Campos H. Coffee, CYP1A2 genotype, and risk of myocardial infarction. JAMA. 2006;295(10):1135–41.
Palatini P, Ceolotto G, Ragazzo F, Dorigatti F, Saladini F, Papparella I, et al. CYP1A2 genotype modifies the association between coffee intake and the risk of hypertension. J Hypertens. 2009;27(8):1594–601.
Guessous I, Dobrinas M, Kutalik Z, Pruijm M, Ehret G, Maillard M, et al. Caffeine intake and CYP1A2 variants associated with high caffeine intake protect non-smokers from hypertension. Hum Mol Genet. 2012;21(14):3283–92.
Soares RN, Schneider A, Valle SC, Schenkel PC. The influence of CYP1A2 genotype in the blood pressure response to caffeine ingestion is affected by physical activity status and caffeine consumption level. Vascul Pharmacol. 2018;106:67–73.
Kohno M, Tajima O, Uezono K, Tabata S, Abe H, Adachi M, et al. Cytochrome P450 1A2 polymorphisms, coffee consumption and impaired glucose metabolism in Japanese men. Endocrinol Metab Syndr. 2013;2:119.
Denden S, Bouden B, Haj Khelil A, Ben Chibani J, Hamdaoui M. Gender and ethnicity modify the association between the CYP1A2 rs762551 polymorphism and habitual coffee intake: evidence from a meta-analysis. Genet Mol Res. 2016;15(2).
Klag MJ, Wang N-Y, Meoni LA, Brancati FL, Cooper LA, Liang K-Y, et al. Coffee intake and risk of hypertension: the Johns Hopkins precursors study. Arch Intern Med. 2002;162(6):657–62.
Winkelmayer WC, Stampfer MJ, Willett WC, Curhan GC. Habitual caffeine intake and the risk of hypertension in women. JAMA. 2005;294(18):2330–5.
Uiterwaal CS, Verschuren WM, Bueno-de-Mesquita HB, Ocké M, Geleijnse JM, Boshuizen HC, et al. Coffee intake and incidence of hypertension. Am J Clin Nutr. 2007;85(3):718–23.
Lopez-Garcia E, Rodriguez-Artalejo F, Rexrode KM, Logroscino G, Hu FB, van Dam RM. Coffee consumption and risk of stroke in women. Circulation. 2009;119(8):1116.
Larsson SC, Virtamo J, Wolk A. Coffee consumption and risk of stroke in women. Stroke. 2011;42(4):908–12.
Larsson SC, Männistö S, Virtanen MJ, Kontto J, Albanes D, Virtamo J. Coffee and tea consumption and risk of stroke subtypes in male smokers. Stroke. 2008;39(6):1681–7.
Liu Y-T, Tantoh DM, Wang L, Nfor ON, Hsu S-Y, Ho C-C, et al. Interaction between coffee drinking and TRIB1 rs17321515 single nucleotide polymorphism on coronary heart disease in a Taiwanese population. Nutrients. 2020;12(5):1301.
Butt MS, Sultan MT. Coffee and its consumption: benefits and risks. Crit Rev Food Sci Nutr. 2011;51(4):363–73.
Guessous I, Eap CB, Bochud M. Blood pressure in relation to coffee and caffeine consumption. Curr Hypertens Rep. 2014;16(9):468.
Djordjevic N, Ghotbi R, Bertilsson L, Jankovic S, Aklillu E. Induction of CYP1A2 by heavy coffee consumption in Serbs and Swedes. Eur J Clin Pharmacol. 2008;64(4):381–5.
Castorena-Torres F, Mendoza-Cantú A, De León MB, Cisneros B, Zapata-Pérez O, López-Carrillo L, et al. CYP1A2 phenotype and genotype in a population from the Carboniferous Region of Coahuila. Mexico Toxicol Lett. 2005;156(3):331–9.
This work was funded by the University of Foresight Education Project, Ministry of Education, Taiwan-Dementia Health Care Professional Training Program (A-01-108-0-0741) and partly by the Ministry of Science and Technology (MOST), Taiwan (MOST 109-2121-M-040-002, 110-2121-M-040-002110-2121-M-040-002).
Ethics approval and consent to participate
All participants signed an informed consent letter. We obtained ethical approval from the Institutional Review Board of Chung Shan Medical University (CS1-20009).
Consent for publication
Availability of data and materials
The data that support the findings of this study are available from Taiwan Biobank but restrictions apply to the availability of these data, which were used under license for the current study, and so are not publicly available. Data are however available from Professor Yung Po Liaw (Email address: Liawyp@csmu.edu.tw, Tel: + 886424730022 ext. 12102; fax: + 886423248179) upon reasonable request and with permission of Taiwan Biobank.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Hou, CC., Tantoh, D.M., Lin, CC. et al. Association between hypertension and coffee drinking based on CYP1A2 rs762551 single nucleotide polymorphism in Taiwanese. Nutr Metab (Lond) 18, 78 (2021). https://doi.org/10.1186/s12986-021-00605-9
- Coffee drinking