Intended for healthcare professionals

CCBYNC Open access
Research

Associations of dairy intake with risk of mortality in women and men: three prospective cohort studies

BMJ 2019; 367 doi: https://doi.org/10.1136/bmj.l6204 (Published 27 November 2019) Cite this as: BMJ 2019;367:l6204
  1. Ming Ding, research associate1,
  2. Jun Li, research associate1 2,
  3. Lu Qi, professor1 3,
  4. Christina Ellervik, assistant professor4 5 6,
  5. Xuehong Zhang, assistant professor7,
  6. JoAnn E Manson, professor2 6 7,
  7. Meir Stampfer, professor2 7,
  8. Jorge E Chavarro, associate professor1 2 7,
  9. Kathryn M Rexrode, associate professor8,
  10. Peter Kraft, professor2,
  11. Daniel Chasman, associate professor6 9 10,
  12. Walter C Willett, professor1 2 7,
  13. Frank B Hu, professor1 2 7
  1. 1Department of Nutrition, Harvard School of Public Health, Boston, MA 02115, USA
  2. 2Department of Epidemiology, Harvard School of Public Health, Boston, MA, USA
  3. 3Department of Epidemiology, School of Public Health and Tropical Medicine, Tulane University, New Orleans, LA, USA
  4. 4Department of Clinical Medicine, Faculty of Health and Medical Sciences, University of Copenhagen, Copenhagen, Denmark
  5. 5Department of Laboratory Medicine, Boston Children’s Hospital, Boston, MA, USA
  6. 6Division of Preventive Medicine, Brigham and Women’s Hospital and Harvard Medical School, Boston, MA, USA
  7. 7Channing Division of Network Medicine, Brigham and Women’s Hospital and Harvard Medical School, Boston, MA, USA
  8. 8Division of Women’s Health, Brigham and Women’s Hospital and Harvard Medical School, Boston, MA, USA
  9. 9Division of Genetics, Brigham and Women's Hospital and Harvard Medical School, Boston, MA, USA
  10. 10Broad Institute of MIT and Harvard, Cambridge, MA, USA
  1. Correspondence to: M Ding mid829{at}mail.harvard.edu
  • Accepted 17 October 2019

Abstract

Objective To examine the association of consumption of dairy foods with risk of total and cause specific mortality in women and men.

Design Three prospective cohort studies with repeated measures of diet and lifestyle factors.

Setting Nurses’ Health Study, Nurses’ Health Study II, and the Health Professionals Follow-up Study, in the United States.

Participants 168 153 women and 49 602 men without cardiovascular disease or cancer at baseline.

Main outcome measure Death confirmed by state vital records, the national death index, or reported by families and the postal system. During up to 32 years of follow-up, 51 438 deaths were documented, including 12 143 cardiovascular deaths and 15 120 cancer deaths. Multivariable analysis further adjusted for family history of cardiovascular disease and cancer, physical activity, overall dietary pattern (alternate healthy eating index 2010), total energy intake, smoking status, alcohol consumption, menopausal status (women only), and postmenopausal hormone use (women only).

Results Compared to the lowest category of total dairy consumption (average 0.8 servings/day), the multivariate pooled hazard ratio for total mortality was 0.98 (95% confidence interval 0.96 to 1.01) for the second category of dairy consumption (average 1.5 servings/day), 1.00 (0.97 to 1.03) for the third (average 2.0 servings/day), 1.02 (0.99 to 1.05) for the fourth (average 2.8 servings/day), and 1.07 (1.04 to 1.10) for highest category (average 4.2 servings/day; P for trend <0.001). For the highest compared to the lowest category of total dairy consumption, the hazard ratio was 1.02 (0.95 to 1.08) for cardiovascular mortality and 1.05 (0.99 to 1.11) for cancer mortality. For subtypes of dairy products, whole milk intake was significantly associated with higher risks of total mortality (hazard ratio per 0.5 additional serving/day 1.11, 1.09 to 1.14), cardiovascular mortality (1.09, 1.03 to 1.15), and cancer mortality (1.11, 1.06 to 1.17). In food substitution analyses, consumption of nuts, legumes, or whole grains instead of dairy foods was associated with a lower mortality, whereas consumption of red and processed meat instead of dairy foods was associated with higher mortality.

