Automatic registration of MRI-based joint models to high-speed biplanar radiographs for precise quantification of in vivo anterior cruciate ligament deformation during gait

Scan-Free and Fully Automatic Tracking of Native Knee Anatomy from Dynamic Stereo-Radiography with Statistical Shape and Intensity Models

Kinematic tracking of native anatomy from stereo-radiography provides a quantitative basis for evaluating human movement. Conventional tracking procedures require significant manual effort and call for acquisition and annotation of subject-specific volumetric medical images. The current work introduces a framework for fully automatic tracking of native knee anatomy from dynamic stereo-radiography which forgoes reliance on volumetric scans. The method consists of three computational steps. First, captured radiographs are annotated with segmentation maps and anatomic landmarks using a convolutional neural network. Next, a non-convex polynomial optimization problem formulated from annotated landmarks is solved to acquire preliminary anatomy and pose estimates. Finally, a global optimization routine is performed for concurrent refinement of anatomy and pose. An objective function is maximized which quantifies similarities between masked radiographs and digitally reconstructed radiographs produced from statistical shape and intensity models. The proposed framework was evaluated against manually tracked trials comprising dynamic activities, and additional frames capturing a static knee phantom. Experiments revealed anatomic surface errors routinely below 1.0 mm in both evaluation cohorts. Median absolute errors of individual bone pose estimates were below 1.0∘ or mm for 15 out of 18 degrees of freedom in both evaluation cohorts. Results indicate that accurate pose estimation of native anatomy from stereo-radiography may be performed with significantly reduced manual effort, and without reliance on volumetric scans.

Fully automatic tracking of native glenohumeral kinematics from stereo-radiography

The current work introduces a system for fully automatic tracking of native glenohumeral kinematics in stereo-radiography sequences. The proposed method first applies convolutional neural networks to obtain segmentation and semantic key point predictions in biplanar radiograph frames. Preliminary bone pose estimates are computed by solving a non-convex optimization problem with semidefinite relaxations to register digitized bone landmarks to semantic key points. Initial poses are then refined by registering computed tomography-based digitally reconstructed radiographs to captured scenes, which are masked by segmentation maps to isolate the shoulder joint. A particular neural net architecture which exploits subject-specific geometry is also introduced to improve segmentation predictions and increase robustness of subsequent pose estimates. The method is evaluated by comparing predicted glenohumeral kinematics to manually tracked values from 17 trials capturing 4 dynamic activities. Median orientation differences between predicted and ground truth poses were 1.7∘ and 8.6∘ for the scapula and humerus, respectively. Joint-level kinematics differences were less than 2∘ in 65%, 13%, and 63% of frames for XYZ orientation DoFs based on Euler angle decompositions. Automation of kinematic tracking can increase scalability of tracking workflows in research, clinical, or surgical applications.

Elevated In Vivo ACL Strain Is Associated With a Straight Knee in Both the Sagittal and the Coronal Planes

BACKGROUND

Noncontact anterior cruciate ligament (ACL) injuries typically occur during deceleration movements such as landing or cutting. However, conflicting data have left the kinematic mechanisms leading to these injuries unclear. Quantifying the influence of sagittal and coronal plane knee kinematics on in vivo ACL strain may help to elucidate noncontact ACL injury mechanisms.

PURPOSE/HYPOTHESIS

The purpose of this study was to measure in vivo sagittal and coronal plane knee kinematics and ACL strain during a single-leg jump. We hypothesized that ACL strain would be modulated primarily by motion in the sagittal plane and that limited coronal plane motion would be measured during this activity.

STUDY DESIGN

Descriptive laboratory study.

METHODS

Seventeen healthy participants (8 male/9 female) underwent magnetic resonance imaging (MRI) followed by high-speed biplanar radiography, obtained as participants performed a single-leg jump. Three-dimensional models of the femur, tibia, and associated ACL attachment site footprints were created from the MRIs and registered to the radiographs to reproduce the position of the knee during the jump. ACL strain, knee flexion/extension angles, and varus/valgus angles were measured throughout the jump. Spearman rank correlations were used to assess relationships between mean ACL strain and kinematic variables.

