Computationally-Optimized Bone Mechanical Modeling from High-Resolution Structural Images

Image-based mechanical modeling of the complex micro-structure of human bone has shown promise as a non-invasive method for characterizing bone strength and fracture risk in vivo. In particular, elastic moduli obtained from image-derived micro-finite element (μFE) simulations have been shown to correlate well with results obtained by mechanical testing of cadaveric bone. However, most existing large-scale finite-element simulation programs require significant computing resources, which hamper their use in common laboratory and clinical environments. In this work, we theoretically derive and computationally evaluate the resources needed to perform such simulations (in terms of computer memory and computation time), which are dependent on the number of finite elements in the image-derived bone model. A detailed description of our approach is provided, which is specifically optimized for μFE modeling of the complex three-dimensional architecture of trabecular bone. Our implementation includes domain decomposition for parallel computing, a novel stopping criterion, and a system for speeding up convergence by pre-iterating on coarser grids. The performance of the system is demonstrated on a dual quad-core Xeon 3.16 GHz CPUs equipped with 40 GB of RAM. Models of distal tibia derived from 3D in-vivo MR images in a patient comprising 200,000 elements required less than 30 seconds to converge (and 40 MB RAM). To illustrate the system's potential for large-scale μFE simulations, axial stiffness was estimated from high-resolution micro-CT images of a voxel array of 90 million elements comprising the human proximal femur in seven hours CPU time. In conclusion, the system described should enable image-based finite-element bone simulations in practical computation times on high-end desktop computers with applications to laboratory studies and clinical imaging

Validation of calcaneus trabecular microstructure measurements by HR-pQCT

OBJECTIVE: Assessment of calcaneus microstructure using high-resolution peripheral quantitative computed tomography (HR-pQCT) might be used to improve fracture risk predictions or to assess responses to pharmacological and physical interventions. To develop a standard clinical protocol for the calcaneus, we validated calcaneus trabecular microstructure measured by HR-pQCT against 'gold-standard' micro-CT measurements. METHODS: Ten human cadaveric feet were scanned in situ using HR-pQCT (isotropic 82μm voxel size) at 100, 150 and 200ms integration times, and at 100ms integration time following removal of the calcaneus from the foot (ex vivo). Dissected portions of these bones were scanned using micro-computed tomography (micro-CT) at an isotropic 17.4μm voxel size. HR-pQCT images were rigidly registered to those obtained with micro-CT and divided into multiple 5mm sided cubes to evaluate and compare morphometric parameters between the modalities. Standard HR-pQCT measurements (derived bone volume fraction (BV/TV(d)); trabecular number, Tb.N; derived trabecular thickness, Tb.Th(d); derived trabecular spacing, Tb.Sp(d)) and corresponding micro-CT voxel-based measurements (BV/TV, Tb.N, Tb.Th, Tb.Sp) were compared. RESULTS: A total of 108 regions of interest were analysed across the 10 specimens. At all integration times HR-pQCT BV/TV(d) was strongly correlated with micro-CT BV/TV (r(2)=0.95-0.98, RMSE=1%), but BV/TV(d) was systematically lower than that measured by micro-CT (mean bias=5%). In contrast, HR-pQCT systematically overestimated Tb.N at all integration times; of the in situ scans, 200ms yielded the lowest mean bias and the strongest correlation with micro-CT (r(2)=0.61, RMSE=0.15mm(-1)). Regional analysis revealed greater accuracy for Tb.N in the superior regions of the calcaneus at all integration times in situ (mean bias=0.44-0.85mm(-1); r(2)=0.70-0.88, p<0.001 versus mean bias=0.63-1.46mm(-1); r(2)<0.10, p≥0.21 for inferior regions). Tb.Sp(d) was underestimated by HR-pQCT compared to micro-CT, but showed similar trends with integration time and the region evaluated as Tb.N. HR-pQCT Tb.Th(d) was also underestimated (mean bias=0.081-0.102mm) and moderately correlated (r(2)=0.55-0.59) with micro-CT Tb.Th, independently from the integration time. Stronger correlations, smaller biases and error were found in the scans of the calcaneus ex vivo compared to in situ. CONCLUSION: Calcaneus trabecular BV/TV(d) and trabecular microstructure, particularly in the superior region of the calcaneus, can be assessed by HR-pQCT. The highest integration time examined, 200ms, compared best with micro-CT. Weaker correlations for microstructure at inferior regions, and also with lower integration times, might limit the use of the proposed protocol, which warrants further investigation in vivo

The SPECTRA Collaboration OMERACT Special Interest Group: Current Research and Future Directions

