Knee cartilage MRI with in situ mechanical loading using prospective motion correction

Time-Resolved Quantification of Patellofemoral Cartilage Deformation in Response to Loading and Unloading via Dynamic MRI With Prospective Motion Correction

BACKGROUND

In vivo cartilage deformation has been studied by static magnetic resonance imaging (MRI) with in situ loading, but knowledge about strain dynamics after load onset and release is scarce.

PURPOSE

To measure the dynamics of patellofemoral cartilage deformation and recovery in response to in situ loading and unloading by using MRI with prospective motion correction.

STUDY TYPE

Prospective.

SUBJECTS

Ten healthy male volunteers (age: [31.4 ± 3.2] years).

FIELD STRENGTH/SEQUENCE

T1-weighted RF-spoiled 2D gradient-echo sequence with a golden angle radial acquisition scheme, augmented with prospective motion correction, at 3 T.

ASSESSMENT

In situ knee loading was realized with a flexion angle of approximately 40° using an MR-compatible pneumatic loading device. The loading paradigm consisted of 2 minutes of unloaded baseline followed by a 5-minute loading bout with 50% body weight and an unloading period of 38 minutes. The cartilage strain was assessed as the mean distance between patellar and femoral bone-cartilage interfaces as a percentage of the initial (pre-load) distance.

STATISTICAL TESTS

Wilcoxon signed-rank tests (significance level: P < 0.05), Pearson correlation coefficient (r).

RESULTS

The cartilage compression and recovery behavior was characterized by a viscoelastic response. The elastic compression ([-12.5 ± 3.1]%) was significantly larger than the viscous compression ([-7.6 ± 1.5]%) and the elastic recovery ([10.5 ± 2.1]%) was significantly larger than the viscous recovery ([6.1 ± 1.8]%). There was a significant residual offset strain ([-3.6 ± 2.3]%) across the cohort. A significant negative correlation between elastic compression and elastic recovery was observed (r = -0.75).

DATA CONCLUSION

The in vivo cartilage compression and recovery time course in response to loading was successfully measured via dynamic MRI with prospective motion correction. The clinical relevance of the strain characteristics needs to be assessed in larger subject and patient cohorts.

LEVEL OF EVIDENCE

2 TECHNICAL EFFICACY: Stage 1.

Articular Cartilage-From Basic Science Structural Imaging to Non-Invasive Clinical Quantitative Molecular Functional Information for AI Classification and Prediction

This review presents the changes that the imaging of articular cartilage has undergone throughout the last decades. It highlights that the expectation is no longer to image the structure and associated functions of articular cartilage but, instead, to devise methods for generating non-invasive, function-depicting images with quantitative information that is useful for detecting the early, pre-clinical stage of diseases such as primary or post-traumatic osteoarthritis (OA/PTOA). In this context, this review summarizes (a) the structure and function of articular cartilage as a molecular imaging target, (b) quantitative MRI for non-invasive assessment of articular cartilage composition, microstructure, and function with the current state of medical diagnostic imaging, (c), non-destructive imaging methods, (c) non-destructive quantitative articular cartilage live-imaging methods, (d) artificial intelligence (AI) classification of degeneration and prediction of OA progression, and (e) our contribution to this field, which is an AI-supported, non-destructive quantitative optical biopsy for early disease detection that operates on a digital tissue architectural fingerprint. Collectively, this review shows that articular cartilage imaging has undergone profound changes in the purpose and expectations for which cartilage imaging is used; the image is becoming an AI-usable biomarker with non-invasive quantitative functional information. This may aid in the development of translational diagnostic applications and preventive or early therapeutic interventions that are yet beyond our reach.

[Influence of the quadriceps muscles on the patellofemoral contact in patients with low flexion patellofemoral instability after MPFL reconstruction]

INTRODUCTION