RESULTS

Mean ACL strain increased with decreasing knee flexion angle (ρ = -0.3; P = .002), and local maxima in ACL strain occurred with the knee in a straight position in both the sagittal and the coronal planes. In addition, limited coronal plane motion (varus/valgus angle) was measured during this activity (mean ± SD, -0.5°± 0.3°). Furthermore, we did not detect a statistically significant relationship between ACL strain and varus/valgus angle (ρ = -0.01; P = .9).

CONCLUSION

ACL strain was maximized when the knee was in a straight position in both the sagittal and coronal planes. Participants remained in <1° of varus/valgus position on average throughout the jump. As a ligament under elevated strain is more vulnerable to injury, landing on a straight knee may be an important risk factor for ACL rupture.

CLINICAL RELEVANCE

These data may improve understanding of risk factors for noncontact ACL injury, which may be useful in designing ACL injury prevention programs.

Use of a Novel Multimodal Imaging Technique to Model In Vivo Quadriceps Force and ACL Strain During Dynamic Activity

BACKGROUND

Quadriceps loading of the anterior cruciate ligament (ACL) may play a role in the noncontact mechanism of ACL injury. Musculoskeletal modeling techniques are used to estimate the intrinsic force of the quadriceps acting at the knee joint.

PURPOSE/HYPOTHESIS

The purpose of this paper was to develop a novel musculoskeletal model of in vivo quadriceps force during dynamic activity. We used the model to estimate quadriceps force in relation to ACL strain during a single-leg jump. We hypothesized that quadriceps loading of the ACL would reach a local maximum before initial ground contact with the knee positioned in extension.

STUDY DESIGN

Descriptive laboratory study.

METHODS

Six male participants underwent magnetic resonance imaging in addition to high-speed biplanar radiography during a single-leg jump. Three-dimensional models of the knee joint, including the femur, tibia, patellofemoral cartilage surfaces, and attachment-site footprints of the patellar tendon, quadriceps tendon, and ACL, were created from the magnetic resonance imaging scans. The bone models were registered to the biplanar radiographs, thereby reproducing the positions of the knee joint at the time of radiographic imaging. The magnitude of quadriceps force was determined for each knee position based on a 3-dimensional balance of the forces and moments of the patellar tendon and the patellofemoral cartilage contact acting on the patella. Knee kinematics and ACL strain were determined for each knee position.

RESULTS

A local maximum in average quadriceps force of approximately 6500 N (8.4× body weight) occurred before initial ground contact. ACL strain increased concurrently with quadriceps force when the knee was positioned in extension.

CONCLUSION

This novel participant-specific modeling technique provides estimates of in vivo quadriceps force during physiologic dynamic loading. A local maximum in quadriceps force before initial ground contact may tension the ACL when the knee is positioned in extension.

CLINICAL RELEVANCE

These data contribute to understanding noncontact ACL injury mechanisms and the potential role of quadriceps activation in these injuries.

[Evaluation of a Bone Coordinate System Constructed Using MR Image Composing]

PURPOSE

To evaluate the accuracy of a bone coordinate system constructed using MR image composing.

METHOD

A femoral coordinate system constructed using image composing of MR images of a whole bovine femur was evaluated using CT images. The MR images were acquired by moving the table and were processed with 3D distortion correction and composing. To evaluate the reproducibility of the measurements, the same operator repeated the construction of the femoral coordinate system. In addition, distortions in the MR images were evaluated in comparison with those in the CT images.

RESULT

The center position of the femoral coordinate system constructed using the MR image composing was 1.6±0.9 mm on the X-axis, 1.5±0.8 mm on the Y-axis, and 0.2±0.3 mm on the Z-axis, and the rotation of each axis was 1° or less. The distortion of the composed MR image was about 0.3%.

CONCLUSION

The femoral coordinate system constructed using MR image composing had the same accuracy as a system constructed with CT images. The effect of MR image composing on the construction of the femoral coordinate system was small.

Design and validation of a semi-automatic bone segmentation algorithm from MRI to improve research efficiency