Objective High-resolution peripheral quantitative computed tomography (HR-pQCT) has the potential to improve radiographic progression determination in clinical trials and longitudinal observational studies. The goal of this work was to describe the current state of research presented at Outcome Measures in Rheumatology (OMERACT) 2016 and ensuing future directions outlined during discussion among attendees. Methods At OMERACT 2016, SPECTRA (Study grouP for xtrEme-Computed Tomography in Rheumatoid Arthritis) introduced efforts to (1) validate the HR-pQCT according to OMERACT guidelines, focusing on rheumatoid arthritis (RA), and (2) find alternatives for automated joint space width (JSW) analysis. The Special Interest Group (SIG) was presented to patient research partners, physicians/researchers, and SIG leaders followed by a 40-min discussion on future directions. Results A consensus definition for RA erosion using HR-pQCT was demonstrated through a systematic literature review and a Delphi exercise. Histopathology and perfusion studies were presented that analyzed the true characteristics of cortical breaks in HR-pQCT images, and to provide criterion validity. Results indicate that readers were able to discriminate between erosion and small vascular channels. Moderate reliability (ICC 0.206–0.871) of direct erosion size measures was shown, which improved (> 0.9) only when experienced readers were considered. Quantification of erosion size was presented for scoring, direct measurement, and volumetric approaches, as well as a reliability exercise for direct measurement. Three methods for JSW measurement were compared, all indicating excellent reproducibility with differences at the extremes (i.e., near-zero and joint edge thickness). Conclusion Initial reports on HR-pQCT are promising; however, to consider its use in clinical trials and longitudinal observational studies, it is imperative to assess the responsiveness of erosion measurement quantification

In vivo assessment of trabecular microarchitecture and bone biomechanical properties by high resolution peripheral quantitative tomography : application to osteoporosis

La microarchitecture osseuse est un des déterminants de la qualité osseuse qui peut maintenant être évaluée in vivo au radius et au tibia distaux avec une résolution isotropique de 82μm par un nouveau scanner à haute résolution (XtremeCT, SCANCO Medical AG). Par ailleurs, l’utilisation d’analyse en éléments finis sur les volumes 3D obtenus permet d’évaluer les propriétés biomécaniques de l’os comme la résistance osseuse. Nous avons montré qu’il s’agissait d’une technique prometteuse pour évaluer la densité, la microarchitecture et les propriétés biomécaniques osseuses au niveau des sites périphériques, notamment parce que ces mesures étaient associées chez la femme avec des fractures ostéoporotiques de toutes sortes. Nous avons également montré que les mêmes mesures étaient tout aussi pertinentes chez l’homme, alors qu’il est moins sujet à l’ostéoporose. Les résultats étaient associés aux fractures ostéoporotiques de toutes sortes, notamment les fractures vertébrales. L’analyse en éléments finis permet donc la mesure in vivo de la résistance osseuse, ce qui pourrait fournir des informations sur la fragilité osseuse et le risque de fracture non accessible par les seules mesures de densité ou de microarchitecture osseuse.Bone microarchitecture is one of the determinants of bone quality that can now be evaluated in vivo at the distal radius and tibia with an isotropic resolution of 82μm with a new high-resolution peripheral scanner (XtremeCT, SCANCO Medical AG). Moreover, the use of finite element analysis on the 3D bone volume acquired allows the assessment of bone biomechanical properties such as bone strength. Our studies show that this technique is promising to assess bone density, microarchitecture and strength at peripheral skeletal sites. Indeed those measures were associated with osteoporotic fractures of all kinds in women. We also demonstrated that those same measures were associated with osteoporotic fractures of all kinds, including vertebral fractures, in men, who are less prone to be affected by osteoporosis. Finite element analysis allows in vivo measurement of bone strength, which might provide additional information about bone fragility and fracture risk that are not assessed by measures of density or microarchitecture

Évaluation de la microarchitecture trabéculaire et des propriétés mécaniques osseuses in vivo chez l’humain par scanner périphérique a haute résolution : application clinique à l’ostéoporose

Bone microarchitecture is one of the determinants of bone quality that can now be evaluated in vivo at the distal radius and tibia with an isotropic resolution of 82μm with a new high-resolution peripheral scanner (XtremeCT, SCANCO Medical AG). Moreover, the use of finite element analysis on the 3D bone volume acquired allows the assessment of bone biomechanical properties such as bone strength. Our studies show that this technique is promising to assess bone density, microarchitecture and strength at peripheral skeletal sites. Indeed those measures were associated with osteoporotic fractures of all kinds in women. We also demonstrated that those same measures were associated with osteoporotic fractures of all kinds, including vertebral fractures, in men, who are less prone to be affected by osteoporosis. Finite element analysis allows in vivo measurement of bone strength, which might provide additional information about bone fragility and fracture risk that are not assessed by measures of density or microarchitecture.La microarchitecture osseuse est un des déterminants de la qualité osseuse qui peut maintenant être évaluée in vivo au radius et au tibia distaux avec une résolution isotropique de 82μm par un nouveau scanner à haute résolution (XtremeCT, SCANCO Medical AG). Par ailleurs, l’utilisation d’analyse en éléments finis sur les volumes 3D obtenus permet d’évaluer les propriétés biomécaniques de l’os comme la résistance osseuse. Nous avons montré qu’il s’agissait d’une technique prometteuse pour évaluer la densité, la microarchitecture et les propriétés biomécaniques osseuses au niveau des sites périphériques, notamment parce que ces mesures étaient associées chez la femme avec des fractures ostéoporotiques de toutes sortes. Nous avons également montré que les mêmes mesures étaient tout aussi pertinentes chez l’homme, alors qu’il est moins sujet à l’ostéoporose. Les résultats étaient associés aux fractures ostéoporotiques de toutes sortes, notamment les fractures vertébrales. L’analyse en éléments finis permet donc la mesure in vivo de la résistance osseuse, ce qui pourrait fournir des informations sur la fragilité osseuse et le risque de fracture non accessible par les seules mesures de densité ou de microarchitecture osseuse

