Biomechanical properties of articular cartilage in different regions and sites of the knee joint: acquisition of osteochondral allografts

Osteochondral allograft (OCA) transplantation involves grafting of natural hyaline cartilage and supporting subchondral bone into the cartilage defect area to restore its biomechanical and tissue structure. However, differences in biomechanical properties and donor-host matching may impair the integration of articular cartilage (AC). This study analyzed the biomechanical properties of the AC in different regions of different sites of the knee joint and provided a novel approach to OCA transplantation. Intact stifle joints from skeletally mature pigs were collected from a local abattoir less than 8 h after slaughter. OCAs were collected from different regions of the joints. The patella and the tibial plateau were divided into medial and lateral regions, while the trochlea and femoral condyle were divided into six regions. The OCAs were analyzed and compared for Young's modulus, the compressive modulus, and cartilage thickness. Young's modulus, cartilage thickness, and compressive modulus of OCA were significantly different in different regions of the joints. A negative correlation was observed between Young's modulus and the proportion of the subchondral bone (r =  - 0.4241, P < 0.0001). Cartilage thickness was positively correlated with Young's modulus (r = 0.4473, P < 0.0001) and the compressive modulus (r = 0.3678, P < 0.0001). During OCA transplantation, OCAs should be transplanted in the same regions, or at the closest possible regions to maintain consistency of the biomechanical properties and cartilage thickness of the donor and recipient, to ensure smooth integration with the surrounding tissue. A 7 mm depth achieved a higher Young's modulus, and may represent the ideal length.

Comparison of MRI T1, T2, and T2* mapping with histology for assessment of intervertebral disc degeneration in an ovine model

An easy, reliable, and time-efficient standardized approach for assessing lumbar intervertebral disc (IVD) degeneration with relaxation times measurements in pre-clinical and clinical studies is lacking. This prospective study aims to determine the most appropriate method for lumbar IVD degeneration (IDD) assessment in sheep by comparing three quantitative MRI sequences (variable-flip-angle T1 mapping, and multi-echo T2 and T2* mapping), correlating them with Pfirrmann grading and histology. Strong intra- and interrater agreements were found for Nucleus pulposus (NP) regions-of-interest (ROI). T1, T2, and T2* mapping correlated with Pfirrmann grading and histological scoring (p < 0.05) except for the most ventral rectangular ROI on T2 maps. Correlations were excellent for all of the T1 ROIs and the T2* NP ROIs. Highly significant differences in T1 values were found between all Pfirrmann grades except between grades I/II and between grades III/IV. Significant differences were identified in the T2 and the T2* values between all grades except between grades I/III. T1, T2, and T2* relaxation times measurements of the NP are an accurate and time-efficient tool to assess lumbar IDD in sheep. Variable-flip-angle T1 mapping may be further considered as a valuable method to investigate IDD and to assess the efficacy of regenerative treatments in longitudinal studies.

Cutoff points of T1 rho/T2 mapping relaxation times distinguishing early-stage and advanced osteoarthritis

INTRODUCTION

The histopathology grading system is the gold standard post-operative method to evaluate cartilage degeneration in knee osteoarthritis (OA). Magnetic resonance imaging (MRI) T1 rho/T2 mapping imaging can be used for preoperative detection. An association between histopathology and T1 rho/T2 mapping relaxation times was suggested in previous research. However, the cutoff point was not determined among different histopathology grades. Our study aimed to determine the cutoff point of T1 rho/T2 mapping.

MATERIAL AND METHODS

T1 rho/T2 mapping images were acquired from 80 samples before total knee replacements. Then the histopathology grading system was applied.

RESULTS

The mean T1 rho/T2 mapping relaxation times of 80 samples were 39.17 ms and 37.98 ms respectively. Significant differences were found in T1 rho/T2 mapping values between early-stage and advanced OA (p < 0.001). The cutoff point for T1 rho was 33 ms with a sensitivity of 94.12 (95% CI: 80-99.3) and a specificity of 91.30 (95% CI: 79.2-97.6). The cutoff point for T2 mapping was suggested as 35.04 ms with a sensitivity of 88.24 (95% CI: 72.5-96.7) and specificity of 97.83 (95% CI: 88.5-99.9). After bootstrap simulation, the 95% CI of the T1 rho/T2 mapping cutoff point was estimated as 29.36 to 36.32 ms and 34.8 to 35.04 ms respectively. The area under the PR curve of T1 rho/T2 mapping was 0.972 (95% CI: 0.925-0.992) and 0.949 (95% CI: 0.877-0.989) respectively.

CONCLUSIONS

The cutoff point of T1 rho relaxation times, which was suggested as 33 ms, could be used to distinguish early-stage and advanced OA.

Detecting Articular Cartilage and Meniscus Deformation Effects Using Magnetization Transfer Ultrashort Echo Time (MT-UTE) Modeling during Mechanical Load Application: Ex Vivo Feasibility Study

OBJECTIVE

Ultrashort echo time (UTE) magnetic resonance imaging (MRI) sequences have improved imaging of short T2 musculoskeletal (MSK) tissues. UTE-MRI combined with magnetization transfer modeling (UTE-MT) has demonstrated robust assessment of MSK tissues. This study aimed to investigate the variation of UTE-MT measures under mechanical loading in tibiofemoral cartilage and meniscus of cadaveric knee joints.

DESIGN

Fourteen knee joints from young (n = 8, 42 ± 12 years old) and elderly (n = 6, 89 ± 4 years old) donors were scanned on a 3-T scanner under 3 loading conditions: load = 300 N (Load1), load = 500 N (Load2), and load = 0 N (Unload). UTE-MT sequences were performed at each loading condition. Macromolecular proton fraction (MMF) was calculated from UTE-MT modeling. Wilcoxon rank sum test was used to examine the MRI data differences between loading conditions.

