Skip to main content

Cooking oil/fat consumption and deaths from cardiometabolic diseases and other causes: prospective analysis of 521,120 individuals

Abstract

Background

Increasing evidence highlights healthy dietary patterns and links daily cooking oil intake with chronic diseases including cardiovascular disease (CVD) and diabetes. However, food-based evidence supporting the consumption of cooking oils in relation to total and cardiometabolic mortality remains largely absent. We aim to prospectively evaluate the relations of cooking oils with death from cardiometabolic (CVD and diabetes) and other causes.

Methods

We identified and prospectively followed 521,120 participants aged 50–71 years from the National Institutes of Health-American Association of Retired Persons Diet and Health Study. Individual cooking oil/fat consumption was assessed by a validated food frequency questionnaire. Hazard ratios (HRs) and 95% confidence intervals (CIs) were estimated for mortality through the end of 2011.

Results

Overall, 129,328 deaths were documented during a median follow-up of 16 years. Intakes of butter and margarine were associated with higher total mortality while intakes of canola oil and olive oil were related to lower total mortality. After multivariate adjustment for major risk factors, the HRs of cardiometabolic mortality for each 1-tablespoon/day increment were 1.08 (95% CI 1.05–1.10) for butter, 1.06 (1.05–1.08) for margarine, 0.99 (0.95–1.03) for corn oil, 0.98 (0.94–1.02) for canola oil, and 0.96 (0.92–0.99) for olive oil. Besides, butter consumption was positively associated with cancer mortality. Substituting corn oil, canola oil, or olive oil for equal amounts of butter and margarine was related to lower all-cause mortality and mortality from certain causes, including CVD, diabetes, cancer, respiratory disease, and Alzheimer’s disease.

Conclusions

Consumption of butter and margarine was associated with higher total and cardiometabolic mortality. Replacing butter and margarine with canola oil, corn oil, or olive oil was related to lower total and cardiometabolic mortality. Our findings support shifting the intake from solid fats to non-hydrogenated vegetable oils for cardiometabolic health and longevity.

Peer Review reports

Background

Cooking oils/fats are known as edible oils from vegetable or animal origin and used for cuisine or salad preparation worldwide. To meet consumption needs, global vegetable oil production was closed to 198 million metric tons in 2018–2019. For animal or artificial cooking fats, 187.1 million Americans used margarine or margarine spread while butter consumption in the USA reached 5.8 pounds per capita in 2019 [1]. Growing controversy focused on the role of cooking oils/fats in the incidence of various chronic disorders including cardiovascular disease (CVD) events. Importantly, the well-known functions of saturated fatty acids (SFAs), monounsaturated fatty acids (MUFAs), or polyunsaturated fatty acids (PUFAs) could not apparently translate to cardiometabolic health effects of cooking oils/fats [2]. Thus, increasing evidence supports shifting away from isolated fatty acids toward food-based patterns for linking dietary cooking oils/fats with all-cause mortality, fatal and nonfatal CVD events [3, 4]. Vegetable cooking oils are regarded as the healthier choice as they contain more unsaturated fatty acids than animal oils. Canola oil and corn oil may ameliorate blood lipid profile and protect against CVD risk factors [5, 6], whereas butter raises total and LDL cholesterol levels [7]. Canola oil- or olive oil-enriched diet could improve glycemic control in patients with type 2 diabetes [6, 8]. However, only a few studies have provided weak evidence of cooking oil/fat consumption in relation to all-cause and cardiometabolic mortality [9,10,11]. The associations of lard, canola oil, and corn oil consumption with cardiometabolic mortality remain lacking.

To fill these gaps, we assessed the long-term associations of 6 typical cooking oils/fats, including butter, margarine, lard, corn oil, canola oil, and olive oil, with all-cause, cardiometabolic, and other major cause-specific mortality in the National Institutes of Health-American Association of Retired Persons (NIH-AARP) Diet and Health Study.

Methods

Study population

The NIH-AARP Diet and Health Study is a large prospective cohort consisting of 617,119 US men and women aged 50–71 years. At baseline in 1995–1996, validated questionnaires were mailed to 3.5 million AARP members from 6 US states (California, Florida, Louisiana, New Jersey, North Carolina, and Pennsylvania) and 2 metropolitan areas (Atlanta, Georgia, and Detroit, Michigan) to collect data on demographics, lifestyle, and dietary characteristics. All participants provided written informed consent. Among 567,169 participants who completed the questionnaires, we excluded participants who were proxy responders, had duplicate records, decided to withdraw, moved or died before entry, and had null person-years of follow-up or extreme total energy intake (< 800 or > 4200 kcal/day for men and < 600 or > 3500 kcal/day for women [12]). Finally, 521,120 persons were selected (Additional file 1: Fig. S1), which was approved by the Institutional Review Board of the National Cancer Institute.

Assessment of diet and cooking oils/fats

Dietary intake was assessed at baseline using a 124-item self-administered food frequency questionnaire (FFQ) developed as the diet history questionnaire (DHQ) and validated by National Cancer Institute [13]. The frequency and portion sizes of food consumption were recorded during the past year. Some questions were asked to collect the frequency of oils/fats used in cooking and added after cooking, such as “How often was oil, butter, or margarine used to fry or saute the vegetables, eggs, or meat you ate?” followed by options from “never” to “4 or more times per day” and “When you ate each of the foods listed in this table, how often was butter or margarine added after cooking or at table?” followed by options from “almost never or never” to “more than half the time.” Participants were also asked to select the types of oils/fats they regularly used, including butter, margarine, lard, corn oil, canola oil, and olive oil. Portion sizes for individual oils/fats were estimated based on the 1994–1996 USDA Continuing Survey of Food Intakes by Individuals (CSFII) [14] and intakes of cooking oils/fats were then calculated by multiplying the frequency of consumption with the corresponding portion size derived from CSFII. Cooking oils/fats included in foods were also accounted for by asking questions such as the frequency of consuming butter or margarine on bread or rolls. Three typical solid fats (butter, margarine, and lard) and 3 commonly consumed vegetable oils (olive oil, canola oil, and corn oil) were calculated and analyzed in our study. The Healthy Eating Index (HEI)-2015 score was established to assess the adherence to an overall healthy dietary pattern according to US Dietary Guideline 2015–2020. The HEI-2015 is composed of 13 components, including total fruit, whole fruit, total vegetables, greens and beans, whole grains, dairy, total protein foods, seafood and plant proteins, fatty acids, refined grains, sodium, added sugars, and saturated fats [15]. Scores are assigned to each component by comparing the density (the amount of dietary component per 1000 kcal) to the relevant standards [15].

