Local plate/rod descriptors of 3D trabecular bone micro-CT images from medial axis topologic analysis

Establishment and Characteristic Analysis of a Dog Model for Autologous Homologous Cranioplasty

Objective

The aim of this study is to establish a large animal (dog) model that can be referred clinically for autologous homologous cranioplasty.

Methods

Our large skull defect dog model was established by emulating the decompressive craniectomy with 22 adult beagle dogs. The autologous bones were taken out from the dogs and divided into two groups, the freeze-drying (FD) group and the single freezing (SF) group. They were then stored in the bone bank at -20°C after being irradiated with 25 KGy. Three months later, the bones were reimplanted. After operation, we closely watch the experimental objects for four more months examining the infection and survival of the bone graft.

Results

Through macroscopic observation, it was found that, among 44 cranial flaps (bilateral) from the rest of the 22 dogs, grade A cranial flaps accounted for 86.4% (19/22) in the SF group and only 31.8% (7/22) in the FD group. Although osteogenic osteoclast, Harvard tube, neovascularization, and angiogenic factors were found through the pathological results, including an electron microscope and calmodulin tracer, it could be verified by using X-CT and micro-CT that early bone resorption could be still found even in grade A bone flap.

Conclusion

By using the common clinical method to preserve the cranial flaps, we established an experimental dog model of autologous cranioplasty for a large area of cranial defect. It was proved that this model could reproduce the infections and bone resorption which typically happened in clinical autologous homologous cranioplasty. As a conclusion, the established model can be used as an effective experimental tool for further research to improve the success rate of autologous homologous cranioplasty.

Diffusion Tensor Imaging of Lumbar Vertebras in Female Adolescent Idiopathic Scoliosis: Initial Findings

OBJECTIVE

The purpose of this study was to characterize diffusion tensor imaging (DTI) features of lumbar vertebras in adolescent idiopathic scoliosis (AIS) patients.

METHODS

Fifty-two AIS patients and 20 healthy volunteers underwent 3-T magnetic resonance scanning including DTI sequence. The fractional anisotropy (FA) and apparent diffusion coefficient (ADC) values on the convex and concave sides of lumbar vertebras were obtained and compared.

RESULTS

The FA and ADC values differed significantly between the convex and concave side of lumbar vertebras in AIS (P < 0.01). The ADC values in AIS differed significantly with healthy volunteers (P < 0.01). The FA values on the convex side of L1 to L2 were significantly lower than L4 to L5 in AIS. The difference of FA values between the concave and convex sides of the apex vertebra correlated significantly with Cobb angle (r = 0.436, P < 0.01).

CONCLUSIONS

The convex and concave sides of lumbar vertebras in AIS patients showed different DTI features.

Vascular imaging with contrast agent in hard and soft tissues using microcomputed-tomography

Vascularization is essential for many tissues and is a main requisite for various tissue-engineering strategies. Different techniques are used for highlighting vasculature, in vivo and ex vivo, in 2-D or 3-D including histological staining, immunohistochemistry, radiography, angiography, microscopy, computed tomography (CT) or micro-CT, both stand-alone and synchrotron system. Vascularization can be studied with or without a contrast agent. This paper presents the results obtained with the latest Skyscan micro-CT (Skyscan 1272, Bruker, Belgium) following barium sulphate injection replacing the bloodstream in comparison with results obtained with a Skyscan In Vivo 1076. Different hard and soft tissues were perfused with contrast agent and were harvested. Samples were analysed using both forms of micro-CT, and improved results were shown using this new micro-CT. This study highlights the vasculature using micro-CT methods. The results obtained with the Skyscan 1272 are clearly defined compared to results obtained with Skyscan 1076. In particular, this instrument highlights the high number of small vessels, which were not seen before at lower resolution. This new micro-CT opens broader possibilities in detection and characterization of the 3-D vascular tree to assess vascular tissue engineering strategies.

The Ellipsoid Factor for Quantification of Rods, Plates, and Intermediate Forms in 3D Geometries

