Introduction

The World Health Organization has reported that obesity has reached epidemic levels, having tripled in the last 50 years. Western-type diets that are high in fat and sugar often lead to obesity and associated metabolic syndrome. Metabolic syndrome is characterized by the presence of a cluster of conditions which includes abdominal obesity, hypertension, dyslipidemia, and hyperglycemia, all of which put a person at a greater risk for cardiovascular disease1. A less well-known risk of excessive fat gain is sarcopenic obesity. Sarcopenic obesity is defined by low muscle mass and function coupled with excessive fat mass2. As muscle mass has been found to be positively associated with longevity3 and other measures of quality of life4,5, it is important to understand the relationship between diet and muscle integrity.

Myosteatosis, or fat infiltration in muscle tissue, is poorly characterized, yet is an important feature of sarcopenic obesity. Fatty infiltrations in muscle tissue in sarcopenic obesity can further exacerbate the symptoms of metabolic syndrome due to the relative loss of muscle mass to fat gain6,7. Studies have found intramuscular lipid accumulation to be associated with dyslipidemia, hyperglycemia, and hypertension in both men and women6,8. Insulin resistance is commonly increased as well, as muscle is the largest organ responsible for glucose regulation9,10,11. Furthermore, proinflammatory cytokines that circulate in the obese state can further disrupt the regenerative capacity of skeletal muscle10. Whereas the muscle undergoes regular remodeling triggered by inflammatory signaling12,13,14, the environment created by local cytokines in obesity facilitates conditions of chronic injury leading to irregular remodeling and pathophysiologic fibrosis15. Furthermore, intramuscular fat infiltration and fibrosis are associated with decreased skeletal muscle contraction which can lead to a loss of muscle strength and physical function4,5,16,17,18,19. These conditions, in turn, are thought to encourage a sedentary lifestyle that further promotes weight gain and obesity.

To study fat infiltration in muscle tissue, a diet-induced obesity model in male Sprague–Dawley rats was previously developed. It has been demonstrated that a high fat/high sucrose (HFS) diet induces obesity and low-grade systemic and local muscle inflammation in this pre-clinical model20,21. The HFS diet led to an increase in fat infiltration and fibrosis in the primarily fast twitch fibred vastus lateralis (VL) muscle in as early as three days following the start of a metabolic challenge20. However, the oxidative capacity of the soleus muscle (~ 100% slow twitch fibred in Sprague Dawley rats)22 seemed to protect against changes in muscle integrity23.

Pre-clinical studies on obesity-induced metabolic disease and associated changes in muscle health have used male animals, as males are typically thought to have a steady metabolic profile, while females undergo periodic cycles of metabolic profiles associated with their pup-bearing function. Up to this point, studies looking at obesity related musculoskeletal degeneration in rats have only been performed in males. As females can differ in the presentation and mechanism behind various diseases24,25, it is important to understand the impact that sex has on musculoskeletal health. Furthermore, studying male and female rats may provide unique opportunities to understand the different relationships that may exist between metabolic responses to an obesogenic diet and corresponding musculoskeletal adaptations and degenerations.

Therefore, the purpose of the present study was to determine if previous findings from the model of diet-induced obesity in male rats can be generalized to both sexes, or whether sex affects the response to the HFS diet, as indicated by markers of metabolic syndrome and changes in muscle integrity. We assessed the VL and soleus muscle in male and female Sprague Dawley rats after a 12-week HFS diet intervention using muscle histology and a quantitative proteomics analysis as our primary outcomes. It was hypothesized that the structural changes in the muscles associated with the HFS diet intervention would be similar in the male and female rats. However, we detected sex-dependent and muscle-specific differences in intramuscular fat infiltration between male and female rats receiving the obesogenic diet.

Methodology

Experimental design

