Skip to main content

ORIGINAL RESEARCH article

Front. Genet., 12 July 2022
Sec. Behavioral and Psychiatric Genetics

Skeletal Muscle Transcriptomic Comparison Between Men and Women in Response to Acute Sprint Exercise

Mingkun Nie&#x;Mingkun Nie1Qingling Liu&#x;Qingling Liu2Cheng Yan,,
Cheng Yan2,3,4*
  • 1School of Physical Education, Xinxiang University, Xinxiang, China
  • 2School of Pharmacy, Xinxiang University, Xinxiang, China
  • 3Key Laboratory of Nano-carbon Modified Film Technology of Henan Province, Xinxiang University, Xinxiang, China
  • 4Diagnostic Laboratory of Animal Diseases, Xinxiang University, Xinxiang, China

Background: Acute sprint exercise is a time-efficient physical activity that improves cardiorespiratory fitness in younger and middle-aged adults. Growing evidence has demonstrated that acute sprint exercise provides equal to or superior health benefits compared with moderate-intensity continuous training, which will dramatically increase aerobic capacity, insulin sensitivity, and muscle capillarization. Although the beneficial effects of acute sprint exercise are well documented, the mechanisms behind how acute sprint exercise prevents disease and benefits health are less understood.

Method: We obtained differentially expressed genes in muscle (vastus lateralis) from men and women before and after an acute sprint exercise. Then, we identified hub genes from the protein–protein interaction (PPI) network of differentially expressed genes (DEGs) and key transcription factors in men and women related to acute sprint exercise. Finally, Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) enrichment analyses are performed on DEGs and sex-biased genes, respectively.

Results: First, we identified 127 sexually dimorphic genes in men (90 upregulated and 37 downregulated) and 75 genes in women (90 upregulated and 37 downregulated) in response to acute sprint exercise. Second, CEBPB, SMAD3, and CDKN1A are identified as the top three hub genes related to men-biased genes. Accordingly, the top three hub genes related to women-biased genes are JUN, ACTB, and SMAD7. In addition, CLOCK, ZNF217, and KDM2B are the top three enriched transcriptional factors in men-biased genes, while XLR, SOX2, JUND, and KLF4 are transcription factors enriched most in women-biased genes. Furthermore, based on GO and KEGG enrichment analyses, we identified potential key pathways in regulating the exercise-related response in men and women, respectively.

Conclusion: In this study, we found the difference in gene expression and enrichment pathways in muscle in men and women in response to acute sprint exercise. These results will shed new light on the mechanism underlying sex-based differences in skeletal muscle remodeling and metabolism related to acute sprint exercise, which may illustrate the mechanisms behind how acute sprint exercise prevents disease and benefits health.

Introduction

Regular physical activity is well accepted as a key intervention for reducing the risk of chronic disease, including mobility disability, heart disease, type 2 diabetes, and several cancers (Hawley, 2004; Warburton et al., 2006). Acute sprint exercise is a time-efficient physical activity that improves cardiorespiratory fitness in younger and middle-aged adults (Willoughby et al., 2016), promotes favorable enhancement in insulin sensitivity (Richards et al., 2010), and has beneficial effects on the enhancement of cognitive abilities with increased neuroprotective factors (Kujach et al., 2019). The previous study has shown that acute sprint exercise provides equal to or superior health benefits compared with moderate-intensity continuous training, which will dramatically increase aerobic capacity, insulin sensitivity, and muscle capillarization (Cocks et al., 2016). Different from other exercise manners, acute sprint exercise is characterized by rapid glycogen degradation and large net ATP breakdown in muscle and is accompanied by the increased secretion of hormones such as insulin, catecholamines, and growth factor (Esbjörnsson-Liljedahl et al., 1985; Esbjörnsson et al., 2009). Gene signatures related to a turnover of skeletal muscle mass are upregulated after sprint exercise (Rundqvist et al., 2019). Although the beneficial effects of acute sprint exercise are well documented, the mechanisms behind how acute sprint exercise prevents disease and benefits health are less understood.

Skeletal muscles have the fourth largest number of sex differentially expressed genes across 29 human tissues, with more than 3,000 differentially expressed genes between male and female skeletal muscles at rest (Welle et al., 2008; Lopes-Ramos et al., 2020). Thus, differentially expressed genes may be the most important cause of sex differences in muscle morphology, function, and plasticity (Pointer et al., 2013). According to a previous study, men usually have significantly larger cross-sectional areas than women, especially in type II fibers (Miller et al., 1993). Conversely, women may be more resistant to muscle fatigue due to a higher proportion of type I fibers (Salomoni et al., 2008; Avin et al., 2010). Furthermore, men and women have sex-based differences in response to sprint interval training. It is reported that after six weeks of sprint interval training, insulin sensitivity was increased only in healthy men but not women (Metcalfe et al., 2012). Sprint interval training reduced 24-h blood glucose concentration in men but not women among adults with overweight (Gillen et al., 2014). Glucose transporter 4 (GLUT4), a key regulator of glucose metabolism in muscle, is higher increased in men compared with women in response to sprint interval training. In summary, previous studies indicate that men showed stronger adaptative responses to sprint interval training than women. However, the mechanism underlying sex-based differences in skeletal muscle remodeling and metabolism in response to acute sprint exercise is still unclear. To the best of our knowledge, the difference of global gene expression in skeletal muscle between men and women in response to once acute sprint exercise (three bouts of 30-s sprint exercise with 20-min rest) has not been previously demonstrated.