The ellipsoid factor (EF) is a method for the local determination of the rod- or plate-like nature of porous or spongy continua. EF at a point within a 3D structure is defined as the difference in axis ratios of the greatest ellipsoid that fits inside the structure and that contains the point of interest, and ranges from -1 for strongly oblate (discus-shaped) ellipsoids, to +1 for strongly prolate (javelin-shaped) ellipsoids. For an ellipsoid with axes a ≤ b ≤ c, EF = a/b - b/c. Here, EF is demonstrated in a Java plugin, "Ellipsoid Factor" for ImageJ, distributed in the BoneJ plugin collection. Ellipsoid Factor utilizes an ellipsoid optimization algorithm, which assumes that maximal ellipsoids are centered on the medial axis, then dilates, rotates, and translates slightly each ellipsoid until it cannot increase in volume any further. EF successfully identifies rods, plates, and intermediate structures within trabecular bone, and summarizes the distribution of geometries with an overall EF mean and SD, EF histogram, and Flinn diagram displaying a/b versus b/c. EF is released to the community for testing, use, and improvement.

Local topological analysis at the distal radius by HR-pQCT: Application to in vivo bone microarchitecture and fracture assessment in the OFELY study

High-resolution peripheral quantitative computed tomography (HR-pQCT) is an in-vivo technique used to analyze the distal radius and tibia. It provides a voxel size of 82μm. In addition to providing the usual microarchitecture parameters, local topological analysis (LTA) depicting rod- and plate-like trabeculae may improve prediction of bone fragility. Thirty-three women with prevalent wrist fractures from the OFELY cohort were compared with age-matched controls. Bone microarchitecture, including the structural model index (SMI), was assessed by HR-pQCT, and micro-finite element analysis (μFE) was computed on trabecular bone images of the distal radius (XtremeCT, Scanco Medical AG). A new LTA method was applied to label each bone voxel as a rod, plate or node. Then the bone volume fraction (BV/TV*), the rod, plate and node ratios over bone volume (RV/BV*, PV/BV*, NV/BV*) or total volume (RV/TV*, PV/TV*, NV/TV*) and the rod to plate ratio (RV/PV*) were calculated. Associations between LTA parameters and wrist fractures were computed in a conditional logistic regression model. Multivariate models were tested to predict the μFE-derived trabecular bone stiffness. RV/TV* (OR=4.41 [1.05-18.62]) and BV/TV* (OR=6.45 [1.06-39.3]), were significantly associated with prevalent wrist fracture, after adjustment for ultra distal radius aBMD. Multivariate linear models including PV/TV* or BV/TV*+RV/PV* predicted trabecular stiffness with the same magnitude as those including SMI. Conversion from plates into rods was significantly associated with bone fragility, with a negative correlation between RV/PV* and trabecular bone stiffness (r=-0.63, p<0.0001). We conclude that our local topological analysis is feasible for a voxel size of 82μm. After further validation, it may improve bone fragility description.

Trabecular architecture and vertebral fragility in osteoporosis

Osteoporosis heightens vertebral fragility owing to the biomechanical effects of diminished bone structure and composition. These biomechanical effects are only partially explained by loss in bone mass, so additional factors that are independent of bone mass are also thought to play an important role in vertebral fragility. Recent advances in imaging equipment, imaging-processing methods, and computational capacity allow researchers to quantify trabecular architecture in the vertebra at the level of the individual trabecular elements and to derive biomechanics-based measures of architecture that are independent of bone mass and density. These advances have shed light on the role of architecture in vertebral fragility. In addition to the adverse biomechanical consequences associated with trabecular thinning and loss of connectivity, a reduction in the number of vertical trabecular plates appears to be particularly harmful to vertebral strength. In the clinic, detailed architecture analysis is primarily applied to peripheral sites such as the distal radius and tibia. Analysis of trabecular architecture at these peripheral sites has shown mixed results for discriminating between patients with and without a vertebral fracture independent of bone mass, but has the potential to provide unique insight into the effects of therapeutic treatments. Overall, it does appear that trabecular architecture has an independent role on vertebral strength. Additional research is required to determine how and where architecture should be measured in vivo and whether assessment of trabecular architecture in a clinical setting improves prospective fracture risk assessment for the vertebra.

