Abstract

Background

Accumulating evidence suggests that consuming coffee may lower the risk of death, but evidence regarding tea consumption in Asians is limited. We examined the association between coffee and tea consumption and mortality in Asian populations.

Methods

We used data from 12 prospective cohort studies including 248 050 men and 280 454 women from the Asia Cohort Consortium conducted in China, Japan, Korea and Singapore. We estimated the study-specific association of coffee, green tea and black tea consumption with mortality using Cox proportional-hazards regression models and the pooled study-specific hazard ratios (HRs) using a random-effects model.

Results

In total, 94 744 deaths were identified during the follow-up, which ranged from an average of 6.5 to 22.7 years. Compared with coffee non-drinkers, men and women who drank at least five cups of coffee per day had a 24% [95% confidence interval (CI) 17%, 29%] and a 28% (95% CI 19%, 37%) lower risk of all-cause mortality, respectively. Similarly, we found inverse associations for coffee consumption with cardiovascular disease (CVD)-specific and cancer-specific mortality among both men and women. Green tea consumption was associated with lower risk of mortality from all causes, CVD and other causes but not from cancer. The association of drinking green tea with CVD-specific mortality was particularly strong, with HRs (95% CIs) of 0.79 (0.68, 0.91) for men and 0.78 (0.68, 0.90) for women who drank at least five cups per day of green tea compared with non-drinkers. The association between black tea consumption and mortality was weak, with no clear trends noted across the categories of consumption.

Conclusions

In Asian populations, coffee consumption is associated with a lower risk of death overall and with lower risks of death from CVD and cancer. Green tea consumption is associated with lower risks of death from all causes and CVD.

Key Messages
  • In a pooled analysis of 12 prospective cohort studies from the Asia Cohort Consortium conducted in China, Japan, Korea and Singapore, coffee consumption was associated with a lower risk of death overall and lower risks of death from cardiovascular disease and cancer.

  • Green tea consumption was associated with lower risks of death from all causes and cardiovascular disease.

  • Green tea consumption did not decrease the risk of cancer-specific deaths.

Introduction

Coffee and tea are the two most commonly consumed beverages in the world. Coffee consumption has increased considerably over the past few decades in Asia. Food-balance data from the United Nations Food and Agricultural Organization (FAO) have shown an increase in coffee supply per capita per day available for human consumption, which increased from 1.93 g in 1970 to 10.65 g in 2013 in Japan and 0.07 g in 1970 to 3.06 g in 2013 in Korea.1 Asia has a long tradition of drinking tea, with green tea being one of the most popular types of tea consumed. China and Japan accounted for >90% of the worldwide green tea production in 2015–2017.2 The potential health benefits associated with drinking coffee and tea have drawn the attention of both researchers and the general public because coffee and tea are rich sources of antioxidants.3,4

Coffee consumption has been inversely associated with several chronic diseases, including type 2 diabetes,5 cardiovascular disease (CVD),6 stroke,7 Parkinson’s disease,8 advanced prostate cancer,9 endometrial cancer10 and liver cancer.11 Moreover, cohort studies in the USA12–15 and Europe16 have consistently found that coffee consumption was associated with lower risk of both all-cause and cause-specific mortality. A recent dose–response meta-analysis of three Japanese cohort studies reported that the largest reduction in the all-cause mortality was noted in the group consuming four cups of coffee per day but that there was no further reduction in the groups that consumed more than four cups per day.17 Comparatively, few cohort studies have explored the association of green tea consumption with all-cause or chronic disease-specific mortality. In Japanese populations, green tea consumption has been associated with a lower risk of gastric cancer in women but not in men.18 Five prospective Asian cohort studies, including three that were included in our analysis, evaluated the association between green tea consumption and mortality, but the results were not consistent for cause-specific mortality.19–23

The objective of this study was to evaluate the association of coffee, green tea and black tea consumption with all-cause and cause-specific mortality in 12 ethnically diverse Asian cohorts and these cohorts included nearly half a million men and women, involved considerable follow-up times and were used to study a broad range of dietary and other lifestyle exposures.

Methods

Study population

The Asia Cohort Consortium is an international collaboration that aims to understand the aetiology of common and rare diseases by combining data from multiple Asian prospective cohort studies.24,25 The members of the Asia Cohort Consortium include established cohorts or new cohorts that are in the development stage across Asia.25 The established cohorts include cohorts in Bangladesh, China, Japan, Korea, Singapore and Taiwan. Among these cohorts, the study presented here included data from 12 prospective cohort studies that agreed to participate in this project and represented the four countries that collected baseline information on coffee or tea consumption and follow-up data on mortality. Of these cohorts, seven were conducted in Japan, two in mainland China, two in Korea and one in Singapore. These studies included the Japan Public Health Center-based Prospective Study (JPHC) I,26 the JPHC II,26 Miyagi Cohort (Miyagi),27 the Ohsaki National Health Insurance Cohort Study (Ohsaki),27 the Life Span Study Cohort (RERF),28 the 3 Prefecture Miyagi,29 the 3 Prefecture Aichi,29 the Seoul Male Cancer Cohort (Seoul Male),30 the Korean Multi-center Cancer Cohort Study (KMCC),31 the Shanghai Men’s Health Study (SMHS),32 the Shanghai Women’s Health Study (SWHS)33 and the Singapore Chinese Health Study (SCHS).34 Among a total of 529 362 participants, we excluded those who were missing data regarding an entry date (n = 858). Thus, 248 050 men and 280 454 women were considered for the subsequent analyses. For the coffee analyses, we excluded those who did not provide information on coffee consumption (n = 198 180), resulting in an analytic sample of 330 324 participants (161 177 men and 169 147 women). For the green tea analyses, we excluded those who did not provide information on green tea consumption (n = 51 312), resulting in an analytic sample of 477 192 participants (225 622 men and 251 570 women). Pooled analysis of the Asia Cohort Consortium (ACC) cohorts was approved by the ethical committee of the National Cancer Center Japan (number 2014–041) and each study was approved by respective ethic committees overseeing the participating studies. Written or oral consent was obtained from all subjects who participated in the study.

Assessment of coffee and tea intake

Each study assessed the participants’ diets at baseline using a self-administered food-frequency questionnaire (FFQ),34–39 and coffee and tea intake was defined as the number of servings (i.e. one cup) consumed per day. Information on coffee consumption was collected in 10 studies (Supplementary Table S1, available as Supplementary data at IJE online). Seven studies (JPHC I, JPHC II, Miyagi, 3 Prefecture Miyagi, Ohsaki, RERF and SCHS) collected data on regular coffee and two studies (KMCC and Seoul Male) collected data on any type of coffee. One study (3 Prefecture Aichi) collected data on both ground coffee and instant coffee. The two Shanghai-based studies (SMHS and SWHS) did not query about coffee consumption due to the low consumption in the population. Information on tea consumption was collected in all 12 studies, with all studies assessing green tea consumption and all but two studies (Seoul Male and KMCC) assessing black tea consumption (Supplementary Table S1, available as Supplementary data at IJE online).

Coffee- and tea-consumption data were obtained from the participants using questions of four to nine predefined categories or by using open-ended questions. JPHC, Miyagi, 3 Prefecture Miyagi, Ohsaki and 3 Prefecture Aichi assessed the intake with similar categories: almost never or does not drink, sometimes drink or one or two, or three or four days per week, one to two cups per day, three to four cups per day and five or more cups per day. Korean and Singapore studies provided information on continuous servings or frequencies per day based on their FFQs; e.g. in the KMCC, continuous servings per day were provided based on frequencies written as none, once per month, two or three times per month, once per week, two or three times per week, four to six times per week, once per day, and two or more times per day. SMHS and SWHS asked participants about the types and amounts of tea they consumed. We categorized total coffee, green tea and black tea consumption data from each cohort using the following categories: almost none, less than one cup per day, one to under three cups per day, three to under five cups per day, and five or more cups per day. For black tea, we collapsed the top two categories into one because of an insufficient number of cases and therefore used four frequency categories: almost never, less than one cup per day, one to under three cups per day, and three or more cups per day. Six studies asked whether participants put sugar or milk/cream in their coffee or tea, but the majority of the studies did not ask the amount used.

Assessment of other environmental factors

Age, height, weight, cigarette-smoking history including intensity and duration, and alcohol-consumption data were collected in all 12 studies using structured questionnaires; body mass index [BMI, weight (kg)/height (m2)] was calculated using the baseline height and weight of the participants. Energy intake was assessed using a structured food-frequency questionnaire in eight of the cohort studies. Educational attainment was also assessed in eight studies. Individual-level data from all participating cohorts were collected and harmonized for the statistical analyses.

Death ascertainment

Data on all-cause and cause-specific mortality were obtained via linkage to death-certificate data or by active follow-ups. Our primary endpoint was death from any cause. Based on the ninth or tenth revision of the International Statistical Classification of Diseases (ICD-9 or 10), outcomes of interest included CVD (using ICD-10 codes I00–I99; ICD-9 codes 390–459) and cancer (ICD-10 codes C00–C97; ICD-9 codes 140–239). Deaths from other causes were defined as all other ICD codes.

Statistical analysis

The means (standard deviation) or proportions of the baseline characteristics based on the coffee or green tea consumption were calculated. Person-years of follow-up were estimated from the baseline entry date until the date of death, loss to follow-up (if applicable) or end of study follow-up, whichever came first. Multivariable Cox proportional-hazards regression models were used to obtain the study-specific hazard ratio (HR) estimates for coffee or tea intake and mortality, adjusting for potential confounding factors, including age at baseline (years; continuous), BMI (kg/m2; <18.5, 18.5–19.9, 20.0–24.9, 25.0–29.9 and ≥30.0), smoking status (never-smoker, former smoker, current smoker with <20 pack-years of smoking and current with ≥20 pack-years of smoking), alcohol intake (non-drinker, 1–155 g/week and ≥156 g/week), education level (less than secondary, secondary and more than secondary school) and energy intake (kcal/day; continuous). Coffee and tea intakes were mutually adjusted in the multivariable analyses. Because information on the education level was not available for some participants in eight studies that assessed this information, participants with missing data were coded with missing indicator variables. For energy intake, missing data were coded using the sex-specific median values of energy intake. The proportions of missing data observed were 2.4% for education level and 3.4% for energy intake among studies that assessed this information.

After calculating the study- and sex-specific HRs for each coffee and green tea intake category, we combined the loge HRs using a random-effects model.40 The study-specific HRs were weighted by the inverse of their variance. Studies that had no participants in a category were not included in the pooled estimate for that category. We estimated P-values for trends using the regression coefficients by treating the category midpoint of each category as an ordinal variable in each individual study and then combined the coefficients using a random-effects model. We tested for and quantified the heterogeneity of the HRs for the highest category across studies using the Q statistics.40

