Local Patterns in 2-Year T1ρ and T2 Changes of Hip Cartilage Are Related to Sex and Functional Data: A Prospective Evaluation on Hip Osteoarthritis Participants.

BACKGROUND: Although T1ρ and T2 have emerged as early indicators for hip osteoarthritis (OA), there is little information regarding longitudinal changes across the cartilage in the early stages of this disease. PURPOSE: To characterize the variability in 2-year hip cartilage T1ρ and T2 changes and investigate associations between these patterns of change and common indicators of hip OA. STUDY TYPE: Prospective. POPULATION: A total of 25 women (age: 51.9 ± 16.3 years old; BMI: 22.6 ± 2.0 kg/m2 ) and 17 men (age: 55.8 ± 14.9 years old; body mass index (BMI): 24.4 ± 3.8 kg/m2 ) who were healthy or with early-to-moderate hip OA. FIELD STRENGTH/SEQUENCE: A 3 T MRI (GE), 3D combined T1ρ /T2 magnetization-prepared angle-modulated partitioned k-space spoiled gradient echo snapshots. ASSESSMENT: Principal component (PC) analysis of Z-score difference maps of 2-year changes in hip cartilage T1ρ and T2 relaxation times, participant hip disability and osteoarthritis outcome scores (HOOS) and functional tests at 2-year follow-up. STATISTICAL TESTS: Shapiro-Wilk test, unpaired t-tests, Kruskal Wallis tests, Pearson or Spearman (ρ) correlations. Significance was set at P < 0.05. RESULTS: Women (-6.40 ± 14.48) had significantly lower T1ρ PC1 scores than men (10.05 ± 26.15). T1ρ PC4 was significantly correlated with HOOSsport , HOOSsymptoms , HOOSpain , HOOSadl , and HOOSqol at 2-year follow-up (ρ: [0.36, 0.50]). T1ρ PC2 and PC4 were significantly correlated with 30-second chair test (ρ = -0.39 and ρ = 0.24, respectively) and side plank (ρ = -0.32 and ρ = 0.21). T1ρ and T2 PC2 were significantly correlated with 40 m walk test (ρ = 0.34 and ρ = 0.31) and 30-second chair rise test (ρ = -0.39 and ρ = -0.32). DATA CONCLUSION: Men exhibited accelerated T1ρ increases across the femoral cartilage compared to women, suggesting sex should be considered when evaluating early hip OA. Participants with poorer HOOS and function exhibited greater T1ρ and T2 increases in superior and anterior femoral cartilage and greater T1ρ increases in the anterior femoral cartilage. These patterns of short-term relaxometry increases could indicate hip OA progression. EVIDENCE LEVEL: 2 TECHNICAL EFFICACY: Stage 3.

Effects of T1p Characteristics of Load-Bearing Hip Cartilage on Bilateral Knee Patellar Cartilage Subregions: Subjects With None to Moderate Radiographic Hip Osteoarthritis.

