Skip to main content

Association of milk consumption with all-cause mortality and cardiovascular outcomes: a UK Biobank based large population cohort study

Abstract

Background

The association of milk consumption with mortality and cardiovascular disease (CVD) outcomes was unclear.

Objective

The present study was performed to reveal the association of full cream, semi-skimmed, skimmed, soy, and other milk with all-cause mortality and CVD outcomes.

Methods

A prospective cohort study was performed using data from the UK Biobank. This study recruited 450,507 participants without CVD at baseline between 2006 and 2010 from UK Biobank and followed them up through 2021. Cox proportional hazard models were adopted to calculate the hazard ratios (HRs) and 95% confidence interval (CI) to understand the correlation between milk consumption and clinical outcomes. Subgroup and sensitivity analyses were further conducted.

Results

Among the participants, 435,486 (96.7%) were milk consumers. Multivariable model indicated that the adjusted HR of association between milk consumption and all-cause mortality was 0.84 (95% CI 0.79 to 0.91; P = 0.000) for semi-skimmed milk; 0.82 (0.76 to 0.88; P = 0.000) for skimmed milk and 0.83 (0.75 to 0.93; P = 0.001) for soy milk. Semi-skimmed, skimmed, and soy milk use were significantly related to lower risks of CVD mortality, CVD event, and stroke.

Conclusion

Compared with non-milk users, semi-skimmed milk, skimmed milk, and soy milk consumption were related to a lower risk of all-cause mortality and CVD outcomes. Among them, skim milk consumption was more beneficial for all-cause mortality, while soy milk consumption was more beneficial for CVD outcomes.

Introduction

Milk contains many essential nutrients, including calcium, phosphorus, and vitamin D. Many dietary guidelines commonly recommend milk in the diet [1,2,3]. Previous studies showed that the association of dairy milk consumption with mortality and cardiovascular disease (CVD) outcomes was highly controversial [1, 4,5,6,7,8]. Several reports showed a neutral association of dairy use with all-cause mortality and worse CVD outcomes [1, 4, 5, 9]. Meanwhile, a population-based cohort study demonstrated that dairy-related supplementation lowers the risk of CVD mortality [7]. Three prospective cohorts of U.S. healthcare professionals indicated that increased full cream milk use was associated with an increased risk of all-cause mortality [6]. However, several studies indicated that milk consumption was related to decreased all-cause mortality and CVD risk [8, 10]. Thus, the relationship of milk consumption with mortality and CVD outcomes is still unclear, and the different milk types' relationship with mortality is also less clear. A previous study reported that compared with those consuming whole milk, participants with reduced fat indicated a lower risk of mortality [11]. However, a cohort study showed that high milk consumption was connected to higher mortality [12]. Another study showed that intake of whole milk was related to a lower risk of CVD [13]. There is no direct evidence on whether semi-skimmed or skimmed milk was more beneficial to health than full-cream milk [14].

Large population-based cohort studies can provide additional information on the correlation of milk consumption with clinical outcomes, given the uncertainty. This prospective cohort study based on UK biobank cohort was conducted to explore the connection between intake of different types of milk, which including full cream, semi-skimmed, skimmed, soy and other milk type, and all-cause mortality and CVD outcomes and investigate the correction factors affecting these associations.

Methods

Study population

The UK Biobank included 502,409 participants from England, Scotland and Wales between March 2006 and July 2010. Inclusion criteria included age 40–69 years and living within a reasonable travel distance (10 miles) from one of the 22 assessment centers. Participants provided biological samples, physical measurements and baseline information in the assessment centers. After 6 participants dropped out of the UK Biobank project, 502,409 were obtained. Furthermore, 1526 participants without information for milk type used, 1283 participants lost to follow-up and 49,093 participants diagnosed with CVD at baseline were excluded. Then, 450,507 participants without CVD at baseline were enrolled in this study (Fig. 1).

Fig. 1
figure 1

Flowchart of participant selection

Exposure assessment

All the enrolled participants attended one of 22 assessment centers across the UK and finished a questionnaire using a touch screen device. One of the questions asked, “What type of milk do you mainly use?” The participants can choose their answer from a list of milk, including skimmed, semi-skimmed, full cream, soy and other milk types.

Assessment of outcomes