MPFL reconstruction represents one of the most important surgical treatment options for recurrent patellar dislocations at low flexion angles associated with low flexion patellofemoral instability. Nevertheless, the role of quadriceps muscles in patients with patellofemoral instability before and after patellofemoral stabilization using MPFL reconstruction has not been fully elucidated. The present study investigates the influence of quadriceps muscles on the patellofemoral contact in patients with low flexion patellofemoral instability (PFI) before and after surgical patellofemoral stabilization using MPFL reconstruction using 3 T MRI datasets in early degrees of flexion (0-30°).

METHODS

In this prospective cohort study, 15 patients with low flexion PFI before and after MPFL reconstruction and 15 subjects with healthy knee joints were studied using dynamic MRI scans. MRI scans were performed in a custom-made pneumatic knee loading device to determine the patellofemoral cartilage contact area (CCA) with and without quadriceps activation (50 N). Comparative measurements were performed using 3D cartilage and bone meshes in 0-30° knee flexion in the patients with patellofemoral instability preoperatively and postoperatively.

RESULTS

The preoperative patellofemoral CCA of patients with low flexion PFI was 67.3 ± 47.3 mm2 in 0° flexion, 118.9 ± 56.6 mm2 in 15° flexion, and 267.6 ± 96.1 mm2 in 30° flexion. With activated quadriceps muscles (50 N), the contact area was 72.4 ± 45.9 mm2 in extension, 112.5 ± 54.9 mm2 in 15° flexion, and 286.1 ± 92.7 mm2 in 30° flexion without statistical significance. Postoperatively determined CCA revealed 159.3 ± 51.4 mm2 , 189.6 ± 62.2 mm2 and 347.3 ± 52.1 mm2 in 0°, 15° and 30° flexion. Quadriceps activation with 50 N showed a contact area in extension of 141.0 ± 63.8 mm2, 206.6 ± 67.7 mm2 in 15° flexion, and 353.5 ± 64.6 mm2 in 30° flexion, also without statistical difference compared with unloaded CCAs. Subjects with healthy knee joints showed an increase of 10.3% in CCA at 30° of flexion (p = 0.003).

CONCLUSION

Although patellofemoral CCA increases significantly after isolated MPFL reconstruction in patients with low flexion patellofemoral instability, there is no significant influence of quadriceps muscles either preoperatively or postoperatively.

Influence of Medial Patellofemoral Ligament Reconstruction on Patellofemoral Contact in Patients With Low-Flexion Patellar Instability: An MRI Study

Background

Medial patellofemoral ligament (MPFL) reconstruction is a well-established procedure for the treatment of patients with patellofemoral instability (PFI) at low flexion angles (0°-30°). Little is known about the effect of MPFL surgery on patellofemoral cartilage contact area (CCA) during the first 30° of knee flexion.

Purpose/Hypothesis

The purpose of this study was to investigate the effect of MPFL reconstruction on CCA using magnetic resonance imaging (MRI). We hypothesized that patients with PFI would have a lower CCA than patients with healthy knees and that CCA would increase after MPFL reconstruction over the course of low knee flexion.

Study Design

Cohort study; Level of evidence, 2.

Methods

In a prospective matched-paired cohort study, the CCA of 13 patients with low-flexion PFI was determined before and after MPFL reconstruction, and the data were compared with those of 13 healthy volunteers (controls). MRI was performed with the knee at 0°, 15°, and 30° of flexion in a custom-designed knee-positioning device. To suppress motion artifacts, motion correction was performed using a Moiré Phase Tracking system via a tracking marker attached to the patella. The CCA was calculated on the basis of semiautomatic cartilage and bone segmentation and registration.

Results

The CCA (mean ± SD) at 0°, 15°, and 30° of flexion for the control participants was 1.38 ± 0.62, 1.91 ± 0.98, and 3.68 ± 0.92 cm², respectively. In patients with PFI, the CCA at 0°, 15°, and 30° of flexion was 0.77 ± 0.49, 1.26 ± 0.60, and 2.89 ± 0.89 cm² preoperatively and 1.65 ± 0.55, 1.97 ± 0.68, and 3.52 ± 0.57 cm² postoperatively. Patients with PFI exhibited a significantly reduced preoperative CCA at all 3 flexion angles when compared with controls (P ≤ .045 for all). Postoperatively, there was a significant increase in CCA at 0° of flexion (P = .001), 15° of flexion (P = .019) and 30° of flexion (P = .026). There were no significant postoperative differences in CCA between patients with PFI and controls at any flexion angle.

