牛乳を飲むほど骨折が増える

乳製品摂取量は大腿骨近位部骨折の発生リスクを高めるが、ヨーグルトとチーズは摂取量が多いほど骨折リスクが低減した。

（2023/11/22 管理者記載）

乳製品摂取量と大腿骨近位部骨折の発生リスクとの関連を調べた用量反応メタ解析の結果、牛乳摂取量の増加は骨折リスクの増大と関連するものの、ヨーグルトとチーズは摂取量が多いほど骨折リスクが低減したことを、米国・メリーランド大学のSuruchi Mishra氏らが報告した。
Journal of Nutritional Science誌2023年9月11日号に報告された。

　これまで、牛乳摂取は骨折の頻度を減少させて死亡リスクも低下させるという報告1）がある一方で、牛乳摂取量が多い人ほど骨折率や死亡率が高いという報告2）もあり、一貫性はない。そこで研究グループは、乳製品の摂取と大腿骨近位部骨折の発生リスクを評価するために用量反応メタ解析を実施した。

　PubMed（MEDLINE）とGoogle Scholarを用いて、1946～2021年12月に英語で発表された、乳製品摂取量と骨折リスクに関する前向きコホート研究を検索し、メタ回帰により用量反応相対リスクを導き出した。解析には13件の研究から成人48万6,950人、骨折1万5,320件を含めた。

　主な結果は以下のとおり。

★牛乳の摂取は、400g/日まで段階的に大腿骨近位部骨折のリスクを増大させ、用量効果として200g/日当たり7％のリスク増大と関連していた（相対リスク[RR]：1.07、95％信頼区間[CI]：1.05～1.10、p＜0.0001）。牛乳摂取量が0g/日群と比較して、400g/日群では大腿骨近位部骨折のリスクが最も高かった（RR：1.15、95％CI：1.09～1.21、p＜0.0001）。

★牛乳摂取量400g/日を超えると用量リスクは低減したが、それでも0g/日群と比較して750g/日まで骨折リスクは上昇した。

★観察された牛乳摂取量のどの範囲においても、0g/日群と比較して大腿骨近位部骨折のリスクが有意に低いという解析結果は得られなかった。

★ヨーグルト摂取に関する5研究の分析では、ヨーグルト250g/日当たり大腿骨近位部骨折リスクが15％低下するという逆相関が得られた（RR：0.85、95％CI：0.82～0.89）。

★同様に、チーズ摂取に関する5研究の分析では、チーズ43g/日（参考：6Pチーズ約2.5個分）当たり大腿骨近位部骨折リスクは19％低下した（RR：0.81、95％CI：0.72～0.92）。

★すべての乳製品の総摂取量と大腿骨近位部骨折の間に明らかな関連性は認められなかった（乳製品の総摂取量250g/日当たりのRR：0.97、95％CI：0.93～1.004、p＝0.079）。

　研究グループは、本研究の限界として「われわれの研究には小児のデータが不足していた。乳製品が小児集団に有益な影響を及ぼす可能性はあるが、小児における大腿骨近位部骨折を評価項目とした研究が不足していることから、乳製品と骨密度に関する文献を今後慎重に評価する必要がある」ことなどを挙げた。

Dairy intake and risk of hip fracture in prospective cohort studies: non-linear algorithmic dose-response analysis in 486 950 adults.

Suruchi Mishra, Kavitha Baruah, Vasanti S. Malik, and Eric L. Ding

Journal of nutritional science. 2023;12;e96. pii: e96.

https://www.ncbi.nlm.nih.gov/pmc/articles/PMC10495826/

Abstract