BACKGROUND: The polyarticular nature of Osteoarthritis (OA) tends to manifest in multi-joints. Associations between cartilage health in connected joints can help identify early degeneration and offer the potential for biomechanical intervention. Such associations between hip and knee cartilages remain understudied. PURPOSE: To investigate T1p associations between hip-femoral and acetabular-cartilage subregions with Intra-limb and Inter-limb patellar cartilage; whole and deep-medial (DM), deep-lateral (DL), superficial-medial (SM), superficial-lateral (SL) subregions. STUDY TYPE: Prospective. SUBJECTS: Twenty-eight subjects (age 55.1 ± 12.8 years, 15 females) with none-to-moderate hip-OA while no radiographic knee-OA. FIELD STRENGTH/SEQUENCE: 3-T, bilateral hip, and knee: 3D-proton-density-fat-saturated (PDFS) Cube and Magnetization-Prepared-Angle-Modulated-Partitioned-k-Space-Spoiled-Gradient-Echo-Snapshots (MAPSS). ASSESSMENT: Ages of subjects were categorized into Group-1 (≤40), Group-2 (41-50), Group-3 (51-60), Group-4 (61-70), Group-5 (71-80), and Group-6 (≥81). Hip T1p maps, co-registered to Cube, underwent an atlas-based algorithm to quantify femoral and acetabular subregional (R2 -R7 ) cartilage T1p . For knee Cube, a combination of V-Net architectures was used to segment the patellar cartilage and subregions (DM, DL, SM, SL). T1p values were computed from co-registered MAPSS. STATISTICAL TESTS: For Intra-and-Inter-limb, 5 optimum predictors out of 13 (Hip subregional T1p , age group, gender) were selected by univariate linear-regression, to predict outcome (patellar T1p ). The top five predictors were stepwise added to six linear mixed-effect (LME) models. In all LME models, we assume the data come from the same subject sharing the same random effect. The best-performing models (LME-modelbest ) selected via ANOVA, were tested with DM, SM, SL, and DL subregional-mean T1p . LME assumptions were verified (normality of residuals, random-effects, and posterior-predictive-checks). RESULTS: LME-modelbest (Intra-limb) had significant negative and positive fixed-effects of femoral-R5 and acetabular-R2 T1p , respectively (conditional-R2 = 0.581). LME-modelbest (Inter-limb) had significant positive fixed-effects of femoral-R3 T1p (conditional-R2 = 0.26). DATA CONCLUSION: Significant positive and negative T1p associations were identified between load-bearing hip cartilage-subregions vs. ipsilateral and contralateral patellar cartilages respectively. The effects were localized on medial subregions of Inter-limb, in particular. EVIDENCE LEVEL: 1 TECHNICAL EFFICACY: Stage 1.

In Vivo Kinematics of the Tibiotalar and Subtalar Joints in Asymptomatic Subjects: A High-Speed Dual Fluoroscopy Study.

Measurements of joint kinematics are essential to understand the pathomechanics of ankle disease and the effects of treatment. Traditional motion capture techniques do not provide measurements of independent tibiotalar and subtalar joint motion. In this study, high-speed dual fluoroscopy images of ten asymptomatic adults were acquired during treadmill walking at 0.5 m/s and 1.0 m/s and a single-leg, balanced heel-rise. Three-dimensional (3D) CT models of each bone and dual fluoroscopy images were used to quantify in vivo kinematics for the tibiotalar and subtalar joints. Dynamic tibiotalar and subtalar mean joint angles often exhibited opposing trends during captured stance. During both speeds of walking, the tibiotalar joint had significantly greater dorsi/plantarflexion (D/P) angular ROM than the subtalar joint while the subtalar joint demonstrated greater inversion/eversion (In/Ev) and internal/external rotation (IR/ER) than the tibiotalar joint. During balanced heel-rise, only D/P and In/Ev were significantly different between the tibiotalar and subtalar joints. Translational ROM in the anterior/posterior (AP) direction was significantly greater in the subtalar than the tibiotalar joint during walking at 0.5 m/s. Overall, our results support the long-held belief that the tibiotalar joint is primarily responsible for D/P, while the subtalar joint facilitates In/Ev and IR/ER. However, the subtalar joint provided considerable D/P rotation, and the tibiotalar joint rotated about all three axes, which, along with translational motion, suggests that each joint undergoes complex, 3D motion.

Automated evaluation of hip abductor muscle quality and size in hip osteoarthritis: Localized muscle regions are strongly associated with overall muscle quality.

Limited information exists regarding abductor muscle quality variation across its length and which locations are most representative of overall muscle quality. This is exacerbated by time-intensive processes for manual muscle segmentation, which limits feasibility of large cohort analyses. The purpose of this study was to develop an automated and localized analysis pipeline that accurately estimates hip abductor muscle quality and size in individuals with mild-to-moderate hip osteoarthritis (OA) and identifies regions of each muscle which provide best estimates of overall muscle quality. Forty-four participants (age 52.7 ± 16.1 years, BMI 23.7 ± 3.4 kg/m2, 14 males) with and without mild-to-moderate radiographic hip OA were recruited for this study. Unilateral hip magnetic resonance (MR) images were acquired on a 3.0 T MR scanner and included axial T1-weighted fast spin echo and 3D axial Iterative Decomposition of water and fat with Echo Asymmetry and Least-squares estimation (IDEAL-IQ) spoiled gradient-recalled echo (SPGR) with multi-peak fat spectrum modeling and single T2* correction. A three dimensional (3D) V-Net convolutional neural network was trained to automatically segment the gluteus medius (GMED), gluteus minimus (GMIN), and tensor fascia lata (TFL) on axial IDEAL-IQ. Agreement between manual and automatic segmentation and associations between axial fat fraction (FF) estimated from IDEAL-IQ and overall muscle FF were evaluated. Dice scores for automatic segmentation were 0.94, 0.87, and 0.91 for GMED, GMIN, and TFL, respectively. GMED, GMIN, and TFL volumetric and FF measures were strongly correlated (r: 0.92-0.99) between automatic and manual segmentations, where all values fell within the 95% limits of agreement of [-9.79 cm3, 17.43 cm3] and [-1.99%, 2.89%], respectively. Axial FF was significantly associated with overall FF with the strongest correlations at 50%, 50%, and 65% the length of the GMED, GMIN, and TFL muscles, respectively (r: 0.93-0.97). An automated and localized analysis can provide efficient and accurate estimates of hip abductor muscle quality and size across muscle length. Specific regions of the muscle may be used to estimate overall muscle quality in an abbreviated evaluation of muscle quality.

