Bone geometry and lower extremity bone stress injuries in male runners

Fractures in Thoroughbred racing and the potential for pre-race identification of horses at risk

Risk rates for and predisposing factors to fractures occurring in Thoroughbred racing that have been published in peer reviewed journals are documented. The potential for currently available techniques to identify horses at increased risk for fracture is discussed on the bases of principles, practicalities, advantages, disadvantages and current data. All are reviewed in light of justifiable decision making and importance of fractures to horseracing's social license.

A hierarchical clustering approach for examining the relationship between pelvis-proximal femur geometry and bone stress injury in runners

Bone stress injury (BSI) risk in runners is multifactorial and not well understood. Unsupervised machine learning approaches can potentially elucidate risk factors for BSI by identifying groups of similar runners within a population which differ in BSI incidence. Here, a hierarchical clustering approach is used to identify groups of collegiate cross country runners based on 2-dimensional frontal plane pelvis and proximal femur geometry, which was extracted from dual-energy X-ray absorptiometry scans and dimensionally reduced by principal component analysis. Seven distinct groups were identified using the cluster tree, with the initial split being highly related to female-male differences. Visual inspection revealed clear differences between groups in pelvis and proximal femur geometry, and groups were found to differ in lower body BSI incidence during the subsequent academic year (Rand index = 0.53; adjusted Rand index = 0.07). Linear models showed between-cluster differences in visually identified geometric measures. Geometric measures were aggregated into a pelvis shape factor based on trends with BSI incidence, and the resulting shape factor was significantly different between clusters (p < 0.001). Lower shape factor values, corresponding with lower pelvis height and ischial span, and greater iliac span and trochanteric span, appeared to be related to increased BSI incidence. This trend was dominated by the effect observed across clusters of male runners, indicating that geometric effects may be more relevant to BSI risk in males, or that other factors masked the relationship in females. More broadly, this work outlines a methodological approach for distilling complex geometric differences into simple metrics that relate to injury risk.

Distal Tibial Bone Properties and Bone Stress Injury Risk in Young Men Undergoing Arduous Physical Training

Trabecular microarchitecture contributes to bone strength, but its role in bone stress injury (BSI) risk in young healthy adults is unclear. Tibial volumetric BMD (vBMD), geometry, and microarchitecture, whole-body areal BMD, lean and fat mass, biochemical markers of bone metabolism, aerobic fitness, and muscle strength and power were measured in 201 British Army male infantry recruits (age 20.7 [4.3] years, BMI 24.0 ± 2.7 kg·m2) in week one of basic training. Tibial scans were performed at the ultra-distal site, 22.5 mm from the distal endplate of the non-dominant leg using High Resolution Peripheral Quantitative Computed Tomography (XtremeCT, Scanco Medical AG, Switzerland). Binary logistic regression analysis was performed to identify associations with lower body BSI confirmed by MRI. 20 recruits (10.0%) were diagnosed with a lower body BSI. Pre-injured participants had lower cortical area, stiffness and estimated failure load (p = 0.029, 0.012 and 0.011 respectively) but tibial vBMD, geometry, and microarchitecture were not associated with BSI incidence when controlling for age, total body mass, lean body mass, height, total 25(OH)D, 2.4-km run time, peak power output and maximum dynamic lift strength. Infantry Regiment (OR 9.3 [95%CI, 2.6, 33.4]) Parachute versus Line Infantry, (p ≤ 0.001) and 2.4-km best effort run time (1.06 [95%CI, 1.02, 1.10], p < 0.033) were significant predictors. Intrinsic risk factors, including ultradistal tibial density, geometry, and microarchitecture, were not associated with lower body BSI during arduous infantry training. The ninefold increased risk of BSI in the Parachute Regiment compared with Line Infantry suggests that injury propensity is primarily a function of training load and risk factors are population-specific.

Variability in tibia-fibular geometry is associated with increased tibial strain from running loads

