Hypoxia treatment and resistance training alters microRNA profiling in rats skeletal muscle

MicroRNAs (miRNAs) may play a crucial regulatory role in the process of muscle atrophy induced by high-altitude hypoxia and its amelioration through resistance training. However, research in this aspect is still lacking. Therefore, this study aimed to employ miRNA microarray analysis to investigate the expression profile of miRNAs in skeletal muscle from an animal model of hypoxia-induced muscle atrophy and resistance training aimed at mitigating muscle atrophy. The study utilized a simulated hypoxic environment (oxygen concentration at 11.2%) to induce muscle atrophy and established a rat model of resistance training using ladder climbing, with a total intervention period of 4 weeks. The miRNA expression profile revealed 9 differentially expressed miRNAs influenced by hypoxia (e.g., miR-341, miR-32-5p, miR-465-5p) and 14 differentially expressed miRNAs influenced by resistance training under hypoxic conditions (e.g., miR-338-5p, miR-203a-3p, miR-92b-3p) (∣log2(FC)∣ ≥ 1.5, p < 0.05). The differentially expressed miRNAs were found to target genes involved in muscle protein synthesis and degradation (such as Utrn, mdm2, eIF4E), biological processes (such as negative regulation of transcription from RNA polymerase II promoter, regulation of transcription, DNA-dependent), and signaling pathways (such as Wnt signaling pathway, MAPK signaling pathway, ubiquitin-mediated proteolysis, mTOR signaling pathway). This study provides a foundation for understanding and further exploring the molecular mechanisms underlying hypoxia-induced rats muscle atrophy and the mitigation of atrophy through resistance training.

of chronic obstructive pulmonary disease 13,14 .Investigations have underscored alterations in miRNA expression within murine skeletal muscle subjected to hypoxic conditions, pinpointing their involvement in the regulation of AMPK-mediated energy metabolism 15 .However, a notable gap in the current body of research pertains to the paucity of knowledge concerning alterations in the miRNA expression profile during hypoxia-induced skeletal muscle atrophy.Several studies have reported that resistance training may elicit alterations in the expression of miRNAs within skeletal muscle [16][17][18] .Particular miRNAs, including miR-1-3p, miR-19b-3p, miR-92a, and others, are believed to participate in the regulatory mechanisms through which resistance training ameliorates agerelated muscle atrophy 19 .Despite this, no investigations have yet documented the shifts in miRNA expression during hypoxia-induced skeletal muscle atrophy or the potential mitigation of such atrophy through resistance training.Nevertheless, considering the biological processes under miRNA regulation, it is conceivable that miRNAs might function as intermediary regulatory factors influencing the impacts of hypoxia and resistance training on differential gene expression within skeletal muscle.
Therefore, this study initially established a rat model of hypoxia-induced skeletal muscle atrophy, attenuated by resistance training, followed by the screening of differentially expressed miRNAs in skeletal muscle.Subsequently, bioinformatics analysis methods were employed to analyze the potential biological functions of these differentially expressed miRNAs, with the aim of providing a foundation for further investigation into the molecular mechanisms underlying hypoxia-induced skeletal muscle atrophy and the mitigation of muscle wasting by resistance training in the future.

Animal model of hypoxia-induced skeletal muscle atrophy and resistance training regimen to mitigate muscle atrophy
The findings from the body composition analysis revealed notable disparities among the rat groups.Specifically, in the hypoxia control group (HC), distinct reductions were observed in the overall muscle content (12.49%), trunk muscle content (14.32%), and leg muscle content (10.06%), all of which were statistically significant when compared to the normoxia control group (NC) (p < 0.05).Conversely, within the hypoxia resistance training group (HR), there were substantial elevations noted in both overall muscle content (15.33%) and leg muscle content (24.42%) in comparison to the hypoxia control group (HC), with statistical significance achieved (p < 0.05).Comparing the hypoxia control group (HC) with the normoxia control group (NC), a substantial decrease of 19.35% in overall fat content (p < 0.01) and a significant reduction of 15.75% in leg fat content (p < 0.05) were evident, while no significant alteration emerged in trunk fat content (p > 0.05).It is noteworthy that the hypoxia resistance training group (HR) exhibited no significant changes in overall fat content, trunk fat content, or leg fat content when compared to the hypoxia control group (HC) (p > 0.05) (Fig. 1).
In comparison to the NC group, the gastrocnemius muscle and extensor digitorum longus muscle in the HC group of rats exhibited a considerable reduction in cross-sectional area (p < 0.05), whereas no significant difference was observed in the soleus muscle (p > 0.05).Conversely, when juxtaposed with the HC group, the HR group of rats displayed a substantial increase in the cross-sectional area of both the gastrocnemius muscle and extensor digitorum longus muscle (p < 0.01), while no significant disparity was detected in the soleus muscle (p > 0.05) (Fig. 2).