Total Ankle Replacement Provides Symmetrical Postoperative Kinematics: A Biplane Fluoroscopy Imaging Study.

BACKGROUND: In vivo measurements of tibiotalar and subtalar joint motion following TAR are unavailable. Using biplane fluoroscopy, we tested the hypothesis that the prosthetic tibiotalar joint and adjacent subtalar joint would demonstrate kinematic and range of motion differences compared to the contralateral untreated limb, and control participants. METHODS: Six patients of 41 identified candidates that all underwent unilateral Zimmer TAR (5.4 ± 1.9 years prior) and 6 control participants were imaged with biplane fluoroscopy during overground walking and a double heel-rise activity. Computed tomography scans were acquired; images were segmented and processed to serve as input for model-based tracking of the biplane fluoroscopy data. Measurements included tibiotalar and subtalar kinematics for the TAR, untreated contralateral, and control limbs. Statistical parametric mapping quantified differences in kinematics throughout overground walking and the double heel-rise activity. RESULTS: Patients with this TAR performed walking and heel-rise activities symmetrically with no significant kinematic differences at the tibiotalar and subtalar joints between limbs. Compared to control participants, patients exhibited reduced dorsi/plantarflexion range of motion that corresponded to decreased peak dorsiflexion, but only in the late stance phase of walking. This reduction in tibiotalar dorsi/plantarflexion range of motion in the TAR group became more apparent with double heel-rise activity. CONCLUSION: Patients with a Zimmer TAR had symmetric kinematics during activities of walking and double heel-rise, but they did exhibit minor compensations in tibiotalar kinematics as compared to controls. CLINICAL RELEVANCE: The lack of significant kinematic compensation at the subtalar joint may explain why secondary subtalar osteoarthritis is reported as being relatively uncommon in patients with some TAR designs.

Subject-specific biomechanical analysis to estimate locations susceptible to osteoarthritis-Finite element modeling and MRI follow-up of ACL reconstructed patients.

The aims of this case-control study were to: (1) Identify cartilage locations and volumes at risk of osteoarthritis (OA) using subject-specific finite element (FE) models; (2) Quantify the relationships between the simulated biomechanical parameters and T2 and T1ρ relaxation times of magnetic resonance imaging (MRI). We created subject-specific FE models for seven patients with anterior cruciate ligament (ACL) reconstruction and six controls based on a previous proof-of-concept study. We identified locations and cartilage volumes susceptible to OA, based on maximum principal stresses and absolute maximum shear strains in cartilage exceeding thresholds of 7 MPa and 32%, respectively. The locations and volumes susceptible to OA were compared qualitatively and quantitatively against 2-year longitudinal changes in T2 and T1ρ relaxation times. The degeneration volumes predicted by the FE models, based on excessive maximum principal stresses, were significantly correlated (r = 0.711, p < 0.001) with the degeneration volumes determined from T2 relaxation times. There was also a significant correlation between the predicted stress values and changes in T2 relaxation time (r = 0.649, p < 0.001). Absolute maximum shear strains and changes in T1ρ relaxation time were not significantly correlated. Five out of seven patients with ACL reconstruction showed excessive maximum principal stresses in either one or both tibial cartilage compartments, in agreement with follow-up information from MRI. Expectedly, for controls, the FE models and follow-up information showed no degenerative signs. Our results suggest that the presented modelling methodology could be applied to prospectively identify ACL reconstructed patients at risk of biomechanically driven OA, particularly by the analysis of maximum principal stresses of cartilage.