3D shape-dependent thinning method for trabecular bone characterization

PURPOSE

Curve and surface thinning are widely-used skeletonization techniques for modeling objects in three dimensions. In the case of trabecular bone analysis, however, neither curve nor surface thinning is really efficient since the internal geometry of the object is usually composed of both rod and plate shapes. The purpose of this paper is to propose an original method called hybrid skeleton which better matches the geometry of the data compared to curve and surface skeletons. In the hybrid skeleton algorithm, 1D curves represent rod-shaped zones whereas 2D surfaces represent plate-shaped elements.

METHODS

The proposed hybrid skeleton algorithm is based on a combination of three methods. (1) A new variant of the method proposed by Bonnassie et al. for the classification of voxels as belonging to plate-like or rod-like structures, where the medial axis (MA) algorithm is replaced by a fast and connected skeletonization algorithm. In addition, the reversibility of the MA algorithm is replaced by an isotropic region-growth method to spread the rod and plate labels back to the original object. (2) A well chosen surface thinning method applied on the plate voxels set. (3) A well chosen curve skeleton thinning method applied on the rod voxels set. The efficiency and the robustness of the proposed algorithm were evaluated using synthesis test vectors. A clinical study was led on micro-CT (computed tomography) images of two different populations of osteoarthritic and osteoporotic trabecular bone samples. The morphological and topological characteristics of the two populations were evaluated using the proposed hybrid skeleton as well as the classification algorithm.

RESULTS

When evaluated on test vectors and compared to Bonnassie's algorithm, the proposed classification algorithm gives a slightly better rate of classification. The hybrid skeleton preserves the shape information of the processed objects. Interesting morphological and topological features as well as volumetric ones were extracted from the skeleton and from the classified volumes, respectively. The extracted features enable the two populations of osteoarthritic and osteoporotic trabecular bone samples to be distinguished.

CONCLUSIONS

Compared to curve-based or surface-based skeletons, the hybrid skeleton better matches the geometry of the data. Each rod is represented by a one-voxel-thick arc and each plate is represented by a one-voxel-thick surface. The hybrid skeleton as well as the proposed classification algorithm introduce relevant parameters linked to the presence of plates in the trabecular bone data, showing that rods and plates contain independent information about trabeculae. The hybrid skeleton offers a new opportunity for precise studies of porous media such as trabecular bone.

High-resolution computed tomography for clinical imaging of bone microarchitecture

BACKGROUND

The role of bone structure, one component of bone quality, has emerged as a contributor to bone strength. The application of high-resolution imaging in evaluating bone structure has evolved from an in vitro technology for small specimens to an emerging clinical research tool for in vivo studies in humans. However, many technical and practical challenges remain to translate these techniques into established clinical outcomes.

QUESTIONS/PURPOSES

We reviewed use of high-resolution CT for evaluating trabecular microarchitecture and cortical ultrastructure of bone specimens ex vivo, extension of these techniques to in vivo human imaging studies, and recent studies involving application of high-resolution CT to characterize bone structure in the context of skeletal disease.

METHODS

We performed the literature review using PubMed and Google Scholar. Keywords included CT, MDCT, micro-CT, high-resolution peripheral CT, bone microarchitecture, and bone quality.

RESULTS

Specimens can be imaged by micro-CT at a resolution starting at 1 μm, but in vivo human imaging is restricted to a voxel size of 82 μm (with actual spatial resolution of ~ 130 μm) due to technical limitations and radiation dose considerations. Presently, this mode is limited to peripheral skeletal regions, such as the wrist and tibia. In contrast, multidetector CT can assess the central skeleton but incurs a higher radiation burden on the subject and provides lower resolution (200-500 μm).

CONCLUSIONS

CT currently provides quantitative measures of bone structure and may be used for estimating bone strength mathematically. The techniques may provide clinically relevant information by enhancing our understanding of fracture risk and establishing the efficacy of antifracture for osteoporosis and other bone metabolic disorders.