Objective: To explore potential factors that could be associated with low bone mineral density (BMD) in female endurance runners.
Methods: Fifty two female endurance runners (1500 m to marathon), aged 18–44 years, took part. Body fat percentage, lumbar spine BMD, and femoral neck BMD were measured using the Hologic QDR 4500w bone densitometer. Data on training, menstrual cycle status, osteoporosis, and health related factors were obtained by questionnaire. Dietary variables were assessed from a prospective seven day dietary record of macronutrients and micronutrients.
Results: The mean (SD) lumbar spine and femoral neck BMD were 1.11 (0.11) and 0.89 (0.12) g/cm2 respectively. A backward elimination regression analysis showed that age, body mass, body fat, distance run, magnesium, and zinc intake were the variables significantly associated with BMD. Lumbar spine BMD (g/cm2) = −1.90 + (0.0486 × age (years)) + (0.342 × log mass (kg)) − (0.000861 × age2 (years)) − (0.00128 × distance (km/week)), with an R2 = 30.1% (SEE = 0.089 (95% confidence interval (CI) 0.05 to 0.23); p<0.001). Femoral neck BMD (g/cm2) = −2.51 − (0.00989 × age (years)) + (0.720 × log mass (kg)) + (0.000951 × magnesium (mg/day)) −(0.0289 × zinc (mg/day)) − (0.00821 × body fat (%)) − (0.00226 × distance (km/week)), with an R2 = 50.2% (SEE = 0.100 (95% CI 0.06 to 0.22); p<0.001). The negative association between skeletal BMD and distance run suggested that participants who ran longer distances had a lower BMD of the lumbar spine and femoral neck. Further, the results indicated a positive association between body mass and BMD, and a negative association between body fat and BMD.
Conclusions: The results suggest a negative association between endurance running distance and lumbar spine and femoral neck BMD, with a positive association between body mass and femoral neck and lumbar spine BMD. However, longitudinal studies are required to assess directly the effect of endurance running and body mass on BMD, and to see if the addition of alternative exercise that would increase lean body mass would have a positive effect on BMD and therefore help to prevent osteoporosis.
Statistics from Altmetric.com
Osteoporosis, which was once considered a disease of the elderly, is now prevalent in young women of any age whose bone mineral density (BMD) has fallen below a critical threshold.1 Various risk factors for osteoporosis have been postulated, such as menstrual cycle disturbances,2 diet,3 and inactivity.4 The research indicates that, although moderate exercise loads may benefit BMD,5 extreme loads may be detrimental to bone health in adults6 and children.7 Indeed, despite large volumes of weight-bearing exercise, female endurance runners have been shown to have lower BMD than their sedentary counterparts.8 However, recent studies in trained female runners have suggested that the periods of prolonged running alone do not affect the balance of bone turnover unless energy balance is simultaneously altered.9 This is apparent in female runners who have a tendency to seek a low body mass by chronic energy deficits,10 reducing luteinising hormone pulsatile activity resulting in a condition known as the female athlete triad.11 Thus there seems to be some conflict between the acquisition of bone by skeletal modelling and remodelling in response to the imposed exercise demands, and the bone mineral losses associated with energy deficits and amenorrhoea. Consequently, an improvement in the understanding of the key risk factors for low BMD in premenopausal, female endurance runners may help to reduce the risk of osteoporosis.12 Our aim was therefore to assess potential factors associated with low BMD at the femoral neck and lumbar spine in female endurance runners, using regression analysis.
Before initiation of the study, the purposes and procedures were explained verbally and in writing to each participant, highlighting the possible risks and benefits. Informed, written consent was obtained from all participants before any testing, in accordance with protocols defined by the Canterbury research ethics committee, who gave approval for the study.
To be accepted on the study, the participants had to fulfil the following criteria: female, white, aged 18–44 years (to decrease the possibility of women going through puberty and the menopause), and currently involved in endurance running. Participants were recruited from local and regional athletics clubs in the United Kingdom. They completed a questionnaire on osteoporosis risk factors and health before the start to evaluate health status. They were excluded from the study if they reported current or previous conditions that might have interfered with bone metabolism—for example, heart disease, long term corticosteroid use, smoking, or alcoholism. Women currently using or who had previously used oral contraceptives were not excluded.
Fifty two female endurance runners (recreational to elite standard) were invited and agreed to participate in the study. The participants exercised on average eight hours a week (2–24 hours), at recreational to elite level for an average of nine years (1–20 years). The mean (SD) distance run per week was 32 (18) km (5–70 km).
Participants completed questionnaires to assess training, dietary, health, and menstrual status (oligo/amenorrhoeic = 0–9 cycle a year; amenorrhoeic = 0–3 cycles a year; eumenorrhoeic ≥ 10 cycles a year). The questionnaires collected information on the history of the menstrual, training, and dietary status, as well as any recent changes (within the last year). A prospective, seven day dietary record was also completed, from which average energy (MJ/day), zinc (mg/day), magnesium (mg/day), phosphorus (mg/day), and calcium (mg/day) intake over this period was assessed (Dietmaster, version 4.0; Swift Computers Ltd, Phoenix, Arizona, USA).
Participants were weighed (wearing minimal clothing) to the nearest 0.1 kg (Balance Beam Scale; Seca, Germany), and height recorded to the nearest centimetre (stadiometer attached to Balance Beam Scale). Bone densitometry was used to assess body fat percentage and BMD. BMD was measured at the lumbar spine (L2–L4) and the femoral neck, and body fat was assessed during a whole body scan, using dual energy x ray absorptiometry (Hologic, QDR 4500w, Hologic Inc, Waltham, Massachusetts, USA). The standard Hologic protocol for positioning the lumbar spine and femoral neck was used. Bone mineral results were expressed as BMD (g/cm2). All scans were conducted and analysed in the osteoporosis unit at the Kent and Canterbury Hospital. The in vivo and in vitro precision of dual energy x ray absorptiometry for BMD in the laboratory was less than 1% for both the lumbar spine and femoral neck.
Statistical analysis was carried out using the Minitab statistical package (version 12.1). A backward, elimination log-linear regression model was used beginning with a saturated model incorporating all possible covariates and then sequentially eliminating the least important variables (based on the size of the t statistic), to assess the impact of all the variables on bone density at the femoral neck and lumbar spine sites. On the basis of proportional models proposed by Nevill and Holder,13 the following multiplicative model was proposed for BMD of the lumbar spine (L2–L4) and femoral neck:
The model can be linearised with a log transformation:and, provided that the error term ε1 has a normal distribution, linear regression can be used to estimate the parameters k, c, d, and e. Note that the proportional error ratio term in equation (1) is dimensionless—that is, a percentage error. Clearly, there is some operational advantage to incorporating age and age2 within an exponential term in equation (1), as well as providing a more plausible asymptotic decline for BMD with age. If a significant age2 parameter “e” was found, a peak in BMD had been identified, suggesting that BMD rises and peaks at age = −d/2e (using elementary differential calculus) and subsequently declines. The model can be developed to incorporate other factors of interest (within the exponential term) such as body fat, dietary macronutrients and micronutrients, and the distance walked a week (km).
To assess for the assumption of normality, an Anderson-Darling normality test and a Box-Cox transformation were used. Box and Cox14 suggested transforming the dependent variable (Y) using either Y‘ = (Yλ−1)/λ if λ = 1 or Y‘ = log(Y) if λ = 0, the value of the parameter λ being chosen to maximise the log-likelihood function. The transformation then provides residuals that best approximate to a normal distribution with constant variance. In the case of BMD of the lumbar spine (L2–L4) and the femoral neck, the optimal transformation parameters were λ = 0.09 and λ = 0.13 respectively, which were not significantly different from λ = 0, confirming that a log transformation of BMD was near optimal.
Pearson’s product-moment correlations were performed between dependent and independent variables to assess the relations.
Alpha was set at 0.05 for all statistical procedures.
Table 1 summarises the participant characteristics for all variables.
Of the 52 participants, the menstrual questionnaire classified 12 (23%) as oligo/amenorrhoeic (only one participant was amenorrhoeic), and 40 (77%) as eumenorrhoeic. Further, 46 (88%) had or were currently using various forms of oral contraceptives over a period of 59 (54) months. The daily energy intake for all participants over the seven day period was 7.5 (2) MJ/day, which placed them below the reference nutrient intake for some vitamins and minerals (table 1). There was no significant correlation (r = 0.04) between energy intake and distance run per week.
The mean (SEM) BMD for all participants at the lumbar spine and femoral neck were 1.11 (0.11) and 0.89 (0.12) g/cm2 respectively. The BMD values for the lumbar spine and femoral neck were on average −0.43 and + 0.45 SD from those of age matched participants in the databases. Of the cohort, two and one of the participants were classified as osteoporotic at the lumbar spine and femoral neck respectively. Osteopenia was diagnosed in nine and six of the participants at the lumbar spine and femoral neck respectively.
The backward elimination log-linear regression analysis showed that log lumbar spine BMD (g/cm2) was significantly associated with age (years) (p = 0.032), age2 (p = 0.021), log body mass (kg) (p = 0.007), and distance run (km/week) (p = 0.035). It was calculated that lumbar spine BMD increased up to the age of (−d)/(2e) = 0.0486/(2 × 0.000861) = 28.2 years, where it peaked and declined thereafter (parameters taken from the fitted regression model). The lumbar spine BMD was associated with body mass only, not body fat percentage, suggesting that it may be the lean component of body mass that is associated with lumbar spine BMD. Distance run was negatively associated with lumbar spine BMD, with participants who ran the longer distances appearing to have a lower BMD (lower by 1% in participants running a further 10 km/week) (fig 1).
Log transformed femoral neck BMD was significantly associated with distance (km/week) (p = 0.006), age (years) (p = 0.002), log body mass (kg) (p = 0.000), magnesium intake (mg/day) (p = 0.006), zinc intake (mg/day) (p = 0.009), and body fat percentage (p = 0.019). The regression model indicated that magnesium intake (mg/day) was positively related, whereas zinc intake (mg/day) was negatively related to femoral neck BMD. Distance run (km/week) was negatively associated with femoral neck BMD, with participants who ran the longer distances appearing to have a lower BMD (lower by 2% in participants running a further 10 km/week). A significant age related decline in femoral neck BMD was observed but the quadratic age2 parameter was not significant, suggesting that a peak in femoral neck BMD might have occurred in late adolescence (the teenage years) rather than early adulthood. Once again, body mass was positively associated with BMD, whereas body fat percentage was negatively associated with BMD, suggesting that it may be the lean component of body mass that was positively associated with BMD at this site (fig 2).
The resulting regression equations were:
R2 = 30.1%; SEE = 0.089 (95% confidence interval (CI) 0.05 to 0.23); p = 0.002
R2 = 50.2%; SEE = 0.100 (95% CI 0.06 to 0.22); p<0.001.
Note that by taking antilogs of the SEEs for lumbar spine and femoral neck BMD regression models, we obtained percentage errors of 9.3% (95% CI 5.1 to 25.8%) and 10.5% (95% CI 6.1 to 24.6%) respectively.
It has been thought that physical activity may help to offset the bone loss in female runners with inadequate diets,17 and it has been shown to exert a positive effect on BMD in postmenopausal women.15 As the frequency, duration, and intensity of the training regimens of endurance athletes far exceed the exercise of most postmenopausal women,15 one may expect the activity to exert a protective effect against bone loss in such athletes. However in this study, after controlling for differences in body size and age, distance run per week was negatively associated with lumbar spine and femoral neck BMD. The regression analysis indicated that participants who ran longer distances had a lower femoral neck and lumbar spine BMD (refer to the negative beta weights in the regression equations; p<0.05) than those who ran shorter distances. Indeed, the association suggested that participants who ran 10 km/week further than other participants had a 1% lower lumbar spine BMD and a 2% lower femoral neck BMD, regardless of the age of training onset. As there is evidence to suggest that low BMD is associated with increased fracture risk,18 such runners may be at an increased risk of fracture. However, these results should be interpreted with caution because the study was cross sectional and does not show a cause and effect relation. Further longitudinal studies are therefore required to support such claims.
Research has also suggested that long distance runners may have more dietary deficiencies,1,19 and that this may affect BMD more than the training per se. However, none of the dietary variables measured here were significantly associated with lumbar spine BMD. This agrees with previous findings,20 which suggests that only abnormal and restrictive eating behaviours are related to low BMD and future fracture risk. However, other studies9,10,19 have suggested that low energy intake can alter luteinising hormone pulsatile activity and decrease BMD levels. Such a mechanism may be site specific, for, although we did not find any direct relation between energy intake and BMD at any site, femoral neck BMD was significantly associated with magnesium and zinc intake, with higher magnesium and lower zinc intakes related to higher femoral neck BMD values. However, the low magnesium and zinc intakes recorded may reflect an undernourished state of the runners (7.5 (2) MJ/day while running 51 (28) km/week) and thus the effect of a low energy intake on BMD.3,19 Indeed, the dietary intakes were measured by a seven day prospective dietary intake record which relies heavily on the honesty and accuracy of the participants reporting, so the results need to be interpreted with some caution until further studies are completed. However, if the micronutrients are looked at individually, low levels of magnesium could downregulate parathyroid hormone secretion, decreasing calcium uptake and altering the bone remodelling process.21 However, the role of zinc in regulating BMD is less clear. Zinc was previously thought to exert a positive influence on BMD and to be essential for bone metabolism.21 Thus the negative influence of zinc on femoral neck BMD in this study was surprising. One possible clue may come from the way zinc affects the immune system by increasing osteoclastic activity.21
Another aspect to consider is the actual forces applied to the limbs during running, as previous research has suggested that they are crucial to BMD.22 Athletes involved in sports and training where forces applied to limbs are in excess of 10 times body weight (gymnastics, weightlifting, and volleyball) have been found to have higher BMD than those involved in sports where forces are only in the range of 5–10 times body weight, such as endurance running.23,24 It seems that, although running includes many more cycles of foot strike and thus force application than these other sports, the lower loads are not as stimulatory to bone accretion. This may explain our findings that long running distances appear to be associated with low BMD, because the longer distances covered would be associated with lower intensity exercise and thus lower forces being applied to the limbs. This supports the theory of Frost,22 in which a minimum force is required to elicit a bone remodelling response, and forces above this threshold would stimulate modelling at bone sites, resulting in a net bone gain. Therefore, magnitude of loading to the bone may be more important than the number of cycles, repetitions, or exercise mode.25,26
One would therefore expect body mass to have a positive effect on BMD. Indeed, this study suggests that body mass has a positive association with lumbar spine and femoral neck BMD. However, the regression analysis indicates that it is body mass, independent of body fat percentage, that is positively associated with BMD. Indeed, lean body mass has been found to correlate with appendicular27 and axial28 bone mass in previous studies. Further, muscle strength and lean mass have been reported to be positively associated with bone mass in a variety of weight bearing and non-weight bearing sports.23,29 Thus the force exerted by the muscle pull on the skeleton during muscular contraction may be a powerful osteogenic stimulus,30 although it has been suggested that, for the oestrogenic stimulus to have an effect on BMD, a regular menstrual cycle may need to be present1,2,31–32 However, we found no significant association between menstrual irregularity and BMD at the lumbar spine or femoral neck. Moreover, there was no significant association between oral contraceptive use and BMD at the lumbar spine or femoral neck from 1–192 months of use. This agrees with the findings of Prior et al33 that use of oral contraceptives did not affect bone density. However, other studies have reported positive34,35 and detrimental28,35 effects of oral contraceptive use on bone mass and fracture risk.36 These negative findings need further rigorous scientific studies because of a number of methodological deficiencies, such as not controlling for confounding variables in the study protocols—for example, smoking and endocrine status.
A major variable shown to influence BMD is age. BMD at the lumbar spine and femoral neck has been reported to peak between the third and fourth decades of life,37 the second and third decades of life,1 or as early as 18 years.38,39 The age at which this peak occurs and the factors affecting it are of utmost importance in reducing fracture risk.40,41 We observed that BMD was not significantly associated with age, with lumbar spine BMD peaking at 28 years, after which it declined with age. However, there was no evidence of an age related peak in femoral neck BMD, with a negative exponential decline apparent across the whole age range. This difference may be due to the lower amount of trabecular bone in the femoral neck (25%) compared with the lumbar spine (35%), promoting earlier accretion of bone mineral in the former.42,43
We observed a negative association between endurance running distance and BMD, and a positive association between body mass and BMD. This suggests that increasing body mass (independently of body fat percentage) may benefit BMD, and that long distance running may be associated with low BMD. As BMD is an important factor in the prediction of fracture risk, recommendations on exercise mode and the resulting effects on body composition may need to be considered. The association between distance run and BMD is a relatively exciting and novel finding which needs to be further investigated through longitudinal research.
Take home message
It may not be the exercise mode alone that affects BMD but the force applied to the limbs during such exercise and the resulting effects on body composition. Therefore planning of exercise regimens may need to consider the type and intensity of exercise undertaken and the subsequent impact on body composition.
If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.