Research Paper Volume 13, Issue 14 pp 18564—18585

Chocolate consumption and all-cause and cause-specific mortality in a US population: a post hoc analysis of the PLCO cancer screening trial

Guo-Chao Zhong1, *, , Tian-Yang Hu2, *, , Peng-Fei Yang3, , Yang Peng4, , Jing-Jing Wu5, , Wei-Ping Sun6, , Long Cheng7, , Chun-Rui Wang8, ,

  • 1 Department of Hepatobiliary Surgery, The Second Affiliated Hospital of Chongqing Medical University, Chongqing, China
  • 2 Department of Cardiology, The Second Affiliated Hospital of Chongqing Medical University, Chongqing, China
  • 3 Department of Nephrology, The Second Affiliated Hospital of Chongqing Medical University, Chongqing, China
  • 4 Department of Geriatrics, The Fifth People’s Hospital of Chengdu, Chengdu, China
  • 5 Department of Nutrition and Food Hygiene, School of Public Health and Management, Chongqing Medical University, Chongqing, China
  • 6 Department of Gastrointestinal Surgery, The Second Affiliated Hospital of Chongqing Medical University, Chongqing, China
  • 7 Department of Breast and Thyroid Surgery, The Second Affiliated Hospital of Chongqing Medical University, Chongqing, China
  • 8 Department of Infectious Diseases, Institute for Viral Hepatitis, The Key Laboratory of Molecular Biology for Infectious Diseases, Chinese Ministry of Education, The Second Affiliated Hospital of Chongqing Medical University, Chongqing, China
* Equal contribution

Received: February 22, 2021       Accepted: July 6, 2021       Published: July 29, 2021      

https://doi.org/10.18632/aging.203302
How to Cite

Copyright: © 2021 Zhong et al. This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Abstract

Few studies with mixed results have examined the association between chocolate consumption and mortality. We aimed to examine this association in a US population. A population-based cohort of 91891 participants aged 55 to 74 years was identified. Chocolate consumption was assessed via a food frequency questionnaire. Cox regression was used to estimate risk estimates. After an average follow-up of 13.5 years, 19586 all-cause deaths were documented. Compared with no regular chocolate consumption, the maximally adjusted hazard ratios of all-cause mortality were 0.89 [95% confidence interval (CI) 0.84–0.94], 0.84 (95% CI 0.79–0.90), 0.86 (95% CI 0.81–0.93), and 0.87 (95% CI 0.82–0.93) for >0–0.5 servings/week, >0.5–1 serving/week, >1–2 servings/week, and >2 servings/week, respectively (Ptrend = 0.009). A somewhat stronger inverse association was observed for mortality from cardiovascular disease and Alzheimer’s disease. A nonlinear dose–response pattern was found for all-cause and cardiovascular mortality (all Pnonlinearity < 0.01), with the lowest risk observed at chocolate consumption of 0.7 servings/week and 0.6 servings/week, respectively. The favorable associations with all-cause and cardiovascular mortality were found to be more pronounced in never smokers than in current or former smokers (all Pinteraction < 0.05). In conclusion, chocolate consumption confers reduced risks of mortality from all causes, cardiovascular disease, and Alzheimer’s disease in this US population.

Introduction

Because of its exquisite taste and flavor, chocolate is a widely consumed confectionery product, especially in Western countries [1]. In 2018/2019, it is estimated that around 7.7 million tons of chocolate confectionery were consumed globally [2]. Mounting studies have investigated health effects of chocolate [3]. Several meta-analyses found that chocolate consumption was inversely associated with risks of stroke, coronary heart disease, and heart failure [46], whereas other studies found that chocolate consumption was not significantly associated with the risk of overall invasive cancer [7] or was positively associated with risks of breast and colorectal cancers [7, 8].

However, evidence on the association of chocolate consumption with mortality is sparse and somewhat mixed. Specifically, higher chocolate consumption was found to be significantly associated with lower mortality from all causes [9] and cardiovascular disease (CVD) [9, 10], whereas no significant association was observed for chocolate consumption and coronary heart disease mortality [11]. These inconclusive results may be due to the differences in study population, study outcome, and/or study method. A study of 470 elderly men observed that higher cocoa intake conferred lower all-cause and cardiovascular mortality [12]. However, findings from this study may not be applicable to chocolate, as cocoa is only a natural product of cocoa beans while chocolate is a solid food product mainly including cocoa, cocoa butter, and sugar [13]; furthermore, this study has a small sample size, and could be subject to small study effects (i.e., small studies are conducted with less methodological rigor and tend to report larger risk estimates compared with large studies) [14].

Currently, evidence on the association between chocolate consumption and health outcomes is weak [15]. Given the need for more data on this topic, we conducted a post hoc analysis of the Prostate, Lung, Colorectal, and Ovarian (PLCO) Cancer Screening Trial to examine the hypothesis that chocolate consumption confers reduced risks of all-cause and cause-specific mortality in a US population.

Results

Participant characteristics

In the whole study population, the mean (standard deviation) age was 65.3 (5.7) years; the proportion of males was 46.3% (Table 1); and the average chocolate consumption (standard deviation) was 1.5 (3.2) servings/week. Compared with chocolate non-consumers, those consuming >2 servings/week were less likely to have a history of diabetes or hypertension, had higher energy intake from diet but lower Healthy Eating Index-2015. Moreover, participants consuming >2 vs. 0 servings/week of chocolate had higher consumption of red meat, processed meat, dairy, and coffee but lower consumption of fruit and vegetable, and higher intakes of saturated fatty acids, added sugar, sodium, magnesium, potassium, and calcium.

Table 1. Baseline characteristics of study population according to chocolate consumptiona.

CharacteristicsEnergy-adjusted chocolate consumption (servings/week)P
0>0–0.5>0.5–1>1–2>2
No. of participants536840760148211090720035<0.001
Age (years)65.8 ± 5.765.3 ± 5.665.3 ± 5.765.0 ± 5.765.2 ± 5.8<0.001
Male2907 (54.2)18217 (44.7)7090 (47.8)5097 (46.7)9232 (46.1)<0.001
Ethnic group
 Non-Hispanic white4372 (81.4)36249 (88.9)13702 (92.4)10236 (93.8)18989 (94.8)<0.001
 Non-Hispanic black433 (8.1)1619 (4.0)381 (2.6)209 (1.9)339 (1.7)
 Hispanic140 (2.6)655 (1.6)208 (1.4)129 (1.2)228 (1.1)
 Othersb423 (7.9)2237 (5.5)530 (3.6)333 (3.1)479 (2.4)
Married or living as married4000 (74.5)31975 (78.4)11695 (78.9)8695 (79.7)15568 (77.7)<0.001
Body mass index (kg/m2)27.1 ± 5.127.0 ± 4.827.4 ± 4.827.2 ± 4.727.3 ± 4.8<0.001
Educational degree
 College below3512 (65.4)25735 (63.1)9479 (64.0)6790 (62.3)12649 (63.1)<0.001
 College graduate884 (16.5)7227 (17.7)2652 (17.9)1948 (17.9)3621 (18.1)
 Postgraduate972 (18.1)7798 (19.1)2690 (18.1)2169 (19.9)3765 (18.8)
Alcohol consumption (g/day)10.6 ± 29.610.3 ± 26.69.8 ± 26.58.7 ± 21.58.1 ± 21.6<0.001
Smoking status
 Current442 (8.2)3472 (8.5)1383 (9.3)982 (9.0)2199 (11.0)<0.001
 Former2418 (45.0)17274 (42.4)6261 (42.2)4417 (40.5)7983 (39.8)
 Never2508 (46.7)20014 (49.1)7177 (48.4)5508 (50.5)9853 (49.2)
Physical activity (min/week)c130.2 ± 133.9124.9 ± 123.8122.2 ± 121.3121.0 ± 120.8117.0 ± 120.5<0.001
Energy intake from diet (kcal/day)1641.8 ± 747.41567.7 ± 658.41748.0 ± 716.51840.3 ± 729.82040.6 ± 789.7<0.001
Healthy Eating Index-201568.8 ± 9.968.2 ± 9.566.5 ± 9.465.8 ± 9.363.1 ± 9.5<0.001
History of diabetes908 (16.9)2651 (6.5)643 (4.3)397 (3.6)687 (3.4)<0.001
History of hypertension1810 (33.7)12628 (31.0)4502 (30.4)3115 (28.6)5860 (29.2)<0.001
Family history of cancer2828 (52.9)22711 (55.9)8410 (56.9)6137 (56.4)11439 (57.2)<0.001
Aspirin user2274 (42.6)17605 (43.4)6494 (44.0)4842 (44.6)8917 (44.7)0.004
Food consumption
 Red meat (g/day)54.1 ± 54.953.2 ± 46.164.9 ± 51.769.0 ± 55.672.7 ± 58.4<0.001
 Processed meat (g/day)16.8 ± 23.114.8 ± 16.917.7 ± 18.718.4 ± 19.019.4 ± 19.9<0.001
 Fruit (g/day)312.2 ± 276.0272.7 ± 219.7266.8 ± 200.5273.0 ± 204.1270.5 ± 209.5<0.001
 Vegetable (g/day)319.1 ± 238.2278.3 ± 187.1280.3 ± 178.1285.3 ± 173.7285.4 ± 175.7<0.001
 Whole grain (servings/day)1.3 ± 1.01.1 ± 0.81.2 ± 0.81.2 ± 0.81.2 ± 0.8<0.001
 Dairy (cups/day)1.2 ± 1.21.3 ± 1.11.4 ± 1.11.5 ± 1.11.6 ± 1.2<0.001
 Coffee (g/day)762.0 ± 792.5843.7 ± 769.7869.5 ± 795.7854.4 ± 804.1846.1 ± 829.8<0.001
 Tea (g/day)285.7 ± 525.1259.2 ± 458.5253.8 ± 446.9262.9 ± 457.9270.6 ± 486.8<0.001
