Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Dairy Products Consumption and Risk of Type 2 Diabetes: Systematic Review and Dose-Response Meta-Analysis

  • Dengfeng Gao ,

    gaomedic@163.com

    Affiliations Department of Cardiology, The Second Affiliated Hospital, Xi'an Jiaotong University School of Medicine, Xi'an, Shaanxi, P.R. China, Key Laboratory of Environment and Genes Related to Diseases (Xi'an Jiaotong University), Ministry of Education, Xi'an, Shaanxi, P. R. China

  • Ning Ning,

    Affiliation Department of Cardiology, The Second Affiliated Hospital, Xi'an Jiaotong University School of Medicine, Xi'an, Shaanxi, P.R. China

  • Congxia Wang,

    Affiliation Department of Cardiology, The Second Affiliated Hospital, Xi'an Jiaotong University School of Medicine, Xi'an, Shaanxi, P.R. China

  • Yuhuan Wang,

    Affiliation Department of Endocrinology, The Second Affiliated Hospital, Xi'an Jiaotong University School of Medicine,Xi'an, Shaanxi, P.R. China

  • Qing Li,

    Affiliation Department of Cardiology, The Second Affiliated Hospital, Xi'an Jiaotong University School of Medicine, Xi'an, Shaanxi, P.R. China

  • Zhe Meng,

    Affiliation Department of Cardiology, The Second Affiliated Hospital, Xi'an Jiaotong University School of Medicine, Xi'an, Shaanxi, P.R. China

  • Yang Liu,

    Affiliation Department of Cardiology, The Second Affiliated Hospital, Xi'an Jiaotong University School of Medicine, Xi'an, Shaanxi, P.R. China

  • Qiang Li

    Affiliation Department of Epidemiology, Xi'an Jiaotong University School of Medicine, Xi'an, Shaanxi, P.R. China

Abstract

Background

The consumption of dairy products may influence the risk of type 2 diabetes mellitus (T2DM), but inconsistent findings have been reported. Moreover, large variation in the types of dairy intake has not yet been fully explored.

Methods and Results

We conducted a systematic review and meta-analysis to clarify the dose–response association of dairy products intake and T2DM risk. We searched PubMed, EMBASE and Scopus for studies of dairy products intake and T2DM risk published up to the end of October 2012. Random-effects models were used to estimate summary relative risk (RR) statistics. Dose-response relations were evaluated using data from different dairy products in each study. We included 14 articles of cohort studies that reported RR estimates and 95% confidence intervals (95% CIs) of T2DM with dairy products intake. We found an inverse linear association of consumption of total dairy products (13 studies), low-fat dairy products (8 studies), cheese (7 studies) and yogurt (7 studies) and risk of T2DM. The pooled RRs were 0.94 (95% CI 0.91–0.97) and 0.88 (0.84–0.93) for 200 g/day total and low-fat dairy consumption, respectively. The pooled RRs were 0.80 (0.69–0.93) and 0.91 (0.82–1.00) for 30 g/d cheese and 50 g/d yogurt consumption, respectively. We also found a nonlinear association of total and low-fat dairy intake and T2DM risk, and the inverse association appeared to be strongest within 200 g/d intake.

Conclusion

A modest increase in daily intake of dairy products such as low fat dairy, cheese and yogurt may contribute to the prevention of T2DM, which needs confirmation in randomized controlled trials.

Introduction

The prevalence of type 2 diabetes mellitus (T2DM) is a growing public-health burden worldwide, particularly in developing countries. The prevalence of T2DM is estimated to reach 552 million worldwide by 2030 [1]. T2DM may cause substantial morbidity and mortality and is associated with enormous economic, health, and societal costs [2], [3]. Moreover, as compared with unaffected people, those with T2DM are at increased risk of other chronic illnesses, including cardiovascular disease; T2DM more than doubles the risk of a heart attack or stroke [4], [5]. Therefore, the identification of modifiable risk factors for primary prevention of T2DM is of considerable public health importance.

T2DM has genetic components but is also directly influenced by modifiable lifestyle factors, including dietary behaviors [6]. Dairy consumption might affect T2DM. Experimental studies indicated that dairy protein, such as whey protein, has insulinotropic and glucose-lowering properties [7].The Multi-Ethnic Study Atherosclerosis [8] and Cardiovascular Health study [9] suggested that fatty acids in dairy might be responsible for lower risk of T2DM. Epidemiological studies of dairy products and T2DM risk have given mixed results [10], [11], [12], [13], [14], [15], [16], [17], [18], [19], [20], [21], [22], [23]. Some cohort studies have reported inverse associations of intake of total and low-fat dairy products, milk and/or yogurt and T2DM risk, but other studies found no association [10], [11], [19], [20], [23]. One meta-analysis of 7 studies reported a significant inverse association of dairy intake and risk of T2DM [24]. However, the large variation in types of dairy consumed has not been fully explored. Furthermore, the dose–response relationship needs to be clarified as well as any gender or geographic differences in the T2DM risk. In addition, possible confounding by other lifestyle factors needs to be explored to firmly establish the potential preventive role of dairy products in T2DM.

We conducted a meta-analysis of population-based cohort studies to investigate dose–response associations of consumption of total, low-fat, and full-fat dairy products as well as different types of dairy products and risk of T2DM.

Methods

Data Sources and Search Strategy

We followed standard criteria for conducting and reporting meta-analyses of observational studies (MOOSE). Two authors (DG and NN) independently did a literature search MEDLINE via PubMed (published from 1966 to March 2013), EMBASE (published from 1980 to March 2013), and Scopus (www.scopus.com) with no restriction on language. To identify studies of milk or dairy product intake and T2DM risk, we used both the medical subject heading (MeSH) terms (“Diabetes Mellitus” AND (milk OR dairy)) and searched the text using the terms (‘diabetes’/exp OR diabetes') AND (‘dairy’/exp OR dairy OR ‘milk’/exp OR milk). We also searched the reference lists of all studies retrieved and published systematic reviews and meta-analysis.

Study Selection

All abstracts retrieved were examined independently by 2 investigators (DG and NN) who then retrieved the full text of potential articles. Disagreements were resolved by consensus, and if necessary, with a third author (CW). We included prospective cohort studies and case-cohort studies assessing the association of consumption of total dairy products or specific types of dairy products and T2DM. To be included in the analyses, articles needed to contain estimates of the relative risk (RR) (such as odds ratios [ORs], hazard ratios [HRs] or risk ratios) with 95% confidence intervals (95% CIs). We excluded animal studies, clinical trials, cross sectional studies, case-control studies, and studies that examined other associations. For the dose–response analysis, a quantitative measure of intake had to be provided. If the article lacked data, we attempted to contact the author.

Data Extraction and Quality Assessment

We extracted the following data from each study: country where the study was conducted, follow-up period, sample size, gender, age, number of cases, dietary assessment method (type, number of food items and whether the food intake had been validated), type of dairy product (e.g., total dairy, milk, cheese), quantity of intake, HRs, RR values, and ORs and 95% CIs for dairy product intake and, when available, the number of cases and participants or person-years for each category of dairy product consumption. Two authors (YL and ZM) independently performed the data extraction. Any disagreements were resolved by discussion.

Two independent reviewers (DG and NN) evaluated the quality of the selected studies by using a modified scoring system that was based on a recently used system (designed with reference to QUATSO [25], MOOSE [26], and STROBE [27]) that allowed for a total score of 0 to 6 points (6 indicating highest quality) [28]. The system allocates one point each for 1) any justification given for the cohort; 2) appropriate inclusion and exclusion criteria used; 3) outcome (diagnosis of T2DM not solely based on self-reporting); 4) intervention (participants' usual dairy consumption assessed with a validated tool); 5) statistical analysis (adjustments made for age, sex, body mass index, and family history of T2DM, total energy intake and physical activity, these being proven risk factors for type 2 diabetes); and 6) any other adjustments performed (such as glycemic load and dietary factors).

Statistical analysis

HRs and RRs were assumed to be approximately the same measure of relative risk. For articles reporting ORs, we estimated the RRs from the ORs using a previously published correction method [29]. To take into account heterogeneity between studies, we used a random-effects models to calculate summary RRs and 95% CIs for the highest versus lowest level of dairy product intake and for the dose–response analysis. The natural logarithm of the RR from each study was weighted by the inverse of its variance and pooled across studies. A two-tailed P<0.05 was considered statistically significant. Articles that reported findings for men and women separately were considered 2 studies when the observed items were combined.

For the dose–response analysis, we used GLST command in Stata software as the method proposed by Greenland and Longnecker [30] and Orsini et al. [31] to compute study-specific slopes (linear trends) and 95% CIs from the natural logs of the RRs and 95% CIs across categories of dairy product intake.

For each study, the median or mean level of dairy product intake for each category was assigned to each corresponding RR. When the median or mean intake per category was not provided, we assigned the midpoint of upper and lower boundaries in each category as the average intake. If the highest or the lowest category was open-ended, we assumed that the open-ended interval length had the same length as the adjacent interval. If the intake was reported in densities (i.e., per 1000 kcal), we recalculated the reported intake as absolute intake using the mean or median energy intake reported in the article [14]. When studies reported the intake in servings and times per day or week, we converted the intake to grams of intake per day using standard units of 244 g (or 244 ml) for milk, 43 g for cheese (2 slices) and 177 g for total dairy products from the serving sizes reported in the US Department of Agriculture Food and Nutrient Database for Dietary Studies [32]. Pooled estimates were expressed in rounded numbers that approximated a normal portion size and fitted within the range of dairy intake of all studies (i.e., 200 g for milk and total, low-fat, and full-fat dairy; 50 g for yogurt; and 30 g for cheese).

To examine a potential nonlinear association between dairy products intake and T2DM risk, we performed a 2-stage, random-effects, dose-response meta-analysis, as recently summarized [33]. In the first stage, we constructed study-specific restricted cubic spline models, with 4 knots at fixed percentiles (5%, 35%, 65%, 95%) of the exposure distribution by using generalized least-squares regression. In the second stage, we combined the 2 regression coefficients and the variance/covariance matrix that had been estimated within each study, using the restricted maximum likelihood method in a multivariate random-effects meta-analysis. The pooled relative risks for specific exposure values were then estimated. A P value for nonlinearity was calculated by testing the null hypothesis that the coefficient of the second spline was equal to zero.

Heterogeneity among studies was assessed by I2, the amount of total variation explained by the between-study variation, and the Q test. We conducted subgroup and random effects univariate and multivariate meta-regression to investigate potential sources of heterogeneity we performed for the primary outcomes. Publication bias was assessed with funnel plots, Begg's test and Egger's test. Stata v12.0 (Stata Corp, College Station, TX) was used for all the statistical analysis.

Results

Study characteristics

We included 15 prospective cohort studies and 1 case–cohort study in our analysis (Figure 1). 6 of the studies [12], [13], [14], [15], [16], [18] were performed in the United States, 6 in Europe [11], [19], [21], [23], [34], [35], 2 in Asia [17], [22] and 2 in Australia [10], [20]. The articles were published between 2005 and 2013 and included 526,998 subjects (including 29,789 T2DM cases). Characteristics of included studies are in Table 1. Figures 2 show assessments by risk of bias. The studies were generally of moderate quality. More than 75% of the studies met 4 of the quality items as reported and 8 studies met all requirements. Interrater reliability for assessing quality items was good (κ = 0.86, P<0.01).

thumbnail
Figure 1. Flow chart for the selection of studies for meta-analysis of the association of dairy products intake and type 2 diabetes (T2DM).

https://doi.org/10.1371/journal.pone.0073965.g001

thumbnail
Figure 2. Methodological quality across included studies.

https://doi.org/10.1371/journal.pone.0073965.g002

thumbnail
Table 1. Characteristics of the cohort studies of dairy products intake and type 2 diabetes mellitus (T2DM).

https://doi.org/10.1371/journal.pone.0076613.t001

Total Dairy Products Intake and T2DM Risk

In all, 13 studies [10], [11], [12], [13], [14], [15], [16], [17], [18], [19], [20], [23], [34] including 457,893 subjects (27,095 cases) were analyzed.

In all, 13 studies(8–18; 21) including 457,893 subjects (27,095 cases) were analyzed.

High versus low intake.

The summary RR for all studies was 0.89 (95% CI 0.81–0.98), with moderate heterogeneity, I2 = 65.4% and Pheterogeneity = 0.000 (Figure 3A).

thumbnail
Figure 3. Forest plot of relative risk (RR) for total dairy products intake and T2DM.

A, highest versus lowest intake. B, dose–response analysis (200 g/d). Weights are from random effects analysis.

https://doi.org/10.1371/journal.pone.0073965.g003

Dose–response analysis.

The summary RR for an increase of 200 g/day was 0.94 (0.91–0.97), with moderate heterogeneity, I2  = 51.6%, Pheterogeneity  = 0.02 (Figure 3B). On subgroup analysis (Table 2), we found an inverse association of total dairy intake and T2DM risk in all strata except European studies and studies not adjusting for family history of T2DM, although in some analyses the associations were not statistically significant. None of the results differed significantly by sex (P = 0.21 for all comparisons). On univariate meta-regression analysis, geographic location, adjustment for family T2DM history, and glycemic load were significant predictors of the heterogeneity (p = 0.05, p = 0.04 and p = 0.04, respectively). But on multivariate meta-regression, we failed to identify the source of heterogeneity. We found no evidence of publication bias by Egger's test (P = 0.37), Begg's test (P = 0.58) or funnel plot(see Appendix Figure 1). We also found a nonlinear association of total dairy product intake and T2DM risk, Pfor nonlinearity <0.001, with most of the risk reduction occurring with intake up to about 200 g/d; higher intake were associated with a further but more modest decrease in risk (Figure 4A).

thumbnail
Figure 4. Dairy products and incidence of T2DM, nonlinear dose–response analysis.

A, total dairy. B, low-fat dairy.

https://doi.org/10.1371/journal.pone.0073965.g004

thumbnail
Table 2. Subgroup analyses of total and low-faty dairy products intake and T2DM, dose–response analysis.

https://doi.org/10.1371/journal.pone.0076613.t002

Low- and Full-fat Dairy Intake and T2DM Risk

8 studies [10], [11], [12], [13], [14], [15], [16], [23] including 260,700 subjects (9,398 cases) were analyzed.

High versus low intake.

Low-fat dairy consumption was inversely associated with T2DM risk, with a pooled RR of 0.81 (0.74–0.89) (Figure 5A). Full-fat dairy consumption was not associated with T2DM risk, with a summary RR of 0.95 (0.85–1.07) (Figure 6A). We found no significant heterogeneity for the associations of low-fat (I2 = 1.8%; Pheterogeneity  = 0.42) or full-fat dairy consumption (I2 = 38.1%; Pheterogeneity  = 0.13).

thumbnail
Figure 5. Forest plot of RR for low-fat dairy products intake and T2DM.

A, highest versus lowest intake. B, dose–response analysis (200 g/d). Weights are from random-effects analysis.

https://doi.org/10.1371/journal.pone.0073965.g005

thumbnail
Figure 6. Forest plot of RR for full-fat dairy products intake and T2DM.

A, highest versus lowest intake. B, dose–response analysis (200 g/d). Weights are from random-effects analysis.

https://doi.org/10.1371/journal.pone.0073965.g006

Dose–response analysis.

The summary RR for a 200-g/day increase in low-fat dairy intake was 0.88 (0.84–0.93), with no evidence of heterogeneity, I2 = 16.3% and Pheterogeneity  = 0.32 (Figure 5B). The summary RR for a 200-g/day increase in full-fat dairy intake was 0.95 (0.88–1.04), with evidence of heterogeneity, I2 = 52.2% and Pheterogeneity  = 0.04 (Figure 6B). We found an inverse association of low-fat dairy intake and T2DM risk for all strata, although in some analyses the associations were not statistically significant (Table 2). On univariate meta-regression analysis, the effect was weaker, although not significantly, for studies with than without adjustment for coffee and meat intake (P = 0.06 for both). But on multivariate meta-regression, we failed to identify the source of heterogeneity. On We found a nonlinear association of low-fat dairy intake and T2DM risk, Pfor nonlinearity  = 0.02, with most of the risk reduction occurring with intake up to about 300 g/day; higher intake (>400 g/day) was not associated with a further decrease in risk (Figure 4B).

Milk Intake and T2DM Risk

High versus low analysis.

9 studies [11], [13], [15], [17], [18], [19], [20], [22], [23] including 327,039 subjects (21,755 cases) were analyzed. For 6 studies [11], [17], [18], [19], [22], [23], data were available on the association of total milk intake and T2DM risk, for 3 studies [13], [15], [20], milk consumption was analyzed as full-fat milk (or whole milk) and low-fat milk (or skim milk) and for 1 study, data was reported as full fat milk only [35]. The summary RR for total milk intake was 0.89 (0.78–1.01), with moderate heterogeneity, I2 = 51.8% and P heterogeneity = 0.043. The summary RR for low-fat and full-fat milk intake was 0.82 (0.69–0.97, I2 = 40% and Pheterogeneity = 0.19) and 1.12 (0.99–1.27, I2 = 0% and Pheterogeneity = 0.79), respectively (Figure 7).

thumbnail
Figure 7. Forest plot of RR for highest versus lowest milk intake and T2DM.

A, total milk. B, low fat milk. C, full fat milk. Weights are from random-effects analysis.

https://doi.org/10.1371/journal.pone.0073965.g007

Dose–response analysis.

8 studies [11], [13], [15], [17], [18], [19], [20], [22], [23] were analyzed. The summary RR for a 200-g/day increase in total milk intake was 0.89 (0.79–1.01), with evidence of moderate heterogeneity, I2 = 66.3% and Pheterogeneity  = 0.005 (data not shown). The summary RR for a 200-g/day increase in full- and low-fat milk intake was 1.27 (0.97–1.67, I2 = 0% and Pheterogeneity = 0.58) and 0.83 (0.70–1.00, I2 = 14% and Pheterogeneity = 0.21), respectively.

Yogurt and Cheese Intake and T2DM Risk

Seven studies [11], [12], [13], [15], [17], [19], [20] including 254,552 subjects (18,532 cases) were analyzed for yogurt intake. Seven studies [11], [13], [15], [17], [19], [20], [23] including 178,429 subjects (14,810 cases) were analyzed for cheese intake.

High versus low intake.

Yogurt and cheese intake were inversely associated with T2DM risk. The pooled RRs were 0.85 (0.75–0.97, I2 = 55% and Pheterogeneity  = 0.02) and 0.82 (0.77–0.87, I2 = 0% and Pheterogeneity  = 0.82), respectively.

Dose–response analysis.

Yogurt and cheese intake were inversely associated with T2DM incidence. The pooled RRs were 0.91 (0.82–1.00, I2 = 74%, Pheterogeneity = 0.001) per 50 g/d and 0.80 (0.69–0.93, I2 = 59%, Pheterogeneity = 0.02) per 30 g/d, respectively.

Other dairy products

Intake of other types of dairy products except ice cream (n = 2 studies) were not significantly associated with T2DM risk (0.84 [0.73–0.95]). The pooled RR for total fermented dairy intake (n = 3 studies) was 0.94 (0.75–1.18) and cream (n = 2 studies) was 0.96 (0.84–1.12).

Discussion

This meta-analysis showed that total dairy intake was associated with a 6% lower risk of T2DM per 200 g/day consumption. When examining different types of dairy products in relation to T2DM risk, we found significant inverse associations of intake of low-fat dairy, low-fat milk, cheese and yogurt and T2DM. We found no association of intake of full-fat dairy as well as total and full-fat milk and T2DM. We also clarified a nonlinear association of both total and low-fat dairy intake and incidence of T2DM.

The hypothesis that dairy products intake protects against T2DM has received much interest among medical professionals and the general population. In intervention studies, the Dietary Approaches to Stop Hypertension (DASH) diet (a dietary pattern focusing on low-fat milk and other dairy products) increased high-density lipoprotein levels, reduced triglycerides levels, reduced blood pressure (both systolic and diastolic), contributed to weight loss, and reduced fasting blood glucose in both men and women as compared with the control diet [36]. In epidemiological studies, the association of dairy products intake and T2DM has been explored with inconsistent results [37].

Our findings for high versus low dairy intake are consistent with results from previous meta-analyses [24], which only included 7 studies. High versus low analyses are limited because true differences in the level and range of intake between studies are not considered and may contribute to heterogeneity in the results. With the accumulated evidence, we were able to enhance the precision of the risk estimates, perform dose–response analyses of different dairy products and explore the shape of the dose–response curve and sources of heterogeneity, thereby increasing the clinical relevance of our findings [38].

In addition, the presence of both linear and nonlinear dose–response relationships of specific dairy products strengthened the findings of an association of dairy products intake and risk of T2DM. In the linear dose-response analysis, we found a 6% and 12% lower risk of T2DM per 200 g/day intake of total and low-fat dairy products, respectively. Furthermore, we discovered a potential nonlinear association of total and low-fat dairy products intake and T2DM. A low threshold of 200 g/day total dairy and 300 g/day low fat dairy may reduce the risk by about 10% or 15% respectively. Intake above that level seems to have further but modest additional benefit for T2DM risk.

Dairy is a major source of dietary calcium and magnesium, 2 minerals that have a role in the development of T2DM, for potential in improving pancreatic B-cell function and insulin sensitivity [39]. Experimental [39], prospective cohort studies [40], [41] and a recent meta-analysis [42] have provided convincing evidence to support the direct effects of calcium and magnesium intake on insulin resistance and T2DM. In this study, we found that the association of dairy intake and T2DM risk remained unchanged after adjusting for diet calcium and/or magnesium (7 studies), so other major components in dairy products could account for the association. Recently, the beneficial physiological effects of dairy protein, such as whey protein, on the control of food intake and glucose metabolism have been reported. Studies have shown the insulinotropic and glucose-lowering properties of whey protein in healthy and T2DM subjects [43]. Furthermore, in addition to milk proteins, trans-palmitoleate, obtained primarily from dairy intake, is associated with reduced incidence of diabetes [9].

Our analysis of high- and low-fat dairy products revealed an inverse association of only low-fat dairy food intake and T2DM risk. This support the present recommendations by health authorities and governments to eat low-fat rather than full-fat dairy products [44]. We think the most prominent relationship was from residual confounding by factors related to a more unhealthy diet or lifestyle. On the other hand, we can not rule out the association between the intake of saturated fatty acid (SFA). Dairy products contributed to 15% of the total dietary SFA intake [45]. Although prospective cohorts demonstrate no significant association between SFA intake and risk of T2DM, some findings from experimental and observational studies have showed that SFA intake was inversely associated with insulin sensitivity [45], [46], [47]. Finally, the likelihood of publication bias effects may cause uncertain results. For analysis of the milk products, only 3 of 14 studies separately evaluated whole vs. low-fat milk, and thus it seems that publication bias could account for the observed difference between low vs. whole fat milk. Furthermore, because cheese, even low-fat cheese, has higher fat and saturated fat than whole milk yet was still associated with lower risk, it appears less likely that the observed difference between whole fat and low fat milk would be due to higher fat or saturated fat content in whole milk. Further confirmatory results of appropriately powered studies are still needed.

Cheese, which has far more fat than whole-fat milk, more than half of which is saturated Evidence suggests that saturated fat intake has an adverse effect on insulin sensitivity and increases the risk of T2DM. In our analysis, we found inverse association between both cheese and yogurt intake and incidence of T2DM. The exact mechanisms responsible for the significant inverse association between cheese and yogurt and T2DM are unknown. It could be partly explained by the fact that both dairy subgroups are a good source for vitamin K2. Vitamin K2 is exclusively synthesized by bacteria and is therefore only present in fermented dairy products such as cheese and yogurt due to the bacterial starter fermentation [48]. Vitamin K2 has recently been linked to a reduced risk of T2DM [49]. Additionally, these dairy subcategories are particularly high in the fat-soluble vitamin D, which has been found to be inversely associated with T2DM [50], [51].

We did not find a consistent pattern of difference or heterogeneity in results by sex or any other study characteristics examined, except for geographic location, which significantly modified the association between total dairy products intake and T2DM risk. We found a significant inverse association among US studies, with no evidence of a protective effect of total dairy food intake in European or Asia studies. This may be a chance finding, because only 3 European studies and 1 Asian study were included in this subgroup analysis or could be due to other factors. As well, differences in the ranges of intake or intake in the referent category could explain these results. Because of the nonlinear association between total dairy food intake and T2DM risk with the strongest reduction at low levels of intake, some studies may have missed an effect because the intake in the referent category may have been already sufficient to reduce risk. For example, in some European studies, intake of total dairy food in the referent category was >200 g/d but was <200 g/d for all US studies. As well, types of dairy food intake may vary between populations. In addition, differences in study size and follow-up time may contribute to the variations. Further cohort studies of specific dairy products and T2DM risk in different populations are needed.

Our meta-analysis contains some limitations. Publication bias is a major concern for analyses that depend on only a few studies. For example, in our analysis, only 4 of the 15 studies separately evaluated full or low-fat milk. So the efficiency of analysis on different milk production was limited. The inverse association we found between dairy products intake and T2DM risk could be due to unmeasured or residual confounding. Higher intake of dairy products, especially low-fat dairy products, is often associated with other lifestyle factors, including increased physical activity, low prevalence of smoking, and overweight/obesity, although different types of dairy products may be differentially associated with some of these confounders. In addition, the results were generally similar in the subgroup analyses when we stratified results by adjustment for confounding factors or other study characteristics, with no heterogeneity between subgroups for total and low-fat dairy product consumption. Only the analysis of total dairy products revealed some indication of heterogeneity, with studies that adjusted for family history of T2DM showing an inverse association with T2DM; studies that did not adjust for family history of T2DM showed a nonsignificant positive association, which suggests potential confounding.

Measurement errors in the assessment of dietary intake are known to bias effect estimates. Our results are based on data from cohort studies, in which dairy intake was mostly assessed by food-frequency questionnaires. In several studies, validation of the food-frequency questionnaires showed good correlations, of ≈0.6–0.7 for milk or (if not assessed) for protein and calcium, which are good indicators for milk intake. However, we cannot exclude that measurement errors might have resulted in attenuated associations. Dietary changes after baseline can also attenuate associations of dietary intake and T2DM risk; however, only 7 of the included studies [12], [13], [14], [17], [18], [22], [23] used repeated assessments of diet, and the results were similar when using only the baseline questionnaire for the analyses (data not shown). Furthermore, dietary intake data were collected between 1984 and 2003. In earlier studies, full-fat dairy was a major contributor to total dairy intake, whereas in later studies intake was more often low-fat dairy and publication year may have explained the study heterogeneity (p = 0.02). Finally, because all the studies were conducted primarily among middle-aged and older people, these results might not be generalizable to dairy intake in earlier life periods, which might have similar or different effects.

In conclusion, our results suggest a inverse association of intake of dairy products, such as low-fat dairy, cheese and yogurt and T2DM risk. Further cohort studies are warranted to investigate the specific types of dairy products in the association, the impact of measurement errors on estimates, any gender-specific recommendations, and biomarkers of dairy intake.

Supporting Information

Figure S1.

Funnel plot of assessing evidence of publication bias. A. For total dairy; B. For low-fat dairy; C. For full-fat dairy.

https://doi.org/10.1371/journal.pone.0073965.s001

(TIF)

Figure S2.

Forest plot of RR for highest versus lowest yogurt and cheese intake and T2DM. A, yogurt. B, cheese. Weights are from random-effects analysis.

https://doi.org/10.1371/journal.pone.0073965.s002

(TIF)

Acknowledgments

We thank Nita G. Forouhi, Dianna J. Magliano and Karen L. Margolis for providing additional data for the meta-analysis. We thank Nicola Orsini for additional help in statistical works. We thank Laura Smales from the http://www.biomedediting.com for manuscript editing.

Author Contributions

Conceived and designed the experiments: DFG YHW CXW. Performed the experiments: DFG NN. Analyzed the data: YL ZM. Contributed reagents/materials/analysis tools: DFG QL QL. Wrote the paper: DFG NN.

References

  1. 1. Guariguata L, Whiting D, Weil C, Unwin N (2011) The International Diabetes Federation diabetes atlas methodology for estimating global and national prevalence of diabetes in adults. Diabetes Res Clin Pract 94: 322–332.
  2. 2. Polonsky KS (2012) The past 200 years in diabetes. N Engl J Med 367: 1332–1340.
  3. 3. Straka RJ, Liu LZ, Girase PS, DeLorenzo A, Chapman RH (2009) Incremental cardiovascular costs and resource use associated with diabetes: an assessment of 29,863 patients in the US managed-care setting. Cardiovasc Diabetol 8: 53.
  4. 4. Woodward M, Zhang X, Barzi F, Pan W, Ueshima H, et al. (2003) The effects of diabetes on the risks of major cardiovascular diseases and death in the Asia-Pacific region. Diabetes Care 26: 360–366.
  5. 5. Huxley R, Barzi F, Woodward M (2006) Excess risk of fatal coronary heart disease associated with diabetes in men and women: meta-analysis of 37 prospective cohort studies. BMJ 332: 73–78.
  6. 6. Lakerveld J, Bot SD, Chinapaw MJ, van Tulder MW, van Oppen P, et al. (2008) Primary prevention of diabetes mellitus type 2 and cardiovascular diseases using a cognitive behavior program aimed at lifestyle changes in people at risk: Design of a randomized controlled trial. BMC Endocr Disord 8: 6.
  7. 7. Jakubowicz D, Froy O (2013) Biochemical and metabolic mechanisms by which dietary whey protein may combat obesity and Type 2 diabetes. J Nutr Biochem 24: 1–5.
  8. 8. Mozaffarian D, de Oliveira Otto MC, Lemaitre RN, Fretts AM, Hotamisligil G, et al. (2013) trans-Palmitoleic acid, other dairy fat biomarkers, and incident diabetes: the Multi-Ethnic Study of Atherosclerosis (MESA). Am J Clin Nutr 97: 854–861.
  9. 9. Mozaffarian D, Cao H, King IB, Lemaitre RN, Song X, et al. (2010) Trans-palmitoleic acid, metabolic risk factors, and new-onset diabetes in U.S. adults: a cohort study. Ann Intern Med 153: 790–799.
  10. 10. Louie JC, Flood VM, Rangan AM, Burlutsky G, Gill TP, et al.. (2012) Higher regular fat dairy consumption is associated with lower incidence of metabolic syndrome but not type 2 diabetes. Nutr Metab Cardiovasc Dis.
  11. 11. Soedamah-Muthu SS, Masset G, Verberne L, Geleijnse JM, Brunner EJ (2012) Consumption of dairy products and associations with incident diabetes, CHD and mortality in the Whitehall II study. Br J Nutr: 1–9.
  12. 12. Margolis KL, Wei F, de Boer IH, Howard BV, Liu S, et al. (2011) A diet high in low-fat dairy products lowers diabetes risk in postmenopausal women. J Nutr 141: 1969–1974.
  13. 13. Liu S, Choi HK, Ford E, Song Y, Klevak A, et al. (2006) A prospective study of dairy intake and the risk of type 2 diabetes in women. Diabetes Care 29: 1579–1584.
  14. 14. Malik VS, Sun Q, van Dam RM, Rimm EB, Willett WC, et al. (2011) Adolescent dairy product consumption and risk of type 2 diabetes in middle-aged women. Am J Clin Nutr 94: 854–861.
  15. 15. Choi HK, Willett WC, Stampfer MJ, Rimm E, Hu FB (2005) Dairy consumption and risk of type 2 diabetes mellitus in men: a prospective study. Arch Intern Med 165: 997–1003.
  16. 16. van Dam RM, Hu FB, Rosenberg L, Krishnan S, Palmer JR (2006) Dietary calcium and magnesium, major food sources, and risk of type 2 diabetes in U.S. black women. Diabetes Care 29: 2238–2243.
  17. 17. Kirii K, Mizoue T, Iso H, Takahashi Y, Kato M, et al. (2009) Calcium, vitamin D and dairy intake in relation to type 2 diabetes risk in a Japanese cohort. Diabetologia 52: 2542–2550.
  18. 18. Pittas AG, Dawson-Hughes B, Li T, Van Dam RM, Willett WC, et al. (2006) Vitamin D and calcium intake in relation to type 2 diabetes in women. Diabetes Care 29: 650–656.
  19. 19. Sluijs I, Forouhi NG, Beulens JW, van der Schouw YT, Agnoli C, et al. (2012) The amount and type of dairy product intake and incident type 2 diabetes: results from the EPIC-InterAct Study. Am J Clin Nutr 96: 382–390.
  20. 20. Grantham NM, Magliano DJ, Hodge A, Jowett J, Meikle P, et al.. (2012) The association between dairy food intake and the incidence of diabetes in Australia: the Australian Diabetes Obesity and Lifestyle Study (AusDiab). Public Health Nutr: 1–7.
  21. 21. Elwood PC, Pickering JE, Fehily AM (2007) Milk and dairy consumption, diabetes and the metabolic syndrome: the Caerphilly prospective study. J Epidemiol Community Health 61: 695–698.
  22. 22. Villegas R, Gao YT, Dai Q, Yang G, Cai H, et al. (2009) Dietary calcium and magnesium intakes and the risk of type 2 diabetes: the Shanghai Women's Health Study. Am J Clin Nutr 89: 1059–1067.
  23. 23. Struijk EA, Heraclides A, Witte DR, Soedamah-Muthu SS, Geleijnse JM, et al.. (2012) Dairy product intake in relation to glucose regulation indices and risk of type 2 diabetes. Nutr Metab Cardiovasc Dis.
  24. 24. Tong X, Dong JY, Wu ZW, Li W, Qin LQ (2011) Dairy consumption and risk of type 2 diabetes mellitus: a meta-analysis of cohort studies. Eur J Clin Nutr 65: 1027–1031.
  25. 25. Wong WC, Cheung CS, Hart GJ (2008) Development of a quality assessment tool for systematic reviews of observational studies (QATSO) of HIV prevalence in men having sex with men and associated risk behaviours. Emerg Themes Epidemiol 5: 23.
  26. 26. Stroup DF, Berlin JA, Morton SC, Olkin I, Williamson GD, et al. (2000) Meta-analysis of observational studies in epidemiology: a proposal for reporting. Meta-analysis Of Observational Studies in Epidemiology (MOOSE) group. JAMA 283: 2008–2012.
  27. 27. von Elm E, Altman DG, Egger M, Pocock SJ, Gotzsche PC, et al. (2008) The Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) statement: guidelines for reporting observational studies. J Clin Epidemiol 61: 344–349.
  28. 28. Carter P, Gray LJ, Troughton J, Khunti K, Davies MJ (2010) Fruit and vegetable intake and incidence of type 2 diabetes mellitus: systematic review and meta-analysis. BMJ 341: c4229.
  29. 29. Zhang J, Yu KF (1998) What's the relative risk? A method of correcting the odds ratio in cohort studies of common outcomes. JAMA 280: 1690–1691.
  30. 30. Greenland S, Longnecker MP (1992) Methods for trend estimation from summarized dose-response data, with applications to meta-analysis. Am J Epidemiol 135: 1301–1309.
  31. 31. Orsini N BR, Greenland S (2006) Generalized least squares for trend estimation of summarized dose-response data. Stata Journal: 17.
  32. 32. Bodner-Montville J AJ, Ingwersen LA, Haggerty ES, Wilkinson Enns C, Perloff BP (2006) USDA Food and Nutrient Database for Dietary Studies: Released on the web. J Food Composition Analysis 19: 7.
  33. 33. Orsini N, Li R, Wolk A, Khudyakov P, Spiegelman D (2012) Meta-analysis for linear and nonlinear dose-response relations: examples, an evaluation of approximations, and software. Am J Epidemiol 175: 66–73.
  34. 34. Ericson U, Sonestedt E, Gullberg B, Hellstrand S, Hindy G, et al. (2013) High intakes of protein and processed meat associate with increased incidence of type 2 diabetes. Br J Nutr 109: 1143–1153.
  35. 35. Montonen J, Jarvinen R, Heliovaara M, Reunanen A, Aromaa A, et al. (2005) Food consumption and the incidence of type II diabetes mellitus. Eur J Clin Nutr 59: 441–448.
  36. 36. Azadbakht L, Mirmiran P, Esmaillzadeh A, Azizi T, Azizi F (2005) Beneficial effects of a Dietary Approaches to Stop Hypertension eating plan on features of the metabolic syndrome. Diabetes Care 28: 2823–2831.
  37. 37. Elwood PC, Pickering JE, Givens DI, Gallacher JE (2010) The consumption of milk and dairy foods and the incidence of vascular disease and diabetes: an overview of the evidence. Lipids 45: 925–939.
  38. 38. Thompson SG (1994) Why sources of heterogeneity in meta-analysis should be investigated. BMJ 309: 1351–1355.
  39. 39. Belin RJ, He K (2007) Magnesium physiology and pathogenic mechanisms that contribute to the development of the metabolic syndrome. Magnes Res 20: 107–129.
  40. 40. Kirii K, Iso H, Date C, Fukui M, Tamakoshi A (2010) Magnesium intake and risk of self-reported type 2 diabetes among Japanese. J Am Coll Nutr 29: 99–106.
  41. 41. Lopez-Ridaura R, Willett WC, Rimm EB, Liu S, Stampfer MJ, et al. (2004) Magnesium intake and risk of type 2 diabetes in men and women. Diabetes Care 27: 134–140.
  42. 42. Dong JY, Xun P, He K, Qin LQ (2011) Magnesium intake and risk of type 2 diabetes: meta-analysis of prospective cohort studies. Diabetes Care 34: 2116–2122.
  43. 43. Sousa GT, Lira FS, Rosa JC, de Oliveira EP, Oyama LM, et al. (2012) Dietary whey protein lessens several risk factors for metabolic diseases: a review. Lipids Health Dis 11: 67.
  44. 44. German JB, Gibson RA, Krauss RM, Nestel P, Lamarche B, et al. (2009) A reappraisal of the impact of dairy foods and milk fat on cardiovascular disease risk. Eur J Nutr 48: 191–203.
  45. 45. Kratz M, Baars T, Guyenet S (2012) The relationship between high-fat dairy consumption and obesity, cardiovascular, and metabolic disease. Eur J Nutr.
  46. 46. Riserus U, Arner P, Brismar K, Vessby B (2002) Treatment with dietary trans10cis12 conjugated linoleic acid causes isomer-specific insulin resistance in obese men with the metabolic syndrome. Diabetes Care 25: 1516–1521.
  47. 47. Niu K, Kobayashi Y, Guan L, Monma H, Guo H, et al.. (2012) Low-fat dairy, but not whole-/high-fat dairy, consumption is related with higher serum adiponectin levels in apparently healthy adults. Eur J Nutr.
  48. 48. Schurgers LJ, Vermeer C (2000) Determination of phylloquinone and menaquinones in food. Effect of food matrix on circulating vitamin K concentrations. Haemostasis 30: 298–307.
  49. 49. Beulens JW, van der AD, Grobbee DE, Sluijs I, Spijkerman AM, et al. (2010) Dietary phylloquinone and menaquinones intakes and risk of type 2 diabetes. Diabetes Care 33: 1699–1705.
  50. 50. Khan H, Kunutsor S, Franco OH, Chowdhury R (2013) Vitamin D, type 2 diabetes and other metabolic outcomes: a systematic review and meta-analysis of prospective studies. Proc Nutr Soc 72: 89–97.
  51. 51. Song Y, Wang L, Pittas AG, Del Gobbo LC, Zhang C, et al. (2013) Blood 25-hydroxy vitamin D levels and incident type 2 diabetes: a meta-analysis of prospective studies. Diabetes Care 36: 1422–1428.