Proximal femur specimens: automated 3D trabecular bone mineral density analysis at multidetector CT--correlation with biomechanical strength measurement

Ultra-low dose hip CT-based automated measurement of volumetric bone mineral density at proximal femoral subregions

BACKGROUND

Forty to fifty percent of women and 13%-22% of men experience an osteoporosis-related fragility fracture in their lifetimes. After the age of 50 years, the risk of hip fracture doubles in every 10 years. x-Ray based DXA is currently clinically used to diagnose osteoporosis and predict fracture risk. However, it provides only 2-D representation of bone and is associated with other technical limitations. Thus, alternative methods are needed.

PURPOSE

To develop and evaluate an ultra-low dose (ULD) hip CT-based automated method for assessment of volumetric bone mineral density (vBMD) at proximal femoral subregions.

METHODS

An automated method was developed to segment the proximal femur in ULD hip CT images and delineate femoral subregions. The computational pipeline consists of deep learning (DL)-based computation of femur likelihood map followed by shape model-based femur segmentation and finite element analysis-based warping of a reference subregion labeling onto individual femur shapes. Finally, vBMD is computed over each subregion in the target image using a calibration phantom scan. A total of 100 participants (50 females) were recruited from the Genetic Epidemiology of COPD (COPDGene) study, and ULD hip CT imaging, equivalent to 18 days of background radiation received by U.S. residents, was performed on each participant. Additional hip CT imaging using a clinical protocol was performed on 12 participants and repeat ULD hip CT was acquired on another five participants. ULD CT images from 80 participants were used to train the DL network; ULD CT images of the remaining 20 participants as well as clinical and repeat ULD CT images were used to evaluate the accuracy, generalizability, and reproducibility of segmentation of femoral subregions. Finally, clinical CT and repeat ULD CT images were used to evaluate accuracy and reproducibility of ULD CT-based automated measurements of femoral vBMD.

RESULTS

Dice scores of accuracy (n = 20), reproducibility (n = 5), and generalizability (n = 12) of ULD CT-based automated subregion segmentation were 0.990, 0.982, and 0.977, respectively, for the femoral head and 0.941, 0.970, and 0.960, respectively, for the femoral neck. ULD CT-based regional vBMD showed Pearson and concordance correlation coefficients of 0.994 and 0.977, respectively, and a root-mean-square coefficient of variation (RMSCV) (%) of 1.39% with the clinical CT-derived reference measure. After 3-digit approximation, each of Pearson and concordance correlation coefficients as well as intraclass correlation coefficient (ICC) between baseline and repeat scans were 0.996 with RMSCV of 0.72%. Results of ULD CT-based bone analysis on 100 participants (age (mean ± SD) 73.6 ± 6.6 years) show that males have significantly greater (p < 0.01) vBMD at the femoral head and trochanteric regions than females, while females have moderately greater vBMD (p = 0.05) at the medial half of the femoral neck than males.

CONCLUSION

Deep learning, combined with shape model and finite element analysis, offers an accurate, reproducible, and generalizable algorithm for automated segmentation of the proximal femur and anatomic femoral subregions using ULD hip CT images. ULD CT-based regional measures of femoral vBMD are accurate and reproducible and demonstrate regional differences between males and females.

Associations Among Hip Structure, Bone Mineral Density, and Strength Vary With External Bone Size in White Women

Bone mineral density (BMD) is heavily relied upon to reflect structural changes affecting hip strength and fracture risk. Strong correlations between BMD and strength are needed to provide confidence that structural changes are reflected in BMD and, in turn, strength. This study investigated how variation in bone structure gives rise to variation in BMD and strength and tested whether these associations differ with external bone size. Cadaveric proximal femurs (n = 30, White women, 36-89+ years) were imaged using nanocomputed tomography (nano-CT) and loaded in a sideways fall configuration to assess bone strength and brittleness. Bone voxels within the nano-CT images were projected onto a plane to create pseudo dual-energy X-ray absorptiometry (pseudo-DXA) images consistent with a clinical DXA scan. A validation study using 19 samples confirmed pseudo-DXA measures correlated significantly with those measured from a commercially available DXA system, including bone mineral content (BMC) (R ²  = 0.95), area (R ²  = 0.58), and BMD (R ²  = 0.92). BMD-strength associations were conducted using multivariate linear regression analyses with the samples divided into narrow and wide groups by pseudo-DXA area. Nearly 80% of the variation in strength was explained by age, body weight, and pseudo-DXA BMD for the narrow subgroup. Including additional structural or density distribution information in regression models only modestly improved the correlations. In contrast, age, body weight, and pseudo-DXA BMD explained only half of the variation in strength for the wide subgroup. Including bone density distribution or structural details did not improve the correlations, but including post-yield deflection (PYD), a measure of bone material brittleness, did increase the coefficient of determination to more than 70% for the wide subgroup. This outcome suggested material level effects play an important role in the strength of wide femoral necks. Thus, the associations among structure, BMD, and strength differed with external bone size, providing evidence that structure-function relationships may be improved by judiciously sorting study cohorts into subgroups. © 2022 The Authors. JBMR Plus published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research.

Patient-Specific Finite Element Modeling of the Whole Lumbar Spine Using Clinical Routine Multi-Detector Computed Tomography (MDCT) Data-A Pilot Study

(1) Background: To study the feasibility of developing finite element (FE) models of the whole lumbar spine using clinical routine multi-detector computed tomography (MDCT) scans to predict failure load (FL) and range of motion (ROM) parameters. (2) Methods: MDCT scans of 12 subjects (6 healthy controls (HC), mean age ± standard deviation (SD): 62.16 ± 10.24 years, and 6 osteoporotic patients (OP), mean age ± SD: 65.83 ± 11.19 years) were included in the current study. Comprehensive FE models of the lumbar spine (5 vertebrae + 4 intervertebral discs (IVDs) + ligaments) were generated (L1−L5) and simulated. The coefficients of correlation (ρ) were calculated to investigate the relationship between FE-based FL and ROM parameters and bone mineral density (BMD) values of L1−L3 derived from MDCT (BMDQCT-L1-3). Finally, Mann−Whitney U tests were performed to analyze differences in FL and ROM parameters between HC and OP cohorts. (3) Results: Mean FE-based FL value of the HC cohort was significantly higher than that of the OP cohort (1471.50 ± 275.69 N (HC) vs. 763.33 ± 166.70 N (OP), p < 0.01). A strong correlation of 0.8 (p < 0.01) was observed between FE-based FL and BMDQCT-L1-L3 values. However, no significant differences were observed between ROM parameters of HC and OP cohorts (p = 0.69 for flexion; p = 0.69 for extension; p = 0.47 for lateral bending; p = 0.13 for twisting). In addition, no statistically significant correlations were observed between ROM parameters and BMDQCT- L1-3. (4) Conclusions: Clinical routine MDCT data can be used for patient-specific FE modeling of the whole lumbar spine. ROM parameters do not seem to be significantly altered between HC and OP. In contrast, FE-derived FL may help identify patients with increased osteoporotic fracture risk in the future.

Prediction of Incidental Osteoporotic Fractures at Vertebral-Specific Level Using 3D Non-Linear Finite Element Parameters Derived from Routine Abdominal MDCT

