Diffusion tensor imaging of hamstring muscles after acute strain injury and throughout recovery in collegiate athletes

OBJECTIVE

To identify the region of interest (ROI) to represent injury and observe between-limb diffusion tensor imaging (DTI) microstructural differences in muscle following hamstring strain injury.

MATERIALS AND METHODS

Participants who sustained a hamstring strain injury prospectively underwent 3T-MRI of bilateral thighs using T1, T2, and diffusion-weighted imaging at time of injury (TOI), return to sport (RTS), and 12 weeks after RTS (12wks). ROIs were using the hyperintense region on a T2-weighted sequence: edema, focused edema, and primary muscle injured excluding edema (no edema). Linear mixed-effects models were used to compare diffusion parameters between ROIs and timepoints and limbs and timepoints.

RESULTS

Twenty-four participants (29 injuries) were included. A significant ROI-by-timepoint interaction was detected for all diffusivity measures. The edema and focused edema ROIs demonstrated increased diffusion at TOI compared to RTS for all diffusivity measures (p-values < 0.006), except λ1 (p-values = 0.058-0.12), and compared to 12wks (p-values < 0.02). In the no edema ROI, differences in diffusivity measures were not observed (p-values > 0.82). At TOI, no edema ROI diffusivity measures were lower than the edema ROI (p-values < 0.001) but not at RTS or 12wks (p-values > 0.69). A significant limb-by-timepoint interaction was detected for all diffusivity measures with increased diffusion in the involved limb at TOI (p-values < 0.001) but not at RTS or 12wks (p-values > 0.42). Significant differences in fractional anisotropy over time or between limbs were not detected.

CONCLUSION

Hyperintensity on T2-weighted imaging used to define the injured region holds promise in describing muscle microstructure following hamstring strain injury by demonstrating between-limb differences at TOI but not at follow-up timepoints.

MRI radiomics for hamstring strain injury identification and return to sport classification: a pilot study

OBJECTIVE

To determine if MRI-based radiomics from hamstring muscles are related to injury and if the features could be used to perform a time to return to sport (RTS) classification. We hypothesize that radiomics from hamstring muscles, especially T2-weighted and diffusion tensor imaging-based features, are related to injury and can be used for RTS classification.

SUBJECTS AND METHODS

MRI data from 32 athletes at the University of Wisconsin-Madison that sustained a hamstring strain injury were collected. Diffusion tensor imaging and T1- and T2-weighted images were processed, and diffusion maps were calculated. Radiomics features were extracted from the four hamstring muscles in each limb and for each MRI modality, individually. Feature selection was performed and multiple support vector classifiers were cross-validated to differentiate between involved and uninvolved limbs and perform binary (≤ or > 25 days) and multiclass (< 14 vs. 14-42 vs. > 42 days) classification of RTS.

RESULT

The combination of radiomics features from all diffusion tensor imaging and T2-weighted images provided the most accurate differentiation between involved and uninvolved limbs (AUC ≈ 0.84 ± 0.16). For the binary RTS classification, the combination of all extracted radiomics offered the most accurate classification (AUC ≈ 0.95 ± 0.15). While for the multiclass RTS classification, the combination of features from all the diffusion tensor imaging maps provided the most accurate classification (weighted one vs. rest AUC ≈ 0.81 ± 0.16).

CONCLUSION

This pilot study demonstrated that radiomics features from hamstring muscles are related to injury and have the potential to predict RTS.

Ultrasound shear wave seeds reduced following hamstring strain injury but not after returning to sport

OBJECTIVES

The purpose of the study was to investigate differences in ultrasound shear wave speed (SWS) between uninjured and injured limbs following hamstring strain injury (HSI) at time of injury (TOI), return to sport (RTS), and 12 weeks after RTS (12wks).

METHODS