Conclusion These data from large cohorts do not support an inverse association between high amount of total dairy consumption and risk of mortality. The health effects of dairy could depend on the comparison foods used to replace dairy. Slightly higher cancer mortality was non-significantly associated with dairy consumption, but warrants further investigation.

Introduction

Dairy products are widely consumed worldwide. They are important sources of protein, vitamin D, and calcium, but several constituents of dairy foods such as saturated fat and cholesterol might have adverse effects on health. Previous observational studies have examined associations of dairy foods with multiple health outcomes. In general, the associations between total dairy and risk of hypertension,12 type 2 diabetes,3 cardiovascular disease,4 colorectal cancer,5 and breast cancer6 were null or moderately inverse. However, higher intakes of dairy products have also been associated with greater risk of prostate cancer7 and ovarian cancer.89

The results have also been inconsistent in prospective observational studies of dairy intake and mortality: positive associations have been seen in some studies,101112 while inverse or null associations have been seen in others.1314151617 In two meta-analyses summarizing these findings, overall dairy intake was not significantly associated with risk of mortality but the heterogeneity between studies was high.1819 Reasons for the heterogeneity across studies included limited sample size in some studies, different types of dairy products, different overall dietary patterns across populations, and inconsistent control for confounding variables.

Therefore, we examined the association of total dairy consumption and dairy subtypes with total and cause specific mortality in three large cohorts, including the Nurses’ Health Study (NHS), Nurses’ Health Study II (NHSII), and Health Professional Follow-up Study (HPFS). The total sample size of our study was 217 755, and we documented 51 438 deaths during 29-32 years of follow-up.

Methods

Study population

The NHS began in 1976, when 121 700 female registered nurses aged 30-55 residing in 11 US states were recruited to complete a baseline questionnaire about their lifestyle and medical history. The NHSII was established in 1989 and consisted of 116 671 younger female registered nurses, aged 25-42 at baseline. The nurses responded to a baseline questionnaire similar to the NHS. The HPFS was initiated in 1986, and was composed of 51 529 male dentists, pharmacists, veterinarians, optometrists, osteopathic physicians, and podiatrists, aged 40-75 at baseline. Participants returned a baseline questionnaire about detailed medical history, lifestyle, and usual diet. Detailed descriptions of the three cohorts are provided elsewhere.202122 In all the three cohorts, questionnaire data were collected at baseline and biennially thereafter, to update information on lifestyle factors and the occurrence of chronic diseases. These cohorts consist of participants with about 95% European ancestry. We excluded participants who reported cardiovascular disease or cancer at baseline for these analyses (1984 for the NHS, 1989 for the NHSII, and 1986 for the HPFS).23 We further excluded participants with extreme and implausible daily energy intakes (<500 or >3500 kcal for women; <800 or >4200 kcal for men; 1 kcal=4.18 kJ=0.00418 MJ). The study protocol was approved by the institutional review boards of the Brigham and Women’s Hospital and Harvard T H Chan School of Public Health, and those of participating registries as required.

Assessment of dairy consumption

In 1984, a 116 item food frequency questionnaire was administered to the NHS participants to obtain information on usual intake of food and beverages. Starting in 1986, an expanded 131 item food frequency questionnaire was administered every four years to update diet. Using a similar food frequency questionnaire, dietary data was collected every four years from the HPFS participants starting in 1986 and from the NHSII participants starting in 1989. In all food frequency questionnaires, participants were asked how often (in nine categories, ranging from “never or less than once per month” to “six or more times per day”), on average, they consumed a standard portion size of each food item during the previous year.

