Next Article in Journal
A Robust, Fully Automatic Detection Method and Calculation Technique of Midline Shift in Intracranial Hemorrhage and Its Clinical Application
Next Article in Special Issue
Assessment of Maxillary Molars Interradicular Septum Morphological Characteristics as Criteria for Ideal Immediate Implant Placement—The Advantages of Cone Beam Computed Tomography Analysis
Previous Article in Journal
Current Advancement in Diagnosing Atrial Fibrillation by Utilizing Wearable Devices and Artificial Intelligence: A Review Study
Previous Article in Special Issue
Myelin-Weighted Imaging Presents Reduced Apparent Myelin Water in Patients with Alzheimer’s Disease
 
 
Correction published on 31 October 2022, see Diagnostics 2022, 12(11), 2635.
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Machine Learning for Opportunistic Screening for Osteoporosis from CT Scans of the Wrist and Forearm

by
Ronnie Sebro
1,2,* and
Cynthia De la Garza-Ramos
1
1
Mayo Clinic Florida, Department of Radiology, Jacksonville, FL 32224, USA
2
Center for Augmented Intelligence, Mayo Clinic Florida, Department of Radiology, Jacksonville, FL 32224, USA
*
Author to whom correspondence should be addressed.
Diagnostics 2022, 12(3), 691; https://doi.org/10.3390/diagnostics12030691
Submission received: 20 February 2022 / Revised: 4 March 2022 / Accepted: 8 March 2022 / Published: 11 March 2022 / Corrected: 31 October 2022
(This article belongs to the Special Issue Advances in Diagnostic Medical Imaging)

Abstract

:
Background: We investigated whether opportunistic screening for osteoporosis can be done from computed tomography (CT) scans of the wrist/forearm using machine learning. Methods: A retrospective study of 196 patients aged 50 years or greater who underwent CT scans of the wrist/forearm and dual-energy X-ray absorptiometry (DEXA) scans within 12 months of each other was performed. Volumetric segmentation of the forearm, carpal, and metacarpal bones was performed to obtain the mean CT attenuation of each bone. The correlations of the CT attenuations of each of the wrist/forearm bones and their correlations to the DEXA measurements were calculated. The study was divided into training/validation (n = 96) and test (n = 100) datasets. The performance of multivariable support vector machines (SVMs) was evaluated in the test dataset and compared to the CT attenuation of the distal third of the radial shaft (radius 33%). Results: There were positive correlations between each of the CT attenuations of the wrist/forearm bones, and with DEXA measurements. A threshold hamate CT attenuation of 170.2 Hounsfield units had a sensitivity of 69.2% and a specificity of 77.1% for identifying patients with osteoporosis. The radial-basis-function (RBF) kernel SVM (AUC = 0.818) was the best for predicting osteoporosis with a higher AUC than other models and better than the radius 33% (AUC = 0.576) (p = 0.020). Conclusions: Opportunistic screening for osteoporosis could be performed using CT scans of the wrist/forearm. Multivariable machine learning techniques, such as SVM with RBF kernels, that use data from multiple bones were more accurate than using the CT attenuation of a single bone.

1. Introduction

Bone mineral density (BMD) decreases with age, with the decrease being more evident and rapid in post-menopausal females [1,2,3]. Decreased BMD increases the risk of frailty fractures including fractures of the spine, forearm, and hips [4,5,6,7,8]. Fractures of the hips are associated with increased 1-year mortality following the fracture of approximately 15–36%, therefore identifying risk factors for frailty fractures and hip fractures is of increased clinical importance [9,10]. Dual-energy X-ray absorptiometry (DEXA) is the gold standard screening test for the evaluation of BMD [11]. DEXA evaluates BMD in the L1–L4 lumbar spine, total hip, and femoral neck and compares these values to those of normal young adults (20–29 years of age) from the National Health and Nutrition Examination Survey (NHANES) III cohort to create BMD T-scores [12,13]. The World Health Organization (WHO) and the International Society for Clinical Densitometry (ISCD) have guidelines for the classification of patients aged 50 years or greater with low BMD as osteoporotic or osteopenic based on DEXA BMD T-scores. Patients aged 50 years or greater with lowest BMD T-scores ≥ −1 have normal BMD, patients aged 50 years or greater with −2.5 ≤ lowest BMD T-scores < −1 have osteopenia, and patients aged 50 years or greater with lowest BMD T-scores ≤ −2.5 have osteoporosis [14]. The University of Sheffield Fracture Risk Assessment Tool (FRAX) utilizes patient age, sex, femoral neck BMD, and other variables to predict the risk of future femoral fracture [15]. Trabecular bone scores (TBS) are based on textural analysis of bone microarchitecture, and TBS in conjunction with BMD scores were shown to improve the overall predictive properties of DEXA measurements for future fractures [16,17].
Opportunistic screening for osteoporosis was performed using data from computed tomography (CT) scans of the chest/thorax [18], abdomen and pelvis [19,20], and cardiac CT [21]. A strong correlation was noted between the CT attenuation of trabecular bone in the thoracolumbar spine and BMD measurements [18,19,20,21]. Because forearm frailty fractures are associated with low BMD, we hypothesized that the CT attenuation of the wrist and forearm bones would correlate with BMD measurements. A recent study showed that the cortical thickness of the second metacarpal is also correlated with DEXA BMD measurements [22]. Therefore, we hypothesized that the CT attenuation in the metacarpals may also be correlated with BMD DEXA measurements. We hypothesized that routinely obtained CT scans of the wrist/forearm bones in patients aged 50 years or greater could be used to provide insight into whether a patient has osteoporosis or osteopenia and whether the patient needs to be subsequently screened using DEXA.
The aims of this study were (1) to evaluate the correlations between the CT attenuation of the wrist/forearm bones (distal third of the radius; ultradistal radius (radius UD) (distal radial metaphysis and epiphysis), distal third of the radius shaft (radius 33%), distal third of the ulna, ultradistal ulna (ulna UD) (distal ulnar metaphysis and diaphysis), distal third of the ulnar shaft (ulna 33%), scaphoid, lunate, triquetrum, pisiform, trapezium, trapezoid, capitate, hamate, and proximal thirds of the first through fifth metacarpals) in men and women; (2) to evaluate the correlations between CT attenuation of the wrist/forearm bones and L1–L4 BMD, L1–L4 BMD T-score, L1–L4 TBS, total hip BMD, total hip BMD T-score, femoral neck BMD, and femoral neck BMD T-score; (3) to evaluate the predictive performance of each wrist/forearm bone regarding predicting (i) osteoporosis, (ii) osteopenia/osteoporosis, (iii) femoral neck BMD ≤ −2.5, and (iv) femoral neck BMD < −1; and (4) to use machine learning to identify the best combination of clinical/demographic variables and the CT attenuation of wrist/forearm bones to predict (i) osteoporosis, (ii) osteopenia/osteoporosis, (iii) femoral neck BMD ≤ −2.5, and (iv) femoral neck BMD < −1.

2. Materials and Methods

The study was compliant with the Health Insurance Portability and Accountability Act of 1996 (HIPAA), and the study protocol was reviewed and approved by the local Institutional Review Board (IRB) at a tertiary care academic medical center. The need for signed informed consent from each patient was waived by the IRB. The work described was carried out in accordance with The Code of Ethics of the World Medical Association (Declaration of Helsinki) for experiments involving humans.

2.1. CT Scanner Protocol and CT Attenuation Measurements