Conclusion

Patients with low-flexion patellar instability showed a significant reduction in patellofemoral CCA at 0°, 15°, and 30° of flexion. MPFL reconstruction increased the contact area significantly at all angles.

Validation of a Finite Element Simulation for Predicting Individual Knee Joint Kinematics

Goal: We introduce an in-vivo validated finite element (FE) simulation approach for predicting individual knee joint kinematics. Our vision is to improve clinicians' understanding of the complex individual anatomy and potential pathologies to improve treatment and restore physiological joint kinematics. Methods: Our 3D FE modeling approach for individual human knee joints is based on segmentation of anatomical structures extracted from routine static magnetic resonance (MR) images. We validate the predictive abilities of our model using static MR images of the knees of eleven healthy volunteers in dedicated knee poses, which are achieved using a customized MR-compatible pneumatic loading device. Results: Our FE simulations reach an average translational accuracy of 2 mm and an average angular accuracy of 1[Formula: see text] compared to the reference knee pose. Conclusions: Reaching high accuracy, our individual FE model can be used in the decision-making process to restore knee joint stability and functionality after various knee injuries.

Biomechanical Effects of Chronic Ankle Instability on the Talar Cartilage Matrix: The Value of T1ρ Relaxation Mapping Without and With Mechanical Loading

BACKGROUND

T1ρ mapping has been proposed for the detection of early cartilage degeneration associated with chronic ankle instability (CAI). However, there are limited data surrounding the influence of ankle loading on T1ρ relaxation.

PURPOSE

To evaluate T1ρ relaxation times of talar cartilage, as an indicator of early degenerative changes, associated with CAI and to investigate the influence of acute axial in situ loading on T1ρ values in CAI patients and healthy controls.

STUDY TYPE

Prospective.

SUBJECTS

A total of 9 patients (age = 21.8 ± 2.5 years, male/female = 2/7) with chronic ankle instability and 18 healthy control subjects (age = 22.8 ± 3.6 years, male/female = 5/13).

FIELD STRENGTH

3 T.

SEQUENCE

3D gradient echo fast low-angle shot (FLASH) sequence augmented with a variable spin-lock preparation period.

ASSESSMENT

Ankle T1ρ mapping was performed without and with axial loading of 500 N. The talar cartilage was segmented in five coronal slices covering the central talocrural joint. Median talar T1ρ values were separately calculated for the medial and lateral facets.

STATISTICAL TESTS

Mann-Whitney U and Wilcoxon signed-rank tests, significance level: P < 0.05.

RESULTS

For the combined cohorts, the statistical analysis yielded significantly lower T1ρ values with loading compared to the no-load measurement for both the lateral (no load: [51.0 ± 4.0] msec, load: [49.5 ± 5.4] msec) as well as the medial compartment (no load: [50.0 ± 5.4] msec, load: [47.8 ± 6.8] msec). In the unloaded scans, the CAI patients showed significantly increased talar T1ρ values ([53.0 ± 7.4] mse ) compared to the healthy control subjects ([48.8 ± 4.1] msec) in the medial compartment.

DATA CONCLUSION

Increased talar T1ρ relaxation times in CAI patients compared to healthy controls suggest that T1ρ relaxation is a sensitive biomarker for CAI-induced early-stage cartilage degeneration. However, the load-induced T1ρ change did not prove to be a viable marker for the altered biomechanical properties of the hyaline talar cartilage.

LEVEL OF EVIDENCE

2 LEVEL OF EFFICACY: Stage 2.

Individual Influence of Trochlear Dysplasia on Patellofemoral Kinematics after Isolated MPFL Reconstruction

INTRODUCTION

The influence of the MPFL graft in cases of patella instability with dysplastic trochlea is a controversial topic. The effect of the MPFL reconstruction as single therapy is under investigation, especially with severely dysplastic trochlea (Dejour types C and D). The purpose of this study was to evaluate the impact of trochlear dysplasia on patellar kinematics in patients suffering from low flexion patellar instability under weight-bearing conditions after isolated MPFL reconstruction.