Multivariate functional principal component analysis identifies waveform features of gait biomechanics related to early-to-moderate hip osteoarthritis.

Clinicians often examine movement patterns to design hip osteoarthritis (OA) interventions, yet traditional biomechanical analyses only report a single timepoint. Multivariate principal component analysis (MFPCA) analyzes the entire waveform (i.e., movement pattern), which clinicians observe to direct treatment. This study investigated hip OA indicators, by (1) employing MFPCA to characterize variance across the hip, knee, and ankle angles in healthy and early-to-moderate hip OA participants; and (2) investigating relationships between these waveform features and hip cartilage health. Bilateral hip magnetic resonance images from 72 participants with Kellgren-Lawrence grades ranging from 0 to 3 were used to calculate mean T 1ρ and T 2 relaxation times in the femoral and acetabular cartilage. MFPCA was performed on lower-limb gait biomechanics and used to identify primary modes of variation, which were related to T 1ρ and T 2 relaxation times. Here, a MFPC = mode of variation = waveform feature. In the femoral cartilage, transverse plane MFPCs 3 and 5 and body mass index (BMI) was related to T 1ρ , while MFPC 2 and BMI were related to T 2 relaxation times. In the acetabular cartilage, sagittal plane MFPC 1 and BMI were related to T 1ρ , while BMI was related to T 2 relaxation times. Greater internal rotation was related to increased T 1ρ and T 2 relaxation times in the femoral cartilage, while the greater extension was related to increased T 1ρ relaxation times in the acetabular cartilage. This study established a data-driven framework to assess relationships between multi-joint biomechanics and quantitative assessments of cartilage health and identified waveform features that could be evaluated in future hip OA intervention studies.

The modified Shriners Hospitals for Children Greenville (mSHCG) multi-segment foot model provides clinically acceptable measurements of ankle and midfoot angles: A dual fluoroscopy study.

BACKGROUND: Several multi-segment foot models have been developed to evaluate foot and ankle motion using skin-marker motion analysis. However, few multi-segment models have been evaluated against a reference standard to establish kinematic accuracy. RESEARCH QUESTION: How accurately do skin-markers estimate foot and ankle motion for the modified Shriners Hospitals for Children Greenville (mSHCG) multi-segment foot model when compared against the reference standard, dual fluoroscopy (DF), during gait, in asymptomatic participants? METHODS: Five participants walked overground as full-body skin-marker trajectory data and DF images of the foot and shank were simultaneously acquired. Using the mSHCG model, ankle and midfoot angles were calculated throughout stance for both motion analysis techniques. Statistical parametric mapping assessed differences in joint angles and marker positions between skin-marker and DF motion analysis techniques. Paired t tests, and linear regression models were used to compare joint angles and range of motion (ROM) calculated from the two techniques. RESULTS: In the coronal plane, the skin-marker model significantly overestimated ROM (p = 0.028). Further, the DF model midfoot ROM was significantly positively related to differences between DF and skin-marker midfoot angles (p = 0.035, adjusted R2 = 0.76). In the sagittal plane, skin-markers underestimated ankle angles by as much as 7.26°, while midfoot angles were overestimated by as much as 9.01°. However, DF and skin-marker joint angles were not significantly different over stance. Skin-markers on the tibia, calcaneus, and fifth metatarsal had significantly different positions than the DF markers along the direction of walking for isolated portions that were less than 10 % of stance. Euclidean distances between DF and skin-markers positions were less than 9.36 mm. SIGNIFICANCE: As the accuracy of the mSHCG model was formerly unknown, the results of this study provide ranges of confidence for key angles calculated by this model.

Compensatory Motion of the Subtalar Joint Following Tibiotalar Arthrodesis: An in Vivo Dual-Fluoroscopy Imaging Study.

