Coffee is among the most consumed non-alcoholic beverages worldwide( Reference Grigg 1 ). Besides the gastronomic aspects, coffee has been suggested to affect human health and disease( Reference Cano-Marquina, Tarin and Cano 2 ). Although in the past it has been considered unhealthful due to its association with some negative outcomes following acute ingestion (i.e. rise in blood pressure and plasma lipids, inhibition of insulin activity), adverse health outcomes of long-term coffee consumption have not been found. Indeed, several epidemiological studies have explored the relationship between coffee intake and cardiovascular risk factors( Reference Jiang, Zhang and Jiang 3 , Reference Steffen, Kuhle and Hensrud 4 ), CVD( Reference Ding, Bhupathiraju and Satija 5 ) and certain types of cancer( Reference Bohn, Blomhoff and Paur 6 ), suggesting that a protective effect may exist. From an epidemiological point of view, long-term coffee drinking exerts different effects from acute intake typically examined under experimental conditions. Long-term experimental studies have demonstrated that habitual coffee drinkers develop tolerance to the acute effects of caffeine and may experience the effects of other coffee compounds that are related with its beneficial role on human health( Reference Corti, Binggeli and Sudano 7 ). Coffee is a complex mixture of several biochemically active components, including caffeine, phenolic compounds and minerals (magnesium and potassium). Many of these compounds may act as antioxidants, thus preventing the state of subclinical inflammation involved in the pathogenesis of diseases such as cancer and CVD( Reference Godos, Pluchinotta and Marventano 8 ). Among them, chlorogenic acids are responsible for ameliorating hypertension, glucose and lipid metabolism, whereas diterpenes demonstrate anti-carcinogenic properties( Reference Ludwig, Clifford and Lean 9 ).
Recently published prospective studies evaluating the association between coffee consumption and mortality in large European cohorts have been mostly conducted in Northern( Reference Ahmed, Levitan and Wolk 10 – Reference Laaksonen, Talala and Martelin 12 ) and Western European countries( Reference Leurs, Schouten and Goldbohm 13 ). Findings obtained are puzzling and inconsistent, suggesting that significant differences between countries and methods of coffee consumption may be relevant for the outcome explored. We recently reported an inverse association between coffee consumption and metabolic impairment and decreased risk of developing hypertension in Polish individuals free of diseases at baseline( Reference Grosso, Stepaniak and Micek 14 , Reference Grosso, Stepaniak and Polak 15 ). However, whether these associations may be translated into a decreased risk of mortality following a prospective evaluation is unknown. To date, no further data from Eastern European countries are available. The aim of the present study was to test the association between coffee consumption and risk of all-cause, CVD and cancer mortality in the HAPIEE (Health, Alcohol and Psychosocial factors In Eastern Europe) study, a large-scale cohort study including adult individuals from the Czech Republic, Russia and Poland.
Methods
Study population
The HAPIEE study is a multicentre prospective cohort study investigating the role of biological, dietary, lifestyle and environmental factors in cardiovascular and other chronic diseases( Reference Peasey, Bobak and Kubinova 16 ). The criteria for sample selection and the methods are reported in detail elsewhere( Reference Peasey, Bobak and Kubinova 16 ). Briefly, a random sample of 28945 individuals (aged 45–69 years) was recruited at the baseline survey conducted in 2002–2005 (response ratio of 59 %) in Novosibirsk (Russia), Krakow (Poland) and six Czech urban centres. The participants, after providing written informed consent, filled in a structured questionnaire and were invited to a clinical examination. Participants who left more than 50 % of food items blank in the FFQ or reported implausible energy (energy intake <2092/>16736 kJ/d (<500/>4000 kcal/d) for females and <3347/>20920 kJ/d (<800/>5000 kcal/d) for males, n 384) were excluded from the analyses. A total of 28561 individuals were included for the final analysis. The study was conducted according to the guidelines laid down in the Declaration of Helsinki and all procedures involving human subjects were approved by the University College London Hospital research ethics committee and by ethical committees in each participating country. Written informed consent was obtained from all subjects.
Demographic, lifestyle and clinical measurements
Sociodemographic and lifestyle characteristics included age, gender, educational and occupational level, smoking and alcohol drinking habits. Physical activity included energy expenditure both at work and leisure time. The overall amount of energy expenditure was calculated and categorized as low, moderate and high activity level. Individuals were categorized according their smoking status as non-smoker and current smoker. Alcohol consumption was categorized as none or low (≤12 g/d) and moderate to high alcohol drinker (>12 g/d). The physical examination included measurement of height, weight and blood pressure using standard procedures( Reference Peasey, Bobak and Kubinova 16 ). BMI was calculated according the formula [weight (kg)]/[height (m)]2. Blood pressure was measured three times at the end of the physical examination and the final value was the mean among the three measurements. Participants were categorized as having hypertension at baseline if they had a systolic and/or a diastolic blood pressure measurement higher than 139 and/or 89 mmHg, respectively, or had taken hypertensive medication within the last 2 weeks. Participants were categorized as having diabetes at baseline if they had fasting plasma glucose >120 mg/dl, previously diagnosed type 2 diabetes or were receiving treatment for diabetes. Participants were categorized as having hypercholesterolaemia if they were previously diagnosed with or using a specific treatment for this lipid abnormality. Presence of CVD or cancer relied on previous diagnosis. Participants were finally asked whether they also had family history of CVD or cancer.
Dietary assessment
Dietary data were collected by using an FFQ based on the tool developed by Willett et al.( Reference Willett, Sampson and Stampfer 17 ) and subsequently adapted in the Whitehall II Study( Reference Brunner, Stallone and Juneja 18 ). The Czech, Russian and Polish FFQ consisted of 136, 147 and 148 food and drink items, respectively (the different numbers of questions are due to country-specific dishes). An instruction manual, which included photographs to facilitate the estimation of portion sizes, was used. Participants were asked how often, on average, they had consumed that amount of the item during the last 3 months, with nine responses ranging from ‘never or less than once per month’ to ‘six or more times per day’. Moreover, participants were asked to include additional foods and frequency of consumption by manual entry.
The average beverage consumption was calculated (in millilitres) by following the standard portion sizes used in the study and then converted into 24 h intake. We categorized daily coffee consumption according to standard cups of coffee (150 ml) as: (i) <1 cup/d; (ii) 1–2 cups/d; (iii) 3–4 cups/d; and (iv) >4 cups/d.
Mortality
The mortality status of participants was obtained through linkage with mortality registers of their respective country. Specifically, mortality data in Russia were obtained from the death register developed by the Institute of Internal Medicine, based on data from medical death certificates, the Novosibirsk office of the State Statistical Bureau and the population registration bureau. In Poland, data from the Local Register of Residents of the City of Krakow and Central Statistical Office were used. In the Czech Republic, data from the National Death Register were used. The follow-up period was until December 2011 for the Czech Republic or December 2010 for Russia and Poland. Participants who were lost to follow-up were censored at their last date of contact (n 527). The underlying cause of death was assigned according to the International Classification of Diseases, 10th Revision (ICD-10). The primary end point in the current analysis was death from any cause. We also conducted analyses according to CVD (ICD-10 codes I00–I99, CHD codes I20–I25, stroke codes I60–I69) and cancer mortality (codes C00–D48).
Statistical analysis
The distribution of background characteristics was compared between different coffee consumption groups using ANOVA or the Mantel–Haenszel test. We calculated person-years of follow-up for each participant from the date of the baseline questionnaire to the date of death from any cause, loss to follow-up or end of the follow-up period, whichever came first. Hazard ratios (HR) and 95 % CI for the association between coffee consumption and all end points, using the lowest category as the reference, were calculated by Cox proportional hazards models. No interactions between countries and categories of coffee consumption were detected, so we calculated the results of the Cox regression in the pooled sample and checked possible differences in an additional analysis by country cohorts. A crude model was adjusted for age (continuous) and gender. A multivariate model was built adjusting for country (Poland, Czech Republic, Russia), educational level (low, medium, high), physical activity (low, medium, high), smoking status (yes/no), history of CVD or cancer (yes/no), health status (presence of hypertension, diabetes, use of cholesterol-lowering medications), family history of CVD (yes/no), family history of cancer (yes/no), menopausal status in women (yes/no), alcohol intake >12 g/d (yes/no), total energy intake (continuous), SFA intake (continuous), PUFA intake (continuous), n-3 PUFA intake (continuous) and use of vitamin supplements (yes/no). The analysis was performed also for CVD and cancer mortality and stratified by gender and smoking status. To evaluate reverse causation, the risk of mortality excluding deaths occurring in the first 2 years from baseline, in individuals with a past history of CVD or cancer, and after sequential exclusion of individuals with CVD risk factors (hypertension, diabetes, use of cholesterol-lowering medications), was also estimated. All reported P values were based on two-sided tests and compared with a significance level of 5 %. Analyses were performed using the statistical software package SPSS Statistics Version 17.0.
Results
A total of 2121 deaths (1435 men and 686 women) occurred during a mean follow-up of 6·1 years. Of these deaths, 913 (43·1 %) were due to CVD and 755 (35·7 %) were due to cancer mortality. Distribution of background characteristics by categories of coffee consumption did not differ significantly between men and women (Table 1). No specific trends were identified, with the exception of the distribution of current smokers (increasing with increased category of coffee consumption) and prevalence of hypertension and CVD (decreasing with increased category of coffee consumption). High alcohol intake was more frequent among men than women (about 18–20 % v. 1–4 %, respectively) and associated with higher coffee consumption, whereas prevalence of hypercholesterolaemia and family history of CVD and cancer were higher in women than in men.
Baseline characteristics by coffee consumption in 28 561 men and women living in the Czech Republic, Poland and Russia in 2002–2005; the HAPIEE (Health, Alcohol and Psychosocial factors In Eastern Europe) study
Frequency and risk of all-cause, CVD and cancer mortality based on coffee consumption are shown in Table 2. Risk of all-cause mortality decreased with increasing coffee consumption in both men and women, both in age-adjusted and multivariate-adjusted models, with no significant difference in the distribution of deaths across categories of exposure between genders (P=0·094; Table 2). A significant risk reduction of 17 % (HR=0·83; 95 % CI 0·71, 0·99) in men and 37 % (HR=0·63; 95 % CI 0·47, 0·84) in women was found for consumption of 3–4 cups/d, whereas further intake led to a non-significant decrease in risk. The risk estimate of all-cause mortality by category of coffee consumption partially reflected the risk of mortality by CVD in both women and men, whereas non-significant decreased risk in women and non-significant increased in men for risk cancer mortality was found. The stratified analysis by country revealed consistent non-significant decreasing mortality risk estimates with increasing coffee consumption across all cohorts and significant decreased risk up to 4 cups/d or less (see online supplementary material, Supplemental Table 1). Moreover, the analysis revealed that the Russian cohort was the only one showing an increased non-significant risk of death by cancer, whereas the others showed a decreased, non-significant risk.
Age- and multivariate-adjusted hazard ratios (HR) and 95 % confidence intervals for all-cause, CVD and cancer mortality by coffee consumption in men and women living in the Czech Republic, Poland and Russia over a period of 6 years; the HAPIEE (Health, Alcohol and Psychosocial factors In Eastern Europe) study
Ref., reference category.
* Adjusted for age (continuous), gender, smoking status (yes/no), BMI (continuous), educational level (low, medium, high), physical activity (low, medium, high), alcohol intake (>12 g/d), hypertension, diabetes, hypercholesterolaemia, history of CVD or cancer, family history of CVD, family history of cancer, total energy intake (continuous), vitamin supplement use, SFA, PUFA and n-3 PUFA intake (continuous) and menopausal status (in women).
A number of analyses to evaluate possible reverse causation are reported in the online supplementary material, Supplemental Table 2. Analysis performed by excluding deaths occurring within the first 2 years from baseline of follow-up showed no significant alteration of the results among both men and women. No substantially different risks were estimated when excluding individuals with a past history of CVD or cancer. Sequential exclusion also of individuals with any CVD risk factors weakened CVD mortality risk in women due to the limited cases in the highest categories of exposure, whereas it led to a non-significant decreased risk of death from cancer in men consumers of 3–4 cups coffee/d.
Results of stratified analyses by smoking status and alcohol intake are presented in Tables 3 and 4, respectively. Distribution of deaths across categories of exposure between smokers and non-smokers was significantly different, with more numerous mortality cases among the categories of higher coffee intake in smokers than in non-smokers (P<0·001). Nevertheless, consumption of 3–4 cups coffee/d was associated with significantly decreased risk of CVD mortality in smoker men (HR=0·66; 95 % CI 0·47, 0·93) and decreased risk of all-cause (HR=0·61; 95 % CI 0·43, 0·87) and cancer mortality (HR=0·59; 95 % CI 0·35, 0·99) in non-smoker women (Table 3). Results were mainly confirmed after stratification by alcohol consumption, with no significant difference between groups (P=0·068; Table 4). However, among men with alcohol intake of ≤12 g/d, consumption of >4 cups coffee/d was associated with a decreased risk of all-cause mortality (HR=0·61; 95 % CI 0·39, 0·94).
Stratified analysis by smoking status of age- and multivariate-adjusted hazard ratios (HR) and 95 % confidence intervals for all-cause, CVD and cancer mortality by coffee consumption in men and women living in the Czech Republic, Poland, and Russia over a period of 6 years; the HAPIEE (Health, Alcohol and Psychosocial factors In Eastern Europe) study
Ref., reference category.
* Adjusted for age (continuous), gender, smoking status (yes/no), BMI (continuous), educational level (low, medium, high), physical activity (low, medium, high), alcohol intake (>12 g/d), hypertension, diabetes, hypercholesterolaemia, history of CVD or cancer, family history of CVD, family history of cancer, total energy intake (continuous), vitamin supplement use, SFA, PUFA and n-3 PUFA intake (continuous), and menopausal status (in women).
Stratified analysis by alcohol intake of multivariate-adjusted hazard ratios (HR) and 95 % confidence intervals for all-cause, CVD and cancer mortality by coffee consumption in men and women living in the Czech Republic, Poland, and Russia over a period of 6 years; the HAPIEE (Health, Alcohol and Psychosocial factors In Eastern Europe) study
Ref., reference category; NA, not applicable.
* Adjusted for age (continuous), gender, smoking status (yes/no), BMI (continuous), educational level (low, medium, high), physical activity (low, medium, high), smoking status (yes/no), hypertension, diabetes, hypercholesterolaemia, history of CVD or cancer, family history of CVD, family history of cancer, total energy intake (continuous), vitamin supplement use, SFA, PUFA and n-3 PUFA intake (continuous) and menopausal status (in women).
Discussion
In the present study we found that daily coffee consumption was significantly correlated with decreased risk of all-cause and CVD death, but the effect was not observed when drinking more than 4 cups/d. Interestingly, we also found that after stratifying the sample by smoking status, coffee consumption was inversely associated with risk of death by cancer in women drinking 3–4 cups coffee/d. Lack of decreased risk of death for higher intakes of coffee has been recently reported, suggesting that the association of coffee drinking with mortality risk may not be linear( Reference Lof, Sandin and Yin 19 ). These findings are in line with a recent dose–response meta-analysis( Reference Crippa, Discacciati and Larsson 20 ). It has been hypothesized that intake of high quantities of coffee may have an acute effect at the vascular level by triggering coronary events and, consequently, increasing the risk of myocardial infarction in specific patient groups, such as the elderly and those at high CVD risk( Reference Greenberg, Dunbar and Schnoll 21 , Reference Rosner, Akesson and Stampfer 22 ). This may depend on the biological effect of caffeine intake that causes an acute increase in arterial wave reflection, which in turn increases the pulsatile load of the heart. Furthermore, it has been reported that excessive intakes may also aggravate cardiac arrhythmias and raise plasma homocysteine( Reference Verhoef, Pasman and Van Vliet 23 ). However, other studies have reported that habitual coffee consumption is not directly associated with myocardial infarction in the general population and that the effects of caffeine in habitual coffee drinkers are attenuated or overall absent due to the development of tolerance( Reference Corti, Binggeli and Sudano 7 , Reference Ammon, Bieck and Mandalaz 24 – Reference Robertson, Wade and Workman 27 ). The contrasting hypotheses on the effects of caffeine on cardiovascular events may reflect the modality of coffee consumption. Indeed, a habitual, rather than occasional intake of coffee may lead to desensitization to the acute effects of caffeine but excessive intake (i.e. more than 4 cups/d) could exacerbate its harmful effects.
The analysis of cancer mortality according coffee consumption revealed non-significant associations towards both increased and decreased risk, but after stratifying the analysis by smoking status, consumption of 3–4 cups coffee/d was associated with a decreased risk of death from cancer of about 40 % compared with non-drinkers. Some other cohort studies specifically explored the association between coffee consumption and cancer mortality, reporting a decreased mortality risk from cancer( Reference Tamakoshi, Lin and Kawado 28 ). Moreover, recent meta-analyses including epidemiological studies showed a protective effect of coffee consumption against liver, colorectal, breast and prostate cancer incidence( Reference Bohn, Blomhoff and Paur 6 ). Coffee has been hypothesized to play a role in the early phases of cancer development and promotion, delaying disease progression and thus prolonging a person’s lifespan. As we found a significant association between coffee consumption and cancer mortality only in non-smoker women, we may hypothesize that smoking status might be a strong confounding factor and that the potentially beneficial effects of coffee are appreciable only when analysis is restricted to individuals free of the confounding effect of smoking. Another hypothesis is that the harmful effects of smoking among smokers counteract the potentially protective effect of coffee. However, some questions remain, for instance the lack of an evident decreased risk among non-smoker men. Knowledge of the type of cancer could help to draw more accurate hypotheses and explanations, but the data presented in the current study lack this information and future studies should take into account our results to better design their analyses.
The polyphenol content of coffee represents the main biological explanation for its benefits on health. During recent years, dietary patterns rich in polyphenols, such as the Mediterranean diet, have been suggested to exert their benefits due to their polyphenol content that may contribute towards their protective effects against diseases with an inflammatory basis( Reference Grosso, Marventano and Yang 29 ). Coffee has lately been demonstrated to be among the main contributors to polyphenol content in the diet of individuals living in Eastern European countries( Reference Grosso, Stepaniak and Topor-Madry 30 – Reference Zujko, Witkowska and Waskiewicz 34 ) as well as among the foods with the highest antioxidant capacity( Reference Pellegrini, Serafini and Colombi 35 ). Among the most important, chlorogenic acids are molecules which belong to the group of hydroxycinnamic acids demonstrating the ability to affect some metabolic pathways, for instance by improving glucose metabolism and decreasing inflammation and endothelial dysfunction( Reference Lopez-Garcia, van Dam and Qi 36 ). This family of molecules has also demonstrated anti-tumour properties by inhibiting key enzymes involved in tumour genesis and metastasis( Reference Jin, Lee and Kang 37 ). Other components, such as the diterpenes kahweol and cafestol, may play a role in the regulation of angiogenesis, apoptosis and inflammation processes and have been proposed as possible anti-tumour agents( Reference Cardenas, Quesada and Medina 38 , Reference Wang, Yoon and Sung 39 ). Together with these effects, an excess of the same compounds, for instance caffeine and diterpenes, may produce harmful effects, such as an acute rise in blood pressure, increase in serum cholesterol and homocysteine levels, and decreased insulin sensitivity, all contributing to raise the risk for CVD( Reference Verhoef, Pasman and Van Vliet 23 , Reference Hartley, Lovallo and Whitsett 40 , Reference Thong and Graham 41 ). The equilibrium between positive and negative effects on health may explain the findings of our and other studies, as an excessive consumption was not associated with the same benefits as moderate daily intake. The antioxidant compounds are thought to counteract the potentially harmful effects of caffeine or diterpenes, but an excessive intake of the latter may disrupt such balance, thus weakening or even annulling the overall beneficial effect of coffee. Moreover, other factors may impact the effects of coffee consumption on health. For instance, some previous studies have emphasized the possible role of associated lifestyle habits, such as smoking status, which largely correlates with mortality. Also in the present study we confirmed that current smoking, associated with coffee drinking, was weakening the effect of coffee on CVD mortality, as a stronger association between coffee and CVD mortality was found among non-smokers. In light also of the aforementioned findings obtained on cancer mortality, an accurate study of smoking status and coffee drinking should be conducted in future investigations aimed to explain the effects of coffee on human health, whatever is the outcome.
It is still debatable if the health status of a person may influence coffee consumption (i.e. those with poorer health conditions drinking less coffee), thus biasing the results. Generally, studies on coffee consumption and mortality report analyses with the exclusion of deaths occurring within the first 2 years from baseline observation with no substantial changing of the results. In our study we performed a similar analysis and the results suggested that health status at baseline did not necessarily explain the aforementioned association. In addition, we also took into account health status by adjusting for a previous history of CVD and cancer, as well as the presence of CVD risk factors, in the final model of multivariable analysis and excluded these cases in an additional sub-analysis in order to reveal possible reverse causation. Besides these strengths, the current study also presented some limitations. First, the data presented refer only to the baseline examination and coffee consumption was self-reported. Therefore, we were not able to control for possible variation of coffee intake during the follow-up period. Second, data regarding the type of coffee and methods of preparation were not collected. Thus, we can only assume homogeneous preparation methods within and between countries, but differences regarding the use of caffeinated and decaffeinated coffee may still exist. Third, the study had relatively short follow-up, thus cases are limited, especially regarding cancer mortality. Fourth, the last category of exposure included fewer individuals than others and this could have affected the statistical power. Nevertheless, for analysis, we preferred to consider this category as separate due to a different biological effect of extreme coffee consumption. Finally, the sample was recruited in urban centres and thus should not be considered nationally representative.
Conclusion
In conclusion, in our cohort coffee consumption was associated with decreased risk of all-cause and CVD mortality. The healthful effect of coffee seems to rely on CVD, but additional analyses focused on non-smokers are needed to further investigate potentially preventive effects also against cancer.
Acknowledgements
Financial support: The study was funded by the Wellcome Trust (grant numbers 064947/Z/01/Z and 081081/Z/06/Z); the US National Institute on Aging (grant number 1R01 AG23522-01); and the MacArthur Foundation Initiative on Social Upheaval and Health (award number 71208). The funders had no role in the design, analysis or writing of this article. Conflict of interest: The authors declare that they have no conflict of interest. Authorship: G.G. designed the study, performed the analysis and drafted the manuscript. A.M. performed the analysis. U.S. contributed in collecting data. D.S. provided critical revision for the data analysis methods. A.P. and M.B. designed the study and provided critical revision. All authors provided final approval for the manuscript. Ethics of human subject participation: This study was conducted according to the guidelines laid down in the Declaration of Helsinki and all procedures involving human subjects were approved by the University College London Hospital research ethics committee and by ethical committees in each participating country. Written informed consent was obtained from all subjects.
Supplementary material
To view supplementary material for this article, please visit http://dx.doi.org/10.1017/S1368980016001749