To investigate whether finite element (FE) analysis of the spine in routine thoracic/abdominal multi-detector computed tomography (MDCT) can predict incidental osteoporotic fractures at vertebral-specific level; Baseline routine thoracic/abdominal MDCT scans of 16 subjects (8(m), mean age: 66.1 ± 8.2 years and 8(f), mean age: 64.3 ± 9.5 years) who sustained incidental osteoporotic vertebral fractures as confirmed in follow-up MDCTs were included in the current study. Thoracic and lumbar vertebrae (T5-L5) were automatically segmented, and bone mineral density (BMD), finite element (FE)-based failure-load, and failure-displacement were determined. These values of individual vertebrae were normalized globally (g), by dividing the absolute value with the average of L1-3 and locally by dividing the absolute value with the average of T5-12 and L1-5 for thoracic and lumbar vertebrae, respectively. Mean-BMD of L1-3 was determined as reference. Receiver operating characteristics (ROC) and area under the curve (AUC) were calculated for different normalized FE (Kload, Kdisplacement,K(load)g, and K(displacement)g) and BMD (KBMD, and K(BMD)g) ratio parameter combinations for identifying incidental fractures. Kload, K(load)g, KBMD, and K(BMD)g showed significantly higher discriminative power compared to standard mean BMD of L1-3 (BMDStandard) (AUC = 0.67 for Kload; 0.64 for K(load)g; 0.64 for KBMD; 0.61 for K(BMD)g vs. 0.54 for BMDStandard). The combination of Kload, Kdisplacement, and KBMD increased the AUC further up to 0.77 (p < 0.001). The combination of FE with BMD measurements derived from routine thoracic/abdominal MDCT allowed an improved prediction of incidental fractures at vertebral-specific level.

Influence of thermodisinfection on microstructure of human femoral heads: duration of heat exposition and compressive strength

Allogeneic bone derived from living donors being necessary to match demand for bone transplantation and thermodisinfection of femoral heads is an established sterilization method. During the thermodisinfection the peripheral bone is exposed to maximum 86 °C for 94 min providing 82.5 °C within the center of the femoral head for at least 15 min. This study examined the compression force of the central and representative peripheral regions of native and thermodisinfected human femoral heads to observe wether different duration and intensity of heat exposure might alter mechanic behaviour. Slices from the equatorial region of human femoral heads were taken from each 14 native and thermodisinfected human femoral heads. The central area revealed a significantly higher compression force for native (p ≤ 0.001) and for thermodisinfected bone (p = 0.002 and p = 0.005) compared with peripheral regions since no relevant differences were found between the peripheral and intermediate areas themselves. A small reduction of compression force for thermodisinfected bone was shown since this did not appear significant due to the small number of specimens. The heat exposure did not alter the pre-existing anatomical changes of the microarchitecture of the native femoral heads from the center towards the peripheral regions. The heterogeneity of microstructure of the femoral head might be of interest concerning clinical applications of bone grafts since the difference between native and thermodisinfected bone appears moderate as shown previously. The different quantity of heat exposure did not reveal any significant influence on compression force which might enable thermodisinfection of preformed bone pieces for surgical indications.

Effect of Statistically Iterative Image Reconstruction on Vertebral Bone Strength Prediction Using Bone Mineral Density and Finite Element Modeling: A Preliminary Study

Statistical iterative reconstruction (SIR) using multidetector computed tomography (MDCT) is a promising alternative to standard filtered back projection (FBP), because of lower noise generation while maintaining image quality. Hence, we investigated the feasibility of SIR in predicting MDCT-based bone mineral density (BMD) and vertebral bone strength from finite element (FE) analysis. The BMD and FE-predicted bone strength derived from MDCT images reconstructed using standard FBP (FFBP) and SIR with (FSIR) and without regularization (FSIRB0) were validated against experimental failure loads (Fexp). Statistical iterative reconstruction produced the best quality images with regard to noise, signal-to-noise ratio, and contrast-to-noise ratio. Fexp significantly correlated with FFBP, FSIR, and FSIRB0. FFBP had a significant correlation with FSIRB0 and FSIR. The BMD derived from FBP, SIRB0, and SIR were significantly correlated. Effects of regularization should be further investigated with FE and BMD analysis to allow for an optimal iterative reconstruction algorithm to be implemented in an in vivo scenario.

MDCT-based Finite Element Analysis of Vertebral Fracture Risk: What Dose is Needed?

PURPOSE

The aim of this study was to compare vertebral failure loads, predicted from finite element (FE) analysis of patients with and without osteoporotic vertebral fractures (OVF) at virtually reduced dose levels, compared to standard-dose exposure from multidetector computed tomography (MDCT) imaging and evaluate whether ultra-low dose derived FE analysis can still differentiate patient groups.

MATERIALS AND METHODS

An institutional review board (IRB) approval was obtained for this retrospective study. A total of 16 patients were evaluated at standard-dose MDCT; eight with and eight without OVF. Images were reconstructed at virtually reduced dose levels (i. e. half, quarter and tenth of the standard dose). Failure load was determined at L1-3 from FE analysis and compared between standard, half, quarter, and tenth doses and used to differentiate between fracture and control groups.

RESULTS

Failure load derived at standard dose (3254 ± 909 N and 3794 ± 984 N) did not significantly differ from half (3390 ± 890 N and 3860 ± 1063 N) and quarter dose (3375 ± 915 N and 3925 ± 990 N) but was significantly higher for one tenth dose (4513 ± 1762 N and 4766 ± 1628 N) for fracture and control groups, respectively. Failure load differed significantly between the two groups at standard, half and quarter doses, but not at tenth dose. Receiver operating characteristic (ROC) curve analysis also demonstrated that standard, half, and quarter doses can significantly differentiate the fracture from the control group.

CONCLUSION

The use of MDCT enables a dose reduction of at least 75% compared to standard-dose for an adequate prediction of vertebral failure load based on non-invasive FE analysis.

Pullout strength of monocortical and bicortical screws in meta -physeal and diaphyseal regions of the canine humerus

Objective: Monocortical screws are commonly employed in locking plate fixation, but specific recommendations for their placement are lacking and use of short monocortical screws in metaphyseal bone may be contra indicated. Objectives of this study were to evaluate axial pullout strength of two different lengths of monocortical screws placed in various regions of the canine humerus compared to bicortical screws, and to derive cortical thickness and bone density values for those regions using quantitative computed tomography analysis (QCT).

Methods: The QCT analysis was performed on 36 cadaveric canine humeri for six regions of interest (ROI). A bicortical, short monocortical, or 50% transcortical 3.5 mm screw was implanted in each ROI and axial pullout testing was performed.

Results: Bicortical screws were stronger than monocortical screws in all ROI except the lateral epicondylar crest. Short monocortical metaphyseal screws were weaker than those placed in other regions. The 50% transcortical screws were stronger than the short monocortical screws in the condyle. A linear relationship between screw length and pull-out strength was observed.

Clinical significance: Cortical thickness and bone density measurements were obtained from multiple regions of the canine humerus using QCT. Use of short monocortical screws may contribute to failure of locking plate fixation of humeral fractures, especially when placed in the condyle. When bicortical screw placement is not possible, maximizing monocortical screw length may optimize fixation stability for distal humeral fractures.

Risk of vertebral compression fractures in multiple myeloma patients: A finite-element study

The purpose of this study was to develop and validate a finite element (FE) model to predict vertebral bone strength in vitro using multidetector computed tomography (MDCT) images in multiple myeloma (MM) patients, to serve as a complementing tool to assess fracture risk. In addition, it also aims to differentiate MM patients with and without vertebral compression fractures (VCFs) by performing FE analysis on vertebra segments (T1-L5) obtained from in vivo routine MDCT imaging scans. MDCT-based FE models were developed from the in vitro vertebrae samples and were then applied to the in vivo vertebrae segments of MM patients (n = 4) after validation. Predicted fracture load using FE models correlated significantly with experimentally measured failure load (r = 0.85, P < 0.001). Interestingly, an erratic behavior was observed in patients with fractures (n = 2) and a more gradual change in FE-predicted strength values in patients without fractures (n = 2). Severe geometric deformations were also observed in models that have already attained fractures. Since BMD is not a reliable parameter for fracture risk prediction in MM subjects, it is necessary to use advanced tools such as FE analysis to predict individual fracture risk. If peaks are observed between adjacent segments in an MM patient, it can be safe to conclude that the spine is experiencing regions of structural instability. Such an FE visualization may have therapeutic consequences to prevent MM associated vertebral fractures.

Focal osteoporosis defects play a key role in hip fracture

BACKGROUND

