Adaptive local thresholding can enhance the accuracy of HR-pQCT-based trabecular bone morphology assessment

Accuracy and precision of segmentation and quantification of wrist bone microarchitecture using photon-counting computed tomography ex vivo

The quantification of bone microarchitecture provides insight into bone health and the effects of disease or treatment, and is therefore highly relevant clinical information. Nonetheless, in vivo quantification of bone microarchitecture is mostly limited to high-resolution peripheral quantitative CT (HR-pQCT). This is a small field of view CT modality of which the gantry size only allows scanning of distal radius and tibia. Photon-counting CT (PCCT) is a novel clinical full-body CT with improved image resolution and quality compared to other clinical CT modalities, yet data on its capabilities in quantifying bone microarchitecture are limited. The aim of this study was to quantify the accuracy of two methods for trabecular bone segmentation on PCCT images as compared to the segmentations on micro-CT (μCT) and to use these segmentations to quantify the accuracy and agreement of trabecular bone morphometry measurements as compared to μCT, as well as the short-term precision. This study analysed multimodal CT data, obtained from eight cadaveric forearms; the data includes two repeated PCCT scans, as well as a single HR-pQCT scan from the forearm, and μCT scans of all individual carpal bones. For each carpal bone, trabecular volumes of interest (VOI) were delineated on the μCT images, and the μCT reference segmentations and VOIs were resampled onto the PCCT and HR-pQCT images. HR-pQCT images were segmented with a global threshold of 320 mgHA/cm3; PCCT images were segmented with either an identical global threshold or with an adaptive thresholding algorithm. Trabecular bone-volume fraction (Tb.BV/TV), trabecular thickness (Tb.Th), trabecular number (Tb.N) and trabecular separation (Tb.Sp) were quantified for all segmented VOIs. Accuracy and agreement were calculated relative to μCT as the gold standard, short-term precision was calculated from the repeated PCCT scan. For PCCT, adaptive threshold segmentation had significantly increased sensitivity compared to global threshold segmentation, along with a lower variance in its sensitivity and specificity. Concerning the microarchitecture quantification, for global threshold segmentation of PCCT images, correlations with μCT were significant, except for Tb.Sp. Correlation coefficients of Tb.BV/TV and Tb.N were not significantly different from those between HR-pQCT and μCT. Adaptive threshold segmentation led to higher correlation coefficients between PCCT and μCT of Tb.Th, Tb.N and Tb.Sp, although correlations of Tb.N remained poor for both PCCT and HR-pQCT. Moreover, adaptive thresholding led to a constant bias of Tb.BV/TV, Tb.Th and Tb.Sp, unlike the bias of HR-pQCT which was proportionally increasing with the size of the measurement. Finally, adaptive threshold segmentation led to a higher short-term precision than global threshold segmentation, with a root-mean-squared coefficient of variation below 0.65 % for all parameters. We conclude that adaptive threshold segmentation is well-suited for the segmentation of PCCT images. Despite measurement error, our results indicate that these segmentations can be used for bone microarchitecture analyses of carpal bones with agreement and short-term precision comparable to HR-pQCT.

Biomechanics of the Human Osteochondral Unit: A Systematic Review

The damping system ensured by the osteochondral (OC) unit is essential to deploy the forces generated within load-bearing joints during locomotion, allowing furthermore low-friction sliding motion between bone segments. The OC unit is a multi-layer structure including articular cartilage, as well as subchondral and trabecular bone. The interplay between the OC tissues is essential in maintaining the joint functionality; altered loading patterns can trigger biological processes that could lead to degenerative joint diseases like osteoarthritis. Currently, no effective treatments are available to avoid degeneration beyond tissues' recovery capabilities. A thorough comprehension on the mechanical behaviour of the OC unit is essential to (i) soundly elucidate its overall response to intra-articular loads for developing diagnostic tools capable of detecting non-physiological strain levels, (ii) properly evaluate the efficacy of innovative treatments in restoring physiological strain levels, and (iii) optimize regenerative medicine approaches as potential and less-invasive alternatives to arthroplasty when irreversible damage has occurred. Therefore, the leading aim of this review was to provide an overview of the state-of-the-art-up to 2022-about the mechanical behaviour of the OC unit. A systematic search is performed, according to PRISMA standards, by focusing on studies that experimentally assess the human lower-limb joints' OC tissues. A multi-criteria decision-making method is proposed to quantitatively evaluate eligible studies, in order to highlight only the insights retrieved through sound and robust approaches. This review revealed that studies on human lower limbs are focusing on the knee and articular cartilage, while hip and trabecular bone studies are declining, and the ankle and subchondral bone are poorly investigated. Compression and indentation are the most common experimental techniques studying the mechanical behaviour of the OC tissues, with indentation also being able to provide information at the micro- and nanoscales. While a certain comparability among studies was highlighted, none of the identified testing protocols are currently recognised as standard for any of the OC tissues. The fibril-network-reinforced poro-viscoelastic constitutive model has become common for describing the response of the articular cartilage, while the models describing the mechanical behaviour of mineralised tissues are usually simpler (i.e., linear elastic, elasto-plastic). Most advanced studies have tested and modelled multiple tissues of the same OC unit but have done so individually rather than through integrated approaches. Therefore, efforts should be made in simultaneously evaluating the comprehensive response of the OC unit to intra-articular loads and the interplay between the OC tissues. In this regard, a multidisciplinary approach combining complementary techniques, e.g., full-field imaging, mechanical testing, and computational approaches, should be implemented and validated. Furthermore, the next challenge entails transferring this assessment to a non-invasive approach, allowing its application in vivo, in order to increase its diagnostic and prognostic potential.

