Next Article in Journal
Peptide YY (PYY) Is Associated with Cardiovascular Risk in Patients with Acute Myocardial Infarction
Next Article in Special Issue
Ectopic Fat Accumulation in Pancreas and Heart
Previous Article in Journal
Lateral Acromioplasty has a Positive Impact on Rotator Cuff Repair in Patients with a Critical Shoulder Angle Greater than 35 Degrees
Previous Article in Special Issue
CD36 and GPR120 Methylation Associates with Orosensory Detection Thresholds for Fat and Bitter in Algerian Young Obese Children
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Endurance Runners with Intramyocellular Lipid Accumulation and High Insulin Sensitivity Have Enhanced Expression of Genes Related to Lipid Metabolism in Muscle

1
Department of Metabolism and Endocrinology Juntendo University Graduate School of Medicine, 2-1-1 Hongo, Bunkyo-ku, Tokyo 113-8421, Japan
2
Sportology Center, Juntendo University, Tokyo 113-8421, Japan
3
Institute of Health and Sports Science & Medicine, Juntendo University, Chiba 270-1695, Japan
4
Department of Exercise Physiology, Graduate School of Health and Sports Science, Juntendo University, Chiba 270-1695, Japan
5
Center for Therapeutic Innovations in Diabetes, Juntendo University, Tokyo 113-8421, Japan
6
Center for Identification of Diabetic Therapeutic Targets, Graduate School of Medicine, Juntendo University, Tokyo 113-8421, Japan
*
Author to whom correspondence should be addressed.
J. Clin. Med. 2020, 9(12), 3951; https://doi.org/10.3390/jcm9123951
Submission received: 25 September 2020 / Revised: 20 November 2020 / Accepted: 3 December 2020 / Published: 6 December 2020

Abstract

:
Context: Endurance-trained athletes have high oxidative capacities, enhanced insulin sensitivities, and high intracellular lipid accumulation in muscle. These characteristics are likely due to altered gene expression levels in muscle. Design and setting: We compared intramyocellular lipid (IMCL), insulin sensitivity, and gene expression levels of the muscle in eight nonobese healthy men (control group) and seven male endurance athletes (athlete group). Their IMCL levels were measured by proton-magnetic resonance spectroscopy, and their insulin sensitivity was evaluated by glucose infusion rate (GIR) during a euglycemic–hyperinsulinemic clamp. Gene expression levels in the vastus lateralis were evaluated by quantitative RT-PCR (qRT-PCR) and microarray analysis. Results: IMCL levels in the tibialis anterior muscle were approximately 2.5 times higher in the athlete group compared to the control group, while the IMCL levels in the soleus muscle and GIR were comparable. In the microarray hierarchical clustering analysis, gene expression patterns were not clearly divided into control and athlete groups. In a gene set enrichment analysis with Gene Ontology gene sets, “RESPONSE TO LIPID” was significantly upregulated in the athlete group compared with the control group. Indeed, qRT-PCR analysis revealed that, compared to the control group, the athlete group had 2–3 times higher expressions of proliferator-activated receptor gamma coactivator-1 alpha (PGC1A), adiponectin receptors (AdipoRs), and fatty acid transporters including fatty acid transporter-1, plasma membrane-associated fatty acid binding protein, and lipoprotein lipase. Conclusions: Endurance runners with higher IMCL levels have higher expression levels of genes related to lipid metabolism such as PGC1A, AdipoRs, and fatty acid transporters in muscle.

1. Introduction

Skeletal muscle has been shown to adapt to various stimuli, including exercise. Indeed, endurance training enhances mitochondrial biogenesis, lipid and glucose transport, lipid oxidation, and oxidative capacity in muscle [1,2,3,4,5,6,7]. In addition, endurance-trained athletes who have high oxidative capacity and enhanced insulin sensitivity also have higher intramyocellular (IMCL) content [8,9]. These characteristics of the muscles in endurance-trained athletes are likely acquired as a long-term effect of training. It is supposed that gene expression levels are chronically altered by endurance training over the long term, which contributes to the altered characteristics in the muscles of endurance-trained athletes.
Several studies have compared global gene expression levels in skeletal muscle using microarrays before and after endurance training for 6 to 20 weeks [10,11,12,13]. However, only a few studies have included microarray analysis of muscles in athletes with long-term endurance training [14,15]. Stepto et al. [14] demonstrated that cyclists with long-term (> 5 years) endurance training have higher expressions of gene clusters related to mitochondrial/oxidative capacity than the control subjects. In addition, VO2 peak was correlated with clusters of mitochondrial, fat, and carbohydrate oxidation genes, respectively. Similarly, Wittwer et al. demonstrated that expression levels of genes related to glycolysis were also enhanced in a professional cyclist with 8 years of endurance training [15]. However, the overall gene profile in the skeletal muscle of endurance runners compared with control subjects has not been clarified yet.
Based on this background, we investigated the differences in gene expression profiles in skeletal muscle of endurance runners and control subjects using microarray analysis. In addition, we performed quantitative RT-PCR (qRT-PCR) to specifically measure gene expression levels related to IMCL accumulation and insulin sensitivity, which were clarified in our previous study [16].

2. Research Design and Methods

2.1. Subjects

Study subjects consisted of 7 endurance runners with average records and 8 healthy men with BMI < 23 kg/m2. The endurance runners’ personal records in 5000 m ranged between 14:33.70 and 15:50.00. They were in good health as determined by medical history, physical examination, and standard blood chemistry analyses. All subjects gave written informed consent to the study, which was approved by the Ethics Committee of Juntendo University.