Mortality

All the participants were followed for address changes via the US Postal Service National Change of Address database, responses to study-related mailings, and direct notifications from cohort members. Deaths were identified by annual linkage to the Social Security Administration Death Master File and confirmed by follow-up searches of the National Death Index Plus. The International Classification of Diseases 9th and 10th Revision codes were used to classify death causes (Additional file 1: Table S1). Follow-up was calculated from the return date of the baseline questionnaire to the time of death or the end of follow-up (31 December 2011), whichever occurred earlier. The complete follow-up rate for mortality exceeds 99% in this cohort study.

Statistical analysis

Intakes of individual cooking oils/fats were expressed as the functions of energy density (g 2000 kcal−1 day−1) using the nutrient density method [16]. We used Cox proportional hazards regression models to estimate hazard ratios (HRs) and 95% confidence intervals (CIs) for mortality. Model 1 was adjusted for age and sex at baseline. Model 2 was further adjusted for race, marital status, education, household income, body mass index (BMI), alcohol, smoking, vigorous physical activity, usual activity at work, perceived health condition, and history of cancer, heart disease, stroke, and diabetes. Our final multivariate model 3 was additionally adjusted for HEI-2015 (with no fat components), total energy intake, and consumption of remaining oils where appropriate. Tests for linear trend were performed by assigning median values to corresponding categories of intake and modeling the values as continuous variables. We also estimated the associations of hypothetical substitution of 1 tablespoon/day olive oil, canola oil, or corn oil for the equivalent amounts of butter and margarine with mortality by simultaneously including individual cooking oils as continuous variables and total cooking oil intake in the same multivariable model (substitution model), which also contained total energy intake, HEI, and other non-dietary covariates. Total oil/fat intake was held constant in this model. By leaving butter or margarine out of this model, regression coefficients of the remaining oils bear the interpretations as the theoretical effects of substituting one of these oils for the equivalent amounts of butter or margarine while holding other oils unchanged. A fixed 1-tablespoon/day increase corresponds to an increment of approximately 14 g/day butter/margarine or 8 g/day vegetable oils [9, 17].

We further separately analyzed the associations for stick margarine and tub/soft margarine, respectively. Subgroup analyses were also conducted according to important potential effect modifiers and P values for interactions were tested by the likelihood-ratio test. Sensitivity analyses were performed by excluding participants who had extreme BMI (< 18.5 or > 40 kg/m2); using the propensity-score adjustment [18] to further control for potential residual confounding from measured variables; further adjusting for history of hypertension and hypercholesteremia and the use of aspirin and multivitamins; excluding the initial 4 years of follow-up; excluding those who had CVD, cancer, or diabetes at baseline; or ending up the follow-up at the year 2004 (midpoint, 8 years of follow-up). We also tested whether the associations were affected by the use of cholesterol-lowering medications among persons who provided this information in the resurvey (n = 293,918).

Statistical analyses were performed with SAS version 9.4 (SAS Institute, Cary, NC, USA). Tests were two-sided and the significance was defined as P < 0.05.

Results

Population characteristics and cooking oil/fat consumption

During an average of 16 years of follow-up (7,307,097 person-years), 129,328 individuals died, including 85,037 men and 44,291 women. At baseline, participants with higher butter consumption were less likely to be married, have prevalent heart disease, and use aspirin, and they had lower protein intake and lower HEI scores. Participants with higher margarine consumption were more likely to have a higher BMI, use aspirin, and have heart disease, stroke, and diabetes, and they had lower household income and lower alcohol intake. The median intakes in the highest tertile among consumers were 13.7 g 2000 kcal− 1 day−1 for butter and 20.6 g 2000 kcal−1 day−1 for margarine, respectively (Table 1). Participant characteristics according to corn oil, canola oil, and olive oil consumption are shown in Additional file 1: Table S2. The Spearman correlations between individual cooking oil/fat consumption are presented in Additional file 1: Table S3.

Table 1 Baseline characteristics of participants from the NIH-AARP Diet and Health Study according to butter and margarine consumption

All-cause mortality

Butter and margarine consumption was strongly associated with higher all-cause mortality in all multivariable-adjusted models. In contrast, intakes of canola oil and olive oil were both inversely associated with all-cause mortality. Corn oil consumption was related to higher all-cause mortality after adjusting for age and sex, but the association became non-significant after adjusting further covariates (Table 2). Compared to non-consumers, the multivariable HRs of all-cause mortality in the highest categories were 1.09 (95% CI, 1.07–1.11) for butter (P-trend < 0.001), 1.07 (95% CI, 1.05–1.09) for margarine (P-trend < 0.001), 0.97 (95% CI, 0.95–0.99) for canola oil (P-trend < 0.001), and 0.96 (95% CI, 0.95–0.98) for olive oil (P-trend < 0.001) (Table 2). Every 1-tablespoon/day increment of butter or margarine consumption was related to 7% and 4% higher all-cause mortality, respectively. In contrast, each 1-tablespoon/day increment of canola oil or olive oil consumption was associated with 2% and 3% of reductions in all-cause mortality, respectively (Fig. 1a).

Table 2 HRs (95% CIs) of all-cause mortality according to cooking oil/fat consumption
Fig. 1
figure 1

Multivariable-adjusted hazard ratios of total and cardiometabolic mortality for 1-tablespoon/day increment in cooking oil/fat consumption. Forest plots show the multivariable HRs of total (a) and cardiometabolic (b) mortality associated with 1-tablespoon/day increment in butter, margarine, corn oil, canola oil, and olive oil consumption. HRs were adjusted for age, sex, BMI, race, education, marital status, household income, smoking, alcohol, vigorous physical activity, usual activity at work, perceived health condition, history of heart disease, stroke, diabetes, and cancer at baseline, Healthy Eating Index-2015, total energy intake, and consumption of remaining oils where appropriate (butter, margarine, lard, corn oil, canola oil, olive oil, and other vegetable oils). Horizontal lines represent 95% CIs

