Low energy availability (EA) is the failure of athletes to consume sufficient energy to cover the energy cost of exercise as well as energy required for optimal metabolic function and health (Loucks et al., 2011). Low EA has been reported in both female (Melin et al., 2015) and male (Viner et al., 2015) endurance athletes and is emerging as one of the most significant factors associated with athlete illness/injury (Tenforde et al., 2017). Original recognition of this syndrome in females, the Female Athlete Triad (Triad), identified the interrelatedness of low EA, menstrual dysfunction, and poor bone health (Joy et al., 2014; Nattiv et al., 1994). Indeed, there is ample evidence of the high prevalence (63%; Melin et al., 2015) and negative effects of low EA on hormonal function (Loucks & Thuma, 2003), bone health (Ihle & Loucks, 2004), and injury risk (Tenforde et al., 2017) in females. On the contrary, until recently, research on low EA in male athletes was lacking almost completely. Indeed, it was not until 2014, when the concept of relative energy deficiency in sport (RED-S) was introduced to include both sexes and a broader spectrum of health and performance-related concerns that the existence of low EA in males was acknowledged (Mountjoy et al., 2014). Since then, several reports have confirmed that male athletes do indeed suffer from the negative consequences of low EA (Tenforde et al., 2016), including low testosterone (TES; Gomez-Merino et al., 2002) and metabolic hormone (Koehler et al., 2016) levels, albeit likely at a lower threshold than females (20–25 kcal·kg−1 fat-free mass [FFM]·day−1; Fagerberg, 2017). Furthermore, it is noteworthy that despite significant reductions (10–40%) in TES levels due to low EA (Tenforde et al., 2016) have been observed, the values have remained within the normal clinical range. This is important as the use of nonathlete clinical cutoff values in assessing suppressed reproductive function in male athletes may lack sensitivity, at least based on current literature.
Despite growing awareness of the causes and outcomes of low EA, several issues related to its detection and practical management remain difficult. This includes the lack of a single sensitive tool to diagnose the phenomenon, particularly in field settings and with application to male athletes. We studied two cohorts of national to world-class distance female and male athletes with two primary purposes: (a) to provide a cross-sectional report on measurements of EA, metabolic and reproductive hormonal function, bone mineral density (BMD), injury/illness rates, and body composition during a high-volume/intensity precompetition training period and (b) to investigate the robustness of risk assessment of the Triad and RED-S diagnostic tools to identify athletes with these symptoms typically associated with low EA.
Methods
Subjects
We used word-of-mouth and online advertisements to recruit national to world-class female and male middle- and long-distance runners and race walkers from the Finnish Athletics Federation (FIN; 12 females and 10 males, at sea level) as well as Athletics Canada, Australia, and U.S. runners/race walkers at altitude in Flagstaff, AZ (FLAG; 27 females and 21 males). Inclusion criteria were as follows: 18–40 years, events from 800 m to 50 km walk, International Association of Athletics Federations (IAAF) score within the last 2 years of ≥880 (FIN: 14:37.17 and 17:36.39 min for 5,000 m for males and females, respectively) or ≥1,050 (FLAG: 13:45.20 min and 16:00.04 min for 5,000 m for males and females, respectively). Participants were informed about the study design before signing an informed consent. The studies were approved by the ethics committees of Helsinki University Hospital and University of Helsinki (BMD measurement) and University of Jyväskylä (other measurements). Study dropouts included 11 participants (FIN: n = 1, poor compliance [insufficient detail]; FLAG: n = 5 [poor compliance] and n = 5 [injuries during the recording week]).
Study Design
In a cross-sectional study design, we examined EA via 7-day energy intake (EI; via dietary recording) and exercise energy expenditure (via training recording), metabolic and reproductive hormonal function, BMD (via dual-energy X-ray absorptiometry), injury and illness rates, and body composition of identified endurance athletes during sea level precompetition training in FIN during spring 2015 (5% Olympians) and during a precompetition altitude training camp in FLAG during spring 2016 (27% Olympians).
Specifically, we divided female and male athletes into groups of apparently suboptimal and normal EA status based on two different criteria: (a) an acute objective measure of reproductive hormone function/menstrual status coupled with assessments using Triad and RED-S tools or (b) a single time-span calculation of EA from food (EI) and training (exercise energy expenditure) diaries over 7 days.
Blood Samples
Fasted blood samples were obtained from an antecubital vein for the analysis of insulin, TES (males), estradiol (E2; females), triiodothyronine (T3), and insulin-like growth factor 1 (IGF-1). FIN blood samples were collected into a 5-ml Vacuette gel serum tube (Greiner Bio-One GmbH, Kremsmünster, Austria) and centrifuged at 3,600 rpm for 10 min to collect serum; this was frozen at −20 °C, before the measurement of IGF-1, E2, TES, insulin by chemiluminometric immunological analysis, and T3 by electrochemiluminescence (IMMULITE 2000 XPi analyzer; Siemens Healthcare Diagnostics Inc., Tarrytown, NY). The assay sensitivities were 0.5 nmol/L (TES), 0.26 nmol/L (IGF-1), 2 mIU/L (insulin), 55 pmol/L (E2), and 0.6 pmol/L (T3). Reliability (coefficient of variation, CV) for between-day measurements was 9.1% (TES), 5.4% (IGF-1), 5.9% (insulin), 6.7% (E2), and 7.2% (T3). Eumenorrheic (EUM) females were measured on Days 3–7 during the follicular phase during a hard training period. Males and amenorrheic (AME) females were measured on an arbitrary day during a typical training week.
FLAG blood samples were collected into 8.5-ml serum separator gel tubes (BD Vacutainer, Franklin Lakes, NJ) and centrifuged at 3,400 rpm. Ferritin, insulin, T3, TES, and E2 were measured via electrochemiluminescence immunoassay method, and IGF-1 was measured via immunochemiluminometric assay. The samples were collected approximately at Day 3 of altitude exposure. The menstrual phase was not controlled but was taken into account during data analysis.
Body Composition and BMD
Body composition and BMD (whole body, femoral neck, and lumbar spine) were measured in an overnight fasted state via dual-energy X-ray absorptiometry (FIN: GE Lunar Prodigy Advance, Garnerville, NY; FLAG: GE Lunar DPX-IQ, Madison, WI), according to standardized protocols for body composition as recommended by Nana et al. (2015) and for BMD as recommended by The International Society for Clinical Densitometry (The International Society of Clinical Densitometry, 2017).
Food and Training Logs
Food and training diaries were recorded during the second week at FLAG and immediately after laboratory tests in FIN over a 7-day period. The principal investigator analyzed all dietary records (FIN: Nutri-Flow® software; Flow-Team, Oulu, Finland; FLAG: Food Processor software; ESHA, Salem, OR). Data collection protocols and methods of estimating EI (a combination of weighed and household measure dietary records) have been described in detail in Heikura et al. (2018), whereas the process of calculating EA following a field-based approach has been explained in Table 1.
Method for Estimating EA and Each of Its Components Based on Food and Training Records (Field-Based Approach)
EA equation | (EI − EEE)/FFM |
---|---|
EI | • Estimate total daily EI (kcal/day) via food records • Usually 4- to 7-day food records provide the best result (in terms of compliance and reflection of real-life dietary habits) • Prefer weighed food records over more vague estimations of portion sizes • Educate the athlete to minimize errors in recording • Analyze food records with a dietary analysis software to get an estimate of daily EI |
EEE and REE | REE • Estimate REE by using the Cunningham equation (Cunningham, 1991) • Divide total daily REE by 24 hr to yield hourly REE EEE • Estimate EEE via training diaries, where information on exercise mode, duration, and intensity (heart rate or ratings of perceived exertion) is included • Assign each training session or its parts a METs (Ainsworth et al., 2000) value that best corresponds to the type and intensity of that activity • Multiply the MET for each training session/part of the by the duration of the session or part of that session to yield MET hours • Multiply MET hours by the hourly REE to yield tEEE • Subtract REE that would have occurred regardless of exercise from tEEE so that only the additional energy cost of exercise is included in the EEE • Use this EEE value as part of the equation provided above |
FFM | FFM can be obtained via dual-energy X-ray absorptiometry or air displacement plethysmography (BodPod) provided “best practice” conditions for body composition analysis are utilized. Alternatively, well-standardized skinfold measures or bioelectrical impedance analysis may be used. |
Interpretation of EA values | The following cutoffs have been established (Loucks et al., 2011): Low EA: EA < 30 kcal·kg−1 FFM·day−1 Moderate EA: EA = 30–45 kcal·kg−1 FFM·day−1 Optimal EA: EA > 45 kcal·kg−1 FFM·day−1 |
Note. EA = energy availability; EI = energy intake; EEE = exercise energy expenditure; FFM = fat-free mass; REE = resting energy expenditure; METs = metabolic equivalents; tEEE = total EEE.
Questionnaires
Female participants completed the Low Energy Availability in Females Questionnaire (LEAF-Q; Melin et al., 2014), which was also used to assess self-reported amenorrhea. To examine the differences between subgroups, cumulative risk scores for Triad (cumulative risk assessment tool; Joy et al., 2014) and RED-S (Mountjoy et al., 2014) were calculated. The Triad tool was applied for males by replacing amenorrhea with low TES: values within the lowest quartile of clinical range (total TES: 9–38 nmol/L); males with low TES scored 1 point, whereas males with normal TES scored 0 point. The use of the lowest quartile for TES values was chosen because of clear evidence that although low EA is associated with significant reductions in TES, values have remained within the normal nonathlete clinical range (Tenforde et al., 2016); therefore, a subclinical threshold of low TES may be more appropriate in identifying low EA in elite endurance athletes. Athletes were grouped according to the Triad scores as low (0–1 points), moderate (2–5 points), or high (≥6 points) Triad risk score based on Joy et al. (2014). RED-S scoring system was applied to include the following conditions: low EA (<30 kcal·kg−1 FFM·day−1), amenorrhea (females) or low TES (males), BMD z score <−1.0, and a history of one or more stress fractures associated with low EA and/or hormonal dysfunction (see above). The presence of each condition was scored as 1 point, whereas the absence of any condition was scored as 0 point. The scores of all conditions were added together to form a RED-S score, and athletes were grouped as having low (0–1 points), moderate (2 points), or high (≥3 points) RED-S risk score. All participants answered an informal questionnaire on injury/illness history during the preceding 12 months (questions assessed the type of illness/injury and days of absence from running training due to each incident) and all-time history of bone-related injuries (number and type of bone-related injuries during the entire athletic career).
Statistical Analysis
Pooled-data statistical were conducted using SPSS Statistics 22 (IBM, Armonk, NY). Data were assessed via Shapiro–Wilk goodness-of-fit test, and nonparametric tests were used with nonnormally distributed data. Females and males were divided into subgroups based on EA (cutoff: 30 kcal·kg−1 FFM·day−1; Loucks et al., 2011) and reproductive function (amenorrhea: absence of ≥3 consecutive menses [Mountjoy et al., 2014] vs. eumenorrhea for females, and low TES vs. normal TES for males).
Differences and/or correlation analysis in body composition, dietary intakes, EA, BMD, injury/illness rates, and blood measures between sexes, subgroups, and cohorts were analyzed with either Student’s t test, one-way analysis of variance, or Pearson’s correlation coefficient (normally distributed data) or Mann–Whitney U test, Kruskal–Wallis test, or Spearman’s correlation coefficient (nonnormally distributed data). BMD was quantified via: Z > −1, normal BMD; Z < −1, a trend for low BMD; and Z < −2, clinically low BMD (Torstveit & Sundgot-Borgen, 2005). Analysis of covariance was applied to FLAG females for menstrual cycle phases, with E2 as a covariate, which was nonsignificant (p > .05). Therefore, we were confident that t tests were reflective of true differences.
Data are presented as means ± SD, with significance set to p < .05. Differences between and within sexes are expressed as Cohen’s effect sizes (Hopkins et al., 2009), with effect size threshold values of <0.2 (trivial), 0.2–0.5 (small), 0.5–0.8 (moderate), and >0.8 (large).
Results
Tables 2 and 3 summarize body composition, training/performance status, and dietary intake from the pooled dataset of 59 athletes, comparing subgroups of males and females based on EA and reproductive status, whereas these subgroups are examined in Tables 4 and 5 in terms of metabolic and health parameters.
Body Composition and Self-Reported Dietary and Training Data in Female and Male Athletes Categorized Into Low (AME and Low TES) and Normal (EUM and Normal TES) Reproductive Function Subgroups
Females | Males | |||||
---|---|---|---|---|---|---|
AME (n = 13) | EUM (n = 22) | ESa | Low TES (n = 10) | Normal TES (n = 14) | ESa | |
Age (years) | 23.8 ± 4.4 | 26.7 ± 2.9** | 0.81 | 27.4 ± 3.4 | 26.9 ± 4.5 | 0.12 |
Height (m) | 1.69 ± 0.04 | 1.68 ± 0.05 | 0.04 | 1.83 ± 0.06 | 1.81 ± 0.06 | 0.31 |
Weight (kg) | 52.9 ± 5.8 | 54.7 ± 5.0 | 0.33 | 70.5 ± 6.7 | 67.8 ± 6.2 | 0.42 |
Body fat (%) | 11.3 ± 4.1 | 12.1 ± 2.9 | 0.22 | 6.6 ± 1.1 | 7.2 ± 1.9 | 0.39 |
Energy (kJ·kg−1·day−1) | 208 ± 40 | 201 ± 36 | 0.17 | 212 ± 39 | 212 ± 22 | 0.01 |
CHO (g·kg−1·day−1) | 6.0 ± 1.3 | 6.0 ± 1.2 | 0.03 | 6.5 ± 2.0 | 6.6 ± 0.9 | 0.05 |
PRO (g·kg−1·day−1) | 2.4 ± 0.6 | 2.0 ± 0.4* | 0.71 | 2.2 ± 0.3 | 2.1 ± 0.3 | 0.39 |
Fat (g·kg−1·day−1) | 1.7 ± 0.5 | 1.8 ± 0.4 | 0.13 | 1.7 ± 0.5 | 1.8 ± 0.4 | 0.16 |
Fiber (g/day) | 54.1 ± 35.4 | 37.7 ± 9.7* | 0.73 | 47.2 ± 14.5 | 48.3 ± 17.9 | 0.07 |
IAAF score | 1,077 ± 45 | 1,102 ± 52 | 0.50 | 1,075 ± 83 | 1,071 ± 84 | 0.05 |
LEAF-Q score | 12.8 ± 4.8 | 8.3 ± 3.7** | 1.08 | N/A | N/A | N/A |
Running (km/week) | 101 ± 28 | 90 ± 30 | 0.40 | 117 ± 41 | 110 ± 24 | 0.22 |
Note. Values are presented as means ± SD. AME = amenorrhea; TES = testosterone; EUM = eumenorrhea; ES = effect size; CHO = carbohydrate; PRO = protein; IAAF = International Association of Athletics Federations; LEAF-Q = Low Energy Availability in Females Questionnaire (Melin et al., 2014); N/A = not available.
aMagnitude of within-sex difference.
*p < .05. **p < .01, significant within-sex difference.
Dietary and Training Data in Female and Male Athletes Categorized Into Low EA and Moderate EA
Females | Males | |||||
---|---|---|---|---|---|---|
Low EA (n = 11) | Moderate EA (n = 24) | ESa | Low EA (n = 6) | Moderate EA (n = 18) | ESa | |
Age (years) | 25.0 ± 4.5 | 25.9 ± 3.5 | 0.22 | 26.9 ± 3.8 | 27.2 ± 4.2 | 0.07 |
Height (m) | 1.69 ± 0.06 | 1.68 ± 0.04 | 0.04 | 1.83 ± 0.08 | 1.81 ± 0.05 | 0.32 |
Weight (kg) | 54.0 ± 6.8 | 54.1 ± 4.6 | 0.02 | 69.8 ± 6.9 | 68.6 ± 6.4 | 0.19 |
Body fat (%) | 11.8 ± 4.6 | 11.8 ± 2.8 | 0.01 | 6.5 ± 0.9 | 7.1 ± 1.7 | 0.44 |
Energy (kJ·kg−1·day−1) | 181 ± 27 | 214 ± 37* | 1.02 | 185 ± 35 | 222 ± 21** | 1.34 |
CHO (g·kg−1·day−1) | 5.4 ± 1.1 | 6.3 ± 1.2 | 0.75 | 5.6 ± 1.9 | 6.9 ± 1.1* | 0.86 |
PRO (g·kg−1·day−1) | 2.1 ± 0.6 | 2.2 ± 0.5 | 0.18 | 2.1 ± 0.2 | 2.1 ± 0.3 | 0.06 |
Fat (g·kg−1·day−1) | 1.4 ± 0.3 | 1.9 ± 0.4** | 1.14 | 1.4 ± 0.5 | 1.9 ± 0.4** | 1.28 |
Fiber (g/day) | 40.7 ± 11.1 | 45.2 ± 27.8 | 0.23 | 45.3 ± 12.4 | 48.7 ± 17.6 | 0.23 |
IAAF score | 1,113 ± 41 | 1,083 ± 52 | 0.65 | 1,057 ± 100 | 1,077 ± 77 | 0.23 |
LEAF-Q score | 11.1 ± 4.8 | 9.4 ± 4.6 | 0.35 | N/A | N/A | N/A |
Running (km/week) | 115 ± 23 | 84 ± 28** | 1.19 | 130 ± 43 | 107 ± 25 | 0.68 |
Note. Values are presented as means ± SD. EA = energy availability; ES = effect size; CHO = carbohydrate; PRO = protein; IAAF = International Association of Athletics Federations; LEAF-Q = Low Energy Availability in Females Questionnaire (Melin et al., 2014); N/A = not available.
aMagnitude of within-sex difference.
*p < .05. **p < .01, significant within-sex difference.
Metabolic and Reproductive Hormone Concentrations, Bone Density, and Injury and Illness History in Female and Male Athletes Categorized Into Low (AME and Low TES) and Normal (EUM and Normal TES) Reproductive Function Subgroups
Females | Males | |||||
---|---|---|---|---|---|---|
AME (n = 13) | EUM (n = 22) | ESa | Low TES (n = 10) | Normal TES (n = 14) | ESa | |
Reproductive hormones | ||||||
E2 (pmol/L) | 27 ± 19 | 111 ± 79** | 1.71 | N/A | N/A | N/A |
total TES (nmol/L) | N/A | N/A | N/A | 15.1 ± 3.0 | 25.0 ± 7.1** | 1.86 |
Metabolic hormones | ||||||
IGF-1 (nmol/L) | 225 ± 75 | 237 ± 67 | 0.17 | 197 ± 50 | 197 ± 52 | 0.00 |
T3 (pmol/L) | 2.9 ± 0.7 | 3.4 ± 0.6* | 0.77 | 3.5 ± 0.7 | 4.0 ± 0.6* | 0.77 |
insulin (pmol/L) | 3.9 ± 4.1 | 4.2 ± 1.5 | 0.11 | 2.7 ± 1.7 | 3.6 ± 1.5 | 0.56 |
Bone characteristics | ||||||
BMD (g·cm2) | 1.164 ± 0.041 | 1.222 ± 0.076* | 0.98 | 1.297 ± 0.083 | 1.282 ± 0.073 | 0.19 |
spine z score | −0.3 ± 0.9 | 0.3 ± 0.8* | 0.73 | 0.3 ± 0.7 | 0.3 ± 0.9 | 0.04 |
left femur z score | 1.0 ± 0.8 | 1.6 ± 1.1 | 0.68 | 1.2 ± 0.7 | 1.5 ± 1.0 | 0.39 |
right femur z score | 1.0 ± 0.8 | 1.5 ± 1.1 | 0.59 | 1.2 ± 0.7 | 1.6 ± 1.0 | 0.52 |
Injuries | ||||||
injuries and illnesses (days) | 58 ± 48 | 32 ± 35 | 0.63 | 16 ± 17 | 14 ± 10 | 0.15 |
bone injuries (days) | 32 ± 41 | 7 ± 18 | 0.85 | 9 ± 19 | 2 ± 8 | 0.52 |
EA (kcal·kg−1 FFM·day−1) | 32 ± 12 | 35 ± 9 | 0.29 | 31 ± 12 | 35 ± 5 | 0.47 |
Note. Values are presented as means ± SD. AME = amenorrhea; TES = testosterone; EUM = eumenorrhea; ES = effect size; E2 = estradiol; N/A = not available; IGF-1 = insulin-like growth factor 1; T3 = triiodothyronine; BMD = bone mineral density; z score = age-matched reference value for BMD; injuries and illnesses = the number of training days lost to injury and illness over the preceding 12 months prior to study; bone injuries = the number of training days lost to bone injuries during the preceding 12 months; EA = energy availability; FFM = fat-free mass.
aMagnitude of within-sex difference.
*p < .05. **p < .01, significant within-sex difference.
Metabolic and Reproductive Hormone Concentrations, Bone Density, and Injury and Illness History in Female and Male Athletes Categorized Into Low EA and Moderate EA
Females | Males | |||||
---|---|---|---|---|---|---|
Low EA (n = 11) | Moderate EA (n = 24) | ESa | Low EA (n = 6) | Moderate EA (n = 18) | ESa | |
Reproductive hormones | ||||||
E2 (pmol/L) | 70 ± 82 | 85 ± 73 | 0.19 | N/A | N/A | N/A |
total TES (nmol/L) | N/A | N/A | N/A | 14.8 ± 3.6 | 22.9 ± 8.0* | 1.40 |
Metabolic hormones | ||||||
IGF-1 (nmol/L) | 205 ± 51 | 245 ± 74 | 0.64 | 201 ± 55 | 196 ± 50 | 0.10 |
T3 (pmol/L) | 3.1 ± 0.6 | 3.3 ± 0.7 | 0.31 | 3.4 ± 0.8 | 3.9 ± 0.6 | 0.71 |
insulin (pmol/L) | 3.1 ± 2.0 | 4.6 ± 2.9 | 0.61 | 2.2 ± 1.8 | 3.6 ± 1.4 | 0.88 |
Bone characteristics | ||||||
BMD (g·cm2) | 1.196 ± 0.082 | 1.202 ± 0.067 | 0.09 | 1.316 ± 0.088 | 1.279 ± 0.072 | 0.46 |
spine z score | 0.1 ± 0.8 | 0.0 ± 0.9 | 0.13 | 0.1 ± 0.8 | 0.4 ± 0.8 | 0.42 |
left femur z score | 1.5 ± 1.0 | 1.3 ± 1.0 | 0.17 | 1.5 ± 0.9 | 1.3 ± 0.9 | 0.22 |
right femur z score | 1.4 ± 1.0 | 1.3 ± 1.0 | 0.17 | 1.6 ± 0.9 | 1.3 ± 0.9 | 0.38 |
Injuries | ||||||
injuries and illnesses (days) | 35 ± 46 | 45 ± 40 | 0.23 | 18 ± 22 | 14 ± 10 | 0.25 |
bone injuries (days) | 9 ± 20 | 19 ± 34 | 0.37 | 15 ± 24 | 2 ± 7 | 0.84 |
EA (kcal·kg−1 FFM·day−1) | 24 ± 6 | 38 ± 8** | 2.00 | 21 ± 6 | 37 ± 4** | 3.20 |
Note. Values are presented as means ± SD. EA = energy availability; ES = effect size; E2, estradiol; TES = testosterone; N/A = not available; IGF-1 = insulin-like growth factor 1; T3 = triiodothyronine; BMD = bone mineral density; z score = age-matched reference value for BMD; injuries and illnesses = the number of training days lost to injury and illness over the preceding 12 months prior to study; bone injuries = the number of training days lost to bone injuries during the preceding 12 months; FFM = fat-free mass.
aMagnitude of within-sex difference.
*p < .05. **p < .01, significant within-sex difference.
Reproductive/Metabolic Function, Bone Health, and Injury/Illness Rates
Overall, 37% of females were AME (FIN: 50% and FLAG: 30%). Forty percent of males (FIN: 55% and FLAG: 33%) were low TES, whereas none were clinically low (<9 nmol/L). FLAG females had higher spine, left, and right femur z scores (p = .024, d = 0.84; p = .022, d = 0.85; p = .041, d = 0.75, respectively); higher IAAF scores (p = .002, d = 1.15); and higher T3 (p < .001, d = 2.59) and insulin (p = .036, d = 0.81) concentrations than FIN females. Number of days missed training due to bone injuries was twofold greater in FIN versus FLAG females (25 ± 41 vs. 12 ± 23 days; p = .33, d = 0.39, respectively) and fivefold greater in FIN versus FLAG males (10 ± 20 vs. 2 ± 8 days; p = .28, d = 0.53, respectively). Similarly, FLAG males had higher values for T3 concentrations (p < .001, d = 2.33), EA (p = .047, d = 0.85), reported intakes of energy (p = .032, d = 0.93) and CHO (p < .001, d = 1.55), and IAAF scores (p < .001, d = 0.91) than FIN males. Conversely, FLAG males had lower whole-body BMD than FIN males (p < .001, d = 2.13). Compared with EUM females, AME females had lower E2, T3, and BMD and higher LEAF-Q scores, whereas low TES in males was associated with low T3 (Table 4). On the contrary, apart from TES values in males, there were no differences in hormone concentrations or BMD (Table 5) between low and moderate EA groups of either sex.
Figure 1 shows the prevalence of 0, 1, and ≥2 career stress fractures in AME and EUM females (1.46 ± 1.27 vs. 1.00 ± 1.38 stress fractures, respectively; p = .26, d = 0.41), and in low versus normal TES males (2.10 ± 2.10 vs. 0.50 ± 0.86 stress fractures, respectively; p = .042, d = 1.01). Effect size calculations showed that AME and low TES athletes had 4.5 times more training absences due to bone injuries compared with others (Table 4).
Correlations
In females, whole-body BMD correlated negatively with LEAF-Q score (r = −.410, p = .015). In males, correlations were found between all-time fracture history and TES (r = −.485, p = .016). In all athletes pooled, both Triad score (r = −.338, p = .009) and RED-S criteria (r = −.396, p = .002) correlated negatively with T3 and positively with all-time fracture history (r = .566, p < .001 and r = .328, p = .011, respectively) and negatively with spine z score (r = −.263, p = .044 and r = −.323, p = .013, respectively).
Triad and RED-S Risk Assessment
Athletes were scored based on available Triad criteria. Triad symptoms were absent in 34% of females, whereas 34%, 26%, and 6% presented with one, two, and three Triad conditions, respectively; this was even higher in AME females. AME females had a significantly higher cumulative risk score for the Triad compared with EUM females (5.2 ± 1.9 vs. 1.7 ± 1.4 points; p < .001, d = 2.10). The same was true for low versus moderate EA females (4.4 ± 2.3 vs. 2.4 ± 2.1 points; p = .027, d = 0.91). Of all females pooled, 31%, 51%, and 17% were classified into low-, moderate-, and high-risk categories, respectively. Low TES males had higher adapted Triad scores compared with normal TES males (2.0 ± 1.1 vs. 0.1 ± 0.4 points; p = .009, d = 1.43). Similarly, males with low EA had higher scores than males with moderate EA (2.8 ± 0.4 vs. 0.6 ± 0.8 points; p < .001, d = 3.00). Of males, 63% were classified into low and 37% into moderate-risk categories. Females had higher scores than males (3.0 ± 2.3 vs. 1.2 ± 1.2; p < .001, d = 1.03). Based on Triad risk categories, T3 levels in females (3.0 ± 0.5 vs. 3.5 ± 0.7 vs. 2.7 ± 0.4 pmol/L for low, moderate, and high risk groups, respectively; p = .021 between moderate and high) and in males (4.0 ± 0.6 vs. 3.4 ± 0.7 pmol/L for low and moderate risk groups, respectively; p = .032) were significantly lower, and the number of all-time fractures in females (0.3 ± 0.5 vs. 1.4 ± 1.4 vs. 2.2 ± 1.3 fractures for low, moderate, and high risk groups, respectively; p = .013 between low and high) and in males (0.3 ± 0.5 vs. 2.7 ± 1.9 fractures for low and moderate risk groups, respectively; p < .001) significantly higher in higher score groups. In addition, although not statistically significant, the number of training days lost due to bone injuries in the preceding year was markedly higher in both females (3 ± 11 vs. 18 ± 34 vs. 34 ± 38 days for low, moderate, and high risk groups, respectively; p = .13) and males (2 ± 8 vs. 10 ± 20 days for low and moderate risk groups, respectively; p = .23) with higher risk scores.
In addition, RED-S scoring system was used to score athletes based on symptoms of low EA. Differences in RED-S scores were significant between AME and EUM females (2.4 ± 0.8 vs. 0.5 ± 0.7; p < .001, d = 2.53) and between low versus moderate EA females (2.2 ± 0.9 vs. 0.7 ± 1.0; p = .001, d = 1.58). Likewise, differences between low and moderate TES males (2.0 ± 1.2 vs. 0.6 ± 0.8; p < .001, d = 2.62) and between low and moderate EA males (2.7 ± 0.8 vs. 0.3 ± 0.5; p < .001, d = 3.45) were significant. The division resulted in 43%, 43%, and 14% of females and 50%, 29%, and 21% of males classified into low, moderate, and high RED-S risk categories, respectively. Assessment of RED-S risk categories did not find significant differences in T3 levels or injury indices in females, whereas in males, T3 concentrations (4.0 ± 0.5 vs. 4.0 ± 0.8 vs. 3.1 ± 0.4 pmol/L for low, moderate, and high risk groups, respectively; p = .018 between low and high) were significantly lower, and the number of all-time fractures (0.6 ± 0.9 vs. 0.7 ± 1.3 vs. 3.2 ± 2.2 fractures for low, moderate, and high risk groups, respectively; p = .021 low vs. high, p = .039 moderate vs. high) significantly higher in higher score groups. In addition, although not statistically significant, the number of training days lost due to bone injuries in the preceding year was markedly higher in both females (7 ± 20 vs. 15 ± 27 vs. 48 ± 50 days for low, moderate, and high risk groups, respectively; p = .08) and males (3 ± 9 vs. 0 ± 0 vs. 18 ± 25 days for low and moderate risk groups, respectively; p = .09) with higher risk scores.
Discussion
This is the first study to examine EA, metabolic and reproductive hormonal function, BMD, injury and illness rates, and dietary characteristics in a large cohort (n = 59) of national to world-class (19% were Olympians) female and male distance athletes during a high-load precompetition training block. We divided athletes into subgroups based on reproductive function (menstrual status and TES levels) as well as a single time-point assessment of EA status (a cutoff point of 30 kcal·kg−1 FFM·day−1). Around 40% of females and males were AME or had low TES, respectively. These athletes had significantly lower sex hormone and T3 concentrations and 4.5 times greater incidence of bone injuries compared with EUM females and males with normal TES levels. We found that the Triad and RED-S tools were able to detect symptoms of impaired endocrine-metabolic function and bone health (all-time fractures) in both sexes with better accuracy than the assessment of EA based on food and training logs.
Seminal work by Loucks et al., within the standardized conditions of the laboratory, showed that low EA is the primary reason for the suppression of the metabolic and reproductive axis in females. These sophisticated studies established that below EA levels of ∼30 kcal·kg−1 FFM·day−1, the likelihood of suppression of female reproductive and metabolic hormones increases (Loucks & Callister, 1993; Loucks & Heath, 1994). Emerging data suggest that males, too, incur reciprocal effects (Tenforde et al., 2016), albeit at a lower EA threshold (between 20 and 25 kcal·kg−1 FFM·day−1; Fagerberg, 2017; Koehler et al., 2016). Such perturbations can have serious effects on bone (Mountjoy et al., 2014; Tenforde et al., 2015) and other aspects of physiological function and performance (Mountjoy et al., 2014).
Although the negative effects of low EA on hormonal function are robustly detected in controlled laboratory settings, our study clearly demonstrates that the real-world measurement of EA, via dietary and exercise recording, is challenging and lacking in sensitivity as diagnostic tool for the presence of low EA. Indeed, EA is assessed indirectly, based on self-reported records of food intake and estimates of energy expenditure. These techniques have inherent problems related to reliability and validity and may significantly skew EA calculations to both under/overestimate the true value during the short time period of assessment (Hills et al., 2014; Schoeller, 1995). For example, self-reporting resulted in an underestimation of EI by 5–21% (Stubbs et al., 2014), whereas the use of metabolic equivalents (METs) overestimated resting energy expenditure by 20% (Byrne et al., 2005) and wearable devices underestimated energy expenditure by ∼100–600 kcal/day (Murakami et al., 2016), albeit in nonathlete populations. Furthermore, there is no reference time frame over which a dietary/exercise assessment should be conducted, nor any standardized protocols by which any of the components should be measured. Although an EA assessment has appeal as apparently low-cost, low-technology, and noninvasive test, in reality, it carries a large burden in terms of athlete compliance and practitioner time. Cross-sectional estimates of low EA in athletes vary from 51% to 63% (Koehler et al., 2013; Melin et al., 2014, 2015). In our study, low EA (based on a cutoff of 30 kcal·kg−1 FFM·day−1) was found in 25% of male and 31% of female athletes but was poorly correlated with factors known to be affected by low EA, including reproductive, metabolic, and bone health (Table 5). Indeed, interestingly, the only difference between EA groups was a finding of lower TES levels in males with low EA. Despite the challenges associated with EA assessment, we are confident that the rigorous control applied to obtain high-quality dietary/training data allowed us to estimate EA to a sufficient detail. Our findings are in line with previous studies (Koehler et al., 2013; Melin et al., 2014) that failed to find an association between dietary EA and physiological markers of low EA in females. However, we also included males in our analysis and indeed showed a significant difference between males who had low versus moderate EA. While interesting, given no differences were found in other markers of EA between male EA groups, whether this finding is of relevance remains to be seen. Furthermore, as in our study, the EA assessment may occur at a different time point (and, in the case of FLAG, a different environment) to the scenarios in which poor match between energy intake and exercise energy expenditure may have caused health issues. Finally, other dietary factors often found in concert with low EA may have an independent or additive contribution to the outcomes (Reed et al., 2011). Indeed, in line with previous studies (Melin et al., 2016) who found significantly higher fiber intakes in oligo-AME versus EUM females, we noted a higher fiber intake reported by our AME females than their EUM counterparts.
Instead, qualitative screening tools (LEAF-Q, RED-S, and Triad risk assessment tools) and/or quantitative measurements of reproductive function (blood reproductive and metabolic hormones, BMD, and resting metabolic rate) may provide a better representation of an individual’s longer term EA, rather than just a snapshot, and may be a more sensitive way to diagnose low EA. Indeed, based on the current cross-sectional analysis, it appears that current reproductive function provides an equally good, but more objective marker of appropriate energy support for overall health: We found lower levels of T3, sex hormones, and BMD (females) in AME females and low TES males compared with others. Furthermore, when we divided athletes into low-, moderate-, and high-risk groups based on available Triad cumulative risk assessment tool (Joy et al., 2014) and a modified version of the RED-S consensus paper (Mountjoy et al., 2014), we found significantly lower T3 concentrations and up to ninefold higher number of all-time fractures in higher risk groups of both female and male athletes. For females, the Triad tool seemed to work better, whereas in males, both tools produced similar outcomes. In addition, although not significantly different, both tools showed up to 10 times higher frequency of lost training days due to bone injuries in the preceding year in the higher risk groups.
The prevalence of amenorrhea (37%) in the current study was in line with some (Barrack et al., 2014; Gibbs et al., 2014) but not all (Melin et al., 2014; Pollock et al., 2010) previous studies, whereas low BMD (17% in females and 0% in males) was significantly less than previously reported in females (Melin et al., 2014; Pollock et al., 2010). Interestingly, our cohort of elite athletes was of higher caliber than previously reported in the literature (27% Olympians in FLAG); thus, various differences in cohort rates of amenorrhea or BMD may not be surprising. While the use of menstrual status is probably one of the best markers of long-term EA, the use of oral contraceptives may mask the equation. In case of the use of oral contraceptives, perhaps other signs of chronic low EA, including low BMD, increased incidence of bone-related injuries, and low metabolic hormone concentrations, could be used to detect low EA in these athletes. Interestingly, low TES males had poorer health, despite being within the reference range. Indeed, the clinical cutoffs may not be applicable for elite athletes in assessing RED-S. The available literature lacks a more systematic analysis of the prevalence of these issues in male athletes.
Our findings of a significantly higher LEAF-Q score (12.8 vs. 8.3 points) in AME versus EUM females support the use of the LEAF-Q to assess increased risk of the Triad (Melin et al., 2014). We detected significant differences in the Triad and RED-S scores between the reproductive subgroups by categorizing athletes based on the Triad (Joy et al., 2014) and RED-S (Mountjoy et al., 2014) tools. This supports the work of Tenforde et al. (2017) who classified 239 AME/oligomenorrheic collegiate females into low- (71%), moderate- (25%), and high-risk (4%) categories based on the Triad tool and found that moderate- and high-risk category athletes were twice and nearly 4 times more likely to sustain a bone injury than low-risk athletes. We noted a higher prevalence of EA-associated problems in the lower caliber Finns. Whether this reflects cultural/genetic differences, the additional resources afforded to higher caliber athletes that might reduce the risk of problems, or the natural selection process by which the healthier athletes progress to higher levels, cannot be judged from our data but merits additional investigation. Future research should focus on finding a practical yet valid tool (set) for the assessment of low EA in female and male athletes to address this and other themes.
Conclusions
A relatively high proportion of national class elite female and male distance athletes showed signs of chronic energy conservation, including amenorrhea (37% of females) and low TES (40% of males). AME females were characterized by lower BMD, sex hormone, and T3 concentrations, and 4.5-fold higher rates of injury, compared with EUM females. Low TES males showed lower T3 concentrations, lower EA, and 4.5-fold higher rates of injury compared with males with normal TES. The assessment of risk scores based on Triad and RED-S criteria revealed lower T3 concentrations and higher number of all-time fractures in both female and male athletes with higher risk scores. Our findings suggest that assessing physiological symptoms of low EA (hormone concentrations, menstrual function, and bone density) or using LEAF-Q, Triad, and RED-S screening tools provides a better assessment of the overall health status of an individual than a snapshot of current EA.
Acknowledgments
The study was designed by I.A. Heikura, T. Stellingwerff, D. Bergland, A.A. Mero, A.L.T. Uusitalo, and L.M. Burke; data were collected and analyzed by I.A. Heikura, T. Stellingwerff, A.L.T. Uusitalo, D. Bergland, and L.M. Burke; data interpretation and manuscript preparation were undertaken by I.A. Heikura, T. Stellingwerff, A.L.T. Uusitalo, D. Bergland, and L.M. Burke. All authors approved the final version of the paper. This study was financially supported by grants from the Canadian Sport Institute, Australian Catholic University, Australian Institute of Sport, Finnish Sport Institute Foundation, and Finnish Sport Research Foundation. The authors wish to thank Risto Puurtinen (University of Jyväskylä), Salli Tommola (University of Jyväskylä), Hypo2 High Performance Sport Center, and Northern Arizona Orthopaedics for their role in the study, and Prof. Jennifer Peat for statistical support during data analysis. Finally, the authors wish to thank all the participants for volunteering for the study. The authors and funding agents do not have any conflicts of interests.
References
Ainsworth, B.E., Haskell, W.L., Whitt, M.C., Irwin, M.L., Swartz, A.M., Strath, S.J., … Leon, A.S. (2000). Compendium of physical activities: An update of activity codes and MET intensities. Medicine & Science in Sports & Exercise, 32(Suppl. 9), 498–504. PubMed ID: 10993420 doi:10.1097/00005768-200009001-00009
Barrack, M.T., Gibbs, J.C., De Souza, M.J., Williams, N.I., Nichols, J.F., Rauh, M.J., & Nattiv, A. (2014). Higher incidence of bone stress injuries with increasing female athlete triad-related risk factors: A prospective multisite study of exercising girls and women. The American Journal of Sports Medicine, 42(4), 949–958. PubMed ID: 24567250 doi:10.1177/0363546513520295
Byrne, N.M., Hills, A.P., Hunter, G.R., Weinsier, R.L., & Schutz, Y. (2005). Metabolic equivalent: One size does not fit all. Journal of Applied Physiology, 99(3), 1112–1119. PubMed ID: 15831804 doi:10.1152/japplphysiol.00023.2004
Cunningham, J.J. (1991). Body composition as a determinant of energy expenditure: A synthetic review and a proposed general prediction equation. The American Journal of Clinical Nutrition, 54(6), 963–969. PubMed ID: 1957828 doi:10.1093/ajcn/54.6.963
Fagerberg, P. (2017). Negative consequences of low energy availability in natural male bodybuilding: A review. International Journal of Sport Nutrition and Exercise Metabolism, 1–31. PubMed ID: 28530498 doi:10.1123/ijsnem.2016-0332
Gibbs, J.C., Nattiv, A., Barrack, M.T., Williams, N.I., Rauh, M.J., Nichols, J.F., & De Souza, M.J. (2014). Low bone density risk is higher in exercising women with multiple triad risk factors. Medicine & Science in Sports & Exercise, 46(1), 167–176. PubMed ID: 23783260 doi:10.1249/MSS.0b013e3182a03b8b
Gomez-Merino, D., Chennaoui, M., Drogou, C., Bonneau, D., & Guezennec, C.Y. (2002). Decrease in serum leptin after prolonged physical activity in men. Medicine & Science in Sports & Exercise, 34(10), 1594–1599. PubMed ID: 12370560 doi:10.1249/01.MSS.0000031097.37179.42
Heikura, I.A., Burke, L.M., Bergland, D., Uusitalo, A.L.T., Mero, A.A., & Stellingwerff, T. (2018). Impact of energy availability, health and sex on hemoglobin mass responses following LHTH altitude training in elite female and male distance athletes. International Journal of Sports Physiology and Performance, 12, 1–26. doi:10.1123/ijspp.2017-0547
Hills, A.P., Mokhtar, N., & Byrne, N.M. (2014). Assessment of physical activity and energy expenditure: An overview of objective measures. Frontiers in Nutrition, 1, 5. PubMed ID: 25988109 doi:10.3389/fnut.2014.00005
Hopkins, W.G., Marshall, S.W., Batterham, A.M., & Hanin, J. (2009). Progressive statistics for studies in sports medicine and exercise science. Medicine & Science in Sports & Exercise, 41(1), 3–13. PubMed ID: 19092709 doi:10.1249/MSS.0b013e31818cb278
Ihle, R., & Loucks, A.B. (2004). Dose-response relationships between energy availability and bone turnover in young exercising women. Journal of Bone and Mineral Research, 19(8), 1231–1240. PubMed ID: 15231009 doi:10.1359/JBMR.040410
The International Society of Clinical Densitometry. (2017, September 15). 2015 ISCD official positions—Adult. Retrieved from https://www.iscd.org/official-positions/2015-iscd-official-positions-adult/
Joy, E., De Souza, M.J., Nattiv, A., Misra, M., Williams, N.I., Mallinson, R.J., … Borgen, J.S. (2014). 2014 female athlete triad coalition consensus statement on treatment and return to play of the female athlete triad. Current Sports Medicine Reports, 13(4), 219–232. PubMed ID: 25014387 doi:10.1249/JSR.0000000000000077
Koehler, K., Achtzehn, S., Braun, H., Mester, J., & Schaenzer, W. (2013). Comparison of self-reported energy availability and metabolic hormones to assess adequacy of dietary energy intake in young elite athletes. Applied Physiology, Nutrition, and Metabolism, 38(7), 725–733. PubMed ID: 23980730 doi:10.1139/apnm-2012-0373
Koehler, K., Hoerner, N.R., Gibbs, J.C., Zinner, C., Braun, H., De Souza, M.J., & Schaenzer, W. (2016). Low energy availability in exercising men is associated with reduced leptin and insulin but not with changes in other metabolic hormones. Journal of Sports Sciences, 34(20), 1921–1929. PubMed ID: 26852783 doi:10.1080/02640414.2016.1142109
Loucks, A.B., & Callister, R. (1993). Induction and prevention of low-T3 syndrome in exercising women. The American Journal of Physiology, 264(5, Pt. 2), R924–R930. PubMed ID: 8498602 doi:10.1152/ajpregu.1993.264.5.R924
Loucks, A.B., & Heath, E.M. (1994). Dietary restriction reduces luteinizing hormone (LH) pulse frequency during waking hours and increases LH pulse amplitude during sleep in young menstruating women. The Journal of Clinical Endocrinology & Metabolism, 78(4), 910–915. PubMed ID: 8157720 doi:10.1210/jcem.78.4.8157720
Loucks, A.B., Kiens, B., & Wright, H.H. (2011). Energy availability in athletes. Journal of Sports Sciences, 29(Suppl. 1), S7–S15. PubMed ID: 21793767 doi:10.1080/02640414.2011.588958
Loucks, A.B., & Thuma, J.R. (2003). Luteinizing hormone pulsatility is disrupted at a threshold of energy availability in regularly menstruating women. The Journal of Clinical Endocrinology & Metabolism, 88(1), 297–311. PubMed ID: 12519869 doi:10.1210/jc.2002-020369
Melin, A., Tornberg, A.B., Skouby, S., Faber, J., Ritz, C., Sjodin, A., & Sundgot-Borgen, J. (2014). The LEAF questionnaire: A screening tool for the identification of female athletes at risk for the female athlete triad. British Journal of Sports Medicine, 48(7), 540–545. PubMed ID: 24563388 doi:10.1136/bjsports-2013-093240
Melin, A., Tornberg, A.B., Skouby, S., Moller, S.S., Faber, J., Sundgot-Borgen, J., & Sjodin, A. (2016). Low-energy density and high fiber intake are dietary concerns in female endurance athletes. Scandinavian Journal of Medicine & Science in Sports, 26(9), 1060–1071. PubMed ID: 26148242 doi:10.1111/sms.12516
Melin, A., Tornberg, A.B., Skouby, S., Moller, S.S., Sundgot-Borgen, J., Faber, J., … Sjodin, A. (2015). Energy availability and the female athlete triad in elite endurance athletes. Scandinavian Journal of Medicine & Science in Sports, 25(5), 610–622. PubMed ID: 24888644 doi:10.1111/sms.12261
Mountjoy, M., Sundgot-Borgen, J., Burke, L., Carter, S., Constantini, N., Lebrun, C., … Ljungqvist, A. (2014). The IOC consensus statement: Beyond the female athlete triad—Relative energy deficiency in sport (RED-S). British Journal of Sports Medicine, 48(7), 491–497. PubMed ID: 24620037 doi:10.1136/bjsports-2014-093502
Murakami, H., Kawakami, R., Nakae, S., Nakata, Y., Ishikawa-Takata, K., Tanaka, S., & Miyachi, M. (2016). Accuracy of wearable devices for estimating total energy expenditure: Comparison with metabolic chamber and doubly labeled water method. JAMA Internal Medicine, 176(5), 702–703. PubMed ID: 26999758 doi:10.1001/jamainternmed.2016.0152
Nana, A., Slater, G.J., Stewart, A.D., & Burke, L.M. (2015). Methodology review: Using dual-energy X-ray absorptiometry (DXA) for the assessment of body composition in athletes and active people. International Journal of Sport Nutrition and Exercise Metabolism, 25(2), 198–215. PubMed ID: 25029265 doi:10.1123/ijsnem.2013-0228
Nattiv, A., Agostini, R., Drinkwater, B., & Yeager, K.K. (1994). The female athlete triad. The inter-relatedness of disordered eating, amenorrhea, and osteoporosis. Clinics in Sports Medicine, 13(2), 405–418. PubMed ID: 8013041
Pollock, N., Grogan, C., Perry, M., Pedlar, C., Cooke, K., Morrissey, D., & Dimitriou, L. (2010). Bone-mineral density and other features of the female athlete triad in elite endurance runners: A longitudinal and cross-sectional observational study. International Journal of Sport Nutrition and Exercise Metabolism, 20(5), 418–426. PubMed ID: 20975110 doi:10.1123/ijsnem.20.5.418
Reed, J.L., Bowell, J.L., Hill, B.R., Williams, B.A., De Souza, M.J., & Williams, N.I. (2011). Exercising women with menstrual disturbances consume low energy dense foods and beverages. Applied Physiology, Nutrition, and Metabolism, 36(3), 382–394. PubMed ID: 21574870 doi:10.1139/h11-030
Schoeller, D.A. (1995). Limitations in the assessment of dietary energy intake by self-report. Metabolism: Clinical and Experimental, 44(2, Suppl. 2), 18–22. PubMed ID: 7869932 doi:10.1016/0026-0495(95)90204-X
Stubbs, R.J., O’Reilly, L.M., Whybrow, S., Fuller, Z., Johnstone, A.M., Livingstone, M.B., … Horgan, G.W. (2014). Measuring the difference between actual and reported food intakes in the context of energy balance under laboratory conditions. The British Journal of Nutrition, 111(11), 2032–2043. PubMed ID: 24635904 doi:10.1017/S0007114514000154
Tenforde, A.S., Barrack, M.T., Nattiv, A., & Fredericson, M. (2016). Parallels with the female athlete triad in male athletes. Sports Medicine, 46(2), 171–182. PubMed ID: 26497148 doi:10.1007/s40279-015-0411-y
Tenforde, A.S., Carlson, J.L., Chang, A., Sainani, K.L., Shultz, R., Kim, J.H., … Fredericson, M. (2017). Association of the female athlete triad risk assessment stratification to the development of bone stress injuries in collegiate athletes. The American Journal of Sports Medicine, 45(2), 302–310. PubMed ID: 28038316 doi:10.1177/0363546516676262
Tenforde, A.S., Fredericson, M., Sayres, L.C., Cutti, P., & Sainani, K.L. (2015). Identifying sex-specific risk factors for low bone mineral density in adolescent runners. The American Journal of Sports Medicine, 43(6), 1494–1504. doi:10.1177/0363546515572142
Torstveit, M.K., & Sundgot-Borgen, J. (2005). Low bone mineral density is two to three times more prevalent in non-athletic premenopausal women than in elite athletes: A comprehensive controlled study. British Journal of Sports Medicine, 39(5), 282–287; discussion 282–7. PubMed ID: 15849292 doi:10.1136/bjsm.2004.012781
Viner, R.T., Harris, M., Berning, J.R., & Meyer, N.L. (2015). Energy availability and dietary patterns of adult male and female competitive cyclists with lower than expected bone mineral density. International Journal of Sport Nutrition and Exercise Metabolism, 25(6), 594–602. PubMed ID: 26131616 doi:10.1123/ijsnem.2015-0073