To evaluate for any potential effect modifiers, we conducted analyses that were stratified by the baseline age (<54 and ≥54 years; the median age of cohorts combined); baseline smoking status (never, former and current for men, never and ever for women); alcohol use (non-drinker and current drinker); BMI level (15–18.4, 18.5–24.9 and ≥25.0 kg/m2); education level (lower than secondary and secondary education or higher); baseline co-morbidities including diabetes, hypertension, cancer and coronary heart disease (yes and no); and energy intake (<median, ≥median kcal/day). We combined the data sets from each participating cohort into one harmonized data set for subgroup analyses because of the small number of deaths in each stratum per study (aggregated analysis). In the aggregated analysis, the study centre was adjusted as a stratification variable. In the sensitivity analysis, we estimated the HRs of coffee and green tea with the all-cause mortality and excluded participants who died within 3 years of follow-up after baseline.

All statistical tests were two-sided and P-values of <0.05 were considered significant. All analyses were carried out using SAS version 9.3 (SAS Institute, Cary, NC) and STATA version 14.0 (Stata Corporation, Texas, USA).

Results

In total, 94 744 deaths were identified during the follow-up periods, which ranged from an average of 6.5 to 22.7 years (Supplementary Table S2, available as Supplementary data at IJE online). Of these, 29 968 deaths from CVD and 32 765 deaths from cancer were documented. Participants’ mean ages ranged from 49.2 to 60.7 years and the mean BMIs ranged from 21.9 to 23.7 kg/m2 in men and 21.9 to 24.1 kg/m2 in women (Table 1).

Table 1

Major characteristics of participating cohorts

Cohort (sex)Baseline cohort sizeAge, years, mean (sd)BMI, kg/m2, mean (sd)Total energy intake, kcal/day, mean (sd)Current smoker, %Current alcohol drinker, %Educational level (more than secondary), %
JPHC I (M)2620 46449.5 (5.9)23.6 (2.8)2188.1 (694.7)53.279.312.7
JPHC I (F)2622 30749.7 (5.9)23.6 (3.1)1430.8 (388.0)5.823.210.5
JPHC II (M)2626 41754.0 (8.8)23.5 (2.9)1748.5 (587.0)51.174.7
JPHC II (F)2629 29354.4 (8.8)23.4 (3.2)1102.0 (291.4)6.419.7
Miyagi Cohort Study (M)2721 50451.7 (7.6)23.5 (2.8)1590.0 (444.2)61.477.014.2
Miyagi Cohort Study (F)2723 36352.2 (7.4)23.7 (3.1)1560.5 (448.3)9.025.812.9
3-Prefecture Miyagi (M)2913 29756.5 (11.0)23.1 (3.3)54.176.5
3-Prefecture Miyagi (F)2916 25757.1 (11.2)23.4 (3.6)9.332.5
Ohsaki (M)2723 00059.5 (10.6)23.3 (3.1)1687.1 (442.1)54.872.17.6
Ohsaki (F)2724 71060.7 (9.9)23.8 (3.5)1281.4 (316.4)8.622.88.0
RERF (M)2819 52352.4 (11.1)21.9 (3.2)78.983.918.8
RERF (F)2829 44851.8 (15.1)22.0 (3.8)13.827.76.5
3-Prefecture Aichi (M)2915 27155.7 (11.0)22.3 (2.8)55.078.9
3-Prefecture Aichi (F)2916 94156.6 (11.4)21.9 (3.1)12.742.9
Seoul male (M)3013 95349.2 (5.2)23.4 (2.4)2020.4 (931.4)50.466.753.7
KMCC (M)31532454.9 (12.7)23.0 (2.9)58.063.6
KMCC (F)31795855.1 (12.0)24.1 (3.4)7.517.9
SMHS (M)3261 34355.3 (9.7)23.7 (3.1)1909.0 (484.7)58.729.359.8
SWHS (F)3374 87452.6 (9.1)24.0 (3.4)1680.8 (402.8)2.41.941.5
SCHS (M)3427 95456.7 (8.0)23.0 (3.2)1749.5 (608.8)36.231.18.0
SCHS (F)3435 30356.3 (8.0)23.2 (3.3)1398.8 (472.0)6.39.02.9
Total528 504*54.3 (10.0)23.3 (3.3)1645.0 (561.1)30.226.625.0
Cohort (sex)Baseline cohort sizeAge, years, mean (sd)BMI, kg/m2, mean (sd)Total energy intake, kcal/day, mean (sd)Current smoker, %Current alcohol drinker, %Educational level (more than secondary), %
JPHC I (M)2620 46449.5 (5.9)23.6 (2.8)2188.1 (694.7)53.279.312.7
JPHC I (F)2622 30749.7 (5.9)23.6 (3.1)1430.8 (388.0)5.823.210.5
JPHC II (M)2626 41754.0 (8.8)23.5 (2.9)1748.5 (587.0)51.174.7
JPHC II (F)2629 29354.4 (8.8)23.4 (3.2)1102.0 (291.4)6.419.7
Miyagi Cohort Study (M)2721 50451.7 (7.6)23.5 (2.8)1590.0 (444.2)61.477.014.2
Miyagi Cohort Study (F)2723 36352.2 (7.4)23.7 (3.1)1560.5 (448.3)9.025.812.9
3-Prefecture Miyagi (M)2913 29756.5 (11.0)23.1 (3.3)54.176.5
3-Prefecture Miyagi (F)2916 25757.1 (11.2)23.4 (3.6)9.332.5
Ohsaki (M)2723 00059.5 (10.6)23.3 (3.1)1687.1 (442.1)54.872.17.6
Ohsaki (F)2724 71060.7 (9.9)23.8 (3.5)1281.4 (316.4)8.622.88.0
RERF (M)2819 52352.4 (11.1)21.9 (3.2)78.983.918.8
RERF (F)2829 44851.8 (15.1)22.0 (3.8)13.827.76.5
3-Prefecture Aichi (M)2915 27155.7 (11.0)22.3 (2.8)55.078.9
3-Prefecture Aichi (F)2916 94156.6 (11.4)21.9 (3.1)12.742.9
Seoul male (M)3013 95349.2 (5.2)23.4 (2.4)2020.4 (931.4)50.466.753.7
KMCC (M)31532454.9 (12.7)23.0 (2.9)58.063.6
KMCC (F)31795855.1 (12.0)24.1 (3.4)7.517.9
SMHS (M)3261 34355.3 (9.7)23.7 (3.1)1909.0 (484.7)58.729.359.8
SWHS (F)3374 87452.6 (9.1)24.0 (3.4)1680.8 (402.8)2.41.941.5
SCHS (M)3427 95456.7 (8.0)23.0 (3.2)1749.5 (608.8)36.231.18.0
SCHS (F)3435 30356.3 (8.0)23.2 (3.3)1398.8 (472.0)6.39.02.9
Total528 504*54.3 (10.0)23.3 (3.3)1645.0 (561.1)30.226.625.0

M, men; F, women; JPHC, Japan Public Health Center-based Prospective Study; Miyagi, Miyagi Cohort Study; Ohsaki, Ohsaki National Health Insurance Cohort Study; RERF, Radiation Effects Research Foundation Life Span Study; Seoul Male, Seoul Male Cohort Study; KMCC, Korean Multi-center Cancer Cohort Study; SMHS, Shanghai Men's Health Study; SWHS, Shanghai Women’s Health Study; SCHS, Singapore Chinese Health Study; BMI, body mass index.

Table 1

Major characteristics of participating cohorts

Cohort (sex)Baseline cohort sizeAge, years, mean (sd)BMI, kg/m2, mean (sd)Total energy intake, kcal/day, mean (sd)Current smoker, %Current alcohol drinker, %Educational level (more than secondary), %
JPHC I (M)2620 46449.5 (5.9)23.6 (2.8)2188.1 (694.7)53.279.312.7
JPHC I (F)2622 30749.7 (5.9)23.6 (3.1)1430.8 (388.0)5.823.210.5
JPHC II (M)2626 41754.0 (8.8)23.5 (2.9)1748.5 (587.0)51.174.7
JPHC II (F)2629 29354.4 (8.8)23.4 (3.2)1102.0 (291.4)6.419.7
Miyagi Cohort Study (M)2721 50451.7 (7.6)23.5 (2.8)1590.0 (444.2)61.477.014.2
Miyagi Cohort Study (F)2723 36352.2 (7.4)23.7 (3.1)1560.5 (448.3)9.025.812.9
3-Prefecture Miyagi (M)2913 29756.5 (11.0)23.1 (3.3)54.176.5
3-Prefecture Miyagi (F)2916 25757.1 (11.2)23.4 (3.6)9.332.5
Ohsaki (M)2723 00059.5 (10.6)23.3 (3.1)1687.1 (442.1)54.872.17.6
Ohsaki (F)2724 71060.7 (9.9)23.8 (3.5)1281.4 (316.4)8.622.88.0
RERF (M)2819 52352.4 (11.1)21.9 (3.2)78.983.918.8
RERF (F)2829 44851.8 (15.1)22.0 (3.8)13.827.76.5
3-Prefecture Aichi (M)2915 27155.7 (11.0)22.3 (2.8)55.078.9
3-Prefecture Aichi (F)2916 94156.6 (11.4)21.9 (3.1)12.742.9
Seoul male (M)3013 95349.2 (5.2)23.4 (2.4)2020.4 (931.4)50.466.753.7
KMCC (M)31532454.9 (12.7)23.0 (2.9)58.063.6
KMCC (F)31795855.1 (12.0)24.1 (3.4)7.517.9
SMHS (M)3261 34355.3 (9.7)23.7 (3.1)1909.0 (484.7)58.729.359.8
SWHS (F)3374 87452.6 (9.1)24.0 (3.4)1680.8 (402.8)2.41.941.5
SCHS (M)3427 95456.7 (8.0)23.0 (3.2)1749.5 (608.8)36.231.18.0
SCHS (F)3435 30356.3 (8.0)23.2 (3.3)1398.8 (472.0)6.39.02.9
Total528 504*54.3 (10.0)23.3 (3.3)1645.0 (561.1)30.226.625.0
Cohort (sex)Baseline cohort sizeAge, years, mean (sd)BMI, kg/m2, mean (sd)Total energy intake, kcal/day, mean (sd)Current smoker, %Current alcohol drinker, %Educational level (more than secondary), %
JPHC I (M)2620 46449.5 (5.9)23.6 (2.8)2188.1 (694.7)53.279.312.7
JPHC I (F)2622 30749.7 (5.9)23.6 (3.1)1430.8 (388.0)5.823.210.5
JPHC II (M)2626 41754.0 (8.8)23.5 (2.9)1748.5 (587.0)51.174.7
JPHC II (F)2629 29354.4 (8.8)23.4 (3.2)1102.0 (291.4)6.419.7
Miyagi Cohort Study (M)2721 50451.7 (7.6)23.5 (2.8)1590.0 (444.2)61.477.014.2
Miyagi Cohort Study (F)2723 36352.2 (7.4)23.7 (3.1)1560.5 (448.3)9.025.812.9
3-Prefecture Miyagi (M)2913 29756.5 (11.0)23.1 (3.3)54.176.5
3-Prefecture Miyagi (F)2916 25757.1 (11.2)23.4 (3.6)9.332.5
Ohsaki (M)2723 00059.5 (10.6)23.3 (3.1)1687.1 (442.1)54.872.17.6
Ohsaki (F)2724 71060.7 (9.9)23.8 (3.5)1281.4 (316.4)8.622.88.0
RERF (M)2819 52352.4 (11.1)21.9 (3.2)78.983.918.8
RERF (F)2829 44851.8 (15.1)22.0 (3.8)13.827.76.5
3-Prefecture Aichi (M)2915 27155.7 (11.0)22.3 (2.8)55.078.9
3-Prefecture Aichi (F)2916 94156.6 (11.4)21.9 (3.1)12.742.9
Seoul male (M)3013 95349.2 (5.2)23.4 (2.4)2020.4 (931.4)50.466.753.7
KMCC (M)31532454.9 (12.7)23.0 (2.9)58.063.6
KMCC (F)31795855.1 (12.0)24.1 (3.4)7.517.9
SMHS (M)3261 34355.3 (9.7)23.7 (3.1)1909.0 (484.7)58.729.359.8
SWHS (F)3374 87452.6 (9.1)24.0 (3.4)1680.8 (402.8)2.41.941.5
SCHS (M)3427 95456.7 (8.0)23.0 (3.2)1749.5 (608.8)36.231.18.0
SCHS (F)3435 30356.3 (8.0)23.2 (3.3)1398.8 (472.0)6.39.02.9
Total528 504*54.3 (10.0)23.3 (3.3)1645.0 (561.1)30.226.625.0

