- 1School of Physical Education, China University of Geosciences (Wuhan), Wuhan, China
- 2Institute of Hydrobiology, Chinese Academy of Sciences, Wuhan, China
This study aims to explore the molecular regulatory mechanisms of acute exercise in the skeletal muscle of mice. Male C57BL/6 mice were randomly assigned to the control group, and the exercise group, which were sacrificed immediately after an acute bout of exercise. The study was conducted to investigate the metabolic and transcriptional profiling in the quadriceps muscles of mice. The results demonstrated the identification of 34 differentially expressed metabolites (DEMs), with 28 upregulated and 6 downregulated, between the two groups. Metabolic pathway analysis revealed that these DEMs were primarily enriched in several, including the citrate cycle, propanoate metabolism, and lysine degradation pathways. In addition, the results showed a total of 245 differentially expressed genes (DEGs), with 155 genes upregulated and 90 genes downregulated. KEGG analysis indicated that these DEGs were mainly enriched in various pathways such as ubiquitin mediated proteolysis and FoxO signaling pathway. Furthermore, the analysis revealed significant enrichment of DEMs and DEGs in signaling pathways such as protein digestion and absorption, ferroptosis signaling pathway. In summary, the identified multiple metabolic pathways and signaling pathways were involved in the exercise-induced physiological regulation of skeletal muscle, such as the TCA cycle, oxidative phosphorylation, protein digestion and absorption, the FoxO signaling pathway, ubiquitin mediated proteolysis, ferroptosis signaling pathway, and the upregulation of KLF-15, FoxO1, MAFbx, and MuRF1 expression could play a critical role in enhancing skeletal muscle proteolysis.
1 Introduction
Skeletal muscle, being the most abundant tissue in the body, is indispensable for movement as the body relies on the contraction and relaxation of skeletal muscle. Furthermore, a key predictor of mortality is muscle mass quality (Baskin et al., 2015). The skeletal muscle plays a crucial role in maintaining the balance of metabolic processes, making it one of the primary sites for glucose disposal (DeFronzo and Tripathy, 2009). Unhealthy lifestyles and habits have led to an alarming increase in the incidence of various chronic diseases, imposing a significant economic burden on healthcare systems worldwide (Kim et al., 2016; Farr and Khosla, 2019; Castellanos et al., 2020). It is well known that exercise, as a stressor, has a significant impact on physical function (Cook et al., 2012; Wang et al., 2021). Moreover, exercise is a major factor in the extensive metabolic and molecular remodeling of skeletal muscles (Egan and Zierath, 2013). The skeletal muscle undergoes various adaptations in response to exercise, including the synchronization of muscle contraction and ATP production, as well as the utilization of energy triggered by mechano- and other metabolic sensors (Egan and Zierath, 2013; Camera et al., 2016; Fan and Evans, 2017).
Studies have shown that regular physical activity enhances exercise performance, induces increased skeletal muscle protein synthesis, inhibits skeletal muscle protein degradation (Qiu et al., 2014), improves mitochondrial function, and subsequently promotes metabolic homeostasis in the body (Shen et al., 2021). However, intensive or prolonged exercise can trigger an oxidative imbalance, which can cause oxidative damage and lead to muscle fatigue (Reid et al., 1992; O'Neill et al., 1996). If the ability to remodel skeletal muscles decreases and their quality deteriorates, it can affect the bodyʼs motor function and increase the risk of falls and injuries. The enhancement of muscle function is achieved through improved the exercise capacity and the optimized metabolic functions. Despite significant efforts to uncover the comprehensive changes in muscle caused by exercise, there are still many regulators of skeletal muscle that have yet to be discovered.
Transcriptomic profiling is a powerful technique that allows for the identification of intricate gene expression patterns (Wang et al., 2009). On the other hand, metabolomics serves as a valuable tool in detecting and elucidating changes in metabolites within a biological system under varying conditions (Dettmer et al., 2007), but few studies have investigated the comprehensive transcriptomic and metabolomic profile of exercise muscular physiology. The objective of this study was to evaluate and analyze the transcriptional and metabolic networks of skeletal muscle that are regulated by exercise in mice. By doing so, we were able to gain insights into the mechanisms through which exercise intervention impacts the overall metabolic balance of skeletal muscle.
2 Material and methods
2.1 Animals
The study used C57BL/6 male mice were 9 weeks old and weighing 22.1 ± 1.3 g. The mice utilized in this study were housed in the animal facilities located at Hubei Provincial Center for Disease Prevention and Control in China, where they received a standard photoperiod of 12-h light/12-h dark cycles. Mice were housed under standard conditions, with access to food and water unrestricted, at a temperature of 21°C ± 1°C. After 1 week of the acclimation period, the animals were randomly assigned to the control or exercise groups. All animal experiments were approved by the local Animal Care and Use Committee and performed following the recognized guidelines for laboratory animal care and management. Ethical approval was granted by the Academic Integrity and Scientific Research Ethics Committee of China University of Geosciences (Wuhan).
2.2 Aerobic capacity test and treadmill exercise protocol
This study utilized a mouse treadmill, following the method previously mentioned by Kruger et al. (2009) and Kemi et al. (2002). To minimize the stress caused by the equipment, the animals were provided with 1 week of exercise adaptation prior to undergoing the exercise test. Each animal underwent an incremental exercise test to exhaustion at least 4 days prior to commencing the experiments to measure their maximal oxygen uptake (
2.3 Tissue samples
After conducting the exercise tests, the animals were administered anesthesia, and then euthanized using cervical dislocation. The muscle samples used in this study referenced the experimental conditions described previously (Song et al., 2023). The posterior thigh of mice was sterilized, and a skin incision was made. The m. gluteus superficialis and the tensor fasciae latae muscle were identified. A sample was taken from the exposed quadriceps femoris muscle. The quadriceps were weighed; mass was normalized to body mass. To maintain its integrity, the specimen was promptly cooled in a pre-chilled diethylpyrocarbonate (DEPC) solution at a temperature of 4°C, effectively eliminating any superficial blood and hair. Then, the muscular tissue was desiccated using filter paper, carefully placed into a cryopreservation tube, and promptly placed in liquid nitrogen for subsequent examination.
2.4 Experimental methods of metabolomics
Retrieve the sample from liquid nitrogen and place it on ice to thaw until it reaches a state where it can be easily cut (all subsequent steps were performed on ice). Chop and mix the samples, and weigh the samples at 20 mg (±1 mg) into the corresponding numbered centrifuge tubes. Add a steel ball using forceps, homogenize for 20 s with a ball mill (30 HZ), and centrifuge the sample at 3,000 r/min for 30 s at 4°C (the time can be increased according to the actual homogenization). After centrifugation, add 400 μL of 70% methanolic water internal standard extract, centrifuged at 2,500 r/min for 5 min, and let stand on ice for 15 min. Centrifuge at 12,000 r/min for 10 min at 4°C, transfer 300 μL of the supernatant to another centrifuge tube with a corresponding number, and let it stand in a −20°C refrigerator for 30 min; Under 4°C conditions, centrifuge again at 12,000 r/min for 3 min, transfer 200 μL of the supernatant into a lined tube inside the corresponding sample bottle for analysis on the machine. The samples were then analysed using ultra-performance liquid chromatography-tandem mass spectrometry (UPLC-MS/MS). Mass spectroscopy data processing was performed using Analyst software (version 1.6.3). Orthogonal projections to latent structures-discriminant analysis (OPLS-DA) were employed to analyze the differentially expressed metabolites (DEMs). DEMs were screened based on a combination of fold change, t-test p-value, and variable importance in the projection (VIP) values from the OPLS-DA model. The screening criteria were set as VIP >1 and p-value <0.05. Use R language to perform OPLS-DA to determine the magnitude of differences between samples. Subsequently, metabolic pathway enrichment analysis of DEMs was performed.
2.5 RNA-seq experimental method
Total RNA was isolated from skeletal muscle using the RNeasy Mini Kit (Qiagen, Hilden, Germany), and the quantity and purity were determined with a UV-Vis spectrophotometer (ND-1000, Nano-Drop Technologies, Rockland, United States). During RNA separation, the RNase-free DNase set (Qiagen, Hilden, Germany) was used to perform on-column DNase digestion, ensuring the elimination of any potential DNA contamination. The transcriptomic sequencing (RNA-seq) was performed using Illumina Novaseq 6000 (Illumina, CA) to construct and sequence libraries. The high-quality data were aligned to a transcript reference obtained from the mouse reference genome (version mm10), through Bowtie2 default parameters (Langmead and Salzberg, 2012). Transcript abundance was normalized and quantified utilizing RSEM (RNA-Seq by Expectation Maximization) (Li and Dewey, 2011). Differentially expressed genes (DEGs) were detected using a threshold of Fold Change ≥2 and FDR <0.01. Subsequently, the obtained DEGs were subjected to Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) functional annotation and enrichment analysis.
2.6 Validation of DEGs by quantitative real-time PCR
To validate the reliability of RNA-seq data, a subset of 5 differentially expressed genes (DEGs) was selected for quantitative real-time PCR (qPCR) using the primers provided in Table 1. The High-capacity cDNA reverse transcription kit (Applied Biosystems) was used according to the manufacturerʼs instructions to synthesize cDNA from the RNA samples. The cDNA was obtained using the T100 PCR thermal cycler (Bio-Rad Laboratories, Munich, Germany) and used for qPCR. The iQ SYBR Green Supermix (Bio-Rad Laboratories, Munich, Germany) was utilized for qPCR to measure the mRNA expression. The amplification products were quantified using an iCycler (Bio-Rad Laboratories, Munich, Germany). The reaction conditions included an initial denaturation at 95°C for 3 min, followed by 42 cycles of denaturation at 95°C for 15 s, annealing at 61°C for 30 s, and elongation at 72°C for 30 s. The mRNA expression levels were normalized to the expression of the housekeeping gene β-actin. The results were analyzed using the 2−Δ(ΔCt) method and presented as fold changes relative to the control group.
2.7 Combined analysis of transcriptome and metabolome
Conjoint Analysis of DEGs and DEMs is conducted based on the relative content data from transcriptomics and metabolomics. Enrichment analysis is performed on the significantly enriched KEGG pathways for both DEGs and DEMs. This analysis identifies the pathways that are jointly enriched by DEGs and DEMs, and these pathways are plotted as bar graphs.
2.8 Statistical analysis
This study presents the mean values and standard errors of the mean, which were determined. The data was tested for normality using the Shapiro-Wilk test, and the results indicated that the data followed a normal distribution. A one-way ANOVA test was conducted to analyze the differences between the groups, and the results were considered statistically significant at p < 0.05. All statistical analyses were performed using SPSS 26.0 software.
3 Results
3.1 Screening for differential metabolites
The OPLS-DA model showed that the two groups of samples could be separated, indicating that the intervention of an acute bout of exercise caused changes in skeletal muscle metabolism in rats (Figure 1). The T score was 20.9% and the orthogonal T score was 30.3% indicating significant differences in metabolites between the experimental group and the control group. The fold change was combined with the VIP values of the OPLS-DA model to filter the DEMs. In the exercise group and the control group, a total of 34 DEMs were identified, comprising 28 upregulated metabolites and 6 downregulated metabolites (Figure 2). These include 15 amino acids and its metabolites, 8 organic acids and its derivatives, 3 nucleotide and its metabolites, 3 fatty acyls, 2 carbohydrates and its metabolites, 1 heterocyclic compound, 1 coenzyme and vitamin, and 1 aldehyde, ketone, and ester.
FIGURE 1. Score plot of the OPLS-DA. The abscissa indicates the predicted principal component; the ordinate indicates the orthogonal principal component; the percentage indicates the explanation rate of this component to the data set. Each point in the graph indicates a sample, and samples of the same group are represented using the same color.
FIGURE 2. Screening of Differentially Expressed Metabolites. (A): Volcano plot of differentially expressed metabolites. Each point in the volcano plot represents a metabolite, with green dots representing downregulated differential metabolites, red dots representing upregulated differential metabolites, and gray representing metabolites that are detected but not significantly different. The abscissa represents the logarithmic value (log2FC) of the relative abundance difference of a metabolite between the two groups of samples. The larger the absolute value of the abscissa, the greater the relative abundance difference of the metabolite between the two groups of samples. Under the filtering conditions of VIP + FC + p-value, the ordinate represents the significance level of the difference (-log10p-value), and the size of the dots represents the VIP value. (B): Clustering heat map of the metabolites. The abscissa is used to display the names of samples, and the ordinate on the right is used to display the names of metabolites. Different colors are filled with different values obtained after normalization for different relative contents. Red indicates high content; green indicates low content.
3.2 Signaling pathway analysis of differential metabolites
Metabolic pathway analysis was conducted on the DEMs between the exercise group and the control group, and a significant enrichment pathway for metabolites was obtained using a hypergeometric distribution algorithm (p < 0.05) (Figure 3). The top 20 pathways in terms of p-value were selected for display from smallest to largest and made into a bubble chart. The DEMs were mainly enriched in propanoate metabolism, protein digestion and absorption, d-amino acid metabolism, lysine degradation, citrate cycle, alanine, aspartate and glutamate metabolism, ferroptosis, pyruvate metabolism and aminoacyl-tRNA biosynthesis.
FIGURE 3. Pathway analysis of differentially expressed metabolites. The abscissa indicates the Rich Factor for each pathway, the ordinate is the pathway name (sorted by p-value), and the color of the dot reflects the p-value size, the redder the enrichment is, the more significant the enrichment is. Dot size relates to the number of differentially expressed metabolites.
3.3 Screening of differential gene expression
The raw data were obtained from 6 sequenced samples of skeletal muscle total RNA samples after sequencing by the transcriptome Illumina Novaseq 6000 sequencing platform. Using Fold Change≥2 and FDR<0.01 as criteria, we identified 245 DEGs in the exercise group compared with the control group. Among them, 155 upregulated genes and 90 downregulated genes accounted for 63% and 37% of the total DEGs, respectively (Figure 4). Furthermore, we have identified 10 key genes that are highly associated with skeletal muscle (Table 2).
FIGURE 4. Identification of differentially expressed genes (A): Volcano plots of differentially expressed genes. In a volcano plot, each dot represents a gene. Blue dots represent the downregulated genes, red dots represent the upregulated genes, and the gray dots represent the non-significant genes (B): Heatmap of differentially expressed genes.
3.4 GO and KEGG analysis of differentially expressed genes
Based on our screened DEGs, the top 15 GO Terms in q-value ranking from GO functional analysis were selected for presentation and bubble plotted (Figure 5). The DEGs were annotated to skeletal muscle cell differentiation, protein folding chaperone, aging, response to muscle stretch, positive regulation of cell migration and response to heat in the biological process functional classification. The 245 differential genes were subjected to KEGG enrichment analysis, and the results showed that a total of 238 signaling pathways were enriched. The top 15 enriched KEGG signaling pathways were selected for display (Figure 6). These include the ubiquitin mediated proteolysis, longevity regulating pathway, cAMP signaling pathway, FoxO signaling pathway, and MAPK signaling pathway.
FIGURE 5. Gene Ontology enrichment analysis of differentially expressed genes. The abscissa is the ratio of the number of genes to all genes for the term, ordinate shows the GO terms, the size of the bubble represents the number of genes enriched in the GO terms, and color showed the q-value of GO terms.
FIGURE 6. Functional annotation of the KEGG pathway for differentially expressed genes. Ordinate: nomenclature of the KEGG signaling pathway; Abscissa: gene numbers that annotated this pathway and the ratio of KEGG in all genes that annotated relative to the total KEGG.
3.5 Integrative KEGG pathway analysis of the transcriptome and metabolome
To further understand the correlation between metabolites and genes involved in the same biological process, we conducted a comprehensive analysis of transcriptomic and metabolomic data, resulting in the identification of 20 KEGG signaling pathways (Figure 7). The KEGG pathways predominantly enriched in differentially expressed metabolites and genes include protein digestion and absorption, biosynthesis of amino acids, cAMP signaling pathway, protein digestion and absorption and ferroptosis, among others.
FIGURE 7. KEGG pathway enrichment of differential genes and metabolites. The abscissa represents signaling pathways, the ordinate represents enriched p-values, expressed as -log (p-value).
3.6 Quantitative real-time PCR
The results of qPCR analysis in the samples are graphically represented in Figure 8. The levels of mRNA in the exercise mice were analyzed and compared to those of control animals. The transcript levels of MAFbx, Atf3, and Sap30 (p < 0.01 vs. CG, n = 3), KLF-15 and UCP3 (p < 0.05 vs. CG, n = 3) were significantly increased in the exercise group. The results of RNA-Seq are consistent with these findings.
FIGURE 8. The effects of acute exercise on the gene expression of muscle. Compared with the control group, the mRNA expression in the exercise group. The relative levels of target gene mRNA expression were normalized against the mRNA expression of internal housekeeping gene (β-actin). The expression of β-actin mRNA was not significantly different between the tested groups allowing a direct comparison (data not shown). Note that columns and error bars represent (mRNA in % to respective control ±SEM), and levels of significance (*p < 0.05, **p < 0.01, n = 3) are indicated.
4 Discussion
Skeletal muscle, a crucial component in controlling movement, plays a fundamental role as the primary regulator of metabolism in various body systems. Its susceptibility to the metabolic regulations imposed by exercise is particularly notable. To the best of our knowledge, this study represents the first investigation that employs transcriptomics and metabolomics to extensively explore the collective effects of acute exercise on skeletal muscle. This study identified a total of 34 metabolites that exhibited differential regulation. These metabolites are primarily associated with several metabolic pathways, including the citrate cycle (TCA cycle), oxidative phosphorylation, and Valine, leucine and isoleucine degradation. This transcriptomics analysis unveiled 245 differentially expressed genes. The KEGG pathways involved in this study include FoxO, ubiquitin mediated proteolysis signaling pathway. By conducting a joint analysis of transcriptomics and metabolomics, 20 enriched KEGG pathways were identified, including biosynthesis of amino acids, ferroptosis, protein digestion and absorption signaling pathway. Based on the aforementioned findings, we focus on a comprehensive analysis to investigate the potential mechanism by which exercise regulates skeletal muscle function.
To begin with, the results of metabolomics in this study demonstrated a significant increase in succinate levels within the TCA cycle following acute exercise, as illustrated in Figure 3. Succinate, being an integral component of the TCA cycle, can accumulate when the flow of succinic acid within the cycle is hindered. The metabolic dysfunction in the TCA cycle can lead to a deterioration in physiological function. Succinate is generated within the mitochondria and its main destiny upon transition to an aerobic environment is to undergo oxidation at the same location. Consequently, in a scenario where aerobic metabolism significantly outweighs anaerobic metabolism, succinate will not accumulate. However, the accumulation of succinate under aerobic conditions can occur if oxygen cost far exceeds the individual’s maximal oxygen uptake (VO2max) and anaerobic energy production is significant (Hochachka and Dressendorfer, 1976). In the analysis of metabolic pathways, succinate is not only closely related to the TCA cycle, but is also involved in various other metabolic pathways, including oxidative phosphorylation (Zhao et al., 2017).
The TCA cycle is an important pathway for intracellular oxidative phosphorylation to produce ATP, and mitochondria are the main site where the TCA cycle takes place. Moderate-intensity exercise can improve mitochondrial function, whereas both high-intensity and exhaustive exercise may disrupt mitochondrial structure and cause mitochondrial disorders (Tondera et al., 2009; Schoepe et al., 2012). In particular, skeletal muscle is an important site for the uptake and utilization of fatty acids, and mitochondria play an important role in cellular energy metabolism, being the main site for energy conversion, the TCA cycle, and oxidative phosphorylation. The mitochondria undergo pathological changes following acute exercise, leading to impaired function (Magalhaes et al., 2013). As a result, the cellular capacity for oxidative phosphorylation is affected. This implies that the disruption of mitochondria by acute exercise may lead to an inadequate energy supply and adversely affect the organism.
Disruptions in mitochondrial function could be related to amino acid pool imbalances and altered amino acid metabolism (Naseer et al., 2020). The fascinating finding in this study was that strenuous exercise resulted in the degradation of amino acids such as valine, leucine, isoleucine, and so on (as depicted in Figure 3), while amino acids such as glycine and serine increased, demonstrating that the body undergoes degradation and production of proteins associated with the injury zone immediately after the onset of the injury. It is believed that succinate is produced as a result of the catabolism of amino acids. The findings of this study demonstrated that the levels of 15 amino acids and their derivatives, such as l-asparagine anhydrous, l-lysine, d-ornithine, pyroglutamic acid, and n, n-dimethylarginine and so on, were upregulated following an acute bout of intensive exercise. These changes in metabolites involve pathways such as d-amino acid metabolism, pyruvate metabolism, propanoate metabolism, and fatty acid degradation and so on. Besides, a source of anaerobic energy may be generated by amino acid catabolism (Hochachka and Dressendorfer, 1976). The noteworthy elevation of amino acids in vivo following acute exercise suggests that the rate of acid production in the skeletal muscle of exercising mice exceeds the rate of elimination, there is an increase in muscle protein degradation, which results in the release of amino acids from muscle tissue into the bloodstream. This phenomenon could be a compensatory mechanism to provide amino acids as a substrate for energy production and gluconeogenesis during exercise. Exercise could stimulate the activity of key enzymes involved in amino acid catabolism (Harris et al., 2005).
The general consensus is that protein synthesis in skeletal muscle is suppressed during exercise, as evidenced by findings in both humans and other species (Rennie et al., 1980; Dohm et al., 1985). The results of transcriptomics in this study show that acute exercise could induce an increase in KLF-15, FoxO1, MAFbx, and MuRF1 expression (Table 2). KLF-15, a key transcription factor, not only inhibits skeletal muscle protein synthesis but also plays a crucial role in enhancing skeletal muscle proteolysis by promoting MuRF1 and MAFbx transcription, and stimulates FoxO1 (Kuo et al., 2013). FoxO1, in particular, is the main regulatory factor among all FoxO subtypes, expressed significantly in skeletal muscle. This study demonstrates that the activation of the FoxO signaling pathway is significantly enhanced, as evidenced by the data presented in Figure 6. In the process of KLF-15 mediated skeletal muscle degradation, FoxO plays a pivotal role as a transcription factor and serves as a central component of multiple degradation pathways. Participating in moderate or vigorous physical activity leads to an intensified generation of reactive oxygen species (ROS) (Vezzoli et al., 2016). The FoxO signaling pathway may be activated as a result of increased generation of ROS (Dodd et al., 2010). The importance of this heightened FoxO signaling lies in its role in increasing FoxO nuclear translocation, therefore leading to an intensified transcription of FoxO target genes, such as MAFbx and MuRF-1, which are directly involved in muscle catabolism (Bodine et al., 2001; Gomes et al., 2001; Crow, 2002), and MAFbx and MuRF-1 play a crucial role in specifically targeting proteins for degradation by the proteasome (Bodine et al., 2001; Gomes et al., 2001; Puthanveetil et al., 2013). In addition, we found that the HSP70, UCP3 and SIRT6 genes were significantly upregulated in skeletal muscle after exercise in this study. The increase in expression of the differential gene HSP70 can induce a protective anti-apoptotic response, alleviate membrane protein denaturation caused by injury stimuli, and protect the protein synthesis pathway in cells from damage (Mosser et al., 2000; Mokranjac et al., 2006). This high expression of UCP3 can reduce the release of ROS and prevent damage to the body (Nabben et al., 2008; Toime and Brand, 2010), while SIRT6 overexpression inhibits the protein degradation pathway related to muscle atrophy (Samant et al., 2017). However, their protective effects may not be sufficient to counteract the physical damage caused by acute exercise.
The ubiquitin-proteasome system (UPS) is the primary pathway for intracellular protein degradation (Meng et al., 2014). The importance of the ubiquitination level of proteins for quality control is further highlighted by the fact that approximately 80% of proteins in the body are degraded by the ubiquitin-proteasome system (Cromm and Crews, 2017). There is evidence that abnormal protein metabolism can lead to excessive activation of the ubiquitin-proteasome system, resulting in protein degradation (Cunha et al., 2012). MAFbx and MuRF-1, which are both downstream signaling proteins of FoxO, are widely recognized as crucial skeletal muscle-specific F-box-type ubiquitin ligases involved in protein ubiquitination (Bodine and Baehr, 2014). This study demonstrates an enhanced process of ubiquitin mediated proteolysis, as depicted in Figure 6. Previous studies have examined the expression of MAFbx and MuRF1 genes in skeletal muscle during exercise. Mascher et al. (2008) found that after one session of knee joint resistance exercise in sedentary young men, the mRNA content of MuRF1 increased by 4.8 times. Louis et al. (2007) had subjects perform running and resistance exercises and found increased expression of MuRF1 and MAFbx through a muscle biopsy.
UPS could regulate the degradation of proteins involved in iron metabolism, thereby modulating ferroptosis (Zhang et al., 2010; Horie et al., 2017). In this study, Figure 7 vividly depicts the crucial discovery that ferroptosis serves as a pivotal biological process, identified through KEGG pathway enrichment analysis of differential genes and metabolites. Ferroptosis, a novel type of cell death has been reported (Dixon et al., 2012). It is a specific mechanism characterized by lipid peroxidation due to iron overload, which leads to cell death (Conrad et al., 2016). Studies have shown that ferroptosis is associated with various cellular processes, such as iron homeostasis, redox homeostasis, and lipid metabolism (Yang et al., 2014; Ingold et al., 2018). UPS, as a key regulatory system for homeostasis of eukaryotic cells (Glickman and Ciechanover, 2002). Ferroptosis impacted muscle fiber differentiation and promoted muscle protein degradation. Interestingly, some ubiquitin ligases are also involved in ferroptosis occurrence (Meng et al., 2022), while myasthenia-specific MAFbx and MuRF1 expression increases with iron accumulation in the skeletal muscle of aging individuals (Bodine and Baehr, 2014; Huang et al., 2021).
Thus, take all together, this study identified multiple metabolic pathways and signaling pathways that were involved in the regulation of skeletal muscle, such as the TCA cycle, oxidative phosphorylation, protein digestion and absorption, the FoxO signaling pathway, ubiquitin mediated proteolysis, and ferroptosis signaling pathway. Additionally, the study demonstrated that acute exercise led to the upregulation of KLF-15, FoxO1, MAFbx, and MuRF1 expression, which could play a critical role in enhancing skeletal muscle proteolysis. However, the study mainly focused on the metabolic and transcription profiling of mouse muscle, and did not involve other potential biomarkers or mechanisms, such as protein spectrometry or post-transcriptional modifications of gene expression.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://www.ncbi.nlm.nih.gov/search/all/?term=PRJNA1002978.
Ethics statement
The animal study was approved by Academic Integrity and Research Ethics Committee of China University of Geosciences (Wuhan). The study was conducted in accordance with the local legislation and institutional requirements.
Author contributions
XY: Conceptualization, Data curation, Formal Analysis, Investigation, Methodology, Resources, Software, Validation, Visualization, Writing–original draft, Writing–review and editing. RL: Conceptualization, Data curation, Formal Analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Software, Validation, Visualization, Writing–original draft, Writing–review and editing. ZQ: Methodology, Writing–review and editing. XC: Methodology, Writing–review and editing. YW: Methodology, Writing–review and editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This work was supported by grants from the Fundamental Research Funds for the Central Universities in China (Grant No. CUG150607).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
Baskin K. K., Winders B. R., Olson E. N. (2015). Muscle as a "mediator" of systemic metabolism. Cell Metab. 21, 237–248. doi:10.1016/j.cmet.2014.12.021
Bodine S. C., Baehr L. M. (2014). Skeletal muscle atrophy and the E3 ubiquitin ligases MuRF1 and MAFbx/atrogin-1. Am. J. Physiology-Endocrinology Metabolism 307, E469–E484. doi:10.1152/ajpendo.00204.2014
Bodine S. C., Latres E., Baumhueter S., Lai V. K., Nunez L., Clarke B. A., et al. (2001). Identification of ubiquitin ligases required for skeletal muscle atrophy. Science 294, 1704–1708. doi:10.1126/science.1065874
Camera D. M., Smiles W. J., Hawley J. A. (2016). Exercise-induced skeletal muscle signaling pathways and human athletic performance. Free Radic. Biol. Med. 98, 131–143. doi:10.1016/j.freeradbiomed.2016.02.007
Castellanos N., Diez G. G., Antúnez-Almagro C., Bailén M., Bressa C., González Soltero R., et al. (2020). A critical mutualism–competition interplay underlies the loss of microbial diversity in sedentary lifestyle. Front. Microbiol. 10, 3142. doi:10.3389/fmicb.2019.03142
Conrad M., Angeli J. P. F., Vandenabeele P., Stockwell B. R. (2016). Regulated necrosis: disease relevance and therapeutic opportunities. Nat. Rev. Drug Discov. 15, 348–366. doi:10.1038/nrd.2015.6
Cook D. B., Stegner A. J., Nagelkirk P. R., Meyer J. D., Togo F., Natelson B. H. J. M. (2012). Responses to exercise differ for chronic fatigue syndrome patients with fibromyalgia. Med. Sci. sports Exerc. 44, 1186–1193. doi:10.1249/MSS.0b013e3182417b9a
Cromm P. M., Crews C. M. (2017). The proteasome in modern drug discovery: second life of a highly valuable drug target. ACS Cent. Sci. 3, 830–838. doi:10.1021/acscentsci.7b00252
Crow M. T. (2002). Hypoxia, BNip3 proteins, and the mitochondrial death pathway in cardiomyocytes. Circ. Res. 91, 183–185. doi:10.1161/01.res.0000030195.38795.cf
Cunha T. F., Bacurau A. V., Moreira J. B., Paixao N. A., Campos J. C., Ferreira J. C., et al. (2012). Exercise training prevents oxidative stress and ubiquitin-proteasome system overactivity and reverse skeletal muscle atrophy in heart failure. PLoS One 7, e41701. doi:10.1371/journal.pone.0041701
Defronzo R. A., Tripathy D. (2009). Skeletal muscle insulin resistance is the primary defect in type 2 diabetes. Diabetes Care 32 (2), S157–S163. doi:10.2337/dc09-S302
Dettmer K., Aronov P. A., Hammock B. D. (2007). Mass spectrometry-based metabolomics. Mass Spectrom. Rev. 26, 51–78. doi:10.1002/mas.20108
Dixon S. J., Lemberg K. M., Lamprecht M. R., Skouta R., Zaitsev E. M., Gleason C. E., et al. (2012). Ferroptosis: an iron-dependent form of nonapoptotic cell death. Cell 149, 1060–1072. doi:10.1016/j.cell.2012.03.042
Dodd S. L., Gagnon B. J., Senf S. M., Hain B. A., Judge A. R. (2010). Ros-mediated activation of NF-kappaB and Foxo during muscle disuse. Muscle Nerve 41, 110–113. doi:10.1002/mus.21526
Dohm G. L., Kasperek G. J., Tapscott E. B., Barakat H. A. (1985). Protein metabolism during endurance exercise. Fed. Proc. 44, 348–352.
Egan B., Zierath J. R. (2013). Exercise metabolism and the molecular regulation of skeletal muscle adaptation. Cell Metab. 17, 162–184. doi:10.1016/j.cmet.2012.12.012
Fan W., Evans R. M. (2017). Exercise mimetics: impact on health and performance. Cell Metab. 25, 242–247. doi:10.1016/j.cmet.2016.10.022
Farr J. N., Khosla S. (2019). Cellular senescence in bone. bone 121, 121–133. doi:10.1016/j.bone.2019.01.015
Glickman M. H., Ciechanover A. J. P. R. (2002). The ubiquitin-proteasome proteolytic pathway: destruction for the sake of construction. Physiol. Rev. 82, 373–428. doi:10.1152/physrev.00027.2001
Gomes M. D., Lecker S. H., Jagoe R. T., Navon A., Goldberg A. L. (2001). Atrogin-1, a muscle-specific F-box protein highly expressed during muscle atrophy. Proc. Natl. Acad. Sci. U. S. A. 98, 14440–14445. doi:10.1073/pnas.251541198
Harris R. A., Joshi M., Jeoung N. H., Obayashi M. (2005). Overview of the molecular and biochemical basis of branched-chain amino acid catabolism. J. Nutr. 135, 1527S–1530S. doi:10.1093/jn/135.6.1527S
Hochachka P. W., Dressendorfer R. H. (1976). Succinate accumulation in man during exercise. Eur. J. Appl. Physiol. Occup. Physiol. 35, 235–242. doi:10.1007/BF00423282
Horie T., Kawamata T., Matsunami M., Ohsumi Y. (2017). Recycling of iron via autophagy is critical for the transition from glycolytic to respiratory growth. J. Biol. Chem. 292, 8533–8543. doi:10.1074/jbc.M116.762963
Huang Y., Wu B., Shen D., Chen J., Yu Z., Chen C. (2021). Family function fully mediates the relationship between social support and perinatal depression in rural Southwest China. Int. J. Biol. Sci. 17, 151. doi:10.1186/s12888-021-03155-9
Ingold I., Berndt C., Schmitt S., Doll S., Poschmann G., Buday K., et al. (2018). Selenium utilization by GPX4 is required to prevent hydroperoxide-induced ferroptosis. Cell 172, 409–422. doi:10.1016/j.cell.2017.11.048
Kemi O. J., Loennechen J. P., Wisloff U., Ellingsen O. (2002). Intensity-controlled treadmill running in mice: cardiac and skeletal muscle hypertrophy. J. Appl. Physiol. 93, 1301–1309. doi:10.1152/japplphysiol.00231.2002
Kim S., Hong J., Jeon R., Kim H.-S. (2016). Adzuki bean ameliorates hepatic lipogenesis and proinflammatory mediator expression in mice fed a high-cholesterol and high-fat diet to induce nonalcoholic fatty liver disease. Nutr. Res. 36, 90–100. doi:10.1016/j.nutres.2015.11.002
Kruger K., Frost S., Most E., Volker K., Pallauf J., Mooren F. C. (2009). Exercise affects tissue lymphocyte apoptosis via redox-sensitive and Fas-dependent signaling pathways. Am. J. Physiol. Regul. Integr. Comp. Physiol. 296, R1518–R1527. doi:10.1152/ajpregu.90994.2008
Kuo T., Harris C. A., Wang J.-C. (2013). Metabolic functions of glucocorticoid receptor in skeletal muscle. Mol. Cell. Endocrinol. 380, 79–88. doi:10.1016/j.mce.2013.03.003
Langmead B., Salzberg S. L. (2012). Fast gapped-read alignment with Bowtie 2. Nat. Methods 9, 357–359. doi:10.1038/nmeth.1923
Li B., Dewey C. N. (2011). RSEM: accurate transcript quantification from RNA-seq data with or without a reference genome. BMC Bioinforma. 12, 323. doi:10.1186/1471-2105-12-323
Louis E., Raue U., Yang Y., Jemiolo B., Trappe S. (2007). Time course of proteolytic, cytokine, and myostatin gene expression after acute exercise in human skeletal muscle. J. Appl. physiology 103, 1744–1751. doi:10.1152/japplphysiol.00679.2007
Magalhaes J., Fraga M., Lumini-Oliveira J., Goncalves I., Costa M., Ferreira R., et al. (2013). Eccentric exercise transiently affects mice skeletal muscle mitochondrial function. Appl. Physiol. Nutr. Metab. 38, 401–409. doi:10.1139/apnm-2012-0226
Mascher H., Tannerstedt J., Brink-Elfegoun T., Ekblom B., Gustafsson T., Blomstrand E. (2008). Repeated resistance exercise training induces different changes in mRNA expression of MAFbx and MuRF-1 in human skeletal muscle. Am. J. Physiology Endocrinol. Metabolism 294, E43–E51. doi:10.1152/ajpendo.00504.2007
Meng J., Zou X., Wu R., Zhong R., Zhu D., Zhang Y. (2014). Accelerated regeneration of the skeletal muscle in RNF13-knockout mice is mediated by macrophage-secreted IL-4/IL-6. Protein and Cell 5, 235–247. doi:10.1007/s13238-014-0025-4
Meng Y., Sun H., Li Y., Zhao S., Su J., Zeng F., et al. (2022). Targeting ferroptosis by ubiquitin system enzymes: a potential therapeutic strategy in cancer. Int. J. Biol. Sci. 18, 5475–5488. doi:10.7150/ijbs.73790
Mokranjac D., Bourenkov G., Hell K., Neupert W., Groll M. (2006). Structure and function of Tim14 and Tim16, the J and J-like components of the mitochondrial protein import motor. EMBO J. 25, 4675–4685. doi:10.1038/sj.emboj.7601334
Mosser D. D., Caron A. W., Bourget L., Meriin A. B., Sherman M. Y., Morimoto R. I., et al. (2000). The chaperone function of hsp70 is required for protection against stress-induced apoptosis. Mol. Cell. Biol. 20, 7146–7159. doi:10.1128/mcb.20.19.7146-7159.2000
Nabben M., Hoeks J., Briedé J. J., Glatz J. F., Moonen-Kornips E., Hesselink M. K., et al. (2008). The effect of UCP3 overexpression on mitochondrial ROS production in skeletal muscle of young versus aged mice. FEBS Lett. 582, 4147–4152. doi:10.1016/j.febslet.2008.11.016
Naseer M. I., Abdulkareem A. A., Jan M. M., Chaudhary A. G., Alharazy S., Alqahtani M. H. (2020). Next generation sequencing reveals novel homozygous frameshift in PUS7 and splice acceptor variants in AASS gene leading to intellectual disability, developmental delay, dysmorphic feature and microcephaly. Saudi J. Biol. Sci. 27, 3125–3131. doi:10.1016/j.sjbs.2020.09.033
O'neill C. A., Stebbins C. L., Bonigut S., Halliwell B., Longhurst J. C. (1996). Production of hydroxyl radicals in contracting skeletal muscle of cats. J. Appl. Physiol. 81, 1197–1206. doi:10.1152/jappl.1996.81.3.1197
Puthanveetil P., Wan A., Rodrigues B. (2013). FoxO1 is crucial for sustaining cardiomyocyte metabolism and cell survival. Cardiovasc Res. 97, 393–403. doi:10.1093/cvr/cvs426
Qiu S. T., Cui D., Lu J., Chen C. Z. (2014). Effect of AMPK/Sirt1 pathway on the muscle mass control induced by exercise. J. Tianjin Univ. Sport 27, 255–258.
Reid M. B., Haack K. E., Franchek K. M., Valberg P. A., Kobzik L., West M. S. (1992). Reactive oxygen in skeletal muscle. I. Intracellular oxidant kinetics and fatigue in vitro. J. Appl. Physiol. 73, 1797–1804. doi:10.1152/jappl.1992.73.5.1797
Rennie M. J., Edwards R. H., Davies C. T., Krywawych S., Halliday D., Waterlow J. C., et al. (1980). Protein and amino acid turnover during and after exercise. Biochem. Soc. Trans. 8, 499–501. doi:10.1042/bst0080499
Samant S. A., Kanwal A., Pillai V. B., Bao R., Gupta M. P. (2017). The histone deacetylase SIRT6 blocks myostatin expression and development of muscle atrophy. Sci. Rep. 7, 11877. doi:10.1038/s41598-017-10838-5
Schoepe M., Schrepper A., Schwarzer M., Osterholt M., Doenst T. (2012). Exercise can induce temporary mitochondrial and contractile dysfunction linked to impaired respiratory chain complex activity. Metabolism-clinical Exp. 61, 117–126. doi:10.1016/j.metabol.2011.05.023
Shen F., Zhao Y., Ding W., Liu K., Ren X., Zhang Q., et al. (2021). Autonomous climbing: an effective exercise mode with beneficial outcomes of aerobic exercise and resistance training. Life Sci. 265, 118786. doi:10.1016/j.lfs.2020.118786
Song Y., Shi X., Gao Z., Li R., Tian J., Cao X., et al. (2023). Acupoint catgut embedding improves lipid metabolism in exercise-induced fatigue rats via the PPAR signaling pathway. Anim. (Basel) 13, 558. doi:10.3390/ani13040558
Toime L. J., Brand M. D. (2010). Uncoupling protein-3 lowers reactive oxygen species production in isolated mitochondria. Free Radic. Biol. Med. 49, 606–611. doi:10.1016/j.freeradbiomed.2010.05.010
Tondera D., Grandemange S., Jourdain A., Karbowski M., Mattenberger Y., Herzig S., et al. (2009). SLP-2 is required for stress-induced mitochondrial hyperfusion. EMBO J. 28, 1589–1600. doi:10.1038/emboj.2009.89
Vezzoli A., Dellanoce C., Mrakic-Sposta S., Montorsi M., Moretti S., Tonini A., et al. (2016). Oxidative stress assessment in response to ultraendurance exercise: thiols redox status and ROS production according to duration of a competitive race. Oxid. Med. Cell Longev. 2016, 6439037. doi:10.1155/2016/6439037
Wang F., Wang X., Liu Y., Zhang Z. (2021). Effects of exercise-induced ROS on the pathophysiological functions of skeletal muscle. Oxidative Med. Cell. Longev. 2021, 3846122–3846125. doi:10.1155/2021/3846122
Wang Z., Gerstein M., Snyder M. (2009). RNA-seq: a revolutionary tool for transcriptomics. Nat. Rev. Genet. 10, 57–63. doi:10.1038/nrg2484
Yang W. S., Sriramaratnam R., Welsch M. E., Shimada K., Skouta R., Viswanathan V. S., et al. (2014). Regulation of ferroptotic cancer cell death by GPX4. Cell 156, 317–331. doi:10.1016/j.cell.2013.12.010
Zhang Y., Mikhael M., Xu D., Li Y., Soe-Lin S., Ning B., et al. (2010). Lysosomal proteolysis is the primary degradation pathway for cytosolic ferritin and cytosolic ferritin degradation is necessary for iron exit. Antioxidants redox Signal. 13, 999–1009. doi:10.1089/ars.2010.3129
Keywords: acute exercise, metabolome, transcriptome, skeletal muscle, quadriceps, mice
Citation: Ye X, Liu R, Qiao Z, Chai X and Wang Y (2023) Integrative profiling of metabolome and transcriptome of skeletal muscle after acute exercise intervention in mice. Front. Physiol. 14:1273342. doi: 10.3389/fphys.2023.1273342
Received: 06 August 2023; Accepted: 25 September 2023;
Published: 06 October 2023.
Edited by:
Simone Luti, University of Florence, ItalyReviewed by:
Elen H. Miyabara, University of São Paulo, BrazilMaurice H. T. Ling, Temasek Polytechnic, Singapore
Copyright © 2023 Ye, Liu, Qiao, Chai and Wang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Renyi Liu, cmVueWkubGl1QGN1Zy5lZHUuY24=