Nutrient intake
 Dietary fiber (g/day)19.4 ± 10.517.1 ± 8.217.8 ± 8.318.5 ± 8.219.4 ± 8.4<0.001
 Saturated fatty acids (g/day)16.6 ± 10.816.8 ± 9.820.2 ± 11.121.9 ± 11.726.2 ± 13.2<0.001
 Added sugar (tsp/day)9.4 ± 8.810.0 ± 7.512.2 ± 8.413.6 ± 8.617.8 ± 10.5<0.001
 Sodium (mg/day)2677.4 ± 1294.62496.7 ± 1103.22767.4 ± 1193.22894.8 ± 1232.13064.3 ± 1292.8<0.001
 Magnesium (mg/day)328.7 ± 146.3304.6 ± 122.1321.2 ± 125.6332.0 ± 125.9351.6 ± 131.7<0.001
 Potassium (mg/day)3287.9 ± 1418.23095.1 ± 1210.13252.6 ± 1233.53348.8 ± 1243.73476.2 ± 1286.0<0.001
 Calcium (mg/day)722.2 ± 419.2695.5 ± 386.7753.7 ± 400.0786.8 ± 404.8843.1 ± 425.6<0.001
aValues are mean ± standard deviation or counts (percentage) as indicated.
b“Others” refers to Asian, Pacific Islander, or American Indian.
cTotal time of moderate to vigorous physical activities per week.

Chocolate consumption and all-cause and cause-specific mortality

During a mean (standard deviation) follow-up length of 13.5 (3.3) years (1238513.2 person-years), a total of 19586 all-cause deaths were observed, of which 5490 (28.0%) and 6175 (31.5%) were classified as deaths from CVD and cancer, respectively. Crude mortality rates of all causes, CVD, and cancer were 158.14, 44.33, and 49.86 per 10000 person-years, respectively. Compared with no regular chocolate consumption, the maximally adjusted HRs of all-cause mortality were 0.89 (95% CI 0.84–0.94), 0.84 (95% CI 0.79–0.90), 0.86 (95% CI 0.81–0.93), and 0.87 (95% CI 0.82–0.93) for >0–0.5 servings/week, >0.5–1 serving/week, >1–2 servings/week, and >2 servings/week, respectively (Ptrend = 0.009) (Table 2). Interestingly, chocolate consumption was also found to be inversely associated with cardiovascular mortality, with somewhat stronger magnitude of the association observed, while no significant association was observed for cancer mortality. In addition, we explored the association of chocolate consumption with the risk of mortality from cerebrovascular, respiratory, and Alzheimer’s diseases. In the fully adjusted model, chocolate consumption was found to be inversely associated with the risk of mortality from Alzheimer’s disease (HR for >2 vs. 0 servings/week 0.69; 95% CI 0.49–0.99) but not cerebrovascular (HR for >2 vs. 0 servings/week 0.97; 95% CI 0.73–1.28) and respiratory (HR for >2 vs. 0 servings/week 0.90; 95% CI 0.74–1.10) diseases (Table 2).

Table 2. Association of chocolate consumption with all-cause and cause-specific mortalitya.

Causes of mortalityEnergy-adjusted chocolate consumption (servings/week)Ptrend
0>0–0.5>0.5–1>1–2>2
All causes
No. of deaths14558617304921874278
Death rateb207.00156.44152.07147.94159.00
Model 1c1.00 (reference)0.81 (0.77–0.86)0.77 (0.73–0.83)0.78 (0.73–0.83)0.82 (0.77–0.87)0.002
Model 2d1.00 (reference)0.89 (0.84–0.94)0.84 (0.79–0.90)0.87 (0.81–0.93)0.88 (0.83–0.94)0.026
Model 3e1.00 (reference)0.89 (0.84–0.94)0.84 (0.79–0.90)0.86 (0.81–0.93)0.87 (0.82–0.93)0.009
Cardiovascular disease
No. of deaths48324088565841159
Death rateb68.7243.7242.6939.5043.08
Model 1c1.00 (reference)0.70 (0.63–0.77)0.67 (0.60–0.75)0.64 (0.57–0.73)0.69 (0.62–0.76)0.001
Model 2d1.00 (reference)0.79 (0.72–0.88)0.76 (0.68–0.85)0.77 (0.68–0.87)0.79 (0.70–0.88)0.027
Model 3e1.00 (reference)0.79 (0.72–0.88)0.76 (0.68–0.86)0.77 (0.68–0.87)0.78 (0.70–0.88)0.020
Cancer
No. of deaths38826599937381397
Death rateb55.2048.2749.5249.9251.92
Model 1c1.00 (reference)0.93 (0.84–1.04)0.94 (0.84–1.06)0.97 (0.86–1.10)1.00 (0.89–1.12)0.105
Model 2d1.00 (reference)0.95 (0.86–1.06)0.95 (0.84–1.07)0.99 (0.87–1.12)0.98 (0.87–1.10)0.595
Model 3e1.00 (reference)0.95 (0.86–1.06)0.95 (0.84–1.07)0.99 (0.87–1.13)0.99 (0.88–1.11)0.431
Cerebrovascular disease
No. of deaths74518177133224
Death rateb10.539.408.839.008.33
Model 1c1.00 (reference)0.94 (0.74–1.20)0.88 (0.67–1.15)0.93 (0.70–1.23)0.84 (0.64–1.09)0.145
Model 2d1.00 (reference)1.04 (0.81–1.33)0.98 (0.74–1.30)1.10 (0.82–1.47)0.98 (0.74–1.29)0.503
Model 3e1.00 (reference)1.05 (0.82–1.35)0.99 (0.75–1.31)1.10 (0.82–1.48)0.97 (0.73–1.28)0.359
Respiratory disease
No. of deaths138899309224468
Death rateb19.6316.3215.4115.1517.39
Model 1c1.00 (reference)0.89 (0.74–1.06)0.82 (0.67–1.00)0.83 (0.67–1.03)0.93 (0.77–1.12)0.284
Model 2d1.00 (reference)0.97 (0.80–1.16)0.87 (0.71–1.06)0.91 (0.73–1.13)0.93 (0.77–1.14)0.904
Model 3e1.00 (reference)0.95 (0.79–1.15)0.85 (0.69–1.04)0.89 (0.72–1.11)0.90 (0.74–1.10)0.658
Alzheimer’s disease
No. of deaths452678256134
Death rateb6.404.854.093.794.98
Model 1c1.00 (reference)0.74 (0.54–1.02)0.61 (0.43–0.88)0.59 (0.40–0.87)0.75 (0.54–1.06)0.088
Model 2d1.00 (reference)0.76 (0.55–1.05)0.62 (0.42–0.89)0.56 (0.38–0.84)0.71 (0.50–1.01)0.031
Model 3e1.00 (reference)0.78 (0.56–1.08)0.63 (0.43–0.91)0.57 (0.38–0.86)0.69 (0.49–0.99)0.041
aValues are hazard ratios (95% confidence intervals).
bCrude death rate per 10000 person-years.
cAdjusted for age (years), sex (male, female), and ethnicity (non-Hispanic white, non-Hispanic black, Hispanic, others).
dAdjusted for model 1 plus educational level (college below, college graduate, postgraduate), marital status (married or living as married, widowed, divorced, separated, never married), study center (10 categories), history of hypertension (yes, no), history of diabetes (yes, no), aspirin use (yes, no), hormone use status (current, former, never) for women, smoking status (current, former, never), alcohol consumption (g/day), body mass index (kg/m2), physical activity (min/week), and energy intake from diet (kcal/day). For all-cause and cancer mortality, the model 2 was further adjusted for family history of cancer (yes, no).
eAdjusted for model 2 plus consumption of red meat (g/day), processed meat (g/day), fruit (g/day), vegetable (g/day), whole grain (servings/day), dairy (cups/day), coffee (g/day), and tea (g/day).

Subgroup analyses

A significant interaction between chocolate consumption and smoking status was observed for mortality from all causes (Pinteraction = 0.006) and CVD (Pinteraction = 0.045) but not cancer (Pinteraction = 0.504) (Table 3). Specifically, the favorable association between chocolate consumption and mortality from all causes and CVD was found to be more pronounced in never smokers than in current or former smokers. No significant interaction was observed for remaining stratification factors, namely age, sex, trial group, BMI, alcohol consumption, milk consumption, and history of hypertension. In addition, we performed a subgroup analysis by smoking status for the association of chocolate consumption with the risk of mortality from Alzheimer’s disease; the results showed that there was no significant difference in the risk of mortality from Alzheimer’s disease between current or past and never smokers (Pinteraction = 0.812).

