Machine learning-augmented and microspectroscopy-informed multiparametric MRI for the non-invasive prediction of articular cartilage composition

Gaining Insight into Updated MR Imaging for Quantitative Assessment of Cartilage Injury in Knee Osteoarthritis

PURPOSE OF THE REVIEW

Knee Osteoarthritis (KOA) entails progressive cartilage degradation, reviewed via MRI for morphology, biochemical composition, and microtissue alterations, discussing clinical advantages, limitations, and research applicability.

RECENT FINDINGS

Compositional MRI, like T2/T2* mapping, T1rho mapping, gagCEST, dGEMRIC, sodium imaging, diffusion-weighted imaging, and diffusion-tensor imaging, provide insights into cartilage injury in KOA. These methods quantitatively measure collagen, glycosaminoglycans, and water content, revealing important information about biochemical compositional and microstructural alterations. Innovative techniques like hybrid multi-dimensional MRI and diffusion-relaxation correlation spectrum imaging show potential in depicting initial cartilage changes at a sub-voxel level. Integration of automated image analysis tools addressed limitations in manual cartilage segmentation, ensuring robust and reproducible assessments of KOA cartilage. Compositional MRI techniques reveal microstructural changes in cartilage. Multi-dimensional MR imaging assesses biochemical alterations in KOA-afflicted cartilage, aiding early degeneration identification. Integrating artificial intelligence enhances cartilage analysis, optimal diagnostic accuracy for early KOA detection and monitoring.

Quantitative MRI methods for the assessment of structure, composition, and function of musculoskeletal tissues in basic research and preclinical applications

Osteoarthritis (OA) is a disabling chronic disease involving the gradual degradation of joint structures causing pain and dysfunction. Magnetic resonance imaging (MRI) has been widely used as a non-invasive tool for assessing OA-related changes. While anatomical MRI is limited to the morphological assessment of the joint structures, quantitative MRI (qMRI) allows for the measurement of biophysical properties of the tissues at the molecular level. Quantitative MRI techniques have been employed to characterize tissues' structural integrity, biochemical content, and mechanical properties. Their applications extend to studying degenerative alterations, early OA detection, and evaluating therapeutic intervention. This article is a review of qMRI techniques for musculoskeletal tissue evaluation, with a particular emphasis on articular cartilage. The goal is to describe the underlying mechanism and primary limitations of the qMRI parameters, their association with the tissue physiological properties and their potential in detecting tissue degeneration leading to the development of OA with a primary focus on basic and preclinical research studies. Additionally, the review highlights some clinical applications of qMRI, discussing the role of texture-based radiomics and machine learning in advancing OA research.

Updates on Compositional MRI Mapping of the Cartilage: Emerging Techniques and Applications

Osteoarthritis (OA) is a widely occurring degenerative joint disease that is severely debilitating and causes significant socioeconomic burdens to society. Magnetic resonance imaging (MRI) is the preferred imaging modality for the morphological evaluation of cartilage due to its excellent soft tissue contrast and high spatial resolution. However, its utilization typically involves subjective qualitative assessment of cartilage. Compositional MRI, which refers to the quantitative characterization of cartilage using a variety of MRI methods, can provide important information regarding underlying compositional and ultrastructural changes that occur during early OA. Cartilage compositional MRI could serve as early imaging biomarkers for the objective evaluation of cartilage and help drive diagnostics, disease characterization, and response to novel therapies. This review will summarize current and ongoing state-of-the-art cartilage compositional MRI techniques and highlight emerging methods for cartilage compositional MRI including MR fingerprinting, compressed sensing, multiexponential relaxometry, improved and robust radio-frequency pulse sequences, and deep learning-based acquisition, reconstruction, and segmentation. The review will also briefly discuss the current challenges and future directions for adopting these emerging cartilage compositional MRI techniques for use in clinical practice and translational OA research studies. EVIDENCE LEVEL: 2 TECHNICAL EFFICACY: Stage 2.

