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

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

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<0.001) and spine TBS (−3.1%, p<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<0.001) and correlated together with 39% of TBS explained by spine BMD (r=0.63, p<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 >−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

(Young Investigator Award) Trabecular Bone Score Helps Classifying Women At Risk Of Fracture Prospectively In The Ofely Study

Introduction: To determine the metabolic effect of teriparatide (TPTD) on bone, 99mTc-MDP skeletal plasma clearance was measured in postmenopausal women with osteoporosis treated with TPTD 20 μg/day. Methods: Ten postmenopausal women with osteoporosis had radionuclide bone scans at baseline, 3, and 18 months after starting TPTD 20 μg/day and after 6 months off therapy. Participants were injected with 600 MBq 99mTc- MDP and whole body bone scans acquired at 10 min, 1, 2, 3, and 4 h. Multiple blood samples were taken between 5 min and 4 h and free 99mTc-MDP measured using ultrafiltration. 99mTc-MDP plasma clearance (Kbone) was evaluated using the Patlak plot method. Regional differences in Kbone were studied by measuring the whole skeleton and subregions. Serum procollagen type I Nterminal propeptide (PINP), bone-specific alkaline phosphatase (BSAP), and urinary N-terminal telopeptide (NTX) were measured at each visit.Discussion: The median increase from baseline in whole skeleton Kbone was 22% (P=0.004) at 3 months and 34% (P= 0.002) at 18 months, decreasing to 0.7% after 6 months off therapy. In subregions, Kbone value increases were statistically significant at 3 months and in all subregions except the pelvis at 18 months. After 6 months off therapy, subregional Kbone values also returned toward baseline. Bone markers increases from baseline were statistically significant at 3 and 18 months (BSAP, 15% and 36%; PINP, 137% and 192%; NTX, 109% and 125%). After 6 months off therapy, PINP and NTX values had declined, though remained above baseline (BSAP, −3%; PINP, 43%; NTX, 56%). Increased Kbone values in the whole body and lower extremities were correlated with increases in most bone markers at 3 and 18 months. Increased skeletal uptake of 99mTc-MDP during treatment with TPTD is indicative of increased bone formation and is supported by increases in bone turnover markers.Conclusion: Changes in Kbone and skeletal uptake measured by radionuclide bone scans in patients taking TPTD are the result of metabolic activity of the drug. These data may provide physicians with useful insights when interpreting bone scan results in this population