1. Introduction
Skeletal muscle represents the primary site for postprandial glucose disposal [
1]. Lipids are important for energy production, with about two-thirds of resting skeletal muscle energy production coming from lipid oxidation. The accumulation of lipids inside muscle cells is a key energy source for physical activity. However, while this may be beneficial for athletes, excess accumulation, particularly in diseases such as diabetes and ageing muscle, can lead to lipotoxicity [
2]. Lipotoxicity results from the ability of fatty acids (FA) to act as detergents, to affect acid-base homeostasis, and produce highly bioactive lipids, such as ceramides. Larger accumulation of lipids inside muscle cells leads to impaired mitochondrial oxidative phosphorylation in ageing and also in subjects with diabetes [
3]. This process can then lead to fat infiltration in muscle which is associated with frailty and muscle loss [
4]. Type 2 diabetes (T2DM) is a chronic metabolic condition where every aspect of the body’s metabolism is altered [
5]. Insulin resistance (IR) is a common feature of T2DM and an early metabolic abnormality [
6]. Insulin resistance leads to impaired glucose absorption, promotes hepatic gluconeogenesis, manifests a serious breakdown in lipid dynamics, and is associated with the metabolic syndrome, especially with its effects on the liver and muscle [
7]. It is known that lipids are stored as triglycerides in the adipose tissues or in lipid droplets in the cytoplasm of non-adipose cells. Insulin resistance promotes hydrolysis of triglycerides, release of free fatty acids, reduced esterification and re-esterification of fatty acids in adipose tissues (AT), reduces glucose uptake in muscle cells, impairs glycogen synthesis, and increases gluconeogenesis, lipogenesis, and synthesis of CRP in the liver [
8,
9]. This promotes ectopic accumulation of lipids in non-adipocyte tissue such as skeletal muscle, abdominal visceral and hepatic fat depots and pancreatic beta cells [
10,
11] along with macrophage infiltration and inflammation [
12]. This fat deposition in non-adipose tissues leads to further IR [
13] and oxidative stress, and has deleterious effects including tissue damage (lipotoxicity) [
14]. In addition, ageing is also associated with larger accumulation of lipids within skeletal muscle leading to muscle loss [
2]. Two pools of lipids are found in skeletal muscle tissues: intra-myocellular lipids (IMCL), stored in spheroid droplets in muscle cells adjacent to mitochondria and extra-myocellular lipids (EMCL), distributed over large regions of muscle fascia. There is also a large interest in understanding the fat composition within these lipid pools. The effects of fatty acid on metabolic signaling and energy metabolism are modulated by degree of unsaturation [
15]. Monounsaturated fatty acids and polyunsaturated fatty acids are synthesized by the desaturation of saturated fatty acids by an oxidative reaction catalyzed by the iron-containing microsomal enzyme stearoyl CoA desaturase. Therefore, this enzyme determines the degree of desaturation, and its activity, accordingly, has profound effects on lipid metabolism [
15]. There is also recent evidence that the accumulation of saturated IMCL is associated with insulin resistance [
16,
17]. It has also been shown that IMCL content in different muscle compartments exhibit differences in muscle cellular adaptation for differences in the type of exercise and muscle fiber types [
18].
In addition to the lipid metabolism, there is also clinical interest in the detection and estimation of metabolites including choline, carnosine, taurine, myo-Insositol and creatine. Creatine supplementation has been shown to increase the muscle creatine and phosphocreatine (PCr) content which is essential for high-intensity exercise performance [
19,
20]. Choline plays an important role in phospholipid synthesis and triglyceride metabolism contributing to structure and function of cell membranes and muscle cells [
21]. Low concentration of choline is associated with several metabolic changes involving anti-fibrotic effects [
22] including muscle wasting [
23], and also plays a role in muscle contraction [
24]. It has also been shown that inositol is involved in excitation – contraction coupling in skeletal muscle [
25]. Taurine is required for proper muscle functioning and has the potential to reduce muscle inflammation and improve muscle strength [
26]. Recently taurine deficiency has been shown to be associated with accelerated ageing [
27]. Inositols via its major isomers (myo-inositol and D-chiro-inositol) participate in both insulin signaling and glucose metabolism [
28] by influencing distinct pathways. For the first time, we also report the feasibility of detecting myo-Inositol in vivo in skeletal muscle.
Over the last few years there has been steady progress in advancing in vivo magnetic resonance spectroscopy (MRS) methodology for quantification of lipids and other metabolites in skeletal muscle [
29]. One-dimensional MRS techniques result in severe overlap of lipid resonances and also hinders detection of wide range of metabolites at lower field strengths [
30,
31]. On the other hand, two-dimensional (2D) localized correlated spectroscopy (L-COSY) techniques have been demonstrated to minimize this problem of spectral overlap by spreading resonances into a second dimension, thereby improving the dispersion of spectral resonances [
30,
31,
32]. The L-COSY techniques permit estimating the IMCL and EMCL content in a specific compartment, and also allow quantitative assessment of lipid unsaturation in vivo. However, single voxel 2D MRS techniques are limited to single muscle compartment with limited spatial coverage and is not efficient. Recently using a multi-voxel 2D COSY technique, Nagarajan et al. reported differences in changes in myolipid metabolism between type I (slow twitch) and type II (fast twitch) muscle fibers, with type I fibers having higher IMCL content relative to type II (fast twitch) muscle fibers in obesity, perhaps related to IR [
33,
34]. The 2D EPSI spectra can resolve the overlapping peaks of various lipid resonances and also the metabolites including Choline, Creatine, Taurine, and myo-Inositol.
Central obesity is a major player in T2DM [
35,
36] and it is therefore important to assess the abdominal fat tissues in various compartments within the abdomen. Similar to IMCL in skeletal muscle, triglycerides in adipose tissue and liver are found to be associated with insulin resistance [
37]. Non-alcoholic fatty liver disease (NAFLD), a common liver disease characterized by hepatic triglyceride accumulation and chronically elevated serum aminotransferase levels, reported to have high prevalence in patients with T2DM [
38]. In addition to liver, there is also emerging evidence on nonalcoholic fatty pancreas disease (NAFPD) for describing fat accumulation in pancreatic cells. Some cross--sectional studies have shown that pancreatic fat content increased in subjects with T2DM [
39]. In adults, NAFPD has been associated with obesity, metabolic syndrome, and insulin resistance / T2DM [
40,
41].
In this study, a major goal was to validate lipid pools (IMCL, EMCL, unsaturation index) and more metabolites using a novel accelerated five-dimensional (5D) echo-planar correlated spectroscopic imaging (EP-COSI) in the calf muscles of T2DM patients and age-matched healthy controls, and to correlate the lipid and metabolite ratios with the abdominal and ectopic fat in liver and pancreas.
2. Materials and Methods
2.1. Human Subjects
Twelve T2DM patients (5 male, 7 females; age = 60.3 ± 8.6 years; BMI = 28.1 ± 6.3 kg/m2), 9 age-matched healthy controls (AMHC) (6 males, 3 females; age = 60.9 ± 7.8 years; BMI = 25.3 ± 2.8 kg/m2), and 10 young healthy controls (YHC) (9 male, 1 female; age = 28.3 ± 1.8 years; BMI = 24.4 ± 3.3 kg/m2) were enrolled in the study. All the subjects provided informed consent in accordance with the institutional review board (IRB) at the University of California at Los Angeles.
The T2DM subjects were recruited from the UCLA Gonda Diabetic Center and from the Southern California area. The AMHC were recruited from the UCLA community, and the general community and family members of the patient subjects. The medical history of T2DM patients was assessed. The exclusion criteria for this study were as follows: history of cardiovascular or neuromuscular disease, HbA1c% > 9.0%, treatment with any medications that alter insulin sensitivity, documented diabetic micro- and macro-angiopathy, severe renal or hepatic disease, malignancy, or chronic inflammatory diseases such as rheumatoid arthritis. Other exclusions included subjects with claustrophobia, metallic implants such as pacemakers, defibrillators and aneurysm clips and certain prostheses contraindicated for MRI. Study inclusion criteria for healthy controls were good health with sedentary lifestyles, without a family history of diabetes, and women with regular menstrual cycles.
2.2. Clinical and Biochemical Data
All T2DM subjects underwent blood sample collection at fasting, and assessment of anthropometry by standardized methodologies. The collected data included age, gender, BMI, and levels of the serum creatinine (Cre), aspartate aminotransferase (AST), alanine aminotransferase (ALT), total bilirubin (TB), alkaline phosphatase (ALP), fasting glucose (Glc), glycated hemoglobin (HbA1c), total cholesterol (TC), triglyceride (TG), high-density lipoprotein (HDL), and low-density lipoprotein (LDL).
2.3. MRI / MR Spectroscopic Imaging
MRI / MRS imaging was performed on a 3T Prisma MR scanner (Siemens Medical Solutions, Erlangen, Germany) running on the VE11C platform.
2.3.1. MR Spectroscopic Imaging of Calf Muscles
MR spectra of calf muscles were acquired using a 15 channel ‘receive’ and a single channel ‘transmit’ knee coil in the feet first supine position. An enhanced version of the accelerated 5D EP-COSI technique was used to obtain spectra from the right calf of each subject with the largest diameter of the calf muscle positioned in the middle of the coil [
42]. The standard 5D EP-COSI sequence which uses (90°–180°-Δt1 -90°-t2) for localization was modified by employing a pair of adiabatic full passage pulses at the refocusing 180
0 RF pulse. Position of the volume of interest (VOI) was determined using the axial, coronal and sagittal images to assure the coverage of soleus, tibialis anterior and gastrocnemius muscle regions. Interference from vascular structures and gross adipose tissue deposits was minimized by appropriate positioning of the VOI. Anatomical MRI studies included the following scans: a three-plane localizer MRI (Repetition time/Echo time (TR/TE) = 14/5 ms, field of view (FOV) = 45 × 45 cm
2, a 320 × 160 data matrix, one average) and a T1-weighted (T1W) axial, coronal and sagittal spin-echo MRI (TR/TE = 600/2.7 ms, FOV = 16 × 16 cm
2, 320 × 256 data matrix, one average). The following parameters were used for acquiring the accelerated 5D EP-COSI acquisition with a maximum echo sampling scheme: TR/TE = 1500/35 ms, average = 1, matrix size = 16 x 16 x 8, voxel resolution = 3.37 cm
3, 64 Δt1 increments, 512 bipolar echo pair, acceleration factor = 8 along k
y, k
z and t
1 dimensions, F
1 and F
2 bandwidths of 1250 Hz and 1190 Hz respectively with scan time of ~25 min. Outer volume saturation bands were included outside the PRESS volume of interest and manual B
0 shimming was performed to minimize B
0 inhomogeneity over the localized VOI to achieve a line width of ~12-14 Hz. A fully sampled non-water-suppressed (NWS) scan was acquired with only the first t
1 increment to be used for eddy current correction and as a reference for coil combination, adding 30 seconds to the total scanning duration.
2.3.2. Abdominal Fat Imaging
Abdominal MRI was performed using a 3D GRE VIBE 6-point Dixon sequence. The acquisition parameters were voxel size = 1.187 x 1.187 x 3 mm
3, matrix size = 320x240, TR=8.85 ms, bandwidth = 1080 Hz/px, flip angle = 500 and shortest possible TEs (equidistant with TE1 = 1.23 ms and an echo time shift of 1.23 ms), a GRAPPA (GeneRalized Autocalibrating Partially Parallel Acquisitions) acceleration factor of 2, and one average. Fifty-two slices were acquired between vertebral levels T12 and L5. Using the water-only and fat-only images extracted from the multi-point Dixon data, fat fraction and water fraction were calculated by the in-line MR image reconstruction software as parametric maps. A gold-standard single-voxel
1H MRS of the liver was performed using a 15-second breath-hold, T2-corrected multi-echo spectroscopic sequence (HISTO, High speed multiple echo acquisition) as described in [
43].
2.4. Data Processing and Analysis
2.4.1. Spectroscopic Data Processing
The acquired undersampled data were reconstructed using a using a Group Sparsity-based compressed sensing algorithm [
42]. The reconstruction was performed over each coil separately, and the individual coil reconstructions were combined as a sum-of-squares. Reconstructed data were further post-processed using a series of steps comprising scaling, spatial reordering, phase correction, resolving averages and oversampling [
42], with a custom MATLAB-based program.
For the 2D spectra recorded in the soleus, tibialis anterior, and gastrocnemius muscles regions of each subject, absolute volumes for the following metabolites were quantified using peak integration. Intramyocellular lipids: IMCL1 and IMCL2 ((5.3 ppm, 2.7 ppm) and (5.3 ppm, 2.0 ppm), extramyocellular lipids: EMCL1 and EMCL2 (5.45 ppm, 2.85 ppm) and (5.45 ppm, 2.15 ppm), choline, Ch_d (3.2 ppm, 3.2 ppm), creatine, Cr_d (3.0 ppm, 3.0 ppm), Cr_3.9 (3.9 ppm, 3.9 ppm), carnosine, Car (8.0 ppm, 8.0 ppm), taurine, Tau (3.4 ppm, 3.4 ppm), myo-inositol, mI (3.5 ppm, 3.5 ppm), triglyceryl fats, TGFR1 (5.4 ppm, 4.3 ppm) and TGFR2 (4.3 ppm, 5.4 ppm), unsaturated fatty acids, Unsat_CH = CH, FAT_5.4 (5.4 ppm, 5.4 ppm), and poly-methylene protons of fat, FAT_1.4 (1.4 ppm, 1.4 ppm). We identified the diagonal peak at [3.5 ppm, 3.5ppm] as myo-Inositol based on earlier ex vivo magic angle spinning spectroscopy of muscle tissues [
44].
Since total muscle creatine (Cr_d) content does not change with T2DM, it was used as an internal reference for the other peaks in the muscle spectrum. The peak ratios were calculated by dividing the absolute peak volume by that of the Cr_d peak. The unsaturation indices (UI) of IMCL and EMCL were estimated using the ratio of the cross-peaks, IMCL1/IMCL2 and EMCL1/EMCL2, respectively [
29,
30,
43]. Average EMCL (Avg EMCL) and IMCL (Avg IMCL) were calculated from EMCL1, EMCL2 and IMCL1, IMCL2 respectively.
2.4.2. Abdominal Fat Segmentation
We used an image analysis software slice-O-matic (Ver. 5.0, Tomovision, Canada) to quantify subcutaneous adipose tissues (SAT), visceral adipose tissues (VAT), and total abdominal fat (TAT). A single trained observer (MKS) guided by an experienced radiologist (ERF/SR) performed the image analysis. The VAT and SAT areas (cm2) were computed automatically from the respective tissue regions in each slice by summing the pixels from the given tissue and multiplying by the pixel surface area. Finally, tissue volume (cm3) for each slice was calculated by multiplying the tissue area (cm2) by slice thickness.
2.4.3. Liver and Pancreatic Fat
Hepatic fat fraction (HFF) was measured by ROIs selected in homogeneous sections of the liver. To calculate pancreatic fat fraction (PFF), ROIs were placed in the head (caput, PFFHead) and body/tail (corpus + cauda, PFFBody+Tail) of the pancreas. Pancreatic regions were segmented into two parts: head and body/tail. To avoid segmentation error, we did not calculate PFF separately in body and tail. The mean pixel value obtained through the software was scaled down by a factor of 1000 to counteract scaling performed during image acquisition on the MRI scanner.
For HISTO, fat fraction values (HFF_MRS) and R2 (1/T2) relaxation rate of water (R2_WAT) were calculated using the vendor-provided software package as described in [
43].
2.5. Statistical Analyses
All statistical analyses were performed with SPSS (V25.0, IBM, Chicago, IL) software. A p value < 0.1 was considered statistically significant. The independent sample t-tests were performed to assess differences in lipids and metabolite ratios between T2DM and AMHC, and between AMHC and YHC groups. Pearson’s correlation was used to determine the association between MRS lipids/metabolites and MRI-based fat estimates, and also, between MRI/MRS parameters and blood chemistry in the T2DM and AMHC groups.
4. Discussion
There is a large interest in understanding the ectopic fat accumulation in progressive metabolic diseases [
44]. Quantification of IMCL and EMCL during ageing and diabetes is of large clinical interest due to its association with inflammation and muscle loss [
4]. The contribution of individual muscle compartments with different combination of fiber types and its association with skeletal muscle metabolism is not known. In that context, there is a large interest in implementation of advanced MRS methodology for investigating various muscle compartments.
Even though conventional one-dimensional MRSI techniques can provide information on distribution of IMCL / EMCL in various compartments [
47], it is a challenge to investigate the composition of the lipids at clinical field strengths. In this context, there is also clinical interest in understanding the role of saturated and unsaturated lipids within these lipid pools in the development of insulin resistance [
16,
17]. To address this issue, we earlier implemented single voxel based localized 2D MRS techniques to separate the overlapping lipid pools and also to resolve various metabolites [
30,
31]. However, single voxel-based techniques have limited spatial coverage involving acquisition from one muscle compartment. In this study, we implemented spatially resolved 2D MRS by accelerated 5D EP-COSI along with body composition imaging to estimate the lipids and metabolites in young and old healthy subjects along with diabetic subjects. We were able to measure the lipid composition within IMCL and EMCL along with various metabolites including choline, carnosine, and taurine. Recent work has shown that circulating Taurine deficiency may be a driver of aging [
27]. In that context, in vivo detection of Taurine in skeletal muscle will open new possibilities of trials with Taurine. We also identified the diagonal peak at [3.5ppm, 3.5ppm] as myo-Inositol based on the earlier ex vivo HRMAS of muscle tissues.
We also evaluated the body composition including subcutaneous fat (SAT), visceral fat (VAT), hepatic fat fraction (HFF) and pancreatic fat fraction (PFF) in age-matched healthy controls and diabetic subjects. We observed lower unsaturation of IMCL and EMCL in GAS, SOL, and TA in T2DM compared to AHMC group though not statistically significant. It is known that the fatty acid synthase enzymes can become defective in diabetes [
15]. We also observed higher choline levels in both T2DM and AHMC group compared to YHC. Choline is an essential micronutrient and plays an important role in several metabolic pathways. Low concentration of Choline is also associated with muscle wasting [
22]. We also observed higher levels of myoinositol in the T2DM group. It is known that inositols are involved in the in excitation-contraction of coupling in skeletal muscle [
24]. The higher choline and inositol levels in the T2DM group might be due to compensatory mechanisms to protect muscle loss.
The correlation analysis revealed an association between IMCL unsaturation index and abdominal fat content. In light of the close relation between IMCL and insulin sensitivity, these results can help to explain the impact of dysfunctional adipose tissue in T2DM. Both EP-COSI and Dixon MRI showed the same trends in differentiating T2DM and healthy controls. However, given the small sample size in this study, we didn’t observe any significant differences between the two groups. In AMHC, visceral fat was correlated with IMCL, which was also found in previous studies in lean and obese children [
48]. We did not see the correlation in the T2DM group.
In a recent study by Sarma et al. adipose tissue distribution including SAT, VAT, HFF and PFF in T2DM patients and these findings were compared with nondiabetic age-matched adults and young healthy subjects [
43]. There was a high correlation between hepatic fat fraction estimated by 6-point Dixon MRI and HISTO-MRS. Increased SAT and VAT reaffirmed that central obesity is connected with the evolution of T2DM favoring insulin resistance. VAT was greater in the T2DM group than AMHC and YHC, but the difference between groups was larger for the YHC. In addition, almost equal percentages of SAT and VAT were observed in T2DM suggesting that the VAT accumulation may play an important role in the etiology of T2DM. The T2DM patients had significantly greater hepatic fat content than both AHMC and YHC subjects, indicating that T2DM is associated with increased fat accumulation in the liver. The HFF percentage was 44.6% higher in individuals with T2DM compared to age-matched healthy controls and 64.4% higher compared to young healthy controls.
A limitation of the present study is that it was a cross sectional investigation using a small number of subjects. While we noted significant differences in VAT, HFF, and PFF of T2DM patients in comparison to AMHC, large cohort studies will be required to determine whether the observed changes are of primary importance in the development of insulin resistance, impaired glucose tolerance, and progression to T2DM. However, this preliminary study using 5D EPCOSI in conjunction with body composition imaging reveals important differences between T2DM and healthy subjects.
Author Contributions
Conceptualization, and designed the experiments, M.A.T.; methodology, M.K.S. and M.A.T.; data curation, M.K.S., and A. S.; data analysis, M.K.S., A.S., and S.A.S.; writing—original draft preparation, M.K.S.; writing—review and editing, M.A.T., A.S., S.A.S., S.S.V., C.H.D., E.R.F., and S.R; funding acquisition, M.A.T. All authors reviewed and approved of the final manuscript.