Fully automated segmentation of callus by micro-CT compared to biomechanics
© The Author(s). 2017
Received: 2 March 2017
Accepted: 1 July 2017
Published: 11 July 2017
A high percentage of closed femur fractures have slight comminution. Using micro-CT (μCT), multiple fragment segmentation is much more difficult than segmentation of unfractured or osteotomied bone. Manual or semi-automated segmentation has been performed to date. However, such segmentation is extremely laborious, time-consuming and error-prone. Our aim was to therefore apply a fully automated segmentation algorithm to determine μCT parameters and examine their association with biomechanics.
The femura of 64 rats taken after randomised inhibitory or neutral medication, in terms of the effect on fracture healing, and controls were closed fractured after a Kirschner wire was inserted. After 21 days, μCT and biomechanical parameters were determined by a fully automated method and correlated (Pearson’s correlation).
The fully automated segmentation algorithm automatically detected bone and simultaneously separated cortical bone from callus without requiring ROI selection for each single bony structure. We found an association of structural callus parameters obtained by μCT to the biomechanical properties. However, results were only explicable by additionally considering the callus location.
A large number of slightly comminuted fractures in combination with therapies that influence the callus qualitatively and/or quantitatively considerably affects the association between μCT and biomechanics. In the future, contrast-enhanced μCT imaging of the callus cartilage might provide more information to improve the non-destructive and non-invasive prediction of callus mechanical properties. As studies evaluating such important drugs increase, fully automated segmentation appears to be clinically important.
Osteoporosis is frequently regarded as a pathology of cancellous bone, which is the reason why biomechanics and trabecular structures within the bone are of high socioeconomic interest. Human bone biopsies were originally collected from different donor sites (e.g. femoral head, vertebral bodies and iliac crest) to develop three-dimensional (3D) methods for the direct quantification of the actual structure type of cancellous bone. Therefore, no assumptions regarding the structure type (e.g. ‘plate-like’ or ‘rod-like’) are needed . These micro-CT (μCT) measures were then correlated to the biomechanical parameters of the biopsies to determine their predictive value for the mechanical properties of bone samples from healthy individuals . As several significant relationships were determined, various fracture healing studies were performed. For example, murine femur fracture callus mechanical properties were predicted by several μCT-derived basic morphometric indices of callus structure (particularly TMD and BV) [3–5].
Parameters of trabecular morphometry, such as trabecular thickness (Tb. Th.), structure model index (SMI) and degree of anisotropy (DA), are rarely reported as a group. Such omissions may occur because they show lower statistical significance combined with a need for more demanding interpretation. However, these measures should be routinely described because of their potential in future preclinical and clinical studies as surrogate measures of callus mechanical properties .
Eighty-five percent of the skeleton consists of cortical bone providing the majority of mechanical strength and stiffness . However, following a fracture, the situation changes. The callus (peri- and endosteal), which mainly consists of trabecular bone, has to maintain relevant components of stability . The size and shape of the callus can vary extensively; therefore, segmentation of this type of tissue is much more difficult than that of unfractured bone .
Therefore, we aimed to apply a fully automated method to identify and separate the callus, the original cortical bone and the marrow without requiring specific ROIs in a closed fracture model. We also aimed to assess the non-volume-dependent μCT parameters after 21 days for association with biomechanical parameters. To this end, we applied medication with the potential of influencing regeneration in diverse ways, namely by causing inhibition (diclofenac and prednisolone) or by having no effect (cefuroxime).
Sixty-four adult male rats (CRL: WI, aged 16 weeks, mean weight ± SD 500 g ± 50 g, taken from an animal experiment) were acclimatised for 2 weeks prior to surgery. The rats were singly fed (temperature 23–25 °C, humidity 55 ± 5%, 12-h light/dark cycle) and allowed access to water and standard laboratory pellets ad libitum. The local animal research committee authorised the animal experiment in accordance with German legislative requirements (reference number of the Regierung von Oberbayern: 55.2-1-54-2531-15-08).
Numbers of animals/group
A material test machine (Wolpert TZZ 707/386, Istron Wolpert GmbH, Darmstadt, Germany) was used to apply three-point bending of the femura. The femura (n = 39) were placed horizontally with the anterior surface upwards. A distance of 15 mm between the bearing and loading bars for each femur was used as recommended by Turner and Burr .
A round-ended intender with a displacement rate of 5 mm/min was directed vertical to the mid-shaft of the femur to apply a bending load until failure (breaking load). The failure criterion was defined as a force reduction of 80%. The termination criterion was defined as a force reduction of >50 N. The test programme Test&Motion (DOLI Elektronik GmbH, München, Germany) was used to determine, from the load-displacement diagram, the highest point (breaking load (N)) and regression (stiffness (N/mm)). Both failure load and stiffness were detected for the right femura.
This approach did not allow for separation of the cortical bone, trabecular bone and air; therefore, thresholding was used for this purpose . Global thresholding was visually determined by two independent examiners (based on slice-wise 2D comparisons between the grey scale and segmented images of all samples and histograms) and, essentially, by the associated histogram (Fig. 3a–f), in which the relative frequency of each density value (mg HA/cm3) or grey value (dimensionless) was displayed by the height of a corresponding line. Local maxima indicated areas representing (larger) objects or background (black). Three areas could be clearly identified in the histogram: a black area (left peak), which represents air or marrow, a red area (middle plateau), which illustrates the callus, and a white area (right maximum), which represents the original cortical bone. All pixels with grey values higher than the threshold were coloured using the same parameters for callus (sigma (0.8), support (1) and threshold (150)) and original cortical bone (sigma (1.5), support (3) and threshold (370)). Pixels with grey values less than 150 were coloured black.
Pearson’s correlation coefficients between groups were used to examine the association of the median results of micro-CT and biomechanics (unpaired data). Additionally, a linear regression model was fitted to the data with biomechanical data as the dependent variable and μCT data as the independent variable. All statistical analyses were performed using GraphPad Prism Version 6.00 (GraphPad Prism Software® San Diego, USA). All tests were conducted on exploratory two-sided 5% significance levels.
All in all, 64 adult male rats were examined. Figure 1 displays 2D axial μCT grey value images (a and e), contoured images according to the methods described by Nyman (b and f) and Morgan (c and g) and corresponding fully automated segmented images (d and h, callus = red, original bone = white, marrow/air = black) [5, 9]. The contouring of the ROI, as applied in our study, is shown in Fig. 3; notably, the contour does not need to be drawn directly adjacent to bone surface. Conventional contouring with a single ROI adjacent to the callus (Fig. 1b, f) led to overestimation of the TV because it includes all tissue within the outer contour in the equation (e.g. bone marrow, air, cortical bone). As a consequence, the volume-dependent parameters (BV/TV and BMD) are underestimated. Two contours (outer callus = green line and outer cortical bone = red line) enclosing the callus according to conventional contouring enabled the evaluation of all parameters within the ROI in an area not corresponding to the fracture (Fig. 1c). No conventional contouring was possible around the fracture gap, especially if slight comminution was present (Fig. 1g, segmented as in c). Figure 1h shows an adequate separation of all islands of the callus and cortical bone via fully automated segmentation, even within the fracture gap.
Fully automated segmentation was performed with three different thresholds for both the cortical bone (Fig. 3a, c, e) and callus (Fig. 3b, d, f). Using incorrect thresholds for the callus or cortical bone leads to misinterpretation of the segmentation. For instance, a low threshold for cortical bone causes false-positive cortical bone segmentation at the expense of the callus. This erroneous segmentation can be observed as misplaced white (cortical) pixels in the callus (Fig. 3a). In Fig. 3c, d, correct thresholding has been applied, and the optimal distinction between the marrow, callus and cortical bone is obtained with the fully automated method.
Median and interquartile range (IQR)
Slight comminution was reported for a high percentage of all usable closed femur fractures created according to Bonnarens et al. . This finding extensively interferes with correct segmentation and requires reliable, easy segmentation procedures. Therefore, we applied a fully automated method to identify and separate all parts of the callus, original cortical bone and marrow without requiring determination of specific ROIs. In addition, we assessed the non-volume-dependent μCT parameters for association with biomechanics.
Contrary to our findings, the ‘semi-automated’ image segmentation method, as described by Morgan et al. in 2009 and again in 2016, precludes the measurement of the endosteal callus. However, the latter contributes to the stability of the entire structure. Morgan et al. also remark that mechanical contributions occur at the level of the ‘entire’ callus [5, 11]. Dickson et al. showed increased biomechanical stability with increased endosteal calcified tissue volume after femur fractures .
Our study was not originally designed to provide a systematic correlation or comparison of methods; we instead aimed to assess the influence of medication on fracture healing (data not shown). However, when interpreting these data, uncertainty was apparent concerning the impact of various parameters. Therefore, we further evaluated the association between μCT and biomechanical data between different animals. Based on the different processing of the two methods, we could not compare the same femura for association of single values between μCT and biomechanics. However, we could compare the medians of the groups, resulting in fewer values with a minor impact. Nevertheless, conclusions can be drawn from the results. The weight range of our rats was very narrow such that results for femura from different rats are comparable .
Unlike the reported positive correlations between BV and biomechanics of unfractured femura  and the very strong positive correlations between BMC and the structural strength using pQCT and DXA (r > 0.95) [6, 17, 18], our results yielded a strongly negative correlation between the load and BV (r = −0.75), BMC (r = −0.73) or BS (r = −0.70). However, the study approach for a regenerating callus is completely different from studies evaluating unfractured bone . In addition, the level of fracture comminution is positively correlated with the amount of callus . Nyman et al. previously showed an inverse correlation between μCT callus BV and strength (the body tries to compensate for the reduced stability) after 28 days, placing these results in a logical context . With regard to TMD, no relevant association could be seen between the measures (r = −0.08). Shefelbine et al. reported a good predictive value of biomechanical measurements in osteotomy models (clear gap) for density measurements. However, in more complex fractures (fragmentation), which are frequently obtained by closed fracture mechanisms according to Bonnarens et al., density measurements lead to poor predictions of fracture callus biomechanics [10, 19].
Again, in contrast to our initial expectation, the SMI had a strongly positive correlation with the breaking load (r = 0.75). During ageing (osteoporosis) and disease, plates are perforated and connecting rods are dissolved as part of a transition from more stable (and dense) plate-like (SMI ≤0) to less stable rod-like (SMI = 3) trabeculae [2, 20]. The SMI of the prednisolone group was considerably lower (−1.30) than the SMI of the other groups. In combination with the lowest load, the latter led to the positive association between the groups. Unlike the biomechanical evidence, the trabecular structure of the prednisolone-affected callus indicates there is a more stable callus. DA, representing the orientation of the trabecular structure responsible for the anisotropic properties of the bone, revealed a weakly negative correlation (−0.45). . Tb. Th. revealed a strongly negative correlation (r = −0.70). The notable, but unexpectedly, lowest load (46.28) compared with the highest Tb. Th. (0.1004) in the prednisolone group contributed to the negative association between the groups. Qualitative evaluation of the colour-coded μCT reconstructions, representing the Tb. Th., showed slightly thicker trabecular-like struts in the prednisolone group compared with those of the other groups (Fig. 4d). Prednisolone calluses seemed to have a lower porosity (note the increased TMD) compared with the other groups. Around the 21st day, the cartilage and original cortical bone tissue are generally remodelled, increasing the void spaces between the built trabeculae . Histologically, a delayed healing of 5 days has previously been described after prednisolone medication, which might explain the thicker trabeculae and the increased thickness compared with that of the other groups . Additionally, the thicker (Tb. Th.) and plate-like (SMI) trabeculae, assumed to have higher stability, do not increase the stability because of the missing outer periosteal callus fracture gap-bridging that is typically observed within this group (Fig. 4c) . Consistent with our findings, increased outer bony bridging has previously been demonstrated to be able to enhance the correlation of μCT parameters with strength . Therefore, bridges across the fracture gap, particularly the outer periosteal callus, can be considered relevant for biomechanical properties, even if they only slightly contribute to the total callus. Our study has shown that the localisation of the callus is important for biomechanical properties; this might explain the unexpected correlations between several biomechanical and structural measures (BV, SMI and Tb. Th.) [22, 24]. Therefore, therapies that enhance bone healing and lead to a bridging callus with increased load might lead to stronger, positive (vice versa for SMI) associations between the modalities and reach statistical significance . This implication is particularly important in a clinical sense because an increasing number of studies are assessing the influence of medication (e.g. rivaroxaban and tadalafil) on fracture healing and bone remodelling, respectively [25, 26].
The breaking load is considered a more precise parameter compared with stiffness. Because of the well-known strong correlation between the methods, the latter can strengthen the results collected based on the breaking load but is not discussed here .
The fully automated segmentation procedure for identifying and separating the callus components, the original cortical bone and the marrow represents a reliable, easy method for assessing non-volume-dependent μCT parameters. Given the increasing number of studies evaluating drugs with the potential of influencing bone regeneration, fully automated segmentation appears to be clinically important. The localisation of the bone, especially the outer bridging callus, might play a significant role in the interpretation of the data.
Soft tissue also contributes to the mechanical properties of the callus, particularly when the callus tissue transforms to cartilage before mineralisation . In the future, non-invasive and non-destructive contrast-enhanced (in vivo) μCT imaging of the cartilage might provide complementary information for improving the prediction of the callus mechanical properties and thus reduce the number of animals needed in studies.
The authors thank Dr. Andres Laib from Scanco Medical (Scanco Medical AG, Brüttisellen, Switzerland) for providing technical support and image analysis. Furthermore, we thank Professor Karl-Heinz Kunzelmann, Department of Restorative Dentistry, LMU-University, Munich, Germany, for providing the μCT 40. Finally, we thank Dr. Eduardo Grande Garcia, Department of Orthopaedic Sports Medicine, Technische Universität München, for the biomechanical support.
This work was supported by the ‘Komission für Klinische Forschung’ at the Faculty of Medicine (Technische Universität München), KKF grant no. C78-08. The funding body had no influence in the design of the study or the collection, analysis, and interpretation of data or in the writing of the manuscript.
Availability of data and materials
The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.
OB and TT are responsible for the study design and conducted the study. OB and PP are responsible for the data collection. OB, CG and KDW are responsible for the data analysis. OB and AH are responsible for the data interpretation. OB drafted the manuscript. OB, CG, PMP and TT revised the manuscript content. All authors read and approved the final manuscript.
All applicable international, national and institutional guidelines for the care and use of animals were followed. The animal experiment was authorised by the local animal research committee and performed in accordance with German legislative requirements at the Institute of Experimental Oncology and Therapy Research, Centre for Preclinical Research at the Technical University of Munich.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Hildebrand T, Laib A, Muller R, Dequeker J, Ruegsegger P. Direct three-dimensional morphometric analysis of human cancellous bone: microstructural data from spine, femur, iliac crest, and calcaneus. J Bone Miner Res. 1999;14:1167–74.View ArticlePubMedGoogle Scholar
- Bevill G, Eswaran SK, Gupta A, Papadopoulos P, Keaveny TM. Influence of bone volume fraction and architecture on computed large-deformation failure mechanisms in human trabecular bone. Bone. 2006;39:1218–25.View ArticlePubMedGoogle Scholar
- Reifenrath J, Angrisani N, Lalk M, Besdo S. Replacement, refinement, and reduction: necessity of standardization and computational models for long bone fracture repair in animals. J Biomed Mater Res A. 2014;102:2884–900.View ArticlePubMedGoogle Scholar
- Butezloff MM, Zamarioli A, Leoni GB, Sousa-Neto MD, Volpon JB. Whole-body vibration improves fracture healing and bone quality in rats with ovariectomy-induced osteoporosis. Acta Cir Bras. 2015;30:727–35.View ArticlePubMedGoogle Scholar
- Morgan EF, Mason ZD, Chien KB, Pfeiffer AJ, Barnes GL, Einhorn TA, et al. Micro-computed tomography assessment of fracture healing: relationships among callus structure, composition, and mechanical function. Bone. 2009;44:335–44.View ArticlePubMedGoogle Scholar
- Jarvinen TL, Sievanen H, Jokihaara J, Einhorn TA. Revival of bone strength: the bottom line. J Bone Miner Res. 2005;20:717–20.View ArticlePubMedGoogle Scholar
- Bagi CM, Hanson N, Andresen C, Pero R, Lariviere R, Turner CH, et al. The use of micro-CT to evaluate cortical bone geometry and strength in nude rats: correlation with mechanical testing, pQCT and DXA. Bone. 2006;38:136–44.View ArticlePubMedGoogle Scholar
- Gabet Y, Muller R, Regev E, Sela J, Shteyer A, Salisbury K, et al. Osteogenic growth peptide modulates fracture callus structural and mechanical properties. Bone. 2004;35:65–73.View ArticlePubMedGoogle Scholar
- Nyman JS, Munoz S, Jadhav S, Mansour A, Yoshii T, Mundy GR, et al. Quantitative measures of femoral fracture repair in rats derived by micro-computed tomography. J Biomech. 2009;42:891–7.View ArticlePubMedGoogle Scholar
- Auregan JC, Coyle RM, Danoff JR, Burky RE, Akelina Y, Rosenwasser MP. The rat model of femur fracture for bone and mineral research: an improved description of expected comminution, quantity of soft callus and incidence of complications. Bone Joint Res. 2013;2:149–54.View ArticlePubMedPubMed CentralGoogle Scholar
- Morgan EF, Pittman J, DeGiacomo A, Cusher D, de Bakker CM, Mroszczyk KA, et al. BMPR1A antagonist differentially affects cartilage and bone formation during fracture healing. J Orthop Res. 2016;34(12):2096–105. doi:10.1002/jor.23233. Epub 2016 Apr 6.
- Bonnarens F, Einhorn TA. Production of a standard closed fracture in laboratory animal bone. J Orthop Res. 1984;2:97–101.View ArticlePubMedGoogle Scholar
- Turner CH, Burr DB. Basic biomechanical measurements of bone: a tutorial. Bone. 1993;14:595–608.View ArticlePubMedGoogle Scholar
- Bouxsein ML, Boyd SK, Christiansen BA, Guldberg RE, Jepsen KJ, Muller R. Guidelines for assessment of bone microstructure in rodents using micro-computed tomography. J Bone Miner Res. 2010;25:1468–86.View ArticlePubMedGoogle Scholar
- Dickson GR, Geddis C, Fazzalari N, Marsh D, Parkinson I. Microcomputed tomography imaging in a rat model of delayed union/non-union fracture. J Orthop Res. 2008;26:729–36.View ArticlePubMedGoogle Scholar
- Voide R, van Lenthe GH, Muller R. Bone morphometry strongly predicts cortical bone stiffness and strength, but not toughness, in inbred mouse models of high and low bone mass. J Bone Miner Res. 2008;23:1194–203.View ArticlePubMedGoogle Scholar
- Hudelmaier M, Kuhn V, Lochmuller EM, Well H, Priemel M, Link TM, et al. Can geometry-based parameters from pQCT and material parameters from quantitative ultrasound (QUS) improve the prediction of radial bone strength over that by bone mass (DXA)? Osteoporos Int. 2004;15:375–81.View ArticlePubMedGoogle Scholar
- Jarvinen TL, Sievanen H, Kannus P, Jarvinen M. Dual-energy X-ray absorptiometry in predicting mechanical characteristics of rat femur. Bone. 1998;22:551–8.View ArticlePubMedGoogle Scholar
- Shefelbine SJ, Simon U, Claes L, Gold A, Gabet Y, Bab I, et al. Prediction of fracture callus mechanical properties using micro-CT images and voxel-based finite element analysis. Bone. 2005;36:480–8.View ArticlePubMedGoogle Scholar
- Wang J, Zhou B, Liu XS, Fields AJ, Sanyal A, Shi X, et al. Trabecular plates and rods determine elastic modulus and yield strength of human trabecular bone. Bone. 2015;72:71–80.View ArticlePubMedGoogle Scholar
- Hildebrand T, Ruegsegger P. Quantification of bone microarchitecture with the structure model index. Comput Methods Biomech Biomed Engin. 1997;1:15–23.View ArticlePubMedGoogle Scholar
- Gerstenfeld LC, Alkhiary YM, Krall EA, Nicholls FH, Stapleton SN, Fitch JL, et al. Three-dimensional reconstruction of fracture callus morphogenesis. J Histochem Cytochem. 2006;54:1215–28.View ArticlePubMedGoogle Scholar
- Newman RJ, Francis MJ, Duthie RB. Nuclear magnetic resonance studies of experimentally induced delayed fracture union. Clin Orthop Relat Res. 1987;(216):253–61.Google Scholar
- Krischak GD, Augat P, Sorg T, Blakytny R, Kinzl L, Claes L, et al. Effects of diclofenac on periosteal callus maturation in osteotomy healing in an animal model. Arch Orthop Trauma Surg. 2007;127:3–9.View ArticlePubMedGoogle Scholar
- Togral G, Arikan M, Korkusuz P, Hesar RH, Eksioglu MF. Positive effect of tadalafil, a phosphodiesterase-5 inhibitor, on fracture healing in rat femur. Eklem Hastalik Cerrahisi. 2015;26:137–44.View ArticlePubMedGoogle Scholar
- Prodinger PM, Burgkart R, Kreutzer K, Liska F, Pilge H, Schmitt A, et al. Does anticoagulant medication alter fracture-healing? A morphological and biomechanical evaluation of the possible effects of rivaroxaban and enoxaparin using a rat closed fracture model. PLoS One. 2016;11:e0159669.View ArticlePubMedPubMed CentralGoogle Scholar
- Leppanen OV, Sievanen H, Jarvinen TL. Biomechanical testing in experimental bone interventions—may the power be with you. J Biomech. 2008;41:1623–31.View ArticlePubMedGoogle Scholar