Hip fractures are mainly caused by accidental falls and trips, which magnify forces in well-defined areas of the proximal femur. Unfortunately, the same areas are at risk of rapid bone loss with ageing, since they are relatively stress-shielded during walking and sitting. Focal osteoporosis in those areas may contribute to fracture, and targeted 3D measurements might enhance hip fracture prediction. In the FEMCO case-control clinical study, Cortical Bone Mapping (CBM) was applied to clinical computed tomography (CT) scans to define 3D cortical and trabecular bone defects in patients with acute hip fracture compared to controls. Direct measurements of trabecular bone volume were then made in biopsies of target regions removed at operation.

METHODS

The sample consisted of CT scans from 313 female and 40 male volunteers (158 with proximal femoral fracture, 145 age-matched controls and 50 fallers without hip fracture). Detailed Cortical Bone Maps (c.5580 measurement points on the unfractured hip) were created before registering each hip to an average femur shape to facilitate statistical parametric mapping (SPM). Areas where cortical and trabecular bone differed from controls were visualised in 3D for location, magnitude and statistical significance. Measures from the novel regions created by the SPM process were then tested for their ability to classify fracture versus control by comparison with traditional CT measures of areal Bone Mineral Density (aBMD). In women we used the surgical classification of fracture location ('femoral neck' or 'trochanteric') to discover whether focal osteoporosis was specific to fracture type. To explore whether the focal areas were osteoporotic by histological criteria, we used micro CT to measure trabecular bone parameters in targeted biopsies taken from the femoral heads of 14 cases.

RESULTS

Hip fracture patients had distinct patterns of focal osteoporosis that determined fracture type, and CBM measures classified fracture type better than aBMD parameters. CBM measures however improved only minimally on aBMD for predicting any hip fracture and depended on the inclusion of trabecular bone measures alongside cortical regions. Focal osteoporosis was confirmed on biopsy as reduced sub-cortical trabecular bone volume.

CONCLUSION

Using 3D imaging methods and targeted bone biopsy, we discovered focal osteoporosis affecting trabecular and cortical bone of the proximal femur, among men and women with hip fracture.

Effects of dose reduction on bone strength prediction using finite element analysis

This study aimed to evaluate the effect of dose reduction, by means of tube exposure reduction, on bone strength prediction from finite-element (FE) analysis. Fresh thoracic mid-vertebrae specimens (n = 11) were imaged, using multi-detector computed tomography (MDCT), at different intensities of X-ray tube exposures (80, 150, 220 and 500 mAs). Bone mineral density (BMD) was estimated from the mid-slice of each specimen from MDCT images. Differences in image quality and geometry of each specimen were measured. FE analysis was performed on all specimens to predict fracture load. Paired t-tests were used to compare the results obtained, using the highest CT dose (500 mAs) as reference. Dose reduction had no significant impact on FE-predicted fracture loads, with significant correlations obtained with reference to 500 mAs, for 80 mAs (R²  = 0.997, p < 0.001), 150 mAs (R² = 0.998, p < 0.001) and 220 mAs (R² = 0.987, p < 0.001). There were no significant differences in volume quantification between the different doses examined. CT imaging radiation dose could be reduced substantially to 64% with no impact on strength estimates obtained from FE analysis. Reduced CT dose will enable early diagnosis and advanced monitoring of osteoporosis and associated fracture risk.

The use of routine non density calibrated clinical computed tomography data as a potentially useful screening tool for identifying patients with osteoporosis

OBJECTIVES

To evaluate whether lumbar vertebral body density CT attenuation values measured in Hounsfield Units (HUs) on routine Computed Tomography (CT) examinations can be reliably measured with limited variability, and to evaluate for a correlation between HUs and bone mineral density as measured by dual energy X-ray absorptiometry (DXA) scan.

METHODS

Retrospective review of a total of 249 routine MDCT examinations, performed to measure HUs at the first non-rib bearing lumbar vertebral body on axial images, cross-referenced to the lateral scout image.

RESULTS

The overall ICC and RC for intra-reader variability on CT HU were 0.987 (95% CI 0.973 - 0.999) and 15.664 (95% CI 11.66-16.97). The overall ICC and RDC for inter-reader variability on CT HU were 0.952 (95% CI 0.892 - 0.999) and 30.20 (95% CI 23.73 - 34.48). The ICC and RC for interscanner variability were 0.98 (95% CI 0.95 - 0.99) and 16.67 (95% CI 13.13 - 22.85). The correlation between the L1 HUs and L1 BMD, L1 t-score, and overall t-score was 0.437, 0.392, and 0.400, respectively.

CONCLUSIONS

CT attenuation values of the first lumbar vertebra can be measured on routine abdomen CTs with limited variability despite multiple readers and scanners. Correlation between HU and BMD as measured by DXA scan was only weakly positive, and by this method measuring the density of a lumbar vertebral body from a routine MDCT scan does not provide the sensitivity or specificity necessary for a screening test. However above a certain measured value (180 HU), patients have a low chance of osteoporosis and therefore may not need additional screening, potentially limiting radiation exposure and cost.

Assessing vertebral fracture risk on volumetric quantitative computed tomography by geometric characterization of trabecular bone structure

The current clinical standard for measuring Bone Mineral Density (BMD) is dual X-ray absorptiometry, however more recently BMD derived from volumetric quantitative computed tomography has been shown to demonstrate a high association with spinal fracture susceptibility. In this study, we propose a method of fracture risk assessment using structural properties of trabecular bone in spinal vertebrae. Experimental data was acquired via axial multi-detector CT (MDCT) from 12 spinal vertebrae specimens using a whole-body 256-row CT scanner with a dedicated calibration phantom. Common image processing methods were used to annotate the trabecular compartment in the vertebral slices creating a circular region of interest (ROI) that excluded cortical bone for each slice. The pixels inside the ROI were converted to values indicative of BMD. High dimensional geometrical features were derived using the scaling index method (SIM) at different radii and scaling factors (SF). The mean BMD values within the ROI were then extracted and used in conjunction with a support vector machine to predict the failure load of the specimens. Prediction performance was measured using the root-mean-square error (RMSE) metric and determined that SIM combined with mean BMD features (RMSE = 0.82 ± 0.37) outperformed MDCT-measured mean BMD (RMSE = 1.11 ± 0.33) (p < 10^-4). These results demonstrate that biomechanical strength prediction in vertebrae can be significantly improved through the use of SIM-derived texture features from trabecular bone.

Fast and easy preoperative estimation of cancellous bone mineral density in patients with proximal femur fractures

INTRODUCTION

Postoperative complications after hip fractures in osteoporotic bone such as implant cutout can be reduced by the use of specially designed implants or additional cement augmentation. It is not yet clear at which degree of osteoporosis, patients will profit from implant augmentation or specially designed implants for geriatric patients. As the surgeon ideally should obtain information on local bone quality at the site of implant anchorage already preoperatively, the aim of the study was to develop an easily applicable radiographic method to estimate bone quality in those patients.

MATERIALS AND METHODS

75 patients with unilateral hip fracture were included. Preoperatively, a CT scan with a calibration device was conducted. Postoperatively, DXA scans were performed. The proposed method measures local cancellous bone mineral density in the contralateral and uninjured femoral head. As a control, 15 young and healthy non-osteoporotic subjects were included. Inter- and intraobserver reliability was investigated for a subgroup of 20 patients.

RESULTS

Study group patients had a mean BMD measured by CT scans of 194.2 mg/cm(3) (SD 40.4). There was a statistically significant correlation with data from DXA scans (r = 0.706, p < 0.001). The control group was significantly younger and showed a significantly higher BMD when compared to the study group (p < 0.001). Reliability evaluation showed no statistically significant difference in inter- and intraobserver measurements. Interclass correlation proved to be very high.

CONCLUSION

The proposed method is an easily applicable, reliable and useful tool to estimate bone quality preoperatively using the contralateral hip as a reference. Obtained data may facilitate the decision-making towards the use of further therapeutic measures to improve implant anchorage in osteoporotic bone such as bone cement augmentation. Thus, our method allows for a more individualized surgical treatment of hip fracture patients adapted to the estimated cancellous bone quality of the patient.

