MR imaging pattern of tibial subchondral bone structure: considerations of meniscal coverage and integrity

OBJECTIVES

To compare regional differences in subchondral trabecular structure using high-resolution MRI in meniscus-covered/meniscus-uncovered tibia in cadaveric knees with intact/torn menisci.

MATERIALS AND METHODS

3D proton density CUBE MRI of 6 cadaveric knees without significant osteoarthritis (OA) was acquired, 0.25-mm resolution. Menisci were evaluated and classified intact or torn. MR data were transferred to ImageJ program to segment tibial 3D volume of interest (VOI). Data was subdivided into meniscus-covered/meniscus-uncovered regions. Segmented VOI was classified into binary data, trabeculae/bone marrow. The trabecular bone data was used to measure MR biomarkers (apparent subchondral plate-connected bone density (adapted from spine MR), apparent trabecular bone volume fraction, apparent mean trabecular thickness, apparent connectivity density, and structure model index (SMI)). Mean value of parameters was analyzed for the effects of meniscal tear/tibial coverage.

RESULTS

Nine torn menisci and 3 intact menisci were present. MR measures of bone varied significantly due to meniscal coverage/tear. Subchondral plate-connected bone density under covered meniscus regions increased from 10.9 to 23.5% with meniscal tear. Values increased in uncovered regions, 19.3% (intact) and 32.4% (torn). This reflects higher density when uncovered (p = 0.048) with meniscal tear (p = 0.007). Similar patterns were found for trabecular bone fraction (coverage p < 0.001, tear p = 0.047), trabecular thickness (coverage p = 0.03), connectivity density (coverage p = 0.002), and SMI (coverage p = 0.015).

CONCLUSION

Quantitative trabecular bone evaluation emphasizes intrinsic structural differences between meniscus-covered/meniscus-uncovered tibias. Results offer insight into bone adaptation with meniscal tear and support the hypothesis that subchondral bone plate-connected bone density could be important in early subchondral bone adaptation.

Subchondral bone microarchitecture analysis in the proximal tibia at 7-T MRI

Background Bone remodels in response to mechanical loads and osteoporosis results from impaired ability of bone to remodel. Bone microarchitecture analysis provides information on bone quality beyond bone mineral density (BMD). Purpose To compare subchondral bone microarchitecture parameters in the medial and lateral tibia plateau in individuals with and without fragility fractures. Material and Methods Twelve female patients (mean age = 58 ± 15 years; six with and six without previous fragility fractures) were examined with dual-energy X-ray absorptiometry (DXA) and 7-T magnetic resonance imaging (MRI) of the proximal tibia. A transverse high-resolution three-dimensional fast low-angle shot sequence was acquired (0.234 × 0.234 × 1 mm). Digital topological analysis (DTA) was applied to the medial and lateral subchondral bone of the proximal tibia. The following DTA-based bone microarchitecture parameters were assessed: apparent bone volume; trabecular thickness; profile-edge-density (trabecular bone erosion parameter); profile-interior-density (intact trabecular rods parameter); plate-to-rod ratio; and erosion index. We compared femoral neck T-scores and bone microarchitecture parameters between patients with and without fragility fracture. Results There was no statistical significant difference in femoral neck T-scores between individuals with and without fracture (-2.4 ± 0.9 vs. -1.8 ± 0.7, P = 0.282). Apparent bone volume in the medial compartment was lower in patients with previous fragility fracture (0.295 ± 0.022 vs. 0.317 ± 0.009; P = 0.016). Profile-edge-density, a trabecular bone erosion parameter, was higher in patients with previous fragility fracture in the medial (0.008 ± 0.003 vs. 0.005 ± 0.001) and lateral compartment (0.008 ± 0.002 vs. 0.005 ± 0.001); both P = 0.025. Other DTA parameters did not differ between groups. Conclusion 7-T MRI and DTA permit detection of subtle changes in subchondral bone quality when differences in BMD are not evident.

The robustness of T2 value as a trabecular structural index at multiple spatial resolutions of 7 Tesla MRI

PURPOSE

To evaluate the robustness of MR transverse relaxation times of trabecular bone from spin-echo and gradient-echo acquisitions at multiple spatial resolutions of 7 T.

METHODS

The effects of MRI resolutions to T₂ and T2* of trabecular bone were numerically evaluated by Monte Carlo simulations. T₂ , T2*, and trabecular structural indices from multislice multi-echo and UTE acquisitions were measured in defatted human distal femoral condyles on a 7 T scanner. Reference structural indices were extracted from high-resolution microcomputed tomography images. For bovine knee trabecular samples with intact bone marrow, T₂ and T2* were measured by degrading spatial resolutions on a 7 T system.

RESULTS