RESULTS

For young donors, MMF increased in all grouped regions of interest (meniscus [M], femoral articular cartilage [FAC], tibial articular cartilage [TAC], articular cartilage regions covered by meniscus [AC-MC], and articular cartilage regions uncovered by meniscus [AC-UC]) when the load increased from 300 to 500 N. The increases in MMF were significant for M (13.3%, P < 0.01) and AC-MC (9.2%, P = 0.04). MMF decreased in all studied regions after unloading, which was significant only for AC-MC (-8.9%, P = 0.01). For elderly donors, MRI parameters did not show significant changes by loading or unloading.

CONCLUSION

This study highlights the potential of the UTE-MT modeling combined with knee loading in differentiating between normal and abnormal knees. Average tissue deformation effects were likely higher and more uniformly distributed in the joints of young donors compared with elderly donors.

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.

No pressure, no diamonds? - Static vs. dynamic compressive in-situ loading to evaluate human articular cartilage functionality by functional MRI

Biomechanical Magnetic Resonance Imaging (MRI) of articular cartilage, i.e. its imaging under loading, is a promising diagnostic tool to assess the tissue's functionality in health and disease. This study aimed to assess the response to static and dynamic loading of histologically intact cartilage samples by functional MRI and pressure-controlled in-situ loading. To this end, 47 cartilage samples were obtained from the medial femoral condyles of total knee arthroplasties (from 24 patients), prepared to standard thickness, and placed in a standard knee joint in a pressure-controlled whole knee-joint compressive loading device. Cartilage samples' responses to static (i.e. constant), dynamic (i.e. alternating), and no loading, i.e. free-swelling conditions, were assessed before (δ₀), and after 30 min (δ₁) and 60 min (δ₂) of loading using serial T1ρ maps acquired on a 3.0T clinical MRI scanner (Achieva, Philips). Alongside texture features, relative changes in T1ρ (Δ₁, Δ₂) were determined for the upper and lower sample halves and the entire sample, analyzed using appropriate statistical tests, and referenced to histological (Mankin scoring) and biomechanical reference measures (tangent stiffness). Histological, biomechanical, and T1ρ sample characteristics at δ₀ were relatively homogenous in all samples. In response to loading, relative increases in T1ρ were strong and significant after dynamic loading (Δ₁ = 10.3 ± 17.0%, Δ₂ = 21.6 ± 21.8%, p = 0.002), while relative increases in T1ρ after static loading and in controls were moderate and not significant. Generally, texture features did not demonstrate clear loading-related associations underlying the spatial relationships of T1ρ. When realizing the clinical translation, this in-situ study suggests that serial T1ρ mapping is best combined with dynamic loading to assess cartilage functionality in humans based on advanced MRI techniques.

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

BACKGROUND

Articular cartilage degeneration is the hallmark change of osteoarthritis, a severely disabling disease with high prevalence and considerable socioeconomic and individual burden. Early, potentially reversible cartilage degeneration is characterized by distinct changes in cartilage composition and ultrastructure, while the tissue's morphology remains largely unaltered. Hence, early degenerative changes may not be diagnosed by clinical standard diagnostic tools.

METHODS

Against this background, this study introduces a novel method to determine the tissue composition non-invasively. Our method involves quantitative MRI parameters (i.e., T₁, T_1ρ, T₂ and [Formula: see text] maps), compositional reference measurements (i.e., microspectroscopically determined local proteoglycan [PG] and collagen [CO] contents) and machine learning techniques (i.e., artificial neural networks [ANNs] and multivariate linear models [MLMs]) on 17 histologically grossly intact human cartilage samples.

RESULTS

Accuracy and precision were higher in ANN-based predictions than in MLM-based predictions and moderate-to-strong correlations were found between measured and predicted compositional parameters.

CONCLUSION

Once trained for the clinical setting, advanced machine learning techniques, in particular ANNs, may be used to non-invasively determine compositional features of cartilage based on quantitative MRI parameters with potential implications for the diagnosis of (early) degeneration and for the monitoring of therapeutic outcomes.

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”.

Identifying the imaging correlates of cartilage functionality based on quantitative MRI mapping - The collagenase exposure model

Cartilage functionality is determined by tissue structure and composition. If altered, cartilage is predisposed to premature degeneration. This pathomimetical study of early osteoarthritis evaluated the dose-dependant effects of collagenase-induced collagen disintegration and proteoglycan depletion on cartilage functionality as assessed by serial T1, T1ρ, T2, and T2* mapping under loading. 30 human femoral osteochondral samples underwent imaging on a clinical 3.0 T MRI scanner (Achieva, Philips) in the unloaded reference configuration (δ₀) and under pressure-controlled quasi-static indentation loading to 15.1 N (δ₁) and to 28.6 N (δ₂). Imaging was performed before and after exposure to low (LC, 0.5 mg/mL; n = 10) or high concentration (HC, 1.5 mg/mL; n = 10) of collagenase. Untreated samples served as controls (n = 10). Loading responses were determined for the entire sample and the directly loaded (i.e. sub-pistonal) and bilaterally adjacent (i.e. peri‑pistonal) regions, referenced histologically, quantified as relative changes, and analysed using adequate parametric and non-parametric statistical tests. Dose-dependant surface disintegration and tissue loss were reflected by distinctly different pre- and post-exposure response-to-loading patterns. While T1 generally decreased with loading, regardless of collagenase exposure, T1ρ increased significantly after HC exposure (p = 0.008). Loading-induced decreases in T2 were significant after LC exposure (p = 0.006), while changes in T2* were ambiguous. In conclusion, aberrant loading-induced changes in T2 and T1ρ reflect moderate and severe matrix changes, respectively, and indicate the close interrelatedness of matrix changes and functionality in cartilage.