Cardiometabolic mortality

The consumption of butter and margarine was positively associated with CVD mortality after multivariate adjustment, whereas olive oil intake was inversely associated with CVD mortality (Table 3). Compared with non-consumers, participants in the highest tertile of olive oil consumption had 5% (HR = 0.95, 95% CI 0.92–0.99; P-trend = 0.001) lower CVD mortality, but those in the highest tertiles of butter and margarine consumption had 8% (HR = 1.08, 95% CI 1.05–1.12; P-trend< 0.001) and 10% (HR = 1.10, 95% CI 1.06–1.14; P-trend< 0.001) higher CVD mortality, respectively. Canola oil consumption was marginally associated with lower CVD mortality (P-trend = 0.052), while corn oil intake was not related to CVD mortality. Similar associations were also observed for heart disease mortality (Additional file 1: Table S4). Besides, borderline trends toward lower and higher stroke mortality were observed for corn oil (P-trend = 0.061) and butter (P-trend = 0.059) consumption, respectively. For diabetes mortality, we detected positive associations for butter and margarine consumption. Compared with non-consumers, participants in the highest tertiles of butter and margarine consumption had 18% (HR = 1.18, 95% CI 1.06–1.32; P-trend = 0.0041) and 12% (HR = 1.12, 95% CI 1.00–1.26; P-trend = 0.047) higher diabetes mortality, respectively. In contrast, olive oil consumption was inversely related to diabetes mortality (HR comparing the highest tertile with non-consumers = 0.87, 95% CI 0.77–0.99; P-trend = 0.019). Overall, each 1-tablespoon/day increment of butter or margarine consumption was associated with 8% and 6% higher cardiometabolic mortality, respectively, while every 1-tablespoon/day increment of olive oil consumption was related to 4% decreased cardiometabolic mortality (Fig. 1b). Restricted-cubic-spline regression yielded similar results (Fig. 2a–e).

Table 3 HRs (95% CIs) of CVD and diabetes mortality according to cooking oil/fat consumption
Fig. 2
figure 2

Cubic spline curves for the association between individual cooking oil intakes and cardiometabolic mortality. Hazard ratios are based on Cox proportional hazards models adjusted for age, sex, BMI, race, education, marital status, household income, smoking, alcohol, vigorous physical activity, usual activity at work, perceived health condition, history of heart disease, stroke, diabetes, and cancer at baseline, Healthy Eating Index-2015, total energy intake, and consumption of remaining oils where appropriate (butter, margarine, lard, corn oil, canola oil, olive oil, and other vegetable oils)

Other non-cardiometabolic mortality

For other cause-specific mortality, butter consumption was associated with higher mortality from cancer, respiratory disease (RD), kidney disease, and chronic liver disease. Intake of margarine was related to higher RD and kidney disease mortality. Nonetheless, canola oil intake turned to be inversely associated with RD and infection mortality, while higher olive oil consumption was related to lower mortality from RD and Alzheimer’s disease (AD). We observed a borderline trend toward higher consumption of corn oil in relation to lower AD mortality (P-trend = 0.061) (Additional file 1: Table S5). Compared with non-consumers, lard consumers had 4% and 13% higher all-cause and RD mortality, respectively (Additional file 1: Table S6).

Substitution for butter and margarine

Substituting 1 tablespoon/day (8 g/day) corn oil, canola oil, or olive oil for equivalent amounts of butter was associated with 5%, 6%, and 7% lower all-cause mortality, respectively. Likewise, substituting 1 tablespoon/day corn oil, canola oil, or olive oil for equivalent amounts of butter was associated with 5%, 6%, and 8% lower cardiometabolic mortality, respectively. For the cause-specific mortality by the replacement of butter, each 1-tablespoon/day increment of canola oil was associated with 7%, 5%, and 11% of reductions in CVD, cancer, and RD mortality, respectively, and each 1-tablespoon/day increment of olive oil was related to 7%, 16%, 21%, and 15% lower mortality from CVD, RD, AD, and diabetes, respectively. Besides, we detected 5% lower cancer mortality when replacing 1 tablespoon/day butter with corn oil (Fig. 3 and Additional file 1: Table S7).

Fig. 3
figure 3

Multivariable-adjusted hazard ratios of total and cause-specific mortality by substitution of 1 tablespoon/day individual cooking oils/fats for butter and margarine. Forest plots show the multivariable HRs of total and cause-specific mortality associated with substitution of 1 tablespoon/day individual cooking oils/fats for equivalent amounts of butter and margarine. HRs were adjusted for age, sex, BMI, race, education, marital status, household income, smoking, alcohol, vigorous physical activity, usual activity at work, perceived health condition, history of heart disease, stroke, diabetes, and cancer at baseline, Healthy Eating Index-2015, and total energy intake. Horizontal lines represent 95% CIs

When replacing margarine, each 1-tablespoon/day increment of corn oil, canola oil, or olive oil was associated with 3%, 5%, and 5% lower all-cause mortality, respectively. Similarly, each 1-tablespoon/day increment of corn oil, canola oil, or olive oil was associated with 4%, 6%, and 7% lower cardiometabolic mortality, respectively. Regarding the cause-specific mortality, replacing with 1 tablespoon/day canola oil was associated with lower CVD mortality (6%) and RD mortality (9%), respectively, and replacing with olive oil was related to lower mortality from CVD (6%), RD (14%), AD (25%), and diabetes (13%), respectively. Each 1-tablespoon/day increment of corn oil was related to 24% lower AD mortality (Fig. 3 and Additional file 1: Table S7).

Subgroup analyses