2.2. Study Design and Measurement of Various Parameters

Regular exercise was prohibited from at least 3 days before the experiment day [16]. Each subject consumed a 3-day isocaloric normal-fat diet. Blood samples were taken after overnight fasting. A biopsy of the vastus lateralis muscle was also performed using a needle with local anesthesia. Muscle samples were immediately placed in the RNA-stabilizing reagent RNAlater (Qiagen, Hilden, Germany) and stored at −80 °C. Extraction of total RNA and DNA from muscle samples was performed as previously described [1,16]. Biochemical analyses of serum samples and total body fat content measurements were performed as described previously [1,16].

2.3. Proton Magnetic Resonance Spectroscopy

IMCL in the tibialis anterior muscle (TA) and soleus muscle (SOL) were measured after overnight fasting using 1H-magnetic resonance spectroscopy (VISART 1.5T EX V4.40; Toshiba, Tokyo), as described previously [1,16,17,18,19,20,21]. Voxels (1.2 × 1.2 × 1.2 cm3) were positioned in the TA or SOL avoiding visible interfascial fat and blood vessels. Imaging parameters were set as follows: repetition time of 1500 ms, echo time of 136 ms, acquisition numbers of 192, and 1024 data points over a 1000-kHz spectral width. After examination, the resonances were line fit using a mixed Lorentzian/Gaussian function. After corrections for T1 and T2 relaxations, the methylene signal intensity (S-fat), with peaks observed at approximately 1.25 parts/million (ppm), was quantified. IMCL was quantified by S-fat with the creatine signal at 3.0 ppm (Cre) as the reference and was calculated as a ratio relative to Cre (S-fat/Cre).

2.4. Hyperinsulinemic Euglycemic Clamp Study

A hyperinsulinemic euglycemic glucose clamp study with a target plasma glucose level of 95 mg/dL and an insulin infusion rate of 100 mU/m2/min was performed with an artificial pancreas (STG22; Nikkiso, Shizuoka, Japan), as reported previously [20]. The steady-state glucose infusion rate (GIR) was observed from 105 to 120 min after the beginning of the study and was used as a marker of peripheral insulin sensitivity [1,16,17,18,19,20,21].

2.5. DNA Microarray Analysis

Microarray analysis was performed with muscle samples from 7 control subjects and 5 endurance athletes, as described previously [1]. The transcript profile of each sample was determined using the Affymetrix Human U133 Plus 2.0 array with standard Affymetrix protocols (Affymetrix, Santa Clara, CA, USA). The data were summarized using Genespring GX software (Agilent Technologies, Santa Clara, CA, USA). The signal intensities for the beta actin (ACTB) and glyceraldehyde 3-phosphate dehydrogenase (GAPDH) genes were used as internal quality controls. To correct for variations between gene chips, the signal data (CEL files) were quantile normalized, with probe set intensities calculated using the MAS5 method. Genes identified as being absent across all samples using Genespring GX were excluded from further analyses. The gene expression profiles have been deposited in the Gene Expression Omnibus database under the accession number GSE155271.

2.6. DNA Microarray Data Analysis

Ward’s hierarchical two-way classification method was performed using differentially expressed genes and the subjects. A false discovery rate of less than 5% was considered significant. A cut-off value of 1.3 for fold change_ratios was used. Differentially expressed genes were also identified using gene set enrichment analysis (GSEA) with the Kolmogorov–Smirnov test performed on each sample group. Analysis was performed with Gene Ontology gene sets (c5.all.v2.5) and Reactome gene sets (c2).

2.7. qRT-PCR

Gene expression levels were examined using qRT-PCR (1). Primers were designed using Primer-BLAST.

2.8. Statistical Analysis

All data are expressed as means ± SD. Differences between groups were evaluated using the unpaired t-test. Statistical significance was set at p < 0.05.

3. Results

3.1. Characteristics of the Subjects

Table 1 shows the anthropometric data in the control and athlete groups. Total cholesterol was significantly higher in the athlete group than in the control group. Maximum oxygen uptake was numerically higher in the athlete group, but this difference was not statistically significant. As shown in Figure 1, IMCL levels in TA were significantly higher in the athlete group (p < 0.05). GIR and IMCL levels in SOL were comparable between the groups.

3.2. Hierarchical Clustering Analysis in the Control and Athlete Groups

We used DNA microarrays to profile the gene expression in skeletal muscle in 13 subjects, 5 from the athlete group and 8 from the control group. Hierarchical clustering analysis based on gene expression demonstrated that gene expression patterns were not clearly divided into control and athlete groups, suggesting that differences in gene expression patterns between controls and athletes identified in this study do not surpass differences in gene expression patterns across individuals (Figure 2).

3.3. Gene Expression Profiles of Muscle in the Control and Athlete Groups

We identified 613 genes or probes whose expressions were significantly upregulated or downregulated in athletes compared with controls (p < 0.05, fold change > |1.3|). We identified 286 downregulated genes and 330 upregulated genes in the athlete group compared to the control group (Figure 3A, Supplementary Tables S1 and S2). We did not find genes with more than a 3-fold change. The upregulated gene with the highest fold change (2.38) was nicotinamide riboside kinase 2 (NMRK2). The downregulated gene with the lowest fold change (−2.68) was pyruvate dehydrogenase kinase 4 (PDK4).

3.4. Gene Set Enrichment Analysis (GSEA) of Skeletal Muscle Tissue in the Control and Athlete Groups