The questionnaire items for dairy products included “skim/low fat milk,” “whole milk,” “ice cream,” “sherbert,” “yogurt,” “cottage/ricotta cheese,” “cream cheese,” “other cheese,” and “cream.” The standard serving sizes were a 240 mL glass for skim/low fat milk and whole milk; one tablespoon (6 g) for cream; 1/2 cup (120 mL) for sherbet or frozen yogurt, ice cream, and cottage and ricotta cheese; and 30 mL for cream cheese and other cheese. Butter consumption was not included because it is nutritionally distinct from other dairy foods. Consumption of total dairy foods was calculated as the sum of intakes of all types of dairy products. We further divided the individual dairy foods into four groups: milk, cheese, yogurt, and other types of dairy foods including ice cream, cream, and sherbert. The validity and reproducibility of the food frequency questionnaire have been described in detail elsewhere.24252627 Briefly, correlation coefficients comparing intakes based on four dietary records, each record lasting one week, collected over a one year period with a single food frequency questionnaire ranged from 0.57 for hard cheese to 0.97 for yogurt.26

Assessment of covariates and effect modifiers

Updated information on age, weight, smoking status, physical activity, family history of cardiovascular disease and cancer, self reported diagnosis of diseases (including hypertension, hypercholesterolemia, cardiovascular disease, and cancer), and cardiovascular drug treatment use (including aspirin, diuretics, calcium blockers, β blockers, other antihypertensive drugs, statins, and other cholesterol lowering drugs) was collected in the biennial follow-up questionnaires in the NHS, NHSII, and HPFS. For female participants, we also ascertained data on menopausal status and postmenopausal hormone use. To measure the overall diet quality using food frequency questionnaire data, we calculated the alternate healthy eating index in the three cohorts.28

Assessment of deaths

Our primary endpoint was death from any cause ascertained by 31 December 2016. In the NHS, NHSII, and HPFS, we performed systematic searches of state vital records and the national death index. This search was supplemented by reports from family members and postal authorities. Using these methods, we were able to ascertain more than 98% of the deaths in each cohort.29 Physicians who were blinded to data on dairy consumption and other risk factors reviewed death certificates and medical records to classify the cause of death according to ICD-8 and ICD-9 (international classification of diseases, 8th and 9th revisions).

Statistical analysis

We followed the participants from the date of the return of the baseline questionnaire until their date of death or the end of follow-up (31 December 2016), whichever came first. We used Cox proportional hazards regression models to examine the associations between total dairy consumption and risk of mortality. The regression models included follow-up time as the time scale and were stratified by age in years. In the multivariable analysis, we further adjusted for family history of cardiovascular disease and cancer, physical activity, overall dietary pattern (alternate healthy eating index 2010), total energy intake, smoking status, alcohol consumption, menopausal status (women only), and postmenopausal hormone use (women only), all of which were updated from follow-up questionnaires. We additionally adjusted for baseline body mass index, hypertension, and hypercholesterolemia in both men and women. We used the cumulative average of dietary intakes from baseline to the censoring events to best represent long term diet and minimize within-person variation,23 and we stopped updating dietary information after self reported diagnosis of intermediate outcomes, including type 2 diabetes, cardiovascular disease, and cancer. We modelled dairy intakes as categorical variables, and tested for trend by modelling the median value within each category as a continuous variable.

All analyses were performed separately in each cohort, and then pooled to obtain the overall effect using a fixed effect model. We used random effects models if significant heterogeneity across cohorts was found. Statistical heterogeneity across studies was assessed by Cochrane Q test, with P<0.1 indicating significant heterogeneity between studies.30 We examined dose-response associations of dairy intake with risk of mortality by pooling individual level data in the three cohorts. We used restricted cubic splines with three knots to flexibly model the non-linear association. To test for a potential non-linear association between dairy consumption and risk of mortality, we used a likelihood ratio test to compare the model with cubic spline terms to the model with a linear term of dairy consumption, with P<0.05 denoting significant non-linearity. A significant P value for non-linearity indicated that the cubic spline model fitted the data better; therefore, we further tested for overall significance of curve using a likelihood ratio test by comparing the model with cubic spline terms to the model only with intercept. Otherwise, a non-significant P value for non-linearity indicated that the linear model fitted the data better.