Bioinformatics analysis of differentially expressed miRNAs
By employing the miRWalk, miRanda, and TargetScan databases for target gene prediction and taking into account the convergence of outcomes across these three databases, a comprehensive compilation of 8,298 target genes emerged for the miRNAs exhibiting differential expression due to the influence of hypoxia on skeletal muscle.Similarly, in the context of hypoxia, the miRNAs displaying altered expression as a consequence of resistance training in skeletal muscle resulted in the projection of 14,275 target genes (Table 2).
GO analysis revealed that among the differentially expressed miRNAs between HC and NC, increased miR-NAs were enriched in biological processes such as small GTPase-mediated signal transduction, positive regulation of cell migration, and regulation of transcription DNA-dependent, totaling 686 terms.Increased miRNAs were also enriched in molecular functions including protein binding, DNA binding, and zinc ion binding, with a total of 205 terms.Additionally, increased miRNAs were enriched in cellular components such as cytoplasm, nucleus, and nucleolus, totaling 169 terms.On the other hand, decreased miRNAs were enriched in biological processes such as positive regulation of cell migration, positive regulation of transcription DNA-dependent, and positive regulation of transcription from RNA polymerase II promoter, totaling 942 terms.They were also enriched in molecular functions including protein binding, zinc ion binding, and molecular_function, totaling 277 terms.Furthermore, decreased miRNAs were enriched in cellular components such as nucleus, plasma membrane, and membrane, totaling 217 terms.
Regarding the differentially expressed miRNAs between HR and HC, increased miRNAs were enriched in biological processes such as regulation of neuronal synaptic plasticity, regulation of transcription DNA-dependent, and dephosphorylation, totaling 197 terms.They were also enriched in molecular functions including protein binding, histone acetyltransferase activity H4 K16 specific, and histone acetyltransferase activity H4 K8 specific, totaling 68 terms.Moreover, increased miRNAs were enriched in cellular components such as nucleus, cytoplasm, and nucleolus, totaling 76 terms.Conversely, decreased miRNAs were enriched in biological processes such as regulation of transcription DNA-dependent, negative regulation of transcription DNA-dependent, and nervous system development, totaling 1416 terms.They were also enriched in molecular functions including zinc ion binding, protein binding, and DNA binding, totaling 413 terms.Additionally, decreased miRNAs were enriched in cellular components such as cytoplasm, nucleus, and Golgi apparatus, totaling 296 terms.The top 15 enriched terms are presented in the Fig. 4.
The KEGG pathway analysis revealed that among the differentially expressed miRNAs between HC and NC, increased miRNAs were enriched in pathways such as Metabolic pathways, Pentose and glucuronate interconversions, and Leukocyte transendothelial migration, totaling 9 pathways.On the other hand, decreased miRNAs were enriched in pathways including the Wnt signaling pathway, Pathways in cancer, and Axon guidance, totaling 25 pathways.In the comparison between HR and HC, increased miRNAs did not exhibit enrichment in any specific pathways.However, decreased miRNAs were enriched in pathways such as the Wnt signaling pathway, Axon guidance, and Focal adhesion, totaling 33 pathways (Fig. 5).All GO and KEGG analysis results can be found in Supplementary Information 1-8.