MATERIAL AND METHODS

Thirteen patients were included in this study, among them were eight patients with mild dysplasia (Dejour type A and B) and five patients with severe dysplasia (Dejour type C and D). By performing a knee MRI with in situ loading, patella kinematics and the patellofemoral cartilage contact area could be measured under the activation of the quadriceps musculature in knee flexion angles of 0°, 15° and 30°. To mitigate MRI motion artefacts, prospective motion correction based on optical tracking was applied. Bone and cartilage segmentation were performed semi-automatically for further data analysis. Cartilage contact area (CCA) and patella tilt were the main outcome measures for this study. Pre- and post-surgery measures were compared for each group.

RESULTS

Data visualized a trending lower patella tilt after MPFL graft installation in both groups and flexion angles of the knee. There were no significant changes in patella tilt at 0° (unloaded pre-surgery: 22.6 ± 15.2; post-surgery: 17.7 ± 14.3; p = 0.110) and unloaded 15° flexion (pre-surgery: 18.9 ± 12.7; post-surgery: 12.2 ± 13.0; p = 0.052) of the knee in patients with mild dysplasia, whereas in patients with severe dysplasia of the trochlea the results happened not to be significant in the same angles with loading of 5 kg (0° flexion pre-surgery: 34.4 ± 12.1; post-surgery: 31.2 ± 16.1; p = 0.5; 15° flexion pre-surgery: 33.3 ± 6.1; post-surgery: 23.4 ± 8.6; p = 0.068). CCA increased in every flexion angle and group, but significant increase was seen only between 0°-15° (unloaded and loaded) in mild dysplasia of the trochlea, where significant increase in Dejour type C and D group was seen with unloaded full extension of the knee (0° flexion) and 30° flexion (unloaded and loaded).

CONCLUSION

This study proves a significant effect of the MPFL graft to cartilage contact area, as well as an improvement of the patella tilt in patients with mild dysplasia of the trochlea. Thus, the MPFL can be used as a single treatment for patient with Dejour type A and B dysplasia. However, in patients with severe dysplasia the MPFL graft alone does not significantly increase CCA.

Kartogenin (KGN)/synthetic melanin nanoparticles (SMNP) loaded theranostic hydrogel scaffold system for multiparametric magnetic resonance imaging guided cartilage regeneration

Cartilage regeneration after injury is still a great challenge in clinics, which suffers from its avascularity and poor proliferative ability. Herein we designed a novel biocompatible cellulose nanocrystal/GelMA (gelatin-methacrylate anhydride)/HAMA (hyaluronic acid-methacrylate anhydride)-blended hydrogel scaffold system, loaded with synthetic melanin nanoparticles (SMNP) and a bioactive drug kartogenin (KGN) for theranostic purpose. We found that the SMNP-KGN/Gel showed favorable mechanical property, thermal stability, and distinct magnetic resonance imaging (MRI) contrast enhancement. Meanwhile, the sustained release of KGN could recruit bone-derived mesenchymal stem cells to proliferate and differentiate into chondrocytes, which promoted cartilage regeneration in vitro and in vivo. The hydrogel degradation and cartilage restoration were simultaneously monitored by multiparametric MRI for 12 weeks, and further confirmed by histological analysis. Together, these results validated the multifunctional hydrogel as a promising tissue engineering platform for noninvasive imaging-guided precision therapy in cartilage regenerative medicine.

Device for Assessing Knee Joint Dynamics During Magnetic Resonance Imaging

BACKGROUND

Knee assessment with and without load using magnetic resonance imaging (MRI) can provide information on knee joint dynamics and improve the diagnosis of knee joint diseases. Performing such studies on a routine MRI-scanner require a load-exerting device during scanning. There is a need for more studies on developing loading devices and evaluating their clinical potential.

PURPOSE

Design and develop a portable and easy-to-use axial loading device to evaluate the knee joint dynamics during the MRI study.