In order to gain insight into the biological effects of gene expression patterns, we performed a GSEA with Gene Ontology gene sets (c5.all.v2.5) (Figure 3B,C, Supplementary Tables S4 and S6) and Reactome gene sets (c2) (Supplementary Table S5). This analysis found 70 gene sets significantly upregulated in the control group (p < 0.05) (Figure 3B, Supplementary Table S3) and 13 gene sets significantly upregulated in the athlete group (p < 0.05) (Figure 3C, Supplementary Table S4). Given that one of the important characteristics of endurance-trained athletes is higher IMCL content [8,9], “RESPONSE TO LIPID” (Biological Process) was among the gene sets significantly upregulated in the athlete group compared to the control group (p = 0.045) (Figure 3C, Supplementary Tables S4 and S6).

3.5. Gene Expression Analysis Using qRT-PCR in the Control and Athlete Groups

Previously, we found higher expressions of fatty acid transporters and fatty acid oxidation-related genes in muscle from nonobese nonathletes with high IMCL accumulation and high insulin sensitivity [16]. Here, we found that the “RESPONSE TO LIPID” gene set was significantly upregulated in the athlete group. Thus, we evaluated the expression of genes related to lipid metabolism.
Gene expression levels in fatty acid transporter 1 (FATP1), plasma membrane-associated fatty acid binding protein (FABPpm), and lipoprotein lipase (LPL), which regulate fatty acid uptake in muscle, were approximately 2 times higher in the athlete group than in the control group (p < 0.05) (Figure 4A). Furthermore, gene expression of proliferator-activated receptor alpha (PPARA), a key transcription factor for fatty acid oxidation; proliferator-activated receptor gamma coactivator-1 alpha (PGC1A), a master regulator of mitochondrial biogenesis; hydroxyacyl-CoA dehydrogenase trifunctional multienzyme complex subunit beta (HADHB), a mitochondrial beta oxidation enzyme; and adiponectin receptor 1 (ADIPOR1) and ADIPOR2, which mediate fatty acid oxidation and glucose uptake by adiponectin [22,23] were also significantly higher in the athlete group (Figure 4B).

4. Discussion

In the present study, we compared the gene expression levels in skeletal muscle and characteristics such as IMCL accumulation and insulin sensitivity between endurance runners and control subjects. IMCL levels in TA were approximately 2.5 times higher in the athlete group compared with the control group, while IMCL levels in SOL and GIR were comparable between the groups. In the microarray hierarchical clustering analysis, gene expression patterns were not clearly divided into control and athlete groups. In the GSEA with Gene Ontology gene sets, “RESPONSE TO LIPID” was significantly upregulated in the athlete group compared to that in the control group. The qRT-PCR analysis revealed that the athlete group had expression levels of several genes related to fatty acid transportation, including FATP1, FABPpm, and LPL, that were approximately 2 times higher than in the control group. In addition, genes related to fatty acid oxidation such as PPARA, PGC1A, HADHB, ADIPOR1, and ADIPOR2, were also higher in the athlete group.
In the present study, IMCL levels in TA were approximately 2.5 times higher in endurance athletes than in healthy controls, while IMCL levels in SOL and GIR were comparable between the groups. However, this is not surprising since IMCL is not always raised in endurance athletes [24,25]. A previous study also suggested that IMCL levels in TA are approximately 2 times higher in endurance runners compared with sprinters, while IMCL levels in SOL was comparable between the groups [21]. SOL is a muscle that is rich in type I fibers and has higher IMCL levels, mitochondria, and capillary density [26,27,28,29]; these characteristics are acquired with endurance training [28]. In fact, IMCL in endurance athletes was raised due to higher proportions of type I fiber [30]. On the other hand, TA contains a lower proportion of type I fibers than SOL [31]. Thus, it is possible that endurance training might induce IMCL accumulation in TA, but to a lesser extent in SOL.
We analyzed the gene expression profiles in skeletal muscle of endurance runners using microarray analysis. A previous study investigating the skeletal muscle gene profile of endurance cyclists with microarray analysis showed enhanced expression of gene clusters related to mitochondrial and oxidative capacity [14]. In the present study, as shown in the hierarchical cluster analysis, gene expression was not clearly divided by controls versus endurance athletes, suggesting that endurance runners are not significantly characterized by gene expression patterns in muscle. It may also be possible that control subjects in the present study had relatively higher maximum oxygen uptake than those of a previous study [14], leading to no clear significant differences in gene expression levels between the groups in the microarray analysis. Nonetheless, in the GSEA, “RESPONSE TO LIPID” was upregulated in the athlete group compared to that in the control group, suggesting that this gene cluster might be important in characterizing endurance runners among healthy young people.
It has been suggested that the oxidative capacity of skeletal muscle may be an important modulator of the association between IMCL accumulation and insulin sensitivity [8]. Previous studies have suggested that fatty acid influx into myocytes increases oxidative capacity in muscle. For example, it has been reported that endurance runners have increased FABPpm and LPL [32,33] expression levels, and free fatty acid uptake in skeletal muscle is enhanced in endurance athletes [34]. In addition, overexpression of FABPpm, FATP-1, and LPL increases fatty acid oxidative capacity and mitochondrial content in muscle [35,36,37]. Moreover, a high-fat diet stimulates mitochondrial fatty acid oxidation and increases muscle expression levels of PGC1A, HADHB, and CPT-1 in rodents [38]. Taken together, higher expression levels of genes including FATP1, FABPpm, and LPL may facilitate IMCL accumulation and promote oxidative capacity in muscle and may characterize the athlete’s paradox.
In the present study, expression levels of PPARA, PGC1A, AdipoRs, and several fatty acid transporters were higher in the athlete group compared with the control group. PPARA is a transcriptional factor activated by fatty acids, leading to the induction of genes involved in fatty acid import and β-oxidation [39,40,41,42]. In addition, in vitro and in vivo studies have demonstrated that PGC1A overexpression results in increased expression of CD36, FABPpm, FATPs, mitochondrial content, and oxidative capacity in muscle [43]. We previously showed that overexpression of FABPpm in C2C12 myotubes results in enhanced expression levels of PPARA, PGC1A, AdipoR1, and AdipoR2, while fatty acid uptake and oxidation were both increased in the myotube [16]. PPARα and PGC-1α activation are partly regulated by adiponectin via AdipoR2 and AdipoR1, respectively [3,22]. Exercise training increases muscle AdipoRs [2,3,4], PGC1A, and PPARA expression [5,6,7]. Thus, it has been suggested that PPARA, PGC1A, AdipoRs, and fatty acid transporters are upregulated by endurance exercise and fatty acid incorporation into myocytes, leading to both higher IMCL levels and higher oxidative capacity in endurance athletes.
One limitation of the present study is that we included only male subjects. It has been reported that fatty acid metabolism in skeletal muscle differs by sex and that females have higher levels of IMCL than males. In addition, the number of subjects in the microarray analysis was small, so further study is required to identify specific features of gene expression patterns in endurance athletes. There was no significant difference in maximum oxygen uptake between the groups. Thus, if we included more unfit subjects as controls, we might have different results in IMCL levels and gene expression profiles. Finally, gene expression levels were evaluated in vastus lateralis muscle, while IMCL levels were measured in TA and SOL. Thus, it is still unclear how IMCL accumulation is associated with muscle gene expression in endurance athletes.
In conclusion, endurance runners had higher IMCL levels in TA but similar IMCL levels in SOL and insulin sensitivity as the heathy controls. These characteristics of endurance runners were associated with elevated expressions of genes related to lipid metabolism, including lipid entry (FATP1, FABPpm, and LPL), lipid oxidation (HADHB), and their upstream regulators (PPARA, PGC1A, AdipoRs) in skeletal muscle.

