Fabric-elasticity relationships of tibial trabecular bone are similar in osteogenesis imperfecta and healthy individuals

Simon, Mathieu<sup>1</sup>; Indermaur, Michael<sup>1</sup>; Schenk, Denis<sup>1</sup>; Hosseinitabatabaei, Seyedmahdi<sup>2,3,4</sup>; Willie, Bettina M.<sup>2,3,4</sup>; Zysset, Philippe<sup>1</sup>

Affiliations

¹ARTORG Centre for Biomedical Engineering Research, University of Bern, Bern, Switzerland

²Research Centre, Shriners Hospital for Children, Montreal, Canada

³Department of Pediatric Surgery, McGill University, Montreal, Canada

⁴Department of Biomedical Engineering, McGill University, Montreal, Canada

Abstract and keywords

Osteogenesis Imperfecta (OI) is an inherited form of bone fragility characterised by impaired synthesis of type I collagen, altered trabecular bone architecture and reduced bone mass. High resolution peripheral computed tomography (HR-pQCT) is a powerful method to investigate bone morphology at peripheral sites including the weight-bearing distal tibia. The resulting 3D reconstructions can be used as a basis of micro-finite element (FE) or homogenized finite element (hFE) models for bone strength estimation. The hFE scheme uses homogenized local bone volume fraction (BV/TV) and anisotropy information (fabric) to compute healthy bone strength within a reasonable computation time using fabric-elasticity relationships. However, it is unclear if these relationships quantified previously for healthy controls are valid for trabecular bone from OI patients. Thus, the aim of this study is to investigate fabric-elasticity relationships in OI trabecular bone compared to healthy controls.

In the present study, the morphology of distal tibiae from 50 adults with OI were compared to 120 healthy controls using second generation HR-pQCT. Six cubic regions of interest (ROIs) were selected per individual in a common anatomical region. A first matching between OI and healthy control group was performed by selecting similar individuals to obtain identical mean and median age and sex distribution. It allowed us to perform a first morphometric analysis and compare the outcome with literature. Then, stiffness tensors of the ROIs were computed using μFE and multiple linear regressions were performed with the Zysset-Curnier orthotropic fabric-elasticity model.

An initial fit was performed on both the OI group and the healthy control group using all extracted ROIs. Then, data was filtered according to a fixed threshold for a defined coefficient of variation (CV) assessing ROI heterogeneity and additional linear regressions were performed on these filtered data sets. These full and filtered data were in turn compared with previous results from μCT reconstructions obtained in other anatomical locations. Finally, the ROIs of both groups were matched according to their BV/TV and degree of anisotropy (DA). Linear regressions were performed using these matched data to detect statistical differences between the two groups.

Compared to healthy controls, we found the OI samples to have significantly lower BV/TV and trabecular number (Tb.N.), significantly higher CV, trabecular separation (Tb.Sp.) and trabecular separation standard deviation (Tb.Sp.SD), but no differences in trabecular thickness (Tb.Th.). These results are in agreement with previous studies. The stiffnesses of highly heterogeneous ROIs were randomly lower with respect to the fabric-elasticity relationships, which reflects the limit of validity of the computational homogenisation methodology. This limitation does not challenge the fabric-elasticity relationship, which extrapolation to heterogeneous ROIs is probably reasonable but can simply not be evaluated with the employed homogenisation methodology. Moreover, due to their low BV/TV, the potential (unknown) errors on these heterogeneous ROIs would have negligible influence on whole bone stiffness in comparison to homogeneous ROIs which are orders of magnitude stiffer.

The filtering of highly heterogeneous ROIs removed these low stiffness ROIs and led to similar correlation coefficients for both OI and healthy groups. Finally, the BV/TV and DA matched data revealed no significant differences in fabric-elasticity parameters between OI and healthy individuals. Moreover, the filtering step did not exclude a particular OI type. Compared to previous studies, the stiffness constants from the 61 μm resolution HR-pQCT ROIs were lower than for the 36 μm resolution μCT ROIs.

In conclusion, OI trabecular bone of the distal tibia was shown to be significantly more heterogeneous and have a lower BV/TV than healthy controls. Despite the reduced linear regression parameters found for HR-pQCT images, the fabric-elasticity relationships between OI and healthy individuals are similar when the trabecular bone ROIs are sufficiently homogeneous to perform the computational stiffness analysis. Accordingly, the elastic properties used for FEA of healthy bones are also valid for OI bones.