In secondary analyses for stick and tub/soft margarine consumption, most of the associations were similar except a significant association of stick but not tub/soft margarine consumption with higher AD mortality (Additional file 1: Table S8). Subgroup analyses showed significant modifications for the associations with all-cause mortality stratified by sex, baseline BMI, smoking status, alcohol drinking, income level, and HEI score (Additional file 1: Table S9). Although positive associations of butter and margarine consumption with all-cause mortality persisted in all the subgroups, the associations for butter consumption were stronger among men than women (P-interaction = 0.022) and among non-obese participants than obese participants (P-interaction< 0.001), while the associations for margarine consumption were more pronounced among non-smokers or former smokers (P-interaction = 0.004) and those with higher income level (P-interaction< 0.001) and higher HEI score (P-interaction = 0.001). Notably, the inverse association of olive oil consumption with all-cause mortality was restricted to alcohol drinkers (P-interaction = 0.003).

Sensitivity analyses

Both significant and non-significant associations did not change materially after exclusion of participants with extreme BMIs, or further adjustment for a propensity score, history of hypertension and hypercholesteremia, or the use of multivitamins and aspirin (Additional file 1: Tables S10-S13). We also detected largely similar associations when we further adjusted for the use of cholesterol-lowering medications; excluded those with CVD, cancer, or diabetes at baseline and those with the first 4 years of follow-up; or censored participants at 8-year follow-up (Additional file 1: Tables S14-S17).

Discussion

In this large prospective cohort, intakes of butter and margarine were associated with higher all-cause, cardiometabolic, and other major cause-specific mortality, whereas intakes of canola oil and olive oil were related to lower mortality. Substituting corn oil, canola oil, or olive oil for butter and margarine was associated with lower all-cause, cardiometabolic, cancer, RD, and AD mortality.

Current dietary recommendations on butter consumption largely depend on the assumed negative effect of high SFA content (> 65%) [19] in relation to higher all-cause mortality and CVD incidence [20, 21] regardless of beneficial ingredients such as vitamins, PUFAs, MUFAs, and ruminant trans-fat. However, a previous multi-center nationwide study reported no significant association of butter consumption with all-cause mortality and CVD incidence [22]. Nonetheless, evidence from intervention trials indicated that the butter-enriched diet elevated total and LDL cholesterol compared with the diet rich in vegetable oils or coconut oil [5, 6]. In the current US prospective study, we observed strong associations with higher total and CVD mortality. These data emphasize the reductions in butter consumption among the US population for the management of cardiovascular health. Interestingly, a previous meta-analysis of 11 studies revealed a weak inverse relationship of butter with the risk of diabetes [9], whereas we found participants consuming > 8.6 g 2000 kcal−1 day−1 had an 18% higher risk of diabetes mortality compared with non-consumers in our study. In accordance with our finding, replacing butter with olive oil was related to a lower risk of diabetes among US women [17]. Our finding of butter consumption in relation to higher cancer mortality was supported by previous studies showing higher incidence of breast cancer and non-Hodgkin lymphoma [23, 24], whereas the other positive associations were first reported. Together, our documented detrimental effects of butter on mortality supported the current US dietary recommendation on decreasing butter intake and highlighted the need for further mechanistic studies on the metabolic effects of butter. In addition, we found higher all-cause and RD mortality among lard consumers. Despite the absent evidence in humans, lard has been frequently used as a high-fat diet to induce obese phenotype and metabolic dysfunction in rodents [25]. More investigations are warranted given the low consumption of lard in the current study.

Margarine contains trans-fat, a well-documented risk factor for arterial calcification and coronary heart disease [26], and has a negative impact on plasma lipid profiles in both healthy individuals and patients with hypercholesterolemia [27]. Our findings were consistent with a recent meta-analysis, showing a positive association of trans-fat with all-cause and CVD mortality [28]. The observed association of margarine intake with modestly higher diabetes mortality was in line with a European multi-center study [29]. In addition, higher incidence of asthma onset was contributed by the intake of margarine [30], supporting our finding of elevated RD mortality. Compared with tub/soft margarine, our secondary analysis showed that stick margarine consumption turned to be much stronger for its positive association with AD mortality, which could be explained by higher trans-fat content (15–21%) [31] and supported by previous evidence suggesting a negative effect of trans-fat on dementia [32]. Taken together, our results suggest the importance of restricting intake of trans-fat containing margarines to decrease the incidence of cardiometabolic diseases.

We observed overall neutral associations of corn oil consumption with mortality only except a marginal inverse association with stroke mortality. Corn oil could ameliorate plasma atherogenic lipids among participants with elevated cholesterol [33]. However, substituting vegetable oils rich in linoleic acid for SFAs [34] was not associated with lower all-cause or CVD mortality [35]. Unlike corn oil, canola oil as a cardioprotective contributor is rich in MUFAs and α-linolenic acid [36]. Canola oil supplementation could ameliorate overall blood lipid profiles and improve glycemic control and inflammation [6]. Consistently, we showed a protective relationship between canola oil consumption and all-cause and heart disease mortality. The observed lower RD and infection mortality might be ascribed to immunomodulatory and anti-inflammatory effects of MUFAs [37].

The cardioprotective association of olive oil has been supported by the Spanish European Prospective into Cancer and Nutrition study [38], the Nurses’ Health Study and the Health Professionals Follow-up Study [39], and the PREvención con DIeta MEDiterránea study of patients at high CVD incidence [40]. Consistent with our finding, a previous US women study also showed a lower risk of diabetes with higher olive oil intake [17]. Mediterranean diet supplemented with olive oil seemed to be effective in reducing the risk of developing diabetes among participants at high cardiovascular risk [41]. A recent meta-analysis also demonstrated that olive oil supplementation could lower HbA1c and fasting plasma glucose levels among diabetic patients [8]. Besides, the association with lower AD mortality might be due to the neuroprotective phenols and oleocanthal in olive oil. Collectively, these data underscored the health benefits of olive oil and provided evidence on habitual olive oil consumption as a key contributor to the healthy Mediterranean diet [40]. Increasing the consumption of olive oil may confer health benefits on cardiometabolic health and reduce mortality.

