Quantifying Complex Micro-Topography of Degenerated Articular Cartilage Surface by Contrast-Enhanced Micro-Computed Tomography and Parametric Analyses

Modulation of early osteoarthritis by tibiofemoral re-alignment in sheep

OBJECTIVE

To investigate whether tibiofemoral alignment influences early knee osteoarthritis (OA). We hypothesized that varus overload exacerbates early degenerative osteochondral changes, and that valgus underload diminishes early OA.

METHOD

Normal, over- and underload were induced by altering alignment via high tibial osteotomy in adult sheep (n = 8 each). Simultaneously, OA was induced by partial medial anterior meniscectomy. At 6 weeks postoperatively, OA was examined in five individual subregions of the medial tibial plateau using Kellgren-Lawrence grading, quantification of macroscopic OA, semiquantitative histopathological OA and immunohistochemical type-II collagen, ADAMTS-5, and MMP-13 scoring, biochemical determination of DNA and proteoglycan contents, and micro-computed tomographic evaluation of the subchondral bone.

RESULTS

Multivariate analyses revealed that OA cartilaginous changes had a temporal priority over subchondral bone changes. Underload inhibited early cartilage degeneration in a characteristic topographic pattern (P ≥ 0.0983 vs. normal), in particular below the meniscal damage, avoided alterations of the subarticular spongiosa (P ≥ 0.162 vs. normal), and prevented the disturbance of otherwise normal osteochondral correlations. Overload induced early alterations of the subchondral bone plate microstructure towards osteopenia, including significantly decreased percent bone volume and increased bone surface-to-volume ratio (all P ≤ 0.0359 vs. normal).

CONCLUSION

The data provide high-resolution evidence that tibiofemoral alignment modulates early OA induced by a medial meniscus injury in adult sheep. Since underload inhibits early OA, these data also support the clinical value of strategies to reduce the load in an affected knee compartment to possibly decelerate structural OA progression.

Objective analysis of facial bone fracture CT images using curvature measurement in a surface mesh model

The human facial skeleton consists of multiple segments and causes difficulty during analytic processes. We developed image analysis software to quantify the amount of injury and validate the smooth curvature of the surface after facial bone reduction surgery. Three-dimensional computed tomography images of facial bone were obtained from 40 patients who had undergone open reduction surgery to treat unilateral zygomaticomaxillary fractures. Analytic software was developed based on the discrete curvature of a triangular mesh model. The discrete curvature values were compared before and after surgery using two regions of interest. For the inferior orbital rim, the weighted average of curvature changed from 0.543 ± 0.034 to 0.458 ± 0.042. For the anterior maxilla, the weighted average of curvature changed from 0.596 ± 0.02 to 0.481 ± 0.031, showing a significant decrement (P < 0.05). The curvature was further compared with the unaffected side using the Bray-Curtis similarity index (BCSI). The BCSI of the inferior orbital rim changed from 0.802 ± 0.041 to 0.904 ± 0.015, and that for the anterior maxilla changed from 0.797 ± 0.029 to 0.84 ± 0.025, demonstrating increased similarity (P < 0.05). In computational biology, adequate analytic software is crucial. The newly developed software demonstrated significant differentiation between pre- and postoperative curvature values. Modification of formulas and software will lead to further advancements.

Relationships between tibial articular cartilage, in vivo external joint moments and static alignment in end-stage knee osteoarthritis: A micro-CT study

Biomechanical factors (e.g., joint loading) have a significant role in the progression of osteoarthritis (OA). However, some relationships between in vivo joint loading indices and tibial cartilage thickness are conflicting. This study investigated relationships between pre-operative in vivo external knee joint moments, joint alignment and regional tibial cartilage thickness using micro-CT in subjects with end-stage knee OA. Tibial plateaus from 25 patients that underwent knee replacement for OA were micro-CT scanned (17 µm/voxel). Prior to surgery, subjects underwent gait analysis to calculate external knee moments. The mechanical axis deviation (MAD) was obtained from pre-operative radiographs. Cartilage thickness (Cart.Th) was analyzed from micro-CT images, in anteromedial, anterolateral, posteromedial and posterolateral subregions of interest. Medial-to-lateral Cart.Th ratios were also explored. Relationships between Cart.Th and joint loading indices were examined using Pearson's correlations. Significant correlations were found between Cart.Th and joint loading indices, positive anteromedially with the first peak knee adduction moment (r = 0.55, p < 0.01) and external rotation moment (ERM; r = 0.52, p < 0.01), and negative with MAD (r = -0.76, p < 0.001). In the lateral regions, these correlations had opposite signs. The medial-to-lateral Cart.Th ratio correlated strongly with ERM (r = 0.63, p = 0.001) and MAD (r = -0.75, p < 0.001). Joint loading indices correlated with regional cartilage thickness values and their medial-to-lateral ratios in end-stage knee OA subjects, with higher regional loads corresponding to thinner cartilage. These relationships have the opposite sign compared to the subchondral bone microarchitecture found in our previous study on the same specimens, which may suggest a complementary bone-cartilage interplay in response to loading.

Tibial cartilage, subchondral bone plate and trabecular bone microarchitecture in varus- and valgus-osteoarthritis versus controls