Previous studies on the relationship between dairy consumption and hip fracture risk have reported inconsistent findings. Therefore, we aimed to conduct an algorithmically driven non-linear dose-response meta-analysis of studies assessing dairy intake and risk of developing incident hip fracture. Meta-analysis from PubMed and Google Scholar searches for articles of prospective studies of dairy intake and risk of hip fracture, supplemented by additional detailed data provided by authors. Meta-regression derived dose-response relative risks, with comprehensive algorithm-driven dose assessment across the entire dairy consumption spectrum for non-linear associations. Review of studies published in English from 1946 through December 2021. A search yielded 13 studies, with 486 950 adults and 15 320 fractures. Non-linear dose models were found to be empirically superior to a linear explanation for the effects of milk. Milk consumption was associated with incrementally higher risk of hip fractures up to an intake of 400 g/d, with a 7 % higher risk of hip fracture per 200 g/d of milk (RR 1⋅07, 95 % CI 1⋅05, 1⋅10; P < 0⋅0001), peaking with 15 % higher risk (RR 1⋅15, 95 % CI 1⋅09, 1⋅21, P < 0⋅0001) at 400 g/d versus 0 g/d. Although there is a dose-risk attenuation above 400 g/d, milk consumption nevertheless continued to exhibit elevated risk of hip fracture, compared to zero intake, up to 750 g/d. Meanwhile, the analysis of five cohort studies of yoghurt intake per 250 g/d found a linear inverse association with fracture risk (RR 0⋅85, 95 % CI 0⋅82, 0⋅89), as did the five studies of cheese intake per 43 g/d (~1 serving/day) (RR 0⋅81, 95 % CI 0⋅72, 0⋅92); these studies did not control for socioeconomic status. However, no apparent association between total dairy intake and hip fracture (RR per 250 g/d of total dairy = 0⋅97, 95 % CI 0⋅93, 1⋅004; P = 0⋅079). There were both non-linear effects and overall elevated risk of hip fracture associated with greater milk intake, while lower risks of hip fracture were reported for higher yoghurt and cheese intakes.

Introduction

Hip fractures are common among elderly individuals and often have devastating sequelae. Dairy consumption has been promoted for fracture prevention because dairy products contain calcium, along with protein and phosphorus, and are sometimes fortified with vitamin D. However, prior systematic reviews have offered a conflicting and incomplete information for the beneficial effect of total dairy intake for hip fracture risk possibly due to the omission of certain studies and the failure to conduct careful dose-response analyses(1–5).

Dairy consumption may carry potential adverse health consequences. Dairy products are the leading source of saturated fat in the US diet(6), and dairy consumption has been associated with risk of cancer of the prostate and other organs(7–11), raising questions about the overall health effects of recommending this product group for fracture prevention.

We therefore undertook an evidence-based review and non-linear dose-response algorithmic meta-analysis of prospective studies to examine associations between dairy intake and hip fracture risk. Notably, we implemented a unique 1000+ model data-driven non-linear analysis algorithm to sequentially compare thousands of knots to determine the optimum non-linear dose-effect association for various dairy products, including milk, yoghurt, cheese and total dairy.

Methods

Search strategy

We searched the published scientific literature for prospective cohort studies of human participants reporting dairy or milk intake and the end point of incident hip fracture. Literature review was conducted in two phases, prior to July 2016, and after July 2016 until the end of 2021. In the first phase, two investigators identified articles published in English from 1946 to July 2016 through a search of MEDLINE (PubMed) and Google Scholar. We also identified additional articles from the reference lists of the extracted articles and reviews, and articles that cited extracted articles. The following terms were used for the MEDLINE search: (dairy[All Fields] OR milk [All Fields] OR human milk[All Fields] OR milk [MeSH Terms] OR yoghurt[All Fields] OR cheese[All Fields]) AND humans [MeSH Terms] AND (bone [All Fields] OR bones [All Fields] OR bone and bones [All Fields] OR bone [All Fields] OR fractures, bone [MeSH Terms] OR fractures [All Fields] OR bone fractures [All Fields] OR fracture [All Fields] OR accidental falls [MeSH Terms] OR dxa[All Fields] OR dexa[All Fields] OR (osteoporosis, postmenopausal[MeSH Terms] OR osteoporosis [All Fields]) OR bone diseases, metabolic[MeSH Terms] OR (bone[All Fields] OR osteopenia[All Fields]) OR (osteomalacia[MeSH Terms] OR osteomalacia[All Fields]) OR (bone density[MeSH Terms] OR bone[All Fields] OR bone density[All Fields] OR (bone[All Fields] AND mineral[All Fields] AND density[All Fields]) OR bone mineral density[All Fields])). Similar search terms were used for Google Scholar. No language restriction was applied.