M, men; F, women; JPHC, Japan Public Health Center-based Prospective Study; Miyagi, Miyagi Cohort Study; Ohsaki, Ohsaki National Health Insurance Cohort Study; RERF, Radiation Effects Research Foundation Life Span Study; Seoul Male, Seoul Male Cohort Study; KMCC, Korean Multi-center Cancer Cohort Study; SMHS, Shanghai Men's Health Study; SWHS, Shanghai Women’s Health Study; SCHS, Singapore Chinese Health Study; BMI, body mass index.

The proportions of coffee drinkers and green tea drinkers ranged from 55.6% to 84.2% and 1.5% to 97.8%, respectively (Supplementary Table S1, available as Supplementary data at IJE online). Frequent coffee and green tea drinkers were more likely to smoke cigarettes than non-drinkers (Supplementary Table S3, available as Supplementary data at IJE online). A lower proportion of men and women who almost never drank coffee vs those who drank more than five cups per day were current smokers (43.1% vs 74.6% and 5.9% vs 23.7%, respectively). Similarly, a lower proportion of men and women who almost never drank green tea vs those who drank more than five cups per day were current smokers (45.9% vs 59.1% and 4.1% vs 9.9%, respectively). Women, but not men, who reported drinking coffee were more likely to drink alcohol, and both women and men who reported drinking green tea were more likely to drink alcohol. Among women, regular coffee drinkers were more educated than non-drinkers, but this trend was not found in men. Participants who drank green tea had completed less education on average than non-drinkers.

Coffee consumption and all-cause and cause-specific mortality

We found that a higher coffee consumption was associated with a lower risk of all-cause mortality as well as lower risks of mortality from CVD and cancer among men and women; the HRs [95% confidence intervals (CIs)] for all-cause mortality when comparing participants who drank five or more cups per day to coffee non-drinkers were 0.76 (0.71, 0.83; P for trend <0.001) among men and 0.72 (0.63, 0.81; P for trend <0.001) among women (Table 2). We observed that the reduction in risk was slightly greater at five or more cups compared with less than one cup per day, but there was a lower risk of mortality among those consuming coffee at both lower and higher levels. The risk of CVD-specific mortality decreased with increasing coffee consumption among both men and women; the HRs (95% CIs) for five or more cups of coffee per day were 0.73 (0.62, 0.85) for men (P for trend <0.001) and 0.73 (0.58, 0.92) for women (P for trend <0.001) compared with those who did not drink coffee. Both men and women who regularly drank coffee had a lower risk of cancer mortality than non-drinkers. The exclusion of the men and women who died during the first 3 years of follow-up resulted in similar HR estimates for all-cause, CVD- and other cause-specific mortalities, but there was an attenuation in the HR estimates for the cancer-specific mortality (Supplementary Table S4, available as Supplementary data at IJE online).

Table 2

Pooled hazard ratios (HRs)a and 95% confidence intervals (CIs) of all-cause and cause-specific mortality according to coffee consumption

Coffee consumption
Almost neverLess than 1 cup/day1 to <3 cups/day3 to <5 cups/day≥5 cups/daybP for trendP for heterogeneityc
Men
 No. of participants35 40947 36962 69611 5374166
 Total person-years434 288603 860786 344133 87448 269
 All-cause-specific mortality (no. of cases)10 21610 60511 5801797706
 HR (95% CI)1.00 (Ref.)0.83 (0.79–0.87)0.78 (0.73–0.83)0.76 (0.67–0.85)0.76 (0.71–0.83)<0.0010.934
 CVD-specific mortality (no. of cases)284328583286430180
 HR (95% CI)1.00 (Ref.)0.80 (0.75–0.85)0.77 (0.68–0.87)0.67 (0.55–0.81)0.73 (0.62–0.85)<0.0010.818
 Cancer-specific mortality (no. of cases)356139944633776291
 HR (95% CI)1.00 (Ref.)0.88 (0.84–0.93)0.86 (0.81–0.92)0.88 (0.79–0.97)0.85 (0.75–0.97)0.0040.544
 Other cause-specific mortality (no. of cases)381237533661591235
 HR (95% CI)1.00 (Ref.)0.81 (0.75–0.87)0.73 (0.67–0.80)0.66 (0.55–0.81)0.72 (0.62–0.84)<0.0010.296
Women
 No. of participants44 70950 86061 96086232995
 Total person-years565 582683 414802 334104 83036 053
 All-cause-specific mortality (no. of cases)946174447253639265
 HR (95% CI)1.00 (Ref.)0.86 (0.82–0.90)0.80 (0.72–0.89)0.65 (0.54–0.78)0.72 (0.63–0.81)<0.0010.686
 CVD-specific mortality (no. of cases)33592385225017275
 HR (95% CI)1.00 (Ref.)0.83 (0.76–0.92)0.75 (0.63–0.89)0.61 (0.48–0.79)0.73 (0.58–0.92)<0.0010.452
 Cancer-specific mortality (no. of cases)260724662764282110
 HR (95% CI)1.00 (Ref.)0.91 (0.86–0.97)0.90 (0.81–1.00)0.75 (0.57–0.98)0.81 (0.63–1.04)0.0390.208
 Other cause-specific mortality (no. of cases)34952593223918580
 HR (95% CI)1.00 (Ref.)0.85 (0.79–0.91)0.75 (0.67–0.84)0.60 (0.50–0.72)0.68 (0.54–0.85)<0.0010.870
Coffee consumption
Almost neverLess than 1 cup/day1 to <3 cups/day3 to <5 cups/day≥5 cups/daybP for trendP for heterogeneityc
Men
 No. of participants35 40947 36962 69611 5374166
 Total person-years434 288603 860786 344133 87448 269
 All-cause-specific mortality (no. of cases)10 21610 60511 5801797706
 HR (95% CI)1.00 (Ref.)0.83 (0.79–0.87)0.78 (0.73–0.83)0.76 (0.67–0.85)0.76 (0.71–0.83)<0.0010.934
 CVD-specific mortality (no. of cases)284328583286430180
 HR (95% CI)1.00 (Ref.)0.80 (0.75–0.85)0.77 (0.68–0.87)0.67 (0.55–0.81)0.73 (0.62–0.85)<0.0010.818
 Cancer-specific mortality (no. of cases)356139944633776291
 HR (95% CI)1.00 (Ref.)0.88 (0.84–0.93)0.86 (0.81–0.92)0.88 (0.79–0.97)0.85 (0.75–0.97)0.0040.544
 Other cause-specific mortality (no. of cases)381237533661591235
 HR (95% CI)1.00 (Ref.)0.81 (0.75–0.87)0.73 (0.67–0.80)0.66 (0.55–0.81)0.72 (0.62–0.84)<0.0010.296
Women
 No. of participants44 70950 86061 96086232995
 Total person-years565 582683 414802 334104 83036 053
 All-cause-specific mortality (no. of cases)946174447253639265
 HR (95% CI)1.00 (Ref.)0.86 (0.82–0.90)0.80 (0.72–0.89)0.65 (0.54–0.78)0.72 (0.63–0.81)<0.0010.686
 CVD-specific mortality (no. of cases)33592385225017275
 HR (95% CI)1.00 (Ref.)0.83 (0.76–0.92)0.75 (0.63–0.89)0.61 (0.48–0.79)0.73 (0.58–0.92)<0.0010.452
 Cancer-specific mortality (no. of cases)260724662764282110
 HR (95% CI)1.00 (Ref.)0.91 (0.86–0.97)0.90 (0.81–1.00)0.75 (0.57–0.98)0.81 (0.63–1.04)0.0390.208
 Other cause-specific mortality (no. of cases)34952593223918580
 HR (95% CI)1.00 (Ref.)0.85 (0.79–0.91)0.75 (0.67–0.84)0.60 (0.50–0.72)0.68 (0.54–0.85)<0.0010.870

CVD, cardiovascular disease.

a

