Introduction

Sarcopenia, first introduced by Rosenberg in 19841, refers to the age-related decline in skeletal muscle mass and quality. In 2016, sarcopenia was classified under the International Classification of Diseases (ICD-10-CM) with the code M62.842. Numerous studies have confirmed its association with aging3,4, atherosclerosis, diabetes, malignancies, and cognitive disorders5, often leading to increased morbidity, mortality, and healthcare costs.

Early research primarily focused on the link between muscle performance decline and muscle mass loss6,7,8. However, substantial reductions in skeletal muscle function with age can occur with minimal muscle mass loss9,10,11,12,13, leading to a growing interest in muscle quality. Skeletal muscle density (SMD) has been widely used in muscle quality assessment, providing a quantitative measure of the density within the muscle’s region of interest (ROI), which can be directly obtained from computed tomography (CT). Numerous studies have shown a strong correlation between low muscle density and adverse health outcomes, yet the factors contributing to muscle density loss remain unclear.

Skeletal muscle fat infiltration, known as myosteatosis, defined as the accumulation of intramuscular and intermuscular fat within skeletal muscle, serves as an indicator of poor muscle quality14. A decrease in computed tomography (CT) attenuation by 1 Hounsfield unit (HU) corresponds to an increase in lipid concentration of 1 g/100 mL15. The primary cause of altered muscle density is skeletal muscle fat infiltration, assessed by CT parameters such as intermuscular adipose tissue (IMAT) and intramuscular adipose tissue content (IMAC). IMAT, defined as adipose tissue accumulation between muscles and muscle fiber bundles, is directly associated with declines in muscle mass, flexibility, and quality. This measure, quantified by CT attenuation ranging from − 190 HU to -30 HU16,17, has been widely used to assess muscle fat infiltration, with Amini et al.18 noting that 191 of the 388 studies reviewed used IMAT as an evaluation metric.

IMAC, introduced by Kitajima in 201019, represents extra-myocyte adipose tissue within muscle fiber bundles. Calculated as the ratio of CT attenuation in the multifidus muscle to CT attenuation in abdominal wall subcutaneous fat, IMAC serves as a standardized SMD indicator, eliminating variations from CT equipment, scan parameters, and individual differences among patients. Most studies on IMAC have focused on disease prognosis, with high IMAC scores often indicating poor outcomes20. However, fat infiltration distribution within skeletal muscle at the third lumbar vertebra (L3) level is uneven. Some individuals exhibit abdominal wall muscle atrophy, while others show multifidus muscle atrophy. Kitajima’s study evaluated only the multifidus muscle, and the IMAC of a single muscle may not adequately reflect extra-myocyte fat within muscle fiber bundles or provide an accurate measure of overall muscle fat content.

It is hypothesized that increases in IMAT and IMAC contribute to reduced SMD. However, it remains unclear which factor has a greater impact on SMD. Considering this background, our study aims to conduct a cross-sectional analysis to identify factors influencing SMD in the abdominal wall and paravertebral muscle groups (SM) at the L3 level in individuals over 50 years and in postmenopausal women, using opportunistic abdominal CT scans of inpatients.

Results

Baseline characteristics

A total of 499 inpatients were recruited for this study. Of these, 112 patients were excluded due to age mismatch, 84 were excluded due to poor health status or the presence of diseases significantly affecting muscle quality, and 11 were excluded because of poor image quality. Figure 1 illustrates the recruitment process. Ultimately, 292 inpatients were enrolled, consisting of 146 males (mean age 62.25 ± 9.99 years) and 146 females (mean age 69.16 ± 9.95 years), with an overall mean age of 65.71 ± 10.54 years.

Fig. 1
figure 1

Flow chart of participant selection for the study, with a total of 292 subjects enrolled.

No significant differences in age or body mass index (BMI) were observed between genders within each group (Tables 1 and 2). Muscle parameters for males and females are shown in Table 3. The inter-observer reliability, measured by the intraclass correlation coefficient (ICC), exceeded 0.75, indicating excellent repeatability (Table 4).