In the second search phase, spanning July 2016 through December 2021, we conducted three parallel searches in effort to be more comprehensive, using two sets of Google Scholar searches, which limits results to <1000 results per search, and using MEDLINE to find additional prospective studies on dairy and risk of hip fractures. The following terms were used for the first Google Scholar search: (Dairy OR milk OR yoghurt OR cheese) AND (Bone OR fracture OR hip fractures OR dxa OR dexa OR osteoporosis OR bone diseases OR osteopenia OR osteomalacia OR bone mineral density). For the second wider Google Scholar search, we used the broader terms: (Dairy OR milk OR yoghurt OR cheese) AND (hip fracture OR fracture). We also conducted a search of MEDLINE between July 2016 through December 2021 using the following search terms: (dairy[All Fields] OR milk [All Fields] OR human milk[All Fields] OR milk [MeSH Terms] OR yoghurt[All Fields] OR cheese[All Fields]) AND humans [MeSH Terms] AND (bone [All Fields] OR bones [All Fields] OR bone and bones [All Fields] OR bone [All Fields] OR fractures, bone [MeSH Terms] OR fractures [All Fields] OR bone fractures [All Fields] OR fracture [All Fields] OR accidental falls [MeSH Terms] OR dxa[All Fields] OR dexa[All Fields] OR (osteoporosis, postmenopausal[MeSH Terms] OR osteoporosis [All Fields]) OR bone diseases, metabolic[MeSH Terms] OR (bone[All Fields] OR osteopenia[All Fields]) OR (osteomalacia[MeSH Terms] OR osteomalacia[All Fields]) OR (bone density[MeSH Terms] OR bone[All Fields] OR bone density[All Fields] OR (bone[All Fields] AND mineral[All Fields] AND density[All Fields]) OR bone mineral density[All Fields])). All search results were downloaded, cross-matched for de-duplication processing before assessment for eligibility within the phase of the search.

Study selection

We first conducted an initial screening of all titles. Abstracts were reviewed if the titles appeared to meet the inclusion criteria. Then, we assessed all potentially relevant studies based on text reviews. We included prospective studies of dairy intake and risk of hip fracture in humans and excluded literature reviews, clinical trials and studies with an ecological, case–control or cross-sectional design. We excluded papers using the criteria listed below. Papers were excluded based on the first applicable exclusion category: animal studies, no dairy consumption, no fracture, review, not prospective, correction, duplicate population or identified in prior search (if second search phase). And for those nearly qualified studies that report any/total fracture but not hip fracture, we emailed authors for additional information, but conditionally excluded such studies unless additional data is provided. We also excluded studies of individuals with pre-existing bone diseases and studies that did not report relative risks (RRs) or hazard ratios (HRs), dairy dose data and the corresponding 95 % confidence intervals (CIs).

Data extraction

Two authors (S.M. and K.B.) independently selected studies and performed the data extraction. To resolve any discrepancies, a third investigator (E.L.D.) was consulted, and any disagreements were settled by consensus of all three authors. The primary exposure variables were milk, yoghurt, cheese, cream, soured milk and yoghurt, total dairy without cream and total dairy with cream consumption. Total dairy with cream included milk, yoghurt, frozen yoghurt and cheese. The outcome of interest was risk of hip fracture. Data extraction was conducted using a standardised data collection form. The following characteristics of the identified articles were recorded: first author, year of publication, cohort name, study population, study design, country, mean/median age, adjustment for body mass index (BMI), adjustment for energy intake (kilocalories), industry funding, percentage of women, follow-up duration, calcium and vitamin D supplement use, and baseline dietary calcium and vitamin D levels, where available, as well as the quality measures noted below.

When data were incomplete or missing, we contacted the authors to obtain additional data. We also attempted to contact authors who we believed were likely to have data on dairy intake and fracture risk based on their prior publications, and obtained additional updated data from the Nurses’ Health Study, Health Professionals Follow-up Study, Swedish Mammography Cohort and the Framingham Offspring Study through personal communication with study authors or research staff. Other authors of cohort studies with potentially relevant data were contacted, but none provided dairy and hip fracture risk data that could allow for their inclusion.