Subgroup analyses demonstrated stronger associations of butter consumption with higher mortality among men than women, which could be due to sex-dependent fatty acid metabolism [42], and higher CVD mortality rates and lower life expectancy in men than women. The interaction between olive oil and alcohol drinking status on mortality might be due to the correlation between olive oil consumption and the Mediterranean pattern which included wine. Future research is warranted to elucidate the observed interactions between smoking, alcohol drinking, BMI, and HEI score and cooking oils/fats.

Strengths of the current study include the large population size, long-term duration with a high follow-up rate (> 99%), and a large number of deaths from various causes. We excluded participants with chronic diseases at baseline or initial 4 years of follow-up to further reduce the possibility of reverse causality and found similar results, indicating the robustness of our findings. The limitations should also be noted. First, the observed associations might be partly due to residual confounding despite comprehensive adjustment for well-known risk factors. Second, although validated, our FFQ might still produce measurement errors. Given the prospective study design, any mismeasurement in cooking oil/fat consumption would likely be random for mortality, resulting in conservative associations. Third, the overall intake level of individual vegetable oils was low in this population with a relatively narrow intake range. Nonetheless, we still detected significant protective associations of these vegetable oils, especially when substituting for solid fats. Fourth, due to the lack of measurement data, we could not further analyze the associations between fats/oils with different cooking methods and mortality, which could make our results more informative. In addition, dietary intakes were only assessed at baseline, and potential dietary changes could occur during the long-term follow-up. Nonetheless, this might not appreciably change our results because we also observed similar associations when censoring participants at a shorter duration (8 years) of follow-up. Moreover, the potential changes in the constituents of specific oils could not be captured. Typically, margarines contain high amount of trans-fat in 1990s. However, trans-fat was phased out in the last decade. In recent years, soft margarines have contained no trans-fat and only hard and tub may still contain. This probably would have weakened the observed positive associations of margarine consumption. Finally, a causal relationship may not be established due to the observational nature.

Conclusion

In summary, consumption of butter and margarine was associated with higher all-cause and cardiometabolic mortality. Intakes of canola oil and olive oil were associated with lower total mortality and corn oil had a neutral association with mortality. From the standpoint of public health, intakes of butter and margarine may be limited while olive oil consumption may be recommended to lower deaths from cardiometabolic diseases. Replacing butter and margarine with corn oil, canola oil, or olive oil may confer health benefits on cardiometabolic health and reduce mortality. Taken together, current dietary recommendations should continue to highlight shifting the intake from solid fats, including butter and margarine, to non-hydrogenated vegetable oils, such as corn oil, canola oil, and olive oil, for the prevention of cardiometabolic diseases and premature deaths.

Availability of data and materials

Because of the sensitive nature of the data collected for this study, requests to access the dataset from qualified researchers trained in human subject confidentiality protocols may be sent to the National Cancer Institute Division of Cancer Epidemiology & Genetics to Linda M. Liao (e-mail: liaolm@mail.nih.gov).

Abbreviations

AARP:

Formerly known as the American Association for Retired Persons

AD:

Alzheimer’s disease

BMI:

Body mass index

CI:

Confidence interval

CVD:

Cardiovascular disease

FFQ:

Food frequency questionnaire

HEI:

Healthy Eating Index

HR:

Hazard ratio

MUFA:

Monounsaturated fatty acid

PUFA:

Polyunsaturated fatty acid

RD:

Respiratory disease

SFA:

Saturated fatty acid