Supplementary Materials

The following are available online at https://www.mdpi.com/2077-0383/9/12/3951/s1, Table S1: Upregulated genes differentially expressed in athletes (p < 0.05 and 1.3-fold change), Table S2: Downregulated genes differentially expressed in athletes (p < 0.05 and -1.3-fold change), Table S3: Significantly upregulated Gene Ontology gene sets in the control group, Table S4: Significantly upregulated Gene Ontology gene sets in the athlete group, Table S5: Significantly upregulated Reactome Pathway gene sets in the control group, Table S6: Gene expression changes in “RESPONSE TO LIPID”.

Author Contributions

S.K. and Y.T. researched the data and contributed to study design, data collection, interpretation of results, wrote and edited the manuscript. K.T., S.-i.I., Y.O., N.S., T.M., and H.N. participated in data collection, data analysis, and contributed to the discussion. R.K. and H.W. contributed to study design, reviewed, and edited the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Culture, Sports, Science and Technology of Japan, High Technology Research Center Grant and Grants-in-Aid for Scientific Research (B) and Ono Medical Foundation and Diabetes Masters Conference and MIZUNO Sports Foundation.

Acknowledgments

We thank Naoko Daimaru, Eriko Magoshi, and Kiyomi Nakamura for their excellent technical assistance.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kakehi, S.; Tamura, Y.; Takeno, K.; Sakurai, Y.; Kawaguchi, M.; Watanabe, T.; Funayama, T.; Sato, F.; Ikeda, S.; Kanazawa, A.; et al. Increased intramyocellular lipid/impaired insulin sensitivity is associated with altered lipid metabolic genes in muscle of high responders to a high-fat diet. Am. J. Physiol. Endocrinol. Metab. 2016, 310, 32–40. [Google Scholar] [CrossRef] [Green Version]
  2. O’Leary, V.B.; Jorett, A.E.; Marchetti, C.M.; Gonzalez, F.; Phillips, S.A.; Ciaraldi, T.P.; Kirwan, J.P. Enhanced adiponectin multimer ratio and skeletal muscle adiponectin receptor expression following exercise training and diet in older insulin-resistant adults. Am. J. Physiol. 2007, 293, 421–427. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Christiansen, T.; Paulsen, S.K.; Bruun, J.M.; Ploug, T.; Pedersen, S.B.; Richelsen, B. Diet-induced weight loss and exercise alone and in combination enhance the expression of adiponectin receptors in adipose tissue and skeletal muscle, but only diet-induced weight loss enhanced circulating adiponectin. J. Clin. Endocrinol. Metab. 2010, 95, 911–919. [Google Scholar] [CrossRef] [PubMed]
  4. Bluher, M.; Bullen, J.W., Jr.; Lee, J.H.; Kralisch, S.; Fasshauer, M.; Kloting, N.; Niebauer, J.; Schon, M.R.; Williams, C.J.; Mantzoros, C.S. Circulating adiponectin and expression of adiponectin receptors in human skeletal muscle: Associations with metabolic parameters and insulin resistance and regulation by physical training. J. Clin. Endocrinol. Metab. 2006, 91, 2310–2316. [Google Scholar] [CrossRef] [Green Version]
  5. Mahoney, D.J.; Parise, G.; Melov, S.; Safdar, A.; Tarnopolsky, M.A. Analysis of global mRNA expression in human skeletal muscle during recovery from endurance exercise. FASEB J. 2005, 19, 1498–1500. [Google Scholar] [CrossRef]
  6. Norrbom, J.; Sundberg, C.J.; Ameln, H.; Kraus, W.E.; Jansson, E.; Gustafsson, T. PGC-1alpha mRNA expression is influenced by metabolic perturbation in exercising human skeletal muscle. J. Appl. Physiol. 2004, 96, 189–194. [Google Scholar] [CrossRef] [Green Version]
  7. Pilegaard, H.; Saltin, B.; Neufer, P.D. Exercise induces transient transcriptional activation of the PGC-1alpha gene in human skeletal muscle. J. Physiol. 2003, 546, 851–858. [Google Scholar] [CrossRef]
  8. Goodpaster, B.H.; He, J.; Watkins, S.; Kelley, D.E. Skeletal muscle lipid content and insulin resistance: Evidence for a paradox in endurance-trained athletes. J. Clin. Endocrinol. Metab. 2001, 86, 5755–5761. [Google Scholar] [CrossRef]
  9. Van Loon, L.J.; Goodpaster, B.H. Increased intramuscular lipid storage in the insulin-resistant and endurance-trained state. Pflugers Arch. 2006, 451, 606–616. [Google Scholar] [CrossRef]
  10. Teran-Garcia, M.; Rankinen, T.; Koza, R.A.; Rao, D.C.; Bouchard, C. Endurance training-induced changes in insulin sensitivity and gene expression. Am. J. Physiol. Endocrinol. Metab. 2005, 288, 1168–1178. [Google Scholar] [CrossRef] [PubMed]
  11. Timmons, J.A.; Larsson, O.; Jansson, E.; Fischer, H.; Gustafsson, T.; Greenhaff, P.L.; Ridden, J.; Rachman, J.; Peyrard-Janvid, M.; Wahlestedt, C.; et al. Human muscle gene expression responses to endurance training provide a novel perspective on Duchenne muscular dystrophy. FASEB J. 2005, 19, 750–760. [Google Scholar] [CrossRef] [PubMed]
  12. Timmons, J.A.; Knudsen, S.; Rankinen, T.; Koch, L.G.; Sarzynski, M.; Jensen, T.; Keller, P.; Scheele, C.; Vollaard, N.B.; Nielsen, S.; et al. Using molecular classification to predict gains in maximal aerobic capacity following endurance exercise training in humans. J. Appl. Physiol. 2010, 108, 1487–1496. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Busso, T.; Fluck, M. A mixed-effects model of the dynamic response of muscle gene transcript expression to endurance exercise. Eur. J. Appl. Physiol. 2013, 113, 1279–1290. [Google Scholar] [CrossRef] [PubMed]
  14. Stepto, N.K.; Coffey, V.G.; Carey, A.L.; Ponnampalam, A.P.; Canny, B.J.; Powell, D.; Hawley, J.A. Global gene expression in skeletal muscle from well-trained strength and endurance athletes. Med. Sci. Sports Exerc. 2009, 41, 546–565. [Google Scholar] [CrossRef] [PubMed]
  15. Wittwer, M.; Billeter, R.; Hoppeler, H.; Fluck, M. Regulatory gene expression in skeletal muscle of highly endurance-trained humans. Acta Physiol. Scand. 2004, 180, 217–227. [Google Scholar] [CrossRef]
  16. Kawaguchi, M.; Tamura, Y.; Kakehi, S.; Takeno, K.; Sakurai, Y.; Watanabe, T.; Funayama, T.; Sato, F.; Ikeda, S.; Ogura, Y.; et al. Association between expression of FABPpm in skeletal muscle and insulin sensitivity in intramyocellular lipid-accumulated nonobese men. J. Clin. Endocrinol. Metab. 2014, 99, 3343–3352. [Google Scholar] [CrossRef] [Green Version]
  17. Sakurai, Y.; Tamura, Y.; Takeno, K.; Kumashiro, N.; Sato, F.; Kakehi, S.; Ikeda, S.; Ogura, Y.; Saga, N.; Naito, H.; et al. Determinants of intramyocellular lipid accumulation after dietary fat loading in non-obese men. J. Diabetes Investig. 2011, 2, 310–317. [Google Scholar] [CrossRef]
  18. Sato, F.; Tamura, Y.; Watada, H.; Kumashiro, N.; Igarashi, Y.; Uchino, H.; Maehara, T.; Kyogoku, S.; Sunayama, S.; Sato, H.; et al. Effects of diet-induced moderate weight reduction on intrahepatic and intramyocellular triglycerides and glucose metabolism in obese subjects. J. Clin. Endocrinol. Metab. 2007, 92, 3326–3329. [Google Scholar] [CrossRef] [Green Version]
  19. Takeno, K.; Tamura, Y.; Kawaguchi, M.; Kakehi, S.; Watanabe, T.; Funayama, T.; Furukawa, Y.; Kaga, H.; Yamamoto, R.; Kim, M.; et al. Relation between insulin sensitivity and metabolic abnormalities in japanese men with BMI of 23–25 kg/m2. J. Clin. Endocrinol. Metab. 2016, 101, 3676–3684. [Google Scholar] [CrossRef] [Green Version]
  20. Tamura, Y.; Tanaka, Y.; Sato, F.; Choi, J.B.; Watada, H.; Niwa, M.; Kinoshita, J.; Ooka, A.; Kumashiro, N.; Igarashi, Y.; et al. Effects of diet and exercise on muscle and liver intracellular lipid contents and insulin sensitivity in type 2 diabetic patients. J. Clin. Endocrinol. Metab. 2005, 90, 3191–3196. [Google Scholar] [CrossRef] [Green Version]
  21. Tamura, Y.; Watada, H.; Igarashi, Y.; Nomiyama, T.; Onishi, T.; Takahashi, K.; Doi, S.; Katamoto, S.; Hirose, T.; Tanaka, Y.; et al. Short-term effects of dietary fat on intramyocellular lipid in sprinters and endurance runners. Metabolism 2008, 57, 373–379. [Google Scholar] [CrossRef] [PubMed]
  22. Kadowaki, T.; Yamauchi, T. Adiponectin and adiponectin receptors. Endocr. Rev. 2005, 26, 439–451. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Yamauchi, T.; Kamon, J.; Ito, Y.; Tsuchida, A.; Yokomizo, T.; Kita, S.; Sugiyama, T.; Miyagishi, M.; Hara, K.; Tsunoda, M.; et al. Cloning of adiponectin receptors that mediate antidiabetic metabolic effects. Nature 2003, 423, 762–769. [Google Scholar] [CrossRef]
  24. Bergman, B.C.; Perreault, L.; Strauss, A.; Bacon, S.; Kerege, A.; Harrison, K.; Brozinick, J.T.; Hunerdosse, D.M.; Playdon, M.C.; Holmes, W.; et al. Intramuscular triglyceride synthesis: Importance in muscle lipid partitioning in humans. Am. J. Physiol. Endocrinol. Metab. 2018, 314, 152–164. [Google Scholar] [CrossRef] [PubMed]
  25. Bruce, C.R.; Anderson, M.J.; Carey, A.L.; Newman, D.G.; Bonen, A.; Kriketos, A.D.; Cooney, G.J.; Hawley, J.A. Muscle oxidative capacity is a better predictor of insulin sensitivity than lipid status. J. Clin. Endocrinol. Metab. 2003, 88, 5444–5451. [Google Scholar] [CrossRef] [Green Version]
  26. Turban, S.; Hajduch, E. Protein kinase C isoforms: Mediators of reactive lipid metabolites in the development of insulin resistance. FEBS Lett. 2011, 585, 269–274. [Google Scholar] [CrossRef] [Green Version]
  27. Kano, Y.; Shimegi, S.; Furukawa, H.; Matsudo, H.; Mizuta, T. Effects of aging on capillary number and luminal size in rat soleus and plantaris muscles. J. Gerontol. A Biol. Sci. Med. Sci. 2002, 57, 422–427. [Google Scholar] [CrossRef] [Green Version]
  28. Poole, D.C.; Mathieu-Costello, O. Relationship between fiber capillarization and mitochondrial volume density in control and trained rat soleus and plantaris muscles. Microcirculation 1996, 3, 175–186. [Google Scholar] [CrossRef]
  29. Sullivan, S.M.; Pittman, R.N. Relationship between mitochondrial volume density and capillarity in hamster muscles. Am. J. Physiol. 1987, 252, 149–155. [Google Scholar] [CrossRef] [Green Version]
  30. Van Loon, L.J.; Koopman, R.; Manders, R.; van der Weegen, W.; van Kranenburg, G.P.; Keizer, H.A. Intramyocellular lipid content in type 2 diabetes patients compared with overweight sedentary men and highly trained endurance athletes. Am. J. Physiol. Endocrinol. Metab. 2004, 287, 558–565. [Google Scholar] [CrossRef]
  31. Polgar, J.; Johnson, M.A.; Weightman, D.; Appleton, D. Data on fibre size in thirty-six human muscles. An autopsy study. J. Neurol. Sci. 1973, 19, 307–318. [Google Scholar] [CrossRef]
  32. Kiens, B.; Roepstorff, C.; Glatz, J.F.; Bonen, A.; Schjerling, P.; Knudsen, J.; Nielsen, J.N. Lipid-binding proteins and lipoprotein lipase activity in human skeletal muscle: Influence of physical activity and gender. J. Appl. Physiol. 2004, 97, 1209–1218. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Talanian, J.L.; Holloway, G.P.; Snook, L.A.; Heigenhauser, G.J.; Bonen, A.; Spriet, L.L. Exercise training increases sarcolemmal and mitochondrial fatty acid transport proteins in human skeletal muscle. Am. J. Physiol. 2010, 299, 180–188. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Turcotte, L.P.; Richter, E.A.; Kiens, B. Increased plasma FFA uptake and oxidation during prolonged exercise in trained vs. untrained humans. Am. J. Physiol. 1992, 262, 791–799. [Google Scholar] [CrossRef]
  35. Nickerson, J.G.; Alkhateeb, H.; Benton, C.R.; Lally, J.; Nickerson, J.; Han, X.X.; Wilson, M.H.; Jain, S.S.; Snook, L.A.; Glatz, J.F.; et al. Greater transport efficiencies of the membrane fatty acid transporters FAT/CD36 and FATP4 compared with FABPpm and FATP1 and differential effects on fatty acid esterification and oxidation in rat skeletal muscle. J. Biol. Chem. 2009, 284, 16522–16530. [Google Scholar] [CrossRef] [Green Version]
  36. Holloway, G.P.; Chou, C.J.; Lally, J.; Stellingwerff, T.; Maher, A.C.; Gavrilova, O.; Haluzik, M.; Alkhateeb, H.; Reitman, M.L.; Bonen, A. Increasing skeletal muscle fatty acid transport protein 1 (FATP1) targets fatty acids to oxidation and does not predispose mice to diet-induced insulin resistance. Diabetologia 2011, 54, 1457–1467. [Google Scholar] [CrossRef] [Green Version]
  37. Levak-Frank, S.; Radner, H.; Walsh, A.; Stollberger, R.; Knipping, G.; Hoefler, G.; Sattler, W.; Weinstock, P.H.; Breslow, J.L.; Zechner, R. Muscle-specific overexpression of lipoprotein lipase causes a severe myopathy characterized by proliferation of mitochondria and peroxisomes in transgenic mice. J. Clin. Invest. 1995, 96, 976–986. [Google Scholar] [CrossRef] [Green Version]
  38. Turner, N.; Bruce, C.R.; Beale, S.M.; Hoehn, K.L.; So, T.; Rolph, M.S.; Cooney, G.J. Excess lipid availability increases mitochondrial fatty acid oxidative capacity in muscle: Evidence against a role for reduced fatty acid oxidation in lipid-induced insulin resistance in rodents. Diabetes 2007, 56, 2085–2092. [Google Scholar] [CrossRef] [Green Version]
  39. Vega, R.B.; Huss, J.M.; Kelly, D.P. The coactivator PGC-1 cooperates with peroxisome proliferator-activated receptor alpha in transcriptional control of nuclear genes encoding mitochondrial fatty acid oxidation enzymes. Mol. Cell. Biol. 2000, 20, 1868–1876. [Google Scholar] [CrossRef] [Green Version]
  40. Duncan, J.G. Peroxisome proliferator activated receptor-alpha (PPARalpha) and PPAR gamma coactivator-1alpha (PGC-1alpha) regulation of cardiac metabolism in diabetes. Pediatr. Cardiol. 2011, 32, 323–328. [Google Scholar] [CrossRef] [Green Version]
  41. Mandard, S.; Muller, M.; Kersten, S. Peroxisome proliferator-activated receptor alpha target genes. Cell. Mol. Life Sci. 2004, 61, 393–416. [Google Scholar] [CrossRef] [PubMed]
  42. Mandard, S.; Patsouris, D. Nuclear control of the inflammatory response in mammals by peroxisome proliferator-activated receptors. PPAR Res. 2013, 2013, 613864. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Supruniuk, E.; Miklosz, A.; Chabowski, A. The implication of PGC-1alpha on fatty acid transport across plasma and mitochondrial membranes in the insulin sensitive tissues. Front. Physiol. 2017, 8, 923. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Absolute intramyocellular lipid (IMCL) and glucose infusion rate (GIR) levels in the control and athlete groups. Absolute values of GIR and IMCL in the control group (Con) and athlete group (Ath). The average in each group is shown as closed circle. * p < 0.05 Intramyocellular lipid, IMCL; tibialis anterior, TA; soleus, SOL; methylene signal intensity, S-fat; creatine signal, Cre.