Segmentation of medical images into different tissue types is essential for many advancements in orthopaedic research; however, manual segmentation techniques can be time- and cost-prohibitive. The purpose of this work was to develop a semi-automatic segmentation algorithm that leverages gradients in spatial intensity to isolate the patella bone from magnetic resonance (MR) images of the knee that does not require a training set. The developed algorithm was validated in a sample of four human participants (in vivo) and three porcine stifle joints (ex vivo) using both magnetic resonance imaging (MRI) and computed tomography (CT). We assessed the repeatability (expressed as mean ± standard deviation) of the semi-automatic segmentation technique on: (1) the same MRI scan twice (Dice similarity coefficient = 0.988 ± 0.002; surface distance = - 0.01 ± 0.001 mm), (2) the scan/re-scan repeatability of the segmentation technique (surface distance = - 0.02 ± 0.03 mm), (3) how the semi-automatic segmentation technique compared to manual MRI segmentation (surface distance = - 0.02 ± 0.08 mm), and (4) how the semi-automatic segmentation technique compared when applied to both MRI and CT images of the same specimens (surface distance = - 0.02 ± 0.06 mm). Mean surface distances perpendicular to the cartilage surface were computed between pairs of patellar bone models. Critically, the semi-automatic segmentation algorithm developed in this work reduced segmentation time by approximately 75%. This method is promising for improving research throughput and potentially for use in generating training data for deep learning algorithms.

Translation and rotation analysis based on stress MRI for the diagnosis of anterior cruciate ligament tears

BACKGROUND

Due to the increasing need for a detailed biomechanical analysis of anterior cruciate ligament (ACL) lesions, the aim of the study was to develop a method of direct measurement of the three-dimensional tibial translation and rotation based on stress MRI.

METHODS

For the purpose of the study, thirty patients with acute ACL rupture and 17 healthy control subjects were selected. Based on clinical examination, they were qualified for MRI examination using the Arthroholder Device prototype to perform anterior tibial translation. Each examination was performed at 30° of knee flexion, initially without tibia translation and then using the force applied to the calf of 80 N. The femur and tibia were separately registered using rigid local SimpleITK landmark refinement; translation and rotation parameters were then calculated using the 3D transformation algorithms. The significance level was set at 0.05.

RESULTS

Initially, the device and method for obtaining the parameters of the 3D translation and rotation were validated. The pooled Standard Deviation for translation parameters was 0.81 mm and for rotation parameters 0.87°. Compared to the control group, statistically significant differences were found in parameters such as Anterior Shift [(median ± interquartile range) 3.89 mm ±6.55 vs. 0.90 mm ±2.78, P=0.002238] and External Rotation (-0.55° ±3.88 vs. -2.87° ±2.40, P=0.005074). Statistically significant correlations were observed in combined groups between Anterior Shift and parameters such as External Rotation (P=0.001611), PCL Tibial Attachment Point (pPCL) Anterior Shift (<0.000001), Rolimeter Measurement (P=0.000016), and Side-to-Side Difference (SSD) (P=0.000383). A significant statistical correlation was also observed between External Rotation and parameters such as Rolimeter (P=0.02261) and SSD (P=0.03458).

CONCLUSIONS

The analysis of the anterior tibia translation using stress MRI and the proposed three-dimensional calculation method allows for a detailed analysis of the tibial translation and rotation parameters. The correlations showed the importance of external rotation during anterior tibial translation.

Automatic tracking of healthy joint kinematics from stereo-radiography sequences

Kinematic tracking of healthy joints in radiography sequences is frequently performed by maximizing similarities between computed perspective projections of 3D computer models and corresponding objects' appearances in radiographic images. Significant human effort associated with manual tracking presents a major bottleneck in biomechanics research methods and limits the scale of target applications. The current work introduces a method for fully-automatic tracking of tibiofemoral and patellofemoral kinematics in stereo-radiography sequences for subjects performing dynamic activities. The proposed method involves the application of convolutional neural networks for annotating radiographs and a multi-stage optimization pipeline for estimating bone pose based on information provided by neural net predictions. Predicted kinematics are evaluated by comparing against manually-tracked trends across 20 distinct trials. Median absolute differences below 1.5 millimeters or degrees for 6 tibiofemoral and 3 patellofemoral degrees of freedom demonstrate the utility of our approach, which improves upon previous semi-automatic methods by enabling end-to-end automation. Implementation of a fully-automatic pipeline for kinematic tracking will benefit evaluation of human movement by enabling large-scale studies of healthy knee kinematics.