References

  1. Consumption of vegetable oils worldwide from 2013/14 to 2018/2019, by oil type (in million metric tons). https://www.statista.com/statistics/263937/vegetable-oils-global-consumption/. Accessed 1 Sept 2020.

  2. Forouhi NG, Krauss RM, Taubes G, Willett W. Dietary fat and cardiometabolic health: evidence, controversies, and consensus for guidance. BMJ. 2018;361:k2139.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Dietary Guidelines Advisory Committee. Scientific Report of the 2015 Dietary Guidelines Advisory Committee (Advisory Report). In: U.S. Department of Health and Human Services (HHS) and the U.S. 3 Department of Agriculture (USDA), editor. 2015.

  4. Mozaffarian D, Ludwig DS. Dietary guidelines in the 21st century—a time for food. JAMA. 2010;304(6):681–2. https://doi.org/10.1001/jama.2010.1116.

    Article  CAS  PubMed  Google Scholar 

  5. Schwingshackl L, Bogensberger B, Benčič A, Knüppel S, Boeing H, Hoffmann G. Effects of oils and solid fats on blood lipids: a systematic review and network meta-analysis. J Lipid Res. 2018;59(9):1771–82. https://doi.org/10.1194/jlr.P085522.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Jenkins DJA, Kendall CWC, Vuksan V, Faulkner D, Augustin LSA, Mitchell S, Ireland C, Srichaikul K, Mirrahimi A, Chiavaroli L, Blanco Mejia S, Nishi S, Sahye-Pudaruth S, Patel D, Bashyam B, Vidgen E, de Souza RJ, Sievenpiper JL, Coveney J, Josse RG, Leiter LA. Effect of lowering the glycemic load with canola oil on glycemic control and cardiovascular risk factors: a randomized controlled trial. Diabetes Care. 2014;37(7):1806–14. https://doi.org/10.2337/dc13-2990.

    Article  CAS  PubMed  Google Scholar 

  7. Engel S, Tholstrup T. Butter increased total and LDL cholesterol compared with olive oil but resulted in higher HDL cholesterol compared with a habitual diet. Am J Clin Nutr. 2015;102(2):309–15. https://doi.org/10.3945/ajcn.115.112227.

    Article  CAS  PubMed  Google Scholar 

  8. Schwingshackl L, Lampousi AM, Portillo MP, Romaguera D, Hoffmann G, Boeing H. Olive oil in the prevention and management of type 2 diabetes mellitus: a systematic review and meta-analysis of cohort studies and intervention trials. Nutr Diab. 2017;7(4):e262. https://doi.org/10.1038/nutd.2017.12.

    Article  CAS  Google Scholar 

  9. Pimpin L, Wu JHY, Haskelberg H, Del Gobbo L, Mozaffarian D. Is butter back? A systematic review and meta-analysis of butter consumption and risk of cardiovascular disease, diabetes, and total mortality. Plos One. 2016;11(6):e0158118. https://doi.org/10.1371/journal.pone.0158118.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Guasch-Ferré M, Hu FB, Martínez-González MA, Fitó M, Bulló M, Estruch R, Ros E, Corella D, Recondo J, Gómez-Gracia E, Fiol M, Lapetra J, Serra-Majem L, Muñoz MA, Pintó X, Lamuela-Raventós RM, Basora J, Buil-Cosiales P, Sorlí JV, Ruiz-Gutiérrez V, Martínez JA, Salas-Salvadó J. Olive oil intake and risk of cardiovascular disease and mortality in the PREDIMED Study. BMC Med. 2014;12(1):78. https://doi.org/10.1186/1741-7015-12-78.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Guasch-Ferre M, Zong G, Willett Walter C, Zock P, Wanders Anne J, Hu Frank B, et al. Associations of monounsaturated fatty acids from plant and animal sources with total and cause-specific mortality in two US prospective cohort studies. Circ Res. 2019;124(8):1266–75. https://doi.org/10.1161/CIRCRESAHA.118.313996.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. WCW. Nutritional epidemiology. New York: Oxford University Press; 1998.

    Google Scholar 

  13. Midthune D, Schatzkin A, Subar AF, Thompson FE, Freedman LS, Carroll RJ, Shumakovich MA, Kipnis V. Validating an FFQ for intake of episodically consumed foods: application to the National Institutes of Health–AARP Diet and Health Study. Public Health Nutr. 2011;14(7):1212–21. https://doi.org/10.1017/S1368980011000632.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Subar AF, Midthune D, Kulldorff M, Brown CC, Thompson FE, Kipnis V, Schatzkin A. Evaluation of alternative approaches to assign nutrient values to food groups in food frequency questionnaires. Am J Epidemiol. 2000;152(3):279–86. https://doi.org/10.1093/aje/152.3.279.

    Article  CAS  PubMed  Google Scholar 

  15. Krebs-Smith SM, Pannucci TE, Subar AF, Kirkpatrick SI, Lerman JL, Tooze JA, et al. Update of the healthy eating index: HEI-2015. J Acad Nutr Diet. 2018;118(9):1591–602. https://doi.org/10.1016/j.jand.2018.05.021.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Willett WC. Chapter 5: food frequency methods; chapter 6: reproducibility and validity of food-frequency questionnaires; 2012.

    Google Scholar 

  17. Guasch-Ferré M, Hruby A, Salas-Salvadó J, Martínez-González MA, Sun Q, Willett WC, Hu FB. Olive oil consumption and risk of type 2 diabetes in US women. Am J Clin Nutr. 2015;102(2):479–86. https://doi.org/10.3945/ajcn.115.112029.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. D’Agostino RB. Propensity score methods for bias reduction in the comparison of a treatment to a non-randomized control group. Stat Med. 1998;17(19):2265–81. https://doi.org/10.1002/(SICI)1097-0258(19981015)17:19<2265::AID-SIM918>3.0.CO;2-B.

    Article  PubMed  Google Scholar 

  19. 2015–2020 Dietary Guidelines for Americans. 8th Edition: U.S. Department of Health and Human Services and U.S. Department of Agriculture December. 2015.

  20. Wang DD, Li Y, Chiuve SE, Stampfer MJ, Manson JE, Rimm EB, et al. Association of specific dietary fats with total and cause-specific mortality. JAMA Intern Med. 2016;176(8):1134–45. https://doi.org/10.1001/jamainternmed.2016.2417.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Xu J, Eilat-Adar S, Loria C, Goldbourt U, Howard BV, Fabsitz RR, Zephier EM, Mattil C, Lee ET. Dietary fat intake and risk of coronary heart disease: the Strong Heart Study. Am J Clin Nutr. 2006;84(4):894–902. https://doi.org/10.1093/ajcn/84.4.894.

    Article  CAS  PubMed  Google Scholar 

  22. Dehghan M, Mente A, Rangarajan S, Sheridan P, Mohan V, Iqbal R, Gupta R, Lear S, Wentzel-Viljoen E, Avezum A, Lopez-Jaramillo P, Mony P, Varma RP, Kumar R, Chifamba J, Alhabib KF, Mohammadifard N, Oguz A, Lanas F, Rozanska D, Bostrom KB, Yusoff K, Tsolkile LP, Dans A, Yusufali AH, Orlandini A, Poirier P, Khatib R, Hu B, Wei L, Yin L, Deeraili A, Yeates K, Yusuf R, Ismail N, Mozaffarian D, Teo K, Anand SS, Yusuf S. Association of dairy intake with cardiovascular disease and mortality in 21 countries from five continents (PURE): a prospective cohort study. Lancet. 2018;392(10161):2288–97. https://doi.org/10.1016/S0140-6736(18)31812-9.

    Article  PubMed  Google Scholar 

  23. Voorrips LE, Brants HAM, Kardinaal AFM, Hiddink GJ, van den Brandt PA, Goldbohm RA. Intake of conjugated linoleic acid, fat, and other fatty acids in relation to postmenopausal breast cancer: the Netherlands Cohort Study on Diet and Cancer. Am J Clin Nutr. 2002;76(4):873–82. https://doi.org/10.1093/ajcn/76.4.873.

    Article  CAS  PubMed  Google Scholar 

  24. Wang J, Li X, Zhang D. Dairy product consumption and risk of non-Hodgkin lymphoma: a meta-analysis. Nutrients. 2016;8(3):120. https://doi.org/10.3390/nu8030120.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Putti R, Migliaccio V, Sica R, Lionetti L. Skeletal muscle mitochondrial bioenergetics and morphology in high fat diet induced obesity and insulin resistance: focus on dietary fat source. Front Physiol. 2016;6:426.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Mozaffarian D, Katan MB, Ascherio A, Stampfer MJ, Willett WC. Trans fatty acids and cardiovascular disease. N Engl J Med. 2006;354(15):1601–13. https://doi.org/10.1056/NEJMra054035.

    Article  CAS  PubMed  Google Scholar 

  27. Lichtenstein AH, Ausman LM, Jalbert SM, Schaefer EJ. Effects of different forms of dietary hydrogenated fats on serum lipoprotein cholesterol levels. N Engl J Med. 1999;340(25):1933–40. https://doi.org/10.1056/NEJM199906243402501.

    Article  CAS  PubMed  Google Scholar 

  28. de Souza RJ, Mente A, Maroleanu A, Cozma AI, Ha V, Kishibe T, et al. Intake of saturated and trans unsaturated fatty acids and risk of all cause mortality, cardiovascular disease, and type 2 diabetes: systematic review and meta-analysis of observational studies. BMJ. 2015;351:h3978.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Buijsse B, Boeing H, Drogan D, Schulze MB, Feskens EJ, Amiano P, Barricarte A, Clavel-Chapelon F, de Lauzon-Guillain B, Fagherazzi G, Fonseca-Nunes A, Franks PW, Huerta JM, Jakobsen MU, Kaaks R, Key TJ, Khaw KT, Masala G, Moskal A, Nilsson PM, Overvad K, Pala V, Panico S, Redondo ML, Ricceri F, Rolandsson O, Sánchez MJ, Sluijs I, Spijkerman AM, Tjonneland A, Tumino R, van der AD, van der Schouw Y, Langenberg C, Sharp SJ, Forouhi NG, Riboli E, Wareham NJ, InterAct Consortium. Consumption of fatty foods and incident type 2 diabetes in populations from eight European countries. Eur J Clin Nutr. 2014;69(4):455–61. https://doi.org/10.1038/ejcn.2014.249.

    Article  CAS  PubMed  Google Scholar 

  30. Bolte G, Winkler G, Hölscher B, Thefeld W, Weiland SK, Heinrich J. Margarine consumption, asthma, and allergy in young adults: results of the German National Health Survey 1998. Ann Epidemiol. 2005;15(3):207–13. https://doi.org/10.1016/j.annepidem.2004.04.004.

    Article  PubMed  Google Scholar 

  31. Liu Q, Rossouw JE, Roberts MB, Liu S, Johnson KC, Shikany JM, Manson JAE, Tinker LF, Eaton CB. Theoretical effects of substituting butter with margarine on risk of cardiovascular disease. Epidemiology. 2017;28(1):145–56. https://doi.org/10.1097/EDE.0000000000000557.

    Article  PubMed  PubMed Central  Google Scholar 

  32. Barnard ND, Bunner AE, Agarwal U. Saturated and trans fats and dementia: a systematic review. Neurobiol Aging. 2014;35(Suppl 2):S65–73. https://doi.org/10.1016/j.neurobiolaging.2014.02.030.

    Article  CAS  PubMed  Google Scholar 

  33. Maki KC, Lawless AL, Kelley KM, Kaden VN, Geiger CJ, Dicklin MR. Corn oil improves the plasma lipoprotein lipid profile compared with extra-virgin olive oil consumption in men and women with elevated cholesterol: results from a randomized controlled feeding trial. J Clin Lipidol. 2015;9(1):49–57. https://doi.org/10.1016/j.jacl.2014.10.006.

    Article  PubMed  Google Scholar 

  34. Sakata Y, Shimokawa H. Saturated fatty acid intake and cardiovascular risk. Eur Heart J. 2013;34(16):1178–80. https://doi.org/10.1093/eurheartj/eht057.

    Article  PubMed  Google Scholar 

  35. Ramsden CE, Zamora D, Majchrzak-Hong S, Faurot KR, Broste SK, Frantz RP, et al. Re-evaluation of the traditional diet-heart hypothesis: analysis of recovered data from Minnesota Coronary Experiment (1968-73). BMJ. 2016;353:i1246.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Lin L, Allemekinders H, Dansby A, Campbell L, Durance-Tod S, Berger A, Jones PJH. Evidence of health benefits of canola oil. Nutr Rev. 2013;71(6):370–85. https://doi.org/10.1111/nure.12033.

    Article  PubMed  Google Scholar 

  37. Wood LG, Scott HA, Garg ML, Gibson PG. Innate immune mechanisms linking non-esterified fatty acids and respiratory disease. Prog Lipid Res. 2009;48(1):27–43. https://doi.org/10.1016/j.plipres.2008.10.001.

    Article  CAS  PubMed  Google Scholar 

  38. Buckland G, Mayén AL, Agudo A, Travier N, Navarro C, Huerta JM, Chirlaque MD, Barricarte A, Ardanaz E, Moreno-Iribas C, Marin P, Quirós JR, Redondo ML, Amiano P, Dorronsoro M, Arriola L, Molina E, Sanchez MJ, Gonzalez CA. Olive oil intake and mortality within the Spanish population (EPIC-Spain). Am J Clin Nutr. 2012;96(1):142–9. https://doi.org/10.3945/ajcn.111.024216.

    Article  CAS  PubMed  Google Scholar 

  39. Guasch-Ferré M, Liu G, Li Y, Sampson L, Manson JE, Salas-Salvadó J, et al. Olive oil consumption and cardiovascular risk in U.S. adults. J Am Coll Cardiol. 2020;75(15):1729–39. https://doi.org/10.1016/j.jacc.2020.02.036.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Salas-Salvadó J, Bulló M, Babio N, Martínez-González MÁ, Ibarrola-Jurado N, Basora J, et al. Reduction in the incidence of type 2 diabetes with the Mediterranean diet: results of the PREDIMED-Reus nutrition intervention randomized trial. Diabetes Care. 2011;34(1):14–9. https://doi.org/10.2337/dc10-1288.

    Article  PubMed  Google Scholar 

  41. Estruch R, Ros E, Salas-Salvadó J, Covas M-I, Corella D, Arós F, Gómez-Gracia E, Ruiz-Gutiérrez V, Fiol M, Lapetra J, Lamuela-Raventos RM, Serra-Majem L, Pintó X, Basora J, Muñoz MA, Sorlí JV, Martínez JA, Fitó M, Gea A, Hernán MA, Martínez-González MA. Primary prevention of cardiovascular disease with a Mediterranean diet supplemented with extra-virgin olive oil or nuts. N Engl J Med. 2018;378(25):e34. https://doi.org/10.1056/NEJMoa1800389.

    Article  CAS  PubMed  Google Scholar 

  42. Decsi T, Kennedy K. Sex-specific differences in essential fatty acid metabolism. Am J Clin Nutr. 2011;94(6 Suppl):1914S–9S. https://doi.org/10.3945/ajcn.110.000893.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

