Digital wrist tomosynthesis (DWT)-based finite element analysis of ultra-distal radius differentiates patients with and without a history of osteoporotic fracture

Despite effective therapies for those at risk of osteoporotic fracture, low adherence to screening guidelines and limited accuracy of bone mineral density (BMD) in predicting fracture risk preclude identification of those at risk. Because of high adherence to routine mammography, bone health screening at the time of mammography using a digital breast tomosynthesis (DBT) scanner has been suggested as a potential solution. BMD and bone microstructure can be measured from the wrist using a DBT scanner. However, the extent to which biomechanical variables can be derived from digital wrist tomosynthesis (DWT) has not been explored. Accordingly, we measured stiffness from a DWT based finite element (DWT-FE) model of the ultra-distal (UD) radius and ulna, and correlate these to reference microcomputed tomography image based FE (μCT-FE) from five cadaveric forearms. Further, this method is implemented to determine in vivo reproducibility of FE derived stiffness of UD radius and demonstrate the in vivo utility of DWT-FE in bone quality assessment by comparing two groups of postmenopausal women with and without a history of an osteoporotic fracture (Fx; n = 15, NFx; n = 51). Stiffness obtained from DWT and μCT had a strong correlation (R2 = 0.87, p < 0.001). In vivo repeatability error was <5 %. The NFx and Fx groups were not significantly different in DXA derived minimum T-scores (p > 0.3), but stiffness of the UD radius was lower for the Fx group (p < 0.007). Logistic regression models of fracture status with stiffness of the nondominant arm as the predictor were significant (p < 0.01). In conclusion this study demonstrates the feasibility of fracture risk assessment in mammography settings using DWT imaging and FE modeling in vivo. Using this approach, bone and breast screening can be performed in a single visit, with the potential to improve both the prevalence of bone health screening and the accuracy of fracture risk assessment.

Stiffness and Strain Properties Derived From Digital Tomosynthesis-Based Digital Volume Correlation Predict Vertebral Strength Independently From Bone Mineral Density

Vertebral fractures are the most common osteoporotic fractures, but their prediction using standard bone mineral density (BMD) measurements from dual energy X-ray absorptiometry (DXA) is limited in accuracy. Stiffness, displacement, and strain distribution properties derived from digital tomosynthesis-based digital volume correlation (DTS-DVC) have been suggested as clinically measurable metrics of vertebral bone quality. However, the extent to which these properties correlate to vertebral strength is unknown. To establish this relationship, two independent experiments, one examining isolated T11 and the other examining L3 vertebrae within the L2-L4 segments from cadaveric donors were utilized. Following DXA and DTS imaging, the specimens were uniaxially compressed to fracture. BMD, bone mineral content (BMC), and bone area were recorded for the anteroposterior and lateromedial views from DXA, stiffness, endplate to endplate displacement and distribution statistics of intravertebral strains were calculated from DTS-DVC and vertebral strength was measured from mechanical tests. Regression models were used to examine the relationships of strength with the other variables. Correlations of BMD with vertebral strength varied between experimental groups (R2adj = 0.19-0.78). DTS-DVC derived properties contributed to vertebral strength independently from BMD measures (increasing R2adj to 0.64-0.95). DTS-DVC derived stiffness was the best single predictor (R2adj = 0.66, p < 0.0001) and added the most to BMD in models of vertebral strength for pooled T11 and L3 specimens (R2adj = 0.95, p < 0.0001). These findings provide biomechanical relevance to DTS-DVC calculated properties of vertebral bone and encourage further efforts in the development of the DTS-DVC approach as a clinical tool.

Intervertebral kinematics during neck motion 6.5 years after fusion and artificial disc replacement

BACKGROUND

Arthroplasty with artificial disc replacement for surgical treatment of cervical spine degeneration was introduced with the notion that motion-preserving approaches would prevent development of adjacent segment disease. Though clinical outcomes favor arthroplasty over the commonly used anterior cervical discectomy with fusion approach, clinical studies confirming the biomechanical basis of these results are lacking. The aim of this study was to compare intervertebral kinematics between arthroplasty and fusion patients 6.5 years post-surgery during physiological motion of the neck.

METHODS

Using a biplane dynamic X-ray system, computed tomography imaging and model based tracking algorithms, three dimensional intervertebral kinematics were measured during neck axial rotation and extension in 14 patients treated for cervical radiculopathy with fusion (n = 8) or arthroplasty (n = 6). The measurements were performed at 2-year (baseline) and 6.5 year post-surgical time points, with the main interest being in the interaction between surgery types and time points. 3 translations and 3 rotations were investigated for the index (C5C6), and upper- (C4C5) and lower adjacent levels (C6C7).

FINDINGS

Surgery-time interaction was significant for axial rotation (P < 0.04) and flexion-extension rotation (P < 0.005) in C4C5 during neck axial rotation, left-right translation (P < 0.04) in C5C6 and anterior-posterior translation in C6C7 (P < 0.04) during neck extension. In contrast with the expectations, axial rotation and flexion-extension decreased in C4C5 during neck rotation and anterior-posterior translation decreased in C6C7 during neck extension for fusion.

INTERPRETATION

The findings do not support the notion that adjacent segment motion increases after fusion.

Measuring the thickness of vertebral endplate and shell using digital tomosynthesis

The vertebral endplate and cortical shell play an important structural role and contribute to the overall strength of the vertebral body, are at highest risk of initial failure, and are involved in degenerative disease of the spine. The ability to accurately measure the thickness of these structures is therefore important, even if difficult due to relatively low resolution clinical imaging. We posit that digital tomosynthesis (DTS) may be a suitable imaging modality for measurement of endplate and cortical shell thickness owing to the ability to reconstruct multiplanar images with good spatial resolution at low radiation dose. In this study, for 25 cadaveric L1 vertebrae, average and standard deviation of endplate and cortical shell thickness were measured using images from DTS and microcomputed tomography (μCT). For endplate thickness measurements, significant correlations between DTS and μCT were found for all variables when comparing thicknesses measured in both the overall endplate volume (R² = 0.25-0.54) and when measurements were limited to a central range of coronal or sagittal slices (R² = 0.24-0.62). When compared to reference values from the overall shell volume, DTS thickness measurements were generally nonsignificant. However, when measurement of cortical shell thickness was limited to a range of central slices, DTS outcomes were significantly correlated with reference values for both sagittal and coronal central regions (R² = 0.21-0.49). DTS may therefore offer a means for measurement of endplate thickness and, within a limited sagittal or coronal measurement volume, for measurement of cortical shell thickness.

The Relationship of Whole Human Vertebral Body Creep to Bone Density and Texture via Clinically Available Imaging Modalities

Creep deformation of human vertebrae accumulates under physiological levels of load and is understood to contribute to the progression toward clinically observable vertebral fracture. However, little information is available in terms of clinically measurable predictors of creep behavior in human vertebrae. In this study, creep tests were performed on 22 human cadaveric T12 vertebrae (13 male, 9 female; age 41-90). Areal and volumetric bone density parameters were measured from the same specimens using dual x-ray absorptiometry and high resolution computed tomography. Image textural analyses (which probe the organization of image intensities within the cancellous bone in low resolution clinical imaging) were performed using digital tomosynthesis (DTS) images. Multiple regression models were constructed to examine the relationship between creep properties and bone density and DTS image textural parameters. For the standard clinical imaging configuration, models including DTS derived image textural parameters alone were generally more explanatory (adjusted R²: 0.14-0.68) than those with bone density parameters forced in the models (adjusted R²: 0.17-0.61). Metrics of textural heterogeneity and anisotropy presented as the most explanatory imaging markers for creep deformation and recovery from creep. These metrics of image texture may help provide, independent from bone mass, important clinically measurable indicators of the time dependent deformation of human vertebrae.

Assessment of Intravertebral Mechanical Strains and Cancellous Bone Texture Under Load Using a Clinically Available Digital Tomosynthesis Modality

Vertebral fractures are the most common osteoporotic fractures, but clinical means for assessment of vertebral bone integrity are limited in accuracy, as they typically use surrogate measures that are indirectly related to mechanics. The objective of this study was to examine the extent to which intravertebral strain distributions and changes in cancellous bone texture generated by a load of physiological magnitude can be characterized using a clinically available imaging modality. We hypothesized that digital tomosynthesis-based digital volume correlation (DTS-DVC) and image texture-based metrics of cancellous bone microstructure can detect development of mechanical strains under load. Isolated cadaveric T11 vertebrae and L2-L4 vertebral segments were DTS imaged in a nonloaded state and under physiological load levels. Axial strain, maximum principal strain, maximum compressive and tensile principal strains, and von Mises equivalent strain were calculated using the DVC technique. The change in textural parameters (line fraction deviation, anisotropy, and fractal parameters) under load was calculated within the cancellous centrum. The effect of load on measured strains and texture variables was tested using mixed model analysis of variance, and relationships of strain and texture variables with donor age, bone density parameters, and bone size were examined using regression models. Magnitudes and heterogeneity of intravertebral strain measures correlated with applied loading and were significantly different from background noise. Image texture parameters were found to change with applied loading, but these changes were not observed in the second experiment testing L2-L4 segments. DTS-DVC-derived strains correlated with age more strongly than did bone mineral density (BMD) for T11.