In this study, we focused on the response to acute sprint exercise, not long-term adaptation exercise training, to increase the understanding of the transcriptomic difference of skeletal muscle in response to acute sprint exercise between men and women. First, we obtained differentially expressed genes in muscle from men and women before and after an acute sprint exercise. Then, we identified hub genes from the protein–protein interaction (PPI) network of DEGs in men and women related to acute sprint exercise. Finally, KEGG and GO enrichment analyses were performed on DEGs and sex-biased genes, respectively. Our study shed light on the understanding of how the health benefits of acute sprint exercise are regulated in men and women. In addition, our results may also indicate the mechanism underlying sex-based differences in skeletal muscle remodeling and metabolism related to acute sprint exercise.

Materials and Methods

Acquisition of Expression Profile

The gene expression profiles related to acute sprint exercise are downloaded from GSE126296 (PMID: 31647849), which contains the expression data of quadriceps femoris muscle (vastus lateralis) from seven males and seven females before and after acute sprint exercise (Rundqvist et al., 2019). As described in the original study, participants have to be in good general health, participate in leisure-time sports but not at an elite level, and be between 20 and 30 years of age. The participants with chronic illness, acute infection, pregnancy, severe asthma, or use of nicotine-containing products are excluded.

Identification of Differentially Expressed Genes

Differentially expressed genes are identified with “limma” package in R software. The screening threshold used in this study is set at FDR <0.05 and |log2 (fold change)| > 0.3. Heat maps and volcano maps were visualized by employing the “ggplot2” and “pheatmap” packages in R software.

Construction of Protein–Protein Interaction Network and Identification of Hub Genes

A PPI network is constructed based on the STRING database. To identify hub genes, we use the radiality algorithm in cytoHubba plug-in in Cytoscape.

Kyoto Encyclopedia of Genes and Genomes, Gene Ontology, and Transcription Factor Enrichment Analyses