We are indebted to the participants in the NIH-AARP Diet and Health Study for their outstanding cooperation. We also thank Sigurd Hermansen and Kerry Grace Morrissey from Westat for study outcome ascertainment and management and Leslie Carroll at Information Management Services for data support and analysis. Cancer incidence data from the Atlanta metropolitan area were collected by the Georgia Center for Cancer Statistics, Department of Epidemiology, Rollins School of Public Health, Emory University, Atlanta, Georgia. Cancer incidence data from California were collected by the California Cancer Registry, California Department of Public Health’s Cancer Surveillance and Research Branch, Sacramento, California. Cancer incidence data from the Detroit metropolitan area were collected by the Michigan Cancer Surveillance Program, Community Health Administration, Lansing, Michigan. The Florida cancer incidence data used in this report were collected by the Florida Cancer Data System (Miami, Florida) under contract with the Florida Department of Health, Tallahassee, Florida. The views expressed herein are solely those of the authors and do not necessarily reflect those of the FCDC or FDOH. Cancer incidence data from Louisiana were collected by the Louisiana Tumor Registry, Louisiana State University Health Sciences Center School of Public Health, New Orleans, Louisiana. Cancer incidence data from New Jersey were collected by the New Jersey State Cancer Registry, The Rutgers Cancer Institute of New Jersey, New Brunswick, New Jersey. Cancer incidence data from North Carolina were collected by the North Carolina Central Cancer Registry, Raleigh, North Carolina. Cancer incidence data from Pennsylvania were supplied by the Division of Health Statistics and Research, Pennsylvania Department of Health, Harrisburg, Pennsylvania. The Pennsylvania Department of Health specifically disclaims responsibility for any analyses, interpretations, or conclusions. Cancer incidence data from Arizona were collected by the Arizona Cancer Registry, Division of Public Health Services, Arizona Department of Health Services, Phoenix, Arizona. Cancer incidence data from Texas were collected by the Texas Cancer Registry, Cancer Epidemiology and Surveillance Branch, Texas Department of State Health Services, Austin, Texas. Cancer incidence data from Nevada were collected by the Nevada Central Cancer Registry, Division of Public and Behavioral Health, State of Nevada Department of Health and Human Services, Carson City, Nevada.