Trabecular Bone Score Helps Classifying Women At Riskof Fracture Prospectively In The Ofely Study

Objectives: Trabecular Bone Score (TBS, Med-Imaps, France) is an index of bone microarchitecture calculated from antero-posterior spine DXA scan and reported to be associated with fracture in prior case-control studies and in a large prospective study with the Prodigy DXA device. Our aim was to assess the ability of TBS to predict incident fracture and improve the classification of fracture prospectively in the OFELY study.Materials/Methods: TBS was assessed in 564 postmenopausal women (66±8 years old) from the OFELY cohort, who had a spine DXA scan (QDR 4500A, Hologic, USA) between year 2000 and 2001. During a mean follow up of 7.8±1.3 years, 94 women sustained a fragility fracture.Results: At the time of baseline DXA scan, women with incident fracture were significantly older (70±9 vs. 65± 8 years), had a lower spine BMD (T-score: −1.9±1.2 vs. −1.3±1.3, p&lt;0.001) and spine TBS (−3.1%, p&lt;0.001) than women without incident fracture. After adjustment for age, BMI and the presence of prevalent fracture, the magnitude of fracture prediction was similar for spine BMD (OR=1.42 [1.11;1.82] per SD decrease [95% CI]) and TBS (OR=1.34 [1.04;1.74]) but the combination of TBS and spine BMD did not improve fracture prediction. Spine BMD and TBS were both correlated with age (respectively r=−0.17 and −0.49, p&lt;0.001) and correlated together with 39% of TBS explained by spine BMD (r=0.63, p&lt;0.001). When using the WHO classification, 38% of the fractures occurred in osteoporotic (fracture rate=29%), 47% in osteopenic (fracture rate=16%) and 15% in women with T-score &gt;−1 (fracture rate=9%). By classifying our population in tertiles of TBS, we found that 47% of the fractures occurred in the lowest tertile of TBS (fracture rate=23%) and 39% of the fracture that occurred in osteopenic women were in the lowest tertile of TBS.Conclusions: Spine BMD and TBS predicted fractures equally well. The addition of TBS to spine BMD added only limited information on fracture risk prediction in our cohort when considering the all range of BMD. Nevertheless combining the osteopenic T-score and the lowest TBS helped defining a subset of osteopenic women at higher risk of fracture.Disclosure of Interest: None declared

Evolution of bone strength with age : a large micro-finite element analysis in women from the ofely cohort

Abstract OC4 European Congress on Osteoporosis & Osteoarthritis (IOF-ECCEO12

Age-related changes in bone strength from HR-pQCT derived microarchitectural parameters with an emphasis on the role of cortical porosity

The high resolution peripheral computed tomography (HR-pQCT) technique has seen recent developments with regard to the assessment of cortical porosity. In this study, we investigated the role of cortical porosity on bone strength in a large cohort of women.The distal radius and distal tibia were scanned by HR-pQCT. We assessed bone strength by estimating the failure load by microfinite element analysis (μFEA), with isotropic and homogeneous material properties. We built a multivariate model to predict it, using a few microarchitecture variables including cortical porosity.Among 857 Caucasian women analyzed with μFEA, we found that cortical and trabecular properties, along with the failure load, impaired slightly with advancing age in premenopausal women, the correlations with age being modest, with |rage| ranging from 0.14 to 0.38. After the onset of the menopause, those relationships with age were stronger for most parameters at both sites, with |rage| ranging from 0.10 to 0.64, notably for cortical porosity and failure load, which were markedly deteriorated with increasing age. Our multivariate model using microarchitecture parameters revealed that cortical porosity played a significant role in bone strength prediction, with semipartial r2=0.22 only at the tibia in postmenopausal women.In conclusion, in our large cohort of women, we observed a small decline of bone strength at the tibia before the onset of menopause. We also found an age-related increase of cortical porosity at both scanned sites in premenopausal women. In postmenopausal women, the relatively high increase of cortical porosity accounted for the decline in bone strength only at the tibia