Quality measures

Study reports were qualitatively examined for means of dietary assessment, diagnosis and fracture assessment, sample size and statistical measures.

Statistical analysis

To determine whether results indicated non-linear or linear relationships between dairy product consumption and fracture risk, the linearity of associations was first examined using spline analyses and dose-response meta-regressions. For non-linear analysis, we comprehensively assessed the dose across the entire dairy consumption spectrum for non-linear associations, for both single-knot and dual-knot splines for comparative model analysis, at individual 25 g intervals, and used data-driven algorithms to select for the optimal single spline. If a single-knot spline was found superior to linear coefficient model, then upon identification of the optimal single spline, double-spline dose models were tested straddling around the single spline to further assess for potential model superiority. Splined variables were created with the use of MKSPLINE in STATA version 13.1 (College Station, TX). Goodness-of-fit tests and chi-square statistics were used to determine the most appropriate knot points and maximal goodness of fit to assess the fit of potential non-linear relationships and to determine the best dose-response inflection point of any potential non-linear association across a range of intakes. The Greenland and Longnecker dose-response generalised least-square meta-regression analysis using GLST command in STATA was used to derive the dose-response RRs, P-values and 95 % CIs(12). For non-linear associations, the shapes of the associations within individual studies were assessed via repeated modelling of different knot points in order to determine potential splines for best model fit, indicating the optimal non-linear spline model. The pooled non-linear dose curves were visualised using Spaghetti Plots, as previously applied(13). Where non-linearity was not indicated, forest plots were produced to portray linear dose-response relationships and corresponding 95 % CIs via random-effects models. Begg's funnel plots and Begg's test are used to assess potential publication bias. To assess statistical heterogeneity, Cochran's Q test was conducted and the I2 statistic was calculated, representing the percentage of total variation attributable to between-study heterogeneity(14). Potential publication bias was assessed with Egger's test. Two-sided P-values <0⋅05 were considered to be statistically significant.

Results

The Phase 1 PubMed and the Google Scholar searches prior to July 2016 yielded 8781 citation titles for review that yielded 11 relevant studies, two of which were identified from the search of references of the retrieved publications and subsequent author contact. The Nurses’ Health Study, Health Professionals Follow-up Study, Swedish Mammography Cohort and Framingham Offspring Study were among the studies that also provided their updated data through personal communication to the senior author. The Phase 2 set of searches of PubMed and Google Scholar from July 2016 to December 2021 yielded 2218 citation titles for review, of which 2 were new relevant and eligible articles, not previously included or captured in personal communications. Altogether, across the two phases of searches, thirteen articles were included in this analysis. The flow of these articles is described in Fig. 1.

Characteristics of included studies

We based our analyses on 13 studies from 16 unique study populations involving 486 950 participants and 15 320 fractures. Additional data from the Swedish Mammography cohort(15) (Michaelsson, personal communication, November 2017) added information on 1568 additional hip fractures beyond the original published study data set(16), for an updated total of 5827 hip fractures from that cohort.

Although all studies assessed dietary intake by food frequency questionnaires, repeat measurements of dietary intake were available only in the Nurses’ Health Study, Health Professionals Follow-up Study and Swedish Mammography Cohorts. All other cohorts collected dietary data only at baseline. Diagnosis of fracture was based on self-report confirmed by medical records, death certificates, radiographic and operative report or events validated through linkage with national patient registries.

Five studies provided data that were not separated by gender(17–22). Two cohorts included only men(16,23). Four cohorts included only women(15,23–25), and two provided data for men and women separately(26,27). One study provided data pooled by gender and also separately for women(28). Seven studies were conducted in the USA(19–25), two in Sweden(15,16), two in Norway(26,27), one in Denmark(22) and one each in Japan(26), France(18) and Europe(17). Of the thirteen studies included in this meta-analysis, three reported receiving industry funding(18,19,21), and one study with partial industry funding for a doctoral student(22). Table 1 shows the baseline characteristics for all included studies.

Total dairy consumption and hip fracture