Machine Learning Prediction of Collagen Fiber Orientation and Proteoglycan Content From Multiparametric Quantitative MRI in Articular Cartilage

BACKGROUND

Machine learning models trained with multiparametric quantitative MRIs (qMRIs) have the potential to provide valuable information about the structural composition of articular cartilage.

PURPOSE

To study the performance and feasibility of machine learning models combined with qMRIs for noninvasive assessment of collagen fiber orientation and proteoglycan content.

STUDY TYPE

Retrospective, animal model.

ANIMAL MODEL

An open-source single slice MRI dataset obtained from 20 samples of 10 Shetland ponies (seven with surgically induced cartilage lesions followed by treatment and three healthy controls) yielded to 1600 data points, including 10% for test and 90% for train validation.

FIELD STRENGTH/SEQUENCE

A 9.4 T MRI scanner/qMRI sequences: T₁ , T₂ , adiabatic T_1ρ and T_2ρ , continuous-wave T_1ρ and relaxation along a fictitious field (T_RAFF ) maps.

ASSESSMENT

Five machine learning regression models were developed: random forest (RF), support vector regression (SVR), gradient boosting (GB), multilayer perceptron (MLP), and Gaussian process regression (GPR). A nested cross-validation was used for performance evaluation. For reference, proteoglycan content and collagen fiber orientation were determined by quantitative histology from digital densitometry (DD) and polarized light microscopy (PLM), respectively.

STATISTICAL TESTS

Normality was tested using Shapiro-Wilk test, and association between predicted and measured values was evaluated using Spearman's Rho test. A P-value of 0.05 was considered as the limit of statistical significance.

RESULTS