Stratified analyses were conducted according to development of cardiovascular disease or cancer, cardiovascular drug treatment use, follow-up time (≤14 years or >14 years), calendar year (≤2000 or >2000), body mass index (≤25 or >25), age (≤75 or >75), smoking status (current smoker or never/past smoker), alternate healthy eating index (median score or less or more than median score), physical activity (median score or less or more than median score), total fat intake (median score or less or more than median score), calcium intake (median score or less or more than median score), and vitamin D intake (median score or less or more than median score). We treated the effect modifiers as time varying variables, and tested for potential effect modification by including an interaction term between the exposure and potential effect modifier in the Cox model and conducting a likelihood ratio test comparing the models with and without the interaction term.

We estimated the associations of substituting one serving per day of nuts and legumes, whole grain, poultry, or red and processed meat for total dairy by including both variables in continuous form in the same multivariate model. The difference in their regression coefficients were used to estimate the hazard ratios for the substitution and the covariance to estimate the 95% confidence intervals.31 All statistical tests were two sided at type I error rate of 0.05 and performed using SAS version 9.2 for UNIX (SAS Institute, Cary, NC). Because we tested multiple types of dairy foods and several death related outcomes, we conservatively corrected for multiple testing using Bonferroni correction and controlled for type I error rate at 0.05/60=0.001.

Patient and public involvement

No patients were involved in setting the research question or the outcome measures, nor were they involved in developing plans for design or implementation of the study. No patients were asked to advise on interpretation or writing up of results. There are no plans to disseminate the results of the research to study participants or the relevant patient community.

Results

Our original analysis included 217 755 participants, of whom 74 805 were from the NHS, 93 348 from the NHSII, and 49 602 from the HPFS; 168 153 women and 49 602 men without cardiovascular disease or cancer at baseline. Table 1 shows the baseline characteristics of the participants by frequency of total dairy consumption. Those who consumed higher amounts of dairy foods drank less alcohol, were more likely to be physically active, were less likely to be current smokers, and had a slightly lower alternate healthy eating index.

Table 1

Baseline characteristics of participants by fifths of total dairy consumption in the Nurses’ Health Study, Nurses’ Health Study II, and Health Professionals Follow-up Study. Data are mean values unless stated otherwise

View this table:

During 29-32 years of follow-up, 51 438 deaths were documented in the three cohorts, including 12 143 cardiovascular deaths and 15 120 cancer deaths. In an age adjusted model, dairy intake was positively associated with risk of total mortality in the NHSII and HPFS, but not NHS (table 2), which might be due to the different distribution of covariates according to dairy intake across the three cohorts as shown in table 1. However, after adjusting for covariates, dairy intake was positively associated with risk of mortality across the three cohorts (P for heterogeneity >0.1). Compared to the lowest category of total dairy consumption, the multivariable pooled hazard ratio for death was 1.07 (95% confidence interval 1.04 to 1.10) for the highest category of dairy consumption (P for linear trend <0.001; table 2). The findings remained significant after correcting for multiple testing. We further investigated the associations of dairy intake with risk of mortality using a dose-response analysis. We saw a non-linear relation between total dairy intake and risk of total mortality (P for non-linearity <0.001; P for overall significance of the curve <0.001; P for linear association <0.001) and risk of cardiovascular mortality (P=0.001; P=0.004; P=0.60; fig 1), such that moderate intake was associated with a slightly lower risk, whereas high intake was associated with an increased risk. Total dairy intake was associated with higher risk of cancer mortality in a linear fashion (P=0.10; P=0.04; P=0.048).

