Digital Tomosynthesis and Fractal Analysis Predict Prevalent Vertebral Fractures in Patients With Multiple Myeloma: A Preliminary In Vivo Study

Textural and geometric measures derived from digital tomosynthesis discriminate women with and without vertebral fracture

Vertebral fractures are a common and debilitating consequence of osteoporosis. Bone mineral density (BMD), measured by dual energy x-ray absorptiometry (DXA), is the clinical standard for assessing overall bone quantity but falls short in accurately predicting vertebral fracture. Fracture risk prediction may be improved by incorporating metrics of microstructural organization from an appropriate imaging modality. Digital tomosynthesis (DTS)-derived textural and microstructural parameters have been previously correlated to vertebral bone strength in vitro, but the in vivo utility has not been explored. Therefore, the current study sought to establish the extent to which DTS-derived measurements of vertebral microstructure and size discriminate patients with and without vertebral fracture. In a cohort of 93 postmenopausal women with or without history of vertebral fracture, DTS-derived microstructural parameters and vertebral width were calculated for T12 and L1 vertebrae, as well as lumbar spine BMD and trabecular bone score (TBS) from DXA images. Fracture patients had lower BMD and TBS, while DTS-derived degree of anisotropy and vertebral width were higher, compared to nonfracture (p < 0.02 to p < 0.003) patients. The addition of DTS-derived parameters (fractal dimension, lacunarity, degree of anisotropy and vertebral width) improved discriminative capability for models of fracture status (AUC = 0.79) compared to BMD alone (AUC = 0.67). For twelve additional participants who were imaged twice, in vivo repeatability errors for DTS parameters were low (0.2 % - 7.3 %). The current results support the complementary use of DTS imaging for assessing bone quality and improving the accuracy of fracture risk assessment beyond that achievable by DXA alone.

Novel Use of Fractal Analysis for Quantifying Polymethylmethacrylate Distribution Patterns in Osteoporotic and Malignant Vertebral Compression Fractures Following Vertebroplasty

Purpose: Fractal analysis is a mathematical tool which allows the evaluation of complex microstructural features within materials that cannot be expressed in traditional geometric terms. The purpose of this study is to quantify the differences in polymethylmethacrylate intravertebral cement spatial distribution patterns following vertebroplasty using fractal analysis through the examination of osteoporotic and malignant compression fractures. Methods: Frontal and lateral post-vertebroplasty radiographs were evaluated from 29 patients with osteoporotic and malignant compression fractures who underwent vertebroplasty. The individually treated vertebra were divided into osteoporotic (n = 35) and malignant groups (n = 41). Images underwent segmentation, thresholding, and binarization prior to fractal analysis. Fractal dimension and lacunarity values were derived from the region of interest in treated vertebrae using the "box-counting" and "gliding-box" techniques respectively using ImageJ. The mean values of both parameters were compared between the 2 groups. Results: The mean fractal dimension was significantly higher in the malignant vertebral compression fracture group (1.53 ± 0.08) compared to the osteoporotic group (1.34 ± 0.17; P < .001). Similarly, mean lacunarity values were significantly higher in the malignant fracture group (0.50 ± 0.09) compared to the osteoporotic group (0.37 ± 0.10; P < .001). Conclusions: Fractal dimension and lacunarity values of cement spatial distribution patterns obtained from the post-vertebroplasty radiographs can differentiate between benign osteoporotic and malignant vertebral compression fractures. This novel technique may be useful for evaluating cement spatial distribution patterns in spine augmentation procedures, although further research is warranted in this area.

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.

Evaluation of mandibular bone changes in multiple myeloma patients on dental panoramic radiographs

OBJECTIVE

The aim of this study is to assess the mandibular bone structures of multiple myeloma (MM) patients on panoramic radiographs using fractal analysis (FA) and radio-morphometric indices.

METHODS

The study group consisted of 67 MM patients, and the control group consisted of 67 individuals without any systemic disease. The mandibular cortical index (MCI) classification, antegonial index (AGI), gonial index (GI), mandibular cortical width (MCW), panoramic mandibular index (PMI) and fractal dimensions (FD) were evaluated on panoramic radiographs. FD measurements were made by identifying 4 related areas. Shapiro-Wilk test was used to determine the normal distribution of the data. Chi-square and independent t tests were used to compare the findings between the two groups.