Bone health assessment via digital wrist tomosynthesis in the mammography setting

Bone fractures attributable to osteoporosis are a significant problem. Though preventative treatment options are available for individuals who are at risk of a fracture, a substantial number of these individuals are not identified due to lack of adherence to bone screening recommendations. The issue is further complicated as standard diagnosis of osteoporosis is based on bone mineral density (BMD) derived from dual energy x-ray absorptiometry (DXA), which, while helpful in identifying many at risk, is limited in fully predicting risk of fracture. It is reasonable to expect that bone screening would become more prevalent and efficacious if offered in coordination with digital breast tomosynthesis (DBT) exams, provided that osteoporosis can be assessed using a DBT modality. Therefore, the objective of the current study was to explore the feasibility of using digital tomosynthesis imaging in a mammography setting. To this end, we measured density, cortical thickness and microstructural properties of the wrist bone, correlated these to reference measurements from microcomputed tomography and DXA, demonstrated the application in vivo in a small group of participants, and determined the repeatability of the measurements. We found that measurements from digital wrist tomosynthesis (DWT) imaging with a DBT scanner were highly repeatable ex vivo (error = 0.05%-9.62%) and in vivo (error = 0.06%-10.2%). In ex vivo trials, DWT derived BMDs were strongly correlated with reference measurements (R = 0.841-0.980), as were cortical thickness measured at lateral and medial cortices (R = 0.991 and R = 0.959, respectively) and the majority of microstructural measures (R = 0.736-0.991). The measurements were quick and tolerated by human patients with no discomfort, and appeared to be different between young and old participants in a preliminary comparison. In conclusion, DWT is feasible in a mammography setting, and informative on bone mass, cortical thickness, and microstructural qualities that are known to deteriorate in osteoporosis. To our knowledge, this study represents the first application of DBT for imaging bone. Future clinical studies are needed to further establish the efficacy for diagnosing osteoporosis and predicting risk of fragility fracture using DWT.

Dynamic foraminal dimensions during neck motion 6.5 years after fusion and artificial disc replacement

Objective

To compare changes in foraminal motion at two time points post-surgery between artificial disc replacement (ADR) and anterior cervical discectomy and fusion (ACDF).

Methods

Eight ACDF and 6 ADR patients (all single-level C5-6) were tested at 2 years (T1) and 6.5 years (T2) post-surgery. The minimum foraminal height (FH.Min) and width (FW.Min) achieved during neck axial rotation and extension, and the range of these dimensions during motion (FH.Rn and FW.Rn, respectively) were measured using a biplane dynamic x-ray system, CT imaging and model-based tracking while patients performed neck axial rotation and extension tasks. Two-way mixed ANOVA was employed for analysis.

Results

In neck extension, significant interactions were found between year post-surgery and type of surgery for FW.Rn at C5-6 (p<0.006) and C6-7 (p<0.005), and for FH.Rn at C6-7 (p<0.01). Post-hoc analysis indicated decreases over time in FW.Rn for ACDF (p<0.01) and increases in FH.Rn for ADR (p<0.03) at the C6-7 adjacent level. At index level, FW.Rn was comparable between ACDF and ADR at T1, but was smaller for ACDF than for ADR at T2 (p<0.002). In axial rotation, differences were found between T1 and T2 but did not depend on type of surgery (p>0.7).

Conclusions

Changes were observed in the range of foraminal geometry at adjacent levels from 2 years to 6.5 years post-surgery that were different between ACDF and ADR for neck extension. These changes are contrary to the notion that motion at adjacent levels continue to increase following ACDF as compared to ADR over the long term.

Correlation of neural foraminal motion after surgical treatment of cervical radiculopathy with long-term patient reported outcomes

Background

Post-surgical changes in adjacent segment motion are considered a factor in further development of degeneration and cervical radiculopathy. The objective was to examine the extent of correlations between physiological motion of cervical foramina and long-term patient reported outcomes (PRO).

Methods

Biplane X-ray imaging and CT-based markerless tracking were used to measure 3D static and dynamic dimensions during neck axial rotation and extension from 18 patients treated for C5-6 radiculopathy with fusion or arthroplasty. Minimum foraminal height (FH.Min) and width (FW.Min), and their range (FH.Range and FW.Range) achieved during a motion task were calculated for adjacent levels (C4-5 and C6-7) at 2.0±0.6 years post-surgery. The modified Japanese Orthopedic Association score (mJOAS), the Neck Disability Index (NDI) including the visual analogue scale (VAS) for neck and arm pain, and the EuroQol EQ-5D score were recorded at 6.5±1.1 years post-surgery. The relationships between 6.5-year outcomes and 2-year foraminal motion were examined using regression.

Results

Worsening patient-reported outcomes were generally associated with lower values of FW.Min (P<0.05 to P<0.008), the associations being stronger for neck extension (r² up to 0.43). Dynamic foraminal measurements from the C6-7 level more significantly and consistently correlated with mJOAS, EQ-5D and NDI Arm Pain VAS (r²=0.27 to 0.43; P<0.03 to P<0.008), whereas those from the C4-5 level correlated with NDI Neck Pain VAS (r²=0.33; P<0.02).

Conclusions

Dynamic 3D foraminal dimensions at 2-year post-surgery, notably FW.Min measured in neck extension at adjacent levels, were associated with PRO at 6.5 years post-surgery. These relationships provide insight into the motion related factors in development of pain and loss of function, and may help develop markers or objective outcome measures.

Digital tomosynthesis based digital volume correlation: A clinically viable noninvasive method for direct measurement of intravertebral displacements using images of the human spine under physiological load

PURPOSE

We have developed a clinically viable method for measurement of direct, patient-specific intravertebral displacements using a novel digital tomosynthesis based digital volume correlation technique. These displacements may be used to calculate vertebral stiffness under loads induced by a patient's body weight; this is particularly significant because, among biomechanical variables, stiffness is the strongest correlate of bone strength. In this proof of concept study, we assessed the feasibility of the method through a preliminary evaluation of the accuracy and precision of the method, identification of a range of physiological load levels for which displacements are measurable, assessment of the relationship of measured displacements with microcomputed tomography based standards, and demonstration of the in vivo application of the technique.

METHODS

Five cadaveric T11 vertebrae were allocated to three groups in order to study (a) the optimization of digital volume correlation algorithm input parameters, (b) accuracy and precision of the method and the ability to measure displacements at a range of physiological load levels, and (c) the correlation between displacements measured using tomosynthesis based digital volume correlation vs. high resolution microcomputed tomography based digital volume correlation and large scale finite element models. Tomosynthesis images of one patient (Female, 60 yr old) were used to calculate displacement maps, and in turn stiffness, using images acquired in both standing and standing-with-weight (8 kg) configurations.

RESULTS

We found that displacements were accurate (2.28 µm total error) and measurable at physiological load levels (above 267 N) with a linear response to applied load. Calculated stiffness among three tested vertebral bodies was within an acceptable range relative to reported values for vertebral stiffness (5651-13260 N/mm). Displacements were in good qualitative and quantitative agreement with both microcomputed tomography based finite element (r²  = 0.762, P < 0.001) and digital volume correlation (r²  = 0.799, P < 0.001) solutions. For one patient tested twice, once standing and once holding weights, results demonstrated excellent qualitative reproducibility of displacement distributions with superior endplate displacements increasing by 22% with added weight.

CONCLUSIONS

The results of this work collectively suggest the feasibility of the method for in vivo measurement of intravertebral displacements and stiffness in humans. These findings suggest that digital volume correlation using digital tomosynthesis imaging may be useful in understanding the mechanical response of bone to disease and may further enhance our ability to assess fracture risk and treatment efficacy for the spine.

Assessment of vertebral wedge strength using cancellous textural properties derived from digital tomosynthesis and density properties from dual energy X-ray absorptiometry and high resolution computed tomography

The purpose of this study was to examine the potential of digital tomosynthesis (DTS) derived cancellous bone textural measures to predict vertebral strength under conditions simulating a wedge fracture. 40 vertebral bodies (T6, T8, T11, and L3 levels) from 5 male and 5 female cadaveric donors were utilized. The specimens were scanned using dual energy X-ray absorptiometry (DXA) and high resolution computed tomography (HRCT) to obtain measures of bone mineral density (BMD) and content (BMC), and DTS to obtain measures of bone texture. Using a custom loading apparatus designed to deliver a nonuniform displacement resulting in a wedge deformity similar to those observed clinically, the specimens were loaded to fracture and their fracture strength was recorded. Mixed model regressions were used to determine the associations between wedge strength and DTS derived textural variables, alone and in the presence of BMD or BMC information. DTS derived fractal, lacunarity and mean intercept length variables correlated with wedge strength, and individually explained up to 53% variability. DTS derived textural variables, notably fractal dimension and lacunarity, contributed to multiple regression models of wedge strength independently from BMC and BMD. The model from a scan orientation transverse to the spine axis and in the anterior-posterior view resulted in highest explanatory capability (R²_adj = 0.91), with a scan orientation parallel to the spine axis and in the lateral view offering an alternative (R²_adj = 0.88). In conclusion, DTS can be used to examine cancellous texture relevant to vertebral wedge strength, and potentially complement BMD in assessment of vertebral fracture risk.