Four out of the five models (RF, GB, MLP, and GPR) yielded high accuracy (R²  = 0.68-0.75 for PLM and 0.62-0.66 for DD), and strong significant correlations between the reference measurements and predicted cartilage matrix properties (Spearman's Rho = 0.72-0.88 for PLM and 0.61-0.83 for DD). GPR algorithm had the highest accuracy (R²  = 0.75 and 0.66) and lowest prediction-error (root mean squared [RMSE] = 1.34 and 2.55) for PLM and DD, respectively.

DATA CONCLUSION

Multiparametric qMRIs in combination with regression models can determine cartilage compositional and structural features, with higher accuracy for collagen fiber orientation than proteoglycan content.

EVIDENCE LEVEL

2 TECHNICAL EFFICACY: Stage 2.

Towards Automatic Cartilage Quantification in Clinical Trials - Continuing from the 2019 IWOAI Knee Segmentation Challenge

Objective

To evaluate whether the deep learning (DL) segmentation methods from the six teams that participated in the IWOAI 2019 Knee Cartilage Segmentation Challenge are appropriate for quantifying cartilage loss in longitudinal clinical trials.

Design

We included 556 subjects from the Osteoarthritis Initiative study with manually read cartilage volume scores for the baseline and 1-year visits. The teams used their methods originally trained for the IWOAI 2019 challenge to segment the 1130 knee MRIs. These scans were anonymized and the teams were blinded to any subject or visit identifiers. Two teams also submitted updated methods. The resulting 9,040 segmentations are available online.The segmentations included tibial, femoral, and patellar compartments. In post-processing, we extracted medial and lateral tibial compartments and geometrically defined central medial and lateral femoral sub-compartments. The primary study outcome was the sensitivity to measure cartilage loss as defined by the standardized response mean (SRM).

Results

For the tibial compartments, several of the DL segmentation methods had SRMs similar to the gold standard manual method. The highest DL SRM was for the lateral tibial compartment at 0.38 (the gold standard had 0.34). For the femoral compartments, the gold standard had higher SRMs than the automatic methods at 0.31/0.30 for medial/lateral compartments.

Conclusion

The lower SRMs for the DL methods in the femoral compartments at 0.2 were possibly due to the simple sub-compartment extraction done during post-processing. The study demonstrated that state-of-the-art DL segmentation methods may be used in standardized longitudinal single-scanner clinical trials for well-defined cartilage compartments.

Contribution of collagen degradation and proteoglycan depletion to cartilage degeneration in primary and secondary osteoarthritis: an in silico study

OBJECTIVES

Current experimental approaches cannot elucidate the effect of maladaptive changes on the main cartilage constituents during the degeneration process in osteoarthritis (OA). In silico approaches, however, allow creating 'virtual knock-out' cases to elucidate these effects in a constituent-specific manner. We used such an approach to study the main mechanisms of cartilage degeneration in different mechanical loadings associated with the following OA etiologies: (1) physiological loading of degenerated cartilage, (2) injurious loading of healthy intact cartilage and (3) physiological loading of cartilage with a focal defect.

METHODS

We used the recently developed Cartilage Adaptive REorientation Degeneration (CARED) framework to simulate cartilage degeneration associated with primary and secondary OA (OA cases (1)-(3)). CARED incorporates numerical description of tissue-level cartilage degeneration mechanisms in OA, namely, collagen degradation, collagen reorientation, fixed charged density loss and tissue hydration increase following mechanical loading. We created 'virtual knock-out' scenarios by deactivating these degenerative processes one at a time in each of the three OA cases.

RESULTS

In the injurious loading of intact and physiological loading of degenerated cartilage, collagen degradation drives degenerative changes through fixed charge density loss and tissue hydration rise. In contrast, the two later mechanisms were more prominent in the focal defect cartilage model.

CONCLUSION

The virtual knock-out models reveal that injurious loading to intact cartilage and physiological loading to degenerated cartilage induce initial degenerative changes in the collagen network, whereas, in the presence of a focal cartilage defect, mechanical loading initially causes proteoglycans (PG) depletion, before changes in the collagen fibril network occur.

Evaluation of Sodium Relaxation Times and Concentrations in the Achilles Tendon Using MRI

Sodium magnetic resonance imaging (MRI) can be used to evaluate the change in the proteoglycan content in Achilles tendons (ATs) of patients with different AT pathologies by measuring the ²³Na signal-to-noise ratio (SNR). As ²³Na SNR alone is difficult to compare between different studies, because of the high influence of hardware configurations and sequence settings on the SNR, we further set out to measure the apparent tissue sodium content (aTSC) in the AT as a better comparable parameter. Ten healthy controls and one patient with tendinopathy in the AT were examined using a clinical 3 Tesla (T) MRI scanner in conjunction with a dual tuned ¹H/²³Na surface coil to measure ²³Na SNR and aTSC in their ATs. ²³Na T₁ and T₂* of the AT were also measured for three controls to correct for different relaxation behavior. The results were as follows: ²³Na SNR = 11.7 ± 2.2, aTSC = 82.2 ± 13.9 mM, ²³Na T₁ = 20.4 ± 2.4 ms, ²³Na T_2s* = 1.4 ± 0.4 ms, and ²³Na T_2l* = 13.9 ± 0.8 ms for the whole AT of healthy controls with significant regional differences. These are the first reported aTSCs and ²³Na relaxation times for the AT using sodium MRI and may serve for future comparability in different studies regarding examinations of diseased ATs with sodium MRI.

Micro- and Macroscale Assessment of Posterior Cruciate Ligament Functionality Based on Advanced MRI Techniques

T2 mapping assesses tissue ultrastructure and composition, yet the association of imaging features and tissue functionality is oftentimes unclear. This study aimed to elucidate this association for the posterior cruciate ligament (PCL) across the micro- and macroscale and as a function of loading. Ten human cadaveric knee joints were imaged using a clinical 3.0T scanner and high-resolution morphologic and T2 mapping sequences. Emulating the posterior drawer test, the joints were imaged in the unloaded (δ0) and loaded (δ1) configurations. For the entire PCL, its subregions, and its osseous insertion sites, loading-induced changes were parameterized as summary statistics and texture variables, i.e., entropy, homogeneity, contrast, and variance. Histology confirmed structural integrity. Statistical analysis was based on parametric and non-parametric tests. Mean PCL length (37.8 ± 1.8 mm [δ0]; 44.0 ± 1.6 mm [δ1] [p < 0.01]), mean T2 (35.5 ± 2.0 ms [δ0]; 37.9 ± 1.3 ms [δ1] [p = 0.01]), and mean contrast values (4.0 ± 0.6 [δ0]; 4.9 ± 0.9 [δ1] [p = 0.01]) increased significantly under loading. Other texture features or ligamentous, osseous, and meniscal structures remained unaltered. Beyond providing normative T2 values across various scales and configurations, this study suggests that ligaments can be imaged morphologically and functionally based on joint loading and advanced MRI acquisition and post-processing techniques to assess ligament integrity and functionality in variable diagnostic contexts.

Fourier Transform Infrared Microspectroscopy Combined with Principal Component Analysis and Artificial Neural Networks for the Study of the Effect of β-Hydroxy-β-Methylbutyrate (HMB) Supplementation on Articular Cartilage

The potential of Fourier Transform infrared microspectroscopy (FTIR microspectroscopy) and multivariate analyses were applied for the classification of the frequency ranges responsible for the distribution changes of the main components of articular cartilage (AC) that occur during dietary β-hydroxy-β-methyl butyrate (HMB) supplementation. The FTIR imaging analysis of histological AC sections originating from 35-day old male piglets showed the change in the collagen and proteoglycan contents of the HMB-supplemented group compared to the control. The relative amount of collagen content in the superficial zone increased by more than 23% and in the middle zone by about 17%, while no changes in the deep zone were observed compared to the control group. Considering proteoglycans content, a significant increase was registered in the middle and deep zones, respectively; 62% and 52% compared to the control. AFM nanoindentation measurements collected from animals administered with HMB displayed an increase in AC tissue stiffness by detecting a higher value of Young’s modulus in all investigated AC zones. We demonstrated that principal component analysis and artificial neural networks could be trained with spectral information to distinguish AC histological sections and the group under study accurately. This work may support the use and effectiveness of FTIR imaging combined with multivariate analyses as a quantitative alternative to traditional collagenous tissue-related histology.

Longitudinal T2 Mapping and Texture Feature Analysis in the Detection and Monitoring of Experimental Post-Traumatic Cartilage Degeneration

Background: Traumatic cartilage injuries predispose articulating joints to focal cartilage defects and, eventually, posttraumatic osteoarthritis. Current clinical-standard imaging modalities such as morphologic MRI fail to reliably detect cartilage trauma and to monitor associated posttraumatic degenerative changes with oftentimes severe prognostic implications. Quantitative MRI techniques such as T2 mapping are promising in detecting and monitoring such changes yet lack sufficient validation in controlled basic research contexts. Material and Methods: 35 macroscopically intact cartilage samples obtained from total joint replacements were exposed to standardized injurious impaction with low (0.49 J, n = 14) or high (0.98 J, n = 14) energy levels and imaged before and immediately, 24 h, and 72 h after impaction by T2 mapping. Contrast, homogeneity, energy, and variance were quantified as features of texture on each T2 map. Unimpacted controls (n = 7) and histologic assessment served as reference. Results: As a function of impaction energy and time, absolute T2 values, contrast, and variance were significantly increased, while homogeneity and energy were significantly decreased. Conclusion: T2 mapping and texture feature analysis are sensitive diagnostic means to detect and monitor traumatic impaction injuries of cartilage and associated posttraumatic degenerative changes and may be used to assess cartilage after trauma to identify “cartilage at risk”.