Table 1 Age comparison of male and female in each age group (\(\:\stackrel{-}{x}\)±s). Significant statistical difference* (p < 0.05).
Table 2 BMI comparison of male and female in each age group (\(\:\stackrel{-}{x}\)±s). BMI, body mass index. Significant statistical difference* (p < 0.05).
Table 3 Muscle parameters of male and female in each age group (\(\:\stackrel{-}{x}\)±s). SMA, skeletal muscle area; SMI, skeletal muscle index; SMD, skeletal muscle density; IMAT, intermuscular adipose tissue; IMAC, intramuscular adipose tissue content; SFA, subcutaneous fat area in abdominal wall; SFD, subcutaneous fat density in abdominal wall.
Table 4 The ICC of muscle parameters between the two observers (\(\:\stackrel{-}{x}\)±s). SMA, skeletal muscle area; SMD, skeletal muscle density; IMAT, intermuscular adipose tissue; SFA, subcutaneous fat area in abdominal wall; SFD, subcutaneous fat density in abdominal wall.

Correlation analysis and influencing factors of SMD

Table 5 presents the correlation analysis results between SMD and various variables, including age, gender, BMI, and other muscle parameters at the L3 level. The analysis showed a strong correlation between SMD and both IMAT and IMAC (p < 0.05), a moderate correlation with gender, age, SMA, skeletal muscle index (SMI), and subcutaneous fat area (SFA) (p < 0.05), and insignificant correlation with BMI (p > 0.05).

Table 5 Correlation between SMD and age, gender, BMI and other muscle parameters. Significant statistical difference* (p < 0.05), p value was obtained from Spearman’s correlation coefficient. R, correlation coefficient; p, p-value.

Table 6 details the impact of SMI, IMAC, IMAT, SFA, gender, and age on SMD. Multiple linear regression analysis identified IMAC, IMAT, and age as significant factors influencing SMD (p < 0.05). The order of influence was IMAT (β = -0.616), IMAC (β = -0.429), and age (β = -0.098). Specifically, as IMAT increased by 1 unit, SMD decreased by 0.968 HU; as IMAC increased by 1 unit, SMD decreased by 17.524 HU; and as age increased by 1 year, SMD decreased by 0.097 HU. SMI, SFA and gender did not significantly influence SMD (p > 0.05).

Table 6 Multiple linear regression analysis about influencing factors of SMD of abdominal wall and paravertebral muscle group at L3 level. Significant statistical difference (p < 0.05) *. B, partial regression coefficient; SE, standard error; β, standardized partial regression coefficient; t, t-statistic; p, p-value.

Discussion

Our analysis of 292 inpatients using abdominal opportunistic CT demonstrated that IMAT had a greater effect on SMD than IMAC and age. Previous studies may not have fully investigated the relative impacts of IMAT and IMAC on SMD.

Most research has focused on evaluating healthy individuals21,22. However, a significant proportion of participants in our study were elderly with underlying medical conditions. Jong Hyuk Lee et al. assessed 2,720 patients undergoing annual physical examinations, finding that 12% had a history of cancer, 27% had hypertension, 14% diabetes, 11% cardiovascular disease, and 5.4% chronic liver or kidney disease23. Relying solely on direct assessments of muscle quality in such populations can lead to inaccuracies. We selected inpatients for their clear and comprehensive medical histories, facilitating the exclusion of conditions affecting muscles, reducing experimental error, and enhancing reliability.

Measurement of muscle quality is also influenced by CT equipment, tube voltage, slice thickness, and iodine contrast use. Lamba et al. recruited 48 patients over 18 years who underwent unenhanced CT scans on both GE and Siemens 64-MDCT scanners within 12 months, observing that Hounsfield unit measurements for unenhanced abdominal soft tissues varied between the two manufacturers11. Lortie et al. retrospectively analyzed the effects of iodine contrast and tube voltage on muscle mass and quality evaluation, finding that iodine contrast significantly increased muscle density and area, while lower tube voltage resulted in higher muscle density and lower muscle area, with the impact on muscle density being greater than on muscle area10. Fuchs et al. noted an 11.64% decrease in muscle density (p < 0.0001) and a 1.11% increase in muscle area (p < 0.0001) when measured on 5 mm slice thickness images compared to 2 mm slices24. Therefore, controlling these conditions is essential to ensure measurement stability and reproducibility, minimizing errors. In this study, we used the same CT scanner with standardized parameters: tube voltage at 120 kV, scan slice thickness at 10 mm, reconstructed slice thickness at 0.625 mm, reconstruction window width of 350 HU, and reconstruction window level at 40 HU.