KEGG, GO, and transcription factor analyses are performed using the Enrichr tool (https://maayanlab.cloud/Enrichr/).

Results

Analysis of Differentially Expressed Genes in Men and Women in Response to Acute Sprint Exercise

A flowchart of analytical approaches used in this study is provided in Figure 1. In order to identify the differentially expressed genes in skeletal muscle between men and women in response to acute sprint exercise, we performed bioinformatics analysis on leg skeletal muscle collected from human subjects before and after three all-out cycle sprints. We set the threshold at FDR <0.05 and | log2 (fold change) | > 0.3. A total of 234 genes (192 upregulated and 42 downregulated) in men were identified as differentially expressed genes in response to acute sprint exercise (Supplementary Table S1). The heat map and volcano plot for differentially expressed genes are shown in Figures 2A and B. Accordingly, 182 differentially expressed genes are identified in women, containing 157 upregulated and 25 downregulated genes (Supplementary Table S2). The heat map and volcano plot for differentially expressed genes are shown in Figures 2C and D. Furthermore, we identified 127 sexually dimorphic genes in men (90 upregulated and 37 downregulated) and 75 genes in women (55 upregulated and 20 downregulated) in response to acute sprint exercise (Figures 2E,F) (Supplementary Tables S3, S4). These results demonstrate that numerous sexually dimorphic genes are differentially expressed in men and women, and the amount of men-biased genes is much greater than that of women-biased genes in response to acute sprint exercise.

FIGURE 1
www.frontiersin.org

FIGURE 1. A flowchart of analytical approaches used in this study. GEO, Gene Expression Omnibus; GO, Gene Ontology; KEGG, Kyoto Encyclopedia of Genes and Genomes.

FIGURE 2
www.frontiersin.org

FIGURE 2. Analysis of differentially expressed genes in response to acute sprint exercise in skeletal muscle from men and women. (A) Heat map of 234 differentially expressed genes in muscle from men. (B) Volcano plot of 234 differentially expressed genes in muscle from men. (C) Heat map of 182 differentially expressed genes in muscle from women. (D) Volcano plot of 182 differentially expressed genes in muscle from women. (E) Venn diagram showing the overlap of differentially expressed genes between men and women. (F) Venn diagrams showing the overlap of the differentially expressed genes that were upregulated and downregulated in men and women. Red dots represent upregulated genes, and blue dots indicate downregulated genes. FC, fold change.

Identifying Hub Genes and Transcription Factors Associated With Sex-Specific Expression Signatures in Response to Acute Sprint Exercise

In order to explore the network of sexually dimorphic genes in response to acute sprint exercise, we constructed the PPI networks based on male-biased genes and female-biased genes, respectively (Figures 3A,B). Hub genes from the PPI network were screened out with cytoHubba application in Cytoscape using the radiality method (Figures 3C,D). The top three hub genes related to men-biased genes are CEBPB, SMAD3, and CDKN1A. Accordingly, the top three hub genes related to women-biased genes are JUN, ACTB, and SMAD7. These genes may encode core proteins with important biological regulatory functions resulting in sex-based differences in skeletal muscle remodeling and metabolism in response to acute sprint exercise. Subsequentially, we performed transcription factor enrichment analysis to provide insights into the underlying change of transcription activity in sex-specific expression signatures. CLOCK, ZNF217, and KDM2B are the top three enriched transcriptional factors in men-biased genes (Figure 3E; Supplementary Table S5), while XLR, SOX2, JUND, and KLF4 are transcription factors enriched most in women-biased genes (Figure 3F; Supplementary Table S6). These results suggest that sexually dimorphic genes may be regulated by different transcription factors in response to acute sprint exercise.

FIGURE 3
www.frontiersin.org

FIGURE 3. Construction of PPI network and identification of differentially expressed genes in men and women. (A) PPI network of differentially expressed genes in men. (B) PPI network of differentially expressed genes in women. (C) Hub gene network of differentially expressed genes in men. (D) Hub gene network of differentially expressed genes in women. (E) Transcription factor enrichment of men-biased genes. (F) Transcription factor enrichment of women-biased genes.

Gene Ontology and Kyoto Encyclopedia of Genes and Genomes Enrichment Analyses of Differentially Expressed Genes in Men and Women in Response to Acute Sprint Exercise

To further explore the functions associated with DEGs, we conducted a series of enrichment analyses. GO enrichment analysis showed that 389 terms were significantly enriched (p-value < 0.05) in men, including eight cellular component terms, 60 molecular function terms, and 321 biological process terms (Figures 4A,C,E). Notable pathways among BP were positive regulation of transcription by RNA polymerase II, positive regulation of transcription, DNA-templated, and cellular response to starvation. In the molecular function category, the top three terms were sequence-specific DNA binding, RNA polymerase II cis-regulatory region sequence-specific DNA binding, and cis-regulatory region sequence-specific DNA binding. Moreover, nucleus, cell projection membrane, and external side of apical plasma membrane were the top three significant terms in the cellular component category. KEGG pathway analysis identified 57 significantly enriched pathways such as MAPK signaling pathway, mTOR signaling pathway, and AGE-RAGE signaling pathway in diabetic complications (Figure 4G), while GO enrichment analysis in women showed that 176 terms were significantly enriched (p-value < 0.05), including seven cellular component terms, 77 molecular function terms, and 392 biological process terms (Figures 4B,D,F). The BPs of the DEGs from the female subjects were primarily enriched in positive regulation of transcription, DNA-templated, positive regulation of transcription by RNA polymerase II, and regulation of transcription by RNA polymerase II. The CC results revealed that DEGs in the female cohorts were mainly enriched in nucleus, intracellular membrane-bounded organelle, and external side of apical plasma membrane. The MF results indicated that DEGs identified in the female cohorts were primarily enriched in sequence-specific DNA binding, RNA polymerase II cis-regulatory region sequence-specific DNA binding, and cis-regulatory region sequence-specific DNA binding. KEGG pathway analysis identified 69 significantly enriched pathways in female, including MAPK signaling pathway, fluid shear stress, and mTOR signaling pathway (Figure 4H). Furthermore, we use the Venn diagram to indicate the overlap of items of GO and KEGG between men and women (Figures 4I–L, Supplementary Tables S5, S6). These results indicate that changes of signaling pathways are different in skeletal muscle between men and women in response to acute sprint exercise.

FIGURE 4
www.frontiersin.org

FIGURE 4. GO and KEGG enrichment analyses of DEGs in men and women. (A) BP enrichment analyses of differentially expressed genes in men. (B) BP enrichment analyses of differentially expressed genes in women. (C) CC enrichment analyses of differentially expressed genes in men. (D) CC enrichment analyses of differentially expressed genes in women. (E) MF enrichment analyses of differentially expressed genes in men. (F) MF enrichment analyses of differentially expressed genes in women. (G) KEGG enrichment analyses of differentially expressed genes in men. (H) KEGG enrichment analyses of differentially expressed genes in women. (I) Venn diagram showing the overlap of BP pathways between men and women. (J) Venn diagram showing the overlap of CC pathways between men and women. (K) Venn diagram showing the overlap of MF pathways between men and women. (L) Venn diagram showing the overlap of KEGG pathways between men and women. GO, Gene Ontology; BP, biological process; CC, cellular component; MF, molecular function; KEGG, Kyoto Encyclopedia of Genes and Genomes; TF, transcriptional factor.

Gene Ontology and Kyoto Encyclopedia of Genes and Genomes Enrichment Analyses of Male-Biased and Female-Biased Genes in Response to Acute Sprint Exercise

To further explore the functions associated with sex-biased genes, we conducted a series of enrichment analyses. GO enrichment analysis showed that 235 terms were significantly enriched (p-value < 0.05) in men, including eight cellular component terms, 40 molecular function terms, and 187 biological process terms (Figures 5A,C,E). Notable pathways among BP were cellular response to amino acid starvation, response to amino acid starvation, and cellular response to starvation. In the molecular function category, the top three terms were 1-phosphatidylinositol-3-kinase regulator activity, phosphatidylinositol 3-kinase regulator activity, and kinase inhibitor activity. Moreover, cell projection membrane, brush border membrane, and filtration diaphragm were the top three significant terms in the cellular component category. KEGG pathway analysis identified 16 significantly enriched pathways such as JAK-STAT signaling pathway, renal cell carcinoma, and prolactin signaling pathway (Figure 5G). GO enrichment analysis in women showed that 300 terms were significantly enriched (p-value < 0.05), including three cellular component terms, 51 molecular function terms, and 246 biological process terms (Figures 5B,D,F). The BPs of the DEGs from the female subjects were primarily enriched cellular response to cadmium ion, cellular zinc ion homeostasis, and response to cadmium ion. The CC results revealed that DEGs in the female cohorts were mainly enriched in extrinsic component of endoplasmic reticulum membrane, nucleus, and apical dendrite. The MF results indicated that DEGs identified in the female cohorts were primarily enriched in thiol-dependent deubiquitinase, deubiquitinase activity, and omega peptidase activity. KEGG pathway analysis identified 36 significantly enriched pathways in females, including mineral absorption, apoptosis, and focal adhesion (Figure 5H). Furthermore, we used the Venn diagram to indicate the overlap of items of GO and KEGG between men- and women-biased genes (Figures 5I–L; Supplementary Tables S7, S8). These results indicate that the signaling pathways enriched by men-biased and women-biased genes are different.

FIGURE 5
www.frontiersin.org

FIGURE 5. GO and KEGG enrichment analyses of sexually dimorphic genes in men and women. (A) BP analyses of sexually dimorphic genes in men. (B) GO pathway analyses of sexually dimorphic genes in women. (C) CC pathway analyses of sexually dimorphic genes in men. (D) CC pathway analyses of sexually dimorphic genes in women. (E) MF enrichment analyses of sexually dimorphic genes in men. (F) MF enrichment analyses of sexually dimorphic genes in women. (G) KEGG enrichment analyses of sexually dimorphic genes in men. (H) KEGG enrichment analyses of sexually dimorphic genes in women. (I) Venn diagram showing the overlap of BP pathways between men-biased and women-biased genes. (J) Venn diagram showing the overlap of CC pathways between men-biased and women-biased genes. (K) Venn diagram showing the overlap of MF pathways between men-biased and women-biased genes. (L) Venn diagram showing the overlap of KEGG pathways between men-biased and women-biased genes. GO, Gene Ontology; BP, biological process; CC, cellular component; MF, molecular function; KEGG, Kyoto Encyclopedia of Genes and Genomes; TF, transcriptional factor.

Discussion

Acute sprint exercise induces numerous adaptations in the body resulting in beneficial effects on health, which may induce fat loss, increase cognitive functions, improve cardiorespiratory fitness, and alleviate metabolic syndrome (Gibala et al., 2006; Vollaard and Metcalfe, 2017; Kujach et al., 2019). Although the effects of sprint exercise-induced adaptations to skeletal muscle have been well-established in various diseases, the difference of global gene expression in skeletal muscle in men and women in response to acute bouts of acute sprint exercise has not been previously demonstrated. In this study, we conducted an analysis of gene expression profiles obtained from muscle biopsies collected from human volunteers before and after acute bouts of acute sprint exercise to discover exercise-related sexually dimorphic genes, which might be part of the mechanism underlying the sex-based differences in skeletal muscle remodeling and metabolism in response to acute sprint exercise.

RNA-seq results show that about 3,000 genes of human in skeletal muscle at rest are identified as differentially expressed genes (Lindholm et al., 2014). However, the sex-specific epigenomic signatures of immediate early genes in response to acute sprint exercise are still unknown. We performed differentially expressed gene analysis to explore the difference in gene expression between men and women in response to acute sprint exercise. A total of 234 genes (192 upregulated and 42 downregulated) in men were identified as differentially expressed genes using bioinformatics analysis, while only 182 differentially expressed genes were identified in women, containing 157 upregulated and 25 downregulated genes. The number of differentially expressed genes in men is dramatically greater than that in women. This result may account for the fact that men have a stronger adaptative response to sprint interval training than women. Cytokine-inducible SH2-containing protein (CISH), known as cytokine-induced STAT inhibitor, is the top upregulated gene in males. A recent study showed that CISH may play an important role in promoting free fatty acid mobilization (Vestergaard et al., 2014; Vendelbo et al., 2015). Our data also suggest that the expression of CISH may be also regulated by sex hormones. After acute sprint exercise, OTUD1 is the top fold-change gene in females. Although previous independent studies showed that OTUD1 was significantly induced by acute sprint exercise, its function related to acute sprint exercise was still unknown (Bolotta et al., 2020). Furthermore, we identified ZNF217 and CLOCK as the top significant enrichment transcription factors in men-biased genes. CLOCK is necessary for the maintenance of the skeletal muscle phenotype and function through regulating MyoD. Skeletal muscles from mice with mutation of Clock and Bmal1 exhibit ∼30% reductions in normalized maximal force (Andrews et al., 2010). Our result is consistent with previous research that CLOCK is a key regulator of muscle fibers, which may be one of the reasons inducing the difference in muscle type between men and women. However, LXR and SOX2 are the top significant enrichment transcription factors in women. LXRβ is the dominant LXR subtype in skeletal muscle regulating lipogenesis and cholesterol efflux (Hessvik et al., 2010). This result indicates that transcription factors induced by acute sprint are different between men and women. Since the biological function of ZNF217 and SOX2 in skeletal muscle is still unclear, we can only speculate that interference of these two factors may play a role in muscle groups of different fiber type composition or metabolism in men or women, respectively.

Our study also showed that men and women shared several GO and KEGG pathways in response to acute sprint exercise. The significantly enriched GO terms for DEGs in men and women can be categorized into several biological functional groups including positive regulation of transcription, cellular response to starvation, and positive regulation of nucleic acid-templated transcription. These results indicate that acute sprint exercise promoted cell proliferation in muscle both in men and women, which is consistent with the previous study (Pearson, 1990; Bazgir et al., 2017; Saito et al., 2020). Moreover, DEGs in men and women significantly enrich a list of KEGG pathways relevant to the MAPK signaling pathway, mTOR signaling pathway, and AGE-RAGE signaling pathway. The AKT/mTOR signaling pathway played an important role in regulating muscle hypertrophy (Bodine et al., 2001). Similar to AKT/mTOR signaling, MAPK also promoted improvements in fuel homeostasis and could prevent skeletal muscle atrophy (Kramer and Goodyear, 1985). What is more, emerging evidence shows that RAGE signaling pathway is critical for skeletal muscle physiology controlling both the activity of muscle precursors during skeletal muscle development and the correct time of muscle regeneration after acute injury (Riuzzi et al., 2018). These data may implicate a mechanism behind muscle growth both in men and women.

Exercise training is associated with profound changes in skeletal muscle, including increased abundance of glucose up-take, activation of AMPK, and mitochondrial biogenesis (Egan and Zierath, 2013). Our results also provide novel insights into the difference in the pathway between men and women in response to acute sprint. Compared with women-biased genes, the enrichment pathways of men-biased genes are mainly about response to starvation or amino acid starvation (GO: cellular response to amino acid starvation). Men always consume more energy than women. Higher enrichment of these pathways may support an increased rate of energy consumption in men. Therefore, men are more urgent to face the “insufficient energy” situation. However, there is a tendency for women to restore the calcium and zinc homeostasis (GO: cellular response to cadmium ion). These results indicate acute sprint may induce a more severe disturbance of ion homeostasis in cellular processes. Moreover, our study also identified novel sex-specific areas, such as cellular response to cadmium ion, that exist in the fine regulatory system of the muscle response to acute sprint exercise at the transcriptional level which merit further exploration. It is well known that gender differences exist in skeletal muscle mass and distribution (Janssen et al., 1985). Females often have less total and lean body mass, a higher body fat percentage, and a smaller muscle fiber cross-sectional area (Staron et al., 1985). KEGG analysis showed that the JAK-STAT signaling pathway is dramatically enriched in the men-biased genes, not in women. The JAK/STAT pathway participates in lipid metabolism in muscle. Activation of the JAK/STAT pathway in skeletal muscle promotes muscle hypertrophy by increasing the proliferation of satellite cells (Moresi et al., 2019). These results indicate that the JAK-STAT pathway may account for gender differences in muscle fiber-type distribution and substrate metabolism in men and women. These findings may also help us to understand how sex differences manifest in gene regulatory networks in muscle.

Sample size significantly affects DEG analysis output in microarray datasets, and more biological samples predict more accurate outcomes (Stretch et al., 2013). The gene expression profile we analyzed is from previously published literature. In the original study, Rundqvist et al. (2019) expressed reservations about the identification of sexually dimorphic changes due to the small sample size. As described in the method, we downloaded the expression profile of muscle from the original dataset including seven males and seven females before and after acute sprint exercise. Then, we identified the differentially expressed genes using the “limma” package. The screening threshold used in this study was FDR<0.05 and |log2 (fold change)| > 0.3. Through this manner, we identified a total of 234 DEGs in men and 182 DEGs in women. We believe that the biological number in our study (7 vs. 7) is sufficient to get the significant DEGs to a certain extent. We also believe that more samples included will increase the statistical power and further improve the conclusion in the future.

This study needs to be improved in the future for the following reasons: first, the protein level of differentially expressed genes remains to be confirmed. Second, further studies are required to understand the role of sex-biased genes in response to acute sprint exercise. Third, even if the health benefits brought by sprint exercise are similar to those of endurance training, it also should be mentioned that acute response to exercise is not always similar to the long-term adaptation exercise training.

Conclusion

In this study, we found sex-biased genes and their related potential key pathways in muscle in response to acute sprint exercise. Thus, these results will shed new light on the mechanism underlying sex-based differences in skeletal muscle remodeling and metabolism related to sprint interval training, which may also illustrate the mechanisms behind how acute sprint exercise prevents disease and benefits health.

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 in the article/Supplementary Material.

Author Contributions

CY designed the research and wrote the manuscript. MN and QL downloaded and analyzed the data. All authors have agreed to the manuscript.

Funding

This work was supported by the Key Scientific and Technological Project of Henan Province (Grant No. 202102310416) and the Key Scientific Research Projects of Henan Colleges and Universities (Grant No. 21A310008).

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.

Acknowledgments

The authors thank participants and staff of Xinxiang University for their contributions.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2022.860815/full#supplementary-material

Abbreviations

AGE, advanced glycation end; BP, biological process; CC, cellular component; CISH, cytokine-inducible SH2-containing protein; DEGs, differentially expressed genes; FDR, false discovery rate; GO, Gene Ontology; KEGG, Kyoto Encyclopedia of Genes and Genomes; MF, molecular function; PPI, protein–protein interaction; RAGE, receptor for advanced glycation end.

References

Andrews, J. L., Zhang, X., McCarthy, J. J., McDearmon, E. L., Hornberger, T. A., Russell, B., et al. (2010). CLOCK and BMAL1 Regulate MyoD and Are Necessary for Maintenance of Skeletal Muscle Phenotype and Function. Proc. Natl. Acad. Sci. U.S.A. 107 (44), 19090–19095. doi:10.1073/pnas.1014523107

PubMed Abstract | CrossRef Full Text | Google Scholar

Avin, K. G., Naughton, M. R., Ford, B. W., Moore, H. E., Monitto-Webber, M. N., Stark, A. M., et al. (2010). Sex Differences in Fatigue Resistance Are Muscle Group Dependent. Med. Sci. Sports Exerc 42 (10), 1943–1950. doi:10.1249/mss.0b013e3181d8f8fa

PubMed Abstract | CrossRef Full Text | Google Scholar

Bazgir, B., Fathi, R., Rezazadeh Valojerdi, M., Mozdziak, P., and Asgari, A. (2017). Satellite Cells Contribution to Exercise Mediated Muscle Hypertrophy and Repair. Cell. J. 18 (4), 473–484. doi:10.22074/cellj.2016.4714

PubMed Abstract | CrossRef Full Text | Google Scholar

Bodine, S. C., Stitt, T. N., Gonzalez, M., Kline, W. O., Stover, G. L., Bauerlein, R., et al. (2001). Akt/mTOR Pathway Is a Crucial Regulator of Skeletal Muscle Hypertrophy and Can Prevent Muscle Atrophy In Vivo. Nat. Cell. Biol. 3 (11), 1014–1019. doi:10.1038/ncb1101-1014

PubMed Abstract | CrossRef Full Text | Google Scholar

Bolotta, A., Filardo, G., Abruzzo, P. M., Astolfi, A., De Sanctis, P., Di Martino, A., et al. (2020). Skeletal Muscle Gene Expression in Long-Term Endurance and Resistance Trained Elderly. Int. J. Mol. Sci. 21 (11). doi:10.3390/ijms21113988

CrossRef Full Text | Google Scholar

Cocks, M., Shaw, C. S., Shepherd, S. O., Fisher, J. P., Ranasinghe, A., Barker, T. A., et al. (2016). Sprint Interval and Moderate-Intensity Continuous Training Have Equal Benefits on Aerobic Capacity, Insulin Sensitivity, Muscle Capillarisation and Endothelial eNOS/NAD(P)Hoxidase Protein Ratio in Obese Men. J. Physiol. 594 (8), 2307–2321. doi:10.1113/jphysiol.2014.285254

PubMed Abstract | CrossRef Full Text | Google Scholar

Egan, B., and Zierath, J. R. (2013). Exercise Metabolism and the Molecular Regulation of Skeletal Muscle Adaptation. Cell. Metab. 17 (2), 162–184. doi:10.1016/j.cmet.2012.12.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Esbjörnsson, M., Norman, B., Suchdev, S., Viru, M., Lindhgren, A., and Jansson, E. (2009). Greater Growth Hormone and Insulin Response in Women Than in Men during Repeated Bouts of Sprint Exercise. Acta Physiol. (Oxf) 197 (2), 107–115. doi:10.1111/j.1748-1716.2009.01994.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Esbjörnsson-Liljedahl, M., Bodin, K., and Jansson, E. (1985). Smaller Muscle ATP Reduction in Women Than in Men by Repeated Bouts of Sprint Exercise. J. Appl. Physiol. (1985) 93 (3), 1075–1083. doi:10.1152/japplphysiol.00732.1999

CrossRef Full Text | Google Scholar

Gibala, M. J., Little, J. P., van Essen, M., Wilkin, G. P., Burgomaster, K. A., Safdar, A., et al. (2006). Short-term Sprint Interval versus Traditional Endurance Training: Similar Initial Adaptations in Human Skeletal Muscle and Exercise Performance. J. Physiol. 575 (Pt 3), 901–911. doi:10.1113/jphysiol.2006.112094

PubMed Abstract | CrossRef Full Text | Google Scholar

Gillen, J. B., Percival, M. E., Skelly, L. E., Martin, B. J., Tan, R. B., Tarnopolsky, M. A., et al. (2014). Three Minutes of All-Out Intermittent Exercise Per Week Increases Skeletal Muscle Oxidative Capacity and Improves Cardiometabolic Health. PLoS One 9 (11), e111489. doi:10.1371/journal.pone.0111489

PubMed Abstract | CrossRef Full Text | Google Scholar

Hawley, J. A. (2004). Exercise as a Therapeutic Intervention for the Prevention and Treatment of Insulin Resistance. Diabetes Metab. Res. Rev. 20 (5), 383–393. doi:10.1002/dmrr.505

PubMed Abstract | CrossRef Full Text | Google Scholar

Hessvik, N. P., Boekschoten, M. V., Baltzersen, M. A., Kersten, S., Xu, X., Andersén, H., et al. (2010). LXRβ Is the Dominant LXR Subtype in Skeletal Muscle Regulating Lipogenesis and Cholesterol Efflux. Am. J. Physiology-Endocrinology Metabolism 298 (3), E602–E613. doi:10.1152/ajpendo.00553.2009

CrossRef Full Text | Google Scholar

Janssen, I., Heymsfield, S. B., Wang, Z. M., and Ross, R. (1985). Skeletal Muscle Mass and Distribution in 468 Men and Women Aged 18-88 Yr. J. Appl. Physiol. (1985) 89 (1), 81–88. doi:10.1152/jappl.2000.89.1.81

PubMed Abstract | CrossRef Full Text | Google Scholar

Kramer, H. F., and Goodyear, L. J. (1985). Exercise, MAPK, and NF-kappaB Signaling in Skeletal Muscle. J. Appl. Physiol. (1985) 103 (1), 388–395. doi:10.1152/japplphysiol.00085.2007

PubMed Abstract | CrossRef Full Text | Google Scholar

Kujach, S., Olek, R. A., Byun, K., Suwabe, K., Sitek, E. J., Ziemann, E., et al. (2019). Acute Sprint Interval Exercise Increases Both Cognitive Functions and Peripheral Neurotrophic Factors in Humans: The Possible Involvement of Lactate. Front. Neurosci. 13, 1455. doi:10.3389/fnins.2019.01455

PubMed Abstract | CrossRef Full Text | Google Scholar

Lindholm, M. E., Huss, M., Solnestam, B. W., Kjellqvist, S., Lundeberg, J., and Sundberg, C. J. (2014). The Human Skeletal Muscle Transcriptome: Sex Differences, Alternative Splicing, and Tissue Homogeneity Assessed with RNA Sequencing. FASEB J. 28 (10), 4571–4581. doi:10.1096/fj.14-255000

PubMed Abstract | CrossRef Full Text | Google Scholar

Lopes-Ramos, C. M., Chen, C.-Y., Kuijjer, M. L., Paulson, J. N., Sonawane, A. R., Fagny, M., et al. (2020). Sex Differences in Gene Expression and Regulatory Networks across 29 Human Tissues. Cell. Rep. 31 (12), 107795. doi:10.1016/j.celrep.2020.107795

PubMed Abstract | CrossRef Full Text | Google Scholar

Metcalfe, R. S., Babraj, J. A., Fawkner, S. G., and Vollaard, N. B. J. (2012). Towards the Minimal Amount of Exercise for Improving Metabolic Health: Beneficial Effects of Reduced-Exertion High-Intensity Interval Training. Eur. J. Appl. Physiol. 112 (7), 2767–2775. doi:10.1007/s00421-011-2254-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Miller, A. E. J., MacDougall, J. D., Tarnopolsky, M. A., and Sale, D. G. (1993). Gender Differences in Strength and Muscle Fiber Characteristics. Eur. J. Appl. Physiol. 66 (3), 254–262. doi:10.1007/bf00235103

CrossRef Full Text | Google Scholar

Moresi, V., Adamo, S., and Berghella, L. (2019). The JAK/STAT Pathway in Skeletal Muscle Pathophysiology. Front. Physiol. 10, 500. doi:10.3389/fphys.2019.00500

PubMed Abstract | CrossRef Full Text | Google Scholar

Pearson, A. M. (1990). Muscle Growth and Exercise. Crit. Rev. Food Sci. Nutr. 29 (3), 167–196. doi:10.1080/10408399009527522

PubMed Abstract | CrossRef Full Text | Google Scholar

Pointer, M. A., Harrison, P. W., Wright, A. E., and Mank, J. E. (2013). Masculinization of Gene Expression Is Associated with Exaggeration of Male Sexual Dimorphism. PLoS Genet. 9 (8), e1003697. doi:10.1371/journal.pgen.1003697

PubMed Abstract | CrossRef Full Text | Google Scholar

Richards, J. C., Johnson, T. K., Kuzma, J. N., Lonac, M. C., Schweder, M. M., Voyles, W. F., et al. (2010). Short-term Sprint Interval Training Increases Insulin Sensitivity in Healthy Adults but Does Not Affect the Thermogenic Response to Beta-Adrenergic Stimulation. J. Physiol. 588 (Pt 15), 2961–2972. doi:10.1113/jphysiol.2010.189886

PubMed Abstract | CrossRef Full Text | Google Scholar

Riuzzi, F., Sorci, G., Sagheddu, R., Chiappalupi, S., Salvadori, L., and Donato, R. (2018). RAGE in the Pathophysiology of Skeletal Muscle. J. Cachexia, Sarcopenia Muscle 9 (7), 1213–1234. doi:10.1002/jcsm.12350

CrossRef Full Text | Google Scholar

Rundqvist, H. C., Montelius, A., Osterlund, T., Norman, B., Esbjornsson, M., and Jansson, E. (2019). Acute Sprint Exercise Transcriptome in Human Skeletal Muscle. PLoS One 14 (10), e0223024. doi:10.1371/journal.pone.0223024

PubMed Abstract | CrossRef Full Text | Google Scholar

Saito, Y., Chikenji, T. S., Matsumura, T., Nakano, M., and Fujimiya, M. (2020). Exercise Enhances Skeletal Muscle Regeneration by Promoting Senescence in Fibro-Adipogenic Progenitors. Nat. Commun. 11 (1), 889. doi:10.1038/s41467-020-14734-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Salomoni, S., Soares, F. A., de Oliveira Nascimento, F. A., and da Rocha, A. F. (2008). Gender Differences in Muscle Fatigue of the Biceps Brachii and Influences of Female Menstrual Cycle in Electromyography Variables. Annu. Int. Conf. IEEE Eng. Med. Biol. Soc. 2008, 2598–2601. doi:10.1109/IEMBS.2008.4649732

PubMed Abstract | CrossRef Full Text | Google Scholar

Staron, R. S., Karapondo, D. L., Kraemer, W. J., Fry, A. C., Gordon, S. E., Falkel, J. E., et al. (1985). Skeletal Muscle Adaptations during Early Phase of Heavy-Resistance Training in Men and Women. J. Appl. Physiol. (1985) 76 (3), 1247–1255. doi:10.1152/jappl.1994.76.3.1247

CrossRef Full Text | Google Scholar

Stretch, C., Khan, S., Asgarian, N., Eisner, R., Vaisipour, S., Damaraju, S., et al. (2013). Effects of Sample Size on Differential Gene Expression, Rank Order and Prediction Accuracy of a Gene Signature. PLoS One 8 (6), e65380. doi:10.1371/journal.pone.0065380

PubMed Abstract | CrossRef Full Text | Google Scholar

Vendelbo, M. H., Christensen, B., Grønbæk, S. B., Høgild, M., Madsen, M., Pedersen, S. B., et al. (2015). GH Signaling in Human Adipose and Muscle Tissue during 'feast and Famine': Amplification of Exercise Stimulation Following Fasting Compared to Glucose Administration. Eur. J. Endocrinol. 173 (3), 283–290. doi:10.1530/eje-14-1157

PubMed Abstract | CrossRef Full Text | Google Scholar

Vestergaard, P. F., Vendelbo, M. H., Pedersen, S. B., Juul, A., Ringgard, S., Møller, N., et al. (2014). GH Signaling in Skeletal Muscle and Adipose Tissue in Healthy Human Subjects: Impact of Gender and Age. Eur. J. Endocrinol. 171 (5), 623–631. doi:10.1530/eje-14-0538

PubMed Abstract | CrossRef Full Text | Google Scholar

Vollaard, N. B. J., and Metcalfe, R. S. (2017). Research into the Health Benefits of Sprint Interval Training Should Focus on Protocols with Fewer and Shorter Sprints. Sports Med. 47 (12), 2443–2451. doi:10.1007/s40279-017-0727-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Warburton, D. E. R., Nicol, C. W., and Bredin, S. S. (2006). Health Benefits of Physical Activity: the Evidence. Can. Med. Assoc. J. 174 (6), 801–809. doi:10.1503/cmaj.051351

PubMed Abstract | CrossRef Full Text | Google Scholar

Welle, S., Tawil, R., and Thornton, C. A. (2008). Sex-related Differences in Gene Expression in Human Skeletal Muscle. PLoS One 3 (1), e1385. doi:10.1371/journal.pone.0001385

PubMed Abstract | CrossRef Full Text | Google Scholar

Willoughby, T. N., Thomas, M. P. L., Schmale, M. S., Copeland, J. L., and Hazell, T. J. (2016). Four Weeks of Running Sprint Interval Training Improves Cardiorespiratory Fitness in Young and Middle-Aged Adults. J. Sports Sci. 34 (13), 1207–1214. doi:10.1080/02640414.2015.1102316

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: muscle, KEGG, acute sprint exercise, hub genes, sex-biased genes

Citation: Nie M, Liu Q and Yan C (2022) Skeletal Muscle Transcriptomic Comparison Between Men and Women in Response to Acute Sprint Exercise. Front. Genet. 13:860815. doi: 10.3389/fgene.2022.860815

Received: 23 January 2022; Accepted: 27 May 2022;
Published: 12 July 2022.

Edited by:

Rinaldo Wellerson Pereira, Catholic University of Brasilia (UCB), Brazil

Reviewed by:

Daniela Caporossi, Foro Italico University of Rome, Italy
Kalyan C. Vinnakota, Ronin Institute, United States

Copyright © 2022 Nie, Liu and Yan. 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: Cheng Yan, eWFuY2hlbmd4eEAxNjMuY29t

These authors have contributed equally to this work

Disclaimer: 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.