Assessment of Myosteatosis on Computed Tomography by Automatic Generation of a Muscle Quality Map Using a Web-Based Toolkit: Feasibility Study

Introduction

Measurement of muscle is one of the fastest growing research areas in medicine, as quantity and quality of muscle in individuals are reportedly associated with morbidity and mortality in various diseases [1,2]. Muscle mass is frequently assessed using medical imaging such as computed tomography (CT) or magnetic resonance imaging (MRI) [3]. However, an issue of discrepancy between the muscle mass measured on imaging and muscle functions, such as strength and mobility, has been raised. Therefore, the assessment of muscle quality by imaging is gaining a focus in the research and clinical diagnoses of relevant diseases. Specifically, the CT density or attenuation can be quantified using a standardized unit (ie, Hounsfield unit [HU]), enabling a standardized muscle quality evaluation on CT.

Recent studies demonstrated that muscle quality is associated with fatty degeneration or fatty infiltration of the muscle (ie, myosteatosis), which may be associated with decreased muscle function and increased disability [4]. Based on the imaging characteristics in CT, fat can be stored as follows: (1) in the intermuscular adipose tissue, which is observed as gross fat between muscle groups; (2) in the intramuscular adipose tissue as extramyocellular lipids, which are denoted as gross fat tissues between muscle fibers in the same muscle group; and (3) in the intramyocellular lipid droplets, which are not visually demonstrated as fat but are indicated by decreased muscle density on CT [5]. It is known that the intermuscular and intramuscular adipose tissues (IMAT) can be depicted as areas with gross fat density of –190 HU to –30 HU, whereas the intramyocellular lipid can be reflected as low-attenuation muscle area (LAMA) with –29 HU to +29 HU within the muscle [6,7]. The healthy muscle without myosteatosis is observed as normal-attenuation muscle area (NAMA) with +30 HU to +150 HU [6].

The muscle density at the L3 vertebral level has been used for muscle quality evaluation on abdominal CT in majority of prior studies because it easily measures muscle density [1,8]. Recent advancements in technology evaluate myosteatosis based on the distribution and amount of IMAT and LAMA. In this paper, we evaluate the feasibility of fully automated quantitative measurement of the skeletal muscle on CT using a web-based toolkit for assessing normal-attenuation muscle and myosteatosis.

Methods

This study, which uses retrospective data, was approved by the institutional review board of Asan Medical Center, and a waiver for written informed consent was obtained.

Generation of a Muscle Quality Map

Among the abdominal CT images, a single axial image at the inferior endplate level of the L3 vertebra at the portal venous phase was selected for image assessment [9]. The selected CT image was saved as a digital imaging and communications in medicine (DICOM) file and was uploaded to our web-based toolkit in a drag-and-drop manner. First, we performed automatic segmentation of the abdominal compartments using a predeveloped deep learning model based on a fully convolutional network [10]. This model was reported to segment abdominal body compartments into the total abdominal muscle area (TAMA), subcutaneous fat area, and visceral fat area with a dice similarity coefficient of 0.97 and mean cross-sectional area error of 2.26% [10].

Second, we performed postprocessing to divide the TAMA into 3 muscle components using the pixel-wise measurement of CT density with the following range of HU for each component [6,7]: NAMA with +30 HU to +150 HU, LAMA with –29 HU to +29 HU, and IMAT with –190 HU to –30 HU. Subsequently, a muscle quality map was generated by combining NAMA, LAMA, and IMAT, which were displayed in different colors on the web-based toolkit (Figure 1A), enabling the quantitative assessment of muscle properties based on histogram analysis. For example, the participants with similar TAMA values can be differentiated into participants with healthy muscles (ie, mostly composed of NAMA) and those with fatty degenerated muscles (ie, comprising several LAMA and IMAT). The technique has been packaged and can be accessed at an iAID Viewer repository [11]. Figure 1 shows the process and examples of muscle quality map generation using the web-based iAID toolkit. A CT image is uploaded to the toolkit in a drag-and-drop manner. Muscle quality map and its measurement values are displayed upon clicking on the “L3 Area Measure” button (Figure 1A). Muscle quality maps and a histogram for a 28-year-old man with higher quality muscle are shown in Figure 1B; and muscle quality maps and a histogram for a 72-year-old man with lower quality muscle are shown in Figure 1C.