RESULTS

The FD values of the study group in ROI1, ROI2, ROI3, ROI4 regions and GI values were found to be statistically significantly lower than the control group (p < 0.001). There was no statistically significant difference between the two groups in terms of AGI, MCW and PMI values (p > 0.05). While C1 was the most common type in the control group, C2 was the most common type in the study group. C2 and C3 were detected more in the study group than in the control group (p < 0.001).

CONCLUSION

Our study showed a consensus with the studies advocating that fractal analysis and radio-morphometric indices are methods that can be used to determine mandibular bone density. The low bone density of MM patients is a condition that physicians should be aware of for interventional dental procedures.

Application of tomosynthesis for vertebral compression fracture diagnosis and bone healing assessment in fracture liaison services

Early identification of vertebral compression fractures (VCFs) is crucial for successful secondary fracture prevention. Tomosynthesis, a low-dose tomographic imaging technique, may facilitate the evaluation and long-term follow-up of VCFs in patients with osteoporosis. Herein, we compared the performances of plain radiography and tomosynthesis for VCF diagnosis and healing assessment in patients enrolled in fracture liaison services in our hospital. Forty-nine patients with new VCFs at the T10–L5 levels were prospectively recruited between August 2018 and May 2020; all patients underwent thoracolumbar plain radiography and tomosynthesis. We evaluated the accuracy of the VCF diagnosis, image quality, and VCFs healing process. Tomosynthesis identified 90 levels of VCF in 49 patients, while plain radiography revealed only 87.8% (79/90) of them. There were 44.9% (22/49) patients with neglected chronic VCFs as seen on tomosynthesis. Tomosynthesis images had improved VCF diagnostic accuracy up to 12.2% and showed significantly more anatomic details than plain radiography. For diagnosis of VCFs, the performance of plain radiographs was poorer than that of tomosynthesis images (plain radiographs: sensitivity 84%, specificity 93.5%, false positive rate 6.5%, and false negative rate 16%; tomosynthesis: sensitivity 93.2%, specificity 100%, false positive rate 0%, and false negative 6.8%), using magnetic resonance imaging (MRI) as gold standard. The Kappa coefficient between Tomosynthesis and MRI is 0.956 while between radiography and MRI is 0.704. Tomosynthesis showed significantly more anatomic details than plain radiography and all the examiners revealed a clear preference for tomosynthesis. Tomosynthesis scored 3.3 times higher on the fracture healing assessment at the 3-month follow-up than plain radiographs. Tomosynthesis is a promising tool for VCF screening and diagnosis in patients with osteoporosis and for monitoring fracture healing status at a low radiation dose and cost.

Structural complexity of the craniofacial trabecular bone in multiple myeloma assessed by fractal analysis

Purpose

This study aimed to evaluate the structural complexity of craniofacial trabecular bone in multiple myeloma by fractal analysis of panoramic and lateral skull radiography, and to compare the fractal dimension values of healthy patients (HPs), pre-treatment patients (PTPs), and patients during bisphosphonate treatment (DTPs).

Materials and Methods

Pairs of digital panoramic and lateral skull radiographs of 84 PTPs and 72 DTPs were selected. After application of exclusion criteria, 43 panoramic and 84 lateral skull radiographs of PTPs, 56 panoramic and 72 lateral skull radiographs of DTPs, and 99 panoramic radiographs of age- and sex-matched HPs were selected. The fractal dimension values from panoramic radiographs were compared among HPs, PTPs, and DTPs and between anatomical locations within patient groups using analysis of variance with the Tukey test. Fractal dimension values from lateral skull radiographs were compared between PTPs and DTPs using the Student t-test. Pearson correlation coefficients were used to assess the relationship between the mandible from panoramic radiographs and the skull from lateral skull radiographs. Intra-examiner agreement was assessed using intraclass correlation coefficients (α=0.05).

Results

The fractal dimension values were not significantly different among HPs, PTPs, and DTPs on panoramic radiographs or between PTPs and DTPs on lateral skull radiographs (P>0.05). The mandibular body presented the highest fractal dimension values (P≤0.05). The fractal dimension values of the mandible and skull in PTPs and DTPs were not correlated.

Conclusion

Fractal analysis was not sensitive for distinguishing craniofacial trabecular bone complexity in multiple myeloma patients using panoramic and lateral skull radiography.

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.

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.