In the defatted trabecular experiment, both T₂ and T2* values showed strong ( |r| > 0.80) correlations with trabecular spacing and number, at a high spatial resolution of 125 µm³ . The correlations for MR image-segmentation-derived structural indices were significantly degraded ( |r| < 0.50) at spatial resolutions of 250 and 500 µm³ . The correlations for T2* rapidly dropped ( |r| < 0.50) at a spatial resolution of 500 µm³ , whereas those for T₂ remained consistently high ( |r| > 0.85). In the bovine trabecular experiments with intact marrow, low-resolution (approximately 1 mm³ , 2 minutes) T₂ values did not shorten ( |r| > 0.95 with respect to approximately 0.4 mm³ , 11 minutes) and maintained consistent correlations ( |r| > 0.70) with respect to trabecular spacing (turbo spin echo, 22.5 minutes).

CONCLUSION

T₂ measurements of trabeculae at 7 T are robust with degrading spatial resolution and may be preferable in assessing trabecular spacing index with reduced scan time, when high-resolution 3D micro-MRI is difficult to obtain.

Characterization of knee osteoarthritis-related changes in trabecular bone using texture parameters at various levels of spatial resolution-a simulation study

Articular cartilage and subchondral bone are the key tissues in osteoarthritis (OA). The role of the cancellous bone increasingly attracts attention in OA research. Because of its fast adaptation to changes in the loading distribution across joints, its quantification is expected to improve the diagnosis and monitoring of OA. In this study, we simulated OA progression-related changes of trabecular structure in a series of digital bone models and then characterized the potential of texture parameters and bone mineral density (BMD) as surrogate measures to quantify trabecular bone structure. Five texture parameters were studied: entropy, global and local inhomogeneity, anisotropy and variogram slope. Their dependence on OA relevant structural changes was investigated for three spatial resolutions typically used in micro computed tomography (CT; 10 μm), high-resolution peripheral quantitative CT (HR-pQCT) (90 μm) and clinical whole-body CT equipment (250 μm). At all resolutions, OA-related changes in trabecular bone architecture can be quantified using a specific (resolution dependent) combination of three texture parameters. BMD alone is inadequate for this purpose but if available reduces the required texture parameter combination to anisotropy and global inhomogeneity. The results are summarized in a comprehensive analysis guide for the detection of structural changes in OA knees. In conclusion, texture parameters can be used to characterize trabecular bone architecture even at spatial resolutions below the dimensions of a single trabecula and are essential for a detailed classification of relevant OA changes that cannot be achieved with a measurement of BMD alone.

Quantitative phenotyping of bone fracture repair: a review

Fracture repair is a complex process that involves the interaction of numerous molecular factors, cell lineages and tissue types. These biological processes allow for an impressive feat of engineering: an elastic soft callus is progressively replaced by a more rigid and mineralized callus. During this reparative phase, the healing bone is exposed to a risk of re-fracture. Bone volume and bone quality are the two major factors determining the strength of the callus. Although both factors are important, often only bone volume is analyzed and reported in preclinical studies. Recent developments in techniques for examining bone quality in the callus will enable the rapid and detailed analysis of its material properties and its microstructure. This review aims to give an overview of the methods available for quantitatively phenotyping the bone callus in preclinical studies such as Raman spectroscopy, nanoindentation, scanning acoustic microscopy, in vivo micro-computed tomography (micro-CT) and high-resolution micro-CT. Consolidated and emerging experimental methods are described with a focus on their applicability, and with examples of their utilization.

A Digital Model to Simulate Effects of Bone Architecture Variations on Texture at Spatial Resolutions of CT, HR-pQCT, and μCT Scanners

The quantification of changes in the trabecular bone structure induced by musculoskeletal diseases like osteoarthritis, osteoporosis, rheumatoid arthritis, and others by means of a texture analysis is a valuable tool which is expected to improve the diagnosis and monitoring of a disease. The reaction of texture parameters on different alterations in the architecture of the fine trabecular network and inherent imaging factors such as spatial resolution or image noise has to be understood in detail to ensure an accurate and reliable determination of the current bone state. Therefore, a digital model for the quantitative analysis of cancellous bone structures was developed. Five parameters were used for texture analysis: entropy, global and local inhomogeneity, local anisotropy, and variogram slope. Various generic structural changes of cancellous bone were simulated for different spatial resolutions. Additionally, the dependence of the texture parameters on tissue mineralization and noise was investigated. The present work explains changes in texture parameter outcomes based on structural changes originating from structure modifications and reveals that a texture analysis could provide useful information for a trabecular bone analysis even at resolutions below the dimensions of single trabeculae.

Prediction of bone strength by μCT and MDCT-based finite-element-models: how much spatial resolution is needed?