Table 3. Subgroup analyses on the association of chocolate consumption with mortality from all causes, cardiovascular disease, and cancera.

Subgroup variableEnergy-adjusted chocolate consumption (servings/week)Pinteraction
0>0–0.5>0.5–1>1–2>2
All-cause mortality
Age (years)
 ≥601.00 (reference)0.88 (0.83–0.93)0.83 (0.77–0.88)0.85 (0.80–0.92)0.87 (0.81–0.93)0.369
 <601.00 (reference)1.03 (0.83–1.28)1.03 (0.81–1.31)1.01 (0.79–1.30)0.96 (0.75–1.21)
Sex
 Male1.00 (reference)0.90 (0.83–0.98)0.86 (0.79–0.95)0.87 (0.79–0.96)0.90 (0.82–0.98)0.879
 Female1.00 (reference)0.87 (0.80–0.94)0.82 (0.75–0.89)0.85 (0.77–0.94)0.84 (0.77–0.92)
Trial group
 Screening1.00 (reference)1.00 (0.94–1.07)0.89 (0.83–0.95)0.95 (0.90–1.02)0.93 (0.87–1.00)0.683
 Control1.00 (reference)0.97 (0.91–1.04)0.88 (0.82–0.94)0.92 (0.87–0.98)0.94 (0.88–1.00)
Body mass index (kg/m2)
 ≥251.00 (reference)0.88 (0.82–0.95)0.86 (0.79–0.93)0.86 (0.79–0.94)0.88 (0.81–0.95)0.755
 <251.00 (reference)0.90 (0.82–0.99)0.82 (0.74–0.92)0.89 (0.79–1.00)0.88 (0.79–0.98)
Alcohol consumption (g/day)b
 No, light or moderate1.00 (reference)0.91 (0.85–0.96)0.85 (0.79–0.91)0.87 (0.81–0.94)0.88 (0.82–0.94)0.306
 Heavy1.00 (reference)0.79 (0.68–0.93)0.80 (0.67–0.95)0.86 (0.71–1.03)0.83 (0.69–0.98)
Smoking status
 Current or former1.00 (reference)0.91 (0.85–0.98)0.92 (0.84–0.99)0.93 (0.85–1.01)0.94 (0.87–1.02)0.006
 Never1.00 (reference)0.87 (0.79–0.95)0.76 (0.69–0.85)0.80 (0.72–0.89)0.81 (0.74–0.90)
Milk consumption (servings/day)
 ≥ median1.00 (reference)0.91 (0.83–0.99)0.84 (0.77–0.93)0.88 (0.79–0.97)0.85 (0.78–0.93)0.118
 < median1.00 (reference)0.87 (0.80–0.94)0.84 (0.77–0.92)0.86 (0.78–0.94)0.91 (0.83–1.00)
History of hypertension
 Yes1.00 (reference)0.92 (0.85–1.01)0.86 (0.77–0.95)0.94 (0.84–1.05)0.93 (0.84–1.03)0.207
 No1.00 (reference)0.86 (0.79–0.92)0.83 (0.76–0.90)0.82 (0.75–0.89)0.83 (0.77–0.90)
Cardiovascular mortality
Age (years)
 ≥601.00 (reference)0.80 (0.72–0.89)0.78 (0.69–0.88)0.77 (0.68–0.88)0.79 (0.70–0.89)0.639
 <601.00 (reference)0.73 (0.50–1.05)0.59 (0.38–0.92)0.68 (0.43–1.09)0.76 (0.50–1.15)
Sex
 Male1.00 (reference)0.81 (0.70–0.94)0.84 (0.71–0.99)0.80 (0.67–0.96)0.83 (0.70–0.98)0.345
 Female1.00 (reference)0.77 (0.67–0.89)0.68 (0.58–0.80)0.73 (0.61–0.87)0.74 (0.63–0.86)
Trial group
 Screening1.00 (reference)0.98 (0.87–1.11)0.88 (0.77–1.00)0.87 (0.77–0.98)0.90 (0.80–1.02)0.627
 Control1.00 (reference)0.91 (0.81–1.02)0.83 (0.73–0.94)0.91 (0.80–1.02)0.89 (0.78–1.01)
Body mass index (kg/m2)
 ≥251.00 (reference)0.81 (0.72–0.91)0.77 (0.67–0.88)0.77 (0.66–0.89)0.81 (0.70–0.93)0.891
 <251.00 (reference)0.78 (0.65–0.93)0.79 (0.64–0.97)0.79 (0.63–0.99)0.76 (0.62–0.93)
Alcohol consumption (g/day)b
 No, light or moderate1.00 (reference)0.91 (0.85–0.96)0.85 (0.79–0.91)0.87 (0.81–0.93)0.88 (0.82–0.94)0.279
 Heavy1.00 (reference)0.80 (0.68–0.93)0.80 (0.67–0.95)0.86 (0.71–1.04)0.83 (0.70–0.99)
Smoking status
 Current or former1.00 (reference)0.79 (0.69–0.90)0.82 (0.70–0.95)0.83 (0.70–0.97)0.86 (0.74–1.00)0.045
 Never1.00 (reference)0.81 (0.70–0.95)0.71 (0.60–0.85)0.71 (0.58–0.86)0.72 (0.60–0.86)
Milk consumption (servings/day)
 ≥ median1.00 (reference)0.88 (0.76–1.03)0.80 (0.68–0.95)0.80 (0.67–0.96)0.84 (0.71–0.99)0.227
 < median1.00 (reference)0.72 (0.63–0.83)0.74 (0.63–0.87)0.75 (0.63–0.89)0.75 (0.64–0.88)
History of hypertension
 Yes1.00 (reference)0.85 (0.73–0.98)0.81 (0.68–0.96)0.85 (0.70–1.02)0.84 (0.71–1.00)0.636
 No1.00 (reference)0.75 (0.65–0.86)0.73 (0.62–0.85)0.70 (0.60–0.83)0.73 (0.63–0.85)
Cancer mortality
Age (years)
 ≥601.00 (reference)0.94 (0.84–1.06)0.92 (0.81–1.04)0.97 (0.85– 1.11)0.98 (0.87–1.11)0.348
 <601.00 (reference)1.05 (0.73–1.50)1.24 (0.84–1.82)1.19 (0.80– 1.77)1.05 (0.71–1.54)
Sex
 Male1.00 (reference)1.07 (0.91–1.26)1.05 (0.88–1.25)1.14 (0.95–1.37)1.10 (0.92–1.31)0.339
 Female1.00 (reference)0.86 (0.74–1.00)0.86 (0.73–1.01)0.87 (0.73–1.04)0.90 (0.76–1.05)
Trial group
 Screening1.00 (reference)0.92 (0.81–1.03)0.93 (0.83–1.06)0.99 (0.89–1.11)1.06 (0.94–1.20)0.265
 Control1.00 (reference)0.95 (0.84–1.06)0.99 (0.88–1.11)0.92 (0.82–1.03)0.95 (0.84–1.08)
Body mass index (kg/m2)
 ≥251.00 (reference)1.01 (0.88–1.16)1.01 (0.87–1.17)1.06 (0.91–1.24)1.01 (0.87–1.17)0.305
 <251.00 (reference)0.86 (0.72–1.02)0.84 (0.69–1.04)0.87 (0.70–1.08)0.95 (0.78–1.16)
Alcohol consumption (g/day)a
 No, light or moderate1.00 (reference)1.02 (0.90–1.15)1.00 (0.88–1.14)1.05 (0.92–1.21)1.05 (0.92–1.20)0.108
 Heavy1.00 (reference)0.69 (0.54–0.88)0.72 (0.54–0.96)0.73 (0.54–1.00)0.71 (0.53–0.95)
Smoking status
 Current or former1.00 (reference)0.99 (0.87–1.13)0.98 (0.84–1.13)1.08 (0.92–1.26)1.04 (0.90–1.20)0.504
 Never1.00 (reference)0.91 (0.76–1.09)0.94 (0.76–1.15)0.88 (0.71–1.10)0.96 (0.78–1.17)
Milk consumption (servings/day)
 ≥ median1.00 (reference)1.02 (0.87–1.20)1.01 (0.84–1.20)1.06 (0.88–1.27)0.98 (0.82–1.17)0.117
 < median1.00 (reference)0.90 (0.78–1.04)0.90 (0.76–1.06)0.94 (0.79–1.12)1.03 (0.88–1.22)
