1. Introduction
The prevalence of obesity and overweight is increasing worldwide. According to the World Health Organization, 39% of adults aged 18 or over were overweight and 13% were obese in 2016. This population is affected by numerous complications and comorbidities [
1,
2,
3,
4,
5,
6] that may result from the accumulation of fat within and around tissues. Accumulation of abdominal fat and particularly visceral fat in obese people [
7,
8], is associated with metabolic syndrome and cardiovascular risks [
2,
6,
8,
9,
10]. Ectopic adipose tissue can also be located within the muscle in obese people [
11,
12] but also in myopathic patients [
13,
14], elderly [
15,
16] or after a period of inactivity [
17] and represents a risk for metabolic disorders [
11,
18] and muscle function [
19]. As such, it is generally considered that fat deposit generates lipotoxicity [
20].
It is essential to be able to assess intramuscular fat (IMF) content in clinical practice. This fat depot is composed of fat located either within the muscle fibers (intramyocellular lipids; IMCL) or between the muscle fibers (extramyocellular lipids; EMCL). The latter can also be found between and around skeletal muscles and constitutes a genuine adipose tissue referred to as intermuscular fat adipose tissue (IMAT) [
19,
21] when evaluated with magnetic resonance imaging (MRI-T1-sequence or mDixon technique) or computed tomography. The location of the lipid depot has clinical significance: IMCL content is known to be positively correlated with insulin resistance [
22,
23] and negatively correlated to muscle function [
24]. EMCL and IMAT have also been associated with dysglycaemia, including insulin sensitivity and glucose tolerance [
21], and with reduced muscle function [
25], but their effects seem distinct from those of IMCL [
21]. In addition, the athlete paradox revealed that chronic exercise training actually increased the muscle triglyceride component of IMCL, but this effect was independent of exercise effects on IMAT/EMCL [
26]. In clinical practice, it is therefore important to distinguish the accumulation of IMCL and EMCL.
Proton magnetic resonance (
1H-MRS) stands out as the more accurate and convenient non-invasive method to quantify IMCL and EMCL [
27]. Unfortunately, this is an expensive method that limits its use in daily clinical practice. Qualitative grayscale analysis, known as muscle echo intensity (EI) appears the less expensive and the most accessible alternative [
27]. Some studies showed that increased EI is negatively correlated to muscle strength [
28,
29,
30]. Initially, EI was used in agribusiness to identify meat quality and the percentage of IMF of pig and beef meat [
31,
32]. However, EI is only a qualitative approach. To overcome this limitation and provide a quantitative evaluation of IMF, one study related IMF content assessed by MRI to EI measured by an ultrasound device in brightness mode (B-mode) [
33]. Prediction equations were generated between muscle EI and IMF. As subcutaneous fat thickness can alter or reflect the absorption of ultrasound waves in deeper tissues [
34,
35], the authors took this parameter into account in their equations, which significantly improved the prediction. However, this study did not differentiate IMCL and EMCL in the quantification of IMF. To date, only one study attempted to quantify IMCL and EMCL from EI in the vastus lateralis (VL) and the biceps femoris (BF). Akima et al. reported significant correlations between EI and EMCL but did not observe significant relationships between EI and IMCL. This could be accounted for by the low spatial resolution of ultrasound imaging or the fact that these authors did not account for the potential confounding effect of subcutaneous fat thickness [
36]. Another limitation of these studies is that they propose prediction equations that are specific to a particular ultrasound device. As proposed by Young et al. [
33], it would be advisable to test other devices because the muscle EI values depend on ultrasound transducer characteristics and image post-processing techniques. Pillen et al. [
34] demonstrated that although ultrasound devices create images that are different on a pixel level, EI can be comparable between ultrasound devices on a macro level, after the application of a conversion factor. That would allow the creation of generic prediction equations relating IMCL, EMCL and IMF to EI.
The purpose of the present study was to compare four ultrasound devices against MRI to quantify IMCL, EMCL and IMF of two plantar flexors [gastrocnemius medialis (GM) and soleus (SOL)] and two knee extensor [VL and vastus medialis (VM)] muscles from EI, corrected for subcutaneous fat thickness and to generate prediction equations.
2. Materials and Methods
2.1. Subjects
A total of thirty subjects (11 women and 19 men) were recruited in the present study. Participants displayed several body mass index (BMI) values and disparate physical activity levels (from sedentary to highly endurance-trained participants), in order to obtain a varied sample of muscle adiposity. The study was conducted in accordance with the Declaration of Helsinki and approved by the ethics committee of Ritsumeikan University (BKC-IRB-2018-062). All the participants were fully informed and gave their written consent before any testing was conducted.
2.2. Design
Participants were asked to come to the laboratory on one occasion, to undergo all the testing procedures. The ultrasound tests were performed first and then the MRS measurements. This allowed us to perform the MRS scanning in the specific muscle area assessed by ultrasound. Body composition was also assessed during this experimental session.
2.3. Body Composition Assessment
Whole-body and right leg fat-free mass and fat mass were assessed by impedancemetry (InBody 770, InBody co., Seoul, Republic of Korea).
2.4. Ultrasound Experimental Measurement and Analysis
A total of four lower limb muscles were investigated in this study (VM, VL, GM, SOL;
Figure 1). In each participant, the right leg was assessed by the same experimenter in a resting, supine position on an examination table. Four ultrasound devices, differing in terms of post-processing techniques and transducer characteristics, were used for analysis: LOGIQ (Logiq S7, GE Health Care, Chicago, IL, USA), SUPERSONIC (Aixplorer, Supersonic Imagine, Aix-en-Provence, France), ALOKA (Noblus, Hitachi–Aloka Medical Japan, Tokyo, Japan) and TELEMED (Echo Blaster 128 CEXT-1Z, Telemed Ltd., Vilnius, Lithuania) devices. Real-time brightness mode (B-mode) was used for all devices with specific parameters. Settings for each ultrasound device are reported in
Table 1. The scanning depth was adjusted for participants with greater subcutaneous fat thickness to gather maximal muscle area. The other parameters were set in order to obtain a clear visualization of the muscle under consideration and its aponeuroses.
To normalize measurements, the leg and thigh lengths were measured. Thigh length was defined as the distance between the antero-superior iliac spine and the superior lateral aspect of the patella. Leg length was defined as the distance between the inferior lateral aspect of the patella and the calcaneus. Ultrasound scans were performed between the lower- and middle-third of the thigh (for VM and VL, respectively) and between the third and fourth quarter of the leg length (for GM and SOL, respectively). These anatomical landmarks were marked on the skin with permanent ink during ultrasound scanning and later with oil capsules taped to the skin during MRS scanning. Ultrasound scanning was performed with a generous amount of ultrasound gel applied to the probe to avoid pressure on the skin and improve imaging quality. Each muscle was scanned three times for each ultrasound device and each participant. A total of forty-eight scans were performed for each participant. Images were exported to a jpeg format.
ImageJ (U.S. National Institutes of Health, Bethesda, MA, USA) was used for image analysis. Muscle EI was assessed using grayscale calculation by the same experimenter. For each image, three values of muscle EI, examined muscle area and subcutaneous fat thickness were calculated. For each muscle and ultrasound device, the average of nine image values (from muscle EI, area and subcutaneous fat thickness) was computed for further analysis.
2.5. MRS Acquisition and Analysis
Non-invasive
1H-MRS measurement of the muscles (VM, VL, GM, SOL;
Figure 1) was performed on a 1.5T MR system (Signa HDxt, GE Medical Systems, Buckinghamshire, England) with an 8ch body array coil. The subject lay in the supine position with their right leg carefully positioned parallel to the main magnetic field. Single voxel MRS measurement was performed using the Point Resolved Spectroscopy (PRESS) sequence (TR/TE 2000/35 ms, 20 × 20 × 20 mm
3, 32 acquisitions). All subjects were instructed to maintain their usual physical activities and dietary habits for a few days before the measurement. Outer volume suppression slabs were used around the voxel of interest to diminish the signal of subcutaneous fat located nearby the voxel. The analyses and estimation of the absolute concentration of IMCL and EMCL were accomplished by LCModel software ver6.3 with Weis customized calculation [
37]. An unsuppressed water signal was used as an internal standard. During data processing T1 and T2 relaxation effects of the water reference were taken into account with the LCModel’s control parameter, atth2o, determined using the following equation: exp (−TE/T2) [1−exp (−TR/T1)] [
37,
38]. The peak chemical shift of IMCL was adjusted as 1.3 ppm and that of EMCL was 1.5 ppm. Concentrations of EMCL and IMCL were computed as millimoles per liter of muscle tissue (mM) and converted into millimoles per kg wet weight as recommended by Weis et al. [
37]. Finally, IMF was calculated as the sum of EMCL and IMCL. Typical example spectra are illustrated in
Figure 2.
2.6. Muscle and Echographer-Specific Equations
Correlations between raw EI and IMCL, EMCL and IMF from MRS were made for all the muscles for each ultrasound device. To account for the potential confounding effect of subcutaneous fat thickness on muscle EI, the correction factor proposed by Young et al. [
33] was also applied to raw EI to generate new equations, corrected for subcutaneous fat thickness, as proposed in the equation below.
where y
2 is the corrected EI, y
1 is the raw EI, x is the subcutaneous fat thickness and cf is the Young correction factor [
33], i.e., 40.5278.
2.7. Statistics
Linear regression models were used to build, examine and compare relationships between muscle EI and fat content (IMF, IMCL, EMCL) for all muscles and ultrasound devices. The strength of the positive and negative correlations were evaluated using the Evans [
39] guidelines: very weak: 0.00–0.019, weak: 0.20–0.39, moderate: 0.40–0.59, strong: 0.60–0.79, very strong: 0.80–1.0. Then, for significative relationships, ANOVAs were conducted to compare the slopes and Y-intercepts between relationships. When the ANOVA revealed a significant difference, a Tukey test was used as a post-hoc test.
Examined muscle area, raw and corrected EI and fat thickness were also compared between devices. Distribution normality and homogeneity of variances of these variables were tested using a Shapiro–Wilk normality test and Bartlett’s test, respectively. Then, one-way ANOVAs with repeated measures were run. When the ANOVA revealed a significant difference, a Tukey test was used as a post-hoc test. Statistical analyses were conducted on GraphPad 7 (GraphPad Software, San Diego, CA, USA). Values are expressed as mean ± standard deviation (SD). The α-level of significance was set at p < 0.05.
4. Discussion
The purpose of the present study was to compare prediction equations of IMCL, EMCL and IMF from EI between four ultrasound devices. IMCL was poorly correlated with muscle EI, whereas EMCL and IMF were moderate to strongly correlated with muscle EI of GM, SOL, VL and VM. All the relations were improved using the Young corrected factor, showing the potential confounding effect of subcutaneous fat thickness on muscle EI measurements for all muscles. Slopes of the relationships were comparable between the four devices, but y-intercepts differed, with few exceptions. Interestingly, these differences disappeared when Young corrected EI were computed against MRS data, giving the possibility to generate generic prediction equations.
4.1. Quantification of IMF, IMCL and EMCL with Ultrasound
As far as we know, only one study tried to quantify EMCL and IMCL using muscle raw EI [
36]. The authors found a moderate correlation between raw EI and EMCL in VL and BF (r = 0.65). In the present study, r values were greater for all the ultrasound devices concerning VL but comparable for other muscles investigated such as GM (excepted for Supersonic), SOL (excepted for Supersonic) and VM. Results from this experiment confirm that raw EI values can be used to quantify EMCL.
Conversely, Akima et al. [
36] failed to find correlations between IMCL and raw EI in VL and BF muscles. The lack of correlation could be explained by the low spatial resolution of ultrasound imaging or by the fact that these authors did not account for the potential confounding effect of subcutaneous fat thickness [
36]. Indeed, Young et al. [
33] found an independent influence of subcutaneous fat thickness on muscle EI, with alteration of reflection or absorption of the ultrasound waves in deeper tissues [
34,
35]. Consequently, Young et al. took the subcutaneous fat thickness into account in their prediction equations. This improved the quantification of IMF from EI values. In the current study, we applied this correction factor for the prediction of EMCL and IMCL. This translated into a modest increase in the strength of the relationship between EMCL and muscle EI but failed to improve the relationship between corrected muscle EI and IMCL. It can be concluded that the lack of correlation between EI and IMCL seems to be explained by the low spatial resolution of ultrasound imaging rather than the confounding effect of subcutaneous fat thickness.
Nevertheless, one should not conclude about the complete inability of ultrasound imaging to quantify IMCL. The EI-IMF relationships were stronger than EMCL-EI relationships, suggesting that ultrasound imaging can quantify at least part of the IMCL included in the IMF content. This good ability of ultrasound imaging to quantify IMF is consistent with previous studies comparing muscle EI to the percentage of IMF quantified from a muscle biopsy sample [
40,
41] or from MRS data [
36]. As previously mentioned, corrected EI from subcutaneous fat thickness improved EI-IMF relationship quality. Indeed, our results showed moderate to very strong correlations for Young corrected EI and IMF relationship in GM, SOL, VL and VM.
4.2. Comparison between US Devices
The ability to quantify muscle fat content from EI may be ultrasound device dependent. The present study is the first to compare different ultrasound devices to assess EMCL and IMF in lower limb muscles. Some anecdotal differences were observed in y-intercept and slopes between devices when computing raw EI values to derive raw EI-IMF and raw EI/EMCL relationships in GM, SOL and VL. This may be ascribed to the fluctuations of ultrasound waves encountering different tissues [
14] and the subsequent image processing, which is specific to each device. Interestingly, the EI-EMCL and EI-IMF relationships did not differ between ultrasound devices when the subcutaneous thickness was considered. This allowed the production of generic equations for the indirect quantification of IMF and EMCL from Young corrected EI data. This result has important practical implications. Indeed, several studies have assessed muscle quality from EI and their relation with muscle strength and function in several populations [
14,
28,
29,
42]. However, this qualitative approach is limited by the fact that EI values are not comparable between ultrasound devices and thus between studies unless a conversion factor is given, which is currently unavailable in the literature. That would require correcting EI values obtained from several devices and settings [
34,
43]. Interestingly, the current results suggest that correcting EI values for the confounding influence of subcutaneous fat thickness allows a generic and quantitative approach to muscle fat content. This opens new avenues for clinicians and scientists interested in the impact of muscle fat content on physical function [
44].
4.3. Limitations
Some limitations of our experimental approach should be acknowledged. First, we investigated two plantar flexors and two knee extensor muscles. The proposed equations are muscle-specific and need to be extended to other muscles, in particular the dorsi-flexors and knee flexors and upper limb muscles which are seldom studied. Second, unlike Young sex-specific equations [
33], our participants’ sampling did not allow us to generate sex-specific equations. Third, the muscle area examined for EI measurement was bigger than the muscle area examined with MRS, because it was impossible to reproduce exactly the location of the MRS examination zone with ultrasound, although attention was paid to evaluating the same location (slice) in the longitudinal direction. For EI measurements, we chose to examine a muscle area that would encompass the MRS area. Given the heterogeneity of EMCL/IMF distribution within muscles, the correlation coefficients between EI and EMCL/IMF measurements would have been stronger, should the EI measurements be performed on the same voxels used for MRS. Fourth, the purpose of this study was not to uniformize the settings and images between US systems, but to compare four ultrasound devices against MRI, implying that settings would differ between systems. This study was thus designed in an ecological way. The settings were chosen for each system with the aim of producing the clearest possible image (as practiced in the clinical context). However, to improve the approach of the current study, normalization to a reference material (
i.e., calibration phantom) is required. That would reduce system-dependent variance, render the data fully independent of instrumental settings and ultimately would allow producing universal equations for the quantification of muscle fat content. Finally, it remains unknown if the generic equations produced in the current study are valid for populations such as severely obese patients, where a very high subcutaneous fat thickness may potentially impair the ability to quantify IMF and EMCL from EI values.