Development and optimization of a high-throughput micro-computed tomography imaging method incorporating a novel analysis technique to evaluate bone mineral density of arthritic joints in a rodent model of collagen induced arthritis

Rheumatoid arthritis (RA) is a chronic autoimmune disease resulting in joint inflammation, pain, and eventual bone loss. Bone loss and remodeling caused by symmetric polyarthritis, the hallmark of RA, is readily detectable by bone mineral density (BMD) measurement using micro-CT. Abnormalities in these measurements over time reflect the underlying pathophysiology of the bone. To evaluate the efficacy of anti-rheumatic agents in animal models of arthritis, we developed a high throughput knee and ankle joint imaging assay to measure BMD as a translational biomarker. A bone sample holder was custom designed for micro-CT scanning, which significantly increased assay throughput. Batch processing 3-dimensional image reconstruction, followed by automated image cropping, significantly reduced image processing time. In addition, we developed a novel, automated image analysis method to measure BMD and bone volume of knee and ankle joints. These improvements significantly increased the throughput of ex vivo bone sample analysis, reducing data turnaround from 5 days to 24 hours for a study with 200 rat hind limbs. Taken together, our data demonstrate that BMD, as quantified by micro-CT, is a robust efficacy biomarker with a high degree of sensitivity. Our innovative approach toward evaluation of BMD using optimized image acquisition and novel image processing techniques in preclinical models of RA enables high throughput assessment of anti-rheumatic agents offering a powerful tool for drug discovery.

Effects of Different Loading Patterns on the Trabecular Bone Morphology of the Proximal Femur Using Adaptive Bone Remodeling

In this study, the changes in the bone density of human femur model as a result of different loadings were investigated. The model initially consisted of a solid shell representing cortical bone encompassing a cubical network of interconnected rods representing trabecular bone. A computationally efficient program was developed that iteratively changed the structure of trabecular bone by keeping the local stress in the structure within a defined stress range. The stress was controlled by either enhancing existing beam elements or removing beams from the initial trabecular frame structure. Analyses were performed for two cases of homogenous isotropic and transversely isotropic beams. Trabecular bone structure was obtained for three load cases: walking, stair climbing and stumbling without falling. The results indicate that trabecular bone tissue material properties do not have a significant effect on the converged structure of trabecular bone. In addition, as the magnitude of the loads increase, the internal structure becomes denser in critical zones. Loading associated with the stumbling results in the highest density; whereas walking, considered as a routine daily activity, results in the least internal density in different regions. Furthermore, bone volume fraction at the critical regions of the converged structure is in good agreement with previously measured data obtained from combinations of dual X-ray absorptiometry (DXA) and computed tomography (CT). The results indicate that the converged bone architecture consisting of rods and plates are consistent with the natural bone morphology of the femur. The proposed model shows a promising means to understand the effects of different individual loading patterns on the bone density.

Axial QCT: clinical applications and new developments

Quantitative computed tomography (QCT) is currently undergoing a renaissance, with an increasing number of studies being published and the definition of both QCT-specific osteoporosis thresholds and treatment criteria. Compared with dual-energy X-ray absorptiometry, the current standard bone mineral density technique, QCT has a number of pertinent advantages, including volumetric measurements, less susceptibility to degenerative spine changes, and higher sensitivity to changes in bone mass. Disadvantages include the higher radiation doses and less experience with fracture prediction and therapy monitoring. Over the last 10 yr, a number of novel applications have been described allowing assessment of bone mineral density and bone quality in larger patient populations, developments that may substantially improve patient care.

Using Anisotropic 3D Minkowski Functionals for Trabecular Bone Characterization and Biomechanical Strength Prediction in Proximal Femur Specimens

The ability of Anisotropic Minkowski Functionals (AMFs) to capture local anisotropy while evaluating topological properties of the underlying gray-level structures has been previously demonstrated. We evaluate the ability of this approach to characterize local structure properties of trabecular bone micro-architecture in ex vivo proximal femur specimens, as visualized on multi-detector CT, for purposes of biomechanical bone strength prediction. To this end, volumetric AMFs were computed locally for each voxel of volumes of interest (VOI) extracted from the femoral head of 146 specimens. The local anisotropy captured by such AMFs was quantified using a fractional anisotropy measure; the magnitude and direction of anisotropy at every pixel was stored in histograms that served as a feature vectors that characterized the VOIs. A linear multi-regression analysis algorithm was used to predict the failure load (FL) from the feature sets; the predicted FL was compared to the true FL determined through biomechanical testing. The prediction performance was measured by the root mean square error (RMSE) for each feature set. The best prediction performance was obtained from the fractional anisotropy histogram of AMF Euler Characteristic (RMSE = 1.01 ± 0.13), which was significantly better than MDCT-derived mean BMD (RMSE = 1.12 ± 0.16, p<0.05). We conclude that such anisotropic Minkowski Functionals can capture valuable information regarding regional trabecular bone quality and contribute to improved bone strength prediction, which is important for improving the clinical assessment of osteoporotic fracture risk.

Cortical bone assessed with clinical computed tomography at the proximal femur

Hip fractures are the most serious of all fragility fractures in older people of both sexes. Trips, stumbles, and falls result in fractures of the femoral neck or trochanter, and the incidence of these two common fractures is increasing worldwide as populations age. Although clinical risk factors and chance are important in causation, the ability of a femur to resist fracture also depends on the size and spatial distribution of the bone, its intrinsic material properties, and the loads applied. Over the past two decades, clinical quantitative computed tomography (QCT) studies of living volunteers have provided insight into how the femur changes with advancing age to leave older men and women at increased risk of hip fractures. In this review, we focus on patterns of cortical bone loss associated with hip fracture, age-related changes in cortical bone, and the effects of drugs used to treat osteoporosis. There are several methodologies available to measure cortical bone in vivo using QCT. Most techniques quantify bone density (g/cm(3)), mass (g), and thickness (mm) in selected, predefined or “traditional” regions of interest such as the “femoral neck” or “total hip” region. A recent alternative approach termed “computational anatomy,” uses parametric methods to identify systematic differences, before displaying statistically significant regions as color-scaled maps of density, mass, or thickness on or within a representative femur model. This review will highlight discoveries made using both traditional and computational anatomy methods, focusing on cortical bone of the proximal femur.

Predicting the Biomechanical Strength of Proximal Femur Specimens with Minkowski Functionals and Support Vector Regression

Regional trabecular bone quality estimation for purposes of femoral bone strength prediction is important for improving the clinical assessment of osteoporotic fracture risk. In this study, we explore the ability of 3D Minkowski Functionals derived from multi-detector computed tomography (MDCT) images of proximal femur specimens in predicting their corresponding biomechanical strength. MDCT scans were acquired for 50 proximal femur specimens harvested from human cadavers. An automated volume of interest (VOI)-fitting algorithm was used to define a consistent volume in the femoral head of each specimen. In these VOIs, the trabecular bone micro-architecture was characterized by statistical moments of its BMD distribution and by topological features derived from Minkowski Functionals. A linear multi-regression analysis and a support vector regression (SVR) algorithm with a linear kernel were used to predict the failure load (FL) from the feature sets; the predicted FL was compared to the true FL determined through biomechanical testing. The prediction performance was measured by the root mean square error (RMSE) for each feature set. The best prediction result was obtained from the Minkowski Functional surface used in combination with SVR, which had the lowest prediction error (RMSE = 0.939 ± 0.345) and which was significantly lower than mean BMD (RMSE = 1.075 ± 0.279, p<0.005). Our results indicate that the biomechanical strength prediction can be significantly improved in proximal femur specimens with Minkowski Functionals extracted from on MDCT images used in conjunction with support vector regression.

Improving bone strength prediction in human proximal femur specimens through geometrical characterization of trabecular bone microarchitecture and support vector regression