Pooled HRs (95% CIs) were calculated using a random-effects model after study-specific HRs (95% CIs) were obtained using Cox proportional-hazards regression models. Multivariate models were adjusted for baseline age (years; continuous), body mass index (BMI, kg/m2; <18.5, 18.5–19.9, 20.0–24.9, 25.0–29.9 and ≥30.0) and smoking status (never-smoker, former smoker, current with <20 pack-years of smoking and current with ≥20 pack-years of smoking), alcohol intake (non-drinker, 1–155 and ≥156 g/week), educational level (less than secondary, secondary and more than secondary school graduate), energy intake (kcal/day; continuous) and green tea consumption (almost never, <1, 1 to <3, 3 to <5 and ≥5 cups/day).

b

The RERF (Radiation Effects Research Foundation Life Span Study), KMCC (Korean Multi-center Cancer Cohort Study) and Seoul male cohort were excluded from the 3 to <5 and ≥5 cups/day categories because there was no one in these categories due to the way that frequency of intake was assessed in these studies (the highest category was >2 cups/day).

c

P for heterogeneity for the highest category was calculated using the Q statistic.

Table 2

Pooled hazard ratios (HRs)a and 95% confidence intervals (CIs) of all-cause and cause-specific mortality according to coffee consumption

Coffee consumption
Almost neverLess than 1 cup/day1 to <3 cups/day3 to <5 cups/day≥5 cups/daybP for trendP for heterogeneityc
Men
 No. of participants35 40947 36962 69611 5374166
 Total person-years434 288603 860786 344133 87448 269
 All-cause-specific mortality (no. of cases)10 21610 60511 5801797706
 HR (95% CI)1.00 (Ref.)0.83 (0.79–0.87)0.78 (0.73–0.83)0.76 (0.67–0.85)0.76 (0.71–0.83)<0.0010.934
 CVD-specific mortality (no. of cases)284328583286430180
 HR (95% CI)1.00 (Ref.)0.80 (0.75–0.85)0.77 (0.68–0.87)0.67 (0.55–0.81)0.73 (0.62–0.85)<0.0010.818
 Cancer-specific mortality (no. of cases)356139944633776291
 HR (95% CI)1.00 (Ref.)0.88 (0.84–0.93)0.86 (0.81–0.92)0.88 (0.79–0.97)0.85 (0.75–0.97)0.0040.544
 Other cause-specific mortality (no. of cases)381237533661591235
 HR (95% CI)1.00 (Ref.)0.81 (0.75–0.87)0.73 (0.67–0.80)0.66 (0.55–0.81)0.72 (0.62–0.84)<0.0010.296
Women
 No. of participants44 70950 86061 96086232995
 Total person-years565 582683 414802 334104 83036 053
 All-cause-specific mortality (no. of cases)946174447253639265
 HR (95% CI)1.00 (Ref.)0.86 (0.82–0.90)0.80 (0.72–0.89)0.65 (0.54–0.78)0.72 (0.63–0.81)<0.0010.686
 CVD-specific mortality (no. of cases)33592385225017275
 HR (95% CI)1.00 (Ref.)0.83 (0.76–0.92)0.75 (0.63–0.89)0.61 (0.48–0.79)0.73 (0.58–0.92)<0.0010.452
 Cancer-specific mortality (no. of cases)260724662764282110
 HR (95% CI)1.00 (Ref.)0.91 (0.86–0.97)0.90 (0.81–1.00)0.75 (0.57–0.98)0.81 (0.63–1.04)0.0390.208
 Other cause-specific mortality (no. of cases)34952593223918580
 HR (95% CI)1.00 (Ref.)0.85 (0.79–0.91)0.75 (0.67–0.84)0.60 (0.50–0.72)0.68 (0.54–0.85)<0.0010.870
Coffee consumption
Almost neverLess than 1 cup/day1 to <3 cups/day3 to <5 cups/day≥5 cups/daybP for trendP for heterogeneityc
Men
 No. of participants35 40947 36962 69611 5374166
 Total person-years434 288603 860786 344133 87448 269
 All-cause-specific mortality (no. of cases)10 21610 60511 5801797706
 HR (95% CI)1.00 (Ref.)0.83 (0.79–0.87)0.78 (0.73–0.83)0.76 (0.67–0.85)0.76 (0.71–0.83)<0.0010.934
 CVD-specific mortality (no. of cases)284328583286430180
 HR (95% CI)1.00 (Ref.)0.80 (0.75–0.85)0.77 (0.68–0.87)0.67 (0.55–0.81)0.73 (0.62–0.85)<0.0010.818
 Cancer-specific mortality (no. of cases)356139944633776291
 HR (95% CI)1.00 (Ref.)0.88 (0.84–0.93)0.86 (0.81–0.92)0.88 (0.79–0.97)0.85 (0.75–0.97)0.0040.544
 Other cause-specific mortality (no. of cases)381237533661591235
 HR (95% CI)1.00 (Ref.)0.81 (0.75–0.87)0.73 (0.67–0.80)0.66 (0.55–0.81)0.72 (0.62–0.84)<0.0010.296
Women
 No. of participants44 70950 86061 96086232995
 Total person-years565 582683 414802 334104 83036 053
 All-cause-specific mortality (no. of cases)946174447253639265
 HR (95% CI)1.00 (Ref.)0.86 (0.82–0.90)0.80 (0.72–0.89)0.65 (0.54–0.78)0.72 (0.63–0.81)<0.0010.686
 CVD-specific mortality (no. of cases)33592385225017275
 HR (95% CI)1.00 (Ref.)0.83 (0.76–0.92)0.75 (0.63–0.89)0.61 (0.48–0.79)0.73 (0.58–0.92)<0.0010.452
 Cancer-specific mortality (no. of cases)260724662764282110
 HR (95% CI)1.00 (Ref.)0.91 (0.86–0.97)0.90 (0.81–1.00)0.75 (0.57–0.98)0.81 (0.63–1.04)0.0390.208
 Other cause-specific mortality (no. of cases)34952593223918580
 HR (95% CI)1.00 (Ref.)0.85 (0.79–0.91)0.75 (0.67–0.84)0.60 (0.50–0.72)0.68 (0.54–0.85)<0.0010.870

CVD, cardiovascular disease.

a

Pooled HRs (95% CIs) were calculated using a random-effects model after study-specific HRs (95% CIs) were obtained using Cox proportional-hazards regression models. Multivariate models were adjusted for baseline age (years; continuous), body mass index (BMI, kg/m2; <18.5, 18.5–19.9, 20.0–24.9, 25.0–29.9 and ≥30.0) and smoking status (never-smoker, former smoker, current with <20 pack-years of smoking and current with ≥20 pack-years of smoking), alcohol intake (non-drinker, 1–155 and ≥156 g/week), educational level (less than secondary, secondary and more than secondary school graduate), energy intake (kcal/day; continuous) and green tea consumption (almost never, <1, 1 to <3, 3 to <5 and ≥5 cups/day).

b

The RERF (Radiation Effects Research Foundation Life Span Study), KMCC (Korean Multi-center Cancer Cohort Study) and Seoul male cohort were excluded from the 3 to <5 and ≥5 cups/day categories because there was no one in these categories due to the way that frequency of intake was assessed in these studies (the highest category was >2 cups/day).

c

P for heterogeneity for the highest category was calculated using the Q statistic.

When we stratified by risk factors, including baseline age, smoking status, alcohol use, BMI, education level and history of co-morbidities, the inverse associations for the all-cause mortality were consistent across the strata in both men (Figure 1) and women (Figure 2).

Hazard ratios (HRs) and 95% confidence intervals (CIs) for at least five cups per day coffee consumption vs non-drinkers in relation to mortality by potential interaction factors among men. Aggregated HRs (95% CIs) were calculated using Cox proportional-hazards regression models. The associations between coffee consumption and mortality were analysed in subgroups of baseline age (<54 and ≥54 years; the median age of cohorts combined), baseline smoking status (never, former and current for men), alcohol use (non-drinker and current drinker), body mass index (BMI) level (15–18.4, 18.5–24.9 and ≥25.0 kg/m2), education level (lower than secondary and secondary education or higher) and baseline co-morbidities including diabetes, hypertension, cancer and coronary heart disease (yes and no). Models were adjusted for baseline age (years; continuous), BMI (kg/m2; <18.5, 18.5–19.9, 20.0–24.9, 25.0–29.9 and ≥30.0), smoking status (never smoked, former smoker, current with <20 pack-years of smoking and current with ≥20 pack-years of smoking), alcohol consumption (non-drinker and current drinker), educational level (less than secondary, secondary and more than secondary school graduate), energy intake (kcal/day; continuous) and green tea consumption (almost never, <1, 1 to <3, 3 to <5 and ≥5 cups/day). Study centre was adjusted as a stratifying variable.
Figure 1

Hazard ratios (HRs) and 95% confidence intervals (CIs) for at least five cups per day coffee consumption vs non-drinkers in relation to mortality by potential interaction factors among men. Aggregated HRs (95% CIs) were calculated using Cox proportional-hazards regression models. The associations between coffee consumption and mortality were analysed in subgroups of baseline age (<54 and ≥54 years; the median age of cohorts combined), baseline smoking status (never, former and current for men), alcohol use (non-drinker and current drinker), body mass index (BMI) level (15–18.4, 18.5–24.9 and ≥25.0 kg/m2), education level (lower than secondary and secondary education or higher) and baseline co-morbidities including diabetes, hypertension, cancer and coronary heart disease (yes and no). Models were adjusted for baseline age (years; continuous), BMI (kg/m2; <18.5, 18.5–19.9, 20.0–24.9, 25.0–29.9 and ≥30.0), smoking status (never smoked, former smoker, current with <20 pack-years of smoking and current with ≥20 pack-years of smoking), alcohol consumption (non-drinker and current drinker), educational level (less than secondary, secondary and more than secondary school graduate), energy intake (kcal/day; continuous) and green tea consumption (almost never, <1, 1 to <3, 3 to <5 and ≥5 cups/day). Study centre was adjusted as a stratifying variable.

Hazard ratios (HRs) and 95% confidence intervals (CIs) for at least five cups per day of coffee consumption vs non-drinkers in relation to mortality by potential interaction factors among women. Aggregated HRs (95% CIs) were calculated using Cox proportional-hazards regression models. Study centre was adjusted as a stratifying variable. The associations between coffee consumption and mortality were analysed in subgroups of baseline age (<54 and ≥54 years; the median age of cohorts combined), baseline smoking status (never and ever for women), alcohol use (non-drinker and current drinker), body mass index (BMI) level (15–18.4, 18.5–24.9 and ≥25.0 kg/m2), education level (lower than secondary and secondary education or higher) and baseline co-morbidities including diabetes, hypertension, cancer and coronary heart disease (yes and no). Models were adjusted for baseline age (years; continuous), BMI (kg/m2; <18.5, 18.5–19.9, 20.0–24.9, 25.0–29.9 and ≥30.0) and smoking status (never smoked, former smoker, current with <20 pack-years of smoking and current with ≥20 pack-years of smoking), alcohol drinking (non-drinker and current drinker), educational level (less than secondary, secondary and more than secondary school graduate), energy intake (kcal/day; continuous) and green tea consumption (almost never, <1, 1 to <3, 3 to <5 and ≥5 cups/day).
Figure 2

Hazard ratios (HRs) and 95% confidence intervals (CIs) for at least five cups per day of coffee consumption vs non-drinkers in relation to mortality by potential interaction factors among women. Aggregated HRs (95% CIs) were calculated using Cox proportional-hazards regression models. Study centre was adjusted as a stratifying variable. The associations between coffee consumption and mortality were analysed in subgroups of baseline age (<54 and ≥54 years; the median age of cohorts combined), baseline smoking status (never and ever for women), alcohol use (non-drinker and current drinker), body mass index (BMI) level (15–18.4, 18.5–24.9 and ≥25.0 kg/m2), education level (lower than secondary and secondary education or higher) and baseline co-morbidities including diabetes, hypertension, cancer and coronary heart disease (yes and no). Models were adjusted for baseline age (years; continuous), BMI (kg/m2; <18.5, 18.5–19.9, 20.0–24.9, 25.0–29.9 and ≥30.0) and smoking status (never smoked, former smoker, current with <20 pack-years of smoking and current with ≥20 pack-years of smoking), alcohol drinking (non-drinker and current drinker), educational level (less than secondary, secondary and more than secondary school graduate), energy intake (kcal/day; continuous) and green tea consumption (almost never, <1, 1 to <3, 3 to <5 and ≥5 cups/day).

Green tea consumption and all-cause and cause-specific mortality

An increased green tea intake was associated with decreased risk of CVD- and other cause-specific mortality but not cancer-specific mortality (Table 3). The HRs (95% CIs) for CVD-specific mortality when comparing five or more cups per day to green tea non-drinkers were 0.79 (0.68, 0.91; P for trend <0.001) among men and 0.78 (0.68–0.90; P for trend 0.014) among women. There was no association between green tea consumption and cancer-specific mortality among men. However, women who consumed five or more cups per day had a 12% higher risk of cancer-specific deaths than non-drinkers. The results were similar when we excluded participants who died during the first 3 years of follow-up (Supplementary Table S5, available as Supplementary data at IJE online) or in subgroup analyses stratified by baseline age, smoking status, alcohol drinking, BMI, education level and history of co-morbidities in men (Supplementary Figure S1, available as Supplementary data at IJE online) and women (Supplementary Figure S2, available as Supplementary data at IJE online). The associations did not vary by energy intake (Supplementary Table S6, available as Supplementary data at IJE online).

Table 3

Pooled hazard ratios (HRs)a and 95% confidence intervals (CIs) of all-cause and cause-specific mortality according to green tea consumption

Green tea consumption
Almost neverLess than 1 cup/day1 to <3 cups/day3 to <5 cups/day≥5 cups/daybP for trendP for heterogeneityc
Men
 No. of participants56 86741 26555 88736 39935 204
 Total person-years606 553480 888636 679396 522437 025
 All-cause-specific mortality (no. of cases)89246649966669529112
 HR (95% CI)1.00 (Ref.)0.91 (0.87–0.95)0.91 (0.86–0.96)0.91 (0.82–1.02)0.91 (0.82–1.00)0.011<0.001
 CVD-specific mortality (no. of cases)26711980273217992437
 HR (95% CI)1.00 (Ref.)0.90 (0.83–0.98)0.85 (0.80–0.91)0.83 (0.70–0.98)0.79 (0.68–0.91)<0.0010.007
 Cancer-specific mortality (no. of cases)32762459367926723630
 HR (95% CI)1.00 (Ref.)0.94 (0.87–1.02)0.99 (0.90–1.08)1.02 (0.91–1.16)1.07 (0.92–1.25)0.234<0.001
 Other cause-specific mortality (no. of cases)29772210325524813045
 HR (95% CI)1.00 (Ref.)0.88 (0.83–0.94)0.89 (0.84–0.95)0.91 (0.84–0.97)0.83 (0.76–0.90)<0.0010.338
Women
 No. of participants108 50231 47036 59632 50842 494
 Total person-years1 425 367391 362549 387393 930563 773
 All-cause-specific mortality (no. of cases)11 1564112689843447188
 HR (95% CI)1.00 (Ref.)0.94 (0.89–0.98)0.91 (0.86–0.95)0.87 (0.78–0.97)0.86 (0.75–1.00)0.032<0.001
 CVD-specific mortality (no. of cases)36391350247913262437
 HR (95% CI)1.00 (Ref.)0.88 (0.82–0.95)0.88 (0.80–0.96)0.80 (0.70–0.90)0.78 (0.68–0.90)0.0140.061
 Cancer-specific mortality (no. of cases)43141381196214822352
 HR (95% CI)1.00 (Ref.)1.01 (0.93–1.11)0.99 (0.90–1.09)1.08 (0.97–1.20)1.12 (1.01–1.24)0.4260.625
 Other cause-specific mortality (no. of cases)32031381245715362399
 HR (95% CI)1.00 (Ref.)0.92 (0.85–0.99)0.87 (0.80–0.95)0.80 (0.68–0.94)0.76 (0.61–0.94)0.001<0.001
Green tea consumption
Almost neverLess than 1 cup/day1 to <3 cups/day3 to <5 cups/day≥5 cups/daybP for trendP for heterogeneityc
Men
 No. of participants56 86741 26555 88736 39935 204
 Total person-years606 553480 888636 679396 522437 025
 All-cause-specific mortality (no. of cases)89246649966669529112
 HR (95% CI)1.00 (Ref.)0.91 (0.87–0.95)0.91 (0.86–0.96)0.91 (0.82–1.02)0.91 (0.82–1.00)0.011<0.001
 CVD-specific mortality (no. of cases)26711980273217992437
 HR (95% CI)1.00 (Ref.)0.90 (0.83–0.98)0.85 (0.80–0.91)0.83 (0.70–0.98)0.79 (0.68–0.91)<0.0010.007
 Cancer-specific mortality (no. of cases)32762459367926723630
 HR (95% CI)1.00 (Ref.)0.94 (0.87–1.02)0.99 (0.90–1.08)1.02 (0.91–1.16)1.07 (0.92–1.25)0.234<0.001
 Other cause-specific mortality (no. of cases)29772210325524813045
 HR (95% CI)1.00 (Ref.)0.88 (0.83–0.94)0.89 (0.84–0.95)0.91 (0.84–0.97)0.83 (0.76–0.90)<0.0010.338
Women
 No. of participants108 50231 47036 59632 50842 494
 Total person-years1 425 367391 362549 387393 930563 773
 All-cause-specific mortality (no. of cases)11 1564112689843447188
 HR (95% CI)1.00 (Ref.)0.94 (0.89–0.98)0.91 (0.86–0.95)0.87 (0.78–0.97)0.86 (0.75–1.00)0.032<0.001
 CVD-specific mortality (no. of cases)36391350247913262437
 HR (95% CI)1.00 (Ref.)0.88 (0.82–0.95)0.88 (0.80–0.96)0.80 (0.70–0.90)0.78 (0.68–0.90)0.0140.061
 Cancer-specific mortality (no. of cases)43141381196214822352
 HR (95% CI)1.00 (Ref.)1.01 (0.93–1.11)0.99 (0.90–1.09)1.08 (0.97–1.20)1.12 (1.01–1.24)0.4260.625
 Other cause-specific mortality (no. of cases)32031381245715362399
 HR (95% CI)1.00 (Ref.)0.92 (0.85–0.99)0.87 (0.80–0.95)0.80 (0.68–0.94)0.76 (0.61–0.94)0.001<0.001

CVD, cardiovascular disease; BMI, body mass index.

a

Pooled HRs (95% CIs) were calculated using a random-effects model after study-specific HRs (95% CIs) were obtained using Cox proportional-hazards regression models. Multivariate models were adjusted for baseline age (years; continuous), body mass index (BMI, kg/m2; <18.5, 18.5–19.9, 20.0–24.9, 25.0–29.9 and ≥30.0) and smoking status (never-smoker, former smoker, current with <20 pack-years of smoking and current with ≥20 pack-years of smoking), alcohol intake (non-drinker, 1–155 and ≥156 g/week), educational level (less than secondary, secondary and more than secondary school graduate), energy intake (kcal/day; continuous) and coffee consumption (almost never, <1, 1 to <3, 3 to <5 and ≥5 cups/day).

b

The RERF (Radiation Effects Research Foundation Life Span Study) was excluded from 3 to <5 cups/day and KMCC (Korean Multi-center Cancer Cohort Study) and Seoul Male cohort were excluded from the 3 to <5 and ≥5 cups/day categories because there was no one in these categories due to the way that frequency of intake was assessed in these studies.

c

P for heterogeneity for the highest category was calculated using the Q statistic.

Table 3

Pooled hazard ratios (HRs)a and 95% confidence intervals (CIs) of all-cause and cause-specific mortality according to green tea consumption

Green tea consumption
Almost neverLess than 1 cup/day1 to <3 cups/day3 to <5 cups/day≥5 cups/daybP for trendP for heterogeneityc
Men
 No. of participants56 86741 26555 88736 39935 204
 Total person-years606 553480 888636 679396 522437 025
 All-cause-specific mortality (no. of cases)89246649966669529112
 HR (95% CI)1.00 (Ref.)0.91 (0.87–0.95)0.91 (0.86–0.96)0.91 (0.82–1.02)0.91 (0.82–1.00)0.011<0.001
 CVD-specific mortality (no. of cases)26711980273217992437
 HR (95% CI)1.00 (Ref.)0.90 (0.83–0.98)0.85 (0.80–0.91)0.83 (0.70–0.98)0.79 (0.68–0.91)<0.0010.007
 Cancer-specific mortality (no. of cases)32762459367926723630
 HR (95% CI)1.00 (Ref.)0.94 (0.87–1.02)0.99 (0.90–1.08)1.02 (0.91–1.16)1.07 (0.92–1.25)0.234<0.001
 Other cause-specific mortality (no. of cases)29772210325524813045
 HR (95% CI)1.00 (Ref.)0.88 (0.83–0.94)0.89 (0.84–0.95)0.91 (0.84–0.97)0.83 (0.76–0.90)<0.0010.338