This preliminary study quantified tibia cartilage thickness (Cart.Th), subchondral bone plate thickness (SBPl.Th) and subchondral trabecular bone (STB) microarchitecture in subjects with varus- or valgus- malaligned knees diagnosed with end-stage knee osteoarthritis (OA) and compared them to controls (non-OA). Tibial plateaus from 25 subjects with knee-OA (undergoing knee arthroplasty) and 15 cadavers (controls) were micro-CT scanned (17 µm/voxel). Joint alignment was classified radiographically for OA subjects (varus-aligned n = 18, valgus-aligned n = 7). Cart.Th, SBPl.Th, STB bone volume fraction (BV/TV) and their medial-to-lateral ratios were analyzed in anteromedial, anterolateral, posteromedial and posterolateral subregions. Varus-OA and valgus-OA were compared to controls. Compared to controls (1.19-1.54 mm), Cart.Th in varus-OA was significantly lower anteromedially (0.58 mm, -59%) and higher laterally (2.19-2.47 mm, +60-63%); in valgus-OA, Cart.Th was significantly higher posteromedially (1.86 mm, +56%). Control medial-to-lateral Cart.Th ratios were around unity (0.8-1.1), in varus-OA significantly below (0.2-0.6) and in valgus-OA slightly above (1.0-1.3) controls. SBPl.Th and BV/TV were significantly higher medially in varus-OA (0.58-0.72 mm and 37-44%, respectively) and laterally in valgus-OA (0.60-0.61 mm and 32-37%), compared to controls (0.26-0.47 mm and 18-37%). In varus-OA, the medial-to-lateral SBPl.Th and BV/TV ratios were above unity (1.4-2.4) and controls (0.8-2.1); in valgus-OA they were closer to unity (0.8-1.1) and below controls. Varus- and valgus-OA tibia differ significantly from controls in Cart.Th, SBPl.Th and STB microarchitecture depending on joint alignment, suggesting structural changes in OA may reflect differences in medial-to-lateral load distribution upon the tibial plateau. Here we identified an inverse relationship between cartilage thickness and underlying subchondral bone, suggesting a whole-joint response in OA to daily stimuli.

Automating three-dimensional osteoarthritis histopathological grading of human osteochondral tissue using machine learning on contrast-enhanced micro-computed tomography

OBJECTIVE

To develop and validate a machine learning (ML) approach for automatic three-dimensional (3D) histopathological grading of osteochondral samples imaged with contrast-enhanced micro-computed tomography (CEμCT).

DESIGN

A total of 79 osteochondral cores from 24 total knee arthroplasty patients and two asymptomatic donors were imaged using CEμCT with phosphotungstic acid -staining. Volumes-of-interest (VOI) in surface (SZ), deep (DZ) and calcified (CZ) zones were extracted depth-wise and subjected to dimensionally reduced Local Binary Pattern -textural feature analysis. Regularized linear and logistic regression (LR) models were trained zone-wise against the manually assessed semi-quantitative histopathological CEμCT grades (diameter = 2 mm samples). Models were validated using nested leave-one-out cross-validation and an independent test set (4 mm samples). The performance was primarily assessed using Mean Squared Error (MSE) and Average Precision (AP, confidence intervals are given in square brackets).

RESULTS

Highest performance on cross-validation was observed for SZ, both on linear regression (MSE = 0.49, 0.69 and 0.71 for SZ, DZ and CZ, respectively) and LR (AP = 0.9 [0.77-0.99], 0.46 [0.28-0.67] and 0.65 [0.41-0.85] for SZ, DZ and CZ, respectively). The test set evaluations yielded increased MSE on all zones. For LR, the performance was also best for the SZ (AP = 0.85 [0.73-0.93], 0.82 [0.70-0.92] and 0.8 [0.67-0.9], for SZ, DZ and CZ, respectively).

CONCLUSION

We present the first ML-based automatic 3D histopathological osteoarthritis (OA) grading method which also adequately perform on grading unseen data, especially in SZ. After further development, the method could potentially be applied by OA researchers since the grading software and all source codes are publicly available.

Full-Field Strain Uncertainties and Residuals at the Cartilage-Bone Interface in Unstained Tissues Using Propagation-Based Phase-Contrast XCT and Digital Volume Correlation

A deeper understanding of the cartilage-bone mechanics is fundamental to unravel onset and progression of osteoarthritis, enabling better diagnosis and treatment. The aim of this study is therefore to explore the capability of X-ray computed (XCT) phase-contrast imaging in a lab-based system to enable digital volume correlation (DVC) measurements of unstained cartilage-bone plugs from healthy adult bovines. DVC strain uncertainties were computed for both articular cartilage and mineralized tissue (calcified cartilage and subchondral bone) in the specimens at increasing propagation distances, ranging from absorption up to four times (4× such effective distance. In addition, a process of dehydration and rehydration was proposed to improve feature recognition in XCT of articular cartilage and mechanical properties of this tissue during the process were assessed via micromechanical probing (indentation), which was also used to determine the effect of long X-ray exposure. Finally, full-field strain from DVC was computed to quantify residual strain distribution at the cartilage-bone interface following unconfined compression test (ex situ). It was found that enhanced gray-scale feature recognition at the cartilage-bone interface was achieved using phase-contrast, resulting in reduced DVC strain uncertainties compared to absorption. Residual strains up to ~7000 µε in the articular cartilage were transferred to subchondral bone via the calcified cartilage and micromechanics revealed the predominant effect of long phase-contrast X-ray exposure in reducing both stiffness and hardness of the articular cartilage. The results of this study will pave the way for further development and refinement of the techniques, improving XCT-based strain measurements in cartilage-bone and other soft-hard tissue interfaces.