This observational, prospective, cross-sectional design included male and female collegiate athletes who sustained an HSI. SWS imaging was performed at TOI, RTS, and 12wks with magnetic resonance imaging. SWS maps were acquired by a musculoskeletal-trained sonographer at the injury location of the injured limb and location-matched on the contralateral limb. The average SWS from three 5 mm diameter Q-boxes on each limb were used for analysis. A linear mixed effects model was performed to determine differences in SWS between limbs across the study time points.

RESULTS

SWS was lower in the injured limb compared to the contralateral limb at TOI (uninjured - injured limb difference: 0.23 [0.05, 0.41] m/s, p = 0.006). No between-limb differences in SWS were observed at RTS (0.15 [-0.05, 0.36] m/s, p = 0.23) or 12wks (-0.11 [-0.41, 0.18] m/s, p = 0.84).

CONCLUSIONS

The SWS in the injured limb of collegiate athletes after HSI was lower compared to the uninjured limb at TOI but not at RTS or 12 weeks after RTS.

CRITICAL RELEVANCE STATEMENT

Hamstring strain injury with structural disruption can be detected by lower injured limb shear wave speed compared to the uninjured limb. Lack of between-limb differences at return to sport may demonstrate changes consistent with healing. Shear wave speed may complement traditional ultrasound or MRI for monitoring muscle injury.

KEY POINTS

• Ultrasound shear wave speed can non-invasively measure tissue elasticity in muscle injury locations. • Injured limb time of injury shear wave speeds were lower versus uninjured limb but not thereafter. • Null return to sport shear wave speed differences may correspond to structural changes associated with healing. • Shear wave speed may provide quantitative measures for monitoring muscle elasticity during recovery.

Association of quantitative diffusion tensor imaging measures with time to return to sport and reinjury incidence following acute hamstring strain injury

Hamstring strain injuries (HSI) are a common occurrence in athletics and complicated by limited prognostic indicators and high rates of reinjury. Assessment of injury characteristics at the time of injury (TOI) may be used to manage athlete expectations for time to return to sport (RTS) and mitigate reinjury risk. Magnetic resonance imaging (MRI) is routinely used in soft tissue injury management, but its prognostic value for HSI is widely debated. Recent advancements in musculoskeletal MRI, such as diffusion tensor imaging (DTI), have allowed for quantitative measures of muscle microstructure assessment. The purpose of this study was to determine the association of TOI MRI-based measures, including the British Athletic Muscle Injury Classification (BAMIC) system, edema volume, and DTI metrics, with time to RTS and reinjury incidence. Negative binomial regressions and generalized estimating equations were used to determine relationships between imaging measures and time to RTS and reinjury, respectively. Twenty-six index injuries were observed, with five recorded reinjuries. A significant association was not detected between BAMIC score and edema volume at TOI with days to RTS (p-values ≥ 0.15) or reinjury (p-values ≥ 0.13). Similarly, a significant association between DTI metrics and days to RTS was not detected (p-values ≥ 0.11). Although diffusivity metrics are expected to increase following injury, decreased values were observed in those who reinjured (mean diffusivity, p = 0.016; radial diffusivity, p = 0.02; principal effective diffusivity eigenvalues, p-values = 0.007-0.057). Additional work to further understand the directional relationship observed between DTI metrics and reinjury status and the influence of external factors is warranted.

Vertical Loading Rate Is Not Associated with Running Injury, Regardless of Calculation Method

INTRODUCTION

Loading rate (LR), the slope of the vertical ground reaction force (vGRF), is commonly used to assess running-related injury risk. However, the relationship between LR and running-related injuries, including bone stress injuries (BSI), is unclear. Inconsistent findings may result from the numerous LR calculation methods that exist and their application across different running speeds.

PURPOSE

This study aimed to assess the influence of calculation method and running speed on LR values and to determine the association of LR during healthy running with subsequent injury.

METHODS