Dynamic foraminal dimensions during neck extension and rotation in fusion and artificial disc replacement: an observational study

BACKGROUND

Changes in the dimensions of the cervical neural foramina (CNF) are considered to be a key factor in nerve root compression and development of cervical radiculopathy. However, to what extent foraminal geometry differs between patients who underwent anterior cervical discectomy and fusion (ACDF) and those who underwent total disc arthroplasty with an artificial disc (AD) during physiological motion is largely unknown.

PURPOSE

The objective of this study is to compare CNF dimensions during physiological neck motion between ACDF and AD.

STUDY DESIGN/SETTING

This is a retrospective comparative analysis of prospectively collected, consecutive, non-randomized series of patients at a single institution.

PATIENT SAMPLE

A total of 16 single-level C5-C6 ACDF (4 males, 12 females; 28-71 years) and 7 single-level C5-C6 cervical arthroplasty patients (3 males, 4 females; 38-57 years), at least 12 months after surgery (23.6±6.8 months) were included.

OUTCOME MEASURES

Patient demographics, preoperative magnetic resonance imaging (MRI)-based measurements of cervical spine degeneration, and 2-year postoperative measurements of dynamic foraminal geometry were the outcome measures.

METHODS

Biplane X-ray images were acquired during axial neck rotation and neck extension. A computed tomography scan was also acquired from C3 to the first thoracic vertebrae. The subaxial cervical vertebrae (C3-C7) were reconstructed into three-dimensional (3D) bone models for use with model-based tracking. Foraminal height (FH) was calculated as the 3D distance between the superior point of the inferior pedicle and the inferior point of the superior pedicle using custom software. Foraminal width (FW) was similarly calculated as the 3D distance between the anterolateral aspect of the superior vertebral body inferior notch and the posterolateral aspect of the inferior vertebral body superior notch. Dynamic foraminal dimensions were quantified as the minimum (FH.Min, FW.Min), the range (FH.Range, FW.Range), and the median (FH.Med, FW.Med) of each trial and then averaged over trials. Mixed model analysis of variance framework was used to examine the differences between ACDF and AD groups. The initial severity of disc degeneration as determined from preoperative MRI images was introduced as covariates in the models.

RESULTS

At the operated level (C5-C6), FH.Med and FH.Range were smaller in ACDF than in AD during axial rotation and neck extension (p<.003 to p<.05). At the superior adjacent level (C4-C5), no significant difference was found. At the inferior adjacent level (C6-C7), FW.Range was greater in ACDF than in AD during axial rotation and extension (p<.05). At the non-adjacent level (C3-C4), FW.Range was greater in ACDF than in AD during extension (p<.008).

CONCLUSIONS

This study demonstrated decreases in foraminal dimensions and their range for ACDF compared with AD at the operated level. In contrast, it demonstrated increases in the range of foraminal dimensions during motion for ACDF compared with AD at the non-operated segments. Together, these data support the notion that increased mobility at the non-operated segments after ACDF may contribute to a greater risk for adjacent segment degeneration. Because of the significant presence of range variables in the findings, the current data also indicate that a dynamic evaluation is likely more appropriate for evaluation of the differences in foramina between ACDF and AD than a static evaluation.

The relationship of whole human vertebral body creep to geometric, microstructural, and material properties

Creep, the time dependent deformation of a structure under load, is an important viscoelastic property of bone and may play a role in the development of permanent deformity of the vertebrae in vivo leading to clinically observable spinal fractures. To date, creep properties and their relationship to geometric, microstructural, and material properties have not been described in isolated human vertebral bodies. In this study, a range of image-based measures of vertebral bone geometry, bone mass, microarchitecture and mineralization were examined in multiple regression models in an effort to understand their contribution to creep behavior. Several variables, such as measures of mineralization heterogeneity, average bone density, and connectivity density persistently appeared as significant effects in multiple regression models (adjusted r²: 0.17-0.56). Although further work is needed to identify additional tissue properties to fully describe the portion of variability not explained by these models, these data are expected to help understand mechanisms underlying creep and improve prediction of vertebral deformities that eventually progress to a clinically observable fracture.

A biomechanical comparison of subscapularis repair techniques in total shoulder arthroplasty: lesser tuberosity osteotomy versus subscapularis peel

Background

The subscapularis peel (SP) and the lesser tuberosity osteotomy (LTO) are 2 common exposure techniques for total shoulder arthroplasty. Although some biomechanical studies have suggested a higher resistance to failure with the LTO, clinical studies have demonstrated no difference in repair failure or tendon healing. We hypothesized that there would be no difference in biomechanically tested repair strength between our SP technique and the previously tested LTO technique.

Methods

Eleven cadaver shoulders were separated into 2 groups: 6 SPs and 5 LTOs. After initial loading for 3000 cycles, the specimens were incrementally loaded to 450 ± 50 N or catastrophic failure. Repair gapping was measured after cyclical loading, and fatigue life was analyzed after incremental loading.

Results

There was no significant difference in mean repair gapping between the SP (2.40 ± 0.36 mm; mean ± standard deviation) and the LTO groups (3.10 ± 2.93 mm; P = .57). There was also no difference in the mean number of cycles to failure (6894 ± 956 vs. 6018 ± 1179; P = .14) and mean load to failure (400 ± 79 N vs. 340 ± 91 N; P = .21) between the SP and LTO techniques. However, there was more variability in bead gapping in the LTO group (P < .01).

Conclusion

No significant differences were found in repair gapping, fatigue failure, and load to failure in comparing the SP and LTO repairs. However, the SP repair demonstrated significantly less variability in repair gapping. These findings suggest that initial fixation biomechanical properties between the 2 constructs are similar in vitro.

A Point-of-Care Raman Spectroscopy-Based Device for the Diagnosis of Gout and Pseudogout: Comparison With the Clinical Standard Microscopy

OBJECTIVE

To demonstrate the usefulness of a novel medical device based on Raman spectroscopy for the rapid point-of-care diagnosis of gout and pseudogout.

METHODS

A shoebox-sized point-of-care Raman spectroscopy (POCRS) device was developed for use in the diagnosis of gout and pseudogout. The device included a disposable syringe microfiltration kit to collect arthropathic crystals from synovial fluid and a customized automated Raman spectroscopy system to chemically identify crystal species. Diagnosis according to the findings of POCRS was compared with the clinical standard diagnosis based on compensated polarized light microscopy (CPLM) of synovial fluid aspirates collected from symptomatic patients (n = 174). Kappa coefficients were used to measure the agreement between POCRS and CPLM findings.

RESULTS

Overall, POCRS and CPLM results were consistent in 89.7% of samples (156 of 174). For the diagnosis of gout, the kappa coefficient for POCRS and CPLM was 0.84 (95% confidence interval [95% CI] 0.75-0.94). For the diagnosis of pseudogout, the kappa coefficient for POCRS and CPLM was 0.61 (95% CI 0.42-0.81).

CONCLUSION

Kappa coefficients indicated that there was excellent agreement between POCRS and CPLM for the diagnosis of gout, with good agreement for the diagnosis of pseudogout. The POCRS device holds the potential to standardize and expedite the time to clinical diagnosis of gout and pseudogout, especially in settings where certified operators trained for CPLM analysis are not available.

The Use of Tomosynthesis in the Global Study of Knee Subchondral Insufficiency Fractures

RATIONALE AND OBJECTIVES

Subchondral insufficiency fractures (SIF), previously termed spontaneous osteonecrosis of the knee, are marked by a sudden onset of severe pain. Other than the size of the lesion, prediction for progression to joint replacement is difficult. The objective was to determine if quantitative analysis of bone texture using digital tomosynthesis imaging would be useful in predicting more rapid progression to joint replacement.

MATERIALS AND METHODS

Tomosynthesis studies of 30 knees with documented SIF were quantified by fractal, mean intercept length (MIL), and line fraction deviation analyses. Fractal dimension, lacunarity, MIL, and line fraction deviation variables measured from these analyses were then correlated to short interval progression to joint replacement surgery.

RESULTS

Higher odds for joint replacement were related to higher values of the standard deviation of slope lacunarity and to morphometric measures (eg, MIL).

CONCLUSIONS

Using digital tomosynthesis images for bone texture assessment may help distinguish condylar bone response in SIF, potentially acting as a clinically relevant predictive tool. In the future, contrasting SIF to the more gradual long-term process of osteoarthritis, there may be a better understanding of the different mechanisms for the two conditions.

Digital tomosynthesis and high resolution computed tomography as clinical tools for vertebral endplate topography measurements: Comparison with microcomputed tomography