Figure 1. Absolute intramyocellular lipid (IMCL) and glucose infusion rate (GIR) levels in the control and athlete groups. Absolute values of GIR and IMCL in the control group (Con) and athlete group (Ath). The average in each group is shown as closed circle. * p < 0.05 Intramyocellular lipid, IMCL; tibialis anterior, TA; soleus, SOL; methylene signal intensity, S-fat; creatine signal, Cre.
Jcm 09 03951 g001
Figure 2. Hierarchical clustering of microarray data in the control and athlete groups. Euclidean distance was used as the metric. The tree was constructed using Ward’s method. Ath1–Ath5, Endurance athletes (n = 5); Con1–Con8 Control subjects (n = 8).
Figure 2. Hierarchical clustering of microarray data in the control and athlete groups. Euclidean distance was used as the metric. The tree was constructed using Ward’s method. Ath1–Ath5, Endurance athletes (n = 5); Con1–Con8 Control subjects (n = 8).
Jcm 09 03951 g002
Figure 3. Microarray findings in the control and athlete groups. (A) Volcano plot comparing gene expression in the control versus athlete group. The x-axis indicates differential expression profiles plotted in fold-induction ratios on the log2 scale. The y-axis indicates the statistical significance of the difference in expression on the log10 scale. Genes upregulated in the athlete group are shown in red. Genes downregulated in the athlete group are shown in blue. (B) The top 20 significantly upregulated Gene Ontology gene sets (c5.all.v2.5) based on GSEA in the control group. (C) The top 20 significantly upregulated Gene Ontology gene sets (c5.all.v2.5) based on GSEA in the athlete group. Gene set enrichment analysis; GSEA, Normalized enrichment score; NES.
Figure 3. Microarray findings in the control and athlete groups. (A) Volcano plot comparing gene expression in the control versus athlete group. The x-axis indicates differential expression profiles plotted in fold-induction ratios on the log2 scale. The y-axis indicates the statistical significance of the difference in expression on the log10 scale. Genes upregulated in the athlete group are shown in red. Genes downregulated in the athlete group are shown in blue. (B) The top 20 significantly upregulated Gene Ontology gene sets (c5.all.v2.5) based on GSEA in the control group. (C) The top 20 significantly upregulated Gene Ontology gene sets (c5.all.v2.5) based on GSEA in the athlete group. Gene set enrichment analysis; GSEA, Normalized enrichment score; NES.
Jcm 09 03951 g003aJcm 09 03951 g003b
Figure 4. Gene expression analysis of skeletal muscle in the control and athlete groups using qRT-PCR. (A) Expression levels of fatty acid transporter in the control group (white bars) and the athlete group (black bars). Data are relative to the expression level in the control group, which was set to 1. (B) Expression levels of genes related to fatty acid β-oxidation in the control group (white bars) and the athlete group (black bars). Data are relative to the expression level in the control group, which was set to 1. * p < 0.05, ** p < 0.01. Values were obtained by normalization to a housekeeping gene (ACTB). Data are presented as means ± SD. CD36: CD36, FATP1: fatty acid transporter protein 1, FABPpm: plasma membrane-associated fatty acid binding protein, LPL: lipoprotein lipase, PPARA: peroxisome proliferator-activated receptor-α, PGC1A: proliferator-activated receptor gamma coactivator-1α, PDHA: pyruvate dehydrogenase-α; HADHB: hydroxyacyl-CoA dehydrogenase-β, SDHB succinate dehydrogenase subunit B, CPT1: carnitine palmitoyltransferase 1, ACSL: long chain acyl-CoA synthetase, ADIPOR1: adiponectin receptor 1, ADIPOR2: adiponectin receptor 2.
Figure 4. Gene expression analysis of skeletal muscle in the control and athlete groups using qRT-PCR. (A) Expression levels of fatty acid transporter in the control group (white bars) and the athlete group (black bars). Data are relative to the expression level in the control group, which was set to 1. (B) Expression levels of genes related to fatty acid β-oxidation in the control group (white bars) and the athlete group (black bars). Data are relative to the expression level in the control group, which was set to 1. * p < 0.05, ** p < 0.01. Values were obtained by normalization to a housekeeping gene (ACTB). Data are presented as means ± SD. CD36: CD36, FATP1: fatty acid transporter protein 1, FABPpm: plasma membrane-associated fatty acid binding protein, LPL: lipoprotein lipase, PPARA: peroxisome proliferator-activated receptor-α, PGC1A: proliferator-activated receptor gamma coactivator-1α, PDHA: pyruvate dehydrogenase-α; HADHB: hydroxyacyl-CoA dehydrogenase-β, SDHB succinate dehydrogenase subunit B, CPT1: carnitine palmitoyltransferase 1, ACSL: long chain acyl-CoA synthetase, ADIPOR1: adiponectin receptor 1, ADIPOR2: adiponectin receptor 2.
Jcm 09 03951 g004
Table 1. Characteristics of the control and athlete groups.
Table 1. Characteristics of the control and athlete groups.
Control Group (n = 8)Athlete Group (n = 7)
Age (years)21.7 ± 1.421.6 ± 0.8
Body mass index (kg/m2)19.7 ± 1.120.5 ± 1.1
Body fat (%)12 ± 3.513.3 ± 3.0
Glucose (mg/dL)83.4 ± 5.185.6 ± 7.1
Insulin (μU/mL)2.8 ± 1.04.3 ± 1.2
Free fatty acids (mmol/L)0.43 ± 0.160.49 ± 0.10
Total cholesterol (mg/dL) 169 ± 17.2208.6 ± 33.3 *
High-density lipoprotein cholesterol (mg/dL)59.4 ± 11.459.4 ± 10.9
Triglycerides (mg/dL) 62.2 ± 31.065.2 ± 20.3
Hemoglobin A1c (%)4.7 ± 0.24.7 ± 0.2
High molecular weight adiponectin (µg/mL)1.81 ± 1.151.89 ± 1.20
Maximum oxygen uptake (mL/kg/min)52.2 ± 6.958.4 ± 3.9
Data are presented as means ± SD. * p < 0.05.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kakehi, S.; Tamura, Y.; Takeno, K.; Ikeda, S.-i.; Ogura, Y.; Saga, N.; Miyatsuka, T.; Naito, H.; Kawamori, R.; Watada, H. Endurance Runners with Intramyocellular Lipid Accumulation and High Insulin Sensitivity Have Enhanced Expression of Genes Related to Lipid Metabolism in Muscle. J. Clin. Med. 2020, 9, 3951. https://doi.org/10.3390/jcm9123951