For CT-based muscle quality assessment, inter-observer, intra-observer, and inter-software variability also play roles. Barbalho et al.25 showed that muscle parameters measured by Slice-O-Matic and OsiriX software exhibited excellent consistency (ICC ≥ 0.954, p < 0.001). Our findings demonstrated excellent inter-observer consistency (ICC > 0.75), consistent with other studies26.

Our study demonstrated a strong negative correlation between IMAC and IMAT with SMD, identifying both as influential factors for SMD. Multiple linear regression analysis revealed standardized regression coefficients of -0.616 for IMAT and − 0.429 for IMAC, indicating that a one-unit increase in IMAT led to a 0.968 HU decrease in SMD, while a one-unit increase in IMAC resulted in a 17.524 HU decrease in SMD. These findings underscore the significantly greater impact of IMAT on SMD compared to IMAC. Brennan et al. demonstrated a significant association (p < 0.01) between lower muscle density in the elderly and higher intracellular lipid concentrations; however, SMD showed a weaker correlation with intracellular lipid concentrations, being more influenced by extra-myocyte lipids and independent of age, gender, race, and obesity27. IMAT represents lipid storage in adipose cells beneath the deep muscle fascia, encompassing visible lipids stored within intermuscular fibers and adipocytes28. Therefore, we infer that muscle density is predominantly influenced by the visible fat component located between muscles and fascicles, with minimal impact from intrafascicular and intramuscular cell fat. This aligns with Brennan’s findings27. Further research is needed to determine whether muscle function and metabolism relate to the ___location of fat infiltration within muscle.

Additionally, our study demonstrated a moderate association between age and SMD, with no significant influence from gender. Aging significantly contributed to muscle quality decline in both males and females. Multiple linear regression analysis revealed a standardized regression coefficient of -0.098 for age, indicating that SMD decreased by 0.097 HU annually. Previous studies consistently reported that IMAT increases with age. Our results also align with Tetsuro Hida29, who observed that aging is associated with reduced muscle area and increased fat deposition in muscle, even after controlling for gender and ethnicity. This suggests that the effect of age on muscle quality is universal, irrespective of spinal disease presence, muscle type, gender, or ethnicity.

Our study excluded patients with conditions that could severely affect muscle quality, such as subcutaneous oedema and hemorrhage, the evaluation of the effect of IMAT and IMAC on muscle density was limited in this pathological state. As in many previous large studies, we used a fixed threshold approach to define IMAT. But in a clinical setting, variations in scanning modes and acquisition schemes could potentially affect muscle measurements and limit the power of our analysis30. Further research in this area, conducted through animal testing studies, is recommended to provide more conclusive evidence. Recently, a novel approach to muscle segmentation is beginning to emerge. The method suggested by Mühlberg et al.31 was used to estimate an individualized threshold for the IMAT by using the attenuation of the subcutaneous adipose tissue. This method is currently used for the segmentation of the pectoral and upper femoral thigh muscles31,32, and extending this approach to large-scale multi-center research applications may enable more accurate assessment.

This study has several limitations. First, we evaluated the entire abdominal wall and paraspinal muscles, which show a good linear correlation with total body muscle33. However, delineating the entire muscle group was time-consuming. Future studies should separately analyze the abdominal wall and individual paraspinal muscles to identify the best indicators of changes in muscle mass and quality. Second, due to the difficulty in distinguishing between adipose tissue within muscle bundles and intramuscular lipids, CT attenuation only reflects the overall composition within the muscle. With advances in magnetic resonance spectroscopy (MRS), intramyocellular lipids (IMCL) warrant further investigation in future studies. Finally, our study was a single-center study, and the sample size of this study was relatively small and lack of younger cohort, necessitating further validation of these conclusions through larger-scale studies.

Conclusions

Our findings indicate that IMAT has a greater effect on SMD than IMAC and age, suggesting a distinct causal mechanism for SMD reduction. This insight may encourage a stronger emphasis on preventing IMAT increase.

Materials and methods

Study Design and participants

This retrospective study recruited 499 inpatients from Tianjin Hospital between July 2019 and July 2022. Clinical information and abdominal CT scans were routinely collected for each participant. Inclusion criteria were: (1) males aged ≥ 50 years and females in a postmenopausal state, and (2) abdominal CT images with optimal clarity, showing no density impact on SM, subcutaneous fat tissue in abdominal wall (SAT), or intra-abdominal fat at the L3 level, regardless of underlying causes such as edema, hemorrhage, or artifacts34,35. Exclusion criteria included conditions severely affecting muscle quality, such as musculoskeletal, neuromuscular, malignant, and chronic wasting diseases (e.g., severe diabetes, chronic obstructive pulmonary disease, chronic liver disease). The study followed the Helsinki Declaration and was approved by the Ethics Committee of Tianjin Hospital (2024YLS061). Due to its retrospective nature, Tianjin Hospital waived the requirement for informed consent.