OBJECTIVES

Finite-element-models (FEM) are a promising technology to predict bone strength and fracture risk. Usually, the highest spatial resolution technically available is used, but this requires excessive computation time and memory in numerical simulations of large volumes. Thus, FEM were compared at decreasing resolutions with respect to local strain distribution and prediction of failure load to (1) validate MDCT-based FEM and to (2) optimize spatial resolution to save computation time.

MATERIALS AND METHODS

20 cylindrical trabecular bone specimens (diameter 12 mm, length 15-20mm) were harvested from elderly formalin-fixed human thoracic spines. All specimens were examined by micro-CT (isotropic resolution 30 μm) and whole-body multi-row-detector computed tomography (MDCT, 250 μm × 250 μm × 500 μm). The resolution of all datasets was lowered in eight steps to ~ 2,000 μm × 2000 μm × 500 μm and FEM were calculated at all resolutions. Failure load was determined by biomechanical testing. Probability density functions of local micro-strains were compared in all datasets and correlations between FEM-based and biomechanically measured failure loads were determined.

RESULTS

The distribution of local micro-strains was similar for micro-CT and MDCT at comparable resolutions and showed a shift toward higher average values with decreasing resolution, corresponding to the increasing apparent trabecular thickness. Small micro-strains (εeff<0.005) could be calculated down to 250 μm × 250 μm × 500 μm. Biomechanically determined failure load showed significant correlations with all FEM, up to r=0.85 and did not significantly change with lower resolution but decreased with high thresholds, due to loss of trabecular connectivity.

CONCLUSION

When choosing connectivity-preserving thresholds, both micro-CT- and MDCT-based finite-element-models well predicted failure load and still accurately revealed the distribution of local micro-strains in spatial resolutions, available in vivo (250 μm × 250 μm × 500 μm), that thus seemed to be the optimal compromise between high accuracy and low computation time.

Cylinders or walls? A new computational model to estimate the MR transverse relaxation rate dependence on trabecular bone architecture

OBJECTIVE

Bone density is distributed in a complex network of interconnecting trabecular plates and rods that are interspersed with bone marrow. A computational model to assess the dependence of the relaxation rate on the geometry of bone can consider the distribution of bone material in the form of two components: cylinders and open walls (walls with gaps). We investigate whether the experimentally known dependence of the transverse relaxation rate on the trabecular bone structure can be usefully interpreted in terms of these two components.

MATERIALS AND METHODS

We established a computer model based on an elementary computational cell. The model includes a variable number of open walls and infinitely long cylinders as well as multiple geometric parameters. The transverse relaxation rate is computed as a function of these parameters. Within the model, increasing the trabecular spacing with a fixed trabecular radius is equivalent to thinning the trabeculae while maintaining constant spacing.

RESULTS

Increasing the number of cylinder and wall gap elements beyond their nearest neighbors does not change the transverse relaxation rate. Although the absolute contribution to the relaxation due to open walls is on average more important than that due to cylinders, the latter drops off rapidly. The change on transverse relaxation rate is larger for changing cylinder geometry than for changing wall geometry, as it can be seen from the effect on the relaxation rate when trabecular spacing is varied, compared to varying the size of wall gaps.

CONCLUSION

Our results provide strong evidence that trabecular thinning, which is associated with increasing age, decreases the relaxation rates. The effect of thinning plates and rods on the transverse relaxation can be understood in terms of simple cylinders and open walls. A reduction in the relaxation rate can be seen as an indication of thinning cylinders, corresponding to reduced bone stability and ultimately, osteoporosis.

Osteochondral repair: evaluation with sweep imaging with fourier transform in an equine model

PURPOSE

To evaluate the status of articular cartilage and bone in an equine model of spontaneous repair by using the sweep imaging with Fourier transform (SWIFT) magnetic resonance (MR) imaging technique.

MATERIALS AND METHODS

Experiments were approved by the Utrecht University Animal Ethics Committee. Six-millimeter-diameter chondral (n = 5) and osteochondral (n = 5, 3-4 mm deep into subchondral bone) defects were created in the intercarpal joints of seven 2-year-old horses and examined with SWIFT at 9.4 T after spontaneous healing for 12 months. Conventional T2 maps and gradient-echo images were obtained for comparison, and histologic assessment of cartilage and micro-computed tomography (CT) of bone were performed for reference. Signal-to-noise ratio (SNR) analysis was performed, and a radiologist evaluated the MR images. Structural bone parameters were derived from SWIFT and micro-CT datasets. Significance of differences was investigated with the Wilcoxon signed rank test and Pearson correlation analysis.

RESULTS