Variation in tibial geometry may alter strain magnitude and distribution during locomotion. We investigated the effect of tibia-fibula geometric variations on tibial strain with running loads applied at various speeds. Participant-specific three-dimensional models of the tibia-fibula were created using lower limb computed tomography scans from 30 cadavers. Finite-element models were developed in FEBio, and running loads from 3, 4 and 5 m s-1 were applied to extract effective strain from the tibial shaft. Linear regression models evaluated the relationship between geometric characteristics and effective strain along the tibial shaft. We found a statistically significant positive relationship between: (i) increased thickness of the midshaft to upper tibia with increased condyle prominence and effective strain at points along the distal anterolateral and proximal posterior regions of the tibial shaft; and (ii) increased midshaft cortical thickness and effective strain at points along the medial aspect of the distal tibial shaft. It is possible that increased thickness in the more proximal region of the tibia causes strain to redistribute to areas that are more susceptible to the applied loads. A thickness imbalance between the upper and distal portions of the tibial shaft could have a negative impact on tibial stress injury risk.

Current risks factors and emerging biomarkers for bone stress injuries in military personnel

INTRODUCTION

Bone stress injuries (BSIs) have plagued the military for over 150 years; they afflict around 5 to 10% of military recruits, more so in women, and continue to place a medical and financial burden on defence. While the tibia generally adapts to the rigours of basic military training, the putative mechanisms for bone maladaptation are still unclear.

METHODS

This paper provides a review of the published literature on current risk factors and emerging biomarkers for BSIs in military personnel; the potential for biochemical markers of bone metabolism to monitor the response to military training; and, the association of novel biochemical 'exerkines' with bone health.

RESULTS

The primary risk factor for BSI in military (and athletic) populations is too much training, too soon. Appropriate physical preparation before training will likely be most protective, but routine biomarkers will not yet identify those at risk. Nutritional interventions will support a bone anabolic response to training, but exposure to stress, sleep loss, and medication is likely harmful to bone. Monitoring physiology using wearables-ovulation, sleep and stress-offer potential to inform prevention strategies.

CONCLUSIONS

The risk factors for BSIs are well described, but their aetiology is very complex particularly in the multi-stressor military environment. Our understanding of the skeletal responses to military training is improving as technology advances, and potential biomarkers are constantly emerging, but sophisticated and integrated approaches to prevention of BSI are warranted.

Geometric variation of the human tibia-fibula: a public dataset of tibia-fibula surface meshes and statistical shape model

Background

Variation in tibia geometry is a risk factor for tibial stress fractures. Geometric variability in bones is often quantified using statistical shape modelling. Statistical shape models (SSM) offer a method to assess three-dimensional variation of structures and identify the source of variation. Although SSM have been used widely to assess long bones, there is limited open-source datasets of this kind. Overall, the creation of SSM can be an expensive process, that requires advanced skills. A publicly available tibia shape model would be beneficial as it enables researchers to improve skills. Further, it could benefit health, sport and medicine with the potential to assess geometries suitable for medical equipment, and aid in clinical diagnosis. This study aimed to: (i) quantify tibial geometry using a SSM; and (ii) provide the SSM and associated code as an open-source dataset.

Methods

Lower limb computed tomography (CT) scans from the right tibia-fibula of 30 cadavers (male n = 20, female n = 10) were obtained from the New Mexico Decedent Image Database. Tibias were segmented and reconstructed into both cortical and trabecular sections. Fibulas were segmented as a singular surface. The segmented bones were used to develop three SSM of the: (i) tibia; (ii) tibia-fibula; and (iii) cortical-trabecular. Principal component analysis was applied to obtain the three SSM, with the principal components that explained 95% of geometric variation retained.

Results

Overall size was the main source of variation in all three models accounting for 90.31%, 84.24% and 85.06%. Other sources of geometric variation in the tibia surface models included overall and midshaft thickness; prominence and size of the condyle plateau, tibial tuberosity, and anterior crest; and axial torsion of the tibial shaft. Further variations in the tibia-fibula model included midshaft thickness of the fibula; fibula head position relative to the tibia; tibia and fibula anterior-posterior curvature; fibula posterior curvature; tibia plateau rotation; and interosseous width. The main sources of variation in the cortical-trabecular model other than general size included variation in the medulla cavity diameter; cortical thickness; anterior-posterior shaft curvature; and the volume of trabecular bone in the proximal and distal ends of the bone.

Conclusion