CT scans were performed on Siemens Somatom Definition Edge and Siemens Flash (Siemens Healthineers, Erlangen, Germany) scanners. CT scans were performed without intravenous contrast at 120 kVp, 200–250 mA, field of view (FOV) 118–120 mm, tilt of 0° in the axial plane at 0.4–0.5 mm, and reconstructed in the coronal and sagittal planes at 1 mm.

2.2. DEXA Scanners

DEXA scans were performed using General Electric (GE) (Waukesha, WI, USA) Luna iDXA Scanners.
The study cohort comprised patients aged 50 years or older, who were evaluated and/or treated between 1 January 2015 and 30 September 2021. Patients were included if they had CT scans of an upper extremity, including the wrist and forearm, and DEXA scans (of the lumbar spine and hips) within 12 months of each other. Patients were excluded if they had prior surgery of the wrist/forearm, hips, or spine with implantation of hardware that would affect CT attenuation measurements.

2.3. Segmentations

The trabecular component of each bone was carefully, manually segmented using 3D-Slicer [23]. Care was taken to avoid bone lesions (hemangiomas and bone islands), osteophytes, and the bony cortex. No specific attenuation threshold was used to segment the trabecular bone. The subcortical regions were segmented manually and then automatically interpolated between CT slices. The entire volume of each bone was segmented, and the mean CT attenuation of each bone was recorded (Figure 1).
We calculated the CT attenuation of the distal third of the radius, radius UD, radius 33%, the distal third of the ulna, ulna UD, ulna 33%, scaphoid, lunate, triquetrum, pisiform, trapezium, trapezoid, capitate, hamate, and proximal thirds of the first through fifth metacarpals. Segmentations were performed by a trained research assistant (3 weeks of training doing all 3D Slicer tutorials; 4 months of experience with 3D Slicer doing similar projects before doing the current project) and were reviewed by a fellowship-trained musculoskeletal radiologist with 10 years of experience.

2.4. Statistical Analysis

A sample size of 100 patients was determined to have 80% power with a type I error rate of 5% to detect a difference between two receiver operator characteristic (ROC) curves if the first ROC curve has an area under the receiver operator characteristic curve (AUC) of 0.80 and the second ROC curve has an AUC of 0.60, and there are 74% controls and 26% cases. Therefore, the cohort was randomly divided into two datasets—a training/validation dataset (96 (49%) patients) and a testing dataset (100 (51%) patients) for multivariable machine learning.
Summary statistics for all clinical and demographic variables were first calculated. T-tests with unequal variances and Fisher’s exact tests were used to compare quantitative and qualitative variables, respectively, between the training/validation and testing datasets.
Pearson’s correlations were used to examine the correlations between the CT attenuation of the wrist/forearm bones and to evaluate the correlations between CT attenuation of the wrist/forearm bones (distal third of the radius and ulna, proximal third of the first through fifth metacarpals, scaphoid, lunate, triquetrum, pisiform, trapezium, trapezoid, capitate, and hamate) stratified by gender in the entire dataset. Hierarchical cluster analysis was used to cluster the correlations between the CT attenuations of the wrist/forearm bones.
Pearson’s correlation coefficients were also used to evaluate the correlations between the CT attenuations of each of the wrist/forearm bones and DEXA measurements (L1–L4 BMD, L1–L4 BMD T-score, L1–L4 TBS, total hip BMD, total hip BMD T-score, femoral neck BMD, and femoral neck BMD T-score) in the entire dataset.
ROC curves were used to identify the optimal CT attenuation threshold for each bone to evaluate the predictive performance of each bone regarding predicting osteoporosis and osteopenia/osteoporosis in the training/validation dataset. We then evaluated these optimized thresholds in the test dataset.
Since we evaluated the predictive properties of several different bones in the wrist/forearm to predict osteoporosis and osteopenia/osteoporosis, we used several machine learning methods to select the best combination of bones and clinical factors (age, gender, height, weight, BMI) that could be used to categorize a patient as (i) osteoporotic, (ii) osteopenic/osteoporotic using the WHO guidelines, (iii) femoral neck BMD T-score ≤ −2.5, and (iv) femoral neck BMD T-score < −1. We evaluated the femoral neck BMD T-scores because these scores are thought to be less affected by degenerative changes in the lumbar spine and hip.
A support vector machine (SVM) is a supervised learning model that is often used for pattern recognition, classification, and regression analysis [24]. C-classification with three different kernels (linear, radial basis function (RBF), and sigmoid) [24] tuned with 10-fold cross-validation was utilized.
Linear   kernel ,   K x , y = x . y
Radial   basis   function ,   K x , y = e x y 2 / 2 σ 2
Sigmoid   kernel ,   K x , y = tan h υ x . y + c
The tuning for each SVM was performed over epsilon ranges from 0 to 1 with 0.1 increments, with cost values ranging from 1 to 6 with increments of 1 using 10-fold cross-validation.
A random forest classifier is an ensemble learning method for classification based on constructing a multitude of decision trees during training [25]. The random forest model was fit and tuned with the number of variables tried at each step starting at 6, a step factor of 1.5, and 10,000 trees used during the tuning step.
Each of the four machine learning models (linear kernel SVM, RBF kernel SVM, sigmoid kernel SVM, random forest classifier) were used to each categorize patients as (i) osteoporotic using WHO guidelines, (ii) osteopenic/osteoporotic using WHO guidelines, (iii) femoral neck BMD T-score ≤ −2.5, and (iv) femoral neck BMD T-score < −1. The optimal models tuned in the training/validation dataset were retained for each analysis. We used these models to predict (i) osteoporotic using WHO guidelines, (ii) osteopenic/osteoporotic using WHO guidelines, (iii) femoral neck BMD T-score ≤ −2.5, and (iv) femoral neck BMD T-score < −1.
All test statistics were two-sided and p-values < 0.05 were considered statistically significant. We compared the ROC curves using DeLong’s test [26] and compared the machine learning models to the CT attenuation of the radius 33% since this region is used to determine osteoporosis and osteopenia on DEXA studies of the forearm/wrist. Statistics were performed using the pROC, e1071, and RandomForest packages in Rv4.1.2 statistical software.

3. Results