Endplate morphology is understood to play an important role in the mechanical behavior of vertebral bone as well as degenerative processes in spinal tissues; however, the utility of clinical imaging modalities in assessment of the vertebral endplate has been limited. The objective of this study was to evaluate the ability of two clinical imaging modalities (digital tomosynthesis, DTS; high resolution computed tomography, HRCT) to assess endplate topography by correlating the measurements to a microcomputed tomography (μCT) standard. DTS, HRCT, and μCT images of 117 cadaveric thoracolumbar vertebrae (T10-L1; 23 male, 19 female; ages 36-100 years) were segmented, and inferior and superior endplate surface topographical distribution parameters were calculated. Both DTS and HRCT showed statistically significant correlations with μCT approaching a moderate level of correlation at the superior endplate for all measured parameters (R(2)Adj=0.19-0.57), including averages, variability, and higher order statistical moments. Correlation of average depths at the inferior endplate was comparable to the superior case for both DTS and HRCT (R(2)Adj=0.14-0.51), while correlations became weak or nonsignificant for higher moments of the topography distribution. DTS was able to capture variations in the endplate topography to a slightly better extent than HRCT, and taken together with the higher speed and lower radiation cost of DTS than HRCT, DTS appears preferable for endplate measurements.

Mechanical loading causes detectable changes in morphometric measures of trabecular structure in human cancellous bone

The relationships between mechanical loads and bone microstructure are of interest to those who seek to predict bone mechanical properties from microstructure or to predict how organization of bone microstructure is driven by mechanical loads. While strains and displacements in the material are inherently responsible for mechanically caused changes in the appearance of the microstructure, it is the morphometric measures of microstructural organization that are often available for assessment of bone quality. Therefore, an understanding of how strain history is reflected in morphometric measures of bone microstructure has practical implications in that it may provide clinically measurable indices of mechanical history in bone and improve interpretation of bone mechanical properties from microstructural information. The objective of the current study was to examine changes in morphometric measures of cancellous bone microstructure in response to varying levels of continuum level strains. The experimental approach included stereologic analysis of microcomputed tomography (μCT) images of human cancellous bone samples obtained at sequentially increasing levels of strain in a custom-made loading apparatus mounted in a μCT scanner. We found that the degree of anisotropy (DA) decreased from baseline to failure and from failure to postfailure. DA partially recovered from postfailure levels upon unloading; however, the final DA was less than at failure and less than at baseline. We also found that average trabecular thickness (Tb.Th.Av) increased with displacements at postfailure and did not recover when unloaded. Average trabecular number decreased when the specimens were unloaded. In addition, the heterogeneity of Tb.Th as measured by intra-specimen standard deviation (Tb.Th.SD) increased and that of trabecular number (Tb.N.SD) decreased with displacements at postfailure. Furthermore, the intraspecimen coefficient of variation of trabecular number decreased at postfailure displacements but did not recover upon unloading. Finally, the coefficient of variation of trabecular separation at unload was less than that at baseline. These measures can be developed into image-based indices to estimate strain history, damage, and residual mechanical properties where direct analysis of stresses and strains, such as through finite element modeling, may not be feasible. It remains to be determined how wide a time interval can be used to estimate strain history before remodeling becomes an overriding effect on the trabecular architecture.

Laser Wavelength Dependence of Background Fluorescence in Raman Spectroscopic Analysis of Synovial Fluid from Symptomatic Joints

Gout is a disease process where the nucleation and growth of crystals in the synovial fluid of joints elicit painful arthritis-like symptoms. Raman spectroscopy is evolving as a potential diagnostic tool in identifying such crystals; however, attainment of sufficient Raman signal while overcoming the background fluorescence remains as a major challenge. The current study focused on assessing whether excitation in 532-700 nm range will provide greater signal intensity than the standard 785 nm while not being impeded by background fluorescence. We characterized the fluorescence spectra, absorption spectra and Raman spectra of synovial fluid from patients who presented "gout-like symptoms" (symptomatic) and controls (asymptomatic). A digestion and filtration method was developed to isolate crystals from synovial fluid while reducing the organic burden. Spectral profile and photobleaching dynamics during Raman spectroscopy were observed under an excitation wavelength range spanning 532 to 785 nm. Absorbance and fluorescence profiles indicated the digestion and filtration worked effectively to extract crystals from symptomatic synovial fluid without introducing additional fluorescence. Raman spectral analyses at 532 nm, 660 nm, 690 nm and 785 nm indicated that both asymptomatic and symptomatic samples had significant levels of fluorescence at excitation wavelengths below 700 nm, which either hindered the collection of Raman signal or necessitated prolonged durations of photobleaching. Raman-based diagnostics were more feasible at the longest excitation wavelength of 785 nm without employing photobleaching. This study further demonstrated that a near-infrared OEM based lower-cost Raman system at 785 nm excitation has sufficient sensitivity to identify crystals isolated from the synovial fluid. In conclusion, while lower excitation wavelengths provide greater signal, the fluorescence necessitates near-infrared wavelengths for Raman analysis of crystal species observed in synovial aspirates.

Heterogeneity of bone mineral density and fatigue failure of human vertebrae

There is increasing interest in using the heterogeneity of tissue properties in a bone for predicting its fracture risk. Heterogeneity of volumetric bone mineral density (BMD) as measured from quantitative computed tomography (QCT) is of particular interest as these measurements are clinically feasible. Previous examinations of the relationship between the BMD heterogeneity and the mechanical behavior of human vertebrae only considered quasistatic strength and were with limited number of samples. McCubbrey et al. (1995) studied the value of regional BMDs for predicting vertebral fatigue life, determined from short-cycle tests at force levels scaled with the estimated strength of the vertebra, but the focus of that work was in best predictor subsets without a specific focus on the heterogeneity of BMD or the positive vs negative direction of the relationships. The previous analysis also did not take into account the censored nature of the fatigue life data. As such, whether BMD heterogeneity is positively or negatively associated with fatigue life and whether this is independent of the average or minimum BMD are not clear. In the present work, we revisited the McCubbrey data for a preliminary examination of the relationship between BMD heterogeneity and fatigue life using survival analysis. The analysis suggests that BMD heterogeneity measured as the intra-vertebral standard deviation of BMDs in a vertebra is negatively associated with short cycle (high-amplitude) fatigue life independent of the average BMD. The results motivate further studies on the role of BMD heterogeneity in fatigue failure and clinical fracture risk of human vertebrae.

The relationships between femoral cortex geometry and tissue mechanical properties

Bone tissue and geometry are constantly modified through modeling and remodeling at the periosteal, endosteal and intracortical envelopes. Results from several studies indicate that femoral bone geometry is a predictor of whole bone strength (e.g. femoral neck strength), however, it is not known whether there is a relationship between bone structural and material properties. Bone geometry can be determined from parameters based on plane X-ray radiogrammetry which are used to evaluate femoral bone quality for implant success. If there is a relationship between these parameters and tissue mechanical properties, this would have implications in the interpretation of such parameters for assessment of fracture risk and in further understanding of bone biology. Following measurement of radiogrammetric parameters from antero-posterior and medio-lateral X-rays (cortical thickness, bone diameter, bone area, moment of inertia, cortical index, Singh index), human femurs were machined into standard test specimens for assessment of tensile fracture toughness (GIc) of the tissue. Results indicated that tensile fracture toughness generally increased with increasing bone size. We also found that fracture toughness of the tissue was significantly related to radiogrammetric indices and that some of these indices explained a greater variability in toughness than porosity, age or gender.

Variability of trabecular microstructure is age-, gender-, race- and anatomic site-dependent and affects stiffness and stress distribution properties of human vertebral cancellous bone

Cancellous bone microstructure is an important determinant of the mechanical integrity of vertebrae. The numerous microstructural parameters that have been studied extensively are generally represented as a single value obtained as an average over a sample. The range of the intra-sample variability of cancellous microstructure and its effect on the mechanical properties of bone are less well-understood. The objectives of this study were to investigate the extent to which human cancellous bone microstructure within a vertebra i) is related to bone modulus and stress distribution properties and ii) changes along with age, gender and locations thoracic 12 (T12) vs lumbar 1 (L1). Vertebrae were collected from 15 male (66±15 years) and 25 female (54±16 years) cadavers. Three dimensional finite element models were constructed using microcomputed tomography images of cylindrical specimens. Linear finite element models were used to estimate apparent modulus and stress in the cylinders during uniaxial compression. The intra-specimen mean, standard deviation (SD) and coefficient of variation (CV) of microstructural variables were calculated. Mixed model statistical analysis of the results demonstrated that increases in the intra-specimen variability of the microstructure contribute to increases in the variability of trabecular stresses and decreases in bone stiffness. These effects were independent from the contribution from intra-specimen average of the microstructure. Further, the effects of microstructural variability on bone stiffness and stress variability were not accounted for by connectivity and anisotropy. Microstructural variability properties (SD, CV) generally increased with age, were greater in females than in males and in T12 than in L1. Significant interactions were found between age, gender, vertebra and race. These interactions suggest that microstructural variability properties varied with age differently between genders, races and vertebral levels. The current results collectively demonstrate that microstructural variability has a significant effect on mechanical properties and tissue stress of human vertebral cancellous bone. Considering microstructural variability could improve the understanding of bone fragility and improve assessment of vertebral fracture risk.

Stiffness of the endplate boundary layer and endplate surface topography are associated with brittleness of human whole vertebral bodies