Computed tomography acquisition

Full abdominal spiral CT imaging was performed on all participants using a 128-slice GE Revolution ES CT scanner (GE Medical Systems, LLC*). Scans covered the region from the upper edge of the T12 vertebra to the lower edge of the L5 vertebra. Scan parameters were set as follows: tube voltage, 120 kV4,36; tube current, 230 mA; table height, 85 cm; pitch, 0.948:1; and slice thickness, 10 mm.

Cross-sectional abdominal images were reconstructed at the L3 level on a GE AW4.7 workstation, aligned parallel to the upper endplate of the L3 vertebra, showing the longest transverse process slice (Fig. 2). Reconstruction parameters were: standard reconstruction algorithm, 0.625 mm reconstruction thickness, window width of 350 HU, window level of 40 HU, and a display field of view (DFOV) of 45 cm.

Fig. 2
figure 2

Selection of the measurement level. With regard to the abdominal coronal and sagittal CT images, transverse images at the L3 level were obtained, which were parallel to the upper endplate of the L3 vertebrae and demonstrated the longest slice of the transverse process.

Muscle quality assessments

Image J 1.53e (Wayne Rasband and contributors, National Institutes of Health, USA) was used for analysis. Muscle segmentation was manually performed using the “Polygon Selections” tool to outline muscle contours at the L3 level. The segmented SM included the psoas, paraspinal (erector spinae and quadratus lumborum), and abdominal wall (external and internal abdominal obliques, transversus abdominis, and rectus abdominis) muscles. SMA (cm2) and SMD (HU) values within the regions of interest were measured using the “Analyze-Measure” function, and muscle area was normalized for height in meters squared (m²) to obtain lumbar SMI (cm²/m²). Within these regions, a threshold of -190 HU to -30 HU was applied to distinguish muscle tissue from fat, allowing IMAT (%) measurement. The same method was used to outline SAT to obtain the cross-sectional area (SFA) and density (SFD) (Fig. 3).

Fig. 3
figure 3

Muscle segmentation and measurement. Female, 53 years old, BMI 27.34, SMA 101.81cm2, SMD 30.42HU, IMAT 4.89%, SFA 166.99cm2, SFD − 99.88HU, IMAC − 0.30. Measurement of SMA (a), SMD (a), IMAT (b) of the SM at L3 level. SFA (c) and SFD (c) was displayed at L3 level by Image J.

IMAC was calculated using the following formula:

IMAC = SMD / SFD.

All measurements were independently performed by two musculoskeletal radiologists with 10 and 20 years of experience, and mean values were recorded. Figure 4 shows a visual representation of the measurement process.

Fig. 4
figure 4

Abdominal wall and the whole paraspinal muscles at the level of L3 with different extent of fat infiltration in postmenopausal females and males over 50 years old.

All the abbreviations of muscle indicators measured at the L3 level are displayed as following:

SM: Abdominal wall and paravertebral muscle groups

SMD: Skeletal muscle density

SMA: Skeletal muscle area

SAT: Subcutaneous fat tissue of abdominal wall

SFA: Subcutaneous fat area of abdominal wall

SFD: Subcutaneous fat density of abdominal wall

IMAC: Intramuscular adipose tissue content

IMAT: Intermuscular adipose tissue

Statistical analysis

All statistical analyses were conducted using IBM SPSS Statistics Version 21 (USA). Continuous variables are expressed as mean ± standard deviation (\(\:\stackrel{-}{x}\)±s), with all continuous variables following a normal distribution. The significance level was set at α = 0.05.

The intraclass correlation coefficient (ICC) was used to assess interobserver reliability, with ICC > 0.75 indicating good agreement. An independent samples t-test compared age and BMI between men and women within the same age group. One-way ANOVA with LSD (chi-square test for group comparisons, p > 0.05) or Dunnett’s T3 (p < 0.05) was used to compare muscle parameters by gender and age group. The relationship between SMD and muscle measurements was analyzed using Spearman’s correlation coefficient. Multiple linear regression analysis identified and compared factors influencing SMD.