Transfer learning from an artificial radiograph-landmark dataset for registration of the anatomic skull model to dual fluoroscopic X-ray images

Registration of 3D anatomic structures to their 2D dual fluoroscopic X-ray images is a widely used motion tracking technique. However, deep learning implementation is often impeded by a paucity of medical images and ground truths. In this study, we proposed a transfer learning strategy for 3D-to-2D registration using deep neural networks trained from an artificial dataset. Digitally reconstructed radiographs (DRRs) and radiographic skull landmarks were automatically created from craniocervical CT data of a female subject. They were used to train a residual network (ResNet) for landmark detection and a cycle generative adversarial network (GAN) to eliminate the style difference between DRRs and actual X-rays. Landmarks on the X-rays experiencing GAN style translation were detected by the ResNet, and were used in triangulation optimization for 3D-to-2D registration of the skull in actual dual-fluoroscope images (with a non-orthogonal setup, point X-ray sources, image distortions, and partially captured skull regions). The registration accuracy was evaluated in multiple scenarios of craniocervical motions. In walking, learning-based registration for the skull had angular/position errors of 3.9 ± 2.1°/4.6 ± 2.2 mm. However, the accuracy was lower during functional neck activity, due to overly small skull regions imaged on the dual fluoroscopic images at end-range positions. The methodology to strategically augment artificial training data can tackle the complicated skull registration scenario, and has potentials to extend to widespread registration scenarios.

Ligament deformation patterns of the craniocervical junction during head axial rotation tracked by biplane fluoroscopes

BACKGROUND

Frequently, treatment decisions for craniocervical injuries and instability are based on imaging findings, but in vivo ligament kinematics were poorly understood. This study was to determine in vivo deformation patterns of primary ligaments in the craniocervical junction (i.e., C0-2), including the cruciform ligament, alar ligaments, and accessory ligaments, during dynamic head axial rotation.

METHODS

The skulls and cervical spines of eight asymptomatic female subjects were dynamically imaged using a biplane fluoroscopic imaging system, when they performed left and right head axial rotations. Using a 3D-to-2D registration technique, the in vivo positions and orientations of cervical segments were determined. An optimization algorithm was implemented to determine ligament wrapping paths, and the resulting ligament deformations were represented by percent elongations. Using paired t-tests, ligament deformations in the end-range position were compared to those in the neutral position.

FINDINGS

No significant differences were observed in segmental motions during left and right head rotations (p > 0.05). In general, slight deformations occurred in each component of the cruciform ligament. For the alar ligaments, the ipsilateral ligament was lengthened from -0.7 ± 13.8% to 16.6 ± 15.7% (p < 0.001*). For the accessory ligaments, the contralateral ligament was lengthened from -2.9 ± 7.5% to 10.1 ± 6.2% (p < 0.001*).

INTERPRETATION

This study reveals that there are distinct deformation patterns in craniocervical junction ligaments during dynamic axial head rotation. These ligament deformation data can enhance our understanding of the synergic function of craniocervical junction ligaments, and guide the treatment of craniocervical instability.

In vivo primary and coupled segmental motions of the healthy female head-neck complex during dynamic head axial rotation

