A novel mechanical parameter to quantify the microarchitecture effect on apparent modulus of trabecular bone:​ a computational analysis of ineffective bone mass

Jin, Yongqiang<sup>1</sup>; Zhang, Teng<sup>1</sup>; Cheung, Jason Pui Yin<sup>1</sup>; Wong, Tak Man<sup>1</sup>; Feng, Xiaoreng<sup>1</sup>; Sun, Tianhao<sup>2</sup>; Zu, Haiyue<sup>3</sup>; Sze, Kam Yim<sup>4</sup>; Lu, William Weijia<sup>1</sup>

Affiliations

¹Department of Orthopaedics and Traumatology, Li Ka Shing Faculty of Medicine, The University of Hong Kong, Hong Kong, China

²Department of Orthopaedics and Traumatology, The University of Hong Kong - Shenzhen Hospital, Shenzhen, China

³Department of Orthopaedics and Traumatology, Faculty of Medicine, The Chinese University of Hong Kong, Hong Kong, China

⁴Department of Mechanical Engineering, Faculty of Engineering, The University of Hong Kong, Hong Kong, China

Abstract and keywords

Background: One of the characteristics of osteoporotic bone is the deterioration of trabecular microarchitecture. Previous studies have shown microarchitecture alone can vary the apparent modulus of trabecular bone significantly independent of bone volume fraction (BV/TV) from morphological and topological perspectives. However, modulus is a mechanical quantity and there is a lack of mechanical explanatory parameters. This study aims to propose a novel mechanical parameter to quantify the microarchitecture effect on the apparent modulus of trabecular bone.

Materials and methods: Fourteen human female cadaveric vertebrae were scanned with a dual-energy X-ray (DXA) equipment followed by a micro-CT (μCT) system at 18 μm isotropic resolution. Four trabecular bone specimens (3.46 × 3.46 × 3.46 mm) were obtained from each vertebral body and converted to voxel-based micro finite element (μFE) models. The apparent modulus (E) of the μFE model was computed using a linear micro finite element analysis (μFEA). The normalized apparent modulus (E*) was computed as E divided by BV/TV. The relationship between E and BV/TV was analyzed by linear, power-law and exponential regressions. Linear regression was performed between E* and BV/TV. Ineffective bone mass (InBM) was defined as the bone mass with a negligible contribution to the load-resistance and represented by elements with von Mises stress less than a certain stress threshold. InBM was quantified as the low von Mises stress ratio (LS_VMR), which is the ratio of the number of InBM elements to the total number of elements in the μFE model. An incremental search technique with coarse and fine search intervals of 10 and 1 MPa, respectively, was adopted to determine the stress threshold for calculating LS_VMR of the μFE model. Correlation between E* and LS_VMR was analyzed using linear and power-law models for each stress threshold. The threshold producing the highest coefficient of determination (R²) in the correlation between E* and LS_VMR was taken as the optimal stress threshold for calculating LS_VMR. Linear regression was performed between E and LS_VMR. Multiple linear regression of E against both BV/TV and LS_VMR was further analyzed.

Results: E significantly (p < .001) correlates to BV/TV whereas E* has no significant (p = .75) correlation with BV/TV. Incremental search suggests 59 MPa to be the optimal stress threshold for calculating LS_VMR. BV/TV alone can explain 59% of the variation in E using power-law regression model (E = 2254.64BV/TV^1.04, R² = 0.59, p < .001). LS_VMR alone can explain 48% of the variation in E using linear regression model (E = 1696.4–1647.1LS_VMR, R² = 0.48, p < .001). With these two predictors taken into consideration, 95% of the variation in E can be explained in a multiple linear regression model (E = 1364.89 + 2184.37BV/TV – 1605.38LS_VMR, adjusted R² = 0.95, p < .001).

Conclusion: LS_VMR can be adopted as the mechanical parameter to quantify the microarchitecture effect on the apparent modulus of trabecular bone.

Keywords: Osteoporosis; Trabecular bone; Bone biomechanics; Fracture risk; μFEA