Healthy preseason running data and subsequent injury records from Division I cross-country athletes ( n = 79) over four seasons (2015-2019) at 2.68 m·s -1 , preferred training pace, and 4.47 m·s -1 were collected. LR at each speed was calculated four ways: 1) maximum and 2) average slope from 20% to 80% of vGRF magnitude at impact peak (IP), 3) average slope from initial contact to IP, and 4) average slope from 3% to 12% of stance time. Linear mixed effects models and generalized estimation equations were used to assess LR associations.

RESULTS

LR values differed depending on speed and calculation method ( P value <0.001). The maximum slope from 20% to 80% of the vGRF at 4.47 m·s -1 produced the highest LR estimate and the average slope from initial contact to IP at 2.68 m·s -1 produced the lowest. Sixty-four injuries (20 BSI) were observed. No significant association was found between LR and all injuries or BSI across any calculation method ( P values ≥0.13).

CONCLUSIONS

Calculation method and running speed result in significantly different LR values. Regardless of calculation method, no association between LR and subsequent injury was identified. Thus, healthy baseline LR may not be useful to prospectively assess running-related injury risk.

Preseason Eccentric Strength Is Not Associated with Hamstring Strain Injury: A Prospective Study in Collegiate Athletes

INTRODUCTION

Established risk factors for hamstring strain injuries (HSI) include older age and prior HSI. However, these are nonmodifiable and have a limited role in injury prevention. Eccentric hamstring strength is a common component of HSI prevention programs, but its association with injury is less clear.

PURPOSE

This study aimed to determine if eccentric hamstring strength was prospectively associated with HSI among collegiate athletes, while controlling for sex, age, and prior HSI. We hypothesized that athletes with lower eccentric hamstring strength or greater between-limb strength asymmetry at preseason would have an increased risk of HSI.

METHODS

Hamstring eccentric strength measures, maximum total force ( FTotal ) and between-limb asymmetry in maximum force ( FAsym ), were measured at preseason on male and female athletes. HSIs were tracked over the subsequent 12 months. Generalized estimating equations were used to identify univariable and multivariable associations between athlete demographics, eccentric hamstring strength, and HSI risk.

RESULTS

Data for 326 athletes (85 female; 30 track, 43 basketball, 160 American football, 93 soccer) were included, and 64 HSIs were observed. Univariable associations between eccentric hamstring strength and subsequent HSI were nonsignificant ( FTotal : odds ratio [OR], 0.99 (95% confidence interval (CI), 0.93-1.05); P = 0.74; FAsym : OR, 1.35 (95% CI, 0.87-2.09); P = 0.23). No relationship between eccentric hamstring strength and HSI ( FAsym : OR, 1.32 (95% CI, 0.84-2.08); P = 0.23) was identified after adjusting for confounders including sex, age, and prior HSI.

CONCLUSIONS

No association between preseason eccentric hamstring strength and risk of subsequent HSI was identified after controlling for known risk factors and sex among collegiate athletes. Eccentric hamstring strengthening may continue to serve as a preventative approach to HSI, but it does not provide additional insight into HSI risk beyond factors such as age and prior HSI.

A hierarchical clustering approach for examining potential risk factors for bone stress injury in runners

Bone stress injuries (BSI) are overuse injuries that commonly occur in runners. BSI risk is multifactorial and not well understood. Unsupervised machine learning approaches can potentially elucidate risk factors for BSI by looking for groups of similar runners within a population that differ in BSI incidence. Here, a hierarchical clustering approach is used to identify groups of collegiate cross country runners (32 females, 21 males) based on healthy pre-season running (4.47 m·s^-1) gait data which were aggregated and dimensionally reduced by principal component analysis. Five distinct groups were identified using the cluster tree. Visual inspection revealed clear differences between groups in kinematics and kinetics, and linear mixed effects models showed between-group differences in metrics potentially related to BSI risk. The groups also differed in BSI incidence during the subsequent academic year (Rand index = 0.49; adjusted Rand index = -0.02). Groups ranged from those including runners spending less time contacting the ground and generating higher peak ground reaction forces and joint moments to those including runners spending more time on the ground with lower loads. The former groups showed higher BSI incidence, indicating that short stance phases and high peak loads may be risk factors for BSI. Since ground contact duration may itself account for differences in peak loading metrics, we hypothesize that the percentage of time a runner is in contact with the ground may be a useful metric to include in machine learning models for predicting BSI risk.