Discussion
This study observed a decrease in skeletal muscle mass and reduced muscle fiber cross-sectional area in rats exposed to hypoxic conditions.Resistance training demonstrated some improvement in hypoxia-induced muscle atrophy.Hypoxia resulted in differential expression of 9 miRNAs in skeletal muscle (miR-341, miR-32-5p, miR-465-5p, miR-326-3p, miR-499-5p, miR-760-3p, miR-203a-3p, miR-205, and miR-1306-3p), while resistance training under hypoxic conditions induced differential expression of 14 miRNAs in skeletal muscle (miR-338-5p, miR-203a-3p, miR-92b-3p, miR-199a-5p, miR-152-3p, miR-497-5p, miR-499-5p, miR-130a-3p, miR-143-3p, miR-101b-3p, miR-222-3p, miR-139-3p, miR-140-5p, and miR-146b-5p).Bioinformatics analysis indicates that The hypoxic environment induced skeletal muscle atrophy in this study, manifested by decreased muscle mass and reduced muscle fiber cross-sectional area, consistent with findings from most previous studies.In conditions of low pressure and hypoxia, rats tend to exhibit slower weight gain.After a 5-week intervention in a simulated low-pressure and low-oxygen environment, akin to an altitude of 6000 m, the body weight, wet weight of the gastrocnemius muscle, and muscle fiber cross-sectional area of rats mirrored those of rats kept under normoxic conditions for 3 weeks 20 .In another study conducted in a decompression chamber to simulate low-pressure and low-oxygen conditions equivalent to an altitude of 4000 m (463 torr), it was noted that rats exposed to hypoxia for 12 weeks underwent a significant decrease in body weight.The soleus muscle and extensor digitorum longus muscle similarly demonstrated marked reductions in size.Moreover, this investigation established that the cross-sectional area of muscle fibers in the soleus muscle and extensor digitorum longus muscle decreased by 22% and 32%, respectively 21 .During a 5-week intervention in a decompression chamber simulating a lowoxygen environment equivalent to an altitude of 5000 m, a noteworthy reduction in the cross-sectional area of gastrocnemius muscle fibers was observed in rats.Remarkably, this reduction occurred concurrently in both types of muscle fibers 22 .
Under hypoxic conditions, resistance training in this study ameliorated hypoxia-induced skeletal muscle atrophy, indicating that appropriate protocols can promote muscle hypertrophy in such environments.Previous research has reported significant improvements in body weight and lean body mass before and after intervention with resistance training under normobaric hypoxia (14.4% oxygen concentration), showing a 6.1% increase in muscle cross-sectional area of the quadriceps, biceps, and hamstrings 23 .Similarly, exercise under both normoxic (FiO2 = 20.9%) and hypoxic (FiO2 = 12.7%) conditions demonstrated that resistance training significantly increased the thickness of the biceps and triceps brachii, stimulating skeletal muscle hypertrophy 24 .However, in this study, the muscle mass of the torso and whole body relative to body weight did not increase, possibly due to  www.nature.com/scientificreports/ the training protocol employed.Specifically, the resistance training involved climbing, which primarily focused on training the muscles of the limbs, resulting in relatively less stimulation of the muscles in the torso.This could also be one of the reasons for the less noticeable change in torso muscle mass.
Within the scope of this study, hypoxia precipitated the differential expression of 9 miRNAs.Some of these miRNAs have previously been linked to hypoxia stimulation or implicated in muscle atrophy.For instance, In rat models exposed to hypoxic conditions (9-11% oxygen concentration for 8 h daily over 4 weeks), miR-499 expression was noted to decrease in the gastrocnemius and biceps femoris muscles 25 .MiR-760-3p involvement in hypoxia-induced cell proliferation and apoptosis has been documented, with its expression diminished in hypoxia-exposed human pulmonary artery smooth muscle cells 26 .Further down the spectrum, miR-203a-3p's role appears tied to skeletal muscle differentiation.In zebrafish embryos, it regulates skeletal muscle differentiation by targeting Dmrt2a 27 .Notably, in a rat model simulating an altitude of 6000 m, miR-203a-3p expression in lung tissue was significantly downregulated after 72 h of intervention 28 .Similarly, within a rat ischemic injury model, miR-205 exhibited significant increases in skeletal muscle expression after ischemic periods and reperfusion 29 .The remaining four differentially expressed miRNAs in this study-miR-341, miR-465-5p, miR-326-3p, and miR-1306-3p-have yet to be linked to low oxygen intervention or muscle atrophy.Their precise roles in mediating low oxygen-induced skeletal muscle atrophy require further exploration.
In hypoxic conditions, resistance training induced distinctive miRNA expression patterns in skeletal muscle, involving 14 miRNAs.Among this cohort of miRNAs, studies have illuminated the role of miR-199a-5p in Duchenne muscular dystrophy (DMD), where it influences the WNT signaling pathway via genes such as FZD4, JAG1, and WNT2, thereby impacting cellular proliferation and myogenic differentiation.Elevated miR-199a-5p expression has also been linked to aberrant muscle fiber degradation 30 .MiR-143-3p upregulation in the blood of amyotrophic lateral sclerosis (ALS) patients marks it as a potential miRNA biomarker for muscle atrophy 31 .miR-146b-5p has emerged as a pivotal player; its elevated expression in skeletal muscle samples of type 2 myotonic dystrophy (DM2) patients suggests its potential as a miRNA biomarker for muscle atrophy 32 .In the context of miR-222-3p and skeletal muscle, its role extends to energy metabolism 33 .Additional miRNAs, including miR-92b-3p, miR-152-3p, and miR-140-5p, have received less attention in the context of skeletal muscle growth, development, and physiology.
Virtually all genes hold the potential to be under miRNA regulation, with each miRNA capable of targeting hundreds of genes 34 .Throughout gene expression, miRNAs have the ability to curtail mRNA transcription into proteins.This is achieved by targeting the 3'UTR of genes.Additionally, miRNAs can bind to mRNA, thereby inducing mRNA degradation, culminating in a reduction in gene expression levels.Upon delving into www.nature.com/scientificreports/ the predicted target genes, it becomes evident that the miRNAs are differentially expressed due to hypoxia, incline target genes linked to atrophy.For example, Utrn (Utrophin) is one of the target genes of miR-1306-3p.This gene, akin in structure and function to the dystrophin gene (which encodes the anti-atrophy protein), is also recognized as a dystrophin-associated protein.Notably, it contributes to the regulatory machinery governing certain facets of atrophy 35,36 .Murine double minute 2 (Mdm2) has emerged as a direct target gene of miR-205.This gene encodes an E3 ubiquitin ligase that resides within the nucleus.Interestingly, Mdm2 showcases heightened expression within a hindlimb suspension-induced muscle atrophy model 37 .Furthermore, Ube1, Rbx1, and other elements also count among the target genes influenced by miR-205.These genes actively participate in the ubiquitin-mediated protein degradation process 38,39 , which potentially extends to the intricate regulation of skeletal muscle mass.Beyond the aforementioned miRNAs, various other differentially expressed miRNAs are likewise associated with target genes central to muscle atrophy regulation.Notably, miR-326-3p target genes include Stk11, Ikbkb, and Pi3k, while miR-465-5p target genes include Aire, Apc8, and Rbx1.These miRNAs could influence muscle atrophy by overseeing the degradation of skeletal muscle proteins.
In the KEGG pathway analysis, it is evident that upregulated miRNAs and downregulated miRNAs target different signaling pathways.It is noteworthy that within these pathways, the upregulated miRNAs target the Wnt signaling pathway, MAPK-TGF beta signaling pathway, while the downregulated miRNAs target the Wnt signaling pathway and GnRH signaling pathway.Additionally, hypoxia-induced resistance training may target signaling pathways such as the Wnt signaling pathway, Insulin signaling pathway, MAPK signaling pathway, TGF beta signaling pathway, Ubiquitin-mediated proteolysis, and mTOR signaling pathway, which have been reported to potentially influence skeletal muscle quality.In the current study, bioinformatics predictions pinpointed 7 miRNAs, including miR-1306-3p, miR-203a-3p, and miR-205, as potential regulators of the Wnt signaling pathway.This points toward their plausible role as regulatory factors implicated in the low oxygen-induced effects on skeletal muscle quality.Studies reported that under hypoxic stimulation, HIF1α negatively regulates skeletal muscle regeneration by inhibiting the classical Wnt signaling pathway 40 .The skeletal muscle TGF-β signaling pathway holds sway over muscle growth and development.TGF-β exerts its regulatory role in muscle atrophy by activating both the classical Smad-dependent pathway and nonclassical routes, such as ERK1/2, JNK1/2, and p38 MAPK.The engagement of Smad3, ERK1/2, and JNK1/2 is integral to TGF-β's capacity to induce muscle atrophy 41 .In a mouse model of skeletal muscle atrophy induced by amyotrophic lateral sclerosis (ALS) (hSOD1G93A mice), the mRNA levels of TGF-β1, total Smad3 protein levels, and p-Smad3 displayed marked increases 42 .Similarly, augmented TFGβ signaling has been discerned in congenital myopathies 43 .Transcriptional and methylation analyses of individuals undergoing acute and chronic resistance exercise interventions revealed that the manipulation of three genes, FOS, SMAD3, and WNT9A, within the TGF-β signaling pathway impacted skeletal muscle synthesis, metabolism, and hypertrophy in response to resistance training 44 .In a rat resistance training model, the TGF-β1-Smad signaling pathway mediated the muscle hypertrophy elicited by resistance training.This mechanism encompasses a substantial reduction in skeletal muscle TGF-β1 expression and phosphorylation of Smad3 COOH-terminal residues (P-Smad3 S423/425) following resistance training.This results in heightened phosphorylation levels of SOH2 (P-Smad2-LS245/250/255) and Smad3 (P-Smad3-LSer208), culminating in a diminished inhibition of satellite cell activation and protein synthesis, ultimately fostering hypertrophy 45 .Furthermore, differentially expressed miRNAs can also target the ubiquitin-mediated protein degradation signaling pathway (such as miR-203a-3p, miR-205, and miR-32-5p) and the mTOR signaling pathway (such as miR-32-5p, miR-326-3p, and miR-760-3p).
This study also has some limitations.Firstly, the utilization of microarray technology in our investigation presents inherent biases and limitations compared to RNAseq methods.While microarray analysis allows for comprehensive profiling of miRNA expression, it is acknowledged that RNAseq possesses certain advantages such as improved sensitivity and the ability to detect novel transcripts.We recognize the necessity to elucidate these methodological differences and their potential impacts on our findings.Moreover, an important consideration pertains to the translatability of our results from rodent models to human physiology.Although our study extensively discusses mechanisms and pathways observed in rodent models, the direct extrapolation of these findings to human biology requires cautious interpretation.We acknowledge the complexities inherent in inter-species variations and biological contexts.Exploring these distinctions in detail is essential to accurately evaluate the significance and applicability of our findings to human conditions.Finally, experimental validation is necessary to confirm the roles of signaling pathways identified by bioinformatics in hypoxia-induced muscle atrophy.