Funding

This research was supported by the Zhejiang Provincial Natural Science Foundation of China (grant number LR18C200001), the National Natural Science Foundation of China (grant no. 81773419), and the Intramural Research Program of the NIH National Cancer Institute.

Author information

Authors and Affiliations

Authors

Contributions

YZ and JJ designed and supervised the conduct of the whole study and obtained funding. PZ, LM, and FW were responsible for the data analysis of the exposures. WH provided statistical assistance. YZ and PZ wrote the first draft of the paper. All authors participated in the study design, generation of hypotheses, interpretation of data, and critical review of the report. The author(s) read and approved the final manuscript.

Corresponding authors

Correspondence to Yu Zhang or Jingjing Jiao.

Ethics declarations

Ethics approval and consent to participate

Written informed consent was obtained from all men and women who participated in the study, which was approved by the Institutional Review Board of the National Cancer Institute (Proposal #201807-0012).

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Figure S1.

Flow of participants in current NIH-AARP prospective cohort. Table S1. Categories for causes of death. Table S2. Baseline characteristics of participants according to corn oil, canola oil and olive oil consumption. Table S3. Spearman correlations between individual cooking oils. Table S4. Multivariable-adjusted HRs (95% CIs) of heart disease and stroke mortality according to individual oil consumption. Table S5. Multivariable-adjusted HRs (95% CIs) of mortality from non-cardiometabolic causes according to individual oil consumption. Table S6. Multivariable-adjusted HRs (95% CIs) of all-cause and cause-specific mortality according to lard consumption. Table S7. Data source of Fig. 3. Multivariable-adjusted HRs (95% CIs) for substituting tablespoon/d canola oil, corn oil, or olive oil for equivalent amounts of butter and margarine. Table S8. Multivariable-adjusted HRs (95% CIs) of all-cause and cause-specific mortality associated with consumption of stick margarine and other margarine. Table S9. Multivariable-adjusted HRs (95% CIs) of all-cause mortality from subgroup analyses. Table S10. Multivariable-adjusted HRs (95% CIs) of all-cause and cause-specific mortality from the sensitivity analysis that excluding those with extreme BMIs. Table S11. Multivariable-adjusted HRs (95% CIs) of all-cause and cause-specific mortality from the sensitivity analysis that further adjusting for a propensity score. Table S12. Multivariable-adjusted HRs (95% CIs) of all-cause and cause-specific mortality from the sensitivity analysis that further adjusted for history of hypertension and hypercholesteremia. Table S13. Multivariable-adjusted HRs (95% CIs) of all-cause and cause-specific mortality from the sensitivity analysis that further adjusted for aspirin and multivitamins use. Table S14. Multivariable-adjusted HRs (95% CIs) of all-cause and cause-specific mortality from the sensitivity analysis that further adjusted for the use of cholesterol-lowering medications (n = 293,918). Table S15. Multivariable-adjusted HRs (95% CIs) of all-cause and cause-specific mortality from the sensitivity analysis that excluding those with cardiovascular disease, cancer, or diabetes at baseline. Table S16. Multivariable-adjusted HRs (95% CIs) of all-cause and cause-specific mortality from the sensitivity analysis that excluding the first 4 years of follow-up. Table S17. Multivariable-adjusted HRs (95% CIs) of all-cause and cause-specific mortality from the sensitivity analysis that followed up for 8 years.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Zhang, Y., Zhuang, P., Wu, F. et al. Cooking oil/fat consumption and deaths from cardiometabolic diseases and other causes: prospective analysis of 521,120 individuals. BMC Med 19, 92 (2021). https://doi.org/10.1186/s12916-021-01961-2

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12916-021-01961-2

Keywords