Variations that could increase the risk of tibial stress injury were observed, these included general tibial thickness, midshaft thickness, tibial length and medulla cavity diameter (indicative of cortical thickness). Further research is needed to better understand the effect of these tibial-fibula shape characteristics on tibial stress and injury risk. This SSM, the associated code, and three use examples for the SSM have been provided in an open-source dataset. The developed tibial surface models and statistical shape model will be made available for use at: https://simtk.org/projects/ssm_tibia.

Tibial-fibular geometry and density variations associated with elevated bone strain and sex disparities in young active adults

Tibial stress fracture is a common injury in runners and military personnel. Elevated bone strain is believed to be associated with the development of stress fractures and is influenced by bone geometry and density. The purpose of this study was to characterize tibial-fibular geometry and density variations in young active adults, and to quantify the influence of these variations on finite element-predicted bone strain. A statistical appearance model characterising tibial-fibular geometry and density was developed from computed tomography scans of 48 young physically active adults. The model was perturbed ±1 and 2 standard deviations along each of the first five principal components to create finite element models. Average male and female finite element models, controlled for scale, were also generated. Muscle and joint forces in running, calculated using inverse dynamics-based static optimization, were applied to the finite element models. The resulting 95th percentile pressure-modified von Mises strain (peak strain) and strained volume (volume of elements above 4000 με) were quantified. Geometry and density variations described by principal components resulted in up to 12.0% differences in peak strain and 95.4% differences in strained volume when compared to the average tibia-fibula model. The average female illustrated 5.5% and 41.3% larger peak strain and strained volume, respectively, when compared to the average male, suggesting that sexual dimorphism in bone geometry may indeed contribute to greater stress fracture risk in females. Our findings identified important features in subject-specific geometry and density associated with elevated bone strain that may have implications for stress fracture risk.

Tibial Bone Geometry Is Associated With Bone Stress Injury During Military Training in Men and Women

Bone stress injuries (BSI) are a common musculoskeletal condition among exercising and military populations and present a major burden to military readiness. The purpose of this investigation was to determine whether baseline measures of bone density, geometry, and strength, as assessed via peripheral quantitative computed tomography (pQCT), are predictive of tibial BSI during Marine Officer Candidates School training. Tibial pQCT scans were conducted prior to the start of physical training (n = 504; Male n = 382; Female n = 122) to measure volumetric bone mineral density (vBMD), geometry, robustness, and estimates of bone strength. Bone parameters were assessed at three tibial sites including the distal metaphysis (4% of tibial length measured from the distal endplate), mid-diaphysis (38% of tibial length measured from the distal endplate), and proximal diaphysis (66% of tibial length measured from the distal endplate). Injury surveillance data was collected throughout training. Four percent (n = 21) of the sample were diagnosed with a BSI at any anatomical site during training, 10 injuries were of the tibia. Baseline bone parameters were then tested for associations with the development of a tibial BSI during training and it was determined that cortical bone measures at diaphyseal (38 and 66%) sites were significant predictors of a prospective tibial BSI. At the mid-diaphysis (38% site), in a simple model and after adjusting for sex, age, and body size, total area [Odds Ratio (OR): 0.987, 0.983], endosteal circumference (OR: 0.853, 0.857), periosteal circumference (OR: 0.863, 0.824), and estimated bending strength (SSI; OR: 0.998, 0.997) were significant predictors of a BSI during training, respectively, such that lower values were associated with an increased likelihood of injury. Similarly, at the proximal diaphysis (66% site), total area (OR: 0.989, 0.985), endosteal circumference (OR: 0.855, 0.854), periosteal circumference (OR: 0.867, 0.823), robustness (OR: 0.007, 0.003), and SSI (OR: 0.998, 0.998) were also significant predictors of BSI in the simple and adjusted models, respectively, such that lower values were associated with an increased likelihood of injury. Results from this investigation support that narrower bones, with reduced circumference, lower total area, and lower estimated strength are associated with increased risk for tibial BSI during military training.

Biomechanical Basis of Predicting and Preventing Lower Limb Stress Fractures During Arduous Training

PURPOSE OF REVIEW

Stress fractures at weight-bearing sites, particularly the tibia, are common in military recruits and athletes. This review presents recent findings from human imaging and biomechanics studies aimed at predicting and preventing stress fractures.

RECENT FINDINGS