Among the eight cohorts that investigated associations between total dairy intake and hip fracture risk, six found no association (Health Professionals Follow-up Study(23), Adventist Health Study-2(20), Three-City Study(18), Framingham Offspring Study(19), EPIC Elderly Network on Aging and Health(17), NHANES-1 Epidemiologic Follow-up study(24)). Two reported decreased risk of hip fracture with increased dairy intake (Feskanich, personal communication, October 2016 and Michaelsson, personal communication, November 2017), albeit only with moderate dairy intake in one of these (Michaelsson, personal communication, November 2017).

Although there were eight study populations for which data on total dairy intake were available (Nurses’ Health Study(23), Health Professionals Follow-up Study(23), Swedish Mammography Cohort(15), Adventist Health Study-2(20), The Three-City Study(18), Framingham Offspring Cohort(19), EPIC Elderly Network on Aging and Health(17), NHANES-1 Epidemiologic Follow-up study(24)) (Feskanich, personal communication, October 2016 and Michaelsson, personal communication, November 2017), we first conducted an assessment of potential non-linearity of the association among the eight studies with multiple dose-response data points per study (19 total data points). These included the (Nurses’ Health Study(23), Health Professionals Follow-up Study(23), Swedish Mammography Cohort(15), Adventist Health Study-2(20), Framingham Offspring Cohort(19), EPIC Elderly Network on Aging and Health(17)) (Feskanich, personal communication, October 2016 and Michaelsson, personal communication, November 2017).