‎View this figure

The quality of the generated muscle quality map was evaluated by a board-certified abdominal radiologist by comparing the original CT image and the muscle quality map image. If the segmentation quality was not acceptable, manual adjustment was performed using a stand-alone software.

Feasibility Study in Healthy Participants

From the electronic database of Asan Medical Center, we retrospectively identified 3928 abdominal CT scans from healthy participants who had no disease or history of previous treatments and had undergone CT for medical checkup or workup for liver donation during the period from January 2008 to June 2016. We used the block randomization technique to randomly select 240 healthy participants (135 men and 105 women) belonging to 6 different age groups, from 20 to 79 years, with 10-year intervals (40 participants per age group).

CT was performed using a 16-channel (LightSpeed 16, GE Healthcare; Somatom Sensation 16, Siemens Medical Solution), a 64-channel (Discovery CT 750 HD, GE Healthcare; Somatom Definition AS, Siemens Medical Solution), or a 128-channel (Somatom Definition Flash, Siemens Medical Solution) CT scanner with the following parameters: 120 kVp, 200–220 mAs (maximum tube current with automated dose modulation), 1.5-5 mm section thickness and intervals, and a pitch of 0.6-1. Contrast agents were administered at a rate of 3-4 mL/s, and CT images (regardless of the CT protocols) were obtained, including the portal venous phase (65-72 seconds after contrast agent injection) in the craniocaudal direction.

Statistical Analysis

To eliminate effects by different sexes, participants within each subgroup were stratified by sex. The mean values of muscle components (TAMA, NAMA, LAMA, and IMAT) were compared among the age groups using 1-way analysis of variance. A P<.05 was considered statistically significant. Correlation between the mean density of TAMA or skeletal muscle area (SMA; ie, NAMA + LAMA) and muscle components related to the adipose tissue (ie, LAMA, IMAT, or myosteatosis area represented by LAMA + IMAT) were calculated using Pearson correlation (absolute magnitude of Pearson correlation coefficient r=0-0.1, negligible correlation; r=0.1-0.4, weak correlation; r=0.4-0.7, moderate correlation; r=0.7-0.9, strong correlation; and r=0.9-1, very strong correlation) [12]. In addition, to investigate the correlation between the adipose tissue within and outside the muscle compartment, the correlation between the aforementioned muscle components related to the adipose tissue—LAMA, IMAT, and myosteatosis area—and the visceral and subcutaneous fat area was also analyzed using Pearson correlation. All statistical analyses were performed using the SPSS, version 21 (IBM Corp).

Results

The characteristics of the 240 participants are presented in Table 1. In male participants, height and weight were significantly different across age groups (P<.001). In female participants, height was significantly different across age groups (P<.001). Body mass index was not significantly different both in male (P=.11) and female participants (P=.13).

Figure 2 illustrates the mean areas of muscle components among different age groups of male and female participants.

The detailed information is presented in Table 2. TAMA showed a significant difference only in male participants (male, P<.001; female, P=.88). The age group of 30s had the highest TAMA (mean 187.9 cm², SD 29.6 cm²) with a gradual decrease in area with aging in male participants. The TAMA in female participants did not differ significantly across the age groups. In contrast, NAMA, LAMA, and IMAT were significantly different among the age groups both in male and female participants (P<.001). In male participants, NAMA showed a gradual decrease with age (peak age: 20 years; mean 153.4 cm², SD 20.0 cm²). In female participants, there was a peak in the 30s age group (mean 92.6 cm², SD 15.1 cm²), followed by a gradual decrease in area with aging. The myosteatosis area showed a gradual increase with aging in female participants.