We investigate the use of different trabecular bone descriptors and advanced machine learning tech niques to complement standard bone mineral density (BMD) measures derived from dual-energy x-ray absorptiometry (DXA) for improving clinical assessment of osteoporotic fracture risk. For this purpose, volumes of interest were extracted from the head, neck, and trochanter of 146 ex vivo proximal femur specimens on multidetector computer tomography. The trabecular bone captured was characterized with (1) statistical moments of the BMD distribution, (2) geometrical features derived from the scaling index method (SIM), and (3) morphometric parameters, such as bone fraction, trabecular thickness, etc. Feature sets comprising DXA BMD and such supplemental features were used to predict the failure load (FL) of the specimens, previously determined through biomechanical testing, with multiregression and support vector regression. Prediction performance was measured by the root mean square error (RMSE); correlation with measured FL was evaluated using the coefficient of determination R2. The best prediction performance was achieved by a combination of DXA BMD and SIM-derived geometric features derived from the femoral head (RMSE: 0.869 ± 0.121, R2: 0.68 ± 0.079), which was significantly better than DXA BMD alone (RMSE: 0.948 ± 0.119, R2: 0.61 ± 0.101) (p < 10-4). For multivariate feature sets, SVR outperformed multiregression (p < 0.05). These results suggest that supplementing standard DXA BMD measurements with sophisticated femoral trabecular bone characterization and supervised learning techniques can significantly improve biomechanical strength prediction in proximal femur specimens.

Trabecular bone structure analysis of the spine using clinical MDCT: can it predict vertebral bone strength?

Recent technical improvements have made it possible to determine trabecular bone structure parameters of the spine using clinical multi-detector computed tomography (MDCT). Therefore, the purpose of this study was to analyze trabecular bone structure parameters obtained from clinical MDCT in relation to high resolution peripheral quantitative computed tomography (HR-pQCT) as a standard of reference and to investigate whether clinical MDCT can predict vertebral bone strength. Fourteen functional spinal segment units between T7 and L3 were harvested from 14 formalin-fixed human cadavers (11 women and 3 men; age 84 ± 10 years). All functional spinal segment units were examined using HR-pQCT (isotropic voxel size of 41 μm(3)) and a clinical whole-body MDCT (interpolated voxel size of 146 × 146 × 300 μm(3)). Trabecular bone structure analyses (histomorphometric and texture measures) were performed in the HR-pQCT as well as MDCT images. Vertebral failure load (FL) of the functional spinal segment units was determined in an uniaxial biomechanical test. The HR-pQCT and MDCT derived trabecular bone structure parameters showed correlations ranging from r = 0.60 to r = 0.90 (p < 0.05). Correlations between trabecular bone structure parameters and FL amounted up to r = 0.86 (p < 0.05) using the HR-pQCT images, and up to r = 0.79 (p < 0.05) using the MDCT images. Correlation coefficients of FL versus trabecular bone structure parameters obtained with HR-pQCT and MDCT were not significantly different (p > 0.05). In this cadaver model, the spatial resolution of clinically available whole-body MDCT scanners was suitable for trabecular bone structure analysis of the spine and to predict vertebral bone strength.

Heterogeneity of bone microstructure in the femoral head in patients with osteoporosis: an ex vivo HR-pQCT study

INTRODUCTION

Trabecular bone in the femoral head has a complicated and heterogeneous structure with few studies having analyzed heterogeneity in this structure quantitatively. We analyze trabecular bone microstructure in the femoral head with osteoporosis (OP) using high resolution peripheral quantitative CT (HR-pQCT) to investigate its regional characteristics.

METHODS

Fifteen femoral heads extracted from female OP patients with femoral neck fracture (85 ± 7, 67-94 years) were scanned by HR-pQCT at 41 μm voxel size. The femoral head was segmented into 15 regions (3 longitudinal regions: superior, center, and inferior, and 5 axial subregions: center, medial, lateral, anterior, posterior). Of these 15 regions, five were excluded due to overlap with the fracture site, leaving a total of 10 regions of cancellous bone microstructures to be quantitatively assessed using the following parameters: bone volume fraction, trabecular thickness, number, separation, connectivity density, structure model index, and degree and orientation of anisotropy. These parameters were compared among each region.

RESULTS

Trabecular bone at the center, superior, and supero-posterior regions of the femoral head had higher bone volume, trabecular number, thickness, narrower bone marrow spaces, higher connectivity and anisotropy, and more plate-like structure. This plate-like structure ran supero-inferiorly and antero-posteriorly at the superior and center regions. Bone volume at the anterior, posterior, and medial regions was almost half of the central and superior regions.

CONCLUSION

Significant heterogeneity of the trabecular bone microstructure in the OP femoral head was showed quantitatively in this study. These data offer new insight into bone microstructural anatomy and may prove to provide useful information on clinical medicine such as hip surgeries.

Impaired bone mineralization accompanied by low vitamin D and secondary hyperparathyroidism in patients with femoral neck fracture

SUMMARY

Although it is well established that a decrease in bone mass increases the risk of osteoporotic fractures, the proportion of fractures attributable to areal bone mineral density (BMD) is rather low. Here, we have identified bone mineralization defects together with low serum 25-hydroxyvitamin D (25-(OH) D) levels as additional factors associated with femoral neck fractures.

INTRODUCTION

Osteoporotic fractures of the femoral neck are associated with increased morbidity and mortality. Although it is well established that a decrease in bone mass increases the risk of osteoporotic fractures, the proportion of fractures attributable to areal BMD is rather low. To identify possible additional factors influencing femur neck fragility, we analyzed patients with femoral neck fracture.

METHODS

We performed a detailed clinical and histomorphometrical evaluation on 103 patients with femoral neck fracture including dual-energy X-ray absorptiometry, laboratory parameters, and histomorphometric and bone mineral density distribution (BMDD) analyses of undecalcified processed biopsies of the femoral head and set them in direct comparison to skeletal healthy control individuals.

RESULTS

Patients with femoral neck fracture displayed significantly lower serum 25-(OH) D levels and increased serum parathyroid hormone (PTH) compared to controls. Histomorphometric analysis revealed not only a decreased bone volume and trabecular thickness in the biopsies of the patients, but also a significant increase of osteoid indices. BMDD analysis showed increased heterogeneity of mineralization in patients with femoral neck fracture. Moreover, patients with femoral neck fracture and serum 25-(OH) D levels below 12 μg/l displayed significantly thinner trabecular bone.

CONCLUSION

Taken together, our data suggest that impaired bone mineralization accompanied by low serum 25-(OH) D levels is of major importance in the etiology of femoral neck fractures. Therefore, balancing serum 25-(OH) D levels and thereby normalizing PTH serum levels may counteract pronounced mineralization defects and might decrease the incidence of femoral neck fractures.

Bone texture analysis is correlated with three-dimensional microarchitecture and mechanical properties of trabecular bone in osteoporotic femurs

Fracture of the proximal femur is a major public health problem in elderly persons. It has recently been suggested that combining texture analysis and bone mineral density measurement improves the failure load prediction in human femurs. In this study, we aimed to compare bone texture analysis with three-dimensional (3D) microarchitecture and mechanical properties of trabecular bone in osteoporotic femurs. Eight femoral heads from osteoporotic patients who fractured their femoral neck provided 31 bone cores. Bone samples were studied using a new high-resolution digital X-ray device (BMA™, D3A Medical Systems) allowing for texture analysis with fractal parameter H (mean), and were examined using micro-computed tomography (microCT) for 3D microarchitecture. Finally, uniaxial compression tests to failure were performed to estimate failure load and apparent modulus of bone samples. The fractal parameter H (mean) was strongly correlated with bone volume fraction (BV/TV) (r = 0.84) and trabecular thickness (Tb.Th) (r = 0.91) (p < 0.01). H (mean) was also markedly correlated with failure load (r = 0.84) and apparent modulus (r = 0.71) of core samples (p < 0.01). Bone volume fraction (BV/TV) and trabecular thickness (Tb.Th) demonstrated significant correlations with failure load (r = 0.85 and 0.72, respectively) and apparent modulus (r = 0.72 and 0.64, respectively) (p < 0.01). Overall, the best predictors of failure load were H (mean), bone volume fraction, and trabecular thickness, with r (2) coefficients of 0.83, 0.76, and 0.80 respectively. This study shows that the fractal parameter H (mean) is correlated with 3D microCT parameters and mechanical properties of femoral head bone samples, which suggests that radiographic texture analysis is a suitable approach for trabecular bone microarchitecture assessment in osteoporotic femurs.