The entire dataset comprised 196 individuals with a median age of 65.0 years (range 50.0–88.0 years), 173 (88.3%) were females, 54 (27.6%) were osteoporotic, 116 (59.2%) were osteopenic, and the other 26 (13.3%) were normal. There were no statistically significant differences between the variables in the training/validation dataset and the test dataset (Table 1).
Figure 2 compares the mean CT attenuation of each bone between diagnoses (osteoporosis, osteopenia, and normal).
There were positive correlations between the CT attenuation of the wrist/forearm bones. In women, the strongest CT attenuation correlations were between the CT attenuations of the ulna and ulna UD (r = 0.74, p < 0.001), ulna and ulna 33% (r = 0.83, p < 0.001), radius and radius UD (r = 0.74, p < 0.001), capitate and trapezoid (r = 0.74, p < 0.001), hamate and capitate (r = 0.81, p < 0.001), and the hamate and trapezoid (r = 0.81, p < 0.001). In men, the strongest CT attenuation correlations were between the CT attenuations of radius and the radius UD (r = 0.82, p < 0.001), radius and radius 33% (r = 0.79, p < 0.001), ulna and ulna UD (r = 0.84, p < 0.001), ulna and ulna 33% (r = 0.84, p < 0.001), scaphoid and lunate (r = 0.84, p < 0.001), triquetrum (r = 0.81, p < 0.001), pisiform (r = 0.85, p < 0.001), capitate (r = 0.80, p < 0.001), lunate and pisiform (r = 0.81, p < 0.001), lunate and trapezium (r = 0.80, p < 0.001), and trapezium and capitate (r = 0.80, p < 0.001) (Table 2).
Figure 3 shows that the radial and ulnar measurements clustered together and were highly correlated, whereas the metacarpal measurements clustered together and were highly correlated. Similarly, the carpal measurements clustered together and were highly correlated.
Similarly, we noted strong positive correlations between the CT attenuation of the wrist/forearm bones and DEXA measurements. The strongest correlations between CT attenuations and DEXA L1–L4 BMD T-scores were with the CT attenuation of the first metacarpal (r = 0.26, p < 0.001), trapezium (r = 0.26, p < 0.001), and scaphoid (r = 0.32, p < 0.001) (Table 3).
The strongest correlations between the femoral neck BMD T-scores and the CT attenuations of the wrist/forearm bones were with the CT attenuations of the ulna (r = 0.40, p < 0.001), scaphoid (r = 0.47, p < 0.001), lunate (r = 0.40, p < 0.001), pisiform (r = 0.44, p < 0.001), trapezoid (r = 0.40, p < 0.001), and capitate (r = 46, p < 0.001). We also found that the strongest correlations between the total hip BMD T-score and CT attenuation of the wrist/forearm bones were with the CT attenuation of the scaphoid (r = 0.48, p < 0.001), pisiform (r = 0.43, p < 0.001), and first metacarpal (r = 0.40, p < 0.001). The strongest correlations with the L1–L4 BMD TBS were with the CT attenuation of the scaphoid (r = 0.24, p < 0.001), lunate (r = 0.23, p < 0.001), trapezium (r = 0.23, p < 0.001), and first metacarpal (r = 0.23, p < 0.001).
The optimal CT attenuation thresholds for each bone and the predictive performance of these optimal thresholds regarding predicting (i) osteoporosis, (ii) osteopenia/osteoporosis, (iii) femoral neck BMD T-score ≤ −2.5, and (iv) and femoral neck BMD T-score <−1 are shown in Table 4.

3.1. Predicting Osteoporosis

3.1.1. Training/Validation Dataset

We found that the CT attenuation of each bone was a significant predictor of osteoporosis in the training/validation dataset (Supplementary Figure S1). The highest AUC was for the hamate (optimal CT threshold 170.166 Hounsfield units (HU), AUC = 0.769), capitate (optimal CT threshold 248.039 HU, AUC = 0.763), and first metacarpal (optimal CT threshold −7.772 HU, AUC = 0.752). The radius 33% had an AUC of 0.705 in the training/validation dataset. The linear kernel SVM (AUC = 0.894) and radial basis function (RBF) kernel SVM (AUC = 0.987) had the highest AUCs of the machine learning models in the training/validation dataset.

3.1.2. Test Dataset

The performances of the CT attenuation thresholds obtained from the training/validation dataset were evaluated in the test dataset. The hamate (AUC = 0.393), capitate (AUC = 0.636), first metacarpal (AUC = 0.639), and radius 33% (AUC = 0.576) showed slightly lower predictive abilities in the test dataset compared to the training dataset. The RBF kernel SVM (AUC = 0.818) had a higher AUC than any of the CT attenuation thresholds for each bone. The RBF kernel SVM was better than the radius 33% model (p = 0.020).

3.2. Predicting Osteopenia/Osteoporosis

3.2.1. Training/Validation Dataset

When predicting osteopenia/osteoporosis in the training/validation dataset, we found that the CT attenuation of each bone was predictive of osteoporosis/osteopenia (Supplementary Figure S2). The first metacarpal (optimal CT threshold 27.779 HU, AUC = 0.823), scaphoid (optimal CT attenuation threshold 250.749 HU, AUC = 0.773), and lunate (optimal CT attenuation threshold 258.091 HU, AUC = 0.768) were the most predictive, whereas the radius UD (AUC = 0.528), third metacarpal (AUC = 0.529), and fourth metacarpal (AUC = 0.579) were the least predictive bones. The radius 33% had an AUC of 0.716 regarding predicting osteopenia/osteoporosis in the training/validation dataset. The RBF kernel SVM (AUC = 0.969) had the highest AUC of all the machine learning methods investigated regarding predicting osteopenia/osteoporosis in the training/validation dataset.

3.2.2. Test Dataset

In the test dataset, we found that the CT attenuations thresholds obtained from the training/validation dataset for the first metacarpal (AUC = 0.445), scaphoid (AUC = 0.651), and lunate (AUC = 0.433) were less predictive in the test dataset. The radius 33% CT attenuation threshold had an AUC of 0.563 in the test dataset. However, the RBF kernel SVM (AUC = 0.805) had the highest AUC and accuracy in the test dataset. The RBF kernel SVM (p = 0.068) was not significantly better than the radius 33% in the test dataset.

3.3. Predicting Femoral Neck BMD T-Score ≤ −2.5

3.3.1. Training/Validation Dataset

When predicting femoral neck BMD T-score ≤ −2.5, we found that the CT attenuation of each of the bones studied was predictive (Supplementary Figure S3). The fifth metacarpal (optimal CT attenuation threshold 24.690 HU, AUC = 0.737), pisiform (optimal CT attenuation threshold 121.626 HU, AUC = 0.736), and first metacarpal (optimal CT attenuation threshold 0.530 HU, AUC = 0.710) were the best bones for predicting femoral neck BMD T-score ≤ −2.5 in the training/validation dataset, whereas the CT attenuation of the radius (AUC = 0.569), ulna UD (AUC = 0.581), and lunate (AUC = 0.600) were the bones that were the worst predictors of femoral neck BMD T-score ≤ −2.5. The RBF kernel SVM (AUC = 0.997) had the highest AUC in the training/validation dataset.

3.3.2. Test Dataset

In the test dataset, we found that the CT attenuation thresholds obtained from the training/validation dataset for the fifth metacarpal (AUC = 0.398), pisiform (AUC = 0.415), first metacarpal (AUC = 0.618), and radius 33% (AUC = 0.426) were less predictive. The RBF kernel SVM had an AUC of 0.770 regarding identifying patients with a femoral neck BMD T-score ≤ −2.5. The RBF kernel SVM (AUC = 0.770, p < 0.001) model was better than the model using the radius 33% CT attenuation threshold.

3.4. Predicting Femoral Neck BMD T-Score < −1

3.4.1. Training/Validation Dataset

We found that a CT attenuation threshold of the trapezium (optimal CT attenuation threshold 165.624 HU, AUC = 0.722) had the best AUC to predict femoral neck BMD T-score <−1 in the training/validation dataset (Supplementary Figure S4). The scaphoid (optimal CT attenuation threshold 229.799 HU, AUC = 0.719) and pisiform (optimal CT attenuation threshold 221.709 HU, AUC = 0.714) were the next best predictors of having a femoral neck BMD T-score < −1. The radius 33% (AUC = 0.605), radius (AUC = 0.603), radius UD (AUC = 0.558), and hamate (AUC = 0.584) were the worst predictors of having a femoral neck BMD < −1 in the training/validation dataset. The RBF kernel SVM (AUC = 0.987) had the highest AUC of all the multivariable machine learning models regarding predicting femoral neck BMD T-score < −1 in the training/validation dataset.

3.4.2. Test Dataset

In the test dataset, we found that the CT attenuations thresholds obtained from the training/validation dataset for the scaphoid (AUC = 0.736), trapezium (AUC = 0.714), pisiform (AUC = 0.437), and radius 33% (AUC = 0.423) were less predictive in the test dataset. The RBF kernel SVM (AUC = 0.818, p < 0.001) had a significantly higher AUC than radius 33% (AUC = 0.423) in the test dataset.