A serial multiparametric quantitative magnetic resonance imaging study to assess proteoglycan depletion of human articular cartilage and its effects on functionality

Water, collagen, and proteoglycans determine articular cartilage functionality. If altered, susceptibility to premature degeneration is increased. This study investigated the effects of enzymatic proteoglycan depletion on cartilage functionality as assessed by advanced Magnetic Resonance Imaging (MRI) techniques under standardized loading. Lateral femoral condylar cartilage-bone samples from patients undergoing knee replacement (n = 29) were serially imaged by Proton Density-weighted and T1, T1ρ, T2, and T2* mapping sequences on a clinical 3.0 T MRI scanner (Achieva, Philips). Using pressure-controlled indentation loading, samples were imaged unloaded and quasi-statically loaded to 15.1 N and 28.6 N, and both before and after exposure to low-concentrated (LT, 0.1 mg/mL, n = 10) or high-concentrated trypsin (HT, 1.0 mg/mL, n = 10). Controls were not treated (n = 9). Responses to loading were assessed for the entire sample and regionally, i.e. sub- and peri-pistonally, and zonally, i.e. upper and lower sample halves. Trypsin effects were quantified as relative changes (Δ), analysed using appropriate statistical tests, and referenced histologically. Histological proteoglycan depletion was reflected by significant sub-pistonal decreases in T1 (p = 0.003) and T2 (p = 0.008) after HT exposure. Loading-induced changes in T1ρ and T2* were not related. In conclusion, proteoglycan depletion alters cartilage functionality and may be assessed using serial T1 and T2 mapping under loading.

An MRI-compatible varus-valgus loading device for whole-knee joint functionality assessment based on compartmental compression: a proof-of-concept study

OBJECTIVE

Beyond static assessment, functional techniques are increasingly applied in magnetic resonance imaging (MRI) studies. Stress MRI techniques bring together MRI and mechanical loading to study knee joint and tissue functionality, yet prototypical axial compressive loading devices are bulky and complex to operate. This study aimed to design and validate an MRI-compatible pressure-controlled varus-valgus loading device that applies loading along the joint line.

METHODS

Following the device's thorough validation, we demonstrated proof of concept by subjecting a structurally intact human cadaveric knee joint to serial imaging in unloaded and loaded configurations, i.e. to varus and valgus loading at 7.5 kPa (= 73.5 N), 15 kPa (= 147.1 N), and 22.5 kPa (= 220.6 N). Following clinical standard (PDw fs) and high-resolution 3D water-selective cartilage (WATSc) sequences, we performed manual segmentations and computations of morphometric cartilage measures. We used CT and radiography (to quantify joint space widths) and histology and biomechanics (to assess tissue quality) as references.

RESULTS

We found (sub)regional decreases in cartilage volume, thickness, and mean joint space widths reflective of areal pressurization of the medial and lateral femorotibial compartments.

DISCUSSION

Once substantiated by larger sample sizes, varus-valgus loading may provide a powerful alternative stress MRI technique.

Functional MRI Mapping of Human Meniscus Functionality and its Relation to Degeneration

Meniscus pathology may promote early osteoarthritis. This study assessed human meniscus functionality (i.e. its response to loading) ex vivo based on quantitative T1, T1ρ, and T2 mapping as a function of histological degeneration and loading. Forty-five meniscus samples of variable degeneration were harvested from the lateral meniscus body region of 45 patients during total knee arthroplasties. Samples underwent serial mapping on a 3.0-T MRI scanner (Achieva, Philips) using a force-controlled and torque-inducing compressive loading device. Samples were measured at three loading positions, i.e. unloaded, loaded to 2 bar (compression force 37 N) and 4 bar (69 N). Histology (Pauli classification) and biomechanics (Elastic Modulus) served as references. Based on histology, samples were trichotomized as grossly intact (n = 14), mildly degenerative (n = 16), and moderate-to-severely degenerative (n = 15) and analyzed using appropriate parametric and non-parametric tests. For T1, we found loading-induced decreases in all samples, irrespective of degeneration. For T1ρ, zonal increases in intact (apex) and decreases in degenerative samples (base) were found, while for T2, changes were ambiguous. In conclusion, force-controlled loading and serial MR imaging reveal response-to-loading patterns in meniscus. Zonal T1ρ response-to-loading patterns are most promising in differentiating degeneration, while T1 and T2 aren’t clearly related to degeneration.and may provide an imaging-based indication of functional tissue properties.

Histological Grade and Magnetic Resonance Imaging Quantitative T1rho/T2 Mapping in Osteoarthritis of the Knee: A Study in 20 Patients

Background

Magnetic resonance imaging (MRI) of osteoarthritis (OA) of the knee is a preoperative method of joint assessment. Histology of the joint is invasive and performed after surgery. T1rho/T2 MRI mapping is a new preoperative method of quantifying joint changes. This study aimed to analyze and compare the histological changes in the joint cartilage with the use of quantitative T1rho/T2 MRI mapping in patients with OA of the knee.

Material/Methods

Twenty patients with OA of the knee (20 knees) underwent preoperative MRI with T1rho mapping, T2 mapping, T1-weighted, and T2-weighted fat-suppressed MRI sequences. The degree of OA of the knee on MRI was graded according to the Osteoarthritis Research Society International (OARSI) criteria and the Kellgren-Lawrence grading system. Histological grading of OA used the OARSI criteria. Four tibiofemoral condyles were assessed histologically, and the degree of cartilage destruction was determined using the OARSI criteria. Two investigators performed cartilage segmentation for T1rho/T2 values.