Bone fragility beyond strength and mineral density: Raman spectroscopy predicts femoral fracture toughness in a murine model of rheumatoid arthritis

Clinical prediction of bone fracture risk primarily relies on measures of bone mineral density (BMD). BMD is strongly correlated with bone strength, but strength is independent of fracture toughness, which refers to the bone's resistance to crack initiation and propagation. In that sense, fracture toughness is more relevant to assessing fragility-related fracture risk, independent of trauma. We hypothesized that bone biochemistry, determined by Raman spectroscopy, predicts bone fracture toughness better than BMD. This hypothesis was tested in tumor necrosis factor-transgenic mice (TNF-tg), which develop inflammatory-erosive arthritis and osteoporosis. The left femurs of TNF-tg and wild type (WT) littermates were measured with Raman spectroscopy and micro-computed tomography. Fracture toughness was assessed by cutting a sharp notch into the anterior surface of the femoral mid-diaphysis and propagating the crack under 3 point bending. Femoral fracture toughness of TNF-tg mice was significantly reduced compared to WT controls (p=0.04). A Raman spectrum-based prediction model of fracture toughness was generated by partial least squares regression (PLSR). Raman spectrum PLSR analysis produced strong predictions of fracture toughness, while BMD was not significantly correlated and produced very weak predictions. Raman spectral components associated with mineralization quality and bone collagen were strongly leveraged in predicting fracture toughness, reiterating the limitations of mineralization density alone.

Combination of texture analysis and bone mineral density improves the prediction of fracture load in human femurs

Twenty-one excised femurs were studied using (1) a high-resolution digital X-ray device to estimate three textural parameters, (2) dual-energy X-ray absorptiometry (DXA) to measure bone mineral density (BMD), and (3) mechanical tests to failure. Textural parameters significantly correlated with BMD (p < 0.05) and bone strength (p < 0.05). Combining texture parameters and BMD significantly improved the fracture load prediction from adjusted r(2) = 0.74 to adjusted r(2) =0.82 (p < 0.05).

INTRODUCTION

The purpose of this study is to determine if the combination of bone texture parameters using a new high-resolution X-ray device and BMD measurement by DXA provided a better prediction of femoral failure load than BMD evaluation alone.

METHODS

The proximal ends of 21 excised femurs were studied using (1) a high-resolution digital X-ray device (BMA, D3A Medical Systems) to estimate three textural parameters: fractal parameter Hmean, co-occurrence, and run-length matrices, (2) DXA to measure BMD, and (3) mechanical tests to failure in a side-impact configuration. Regions of interest in the femoral neck, intertrochanteric region, and greater trochanter were selected for DXA and bone texture analysis. Every specimen was scanned twice with repositioning before mechanical testing to assess reproducibility using intraclass correlation coefficient (ICC) with 95% confidence interval. The prediction of femoral failure load was evaluated using multiple regression analysis.

RESULTS

Thirteen femoral neck and 8 intertrochanteric fractures were observed with a mean failure load of 2,612 N (SD, 1,382 N). Fractal parameter Hmean, co-occurrence, and run-length matrices each significantly correlated with site-matched BMD (p < 0.05) and bone strength (p < 0.05). The ICC of the textural parameters varied between 0.65 and 0.90. Combining bone texture parameters and BMD significantly improved the fracture load prediction from adjusted r(2) =0.74 to adjusted r(2) = 0.82 (p < 0.05).

CONCLUSION

In these excised femurs, the combination of bone texture parameters with BMD demonstrated a better performance in the failure load prediction than that of BMD alone.

Generation of an atlas of the proximal femur and its application to trabecular bone analysis

Automatic placement of anatomically corresponding volumes of interest and comparison of parameters against a standard of reference are essential components in studies of trabecular bone. Only recently, in vivo MR images of the proximal femur, an important fracture site, could be acquired with high-spatial resolution. The purpose of this MRI trabecular bone study was two-fold: (1) to generate an atlas of the proximal femur to automatically place anatomically corresponding volumes of interest in a population study and (2) to demonstrate how mean models of geodesic topological analysis parameters can be generated to be used as potential standard of reference. Ten females were used to generate the atlas and geodesic topological analysis models, and 10 females were used to demonstrate the atlas-based trabecular bone analysis. All alignments were based on three-dimensional (3D) multiresolution affine transformations followed by 3D multiresolution free-form deformations. Mean distances less than 1 mm between aligned femora, and sharp edges in the atlas and in fused gray-level images of registered femora indicated that the anatomical variability was well accommodated and explained by the free-form deformations.

The combination of structural parameters and areal bone mineral density improves relation to proximal femur strength: an in vitro study with high-resolution peripheral quantitative computed tomography

The aim of this study was to assess structural indices from high-resolution peripheral quantitative computed tomography (HR-pQCT) images of the human proximal femur along with areal bone mineral density (aBMD) and compare the relationship of these parameters to bone strength in vitro. Thirty-one human proximal femur specimens (8 men and 23 women, median age 74 years, range 50-89) were examined with HR-pQCT at four regions of interest (femoral head, neck, major and minor trochanter) with 82 μm and in a subgroup (n = 17) with 41 μm resolution. Separate analyses of cortical and trabecular geometry, volumetric BMD (vBMD), and microarchitectural parameters were obtained. In addition, aBMD by dual-energy X-ray absorptiometry (DXA) was performed at conventional hip regions and maximal compressive strength (MCS) was determined in a side-impact biomechanical test. Twelve cervical and 19 trochanteric fractures were confirmed. Geometry, vBMD, microarchitecture, and aBMD correlated significantly with MCS, with Spearman's correlation coefficients up to 0.77, 0.89, 0.90, and 0.85 (P < 0.001), respectively. No differences in these correlations were found using 41 μm compared to 82 μm resolution. In multiple regression analysis of MCS, a combined model (age- and sex-adjusted) with aBMD and structural parameters significantly increased R (2) values (up to 0.90) compared to a model holding aBMD alone (R (2) up to 0.78) (P < 0.05). Structural parameters and aBMD are equally related to MCS, and both cortical and trabecular structural parameters obtained from HR-pQCT images hold information on bone strength complementary to that of aBMD.

Prediction of biomechanical properties of trabecular bone in MR images with geometric features and support vector regression

Whole knee joint MR image datasets were used to compare the performance of geometric trabecular bone features and advanced machine learning techniques in predicting biomechanical strength properties measured on the corresponding ex vivo specimens. Changes of trabecular bone structure throughout the proximal tibia are indicative of several musculoskeletal disorders involving changes in the bone quality and the surrounding soft tissue. Recent studies have shown that MR imaging also allows non-invasive 3-D characterization of bone microstructure. Sophisticated features like the scaling index method (SIM) can estimate local structural and geometric properties of the trabecular bone and may improve the ability of MR imaging to determine local bone quality in vivo. A set of 67 bone cubes was extracted from knee specimens and their biomechanical strength estimated by the yield stress (YS) [in MPa] was determined through mechanical testing. The regional apparent bone volume fraction (BVF) and SIM derived features were calculated for each bone cube. A linear multiregression analysis (MultiReg) and a optimized support vector regression (SVR) algorithm were used to predict the YS from the image features. The prediction accuracy was measured by the root mean square error (RMSE) for each image feature on independent test sets. The best prediction result with the lowest prediction error of RMSE = 1.021 MPa was obtained with a combination of BVF and SIM features and by using SVR. The prediction accuracy with only SIM features and SVR (RMSE = 1.023 MPa) was still significantly better than BVF alone and MultiReg (RMSE = 1.073 MPa). The current study demonstrates that the combination of sophisticated bone structure features and supervised learning techniques can improve MR-based determination of trabecular bone quality.

Three-dimensional micro-level computational study of Wolff's law via trabecular bone remodeling in the human proximal femur using design space topology optimization

The law of bone remodeling, commonly referred to as Wolff's Law, asserts that the internal trabecular bone adapts to external loadings, reorienting with the principal stress trajectories to maximize mechanical efficiency creating a naturally optimum structure. The goal of the current study was to utilize an advanced structural optimization algorithm, called design space optimization (DSO), to perform a micro-level three-dimensional finite element bone remodeling simulation on the human proximal femur and analyse the results to determine the validity of Wolff's hypothesis. DSO optimizes the layout of material by iteratively distributing it into the areas of highest loading, while simultaneously changing the design domain to increase computational efficiency. The result is a "fully stressed" structure with minimized compliance and increased stiffness. The large-scale computational simulation utilized a 175 μm mesh resolution and the routine daily loading activities of walking and stair climbing. The resulting anisotropic trabecular architecture was compared to both Wolff's trajectory hypothesis and natural femur samples from literature using a variety of visualization techniques, including radiography and computed tomography (CT). The results qualitatively revealed several anisotropic trabecular regions, that were comparable to the natural human femurs. Quantitatively, the various regional bone volume fractions from the computational results were consistent with quantitative CT analyses. The global strain energy proceeded to become more uniform during optimization; implying increased mechanical efficiency was achieved. The realistic simulated trabecular geometry suggests that the DSO method can accurately predict bone adaptation due to mechanical loading and that the proximal femur is an optimum structure as the Wolff hypothesized.

Densitometric and geometric measurement of the proximal femur in elderly women with and without osteoporotic vertebral fractures by volumetric quantitative multi-slice CT

There is a lack of research on volumetric multi-slice CT (MSCT) application in hip densitometric assessment and geometric measures in elderly women with osteoporotic vertebral fractures. A total of 237 elderly women were divided into three groups based on BMD values of the lumbar spine (AP-SPINE) and/or the femoral neck (NECK) by dual energy X-ray absorptiometry (DXA): osteoporosis with (OP_FX, 53 cases) or without vertebral fracture (OP_NONFX, 94 cases), or normal BMD (CONTROL, 90 cases). Volumetric BMD of trabecular bone (TRAB), integral bone (INTGL) and cortical bone (CORT) with neck axis length (NAL) and minimum cross-section area (mCSA) measures of the left femoral neck were calculated, respectively, by using OsteoCAD software based on MSCT scans of the abdominal-pelvic region of all participants, then the index of femoral neck strength (FNSI) was estimated. The values of TRAB, CORT and INTGL of OP_FX were significantly lower than those of OP_NONFX, with the decrease in 6.8-21.8%, as well as being lower than those in CONTROL, whereas no significant differences in the values of AP-SPINE and NECK were found between OP_FX and OP_NONFX. No significant difference of the value of mCSA was found among these three groups. The NAL value of OP_NONFX was larger than that of CONTROL. FNSI of femoral neck in OP_FX (0.42 ± 0.15 g(2)/cm(4)) was significantly lower than OP_NONFX (0.50 ± 0.14 g(2)/cm(4)) (p < 0.05). vQCT measurement seemed to be more effective than DXA in evaluating hip densitometric changes and discriminating osteoporotic elderly subjects with fractured vertebrae from the non-fractured in a group of Chinese women.

Development of a novel method for surgical implant design optimization through noninvasive assessment of local bone properties

A method was developed to improve the design of locking implants by finding the optimal paths for the anchoring elements, based on a high resolution pQCT assessment of local bone mineral density (BMD) distribution and bone micro-architecture (BMA). The method consists of three steps: (1) partial fixation of the implant to the bone and creation of a reference system, (2) implant removal and pQCT scan of the bone, and (3) determination of BMD and BMA of all implant-anchoring locations along the actual and alternative directions. Using a PHILOS plate, the method uncertainty was tested on an artificial humerus bone model. A cadaveric humerus was used to quantify how the uncertainty of the method affects the assessment of bone parameters. BMD and BMA were determined along four possible alternative screw paths as possible criteria for implant optimization. The method is biased by a 0.87 ± 0.12 mm systematic uncertainty and by a 0.44 ± 0.09 mm random uncertainty in locating the virtual screw position. This study shows that this method can be used to find alternative directions for the anchoring elements, which may possess better bone properties. This modification will thus produce an optimized implant design.

Automated 3D trabecular bone structure analysis of the proximal femur--prediction of biomechanical strength by CT and DXA

SUMMARY

The standard diagnostic technique for assessing osteoporosis is dual X-ray absorptiometry (DXA) measuring bone mass parameters. In this study, a combination of DXA and trabecular structure parameters (acquired by computed tomography [CT]) most accurately predicted the biomechanical strength of the proximal femur and allowed for a better prediction than DXA alone.

INTRODUCTION

An automated 3D segmentation algorithm was applied to determine specific structure parameters of the trabecular bone in CT images of the proximal femur. This was done to evaluate the ability of these parameters for predicting biomechanical femoral bone strength in comparison with bone mineral content (BMC) and bone mineral density (BMD) acquired by DXA as standard diagnostic technique.

METHODS

One hundred eighty-seven proximal femur specimens were harvested from formalin-fixed human cadavers. BMC and BMD were determined by DXA. Structure parameters of the trabecular bone (i.e., morphometry, fuzzy logic, Minkowski functionals, and the scaling index method [SIM]) were computed from CT images. Absolute femoral bone strength was assessed with a biomechanical side-impact test measuring failure load (FL). Adjusted FL parameters for appraisal of relative bone strength were calculated by dividing FL by influencing variables such as body height, weight, or femoral head diameter.

RESULTS

The best single parameter predicting FL and adjusted FL parameters was apparent trabecular separation (morphometry) or DXA-derived BMC or BMD with correlations up to r = 0.802. In combination with DXA, structure parameters (most notably the SIM and morphometry) added in linear regression models significant information in predicting FL and all adjusted FL parameters (up to R(adj) = 0.872) and allowed for a significant better prediction than DXA alone.

CONCLUSION

A combination of bone mass (DXA) and structure parameters of the trabecular bone (linear and nonlinear, global and local) most accurately predicted absolute and relative femoral bone strength.

Prediction of femoral fracture load: cross-sectional study of texture analysis and geometric measurements on plain radiographs versus bone mineral density

PURPOSE

To use standard radiographs to determine which combination of co-occurrence textural parameters, geometric measurements, and cortical thickness measurements from femur radiographs provided the best estimate of femoral failure load and to compare these with total hip dual-energy x-ray absorptiometry bone mineral density (BMD) evaluation.

MATERIALS AND METHODS

Digital radiographs of 40 pairs of excised femurs (24 women, 16 men; mean age, 82 years + or - 12 [standard deviation]) were obtained. Regions of interest in the femoral neck, greater trochanter, intertrochanteric area, and femoral head were then selected. Three textural parameters derived from a co-occurrence matrix were estimated with imaging software. Neck-shaft angle, femoral neck axis length, calcar femorale thickness, and internal and external femoral shaft thickness were assessed. The femurs were randomly allocated to single-stance (femoral neck fracture) or side-impact (intertrochanteric fracture) configurations for failure load measurement.

RESULTS

Textural parameters correlated significantly with site-matched BMD. Stepwise regression analysis was performed, and total hip BMD explained 73% and 78% of the failure load in single-stance and side-impact configurations, respectively. Combining internal femoral shaft thickness with one or two textural parameters explained 72%-79% of failure load variance in the single-stance configuration and 63%-76% of failure load variance in the side-impact configuration.

CONCLUSION

In these excised femurs, combining textural parameters with cortical thickness measurements had a performance comparable to that of BMD alone in the explanation of femoral failure load.

Changing structure of the femoral neck across the adult female lifespan