Bone remodeling and responsiveness to mechanical stimuli in individuals with type 1 diabetes mellitus

Type 1 diabetes mellitus (T1DM) has been linked to increased osteocyte apoptosis, local accumulation of mineralized lacunar spaces, and microdamage suggesting an impairment of the mechanoregulation network in affected individuals. Diabetic neuropathy might exacerbate this dysfunction through direct effects on bone turnover, and indirect effects on balance, muscle strength, and gait. However, the in vivo effects of impaired bone mechanoregulation on bone remodeling in humans remain underexplored. This longitudinal cohort study assessed consenting participants with T1DM and varying degree of distal symmetric sensorimotor polyneuropathy (T1DM, n = 20, median age 46.5 yr, eight female) and controls (CTRL; n = 9, median age 59.0 yr, four female) at baseline and 4-yr follow-up. Nerve conduction in participants with T1DM was tested using DPNCheck and bone remodeling was quantified with longitudinal high-resolution peripheral quantitative-computed tomography (HR-pQCT, 82 μm) at the standard distal sites. Local trabecular bone formation (Tb.F) and resorption (Tb.R) sites were captured by implementing 3D rigid image registration of HR-pQCT images, and the mechanical environment across the bone microarchitecture at these sites was simulated using micro-finite element analysis. We calculated odds ratios to determine the likelihood of bone formation (ORF) and resorption (ORR) with increasing/decreasing strain in percent as markers for mechanoregulation. At the distal radius, Tb.F was 47% lower and Tb.R was 59% lower in T1DM participants compared with CTRL (P < .05). Tb.F correlated positively with nerve conduction amplitude (R = 0.69, P < .05) in participants with T1DM and negatively with glycated hemoglobin (HbA1c) (R = -0.45, P < .05). Additionally, ORF was 34% lower and ORR was 18% lower in T1DM compared with CTRL (P < .05). Our findings represent in vivo evidence suggesting that bone remodeling in individuals with T1DM is in a state of low responsiveness to mechanical stimuli, resulting in impaired bone formation and resorption rates; these correlate to the degree of neuropathy and level of diabetes control.

A Laplace-Hamming Binarization Approach for Second-Generation HR-pQCT Rescues Fine Feature Segmentation

Although second-generation high-resolution peripheral quantitative computed tomography (XCTII) provides the highest-resolution in vivo bone microstructure assessment, the manufacturer's standard image processing protocol omits fine features in both trabecular and cortical compartments. To optimize fine structure segmentation, we implemented a binarization approach based on a Laplace-Hamming (LH) segmentation and documented the reproducibility and accuracy of XCTII structure segmentation using both the standard Gaussian-based binarization and the proposed LH segmentation approach. To evaluate reproducibility, 20 volunteers (9 women, 11 men; aged 23-75 years) were recruited, and three repeat scans of the radii and tibias were acquired using the manufacturer's standard in vivo protocol. To evaluate accuracy, cadaveric structure phantoms (14 radii, 6 tibias) were scanned on XCTII using the same standard in vivo protocol and on μCT at 24.5 μm resolution. XCTII images were analyzed twice-first, with the manufacturer's standard patient evaluation protocol and, second, with the proposed LH segmentation approach. The LH approach rescued fine features evident in the grayscale images but omitted or overrepresented (thickened) by the standard approach. The LH approach significantly reduced error in trabecular volume fraction (BV/TV) and thickness (Tb.Th) compared with the standard approach; however, higher error was introduced for trabecular separation (Tb.Sp). The LH approach improved the correlation between XCTII and μCT for cortical porosity (Ct.Po) and significantly reduced error in cortical pore diameter (Ct.Po.Dm) compared with the standard approach. The LH approach resulted in improved precision compared with the standard approach for BV/TV, Tb.Th, Ct.Po, and Ct.Po.Dm at the radius and for Ct.Po at the tibia. Our results suggest that the proposed LH approach produces substantially improved binary masks, reduces proportional bias, and provides greater accuracy and reproducibility in important outcome metrics, all due to more accurate segmentation of the fine features in both trabecular and cortical compartments. © 2023 The Authors. Journal of Bone and Mineral Research published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research (ASBMR).

Association between radiographic hand osteoarthritis and bone microarchitecture in a population-based sample

Background

Subchondral bone plays an important role in the pathogenesis of radiographic osteoarthritis (OA). However, the bony changes that occur in hand OA (HOA) are much less understood. This study aimed to describe the association between radiographic HOA and high-resolution peripheral quantitative computed tomography (HRpQCT) measures of the hand and radius in a population-based sample.

Methods

A total of 201 participants (mean age 72, 46% female) from the Tasmanian Older Adult Cohort (TASOAC) study underwent HRpQCT assessment of the 2nd distal and proximal interphalangeal (DIP, PIP), 1st carpometacarpal (CMC) joint, and distal radius. Radiographic HOA was assessed at the 2nd DIP, PIP joints, and the 1st CMC joint using the OARSI atlas.

Results

Proximal osteophyte and joint space narrowing (JSN) scores were consistently more strongly associated with HRpQCT measures compared to the distal site with positive associations for indices of bone size (total and trabecular bone area and cortical perimeter but inconsistent for cortical area) and negative associations for volumetric bone mineral density (vBMD). There was a decrease in trabecular number and bone volume fraction with increasing osteophyte and JSN score as well as an increase in trabecular separation and inhomogeneity. Osteophyte and JSN scores in the hand were not associated with HRpQCT measures at the distal radius.

Conclusions

This hypothesis generating data suggests that bone size and trabecular disorganization increase with both osteophyte formation and JSN (proximal more than distal), while local vBMD decreases. This process appears to be primarily at the site of pathology rather than nearby unaffected bone.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13075-022-02907-6.