Results

Histology of the four knee joint condyles confirmed mild to severe OA. The histology of the cartilage thickness (P<0.001) and the MRI findings of the distal medial condyle (P<0.00) were significantly different from the other three knee joint condyles. The T2 and T1rho values of each condyle were significantly correlated with the histological grade (II–IV) of the joint condyles, including the cartilage volume, cartilage defects, thickness, and bone lesions (P<0.05).

Conclusions

In 20 patients with OA of the knee, preoperative T2/T1rho MRI identified Grade II–IV OA changes in the joint.

Human articular cartilage mechanosensitivity is related to histological degeneration - a functional MRI study

OBJECTIVE

To investigate changes in response to sequential pressure-controlled loading and unloading in human articular cartilage of variable histological degeneration using serial T1ρ mapping.

METHOD

We obtained 42 cartilage samples of variable degeneration from the medial femoral condyles of 42 patients undergoing total knee replacement. Samples were placed in a standardized artificial knee joint within an MRI-compatible whole knee-joint compressive loading device and imaged before (δ₀), during (δ_ld1, δ_ld2, δ_ld3, δ_ld4, δ_ld5) and after (δ_rl1, δ_rl2, δ_rl3, δ_rl4, δ_rl5) pressure-controlled loading to 0.663 ± 0.021 kN (94% body weight) using serial T1ρ mapping (spin-lock multigradient echo sequence; 3.0T MRI system [Achieva, Philips]). Reference assessment included histology (Mankin scoring) and conventional biomechanics (Tangent stiffness). We dichotomized sample into intact (n = 21) and degenerative (n = 21) based on histology and analyzed data using Mann Whitney, Kruskal Wallis, one-way ANOVA tests and Spearman's correlation, respectively.

RESULTS

At δ₀, we found no significant differences between intact and degenerative samples, while the response-to-loading patterns were distinctly different. In intact samples, T1ρ increases were consistent and non-significant, while in degenerative samples, T1ρ increases were significantly higher (P = 0.004, δ0 vs δld1, δ0 vs δld3), yet undulating and variable. With unloading, T1ρ increases subsided, yet were persistently elevated beyond δ0.

CONCLUSION

Cartilage mechanosensitivity is related to histological degeneration and assessable by serial T1ρ mapping. Unloaded, T1ρ characteristics are not significantly different in intact vs degenerative cartilage, while load bearing is organized in intact cartilage and disorganized in degenerative cartilage.

Magnetic resonance imaging (MRI) studies of knee joint under mechanical loading: Review

Osteoarthritis (OA) is a very common disease that affects the human knee joint, particularly the articular cartilage and meniscus components which are regularly under compressive mechanical loads. Early-stage OA diagnosis is essential as it allows for timely intervention. The primary non-invasive approaches currently available for OA diagnosis include magnetic resonance imaging (MRI), which provides excellent soft tissue contrast at high spatial resolution. MRI-based knee investigation is usually performed on joints at rest or in a non-weight-bearing condition that does not mimic the actual physiological condition of the joint. This discrepancy may lead to missed detections of early-stage OA or of minor lesions. The mechanical properties of degenerated musculoskeletal (MSK) tissues may vary markedly before any significant morphological or structural changes detectable by MRI. Recognizing distinct deformation characteristics of these tissues under known mechanical loads may reveal crucial joint lesions or mechanical malfunctions which result from early-stage OA. This review article summarizes the large number of MRI-based investigations on knee joints under mechanical loading which have been reported in the literature including the corresponding MRI measures, the MRI-compatible devices employed, and potential challenges due to the limitations of clinical MRI sequences.

A multi-purpose force-controlled loading device for cartilage and meniscus functionality assessment using advanced MRI techniques

Response to loading of soft tissues as assessed by advanced magnetic resonance imaging (MRI) techniques is a promising approach to evaluate tissue functionality beyond (statically obtained) structural and compositional features. As cartilage and meniscus pathologies are closely intertwined in osteoarthritis (OA) and beyond, both tissues should ideally be studied to elucidate further the underlying mechanisms involved in load transmission and its failure leading to OA. Hence, we devised, constructed and validated a dedicated MRI-compatible pneumatic force-controlled loading device to study cartilage and meniscus functionality in a standardized and reproducible manner and in reference to alternative tissue evaluation methods. Mechanical reference measurements using digital force sensors confirmed the reproducible application of forces in the range of 0-76N. To demonstrate the device's utility in a basic research context, MRI measurements of human articular cartilage (obtained from the lateral femoral condyle, n = 5) and meniscus (obtained from lateral meniscus body, n = 5) were performed in the unloaded (δ₀) and loaded configurations (δ₁: [cartilage] 0.75 bar corresponding to 15.1 N, [meniscus] 2 bar corresponding to 37.1 N; δ₂: [cartilage] 1.5 bar corresponding to 28.6 N, [meniscus] 4 bar corresponding to 69.1 N). Cartilage samples were directly indented, while meniscus samples were subject to torque-induced compression using a dedicated lever compression device. Morphological MR Imaging using Proton Density-weighted sequences and quantitative MR Imaging using T2 and T1ρ mapping were performed serially and at high resolution. For reference, samples underwent subsequent biomechanical and histological reference evaluation. In conclusion, the force-controlled loading device has been validated for the non-invasive response-to-loading assessment of human cartilage and meniscus samples by advanced MRI techniques. Hereby, both tissues may be functionally evaluated in combination, beyond mere static analysis and in reference to histological and biomechanical measures.