The anatomic distribution of cortical and cancellous bone in the femoral neck may be critical in determining resistance to fracture. We investigated the effects of aging on femoral neck bone in women. In this cross-sectional study, we used clinical multidetector computed tomography (MDCT) of the hips to investigate aging effects in 100 female volunteers aged 20 to 90 years. We developed a clinically efficient protocol to measure cortical thickness (C.Th) and cortical, trabecular, and integral bone mineral density (CtBMD, TrBMD, and iBMD in mg/cm(3)) in anatomic quadrants of the femoral neck. We used a nested ANOVA to evaluate their associations with height, weight, location in the femoral neck, and age of the subject. Age was the principal determinant of both cortical thickness and BMD. Age had significantly different effects within the anatomic quadrants; compared with young women, elderly subjects had relative preservation of the inferoanterior (IA) quadrant but strikingly reduced C.Th and BMD superiorly. A model including height, weight, and region of interest (and their interactions) explained 83% of the measurement variance (p < .0001). There were marked C.Th and BMD differences between age 25 and age 85 in the already thin superior quadrants. At 25 years the predicted C.Th of the superoposterior quadrant was 1.63 mm, whereas at 85 years it was 0.33 mm [-1.33 mm, 95% confidence interval (CI) of difference over 60 years -1.69 to -0.95]. By contrast, at 25 years mean C.Th of the IA quadrant was 3.9 mm, whereas at 85 years it was 3.3 mm (-0.6 mm, 95% CI -0.83 to -0.10). CtBMD of the IA region was equivalent at 25 and 85 years. In conclusion, elderly women had relative preservation of IA femoral neck bone over seven decades compared with young women but markedly lower C.Th and BMD in the other three quadrants. The IA quadrant transmits mechanical load from walking. Mechanical theory and laboratory tests on cadaveric femurs suggest that localized bone loss may increase the risk of fracture in elderly fallers. It remains to be determined whether this MDCT technique can provide better prediction of hip fracture than conventional clinical dual X-ray absorptiometry (DXA).

Development and testing of texture discriminators for the analysis of trabecular bone in proximal femur radiographs

PURPOSE

Texture analysis of femur radiographs may serve as a potential low cost technique to predict osteoporotic fracture risk and has received considerable attention in the past years. A further application of this technique may be the measurement of the quality of specific bone compartments to provide useful information for treatment of bone fractures. Two challenges of texture analysis are the selection of the best suitable texture measure and reproducible placement of regions of interest (ROIs). The goal of this in vitro study was to automatically place ROIs in radiographs of proximal femur specimens and to calculate correlations between various different texture analysis methods and the femurs' anchorage strength.

METHODS

Radiographs were obtained from 14 femoral specimens and bone mineral density (BMD) was measured in the femoral neck. Biomechanical testing was performed to assess the anchorage strength in terms of failure load, breakaway torque, and number of cycles. Images were segmented using a framework that is based on the usage of level sets and statistical in-shape models. Five ROIs were automatically placed in the head, upper and lower neck, trochanteric, and shaft compartment in an atlas subject. All other subjects were registered rigidly, affinely, and nonlinearly, and the resulting transformation was used to map the five ROIs onto the individual femora.

RESULTS

In each ROI, texture features were extracted using gray level co-occurence matrices (GLCM), third-order GLCM, morphological gradients (MGs), Minkowski dimensions (MDs), Minkowski functionals (MFs), Gaussian Markov random fields, and scaling index method (SIM). Coefficients of determination for each texture feature with parameters of anchorage strength were computed. In a stepwise multiregression analysis, the most predictive parameters were identified in different models. Texture features were highly correlated with anchorage strength estimated by the failure load of up to R2=0.61 (MF and MG features, p<0.01) and were partially independent of BMD. The correlations were dependent on the choice of the ROI and the texture measure. The best predictive multiregression model for failure load R2adj=0.86 (p<0.001) included a set of recently developed texture methods (MF and SIM) but excluded bone mineral density and commonly used texture measures.

CONCLUSIONS

The results suggest that texture information contained in trabecular bone structure visualized on radiographs may predict whether an implant anchorage can be used and may determine the local bone quality from preoperative radiographs.

Assessment of volumetric bone mineral density of the femoral neck in postmenopausal women with and without vertebral fractures using quantitative multi-slice CT

OBJECTIVE

To demonstrate the validity and reliability of volumetric quantitative computed tomography (vQCT) with multi-slice computed tomography (MSCT) and dual energy X-ray absorptiometry (DXA) for hip bone mineral density (BMD) measurements, and to compare the differences between the two techniques in discriminating postmenopausal women with osteoporosis-related vertebral fractures from those without.

METHODS

Ninety subjects were enrolled and divided into three groups based on the BMD values of the lumbar spine and/or the femoral neck by DXA. Groups 1 and 2 consisted of postmenopausal women with BMD changes <-2SD, with and without radiographically confirmed vertebral fracture (n=11 and 33, respectively). Group 3 comprised normal controls with BMD changes > or =-1SD (n=46). Post-MSCT (GE, LightSpeed16) scan reconstructed images of the abdominal-pelvic region, 1.25 mm thick per slice, were processed by OsteoCAD software to calculate the following parameters: volumetric BMD values of trabecular bone (TRAB), cortical bone (CORT), and integral bone (INTGL) of the left femoral neck, femoral neck axis length (NAL), and minimum cross-section area (mCSA). DXA BMD measurements of the lumbar spine (AP-SPINE) and the left femoral neck (NECK) also were performed for each subject.

RESULTS

The values of all seven parameters were significantly lower in subjects of Groups 1 and 2 than in normal postmenopausal women (P<0.05, respectively). Comparing Groups 1 and 2, 3D-TRAB and 3D-INTGL were significantly lower in postmenopausal women with vertebral fracture(s) [(109.8+/-9.61) and (243.3+/-33.0) mg/cm3, respectively] than in those without [(148.9+/-7.47) and (285.4+/-17.8) mg/cm(3), respectively] (P<0.05, respectively), but no significant differences were evident in AP-SPINE or NECK BMD.

CONCLUSION

the femoral neck-derived volumetric BMD parameters using vQCT appeared better than the DXA-derived ones in discriminating osteoporotic postmenopausal women with vertebral fractures from those without. vQCT might be useful to evaluate the effect of osteoporotic vertebral fracture status on changes in bone mass in the femoral neck.

Femoral bone strength and its relation to cortical and trabecular changes after treatment with PTH, alendronate, and their combination as assessed by finite element analysis of quantitative CT scans

The "PTH and Alendronate" or "PaTH" study compared the effects of PTH(1-84) and/or alendronate (ALN) in 238 postmenopausal, osteoporotic women. We performed finite element analysis on the QCT scans of 162 of these subjects to provide insight into femoral strength changes associated with these treatments and the relative roles of changes in the cortical and trabecular compartments on such strength changes. Patients were assigned to either PTH, ALN, or their combination (CMB) in year 1 and were switched to either ALN or placebo (PLB) treatment in year 2: PTH-PLB, PTH-ALN, CMB-ALN, and ALN-ALN (year 1-year 2) treatments. Femoral strength was simulated for a sideways fall using nonlinear finite element analysis of the quantitative CT exams. At year 1, the strength change from baseline was statistically significant for PTH (mean, 2.08%) and ALN (3.60%), and at year 2, significant changes were seen for the PTH-ALN (7.74%), CMB-ALN (4.18%), and ALN-ALN (4.83%) treatment groups but not for PTH-PLB (1.17%). Strength increases were primarily caused by changes in the trabecular density regardless of treatment group, but changes in cortical density and mass also played a significant role, the degree of which depended on treatment. For PTH treatment at year 1 and for ALN-ALN treatment at year 2, there were significant negative and positive strength effects, respectively, associated with a change in external bone geometry. Average changes in strength per treatment group were somewhat consistent with average changes in total hip areal BMD as measured by DXA, except for the PTH group at year 1. The relation between change in femoral strength and change in areal BMD was weak (r(2) = 0.14, pooled, year 2). We conclude that femoral strength changes with these various treatments were dominated by trabecular changes, and although changes in the cortical bone and overall bone geometry did contribute to femoral strength changes, the extent of these latter effects depended on the type of treatment.