Retrospective Study of Bone Mineral Density through Phantom-Less Quantitative CT ()
1. Introduction
Osteoporosis is a bone disease characterized by low bone mineral density (BMD) and micro-architectural deterioration of bone tissue [1]. BMD measurement by using dual-energy X-ray absorptiometry (DXA) is defined by the World Health Organization (WHO) Working Group as the gold standard and is clinically used to diagnose osteoporosis [2]. BMD is clearly one of the major determinants of bone strength and fracture risk [3]-[5]. Quantitative computed tomography (QCT) has been shown to provide a similar accurate fracture prediction. In QCT, there are two methods to measure the BMD using calibration phantoms. One is called the simultaneous calibration method, in which the subject is scanned on top of a calibration phantom, which contains several inserts of known values of hydroxyapatite (HA) or similar substances characteristic of bone material. Based on the CT values measured in these inserts, CT values measured in the bone can be converted to BMD. This method has the advantage of minimizing the difference in X-ray beam properties between the bone of interest and the inserts of the calibration phantom. Another is called the asynchronous calibration method, where a calibration phantom is measured separately from the subject. This method is suitable for current modern CTs, which have stable performance. This technique requires additional phantom imaging and therefore limits its clinical applicability. All of these phantom-based calibration methods are rarely used in a clinical setting, since they require additional time and clinical resources and produce image artifacts. Instead of using these calibration phantoms, certain tissues called internal calibration materials within the body, such as muscle, blood, and subcutaneous adipose tissue, can be used. It is assumed that the composition of these internal calibration materials does not vary among subjects. The development and validation of internal calibration techniques were defined as a priority by the International Society of Clinical Densitometry [6] to move the field towards clinically implemented skeletal health assessment. Internal density calibration techniques have been in development for several years. However, these techniques, either commercially available or custom-made, rely on the basic assumptions of ground truth values for each calibration tissue. Weaver et al. determined the ground truth reference values of −69 mg/cc for fat and 77 mg/cc for muscle and then used these reference values for the internal density calibrations [7]. Similarly, using clinically available CT scanners, the regions of interest (ROI) of subcutaneous fat and paraspinal muscle are selected and used to measure BMD [8] [9]. Furthermore, ROIs of air and blood or fat have also been paired with developed reference values of equivalent BMD to measure bone’s BMD [10]. These phantom-less methods have been applied to abdominal CT images to determine the BMD of lumbar spines [7] [11]. In our study, we applied the phantom-less method to clinical chest and abdominal CT images to determine the BMD of thoracic and lumbar spines. The first step of this study was to determine the reference values of aortic blood in the chest based on a linear model, which can be used to determine BMD values of the thoracic spine. The second step was to determine the thoracic BMD values based on the internal reference tissues of aortic blood and anterior subcutaneous fat located in the chest and determine the lumbar BMD values based on the internal reference tissues of psoas muscle and anterior subcutaneous fat located in the abdomen. The third step was to analyze the correlation between these thoracic and lumbar BMD values.
2. Materials and Methods
2.1. Study Subjects
(a) (b)
Figure 1. The mean HU values in three regions of interest (ROI) were measured including right psoas muscle and anterior subcutaneous fat in a soft tissue window for internal calibration, and lumbar vertebrae in the bone window.
(a) (b)
Figure 2. The mean CT values in three regions of interest (ROI) were measured including aortic blood and anterior subcutaneous fat in a soft tissue window for internal calibration, and thoracic vertebrae in the bone window.
CT images of the chest and abdomen without contrast of 400 females were studied retrospectively. Their ages ranged from 20 to 80 years old with a median age of 59 years old. Exclusions were made for conditions affecting BMD, such as chronic kidney disease, vitamin D deficiency, long-term steroid use, hyperparathyroidism and hyperthyroidism, cirrhosis, acute leukemia, and diffuse spinal metastases. These clinical images were acquired by GE scanners (Revolution, GE Healthcare, Waukesha, WI, USA). The chest and abdomen were scanned using a standard helical imaging protocol (120 kVp, 5.0 mm slice thickness, and standard reconstruction kernel) with exposure set by automatic exposure control. Both thoracic and lumbar vertebral Hounsfield unit (HU) values were measured in the mid-vertebral trabecular region from the T11-L4 vertebrae. The average HU value of the thoracic spine and the average HU value of the lumbar spine were then calculated. Additional measurements were made including the HU values of the anterior subcutaneous fat of the chest and abdomen, right psoas muscle, and aortic blood in the chest (Figure 1 and Figure 2). To enable phantom-less calibration for QCT, aortic blood and psoas muscle were selected as internal calibration references. These tissues were chosen based on their anatomical proximity to the spine [8] [10], consistent composition [6] [8], and stable attenuation properties in CT imaging, as supported by previous studies [6] [9] [10] [12].
2.2. Determination of the Reference Value for Aortic Blood
In this study, the relationship between BMD of soft tissues, such as subcutaneous fat, psoas muscle and aortic blood, and their measured CT HU values are assumed as a linear relationship:
(1)
where the constants a and b are obtained by phantom-less calibration with the ground truth reference values of
for subcutaneous fat and
for paraspinal muscle:
(2)
and
(3)
where <…> stands for an ensemble average over N subjects. After determining the constants a and b, the mean BMD value of aortic blood can be calculated by using the HU values of aortic blood, which is
.
2.3. Calculation of the BMD Values of Thoracic Spines
Then, the values of thoracic BMD can be determined by the HU values of the thoracic spine vertebrae, aortic blood and anterior subcutaneous fat on chest CT images:
Similarly, the values of lumbar BMD can be determined by the HU values of the lumbar spine vertebrae, psoas muscle and anterior subcutaneous fat on abdominal CT images:
2.4. Statistical Analysis
The Pearson correlation coefficient for the BMD values between thoracic and lumbar spines was calculated. The Pearson correlation coefficient is given by:
(4)
where
is the
,
is the
of the subject I,
is the mean of
,
is the mean of
, and
is the total number of subjects. Then, linear regression models and their R-squared value were also calculated. The R-squared value shows how well the data fits the regression model.
The receiver operating characteristic (ROC) of thoracic BMD values was studied with lumbar BMD values considered the gold standard for osteoporosis diagnosis. The area under the ROC curve was used to evaluate the suitability of thoracic BMD values for osteoporosis diagnosis.
3. Results
For the 400 female subjects examined, their ages ranged from 20 to 80 years old. Both thoracic and lumbar vertebral CT HU values were measured in the mid-vertebral trabecular region. Thoracic BMD values ranged from 60 - 350 mg/cm3, while lumbar BMD values ranged from 60 - 350 mg/cm3 in most subjects. Figure 3 and Figure 4 show the changes in BMD values of the thoracic and lumbar spine with age. These results indicate that spinal BMD gradually decreases with increasing age. The Pearson correlation coefficient between thoracic and lumbar BMD is 0.953 (p < 0.0001), and linear regression line with zero intercept shows R-squared value = 0.9893, as shown in Figure 5. These results show that there is a very high correlation between thoracic and lumbar BMD.
For subjects aged 20 - 40, lumbar BMD exhibited a mean of 233 mg/cm3 and standard deviation of 55 mg/cm3. A threshold lumbar BMD value of 123 mg/cm3 corresponded to a T-score of −2 and was closely aligned with ACR guidelines (BMD > 120 mg/cm3) [13]. If we use this threshold of lumbar BMD value as a criterion for osteoporosis diagnosis, the area under the ROC curve of thoracic BMD values would be 0.969, as shown in Figure 6. These results indicate that thoracic BMD values can also potentially replace lumbar BMD values for osteoporosis diagnosis. In routine clinical settings, especially for osteoporosis screening, high sensitivity is desired, which could be achieved by increasing the BMD threshold. For example, according to the above ROC curve, if the threshold of thoracic BMD value for osteoporosis diagnosis/screening is set to 150 mg/cm3, the sensitivity for osteoporosis diagnosis can be 0.91, and the specificity is 0.89.
Figure 7 illustrates the variations in lumbar BMD among different racial groups. The median lumbar BMD values were 197 mg/cm3, 188 mg/cm3, 185 mg/cm3 and 167 mg/cm3 for Hispanic, Asian, Black and White. The results show that the median BMD values decrease sequentially from Hispanic to Asian to Black and finally to White, with Hispanics exhibiting the highest BMD and Whites the lowest.
Figure 3. The results show that lumbar BMD values decreased with the subjects’ ages, r = −2.1895, p < 0.0001.
Figure 4. The results show that thoracic BMD values also decreased with the subjects’ ages, r = −2.012, p < 0.0001.
Figure 5. The Pearson correlation coefficient between thoracic and lumbar BMD is 0.953 (p < 0.0001), and linear regression line with zero intercept shows R-squared value = 0.9893.
Figure 6. The area under the ROC curve of thoracic BMD is 0.969.
Figure 7. Lumbar BMD variations between different racial groups.
4. Discussion
Quantifying BMD is important to monitor and evaluate bone health status. Currently, dual-energy X-ray absorptiometry (DXA) is the gold standard for assessing BMD. However, it has a few drawbacks including spuriously elevated BMD measurements in patients with degenerative disease, compression fractures, and vascular calcifications, which subsequently underestimates the risk of fractures. Measuring vertebral HU values on CT scans has been proposed as an alternate method of assessing BMD and bone quality [14]. The HU is a quantitative measurement of X-ray mass attenuation used by radiologists in the interpretation of CT images. Different tissue densities correlate with different HU values. The use of HU to objectively measure tissue densities has aided radiologists in the interpretation of CT images and diagnosis of disease and may also be used to assess bone quality and BMD [14] [15]. One study showed that thoracic and lumbar quantitative CTs provide a more sensitive method for detecting bone mineral loss when compared to DXA [16]. Furthermore, thoracic CT without contrast is commonly utilized for lung nodule and cancer screening. If this imaging modality can simultaneously assess bone health, it would significantly enhance the clinical value of these CT examinations, providing dual benefits for patients by detection of osteoporosis as well as the early detection of lung nodules and cancer.
Our study reveals variations in BMD among different racial groups, likely attributed to differences in diet, genetics, and lifestyle factors [17].
There are several reasons for choosing aortic blood and psoas muscle in our study: 1) Proximity to the spine: The aortic blood and psoas muscle are located adjacent to the lumbar spine, minimizing the impact of variations in beam hardening and scatter. This proximity ensures that the calibration accurately reflects the imaging conditions specific to the region of interest [8] [10]; 2) Consistent composition across subjects: Both aortic blood and psoas muscle exhibit relatively stable and predictable compositions, resulting in consistent HU values across individuals and imaging protocols. This consistency makes them suitable for calibration without the need for external phantoms [6] [8]; 3) Established use in literature: Previous studies have demonstrated the feasibility of using internal soft tissues for phantom-less calibration. For instance, Lee et al. [10] reported that aortic blood and psoas muscle provide reliable internal calibration metrics for BMD assessment, achieving high precision and reproducibility. Similarly, Therkildsen et al. highlighted the potential of internal calibration in routine clinical practice [9]; 4) Practical feasibility in routine clinical imaging: These tissues are routinely visualized in thoracic and abdominal CT scans, allowing opportunistic bone health assessment without additional imaging or radiation exposure [7] [12].
Current study focuses on primary osteoporosis, a condition where bone loss occurs due to the natural aging process. To accurately study primary osteoporosis, it is essential to select subjects who are as healthy as possible and free from conditions that might lead to secondary osteoporosis. Secondary osteoporosis is caused by various underlying diseases or conditions, which introduce significant changes in BMD independent of the aging process. The following explains the rationale for excluding these conditions: 1) Chronic Kidney Disease (CKD): CKD is a common cause of secondary osteoporosis due to renal osteodystrophy. It involves abnormalities in bone turnover, mineralization, and volume resulting from impaired kidney function. CKD often leads to secondary hyperparathyroidism, altered vitamin D metabolism, and calcium-phosphate imbalances, all of which severely affect bone health and BMD [18]. Including patients with CKD could introduce pathological bone loss unrelated to primary osteoporosis; 2) Vitamin D Deficiency: Vitamin D plays a pivotal role in calcium absorption and bone mineralization. Severe deficiency leads to conditions like osteomalacia in adults, characterized by reduced BMD. Such variations, driven by a treatable deficiency rather than the natural aging process, could skew the study results [19]; 3) Long-Term Steroid Use: Glucocorticoids are well-known for inducing osteoporosis by decreasing bone formation and increasing bone resorption. Including patients on long-term steroids could confound BMD measurements, as the bone loss in such cases is pharmacologically induced rather than age-related; 4) Endocrine Disorders: Both hyperparathyroidism and hyperthyroidism significantly affect BMD and introduce systemic influences that could compromise the study’s focus on primary osteoporosis; 5) Cirrhosis: Chronic liver disease contributes to hepatic osteodystrophy, which involves reduced bone formation and decreased BMD due to malnutrition, hypogonadism, and vitamin D metabolism abnormalities. These factors would confound the assessment of age-related bone loss; 6) Acute Leukemia and Diffuse Spinal Metastases: These conditions directly impact bone marrow and structural integrity, causing localized or diffuse changes in BMD unrelated to physiological aging. By excluding these diseases, which are common causes of secondary osteoporosis, the study ensures that the observed BMD variations are primarily due to the aging process. This approach helps isolate the natural progression of bone loss in healthy individuals, thereby improving the validity and reliability of the findings.
Despite the strong results validating the internal calibration techniques, it is still uncertain if these reference values are applicable across different subjects and across imaging protocols for opportunistic screening. Therefore, a broader approach for bone mineral density calibration is necessary to reduce the assumptions made by relying on these ground truth reference values.
5. Limitations
There are a few limitations in this study. First, we did not have lumbar BMD measurements by DXA, which is the consensus clinical standard to diagnose osteoporosis. In addition, the definition of the ROIs in spinal measurements is somewhat arbitrary. Second, our conclusions are limited to non-contrast scans. If intravenous contrast is used for chest and abdominal CT, some notable measurement bias in the thoracic and lumbar BMD could be observed [20]. Further investigations are needed to establish guidance for reliable use of thoracic and lumbar BMD in cases where CT scans include contrast agents.
6. Conclusion
A significant correlation was observed between thoracic and lumbar BMD values. The correlation coefficient between thoracic and lumbar BMD was about 0.953. Both thoracic and lumbar BMD values decrease with the subjects’ age. In addition, if the threshold lumbar BMD value of 123 mg/cm3 is used as a criterion for osteoporosis diagnosis, the area under the ROC curve of thoracic BMD values is 0.969. These findings support the use of CT as a valuable tool in screening for osteoporosis based on thoracic and lumbar BMDs.