Conclusions
Resistance training has shown potential in ameliorating hypoxia-induced skeletal muscle atrophy in rats, with miRNAs emerging as key regulators in this context.Differentially expressed miRNAs exert regulatory influence over both the progression of hypoxia-induced muscle atrophy and its amelioration through resistance training in rats, potentially operating through a multitude of signaling pathways.

Animals and grouping
This study was performed in accordance with relevant guidelines and regulations.Animal experiments were carried out in accordance with the ARRIVE guidelines and approved by the Ethics Committee for Sports Science Experiments at Beijing Sport University (approval number: 2023043A).Male Sprague Dawley (SD) rats, aged eight weeks, were procured from Beijing Vital River Laboratory Animal Technology Co., Ltd.Following a one-week acclimatization period, body composition was evaluated using dual-energy X-ray absorptiometry (DEXA).Subsequently, the rats were allocated randomly into four groups, each comprising 10 rats: the normoxic control group (NC), the normoxic resistance exercise group (NR), the hypoxic control group (HC), and www.nature.com/scientificreports/ the hypoxic resistance exercise group (HR).The rats were individually housed in cages under tightly controlled environmental conditions, maintaining a temperature range of 22-25 °C, relative humidity of 50-80%, and a 12-h light cycle from 7:00 AM to 7:00 PM.They were provided with standard rodent maintenance feed and had unrestricted access to food.Daily recordings of body weight and food intake were recorded.