While previous studies have greatly improved our knowledge on the motion capability of the cervical spine, few reported on the kinematics of the entire head-neck complex (C0-T1) during dynamic activities of the head in the upright posture. This study investigated in vivo kinematics of the entire head-neck complex (C0-T1) of eight female asymptomatic subjects during dynamic left-right head axial rotation using a dual fluoroscopic imaging system and 3D-to-2D registration techniques. During one-sided head rotation (i.e., left or right head rotation), the primary rotation of the overall head-neck complex (C0-T1) reached 55.5 ± 10.8°, the upper cervical spine region (C0-2) had a primary axial rotation of 39.7 ± 9.6° (71.3 ± 8.5% of the overall C0-T1 axial rotation), and the lower cervical spine region (C2-T1) had a primary rotation of 10.0 ± 3.7° (18.6 ± 7.2% of the overall C0-T1 axial rotation). Coupled bending rotations occurred in the upper and lower cervical spine regions in similar magnitude but opposite directions (upper: contralateral bending of 18.2 ± 5.9° versus lower: ipsilateral bending of 21.4 ± 5.1°), resulting in a compensatory cervical lateral curvature that balances the head to rotate horizontally. Furthermore, upper cervical segments (C0-1 or C1-2) provided main mobility in different rotational degrees of freedom needed for head axial rotations. Additionally, we quantitatively described both coupled segmental motions (flexion-extension and lateral bending) by correlation with the overall primary axial rotation of the head-neck complex. This investigation offers comprehensive baseline data regarding primary and coupled motions of craniocervical segments during head axial rotation.

Patellar Tendon Orientation and Strain Are Predictors of ACL Strain In Vivo During a Single-Leg Jump

BACKGROUND

There is little in vivo data that describe the relationships between patellar tendon orientation, patellar tendon strain, and anterior cruciate ligament (ACL) strain during dynamic activities. Quantifying how the quadriceps load the ACL via the patellar tendon is important for understanding ACL injury mechanisms.

HYPOTHESIS

We hypothesized that flexion angle, patellar tendon orientation, and patellar tendon strain influence ACL strain during a single-leg jump. Specifically, we hypothesized that patellar tendon and ACL strains would increase concurrently when the knee is positioned near extension during the jump.

STUDY DESIGN

Descriptive laboratory study.

METHODS

Models of the femur, tibia, ACL, patellar tendon, and quadriceps tendon attachment sites of 8 male participants were generated from magnetic resonance imaging (MRI). High-speed biplanar radiographs during a single-leg jump were obtained. The bone models were registered to the radiographs, thereby reproducing the in vivo positions of the bones, ligament, and tendon attachment sites. Flexion angle, patellar tendon orientation, patellar tendon strain, and ACL strain were measured from the registered models. ACL and patellar tendon strains were approximated by normalizing their length at each knee position to their length at the time of MRI. Two separate bivariate linear regression models were used to assess relationships between flexion angle and patellar tendon orientation and between ACL strain and patellar tendon strain. A multivariate linear regression model was used to assess whether flexion angle and patellar tendon strain were significant predictors of ACL strain during the inflight and landing portions of the jump.

RESULTS

Both flexion angle and patellar tendon strain were significant predictors (P < .05) of ACL strain. These results indicate that elevated ACL and patellar tendon strains were observed concurrently when the knee was positioned near extension.

CONCLUSION

Concurrent increases in patellar tendon and ACL strains indicate that the quadriceps load the ACL via the patellar tendon when the knee is positioned near extension.

CLINICAL RELEVANCE

Increased ACL strain when the knee is positioned near extension before landing may be due to quadriceps contraction. Thus, landing with unanticipated timing on an extended knee may increase vulnerability to ACL injury as a taut ligament is more likely to fail.

The effect of cognitive loading on the lower extremity movement coordination variability in patients with anterior cruciate ligament reconstruction

BACKGROUND

The altered coordination variability was related to anterior cruciate ligament (ACL) re-injury after ACL reconstruction (ACL-R). As motor performance is affected by the cognitive loading, understanding the interaction of cognition and coordination variability is crucial for addressing secondary injury prevention and restoring function in rehabilitation programs.

RESEARCH QUESTION

To investigate the lower extremity movement coordination variability asymmetry in individuals following ACL-R and determine the effects of cognitive loading on the coordination variability.

METHODS

Twenty-five males who received unilateral ACL-R using hamstring tendon autograft (7.4 ± 1.3 months past reconstruction). Participants performed walking without (single-task condition) and with the concurrent cognitive task (dual-task condition). The coordination variability in hip-knee coupled motion in different gait phases was calculated using vector coding technique.

RESULTS