4. Discussion

We found strong positive correlations between the CT attenuation of the wrist/forearm bones. There were stronger correlations between the CT attenuations of the radius and ulna measurements, stronger correlations between the carpal CT attenuations, and stronger correlations between the metacarpal CT attenuations. The CT attenuations of the wrist/forearm bones were largely positively correlated with the DEXA measurements. The RBF kernel SVM was best for predicting osteoporosis with a higher AUC (AUC = 0.818) than other models and statistically better than radius 33%. The RBF kernel SVM also had the highest AUC (AUC = 0.805) regarding predicting patients with osteoporosis/osteopenia. We found that the RBF kernel SVM was best for predicting a femoral neck BMD T-score ≤ −2.5 and was statistically better than the CT attenuation of radius 33%. The data also showed that the RBF kernel SVM was the best for predicting a femoral neck BMD T-score < −1 and statistically better than the CT attenuation of radius 33%.
These results have significant clinical implications. We showed that opportunistic screening for osteoporosis and osteopenia can be performed using routine CT scans of the wrist/forearm obtained for clinical care. We also showed that the accuracy of a machine learning model using the CT attenuation of multiple bones and clinical/demographic variables exceeded that of a single bone. We provided CT attenuation thresholds for each wrist/forearm bone that could be used to identify patients who should go on to get screened for osteoporosis or osteopenia/osteoporosis with a formal DEXA study. This has the potential to identify patients earlier in their course of BMD loss and get patients to clinical care earlier.
The radius transmits the majority of the force from the wrist to the elbow and is therefore often fractured when individuals with diminished BMD fall on an outstretched hand [6,7,27]. Forearm fractures, in particular, distal radius fractures, are often frailty fractures related to diminished BMD. Although the radius is often evaluated on DEXA studies for the prediction of osteopenia/osteoporosis and osteoporosis, we found that the CT attenuation of the radius was not the single best bone for predicting osteopenia/osteoporosis and osteoporosis in the analysis. Bone homeostasis is regulated by the activities of three cell types, namely, the osteoclasts, osteocytes, and osteoblasts, and is generally kept in dynamic equilibrium to maintain bone mass [28]. Our data suggested that all bones may be differentially affected when there is diminished bone mass because of the lack of perfect correlation between the CT attenuation of bones and the likely BMD of each bone [29]. Each bone also likely has a different trabecular bone structure in order to best suit its mechanical/structural demands, and this may affect the mean CT attenuation between bones.
The CT attenuation of the hamate had the best performance to differentiate patients with osteoporosis from those without osteoporosis. The patients in this study were aged 50 years or older, and likely had degenerative changes of the wrist. Degenerative changes of the carpus often affect the scaphoid, trapezium, and trapezoid (triscaphe and first carpometacarpal joint degenerative changes), scapholunate, and the pisotriquetral articulations. These degenerative changes are manifested by increased sclerosis and eburnation, and this increased sclerosis likely affects the CT attenuation of each bone. We hypothesized that the hamate was least likely to be affected by degenerative changes and this was likely the reason it had the best performance. One strength of the multivariable SVM and random forest models was that these models utilized the CT attenuation of all the bones and, therefore, had better performance than the CT attenuation of any individual bone.
A prior study noted that the second metacarpal cortical index can be used for opportunistic screening for osteoporosis [22]. We showed that the CT attenuation of the second metacarpal is predictive of osteoporosis and osteopenia/osteoporosis. We showed that the CT attenuation of the other metacarpals was also predictive of osteoporosis and osteopenia/osteoporosis. However, the CT attenuation of the metacarpals was less accurate than the multivariable machine learning models for the diagnosis of osteoporosis and osteopenia/osteoporosis.
One limitation of this study was the small sample size. We found that most patients with CT scans of the wrist/forearm at our institution did not have concurrent DEXA studies, which suggests that there is under-screening/inadequate screening for osteopenia and osteoporosis. Our method may help to identify patients who had a CT scan of the wrist/forearm who should go on to have a formal DEXA study to screen for osteopenia and osteoporosis in the near future. Patients all had their CT scans of the wrist/forearm performed using Siemens scanners, which is a limitation; however, a recently published article suggests that there is minimal bias in CT attenuation measurements between CT manufacturers [30]. This study was a retrospective study at a single multi-center tertiary care academic institution. Another limitation is that most patients were white; therefore, it is unclear how our results will port to other races/ethnicities. While we showed methods regarding how to categorize patients as osteoporotic or osteopenic/osteoporotic using the CT attenuation of the wrist/forearm bones based on DEXA studies, further work is required to evaluate how well these measurements predict future frailty fractures.

5. Conclusions

In summary, opportunistic screening for osteoporosis and osteopenia/osteoporosis could be performed using the CT attenuation of each of the bones from CT scans of the wrist/forearm. We used machine learning to show that using the CT attenuation of multiple bones was more accurate than using the CT attenuation of a single bone. DEXA scans currently evaluate only the lumbar spine and the hips to assess global bone mineral density. CT attenuation data from routine CT scans of the wrist and forearm can be used to identify patients at risk for osteoporosis who should go on to have formal screening for osteoporosis using DEXA scans.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/diagnostics12030691/s1, Figure S1: Performance of the CT attenuation of each bone and multivariable machine learning models to predict osteoporosis; Figure S2: Performance of the CT attenuation of each bone and multivariable machine learning models to predict osteopenia/osteoporosis; Figure S3: Performance of the CT attenuation of each bone and multivariable machine learning models to predict femoral neck BMD T-score ≤ −2.5; Figure S4: Performance of the CT attenuation of each bone and multivariable machine learning models to predict femoral neck BMD T-score < −1.

Author Contributions

Conceptualization, R.S.; methodology, R.S.; software, R.S.; validation, R.S. and C.d.l.G.-R.; formal analysis, R.S.; investigation, R.S.; resources, R.S.; data curation, C.d.l.G.-R.; writing—original draft preparation, R.S.; writing—review and editing, R.S. and C.d.l.G.-R.; visualization, R.S.; supervision, R.S.; project administration, R.S.; funding acquisition, R.S. All authors read and agreed to the published version of the manuscript.

Funding

R.S. was funded by R21 NIH/NIMH MH093415.

Institutional Review Board Statement

