Low-dose MDCT: evaluation of the impact of systematic tube current reduction and sparse sampling on quantitative paraspinal muscle assessment
Introduction
Cachexia and sarcopenia are katabolic diseases which have great impact on muscle structure and composition (1-3). In course of wasting diseases the assessment and monitoring of overall muscle composition plays a pivotal role for initiation and adaption of further treatment steps ranging from nutritional supplements to adjustment of chemotherapy (3). In all diseases affecting the musculoskeletal system monitoring the composition of large or representative muscle groups is a matter of special interest to track the individual disease progression (4). In literature analogous fatty infiltration patterns of the paravertebral musculature with the thigh have been described before (5). As the paravertebral musculature seems to underlie the same katabolic patterns in muscle decay as the large muscle groups, it can be considered suitable for overall muscle status assessment. Furthermore, the psoas muscle has been identified as an important target muscle for tumor related sarcopenia/cachexia (6). Weerink et al. reported that the presence of low psoas mass prior to surgery, as an indicator for sarcopenia, and a strong predictor for the development of postoperative complications (6). In further studies the psoas muscle was used as target muscle for prediction of complications and survival in patients with gastrectomy, cystectomy, and liver surgery (7-9).
There is a plentitude of imaging methods in clinical use from dual energy X-ray absorptiometry (DXA) and multi-detector computed tomography (MDCT) to single-voxel proton magnetic resonance spectroscopy (MRS) and chemical shift encoding-based water-fat magnetic resonance imaging (MRI) (10-15). Although, MRS and MRI enable the investigator to extract surrogate parameters of the muscle composition like the fat fraction (FF) and even for the identification of the chemical structure of fatty acids and their magnitude with high congruence to histology, they require advanced knowledge during image acquisition and are rather high priced (15). MDCT is broadly available in even smaller hospitals’ allocation and used in monitoring of oncological patients, liver fat and subcutaneous-/ visceral fat assessment (16,17). The number of MDCT scans and corresponding radiation dose a subject received, was shown to be directly correlated with estimated individual cancer risk (18-20). One possibility to tackle this modality specific issue is to reduce the radiation dose received by each subject undergoing an MDCT scan.
Radiation dose reduction can generally be achieved by different approaches including lowering of tube current or by using sparse sampling, which is a technique that allows for acquisition of fewer projections in MDCT scans (21,22). However, the extent of tube current attenuation and decrease of projection numbers is limited with regard to preservation of scan quality (23,24). Sparse sampling is a technologically advanced method in which radiation exposure is lowered as a consequence of decreased projection numbers, whereas the needed energy delivery is stable (23).
In this study the limits of virtual tube current reduction of MDCT were compared to the sparse sampling technique and the validity of the acquired images for the assessment of paraspinal muscle composition was investigated.
Methods
Subjects
Eighteen patients were included in this retrospective study (6 males, age: 71.5±15.9 years, and 12 females, age: 71.0±8.9 years). The included patients underwent non-contrast MDCT imaging of the spine at our department due to suspected degenerative spine diseases or trauma. Our hospital’s picture archiving and communication system (PACS, IDS7; Sectra AB, Linköping, Sweden) was used for patient identification.
Identification criteria were: (I) existing MDCT scan of the lumbar vertebral spine, (II) no medical history of metabolic or muscular diseases, (III) no relevant, image quality reducing artifacts, (IV) no previous surgery with instrumentation at the spine, and (5) no presence of any implants in the field of view.
This study was approved by the local institutional review board (62/18S) and was conducted in accordance with the Declaration of Helsinki. Written and oral informed consent for future anonymized data processing was given by each participant when signing the informed consent before the CT scan.
MDCT
Image acquisition was performed in supine position using a 128-slice MDCT scanner (Ingenuity Core 128; Philips Healthcare, Cleveland, OH). An initial scout scan was used for planning purposes. Then a helical scan was performed with implicit tube current modulation. Scan parameters are shown in Table 1.
Full table
Virtual tube current reduction
Tube current reduction was performed virtually. For this purpose, a simulation algorithm based on raw imaging data was used, which allowed for generating MDCT scans with virtually lowered tube currents (25-27). The original scan with full dose without dose reduction or sparse sampling was defined as D100P100. The corresponding simulations had 50% (D50P100), 10% (D10P100), 5% (D5P100), and 3% (D3P100) of the original X-ray tube current.
Sparse sampling and statistical iterative reconstruction (SIR)
Sparse sampling was simulated by reading only a reduced amount of projection angles and by deleting the remaining projections in the sonogram (27,28). Sparse-sampled images were generated as if MDCT was performed with only 50% (D100P50), 10% (D100P10), 5% (D100P5), and 3% (D100P3) of the original projection data.
SIR
For both virtual tube current reduction as well as sparse sampling the same in-house developed SIR algorithm that was based on ordered-subset separable paraboloidal surrogate combining a momentum accelerating approach was used (29,30). In detail, a Gaussian noise model was applied and the likelihood term for SIR was computed with log-converted projection data. To enhance convergence and to further depress image noise while achieving adequate bone/soft tissue contrast, a regularization term based on a Huber penalty was applied. The calibration data served for calculating linear attenuation coefficients of resulting imaging data, which were then translated to Hounsfield units by using air and water information.
Segmentation of the psoas muscle and the erector spinae and muscle density (MD) measurements
The psoas muscle and the erector spinae muscles were segmented bilaterally at the level of L4 over a craniocaudal distance of 10 mm in the axial reformations of the original dose images (Figure 1). Segmentation was performed by using the free open-source software Medical Imaging Interaction Toolkit (MITK, developed by the Division of Medical and Biological Informatics, German Cancer Research Center, Heidelberg, Germany) by a radiologist with 4 years of experience. Segmentations were superimposed on all low-dose scans. MD was extracted in Hounsfield Units (HU) of each muscle in each scan. Volume weighted mean MD of the right and left psoas and erector spinae muscle was calculated, respectively. A second reader performed segmentation of all included subjects to evaluate the inter-rater reliability of the results.
Signal-to-noise ratio (SNR) calculation
The SNR from each virtually reduced and each sparse sampling series was calculated as follows: SNR = mean in HU/ standard deviation in HU.
Statistical analysis
For the statistical analyses SPSS (version 22.0; IBM SPSS Statistics for Windows, Armonk, NY, USA) was used. Statistical significance was considered at p < 0.05 (two-sided) in all conducted tests.
The Kolmogorov–Smirnov test indicated normally distributed data.
MD values derived from virtually lowered tube current and sparse sampling were compared with each other and with those derived from original dose images as gold standard by using paired t-tests. To assess the reliability of the acquired results, intraclass correlation coefficient (ICC) was calculated.
Results
Mean MD was 6.41±13.55 HU for the erector spinae muscles derived from original dose images (D100P100). Sparse sampling derived mean MD showed no significant changes (P=0.57) down to 5% of the original projections (D100P5). For the psoas muscle, mean MD amounted 39.86±6.14 HU. No significant changes in mean MD of the psoas muscle were observed down to D100P5 (P=0.22) (Table 2 and Figure 2). Representative images of virtually lowered tube current and sparse sampling are shown in Figure 2.
Full table
The virtual tube lowering showed no significant (P>0.05) changes for the mean MD of the erector spinae at 5% of the original dose D5P100 (P=0.19) as compared to D100P100. However, mean MD derived from D50P100, D10P100 and D3P100 were significantly different as compared to D100P100 (P<0.01). In the psoas muscle, all virtually reduced tube current series showed significantly different mean MD as compared to the original dose.
The comparison of corresponding sparse sampling and virtually tube reduced series (D100P50/D50P100, D100P10/D10P100, D100P5/D5P100 and D100P3/D3P100) is shown in Table 3. The 50% and 10% reduced series of the psoas and the 5% reduced series of the erector spinae MD did not differ significantly from each other (P>0.05) (Figure 3).
Full table
The Signal-to-noise ratio (SNR) for each virtually tube reduced series is shown in Table 4. In the erector spinae the SNR did not reveal significant differences for the sparse sampling series down to 5% (D100P5) of the original dose (P=0.08). A high ICC was given between the two readers in both muscle compartments (ICC =0.98, P<0.01).
Full table
Discussion
In this study, we compared different approaches to effectively minimize radiation exposure in MDCT scans using virtual tube current reduction and sparse sampling. With sparse sampling, non-significant changes in measured MD could be detected with decreased projection numbers down to 5% of the original data. In contrast, MD extracted from virtually reduced tube current images was considerably affected already at 50% of the original dose.
In clinical settings, sparse sampling could be an option in the future with regard to body composition assessment by measuring the MD of representative muscle groups with drastically minimized radiation exposure. Sollmann et al. already showed the feasibility of sparse sampling not only for BMD determination, but also in fracture diagnostic and imaging of degenerative changes of the spine (25,26,31). However, available MDCT scanners are not yet capable of applying sparse sampling, because of the continuous X-ray delivery during the scanning procedure, which is a shortcoming that might be tackled by future scanner generations.
The presented study showed not only a future perspective of sparse sampling, but also pointed to the possibility of applying much lower tube currents in combination with iterative image reconstruction to assess katabolic muscle changes in course of sarcopenia and cachexia. The presented findings can be directly transferred to clinical practice as a lowering of tube current can be achieved without technical effort. An aspect which should be kept in mind considering the application of tube current reduction is the increase of image noise and artifacts accompanying the scan (32). The reason for overall better results for sparse-sampled data most likely relates to the inherently different method of sparse sampling compared to tube current reduction. Specifically, while the approach of tube current reduction leads to considerable increases in image noise, that can negatively affect diagnostic use, sparse sampling implies maintained energy delivery for the individual projection image, but lowered overall radiation exposure due to the decrease in total projection numbers (33,34). Consequently, preserved image quality for the individual projection can be achieved whilst circumventing the influence of electronic readout noise, leading to largely preserved structural image information.
In the past, multiple approaches for dose reduction in MDCT-based body composition imaging were performed. Amongst others, Chang et al. showed the possibility of tube current reduction up to 50% of original doses in liver imaging with a preservation of image quality and contrast-to-noise ratio and stable quantitative measures (35,36). Furthermore, Yamada et al. proved the feasibility of lowering radiation doses to 30% of the original dose in imaging visceral and subcutaneous adipose tissue, which constitute the largest share of human fat deposition localizations (37). Taking the results of these groups into consideration, Elhamiasl and Nuyts described a phantom-based technical approach to identifying the lowest exposure scan which still produces sufficient information for clinical applications analyzing the noise power spectrum of each dose-reduced image (38). In the proposed model presented in this study several parameters influencing image quality (like electrical noise, crosstalk between detector channels and the beam hardening effect amongst others) have been taken into account to determine the lowest possible radiation dose that still produces sufficient information for proton therapy treatment planning (38).
To best of our knowledge, our presented results are the first that affirm the insights these prior studies gave into dose reduction in body composition imaging on basis of muscle density measurements. Within the limits of disease progression monitoring in entities like cachexia and sarcopenia the revealed potential but also the limits of sparse sampling based projection number reduction could be useful for clinical purposes in the future.
With that said, there are certain limitations that come along with the presented study. First, the diagnostic value of the dose or projection reduced MDCT scans with regard to other pathologies is mitigated. Like Jensen et al. showed the detection and evaluation of liver metastases is compromised in modest radiation dose reduction (39). Second, the presented cohort is rather small comprising no patients suffering from cachexia or from an oncological disease, which is a shortcoming keeping in mind the potential future application settings. Last, sparse sampling is not applicable in current clinical settings due to technical scanner properties, although it showed promising results and will pose an option for upcoming MDCT scanner generations.
Conclusions
The reported insights into virtual tube current reduction and sparse sampling based muscle density evaluation of the psoas and erector spinae muscles suggest the possibility of considerable dose reduction in MDCT scans with image quality preservation. Additionally, the results of the current study showed the potential of sparse sampling for future clinical applications in body composition imaging.
Acknowledgments
Funding: The present work was supported by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) - Project 432290010 (to TB, PBN and JSK).
Footnote
Conflicts of Interest: All authors have completed the ICMJE uniform disclosure form (available at http://dx.doi.org/10.21037/qims-20-1220). The authors have no conflicts of interest to declare.
Ethical Statement: This study was approved by the Institutional Review Board and was conducted in accordance with the Declaration of Helsinki. Written informed consent was obtained from all study participants.
Open Access Statement: This is an Open Access article distributed in accordance with the Creative Commons Attribution-NonCommercial-NoDerivs 4.0 International License (CC BY-NC-ND 4.0), which permits the non-commercial replication and distribution of the article with the strict proviso that no changes or edits are made and the original work is properly cited (including links to both the formal publication through the relevant DOI and the license). See: https://creativecommons.org/licenses/by-nc-nd/4.0/.
References
- Batsis JA, Villareal DT. Sarcopenic obesity in older adults: aetiology, epidemiology and treatment strategies. Nat Rev Endocrinol 2018;14:513-37. [Crossref] [PubMed]
- Filippin LI, Teixeira VN, da Silva MP, Miraglia F, da Silva FS. Sarcopenia: a predictor of mortality and the need for early diagnosis and intervention. Aging Clin Exp Res 2015;27:249-54. [Crossref] [PubMed]
- Ni J, Zhang L. Cancer Cachexia: Definition, Staging, and Emerging Treatments. Cancer Manag Res 2020;12:5597-605. [Crossref] [PubMed]
- Spitali P, Hettne K, Tsonaka R, Charrout M, van den Bergen J, Koeks Z, Kan HE, Hooijmans MT, Roos A, Straub V, Muntoni F, Al-Khalili-Szigyarto C, Koel-Simmelink MJA, Teunissen CE, Lochmüller H, Niks EH, Aartsma-Rus A. Tracking disease progression non-invasively in Duchenne and Becker muscular dystrophies. J Cachexia Sarcopenia Muscle 2018;9:715-26. [Crossref] [PubMed]
- Burian E, Inhuber S, Schlaeger S, Dieckmeyer M, Klupp E, Franz D, Weidlich D, Sollmann N, Löffler M, Schwirtz A, Rummeny EJ, Zimmer C, Kirschke JS, Karampinos DC, Baum T. Association of thigh and paraspinal muscle composition in young adults using chemical shift encoding-based water-fat MRI. Quant Imaging Med Surg 2020;10:128-36. [Crossref] [PubMed]
- Weerink LBM, van der Hoorn A, van Leeuwen BL, de Bock GH. Low skeletal muscle mass and postoperative morbidity in surgical oncology: a systematic review and meta-analysis. J Cachexia Sarcopenia Muscle 2020;11:636-49. [Crossref] [PubMed]
- Cornet M, Lim C, Salloum C, Lazzati A, Compagnon P, Pascal G, Azoulay D. Prognostic value of sarcopenia in liver surgery. J Visc Surg 2015;152:297-304. [Crossref] [PubMed]
- Yaguchi Y, Kumata Y, Horikawa M, Kiyokawa T, Iinuma H, Inaba T, Fukushima R. Clinical Significance of Area of Psoas Major Muscle on Computed Tomography after Gastrectomy in Gastric Cancer Patients. Ann Nutr Metab 2017;71:145-9. [Crossref] [PubMed]
- Ahmadi H, Montie JE, Weizer AZ, Morgan T, Montgomery JS, Lee CT. Patient Psoas Muscle Mass as a Predictor of Complications and Survival After Radical Cystectomy. Curr Urol Rep 2015;16:79. [Crossref] [PubMed]
- Löffler MT, Sollmann N, Mei K, Valentinitsch A, Noël PB, Kirschke JS, Baum T. X-ray-based quantitative osteoporosis imaging at the spine. Osteoporos Int 2020;31:233-50. [Crossref] [PubMed]
- Baum T, Cordes C, Dieckmeyer M, Ruschke S, Franz D, Hauner H, Kirschke JS, Karampinos DC. MR-based assessment of body fat distribution and characteristics. Eur J Radiol 2016;85:1512-8. [Crossref] [PubMed]
- Burian E, Syväri J, Holzapfel C, Drabsch T, Kirschke JS, Rummeny EJ, Zimmer C, Hauner H, Karampinos DC, Baum T, Franz D. Gender- and Age-Related Changes in Trunk Muscle Composition Using Chemical Shift Encoding-Based Water-Fat MRI. Nutrients 2018;10:1972. [Crossref] [PubMed]
- Sollmann N, Dieckmeyer M, Schlaeger S, Rohrmeier A, Syvaeri J, Diefenbach MN, Weidlich D, Ruschke S, Klupp E, Franz D, Rummeny EJ, Zimmer C, Kirschke JS, Karampinos DC, Baum T. Associations Between Lumbar Vertebral Bone Marrow and Paraspinal Muscle Fat Compositions-An Investigation by Chemical Shift Encoding-Based Water-Fat MRI. Front Endocrinol (Lausanne) 2018;9:563. [Crossref] [PubMed]
- Franz D, Weidlich D, Freitag F, Holzapfel C, Drabsch T, Baum T, Eggers H, Witte A, Rummeny EJ, Hauner H, Karampinos DC. Association of proton density fat fraction in adipose tissue with imaging-based and anthropometric obesity markers in adults. Int J Obes (Lond) 2018;42:175-82. [Crossref] [PubMed]
- Karampinos DC, Ruschke S, Dieckmeyer M, Diefenbach M, Franz D, Gersing AS, Krug R, Baum T. Quantitative MRI and spectroscopy of bone marrow. Quantitative MRI and spectroscopy of bone marrow. J Magn Reson Imaging 2018;47:332-53. [Crossref] [PubMed]
- Graffy PM, Sandfort V, Summers RM, Pickhardt PJ. Automated Liver Fat Quantification at Nonenhanced Abdominal CT for Population-based Steatosis Assessment. Radiology 2019;293:334-42. [Crossref] [PubMed]
- Gao Y, Wang YC, Lu CQ, Zeng C, Chang D, Ju S. Correlations between the abdominal fat-related parameters and severity of coronary artery disease assessed by computed tomography. Quant Imaging Med Surg 2018;8:579-87. [Crossref] [PubMed]
- Richards PJ, George J, Metelko M, Brown M. Spine computed tomography doses and cancer induction. Spine (Phila Pa 1976) 2010;35:430-3. [Crossref] [PubMed]
- Richards PJ, George J. Diagnostic CT radiation and cancer induction. Skeletal Radiol 2010;39:421-4. [Crossref] [PubMed]
- Smith-Bindman R, Lipson J, Marcus R, Kim KP, Mahesh M, Gould R, Berrington de González A, Miglioretti DL. Radiation dose associated with common computed tomography examinations and the associated lifetime attributable risk of cancer. Arch Intern Med 2009;169:2078-86. [Crossref] [PubMed]
- Omoumi P, Verdun FR, Ben Salah Y, Vande Berg BC, Lecouvet FE, Malghem J, Ott JG, Meuli R, Becce F. Low-dose multidetector computed tomography of the cervical spine: optimization of iterative reconstruction strength levels. Acta Radiol 2014;55:335-44. [Crossref] [PubMed]
- Becce F, Ben Salah Y, Verdun FR, Vande Berg BC, Lecouvet FE, Meuli R, Omoumi P. Computed tomography of the cervical spine: comparison of image quality between a standard-dose and a low-dose protocol using filtered back-projection and iterative reconstruction. Skeletal Radiol 2013;42:937-45. [Crossref] [PubMed]
- Abbas S, Lee T, Shin S, Lee R, Cho S. Effects of sparse sampling schemes on image quality in low-dose CT. Med Phys 2013;40:111915 [Crossref] [PubMed]
- Willemink MJ, Noel PB. The evolution of image reconstruction for CT-from filtered back projection to artificial intelligence. Eur Radiol 2019;29:2185-95. [Crossref] [PubMed]
- Sollmann N, Mei K, Hedderich DM, Maegerlein C, Kopp FK, Löffler MT, Zimmer C, Rummeny EJ, Kirschke JS, Baum T, Noël PB. Multi-detector CT imaging: impact of virtual tube current reduction and sparse sampling on detection of vertebral fractures. Eur Radiol 2019;29:3606-16. [Crossref] [PubMed]
- Sollmann N, Mei K, Schwaiger BJ, Gersing AS, Kopp FK, Bippus R, Maegerlein C, Zimmer C, Rummeny EJ, Kirschke JS, Noël PB, Baum T. Effects of virtual tube current reduction and sparse sampling on MDCT-based femoral BMD measurements. Osteoporos Int 2018;29:2685-92. [Crossref] [PubMed]
- Muenzel D, Koehler T, Brown K, Zabić S, Fingerle AA, Waldt S, Bendik E, Zahel T, Schneider A, Dobritz M, Rummeny EJ, Noël PB. Validation of a low dose simulation technique for computed tomography images. PLoS One 2014;9:e107843 [Crossref] [PubMed]
- Zabić S, Wang Q, Morton T, Brown KM. A low dose simulation tool for CT systems with energy integrating detectors. Med Phys 2013;40:031102 [Crossref] [PubMed]
- Rayudu NM, Anitha DP, Mei K, Zoffl F, Kopp FK, Sollmann N, Löffler MT, Kirschke JS, Noël PB, Subburaj K, Baum T. Low-dose and sparse sampling MDCT-based femoral bone strength prediction using finite element analysis. Arch Osteoporos 2020;15:17. [Crossref] [PubMed]
- Kim D, Ramani S, Fessler JA. Combining ordered subsets and momentum for accelerated X-ray CT image reconstruction. IEEE Trans Med Imaging 2015;34:167-78. [Crossref] [PubMed]
- Sollmann N, Mei K, Riederer I, Probst M, Löffler MT, Kirschke JS, Noël PB, Baum T. Low-dose MDCT: evaluation of the impact of systematic tube current reduction and sparse sampling on the detection of degenerative spine diseases. Eur Radiol 2021;31:2590-600. [Crossref] [PubMed]
- Duan X, Wang J, Leng S, Schmidt B, Allmendinger T, Grant K, Flohr T, McCollough CH. Electronic noise in CT detectors: Impact on image noise and artifacts. AJR Am J Roentgenol 2013;201:W626-32 [Crossref] [PubMed]
- Willemink MJ, Leiner T, de Jong PA, de Heer LM, Nievelstein RA, Schilham AM, Budde RP. Iterative reconstruction techniques for computed tomography part 2: initial results in dose reduction and image quality. Eur Radiol 2013;23:1632-42. [Crossref] [PubMed]
- Noël PB, Renger B, Fiebich M, Münzel D, Fingerle AA, Rummeny EJ, Dobritz M. Does iterative reconstruction lower CT radiation dose: evaluation of 15,000 examinations. PLoS One 2013;8:e81141 [Crossref] [PubMed]
- Zhang X, Chen J, Yu N, Ren Z, Tian Q, Tian X, He T, Guo C. Improving image quality with model-based iterative reconstruction at quarter of nominal dose in upper abdominal CT. Br J Radiol 2019;92:20180137 [PubMed]
- Chang W, Lee JM, Lee K, Yoon JH, Yu MH, Han JK, Choi BI. Assessment of a model-based, iterative reconstruction algorithm (MBIR) regarding image quality and dose reduction in liver computed tomography. Invest Radiol 2013;48:598-606. [Crossref] [PubMed]
- Yamada Y, Jinzaki M, Niijima Y, Hashimoto M, Yamada M, Abe T, Kuribayashi S. CT Dose Reduction for Visceral Adipose Tissue Measurement: Effects of Model-Based and Adaptive Statistical Iterative Reconstructions and Filtered Back Projection. AJR Am J Roentgenol 2015;204:W677-83 [Crossref] [PubMed]
- Elhamiasl M, Nuyts J. Low-dose x-ray CT simulation from an available higher-dose scan. Phys Med Biol 2020;65:135010 [Crossref] [PubMed]
- Jensen CT, Wagner-Bartak NA, Vu LN, Liu X, Raval B, Martinez D, Wei W, Cheng Y, Samei E, Gupta S. Detection of Colorectal Hepatic Metastases Is Superior at Standard Radiation Dose CT versus Reduced Dose CT. Radiology 2019;290:400-9. [Crossref] [PubMed]