SWIFT was able to depict the different outcomes of spontaneous healing of focal chondral versus osteochondral defects. SWIFT produced constant signal intensity throughout cartilage, whereas T2 mapping showed elevated T2 values (P = .06) in repair tissue (mean T2 in superficial region of interest in an osteochondral lesion = 50.0 msec ± 10.2) in comparison to adjacent intact cartilage (mean T2 = 32.7 msec ± 4.2). The relative SNR in the subchondral plate with SWIFT (0.91) was more than four times higher than that with conventional fast spin-echo (0.12) and gradient-echo (0.19) MR imaging. The correlation between bone volume-to-tissue volume fractions determined with SWIFT and micro-CT was significant (r = 0.83, P < .01).

CONCLUSION

SWIFT enabled assessment of spontaneous osteochondral repair in an equine model.

Micro-finite element analysis applied to high-resolution MRI reveals improved bone mechanical competence in the distal femur of female pre-professional dancers

Micro-finite element analysis applied to high-resolution (0.234-mm length scale) MRI reveals greater whole and cancellous bone stiffness, but not greater cortical bone stiffness, in the distal femur of female dancers compared to controls. Greater whole bone stiffness appears to be mediated by cancellous, rather than cortical bone adaptation.

INTRODUCTION

The purpose of this study was to compare bone mechanical competence (stiffness) in the distal femur of female dancers compared to healthy, relatively inactive female controls.

METHODS

This study had institutional review board approval. We recruited nine female modern dancers (25.7±5.8 years, 1.63±0.06 m, 57.1±4.6 kg) and ten relatively inactive, healthy female controls matched for age, height, and weight (32.1±4.8 years, 1.6±0.04 m, 55.8±5.9 kg). We scanned the distal femur using a 7-T MRI scanner and a three-dimensional fast low-angle shot sequence (TR/TE=31 ms/5.1 ms, 0.234 mm×0.234 mm×1 mm, 80 slices). We applied micro-finite element analysis to 10-mm-thick volumes of interest at the distal femoral diaphysis, metaphysis, and epiphysis to compute stiffness and cross-sectional area of whole, cortical, and cancellous bone, as well as cortical thickness. We applied two-tailed t-tests and ANCOVA to compare groups.

RESULTS

Dancers demonstrated greater whole and cancellous bone stiffness and cross-sectional area at all locations (p<0.05). Cortical bone stiffness, cross-sectional area, and thickness did not differ between groups (>0.08). At all locations, the percent of intact whole bone stiffness for cortical bone alone was lower in dancers (p<0.05). Adjustment for cancellous bone cross-sectional area eliminated significant differences in whole bone stiffness between groups (p>0.07), but adjustment for cortical bone cross-sectional area did not (p<0.03).

CONCLUSIONS

Modern dancers have greater whole and cancellous bone stiffness in the distal femur compared to controls. Elevated whole bone stiffness in dancers may be mediated via cancellous, rather than cortical bone adaptation.

In vivo estimation of bone stiffness at the distal femur and proximal tibia using ultra-high-field 7-Tesla magnetic resonance imaging and micro-finite element analysis

The goal of this study was to demonstrate the feasibility of using 7-Tesla (7T) magnetic resonance imaging (MRI) and micro-finite element analysis (µFEA) to evaluate mechanical and structural properties of whole, cortical, and trabecular bone at the distal femur and proximal tibia in vivo. 14 healthy subjects were recruited (age 40.7 ± 15.7 years). The right knee was scanned on a 7T MRI scanner using a 28 channel-receive knee coil and a three-dimensional fast low-angle shot sequence (TR/TE 20 ms/5.02 ms, 0.234 mm × 0.234 mm × 1 mm, 80 axial images, 7 min 9 s). Bone was analyzed at the distal femoral metaphysis, femoral condyles, and tibial plateau. Whole, cortical, and trabecular bone stiffness was computed using µFEA. Bone volume fraction (BVF), bone areas, and cortical thickness were measured. Trabecular bone stiffness (933.7 ± 433.3 MPa) was greater than cortical bone stiffness (216 ± 152 MPa) at all three locations (P < 0.05). Across locations, there were no differences in bone stiffness (whole, cortical, or trabecular). Whole, cortical, and trabecular bone stiffness correlated with BVF (R ≥ 0.69, P < 0.05) and inversely correlated with corresponding whole, cortical, and trabecular areas (R ≤ -0.54, P < 0.05), but not with cortical thickness (R < -0.11, P > 0.05). Whole, cortical, and trabecular stiffness correlated with body mass index (R ≥ 0.62, P < 0.05). In conclusion, at the distal femur and proximal tibia, trabecular bone contributes 66-74% of whole bone stiffness. 7T MRI and µFEA may be used as a method to provide insight into how structural properties of cortical or trabecular bone affect bone mechanical competence in vivo.