- Research article
- Open access
- Published:
Osteoporosis imaging: effects of bone preservation on MDCT-based trabecular bone microstructure parameters and finite element models
BMC Medical Imaging volume 15, Article number: 22 (2015)
Abstract
Background
Osteoporosis is defined as a skeletal disorder characterized by compromised bone strength due to a reduction of bone mass and deterioration of bone microstructure predisposing an individual to an increased risk of fracture. Trabecular bone microstructure analysis and finite element models (FEM) have shown to improve the prediction of bone strength beyond bone mineral density (BMD) measurements. These computational methods have been developed and validated in specimens preserved in formalin solution or by freezing. However, little is known about the effects of preservation on trabecular bone microstructure and FEM. The purpose of this observational study was to investigate the effects of preservation on trabecular bone microstructure and FEM in human vertebrae.
Methods
Four thoracic vertebrae were harvested from each of three fresh human cadavers (n = 12). Multi-detector computed tomography (MDCT) images were obtained at baseline, 3 and 6 month follow-up. In the intervals between MDCT imaging, two vertebrae from each donor were formalin-fixed and frozen, respectively. BMD, trabecular bone microstructure parameters (histomorphometry and fractal dimension), and FEM-based apparent compressive modulus (ACM) were determined in the MDCT images and validated by mechanical testing to failure of the vertebrae after 6 months.
Results
Changes of BMD, trabecular bone microstructure parameters, and FEM-based ACM in formalin-fixed and frozen vertebrae over 6 months ranged between 1.0–5.6 % and 1.3–6.1 %, respectively, and were not statistically significant (p > 0.05). BMD, trabecular bone microstructure parameters, and FEM-based ACM as assessed at baseline, 3 and 6 month follow-up correlated significantly with mechanically determined failure load (r = 0.89–0.99; p < 0.05). The correlation coefficients r were not significantly different for the two preservation methods (p > 0.05).
Conclusions
Formalin fixation and freezing up to six months showed no significant effects on trabecular bone microstructure and FEM-based ACM in human vertebrae and may both be used in corresponding in-vitro experiments in the context of osteoporosis.
Background
Osteoporosis is defined as a skeletal disorder characterized by compromised bone strength due to a reduction of bone mass and deterioration of bone microstructure predisposing an individual to an increased risk of fracture [1]. Osteoporotic vertebral and hip fractures are associated with an increased mortality [2]. Due to the aging population, the prevalence of osteoporosis and consecutively the incidence of osteoporotic fractures is expected to increase [3]. Therefore, osteoporosis is classified as public health problem.
The World Health Organisation (WHO) based the diagnosis of osteoporosis on the measurement of bone mineral density (BMD) at the spine and hip using dual-energy X-ray absorptiometry (DXA) [4]. Alternatively, BMD can be assessed by using quantitative computed tomography (QCT). QCT allows for the assessment of volumetric BMD (vBMD), in contrast to DXA assessing areal BMD (aBMD) [5]. Importantly, BMD values often underestimate fracture risk, since osteoporotic fractures frequently occur in patients with non-pathological BMD values [6, 7]. Therefore, considerable research effort has been undertaken to improve fracture risk prediction by using high-resolution imaging techniques including high-resolution peripheral quantitative computed tomography (hr-pQCT), multi-detector computed tomography (MDCT), and magnetic resonance imaging (MRI) [8]. Trabecular bone microstructure parameters and finite element models (FEM) have been computed in the acquired images which significantly improved prediction of bone strength beyond BMD [9]. These computational methods have been developed in human specimens in-vitro and validated with mechanical bone strength measurements as gold standard.
After harvesting, the specimens are usually preserved in formalin solution or are frozen until imaging and mechanical testing is performed. Previous studies have reported that human cortical and trabecular bone samples showed no significant differences in their mechanical properties and tissue parameters (including mineral and lipid content and composition) after freezing or formalin fixation up to several weeks [10, 11]. However, weakened viscoelastic and plastic properties of bovine, murine, and human bone by formalin fixation up to six months were demonstrated as compared to freezing [12–14]. Lochmüller et al. reported that DXA-based BMD measurements in human cadavers within 48 h of death and after 10 months of formalin fixation were not significantly different [15].
Little is known about the effects of preservation on trabecular bone microstructure and FEM in human bone specimens. Therefore, the purpose of our study was to investigate the effects of preservation (formalin fixation and freezing) on trabecular bone microstructure and FEM in intact, human vertebrae as determined by MDCT imaging at baseline, 3 and 6 month follow-up and validated by mechanical testing to failure of the vertebrae after 6 months.
Methods
Specimens
Donors with a history of pathological bone changes other than osteoporosis (i.e., bone metastases, hematological, or metabolic bone disorders) were excluded at the outset. Four thoracic vertebrae between the thoracic vertebra 5 and 12 were harvested from each of three fresh human cadavers (n = 12). The donors consisted of one osteoporotic woman aged 74 years and two non-osteoporotic men aged 46 and 62 years, respectively. They had dedicated their body for educational and research purposes to the local Institute of Anatomy in compliance with local institutional and legislative requirements. The study protocol was reviewed and approved by the local Institutional Review Boards (Ethikkommission der Fakultaet fuer Medizin der Technischen Universitaet Muenchen). The surrounding muscle, fat tissue, and intervertebral discs were completely removed from the vertebrae. Each vertebra was embedded in resin (Rencast Isocyanat and Polyol, Huntsman Group, Bad Säckingen, Germany) up to 2 mm above respectively below their vertebral endplates for the purpose of mechanical testing. The resin fixation was performed with parallel alignment of the upper and lower endplate of the vertebrae with the outer surface of the resin chock to guarantee strict axial loading conditions of the vertebrae during the uniaxial mechanical test.
MDCT imaging was performed at baseline, 3 and 6 month follow-up. In the intervals between the MDCT acquisitions, two vertebrae from each donor were stored in a 3.5 % formalin solution, while the other two vertebrae were stored in sealed plastic bags in a freezer at −40 °C. The vertebrae in the freezer were defrosted for 18 h at 20 °C before 3 and 6 month follow-up MDCT imaging, respectively. All vertebrae were degassed in sodium chloride solution at least 3 h before MDCT imaging to prevent air artifacts. The vertebrae were sealed in vacuum plastic boxes filled with sodium chloride solution during MDCT imaging.
Imaging
MDCT images of the vertebrae at baseline, 3 and 6 month follow-up were acquired by using a clinical whole-body 256-row CT scanner (iCT, Philips Medical Care, Best, Netherlands). Scan parameters were a tube voltage of 120 kVp, a tube load of 585 mAs, an image matrix of 1024 × 1024 pixels, and a field of view of 150 mm. Transverse sections were reconstructed with a high-resolution bone kernel (YE). The interpolated voxel size was of 146 × 146 × 300 μm3, while the real spatial resolution, as determined at ρ50 of the modulation-transfer-function, was 250 × 250 × 600 μm3. A dedicated calibration phantom (Mindways Osteoporosis Phantom, San Francisco, CA, USA) was placed in the scanner mat beneath the vertebrae.
Assessment of BMD and trabecular bone microstructure
MDCT images obtained at baseline, 3 and 6 month follow-up were transferred to a remote LINUX workstation and loaded into an in-house developed program based on IDL (Interactive Data Language, Research Systems, Bolder, CO, USA). Firstly, the 15 most central slices displaying the vertebra equidistant to its endplates were identified. Then, 15 circular regions of interest (ROIs) were manually placed in the ventral half of the vertebral body in the selected slices of the MDCT images similar to QCT-based BMD measurements [5]. The circular ROIs had a diameter of 10 mm (Fig. 1). ROIs’ pixel attenuations in [HU; Hounsfield Units] were converted into BMD values in [mg/cm3 calcium hydroxyapatite] by using the calibration phantom. Afterwards, MDCT images were binarized to calculate trabecular bone microstructure parameters. An optimized global threshold was applied to all MDCT images. Similar to previous studies, 200 mg/cm3 calcium hydroxyapatite was identified as optimized global threshold [16, 17]. Four morphometric parameters were calculated in the ROIs in analogy to standard histomorphometry using the mean intercept length method [18]: bone volume divided by total volume (BV/TV), trabecular number (TbN; [mm−1]), trabecular separation (TbSp; [mm]), and trabecular thickness (TbTh; [mm]). Parameters were labeled as apparent (app.) values, since given the limited spatial resolution they cannot depict the true trabecular microstructure. Furthermore, fractal dimension (FD) as texture measurement of the trabecular bone microstructure was determined in the MDCT images using a box counting algorithm as previously described [16]. The reproducibility error expressed as the root mean square error coefficient of variation amounted to 1.2 % for BMD and ranged between 0.5 % and 2.0 % as outlined in a previous study [16].
FEM
Finite element models (FEM) were computed in the baseline, 3 and 6 month follow-up MDCT images to assess apparent compressive modulus (ACM) of each vertebral body in the superior-inferior direction. Three-dimensional models of the vertebrae were created from the MDCT images by identifying the contour of the vertebrae. The in-plane MDCT resolution was selected as mesh refinement. The uniform hexahedral meshes were generated by using ANSYS Workbench (ANSYS, Canonsburg, PA, USA). The material properties of each element were assigned by using a mapping procedure. Firstly, the elements’ values in [HU] were converted into BMD values ρ BMD in [g/cm3 calcium hydroxyapatite] by using the calibration phantom. Secondly, the elements’ information (position and ρ BMD ) were saved in a text file. Thirdly, a subroutine written in APDL (ANSYS Parametric Design Language) was used to read the text file and assign the material properties to each element (Fig. 2). The equation ρ ash = 1.22 ρ BMD + 0.0526 g/cm3 was used for the conversion of ρ BMD into ρ ash [19, 20]. The isotropic elastic modulus E in [N/mm2] was determined for each element by using the established relationships between E and ρ ash as reported previously [21–23]: E = 33900ρ 2.20 ash ; ρ ash ≤ 0.27, E = 5307ρ ash + 469; 0.27 < ρ ash < 0.6, and E = 10200 ρ 2.01 ash ; ρ ash ≥ 0.6. Each element was assigned a Poisson’s ratio of ν = 0.3 [23]. Finally, the apparent compressive modulus (ACM) of the FEMs in [N/mm2] was obtained by applying a displacement force on one vertebral endplate and fixing the opposite one.
Mechanical testing
After 6 month follow-up MDCT imaging, the resin embedded vertebrae were fixed in a mechanical testing system (Wolpert Werkstoffprüfmaschinen AG, Schaffhausen, Switzerland). The mechanical testing was performed similar to previous studies [16, 24, 25]. Firstly, ten pre-conditioning cycles with uniaxial tension-compression up to a load between 10 N and 400 N with a rate of 5 mm/min were applied. Then, a monotonic, uniaxial compression was performed at the same rate. The load–displacement curve was recorded and vertebral failure load (FL) was defined as the first peak of the load–displacement curve with a subsequent drop of >10 %.
Statistical analysis
The statistical analyses were performed with SPSS (SPSS, Chicago, IL, USA). All tests were done using a two-sided 0.05 level of significance. Mean and standard deviation (SD) of FEM-based ACM, BMD, and trabecular bone microstructure parameters were calculated at each time point separately for the formalin-fixed and frozen vertebrae. The Kolmogorov-Smirnov test showed for most parameters a significant difference from a normal distribution (p < 0.05). Therefore, changes of FEM-based ACM, BMD, and trabecular bone microstructure parameters over 6 months were assessed by using the Friedmann test separately for the formalin-fixed and frozen vertebrae. The root-mean-square coefficients of variation (RMSCV) in [%] were calculated to express the changes of each parameter over time [26]. Correlations between FEM-based ACM, BMD, and trabecular bone microstructure parameters with FL were evaluated with the Spearman’s rank correlation coefficient (r). Significant differences between correlation coefficients were assessed by using the Fisher Z transformation.
Results
Failure load values of all specimens are shown in Table 1. Mean ± SD of FEM-based ACM, BMD, and trabecular bone microstructure parameters of the formalin-fixed and frozen vertebrae at baseline, 3 and 6 month follow-up are listed in Table 2. Changes of the computed parameters expressed as RMSCV ranged between 1.0–5.6 % and 1.3–6.1 % in formalin-fixed and frozen vertebrae, respectively (Table 2). Neither formalin fixation nor freezing significantly changed the computed parameters over six months (p > 0.05; Table 2).
FEM-based ACM, BMD, and trabecular bone microstructure parameters at baseline, 3 and 6 month follow-up showed significant correlations with FL for both, formalin-fixed and frozen vertebrae (p < 0.05; Table 3). Correlations coefficients r ranged between 0.89 and 0.99, and were not significantly different for the two preservation methods as compared by using the Fisher Z transformation (p > 0.05).
Discussion
Human bone specimens have been frequently used to assess mechanical features of new orthopedic implants and to validate new computational methods for the improvement of fracture risk prediction in the context of osteoporosis. Fresh bone specimens would represent the best conditions to guarantee the original structural and mechanical properties. However, the availability of fresh specimens is limited and the setup of many studies requires some type of preservation of the bone tissue due to time constraints. Freezing has the advantages of not significant altering the mechanical properties of human bone specimens [11, 14]. The disadvantage of freezing is the risk of infection of investigators working on bone specimens from a variety of pathogens including HIV and the hepatitis virus [27]. Therefore, bone specimens are often embalmed in formalin solution to minimize the risk of infection. Furthermore, bone specimens available from pathology dissections include high numbers of patients with severe diseases. In contrast, specimens from courses of macroscopic dissections are usually embalmed in formalin solution and could be used to constitute more representative study samples [15]. These cadaver bodies are generally stored for one year or more after they are embalmed. However, formalin fixation may alter the mechanical properties, BMD, and trabecular bone microstructure of the specimens which are particularly important in the context of osteoporosis.
Controversial findings have been reported with regard to the changes of the mechanical properties due to formalin fixation. Haaren et al. reported that long-term preservation by freezing or formalin fixation up to one year did not alter the mechanical properties of cortical bone in goats [28]. In contrast, Wilke et al. reported that formalin fixation strongly influences the mechanical properties of calf spines [29]. Consistently, further studies reported weakened viscoelastic and plastic properties of bovine, murine, and human bone by formalin fixation up to six months as compared to freezing [12–14].
Changes of DXA-based BMD measurements at the lumbar spine and proximal femur due to formalin fixation were assessed by Lochmüller et al. in seven intact human cadavers [15]. They measured BMD within 48 h of death and after 10 months of formalin fixation, and observed no significant deviation in BMD values. Edmondston et al. investigated the correlation of DXA-based BMD and mechanically determined failure load in ten fresh and ten formalin-fixed sheep lumbar spines [30]. The slopes of the regression for BMD and failure load of both groups were not significantly different.
In the line of these studies, we investigated the effects of preservation on QCT-based BMD and advanced computational methods of osteoporosis research, i.e., trabecular bone microstructure parameters and FEM. In consistency with Lochmüller et al., we observed no significant changes of QCT-based BMD over six months for both formalin-fixed and frozen vertebrae [15]. Furthermore, we demonstrated for the first time that trabecular bone microstructure parameters and FEM-based ACM are not significantly altered due to formalin fixation or freezing over six months. Furthermore, the computed parameters correlated well with mechanically determined bone strength independent of the preservation method. These findings are consistent with the association of BMD and failure load in fresh and formalin-fixed sheep lumbar spines as reported by Edmondston et al. [30].
The strength of our study was use of intact, human vertebrae, since this is the most important in-vitro scenario in the context of osteoporosis research. Previous studies have been often limited by investigating trabecular or cortical bone samples only which were sometimes not even harvested from human donors [13, 28, 30]. The limitation of our study was the relatively small sample size, i.e., all the vertebrae were harvested from three donors only.
Conclusions
Formalin fixation and freezing up to six months showed not significant effects on QCT-based BMD, trabecular bone microstructure, and FEM in intact, human vertebrae. Therefore, both preservation methods may be used in corresponding in-vitro experiments in the context of osteoporosis.
References
NIH. NIH Consensus Development Panel on Osteoporosis Prevention, Diagnosis, and Therapy, March 7–29, 2000: highlights of the conference. South Med J. 2001;94:569–73.
Ioannidis G, Papaioannou A, Hopman WM, Akhtar-Danesh N, Anastassiades T, Pickard L, et al. Relation between fractures and mortality: results from the Canadian Multicentre Osteoporosis Study. CMAJ. 2009;181:265–71.
Burge R, Dawson-Hughes B, Solomon DH, Wong JB, King A, Tosteson A. Incidence and economic burden of osteoporosis-related fractures in the United States, 2005–2025. J Bone Miner Res. 2007;22:465–75.
WHO Study Group. Assessment of fracture risk and its application to screening for postmenopausal osteoporosis. Report of a WHO Study Group. World Health Organ Tech Rep Ser. 1994;843:1–129.
Adams JE. Quantitative computed tomography. Eur J Radiol. 2009;71:415–24.
Schuit SC, van der Klift M, Weel AE, de Laet CE, Burger H, Seeman E, et al. Fracture incidence and association with bone mineral density in elderly men and women: the Rotterdam Study. Bone. 2004;34:195–202.
Siris ES, Chen YT, Abbott TA, Barrett-Connor E, Miller PD, Wehren LE, et al. Bone mineral density thresholds for pharmacological intervention to prevent fractures. Arch Intern Med. 2004;164:1108–12.
Link TM. Osteoporosis imaging: state of the art and advanced imaging. Radiology. 2012;263:3–17.
Baum T, Karampinos DC, Liebl H, Rummeny EJ, Waldt S, Bauer JS. High-resolution bone imaging for osteoporosis diagnostics and therapy monitoring using clinical MDCT and MRI. Curr Med Chem. 2013;20:4844–52.
Boskey AL, Cohen ML, Bullough PG. Hard tissue biochemistry: a comparison of fresh-frozen and formalin-fixed tissue samples. Calcif Tissue Int. 1982;34:328–31.
Ohman C, Dall’Ara E, Baleani M, Van Sint JS, Viceconti M. The effects of embalming using a 4 % formalin solution on the compressive mechanical properties of human cortical bone. Clin Biomech (Bristol, Avon). 2008;23:1294–8.
Unger S, Blauth M, Schmoelz W. Effects of three different preservation methods on the mechanical properties of human and bovine cortical bone. Bone. 2010;47:1048–53.
Nazarian A, Hermannsson BJ, Muller J, Zurakowski D, Snyder BD. Effects of tissue preservation on murine bone mechanical properties. J Biomech. 2009;42:82–6.
Linde F, Sorensen HC. The effect of different storage methods on the mechanical properties of trabecular bone. J Biomech. 1993;26:1249–52.
Lochmuller EM, Krefting N, Burklein D, Eckstein F. Effect of fixation, soft-tissues, and scan projection on bone mineral measurements with dual energy X-ray absorptiometry (DXA). Calcif Tissue Int. 2001;68:140–5.
Baum T, Grabeldinger M, Rath C, Grande GE, Burgkart R, Patsch JM, et al. Trabecular bone structure analysis of the spine using clinical MDCT: can it predict vertebral bone strength? J Bone Miner Metab. 2014;32:56–64.
Baum T, Carballido-Gamio J, Huber MB, Muller D, Monetti R, Rath C, et al. Automated 3D trabecular bone structure analysis of the proximal femur--prediction of biomechanical strength by CT and DXA. Osteoporos Int. 2010;21:1553–64.
Parfitt AM, Drezner MK, Glorieux FH, Kanis JA, Malluche H, Meunier PJ, et al. Bone histomorphometry: standardization of nomenclature, symbols, and units. Report of the ASBMR Histomorphometry Nomenclature Committee. J Bone Miner Res. 1987;2:595–610.
Keyak JH, Falkinstein Y. Comparison of in situ and in vitro CT scan-based finite element model predictions of proximal femoral fracture load. Med Eng Phys. 2003;25:781–7.
Les CM, Keyak JH, Stover SM, Taylor KT, Kaneps AJ. Estimation of material properties in the equine metacarpus with use of quantitative computed tomography. J Orthop Res. 1994;12:822–33.
Crawford RP, Cann CE, Keaveny TM. Finite element models predict in vitro vertebral body compressive strength better than quantitative computed tomography. Bone. 2003;33:744–50.
Keller TS. Predicting the compressive mechanical behavior of bone. J Biomech. 1994;27:1159–68.
Kopperdahl DL, Morgan EF, Keaveny TM. Quantitative computed tomography estimates of the mechanical properties of human vertebral trabecular bone. J Orthop Res. 2002;20:801–5.
Chevalier Y, Charlebois M, Pahra D, Varga P, Heini P, Schneider E, et al. A patient-specific finite element methodology to predict damage accumulation in vertebral bodies under axial compression, sagittal flexion and combined loads. Comput Methods Biomech Biomed Engin. 2008;11:477–87.
Dall’Ara E, Pahr D, Varga P, Kainberger F, Zysset P. QCT-based finite element models predict human vertebral strength in vitro significantly better than simulated DEXA. Osteoporos Int. 2012;23:563–72.
Gluer CC, Blake G, Lu Y, Blunt BA, Jergas M, Genant HK. Accurate assessment of precision errors: how to measure the reproducibility of bone densitometry techniques. Osteoporos Int. 1995;5:262–70.
Duma SM, Rudd RW, Crandall JR. A protocol system for testing biohazardous materials in an impact biomechanics research facility. Am Ind Hyg Assoc J. 1999;60:629–34.
van Haaren EH, van der Zwaard BC, van der Veen AJ, Heyligers IC, Wuisman PI, Smit TH. Effect of long-term preservation on the mechanical properties of cortical bone in goats. Acta Orthop. 2008;79:708–16.
Wilke HJ, Krischak S, Claes LE. Formalin fixation strongly influences biomechanical properties of the spine. J Biomech. 1996;29:1629–31.
Edmondston SJ, Singer KP, Day RE, Breidahl PD, Price RI. Formalin fixation effects on vertebral bone density and failure mechanics: an in-vitro study of human and sheep vertebrae. Clin Biomech (Bristol, Avon). 1994;9:175–9.
Acknowledgements
This work was supported by grants of the Deutsche Forschungsgemeinschaft (DFG BA 4085/2-1 and BA 4906/1-1), by the Commission for Clinical Research, Technische Universität München (TUM), TUM School of Medicine, Munich, Germany (Project No. 8762152), and by the Technische Universität München within the funding program Open Access Publishing.
Author information
Authors and Affiliations
Corresponding author
Additional information
Competing interests
The authors declare that they have no competing interests.
Authors’ contributions
Study conception and design: TB, RB, JSB. Acquisition of data: TB, EGG, RB, OG, HL, PMJ, MG, TZ, EJR, SW, JSB. Analysis and interpretation of data: TB, EGG, RB, OG, HL, PMJ, MG, TZ, EJR, SW, JSB. Drafting of Article: TB. Review/revision: TB, EGG, RB, OG, HL, PMJ, MG, TZ, EJR, SW, JSB. Final Approval: TB, EGG, RB, OG, HL, PMJ, MG, TZ, EJR, SW, JSB. All authors read and approved the final manuscript.
Rights and permissions
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. 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.
About this article
Cite this article
Baum, T., Grande Garcia, E., Burgkart, R. et al. Osteoporosis imaging: effects of bone preservation on MDCT-based trabecular bone microstructure parameters and finite element models. BMC Med Imaging 15, 22 (2015). https://doi.org/10.1186/s12880-015-0066-z
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12880-015-0066-z