Women
 No. of participants108 50231 47036 59632 50842 494
 Total person-years1 425 367391 362549 387393 930563 773
 All-cause-specific mortality (no. of cases)11 1564112689843447188
 HR (95% CI)1.00 (Ref.)0.94 (0.89–0.98)0.91 (0.86–0.95)0.87 (0.78–0.97)0.86 (0.75–1.00)0.032<0.001
 CVD-specific mortality (no. of cases)36391350247913262437
 HR (95% CI)1.00 (Ref.)0.88 (0.82–0.95)0.88 (0.80–0.96)0.80 (0.70–0.90)0.78 (0.68–0.90)0.0140.061
 Cancer-specific mortality (no. of cases)43141381196214822352
 HR (95% CI)1.00 (Ref.)1.01 (0.93–1.11)0.99 (0.90–1.09)1.08 (0.97–1.20)1.12 (1.01–1.24)0.4260.625
 Other cause-specific mortality (no. of cases)32031381245715362399
 HR (95% CI)1.00 (Ref.)0.92 (0.85–0.99)0.87 (0.80–0.95)0.80 (0.68–0.94)0.76 (0.61–0.94)0.001<0.001
Green tea consumption
Almost neverLess than 1 cup/day1 to <3 cups/day3 to <5 cups/day≥5 cups/daybP for trendP for heterogeneityc
Men
 No. of participants56 86741 26555 88736 39935 204
 Total person-years606 553480 888636 679396 522437 025
 All-cause-specific mortality (no. of cases)89246649966669529112
 HR (95% CI)1.00 (Ref.)0.91 (0.87–0.95)0.91 (0.86–0.96)0.91 (0.82–1.02)0.91 (0.82–1.00)0.011<0.001
 CVD-specific mortality (no. of cases)26711980273217992437
 HR (95% CI)1.00 (Ref.)0.90 (0.83–0.98)0.85 (0.80–0.91)0.83 (0.70–0.98)0.79 (0.68–0.91)<0.0010.007
 Cancer-specific mortality (no. of cases)32762459367926723630
 HR (95% CI)1.00 (Ref.)0.94 (0.87–1.02)0.99 (0.90–1.08)1.02 (0.91–1.16)1.07 (0.92–1.25)0.234<0.001
 Other cause-specific mortality (no. of cases)29772210325524813045
 HR (95% CI)1.00 (Ref.)0.88 (0.83–0.94)0.89 (0.84–0.95)0.91 (0.84–0.97)0.83 (0.76–0.90)<0.0010.338
Women
 No. of participants108 50231 47036 59632 50842 494
 Total person-years1 425 367391 362549 387393 930563 773
 All-cause-specific mortality (no. of cases)11 1564112689843447188
 HR (95% CI)1.00 (Ref.)0.94 (0.89–0.98)0.91 (0.86–0.95)0.87 (0.78–0.97)0.86 (0.75–1.00)0.032<0.001
 CVD-specific mortality (no. of cases)36391350247913262437
 HR (95% CI)1.00 (Ref.)0.88 (0.82–0.95)0.88 (0.80–0.96)0.80 (0.70–0.90)0.78 (0.68–0.90)0.0140.061
 Cancer-specific mortality (no. of cases)43141381196214822352
 HR (95% CI)1.00 (Ref.)1.01 (0.93–1.11)0.99 (0.90–1.09)1.08 (0.97–1.20)1.12 (1.01–1.24)0.4260.625
 Other cause-specific mortality (no. of cases)32031381245715362399
 HR (95% CI)1.00 (Ref.)0.92 (0.85–0.99)0.87 (0.80–0.95)0.80 (0.68–0.94)0.76 (0.61–0.94)0.001<0.001

CVD, cardiovascular disease; BMI, body mass index.

a

Pooled HRs (95% CIs) were calculated using a random-effects model after study-specific HRs (95% CIs) were obtained using Cox proportional-hazards regression models. Multivariate models were adjusted for baseline age (years; continuous), body mass index (BMI, kg/m2; <18.5, 18.5–19.9, 20.0–24.9, 25.0–29.9 and ≥30.0) and smoking status (never-smoker, former smoker, current with <20 pack-years of smoking and current with ≥20 pack-years of smoking), alcohol intake (non-drinker, 1–155 and ≥156 g/week), educational level (less than secondary, secondary and more than secondary school graduate), energy intake (kcal/day; continuous) and coffee consumption (almost never, <1, 1 to <3, 3 to <5 and ≥5 cups/day).

b

The RERF (Radiation Effects Research Foundation Life Span Study) was excluded from 3 to <5 cups/day and KMCC (Korean Multi-center Cancer Cohort Study) and Seoul Male cohort were excluded from the 3 to <5 and ≥5 cups/day categories because there was no one in these categories due to the way that frequency of intake was assessed in these studies.

c

P for heterogeneity for the highest category was calculated using the Q statistic.

Black tea consumption and all-cause and cause-specific mortality

We evaluated the association of black tea with all-cause and cause-specific mortality (10 studies; 205 783 men and 226 963 women). Because of the limited number of participants who drank black tea, we examined this association in the aggregated data set only. The HRs (95% CIs) for those who consumed fewer than three cups per day, compared with black tea non-drinkers, were generally below 1.0 for both men and women (Supplementary Table S7, available as Supplementary data at IJE online).

Discussion

In the Asia Cohort Consortium, we found that coffee consumption was associated with a lower risk of death from all causes, CVD and cancer. Consuming more than five cups of coffee was associated with a 24% and 28% lower risk of death in men and women, respectively. Green tea consumption was also inversely associated with the risk of death from all causes and CVD among both men and women but was associated with an increased risk of cancer death among women. In our study, black tea was less commonly consumed and the associations of black tea with mortality were inconsistent.

Our study found inverse associations between coffee consumption and all-cause and CVD-specific mortality, which corresponded to the results of previous prospective studies. The National Institutes of Health–AARP Diet and Health Study found that coffee consumption was inversely associated with the risk of death from heart disease, respiratory disease, stroke, injuries and accidents, diabetes and infections.13 A pooled analysis of three large US cohort studies found a lower risk of all-cause mortality among those who consumed one to five cups of coffee per day but not among those who consumed more than five cups per day unless the analysis was restricted to never-smokers.12 This pooled analysis also found an inverse association with the CVD-specific mortality. The European Prospective Investigation into Cancer and Nutrition study,16 Multiethnic Cohort study,15 Prostate, Lung, Colorectal and Ovarian Cancer Screening Trial cohort41 and the UK Biobank study42 observed that drinking coffee was associated with a lower all-cause mortality. An increased risk of mortality from ischaemic heart disease with increased coffee consumption was observed in men, whereas a decreased risk of mortality from ischaemic heart disease was observed in women in the Netherlands Cohort Study.43 The JPHC study,44 Japan Collaborative Cohort Study,45 Miyagi Cohort Study46 and SCHS47 reported inverse associations for coffee consumption with all-cause and CVD mortality. The data from three of these cohorts were included in the current analysis.

We found that coffee consumption lowered the risk of cancer-specific mortality. Although several Western cohort studies have found no association12,13,16,41 or a weak inverse association15 for overall cancer death, stronger inverse associations have been reported for coffee consumption and deaths from colorectal cancer and liver cancer.48 Our study warrants further investigation on whether coffee consumption is associated with the incidences of specific types of cancer in Asian populations.

Our study suggests that there is an inverse association between green tea consumption and the risk of death from all causes and CVD. The evidence for an association between green tea consumption and mortality is largely restricted to Japanese and Chinese studies. The Japan Collaborative Cohort Study for Evaluation of Cancer Risk (JACC Study) of 76 979 individuals aged 40–79 years found a 58% lower risk of mortality from coronary heart disease (95% CI 0.17–0.88) among women who drank six or more cups per day of green tea compared with those who did not drink green tea.45 The Prospective Shizuoka Elderly Cohort of 12 251 elderly residents aged 65–84 years in Shizuoka, Japan, reported an inverse association for CVD mortality (HR, 0.24, 95% CI 0.14, 0.40 for those who consumed seven or more cups per day vs those consuming less than one cup per day) but no association for cancer mortality (HR, 0.82, 95% CI 0.45, 1.50 for those consuming seven or more cups per day vs those consuming less than one cup per day).20 The JPHC study indicated that green tea consumption lowered the risk of stroke, as ascertained by computed tomography scanning and/or magnetic resonance imaging.49 The SMHS and SWHS studies found that green tea consumption was inversely associated with the risk of all-cause and CVD-specific mortality.23 A randomized double-blinded trial of 240 Chinese men and women aged ≥18 years demonstrated that a daily capsule of green tea extract (375 mg) decreased the circulating levels of total and low-density lipoprotein (LDL) cholesterol and increased the levels of high-density lipoprotein cholesterol.50 Another Japanese double-blinded randomized trial of 240 men and women aged 25–55 years demonstrated that supplementation with green tea beverages containing 583 mg of catechins (intervention group) vs one containing only 96 mg of catechins (control group) led to greater reductions in body fat, systolic blood pressure and LDL cholesterol.51

In the present study, for green tea consumption, we observed no association among men and an increased risk of death from cancer among women. Most previous cohort studies found that green tea consumption was not associated with cancer-specific mortality.20–22,47 Previous cohort studies did not report associations between green tea consumption and cancers of the stomach,52 breast,53 prostate54 or colon.55 Although the reason why we observed an increased risk at five or more cups per day of green tea is not clear, the potential mechanism may include the thermal injury to the mucosa and inflammation from drinking tea that has a high temperature.56,57 For black tea consumption, a few cohort studies that assessed the association between the overall tea consumption and mortality in Western countries, where black tea is more common, showed an increased risk of total mortality58 or showed no association,59,60 whereas our study suggested a weak inverse, albeit inconsistent, association with mortality.

Several potential biologically active compounds found in coffee, including chlorogenic acid, caffeine, diterpenes, trigonelline, melanoidins, potassium and magnesium, may explain the potential health benefits of coffee.4 These compounds have been shown to have antioxidant, anti-inflammatory and antidiabetic activities in experimental and clinical trials. Coffee polyphenols have been shown to scavenge free radicals61,62 and modulate the postprandial glycaemic response in mouse models63 and in humans.64 Antioxidants in green tea may similarly be responsible for the potential health benefits related to green tea consumption. Green tea is rich in catechins, the major form of which is epigallocatechin-3-gallate, which increases the resistance to oxidative damage. Green tea catechins have been shown to reduce LDL oxidation65 and adipose insulin resistance.66