STUDY TYPE

Prospective study.

SUBJECTS

Nine healthy subjects.

FIELD STRENGTH/SEQUENCE

A 0.25 T standing-open MRI and 3.0 T MRI. PD-T₂ -weighted FSE, 3D-fast-spoiled-gradient-echo, FS-PD, and CartiGram sequences.

ASSESSMENT

Design and development of loading device, calibration of loads, MR safety assessment (using projectile angular displacement, torque, and temperature tests). Scoring system for ease of doing. Qualitative (by radiologist) and quantitative (using structural similarity index measure [SSIM]) image-artifact assessment. Evaluation of repeatability, comparison with various standing stances load, and loading effect on knee MR parameters (tibiofemoral bone gap [TFBG], femoral cartilage thickness [FCT], tibial cartilage thickness [TCT], femoral cartilage T₂ -value [FCT2], and tibia cartilage T₂ -value [TCT2]). The relative percentage change (RPC) in parameters due to the device load was computed.

STATISTICAL TEST

Pearson's correlation coefficient (r).

RESULTS

The developed device is conditional-MR safe (details in the manuscript and supplementary materials), 15 × 15 × 45 cm³ dimension, and <3 kg. The ease of using the device was 4.9/5. The device introduced no visible image artifacts, and SSIM of 0.9889 ± 0.0153 was observed. The TFBG intraobserver variability (absolute difference) was <0.1 mm. Interobserver variability of all regions of interest was <0.1 mm. The load exerted by the device was close to the load during standing on both legs in 0.25 T scanner with r > 0.9. Loading resulted in RPC of 1.5%-11.0%, 7.9%-8.5%, and -1.5% to 13.0% in the TFBG, FCT, and TCT, respectively. FCT2 and TCT2 were reduced in range of 1.5-2.7 msec and 0.5-2.3 msec due to load.

DATA CONCLUSION

The proposed device is conditionally MR safe, low cost (material cost < INR 6000), portable, and effective in loading the knee joint with up to 50% of body weight.

EVIDENCE LEVEL

1 TECHNICAL EFFICACY: Stage 1.

Retrospective Camera-Based Respiratory Gating in Clinical Whole-Heart 4D Flow MRI

Background

Respiratory gating is generally recommended in 4D flow MRI of the heart to avoid blurring and motion artifacts. Recently, a novel automated contact‐less camera‐based respiratory motion sensor has been introduced.

Purpose

To compare camera‐based respiratory gating (CAM) with liver‐lung‐navigator‐based gating (NAV) and no gating (NO) for whole‐heart 4D flow MRI.

Study Type

Retrospective.

Subjects

Thirty two patients with a spectrum of cardiovascular diseases.

Field Strength/Sequence

A 3T, 3D‐cine spoiled‐gradient‐echo‐T1‐weighted‐sequence with flow‐encoding in three spatial directions.

Assessment

Respiratory phases were derived and compared against each other by cross‐correlation. Three radiologists/cardiologist scored images reconstructed with camera‐based, navigator‐based, and no respiratory gating with a 4‐point Likert scale (qualitative analysis). Quantitative image quality analysis, in form of signal‐to‐noise ratio (SNR) and liver‐lung‐edge (LLE) for sharpness and quantitative flow analysis of the valves were performed semi‐automatically.

Statistical Tests

One‐way repeated measured analysis of variance (ANOVA) with Wilks's lambda testing and follow‐up pairwise comparisons. Significance level of P ≤ 0.05. Krippendorff's‐alpha‐test for inter‐rater reliability.

Results

The respiratory signal analysis revealed that CAM and NAV phases were highly correlated (C = 0.93 ± 0.09, P < 0.01). Image scoring showed poor inter‐rater reliability and no significant differences were observed (P ≥ 0.16). The image quality comparison showed that NAV and CAM were superior to NO with higher SNR (P = 0.02) and smaller LLE (P < 0.01). The quantitative flow analysis showed significant differences between the three respiratory‐gated reconstructions in the tricuspid and pulmonary valves (P ≤ 0.05), but not in the mitral and aortic valves (P > 0.05). Pairwise comparisons showed that reconstructions without respiratory gating were different in flow measurements to either CAM or NAV or both, but no differences were found between CAM and NAV reconstructions.