Spatial frequency analysis detects altered tissue organization following hamstring strain injury at time of injury but not return to sport

Background

Hamstring strain injury (HSI) diagnosis is often corroborated using ultrasound. Spatial frequency analysis (SFA) is a quantitative ultrasound method that has proven useful in characterizing altered tissue organization. The purpose of this study was to determine changes in muscular tissue organization using SFA following HSI.

Methods

Ultrasound B-mode images were captured at time of injury (TOI) and return to sport (RTS) in collegiate athletes who sustained an HSI. Spatial frequency parameters extracted from two-dimensional Fourier Transforms in user-defined regions of interest (ROI) were analyzed. Separate ROIs encompassed injured and adjacent tissue within the same image of the injured limb and mirrored locations in the contralateral limb at TOI. The ROIs for RTS images were drawn to correspond to the injury-matched location determined from TOI imaging. Peak spatial frequency radius (PSFR) and the fascicular banded pattern relative to image background (Mmax%) were compared between injured and adjacent portions within the same image with separate paired t-tests. Within-image differences of SFA parameters in the injured limb were calculated and compared between TOI and RTS with Wilcoxon rank sum tests.

Results

Within the injured limb at TOI, PSFR differences in injured and healthy regions did not strictly meet statistical significance (p = 0.06), while Mmax% was different between regions (p < 0.001). No differences were observed between regions in the contralateral limb at TOI (PSFR, p = 0.16; Mmax%, p = 0.30). Significant within-image differences in PSFR (p = 0.03) and Mmax% (p = 0.04) at RTS were detected relative to TOI.

Conclusions

These findings are a first step in determining the usefulness of SFA in muscle injury characterization and provide quantitative assessment of both fascicular disruption and edema presence in acute HSI.

Lower step rate is associated with a higher risk of bone stress injury: a prospective study of collegiate cross country runners

OBJECTIVES

To determine if running biomechanics and bone mineral density (BMD) were independently associated with bone stress injury (BSI) in a cohort of National Collegiate Athletic Association Division I cross country runners.

METHODS

This was a prospective, observational study of 54 healthy collegiate cross country runners over three consecutive seasons. Whole body kinematics, ground reaction forces (GRFs) and BMD measures were collected during the preseason over 3 years via motion capture on an instrumented treadmill and total body densitometer scans. All medically diagnosed BSIs up to 12 months following preseason data collection were recorded. Generalised estimating equations were used to identify independent risk factors of BSI.

RESULTS

Univariably, step rate, centre of mass vertical excursion, peak vertical GRF and vertical GRF impulse were associated with BSI incidence. After adjusting for history of BSI and sex in a multivariable model, a higher step rate was independently associated with a decreased risk of BSI. BSI risk decreased by 5% (relative risk (RR): 0.95; 95% CI 0.91 to 0.98) with each one step/min increase in step rate. BMD z-score was not a statistically significant risk predictor in the final multivariable model (RR: 0.93, 95% CI 0.85 to 1.03). No other biomechanical variables were found to be associated with BSI risk.

CONCLUSION

Low step rate is an important risk factor for BSI among collegiate cross country runners and should be considered when developing comprehensive programmes to mitigate BSI risk in distance runners.

Foot Angle and Loading Rate during Running Demonstrate a Nonlinear Relationship

Vertical loading rates are typically found to be lower in forefoot compared to rearfoot strikers, promoting the idea that forefoot striking is desirable and may reduce running injury risk. However, prior work using linear models has shown that foot inclination angle (FIA) at initial contact is a poor predictor of vertical loading rate, suggesting a more complex association exists.