Table 2

Hazard ratios of total, cardiovascular, and cancer mortality according to fifths of total dairy intake in the Nurses’ Health Study, Nurses’ Health Study II, and Health Professionals Follow-up Study

View this table:
Fig 1
Fig 1

Dose-response associations of total dairy intake with risks of total mortality and mortality due to cardiovascular disease and cancer, based on pooled analyses of cohorts from the Nurses’ Health Study, Nurses’ Health Study II, and Health Professionals Follow-up Study. Multivariable adjusted model was further adjusted for family history of cancer (yes, no), family history of cardiovascular disease (yes, no), baseline disease status (hypertension, hypercholesterolemia), baseline body mass index (<20.9, 21.0-22.9, 23.0-24.9, 25.0-29.9, 30.0-34.9, ≥35.0), physical activity (<3.0, 3.0-8.9, 9.0-17.9, 18.0-26.9, ≥27.0 metabolic equivalent hours per week), overall dietary pattern (alternate healthy eating index score, in fifths), total energy intake (fifths), smoking status (Nurses’ Health Study and Health Professionals Follow-up Study: never, former (categorized into 1-4, 5-14, 15-24, 25-34, 35-44, ≥45 cigarettes/day, or unknown), current (categorized in to 1-4, 5-14, 15-24, 25-34, 35-44, ≥45 cigarettes/day, or unknown), alcohol consumption (0, ≥0-5, ≥5-10, ≥10-15, ≥15 g/day), postmenopausal status (yes, no; women only), and current postmenopausal hormone use (yes, no; women only)

We further examined associations of subtypes of dairy foods with risk of total mortality and mortality due to cardiovascular disease and cancer (table 3, table S1). We categorized individual dairy products into four groups, and used different cutoff values for different types of dairy foods owing to low amounts of consumption for some specific dairy products. Intake of skimmed or low fat milk was associated with slightly higher risk of total mortality and cardiovascular mortality, and with a lower risk of colorectal cancer. Whole milk intake was significantly associated with higher risk of total mortality, cardiovascular mortality, and cancer mortality, including lung cancer, ovarian cancer, and prostate cancer. The associations of whole milk intake and risk of total mortality and cancer mortality remained significant after correcting for multiple testing (P for trend <0.001). Cheese and yogurt intake was not associated with risk of total or cause specific mortality. Further adjustment for smoking pack years and other dietary factors such as red meat did not appreciably change the results.

Table 3

Hazard ratios of total mortality and mortality due to cardiovascular disease and cancer according to consumptions of dairy food subtypes, based on pooled analyses of cohorts from the Nurses’ Health Study, Nurses’ Health Study II, and Health Professionals Follow-up Study*†

View this table:

We estimated that substituting one serving per day of total dairy for nuts and legumes was associated with a 14% lower risk of total mortality (hazard ratio 0.86, 95% confidence interval 0.78 to 0.95; fig 2). Substituting one dairy serving per day for whole grains was associated with an 11% lower risk of total mortality (0.89, 0.84 to 0.93). However, a substitution for red and processed meat was also associated with a 5% higher risk of total mortality (1.05, 1.01 to 1.09).

Fig 2
Fig 2

Hazard ratios for total mortality and mortality due to cardiovascular disease and cancer associated with replacement of total dairy with other foods in the Nurses’ Health Study, Nurses’ Health Study II, and Health Professionals Follow-up Study. Error bars=95% confidence intervals. Multivariable adjusted model was further adjusted for family history of cancer (yes, no), family history of cardiovascular disease (yes, no), baseline disease status (hypertension, hypercholesterolemia), baseline body mass index (<20.9, 21.0-22.9, 23.0-24.9, 25.0-29.9, 30.0-34.9, ≥35.0), physical activity (<3.0, 3.0-8.9, 9.0-17.9, 18.0-26.9, ≥27.0 metabolic equivalent hours per week), overall dietary pattern (alternate healthy eating index score, in fifths), total energy intake (fifths), smoking status (Nurses’ Health Study and Health Professionals Follow-up Study: never, former (categorized into 1-4, 5-14, 15-24, 25-34, 35-44, ≥45 cigarettes/day, or unknown), current (categorized in to 1-4, 5-14, 15-24, 25-34, 35-44, ≥45 cigarettes/day, or unknown), alcohol consumption (0, ≥0-5, ≥5-10, ≥10-15, ≥15 g/day), postmenopausal status (yes, no; women only), and current postmenopausal hormone use (yes, no; women only)

In sensitivity analyses, we used time varying dairy intake updated every four years as exposure to provide more weight to the most recent diet. We found that dairy intake was not significantly associated with total mortality and mortality due to cancer and cardiovascular disease (table S2). To minimize reverse causation, we excluded participants who died during the first four years of follow-up, and results did not change (table S3). Significant interactions were found between dairy consumption and risk of mortality by body mass index (P value for interaction=0.03; table S4). A positive association was found among non-obese participants and no association was found among obese participants. We saw no significant differences in the associations between dairy consumption and risk of total mortality when stratified by disease status; cardiovascular drug treatment use; follow-up time; calendar year; age; smoking; diet quality; physical activity; and intakes of fat, calcium, and vitamin D (table S4).

Discussion

Principal findings

The present analysis of three large cohort studies involving 217 755 participants and 51 438 deaths did not support an inverse association between dairy intake and mortality. The robustness of our findings was confirmed by using time varying intake of dairy that was updated every four years as the exposure, excluding deaths in the first four years of follow-up, and the results were unchanged. A dose-response analysis suggested that two servings per day of dairy consumption was associated with lowest total and cardiovascular disease mortality. Substitution analyses suggested that replacing dairy with nuts or legumes and whole grains could lower mortality risk, but replacing total dairy with red and processed meat could increase mortality risk.

Comparison with other studies

The association between dairy intake and risk of mortality has been extensively investigated. Some early studies showed positive associations between dairy intake and risk of mortality,101112 while others showed inverse or null associations.1314151617 By summarizing these studies, one meta-analysis showed that overall dairy intake was not associated with risk of mortality.1819 More prospective cohort studies with large sample sizes have since been conducted worldwide but with conflicting conclusions. For example, in two Swedish cohorts including 106 772 participants, dairy intake was associated with higher risk of mortality,11 but in a Danish cohort of 74 241 participants and a Dutch cohort of 34 409 participants, dairy intake was not associated with higher risk of mortality.3233 In addition, in a Japanese cohort of 94 980 participants34 and in the Prospective Urban Rural Epidemiology study including 136 384 participants from 21 countries and five continents, dairy intake was associated with lower risk of mortality.3536

In our study, we found that total dairy intake was associated with slightly higher mortality, and the positive association remained even after correcting for multiple testing. The divergent findings across studies could be due to multiple reasons. For example, specific types of dairy foods varied across studies, the degree of residual confounding could have differed across studies because the covariates adjusted for were varied, and the effects of dairy foods could have depended on the nutritional quality of the background diet. However, the most important reason for the inconsistent findings could be the varied amount of dairy intake across populations—northern European countries consume the largest amount of dairy products worldwide whereas Asian countries such as Japan and China consume the least amount owing to high prevalence of lactose intolerance.3738

In our study, whole milk intake was associated with higher risk of cancer mortality, including ovarian cancer and prostate cancer mortality. The positive associations of whole milk intake and cancer remained significant after correcting for multiple testing. Previous studies also showed that whole milk intake was associated with higher risks of prostate and ovarian cancer.893940 The effects of milk consumption on plasma insulin-like growth factor 1,4142 which predicts higher risks of prostate and breast cancer,43 provides a plausible mechanism. Galactose, a lactose metabolite contained in milk, has been thought to increase risk of ovarian cancer, but no relation was seen in a pooled analysis.44

We also observed that whole milk intake was associated with higher risk of cardiovascular mortality. This finding could be explained by saturated fatty acids from whole milk affecting the blood lipid profile and promoting atherosclerosis. Results from short term intervention studies indicated that a diet higher in saturated fatty acids from whole milk and butter increased low density lipoprotein cholesterol when substituted for carbohydrates or unsaturated fatty acids.34 Furthermore, the Dietary Approaches to Stop Hypertension dietary pattern lowered systolic blood pressure, and the high content of low fat dairy products in that diet might have contributed to the effect.45 Observational studies also showed that saturated fat intake was associated with higher risk of coronary heart disease,46 although the effects depended on the substitution of unsaturated fats.47 One cohort study found that the replacement of dairy fat with polyunsaturated fat was significantly associated with lower risk of cardiovascular disease.48 Moreover, whole milk could be a marker of unhealthy diet and lifestyle.

Although we adjusted for a wide range of diet and lifestyle factors, residual confounding remains a possibility. However, further adjustment for smoking pack years and other dietary factors such as red meat did not appreciably change the results. As for skimmed or reduced fat milk, participants who developed hypertension, hypercholesterolemia, or cardiovascular disease could change to adopt a healthier lifestyle by consuming more skimmed or reduced milk instead of whole milk. Although we stopped updating dietary intakes after the development of cardiovascular disease and cancer, reverse causation could remain and potentially explain a slightly positive association observed between skim milk intake and risk of cardiovascular mortality. Reverse causation might also be the reason for non-significant association between time-varying dairy intake and risk of total mortality in the sensitivity analysis.

Our study found that intake of skimmed or low fat milk was associated with lower risk of colorectal cancer mortality, and calcium contained in dairy foods might have an anti-carcinogenic effect. In previous meta-analyses and pooled analyses, milk consumption was inversely associated with risk of colorectal cancer, possibly because of its high calcium content.54950 Many prospective cohort studies showed that calcium intake was associated with lower risk of colorectal cancer,405152 and clinical trials also showed that calcium supplementation reduced the recurrence of colorectal adenoma.53 Animal studies showed that calcium can suppress colonic epithelial cell proliferation, induce apoptosis in colonic epithelium cells, and inhibit mutations of the colonic k ras gene.5455

Our finding that dairy intake was not associated with lower risk of mortality was also consistent with mendelian randomisation studies using lactase persistence genotype rs4988235 as the instrumental variable. Two studies conducted in a Danish population of 97 811 participants found that the T allele of rs4988235 was not associated with lower risks of type 2 diabetes, total mortality, and mortality due to cardiovascular disease and cancer.5657 Two studies conducted within the CHARGE consortium involving more than 171 000 participants showed that the T allele of rs4988235 was associated with slightly higher body mass index and systolic blood pressure, two important risk factors of mortality.5859

We conducted substitution analyses to predict the effects of replacing dairy products with other foods on mortality. We found that replacing dairy with nuts or legumes, poultry, and whole grains was associated with lower mortality. However, replacing total dairy with red and processed meat was associated with increased mortality. These data suggest that dietary recommendations on dairy intake should consider healthier alternatives. Moreover, the health effects of dairy intake could depend on the background diet. For example, in populations with high carbohydrate intake and poor quality diets, dairy foods and other protein sources might add important nutritional value.

Strengths and limitations of study

Our study has several strengths. Firstly, dairy intake was measured repeatedly in all three cohorts, and we calculated cumulative averages for dairy intake to minimize the random measurement error caused by within-person variation. Furthermore, we stopped updating dairy intake on the diagnosis of potential intermediate conditions, which minimized the potential for reverse causation due to change of lifestyle factors because of chronic diseases. Secondly, our study included 217 755 participants with 51 438 deaths. The large sample size and number of deaths provided sufficient power for the analyses.

Several limitations also need to be considered. Firstly, the NHS, NHSII, and HPFS are predominantly white healthcare professionals, which could limit the generalizability of the findings to other nationalities and races. However, their occupations are a distinct advantage that allows us to collect high quality data using self reported questionnaires and enhance the internal validity of the study by reducing confounding. Secondly, given the observational study design, we could not directly establish a causal relation between dairy intake and mortality.

Conclusions and policy implications

In conclusion, our data did not support an inverse association between high amount of dairy consumption and risk of mortality. A dose-response analysis suggested that two servings per day of dairy consumption was associated with lowest cardiovascular mortality, but higher intake was associated with a slightly higher mortality, especially cancer mortality. Substitution analyses suggested that the health effects of dairy products could depend on the comparison foods used to replace dairy.

What is already known on this topic

  • The relation between dairy intake and various health outcomes including type 2 diabetes, cardiovascular disease, and cancer have been extensively examined, and most studies have shown no appreciable beneficial or adverse associations

  • However, evidence on the association between dairy intake and mortality from prospective cohort studies is more limited

What this paper adds

  • Total dairy intake was not associated with lower risk of total mortality

  • The health effects of dairy could depend on the comparison foods used to replace dairy

  • Slightly higher cancer mortality was non-significantly associated with dairy consumption, but warrants further investigation

Acknowledgments

We thank the participants and staff of the NHS, NHSII, and HPFS cohorts for their valuable contributions as well as the following state cancer registries for their help: Alabama, Arizona, Arkansas, California, Colorado, Connecticut, Delaware, Florida, Georgia, Idaho, Illinois, Indiana, Iowa, Kentucky, Louisiana, Maine, Maryland, Massachusetts, Michigan, Nevada, New Hampshire, New Jersey, New York, North Carolina, North Dakota, Ohio, Oklahoma, Oregon, Pennsylvania, Rhode Island, South Carolina, Tennessee, Texas, Virginia, Washington, and Wyoming.

Footnotes

  • Contributors: MD, JL, LQ, CE, XZ, JEM, MS, JEC, KMR, PK, DC, WCW, and FBH designed the study and collected data. MD mainly conducted analyses and wrote the manuscript. FBH supervised the data analysis and reviewed and edited the manuscript. All authors contributed substantially to the interpretation of data and the drafting or critical revision of the manuscript for important intellectual content. The authors assume full responsibility for analyses and interpretation of these data. MD and FBH are the guarantors. The corresponding author attests that all listed authors meet authorship criteria and that no others meeting the criteria have been omitted.

  • Funding: The NHS, NHSII, and HPFS are supported by grants (UM1 CA186107, P01 CA87969, R01 CA49449, R01 HL034594, R01 HL088521, UM1 CA176726, R01 CA67262, UM1 CA167552, R01 HL35464, R01 HL60712) from the National Institutes of Health. The funding sources did not participate in the design or conduct of the study; collection, management, analysis or interpretation of the data; or preparation, review, or approval of the manuscript.

  • Competing interests: All authors have completed the ICMJE uniform disclosure form at www.icmje.org/coi_disclosure.pdf and declare: support from the National Institutes of Health for the submitted work; FBH reports support from grants HL60712, HL118264, and DK112940 from the National Institutes of Health, research support from the California Walnut Commission, honorariums for lectures from Metagenics and Standard Process, and honorariums from Diet Quality Photo Navigation, outside the submitted work; the remaining authors report no other relationships or activities that could appear to have influenced the submitted work.

  • Ethical approval: The study protocol was approved by the institutional review boards of the Brigham and Women’s Hospital and the Harvard T H Chan School of Public Health. The completion of the self administered questionnaire was considered to imply informed consent.

  • Data sharing: No additional data available.

  • The lead authors affirm that the manuscript is an honest, accurate, and transparent account of the study being reported; that no important aspects of the study have been omitted; and that any discrepancies from the study as planned (and, if relevant, registered) have been explained.

This is an Open Access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/4.0/.

References