Data Conclusion

Camera‐based respiratory gating performed as well as conventional liver‐lung‐navigator‐based respiratory gating. Quantitative image quality analysis showed that both techniques were equivalent and superior to no‐gating‐reconstructions. Quantitative flow analysis revealed local flow differences (tricuspid/pulmonary valves) in images of no‐gating‐reconstructions, but no differences were found between images reconstructed with camera‐based and navigator‐based respiratory gating.

Level of Evidence

3

Technical Efficacy

Stage 2

Quantification of patellofemoral cartilage deformation and contact area changes in response to static loading via high-resolution MRI with prospective motion correction

BACKGROUND

Higher-resolution MRI of the patellofemoral cartilage under loading is hampered by subject motion since knee flexion is required during the scan.

PURPOSE

To demonstrate robust quantification of cartilage compression and contact area changes in response to in situ loading by means of MRI with prospective motion correction and regularized image postprocessing.

STUDY TYPE

Cohort study.

SUBJECTS

Fifteen healthy male subjects.

FIELD STRENGTH

3 T.

SEQUENCE

Spoiled 3D gradient-echo sequence augmented with prospective motion correction based on optical tracking. Measurements were performed with three different loads (0/200/400 N).

ASSESSMENT

Bone and cartilage segmentation was performed manually and regularized with a deep-learning approach. Average patellar and femoral cartilage thickness and contact area were calculated for the three loading situations. Reproducibility was assessed via repeated measurements in one subject.

STATISTICAL TESTS

Comparison of the three loading situations was performed by Wilcoxon signed-rank tests.

RESULTS

Regularization using a deep convolutional neural network reduced the variance of the quantified relative load-induced changes of cartilage thickness and contact area compared to purely manual segmentation (average reduction of standard deviation by ∼50%) and repeated measurements performed on the same subject demonstrated high reproducibility of the method. For the three loading situations (0/200/400 N), the patellofemoral cartilage contact area as well as the mean patellar and femoral cartilage thickness were significantly different from each other (P < 0.05). While the patellofemoral cartilage contact area increased under loading (by 14.5/19.0% for loads of 200/400 N), patellar and femoral cartilage thickness exhibited a load-dependent thickness decrease (patella: -4.4/-7.4%, femur: -3.4/-7.1% for loads of 200/400 N).

DATA CONCLUSION

MRI with prospective motion correction enables quantitative evaluation of patellofemoral cartilage deformation and contact area changes in response to in situ loading. Regularizing the manual segmentations using a neural network enables robust quantification of the load-induced changes.

LEVEL OF EVIDENCE

2 Technical Efficacy: Stage 2 J. Magn. Reson. Imaging 2019;50:1561-1570.

Optical tracking with two markers for robust prospective motion correction for brain imaging

OBJECTIVE

Prospective motion correction (PMC) during brain imaging using camera-based tracking of a skin-attached marker may suffer from problems including loss of marker visibility due to the coil and false correction due to non-rigid-body facial motion, such as frowning or squinting. A modified PMC system is introduced to mitigate these problems and increase the robustness of motion correction.

MATERIALS AND METHODS

The method relies on simultaneously tracking two markers, each providing six degrees of freedom, that are placed on the forehead. This allows us to track head motion when one marker is obscured and detect skin movements to prevent false corrections. Experiments were performed to compare the performance of the two-marker motion correction technique to the previous single-marker approach.

RESULTS

Experiments validate the theory developed for adaptive marker tracking and skin movement detection, and demonstrate improved image quality during obstruction of the line-of-sight of one marker when subjects squint or when subjects squint and move simultaneously.

CONCLUSION

The proposed methods eliminate two common failure modes of PMC and substantially improve the robustness of PMC, and they can be applied to other optical tracking systems capable of tracking multiple markers. The methods presented can be adapted to the use of more than two markers.