History of hypertension
 Yes1.00 (reference)1.09 (0.90–1.31)1.02 (0.82–1.26)1.21 (0.97–1.51)1.13 (0.92–1.39)0.204
 No1.00 (reference)0.89 (0.78–1.02)0.91 (0.78–1.050.90 (0.77–1.05)0.92 (0.80–1.06)
aValues are hazard ratios (95% confidence intervals). Hazard ratios were adjusted for age (years), sex (male, female), ethnicity (non-Hispanic white, non-Hispanic black, Hispanic, others), educational level (college below, college graduate, postgraduate), marital status (married or living as married, widowed, divorced, separated, never married), study center (10 categories), history of hypertension (yes, no), history of diabetes (yes, no), aspirin use (yes, no), hormone use status (current, former, never) for women, smoking status (current, former, never), alcohol consumption (g/day), body mass index (kg/m2), physical activity (min/week), energy intake from diet (kcal/day), and consumption of red meat (g/day), processed meat (g/day), fruit (g/day), vegetable (g/day), whole grain (servings/day), dairy (cups/day), coffee (g/day), and tea (g/day). For all-cause and cancer mortality, the model was further adjusted for family history of cancer (yes, no). In subgroup analyses stratified by sex, smoking status, and history of hypertension, hazard ratios were not adjusted for the stratification factor.
bFor men, light, moderate, and heavy alcohol consumption referred to ≤ 6 g/day, > 6–28 g/day, and > 28 g/day, respectively; for women, light, moderate, and heavy alcohol consumption referred to ≤ 6 g/day, > 6–14 g/day, and > 14 g/day, respectively.

Dose–response analyses

In the entire study population, chocolate consumption was found to be associated with lower risks of death from all causes (Pnonlinearity <0.001) and CVD (Pnonlinearity <0.001) in a nonlinear dose–response pattern, with the lowest risk observed at chocolate consumption of 0.7 servings/week and 0.6 servings/week, respectively (Figure 1). Given the aforementioned significant interaction between chocolate consumption and smoking status, we further conducted smoking status-specific dose–response analyses (Supplementary Figures 1 and 2 for never and current or former smokers, respectively). The nonlinear dose–response associations of chocolate consumption with risks of death from all causes and CVD were only found among never smoker, with the lowest risk observed at chocolate consumption of 0.6 servings/week and 0.7 servings/week, respectively. No significant dose–response association was found for cancer mortality in the entire study population as well as in never or current or former smokers.

Nonlinear dose–response analyses on energy-adjusted chocolate consumption and mortality from all causes, cardiovascular disease, and cancer in the whole study population. The reference level was set at 0 servings/week. Hazard ratio was adjusted for age, sex, ethnicity, educational level, marital status, study center, history of hypertension, history of diabetes, aspirin use, hormone use status for women, smoking status, alcohol consumption, body mass index, physical activity, energy intake from diet, and consumption of red meat, processed meat, fruit, vegetable, whole grain, dairy, coffee, and tea. For all-cause and cancer mortality, the hazard ratio was further adjusted for family history of cancer. The red solid line represents the fitted nonlinear trend, and the purple short-dash line represents corresponding 95% confidence interval.

Figure 1. Nonlinear dose–response analyses on energy-adjusted chocolate consumption and mortality from all causes, cardiovascular disease, and cancer in the whole study population. The reference level was set at 0 servings/week. Hazard ratio was adjusted for age, sex, ethnicity, educational level, marital status, study center, history of hypertension, history of diabetes, aspirin use, hormone use status for women, smoking status, alcohol consumption, body mass index, physical activity, energy intake from diet, and consumption of red meat, processed meat, fruit, vegetable, whole grain, dairy, coffee, and tea. For all-cause and cancer mortality, the hazard ratio was further adjusted for family history of cancer. The red solid line represents the fitted nonlinear trend, and the purple short-dash line represents corresponding 95% confidence interval.

Sensitivity analyses

The initial associations of chocolate consumption with risks of death from all causes, CVD, and cancer did not change materially in a large range of sensitivity analyses (Supplementary Table 1).

Discussion

In this post hoc analysis, chocolate consumption was found to be inversely associated with all-cause and cardiovascular mortality in a nonlinear dose–response manner, with the lowest risk observed at chocolate consumption of 0.7 servings/week and 0.6 servings/week, respectively, whereas no significant association was found for cancer mortality. Subgroup analyses further found that the inverse association of chocolate consumption with all-cause and cardiovascular mortality was more pronounced in never smokers. In addition, chocolate consumption was found to be inversely associated with the risk of death from Alzheimer’s disease but not cerebrovascular and respiratory diseases.

To date, only a small number of studies have determined the potential association between chocolate consumption and mortality [9, 16]. A cohort study of 1169 acute myocardial infarction survivors [16] observed that higher chocolate consumption conferred lower mortality from CVD but not all causes. Subsequently, a prospective study of 84709 postmenopausal women revealed an inverse association between chocolate consumption and mortality from all causes, CVD, and Alzheimer’s disease [9]. However, whether these observations can be extended to the general population is unknown. Based on prospective data from the PLCO Cancer Screening Trial, we found that chocolate consumption was inversely associated with risks of death from all causes, cardiovascular disease, and Alzheimer’s disease, which is consistent with the results from the study by Sun et al. [9]. Our study presents the following important information. First, our study provided data on the association of chocolate consumption with cancer mortality. Second, our study identified the shape and the nadir of dose–response curve for the chocolate-mortality association, which are essential for recommending optimal chocolate consumption to the public. Third, our study found that smoking behavior could significantly modify the observed inverse association of chocolate consumption with all-cause and cardiovascular mortality.

Previously, a prospective study of healthy men and women found that increasing chocolate consumption was monotonically associated with decreased cardiovascular mortality (Plinearity = 0.011) [10]. In contrast, our dose–response analysis found that the association of chocolate consumption with cardiovascular mortality was nonlinear, with the lowest risk observed at 0.6 servings/week, indicating that individuals consuming low amounts of chocolate could benefit most. This finding can be explained by the fact that chocolate is rich in sugar, carbohydrate, and fat [17]. High intakes of these macronutrients are associated with increased mortality [1820], which possibly weakens or even negates the favorable association between chocolate and mortality. In fact, several prospective studies also found that low to moderate but not high chocolate consumption was favorably associated with risks of atrial fibrillation [21] and heart failure [22, 23].

In our study, the information on chocolate types, such as milk chocolate and dark chocolate, was not collected by the DHQ. However, it should be noted that different types of chocolate have different nutritional ingredients [17]. For example, dark chocolate has a significantly higher proportion of cocoa solids than milk chocolate [24], resulting in that dark chocolate contains more polyphenols, a class of bioactive substances mainly found in vegetables and fruits [17]. In addition, milk is found to be capable of reducing the absorption of antioxidants from chocolate [25]. Thus, it is biologically possible that dark chocolate has stronger physiological effects in human body than milk chocolate. Indeed, a crossover study of 16 young men found that compared with participants consuming milk chocolate, those consuming dark chocolate felt more satiated and had lower desire to eat something sweet and lower energy intake [26]. Dark chocolate was also found to have a higher ability to increase total antioxidant capacity of human plasma than milk chocolate [25]. These facts raise a possibility that dark chocolate has more benefits in improving health outcomes than milk chocolate. Unfortunately, the unavailability of chocolate type data in our study hindered us to determine whether there is a difference in reducing mortality between dark chocolate and milk chocolate. Nonetheless, in the US, milk chocolate is more popular than dark chocolate at the time of the study [27], indicating that the observed associations with all-cause and cardiovascular mortality could be largely driven by milk chocolate. In addition, after a thorough literature review, we failed to find an observational study that has examined the potential impacts of chocolate types on health outcomes [4, 6, 15], reminding us that this seems to be a common limitation of studies in this filed. Hence, future studies should pay much attention to the assessment of chocolate types, and investigate the associations with health outcomes by specific chocolate types.

Interestingly, we observed that the inverse association of chocolate consumption with all-cause and cardiovascular mortality was more pronounced in never smokers than in current or former smokers. Similarly, Mink and colleagues also observed that the favorable association between dietary intake of flavonol, a bioactive ingredient in chocolate, and all-cause mortality was more pronounced in never smokers than in ever smokers [11]. These observations suggest that smoking behavior possibly has interactions with chocolate consumption in biological pathways. It is well established that smoking leads to increased levels of systemic inflammation and oxidative stress in human body [28, 29], resulting in that smokers have a higher risk of death than never smokers [30]. Thus, it is highly possible that the favorable association of chocolate consumption with all-cause and cardiovascular mortality has been weakened or even overcome by the above-mentioned harmful effects of smoking behavior. In addition, as mentioned above, milk has been reported to be able to interfere with the absorption of antioxidants from chocolate [25]; hence, we performed a subgroup analysis to determine whether the association of chocolate consumption with mortality could be modified by this factor. However, we did not observe an expected interaction after stratifying for milk consumption. In fact, Mostofsky and colleagues also found that the association between chocolate consumption and the risk of heart failure could not be modified by milk consumption [23]. The exact reason behind this phenomenon is unclear, but may be due to the fact that chocolate and milk are not consumed simultaneously in the daily life.

The inverse association of chocolate consumption with all-cause and cardiovascular mortality is biologically plausible. Several randomized controlled trials have found that dark chocolate consumption efficiently decreases blood pressure [3134], increases insulin sensitivity [32, 34], improves cholesterol profile [33, 35, 36] and endothelial function [37], and inhibits platelet aggregation [38]; these cardiometabolic effects of chocolate have been suggested to mediate its favorable association with CVD risk [39, 40]. Chocolate is a significant source of flavonoids and is especially rich in proanthocyanidins, catechins, and epicatachins [41]. In vitro and animal studies have indicated that flavonoids have antioxidant and anti-inflammatory effects and favorable impacts on glucose-insulin homeostasis [41]. Also, dietary flavonoid intake has been identified to be associated with lower all-cause and cardiovascular mortality [42]. Therefore, the above-mentioned cardiometabolic effects of chocolate may be explained by the biological effects of flavonoids in nature, at least partly. In addition, we observed an inverse association between chocolate consumption and the risk of death from Alzheimer’s disease. In fact, previous observational studies had found that chocolate consumption was associated with a decreased risk of cognitive decline [43, 44]. The observed association of chocolate consumption with the risk of death from Alzheimer’s disease could be attributable to the neuroprotective effects of cocoa [45]. An in vitro Alzheimer’s disease model study suggested that cocoa exerted its neuroprotective effects through activating BDNF signaling pathway [46]. In addition, chocolate contains a certain amount of caffeine, which was found to have protective effects on mortality from all causes and CVD [47] and have some beneficial effects on the development of Alzheimer’s disease [48]; thus, the roles of caffeine in the observed inverse association should not be ignored. More studies are needed to determine which ingredients in chocolate are actually responsible for the inverse association we observed.

Our study has several limitations that should be acknowledged and discussed. First, chocolate consumption was only assessed at baseline in the PLCO Cancer Screening Trial by the DHQ. The assessment of dietary exposure at one time point possibly results in non-differential bias, considering that dietary habits could change over time. Nonetheless, a classic assumption in nutritional epidemiology is that an adult’s dietary habits would not change dramatically during several years. Moreover, it has been suggested that the method only using baseline diet generally yields a weaker association than that using the cumulative averages [49]. Second, in the present study, we used death certificates to obtain the underlying causes of death. However, it should be highlighted that the cause of death from death certificates is misclassified in some conditions, especially deaths due to respiratory or digestive diseases [50]. Hence, our results on chocolate consumption and cause-specific mortality might be subject to misclassification bias. Third, although we had controlled for a wide range of possible confounders, our results could be susceptible to residual confounding because of the observational design of this study. Fourth, given that CVD mortality was a primary outcome in our study, we excluded subjects with a history of acute myocardial infarction or stroke at baseline; however, as the PLCO Cancer Screening Trial did not collect data on coronary artery disease, thus we could not exclude subjects with this disease. Hence, our results might be subject to reverse causation and confounded by other CVDs (e.g., coronary artery disease). Fifth, chocolate consumption was assessed by a self-administrated food frequency questionnaire in our study. Thus, the assessment of chocolate consumption might be subject to measurement errors. Nonetheless, this bias is nondifferential, given that it is not expected to be associated with the future risk of death, and thus tends to attenuate the association of our interest. Moreover, chocolate consumption might have been underestimated in our study as it might be added to other foods. Finally, our findings derived from a US population aged 55 to 74 years with a relatively high chocolate consumption, and thus might not be extended to other age groups or other populations with lower chocolate consumption.

In conclusion, in this US population aged 55 to 74 years, chocolate consumption is associated with lower risks of death from all causes, cardiovascular disease, and Alzheimer’s disease. Due to the observational design of our study, these findings do not imply causation and needed to be further confirmed in other populations and settings. If confirmed, eating chocolate may be a good choice for improving longevity. Future studies are warranted to clarify the potential influence of chocolate types on the observed association.

Materials and Methods

This study was reported following the Strengthening the Reporting of Observational Studies in Epidemiology statement [51].

Study population

A total of 76682 men and 78215 women aged 55 to 74 years were enrolled to the PLCO Cancer Screening Trial between November 1993 and September 2001 in ten screening centers across the US (Washington, Pittsburgh, Honolulu, Denver, Marshfield, Minneapolis, Birmingham, Salt Lake City, Detroit, and St Louis). The PLCO Cancer Screening Trial was a multicenter randomized controlled trial designed to determine whether screening for prostate, lung, colorectal, and ovarian cancers could reduce mortality from these cancers, and its study design and implementation were described elsewhere [52]. The PLCO Cancer Screening Trial was approved by the Institutional Review Boards of the US National Cancer Institute and each study center, and written informed consent was obtained from all individuals. The study was conducted in accordance with the Declaration of Helsinki.

The following individuals were further excluded from our study: (1) individuals with an invalid diet history questionnaire (DHQ), which refers to missing the date of DHQ completion, missing ≥8 DHQ items, death date prior to DHQ completion date, or the presence of extreme values of energy intake (i.e., top 1% or bottom 1%) (n = 41444); (2) individuals diagnosed with any cancer before baseline questionnaire or DHQ completion (n = 9684); (3) individuals having a history of acute myocardial infarction or stroke at baseline (n = 9932); and (4) individuals failing to return the baseline questionnaire (n = 1946). After exclusions, a total of 91891 individuals were included (Figure 2). Importantly, no marked differences were observed in sociodemographic characteristics and medical histories between included and excluded populations (all standardized differences <0.1; Supplementary Table 2), suggesting a small possibility of selection bias due to the exclusion of a large number of individuals in our study. For all eligible individuals, follow-up length was computed from the date of DHQ completion to the date of death, loss to follow-up, the date of mortality, or the end of follow-up (i.e., December 31, 2015), whichever occurred first (Figure 3).

The study flow chart of identifying eligible participants. The percentage in the figure was calculated as the number of participants excluded by each exclusion criterion divided by the total number of participants in the PLCO Cancer Screening Trial (i.e., 154897). PLCO, Prostate, Lung, Colorectal, and Ovarian; DHQ, diet history questionnaire.

Figure 2. The study flow chart of identifying eligible participants. The percentage in the figure was calculated as the number of participants excluded by each exclusion criterion divided by the total number of participants in the PLCO Cancer Screening Trial (i.e., 154897). PLCO, Prostate, Lung, Colorectal, and Ovarian; DHQ, diet history questionnaire.

The timeline and follow-up scheme of the present study. DHQ, diet history questionnaire.

Figure 3. The timeline and follow-up scheme of the present study. DHQ, diet history questionnaire.

Assessment of chocolate consumption

Chocolate consumption was evaluated using the above-mentioned DHQ in the present study. Participants were asked to answer “How often did you eat chocolate during the past year?”. There were 11 predefined answers: never, 1–6 times/year, 7–11 times/year, 1 time/month, 2–3 times/month, 1 time/week, 2 times/week, 3–4 times/week, 5–6 times/week, 1 time/day, and ≥2 times/day. Chocolate consumers were further asked to answer “Each time you ate chocolate, how much did you usually eat?”. The predefined answers were “Less than 1 average bar or less than 28.35 g”, “1 average bar or 28.35 to 56.70 g”, and “More than 1 average bar or more than 56.70 g”. Chocolate consumption (servings/week, 1 serving size of chocolate was defined as 28.35g chocolate [53]) was estimated by multiplying the frequency of chocolate consumption by portion size. Notably, chocolate consumption involved in all analyses was adjusted for energy intake with the residual method for removing extraneous variation of chocolate consumption resulting from energy intake [54].

Outcome assessment

In the PLCO Cancer Screening Trial, a mailed annual study update form was used to confirm participants’ vital status. For participants failing to return this form, they were repeatedly contacted through telephone or e-mail. Additionally, the information on death was adjudicated by periodic linkage to the National Death Index for increasing its completeness. The International Classification of Diseases, ninth Revision (ICD-9), which was frequently used when the trial began, was used to define the underlying causes of mortality obtained from death certificates: CVD (ICD codes 390–459) and cancer (ICD codes 140–209).

Assessment of covariates

A self-administrated baseline questionnaire was used to collect the following data: sex, marital status, ethnicity, height, body weight, educational degree, family history of cancer, history of diabetes, history of hypertension, aspirin use, and smoking status. Body mass index (BMI) was computed as body weight (kg) divided by height squared (m2). A DHQ (version 1.0, National Cancer Institute, 2007) was used to assess the remaining covariates: age at DHQ completion, alcohol consumption, energy intake from diet, food consumption, and nutrient intake. The DHQ was a self-administered 137-item food frequency questionnaire, which was developed to evaluate the frequency and portion size of food consumption and nutrient intake over the past 12 months. Notably, the validity of the DHQ had been confirmed against four 24-h dietary recalls in The Eating at America's Table Study [55]. Daily food consumption was estimated by multiplying the food frequency by portion size; daily nutrient intake was estimated based on the Nutrition Data Systems for Research [56] and the USDA's 1994-96 Continuing Survey of Food Intakes by Individuals [57]. Healthy Eating Index-2015, an indicator of diet quality, with higher scores indicating higher diet quality, was computed as previously described [58]. Physical activity level was approximated based on the frequency and duration of moderate and strenuous activities, which were collected with a self-administrated supplemental questionnaire. The supplemental questionnaire was introduced in 2006 in the PLCO Cancer Screening Trial; in our study, the mean (standard deviation) time from trial entry to the completion of this questionnaire was 3309 (680) days (Figure 3).

Statistical analysis

As there were several variables with missing data (see Supplementary Table 3), for reducing selection bias and increasing statistical power, we assumed that these data were missing at random, and then used multiple imputation with chained equations to impute missing data (the number of imputations = 25) [59]. All variables involved in statistical analyses were used to produce the imputed data sets.

Cox proportional hazards regression was performed to estimate hazard ratios (HRs) and 95% confidence intervals (CIs) for the association of chocolate consumption with all-cause and cause-specific mortality, with person-year as time variable. No evidence for the violation of proportional hazards assumption was detected using Schoenfeld residuals. Chocolate consumption was classified into five categories based on its distribution in our study population and published articles on chocolate consumption and health outcomes [6, 23, 60] (0, >0–0.5, >0.5–1, >1–2, >2 servings/week), with 0 servings/week as the reference group. A P for quadratic trend was obtained by assigning the median value of each category to each individual in the category, and then tested the statistical significance of its squared value in regression analyses. Covariate included in multivariable analyses was selected based on the change-in-estimate approach [61] and the existing literature. Specifically, model 1 was adjusted for age, sex, and ethnicity; model 2 was further adjusted for educational degree, marital status, study center, history of hypertension, history of diabetes, aspirin use, hormone use status for women, smoking status, alcohol consumption, BMI, physical activity, and energy intake from diet; and model 3 was further adjusted for consumption of red meat, processed meat, fruit, vegetable, whole grain, dairy, coffee, and tea. Of note, for all-cause and cancer mortality, model 2 was additionally adjusted for family history of cancer.

Prespecified subgroup analyses were performed after stratifying for age (≥ 65 vs. < 65 years), sex (male vs. female), trial group (screening group vs. control group), BMI (≥ 25 vs. < 25 kg/m2) [62, 63], alcohol consumption (no, light, or moderate vs. heavy), smoking status (current or former vs. never), milk consumption (≥ median vs. < median), and history of hypertension (yes vs. no). In our study, for men, light, moderate, and heavy alcohol consumption referred to ≤ 6 g/day, > 6–28 g/day, and > 28 g/day, respectively; for women, light, moderate, and heavy alcohol consumption referred to ≤ 6 g/day, > 6–14 g/day, and > 14 g/day, respectively [64]. A Pinteraction was obtained via a likelihood ratio test, in which models with and without interaction terms were compared, before conducting the above-mentioned subgroup analyses to avoid spurious subgroup differences.

Restricted cubic spline functions [65] with four knots located at the 5th, 35th, 65th, and 95th percentiles were used to describe the shape and the nadir of the dose–response curves for the association of chocolate consumption with all-cause and cause-specific mortality, with the reference level set at 0 servings/week. Of note, the choice of the number and location of knots was based on the recommendations by Harrell [66] and Akaike's information criterion [67]. Importantly, for minimizing the potential influence of extreme values, individuals with the top 2.5% of chocolate consumption (i.e., ≥ 9.1 servings/week) were excluded from the dose–response analysis. A Pnonlinearity was obtained by testing the null hypothesis that regression coefficients of the second and third splines were equal to 0 [65].

A series of sensitivity analyses were conducted to assess the stability of our results: (1) excluded individuals with extreme values of calorie intake, which were defined as <800 or >4000 kcal/day for men and <500 or >3500 kcal/day for women [68]; (2) excluded deaths observed within the first five years of follow-up to evaluate the potential influence of reverse causation; (3) adjusted for propensity score on crude model (all covariates in model 3 were used to compute propensity score for each participant with a logistic regression); (4) additionally adjusted for Healthy Eating Index-2015 to evaluate the potential influence of diet quality; (5) additionally adjusted for intakes of sodium, added sugars, and saturated fatty acids to indirectly evaluate the potential influence of adherence to a healthy eating pattern, as the Dietary Guidelines for Americans recommended that a healthy eating pattern should limit the intakes of these nutrients [69].

Continuous variables were shown as mean (standard deviation), and categorical variables were shown as counts (percentage). The ANOVA test and the χ2 test were employed to compare the differences of continuous and categorical variables, respectively, across the categories of chocolate consumption. Statistical analyses were performed using Stata version 12.0 (StataCorp, College Station, TX), and the corresponding results were considered statistically significant when a two-tailed P value was less than 0.05.

Author Contributions

GCZ, TYH, LC, and CRW developed the hypothesis, study design and concept, and all authors made useful suggestions for study design and concept. GCZ acquired original data. GCZ and TYH drafted study protocol. PFY, YP, JJW, WPS, and CRW were responsible for statistical analyses, and all author interpreted the results together. GCZ drafted the initial manuscript, and all author made critical revisions for the initial manuscript. CRW were responsible for manuscript submission. All authors were involved in answering reviewers’ questions. GCZ and CRW had full access to original data. LC and CRW acted as guarantors for the integrity of the data and the accuracy of statistical analysis. All authors approved the final version of the article, including the authorship list. The corresponding author attests that all listed authors meet authorship criteria and that no others meeting the criteria have been omitted.

Acknowledgments

The authors thank Prof. Yong Zhang (School of Public Health and Management, Chongqing Medical University, Chongqing, China) for his assistance in the statistical analysis. The authors sincerely appreciate the National Cancer Institute for access to NCI’s data collected by the Prostate, Lung, Colorectal and Ovarian Cancer Screening Trial. The statements contained herein are solely those of the authors and do not represent or imply concurrence or endorsement by NCI.

Conflicts of Interest

The authors declare no conflicts of interest related to this study.

Funding

This study did not receive any external funding.

References

  • 1. Santander Muñoz M, Rodríguez Cortina J, Vaillant FE, Escobar Parra S. An overview of the physical and biochemical transformation of cocoa seeds to beans and to chocolate: Flavor formation. Crit Rev Food Sci Nutr. 2020; 60:1593–613. https://doi.org/10.1080/10408398.2019.1581726 [PubMed]
  • 2. International Cocoa Organization. Consumption of retail chocolate confectionery worldwide 2012/13-2018/19. Hamburg, Germany: Statista Research Department; 2015 (Accessed April 21, 2020 https://www.statista.com/statistics/238849/global-chocolate-consumption/).
  • 3. Morze J, Schwedhelm C, Bencic A, Hoffmann G, Boeing H, Przybylowicz K, Schwingshackl L. Chocolate and risk of chronic disease: a systematic review and dose-response meta-analysis. Eur J Nutr. 2020; 59:389–97. https://doi.org/10.1007/s00394-019-01914-9 [PubMed]
  • 4. Yuan S, Li X, Jin Y, Lu J. Chocolate Consumption and Risk of Coronary Heart Disease, Stroke, and Diabetes: A Meta-Analysis of Prospective Studies. Nutrients. 2017; 9:688. https://doi.org/10.3390/nu9070688 [PubMed]
  • 5. Gong F, Yao S, Wan J, Gan X. Chocolate Consumption and Risk of Heart Failure: A Meta-Analysis of Prospective Studies. Nutrients. 2017; 9:402. https://doi.org/10.3390/nu9040402 [PubMed]
  • 6. Ren Y, Liu Y, Sun XZ, Wang BY, Zhao Y, Liu DC, Zhang DD, Liu XJ, Zhang RY, Sun HH, Liu FY, Chen X, Cheng C, et al. Chocolate consumption and risk of cardiovascular diseases: a meta-analysis of prospective studies. Heart. 2019; 105:49–55. https://doi.org/10.1136/heartjnl-2018-313131 [PubMed]
  • 7. Greenberg JA, Neuhouser ML, Tinker LF, Lane DS, Paskett ED, Van Horn LV, Wassertheil-Smoller S, Shikany JM, Qi L, Sealy-Jefferson S, Manson JE. Chocolate Candy and Incident Invasive Cancer Risk in the Women's Health Initiative: An Observational Prospective Analysis. J Acad Nutr Diet. 2021; 121:314–26.e4. https://doi.org/10.1016/j.jand.2020.06.014 [PubMed]
  • 8. Tavani A, Giordano L, Gallus S, Talamini R, Franceschi S, Giacosa A, Montella M, La Vecchia C. Consumption of sweet foods and breast cancer risk in Italy. Ann Oncol. 2006; 17:341–45. https://doi.org/10.1093/annonc/mdj051 [PubMed]
  • 9. Sun Y, Liu B, Snetselaar L, Wallace R, Shadyab A, Chen GC, Shikany J, Manson J, Bao W. Chocolate Consumption in Relation to All-Cause and Cause-Specific Mortality in Women: The Women's Health Initiative. Current Developments in Nutrition. 2020; 4:1493.
  • 10. Kwok CS, Boekholdt SM, Lentjes MA, Loke YK, Luben RN, Yeong JK, Wareham NJ, Myint PK, Khaw KT. Habitual chocolate consumption and risk of cardiovascular disease among healthy men and women. Heart. 2015; 101:1279–87. https://doi.org/10.1136/heartjnl-2014-307050 [PubMed]
  • 11. Mink PJ, Scrafford CG, Barraj LM, Harnack L, Hong CP, Nettleton JA, Jacobs DR Jr. Flavonoid intake and cardiovascular disease mortality: a prospective study in postmenopausal women. Am J Clin Nutr. 2007; 85:895–909. https://doi.org/10.1093/ajcn/85.3.895 [PubMed]
  • 12. Buijsse B, Feskens EJ, Kok FJ, Kromhout D. Cocoa intake, blood pressure, and cardiovascular mortality: the Zutphen Elderly Study. Arch Intern Med. 2006; 166:411–17. https://doi.org/10.1001/archinte.166.4.411 [PubMed]
  • 13. Corti R, Flammer AJ, Hollenberg NK, Lüscher TF. Cocoa and cardiovascular health. Circulation. 2009; 119:1433–41. https://doi.org/10.1161/CIRCULATIONAHA.108.827022 [PubMed]
  • 14. Nüesch E, Trelle S, Reichenbach S, Rutjes AW, Tschannen B, Altman DG, Egger M, Jüni P. Small study effects in meta-analyses of osteoarthritis trials: meta-epidemiological study. BMJ. 2010; 341:c3515. https://doi.org/10.1136/bmj.c3515 [PubMed]
  • 15. Veronese N, Demurtas J, Celotto S, Caruso MG, Maggi S, Bolzetta F, Firth J, Smith L, Schofield P, Koyanagi A, Yang L, Solmi M, Stubbs B. Is chocolate consumption associated with health outcomes? An umbrella review of systematic reviews and meta-analyses. Clin Nutr. 2019; 38:1101–08. https://doi.org/10.1016/j.clnu.2018.05.019 [PubMed]
  • 16. Janszky I, Mukamal KJ, Ljung R, Ahnve S, Ahlbom A, Hallqvist J. Chocolate consumption and mortality following a first acute myocardial infarction: the Stockholm Heart Epidemiology Program. J Intern Med. 2009; 266:248–57. https://doi.org/10.1111/j.1365-2796.2009.02088.x [PubMed]
  • 17. Montagna MT, Diella G, Triggiano F, Caponio GR, De Giglio O, Caggiano G, Di Ciaula A, Portincasa P. Chocolate, "Food of the Gods": History, Science, and Human Health. Int J Environ Res Public Health. 2019; 16:4960. https://doi.org/10.3390/ijerph16244960 [PubMed]
  • 18. Yang Q, Zhang Z, Gregg EW, Flanders WD, Merritt R, Hu FB. Added sugar intake and cardiovascular diseases mortality among US adults. JAMA Intern Med. 2014; 174:516–24. https://doi.org/10.1001/jamainternmed.2013.13563 [PubMed]
  • 19. Seidelmann SB, Claggett B, Cheng S, Henglin M, Shah A, Steffen LM, Folsom AR, Rimm EB, Willett WC, Solomon SD. Dietary carbohydrate intake and mortality: a prospective cohort study and meta-analysis. Lancet Public Health. 2018; 3:e419–28. https://doi.org/10.1016/S2468-2667(18)30135-X [PubMed]
  • 20. Ho FK, Gray SR, Welsh P, Petermann-Rocha F, Foster H, Waddell H, Anderson J, Lyall D, Sattar N, Gill JMR, Mathers JC, Pell JP, Celis-Morales C. Associations of fat and carbohydrate intake with cardiovascular disease and mortality: prospective cohort study of UK Biobank participants. BMJ. 2020; 368:m688. https://doi.org/10.1136/bmj.m688 [PubMed]
  • 21. Mostofsky E, Berg Johansen M, Tjønneland A, Chahal HS, Mittleman MA, Overvad K. Chocolate intake and risk of clinically apparent atrial fibrillation: the Danish Diet, Cancer, and Health Study. Heart. 2017; 103:1163–67. https://doi.org/10.1136/heartjnl-2016-310357 [PubMed]
  • 22. Petrone AB, Gaziano JM, Djoussé L. Chocolate consumption and risk of heart failure in the Physicians' Health Study. Eur J Heart Fail. 2014; 16:1372–76. https://doi.org/10.1002/ejhf.180 [PubMed]
  • 23. Mostofsky E, Levitan EB, Wolk A, Mittleman MA. Chocolate intake and incidence of heart failure: a population-based prospective study of middle-aged and elderly women. Circ Heart Fail. 2010; 3:612–16. https://doi.org/10.1161/CIRCHEARTFAILURE.110.944025 [PubMed]
  • 24. Khodorowsky K, Robert H. The little book of chocolate. France: Flammarion-Père Castor; 2001.
  • 25. Serafini M, Bugianesi R, Maiani G, Valtuena S, De Santis S, Crozier A. Plasma antioxidants from chocolate. Nature. 2003; 424:1013. https://doi.org/10.1038/4241013a [PubMed]
  • 26. Sørensen LB, Astrup A. Eating dark and milk chocolate: a randomized crossover study of effects on appetite and energy intake. Nutr Diabetes. 2011; 1:e21. https://doi.org/10.1038/nutd.2011.17 [PubMed]
  • 27. Katz DL, Doughty K, Ali A. Cocoa and chocolate in human health and disease. Antioxid Redox Signal. 2011; 15:2779–811. https://doi.org/10.1089/ars.2010.3697 [PubMed]
  • 28. Lee J, Taneja V, Vassallo R. Cigarette smoking and inflammation: cellular and molecular mechanisms. J Dent Res. 2012; 91:142–49. https://doi.org/10.1177/0022034511421200 [PubMed]
  • 29. Bernhard D, Wang XL. Smoking, oxidative stress and cardiovascular diseases--do anti-oxidative therapies fail? Curr Med Chem. 2007; 14:1703–12. https://doi.org/10.2174/092986707781058959 [PubMed]
  • 30. Jones MR, Tellez-Plaza M, Navas-Acien A. Smoking, menthol cigarettes and all-cause, cancer and cardiovascular mortality: evidence from the National Health and Nutrition Examination Survey (NHANES) and a meta-analysis. PLoS One. 2013; 8:e77941. https://doi.org/10.1371/journal.pone.0077941 [PubMed]
  • 31. Taubert D, Roesen R, Lehmann C, Jung N, Schömig E. Effects of low habitual cocoa intake on blood pressure and bioactive nitric oxide: a randomized controlled trial. JAMA. 2007; 298:49–60. https://doi.org/10.1001/jama.298.1.49 [PubMed]
  • 32. Grassi D, Desideri G, Necozione S, Lippi C, Casale R, Properzi G, Blumberg JB, Ferri C. Blood pressure is reduced and insulin sensitivity increased in glucose-intolerant, hypertensive subjects after 15 days of consuming high-polyphenol dark chocolate. J Nutr. 2008; 138:1671–76. https://doi.org/10.1093/jn/138.9.1671 [PubMed]
  • 33. Grassi D, Necozione S, Lippi C, Croce G, Valeri L, Pasqualetti P, Desideri G, Blumberg JB, Ferri C. Cocoa reduces blood pressure and insulin resistance and improves endothelium-dependent vasodilation in hypertensives. Hypertension. 2005; 46:398–405. https://doi.org/10.1161/01.HYP.0000174990.46027.70 [PubMed]
  • 34. Grassi D, Lippi C, Necozione S, Desideri G, Ferri C. Short-term administration of dark chocolate is followed by a significant increase in insulin sensitivity and a decrease in blood pressure in healthy persons. Am J Clin Nutr. 2005; 81:611–14. https://doi.org/10.1093/ajcn/81.3.611 [PubMed]
  • 35. Leyva-Soto A, Chavez-Santoscoy RA, Lara-Jacobo LR, Chavez-Santoscoy AV, Gonzalez-Cobian LN. Daily Consumption of Chocolate Rich in Flavonoids Decreases Cellular Genotoxicity and Improves Biochemical Parameters of Lipid and Glucose Metabolism. Molecules. 2018; 23:2220. https://doi.org/10.3390/molecules23092220 [PubMed]
  • 36. Mellor DD, Sathyapalan T, Kilpatrick ES, Beckett S, Atkin SL. High-cocoa polyphenol-rich chocolate improves HDL cholesterol in Type 2 diabetes patients. Diabet Med. 2010; 27:1318–21. https://doi.org/10.1111/j.1464-5491.2010.03108.x [PubMed]
  • 37. Engler MB, Engler MM, Chen CY, Malloy MJ, Browne A, Chiu EY, Kwak HK, Milbury P, Paul SM, Blumberg J, Mietus-Snyder ML. Flavonoid-rich dark chocolate improves endothelial function and increases plasma epicatechin concentrations in healthy adults. J Am Coll Nutr. 2004; 23:197–204. https://doi.org/10.1080/07315724.2004.10719361 [PubMed]
  • 38. Innes AJ, Kennedy G, McLaren M, Bancroft AJ, Belch JJ. Dark chocolate inhibits platelet aggregation in healthy volunteers. Platelets. 2003; 14:325–27. https://doi.org/10.1080/0953710031000123681 [PubMed]
  • 39. Buijsse B, Weikert C, Drogan D, Bergmann M, Boeing H. Chocolate consumption in relation to blood pressure and risk of cardiovascular disease in German adults. Eur Heart J. 2010; 31:1616–23. https://doi.org/10.1093/eurheartj/ehq068 [PubMed]
  • 40. Ding EL, Hutfless SM, Ding X, Girotra S. Chocolate and prevention of cardiovascular disease: a systematic review. Nutr Metab (Lond). 2006; 3:2. https://doi.org/10.1186/1743-7075-3-2 [PubMed]
  • 41. Mozaffarian D, Wu JHY. Flavonoids, Dairy Foods, and Cardiovascular and Metabolic Health: A Review of Emerging Biologic Pathways. Circ Res. 2018; 122:369–84. https://doi.org/10.1161/CIRCRESAHA.117.309008 [PubMed]
  • 42. Grosso G, Micek A, Godos J, Pajak A, Sciacca S, Galvano F, Giovannucci EL. Dietary Flavonoid and Lignan Intake and Mortality in Prospective Cohort Studies: Systematic Review and Dose-Response Meta-Analysis. Am J Epidemiol. 2017; 185:1304–16. https://doi.org/10.1093/aje/kww207 [PubMed]
  • 43. Crichton GE, Elias MF, Alkerwi A. Chocolate intake is associated with better cognitive function: The Maine-Syracuse Longitudinal Study. Appetite. 2016; 100:126–32. https://doi.org/10.1016/j.appet.2016.02.010 [PubMed]
  • 44. Moreira A, Diógenes MJ, de Mendonça A, Lunet N, Barros H. Chocolate Consumption is Associated with a Lower Risk of Cognitive Decline. J Alzheimers Dis. 2016; 53:85–93. https://doi.org/10.3233/JAD-160142 [PubMed]
  • 45. Camandola S, Plick N, Mattson MP. Impact of Coffee and Cacao Purine Metabolites on Neuroplasticity and Neurodegenerative Disease. Neurochem Res. 2019; 44:214–27. https://doi.org/10.1007/s11064-018-2492-0 [PubMed]
  • 46. Cimini A, Gentile R, D'Angelo B, Benedetti E, Cristiano L, Avantaggiati ML, Giordano A, Ferri C, Desideri G. Cocoa powder triggers neuroprotective and preventive effects in a human Alzheimer's disease model by modulating BDNF signaling pathway. J Cell Biochem. 2013; 114:2209–20. https://doi.org/10.1002/jcb.24548 [PubMed]
  • 47. Grosso G, Micek A, Godos J, Sciacca S, Pajak A, Martínez-González MA, Giovannucci EL, Galvano F. Coffee consumption and risk of all-cause, cardiovascular, and cancer mortality in smokers and non-smokers: a dose-response meta-analysis. Eur J Epidemiol. 2016; 31:1191–205. https://doi.org/10.1007/s10654-016-0202-2 [PubMed]
  • 48. Londzin P, Zamora M, Kąkol B, Taborek A, Folwarczna J. Potential of Caffeine in Alzheimer's Disease-A Review of Experimental Studies. Nutrients. 2021; 13:537. https://doi.org/10.3390/nu13020537 [PubMed]
  • 49. Hu FB, Stampfer MJ, Rimm E, Ascherio A, Rosner BA, Spiegelman D, Willett WC. Dietary fat and coronary heart disease: a comparison of approaches for adjusting for total energy intake and modeling repeated dietary measurements. Am J Epidemiol. 1999; 149:531–40. https://doi.org/10.1093/oxfordjournals.aje.a009849 [PubMed]
  • 50. Kircher T, Nelson J, Burdo H. The autopsy as a measure of accuracy of the death certificate. N Engl J Med. 1985; 313:1263–69. https://doi.org/10.1056/NEJM198511143132005 [PubMed]
  • 51. von Elm E, Altman DG, Egger M, Pocock SJ, Gøtzsche PC, Vandenbroucke JP, and STROBE Initiative. The Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) statement: guidelines for reporting observational studies. Prev Med. 2007; 45:247–51. https://doi.org/10.1016/j.ypmed.2007.08.012 [PubMed]
  • 52. Prorok PC, Andriole GL, Bresalier RS, Buys SS, Chia D, Crawford ED, Fogel R, Gelmann EP, Gilbert F, Hasson MA, Hayes RB, Johnson CC, Mandel JS, et al, and Prostate, Lung, Colorectal and Ovarian Cancer Screening Trial Project Team. Design of the Prostate, Lung, Colorectal and Ovarian (PLCO) Cancer Screening Trial. Control Clin Trials. 2000 (Suppl 6); 21:273S–309S. https://doi.org/10.1016/s0197-2456(00)00098-2 [PubMed]
  • 53. Greenberg JA, Manson JE, Neuhouser ML, Tinker L, Eaton C, Johnson KC, Shikany JM. Chocolate intake and heart disease and stroke in the Women's Health Initiative: a prospective analysis. Am J Clin Nutr. 2018; 108:41–48. https://doi.org/10.1093/ajcn/nqy073 [PubMed]
  • 54. Willett WC, Howe GR, Kushi LH. Adjustment for total energy intake in epidemiologic studies. Am J Clin Nutr. 1997 (Suppl 4); 65:1220S–28S. https://doi.org/10.1093/ajcn/65.4.1220S [PubMed]
  • 55. Subar AF, Thompson FE, Kipnis V, Midthune D, Hurwitz P, McNutt S, McIntosh A, Rosenfeld S. Comparative validation of the Block, Willett, and National Cancer Institute food frequency questionnaires : the Eating at America's Table Study. Am J Epidemiol. 2001; 154:1089–99. https://doi.org/10.1093/aje/154.12.1089 [PubMed]
  • 56. Loth KA. Nutritional Data Systems for Research. Singapore: Springer Singapore; 2015. Accessed April 21, 2020. https://doi.org/10.1007/978-981-287-087-2_6-1.
  • 57. Tippett KS, Cypel YS. Design and Operation: The Continuing Survey of Food Intakes by Individuals and the Diet and Health Knowledge Survey, 1994–96. USA: United States Department of Agriculture and Agricultural Research Service; 1997. (Accessed April 21, 2020, at https://www.ars.usda.gov/ARSUserFiles/80400530/pdf/Design.pdf).
  • 58. Krebs-Smith SM, Pannucci TE, Subar AF, Kirkpatrick SI, Lerman JL, Tooze JA, Wilson MM, Reedy J. Update of the Healthy Eating Index: HEI-2015. J Acad Nutr Diet. 2018; 118:1591–602. https://doi.org/10.1016/j.jand.2018.05.021 [PubMed]
  • 59. Spratt M, Carpenter J, Sterne JA, Carlin JB, Heron J, Henderson J, Tilling K. Strategies for multiple imputation in longitudinal studies. Am J Epidemiol. 2010; 172:478–87. https://doi.org/10.1093/aje/kwq137 [PubMed]
  • 60. Steinhaus DA, Mostofsky E, Levitan EB, Dorans KS, Håkansson N, Wolk A, Mittleman MA. Chocolate intake and incidence of heart failure: Findings from the Cohort of Swedish Men. Am Heart J. 2017; 183:18–23. https://doi.org/10.1016/j.ahj.2016.10.002 [PubMed]
  • 61. Maldonado G, Greenland S. Simulation study of confounder-selection strategies. Am J Epidemiol. 1993; 138:923–36. https://doi.org/10.1093/oxfordjournals.aje.a116813 [PubMed]
  • 62. World Health Organization. Report of a WHO consultation on obesity. Obesity: preventing and managing the global epidemic. Geneva, Switzerland: WHO; 1998.
  • 63. Lebwohl B, Cao Y, Zong G, Hu FB, Green PHR, Neugut AI, Rimm EB, Sampson L, Dougherty LW, Giovannucci E, Willett WC, Sun Q, Chan AT. Long term gluten consumption in adults without celiac disease and risk of coronary heart disease: prospective cohort study. BMJ. 2017; 357:j1892. https://doi.org/10.1136/bmj.j1892 [PubMed]
  • 64. Xi B, Veeranki SP, Zhao M, Ma C, Yan Y, Mi J. Relationship of Alcohol Consumption to All-Cause, Cardiovascular, and Cancer-Related Mortality in U.S. Adults. J Am Coll Cardiol. 2017; 70:913–22. https://doi.org/10.1016/j.jacc.2017.06.054 [PubMed]
  • 65. Desquilbet L, Mariotti F. Dose-response analyses using restricted cubic spline functions in public health research. Stat Med. 2010; 29:1037–57. https://doi.org/10.1002/sim.3841 [PubMed]
  • 66. Harrell FE Jr. Regression modeling strategies: with applications to linear models, logistic and ordinal regression, and survival analysis. Springer; 2015. http://citeseerx.ist.psu.edu/viewdoc/download?doi=10.1.1.468.8845&rep=rep1&type=pdf.
  • 67. Kunzmann AT, Coleman HG, Huang WY, Berndt SI. The association of lifetime alcohol use with mortality and cancer risk in older adults: A cohort study. PLoS Med. 2018; 15:e1002585. https://doi.org/10.1371/journal.pmed.1002585 [PubMed]
  • 68. Willett W. Nutritional epidemiology. Oxford: Oxford University Press; 2012.
  • 69. Dietary Guidelines Advisory Committee. Nutrition and Your Health: 2015–2020 Dietary Guidelines for Americans. Washington, DC: US Government Printing Office; 2015.