The study was reviewed and approved by the local IRB (IRB #21-000724).

Informed Consent Statement

The study protocol was reviewed and approved by the IRB board at the senior author’s institution. The need for signed informed consent from each patient was waived because of the retrospective nature of the study.

Data Availability Statement

Data available upon request from the corresponding author and IRB approval.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Warming, L.; Hassager, C.; Christiansen, C. Changes in bone mineral density with age in men and women: A longitudinal study. Osteoporos. Int. 2002, 13, 105–112. [Google Scholar] [CrossRef] [PubMed]
  2. Padilla Colón, C.J.; Molina-Vicenty, I.L.; Frontera-Rodríguez, M.; García-Ferré, A.; Rivera, B.P.; Cintrón-Vélez, G.; Frontera-Rodríguez, S. Muscle and Bone Mass Loss in the Elderly Population: Advances in diagnosis and treatment. J. Biomed. 2018, 3, 40–49. [Google Scholar] [CrossRef] [PubMed]
  3. Ji, M.X.; Yu, Q. Primary osteoporosis in postmenopausal women. Chronic. Dis. Transl. Med. 2015, 1, 9–13. [Google Scholar] [CrossRef] [PubMed]
  4. Melton, L.J., 3rd; Atkinson, E.J.; O’Connor, M.K.; O’Fallon, W.M.; Riggs, B.L. Bone density and fracture risk in men. J. Bone Miner. Res. 1998, 13, 1915–1923. [Google Scholar] [CrossRef]
  5. The European Prospective Osteoporosis Study (EPOS) Group. The relationship between bone density and incident vertebral fracture in men and women. J. Bone Miner. Res 2002, 17, 2214–2221. [Google Scholar] [CrossRef]
  6. Krolner, B.; Tondevold, E.; Toft, B.; Berthelsen, B.; Nielsen, S.P. Bone mass of the axial and the appendicular skeleton in women with Colles’ fracture: Its relation to physical activity. Clin. Physiol. 1982, 2, 147–157. [Google Scholar] [CrossRef]
  7. Mallmin, H.; Ljunghall, S. Distal radius fracture is an early sign of general osteoporosis: Bone mass measurements in a population-based study. Osteoporos. Int. 1994, 4, 357–361. [Google Scholar] [CrossRef]
  8. Marks, R. Hip fracture epidemiological trends, outcomes, and risk factors, 1970–2009. Int. J. Gen. Med. 2010, 3, 1–17. [Google Scholar] [CrossRef]
  9. Morri, M.; Ambrosi, E.; Chiari, P.; Orlandi Magli, A.; Gazineo, D.; D’Alessandro, F.; Forni, C. One-year mortality after hip fracture surgery and prognostic factors: A prospective cohort study. Sci. Rep. 2019, 9, 18718. [Google Scholar] [CrossRef]
  10. Mariconda, M.; Costa, G.G.; Cerbasi, S.; Recano, P.; Aitanti, E.; Gambacorta, M.; Misasi, M. The determinants of mortality and morbidity during the year following fracture of the hip: A prospective study. Bone Joint J. 2015, 97, 383–390. [Google Scholar] [CrossRef]
  11. Lewiecki, E.; Lane, N. Common mistakes in the clinical use of bone mineral density testing. Nat. Rev. Rheumatol. 2008, 4, 667–674. [Google Scholar] [CrossRef]
  12. Blake, G.M.; Fogelman, I. The role of DXA bone density scans in the diagnosis and treatment of osteoporosis. Postgrad. Med. J. 2007, 83, 509–517. [Google Scholar] [CrossRef]
  13. Lorente-Ramos, R.; Azpeitia-Armán, J.; Muñoz-Hernández, A.; García-Gómez, J.M.; Díez-Martínez, P.; Grande-Bárez, M. Dual-energy X-ray absorptiometry in the diagnosis of osteoporosis: A practical guide. AJR Am. J. Roentgenol. 2011, 196, 897–904. [Google Scholar] [CrossRef]
  14. Kanis, J.A.; McCloskey, E.V.; Johansson, H.; Oden, A.; Melton, L.J., 3rd; Khaltaev, N. A reference standard for the description of osteoporosis. Bone 2008, 42, 467–475. [Google Scholar] [CrossRef]
  15. Unnanuntana, A.; Gladnick, B.P.; Donnelly, E.; Lane, J.M. The assessment of fracture risk. J. Bone Joint Surg. Am. 2010, 92, 743–753. [Google Scholar] [CrossRef]
  16. Hans, D.; Goertzen, A.L.; Krieg, M.A.; Leslie, W.D. Bone microarchitecture assessed by TBS predicts osteoporotic fractures independent of bone density: The Manitoba study. J. Bone Miner. Res. 2011, 26, 2762–2769. [Google Scholar] [CrossRef] [PubMed]
  17. Winzenrieth, R.; Dufour, R.; Pothuaud, L.; Hans, D. A retrospective case-control study assessing the role of trabecular bone score in postmenopausal Caucasian women with osteopenia: Analyzing the odds of vertebral fracture. Calcif. Tissue Int. 2010, 86, 104–109. [Google Scholar] [CrossRef]
  18. Park, S.H.; Jeong, Y.M.; Lee, H.Y.; Kim, E.Y.; Kim, J.H.; Park, H.K.; Ahn, H.K. Opportunistic use of chest CT for screening osteoporosis and predicting the risk of incidental fracture in breast cancer patients: A retrospective longitudinal study. PLoS ONE 2020, 15, e0240084. [Google Scholar] [CrossRef]
  19. Lee, S.J.; Graffy, P.M.; Zea, R.D.; Ziemlewicz, T.J.; Pickhardt, P.J. Future Osteoporotic Fracture Risk Related to Lumbar Vertebral Trabecular Attenuation Measured at Routine Body CT. J. Bone Miner. Res. 2018, 33, 860–867. [Google Scholar] [CrossRef]
  20. Pickhardt, P.J.; Pooler, B.D.; Lauder, T.; del Rio, A.M.; Bruce, R.J.; Binkley, N. Opportunistic screening for osteoporosis using abdominal computed tomography scans obtained for other indications. Ann. Intern. Med. 2013, 158, 588–595. [Google Scholar] [CrossRef]
  21. Therkildsen, J.; Nissen, L.; Jørgensen, H.S.; Thygesen, J.; Ivarsen, P.; Frost, L.; Isaksen, C.; Langdahl, B.L.; Hauge, E.-M.; Boettcher, M.; et al. Thoracic bone mineral density derived from cardiac CT is associated with greater fracture rate. Radiology 2020, 296, 499–508. [Google Scholar] [CrossRef]
  22. Patel, B.; Aqil, A.; Riaz, O.; Jeffers, R.; Dickson, D. The 2nd Metacarpal Cortical Index as a Simple Screening Tool for Osteopenia. J. Bone Metab. 2020, 27, 261–266. [Google Scholar] [CrossRef]
  23. Fedorov, A.; Beichel, R.; Kalpathy-Cramer, J.; Finet, J.; Fillion-Robin, J.-C.; Pujol, S.; Bauer, C.; Jennings, D.; Fennessy, F.M.; Sonka, M.; et al. 3D Slicer as an Image Computing Platform for the Quantitative Imaging Network. Magn. Reson. Imaging 2012, 30, 1323–1341. [Google Scholar] [CrossRef]
  24. Mehta, S.D.; Sebro, R. Computer-Aided Detection of Incidental Lumbar Spine Fractures from Routine Dual-Energy X-ray Absorptiometry (DEXA) Studies Using a Support Vector Machine (SVM) Classifier. J. Digit. Imaging 2020, 33, 204–210. [Google Scholar] [CrossRef]
  25. Mehta, S.D.; Sebro, R. Random forest classifiers aid in the detection of incidental osteoblastic osseous metastases in DEXA studies. Int. J. Comput. Assist. Radiol. Surg. 2019, 14, 903–909. [Google Scholar] [CrossRef] [PubMed]
  26. DeLong, E.R.; DeLong, D.M.; Clarke-Pearson, D.L. Comparing the areas under two or more correlated receiver operating characteristic curves: A nonparametric approach. Biometrics 1988, 44, 837–845. [Google Scholar] [CrossRef]
  27. Markolf, K.L.; Dunbar, A.M.; Hannani, K. Mechanisms of load transfer in the cadaver forearm: Role of the interosseous membrane. J. Hand Surg. Am. 2000, 25, 674–682. [Google Scholar] [CrossRef]
  28. Al-Bari, A.A.; Al Mamun, A. Current advances in regulation of bone homeostasis. FASEB Bioadv. 2020, 2, 668–679. [Google Scholar] [CrossRef]
  29. Sebro, R.; Ashok, S.S. A Statistical Approach Regarding the Diagnosis of Osteoporosis and Osteopenia From DXA: Are We Underdiagnosing Osteoporosis? JBMR Plus 2021, 5, e10444. [Google Scholar] [CrossRef]
  30. Jacobsen, M.C.; Schellingerhout, D.; Wood, C.A.; Tamm, E.P.; Godoy, M.C.; Sun, J.; Cody, D.D. Intermanufacturer Comparison of Dual-Energy CT Iodine Quantification and Monochromatic Attenuation: A Phantom Study. Radiology 2018, 287, 224–234. [Google Scholar] [CrossRef]
Figure 1. Semi-automated volumetric segmentation mask of the bones of the wrist and forearm. Coronal mask of the wrist/forearm demonstrating segmentation of the distal radius, ulna, carpal bones, and bases of the metacarpals.
Figure 1. Semi-automated volumetric segmentation mask of the bones of the wrist and forearm. Coronal mask of the wrist/forearm demonstrating segmentation of the distal radius, ulna, carpal bones, and bases of the metacarpals.
Diagnostics 12 00691 g001
Figure 2. Box and whisker plots of the CT attenuation of each bone of the wrist and forearm by diagnosis: 2—osteoporosis, 1—osteopenia, and 0—normal.
Figure 2. Box and whisker plots of the CT attenuation of each bone of the wrist and forearm by diagnosis: 2—osteoporosis, 1—osteopenia, and 0—normal.
Diagnostics 12 00691 g002
Figure 3. Hierarchical cluster analysis of the correlations between the CT attenuation of each bone of the wrist and forearm. Radius—distal third of the radius; radius UD—ultradistal radius (radius epiphysis and metaphysis); radius 33%—distal third of the radial shaft; ulna—distal third of the ulna; ulna UD—distal ulna (ulnar epiphysis and metaphysis); ulna 33%—distal third of the ulnar shaft; 1st metacarpal—proximal third of the first metacarpal; 2nd metacarpal—proximal third of the second metacarpal; 3rd metacarpal—proximal third of the third metacarpal; 4th metacarpal—proximal third of the fourth metacarpal; 5th metacarpal—proximal third of the fifth metacarpal.
Figure 3. Hierarchical cluster analysis of the correlations between the CT attenuation of each bone of the wrist and forearm. Radius—distal third of the radius; radius UD—ultradistal radius (radius epiphysis and metaphysis); radius 33%—distal third of the radial shaft; ulna—distal third of the ulna; ulna UD—distal ulna (ulnar epiphysis and metaphysis); ulna 33%—distal third of the ulnar shaft; 1st metacarpal—proximal third of the first metacarpal; 2nd metacarpal—proximal third of the second metacarpal; 3rd metacarpal—proximal third of the third metacarpal; 4th metacarpal—proximal third of the fourth metacarpal; 5th metacarpal—proximal third of the fifth metacarpal.
Diagnostics 12 00691 g003
Table 1. Comparison between training/validation and test datasets.
Table 1. Comparison between training/validation and test datasets.
VariableAll
(N = 196)
Training/
Validation Dataset (N = 96)
Test Dataset
(N = 100)
p-Value
Age64.9 (8.7)64.5 (8.8)65.2 (8.7)0.593
Race/ethnicity 0.248
American Indian/Alaskan native1 (0.5%)1 (1.0%)0 (0.0%)
Asian2 (1.0%)0 (0.0%)2 (2.0%)
Black/African-American1 (0.5%)0 (0.0%)1 (1.0%)
Hispanic2 (1.0%)0 (0.0%)2 (2.0%)
Other1 (0.5%)0 (0.0%)1 (1.0%)
White189 (96.4%)95 (99.0%)94 (94.0%)
Gender, male23 (11.7%)12 (12.5%)11 (11.0%)0.826
Height (m)1.65 (0.08)1.65 (0.08)1.64 (0.08)0.47
Weight (kg)76.2 (19.2)76.5 (19.1)75.8 (19.4)0.8
BMI (kg/m2)27.9 (6.4)27.9 (6.4)27.9 (6.4)0.991
Diagnosis 1
Osteoporosis54 (27.6%)26 (27.1%)28 (28.0%)
Osteopenia116 (59.2%)57 (59.4%)59 (59.0%)
Normal26 (13.3%)13 (13.5%)13 (13.0%)
BMI—body mass index; osteoporosis—minimum BMD T-score ≤ −2.5; osteopenia—−2.5 < minimum BMD T-score < −1; normal—minimum BMD T-score ≥ −1; ***—p-value < 1 × 10−3; **—1 × 10−3 < p-value < 1 × 10−2; *—1 × 10−2 < p-value < 5 × 10−2.
Table 2. Correlations between the CT attenuations at different bony sites in males and females.
Table 2. Correlations between the CT attenuations at different bony sites in males and females.
Males Top DiagonalRadiusRadius UDRadius 33%UlnaUlna UDUlna 33%ScaphoidLunateTriquetrumPisiformTrapeziumTrapezoidCapitateHamate1 MC2 MC3 MC4 MC5 MC
Females bottom diagonal
Radius-0.82 ***0.79 ***0.55 **0.46 *0.320.370.360.310.340.400.56 **0.50 *0.300.46 *0.120.08−0.190.05
Radius UD0.74 ***-0.45 *0.51 *0.55 **0.220.54 **0.52 *0.56 **0.47 *0.57 **0.61 **0.59 **0.390.64 **0.240.270.030.12
Radius 33%0.73 ***0.40 ***-0.57 **0.350.56 **0.250.180.040.270.240.410.310.210.300.00−0.09−0.030.04
Ulna0.61 ***0.39 ***0.64 ***-0.84 ***0.84 ***0.62 **0.57 **0.370.65 ***0.62 **0.62 **0.65 ***0.49 *0.43 *0.320.090.270.10
Ulna UD0.46 ***0.43 ***0.35 ***0.74 ***-0.62 **0.64 **0.70 ***0.46 *0.62 **0.67 ***0.68 ***0.67 ***0.390.57 **0.320.240.300.24
Ulna 33%0.46 ***0.24 **0.69 ***0.83 ***0.52 ***-0.54 **0.43 *0.220.59 **0.50 *0.360.410.42 *0.350.330.060.49 *0.23
Scaphoid0.38 ***0.28 ***0.25 ***0.42 ***0.45 ***0.30 ***-0.84 ***0.81 ***0.85 ***0.78 ***0.71 ***0.80 ***0.66 ***0.46 *0.55 **0.45 *0.320.27
Lunate0.27 ***0.15 *0.23 **0.42 ***0.42 ***0.30 ***0.60 ***-0.76 ***0.81 ***0.80 ***0.75 ***0.78 ***0.58 **0.59 **0.55 **0.50 *0.330.42 *
Triquetrum0.37 ***0.33 ***0.22 **0.41 ***0.51 ***0.26 ***0.66 ***0.64 ***-0.69 ***0.71 ***0.59 **0.77 ***0.61 **0.50 *0.59 **0.55 **0.290.45 *
Pisiform0.28 ***0.25 ***0.20 ***0.37 ***0.46 ***0.34 ***0.59 ***0.55 ***0.59 ***-0.74 ***0.72 ***0.69 ***0.56 **0.54 **0.51 *0.310.370.33
Trapezium0.31 ***0.23**0.19*0.38 ***0.41 ***0.31 ***0.66 ***0.57 ***0.61 ***0.51 ***-0.74 ***0.80 ***0.53 *0.71 ***0.59 **0.68 ***0.51 *0.41
Trapezoid0.40 ***0.29 ***0.26 ***0.40 ***0.38 ***0.24 **0.66 ***0.57 ***0.58 ***0.56 ***0.59 ***-0.74 ***0.380.53 **0.50 *0.390.140.19
Capitate0.39 ***0.31 ***0.24 **0.45 ***0.47 ***0.33 ***0.71 ***0.72 ***0.73 ***0.62 ***0.66 ***0.74 ***-0.62 **0.55 **0.410.45 *0.160.32
Hamate0.38 ***0.34 ***0.25 ***0.46 ***0.50 ***0.31 ***0.66 ***0.65 ***0.69 ***0.62 ***0.62 ***0.69 ***0.81 ***-0.280.290.290.180.17
1 MC0.41 ***0.39 ***0.27 ***0.39 ***0.38 ***0.31 ***0.44 ***0.24 **0.34 ***0.38 ***0.48 ***0.35 ***0.35 ***0.38 ***-0.46 *0.52 *0.51 *0.64 **
2 MC0.32 ***0.27 ***0.25 ***0.37 ***0.29 ***0.33 ***0.35 ***0.22 **0.22 **0.34 ***0.31 ***0.42 ***0.37 ***0.37 ***0.48 ***-0.72 ***0.63 **0.68 ***
3 MC0.23 ***0.20 **0.26 ***0.22 **0.18 *0.23 **0.24 **0.060.130.24 ***0.15 *0.28 ***0.20 **0.26 ***0.36 ***0.63 ***-0.60 **0.68 ***
4 MC0.27 ***0.32 ***0.23 **0.19 *0.22 **0.19*0.30 ***0.150.24**0.21**0.19*0.25 ***0.24**0.34 ***0.45 ***0.55 ***0.68 ***-0.65 ***
5 MC0.34 ***0.31 ***0.30 ***0.28 ***0.28 ***0.26 ***0.24 **0.140.22 **0.27 ***0.140.20 **0.20 **0.26 ***0.37 ***0.49 ***0.54 ***0.70 ***-
Radius—distal third of the radius; radius UD—ultradistal radius (radius epiphysis/metaphysis); radius 33%—distal third of the radial shaft; ulna—distal third of the ulna; ulna UD—distal ulna (ulnar epiphysis/metaphysis); ulna 33%—distal third of the ulnar shaft; 1 MC—proximal third of the first metacarpal; 2 MC—proximal third of the second metacarpal; 3 MC—proximal third of the third metacarpal; 4 MC—proximal third of the fourth metacarpal; 5 MC—proximal third of the fifth metacarpal; ***—p-value < 1 × 10−3; **—1 × 10−3 < p-value < 1 × 10−2; *—1 × 10−2 < p-value < 5 × 10−2.
Table 3. Correlations between the CT attenuation at different bony sites and DEXA measurements.
Table 3. Correlations between the CT attenuation at different bony sites and DEXA measurements.
AllL1–L4 BMDL1–L4 BMD
T-Score
L1–L4 TBS
(N = 133)
Femoral Neck BMDFemoral Neck BMD T-ScoreTotal Hip BMDTotal Hip BMD T-Score
Radius0.040.030.120.26 ***0.24 ***0.29 ***0.26 ***
Radius UD0.00−0.020.080.23 ***0.17 *0.18 *0.14 *
Radius 33%0.110.130.110.23 **0.25 ***0.31 ***0.29 ***
Ulna0.23 **0.24 ***0.18 *0.40 ***0.40 ***0.43 ***0.38 ***
Ulna UD0.27 ***0.17*0.22*0.40 ***0.39 ***0.42 ***0.38 ***
Ulna 33%0.24 ***0.24 ***0.19 *0.34 ***0.37 ***0.37 ***0.35 ***
Scaphoid0.26 ***0.32 ***0.24 **0.39 ***0.47 ***0.48 ***0.48 ***
Lunate0.22 **0.23 **0.23 **0.36 ***0.40 ***0.40 ***0.38 ***
Triquetrum0.14 *0.17 *0.18 *0.33 ***0.38 ***0.40 ***0.39 ***
Pisiform0.18 *0.25 ***0.17 *0.40 ***0.44 ***0.47 ***0.43 ***
Trapezium0.28 ***0.26 ***0.23 **0.32 ***0.38 ***0.40 ***0.37 ***
Trapezoid0.15 *0.18 *0.150.35 ***0.40 ***0.38 ***0.34 ***
Capitate0.19 **0.24 ***0.21 *0.39 ***0.46 ***0.43 ***0.42 ***
Hamate0.21 **0.22 **0.140.30 ***0.35 ***0.35 ***0.34 ***
1 MC0.24 ***0.26 ***0.170.36 ***0.38 ***0.43 ***0.40 ***
2 MC0.24 ***0.19 **0.23 **0.23 **0.31 ***0.37 ***0.34 ***
3 MC0.18 *0.040.21 *0.15 *0.17 *0.21 **0.19 **
4 MC0.22 **0.100.150.17 *0.23 **0.27 ***0.27 ***
5 MC0.27 ***0.16 *0.160.29 ***0.32 ***0.38 ***0.37 ***
Radius—distal third of the radius; radius UD—ultradistal radius (radius epiphysis/metaphysis); radius 33%—distal third of the radial shaft; ulna—distal third of the ulna; ulna UD—distal ulna (ulnar epiphysis/metaphysis); ulna 33%—distal third of the ulnar shaft; 1 MC—proximal third of the first metacarpal; 2 MC—proximal third of the second metacarpal; 3 MC—proximal third of the third metacarpal; 4 MC—proximal third of the fourth metacarpal; 5 MC—proximal third of the fifth metacarpal; ***—p-value < 1 × 10−3; **—1 × 10−3 < p-value < 1 × 10−2; *—1 × 10−2 < p-value < 5 × 10−2.
Table 4. Performance of the CT attenuation of each bone and multivariable machine learning models to predict osteoporosis and osteopenia/osteoporosis.
Table 4. Performance of the CT attenuation of each bone and multivariable machine learning models to predict osteoporosis and osteopenia/osteoporosis.
Test Dataset
OsteoporosisTraining/
Validation
Dataset CT
Attenuation Threshold
AUCSensitivitySpecificityAUCAccuracyPositive Predictive Value (PPV)Negative Predictive Value (NPV)
Radius90.1790.7080.5000.6390.5690.6000.3500.767
Radius UD154.9980.7250.6070.6250.6160.6200.3860.804
Radius 33%−13.7170.7050.5000.6530.5760.6100.3590.770
Ulna67.1210.7190.7500.6670.7080.6900.4670.873
Ulna UD98.4460.7320.5000.8060.6530.7200.5000.806
Ulna 33%3.8720.6690.7500.6110.6810.6500.4290.863
Scaphoid247.5920.7630.5710.5830.5770.5800.3480.778
Lunate248.3870.7620.001.000.3650.720-0.720
Triquetrum207.8820.7300.001.000.3900.720-0.720
Pisiform162.2980.7530.7140.6530.6840.6700.4440.855
Trapezium141.8240.7340.001.000.3830.720-0.720
Trapezoid231.0700.6990.5000.7220.6110.6600.4120.788
Capitate248.0390.7630.5360.7360.6360.6800.4410.803
Hamate170.1660.7690.001.000.3930.720-0.720
1 MC−7.7720.7520.5000.7780.6390.7000.4670.800
2 MC16.0230.6860.001.000.4150.720-0.720
3 MC61.5550.5650.001.000.4660.720-0.720
4 MC50.8370.6000.001.000.4150.720-0.720
5 MC−34.8600.5660.001.000.4080.720-0.720
Linear kernel SVM 0.8940.8830.4350.6800.7800.8400.526
Radial basis function kernel SVM 0.9870.5840.9570.8180.6700.9780.407
Sigmoid kernel SVM 0.6270.8440.7390.8180.8200.9150.586
Random Forest classifier 0.5020.9870.0870.5370.7800.7840.667
Osteopenia/OsteoporosisTraining/
Validation
Dataset CT
Attenuation
Threshold
AUCSensitivitySpecificityAUCAccuracyPositive Predictive Value (PPV)Negative Predictive Value (NPV)
Radius149.1990.6350.2620.7780.5200.3290.8890.135
Radius UD160.4960.5280.001.000.4720.129-0.129
Radius 33%10.9420.7160.4590.6670.5630.4860.9030.154
Ulna117.2590.7360.001.000.4320.129-0.129
Ulna UD162.0880.6430.7050.5560.6300.6860.9150.217
Ulna 33%73.3650.7080.001.000.4540.129-0.129
Scaphoid250.7490.7730.5250.7780.6510.5570.9410.194
Lunate258.0910.7680.001.000.4330.129-0.129
Triquetrum213.9980.6100.001.000.3920.129-0.129
Pisiform220.0410.7540.001.000.4230.129-0.129
Trapezium183.7380.7170.001.000.3100.129-0.129
Trapezoid269.5940.7260.6560.7780.7170.6710.9520.250
Capitate294.0580.7550.6230.8890.7560.6570.9740.258
Hamate171.5030.6730.001.000.4230.129-0.129
1 MC27.7790.8230.001.000.4450.129-0.129
2 MC30.5840.7520.7210.8890.8050.7430.9780.320
3 MC31.1970.5290.001.000.4090.129-0.129
4 MC55.3760.5790.7700.5560.6630.7430.9220.263
5 MC52.1120.6150.001.000.4070.390-0.390
Linear kernel SVM 0.8560.4430.8890.6740.6200.8710.507
Radial basis function kernel SVM 0.9690.8850.6670.8050.8000.8060.788
Sigmoid kernel SVM 0.5420.6070.7780.7160.6700.8040.556
Random Forest classifier 0.5110.9670.2220.5950.6800.6630.818
Femoral Neck BMD ≤ −2.5Training/
Validation
Dataset CT
Attenuation
Threshold
AUCSensitivitySpecificityAUCAccuracyPositive Predictive Value (PPV)Negative Predictive Value (NPV)
Radius132.4950.5690.001.000.3940.810-0.810
Radius UD184.1540.6180.7890.5310.6600.5800.2830.915
Radius 33%20.9080.6250.001.000.4260.810-0.810
Ulna67.1210.6030.7890.5560.6730.6000.2940.918
Ulna UD82.7300.5810.5260.7900.6580.7400.3700.877
Ulna 33%35.5200.6210.001.000.3750.810-0.810
Scaphoid202.9160.6570.6320.6790.6550.6700.3160.887
Lunate224.8380.6840.5260.8640.6950.8000.4760.886
Triquetrum208.3340.6670.6320.7280.6800.7100.3530.894
Pisiform121.6260.7360.001.000.4150.810-0.810
Trapezium149.5970.6270.6320.6910.6610.6800.3240.889
Trapezoid207.9530.6630.6320.6790.6550.6700.3160.887
Capitate248.0390.6470.7370.6670.7020.6800.3410.915
Hamate185.7430.6000.8420.5680.7050.6200.3140.939
1 MC0.5300.7100.5790.6420.6100.6300.2750.867
2 MC−7.2730.6810.5260.6300.5780.6100.2500.850
3 MC−47.2510.6090.8950.1360.5150.2800.1950.846
4 MC−13.1460.6720.001.000.4580.810-0.810
5 MC24.6900.7370.001.000.3980.810-0.810
Linear kernel SVM 0.9150.9470.5930.7950.6600.5350.980
Radial basis function kernel SVM 0.9970.5790.8640.7700.8100.5000.897
Sigmoid kernel SVM 0.7360.9470.5310.7490.6100.3210.977
Random Forest classifier 0.4890.4210.9010.6610.8100.5000.869
Femoral Neck BMD < −1Training/
Validation
Dataset CT
Attenuation
Threshold
AUCSensitivitySpecificityAUCAccuracyPositive Predictive Value (PPV)Negative Predictive Value (NPV)
Radius130.3360.6030.001.000.4150.270-0.270
Radius UD163.2090.5580.001.000.4920.270-0.270
Radius 33%10.9420.6050.001.000.4230.270-0.270
Ulna94.0090.6470.7400.6520.6960.7200.8570.486
Ulna UD185.5440.6840.001.000.3630.270-0.270
Ulna 33%27.4060.6180.7270.7390.7330.7300.8830.500
Scaphoid229.7990.7190.5580.9130.7360.6600.9530.439
Lunate268.1930.7070.001.000.3310.270-0.270
Triquetrum287.3660.6410.8310.5650.6980.7600.8360.556
Pisiform221.7090.7140.001.000.4370.270-0.270
Trapezium165.6240.7220.5580.8700.7140.6400.9110.418
Trapezoid236.0410.6930.6100.6960.6530.6400.8490.404
Capitate257.4990.6930.5450.8700.7080.7900.8420.625
Hamate160.0720.5840.001.000.2990.270-0.270
1 MC26.3900.7100.7140.6090.6610.6800.8250.432
2 MC9.5760.7000.6230.8700.7460.6800.9180.451
3 MC54.5740.4910.001.000.4240.270-0.270
4 MC5.1990.6160.001.000.4270.270-0.270
5 MC1.2940.6740.5970.6960.6470.6300.8460.396
Linear kernel SVM 0.8950.4680.8260.6780.5500.9000.317
Radial basis function kernel SVM 0.9870.5840.9570.8180.6700.9780.407
Sigmoid kernel SVM 0.6270.8440.7390.8180.8200.9150.586
Random Forest classifier 05020.9870.0430.5150.7700.7760.500
Radius—distal third of the radius; Radius UD—ultradistal radius (radius epiphysis/metaphysis); Radius 33%—distal third of the radial shaft; Ulna—distal third of the ulna; Ulna UD—distal ulna (ulnar epiphysis/metaphysis); Ulna 33%—distal third of the ulnar shaft; 1 MC—proximal third of the first metacarpal; 2 MC—proximal third of the second metacarpal; 3 MC—proximal third of the third metacarpal; 4 MC—proximal third of the fourth metacarpal; 5 MC—proximal third of the fifth metacarpal; -—Undefined.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Sebro, R.; De la Garza-Ramos, C. Machine Learning for Opportunistic Screening for Osteoporosis from CT Scans of the Wrist and Forearm. Diagnostics 2022, 12, 691. https://doi.org/10.3390/diagnostics12030691

AMA Style

Sebro R, De la Garza-Ramos C. Machine Learning for Opportunistic Screening for Osteoporosis from CT Scans of the Wrist and Forearm. Diagnostics. 2022; 12(3):691. https://doi.org/10.3390/diagnostics12030691

Chicago/Turabian Style

Sebro, Ronnie, and Cynthia De la Garza-Ramos. 2022. "Machine Learning for Opportunistic Screening for Osteoporosis from CT Scans of the Wrist and Forearm" Diagnostics 12, no. 3: 691. https://doi.org/10.3390/diagnostics12030691

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop