Abstract
Background Debate about whether high intake of dietary saturated fatty acids (SFA) causes coronary heart disease (the lipid hypothesis) is ongoing.
Methods Using worldwide Global Burden of Disease (GBD) data on cardiovascular disease deaths/100k/year, ages 15-69 years old in male and female cohorts (CVD) and dietary and other CVD risk factors, we formatted and population weighted data from 195 countries. The formatted rows of data (n=7846 cohorts) each represented about 1 million people, totaling about 7.8 billion people in 2020. We correlated CVD with dietary and other risk factors worldwide and in appropriate subsets. Outcome measures included CVD versus dietary and other risk factor correlations worldwide and in subsets.
Findings After empirical data exploration, we defined a “fat-soluble vitamin variable” Kcal/day (FSVV) as, “FSVV=processed meat + red meat + fish + milk + poultry + eggs + (SFA + polyunsaturated fatty acids (PUFA) + trans fatty acids (TFA)) * 0.46 (all in Kcal/day).” Low density lipoprotein cholesterol mmol/L correlated strongly positively with FSVV worldwide (r=0.780, 95% CI 0.771 to 0.788, p<0001, n=7846 cohorts), so we considered FSVV our marker variable to test both the lipid hypothesis and our fat-soluble vitamin hypothesis. LDL-C correlated negatively with CVD worldwide (r= -0.279, 95% CI -0.299 to -0.258, p<0.0001 as did FSVV versus CVD (r= -0.329, 95% CI -0.349 to -0.309, p<0.0001). However, FSVV correlated positively with CVD in the highest FSVV cohorts (mean male/female FSVV ≥567.27 Kcal/day: r=0.523, 95% CI 0.476 to 0.567, p<0001, n=974 cohorts).
Interpretation Since FSVV correlated positively with CVD only in high FSVV cohorts, the data supported the lipid hypothesis only in GBD cohorts with high FSVV intake. Since both FSVV and LDL-C correlated negatively with CVD worldwide, only the fat-soluble vitamin hypothesis was supported worldwide. This GBD cohort data analysis methodology could be used to help develop food policy and education strategies for improving public health.
Funding none
Evidence before this study Relating to dietary and other risk factors for early death from cardiovascular diseases worldwide, the EAT-Lancet Commission proposed the “Planetary Health Diet” to reduce non-communicable diseases and mitigate climate change. The Planetary Health Diet has been controversial and did not achieve the goal of leading a worldwide “Great Food Transformation” towards a more plant-based diet. The comparative risk assessment (CRA) systematic literature review-based methodology used in drafting the Planetary Health Diet has no significant competitors as published methodologies used to parse health effects of worldwide dietary and other risk factors. No consensus exists about an optimal range of diets for optimum human health, including cardiovascular disease prevention/treatment.
Added value of this study□ Cardiovascular disease deaths/100k/year in 15-69-year-old males and females (CVD) worldwide correlated with dietary risk factors had two patterns based on the amounts of kilocalories/day (Kcal/day) of animal foods and added fats. We differentiated the two patterns by defining a “fat-soluble vitamin variable” Kcal/day (FSVV=processed meat + red meat + fish + milk + poultry + eggs + (added saturated fatty acids (SFA)+polyunsaturated fatty acids (PUFA)+trans fatty acids (TFA)) * 0.46 (all in Kcal/day). Worldwide, all nine risk factors comprising FSVV correlated negatively with CVD, but in high FSVV cohorts (FSVV≥567.27 Kcal/day), all nine risk factors in the FSVV correlated positively with CVD.
Implications of all the available evidence Animal foods and added fats afforded significant protection from CVD worldwide. The data suggested that the fat-soluble vitamins in animal foods and added fat facilitating gut absorption may have been partially or entirely responsible for the worldwide protection from CVD. However, in high FSVV intake cohorts, SFA in food and in extracted added fats facilitating atherosclerosis may have partially or entirely accounted for the positive correlation of FSVV and CVD. The findings suggest that a yet not precisely defined moderate amount of animal food and added fat intake contributed to lower cardiovascular risk relative to very high or very low intakes.
Introduction
The scientific validity of the Dietary Guidelines for Americans for 2015-2020, including guidelines based on the lipid hypothesis, was challenged by Journalist Nina Teicholz in the BMJ.1 The Center for Science in the Public Interest called for the BMJ to retract the article. Peer reviewers the BMJ selected to adjudicate the far-reaching dispute concluded that, “Teicholz’s criticisms of the methods used by Dietary Guidelines for Americans Committee are within the realm of scientific debate.”2 Currently, no methodology for relating cardiovascular disease events to food intake has been generally accepted as rigorous, replicable, and scientifically valid.
This paper will use population weighted, formatted, worldwide global burden of disease (GBD) data from the Institute of Health Metrics and Evaluation (IHME) to assess the impact of diet, metabolic risk factors, tobacco use, air pollution, and other risk factors on cardiovascular disease deaths/100k/year in male and female cohorts 15-69-years-old (CVD).
We were open to the possible support of the data for the lipid hypothesis but also for our newly defined, “fat-soluble vitamin hypothesis,” which an initial exploration of the GBD data prompted. We defined it as “insufficient dietary intake of fat-soluble vitamin containing animal foods and added fat for gut absorption increases cardiovascular disease risk.” We will elaborate in the results and the discussion.
Methods
As volunteer collaborators with the IHME, we acquired and utilised raw GBD worldwide ecological data (GBD: 195 countries and 365 subnational locations, n=1120 male and female cohorts). Data consisted of the rate of cardiovascular disease early deaths, metabolic risk factors, dietary risk factors and covariates, and other risk factors of male and female cohorts 15-49 years old and 50-69 years old from each year 1990-2017. GBD worldwide citations of over 12,000 surveys constituting ecological data inputs for this analysis are available online from IHME.3 The main characteristics of IHME GBD data sources, the protocol for the GBD study, and all risk factor values have been published by IHME GBD data researchers and discussed elsewhere.4–8 These include detailed descriptions of categories of input data, potentially important biases, and the methodologies of analysis. We did not clean or pre-process any of the GBD data. GBD cohort risk factor and health outcome data from the IHME had no missing records.
Food risk factors came from surveys that IHME utilised as gram/day (g/d) consumed on average. GBD dietary covariate data originally came from Food and Agriculture Organization surveys of animal and plant food commodities available per capita in countries worldwide (potatoes, corn, rice, sweet potatoes, poultry, and eggs)—as opposed to consumed on average.9
Supplementary Table 1 lists the relevant GBD dietary and other risk factors, covariates, and other available variables with definitions of those risk factor exposures.7 Once IHME staff time constraints due to COVID-19 data collection and modeling are over, the updated GBD 2020 data with all the variables as the GBD 2017 data that we used for this analysis may be obtained by volunteer researchers collaborating with IHME.10
Study design and population
For cardiovascular disease early deaths and all risk factors, we averaged the values for ages 15-49 years old together with 50-69 years old for each male and female cohort for each year. Finally, for each male and female cohort, data from all 28 years (1990-2017) of the means of the rate of CVD and risk factor exposures were averaged using the computer software program R.
To weigh the data according to population, internet searches (mostly Wikipedia) yielded the most recent population estimates for countries and subnational states, provinces, and regions. The 1120 GBD cohorts available were population weighted by a software program in R, resulting in an analysis dataset with 7846 population weighted cohorts representing about 7.8 billion people in 2020. Each male or female cohort in the population-weighted analysis dataset represented approximately 1 million people (range: < 100,000 to 1.5 million). World population data from the World Bank or the Organization for Economic Co-operation and Development could not be used because they did not include all 195 countries or any subnational data.
Supplementary Table 2 details how omega-3 fatty acid g/d was converted to fish g/d using data on the omega-3 fatty acid content of frequently eaten fish from the National Institutes of Health Office of Dietary Supplements (USA).11 As shown in Supplementary Table 3, we converted all of the animal and plant food data, including alcohol and sugary beverage consumption, from g/d to kilocalories/day (Kcal/day). For the g/d to Kcal/day conversions, we used the Nutritionix track app,12 which tracks types and quantities of foods consumed. Saturated fatty acids (SFA: 0-1 portion of the entire diet Kcal/day) was not available with GBD 2017 data, so GBD SFA risk factor data from GBD 2016 was used. Polyunsaturated fatty acid (PUFA) and trans fatty acid (TFA) GBD risk factor data from 2017 (0-1 portion of the entire diet Kcal/day) were also utilised, but monounsaturated fat data were not available. These fatty acid data expressed as 0-1 portion of the entire diet were converted to Kcal/day by multiplying by the total Kcal/day available (a covariate from the Food and Agriculture Organization9) per capita for each cohort.
Outcome variable and covariates
CVD, the principal outcome variable of this analysis, was a combination variable consisting of the deaths/100k/year of male and female cohorts from (1) ischemic heart disease, (2) stroke, (3) hypertensive heart disease, (4) rheumatic heart disease, (5) non-rheumatic valvular disease, (6) subarachnoid haematoma, (7) myocarditis, (8) alcoholic cardiomyopathy, (9) endocarditis, aortic aneurysm, and (10) atrial fibrillation. Given the multiple dietary and other influences on the causation of CVD, we did not attempt to differentiate between dietary risk factors and covariates.
Statistical methods
To determine the strengths of the risk factor correlations with CVD of population weighted worldwide cohorts (n=7846 cohorts) or subgroups of cohorts (e.g., highest and lowest CVD cohorts, etc.), we utilised Pearson correlation coefficients: r, 95% confidence intervals (CIs), and p values. To identify confounding in the univariate correlations, we used partial correlations, holding constant the suspected confounder variables.
Sensitivity analyses
We performed subset analyses on the 1000 cohorts with the lowest and highest levels of CVD.
We used SAS and SAS OnDemand for Academics software 9.4 (SAS Institute, Cary, NC) for the data analysis.
Results
Table 1 shows the 44 dietary, metabolic, and other worldwide risk factors potentially contributing to CVD (n=7846 cohorts). See Supplementary Table 1 for definitions of the risk factor and covariate variables. In Table 1, low density lipoprotein cholesterol mmol/L (LDL-C) negatively correlated with CVD worldwide (r=-0.279, 95% CI -0.299 to -0.258, p<0001,), meaning the higher the LDL-C the lower the CVD. Figure 1 gives the worldwide plot of LDL-C versus CVD, showing widely scattered data points with a negative slope of the least square regression line.
Table 1 shows that CVD negatively correlated with all the animal foods and all the fatty acids (processed meats, red meats, fish, milk, poultry, eggs, SFA, PUFA, and TFA). Based on these unexpected finding, we formulated our “fat-soluble vitamin hypothesis”—insufficient dietary intake of fat-soluble vitamin containing animal foods and added fat for gut absorption increases cardiovascular disease risk. To test this hypothesis, we summed the Kcal/day of animal foods and added fats to produce a fat-soluble vitamins variable (Kcal/day) (FSVV =processed meats Kcal/day + red meats Kcal/day + fish Kcal/day + milk Kcal/day + poultry Kcal/day + eggs Kcal/day + added SFA Kcal/day + added PUFA Kcal/day + added TFA Kcal/day). In determining the portion of SFA, PUFA, and TFA added in addition to the fatty acids in the animal and plant foods, an adjustment factor (fatty acids * 0.46), adapted from the website “Our World in Data,” differentiated the fatty acids in individual foods (54% of the total) from the added fatty acids in oils, butter, lard, etc. (46% of the total).13 This adjustment prevented double counting fatty acids (Kcal/day) and allowed a more accurate determination of FSVV.
FSVV (with abundant SFA in animal foods and added extracted fat) correlated strongly with LDL-C worldwide (r=0.780, 95% CI 0.771 to 0.788, p<0001, n=7846 cohorts). Consequently, we considered that FSVV could serve as our marker for the lipid hypothesis as well as the fat-soluble vitamin hypothesis.
Table 2 shows a subgroup of 500 pairs of male and female cohorts with the lowest CVD values, totaling 1000 total cohorts (an arbitrary number). Notably, the mean FSVV in this low CVD subset (FSVV=531.38 Kcal/day) was much higher than the mean FSVV worldwide (FSVV=285.36 Kcal/day, Table 1) or the highest CVD subset mean FSVV (FSVV=235.66 Kcal/day, n=1000 cohorts, Table 2).
To test the lipid hypothesis with FSVV as the marker variable, we calculated the median FSVV of the low CVD 1000 cohorts (median FSVV=567.27 Kcal/day) and formed a subset of the highest FSVV cohorts (FSVV≥567.27 Kcal/day, n=974 cohorts, Table 3). In Table 3, we see a much higher mean FSVV (FSVV≥692.58 Kcal/day) than in Table 2 and a strong positive correlation of FSVV with CVD (r=0.523, 95% CI 0.476 to 0.567, p<0001, n=974 cohorts). In contrast, LDL-C and CVD were negatively correlated in this high FSVV subset (r=-0.254, 95% CI -0.194 to -0.312, p<0001, n=974 cohorts). Figure 2 graphs the contrast between worldwide FSVV negatively correlating with CVD versus FSVV positively correlating with CVD in the highest FSVV cohorts (FSVV≥567.27 Kcal/day).
Table 4 shows FSVV’s component risk factors (processed meat, red meat, fish, milk, poultry, eggs, and added (SFA+PUFA+TFA)) for the 34 countries in Table 2 in descending order by CVD. Table 5 gives all the CVD risk factors for seven representative countries out of the 34 countries in Table 4 (CVD≤ 293.97, n=1000 cohorts).
The average alcohol intake was 81.0 Kcal/day worldwide. Alcohol Kcal/day correlated negatively with CVD worldwide (r=-0.061, 95% CI -0.083 to -0.039, p<0001, Table 1). However, in high FSVV cohorts (FSVV≥567.27 Kcal/day), mean alcohol=155.02 Kcal/day), alcohol positively correlated with CVD (r=0.475, 95% CI 0.425 to 0.522, p<0001, n=974 cohorts, Table 3).
Vitamin A deficiency incidence/100k in children ≤ 5 years old, the only fat-soluble vitamin in the GBD database, correlated positively with CVD worldwide (r=0.210, 95% CI 0.189 to 0.231, p<0001, Table 1).
Dietary fiber was not significantly correlated with CVD worldwide (r=0.019, 95% CI -0.003 to 0.041, p=0.09, Table 1). Once dietary fiber was adjusted for FSVV, alcohol, sugary beverages, and potatoes by partial correlation analysis; dietary fiber correlated negatively with CVD (r=- 0.052, 95% CI -0.074 to -0.030, p<0.0001).
Physical activity positively correlated with CVD worldwide (r=0.160, 95% CI 0.139 to 0.182, p<0.0001, Table 1). However, confounding factors existed:
countries with higher FSVV levels also had less physical activity (FSVV versus physical activity worldwide: r=-0.366, 95% CI: -0.385 to -0.347, p< 0.0001),
females had less physical activity than males (mean physical activity males=5061 METs/week, mean physical activity females=4356 METs/week), along with less CVD (mean CVD males=656.2, mean CVD females=428.6), and
males smoked tobacco more than females (mean smoking (0-1) males=0.339 or 33.9%, mean smoking (0-1) females=0.070 or 7%), and males had more physical activity (as #2 above).
Once physical activity was adjusted for FSVV, smoking, and sex by partial correlation analysis, it correlated negatively with CVD (r=-0.127, 95% CI -0.149 to -0.104, p<0.0001).
In cohorts mostly in developed countries, with FSVV ≥567.27 Kcal/day, discontinuing breast feeding correlated positively with CVD (r=0.268, 95% CI 0.208 to 0.325, p<0.0001, Table 3).
Secondhand smoking negatively correlated with CVD worldwide (r=-0.225, 95% CI -0.246 to - 0.204, p<0.0001, Table 1). After adjusting secondhand smoking for (1) smoking, (2) sublingual tobacco use, (3) household smoke, (4) ambient air pollution, and (5) sex by partial correlation analysis; secondhand smoke correlated positively with CVD (r=0.048, 95% CI 0.025 to 0.070, p<0.0001).
The metabolic risk factors, BMI, FPG, and LDL-C, all correlated negatively with CVD worldwide (BMI r=-0.240, 95% CI -0.261 to -0.219, p<0.0001, FPG r=-0.178, 95% CI -0.200 to -0.157, p<0.0001, and LDL-C r=-0.279, 95% CI -0.299 to -0.258, p<0.0001, respectively, Table 1). Once they were adjusted for FSVV, alcohol, sugary beverages, and potatoes by partial correlation analysis; they all correlated positively with CVD (BMI r=0.028, 95% CI 0.006 to 0.050, p=0.0144, FPG r=0.076, 95% CI 0.054 to 0.098, p<0.0001, LDL-C r=0.170, 95% CI 0.148 to 0.191, p<0.0001, respectively). In the cohorts with FSVV≥567.27 Kcal/day, BMI and FPG positively correlated with CVD (BMI r=0.484, 95% CI 0.434 to 0.530, p<0.0001 and FPG r=0.091, 95% CI 0.028 to 0.153, p=0.0046, respectively).
Discussion
Since this GBD data-based worldwide risk factor—health outcome analysis methodology has no precedent, there is no medical literature of previous such studies. The best comparator for this CVD versus risk factor analysis with IHME GBD data would be the comparative risk assessment (CRA) systematic literature review-based methodology used in drafting the “Planetary Health Diet.” Proposed by the EAT-Lancet Commission, the Planetary Health Diet was meant to spark a global transformation in human diets to improve health and mitigate climate change. The EAT- Lancet Commission authors said, “Transformation to healthy diets by 2050 will require substantial dietary shifts, including a greater than 50% reduction in global consumption of unhealthy foods such as red meat and sugar, and a greater than 100% increase in the consumption of healthy foods such as nuts, fruits, vegetables and legumes.”14 Reducing diet-caused non-communicable diseases (e.g., obesity, type 2 diabetes, CVD, cancers, etc.) was a major goal. WHO data show that about 44% of worldwide non-communicable disease deaths were CVD related,15 so these CVD analysis findings from GBD data can be fairly contrasted with the Planetary Health Diet proposal from CRA systematic literature review and expert opinion methodology.
When introduced by the EAT-Lancet Commission, the Planetary Health Diet sparked controversy. Dr. Francisco Zagmutt and colleagues disputed the rigor of the analytics, “A truly effective global solution to the problem of human nutrition and environmental impact must be replicable, transparent, and supported with correct quantification of its impact. Unfortunately, the report did not meet these criteria.”16 The World Health Organization (WHO) withdrew sponsorship over the exclusivity of one dietary approach for the entire planet and the negative impact of reducing cattle and pig ranching in developing countries.17
The methodologies and findings of this GBD data-based analysis of diet and other risk factors related to global CVD were markedly different than those of the EAT-Lancet Commission. Some of the consequential differences were as follows:
The Planetary Health Diet recommendations were based on systematic literature reviews primarily from researchers and participants in developed countries, largely missing at least 70% of the world’s population. This IHME GBD data-based population weighted analysis proportionally included people in developing countries.
To Dr. Zagmutt’s point about the lack of rigor of the Planetary Health Diet analytics, this IHME GBD data analysis used a single population weighted analysis dataset, spanning 28 years, making it replicable (with 2020 GBD data), transparent, and supported with reasonable quantification of the impact of global human nutrition on CVD.
The overall amounts of animal foods (Kcal/day) recommended by the Planetary Health Diet (mean animal foods: 304 Kcal/day)18 do not differ substantially from the overall animal foods (Kcal/day) in the GBD lowest CVD cohorts (338 Kcal/day, n=1000 cohorts, CVD < 292.97 deaths/year) in this analysis (Table 2). However, the Planetary Health Diet recommended 40 Kcal/day of fish18 while the GBD data show that the world fish consumption from 1990-2017 averaged 10.0 Kcal/day (Table 1). Levels of overfishing are unsustainable now, so quadrupling fish consumption worldwide is impossible. Likewise, the Planetary Health Diet recommended average milk consumption to be 153 Kcal/day18 while the GBD data show the world average milk consumption per person from 1990-2017 of 25.0 Kcal/day (Table 1). Increasing the world’s dairy cows six-fold, most of whom would eventually be slaughtered for meat, would seem to contradict major goals of the EAT-Lancet Commission—reducing methane emissions from cows and reducing red meat consumption. Instead, increasing all livestock production in developing countries would be projected to improve CVD outcomes (Table 1) as would decreasing cattle, pig, and poultry raising, but not fish or eggs, in high FSVV cohorts (Table 3).
Whereas the Planetary Health Diet has a single set of recommended ranges of quantities of plant and animal foods for all humanity, the GBD data-based analysis methodology can apply to people in developing countries as well as developed countries.
Fat-soluble vitamins include vitamins A, D, E, and K (K1=phylloquinones and K2= menaquinones). Vitamin K2 (menaquinones from animal foods and fermented plant products) inhibits arterial calcification, a major factor in the pathogenesis of atherosclerosis and cardiovascular disease.19 Vitamin A and its beta-carotene precursor are important in preventing congenital heart disease.20, 21 Vitamin E protects the vasculature against oxidative stress, including lipid peroxidation and the production of atherogenic forms of LDL-C—a part of the pathogenesis of cardiovascular disease.22 Regarding vitamin D, a meta-analysis of prospective studies showed an inverse correlation between cardiovascular disease and circulating 25(OH)-vitamin D ranging from 20 to 60 nmol/L.23
The four fat-soluble vitamins come primarily from animal foods. Vitamin K2 levels are high in aged cheeses, other fermented dairy products,24 organ meats (e.g., bovine and pork liver), chicken, and egg yolks.25 The Japanese had the world’s lowest CVD from 1990-2017 (mean CVD=169.2 deaths/100k/year). Paradoxically, crucial CVD risk factors were high in the Japanese (e.g., smoking: mean male/female smoking=26.8%, salt intake (mean sodium= 6.01 g/d, systolic blood pressure (mean SBP= 137.6 mm/Hg, Table 5). Relative to fermented cheeses, the content of vitamin K2 in fish was low to negligible (Japanese fish consumption mean=260.5 Kcal/day, Table 4).25 Instead, Japanese consumed high quantities of vitamin K2 in the form of fermented soy products like natto and miso.26, 27 People that also consumed high amounts of fermented legumes were the Taiwanese (soybean Douchi, soybean Meitauza, and soybean curd28), South Koreans (Cheonggukjang, a fermented soybean stew28), and Thais (soybean Thua nao29). The Taiwanese, South Koreans, and Thais were also listed among the 34 low CVD countries (Table 4).
In high FSVV cohorts, (FSVV≥ 567.27 Kcal/day), why LDL-C correlated negatively with FSVV (r=-0.254 95% CI -0.312 to -0.194, p<0.0001, n=974 cohorts) is not clear. However, multiple studies have documented that LDL-C may be influenced by genetic, racial, or ethnic differences that do not necessarily correspond with risk of CVD.30, 31 Additionally, LDL-C lowering medications in wide use especially in high FSVV intake developed countries may also confound the relationship of LDL-C with both LSVV and CVD risk in this subset cohort analysis.
Since there are high FSVV and low FSVV intake individual people widely distributed worldwide, prospective observational studies of individuals in high FSVV areas, rather than GBD cohort studies like this one, could better demonstrate positive correlations between LDL-C (or total cholesterol) and CVD. This may explain why the Seven Countries Study32 of Ancel Keys and the Framingham Heart Study,33 both in high FSVV regions, supported the lipid hypothesis. However, the Prospective Urban Rural Epidemiology (PURE) Study, with individual data on 135,335 patients from three high-income, 11 middle-income, and four low-income countries was clearly at variance with the lipid hypothesis.34 The PURE Study found that higher saturated fat intake was associated with lower risk of stroke (quintile 5 vs quintile 1, HR 0.79, (95% CI 0.64 to 0.98), ptrend=0.0498).
Defining the optimal FSVV range that will reduce CVD risk while not increasing other diseases (e.g., cancers and obesity) will require more research. Supplements of fat-soluble vitamins to reduce the animal food and added fat intake required to minimise CVD should be studied, especially in developing countries with less access to animal foods.
Given the negative correlations of CVD with both fish and eggs in both tables 1 and 3, the reduction of FSVV in high FSVV countries might best come from meat, poultry, and added fats. Although milk (including other dairy products) correlated positively with CVD in Table 3 (r=0.185, 95% CI 0.123 to 0.245, p<0001), this probably related to Japan and Peru having the lowest and fifth lowest CVD (Japan: mean CVD=169.2, mean milk=29.0 Kcal/day, Peru: mean CVD=197.7, mean milk=15.3 Kcal/day). Without Japan and Peru in the analysis (n=834 cohorts), milk correlated negatively with CVD (r=-0.071, 95% CI -0.138 to -0.030, p=0.0405). Of the 34 low CVD countries in Table 4, 20 had high milk intakes (milk Kcal/day range: 68.8 Kcal/day-124.1Kcal/day versus mean milk worldwide=25.0 Kcal/day, Table 1) and these high milk intake cohorts all had relatively high per capita intakes of cheese.35, 36 Fermented, full-fat dairy products have high levels of fat-soluble vitamin K2,24, 25 which have been found to be important in preventing calcification of arteries and atherosclerosis.37
Potatoes Kcal/day available positively correlated with CVD both worldwide (r=0.050, 95% CI 0.028 to 0.073, p< 0.0001, Table 1) and in cohorts with FSVV≥ 567.27 Kcal/day (r=0.109, 95% CI 0.046 to 0.170, p< 0.0001, Table 3). Notably, half, or more of the potatoes consumed worldwide are in the form of highly processed food products.38 Recent large prospective observational studies have found higher consumption of ultra-processed foods associated with an increased risk of cardiovascular disease incidence and CVD mortality.39, 40 Data from 79 high- and middle-income countries showed that ultra-processed products dominate the food supplies of high-income countries and that their consumption is now rapidly increasing in middle-income countries.41
Dietary sodium negatively correlated with CVD worldwide (r=-0.213, 95% CI -0.234 to -0.192, p<0.0001, Table 1) and was not significantly correlated with CVD in high FSVV cohorts (FSVV≥ 567.27 Kcal/day: r=-0.040, 95% CI 0.023 to 0.170, p=0.21, Table 3) This was contrary to an analysis of sodium intake versus cardiovascular deaths in 66 countries by Mozaffarian, et. al who attributed 9.5% of cardiovascular deaths worldwide to high sodium intake.42 However, O’Donnell and colleagues found the relationship of sodium to cardiovascular disease to be J shaped curve and suggested that the lowest cardiovascular risk is with moderate sodium intake in the 3-5 g/d range.43 This analysis is consistent with O’Donnell’s study.
The moderately strong positive correlation of child ≤5 years old severe underweight (weight >2 SD below the mean for height) with CVD worldwide (r=0.306, 95% CI 0.285 to 0.325, p<0.0001, Table 1) suggested a relationship of infant/child malnutrition with later CVD. Babies surviving the so called “Dutch famine” toward the end of the second world war (1944-45) have been shown to have higher subsequent heart disease incidence than earlier or later cohorts in Holland.44 This suggested that in utero and early infancy severe malnutrition may subsequently contribute substantially to CVD, especially in developing countries.
Prematurely stopping breast feeding positively correlated with CVD in the countries with FSVV ≥567.27 Kcal/day (r=0.268 95% CI: 0.208 to 0.325, p=0.0027, Table 3). A review of observational studies of breast feeding related to metabolic risk factors for CVD in developed countries suggested that breast feeding was associated with increased insulin sensitivity and decreased systolic blood pressure in later life. Breast feeding also had metabolic benefits for the mother.45
According to the WHO, ambient air pollution (particulate matter ≤2.5 micrometers diameter (PM 2.5)) causes cardiovascular and respiratory diseases and cancers.46 The WHO considers >10 μg/m3 of PM2.5 particles a health hazard. In this GBD data analysis, air pollution (PM2.5) in countries with FSVV≥567.27 Kcal/day and mean CVD=273.7 deaths/100k/year averaged 12.7 μg/m3. Worldwide, comparable numbers for air pollution (PM 2·5) were mean=44.7 μg/m3 and mean CVD=543.7 (Table 1).
These IHME GBD data support the link of smoking tobacco to CVD (r=0.298, 95% CI 0.278 to 0.318, p<0.0001, Table 1). The mean incidence of smoking worldwide=20.5% (Table 1). In the 34 countries with the lowest CVD, the mean smoking incidence=23.3% (n=974 cohorts, Table 3). Smoking was clearly a major risk factor for CVD, but these data suggest that diet was more influential. Given the multiple dietary and other contributors to CVD, the CVD population attributable fraction (PAF) for smoking was probably far below the WHO estimate that tobacco accounts for 20% of coronary artery disease mortality.47
The positive correlation of blood lead levels with CVD worldwide (r=0.180, 95% CI 0.159 to 0.201, p<0.0001, Table 1) was consistent with reports that lead toxicity leads to hypertension.48 Also, environmental lead exposure has been linked with all-cause mortality, cardiovascular disease mortality, and ischaemic heart disease mortality in the US.49
This analysis showed that kidney disease correlated positively with CVD (r=0.194, 95% CI 0.173 to 0.215, p<0.0001, Table 1). A systematic literature review of kidney disease revealed that chronic kidney disease resulted in an estimated 1.2 million deaths in 2017, of which a large portion were from CVD.50 Other CVD risk factors may also lead to kidney disease (e.g., systolic hypertension, types 1 and 2 diabetes).
Limitations
Our study was subject to all the limitations discussed in previous GBD publications.51, 52 These included gaps, biases, and inconsistencies in data sources as well as limitations in the methods of data processing and estimation. Having comprehensive data on dietary inputs is key to more accurate and reliable analyses. These GBD data on animal foods, plant foods, alcohol, sugary beverages, and fatty acids were not comprehensive and comprised only 1191.4 Kcal/day per person on average worldwide. Subnational data on all risk factors were available on only four countries. Because the data formatting and statistical methodology were new, this was necessarily a post hoc analysis and no pre-analysis protocol was possible. This GBD data analysis should be repeated with the most recently released GBD 2019 data when those data become available to volunteer collaborators.
Generalisability
This new the GBD statistical analysis methodology for finding correlations between dietary and other risk factors and health outcomes (e.g., CVD) can be applied to any of the risk factors and health outcomes available with the new IHME GBD 2020 data. With further experience with the methodology, we can determine the generalisability.
Conclusion
The lipid hypothesis (e.g., dietary SFA from animal foods and added fats causes CVD) was consistent with the strong positive correlation of FSVV with CVD in high FSVV cohorts. However, the lipid hypothesis was not supported worldwide. Global data analysis supported the fat-soluble vitamin hypothesis because FSVV negatively correlated with CVD worldwide. This methodology of analysing IHME GBD data might well have advantages over the CRA systematic literature review and expert opinion method in developing policy recommendations, clinical practice guidelines, and public health recommendations. The EAT-Lancet Commission might consider using IHME GBD data to guide another version of a planetary health diet.
Data Availability
Our dataset, SAS codes, and Excel files are posted on the Mendeley data repository.
Article information
Contributors
DKC had full access to the GBD data in the study and takes responsibility for the integrity of the data and the accuracy of the data analysis. DKC conceived and designed the study, acquired and analysed the GBD data from IHME, interpreted the study findings, drafted the manuscript, critically reviewed and edited the manuscript and tables, and approved the final version for publication.
CW designed software programs in R to format and population weight the data, aided with the SAS statistical analysis, critically reviewed the manuscript, and approved the final version for publication.
Competing interest disclosures
None reported.
Competing interests’ statement
Both authors have completed the ICMJE uniform disclosure form at www.icmje.org/coi_disclosure.pdf and declare: no support from any Organization for the submitted work; no financial relationships with any organizations that might have an interest in the submitted work in the previous three years; no other relationships or activities that could appear to have influenced the submitted work.
Funding/Support
This research received no grant from any funding agency in the public, commercial or not-for-profit sectors. The Bill and Melinda Gates Foundation funded the acquisition of the data by the IHME for this analysis.
Role of the Funder/Sponsor
The IHME provided the data for this analysis but neither the IHME nor Bill and Melinda Gates Foundation had any role in the analysis of the data.
Data Sharing Statement
The formatted analysis dataset for this analysis, SAS codes, and Excel files are posted on the Mendeley data repository https://data.mendeley.com/datasets/g6b39zxck4/6)
Additional Contributions
We thank Scott Glenn and Brent Bell from IHME who supplied us with the GBD risk factor exposure data.
Checklist
This report followed the STrengthening the Reporting of OBservational studies in Epidemiology (STROBE) guidelines for reporting global health estimates.53
Footnotes
Editing for clarity and word count was done on the abstract. The text was edited and updated. A comparison with the EAT-Lancet Planetary Health Diet was done.