Prediction of Incident Hip Fracture with the Estimated Femoral Strength by Finite Element Analysis of DXA Scans in the Study of Osteoporotic Fractures

A bone fractures only when loaded beyond its strength. The purpose of this study was to determine the association of femoral strength, as estimated by finite element (FE) analysis of DXA scans, with incident hip fracture in comparison to hip BMD, FRAX(®) and hip structure analysis (HSA) variables. This prospective case-cohort study included a random sample of 1941 women and 668 incident hip fracture cases (295 in the random sample) during a mean±SD follow-up of 12.8±5.7 yrs from the Study of Osteoporotic Fractures (n=7860 community-dwelling women ≥67 yr of age). We analyzed the baseline DXA scans (Holgoic 1000) of the hip using a validated plane-stress, linear-elastic finite element (FE) model of the proximal femur and estimated the femoral strength during a simulated sideways fall. Cox regression accounting for the case-cohort design assessed the association of estimated femoral strength with hip fracture. The age-BMI-adjusted hazard ratio (HR) per SD decrease for estimated strength (2.21, 95% CI 1.95–2.50) was greater than that for TH BMD (1.86, 95% CI 1.67–2.08; p<0.05), FN BMD (2.04, 95% CI 1.79–2.32; p>0.05), FRAX(®) scores (range 1.32–1.68; p<0.0005) and many HSA variables (range 1.13–2.43; p<0.005), and the association was still significant (p<0.05) after further adjustment for hip BMD or FRAX(®) scores. The association of estimated strength with incident hip fracture was strong (Harrell's C index 0.770), significantly better than TH BMD (0.759, p<0.05) and FRAX(®) scores (0.711–0.743, p<0.0001) but not FN BMD (0.762, p>0.05) Similar findings were obtained for intra- and extra-capsular fractures. In conclusion, the estimated femoral strength from FE analysis of DXA scans is an independent predictor and performs at least as well as FN BMD in predicting incident hip fracture in postmenopausal women

ST-V-Net: incorporating shape prior into convolutional neural networks for proximal femur segmentation

We aim to develop a deep-learning-based method for automatic proximal femur segmentation in quantitative computed tomography (QCT) images. We proposed a spatial transformation V-Net (ST-V-Net), which contains a V-Net and a spatial transform network (STN) to extract the proximal femur from QCT images. The STN incorporates a shape prior into the segmentation network as a constraint and guidance for model training, which improves model performance and accelerates model convergence. Meanwhile, a multi-stage training strategy is adopted to fine-tune the weights of the ST-V-Net. We performed experiments using a QCT dataset which included 397 QCT subjects. During the experiments for the entire cohort and then for male and female subjects separately, 90% of the subjects were used in ten-fold stratified cross-validation for training and the rest of the subjects were used to evaluate the performance of models. In the entire cohort, the proposed model achieved a Dice similarity coefficient (DSC) of 0.9888, a sensitivity of 0.9966 and a specificity of 0.9988. Compared with V-Net, the Hausdorff distance was reduced from 9.144 to 5.917 mm, and the average surface distance was reduced from 0.012 to 0.009 mm using the proposed ST-V-Net. Quantitative evaluation demonstrated excellent performance of the proposed ST-V-Net for automatic proximal femur segmentation in QCT images. In addition, the proposed ST-V-Net sheds light on incorporating shape prior to segmentation to further improve the model performance

Validity and sensitivity of a human cranial finite element model: Implications for comparative studies of biting performance

Finite element analysis (FEA) is a modelling technique increasingly used in anatomical studies investigating skeletal form and function. In the case of the cranium this approach has been applied to both living and fossil taxa to (for example) investigate how form relates to function or infer diet or behaviour. However, FE models of complex musculoskeletal structures always rely on simplified representations because it is impossible completely to image and represent every detail of skeletal morphology, variations in material properties and the complexities of loading at all spatial and temporal scales. The effects of necessary simplifications merit investigation. To this end, this study focuses on one aspect, model geometry, which is particularly pertinent to fossil material where taphonomic processes often destroy the finer details of anatomy or in models built from clinical CTs where the resolution is limited and anatomical details are lost. We manipulated the details of a finite element (FE) model of an adult human male cranium and examined the impact on model performance. First, using digital speckle interferometry, we directly measured strains from the infraorbital region and frontal process of the maxilla of the physical cranium under simplified loading conditions, simulating incisor biting. These measured strains were then compared with predicted values from FE models with simplified geometries that included modifications to model resolution, and how cancellous bone and the thin bones of the circum-nasal and maxillary regions were represented. Distributions of regions of relatively high and low principal strains and principal strain vector magnitudes and directions, predicted by the most detailed FE model, are generally similar to those achieved in vitro. Representing cancellous bone as solid cortical bone lowers strain magnitudes substantially but the mode of deformation of the FE model is relatively constant. In contrast, omitting thin plates of bone in the circum-nasal region affects both mode and magnitude of deformation. Our findings provide a useful frame of reference with regard to the effects of simplifications on the performance of FE models of the cranium and call for caution in the interpretation and comparison of FEA results

A New Hip Fracture Risk Index Derived from FEA-Computed Proximal Femur Fracture Loads and Energies-to-Failure

Hip fracture risk assessment is an important but challenging task. Quantitative CT-based patient specific finite element analysis (FEA) computes the force (fracture load) to break the proximal femur in a particular loading condition. It provides different structural information about the proximal femur that can influence a subject overall fracture risk. To obtain a more robust measure of fracture risk, we used principal component analysis (PCA) to develop a global FEA computed fracture risk index that incorporates the FEA-computed yield and ultimate failure loads and energies to failure in four loading conditions (single-limb stance and impact from a fall onto the posterior, posterolateral, and lateral aspects of the greater trochanter) of 110 hip fracture subjects and 235 age and sex matched control subjects from the AGES-Reykjavik study. We found that the first PC (PC1) of the FE parameters was the only significant predictor of hip fracture. Using a logistic regression model, we determined if prediction performance for hip fracture using PC1 differed from that using FE parameters combined by stratified random resampling with respect to hip fracture status. The results showed that the average of the area under the receive operating characteristic curve (AUC) using PC1 was always higher than that using all FE parameters combined in the male subjects. The AUC of PC1 and AUC of the FE parameters combined were not significantly different than that in the female subjects or in all subjectsComment: 27 pages, 4 figure

The interindividual variation in femoral neck width is associated with the acquisition of predictable sets of morphological and tissue‐quality traits and differential bone loss patterns

A better understanding of femoral neck structure and age‐related bone loss will benefit research aimed at reducing fracture risk. We used the natural variation in robustness (bone width relative to length) to analyze how adaptive processes covary traits in association with robustness, and whether the variation in robustness affects age‐related bone loss patterns. Femoral necks from 49 female cadavers (29–93 years of age) were evaluated for morphological and tissue‐level traits using radiography, peripheral quantitative computed tomography, micro–computed tomography, and ash‐content analysis. Femoral neck robustness was normally distributed and varied widely with a coefficient of variation of 14.9%. Age‐adjusted partial regression analysis revealed significant negative correlations ( p   0.2). The results indicated that slender femora were constructed with a different set of traits compared to robust femora, and that the natural variation in robustness was a determinant of age‐related bone loss patterns. Clinical diagnoses and treatments may benefit from a better understanding of these robustness‐specific structural and aging patterns. © 2012 American Society for Bone and Mineral Research.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/92024/1/1614_ftp.pd

Multi-view information fusion using multi-view variational autoencoders to predict proximal femoral strength

The aim of this paper is to design a deep learning-based model to predict proximal femoral strength using multi-view information fusion. Method: We developed new models using multi-view variational autoencoder (MVAE) for feature representation learning and a product of expert (PoE) model for multi-view information fusion. We applied the proposed models to an in-house Louisiana Osteoporosis Study (LOS) cohort with 931 male subjects, including 345 African Americans and 586 Caucasians. With an analytical solution of the product of Gaussian distribution, we adopted variational inference to train the designed MVAE-PoE model to perform common latent feature extraction. We performed genome-wide association studies (GWAS) to select 256 genetic variants with the lowest p-values for each proximal femoral strength and integrated whole genome sequence (WGS) features and DXA-derived imaging features to predict proximal femoral strength. Results: The best prediction model for fall fracture load was acquired by integrating WGS features and DXA-derived imaging features. The designed models achieved the mean absolute percentage error of 18.04%, 6.84% and 7.95% for predicting proximal femoral fracture loads using linear models of fall loading, nonlinear models of fall loading, and nonlinear models of stance loading, respectively. Compared to existing multi-view information fusion methods, the proposed MVAE-PoE achieved the best performance. Conclusion: The proposed models are capable of predicting proximal femoral strength using WGS features and DXA-derived imaging features. Though this tool is not a substitute for FEA using QCT images, it would make improved assessment of hip fracture risk more widely available while avoiding the increased radiation dosage and clinical costs from QCT.Comment: 16 pages, 3 figure

ST-V-Net: Incorporating Shape Prior Into Convolutional Neural Netwoks For Proximal Femur Segmentation

We aim to develop a deep-learning-based method for automatic proximal femur segmentation in quantitative computed tomography (QCT) images. We proposed a spatial transformation V-Net (ST-V-Net), which contains a V-Net and a spatial transform network (STN) to extract the proximal femur from QCT images. The STN incorporates a shape prior into the segmentation network as a constraint and guidance for model training, which improves model performance and accelerates model convergence. Meanwhile, a multi-stage training strategy is adopted to fine-tune the weights of the ST-V-Net. We performed experiments using a QCT dataset which included 397 QCT subjects. During the experiments for the entire cohort and then for male and female subjects separately, 90% of the subjects were used in ten-fold stratified cross-validation for training and the rest of the subjects were used to evaluate the performance of models. In the entire cohort, the proposed model achieved a Dice similarity coefficient (DSC) of 0.9888, a sensitivity of 0.9966 and a specificity of 0.9988. Compared with V-Net, the Hausdorff distance was reduced from 9.144 to 5.917 mm, and the average surface distance was reduced from 0.012 to 0.009 mm using the proposed ST-V-Net. Quantitative evaluation demonstrated excellent performance of the proposed ST-V-Net for automatic proximal femur segmentation in QCT images. In addition, the proposed ST-V-Net sheds light on incorporating shape prior to segmentation to further improve the model performance

A Deep Learning-Based Method for Automatic Segmentation of Proximal Femur from Quantitative Computed Tomography Images

Purpose: Proximal femur image analyses based on quantitative computed tomography (QCT) provide a method to quantify the bone density and evaluate osteoporosis and risk of fracture. We aim to develop a deep-learning-based method for automatic proximal femur segmentation. Methods and Materials: We developed a 3D image segmentation method based on V-Net, an end-to-end fully convolutional neural network (CNN), to extract the proximal femur QCT images automatically. The proposed V-net methodology adopts a compound loss function, which includes a Dice loss and a L2 regularizer. We performed experiments to evaluate the effectiveness of the proposed segmentation method. In the experiments, a QCT dataset which included 397 QCT subjects was used. For the QCT image of each subject, the ground truth for the proximal femur was delineated by a well-trained scientist. During the experiments for the entire cohort then for male and female subjects separately, 90% of the subjects were used in 10-fold cross-validation for training and internal validation, and to select the optimal parameters of the proposed models; the rest of the subjects were used to evaluate the performance of models. Results: Visual comparison demonstrated high agreement between the model prediction and ground truth contours of the proximal femur portion of the QCT images. In the entire cohort, the proposed model achieved a Dice score of 0.9815, a sensitivity of 0.9852 and a specificity of 0.9992. In addition, an R2 score of 0.9956 (p<0.001) was obtained when comparing the volumes measured by our model prediction with the ground truth. Conclusion: This method shows a great promise for clinical application to QCT and QCT-based finite element analysis of the proximal femur for evaluating osteoporosis and hip fracture risk

Biomechanical stimulation of osteoblast gene expression requires phosphorylation of the RUNX2 transcription factor

Bone can adapt its structure in response to mechanical stimuli. At the cellular level, this involves changes in chromatin organization, gene expression, and differentiation, but the underlying mechanisms are poorly understood. Here we report on the involvement of RUNX2, a bone‐related transcription factor, in this process. Fluid flow shear stress loading of preosteoblasts stimulated translocation of extracellular signal‐regulated kinase (ERK)/mitogen‐activated protein kinase (MAPK) to the nucleus where it phosphorylated RUNX2 on the chromatin of target genes, and increased histone acetylation and gene expression. MAPK signaling and two RUNX2 phosphoacceptor sites, S301 and S319, were critical for this response. Similarly, in vivo loading of mouse ulnae dramatically increased ERK and RUNX2 phosphorylation as well as expression of osteoblast‐related genes. These findings establish ERK/MAPK‐mediated phosphorylation of RUNX2 as a critical step in the response of preosteoblasts to dynamic loading and define a novel mechanism to explain how mechanical signals induce gene expression in bone. © 2012 American Society for Bone and Mineral Research.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/91365/1/1574_ftp.pd

Mapping anisotropy improves QCT-based finite element estimation of hip strength in pooled stance and side-fall load configurations

Hip fractures are one of the most severe consequences of osteoporosis. Compared to the clinical standard of DXA-based aBMD at the femoral neck, QCT-based FEA delivers a better surrogate of femoral strength and gains acceptance for the calculation of hip fracture risk when a CT reconstruction is available. Isotropic, homogenised voxel-based, finite element (hvFE) models are widely used to estimate femoral strength in cross-sectional and longitudinal clinical studies. However, fabric anisotropy is a classical feature of the architecture of the proximal femur and the second determinant of the homogenised mechanical properties of trabecular bone. Due to the limited resolution, fabric anisotropy cannot be derived from clinical CT reconstructions. Alternatively, fabric anisotropy can be extracted from HR-pQCT images of cadaveric femora. In this study, fabric anisotropy from HR-pQCT images was mapped onto QCT-based hvFE models of 71 human proximal femora for which both HR-pQCT and QCT images were available. Stiffness and ultimate load computed from anisotropic hvFE models were compared with previous biomechanical tests in both stance and side-fall configurations. The influence of using the femur-specific versus a mean fabric distribution on the hvFE predictions was assessed. Femur-specific and mean fabric enhance the prediction of experimental ultimate force for the pooled, i.e. stance and side-fall, (isotropic: r2=0.81, femur-specific fabric: r2=0.88, mean fabric: r2=0.86,p<0.001) but not for the individual configurations. Fabric anisotropy significantly improves bone strength prediction for the pooled configurations, and mapped fabric provides a comparable prediction to true fabric. The mapping of fabric anisotropy is therefore expected to help generate more accurate QCT-based hvFE models of the proximal femur for personalised or multiple load configurations