Towards Patient-Specific Computational Modelling of Articular Cartilage on the Basis of Advanced Multiparametric MRI Techniques

Cartilage degeneration is associated with tissue softening and represents the hallmark change of osteoarthritis. Advanced quantitative Magnetic Resonance Imaging (qMRI) techniques allow the assessment of subtle tissue changes not only of structure and morphology but also of composition. Yet, the relation between qMRI parameters on the one hand and microstructure, composition and the resulting functional tissue properties on the other hand remain to be defined. To this end, a Finite-Element framework was developed based on an anisotropic constitutive model of cartilage informed by sample-specific multiparametric qMRI maps, obtained for eight osteochondral samples on a clinical 3.0 T MRI scanner. For reference, the same samples were subjected to confined compression tests to evaluate stiffness and compressibility. Moreover, the Mankin score as an indicator of histological tissue degeneration was determined. The constitutive model was optimized against the resulting stress responses and informed solely by the sample-specific qMRI parameter maps. Thereby, the biomechanical properties of individual samples could be captured with good-to-excellent accuracy (mean R² [square of Pearson's correlation coefficient]: 0.966, range [min, max]: 0.904, 0.993; mean Ω [relative approximated error]: 33%, range [min, max]: 20%, 47%). Thus, advanced qMRI techniques may be complemented by the developed computational model of cartilage to comprehensively evaluate the functional dimension of non-invasively obtained imaging biomarkers. Thereby, cartilage degeneration can be perspectively evaluated in the context of imaging and biomechanics.

Detection of Early-Stage Degeneration in Human Articular Cartilage by Multiparametric MR Imaging Mapping of Tissue Functionality

To assess human articular cartilage tissue functionality by serial multiparametric quantitative MRI (qMRI) mapping as a function of histological degeneration. Forty-nine cartilage samples obtained during total knee replacement surgeries were placed in a standardized artificial knee joint within an MRI-compatible compressive loading device and imaged in situ and at three loading positions, i.e. unloaded, at 2.5 mm displacement (20% body weight [BW]) and at 5 mm displacement (110% BW). Using a clinical 3.0 T MRI system (Achieva, Philips), serial T1, T1ρ, T2 and T2* maps were generated for each sample and loading position. Histology (Mankin scoring) and biomechanics (Young’s modulus) served as references. Samples were dichotomized as intact (int, n = 27) or early degenerative (deg, n = 22) based on histology and analyzed using repeated-measures ANOVA and unpaired Student’s t-tests after log-transformation. For T1ρ, T2 and T2*, significant loading-induced differences were found in deg (in contrast to int) samples, while for T1 significant decreases in all zones were observed, irrespective of degeneration. In conclusion, cartilage functionality may be visualized using serial qMRI parameter mapping and the response-to-loading patterns are associated with histological degeneration. Hence, loading-induced changes in qMRI parameter maps provide promising surrogate parameters of tissue functionality and status in health and disease.

Linking hyperelastic theoretical models and experimental data of vaginal tissue through histological data

Mechanical characterization of living tissues and computer-based simulations related to medical issues, has become increasingly important to improve diagnostic processes and treatments evaluation. This work proposes a link between the mechanical testing and the material model predictions through histological data of vaginal tissue. Histological data was used to link tensile testing experiments with material-dependent parameters; the approach was adequate to capture the nonlinear response of ovine vaginal tissue over a large strain range. The experimental data obtained on a previous study, has two main components: tensile testing and histological analysis of the ovine vaginal tissue. Uniaxial tensile test data and histological data were collected from three sheep groups: virgins, pregnant and parous. The distal part of vaginal wall was selected since it is prone to tears induced by vaginal delivery. The HGO (Holzapfel-Gasser-Ogden) model parameters were fitted using a stochastic approach, namely the Simple Genetic Algorithm (SGA). The SGA was able to fit the experimental data successfully (R² > 0.986). The dimensionless coefficient ξ, was highly correlated with histological data. The ratio was seen to increase linearly with increasing collagen content. Coefficient ξ brings a new way of interpreting and understanding experimental data; it connects the nonlinear mechanical behaviour (tensile test) with tissue's morphology (histology). It can be used as an 'inverse' (approximate) method to estimate the mechanical properties without direct experimental measurements, through basic histology. In this context, the proposed methodology appears very promising in estimating the response of the tissue via histological information.

Cartilage-on-cartilage contact: effect of compressive loading on tissue deformations and structural integrity of bovine articular cartilage

OBJECTIVE

This study aims to characterize the deformations in articular cartilage under compressive loading and link these to changes in the extracellular matrix constituents described by magnetic resonance imaging (MRI) relaxation times in an experimental model mimicking in vivo cartilage-on-cartilage contact.

DESIGN

Quantitative MRI images, T₁, T₂ and T_1ρ relaxation times, were acquired at 9.4T from bovine femoral osteochondral explants before and immediately after loading. Two-dimensional intra-tissue displacement and strain fields under cyclic compressive loading (350N) were measured using the displacement encoding with stimulated echoes (DENSE) method. Changes in relaxation times in response to loading were evaluated against the deformation fields.

RESULTS

Deformation fields showed consistent patterns among all specimens, with maximal strains at the articular surface that decrease with tissue depth. Axial and transverse strains were maximal around the center of the contact region, whereas shear strains were minimal around the contact center but increased towards contact edges. A decrease in T₂ and T_1ρ was observed immediately after loading whereas the opposite was observed for T₁. No correlations between cartilage deformation patterns and changes in relaxation times were observed.

CONCLUSIONS

Displacement encoding combined with relaxometry by MRI can noninvasively monitor the cartilage biomechanical and biochemical properties associated with loading. The deformation fields reveal complex patterns reflecting the depth-dependent mechanical properties, but intra-tissue deformation under compressive loading does not correlate with structural and compositional changes. The compacting effect of cyclic compression on the cartilage tissue was revealed by the change in relaxation time immediately after loading.

Revisiting Radiographic L5-S1 Parallelism Using MRI T1 Mapping

Purpose:

Thirty years ago, we reported that parallel aspect of the L5-S1 disc on a lateral view of the spine might be considered to be an initial stage of disk degeneration. The current study represents an attempt to increase the validity of parallel sign on conventional radiograph using MR real T1 mapping.

Methods:

Forty-four young asymptomatic volunteers (mean age 21.6 ± 2.3) underwent lumbar spine MRI, twice the same day, morning and afternoon. Dedicated sequences using the inversion-recovery technique were used to calculate the T1 relaxation time. A region of interest (ROI) representing the nucleus pulposus was defined in each disk. The volunteers were stratified according to the presence or absence of a parallel morphology of L5-S1. Correlation between endplates angles, sacral slopes and T1 values were then evaluated.

Results:

L5-S1 space looks parallel for angles <10° (mean value 6.9° ± 1.4°). Sacral slope was lower in parallel disks (31.7 ± 4.9° vs. 40.1 ± 5.6°), showing a significant difference of 8.4° (p < 0.05). The T1 relaxation values show a significant difference between the two groups (p < 0.05) with a difference of 96 ms for the morning (1090.9 ± 33.3 ms for the parallel group and 1186.9 ± 41.2 ms for the non-parallel) and 121.9 ms for the afternoon (respectively 1004.7 ± 22.2 ms and 1126.6 ± 12.9 ms).

Conclusion:

The difference between the two groups suggests that parallel morphology of the L5-S1 disk is associated with lower water content.

Ultrashort time-to-echo quantitative magnetic resonance imaging of the triangular fibrocartilage: differences in position

PURPOSE

To compare T2* values of the triangular fibrocartilage (TFC) obtained by ultrashort time-to-echo (UTE) techniques at the neutral position, ulnar flexion of the wrist, and pronation of the forearm.

MATERIALS AND METHODS

MR imaging was performed in ten healthy volunteers with a 3-T MR system by using an eight-channel knee coil. Coronal wrist T2* maps from three-dimensional cone UTE pulse sequences were obtained at the neutral, ulnar flexion, and pronation positions (TR: 19 ms, TE: 0.032 ms/4 ms/8 ms/12 ms, FOV: 18 cm, matrix: 430 × 430, section thickness: 1.5 mm, scan time: 8 min 31 s). UTE-T2* maps were calculated on a pixel-by-pixel basis for all structures of the wrist visualized in the coronal planes. The entire region of interest (ROI) for TFC was manually delineated, and the average T2* value was calculated for each ROI by three radiologists. The Kruskal-Wallis test, Wilcoxon signed-rank test, or intraclass correlation coefficients (ICC) were used for statistics.

RESULTS

The difference in the average T2* value among the three groups according to the forearm/wrist position was significant (p < 0.001). The T2* value of the TFC at pronation (mean ± 2 SD: 7.92 ± 1.37 ms) was significantly lower than those at the neutral (10.08 ± 1.90 ms) and ulnar flexion positions (9.15 ± 1.03 ms) (p < 0.017). The ICC showed a substantial interobserver agreement in the T2* value measurements of the TFC (ICC = 0.986).

CONCLUSION

T2* relaxation time measurement of the TFC using UTE may be useful for assessing the loading effect by the forearm/wrist position.

KEY POINTS

• The T2* value of the TFC may reflect the biomechanics of the wrist joint. • Acute loading at pronation results in a decrease in the T2* value of the TFC. • Quantitative wrist UTE MRI was successfully performed in vivo.

Multiparametric MRI and Computational Modelling in the Assessment of Human Articular Cartilage Properties: A Comprehensive Approach

Quantitative magnetic resonance imaging (qMRI) is a promising approach to detect early cartilage degeneration. However, there is no consensus on which cartilage component contributes to the tissue's qMRI signal properties. T1, T1ρ, and T2^⁎ maps of cartilage samples (n = 8) were generated on a clinical 3.0-T MRI system. All samples underwent histological assessment to ensure structural integrity. For cross-referencing, a discretized numerical model capturing distinct compositional and structural tissue properties, that is, fluid fraction (FF), proteoglycan (PG) and collagen (CO) content and collagen fiber orientation (CFO), was implemented. In a pixel-wise and region-specific manner (central versus peripheral region), qMRI parameter values and modelled tissue parameters were correlated and quantified in terms of Spearman's correlation coefficient ρs. Significant correlations were found between modelled compositional parameters and T1 and T2^⁎, in particular in the central region (T1: ρs ≥ 0.7 [FF, CFO], ρs ≤ −0.8 [CO, PG]; T2^⁎: ρs ≥ 0.67 [FF, CFO], ρs ≤ −0.71 [CO, PG]). For T1ρ, correlations were considerably weaker and fewer (0.16 ≤ ρs ≤ −0.15). QMRI parameters are characterized in their biophysical properties and their sensitivity and specificity profiles in a basic scientific context. Although none of these is specific towards any particular cartilage constituent, T1 and T2^⁎ reflect actual tissue compositional features more closely than T1ρ.

Fatigue of soft fibrous tissues: Multi-scale mechanics and constitutive modeling

In recent experimental studies a possible damage mechanism of collagenous tissues mainly caused by fatigue was disclosed. In this contribution, a multi-scale constitutive model ranging from the tropocollagen (TC) molecule level up to bundles of collagen fibers is proposed and utilized to predict the elastic and inelastic long-term tissue response. Material failure of collagen fibrils is elucidated by a permanent opening of the triple helical collagen molecule conformation, triggered either by overstretching or reaction kinetics of non-covalent bonds. This kinetics is described within a probabilistic framework of adhesive detachments of molecular linkages providing collagen fiber integrity. Both intramolecular and interfibrillar linkages are considered. The final constitutive equations are validated against recent experimental data available in literature for both uniaxial tension to failure and the evolution of fatigue in subsequent loading cycles. All material parameters of the proposed model have a clear physical interpretation.

STATEMENT OF SIGNIFICANCE

Irreversible changes take place at different length scales of soft fibrous tissues under supra-physiological loading and alter their macroscopic mechanical properties. Understanding the evolution of those histologic pathologies under loading and incorporating them into a continuum mechanical framework appears to be crucial in order to predict long-term evolution of various diseases and to support the development of tissue engineering.

Non-invasive T1ρ mapping of the human cartilage response to loading and unloading

OBJECTIVE

To define the physiological response to sequential loading and unloading in histologically intact human articular cartilage using serial T1ρ mapping, as T1ρ is considered to indicate the tissue's macromolecular content.

METHOD

18 macroscopically intact cartilage-bone samples were obtained from the central lateral femoral condyles of 18 patients undergoing total knee replacement. Serial T1ρ mapping was performed on a clinical 3.0-T MRI system using a modified prostate coil. Spin-lock multiple gradient-echo sequences prior to, during and after standardized indentation loading (displacement controlled, strain 20%) were used to obtain seven serial T1ρ maps: unloaded (δ₀), quasi-statically loaded (indentation₁-indentation₃) and under subsequent relaxation (relaxation₁-relaxation₃). After manual segmentation, zonal and regional regions-of-interest were defined. ROI-specific relative changes were calculated and statistically assessed using paired t-tests. Histological (Mankin classification) and biomechanical (unconfined compression) evaluations served as references.

RESULTS

All samples were histologically and biomechanically grossly intact (Mankin sum: 1.8 ± 1.2; Young's Modulus: 0.7 ± 0.4 MPa). Upon loading, T1ρ consistently increased throughout the entire sample thickness, primarily subpistonally (indentation₁ [M ± SD]: 9.5 ± 7.8% [sub-pistonal area, SPA] vs 4.2 ± 5.8% [peri-pistonal area, PPA]; P < 0.001). T1ρ further increased with ongoing loading (indentation₃: 14.1 ± 8.1 [SPA] vs 7.7 ± 5.9% [PPA]; P < 0.001). Even upon unloading (i.e., relaxation), T1ρ persistently increased in time.

CONCLUSION

Serial T1ρ-mapping reveals distinct and complex zonal and regional changes in articular cartilage as a function of loading and unloading. Thereby, longitudinal adaptive processes in hyaline cartilage become evident, which may be used for the tissue's non-invasive functional characterization by T1ρ.

[Chondral and osteochondral defects : Representation by imaging methods]

Morphological imaging of cartilage at high resolution allows the differentiation of chondral and osteochondral lesions. Nowadays, magnetic resonance imaging is the principal diagnostic tool in the assessment of cartilage structure and composition. Conventional radiography, computed tomography, ultrasound or optical coherence tomography are adjunct diagnostic modalities in the assessment of cartilage pathologies. The present article discusses the up-to-date diagnostic practice of cartilage imaging in terms of its scientific basis and current clinical status, requirements, techniques and image interpretation. Innovations in the field such as functional MRI are discussed as well due to their mid- to long-term clinical perspective.

T2 MR imaging vs. computational modeling of human articular cartilage tissue functionality

The detection of early stages of cartilage degeneration remains diagnostically challenging. One promising non-invasive approach is to functionally assess the tissue response to loading by serial magnetic resonance (MR) imaging in terms of T2 mapping under simultaneous mechanical loading. As yet, however, it is not clear which cartilage component contributes to the tissue functionality as assessed by quantitative T2 mapping. To this end, quantitative T2 maps of histologically intact cartilage samples (n=8) were generated using a clinical 3.0-T MR imaging system. Using displacement-controlled quasi-static indentation loading, serial T2 mapping was performed at three defined strain levels and loading-induced relative changes were determined in distinct regions-of-interest. Samples underwent conventional biomechanical testing (by unconfined compression) as well as histological assessment (by Mankin scoring) for reference purposes. Moreover, an anisotropic hyperelastic constitutive model of cartilage was implemented into a finite element (FE) code for cross-referencing. In efforts to simulate the evolution of compositional and structural intra-tissue changes under quasi-static loading, the indentation-induced changes in quantitative T2 maps were referenced to underlying changes in cartilage composition and structure. These changes were parameterized as cartilage fluid, proteoglycan and collagen content as well as collagen orientation. On a pixel-wise basis, each individual component correlation with T2 relaxation times was determined by Spearman's ρ_s and significant correlations were found between T2 relaxation times and all four tissue parameters for all indentation strain levels. Thus, the biological changes in functional MR Imaging parameters such as T2 can further be characterized to strengthen the scientific basis of functional MRI techniques with regards to their perspective clinical applications.

Functional in situ assessment of human articular cartilage using MRI: a whole-knee joint loading device

The response to loading of human articular cartilage as assessed by magnetic resonance imaging (MRI) remains to be defined in relation to histology and biomechanics. Therefore, an MRI-compatible whole-knee joint loading device for the functional in situ assessment of cartilage was developed and validated in this study. A formalin-fixed human knee was scanned by computed tomography in its native configuration and digitally processed to create femoral and tibial bone models. The bone models were covered by artificial femoral and tibial articular cartilage layers in their native configuration using cartilage-mimicking polyvinyl siloxane. A standardized defect of 8 mm diameter was created within the artificial cartilage layer at the central medial femoral condyle, into which native cartilage samples of similar dimensions were placed. After describing its design and specifications, the comprehensive validation of the device was performed using a hydraulic force gauge and digital electronic pressure-sensitive sensors. Displacement-controlled quasi-static uniaxial loading to 2.5 mm [Formula: see text] and 5.0 mm [Formula: see text] of the mobile tibia versus the immobile femur resulted in forces of [Formula: see text] N [Formula: see text] and [Formula: see text] N [Formula: see text] (on the entire joint) and local pressures of [Formula: see text] MPa [Formula: see text] and [Formula: see text] MPa [Formula: see text] (at the site of the cartilage sample). Upon confirming the MRI compatibility of the set-up, the response to loading of macroscopically intact human articular cartilage samples ([Formula: see text]) was assessed on a clinical 3.0-T MR imaging system using clinical standard proton-density turbo-spin echo sequences and T2-weighted multi-spin echo sequences. Serial imaging was performed at the unloaded state [Formula: see text] and at consecutive loading positions (i.e. at [Formula: see text] and [Formula: see text]. Biomechanical unconfined compression testing (Young's modulus) and histological assessment (Mankin score) served as the standards of reference. All samples were histologically intact (Mankin score, [Formula: see text]) and biomechanically reasonably homogeneous (Young's modulus, [Formula: see text] MPa). They could be visualized in their entirety by MRI and significant decreases in sample height [[Formula: see text]: [Formula: see text] mm; [Formula: see text]: [Formula: see text] mm; [Formula: see text]: [Formula: see text] mm; [Formula: see text] (repeated-measures ANOVA)] as well as pronounced T2 signal decay indicative of tissue pressurization were found as a function of compressive loading. In conclusion, our compression device has been validated for the noninvasive response-to-loading assessment of human articular cartilage by MRI in a close-to-physiological experimental setting. Thus, in a basic research context cartilage may be functionally evaluated beyond mere static analysis and in reference to histology and biomechanics.

Real T1 relaxation time measurement and diurnal variation analysis of intervertebral discs in a healthy population of 50 volunteers

PURPOSE

To measure the real T1 relaxation time of the lumbar intervertebral discs in a young and healthy population, using different inversion recovery times, and assess diurnal variation.

MATERIAL AND METHODS

Intervertebral discs from D12 to S1 of 50 healthy volunteers from 18 to 25 years old were evaluated twice the same day, in the morning and in the late afternoon. Dedicated MRI sequences with different inversion recovery times (from 100 to 2500ms) were used to calculate the real T1 relaxation time. Three regions of interest (ROIs) were defined in each disc, the middle representing the nucleus pulposus (NP) and the outer parts the annulus fibrosus (AF) anterior and posterior. Diurnal variation and differences between each disc level were analyzed.

RESULTS

T1 mean values in the NP were 1142±12ms in the morning and 1085±13ms in the afternoon, showing a highly significant decrease of 57ms (p<0.001). A highly significant difference between the levels of the spine was found. The mean T1 of the anterior part of the AF was 577±9ms in the morning and 554±8ms in the afternoon. For the posterior part, the mean values were 633±8ms in the morning and 581±7ms in the evening. It shows a highly significant decrease of 23ms for the anterior part and 51ms for the posterior part (all p<0.001).

CONCLUSION

T1 mapping is a promising method of intervertebral disc evaluation. Significant diurnal variation and difference between levels of the lumbar spine were demonstrated. A potential use for longitudinal study in post-operative follow up or sport medicine needs to be evaluated.

Ex vivo quantitative multiparametric MRI mapping of human meniscus degeneration

OBJECTIVES

To evaluate the diagnostic performance of T1, T1ρ, T2, T2*, and UTE-T2* (ultrashort-echo time-enhanced T2*) mapping in the refined graduation of human meniscus degeneration with histology serving as standard-of-reference.

MATERIALS AND METHODS

This IRB-approved intra-individual comparative ex vivo study was performed on 24 lateral meniscus body samples obtained from 24 patients undergoing total knee replacement. Samples were assessed on a 3.0-T MRI scanner using inversion-recovery (T1), spin-lock multi-gradient-echo (T1ρ), multi-spin-echo (T2) and multi-gradient-echo (T2* and UTE-T2*) sequences to determine relaxation times of quantitative MRI (qMRI) parameters. Relaxation times were calculated on the respective maps, averaged to the entire meniscus and to its zones. Histologically, samples were analyzed on a four-point score according to Williams (0-III). QMRI results and Williams (sub)scores were correlated using Spearman's ρ, while Williams grade-dependent differences were assessed using Kruskal-Wallis and Dunn's tests. Sensitivities and specificities in the detection of intact (Williams grade [WG]-0) and severely degenerate meniscus (WG-II-III) were calculated.

RESULTS

Except for T2*, significant increases in qMRI parameters with increasing Williams grades were observed. T1, T1ρ, T2, and UTE-T2* exhibited high sensitivity and variable specificity rates. Significant marked-to-strong correlations were observed for these parameters with each other, with histological WGs and the subscores tissue integrity and cellularity.

CONCLUSIONS

QMRI mapping holds promise in the objective evaluation of human meniscus. Although sufficient discriminatory power of T1, T1ρ, T2, and UTE-T2* was only demonstrated for the histological extremes, these data may aid in the future MRI-based parameterization and quantification of human meniscus degeneration.