Keywords: Bone; Elasticity; Fabric; Osteogenesis imperfecta

Links

DOI: 10.1016/j.bone.2021.116282

PubMed: 34896360

WoS: 000773552300007

History

Received: 2021-07-28

Revised: 2021-11-28

Accepted: 2021-11-30

Online: 2021-12-8

Cited Works (16)

Year	Entry
2012	Pialat JB, Burghardt AJ, Sode M, Link TM, Majumdar S. Visual grading of motion induced image degradation in high resolution peripheral computed tomography: impact of image quality on measures of bone density and micro-architecture. Bone. January 2012;50(1):111-118.
2009	Pahr DH, Zysset PK. A comparison of enhanced continuum FE with micro FE models of human vertebral bodies. J Biomech. March 11, 2009;42(4):455-462.
2020	Fleps I, Bahaloo H, Zysset PK, Ferguson SJ, Pálsson H, Helgason B. Empirical relationships between bone density and ultimate strength: a literature review. J Mech Behav Biomed Mater. October 2020;110:103866.
2011	Varga P, Dall’Ara E, Pahr DH, Pretterklieber M, Zysset PK. Validation of an HR-pQCT-based homogenized finite element approach using mechanical testing of ultra-distal radius sections. Biomech Model Mechanobiol. July 2011;10(4):431-444.
2008	Boutroy S, Van Rietbergen B, Sornay-Rendu E, Munoz F, Bouxsein ML, Delmas PD. Finite element analysis based on in vivo HR-pQCT images of the distal radius is associated with wrist fracture in postmenopausal women. J Bone Miner Res. March 2008;23(3):392-399.
1995	Zysset PK, Curnier A. An alternative model for anisotropic elasticity based on fabric tensors. Mech Mater. November 1995;21(4):243-250.
1998	Zysset PK, Goulet RW, Hollister SJ. A global relationship between trabecular bone morphology and homogenized elastic properties. J Biomech Eng. October 1998;120(5):640-646.
2005	Boutroy S, Bouxsein ML, Munoz F, Delmas PD. In vivo assessment of trabecular bone microarchitecture by high-resolution peripheral quantitative computed tomography. J Clin Endocrinol Metab. December 2005;90(12):6508-6515.
2006	Pasco JA, Seeman E, Henry MJ, Merriman EN, Nicholson GC, Kotowicz MA. The population burden of fractures originates in women with osteopenia, not osteoporosis. Osteoporos Int. September 2006;17(9):1404-1409.
2018	Rolvien T, Stürznickel J, Schmidt FN, Butscheidt S, Schmidt T, Busse B, Mundlos S, Schinke T, Kornak U, Amling M, Oheim R. Comparison of bone microarchitecture between adult osteogenesis imperfecta and early-onset osteoporosis. Calcif Tiss Int. November 2018;103(5):512-521.
2017	Bonaretti S, Majumdar S, Lang TF, Khosla S, Burghardt AJ. The comparability of HR-pQCT bone measurements is improved by scanning anatomically standardized regions. Osteoporos Int. July 2017;28(7):2115-2128.
2017	Ghasem-Zadeh A, Burghardt A, Wang X-F, Iuliano S, Bonaretti S, Bui M, Zebaze R, Seeman E. Quantifying sex, race, and age specific differences in bone microstructure requires measurement of anatomically equivalent regions. Bone. August 2017;101:206-213.
1999	Boyde A, Travers R, Glorieux FH, Jones SJ. The mineralization density of iliac crest bone from children with osteogenesis imperfecta. Calcif Tiss Int. March 1999;64(3):185-190.
2020	Stadelmann MA, Schenk DE, Maquer G, Lenherr C, Buck FM, Bosshardt DD, Hoppe S, Theumann N, Alkalay RN, Zysset PK. Conventional finite element models estimate the strength of metastatic human vertebrae despite alterations of the bone's tissue and structure. Bone. December 2020;141:115598.
2020	Schenk D, Mathis A, Lippuner K, Zysset P. In vivo repeatability of homogenized finite element analysis based on multiple HR-pQCT sections for assessment of distal radius and tibia strength. Bone. December 2020;141:115575.
1985	Cowin SC. The relationship between the elasticity tensor and the fabric tensor. Mech Mater. 1985;4(2):137-147.