The injured leg demonstrated greater coordination variability in hip flexion/knee flexion (HF/KF) during mid-stance phase (P = 0.012) than the uninjured leg in both conditions. No significant differences were observed in other phases of HF/KF variability or other measures in all phases between the injured and uninjured legs. Both legs increased the HF/KF coordination variability during loading response phase in dual-task condition than that in single-task condition (P < 0.001).

SIGNIFICANCE

Individuals following ACL-R demonstrated coordination variability asymmetry of sagittal plane hip-knee coupled motion. The dual cognitive task increased the coordination variability of hip flexion/knee flexion during loading response phase in individuals following ACL-R.

In Vivo Mechanical Function of the Distal Radial Ulnar Ligaments During Rotation of the Wrist

PURPOSE

The purpose of this study was to investigate changes in length of the volar and dorsal radioulnar ligaments (VRULs and DRULs), and the distal radioulnar joint (DRUJ) space during unweighted and weighted rotation of the wrist using magnetic resonance imaging and biplanar fluoroscopy.

METHODS

Fourteen wrists in 7 normal adult volunteers were imaged to define the 3-dimensional geometry of the DRUJ and the insertion sites of the superficial and deep bundles of the VRULs and DRULs. Subjects were imaged at 10 positions of forearm rotation ranging from full pronation to full supination, with or without a 5-pound weight. Lengths of the superficial and deep VRUL and DRUL bundles and DRUJ space were measured (in millimeters) at each position to evaluate ligament function and DRUJ stability.

RESULTS

In the unweighted and weighted trials, maximal elongation of both deep and superficial VRUL bundles occurred in supination and maximal lengths of the deep and superficial DRUL bundles occurred in pronation. Maximum DRUJ space occurred during pronation and a minimum occurred in 30° of supination. In weighted trials, there was a significant increase in deep and superficial VRUL bundle length at positions between 30° of pronation and 30° of supination; however, there was no effect of weight on DRULs length. In weighted trials, there was a significant increase in DRUJ space at positions between full pronation and 15° of supination.

CONCLUSIONS

This study demonstrates elongation of the VRULs in supination and the DRULs in pronation. There was no evidence of reciprocal loading of superficial/deep ligament bundles on either the dorsal or the volar aspects of the DRUJ. The effect of loading the wrist during rotation was apparent primarily in the VRULs, but not the DRULs. The DRUJ space was lowest at approximately 30° of supination.

CLINICAL RELEVANCE

These results add information to the literature regarding the complicated biomechanics of the triangular fibrocartilage complex and DRUJ. Future work should evaluate changes in biomechanics caused by triangular fibrocartilage complex tears to determine how tear severity and location relate to clinical symptoms.

Reconsidering Reciprocal Length Patterns of the Anteromedial and Posterolateral Bundles of the Anterior Cruciate Ligament During In Vivo Gait

BACKGROUND

Some cadaveric studies have indicated that the anterior cruciate ligament (ACL) consists of anteromedial and posterolateral bundles that display reciprocal function with regard to knee flexion. However, several in vivo imaging studies have suggested that these bundles elongate in parallel with regard to flexion. Furthermore, the most appropriate description of the functional anatomy of the ACL is still debated, with the ACL being described as consisting of 2 or 3 bundles or as a continuum of fibers.

HYPOTHESIS

As long as their origination and termination locations are defined within the ACL attachment site footprints, ACL bundles elongate in parallel with knee extension during gait.

STUDY DESIGN

Descriptive laboratory study.

METHODS

High-speed biplanar radiographs of the right knee joint were obtained during gait in 6 healthy male participants (mean ± SD: body mass index, 25.5 ± 1.2 kg/m²; age, 29.2 ± 3.8 years) with no history of lower extremity injury or surgery. Three-dimensional models of the right femur, tibia, and ACL attachment sites were created from magnetic resonance images. The bone models were registered to the biplanar radiographs, thereby reproducing the in vivo positions of the knee joint. For each knee position, the distances between the centroids of the ACL attachment sites were used to represent ACL length. The lengths of 1000 virtual bundles were measured for each participant by randomly sampling locations on the attachment site surfaces and measuring the distances between each pair of locations. Spearman rho rank correlations were performed between the virtual bundle lengths and ACL length.

RESULTS