Stress magnitude and variability as estimated from large scale finite element (FE) analyses have been associated with compressive strength of human vertebral cancellous cores but these relationships have not been explored for whole vertebral bodies. In this study, the objectives were to investigate the relationship of FE-calculated stress distribution parameters with experimentally determined strength, stiffness, and displacement based ductility measures in human whole vertebral bodies, investigate the effect of endplate loading conditions on vertebral stiffness, strength, and ductility and test the hypothesis that endplate topography affects vertebral ductility and stress distributions. Eighteen vertebral bodies (T6-L3 levels; 4 female and 5 male cadavers, aged 40-98 years) were scanned using a flat-panel CT system and followed with axial compression testing with Wood's metal as filler material to maintain flat boundaries between load plates and specimens. FE models were constructed using reconstructed CT images and filler material was added digitally. Two different FE models with different filler material modulus simulating Wood's metal and intervertebral disc (W-layer and D-layer models) were used. Element material modulus to cancellous bone was based on image gray value. Average, standard deviation, and coefficient of variation of von Mises stress in vertebral bone for W-layer and D-layer models and also the ratios of FE parameters from the two models (W/D) were calculated. Inferior and superior endplate surface topographical distribution parameters were calculated. Experimental stiffness, maximum load and work to fracture had the highest correlation with FE-calculated stiffness while experimental ductility measures had highest correlations with FE-calculated average von Mises stress and W-layer to D-layer stiffness ratio. Endplate topography of the vertebra was also associated with its structural ductility and the distribution parameter that best explained this association was kurtosis of inferior endplate topography. Our results indicate that endplate topography variations may provide insight into the mechanisms responsible for vertebral fractures.

Increased microstructural variability is associated with decreased structural strength but with increased measures of structural ductility in human vertebrae

The lack of accuracy in the prediction of vertebral fracture risk from average density measurements, all external factors being equal, may not just be because bone mineral density (BMD) is less than a perfect surrogate for bone strength but also because strength alone may not be sufficient to fully characterize the structural failure of a vertebra. Apart from bone quantity, the regional variation of cancellous architecture would have a role in governing the mechanical properties of vertebrae. In this study, we estimated various microstructural parameters of the vertebral cancellous centrum based on stereological analysis. An earlier study indicated that within-vertebra variability, measured as the coefficient of variation (COV) of bone volume fraction (BV/TV) or as COV of finite element-estimated apparent modulus (E(FE)) correlated well with vertebral strength. Therefore, as an extension to our earlier study, we investigated (i) whether the relationships of vertebral strength found with COV of BV/TV and COV of E(FE) could be extended to the COV of other microstructural parameters and microcomputed tomography-estimated BMD and (ii) whether COV of microstructural parameters were associated with structural ductility measures. COV-based measures were more strongly associated with vertebral strength and ductility measures than average microstructural measures. Moreover, our results support a hypothesis that decreased microstructural variability, while associated with increased strength, may result in decreased structural toughness and ductility. The current findings suggest that variability-based measures could provide an improvement, as a supplement to clinical BMD, in screening for fracture risk through an improved prediction of bone strength and ductility. Further understanding of the biological mechanisms underlying microstructural variability may help develop new treatment strategies for improved structural ductility.

Cancellous bone properties and matrix content of TGF-beta2 and IGF-I in human tibia: a pilot study

Transforming and insulin-like growth factors are important in regulating bone mass. Thus, one would anticipate correlations between matrix concentrations of growth factors and functional properties of bone. We therefore investigated the relationships of (1) TGF-beta2 and (2) IGF-I matrix concentrations with the trabecular microstructure, stress distribution, and mechanical properties of tibial cancellous bone from six male human cadavers. Trabecular stress amplification (VMExp/sigma(app)) and variability (VMCOV) were calculated using microcomputed tomography (muCT)-based finite element simulations. Bone volume fraction (BV/TV), surface/volume ratio (BS/BV), trabecular thickness (Tb.Th), number (Tb.N) and separation (Tb.Sp), connectivity (Eu.N), and anisotropy (DA) were measured using 3-D morphometry. Bone stiffness and strength were measured by mechanical testing. Matrix concentrations of TGF-beta2 and IGF-I were measured by ELISA. We found higher matrix concentrations of TGF-beta2 were associated with higher Tb.Sp and VMExp/sigma(app) for pooled data and within subjects. Similarly, a higher matrix concentration of IGF-I was associated with lower stiffness, strength, BV/TV and Tb.Th and with higher BS/BV, Tb.Sp, VMExp/sigma(app) and VMCOV for pooled data and within subjects. IGF-I and Tb.N were negatively associated within subjects. It appears variations of the stress distribution in cancellous bone correlate with the variation of the concentrations of TGF-beta2 and IGF-I in bone matrix: increased local matrix concentrations of growth factors are associated with poor biomechanical and architectural properties of tibial cancellous bone.

Analytical Approach to Recovering Bone Porosity From Effective Complex Shear Modulus

This work deals with the study of the analytical relations between porosity of cancellous bone and its mechanical properties. The Stieltjes representation of the effective shear complex modulus of cancellous bone is exploited to recover porosity. The microstructural information is contained in the spectral measure in this analytical representation. The spectral function can be recovered from the effective measurements over a range of frequencies. The problem of reconstruction of the spectral measure is very ill-posed. Regularized algorithm is derived to ensure stability of the results. The proposed method does not use any specific assumptions about the microgeometry of bone. The approach does not rely on correlation analysis, it uses analytical relationships. For validation purposes, complex shear modulus over a range of frequencies was calculated by the finite element method using micro-computed tomography (micro-CT) images of human cancellous bone. The calculated values were used in numerical algorithm to recover bone porosity. At the microlevel, bone was modeled as a heterogeneous medium composed of trabeculae tissue and bone marrow treated as transversely isotropic elastic and isotropic viscoelastic materials, respectively. Recovered porosity values are in excellent agreement with true porosity found from the corresponding micro-CT images.

Analysis of crystals leading to joint arthropathies by Raman spectroscopy: comparison with compensated polarized imaging

The current study assessed the feasibility of the application of Raman spectroscopy toward the diagnosis of gout and pseudogout. First, the lowest concentrations of monosodium urate monohydrate (MSUM) and calcium pyrophosphate dihydrate (CPPD) crystals detectable by Raman spectroscopy were investigated by mixing known amounts of synthetic crystals with synovial fluid in the concentration range of 1 to 100 microg/mL. Second, a digestion protocol was developed for clinical samples to improve crystal extraction. The ensuing centrifugation of the digest congregated crystals at a well-defined point and allowed for point-and-shoot Raman analysis without having to conduct an extensive search for individual crystals. Finally, synovial fluid samples obtained from patients (n = 35) were cross-analyzed by polarized light microscopy (PLM) and the Raman method to compare and contrast the diagnoses of the two methods. It was found that Raman spectroscopy can detect MSUM and CPPD crystals with good sensitivity and specificity at concentrations as low as 5 microg/mL and 2.5 microg/mL, respectively, using the current method. This detection limit of Raman analysis is lower than that reported for PLM. Raman and PLM diagnoses of clinical samples agreed in 32 out of 35 samples in the entire sample pool. However, the rate of disagreement between PLM-based and Raman-based diagnoses was noteworthy within the subset of diseased samples (3 out of 10), indicating that PLM has limitations and that the confirmation by a secondary method is essential for a reliable outcome. The proposed protocol of sample preparation and Raman analysis ascribes baseline feasibility to the diagnosis of gout and pseudogout by Raman spectroscopy, thus justifying further studies using a larger clinical sample set for obtaining sensitivity and specificity.

Human cancellous bone from T12-L1 vertebrae has unique microstructural and trabecular shear stress properties

Increase of trabecular stress variability with loss of bone mass has been implicated as a mechanism for increased cancellous bone fragility with age and disease. In the current study, a previous observation that trabecular shear stress estimates vary along the human spine such that the cancellous tissue from the thoracic 12 (T12)-lumbar 1 (L1) junction experiences the highest trabecular stresses for a given load was tested as a formal hypothesis using multiple human spines. Thoracic 4, T5, T7, T9, T10, T12, L1, L2, L4 and L5 vertebrae from 10 human cadaver spines were examined. One specimen in the central anterior region was cored in the supero-inferior (SI) direction and another in the postero-lateral region was cored in the transverse (TR) direction from each vertebra. Micro-CT-based large-scale finite element models were constructed for each specimen and compression in the long axis of the cylindrical specimens was simulated. Cancellous bone modulus and the mean, the standard deviation, variability and amplification of trabecular von Mises stresses were computed. Bone volume fraction, trabecular number, trabecular thickness, trabecular separation, connectivity density and degree of anisotropy were calculated using 3D stereology. The results were analyzed using a mixed model in which spine level was modeled using a quadratic polynomial. The maximum of trabecular shear stress amplification and minimum of bone volume fraction were found in the cancellous tissue from the T12-L1 location when results from the samples of the same vertebra were averaged. When groups were separated, microstructure and trabecular stresses varied with spine level, extrema being at the T12-L1 levels, for the TR specimens only. SI/TR ratio of measured parameters also had quadratic relationships with spine level, the extrema being located at T12-L1 levels for most parameters. For microstructural parameters, these ratios approached to a value of one at the T12-L1 level, suggesting that T12-L1 vertebrae have more uniform cancellous tissue properties than other levels. The mean intercept length in the secondary principal direction of trabecular orientation could account for the variation of all mechanical parameters with spine level. Our results support that cancellous tissue from T12-L1 levels is unique and may explain, in part, the higher incidence of vertebral fractures at these levels.