‎View this figure

All muscle components related to the adipose tissue (ie, LAMA, IMAT, and myosteatosis area) had a moderate-to-strong negative correlation with the mean density of TAMA (r=–0.629 to –0.884) and SMA (r=–0.647 to –0.898) (Table 3). They had a moderate positive correlation (r=0.474 to 0.686) with visceral fat area and weak-to-moderate positive correlation (r= 0.274 to 0.459) with subcutaneous fat area in both male and female participants (Table 3). In particular, LAMA showed higher correlation with visceral and subcutaneous fat compartments than that shown by IMAT. All the correlations in Table 3 are statistically significant (P<.05).

Discussion

Principal Results

We developed a web-based toolkit to generate a pixel-based automatic categorization of muscle components (ie, muscle quality map) within the generated segmented muscle compartment using predeveloped deep learning models. These muscle quality maps can illustrate spatial distribution of fat infiltration and provide insights into the muscle quality in individuals. In our study, NAMA gradually decreased with age; and LAMA and IMAT gradually increased with age. These results indicate that fat infiltration or fatty degeneration increases with aging. Therefore, NAMA might be an effective imaging biomarker for evaluation of muscle in individuals. Additionally, LAMA and IMAT might be used as biomarkers of myosteatosis and related diseases, such as metabolic syndrome [1,5].

Comparison With Previous Work

Sarcopenia, which is defined as the loss of muscle mass or function, is associated with increased morbidity and mortality in various diseases [1,2]. Recently, the revised consensus of the European Working Group on Sarcopenia in Older People emphasized muscle strength and quality as a key characteristic of sarcopenia in the definition and diagnostic criteria of sarcopenia [8]. Accordingly, both qualitative and quantitative evaluations of muscle are highly recommended in sarcopenia diagnosis. Among the diagnostic tests for muscle quantity measurements, muscle quality evaluation can be performed only using cross-sectional imaging such as CT and MRI. The dual-energy x-ray absorptiometry and bioelectrical impedance analysis cannot differentiate healthy muscles from fatty degenerated muscles. Currently, the CT density or attenuation is well-calibrated and standardized (ie, HU of zero for water) across the CT acquisition protocols and imaging machines, whereas the signal intensity of MRI might differ between protocols and machines. Therefore, the muscle quality map on CT might be the optimal option to evaluate myosteatosis.

There are several software programs (both open source and licensed) that analyze muscle and fat composition in CT [13]. However, a software providing muscle quality maps at the L3 vertebral level is not yet available. We developed our own software for performing a fully automated segmentation of body compartments and generating a muscle quality map at the L3 vertebral level [11]. This software is publicly available for academic purposes. Considering the increasing recognition of myosteatosis and sarcopenia as prognostic biomarkers for various diseases in older patients, measuring muscle quality and muscle mass easily is highly recommended for researchers and clinicians.

Limitations

Our study has some limitations. First, we tested the muscle quality map in a relatively small number of participants (N=240). We randomly assigned 40 participants into each age group; thus, our study is an experimental study rather than an epidemiologic study and may not reflect the real-world trends of muscle quality with aging. There was an imbalance between male and female participants, which may preclude the real-world data. In addition, we could not provide any criteria or cutoff value to diagnose myosteatosis. Therefore, an epidemiological study using large populations might be required. Second, the segmented areas of each muscle component might be affected by the different CT protocols and scanners. However, this is beyond the scope of the current study and needs to be evaluated in a separate study. Third, the general applicability of our toolkit should be validated in an external dataset, particularly because the predeveloped segmentation model was trained at our institution and could thus lead to overfitting.

Conclusion

In conclusion, the automated web-based toolkit is feasible and enables a quantitative CT assessment of normal-attenuation muscle and myosteatosis, which can be effective quantitative biomarkers for evaluating structural and functional changes brought on by aging in the skeletal muscle.