The virtual bundle lengths were highly correlated with the length of the ACL, defined as the distance between the centroids of the attachment sites (rho = 0.91 ± 0.1, across participants; P < 5 × 10^-5). The lengths of the bundles that originated and terminated in the anterior and medial aspects of the ACL were positively correlated (rho = 0.81 ± 0.1; P < 5 × 10^-5) with the lengths of the bundles that originated and terminated in the posterior and lateral aspects of the ACL.

CONCLUSION

As long as their origination and termination points are specified within the footprint of the attachment sites, ACL bundles elongate in parallel as the knee is extended.

CLINICAL RELEVANCE

These data elucidate ACL functional anatomy and may help guide ACL reconstruction techniques.

In Vivo Anterior Cruciate Ligament Deformation During a Single-Legged Jump Measured by Magnetic Resonance Imaging and High-Speed Biplanar Radiography

BACKGROUND

The in vivo mechanics of the anterior cruciate ligament (ACL) and its bundles during dynamic activities are not completely understood. An improved understanding of how the ACL stabilizes the knee is likely to aid in the identification and prevention of injurious maneuvers.

PURPOSE/HYPOTHESIS

The purpose was to measure in vivo ACL strain during a single-legged jump through use of magnetic resonance imaging (MRI) and high-speed biplanar radiography. We hypothesized that ACL strain would increase with the knee near extension, and a peak in ACL strain would occur just before landing from the jump, potentially due to quadriceps contraction in anticipation of landing.

STUDY DESIGN

Descriptive laboratory study.

METHODS

Models of the femur, tibia, and ACL attachment sites of 8 male participants were generated from MRI scans through use of solid modeling. High-speed biplanar radiographs were obtained from these participants as they performed a single-legged jump. The bone models were registered to the biplanar radiographs, thereby reproducing the in vivo positions of the joint throughout the jump. ACL and bundle elongations were defined as the centroid to centroid distances between attachment sites for each knee position. ACL strain was defined as ACL length normalized to its length measured in the position of the knee at the time of MRI.

RESULTS

Peaks in ACL strain were observed before toe-off and 55 ± 35 milliseconds before initial ground contact. These peaks were associated with the knee positioned at low flexion angles. Mean ACL strain was inversely related to mean flexion angle (rho = -0.73, P < .001), such that ACL strain generally increased with knee extension throughout the jumping motion. ACL bundle lengths were significantly (rho > 0.85, P < .001) correlated with overall ACL length.

CONCLUSION

These findings provide insight into how landing in extension can increase the risk of ACL injury. Specifically, this study shows that peak ACL strain can occur just before landing from a single-legged jump. Thus, when an individual lands on an extended knee, the ACL is relatively taut, which may make it particularly vulnerable to injury, especially in the presence of a movement perturbation or unanticipated change in landing strategy.

CLINICAL RELEVANCE

This study provides a novel measurement of dynamic ACL strain during an athletic maneuver and lends insight into how landing in extension can increase the likelihood of ACL failure.

In vivo attachment site to attachment site length and strain of the ACL and its bundles during the full gait cycle measured by MRI and high-speed biplanar radiography

The purpose of this study was to measure in vivo attachment site to attachment site lengths and strains of the anterior cruciate ligament (ACL) and its bundles throughout a full cycle of treadmill gait. To obtain these measurements, models of the femur, tibia, and associated ACL attachment sites were created from magnetic resonance (MR) images in 10 healthy subjects. ACL attachment sites were subdivided into anteromedial (AM) and posterolateral (PL) bundles. High-speed biplanar radiographs were obtained as subjects ambulated at 1 m/s. The bone models were registered to the radiographs, thereby reproducing the in vivo positions of the bones and ACL attachment sites throughout gait. The lengths of the ACL and both bundles were estimated as straight line distances between attachment sites for each knee position. Increased attachment to attachment ACL length and strain were observed during midstance (length = 28.5 ± 2.6 mm, strain = 5 ± 4%, mean ± standard deviation), and heel strike (length = 30.5 ± 3.0 mm, strain = 12 ± 5%) when the knee was positioned at low flexion angles. Significant inverse correlations were observed between mean attachment to attachment ACL lengths and flexion (rho = -0.87, p < 0.001), as well as both bundle lengths and flexion (rho = -0.86, p < 0.001 and rho = -0.82, p < 0.001, respectively). AM and PL bundle attachment to attachment lengths were highly correlated throughout treadmill gait (rho = 0.90, p < 0.001). These data can provide valuable information to inform design criteria for ACL grafts used in reconstructive surgery, and may be useful in the design of rehabilitation and injury prevention protocols.

In vivo assessment of the interaction of patellar tendon tibial shaft angle and anterior cruciate ligament elongation during flexion

A potential cause of non-contact anterior cruciate ligament (ACL) injury is landing on an extended knee. In line with this hypothesis, studies have shown that the ACL is elongated with decreasing knee flexion angle. Furthermore, at low flexion angles the patellar tendon is oriented to increase the anterior shear component of force acting on the tibia. This indicates that knee extension represents a position in which the ACL is taut, and thus may have an increased propensity for injury, particularly in the presence of excessive force acting via the patellar tendon. However, there is very little in vivo data to describe how patellar tendon orientation and ACL elongation interact during flexion. Therefore, this study measured the patellar tendon tibial shaft angle (indicative of the relative magnitude of the shear component of force acting via the patellar tendon) and ACL length in vivo as subjects performed a quasi-static lunge at varying knee flexion angles. Spearman rho rank correlations within each individual revealed that flexion angles were inversely correlated to both ACL length (rho = -0.94 ± 0.07, mean ± standard deviation, p < 0.05) and patellar tendon tibial shaft angle (rho = -0.99 ± 0.01, p < 0.05). These findings indicate that when the knee is extended, the ACL is both elongated and the patellar tendon tibial shaft angle is increased, resulting in a relative increase in anterior shear force on the tibia acting via the patellar tendon. Therefore, these data support the hypothesis that landing with the knee in extension is a high risk scenario for ACL injury.

A Comparison of Knee Abduction Angles Measured by a 3D Anatomic Coordinate System Versus Videographic Analysis: Implications for Anterior Cruciate Ligament Injury

Background

Knee positions involved in noncontact anterior cruciate ligament (ACL) injury have been studied via analysis of injury videos. Positions of high ACL strain have been identified in vivo. These methods have supported different hypotheses regarding the role of knee abduction in ACL injury.

Purpose/Hypothesis

The purpose of this study was to compare knee abduction angles measured by 2 methods: using a 3-dimensional (3D) coordinate system based on anatomic features of the bones versus simulated 2-dimensional (2D) videographic analysis. We hypothesized that knee abduction angles measured in a 2D videographic analysis would differ from those measured from 3D bone anatomic features and that videographic knee abduction angles would depend on flexion angle and on the position of the camera relative to the patient.

Study Design

Descriptive laboratory study.

Methods

Models of the femur and tibia were created from magnetic resonance images of 8 healthy male participants. The models were positioned to match biplanar fluoroscopic images obtained as participants posed in lunges of varying flexion angles (FLAs). Knee abduction angle was calculated from the positioned models in 2 ways: (1) varus-valgus angle (VVA), defined as the angle between the long axis of the tibia and the femoral transepicondylar axis by use of a 3D anatomic coordinate system; and (2) coronal plane angle (CPA), defined as the angle between the long axis of the tibia and the long axis of the femur projected onto the tibial coronal plane to simulate a 2D videographic analysis. We then simulated how changing the position of the camera relative to the participant would affect knee abduction angles.

Results

During flexion, when CPA was calculated from a purely anterior or posterior view of the joint-an ideal scenario for measuring knee abduction from 2D videographic analysis-CPA was significantly different from VVA (P < .0001). CPA also varied substantially with the position of the camera relative to the participant.

Conclusion

How closely CPA (derived from 2D videographic analysis) relates to VVA (derived from a 3D anatomic coordinate system) depends on FLA and camera orientation.

Clinical Relevance

This study provides a novel comparison of knee abduction angles measured from 2D videographic analysis and those measured within a 3D anatomic coordinate system. Consideration of these findings is important when interpreting 2D videographic data regarding knee abduction angle in ACL injury.