Intervention protocol and sample collection
Following a week of adaptive training, the formal experimental intervention commenced.Adaptation phase was conducted under normoxic conditions.During the adaptation phase, rats underwent ladder climbing training without the addition of weights.For the formal intervention, rats in the HC and HR groups were exposed to a controlled oxygen concentration of 11.2%, simulating hypoxic conditions equivalent to an altitude of 4500 m.This hypoxic exposure was maintained for four weeks.In NR and HR, resistance exercise was conducted using a ladder-climbing simulation apparatus, with ladder dimensions measuring 110 cm × 18 cm and a step spacing of 2 cm.The resistance training protocol adhered to previous research practices and included the following steps: starting with a load equivalent to 30% of their body weight for the initial training session, followed by a gradual increase of 15% of their body weight for each subsequent session until reaching a load of 130% of their body weight.This load was then maintained for the remainder of the intervention period.Training sessions occurred every alternate day, comprising three sets of five repetitions per session, with a rest interval of 3-5 min between sets.This training regimen was conducted consistently over four weeks.The rats in the NR group underwent resistance training in a normoxic environment, while those in the HR group underwent resistance training in a hypoxic environment (oxygen concentration of 11.2%).Sample collection occurred 72 h after the intervention.Rats were humanely euthanized under anesthesia, and blood was drawn from the abdominal aorta.Skeletal muscles, including the gastrocnemius (GAS), extensor digitorum longus (EDL), and soleus (SOL), were meticulously isolated.Muscle samples were cleansed in physiological saline to eliminate residual blood, excess moisture was absorbed with filter paper, and the weight of wet muscle was measured.Skeletal muscle tissues were then segregated and treated as follows: muscles designated for histological examination were fixed in a 4% paraformaldehyde solution; muscles intended for microarray analysis were preserved in RNAlater solution, initially stored at 4 °C, then shifted to -20 °C, and finally transferred to an ultralow temperature freezer at − 80 °C.

Body composition analysis
Rat body composition was assessed using a dual-energy X-ray absorptiometry (DXA) scanner (Lunar iDXA, GE Healthcare, USA) both before the adaptation phase and upon completion of the hypoxic training intervention.Rats were anesthetized with a sodium pentobarbital solution at a concentration of 3 mg/ml, administered at a rate of 1 ml/100 g of body weight.For the groups subjected to hypoxic intervention, anesthesia was administered within a low-oxygen environment, while the remaining rats were anesthetized under normoxic conditions.Body composition measurements were conducted on the anesthetized rats.

Hematoxylin-eosin staining
Following fixation, skeletal muscle tissue was cut into cubes measuring 0.5 cm × 0.5 cm × 0.5 cm.These cubes were positioned to ensure that the muscle fibers were perpendicular to the section before being embedded.Paraffin sections with a thickness of 5 μm were subsequently prepared.To examine histological alterations in the skeletal muscle tissue, hematoxylin and eosin staining were applied.From each slide, three fields of view were selected, and within each field of view, five images were captured.The cross-sectional area of the skeletal muscle was quantified utilizing ImageJ software.

RNA isolation and small RNA profiling
RNA extraction was carried out from 50 mg of gastrocnemius muscle samples using the miRNeasy Mini Kit (Qiagen, Germany).The concentration and purity of the extracted RNA were determined using a Nanodrop 2000 spectrophotometer, ensuring an OD260/OD280 ratio within the range of 2.0 to 2.2.A 1.5% formaldehyde denaturing gel electrophoresis was performed at 120 V for 15 min, followed by gel imaging to evaluate the integrity of the 28S and 18S ribosomal RNA bands.
For miRNA expression profiling in the gastrocnemius muscle, the Affymetrix miRNA 4.0 microarray (Affymetrix, USA) was employed.This microarray comprises 30,424 miRNA probe sets, encompassing 728 mature rat miRNAs and 490 pre-miRNAs.The labeling of samples followed the standard procedure (Flash Tag Biotin HSR RNA Labeling Kit, Affymetrix, USA).Subsequently, the samples were subjected to microarray hybridization (GeneChip2 Hybridization Oven 640, Affymetrix, USA), followed by washing (GeneChip2 Fluidics Station 450, USA) and finally scanning (GeneChip2 Scanner 3000 7G, Affymetrix, USA).Signal intensity values were evaluated using Affymetrix Expression Console™ software (version 1.3.1).The experimental outcomes were analyzed and processed, and the sample expression data were subjected to normalization procedures.