AMA Style

Kakehi S, Tamura Y, Takeno K, Ikeda S-i, Ogura Y, Saga N, Miyatsuka T, Naito H, Kawamori R, Watada H. Endurance Runners with Intramyocellular Lipid Accumulation and High Insulin Sensitivity Have Enhanced Expression of Genes Related to Lipid Metabolism in Muscle. Journal of Clinical Medicine. 2020; 9(12):3951. https://doi.org/10.3390/jcm9123951

Chicago/Turabian Style

Kakehi, Saori, Yoshifumi Tamura, Kageumi Takeno, Shin-ichi Ikeda, Yuji Ogura, Norio Saga, Takeshi Miyatsuka, Hisashi Naito, Ryuzo Kawamori, and Hirotaka Watada. 2020. "Endurance Runners with Intramyocellular Lipid Accumulation and High Insulin Sensitivity Have Enhanced Expression of Genes Related to Lipid Metabolism in Muscle" Journal of Clinical Medicine 9, no. 12: 3951. https://doi.org/10.3390/jcm9123951

APA Style

Kakehi, S., Tamura, Y., Takeno, K., Ikeda, S. -i., Ogura, Y., Saga, N., Miyatsuka, T., Naito, H., Kawamori, R., & Watada, H. (2020). Endurance Runners with Intramyocellular Lipid Accumulation and High Insulin Sensitivity Have Enhanced Expression of Genes Related to Lipid Metabolism in Muscle. Journal of Clinical Medicine, 9(12), 3951. https://doi.org/10.3390/jcm9123951

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop