Open-source pipeline for automatic segmentation and microstructural analysis of murine knee subchondral bone

μCT images are commonly analysed to assess changes in bone density and microstructure in preclinical murine models. Several platforms provide automated analysis of bone microstructural parameters from volumetric regions of interest (ROI). However, segmentation of the regions of subchondral bone to create the volumetric ROIs remains a manual and time-consuming task. This study aimed to develop an automated end-to-end pipeline, combining segmentation and microstructural analysis, to evaluate subchondral bone in the mouse proximal knee.

A segmented dataset of μCT scans from 62 knees (healthy and arthritic) from 10-week male C57BL/6 mice was used to train a U-Net type architecture to automate segmentation of the subchondral trabecular bone. These segmentations were used in tandem with the original scans as input for microstructural analysis along with thresholded trabecular bone. Manually and U-Net segmented ROIs were fed into two available pipelines for microstructural analysis: the ITKBoneMorphometry library and CTan (SKYSCAN). Outcome parameters were compared between pipelines, including: bone volume (BV), total volume (TV), BV/TV, trabecular number (TbN), trabecular thickness (TbTh), trabecular separation (TbSp), and bone surface density (BSBV).

There was good agreement for all bone measures comparing the manual and U-Net pipelines utilizing ITK (R = 0.88-0.98) and CTAn (R = 0.91-0.98). ITK and CTAn showed good agreement for BV, TV, BV/TV, TbTh and BSBV (R = 0.9-0.98). However, limited agreement was seen between TbN (R = 0.73) and TbSb (R = 0.59) due to methodological differences in how spacing is evaluated. Microstructural parameters generated from manual and automatic segmentations showed high correlation across all measures. Using the CTAn pipeline yielded strong R² values (0.83-0.96) and very strong agreement based on ICC (0.90-0.98). The ITK pipeline yielded similarly high R² values (0.91-0.96, except for TbN (0.77)), and ICC values (0.88-0.98). The automated segmentations yield lower average values for BV, TV and BV/TV (ranging from 14 % to 6.3 %), but differences were not found to be influenced by the mean ROI values.

This integrated pipeline seamlessly automated both segmentation and quantification of the proximal tibia subchondral bone microstructure. This automated pipeline allows the analysis of large volumes of data, and its open-source nature may enable the standardization of microstructural analysis of trabecular bone across different research groups.