BACKGROUND: Tibiotalar arthrodesis is a common treatment for end-stage tibiotalar osteoarthritis, and is associated with a long-term risk of concomitant subtalar osteoarthritis. It has been clinically hypothesized that subtalar osteoarthritis following tibiotalar arthrodesis is the product of compensatory subtalar joint hypermobility. However, in vivo measurements of subtalar joint motion following tibiotalar arthrodesis have not been quantified. Using dual-fluoroscopy motion capture, we tested the hypothesis that the subtalar joint of the limb with a tibiotalar arthrodesis would demonstrate differences in kinematics and increased range of motion compared with the subtalar joint of the contralateral, asymptomatic, untreated ankle. METHODS: Ten asymptomatic patients who had undergone unilateral tibiotalar arthrodesis at a mean (and standard deviation) of 4.0 ± 1.8 years previously were evaluated during overground walking and a double heel-rise task. The evaluation involved markerless tracking with use of dual fluoroscopy integrated with 3-dimensional computed tomography, which allowed for dynamic measurements of subtalar and tibiotalar dorsiflexion-plantar flexion, inversion-eversion, and internal-external rotation. Range of motion, stance time, swing time, step length, and step width were also measured. RESULTS: During the early stance phase of walking, the subtalar joint of the limb that had been treated with arthrodesis was plantar flexed (-4.7° ± 3.3°), whereas the subtalar joint of the untreated limb was dorsiflexed (4.6° ± 2.2°). Also, during the early stance phase of walking, eversion of the subtalar joint of the surgically treated limb (0.2° ± 2.3°) was less than that of the untreated limb (4.5° ± 3.2°). During double heel-rise, the treated limb exhibited increased peak subtalar plantar flexion (-7.1° ± 4.1°) compared with the untreated limb (0.2° ± 1.8°). CONCLUSIONS: A significant increase in subtalar joint plantar flexion was found to be a primary compensation during overground walking and a double heel-rise activity following tibiotalar arthrodesis. CLINICAL RELEVANCE: Significant subtalar joint plantar flexion compensations appear to occur following tibiotalar arthrodesis. We found an increase in subtalar plantar flexion and considered the potential relationship of this finding with the increased rate of subtalar osteoarthritis that occurs following ankle arthrodesis.

Subject-Specific Axes of Rotation Based on Talar Morphology Do Not Improve Predictions of Tibiotalar and Subtalar Joint Kinematics.

Use of subject-specific axes of rotation may improve predictions generated by kinematic models, especially for joints with complex anatomy, such as the tibiotalar and subtalar joints of the ankle. The objective of this study was twofold. First, we compared the axes of rotation between generic and subject-specific ankle models for ten control subjects. Second, we quantified the accuracy of generic and subject-specific models for predicting tibiotalar and subtalar joint motion during level walking using inverse kinematics. Here, tibiotalar and subtalar joint kinematics measured in vivo by dual-fluoroscopy served as the reference standard. The generic model was based on a cadaver study, while the subject-specific models were derived from each subject's talus reconstructed from computed tomography images. The subject-specific and generic axes of rotation were significantly different. The average angle between the modeled axes was 12.9° ± 4.3° and 24.4° ± 5.9° at the tibiotalar and subtalar joints, respectively. However, predictions from both models did not agree well with dynamic dual-fluoroscopy data, where errors ranged from 1.0° to 8.9° and 0.6° to 7.6° for the generic and subject-specific models, respectively. Our results suggest that methods that rely on talar morphology to define subject-specific axes may be inadequate for accurately predicting tibiotalar and subtalar joint kinematics.

Application of High-Speed Dual Fluoroscopy to Study In Vivo Tibiotalar and Subtalar Kinematics in Patients With Chronic Ankle Instability and Asymptomatic Control Subjects During Dynamic Activities.