PURPOSE

To determine if a nonlinear model superiorly describes the relationship between FIA and average vertical loading rate (AVLR). Secondary analyses assessed the influence of sex and sport on the association between FIA and AVLR.

METHODS

Whole body kinematics and vertical ground reaction forces were collected for 170 healthy National Collegiate Athletic Association Division I athletes (97 males; 81 cross-country runners) during treadmill running at 2.68, 3.35, and 4.47 m·s. Foot inclination angle and AVLR were calculated for 15 strides and averaged across strides for each limb. Polynomial mixed effects models assessed linear and nonlinear trends in the relationship between FIA and AVLR across the entire sample and accounting for sex and sport participation.

RESULTS

Average vertical loading rate was lowest at the extremes of FIA (i.e., -15°, 20°), whereas greater AVLR were observed between 5° and 10°. The cubic model resulted in a significantly better fit than the linear model (P < 0.001). Average vertical loading rate was also more variable among FIA associated with rearfoot and midfoot strike than forefoot strike. Adding sex to the model did not influence model fit; though, controlling for sport minimally improved model fit.

CONCLUSIONS

The relationship between FIA and AVLR is best represented by a cubic model. Consequently, FIA should be treated as a continuous variable. Reducing FIA into categories may misrepresent the relationship between FIA and other gait variables.

Increasing running step rate reduces patellofemoral joint forces

PURPOSE

Increasing step rate has been shown to elicit changes in joint kinematics and kinetics during running, and it has been suggested as a possible rehabilitation strategy for runners with patellofemoral pain. The purpose of this study was to determine how altering step rate affects internal muscle forces and patellofemoral joint loads, and then to determine what kinematic and kinetic factors best predict changes in joint loading.

METHODS

We recorded whole body kinematics of 30 healthy adults running on an instrumented treadmill at three step rate conditions (90%, 100%, and 110% of preferred step rate). We then used a 3-D lower extremity musculoskeletal model to estimate muscle, patellar tendon, and patellofemoral joint forces throughout the running gait cycles. In addition, linear regression analysis allowed us to ascertain the relative influence of limb posture and external loads on patellofemoral joint force.

RESULTS

Increasing step rate to 110% of the preferred reduced peak patellofemoral joint force by 14%. Peak muscle forces were also altered as a result of the increased step rate with hip, knee, and ankle extensor forces, and hip abductor forces all reduced in midstance. Compared with the 90% step rate condition, there was a concomitant increase in peak rectus femoris and hamstring loads during early and late swing, respectively, at higher step rates. Peak stance phase knee flexion decreased with increasing step rate and was found to be the most important predictor of the reduction in patellofemoral joint loading.

CONCLUSION

Increasing step rate is an effective strategy to reduce patellofemoral joint forces and could be effective in modulating biomechanical factors that can contribute to patellofemoral pain.

Web-based therapeutic exercise resource center as a treatment for knee osteoarthritis: a prospective cohort pilot study

BACKGROUND

Although beneficial effects of exercise in the management of knee osteoarthritis (OA) have been established, only 14 -18% of patients with knee OA receive an exercise from their primary care provider. Patients with knee OA cite lack of physician exercise advice as a major reason why they do not exercise to improve their condition. The purpose of this pilot study was to investigate use of a web-based Therapeutic Exercise Resource Center (TERC) as a tool to prescribe strength, flexibility and aerobic exercise as part of knee OA treatment. It was hypothesized that significant change in clinical outcome scores would result from patients' use of the TERC.

METHODS