The study’s main outcomes were all-cause mortality, CVD mortality and CVD events. The secondary outcomes were stroke and myocardial infarction incidence and mortality. The date and cause of death were determined by linking to the Death Registry of the National Health Service (NHS) Information Centre for the UK and Welsh participants and the Death Registry of the Scottish NHS Central Registry for Scottish participants. Additionally, dates and reasons for hospitalization were identified by linking to Scottish morbidity records for Scottish participants and health event statistics for England and Wales participants. Detailed information can be found at http://content.digital.nhs.uk/services. This date was used as the end date for follow-up unless death or hospital admission occurred first for CVD outcomes. The CVD events were considered as hospitalization or death due to the following ICD-10 codes according to the hospital or death record: CVD death codes: I00–I99, CVD codes: I20–I25 and I60–I64, myocardial infarction codes: I21, I22, I23, I24.1 or I25.2 and stroke codes: I60–I64 [15,16,17,18].

Other variables

The baseline questionnaires were used to evaluate some potential confounding variables: sociodemographic factors (age, sex, ethnic background and household income), socioeconomic status (Townsend Deprivation Index), lifestyle habits (smoking status, alcohol consumption, tea intake, processed meat intake obesity, dietary intake (fruits and vegetables). Furthermore, medication use (blood pressure drugs and cholesterol-lowering drugs use), vitamin supplements (vitamin A, vitamin B, vitamin C, vitamin D, vitamin E, multivitamins, or folic acid), minerals supplements (calcium, iron, zinc or selenium) and comorbidities (hypertension, diabetes, high cholesterol and long-term illness) were evaluated. The Townsend Deprivation Index was provided in the UK Biobank. The information for medical history (hypertension, diabetes, high cholesterol and long-term illness) was obtained through self-report at baseline. The body mass index (BMI) was calculated as weight in kilograms divided by height in meters squared. In this study, obesity was considered as BMI ≥ 30, which consists with the European guidelines for obesity management in adults [19]. Details of these evaluations could be obtained from UK Biobank (www.ukbiobank.ac.uk).

Statistical analysis

Lilliefors test was adopted to detect if the data have a normal distribution [20, 21]. Continuous variables were indicated as mean ± standard deviation for normally distributed data; otherwise, they were shown as median and interquartile range. All categorical variables were presented as counts with percentages. Cox regression models were adopted to analyze the connection between milk consumption and all-cause mortality and CVD outcomes (CVD mortality, myocardial infarction mortality, stroke mortality, CVD events, myocardial infarction and stroke). Two models were constructed, and the basic model was adjusted for baseline age (years) and sex (male or female). The multivariable model was further adjusted for Townsend Deprivation Index, ethnic background (white or others), household income (< £18,000, £18,000–£30,999, £31,000–£51,999, £52,000–£100,000, or > £100,000), obesity (yes or no), fruit consumption, vegetable consumption, smoking status (current, never or previous). In addition, alcohol intake (< 1, 1–2, 3–4, or > 4 times/week), tea intake (< 2.0, 2.0–3.9, or ≥ 4.0 servings/day), processed meat intake (< 2, 2–4, 5–6, or > 6 times/week), vitamin use (yes or no), minerals use (yes or no), blood pressure drugs use (yes or no), cholesterol-lowering drugs use (yes or no), hypertension (yes or no), diabetes (yes or no), high cholesterol (yes or no) and longstanding illness (yes or no).

Additionally, several subgroup analyses were conducted by age (≥ 60 vs < 60 years), ethnic background (white vs others), sex (male vs female), current smoking status (yes vs no), diabetes (yes vs no), high cholesterol (yes vs no), hypertension (yes vs no), obesity (yes vs no), longstanding illness (yes vs no), blood pressure drugs use (yes vs no), cholesterol-lowering drugs use (yes vs no) and minerals use (yes vs no).

Sensitivity analysis

Several sensitivity analyses were conducted to determine the stability of cox regression model results. First, the participants who had undergone an outcome event during the first 2 years of follow-up were removed. Second, the participants who use vitamin were excluded. Third, all missing covariate data were imputed using multiple imputations. Fourth, the participants diagnosed with cancer at baseline were removed. All results were indicated as HR and 95% CI. R version 4.1.2 (www.r-project.org) was adopted for the analysis, and a two-sided P value < 0.05 was set as statistically significant.

Results

Characteristics of participants at baseline

The median age of participants was 57 years (interquartile range 50.00–63.00); 250,478 (55.6%) were women, 423,954 (94.1%) were white, and 435,486 (96.7%) were milk consumers. Full cream, semi-skimmed, skimmed, soy and other milk were used by 31,418, 290,422, 89,940, 17,853, and 5,853 participants, respectively. Semi-skimmed milk (64.5%) was the most commonly consumed, followed by skimmed (20.0%), full cream (7.0%), soy (4.0%) and others (1.3%). Skimmed milk consumers were more likely to be women, white, obese, former smokers, blood pressure drug users, cholesterol-lowering drug users, and in poor health. The baseline characteristics of enrolled participants are indicated in Table 1.

Table 1 Baseline features of participants

Association between milk consumption and all-cause mortality

The median follow-up time was 12.1 years (IQR, 11.4 to 12.8 years). During the follow-up period, 21,119 deaths were observed, including 3,679 CVD, 810 myocardial infarction, and 718 stroke deaths. The basic model with adjusted age and sex found that semi-skimmed, skimmed, and soy milk consumption was connected to a lower risk of all-cause mortality, CVD mortality, myocardial infarction mortality and CVD events.

After adjusting for age, sex, and Townsend Deprivation Index, etc., the adjusted hazard ratios of association between milk consumption and all-cause mortality were 0.84 (95% CI 0.79 to 0.91; P = 0.000) for semi-skimmed milk; 0.82 (0.76 to 0.88; P = 0.000) for skimmed milk and 0.83 (0.75 to 0.93; P = 0.001) for soy milk. However, no significant correlation was observed between milk consumption and death for full cream milk (HR 1.05; 95% CI 0.97 to 1.14; P = 0.209) and other milk (HR 1.09; 95% CI 0.96 to 1.24; P = 0.163).

Milk consumption and CVD outcomes

A total of 810 myocardial infarction deaths, 6443 myocardial infarction events, 718 stroke deaths and 6291 stroke events were recorded. As shown in Table 2, that semi-skimmed milk consumption was significantly associated with lower risks of CVD mortality (HR 0.75, 95% CI 0.63 to 0.88; P = 0.000), CVD events (HR 0.89, 95% CI 0.81 to 0.97; P = 0.007), stroke mortality (HR 0.70, 95% CI 0.49 to 1.00; P = 0.048) and stroke (HR 0.82, 95% CI 0.72 to 0.93; P = 0.002). Similarly, skimmed milk was found related to lower risks of CVD mortality (HR 0.74, 95% CI 0.62 to 0.88; P = 0.001), CVD events (HR 0.87, 95% CI 0.79 to 0.96; P = 0.003), and stroke (HR 0.80, 95% CI 0.70 to 0.92; P = 0.001). Soy milk was observed to be connected to lower risks of CVD mortality (HR 0.67, 95% CI 0.51 to 0.86; P = 0.002) CVD events (HR 0.82, 95% CI 0.72 to 0.93; P = 0.002), and stroke (HR 0.76, 95% CI 0.63 to 0.92; P = 0.004). However, no significant associations were found between full cream milk or other milk consumption and CVD outcomes (P > 0.05).

Table 2 Associations of milk type used with risk of all cause mortality and cardiovascular outcomes. Values were numbers (percentages) unless stated otherwise

Subgroup analysis

Stratified analysis was performed to observe the correlation between milk consumption and risk of all-cause mortality and CVD outcomes via multivariable model. Full cream milk consumption was connected to an increased risk of all-cause mortality (HR 1.20, 95% CI 1.05 to 1.37; P = 0.007) in participants who were women or with no hypertension (HR 1.13, 95% CI 1.02 to 1.26; P = 0.019). Moreover, other milk consumption was related to a higher risk of all-cause mortality in participants with no hypertension (HR 1.24, 95% CI 1.06 to 1.45; P = 0.008), not obese (HR 1.22, 95% CI 1.05 to 1.41; P = 0.008), not using blood pressure drugs (HR 1.17, 95% CI 1.01 to 1.36; P = 0.043) (Fig. 2).

Fig. 2
figure 2

Association of milk consumption and the risk of all cause mortality stratified by potential risk factors. The results were adjusted for baseline age (years) and sex (male or female), Townsend Deprivation Index, ethnic background (white or others), household income (< £18,000, £18,000–£30,999, £31,000–£51,999, £52,000–£100,000, or > £100,000), obesity (yes or no), fruit consumption, vegetable consumption, smoking status (current, never or previous), alcohol intake (< 1, 1–2, 3–4, or > 4 times/week), tea intake (< 2.0, 2.0–3.9, or ≥ 4.0 servings/day), processed meat intake (< 2, 2–4, 5–6, or > 6 times/week), vitamin use (yes or no), minerals use (yes or no), blood pressure drugs use (yes or no), cholesterol lowering drugs use (yes or no), hypertension (yes or no), diabetes (yes or no), high cholesterol (yes or no) and longstanding illness (yes or no)

The association between semi-skimmed milk intake and all cause mortality was stronger among participants with hypertension (P for interaction = 0.021), and connection of skimmed milk consumption with all-cause mortality was stronger among smokers (P for interaction = 0.010), hypertensive participants (P for interaction = 0.004) and blood pressure drug users (P for interaction = 0.030) for all-cause mortality. Additionally, the relationship between soy milk consumption and all-cause mortality was stronger among smokers (P for interaction = 0.002).

The association between semi-skimmed milk consumption and CVD outcomes was stronger among smokers (P for interaction = 0.047) and non-minerals users (P for interaction = 0.032) for CVD mortality and hypertensive participants (P for interaction = 0.034) for CVD events. Skimmed milk consumption and CVD outcomes were stronger among participants not on mineral supplements (P for interaction = 0.003) for CVD mortality and hypertension (P for interaction = 0.043) for CVD events. The connection between soy milk consumption and CVD outcomes was stronger among participants not using mineral supplements (P for interaction = 0.020) for CVD mortality and younger than 60 years old (P for interaction = 0.040 and not using mineral supplements (P for interaction = 0.027) for CVD events. We did not observe other significant interactions (all P for interaction ≥ 0.05). The associations of milk consumption with all-cause mortality and CVD outcomes were not modified by any other risk factors (Figs. 2 and 3).

Fig. 3
figure 3

Associations of milk consumption and the risk of (A) cardiovascular mortality and (B) cardiovascular events stratified by potential risk factors. The results were adjusted for baseline age (years) and sex (male or female), Townsend Deprivation Index, ethnic background (white or others), household income (< £18,000, £18,000–£30,999, £31,000–£51,999, £52,000–£100,000, or > £100,000), obesity (yes or no), fruit consumption, vegetable consumption, smoking status (current, never or previous), alcohol intake (< 1, 1–2, 3–4, or > 4 times/week), tea intake (< 2.0, 2.0–3.9, or ≥ 4.0 servings/ day), processed meat intake (< 2, 2–4, 5–6, or > 6 times/ week), vitamin use (yes or no), minerals use (yes or no), blood pressure drugs use (yes or no), cholesterol lowering drugs use (yes or no), hypertension (yes or no), diabetes (yes or no), high cholesterol (yes or no) and longstanding illness (yes or no)

Sensitivity analysis

The relationship of milk consumption with decreased risk of all-cause mortality did not change significantly after removing participants who developed CVD events or died during the first two years of follow-up (Additional file 1: Table S1). When the participants who took vitamin were excluded, the results did not change appreciably, but full cream, semi-skimmed, skimmed, and soy milk were related to lower risks of stroke mortality (Additional file 1: Table S2). Then, the result remained stable when all the missing covariate data were imputed using multiple imputations (Additional file 1: Table S3). Finally, when the participants diagnosed with cancer at baseline were removed, the result remained stable (Additional file 1: Table S4).

Discussion

Milk is an important food source that positively affects nutrient and energy intake, but heavy metals can counteract these benefits and affect human health. For example, the toxic trace elements lead (Pb) and cadmium (Cd) [22] can disturb the trace mineral composition of milk and negatively affect its nutritional quality [23]. Previous reports showed the presence of heavy metals in milk collected from different regions of the world [24, 25]. According to previous studies, the concentrations of Pb and Cd in milk were generally higher in developing countries but lower in developed countries, heavy metals such as Pb and Cd. Both of Pb and Cd can cross the placental barrier, which may affect the fetal brain differentiation, leading to neurotoxic effects, including memory loss and language impairment [26, 27].

Milk is rich in protein, amino acids, minerals and trace elements. Milk is a good nutrient for human health [11]. In this prospective cohort study with 450,507 participants recruited from UK Biobank between 2006 and 2010, it was observed that compared with non-milk users, semi-skimmed, skimmed and soya milk consumption was related to a 16%, 18% and 17% lower risk of all-cause mortality, respectively. Additionally, semi-skimmed skimmed and soya milk consumption was significantly correlated with decreased CVD mortality, CVD events and stroke risks. The protective association of skimmed milk consumption against all-cause mortality was stronger in smokers and blood pressure medication users. Association of semi-skimmed, skimmed and soya milk consumption with CVD mortality was stronger among participants not on mineral supplements. Correlation of semi-skimmed and skimmed milk with CVD events was stronger among participants with hypertension.

The association of milk consumption with all-cause mortality and CVD outcomes has been controversial [7, 13, 28]. Xu et al. reported that participants with low fat milk intake indicated a decreased risk of mortality compared to participants consuming whole milk [11]. Similarly, Wang et al. found that skimmed milk was related to lower risk of mortality compared to whole milk [8], while another study using two large Swedish cohorts showed that high milk consumption was connected to higher mortality [12]. In the present study, we observed that full cream milk was not related to all-cause mortality and CVD outcomes, while semi-skimmed and skimmed milk were connected to lower all-cause mortality risk, CVD mortality and incidence of CVD. Skimmed milk showed better protection than semi-skimmed milk. Studies have shown that increased intake of low-fat dairy was associated with lower CVD mortality, but not high-fat dairy intake [29]. Consumption of low-fat dairy products was related to a reduction in hypertension risk and serum LDL-C [29], which may account for the potential CVD benefits of low-fat milk. In this study, low-fat milk consumers were more likely to be obese (41.5%), possibly because obese people tended to choose skim milk. A previous randomized controlled trial showed no significant effect of low-fat milk or whole cheese consumption on ambulatory blood pressure compared with a dairy-free diet [30]. However, compared with a low-fat milk diet, consumption of regular-fat cheese increased LDL particle size, which may be associated with CVD events [31].

The connection between skimmed milk consumption and CVD outcomes in this study was stronger among participants with hypertension for CVD events. A previous study indicated that hypertension is a risk factor for CVD [32]. It was also related to hyperlipidemia, with many common risk factors [33]. Skimmed milk may provide more benefits for hypertensive patients.

Previous studies have shown that dietary habits were one of the major determinants of CVD, responsible for approximately half of deaths from cardiometabolic disease [34]. Skimmed milk may show better results against CVD outcomes compared to full cream milk. At the same time, milk with different fat levels exhibits different acidity and stability, casein content, total solids content and complex interactions [8]; thus, skimmed milk benefits may not be solely attributable to lower milk fat intake [35].

A previous study indicated that skimmed milk presented higher inhibitory activities for oxidative damage of deoxyribose [36]. Additionally, another report showed that soya milk consumption could lead to decreased oxidative damage for DNA bases [37]. Moreover, skimmed milk proteins exhibited ability of good free-radical quenching [38]. Antioxidant effect of milk may be protective against cardiovascular disease, while the protective association of semi-skimmed, skimmed, and soy milk with lower risks of all-cause mortality and CVD outcomes needs to be further explored.

In this study, soy milk consumption was related to 17%, 33%, 18% and 24% lower risk of all-cause mortality, CVD mortality, CVD events and stroke, respectively. Soy milk showed a better reduction in CVD outcomes than semi-skimmed and skimmed milk. Previous clinical trials indicated that soy milk consumption positively affected several cardiovascular events [39, 40]. Positive health effects of whole soy products, such as soy milk, have also been documented [41, 42]; a previous study indicated that soy milk consumption could decrease triglycerides and systolic blood pressure levels. Furthermore, soy peptides could act as inhibitory peptides for angiotensin-converting enzymes, leading to a hypotensive effect [43]. The renin-angiotensin system plays a crucial role in the regulation of renal, cardiac and vascular physiology. Activation of the renin-angiotensin system is a major cause of many common pathological conditions, including hypertension, heart failure, and renal disease. Metabolites of angiotensin I and II were initially considered biologically inactive, but there is now increasing evidence that they play key roles in cardiovascular physiology and pathophysiology [44].

Strengths and limitations

This study had several strengths, including the large sample size, minimal loss to follow-up and population-based cohort study design. However, there are several limitations. First, we adjusted the potential confounders; however, residual confounding can not be completely ruled out. Second, this study is observational, and no causality was obtained. Third, this study did not record detailed information on the dose and duration of milk use, so dose–response associations between milk consumption cannot be assessed. Fourth, due to the low response rate, potential selection biases may limit the generalizability of the conclusions in the wider UK population.

Conclusions

A total of 450,507 participants without CVD at baseline were enrolled in this study, and 435,486 (96.7%) were milk consumers. Compared with non-milk users, semi-skimmed, skimmed, and soy milk consumption was related to a lower risk of all-cause mortality and CVD outcomes. Skim milk consumption was more beneficial for all-cause mortality, while soy milk consumption was more beneficial for CVD events. This study recommends skimmed or soy milk over full cream milk to improve clinical outcomes. Additionally, the difference between skimmed and full cream milk should be discovered for the actual reason for this impact. For this purpose, ingredient analyses such as mass spectrum need to be conducted, and well-designed biological experiments should be performed based on the results.

Availability of data and materials

The datasets used and/or analyzed during the current study are available from the corresponding author upon reasonable request.

References

  1. Pala V, Sieri S, Chiodini P, et al. Associations of dairy product consumption with mortality in the European Prospective Investigation into Cancer and Nutrition (EPIC)-Italy cohort. Am J Clin Nutr. 2019;110:1220–30.

    Article  PubMed  Google Scholar 

  2. Drouin-Chartier JP, Cote JA, Labonte ME, et al. Comprehensive review of the impact of dairy foods and dairy fat on cardiometabolic risk. Adv Nutr. 2016;7:1041–51.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Wang SS, Lay S, Yu HN, Shen SR. Dietary guidelines for chinese residents (2016): comments and comparisons. J Zhejiang Univ Sci B. 2016;17:649–56.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Fontecha J, Calvo MV, Juarez M, Gil A, Martinez-Vizcaino V. Milk and dairy product consumption and cardiovascular diseases: an overview of systematic reviews and meta-analyses. Adv Nutr. 2019;10:S164–89.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Tognon G, Nilsson LM, Shungin D, et al. Nonfermented milk and other dairy products: associations with all-cause mortality. Am J Clin Nutr. 2017;105:1502–11.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Ding M, Li J, Qi L, et al. Associations of dairy intake with risk of mortality in women and men: three prospective cohort studies. BMJ. 2019;367: l6204.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Mazidi M, Mikhailidis DP, Sattar N, Howard G, Graham I, Banach M. Consumption of dairy product and its association with total and cause specific mortality—a population-based cohort study and meta-analysis. Clin Nutr. 2019;38:2833–45.

    Article  PubMed  Google Scholar 

  8. Wang S, Liu Y, Cai H, et al. Decreased risk of all-cause and heart-specific mortality is associated with low-fat or skimmed milk consumption compared with whole milk intake: a cohort study. Clin Nutr. 2021;40:5568–75.

    Article  CAS  PubMed  Google Scholar 

  9. Guo J, Astrup A, Lovegrove JA, Gijsbers L, Givens DI, Soedamah-Muthu SS. Milk and dairy consumption and risk of cardiovascular diseases and all-cause mortality: dose-response meta-analysis of prospective cohort studies. Eur J Epidemiol. 2017;32:269–87.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Wang XJ, Jiang CQ, Zhang WS, et al. Milk consumption and risk of mortality from all-cause, cardiovascular disease and cancer in older people. Clin Nutr. 2020;39:3442–51.

    Article  CAS  PubMed  Google Scholar 

  11. Xu X, Kabir A, Barr ML, Schutte AE. Different types of long-term milk consumption and mortality in adults with cardiovascular disease: a population-based study in 7236 Australian adults over 84 years. Nutrients. 2022;14:704.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Michaelsson K, Wolk A, Langenskiold S, et al. Milk intake and risk of mortality and fractures in women and men: cohort studies. BMJ. 2014;349: g6015.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Talaei M, Hosseini N, van Dam RM, et al. Whole milk consumption and risk of cardiovascular disease and mortality: Isfahan cohort study. Eur J Nutr. 2019;58:163–71.

    Article  PubMed  Google Scholar 

  14. Louie J, Rangan AM. No need to change dairy food dietary guidelines yet. Lancet. 2018;392:2242–4.

    Article  PubMed  Google Scholar 

  15. Li ZH, Gao X, Chung VC, et al. Associations of regular glucosamine use with all-cause and cause-specific mortality: a large prospective cohort study. Ann Rheum Dis. 2020;79:829–36.

    Article  CAS  PubMed  Google Scholar 

  16. Ma H, Li X, Sun D, et al. Association of habitual glucosamine use with risk of cardiovascular disease: prospective study in UK Biobank. BMJ. 2019;365: l1628.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Liu D, Li ZH, Shen D, et al. Association of sugar-sweetened, artificially sweetened, and unsweetened coffee consumption with all-cause and cause-specific mortality : a large prospective cohort study. Ann Intern Med. 2022. https://doi.org/10.7326/L22-0369.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Li ZH, Zhong WF, Liu S, et al. Associations of habitual fish oil supplementation with cardiovascular outcomes and all cause mortality: evidence from a large population based cohort study. BMJ. 2020;368: m456.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Yumuk V, Tsigos C, Fried M, et al. European guidelines for obesity management in adults. Obes Facts. 2015;8:402–24.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Abdi H, Molin P. Lilliefors/Van Soest’s test of normality. Thousand Oaks: Sage; 2007.

    Google Scholar 

  21. Dallal GE, Wilkinson L. An analytic approximation to the distribution of Lilliefors’s test statistic for normality. Am Stat. 1986;40(4):294–6.

    Google Scholar 

  22. Varol M, Sunbul MR. Macroelements and toxic trace elements in muscle and liver of fish species from the largest three reservoirs in Turkey and human risk assessment based on the worst-case scenarios. Environ Res. 2020;184: 109298.

    Article  CAS  PubMed  Google Scholar 

  23. Patra RC, Swarup D, Kumar P, Nandi D, Naresh R, Ali SL. Milk trace elements in lactating cows environmentally exposed to higher level of lead and cadmium around different industrial units. Sci Total Environ. 2008;404:36–43.

    Article  CAS  PubMed  Google Scholar 

  24. Licata P, Trombetta D, Cristani M, et al. Levels of “toxic” and “essential” metals in samples of bovine milk from various dairy farms in Calabria, Italy. Environ Int. 2004;30:1–6.

    Article  CAS  PubMed  Google Scholar 

  25. Boudebbouz A, Boudalia S, Bousbia A, Habila S, Boussadia MI, Gueroui Y. Heavy metals levels in raw cow milk and health risk assessment across the globe: a systematic review. Sci Total Environ. 2021;751: 141830.

    Article  CAS  PubMed  Google Scholar 

  26. Khalil N, Morrow LA, Needleman H, Talbott EO, Wilson JW, Cauley JA. Association of cumulative lead and neurocognitive function in an occupational cohort. Neuropsychology. 2009;23:10–9.

    Article  PubMed  Google Scholar 

  27. Rehman K, Fatima F, Waheed I, Akash M. Prevalence of exposure of heavy metals and their impact on health consequences. J Cell Biochem. 2018;119:157–84.

    Article  CAS  PubMed  Google Scholar 

  28. Stasinopoulos LC, Zhou A, Hypponen E. Association of supplemental calcium and dairy milk intake with all-cause and cause-specific mortality in the UK Biobank: a prospective cohort study. Br J Nutr. 2020;123:574–82.

    Article  CAS  PubMed  Google Scholar 

  29. Huo YKS, Bongard V, Simon C, et al. Low-fat and high-fat dairy products are differently related to blood lipids and cardiovascular risk score. Eur J Prev Cardiol. 2014;21:1557–67.

    Article  Google Scholar 

  30. Rancourt-Bouchard M, Gigleux I, Guay V, et al. Effects of regular-fat and low-fat dairy consumption on daytime ambulatory blood pressure and other cardiometabolic risk factors: a randomized controlled feeding trial. Am J Clin Nutr. 2020;111:42–51.

    Article  PubMed  Google Scholar 

  31. Grammer TB, Kleber ME, Marz W, et al. Low-density lipoprotein particle diameter and mortality: the Ludwigshafen risk and cardiovascular health study. Eur Heart J. 2015;36:31–8.

    Article  CAS  PubMed  Google Scholar 

  32. Kokubo Y, Matsumoto C. Hypertension is a risk factor for several types of heart disease: review of prospective studies. Adv Exp Med Biol. 2017;956:419–26.

    Article  PubMed  Google Scholar 

  33. Ruixing Y, Jinzhen W, Weixiong L, Yuming C, Dezhai Y, Shangling P. The environmental and genetic evidence for the association of hyperlipidemia and hypertension. J Hypertens. 2009;27:251–8.

    Article  PubMed  Google Scholar 

  34. Micha R, Penalvo JL, Cudhea F, Imamura F, Rehm CD, Mozaffarian D. Association between dietary factors and mortality from heart disease, stroke, and type 2 diabetes in the United States. JAMA. 2017;317:912–24.

    Article  PubMed  PubMed Central  Google Scholar 

  35. O’Sullivan TA, Schmidt KA, Kratz M. Whole-fat or reduced-fat dairy product intake, adiposity, and cardiometabolic health in children: a systematic review. Adv Nutr. 2020;11:928–50.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Chiang SH, Wang SY, Chang CY, Chen CW. Bovine colostrum whey protein hydrolysate inhibits cell DNA damage and LDL oxidation in vitro. Molecules. 2017;22:456.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Mitchell JH, Collins AR. Effects of a soy milk supplement on plasma cholesterol levels and oxidative DNA damage in men–a pilot study. Eur J Nutr. 1999;38:143–8.

    Article  CAS  PubMed  Google Scholar 

  38. Prakash PK, Eligar SM, Prakruthi M, Jyothi LA. Comparative assessment of antioxidant and immunomodulatory properties of skimmed milk protein hydrolysates and their incorporation in beverage mix. J Sci Food Agric. 2022;102:6414–22.

    Article  CAS  PubMed  Google Scholar 

  39. Miraghajani MS, Esmaillzadeh A, Najafabadi MM, Mirlohi M, Azadbakht L. Soy milk consumption, inflammation, coagulation, and oxidative stress among type 2 diabetic patients with nephropathy. Diabetes Care. 2012;35:1981–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Miraghajani MS, Najafabadi MM, Surkan PJ, Esmaillzadeh A, Mirlohi M, Azadbakht L. Soy milk consumption and blood pressure among type 2 diabetic patients with nephropathy. J Ren Nutr. 2013;23:277–82.

    Article  CAS  PubMed  Google Scholar 

  41. Reinwald S, Akabas SR, Weaver CM. Whole versus the piecemeal approach to evaluating soy. J Nutr. 2010;140:2335S-2343S.

    Article  CAS  PubMed  Google Scholar 

  42. Torres YTN, Palacios-Gonzalez B, Noriega-Lopez L. Tovar-Palacio AR [Glycemic, insulinemic index, glycemic load of soy beverage with low and high content of carbohydrates]. Rev Invest Clin. 2006;58:487–97.

    Google Scholar 

  43. De Leo F, Panarese S, Gallerani R, Ceci LR. Angiotensin converting enzyme (ACE) inhibitory peptides: production and implementation of functional food. Curr Pharm Des. 2009;15:3622–43.

    Article  PubMed  Google Scholar 

  44. Paz OM, Riquelme JA, Garcia L, et al. Counter-regulatory renin-angiotensin system in cardiovascular disease. Nat Rev Cardiol. 2020;17:116–29.

    Article  Google Scholar 

Download references

Acknowledgements

This research has been conducted using the UK Biobank Resource (https://www.ukbiobank.ac.uk) under Application Number 80610.

Funding

This work was supported by the National Natural Science Foundation of China (Grant nos. 82072441, 81871783, and 81672176) and Hunan Provincial Natural Science Outstanding Youth Fund (2022JJ10095). The study funders/sponsors had no role in the design and conduct of the study; collection, management, analysis, and interpretation of the data; preparation, review, or approval of the manuscript; and decision to submit the manuscript for publication.

Author information

Authors and Affiliations

Authors

Contributions

TL designed the study and performed the analysis. JZ and RW drafted the manuscript. ZW, ZL, WW and TL contributed to the revision of the manuscript. All authors have read and approved the final manuscript.

Corresponding authors

Correspondence to Rongjun Wan or Tang Liu.

Ethics declarations

Ethics approval and consent to participate

This research was conducted under the UK Biobank application 80610. The UK Biobank study was approved by the National Health and Social Care Information Management Board and the North West Multicentre Research Ethics Committee (11/NW/0382).

Consent for publication

All studies included in this study got informed consent from each study participant and each study was approved by the ethics committee or institutional review board.

Competing interests

The authors declare that they do not have any 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: Table S1.

Associations of milk consumption with the risk of all-cause mortality and CVD outcomes after excluding participants who experienced an outcome event during the first two years of follow-up (n = 448,657). Table S2. Associations of milk consumption with the risk of all-cause mortality and CVD outcomes after excluding participants who use vitamin (n = 307,594). Table S3. Associations of milk consumption with the risk of all-cause mortality and CVD outcomes with all missing covariate data imputed using multiple imputation (n = 450,507). Table S4. Associations of milk consumption with the risk of all-cause mortality and CVD outcomes after excluding participants who were diagnosed with cancer at baseline (n = 413,422).

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

Zhou, J., Wu, Z., Lin, Z. et al. Association of milk consumption with all-cause mortality and cardiovascular outcomes: a UK Biobank based large population cohort study. J Transl Med 21, 130 (2023). https://doi.org/10.1186/s12967-023-03980-4

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12967-023-03980-4

Keywords