Links

DOI: 10.1016/j.bone.2020.115314

PubMed: 32156663

WoS: 000530069000007

History

Received: 2019-12-02

Revised: 2020-03-03

Accepted: 2020-03-06

Online: 2020-03-08

Cited Works (33)

Year	Entry
1993	Goldstein SA, Goulet R, McCubbrey D. Measurement and significance of three-dimensional architecture to the mechanical integrity of trabecular bone. Calcif Tiss Int. February 1993;53(suppl 1):S127-S133.
1994	Keaveny TM, Guo XE, Wachtel EF, McMahon TA, Hayes WC. Trabecular bone exhibits fully linear elastic behavior and yields at low strains. J Biomech. 1994;27(9):1127-1136.
2001	Wehrli FW, Gomberg BR, Saha PK, Song HK, Hwang SN, Snyder PJ. Digital topological analysis of in vivo magnetic resonance microimages of trabecular bone reveals structural implications of osteoporosis. J Bone Miner Res. October 2001;16(8):1520-1531.
2000	Saha PK, Gomberg BR, Wehrli FW. Three-dimensional digital topological characterization of cancellous bone architecture. Int J Imaging Syst Technol. 2000;11(1):81-90.
1985	Kleerekoper M, Villanueva AR, Stanciu J, Rao DS, Parfitt AM. The role of three-dimensional trabecular microstructure in the pathogenesis of vertebral compression fractures. Calcif Tiss Int. 1985;37(6):594-597.
2004	Schuit SCE, van der Klift M, Weel AEAM, de Laet CEDH, Burger H, Seeman E, Hofman A, Uitterlinden AG, van Leeuwen JPTM, Pols HAP. Fracture incidence and association with bone mineral density in elderly men and women: the Rotterdam Study. Bone. January 2004;34(1):195-202.
2019	Torres AM, Trikanad AA, Aubin CA, Lambers FM, Luna M, Rimnac CM, Zavattieri P, Hernandez CJ. Bone-inspired microarchitectures achieve enhanced fatigue life. Proc Natl Acad Sci USA. December 3, 2019;116(49):24457-24462.
1993	Consensus development conference: diagnosis, prophylaxis, and treatment of osteoporosis. Am J Med. June 1993;94(6):646-650.
1987	Mosekilde L, Mosekilde L, Danielsen CC. Biomechanical competence of vertebral trabecular bone in relation to ash density and age in normal individuals. Bone. 1987;8(2):79-85.
2007	Bessho M, Ohnishi I, Matsuyama J, Matsumoto T, Imai K, Nakamura K. Prediction of strength and strain of the proximal femur by a CT-based finite element method. J Biomech. 2007;40(8):1745-1753.
2008	Kanis JA, McCloskey EV, Johansson H, Oden A, Melton LJ III, Khaltaev N. A reference standard for the description of osteoporosis. Bone. March 2008;42(3):467-475.
1997	Hildebrand T, Rüegsegger P. Quantification of bone microarchitecture with the structure model index. Comput Methods Prog Biomed. 1997;1(1):15-23.
2013	Nishiyama KK, Macdonald HM, Hanley DA, Boyd SK. Women with previous fragility fractures can be classified based on bone microarchitecture and finite element analysis measured with HR-pQCT. Osteoporos Int. May 2013;24(5):1733-1740.
1999	Hildebrand T, Laib A, Müller R, Dequeker J, Rüegsegger P. Direct three-dimensional morphometric analysis of human cancellous bone: microstructural data from spine, femur, iliac crest, and calcaneus. J Bone Miner Res. July 1999;14(7):1167-1174.
1988	Rice JC, Cowin SC, Bowman JA. On the dependence of the elasticity and strength of cancellous bone on apparent density. J Biomech. 1988;21(2):155-168.
2014	Kawalilak CE, Johnston JD, Olszynski WP, Kontulainen SA. Characterizing microarchitectural changes at the distal radius and tibia in postmenopausal women using HR-pQCT. Osteoporos Int. August 2014;25(8):2057-2066.
2000	Niebur GL, Feldstein MJ, Yuen JC, Chen TJ, Keaveny TM. High-resolution finite element models with tissue strength asymmetry accurately predict failure of trabecular bone. J Biomech. December 2000;33(12):1575-1583.
2006	Stauber M, Müller R. Volumetric spatial decomposition of trabecular bone into rods and plates: a new method for local bone morphometry. Bone. 2006;38(4):475-484.
2008	Liu XS, Sajda P, Saha PK, Wehrli FW, Bevill G, Keaveny TM, Guo XE. Complete volumetric decomposition of individual trabecular plates and rods and its morphological correlations with anisotropic elastic moduli in human trabecular bone. J Bone Miner Res. February 2008;23(2):223-235.
2009	Liu XS, Bevill G, Keaveny TM, Sajda P, Guo XE. Micromechanical analyses of vertebral trabecular bone based on individual trabeculae segmentation of plates and rods. J Biomech. February 9, 2009;42(3):249-256.
2000	Legrand E, Chappard D, Pascaretti C, Duquenne M, Krebs S, Rohmer V, Basle M, Audran M. Trabecular bone microarchitecture, bone mineral density, and vertebral fractures in male osteoporosis. J Bone Miner Res. January 2000;15(1):13-19.
2006	Liu XS, Sajda P, Saha PK, Wehrli FW, Guo XE. Quantification of the roles of trabecular microarchitecture and trabecular type in determining the elastic modulus of human trabecular bone. J Bone Miner Res. October 2006;21(10):1608-1617.
1998	Kopperdahl DL, Keaveny TM. Yield strain behavior of trabecular bone. J Biomech. July 1998;31(7):601-608.
2006	Stauber M, Müller R. Age-related changes in trabecular bone microstructures: global and local morphometry. Osteoporos Int. April 2006;17(4):616-626.
2006	Stauber M, Rapillard L, van Lenthe GH, Zysset P, Müller R. Importance of individual rods and plates in the assessment of bone quality and their contribution to bone stiffness. J Bone Miner Res. April 2006;21(4):586-595.
1992	Weinans H, Huiskes R, Grootenboer HJ. The behavior of adaptive bone-remodeling simulation models. J Biomech. December 1992;25(12):1425-1441.
1990	Jensen KS, Mosekilde L, Mosekilde L. A model of vertebral trabecular bone architecture and its mechanical properties. Bone. 1990;11(6):417-423.
2003	Crawford RP, Cann CE, Keaveny TM. Finite element models predict in vitro vertebral body compressive strength better than quantitative computed tomography. Bone. 2003;33(4):744-750.
2000	Kanis JA, Johnell O, Oden A, Sernbo I, Redlund-Johnell I, Dawson A, De Laet C, Jonsson B. Long-term risk of osteoporotic fracture in Malmö. Osteoporos Int. September 2000;11(8):669-674.
1994	Greenspan SL, Myers ER, Maitland LA, Resnick NM, Hayes WC. Fall severity and bone mineral density as risk factors for hip fracture in ambulatory elderly. JAMA. January 12, 1994;271(2):128-133.
1993	Recker RR. Architecture and vertebral fracture. Calcif Tiss Int. 1993;53(suppl 1):S139-S142.
2006	Johnell O, Kanis JA. An estimate of the worldwide prevalence and disability associated with osteoporotic fractures. Osteoporos Int. December 2006;17(12):1726-1733.
1994	Mullender MG, Huiskes R, Weinans H. A physiological approach to the simulation of bone remodeling as a self-organizational control process. J Biomech. November 1994;27(11):1389-1394.

Cited By (1)

Year	Entry
2021	Sabet FA, Koric S, Idkaidek A, Jasiuk I. High-performance computing comparison of implicit and explicit nonlinear finite element simulations of trabecular bone. Comput Methods Prog Biomed. March 2021;200:105870.