Differential expression miRNA bioinformatics analysis
For predicting target genes, the miRWalk platform (http:// mirwa lk.umm.uni-heide lberg.de/) was utilized, incorporating three target gene prediction databases: miRWalk, miRanda, and TargetScan.Genes that were consistently predicted by all three target gene databases were chosen as the target genes for subsequent analysis.The miRWalk platform was also utilized for conducting Gene Ontology (GO) enrichment analysis and KEGG pathway analysis of the miRNAs [46][47][48] .Additionally, Cytoscape version 3.10.0was employed to construct an interactive regulatory network illustrating the interactions between miRNAs and KEGG signaling pathways.

Data analysis
The data were analyzed using SPSS v22.0 software (SPSS Inc., Chicago, IL, USA).Results are presented as mean ± standard deviation (mean ± SD).To compare multiple groups, a two-way analysis of variance (ANOVA) was conducted.In case of a significant interaction effect between the two factors, simple effects tests were performed.For post hoc comparisons in the absence of an interaction effect, the Tukey method was employed instead of the least significant difference (LSD) method to minimize potential multiple testing bias.Statistical significance was set at P < 0.05.Statistical analysis was performed using a t-test to evaluate the alterations in miRNA expression between the two groups.The p-values were adjusted using the Benjamini-Hochberg correction method, with an FDR < 0.05 threshold applied.Additionally, |log2(FC)| ≥ 1.5 was used as the criterion for filtering differentially expressed miRNAs.

Institutional review board statement
This study was performed in accordance with relevant guidelines and regulations.Animal experiments were carried out in accordance with the ARRIVE guidelines and approved by the Ethics Committee for Sports Science Experiments at Beijing Sport University (approval number: 2023043A).

Figure 1 .
Figure 1.Changes in Rat Body Composition After Intervention (* and ** indicate p < 0.05 and p < 0.01 compared to the NC group, respectively; # and ## indicate p < 0.05 and p < 0.01 compared to the HC group, respectively).

Figure 2 .Figure 3 .
Figure 2. Cross-Sectional Area of Gastrocnemius Muscle, Soleus Muscle, and Extensor Digitorum Longus Muscle in Different Groups.(A) HE staining and cross-sectional area of muscle fibers in GAS.(B) HE staining and cross-sectional area (CSA) of muscle fibers in the SOL.(C) HE staining and cross-sectional area of muscle fibers in EDL (GAS: gastrocnemius muscle, SOL: soleus muscle, EDL: extensor digitorum longus muscle.Scale bars = 25 µm.* and ** indicate p < 0.05 and p < 0.01 compared to the NC group, respectively; # and ## indicate p < 0.05 and p < 0.01 compared to the HC group, respectively).

Figure 4 .
Figure 4. GO Analysis of Differentially Expressed miRNAs in Skeletal Muscle.(A) GO analysis of downregulated miRNA expression af-fected by hypoxia.(B) GO analysis of up-regulated miRNA expression affected by hypoxia.(C) GO analysis of down-regulated miRNAs affected by resistance training under hypoxia.(D) GO analysis of up-regulated miRNAs affected by resistance training under hypoxia.

Table 1 .
Differentially Expressed miRNAs in Rat Skeletal Muscle.

Table 2 .
Predicted Target Genes of Differentially Expressed miRNAs in Skeletal Muscle.