Sixty five individuals diagnosed with mild/moderate knee OA based on symptoms and radiographs were enrolled through outpatient physician clinics. Using exercise animations to facilitate proper technique, the TERC assigned and progressed patients through multiple levels of exercise intensity based on exercise history, co-morbidities and a validated measure of pain and function. Subjects completed a modified short form WOMAC (mSF-WOMAC), World Health Organization Quality of Life (WHO-QOL) and Knee Self-Efficacy Scale (K-SES) at baseline and completion of the 8 week program, and a user satisfaction survey. Outcomes were compared over time using paired t-tests and effect sizes calculated using partial point biserial (pr).

RESULTS

Fifty two participants completed the 8 week program with average duration of knee pain 8.0 ± 11.0 yrs (25 females; 61.0 ± 9.4 yrs; body mass index, 28.8 ± 6.3 kg/m2). During the study period, all outcome measures improved: mSF-WOMAC scores decreased (better pain and function) (p<.001; large effect, pr=0.70); WHO-QOL physical scores increased (p=.015; medium effect, pr=0.33); and K-SES scores increased (p<.001; large effect, pr=0.54). No significant differences were found in study outcomes as a function of gender, age, BMI or symptom duration. Patients reported very positive evaluation of the TERC (94% indicated the website was easy to use; 90% specified the exercise animations were especially helpful).

CONCLUSION

This pilot study demonstrated the web-based TERC to be feasible and efficacious in improving clinical outcomes for patients with mild/moderate knee OA and supports future studies to compare TERC to current standard of care, such as educational brochures.

Changes in muscle activation patterns when running step rate is increased

Running with a step rate 5-10% greater than one's preferred can substantially reduce lower extremity joint moments and powers, and has been suggested as a possible strategy to aid in running injury management. The purpose of this study was to examine how neuromuscular activity changes with an increase in step rate during running. Forty-five injury-free, recreational runners participated in this study. Three-dimensional motion, ground reaction forces, and electromyography (EMG) of 8 muscles (rectus femoris, vastus lateralis, medial gastrocnemius, tibialis anterior, medial and lateral hamstrings, and gluteus medius and maximus) were recorded as each subject ran at their preferred speed for three different step rate conditions: preferred, +5% and +10% of preferred. Outcome measures included mean normalized EMG activity for each muscle at specific periods during the gait cycle. Muscle activities were found to predominantly increase during late swing, with no significant change in activities during the loading response. This increased muscle activity in anticipation of foot-ground contact likely alters the landing posture of the limb and the subsequent negative work performed by the joints during stance phase. Further, the increased activity observed in the gluteus maximus and medius suggests running with a greater step rate may have therapeutic benefits to those with anterior knee pain.

Effects of step rate manipulation on joint mechanics during running

PURPOSE

the objective of this study was to characterize the biomechanical effects of step rate modification during running on the hip, knee, and ankle joints so as to evaluate a potential strategy to reduce lower extremity loading and risk for injury.

METHODS

three-dimensional kinematics and kinetics were recorded from 45 healthy recreational runners during treadmill running at constant speed under various step rate conditions (preferred, ± 5%, and ± 10%). We tested our primary hypothesis that a reduction in energy absorption by the lower extremity joints during the loading response would occur, primarily at the knee, when step rate was increased.

RESULTS

less mechanical energy was absorbed at the knee (P < 0.01) during the +5% and +10% step rate conditions, whereas the hip (P < 0.01) absorbed less energy during the +10% condition only. All joints displayed substantially (P < 0.01) more energy absorption when preferred step rate was reduced by 10%. Step length (P < 0.01), center of mass vertical excursion (P < 0.01), braking impulse (P < 0.01), and peak knee flexion angle (P < 0.01) were observed to decrease with increasing step rate. When step rate was increased 10% above preferred, peak hip adduction angle (P < 0.01) and peak hip adduction (P < 0.01) and internal rotation (P < 0.01) moments were found to decrease.

CONCLUSION

we conclude that subtle increases in step rate can substantially reduce the loading to the hip and knee joints during running and may prove beneficial in the prevention and treatment of common running-related injuries.