Across the range of doses, no spline knot point yielded significant evidence of a non-linear association, with 100 g/d of total dairy being the closest potential spline point, which was not significant for non-linearity (P = 0⋅19). Thus, we analysed all available studies of total dairy and hip fracture in linear dose response fashion. Pooling of the studies in a linear dose increment, expressed per 250 g/d, yielded no significant association between total dairy consumption and hip fracture risk (RR 0⋅97, 95 % CI 0⋅93, 1⋅00). There was no significant heterogeneity, with I2 = 34⋅5 % and overlapping null (95 % CI 0, 71 % and P = 0⋅15). The forest plot is shown in Fig. 2. No publication bias trend was observed (Egger's Test P = 0⋅66).

Milk intake and hip fracture

For milk intake, one study (Michaelsson, personal communication, November 2018)(15) showed elevated risk of hip fracture with increased milk intake, seven showed no association(16,18,19,21,22,25,28), and three studies (Nurses’ Health Study(23), Health Professionals Follow-up Study(23) and middle-aged Norwegian cohort(26)) (Feskanich, personal communication, October 2016) showed lower risk with increased milk intake.

We conducted an assessment of potential non-linearity of the association among the studies with multiple dose-response data points per study (fifty-nine data points total). Across the ranges of different doses algorithmically tested, a strong spline knot point yielded significant strongest evidence of a non-linear association, with 400 g/d being the most robust dose-risk inflection point (P for non-linearity < 0⋅0001); double-knot spline models straddling 400 g/d were not superior to a single knot at 400 g/d. Thus, we analysed the fourteen populations with milk and hip fracture data in a non-linear dose response fashion. Up to a daily dose of 400 g/d, the results showed that milk intake was associated with increased fracture risk, with a dose effect of 7 % increased risk per 200 g/d of milk (RR 1⋅07, 95 % CI 1⋅05, 1⋅10; P < 0⋅0001). At 400 g/d of milk intake, compared to 0 g/d, there was a 15 % higher risk of hip fracture (RR 1⋅15, 95 % CI 1⋅09, 1⋅21, P < 0⋅0001). However, above 400 g/d of milk, although there is a slight attenuation in the higher risk with a dose trend effect of a −4⋅5 % per 200 g/d of milk (RR 0⋅955, 95 % CI 0⋅92, 0⋅99, P = 0⋅007) compared to the risk peak of 400 g/d – milk consumption nevertheless continued to exhibit elevated risk of hip fracture with intakes up to 750 g/d. In no range of observed milk intake did the analysis show any significantly lower risk of hip fracture compared to 0 g/d. The test of heterogeneity yielded a between-study variance of 0 in all models, indicating no significant heterogeneity. Results are shown in the Spaghetti Plot of Fig. 3.

Yoghurt, fermented milk, cheese and hip fracture

A smaller number of studies examined yoghurt, fermented milk or cheese intake in relation to hip fractures, and none of these studies controlled for socioeconomic status (SES). Because of the limited number of studies, we did not perform a non-linear analysis, which requires more datapoints for reliable computation; and none were apparent via tested knots. Thus, pooling the five studies(15,16,18,19,21) (Michaelsson personal communication, November 2017) of yoghurt and fermented milk intake with hip fracture, the linear analysis per 250 g/d yielded an inverse association of 15 % lower hip fracture risk (RR 0⋅85, 95 % CI 0⋅82, 0⋅89). The model did not reveal significant heterogeneity, with between-study variance of 0, P = 0⋅60 (Fig. 4). No publication bias trend was observed (Egger's Test P = 0⋅16). Similarly, pooling the five studies of cheese intake and hip fracture in linear dose analysis found that cheese intake per 43 g/d (~1 serving/day) was associated with an inverse association of 19 % lower hip fracture risk (RR 0⋅81, 95 % CI 0⋅72, 0⋅92), while one study without precise cheese dose (any versus no cheese) found no association for cheese and hip fracture(22). Heterogeneity was not significant, with I2 = 52 (95 % CI 0–82, P = 0⋅08) (Fig. 5). No publication bias trend was observed (Egger's Test P = 0⋅37).

Discussion

In the present analysis, there was no significant relationship between overall dairy intake and hip fracture risk. These findings, derived from the most complete data set and most comprehensive and algorithmically driven meta-analysis of prospective cohorts to date, contradict the findings of earlier meta-analyses for milk(1,2,4–6).

Among dairy subfractions, the largest body of data was for milk intake, which was directly and significantly associated with risk of hip fracture. This risk was moderated at intakes above 400 g/d, but at no level of milk intake was there a suggestion of a significant protective association between milk and fracture risk.

There were noteworthy differences in findings between cohorts in the US and Swedish Mammography Cohort. However, even with sensitivity exclusion of the large Swedish Mammography cohort, milk still did not exhibit any inverse association with fractures. Moreover, the cohort of Swedish men and two cohorts in Norway also were in agreement with the overall findings on milk and fracture. Nevertheless, it is important to consider potential factors that could explain the observed variability. Although all studies included in this review used similar designs and diagnostic criteria, there were differences in age, body weight and vitamin D exposure, which may have influenced their findings.

Age

The three cohorts that reported protective effects of milk intake on hip fracture risk (middle-aged Norwegian cohorts, Nurses’ Health Study and Health Professionals Follow-up Study) had participants with mean baseline ages ranging between 47 and 58 years, whereas mean baseline age of participants in all other cohorts ranged between 53⋅5 and 77 years.

Body weight

In the Nurses’ Health Study and Health Professionals Follow-up Study, the inverse associations between dairy consumption and hip fracture risk were mainly limited to those with larger BMI(23). Obesity is much more common in the US than in other countries(31), affecting approximately 40 % among US adults >20 years old, compared with approximately 19 % of adults in Sweden(31,32).

Vitamin D

Vitamin D enhances calcium absorption(33,34). It is noteworthy that the US-based Nurses’ Health Study, which reported baseline intakes of 9⋅7 μg vitamin D and 1060 mg calcium per day, found a protective effect of milk intake, while the Swedish Mammography Cohort, which reported baseline intakes of 4⋅6 μg vitamin D and 916⋅5 mg calcium per day, found increased hip fracture risk with increased milk intake.

In the US, all milk is fortified with vitamin D and is the single biggest contributor to vitamin D intakes (58 % in men, 39 % in women)(35). In Sweden, only low-fat milk are fortified(36) (although they remain significant contributors to vitamin D intake)(37). In Norway, milk is not routinely fortified with vitamin D(38). Milk's role in vitamin D status is not entirely straightforward, however. Despite vitamin D fortification, dairy calcium can suppress vitamin D activation, depressing circulating levels of 1,25-dihydroxyvitamin D(8).

These observations suggest that discordant findings among studies may be in part attributable to differences in age, participant weight and vitamin D fortification. Although differences in fat content of milk may seem to be an additional possible explanatory factor for between-study differences, analyses of fracture risk by fat content did not yield different effects (Michaelsson, personal communication, November 2017)(15,16).

The mechanisms behind potentially higher risk of hip fracture with milk intake are not entirely clear. d-galactose, a breakdown product of lactose, has been linked to oxidative stress and inflammation, potentially increasing the risk of fracture and mortality(15,16). Higher levels of d-galactose in milk, compared with lower or non-existent levels in cheese, may partly explain positive associations of milk and inverse association of cheese with hip fracture risk reported in the Swedish Mammography Cohort(15). Furthermore, milk intake during adolescence modestly increases height, a risk factor for fracture(39), which a study of teenage milk intake yielded higher hip fracture risk in men, and no effect in women in later adulthood(39). Similarly, a Canadian study also found no association between childhood milk consumption and fracture risk(40).

Fewer studies examined relationships between yoghurt and cheese intake and hip fracture risk, but overall these studies indicated general inverse relationships for yoghurt and hip fracture, and possibly for cheese, though unclear. A seemingly protective association with cheese may be surprising, particularly given that cheese is high in sodium, a micronutrient strongly associated with loss of bone density(41). The interpretation of this finding is complicated by the fact that the included studies did not adjust for SES except those by Michaelson, Holvik et al. and Feart et al. In the US and Europe, yoghurt and cheese consumption is positively associated with SES, while SES appears to have a neutral or negative association with milk consumption(42–45). Indeed, the one Danish prospective study(22) that did adjust for education level did not find any inverse association for milk and hip fracture risk.

Strengths

Our study has notable strengths. It includes a larger number of prospective studies that were available for previous meta-analyses of dairy-fracture risk(1,23,27,46–48), and we also included de novo updated estimates from personal communications with Swedish collaborators (Michaelson, personal communication, November 2017). Furthermore, our study included a detailed algorithmically driven dose-response meta-regression incorporating non-linear splines and accounting for differing reference intake levels, which together allowed us to algorithmically identify significant non-linear associations for milk intake that had not been previously identified. Furthermore, our results of null associations for total dairy with hip fracture risk are corroborated by a more recent study(49).

Limitations

This study also has limitations. Data on dairy consumption were collected using food frequency questionnaires, which may have minor variation in the interpretation of one glass of milk. Data on incident hip fractures were obtained through self-reports in some of the studies included. A large proportion of people are discharged to a nursing facilities after suffering from a hip fracture and mortality after hip fracture is high, which may have contributed to lower reports of fracture events and loss to follow-up. Although SES is sometimes associated with cheese and yoghurt intake, however, SES was explored in the Nurses’ Health Study and HPFS, which have little SES confounding due to the relatively homogeneous populations of health professionals. That said, there may be other non-linear associations that were not detected due to limited available studies and data points. While the Swedish Mammography Cohort was an influential milk study, it was the most comprehensive study of its kind, and its prominence and large size (and simultaneous publication with the Cohort of Swedish Men in the same BMJ paper) indicates that it was unlikely to have be suffered file-drawer publication bias that smaller studies tend to suffer, but which large studies don't. Plus, an analysis without the Swedish Mortality Cohort (SMC) still exhibited non-linearity, albeit smaller slopes. Furthermore, the SMC study's exclusion did not materially change the results for total dairy, yoghurt or cheese, with risk of hip fracture. Moreover, additional follow-up data from the SMC with ~1500 more hip fracture cases only further strengthened the findings. Finally, our study lacked data on dairy intake and hip fracture risk in children; it is possible dairy may have beneficial effects on paediatric populations, but this requires a careful future evaluation of the dairy-bone density literature, given the lack of studies with hip fracture end points in children.

Conclusions

Overall, increasing milk intake is associated with significantly increased risk of hip fracture, while intake and risk have inverse relationships for yoghurt and cheese intake. Because of the divergence in dairy products, total dairy intake is not associated importantly with risk of hip fracture.