Our study has several limitations. First, because decaffeinated coffee was not commonly consumed in Asia at the time at which these cohorts were initiated, the caffeine type was generally not assessed. Thus, we were unable to distinguish decaffeinated coffee from total coffee consumption. Nevertheless, cohorts that have evaluated the associations of decaffeinated and caffeinated coffee with mortality have consistently observed similar results for both types of coffee,13,15,41, suggesting that compounds other than caffeine are important. Second, coffee consumption was only assessed at baseline; therefore, changes in consumption habits over time were not considered. Third, we cannot rule out the possibility that residual confounding factors, including smoking and alcohol consumption, may be present. A positive association for green tea consumption among women could be partly explained by residual confounding factors. In our study, coffee drinkers were more likely to smoke and green tea drinkers were more likely to drink alcohol. Thus, it is possible that there is a residual confounding effect by these risk factors, or other unmeasured or poorly measured risk factors may have affected the HR estimates. Additionally, the possibility of reverse causation could be another limiting factor. However, the exclusion of participants who died during the first 3 years of the follow-up period did not change the results. Additionally, we found similar associations with the presence or absence of the baseline co-morbidities. Additionally, we did not have detailed information on additives (e.g. sugar or cream) or brewing methods for either coffee or tea and therefore we were not able to examine the effects of additives or types of brewing on mortality. For the black tea analysis, we were not able to include all of the studies or able to examine the association above two cups per day because of only a few cases in that category. The strengths of this study include its use of multiple large Asian cohorts to evaluate the associations of coffee and tea consumption with mortality. Additionally, the adjustment for potential confounding factors and robust associations in subgroup and sensitivity analyses suggests that our findings may not be fully explained by confounding factors or reverse causation. The contribution of our study to current knowledge lies primarily in the convincing finding that coffee and green tea consumption are associated with a lower risk of mortality across a diverse sample of Asian populations. Although a few Asian studies have published analyses of coffee or green tea, other Asian studies have not analysed or published their results, owing, in part, to a limited statistical power. Our pooled analysis of 10 or more studies across diverse populations in four countries enabled us to examine a wide range of tea and coffee intake, and to explore the potential effect modifications by age, smoking and other important risk factors for mortality.

In conclusion, drinking as little as one cup of coffee per day was associated with a lower risk of death during the follow-up period among Asian populations. Green tea consumption was associated with a lower risk of CVD death, but the association between black tea consumption and mortality was unclear. These findings warrant further observational and intervention studies in Asian populations to explore the effects of the specific types of coffee and tea.

Supplementary data

Supplementary data are available at IJE online.

Ethics approval

The present study was reviewed and approved by the Institutional Review Boards of their corresponding institutions.

Funding

The cohorts participating in the pooled analysis were supported by the following grants: Japan Public Health Center-based Prospective Study (JPHC Study), National Cancer Center Research and Development Fund (since 2011) and a Grant-in-Aid for Cancer Research from the Ministry of Health, Labour and Welfare of Japan (from 1989 to 2010); Miyagi Cohort Study, National Cancer Center Research and Development Fund; 3-prefecture Miyagi, National Cancer Center Research and Development Fund; Ohsaki National Health Insurance Cohort Study, National Cancer Center Research and Development Fund; Radiation Effects Research Foundation Life Span Study, The Japanese Ministry of Health, Labour and Welfare and the US Department of Energy; 3-prefecture AICHI, National Cancer Center Research and Development Fund (30-A-15); Shanghai Men's Health Study, the US National Cancer Institute (UM1CA173640); Shanghai Women’s Health Study, the US National Cancer Institute (R37 CA70867 and UM1CA182910); and Singapore Chinese Health Study (SCHS), the US National Cancer Institute (R01 CA144034 and UM1 CA182876). This work was also supported, in part, by the Intramural Research Program of the National Institutes of Health.

Data availability

Investigators are granted access to the Asia Cohort Consortium data upon reasonable request and with the approval of ACC members and the Institutional Review Board. We cannot publicly provide individual data due to participant privacy based on the ethical guidelines in Japan. Additionally, the informed consent that the participating studies obtained did not include a provision for publicly sharing data.

Acknowledgements

The authors are grateful to all the participants and staff of the Asia Cohort Consortium Center at the National Cancer Center Japan.

Author contributions

J.E.L. and R.S. had the idea for the study and designed the study. S.S., J.E.L., E.L. and R.S. wrote the first draft. X.-O.S., S.T., N.S., I.T., S.K., Y.S., Y.T., A.S., K.O., I.O., H.I., M-H.S., Y.-O.A., S.K.P., A.S., Y.-B.X., H.C., W.-P.K., J.-M.Y., K.-Y.Y., K.S.C., H.A., W.Z., M.I., D.K. and K.M. collected the data. S.K.A., M.S.R., E.S., M.R.I. and M.I. coordinated the ACC data collection and analysis. S.S. analysed the data. All authors contributed to data interpretation and approved the final manuscript.

Conflict of interest

None declared.

References

1

Food and Agriculture Organization of the United Nations. FAO Statistical Databases (FAOSTAT). Food Supply—Crops Primary Equivalent (Food Balance).

2019
. http://www.fao.org/faostat/en/#data (26 July 2019, date last accessed).

2

The FAO Intergovernmental Group on Tea. Twenty-Third Session of the Intergovernmental Group on Tea. Current Market Situation and Medium Term Outlook.

2018
. http://www.fao.org/economic/est/est-commodities/tea/tea-meetings/tea23/en (26 July 2019, date last accessed).

3

Yang
CS
,
Hong
J.
Prevention of chronic diseases by tea: possible mechanisms and human relevance
.
Annu Rev Nutr
2013
;
33
:
161
81
.

4

Ludwig
IA
,
Clifford
MN
,
Lean
ME
,
Ashihara
H
,
Crozier
A.
Coffee: biochemistry and potential impact on health
.
Food Funct
2014
;
5
:
1695
717
.

5

van Dam
RM
,
Hu
FB.
Coffee consumption and risk of type 2 diabetes: a systematic review
.
JAMA
2005
;
294
:
97
104
.

6

Ding
M
,
Bhupathiraju
SN
,
Satija
A
,
van Dam
RM
,
Hu
FB.
Long-term coffee consumption and risk of cardiovascular disease: a systematic review and a dose-response meta-analysis of prospective cohort studies
.
Circulation
2014
;
129
:
643
59
.

7

Larsson
SC
,
Orsini
N.
Coffee consumption and risk of stroke: a dose-response meta-analysis of prospective studies
.
Am J Epidemiol
2011
;
174
:
993
1001
.

8

Hernan
MA
,
Takkouche
B
,
Caamano-Isorna
F
,
Gestal-Otero
JJ.
A meta-analysis of coffee drinking, cigarette smoking, and the risk of Parkinson's disease
.
Ann Neurol
2002
;
52
:
276
84
.

9

Lu
Y
,
Zhai
L
,
Zeng
J
et al.
Coffee consumption and prostate cancer risk: an updated meta-analysis
.
Cancer Causes Control
2014
;
25
:
591
604
.

10

Je
Y
,
Giovannucci
E.
Coffee consumption and risk of endometrial cancer: findings from a large up-to-date meta-analysis
.
Int J Cancer
2012
;
131
:
1700
10
.

11

Larsson
SC
,
Wolk
A.
Coffee consumption and risk of liver cancer: a meta-analysis
.
Gastroenterology
2007
;
132
:
1740
45
.

12

Ding
M
,
Satija
A
,
Bhupathiraju
SN
et al.
Association of coffee consumption with total and cause-specific mortality in 3 large prospective cohorts
.
Circulation
2015
;
132
:
2305
15
.

13

Freedman
ND
,
Park
Y
,
Abnet
CC
,
Hollenbeck
AR
,
Sinha
R.
Association of coffee drinking with total and cause-specific mortality
.
N Engl J Med
2012
;
366
:
1891
904
.

14

Tsujimoto
T
,
Kajio
H
,
Sugiyama
T.
Association between caffeine intake and all-cause and cause-specific mortality: a population-based prospective cohort study
.
Mayo Clin Proc
2017
;
92
:
1190
202
.

15

Park
SY
,
Freedman
ND
,
Haiman
CA
,
Le Marchand
L
,
Wilkens
LR
,
Setiawan
VW.
Association of coffee consumption with total and cause-specific mortality among nonwhite populations
.
Ann Intern Med
2017
;
167
:
228
35
.

16

Gunter
MJ
,
Murphy
N
,
Cross
AJ
et al.
Coffee drinking and mortality in 10 European countries: a multinational cohort study
.
Ann Intern Med
2017
;
167
:
236
47
.

17

Crippa
A
,
Discacciati
A
,
Larsson
SC
,
Wolk
A
,
Orsini
N.
Coffee consumption and mortality from all causes, cardiovascular disease, and cancer: a dose-response meta-analysis
.
Am J Epidemiol
2014
;
180
:
763
75
.

18

Sasazuki
S
,
Tamakoshi
A
,
Matsuo
K
et al. ; Research Group for the Development and Evaluation of Cancer Prevention Strategies in Japan.
Green tea consumption and gastric cancer risk: an evaluation based on a systematic review of epidemiologic evidence among the Japanese population
.
Jpn J Clin Oncol
2012
;
42
:
335
46
.

19

Liu
J
,
Liu
S
,
Zhou
H
et al.
Association of green tea consumption with mortality from all-cause, cardiovascular disease and cancer in a Chinese cohort of 165,000 adult men
.
Eur J Epidemiol
2016
;
31
:
853
65
.

20

Suzuki
E
,
Yorifuji
T
,
Takao
S
et al.
Green tea consumption and mortality among Japanese elderly people: the prospective Shizuoka elderly cohort
.
Ann Epidemiol
2009
;
19
:
732
39
.

21

Kuriyama
S
,
Shimazu
T
,
Ohmori
K
et al.
Green tea consumption and mortality due to cardiovascular disease, cancer, and all causes in Japan: the Ohsaki study
.
JAMA
2006
;
296
:
1255
65
.

22

Saito
E
,
Inoue
M
,
Sawada
N
et al.
Association of green tea consumption with mortality due to all causes and major causes of death in a Japanese population: the Japan Public Health Center-based Prospective Study (JPHC Study)
.
Ann Epidemiol
2015
;
25
:
512
18.e3
.

23

Zhao
LG
,
Li
HL
,
Sun
JW
et al.
Green tea consumption and cause-specific mortality: results from two prospective cohort studies in China
.
J Epidemiol
2017
;
27
:
36
41
.

24