Age-related changes in porosity and mineralization and in-service damage accumulation

It has been proposed that bone damageability (i.e. bone's susceptibility to formation of damage) increases with the elevation or suppression of bone turnover. Suppression of turnover via bisphosphonates increases local bone mineralization, which theoretically should increase the susceptibility of bone to microcrack formation. Elevation of bone turnover has also been proposed to increase bone microdamage through an increase in bone intracortical porosity and local stresses and strains. The goal of this paper was to investigate the above proposals, i.e., whether or not increases to mineral content and porosity increase bone in-service damageability. To do this, we measured in vivo diffuse damage area (Df.Dm.Ar, %) and microcrack density (Cr.Dn) (cracks/mm(2)) in the same specimen from human cortical bone of the midshaft of the proximal femur obtained from cadavers with an age range of eight decades and examined their relationships with porosity, mineralization and age. Results of this study showed that Cr.Dn and Df.Dm.Ar increased with a decrease in bulk mineralization. This finding does not appear to support the proposal that damage accumulation increases with low bone turnover that results in increases mineralization. It was proposed however that the negative correlation between damage accumulation and mineralization may be attributed to highly mineralized regions of bone existing with under-mineralized regions resulting in an overall decrease in average bone mineralization. It was also found that microdamage accumulates with increasing porosity which does appear to support the proposal that elevated bone turnover that results in increased porosity can accelerate microdamage accumulation. Finally, it was shown that linear microcracks and Df.Dm.Ar accumulate with age differently, but because they correlate with each other, one may be the precursor for the other.

Trabecular shear stress amplification and variability in human vertebral cancellous bone: relationship with age, gender, spine level and trabecular architecture

Trabecular shear stress magnitude and variability have been implicated in damage formation and reduced bone strength associated with bone loss for human vertebral bone. This study addresses the issue of whether these parameters change with age, gender or anatomical location, and if so whether this is independent of bone mass. Additionally, 3D-stereology-based architectural parameters were examined in order to establish the relationship between stress distribution parameters and trabecular architecture. Eighty cancellous bone specimens were cored from the anterior region of thoracic 12 and donor-matched lumbar 1 vertebrae from a randomly selected population of 40 cadavers. The specimens were scanned at 21-microm voxel size using microcomputed tomography (microCT) and reconstructed at 50microm. Bone volume fraction (BV/TV), trabecular number (Tb.N), trabecular thickness (Tb.Th), trabecular separation (Tb.Sp), bone surface-to-volume ratio (BS/BV), degree of anisotropy (MIL1/MIL3), and connectivity density (-#Euler/Vol) were calculated directly from micro-CT images. Large-scale finite element models were constructed and superoinferior compressive loading was simulated. Apparent cancellous modulus (EFEM) was calculated. The average trabecular von Mises stress generated per uniaxial apparent stress (sigma (-)VM / sigmaapp) and coefficient of variation of trabecular von Mises stresses (COV) were calculated as measures of the magnitude and variability of shear stresses in the trabeculae. Mixed-models and regression were used for analysis. sigma(-)VM / sigmaapp and COV were not different between genders and vertebrae. Both sigma(-)VM / sigmaapp and COV increased with age accompanied by a decrease in BV/TV. Strong relationship of sigma(-)VM / sigmaapp with BV/TV was found whereas COV was strongly related to EFEM/(BV/TV). The results from T12 and L1 were not different and highly correlated with each other. The relationship of sigma(-)VM / sigmaapp with COV was observed to be different between males and females. This difference could not be explained by architectural parameters considered in this study. Our results support the relevance of trabecular shear stress amplification and variability in age-related vertebral bone fragility. The relationships found are expected to help understand the micro-mechanisms by which cancellous bone mass and mechanical properties are modulated through a collection of local stress parameters.

The effect of yield damage on the viscoelastic properties of cortical bone tissue as measured by dynamic mechanical analysis

We have previously shown, using Dynamic Mechanical Analysis (DMA), that the presence of a defect in cortical bone tissue affects the apparent viscoelastic properties of that bone. However, mechanically induced damage is more complex than a machined defect making it difficult to predict its effect on bone viscoelasticity. We performed DMA measurements before and after introduction of yield damage into cortical bone beams from sheep radii. The specimens were placed in a DMA machine and baseline measurements of storage modulus (E1) and loss factor (tandelta) were performed using a 3-point bending configuration for a frequency range of 1-10 Hz. Measurements were done in all four bending directions (cranial, caudal, medial, and lateral) in random order. After subjecting the specimens to monotonic yield damage in a servohydraulic testing machine with the load applied to the cranial surface, oscillatory tests were repeated. To supplement results from the current experiment, additional analyses were performed on data from experiments where bone was either cut or fatigue-loaded between viscoelasticity measurements. Introduction of mechanical damage increased tan delta and frequency sensitivity of E1, consistent with the assertion that increased energy dissipation in damaged bone might contribute to its increased resistance to fatigue and fracture.

The Effect of Regional Variations of the Trabecular Bone Properties on the Compressive Strength of Human Vertebral Bodies

Cancellous centrum is a major component of the vertebral body and significantly contributes to its structural strength and fracture risk. We hypothesized that the variability of cancellous bone properties in the centrum is associated with vertebral strength. Microcomputed tomography (micro-CT)-based gray level density (GLD), bone volume fraction (BV/TV), and finite element modulus (E) were examined for different regions of the trabecular centrum and correlated with vertebral body strength determined experimentally. Two sets of images in the cancellous centrum were digitally prepared from micro-CT images of eight human vertebral bodies (T10-L5). One set included a cubic volume (1 per vertebral centrum, n = 8) in which the largest amount of cancellous material from the centrum was included but all the shell materials were excluded. The other set included cylindrical volumes (6 per vertebral centrum, n = 48) from the anterior (4 regions: front, center, left, and right of the midline of vertebra) and the posterior (2 regions: left and right) regions of the centrum. Significant positive correlations of vertebral strength with GLD (r (2) = 0.57, p = 0.03) and E (r (2) = 0.63, p = 0.02) of the whole centrum and with GLD (r (2) = 0.65, p = 0.02), BV/TV (r (2) = 0.72, p = 0.01) and E (r (2) = 0.85, p = 0.001) of the central region of the vertebral centrum were found. Vertebral strength decreased with increasing coefficient of variation of GLD, BV/TV, and E calculated from subregions of the vertebral centrum. The values of GLD, BV/TV, and E in centrum were significantly smaller for the anterior region than for the posterior region. Overall, these findings supported the significant role of regional variability of centrum properties in determining the whole vertebral strength.

Biomechanical analysis of differing pedicle screw insertion angles

BACKGROUND

Pedicle screw fixation to stabilize lumbar spinal fusion has become the gold standard for posterior stabilization. A significant percentage of surgical candidates are classified as obese or morbidly obese. For these patients, the depth of the incisions and soft tissue makes it extremely difficult to insert pedicle screws along the pedicle axis. As such, the pedicle screws can only be inserted in a much more sagittal axis. However, biomechanical stability of the angled screw insertion has been controversial. We hypothesized that the straight or parallel screw was a more stable construct compared to the angled or axially inserted screw when subjected to caudal cyclic loading.

METHODS

We obtained 12 fresh frozen lumbar vertebrae from L3 to L5 from five cadavers. Schantz screws (6.0 mm) were inserted into each pedicle, one angled and along the axis of the pedicle and the other parallel to the spinous process. Fluoroscopic imaging was used to guide insertion. Each screw was then subjected to caudal cyclic loads of 50 N for 2000 cycles at 2 Hz. Analysis of initial damage, initial rate of damage, and total damage during cyclic loading was undertaken.

FINDINGS

Average total fatigue damage for straight screws measured 0.398+/-0.38 mm, and 0.689+/-0.96 mm for angled screws. Statistical analysis for total fatigue damage ratio of angled to straight screws revealed that a significant stability was achieved in straight-screw construct (P<0.03).

INTERPRETATION

This study showed that straight screw insertion results in a more stable pedicle-screw construct. The angled screw insertion technique resulted in more scattered values of damage indicating that the outcome from the angled screw fixation is less predictable. This validates the use of this technique to implant pedicle screws across the axis of the pedicle (parallel to the mid sagittal line) rather than along the axis, and has broad implications in instrumented posterior lumbar spinal surgery.

Evaluation of filler materials used for uniform load distribution at boundaries during structural biomechanical testing of whole vertebrae

This study was designed to compare the compressive mechanical properties of filler materials, Wood's metal, dental stone, and polymethylmethacrylate (PMMA), which are widely used for performing structural testing of whole vertebrae. The effect of strain rate and specimen size on the mechanical properties of the filler materials was examined using standardized specimens and mechanical testing. Because Wood's metal can be reused after remelting, the effect of remelting on the mechanical properties was tested by comparing them before and after remelting. Finite element (FE) models were built to simulate the effect of filler material size and properties on the stiffness of vertebral body construct in compression. Modulus, yield strain, and yield strength were not different between batches (melt-remelt) of Wood's metal. Strain rate had no effect on the modulus of Wood's metal, however, Young's modulus decreased with increasing strain rate in dental stone whereas increased in PMMA. Both Wood's metal and dental stone were significantly stiffer than PMMA (12.7 +/- 1.8 GPa, 10.4 +/- 3.4 GPa, and 2.9 +/- 0.4 GPa, respectively). PMMA had greater yield strength than Wood's metal (62.9 +/- 8.7 MPa and 26.2 +/- 2.6 MPa). All materials exhibited size-dependent modulus values. The FE results indicated that filler materials, if not accounted for, could cause more than 9% variation in vertebral body stiffness. We conclude that Wood's metal is a superior moldable bonding material for biomechanical testing of whole bones, especially whole vertebrae, compared to the other candidate materials.

Effect of fixation and embedding on Raman spectroscopic analysis of bone tissue

Raman spectroscopy provides valuable information on the physicochemical properties of hard tissues. While the technique can analyze tissues in their native state, analysis of fixed, embedded, and sectioned specimens may be necessary on certain occasions. The information on the effects of fixatives and embedding media on Raman spectral properties is limited. We examined the effect of ethanol and glycerol as fixatives and a variety of embedding media (Araldite, Eponate, Technovit, glycol methacrylate, polymethyl methacrylate, and LR white) on Raman spectral properties (mineralization, crystallinity, and carbonation) measured from the cortical bone of mouse humeri. Humeri were fixed in ethanol or glycerol, followed by embedding in one of the media. Nonfixed, freeze-dried, and fixed but not embedded sections were also examined. Periosteal, endosteal, and midosteal regions of the intracortical envelope were analyzed. Raman spectra of fixative solutions and embedding media were also recorded separately in order to examine the specifics of overlap between spectra. We found significant effects of fixation, embedding, and anatomical location on Raman spectral properties. The interference of ethanol with tissue seemed to be relatively less pronounced than that of glycerol. However, there was no single combination of fixation and embedding that left Raman spectral parameters unaltered. We conclude that careful selection of a fixation and embedding combination should be made based on the parameter of interest and the type of tissue. It may be necessary to process multiple samples from the tissue, each using a combination appropriate for the Raman parameter in question.

Do sacrificial bonds affect the viscoelastic and fracture properties of bone?

Sacrificial bonds have been suggested as a toughening mechanism for bone tissue. Ionic bridges formed by divalent calcium ions between collagen molecules have been proposed as candidates for sacrificial bonds. If this mechanism is active at the macroscopic level, we should observe changes in mechanical properties of bone when calcium ions are maintained or removed from the tissue. To test this hypothesis, we measured viscoelastic and monotonic mechanical properties of cortical bone subjected to differing ionic environments. Storage modulus of bone could be changed up to 3.8% by the presence or absence of Na+ or Ca++ in the environment in a reversible fashion when bones were monitored continuously during treatments. A long-term one-time treatment increased the viscoelastic properties of bone soaked in Na+ solutions whereas the viscoelastic properties of bones soaked in Ca++ solutions were maintained. However, the strength and toughness of bone specimens soaked and fractured in treatment solutions were not improved. The presence of Ca++ affected the mechanical behavior of mineralized bone tissue at the macro scale. These effects were reversible, consistent with the original proposal. However, these effects may not necessarily indicate an increase in strength or toughness of the tissue at the macro scale.

Comparison of the Linear Finite Element Prediction of Deformation and Strain of Human Cancellous Bone to 3D Digital Volume Correlation Measurements

The mechanical properties of cancellous bone and the biological response of the tissue to mechanical loading are related to deformation and strain in the trabeculae during function. Due to the small size of trabeculae, their motion is difficult to measure. To avoid the need to measure trabecular motions during loading the finite element method has been used to estimate trabecular level mechanical deformation. This analytical approach has been empirically successful in that the analytical models are solvable and their results correlate with the macroscopically measured stiffness and strength of bones. The present work is a direct comparison of finite element predictions to measurements of the deformation and strain at near trabecular level. Using the method of digital volume correlation, we measured the deformation and calculated the strain at a resolution approaching the trabecular level for cancellous bone specimens loaded in uniaxial compression. Smoothed results from linearly elastic finite element models of the same mechanical tests were correlated to the empirical three-dimensional (3D) deformation in the direction of loading with a coefficient of determination as high as 97% and a slope of the prediction near one. However, real deformations in the directions perpendicular to the loading direction were not as well predicted by the analytical models. Our results show, that the finite element modeling of the internal deformation and strain in cancellous bone can be accurate in one direction but that this does not ensure accuracy for all deformations and strains.

Matrix concentration of insulin-like growth factor I (IGF-I) is negatively associated with biomechanical properties of human tibial cancellous bone within individual subjects

Insulin-like growth factor-I (IGF-I), abundant in bone matrix, is believed to play an important role during bone development and remodeling. To our knowledge, however, few studies have addressed the relationship between the concentration of IGF-I in bone matrix and the biomechanical properties of bone tissue. In this study, forty-five cylindrical specimens of cancellous bone were harvested from six human tibiae and scanned using micro-computed tomography (microCT). The bone volume fraction (BV/TV) was calculated from three-dimensional (3D) microCT images. Mechanical tests were then performed on a servohydraulic testing system to determine the strength and stiffness of cancellous bone. Following mechanical testing, the concentration of IGF-I in bone matrix was measured by using an enzyme-linked immunoabsorbent assay (ELISA). Within each subject, the concentration of IGF-I in bone matrix had significant (P<0.01) negative correlations with the bone volume fraction, strength, and stiffness of cancellous bone. In particular, the anterior quadrant of the proximal tibia was significantly (P<0.02) greater in IGF-I matrix concentration and marginally significantly lower in strength (P=0.053) and stiffness (P=0.059) than the posterior quadrant. The negative correlations between the cancellous bone matrix concentration of IGF-I and cancellous bone biomechanical properties within subjects found in this study may help us understand the variation of the biomechanical properties of cancellous bone in proximal human tibiae.

Effect of Microcomputed Tomography Voxel Size on the Finite Element Model Accuracy for Human Cancellous Bone

The level of structural detail that can be acquired and incorporated in a finite element (FE) analysis might greatly influence the results of microcomputed tomography (μCT)-based FE simulations, especially when relatively large bones, such as whole vertebrae, are of concern. We evaluated the effect of scanning and reconstruction voxel size on the μCT-based FE analyses of human cancellous tissue samples for fixed- and free-end boundary conditions using different combinations of scan/reconstruction voxel size. We found that the bone volume fraction (BV/TV) did not differ considerably between images scanned at 21 and 50 μm and reconstructed at 21, 50, or 110 μm (−0.5% to 7.8% change from the 21/21 μm case). For the images scanned and reconstructed at 110 μm, however, there was a large increase in BV/TV compared to the 21/21 μm case (58.7%). Fixed-end boundary conditions resulted in 1.8% [coefficient of variation (COV)] to 14.6% (E) difference from the free-end case. Dependence of model output parameters on scanning and reconstruction voxel size was similar between free- and fixed-end simulations. Up to 26%, 30%, 17.8%, and 32.3% difference in modulus (E), and average (VMExp), standard deviation (VMSD) and coefficient of variation (COV) of von Mises stresses, respectively, was observed between the 21/21 μm case and other scan/reconstruction combinations within the same (free or fixed) simulation group. Observed differences were largely attributable to scanning resolution, although reconstruction resolution also contributed significantly at the largest voxel sizes. All 21/21 μm results (taken as the gold standard) could be predicted from the 21/50 radj2=0.91-0.99;p<0.001, 21/110 radj2=0.58-0.99;p<0.02 and 50/50 results radj2=0.61-0.97;p<0.02. While BV/TV, VMSD, and VMExp/σz from the 21/21 could be predicted by those from the 50/110 radj2=0.63-0.93;p<0.02 and 110/110 radj2=0.41-0.77;p<0.05 simulations as well, prediction of E, VMExp, and COV became marginally significant 0.04<p<0.13 at 50/110 and nonsignificant at 110/110 0.21<p<0.70. In conclusion, calculation of cancellous bone modulus, mean trabecular stress, and other parameters are subject to large errors at 110/110 μm voxel size. However, enough microstructural details for studying bone volume fraction, trabecular shear stress scatter, and trabecular shear stress amplification VMExp/σz can be resolved using a 21/110 μm, 50/110 μm, and 110/110 μm voxels for both free- and fixed-end constraints.

The effect of microcomputed tomography scanning and reconstruction voxel size on the accuracy of stereological measurements in human cancellous bone

Stereological parameters have been used as an approximation for the architecture of trabecular bone. Structural indices such as bone volume fraction (BV/TV), trabecular number (Tb.N), trabecular thickness (Tb.Th), trabecular separation (Tb.Sp), bone surface-to-volume ratio (BS/BV), degree of anisotropy (MIL1/MIL3), and connectivity density (-Euler/Vol) have been widely studied to investigate pathological conditions in bone. Due to its high resolution and nondestructiveness, microcomputed tomography (micro-CT) has been utilized to take precise three-dimensional (3D) images of trabecular microstructures. However, spatial limitations for applying micro-CT-based analyses to large specimens, such as whole vertebral bodies, require using larger scanning and reconstruction voxel sizes. In this study, combinations of three different scanning and reconstruction voxel size were used to represent best possible voxel size (21 microm; best in our scanner for the specimen size used) relative to other voxel sizes used in this study, commonly used intermediate voxel sizes (50 microm), and those applicable to scans of whole human vertebral bodies (110 microm) in order to examine the effect of scanning and reconstruction voxel size on stereological measures for human cancellous bone. The error in stereological parameters calculated using combinations of large voxel sizes compared to the gold standard (best possible case) ranged from 0.1% to 102%. The signed magnitude of the error in other cases relative to the gold standard was a function of either scanning or reconstruction voxel size or both (r2=0.55-0.95). For most of the structural indices, the results from analysis of images with larger voxel sizes were correlated with those from the gold standard (r2=0.55-0.99) except for Tb.N at 110/110 microm, MIL1/MIL3 at larger than 110 microm reconstruction voxel size, and -Euler/Vol at any combination of voxel sizes. Overall, it was observed that resampling a high resolution image at lower resolutions (corresponding to increasing reconstruction voxel size in this study) had different effects on the calculated parameters than scanning at the same low resolution (corresponding to increasing scanning voxel size in this study). Our results show that investigations of image resolution should include actual scans at the resolution of interest rather than simply coarsening of high-resolution images as is customarily done.

Fracture mechanics of cortical bone tissue: a hierarchical perspective

The performance of bone tissue in the presence of flaws is a highly remarkable one. Bone tissue is the outcome of an adaptive evolutionary process; thus, insight into the mechanisms by which it fails would provide valuable information not only for development of mechanically superior biomimetic materials but also for development of treatment modalities to prevent debilitating bone fractures. Clinically, fractures of skeletal organs occur as a result of aging, disease, overuse, and trauma. Fracture mechanics, a sub-discipline of solid mechanics that investigates the performance of cracked materials, has been employed extensively in characterizing the mechanisms by which bone tissue fractures. At present the fracture mechanisms at the macroscale are better characterized than at the microscale. On the other hand, a mechanistic understanding of damage evolution at the submicroscopic scale is largely limited to postulations with little experimental insight. The challenge of skeletal fragility will be dealt with more efficiently with deeper understanding of the fracture process at each hierarchical size scale. The most recent review on this subject matter was a decade ago, and there have been numerous developments in the fracture mechanics of bone since then. This review recaps the existing literature with an emphasis on the hierarchical nature of the fracture process in bone, entailing the supramolecular, microscopic, and macroscopic scales.

Effects of end boundary conditions and specimen geometry on the viscoelastic properties of cancellous bone measured by dynamic mechanical analysis

The viscoelastic properties of cancellous bone can be measured nondestructively in compression testing using a dynamic mechanical analyzer. In this study, we examined the effects of end boundary conditions and specimen geometry on the viscoelastic properties of cancellous bone measured by dynamic mechanical analysis. During dynamic compression testing, the cancellous bone specimens may be mechanically fixed (e.g., glued) to the loading platens or they may be free to expand across the platen surface. When specimens of cancellous bone were tested between platens with gluing, the dependence of loss tangent on frequency was not consistent with previously observed strain-rate-dependent mechanical behavior of cancellous bone. When long specimens of cancellous bone (length = 10 mm, diameter = 8 mm) were tested without gluing, the relationship between loss tangent and frequency depended on the level of load applied. For short specimens (length = 5 mm, diameter = 8 mm) tested without gluing, however, the frequency dependence of loss tangent agreed with existing data reported for the strain-rate-dependent behavior of cancellous bone and also with the frequency dependence of cortical bone viscoelasticity. Therefore, we recommend that short cancellous bone cylinders with a length of 5 mm and a diameter of 8 mm should be used without gluing in the dynamic mechanical analysis of cancellous bone. This is consistent with the American Society for Testing and Materials testing recommendations for plastics, but different from current practice for unimodal mechanical testing of cancellous bone.

The dependence between the strength and stiffness of cancellous and cortical bone tissue for tension and compression: extension of a unifying principle

A strong positive correlation between the apparent ultimate strength and stiffness of bone tissue that can be expressed by a unified relationship has been observed for cortical bone in tension and low-density cancellous bone in compression. For practical purposes, the existence of a relationship between strength and stiffness is significant in that bone stiffness can be measured in vivo using non-invasive methods. It is generally accepted that bone strength is greater in compression than in tension whereas there is no substantial evidence that bone stiffness in compression is different from that in tension. This might suggest that compressive strength would relate to the stiffness, if at all, in a way that is different from tensile strength. In order to examine similarities and differences in the way strength is associated with stiffness between modes of loading and tissue type, we tested equine cortical bone and bovine cancellous bone in compression and examined these data together with previously reported data from compression testing of human cancellous bone as well as tensile testing of cortical bone from various sources. We have found for cortical bone that (i) the sensitivity of strength to stiffness is the same for tension and compression (p>0.75, ANCOVA), and (ii) the difference between the magnitudes of compressive and tensile strength for cortical bone is the result of an additive, rather than a multiplicative factor (52.1 MPa after adjusting to 1 microstrain/s, p<0.0001, ANOVA). High-density bovine tibial cancellous bone, on the other hand, has a steeper slope for its compressive strength-stiffness relationship than that for cortical bone and human cancellous bone, resulting in a transitional relationship between compressive strength and stiffness for a range of bone types and densities. Based on the current results and previous work, it is suggested that the offset strength in the compressive strength-stiffness relationship may be a direct manifestation of the difference between the compressive and tensile strengths of the bone material that constitutes the building blocks of the bone structure. Deviation of high-density cancellous bone compressive behavior from the other bone types and densities is attributed to stress distribution differences between the bone types.

Apparent viscoelastic anisotropy as measured from nondestructive oscillatory tests can reflect the presence of a flaw in cortical bone

There is evidence that damage, viscoelastic stiffness properties, and postyield mechanical properties are related in bone tissue. Our objective was to test whether presence of a flaw would have an influence on the apparent viscoelastic properties of bone. Examining the effect of flaw orientation on apparent viscoelastic properties and utilization of dynamic mechanical analysis (DMA) as a nondestructive means for detection of damage were our secondary objectives. Cortical bone beams (2 x 2 x 19 mm) machined from the cranial cortex of the radii of six Warhill sheep were used. The specimens were placed in a DMA machine and baseline measurements of storage modulus (E1) and loss factor (tan delta), once for loads in the craniocaudal and once in the mediolateral directions, were performed using a three-point bending configuration for a frequency range of 1-10 Hz. Craniocaudal/mediolateral measurement ratio was calculated as a measure of anisotropy for tan delta and E1. After cutting a thin through-thickness macroscopic notch on the caudal surface at the center of each beam, oscillatory tests were repeated. Two-way repeated measures analysis of variance followed by Tukey's test was used with group (craniocaudal, mediolateral, notched craniocaudal, and notched mediolateral measurements) and frequency as factors. Regression analysis and analysis of covariance were used for examining the relationship between viscoelastic parameters and frequency. Tan delta and E1 were not different between craniocaudal and mediolateral measurements before the flaw was introduced (p > 0.8 and p = 1, respectively). In the presence of the flaw, tan delta was significantly increased (p < 0.003) whereas E1 was significantly reduced (p < 0.001) for craniocaudal measurements. Tan delta and E1 were nearly isotropic in the tested directions before the introduction of a flaw into the bone tissue. Introduction of a flaw resulted in increased tan delta and E1 anisotropy. Presence of a notch resulted in a significant increase in tan delta anisotropy with increasing frequency. In conclusion, we have demonstrated that cortical bone tissue exhibits a different apparent viscoelastic behavior in the presence of a flaw and depending on the flaw's orientation. Our finding that the presence of a notch and its orientation can be detected by nondestructive DMA suggests that in vivo techniques may be developed for detection of cortical bone damage.

A rate-dependent microcrack-bridging model that can explain the strain rate dependency of cortical bone apparent yield strength

Although there are empirical correlations between strain rate, cortical and cancellous bone apparent stiffness, apparent yield strength, apparent ultimate strength and cortical bone fracture toughness, a mechanistic description for these phenomena is lacking. Microcracking is a major mechanism in cortical and cancellous bone failure, however, microdamage content alone cannot explain the strain rate dependence of bone strength without considering time-dependent behavior of the crack. Using a rate-dependent model of a fiber-bridged microcrack and data from the literature, we demonstrate that the experimental apparent yield strength of bone can be predicted directly from measurements of apparent moduli of elasticity of bone constituents and failure strain of the collagenous matrix. Yield strength predictions for estrogen depleted bone were made using the model and data from ovariectomized sheep. It was predicted that the yield strength of estrogen-deficient bone is comparable to that of normal bone within strain rates associated with physiological activities. For high strain rates, however, the strength of estrogen-depleted bone was predicted to be much weaker than normals suggesting a higher fracture risk due to impact from falls, for individuals with estrogen-depleted bones such as in post-menopausal osteoporosis.