Peripheral quantitative computed tomography (pQCT) provides evidence that cortical bone geometry (tibial width and area) is associated with tibial stress fracture risk during weight-bearing exercise. The contribution of bone trabecular microarchitecture, cortical porosity, and bone material properties in the pathophysiology of stress fractures is less clear, but high-resolution pQCT and new techniques such as impact microindentation may improve our understanding of the role of microarchitecture and material properties in stress fracture prediction. Military studies demonstrate osteogenic outcomes from high impact, repetitive tibial loading during training. Kinetic and kinematic characteristics may influence stress fracture risk, but there is no evidence that interventions to modify biomechanics can reduce the incidence of stress fracture. Strategies to promote adaptive bone formation, in combination with improved techniques to assess bone strength, present exciting opportunities for future research to prevent stress fractures.

Training Load Capacity, Cumulative Risk, and Bone Stress Injuries: A Narrative Review of a Holistic Approach

Bone stress injuries (BSIs) are a common orthopedic injury with short-term, and potentially long-term, effects. Training load capacity, influenced by risk factors, plays a critical role in the occurrence of BSIs. Many factors determine how one's body responds to repetitive loads that have the potential to increase the risk of a BSI. As a scientific community, we have identified numerous isolated BSI risk factors. However, we have not adequately analyzed the integrative, holistic, and cumulative nature of the risk factors, which is essential to determine an individual's specific capacity. In this narrative review, we advocate for a personalized approach to monitor training load so that individuals can optimize their health and performance. We define “cumulative risk profile” as a subjective clinical determination of the number of risk factors with thoughtful consideration of their interaction and propose that athletes have their own cumulative risk profile that influences their capacity to withstand specific training loads. In our narrative review, we outline BSI risk factors, discuss the relationship between BSIs and training load, highlight the importance of individualizing training load, and emphasize the use of a holistic assessment as a training load guide.

Effects of Low Energy Availability on Bone Health in Endurance Athletes and High-Impact Exercise as A Potential Countermeasure: A Narrative Review

Endurance athletes expend large amounts of energy in prolonged high-intensity exercise and, due to the weight-sensitive nature of most endurance sports, often practice periods of dietary restriction. The Female Athlete Triad and Relative Energy Deficiency in Sport models consider endurance athletes at high-risk for suffering from low energy availability and associated health complications, including an increased chance of bone stress injury. Several studies have examined the effects of low energy availability on various parameters of bone structure and markers of bone (re)modelling; however, there are differences in findings and research methods and critical summaries are lacking. It is difficult for athletes to reduce energy expenditure or increase energy intake (to restore energy availability) in an environment where performance is a priority. Development of an alternative tool to help protect bone health would be beneficial. High-impact exercise can be highly osteogenic and energy efficient; however, at present, it is rarely utilized to promote bone health in endurance athletes. Therefore, with a view to reducing the prevalence of bone stress injury, the objectives of this review are to evaluate the effects of low energy availability on bone health in endurance athletes and explore whether a high-impact exercise intervention may help to prevent those effects from occurring.

High Cortico-Trabecular Transitional Zone Porosity and Reduced Trabecular Density in Men and Women with Stress Fractures

To determine whether stress fractures are associated with bone microstructural deterioration we quantified distal radial and the unfractured distal tibia using high resolution peripheral quantitative computed tomography in 26 cases with lower limb stress fractures (15 males, 11 females; mean age 37.1 ± 3.1 years) and 62 age-matched healthy controls (24 males, 38 females; mean age 35.0 ± 1.6 years). Relative to controls, in men, at the distal radius, cases had smaller cortical cross sectional area (CSA) (p = 0.012), higher porosity of the outer transitional zone (OTZ) (p = 0.006), inner transitional zone (ITZ) (p = 0.043) and the compact-appearing cortex (CC) (p = 0.023) while trabecular vBMD was lower (p = 0.002). At the distal tibia, cases also had a smaller cortical CSA (p = 0.008). Cortical porosity was not higher, but trabecular vBMD was lower (p = 0.001). Relative to controls, in women, cases had higher distal radial porosity of the OTZ (p = 0.028), ITZ (p = 0.030) not CC (p = 0.054). Trabecular vBMD was lower (p = 0.041). Distal tibial porosity was higher in the OTZ (p = 0.035), ITZ (p = 0.009), not CC. Stress fractures are associated with compromised cortical and trabecular microstructure.