The current study is reported in accordance to ARRIVE guidelines. Twenty-four male and twenty female Sprague–Dawley rats were housed in pairs on a 12-h dark/light cycle. All animals were purchased from the University of Calgary’s Life and Environmental Science Animal Research Centre where they were bred and raised according to Canadian guidelines. Within each sex, rats were either allocated to a HFS diet group (39.2% of total energy as fat, 43.5% sucrose, 15.8% protein; custom Diet #102412, Dyets, Inc, n = 12 males, n = 12 females) or to a chow control diet (Chow, 13.4% of total energy as fat, 3.8% sucrose, 29.8% protein, Lab Diet 5001, n = 12 males, n = 8 females). All experiments were approved by the University of Calgary Life and Environmental Sciences Animal Care Committee, and experiments were conducted in accordance with the Canadian Council on Animal Care standards. The animals were 12 weeks old at the start of the ad libitum diet intervention. All animals remained on the diet for 12 weeks and were euthanized at 24 weeks of age. Animals were anesthetized with isoflurane and then euthanized by a cardiac blood draw and aortic snip. Prior to euthanasia, body mass and composition were measured using Dual Energy X-ray Absorptiometry (DXA) and analyzed with software for small animal analysis (GE Lunar DXA).

Muscle characterization

The VL and soleus muscles were harvested and flash-frozen in liquid nitrogen and stored at minus 80 °C until analysis. Myosin heavy chain (MHC) isoforms from both muscles were identified using SDS-PAGE with a 7.5% acrylamide separating and 4.5% stacking gel as previously described by Joumaa et al26. The electrophoresis was performed using a Bio-Rad (USA) Mini-Protean unit with 0.75 mm spacer plates at a constant voltage of 73 V for 18 h followed by a constant current of 1.5 mA/gel for 23.5 h at 4 °C. Gels were stained with Coomassie brilliant blue and imaged using a GS-800 Calibrated Densitometer (Bio-Rad, USA). Quantification of the band optical density was performed using ImageJ27, and proportions of MHC IIa, IIx, IIb, and I were determined in terms of total MHC content of the muscle.

Cytokine, adipokine, and blood lipid profile measurements

Following an overnight fast, blood was collected immediately following anesthesia via cardiac puncture. Samples were centrifuged at 7000 rpm for 15 min, and two serum aliquots stored at -80°C until analysis. One sample was used to quantify twenty-seven serum cytokines and adipokines using a Rat Cytokine/Chemokine 27-Plex Discovery Assay Array (Eotaxin, EGF, Fractalkine, IL-1α, IL-1 β, IL-2, IL-4, IL-5, IL-6, IL-10, IL-12(p70), IL-13, IL-17A, IL-18, IP-10/CXCL10, GRO/KC, IFN-γ, TNF-α, G-CSF, GM-CSF, MCP-1, leptin, LIX, MIP-1α, MIP-2, RANTES, vascular endothelial growth factor (VEGF), Eve Technologies, Calgary, AB). The other serum sample was analyzed for lipid profile (total cholesterol, LDL-cholesterol, HDL-cholesterol, and triglycerides) using colorimetric assays (Calgary Lab Services, Calgary, AB, Canada).

Insulin sensitivity analysis

One week before sacrifice, animals underwent an insulin sensitivity test. Animals were fasted for 14 h prior to blood collection. Blood was collected via tail nick in a chilled tube for insulin analysis. Blood glucose was measured immediately with a blood glucose meter (OneTouch Verio and Blood Glucose Monitoring System, Lifescan, Switzerland). Insulin resistance was calculated using the HOMA-IR index which uses the fasted blood glucose concentration and insulin levels (homeostatic model assessment for insulin resistance—HOMA-IR28. A value higher than 2.5 is indicative of increased insulin resistance in humans.

Histological staining procedures

Slide preparation

A 3 mm mid-belly slice of the previously flash-frozen muscle was cut, mounted on OCT compound and sectioned at 10 μm using a Leica cryostat at − 20 °C, dried at room temperature, and then stored in − 20 °C until used. Oil Red O (ORO) staining for intramuscular fat, Picrosirius red staining for collagen/fibrosis, and CD68+ immunohistochemistry was completed on one slide per staining procedure.

Oil red O staining

As previously described by Collins et al. (2016), ORO stock solution was made by adding 2.5 g Oil Red O (Sigma-Aldrich, Oakville, ON, Canada) to 400 mL of 99% isopropanol (Fisher Scientific, Ottawa, ON, Canada). ORO working solution was then prepared from a ratio of three parts ORO stock: two parts distilled water and was filtered using a course filter, then was refrigerated at 4 °C for 10 min before use. Sections were first washed quickly in distilled water and then quickly dipped in 60% isopropanol and stained for 20 min in the ORO working solution. Slides were then rinsed first in 60% isopropanol and then in distilled water, and finally cover-slipped with glycerol jelly.

Picrosirius red staining

Slides were dried for 2 weeks in a 38 °C oven, fixed overnight in neutral buffered formalin (NBF) and then immersed in Bouin’s fixative solution for 1 h at 60 °C in a water bath, and then washed in distilled water. The Picrosirius red solution was made by adding 0.1 g of Sirius red stain (Sigma-Aldrich), to 100 mL of a saturated aqueous picric acid (1.3% in water, Sigma-Aldrich). Sections were left in this solution to stain for 1 h at room temperature and then rinsed in 2 changes of 1% acetic acid, dehydrated in three changes of 100% ethanol, cleared in 2 changes of xylene, and mounted with Cytoseal 60 mounting media (Collins et al. 2016).

Stain quantification

ORO and Picrosirius red stained sections were imaged at 40× magnification with an Olympus DP73 camera mounted on an Olympus BX53, and the corresponding images were analyzed using the software ImageJ (Rasband, WS, ImageJ, U.S. National Institutes of Health, Bethesda, Maryland, USA). Images were taken of the entire muscle section and were quantified by one reviewer who was unaware of which diet group each muscle section belonged to. The relative percentage of intramuscular fat and fibrosis was determined per cross sectional area of the muscle.

Immunofluorescence

Immunofluorescence was performed as described previously29. In brief, slides were first fixed in cold acetone for 30 min, quenched for 20 min in 3% hydrogen peroxide, and then blocked using 2.5% normal horse serum (Vector Laboratories, Burlingame, CA). Sections were then incubated overnight at 4 °C with the primary antibody, mouse anti-rat CD68+, a marker for inflammatory macrophages, at a dilution ratio of 1:200. The secondary antibody, IMPRES anti-mouse IgG peroxidase (IMPRESS kit, Biolynx Inc) was applied the next morning for 30 min at room temperature (IMPRESS kit, Vector Laboratories). Slides were then incubated for 20 min in a Tyramide Plus Cyanine 3 Amplification system from AKOYA Biosciences. To visualize nuclei that were present in the muscle section, DAPI was applied, and sections were mounted in Vectashield mounting media for fluorescence (Vector Laboratories). All washes were performed in Phosphate buffered saline (PBS). Ten images per soleus cross-section and 20 images per VL cross-sections were obtained using an Olympus BX51 fluorescent microscope, with appropriate filters for cy3 and DAPI, and attached to an Olympus DP28 camera.

Statistical analysis

All statistical analyses were performed using R statistical software (v. 4.2.1, R Core team). Non-parametric Mann–Whitney U tests were performed on data pertaining to body composition and markers of metabolic syndrome to determine a difference between the chow and HFS animals for both sexes. A two-way ANOVA was performed to determine the main effects of Diet and Sex and their interaction, Diet × Sex, on MHC isoform distributions and histological measures. If a significant interaction effect was found, a one-way ANOVA with Tukey’s post-hoc test was performed to determine differences between groups. Sex-dependent effects were not looked at for the serum cytokines because of unusual results from the male-chow animals. Instead, Mann–Whitney U tests were performed to determine differences between the chow and HFS animals for both sexes. Least-square linear regressions were used to determine if there was a relationship between body fat percentage and intramuscular fat percentage in both the male and female VL and soleus.

Proteomics

Shotgun proteomic analysis

Shotgun proteomics were carried out as described previously by30,31,32 and was performed on the previously described flash-frozen VL and soleus muscle tissue that was stored at -80 °C after sacrifice of the animals. Muscle samples were added to a lysing tube with 500 uL of lysis buffer (2% SDS, 200 mM ammonium bicarbonate, protease inhibitor tablets) and 4–5 small stainless steel homogenizing beads. Samples were homogenised at 30 Hz for 20 min in the homogenizer and then sonicated with a Fisherbrand™ Model 120 Sonic Dismembrator. After sonication, samples were centrifuged at 14,000×g for 15 min at 4 °C to pellet unwanted debris and the supernatant fraction transferred into a protein lobind Eppendorf tube. Total protein concentrations for each sample were then determined using a ThermoFisher BCA kit.

Samples for the proteomic analysis were prepared using the filter-assisted separation of peptides (FASP) method. Initially, 100 μg of protein was precipitated by adding trichloroacetic acid (TCA) followed by an incubation on ice. Samples were then centrifuged at 14,000 ×g for 15 min at 4 °C, washed 3 times in ice cold acetone, and stored at − 20 °C. Samples were resuspended in 8 M tris-urea solution by shaking and then denatured with the addition of 10 mM DTT at 37 °C for 30 min. 50 mM iodoacetamide was added in the dark at room temperature to complete carbamidomethyl modification of the cysteines. Samples were moved to the top of a 30 kDa filter along with 100 uL of wash solution and were then centrifuged at 14,000×g for 15 min. Samples were then washed 3 times with 8 M tris-urea and 3 times with 50 mM ammonium bicarbonate. The samples were trypsinized at 37 °C overnight at a ratio of 1:10 trypsin: total protein, and subsequently, eluted off the filter membrane by washing 3 times with 50 mM ammonium bicarbonate. Finally, samples were subjected to a c18 clean up with Waters solid-phase extraction (SPE) column, according to the manufacturer’s directions, before injection on the orbitrap mass spectrometer in DDA mode.

High performance liquid chromatography (HPLC) and mass spectrometry

Samples were then subjected to liquid chromatography and tandem mass spectrometry at the Southern Alberta Mass Spectrometry Facility at the University of Calgary. Using a process previously described by31,32,33,34, tryptic peptides were analyzed on an Orbitrap Fusion Lumos Tribrid mass spectrometer (Thermo Scientific) operated with Xcalibur (version 4.4.16.14) and coupled to a Thermo Scientific Easy-nLC (nanoflow Liquid Chromatography) 1200 system. Tryptic peptides (2 μg) were loaded onto a C18 trap (75 µm × 2 cm; Acclaim PepMap 100, P/N 164946; ThermoScientific) at a flow rate of 2 μL/min of solvent A (0.1% formic acid in LC–MS grade water). Peptides were eluted using a 120 min gradient from 5 to 40% (5% to 28% in 105 min followed by an increase to 40% B in 15 min) of solvent B (0.1% formic acid in 80% LC–MS grade acetonitrile) at a flow rate of 0.3 μL/min and separated on a C18 analytical column (75 μm × 50 cm; PepMap RSLC C18; P/N ES803; ThermoScientific). Peptides were then electrosprayed using 2.1 kV voltage into the ion transfer tube (300°C) of the Orbitrap Lumos operating in positive mode. The Orbitrap first performed a full MS scan at a resolution of 120,000 FWHM to detect the precursor ion having a m/z between 375 and 1,575 and a + 2 to + 7 charge. The Orbitrap AGC (Auto Gain Control) and the maximum injection time were set at 4 × 105 and 50 ms, respectively. The Orbitrap was operated using the top speed mode with a 3 s cycle time for precursor selection. The most intense precursor ions presenting a peptidic isotopic profile and having an intensity threshold of at least 5,000 were isolated using the quadrupole and fragmented with HCD (30% collision energy) in the ion routing multipole. The fragment ions (MS2) were analyzed in the ion trap at a rapid scan rate. The AGC and the maximum injection time were set at 1e4 and 35 ms, respectively, for the ion trap. Dynamic exclusion was enabled for 45 s to avoid the acquisition of same precursor ion having a similar m/z (plus or minus 10 ppm).

Proteomic data and bioinformatics analysis

Spectra data were matched to peptide sequences from a Rat FASTA reference file obtained from Uniprot on Sept 23rd, 2020 using MaxQuant (v.2.4.9.0). Variable modifications included Oxidation on Methionine, Acetylation on any N-terminus, and Deamidation on Asparagine and Glutamine. Samples were re-quantified and matched between runs. The first search peptide tolerance was set to 10, the minimum peptide length was set to 5, and the maximum mass was set to 6,600 Da, and the minimum score for the modified peptides was set to 20. All other settings on the MaxQuant were set to default.

Statistical analysis was performed on MSstatsShiny UI (v0.1.0; https://msstatsshiny.com/app/MSstatsShiny) running MSstats (v.4.2.0). The MaxQuant files “ProteinGroups” and “evidence” were used for the analysis, while the annotation file was created according to directions provided by MSstats35. Unique peptides were included in the analysis. For data processing, log2 was defined for data transformation and equalize medians was the method chosen for normalization. The options “Use all features”, “assume all NA as censored”, “Do not apply cutoff to censor missing values”, and “Model based imputation” were enabled. Statistical inference was performed using non-paired t-tests, and the data were imported to excel for final assessment. Significant enriched proteins for each group were found by separating upregulated and downregulated proteins on a log twofold change, where positive values were upregulated, and negative were downregulated. A cut-off using a p-value of 0.05 was used to determine significance. Metascape analysis was used to identify enriched pathways36. Protein–protein interactions were encoded into networks using the metascape website (https://metascape.org/), and enriched pathways were plotted as heatmaps.

Results

Body composition and markers for metabolic syndrome

Animals of both sexes, fed the HFS diet demonstrated a statistically higher percentage of body fat and more body mass than did the chow control animals (Fig. 1A,B, p < 0.001). With respect to the blood lipid profile, both male and female animals had a significant increase in blood triglycerides while on the HFS diet compared to the control group animals on the chow diet (Fig. 1C, p  = 0.008 and p = 0.003, respectively). However, only the female animals on the HFS diet had a significant increase in total cholesterol and HDL cholesterol compared to the chow-fed animals (Fig. 1D,E, p  = 0.022 and p = 0.047, respectively), while the males exhibited no change in total and HDL cholesterol (Fig. 1D,E, p  = 0.95 and p = 0.62, respectively). Both male and female HFS diet animals had statistically higher HOMA-IR indexes than their sex-matched chow-fed controls (Fig. 1F, p  = 0.033 and p = 0.004, respectively). However, there was considerable variability in the development of insulin resistance in response to the HFS diet for both male and female animals (Fig. 1F).

Fig. 1
figure 1

Markers of metabolic syndrome for both male and female animals at the end of the diet intervention protocol for rats in the chow and HFS-fed groups. All indicators of metabolic syndrome in female animals were found to be increased in the HFS group compared to the chow-fed controls. For male animals, all values were elevated except for total cholesterol and HDL cholesterol in the HFS animals compared to the controls. *Indicates statistical significance to the level of p < 0.05 compared to the chow-fed animals. Data are shown in box and whisker plots where the line across the box represents the median, and the whiskers extending from the bottom and top of the box represent the minimum value to the lower quartile and the upper quartile to the maximum value respectively.

MHC isoform distribution

Vastus lateralis

The two-way ANOVA analysis of the MHC IIx showed a statistically significant main effect of Sex (F(1,29) = 5.79, p = 0.023) and the interaction of Diet × Sex (F(1,29) = 5.54, p = 0.025.) indicating the percentage of MHC IIx may be influenced by both diet and sex (Fig. 2A). Tukey’s post-hoc evaluation revealed that the percentage of MHC IIx was significantly higher in the male chow-fed animal than in the female chow-fed animals (p = 0.011; Fig. 2A). The two-way ANOVA indicated that there was a Sex × Diet interaction effect for the percentage of MHC IIb (F(29,1) = 5.78, p = 0.023). However, follow up post-hoc tests did not show any significant differences between the groups.

Fig. 2
figure 2

Myosin heavy chain (MHC) isoform distribution in VL and soleus muscles of male and female rats fed both a chow-control diet and a HFS diet. *Indicates statistical significance to the level of p < 0.05 compared to the chow-fed male animals. #Indicates statistical significance to the level of p < 0.05 compared to the HFS-fed male animals. Error bars represent standard error of the mean.

Soleus

For the soleus muscle, the two-way ANOVA revealed a significant effect of Sex for the percentage of both MHC type I (F(1,28) = 38.42, p < 0.001) and MHC IIa (F(1,28) = 38.42, p < 0.001). Diet was not found to affect MHC isoform distribution in the soleus muscle. Follow-up Tukey’s post-hoc tests revealed that the percentage of MHC I was significantly higher in all male compared to the female animals and MHC IIa was significantly higher in all female animals compared to the male animals, irrespective of their diet (p = 0.011, Fig. 2B).

Vastus lateralis and soleus morphological changes

Vastus lateralis

The two-way ANOVA analysis showed a statistically significant main effect of Diet (F(1,41) = 33.932, p < 0.001) but not Sex (F(1,41) = 0.015, p = 0.902) on the intramuscular fat infiltration, represented by the % ORO staining, in the VL muscle (Fig. 3A). A Tukey’s post-hoc evaluation revealed that both the male and female HFS-fed animals had significantly higher intramuscular fat infiltration than their chow-fed counterparts (Fig. 3B) (p < 0.05). Similarly, the two-way ANOVA indicated a significant main effect of Diet (F(1,39) = 7.478, p = 0.009) but not Sex (F(1,39) = 2.285, p = 0.139) on the collagen content in the VL muscle, as represented by the % Picrosirius Red staining (Fig. 3D). However, post-hoc analysis did not show any meaningful differences between groups (Fig. 3E). In terms of a macrophage response to the HFS diet in the VL muscle, a two-way ANOVA determined that both Diet (F(1,39) = 6.855, p = 0.013) and Sex (F(1,39) = 22.392, p < 0.001) had a significant main effect on macrophage infiltration, as represented by macrophage count per image frame (0.15 mm2) (Fig. 3G). Tukey’s post-hoc analysis revealed that macrophage count was significantly higher in chow-fed males compared to chow-fed females (p = 0.043), and higher in HFS-males compared to HFS-females (p = 0.002, Fig. 3H). There was no difference in macrophage number between chow-fed females and HFS-fed females (p > 0.05, Fig. 3H).

Fig. 3
figure 3

Histological assessments of the VL and soleus muscle in control diet (chow) and HFS diet group rats at 24-week-old and following 12 weeks of HFS diet exposure for the experimental male and female Sprague–Dawley rats. (A) Oil Red O staining for intramuscular lipids in VL and soleus muscle sections in male and female animals taken at 100× magnification. (B, C) raw values for % Oil Red O staining per cross-sectional mid-belly muscle area represented by box and whisker plots. (B) intramuscular fat in the VL was increased after 12-weeks in the HFS diet compared to the chow control animals for both male and females (C) intramuscular fat in the soleus only increased in the female HFS animals following the 12-week metabolic challenge. (D) Picrosirius Red staining for collagen in both VL and soleus muscle sections in male and female animals taken at 100× magnification. (E, F) raw values for % Picrosirius Red staining per cross-sectional muscle area represented by box and whisker plots. No statistically significant difference was found in % Picrosirius Red staining between groups for both the VL and soleus muscle. (G) Immunohistochemistry staining for CD68 + cells in VL and soleus muscle sections, imaged at 200× magnification. Yellow/red CD68 + stain overlapping with blue DAPI stain, indicates the presence of a macrophage. (H, I) raw values for CD68 + macrophage count per image frame (area of 0.15mm2) for each animal, where 15–20 (VL) or 6–8 (soleus) images were randomly selected and evaluated for a given mid-belly muscle cross-section. (H) Shows that the macrophage counts in the VL muscle increased in chow-fed males compared to chow-fed females and increased the HFS-fed males compared to HFS-fed females (p < 0.05). No change in macrophage count was observed between chow-fed females and HFS-fed females (I) Shows that the macrophage response in the soleus muscle increased in the chow-fed males compared to the chow-fed females and increased in the HFS-females compared to the HFS-fed males to the HFS diet. No change was observed in macrophage count between the chow-fed male and HFS-fed male groups. Groups without a common letter are significantly different (p < 0.05).

Soleus

For the soleus muscle, the two-way ANOVA revealed a statistically significant main effect of Diet (F(1,40) = 19.52, p < 0.001), Sex (F(1,40) = 15.87, p < 0.001), and Diet × Sex interaction (F(1,40) = 12.82, p < 0.001) indicating that the level of the effect of the HFS diet on intramuscular fat infiltration in the soleus muscle is dependent on the sex of the animal (Fig. 3A,C). The Tukey’s post-hoc test showed that the percent ORO staining in the HFS-female group was higher compared to all other groups (p < 0.05). No main effect differences were found in the collagen content in the soleus muscle (p > 0.05) (Fig. 3D,F). According to the two-way ANOVA, macrophage count in the soleus muscle varied in response to Diet (F(1,35) = 18.88, p < 0.001) and a Diet × Sex interaction (F(1,35) = 30.971, p < 0.001) indicating that the amount of macrophage infiltration is dependent on both Diet and Sex (Fig. 3G,I). Tukey’s post-hoc analysis revealed that the macrophage count in the chow-males was significantly higher than the chow-fed females (p < 0.05, Fig. 3I). However, the HFS-fed females had a significantly higher macrophage count compared to all other groups (p < 0.05, Fig. 3I). There was no change in macrophage number observed between the chow-fed male and HFS-fed male groups (p > 0.05, Fig. 3I).

Relationship between body fat and intramuscular fat

For both male and female animals there was a significant positive relationship between body fat percentage and intramuscular fat in the VL muscle (Fig. 4A,C; males: R2 = 0.71, p < 0.001; females R2 = 0.61, p < 0.001). For the soleus muscle, there was a significant relationship between body fat percentage and intramuscular fat for both females and males (Fig. 4D,B; females R2 = 0.52, p < 0.001; males R2 = 0.20, p = 0.03).

Fig. 4
figure 4

Relationships between body fat percentage and intramuscular fat percentage in the VL and soleus muscle in male and female rats. The corresponding R2 values are shown on the figures. Medium to strong relationships were found between body fat percentage and intramuscular fat infiltration for male VL and female VL and soleus (p < 0.05).

Serum cytokines

Both male and female HFS diet animals had a statistically higher percentage of serum leptin levels compared to the chow-fed controls (Fig. 5A). For the male animals, all other serum cytokine levels were higher in the chow animals than in the HFS animals (Fig. 5B–D, p < 0.05). For the female animals, the HFS diet group exhibited similar serum cytokine levels to the controls for MCP-1 and IL-6 (Fig. 5B,D, p > 0.05). However, TNFα levels in female animals consuming the HFS diet were higher compared to the chow-fed animals (Fig. 5C, p = 0.05). Additional serum cytokine data can be found in Supplementary Table S3.

Fig. 5
figure 5

Serum cytokine and adipokine levels for 24-week-old male and female rats after 12 weeks on a HFS diet and the chow control diet. *indicates statistical significance to the level of p < 0.05 compared to the chow-fed control animals. Data are shown in box and whisker plots.

Quantitative proteomics analysis of muscle proteins

The proteomes of the male and female soleus and VL muscles were analyzed to determine whether their response to the HFS diet differed (Fig. 6A). In the proteomics analysis for the VL muscle (Supplementary Table S1), 95 unique proteins were identified that were upregulated in the female HFS-fed animals and 104 unique proteins that were upregulated in the male HFS-fed animals compared to the sex-matched chow control animals. Using the online meta-analysis tool, Metascape36, several enriched pathways between HFS-fed females and HFS-fed males were identified. Male animals fed the HFS diet showed enrichment in the “innate immune system” pathway in the VL muscle compared to chow control animals (Fig. 6B). In the proteomics analysis for the soleus muscle, 463 unique proteins were identified that were upregulated in the female HFS-fed animals and 90 unique proteins that were upregulated in the male HFS-fed animals compared to the sex-matched chow control animals (Supplementary Table S2). Using Metascape, several enriched pathways were identified including innate immune system and cytokine signalling in the immune system of the soleus muscle (Fig. 6C). The quantitative proteomics analysis indicates key biological differences between male and female VL and soleus muscles.

Fig. 6
figure 6

Quantitative proteomics of the VL and soleus muscle of male and female Sprague–Dawley rats that have been exposed to the HFS diet for 12 weeks. (A) Proteomics workflow and analysis (N = 3, HFS males and N = 3, HFS females). Metascape software was used to identify pathway enrichment. A list of all detected proteins is provided in supplementary tables S1 and S2. (B) Metascape analysis of different pathways between HFS diet animals and Chow-fed control animals for the VL muscle. (C) Metascape analysis of different pathways between HFS diet and Chow-fed control animals for the soleus muscle. Blue writing indicates upregulated pathways.

Discussion

The link between obesity and the decline in muscle integrity, classified by the increase in intramuscular lipid content and fibrosis, has been previously well established19,20,21,23,37,38,39. However, previous studies failed to consider female models of diet-induced obesity. The aim of the current study was to determine if previous findings from the model of diet-induced obesity in male rats can be generalized to both sexes, or whether sex affects the response to the HFS diet, as indicated by markers of metabolic syndrome and changes to muscle integrity. An important finding from this study was that female Sprague–Dawley rats became obese and showed symptoms consistent with metabolic syndrome after 12-weeks on the HFS diet. Similar to the male animals, female animals exposed to the HFS diet displayed a significant increase in body fat percentage, elevated blood triglycerides, and increased insulin resistance (Fig. 1). However, female animals on the HFS diet also experienced an increase in HDL cholesterol levels where males did not. It has been well established in previous literature that estrogen is associated with an increase in HDL cholesterol40. Exposure to the HFS diet also resulted in various effects of both diet and sex on intramuscular fat infiltration, fibrosis, and macrophage invasion in the VL and soleus muscles. The most notable finding was the significant increase in intramuscular fat in the soleus muscle for the female animals on the HFS diet, compared to all the other experimental groups.

Overall, the MHC isoform distribution results did not indicate that there were MHC isoform distribution changes in response to the HFS diet challenge. Although the two-way ANOVA analysis revealed an interaction effect between Diet and Sex in the VL muscle, post-hoc analysis only identified the male chow-fed animals to be different from the female chow-fed animals (Fig. 2A). It has been shown that male mice exposed to a HFS diet had a fibre type change from type I to type I/IIa hybrid. However, this result was obtained using histology and SDS-PAGE on single muscle fibres and after mice had been fed an obesogenic diet for approximately a year41. It is possible that co-expression of MHC isoforms within fibres was missed in our study because SDS-PAGE was performed on whole muscle samples. However, MHC transitions are slow processes and 12-weeks on the HFS diet was potentially too short a time to observe a change in MHC expression. The soleus muscle showed distinct differences in composition between the male and female animals in both diet groups (Fig. 2B). It has been reported previously that the soleus muscle in adult Sprague Dawley rats expresses nearly 100% MHC I22. We confirmed this finding as the male animals in the current study had 100% MHC I in the soleus muscle. However, the female animals showed a mixed distribution of MHC I and MHC IIa isoforms. MHC IIa is the most oxidative of the MHC II (fast) isoforms and it makes up less than 10% of the soleus muscle MHC content in female Sprague Dawley rats, indicating that the soleus in the female rats still has a high oxidative capacity.

Leptin was found to increase in both HFS diet male and female animals compared to the chow-fed controls (Fig. 5A). Leptin is a peptide hormone predominantly produced by adipocytes that influences appetite and satiety to regulate energy homeostasis42. Consequently, leptin levels were associated with increases in body fat percentage in both male and female animals following the HFS diet intervention period. The proinflammatory cytokines MCP-1 and IL-6 appeared to be restored to homeostatic levels for the HFS-females following the 12-week metabolic challenge (Fig. 5B,D). After 12-weeks on the HFS diet, obesity may have become a chronic condition, and a rise in inflammatory cytokines in response to the HFS diet may only occur as an acute response to the obesity-inducing diet. This has been suggested by a previous study where MCP-1 and TNFα showed an initial increase to a HFS diet but then returned to levels similar to those of the chow control group after 4 weeks20. In the present case however, TNFα continued to be elevated after 12 weeks on the HFS diet (Fig. 5C). This sustained increase is likely because the HFS diet has been shown to compromise the intestinal barrier thereby allowing lipopolysaccharides (LPS) to translocate into systemic circulation and induce metabolic endotoxemia43,44. LPS is a metabolic by-product of Gram-negative bacteria and is associated with increases in TNFα45.

It is important to note that there was great variability in the response of the animals to the HFS diet. Sprague–Dawley rats are an outbred strain which means they have been bred to maximize genetic diversity and heterozygosity within a population. Although this likely introduced some of the variability in our results, this scenario is more representative of the human population than an inbred rat strain would be, as humans also have heterogeneous responses to diet46. However, small university-housed colonies, such as the one used in the current study, can experience genetic drifts, and show genetic divergence from other small colonies47. This may offer a possible explanation of the unusual serum cytokine results in the male animals where the chow-fed control group displayed significantly higher levels than the HFS group for almost all inflammatory markers measured (Fig. 5). This finding is in stark contrast to results reported from a similar study that observed either an increase or no change in serum inflammatory cytokines in male Sprague–Dawley rats who underwent a 4-week metabolic challenge20.

The morphological changes of the VL muscle observed following a 12-week metabolic challenge are similar to those observed in previous studies in male Sprague–Dawley rats that found that the HFS diet was associated with an increase in fat infiltration in the VL muscle21,23. It has been shown previously that intramuscular fat in the VL increases proportionally with the increase in body fat percentage21. Here, we confirm this result for the female VL muscle in response to the HFS diet, as there was a significant correlation between body fat and VL intramuscular fat (Fig. 4C). Surprisingly, collagen content was not significantly changed in the VL muscle in response to the HFS diet for either sex. However, this may be attributed to the length of the metabolic challenge, as increases in collagen deposition/fibrosis have been observed following 28 weeks on a HFS diet21.

Fat is an active endocrine organ that can release inflammatory cytokines that attract proinflammatory macrophages14,48. Macrophage infiltration into the VL and soleus muscle in the HFS animals confirms previous findings that it seems to follow intramuscular fat infiltration, giving an indication of the local inflammatory environment. The proteomics analysis confirms the inflammation, as the innate immune system is triggered when immune cells, such as macrophages, detect tissue damage (Fig. 6). By comparing the male and female HFS diet groups of animals, we were able to determine which sex had a more pronounced inflammatory response to intramuscular fat infiltration and what specific proteins are responsible for these differences (the complete list is shown in Supplementary Tables S1-S2). We observed that males had a larger inflammatory response in the VL muscle, and females had a stronger inflammatory response to the infiltrated fat in the soleus muscle, similar to what was observed for the macrophage infiltration (Fig. 3H and I). It has been observed in certain species that females in general have a more robust and persistent immune response than males49,50. Specifically, it has been shown that testosterone decreases certain immune functions such as cytokine production51,52. The chronic inflammation involved in the immune response can have various deleterious affects on healthy tissue, putting females at a higher risk for inflammatory diseases and autoimmune conditions53.

One of the most exciting findings of the present study is that there was a large fat infiltration in the soleus muscle of the female HFS diet group animals (Fig. 3C) compared to the males. This result detected in females is contrary to findings recorded in previous studies that suggest that the oxidative capacity of the soleus may be protective against intramuscular fat infiltration and fibrosis in male Sprague–Dawley rats23. A possible explanation may be offered by how males and females store fat. It is thought that estrogen mediates fat deposition preferentially away from visceral fat deposits and towards subcutaneous deposits in female rats54. Directing fat away from visceral deposits may make it more likely to be deposited ectopically in muscle tissues, as our data suggest that intramuscular fat infiltration in female soleus muscle also increases with body fat percentage. From an evolutionary perspective, because visceral fat contributes more to metabolic dysregulation55, it may have been useful for females to have other mechanisms for storing fat than in the visceral region, thereby reducing the risk of disease and allowing them to care for their offspring and nurture fetuses while pregnant.

A limitation to the current study is that only mid-belly slices of both the VL and soleus were analysed. This may not be representative of the fat distribution throughout the muscle. It has been found that there are regional differences in intramuscular fat and that mid-belly fat may not necessarily be representative of whole muscle fat content56,57,58. Another limitation stems from the serum cytokine profile collected in this study. Serum cytokines are systemic markers of inflammation while local inflammation may differ from the systemic state and may give more insight about muscle adaptation in the presence of a metabolic challenge. Results from this study call for further investigation into different muscle types and how they may be affected by obesity, as we have shown that not all muscles respond the same in the presence of a HFS diet exposure. Future studies should be aimed at quantifying the effects of a HFS diet on female sex hormones to determine the potential role of sex hormones in the regulation of fat storage.

Like humans, Sprague–Dawley rats are heterogeneous in their response to an obesity-inducing diet. In this study we illustrate similarities and differences between male and female Sprague–Dawley rats in their response to a 12-week metabolic challenge. In summary, it was found that female Sprague–Dawley rats are susceptible to obesity and show similar symptoms of metabolic disease as male rats when exposed to a HFS diet. The female animals also displayed a similar response to the diet with respect to intramuscular fat infiltration in the VL muscle. However, females receiving the HFS diet were more susceptible to intramuscular fat infiltration in the high oxidative capacity soleus muscle than males. Including females in this research may allow for a better understanding of the mechanistic relationship between diets, obesity, metabolic syndrome, and the effects of these factors on adaptation and degeneration of musculoskeletal tissues, by studying the differences in metabolic response and tissue adaptation/degeneration between the sexes.