Zheng
W
,
McLerran
DF
,
Rolland
B
et al.
Association between body-mass index and risk of death in more than 1 million Asians
.
N Engl J Med
2011
;
364
:
719
29
.

26

Watanabe
S
,
Tsugane
S
,
Sobue
T
,
Konishi
M
,
Baba
S.
Study design and organization of the JPHC study: Japan Public Health Center-based Prospective Study on cancer and cardiovascular diseases
.
J Epidemiol
2001
;
11
:
S3
7
.

27

Tsuji
I
,
Nishino
Y
,
Ohkubo
T
et al.
A prospective cohort study on National Health Insurance beneficiaries in Ohsaki, Miyagi Prefecture, Japan: study design, profiles of the subjects and medical cost during the first year
.
J Epidemiol
1998
;
8
:
258
63
.

28

Ozasa
K
,
Shimizu
Y
,
Suyama
A
et al.
Studies of the mortality of atomic bomb survivors, Report 14, 1950-2003: an overview of cancer and noncancer diseases
.
Radiat Res
2012
;
177
:
229
43
.

29

Sado
J
,
Kitamura
T
,
Kitamura
Y
et al. ; Three-Prefecture Cohort Study Group.
Rationale, design, and profile of the Three-Prefecture Cohort in Japan: a 15-year follow-up
.
J Epidemiol
2017
;
27
:
193
99
.

30

Rhee
CW
,
Kim
JY
,
Park
BJ
,
Li
ZM
,
Ahn
Y-O.
Impact of individual and combined health behaviors on all causes of premature mortality among middle aged men in Korea: the Seoul Male Cohort Study
.
J Prev Med Public Health
2012
;
45
:
14
20
.

31

Yoo
KY
,
Shin
HR
,
Chang
SH
et al.
Korean Multi-center Cancer Cohort Study including a Biological Materials Bank (KMCC-I)
.
Asian Pac J Cancer Prev
2002
;
3
:
85
92
.

32

Shu
XO
,
Li
H
,
Yang
G
et al.
Cohort profile: The Shanghai Men's Health Study
.
Int J Epidemiol
2015
;
44
:
810
18
.

33

Zheng
W
,
Chow
WH
,
Yang
G
et al.
The Shanghai Women's Health Study: rationale, study design, and baseline characteristics
.
Am J Epidemiol
2005
;
162
:
1123
31
.

34

Hankin
JH
,
Stram
DO
,
Arakawa
K
et al.
Singapore Chinese Health Study: development, validation, and calibration of the quantitative food frequency questionnaire
.
Nutr Cancer
2001
;
39
:
187
95
.

35

Shu
XO
,
Yang
G
,
Jin
F
et al.
Validity and reproducibility of the food frequency questionnaire used in the Shanghai Women's Health Study
.
Eur J Clin Nutr
2004
;
58
:
17
23
.

36

Kim
MK
,
Lee
SS
,
Ahn
YO.
Reproducibility and validity of a self-administered semiquantitative food frequency questionnaire
.
Korean J Community Nutr
1996
;
1
:
376
94
.

37

Tsubono
Y
,
Kobayashi
M
,
Sasaki
S
,
Tsugane
S
; JPHC.
Validity and reproducibility of a self-administered food frequency questionnaire used in the baseline survey of the JPHC Study Cohort I
.
J Epidemiol
2003
;
13
:
S125
33
.

38

Ogawa
K
,
Tsubono
Y
,
Nishino
Y
et al.
Validation of a food-frequency questionnaire for cohort studies in rural Japan
.
Public Health Nutr
2003
;
6
:
147
57
.

39

Sauvaget
C
,
Allen
N
,
Hayashi
M
,
Spencer
E
,
Nagano
J.
Validation of a food frequency questionnaire in the Hiroshima/Nagasaki Life Span Study
.
J Epidemiol
2002
;
12
:
394
401
.

40

Venkateshan
SP
,
Sidhu
S
,
Malhotra
S
,
Pandhi
P.
Meta-analysis in clinical trials
.
Control Clin Trials
2009
;
83
:
1
88
.

41

Loftfield
E
,
Freedman
ND
,
Graubard
BI.
Association of coffee consumption with overall and cause-specific mortality in a large US prospective cohort study
.
Am J Epidemiol
2015
;
182
:
1010
22
.

42

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
:
1086
97
.

43

Leurs
LJ
,
Schouten
LJ
,
Goldbohm
RA
,
van den Brandt
PA.
Total fluid and specific beverage intake and mortality due to IHD and stroke in the Netherlands Cohort Study
.
Br J Nutr
2010
;
104
:
1212
21
.

44

Saito
E
,
Inoue
M
,
Sawada
N
et al.
Association of coffee intake with total and cause-specific mortality in a Japanese population: the Japan Public Health Center-based Prospective Study
.
Am J Clin Nutr
2015
;
101
:
1029
37
.

45

Mineharu
Y
,
Koizumi
A
,
Wada
Y
et al. ; JACC study Group.
Coffee, green tea, black tea and oolong tea consumption and risk of mortality from cardiovascular disease in Japanese men and women
.
J Epidemiol Community Health
2011
;
65
:
230
40
.

46

Sugiyama
K
,
Kuriyama
S
,
Akhter
M
et al.
Coffee consumption and mortality due to all causes, cardiovascular disease, and cancer in Japanese women
.
J Nutr
2010
;
140
:
1007
13
.

47

Odegaard
AO
,
Koh
WP
,
Yuan
JM
,
Pereira
MA.
Beverage habits and mortality in Chinese adults
.
J Nutr
2015
;
145
:
595
604
.

48

Poole
R
,
Kennedy
OJ
,
Roderick
P
,
Fallowfield
JA
,
Hayes
PC
,
Parkes
J.
Coffee consumption and health: umbrella review of meta-analyses of multiple health outcomes
.
BMJ
2017
;
359
:
j5024
.

49

Kokubo
Y
,
Iso
H
,
Saito
I
et al.
The impact of green tea and coffee consumption on the reduced risk of stroke incidence in Japanese population: the Japan public health center-based study cohort
.
Stroke
2013
;
44
:
1369
74
.

50

Maron
DJ
,
Lu
GP
,
Cai
NS
et al.
Cholesterol-lowering effect of a theaflavin-enriched green tea extract: a randomized controlled trial
.
Arch Intern Med
2003
;
163
:
1448
53
.

51

Nagao
T
,
Hase
T
,
Tokimitsu
I.
A green tea extract high in catechins reduces body fat and cardiovascular risks in humans
.
Obesity (Silver Spring)
2007
;
15
:
1473
83
.

52

Myung
SK
,
Bae
WK
,
Oh
SM
et al.
Green tea consumption and risk of stomach cancer: a meta-analysis of epidemiologic studies
.
Int J Cancer
2009
;
124
:
670
77
.

53

Sun
CL
,
Yuan
JM
,
Koh
WP
,
Yu
MC.
Green tea, black tea and breast cancer risk: a meta-analysis of epidemiological studies
.
Carcinogenesis
2006
;
27
:
1310
15
.

54

Zheng
J
,
Yang
B
,
Huang
T
,
Yu
Y
,
Yang
J
,
Li
D.
Green tea and black tea consumption and prostate cancer risk: an exploratory meta-analysis of observational studies
.
Nutr Cancer
2011
;
63
:
663
72
.

55

Wang
ZH
,
Gao
QY
,
Fang
JY.
Green tea and incidence of colorectal cancer: evidence from prospective cohort studies
.
Nutr Cancer
2012
;
64
:
1143
52
.

56

Yu
C
,
Tang
H
,
Guo
Y
, China Kadoorie Biobank Collaborative Group et al.
Hot tea consumption and its interactions with alcohol and tobacco use on the risk for esophageal cancer: a population-based cohort study
.
Ann Intern Med
2018
;
168
:
489
97
.

57

Mao
XQ
,
Jia
XF
,
Zhou
G
et al.
Green tea drinking habits and gastric cancer in southwest China
.
Asian Pac J Cancer Prev
2011
;
12
:
2179
82
.

58

Hertog
MG
,
Sweetnam
PM
,
Fehily
AM
,
Elwood
PC
,
Kromhout
D.
Antioxidant flavonols and ischemic heart disease in a Welsh population of men: the Caerphilly Study
.
Am J Clin Nutr
1997
;
65
:
1489
94
.

59

de Koning Gans
JM
,
Uiterwaal
CS
,
van der Schouw
YT
et al.
Tea and coffee consumption and cardiovascular morbidity and mortality
.
Arterioscler Thromb Vasc Biol
2010
;
30
:
1665
71
.

60

Andersen
LF
,
Jacobs
DR
Jr
,
Carlsen
MH
,
Blomhoff
R.
Consumption of coffee is associated with reduced risk of death attributed to inflammatory and cardiovascular diseases in the Iowa Women's Health Study
.
Am J Clin Nutr
2006
;
83
:
1039
46
.

61

Chen
JH
,
Ho
C-T.
Antioxidant activities of caffeic acid and its related hydroxycinnamic acid compounds
.
J Agric Food Chem
1997
;
45
:
2374
78
.

62

Nardini
M
,
D'Aquino
M
,
Tomassi
G
,
Gentili
V
,
Di Felice
M
,
Scaccini
C.
Inhibition of human low-density lipoprotein oxidation by caffeic acid and other hydroxycinnamic acid derivatives
.
Free Radic Biol Med
1995
;
19
:
541
52
.

63

Murase
T
,
Yokoi
Y
,
Misawa
K
et al.
Coffee polyphenols modulate whole-body substrate oxidation and suppress postprandial hyperglycaemia, hyperinsulinaemia and hyperlipidaemia
.
Br J Nutr
2012
;
107
:
1757
65
.

64

Johnston
KL
,
Clifford
MN
,
Morgan
LM.
Coffee acutely modifies gastrointestinal hormone secretion and glucose tolerance in humans: glycemic effects of chlorogenic acid and caffeine
.
Am J Clin Nutr
2003
;
78
:
728
33
.

65

Yang
TT
,
Koo
MW.
Inhibitory effect of Chinese green tea on endothelial cell-induced LDL oxidation
.
Atherosclerosis
2000
;
148
:
67
73
.

66

Yan
J
,
Zhao
Y
,
Suo
S
,
Liu
Y
,
Zhao
B.
Green tea catechins ameliorate adipose insulin resistance by improving oxidative stress
.
Free Radic Biol Med
2012
;
52
:
1648
57
.

This article is published and distributed under the terms of the Oxford University Press, Standard Journals Publication Model (https://academic.oup.com/journals/pages/open_access/funder_policies/chorus/standard_publication_model)

Supplementary data