BACKGROUND: Abnormal angular and translational (ie, kinematic) motion at the tibiotalar and subtalar joints is believed to cause osteoarthritis in patients with chronic ankle instability (CAI). METHODS: In this preliminary study the investigators quantified and compared in vivo tibiotalar and subtalar kinematics in 4 patients with CAI (3 women) and 10 control subjects (5 men) using dual fluoroscopy during a balanced, single-leg heel-rise and treadmill walking at 0.5 and 1.0 m/s. RESULTS: During balanced heel-rise, 69%, 54%, and 66% of mean CAI tibiotalar internal rotation/external rotation (IR/ER), subtalar inversion/eversion, and subtalar IR/ER angles, respectively, were outside the 95% confidence intervals of control subjects. During 0.5-m/s gait, 50% and 60% of mean CAI tibiotalar dorsi/plantarflexion and subtalar IR/ER angles, respectively, were outside the 95% confidence intervals of control subjects. During 1.0-m/s gait, 62%, 65%, and 73% of mean CAI subtalar dorsi/plantarflexion, inversion/eversion, and IR/ER, respectively, were outside the 95% confidence intervals of control subjects. Patients with CAI exhibited less tibiotalar and subtalar translational motion during gait; no clear differences in translations were noted during balanced heel-rise. CONCLUSION: Overall, the balanced heel-rise activity exposed more tibiotalar and subtalar kinematic variation between patients with CAI and control subjects. Therefore, weight-bearing activities involving large range of motion, balance, and stability may be best for studying kinematic adaptations in patients with CAI. CLINICAL RELEVANCE: These preliminary results suggest that patients with CAI require more tibiotalar external rotation, subtalar eversion, and subtalar external rotation during weight-bearing stability exercises, all with less overall joint translation.

Predicting tibiotalar and subtalar joint angles from skin-marker data with dual-fluoroscopy as a reference standard.

Evidence suggests that the tibiotalar and subtalar joints provide near six degree-of-freedom (DOF) motion. Yet, kinematic models frequently assume one DOF at each of these joints. In this study, we quantified the accuracy of kinematic models to predict joint angles at the tibiotalar and subtalar joints from skin-marker data. Models included 1 or 3 DOF at each joint. Ten asymptomatic subjects, screened for deformities, performed 1.0m/s treadmill walking and a balanced, single-leg heel-rise. Tibiotalar and subtalar joint angles calculated by inverse kinematics for the 1 and 3 DOF models were compared to those measured directly in vivo using dual-fluoroscopy. Results demonstrated that, for each activity, the average error in tibiotalar joint angles predicted by the 1 DOF model were significantly smaller than those predicted by the 3 DOF model for inversion/eversion and internal/external rotation. In contrast, neither model consistently demonstrated smaller errors when predicting subtalar joint angles. Additionally, neither model could accurately predict discrete angles for the tibiotalar and subtalar joints on a per-subject basis. Differences between model predictions and dual-fluoroscopy measurements were highly variable across subjects, with joint angle errors in at least one rotation direction surpassing 10° for 9 out of 10 subjects. Our results suggest that both the 1 and 3 DOF models can predict trends in tibiotalar joint angles on a limited basis. However, as currently implemented, neither model can predict discrete tibiotalar or subtalar joint angles for individual subjects. Inclusion of subject-specific attributes may improve the accuracy of these models.

Accuracy and feasibility of high-speed dual fluoroscopy and model-based tracking to measure in vivo ankle arthrokinematics.

The relationship between altered tibiotalar and subtalar kinematics and development of ankle osteoarthritis is unknown, as skin marker motion analysis cannot measure articulations of each joint independently. Here, we quantified the accuracy and demonstrated the feasibility of high-speed dual fluoroscopy (DF) to measure and visualize the three-dimensional articulation (i.e., arthrokinematics) of the tibiotalar and subtalar joints. Metal beads were implanted in the tibia, talus and calcaneus of two cadavers. Three-dimensional surface models of the cadaver and volunteer bones were reconstructed from computed tomography images. A custom DF system was positioned adjacent to an instrumented treadmill. DF images of the cadavers were acquired during maximal rotation about three axes (dorsal-plantar flexion, inversion-eversion, internal-external rotation) and simulated gait (treadmill at 0.5 and 1.0 m/s). Positions of implanted beads were tracked using dynamic radiostereometric analysis (DRSA). Bead locations were also calculated using model-based markerless tracking (MBT) and compared, along with joint angles and translations, to DRSA results. The mean positional difference between DRSA and MBT for all frames defined bias; standard deviation of the difference defined precision. The volunteer was imaged with DF during treadmill gait. From these movements, joint kinematics and tibiotalar and subtalar bone-to-bone distance were calculated. The mean positional and rotational bias (±standard deviation) of MBT was 0.03±0.35 mm and 0.25±0.81°, respectively. Mean translational and rotational precision was 0.30±0.12 mm and 0.63±0.28°, respectively. With excellent measurement accuracy, DF and MBT may elucidate the kinematic pathways responsible for osteoarthritis of the tibiotalar and subtalar joints in living subjects.