Constance de Meeûs d'Argenteuil1†
Berit Boshuizen1,2†
Carmen Vidal Moreno de Vega1
Luc Leybaert3
Lorie de Maré1
Klara Goethals4
Ward De Spiegelaere5
Maarten Oosterlinck6
Cathérine Delesalle1*- 1Department of Translational Physiology, Infectiology and Public Health, Research Group of Comparative Physiology, Faculty of Veterinary Medicine, Research Group of Comparative Physiology, Ghent University, Merelbeke, Belgium
- 2Wolvega Equine Hospital, Oldeholtpade, Netherlands
- 3Department of Basic and Applied Medical Sciences, Faculty of Medicine and Health Sciences, Ghent University, Ghent, Belgium
- 4Department of Veterinary and Biosciences, Faculty of Veterinary Medicine, Research Group Biometrics, Ghent University, Merelbeke, Belgium
- 5Department of Morphology, Imaging, Orthopedics, Rehabilitation and Nutrition, Faculty of Veterinary Medicine, Ghent University, Merelbeke, Belgium
- 6Department of Large Animal Surgery, Anaesthesia and Orthopaedics, Faculty of Veterinary Medicine, Ghent University, Merelbeke, Belgium
Training-induced follow-up of multiple muscle plasticity parameters in postural stability vs. locomotion muscles provides an integrative physiological view on shifts in the muscular metabolic machinery. It can be expected that not all muscle plasticity parameters show the same expression time profile across muscles. This knowledge is important to underpin results of metabolomic studies. Twelve non-competing Standardbred mares were subjected to standardized harness training. Muscle biopsies were taken on a non-training day before and after 8 weeks. Shifts in muscle fiber type composition and muscle fiber cross-sectional area (CSA) were compared in the m. pectoralis, the m. vastus lateralis, and the m. semitendinosus. In the m. vastus lateralis, which showed most pronounced training-induced plasticity, two additional muscle plasticity parameters (capillarization and mitochondrial density) were assessed. In the m. semitendinosus, additionally the mean minimum Feret's diameter was assessed. There was a significant difference in baseline profiles. The m. semitendinosus contained less type I and more type IIX fibers compatible with the most pronounced anaerobic profile. Though no baseline fiber type-specific and overall mean CSA differences could be detected, there was a clear post-training decrease in fiber type specific CSA, most pronounced for the m. vastus lateralis, and this was accompanied by a clear increase in capillary supply. No shifts in mitochondrial density were detected. The m. semitendinosus showed a decrease in fiber type specific CSA of type IIAX fibers and a decrease of type I fiber Feret's diameter as well as mean minimum Feret's diameter. The training-induced increased capillary supply in conjunction with a significant decrease in muscle fiber CSA suggests that the muscular machinery models itself toward an optimal smaller individual muscle fiber structure to receive and process fuels that can be swiftly delivered by the circulatory system. These results are interesting in view of the recently identified important fuel candidates such as branched-chain amino acids, aromatic amino acids, and gut microbiome-related xenobiotics, which need a rapid gut–muscle gateway to reach these fibers and are less challenging for the mitochondrial system. More research is needed with that respect. Results also show important differences between muscle groups with respect to baseline and training-specific modulation.
Introduction
Skeletal muscles are remarkable in adapting to certain stimuli, such as exercise (1). The muscular compartment represents ±40 up to 57% of total body mass, depending on animal species and breed, which implies that muscular plasticity has a crucial impact on functionality, metabolism and performance capacity of the whole-body system (1–3). It also implies that shifts that take place at the level of certain specific features of the muscular compartment can be viewed as a marker for occurrence of metabolic shifts in answer to training. Several of these markers have been studied in the past, either alone or combined, such as shifts in muscle fiber type composition, muscle fiber cross-sectional area (CSA), mitochondrial density and capillarization (4–7).
Muscle fibers are generally classified as either type I fibers (slow oxidative), type IIA fibers (fast oxidative-glycolytic) or type IIX fibers (fast glycolytic), co-existing with a wide array of transitional hybrid fiber types, such as type I/IIA and type IIA/IIX (8). In horses, the type IIA fiber is the predominant type, followed by the type IIX fiber and the type I fiber (9–12). In general, training can shift muscle fibers from type IIX to the transitional hybrid fiber type IIAX, subsequently to IIA and finally to type I and vice versa (13, 14). Besides shifts in muscle fiber type composition, fiber type specific CSA and mean CSA shifts across all fiber types occur in answer to training. Both fiber type composition and fiber size provide important information on main performance capacity, such as force generation, muscle fiber conduction velocity, and endurance capacity in different species, such as horses and humans (15, 16). It has been shown, for example, that muscle fiber type composition and CSA are directly related to the performance capacity of human endurance athletes (12, 17). Muscle fiber CSA positively correlates with glycogen storage capacity (18), and thus an anaerobic fast-twitch muscle fiber profile (19, 20). Aerobic exercise, on its turn, induces muscle hyperplasia and decreases the muscle CSA, resulting in a smaller distance for oxygen diffusion and swift fuel transport into muscle fibers (21, 22).
The change in muscular capillary supply in answer to training is deemed to be an important plasticity parameter (19). The baseline muscular capillarization level depends on different parameters such as muscle fiber type, muscle fiber CSA, oxidative capacity, and metabolic environment (23, 24). Physiologically, it is interesting to notice that in both humans (25–28) and horses (4, 24, 29–31), baseline muscular capillarization level increases in answer to training, no matter what type of training is imposed (either aerobic endurance vs. anaerobic resistance/sprint type of training).
Finally, also mitochondrial plasticity occurs following training (32, 33). However, unlike capillarization, the observed mitochondrial plasticity seems to depend on the type of imposed exercise (34). On top of that, both an acute (following a few high-intensity exercise sessions) (35) and a long-term (following several weeks to months of training) (36–38) mitochondrial plasticity can occur. Mitochondria are pivotal for oxidative metabolism and adapt to different types of exercise by both morphological and functional modulations (36, 39, 40). Morphological adaptations comprise changes in mitochondrial density, and at the ultrastructural level, changes in the amount and length of mitochondrial cristae and density of mitochondrial DNA. Functional adaptations encompass changes in mitochondrial enzymatic activity [e.g., citrate synthase (CS), succinate dehydrogenase activity, monocarboxylate transporter 1 (MCT1), basigin, 3-hydroxyacyl-CoA dehydrogenase (HAD), hexokinase (HK) and lactate dehydrogenase (LDH)] or changes in mitochondrial respiration (41–43). Endurance exercise, for example, increases mitochondrial density, whereas strength training decreases mitochondrial density by inducing muscle fiber hypertrophy and mitochondrial fragmentation (44, 45). In horses, only a limited number of studies on the effects of different types of exercise on muscle mitochondrial density are available (4, 23, 46–48).
Skeletal muscles contain two distinct mitochondrial subpopulations: intermyofibrillar mitochondria (IMF), distributed between myofibrils and subsarcolemmal (SS) mitochondria situated underneath the sarcolemmal membrane (36, 43). Both specific populations are reported to react differently in answer to training (42, 43, 49). Oxidative or slow-twitch fibers have a high mitochondrial density, predominantly in the subsarcolemmal region, whereas fast-twitch fibers have a low mitochondrial density that appears more fragmented throughout the muscle fiber (50). Mitochondrial populations form a dynamic network within the muscle fibers, allowing them to interconnect via fusion and fission (42, 43, 49). Inter-mitochondrial connections expand in answer to aerobic exercise, due to increased fusion between mitochondria of either subtype, i.e., SS or IMF (51).
Longitudinal follow-up of all four muscle plasticity parameters (shifts in muscle fiber type composition, muscle fiber CSA, capillarization, and mitochondrial density) provides a more integrative physiological view on what changes to expect at the level of the muscular metabolic machinery, when compared to focusing on a subset of these plasticity parameters. It can be expected that not all four muscle plasticity parameters show the same expression time profile and thus differences can be expected between these plasticity parameters when it comes to early identification of training adaptation (4, 6).
A study by Tyler et al. reports the longitudinal follow-up of all aforementioned muscle plasticity parameters (fiber type composition, CSA, capillarization, and mitochondrial density) in answer to respectively 7 weeks of endurance training and 8 weeks of high-intensity exercise in Standardbreds, focusing on one muscle: the m. gluteus medius (4). In that study, both capillarization and mitochondrial density showed significant changes in answer to both types of training, whereas no significant changes could be found with respect to shifts in fiber type composition and CSA. Therefore, not all muscle plasticity parameters seem equally discriminative or sensitive when it comes to follow-up of training adaptation. Another study reports the absence of change in mitochondrial density in the m. gluteus medius of Quarter Horses subjected to respectively 9 (46) and 14 weeks (52) of submaximal exercise. However, these researchers did find training-induced increased mitochondrial respiratory activity, expressed by upregulation of OXPHOS complex II activity, and electron transport chain capacity, underlining the fact that mitochondrial functional shifts represent an additional set of possible muscular plasticity parameters (46, 53, 54).
Many training studies have been combining two or three out of the four aforementioned muscle plasticity parameters (5–7), most of them, focusing on the m. gluteus medius, known as an archetype locomotion-function type of muscle group (55, 56).
However, to our knowledge, at this point, no training studies are available involving longitudinal follow-up of all four muscle plasticity parameters in multiple muscles with differing main mechanical function (muscles mainly involved in postural stability vs. muscles involved in locomotion). It can be expected that each type of training modulates a different core set of muscles (55, 56). By involving both “main posture vs. main locomotion directed muscle groups” at the same time, a much more integrative and subtle physiological insight is obtained into how the muscular compartment adapts in response to a certain training type. Depending on the main mechanical function of a certain muscle group, different metabolic challenges are imposed by the applied training type. Obtaining such detailed view on muscular adaptation is crucial for further optimization of existing training protocols and understanding of muscular pathologies.
The aim of the current study was to map out and compare the evolution of four important muscle plasticity parameters in a core set of muscles in answer to training, to obtain a better view on flexibility of these parameters and their expression profile in muscles mainly involved in postural stability vs. muscles involved in locomotion. For this purpose, 12 Standardbred horses were subjected to a harness training protocol of 8 weeks duration.
Materials and Methods
Study Setup
Twelve untrained, non-competing Standardbred mares aged between 3 and 5 years were subjected to a standardized training trial (N = 12). Horses were comparable in body condition score (5.5 ± 0.70) and body weight (446.5 ± 20.94 kg) and housed at the same training facility in individual boxes with straw bedding. Turn-out was provided on paddocks 2 h a day. Horses were fed the same concentrate feed twice a day, at 8 a.m. and 8 p.m. and had ad libitum access to tab water and roughage. Monitoring for vital signs (pulse, temperature, respiratory rate, capillary refilling time, appetite, and stool production) was performed twice a day and horses were allowed 2 weeks of acclimation prior to start of the training trial. All horses remained in good health throughout the study. No adverse events occurred.
Training Regime
All horses (N = 12) were harness trained during eight consecutive weeks by the same experienced driver 4 days a week on the same oval-shaped sand race track, between 9 and 12 a.m., during the cooler European months. Horses were equipped with a GPS tracker and heart rate monitoring system (Polar® Equine H7, Polar Electro Oy, Finland). Each training session consisted of a warming-up period of 10 min jogging (±20 km/h), followed by either 30 min of aerobic training (3 days a week, ± 25 km/h, mean HR < 150 BPM) or interval training (day 4) consisting of 3 × 3 min at high speed (±35 km/h). The track was watered and raked on each occasion, to maximally assure equal study conditions throughout the trial. Muscle biopsies were taken on a resting day before and after the 8-week training period, at the level of the m. pectoralis (expected postural predominant), the m. semitendinosus (expected pure locomotor), and the m. vastus lateralis of the quadriceps femoris (expected locomotor predominant). This study was approved by the Animal Ethics Committee of the Ghent University (EC 2016/40).
Muscle Biopsies
Briefly, the horses were sedated with detomidine (10 μg/kg bwt) (Detogesic, Vetcare, Finland) and butorphanol (20 μg/kg bwt) (Butomidor, Richter Pharma AG, Wels, Austria). The area was clipped, shaved, and subsequently disinfected. Local anesthetic ointment was applied (Emla 5%, Astra-Zeneca, Rueil-Malmaison, France). After 10 min, local anesthetic (Lidocaine Hydrochloride, Braun, Germany) was injected subcutaneously and a small stab incision was made with a surgical blade number 11. Thereafter, a 14G Bergström needle was inserted into the muscle, until a depth of 4 cm was reached on each occasion. Samples were placed in Tissue-Tek O.C.T compound (Sakura Finetek, Torrance, California) and were immediately snap-frozen in liquid nitrogen-cooled isopentane (Sigma-Aldrich, Dorset, UK) and stored at −80°C until further processing.
Muscle Fiber Typing
Cryosections of 8 μm were created from the Tissue-Tek embedded samples and were collected onto coated glass slides (Thermo Scientific SuperFrost Plus Adhesion slides, Fisher scientific, Belgium) and stored at −20°C until further processing. Briefly, the sections were air-dried and then blocked for 120 min in 1% bovine serum albumin (BSA) in phosphate buffered saline (PBS) solution at room temperature. The slides were than shortly washed with permeabilization solution (BSA and 0.2% triton in PBS solution). Thereafter, the slides were incubated overnight at −4°C with primary antibodies for the different myosin heavy chains (57) dissolved in 0.5% BSA in PBS, for type I, type IIA, type IIX, and sarcolemma [BA-D5 (type I), dilution 1:400; DSHB, RRID:AB_2235587; SC-71 (type IIA), dilution 1:100; DSHB, RRID:AB_2147165; 6H1 (type IIX), dilution 1:50; DSHB, RRID:AB_1157897 and laminin (sarcolemma), dilution 1:500; Thermo Fisher Scientific Cat: PA1-36119, RRID:AB_2133620]. No controls were applied. After rinsing the slides five consecutive times during 5 min in PBS, and washing them shortly with permeabilization solution, they were incubated with the secondary antibodies dissolved in 0.5% BSA in PBS for 1 h at room temperature for type I, IIA, IIX (dilution 1:300), and sarcolemma (dilution 1:600) [Alexa fluor 488 goat anti-mouse IgG2b (type I), Thermo Fisher Scientific Cat: A-21141, RRID:AB_2535778; Alexa fluor 350 goat anti-mouse IgG1 (type IIA), Thermo Fisher Scientific Cat: A21120, RRID:AB_2535763; Alexa fluor 594 goat anti-mouse IgM (type IIX), Thermo Fisher Scientific Cat: A-21044, RRID:AB_2535713; Alexa fluor 568 goat anti-rabbit IgG (sarcolemma), Thermo Fisher Scientific Cat: A-11011, RRID:AB_143157]. After rinsing the slides five consecutive times during 5 min in PBS, they were air dried and fluorescent mounting medium (Dako, Agilent, S3023) was applied on the slides with a cover glass.
For the capillaries and the mitochondria, the slides were incubated overnight at −4°C with the primary antibodies targeting the Von Willebrand factor (capillaries) (polyclonal rabbit anti-human Von Willebrand factor, dilution 1:1600; Agilent Cat# A0082, RRID:AB_2315602) and COX IV (mitochondria) (anti-COX IV antibody, dilution 1:100; Abcam Cat# ab14744, RRID:AB_301443) diluted in 0.5% BSA in PBS. After rinsing the slides and washing them with permeabilization solution, the sections were incubated for 1 h at room temperature with the secondary antibodies diluted in 0.5% BSA in PBS [Alexa fluor 568 goat anti-rabbit IgG (capillaries), dilution 1:600; Thermo Fisher Scientific Cat# A-11011, RRID:AB_143157 and Alexa fluor 647 goat anti-mouse IgG2a (mitochondria), dilution 1:500; Thermo Fisher Scientific Cat# A-21241, RRID:AB_2535810]. The sections were visualized with a Zeiss Palm Micro Beam fluorescence microscope and pictures were taken with the Zen Blue Pro Software (Zeiss). Fibers were classified as either type I (green), type IIA (blue), type IIX (red), or as hybrid type when staining for more than one myosin heavy chain was present. Total fiber count, fiber type percentages, and CSAs were determined as well as capillary count and mitochondrial intensity using Image Pro v.10 analyzer software (Media Cybernetics, Inc., Rockville, USA). All sections were evaluated in duplicate, with a minimum of 250 fibers per section. The mean CSA of each fiber type was measured, called from now on fiber type specific CSA, as well as the overall mean CSA across all muscle fiber types. The mean CSA of all fibers was calculated for all three muscles in two different ways: one calculating the average of all the CSA obtained for each fiber type that we will refer to as “mean mathematical CSA”, and another method that relies on a count per frame approach to which we will refer as “mean frame CSA”. For the m. semitendinosus, we also calculated for each fiber type and for all fibers the mean minimum Feret's diameter, which is the closest possible distance between the two parallel tangents of a muscle fiber. Capillaries surrounding type I and type IIA fibers were counted (CAF), and capillary density (CD) (number of capillaries/mm2) and capillaries per fiber ratio (C/F) (depicting the number of fibers in contact with each capillary) were calculated. CAF and C/F were counted over the whole section, and CD was done by counting and summing four fields of each 250,000 μm2, chosen at random, but covering the whole section. For assessment of mitochondrial density, each muscle fiber image was divided into two segments: one from the sarcolemma going 15 μm deeper into the fiber (subsarcolemmal mitochondria) and one from 15 μm deep, to the center of the fiber (intermyofibrillar mitochondria). Mitochondrial density was expressed as the intensity of fluorescence in type I and type IIA fibers and intensity in the SS and the IMF region and the fold change was calculated by dividing mitochondrial density of type I by mitochondrial density of type IIA. To assess mitochondrial density, a minimum of 75 fibers for each fiber type were chosen at random in different fields, but covering the whole section. To look at the effect of training, a ratio between SS mitochondria and IMF mitochondria was calculated as well as the difference in mitochondrial density between type I and type IIA fibers before and after training.
Statistical Analysis
Different fiber types were counted and classified per type and relative percentages of each type were calculated for each of the muscles. The mean CSA using both methods as well as mean fiber CSA of all types of fibers separately were determined in the three muscles (in μm2).
Statistical analysis was performed in R (104). Normality of the data was graphically checked using histograms and Shapiro–Wilk's test. Because normality of the data could not be assumed, data were analyzed using nonparametric tests. The results are given as median (minimum–maximum). Significance was set at p < 0.05.
The effect of muscle type (m. pectoralis vs. m. vastus lateralis vs. m. semitendinosus) on the percentages of fiber types, fiber type specific CSA, mean frame CSA, and mean mathematical CSA was analyzed using a Friedman test. If the effect was significant, pairwise comparisons were tested using a Wilcoxon signed rank test at a significance level of 0.017 to account for multiple comparisons (Bonferroni correction for multiple pairwise comparisons: α = 0.05/3).
The effect of training on muscle fiber type composition, fiber type specific CSA, mean frame CSA and mean mathematical CSA per studied muscle, and the fiber type specific minimum Feret's diameter, mean minimum Feret's diameter, capillary, supply and mitochondrial density was analyzed using a Wilcoxon signed rank test.
Results
Comparison of Baseline Muscle Fiber Type Composition, Fiber Type Specific CSA, Mean Mathematical CSA, and Mean Frame CSA in the m. pectoralis, m. vastus lateralis, and m. semitendinosus
In the current study, apart from the “pure” fiber types, only hybrid type IIA/IIX (IIAX) fibers were included in the analysis, since the hybrid type I/IIA was only sporadically found (<0.5%).
Only the baseline profiles of fiber types I (p = 0.0057) and IIX (p = 0.0001) showed significant differences between the studied muscles (Table 1).
Table 1. Baseline muscle fiber type composition of the m. pectoralis, the m. vastus lateralis, and the m. semitendinosus.
There was a significantly higher percentage of fiber type I content for both the m. pectoralis (p = 0.0017) and the m. vastus lateralis (p = 0.0081) compared to the m. semitendinosus. A significantly higher percentage of type IIX fibers was present in the m. semitendinosus compared to the m. pectoralis (p = 0.0005) and the m. vastus lateralis (p = 0.0005). There were no significant differences in percentage of type IIA and IIAX fibers between the three different muscles (Table 1).
No significant difference in fiber type specific CSA of type I, type IIA, and type IIAX could be found between the three muscles before training. Additionally, there was no significant difference observed in baseline mean frame and mean mathematical CSA between the three muscles.
Training-Induced Shifts in Muscle Fiber Type Composition, Fiber Type Specific CSA, Mean Frame CSA, and Mean Mathematical CSA in the m. pectoralis, the m. vastus lateralis, and the m. semitendinosus
No significant shifts in fiber type composition could be detected in the m. pectoralis and the m. vastus lateralis for either type I, type IIA, type IIAX, or type IIX; however, the m. semitendinosus showed a significant increase in type I fibers in answer to training [from 13.6 (0.1–30.4) to 26.9% (6.9–34.2)] (p = 0.0034).
Fiber type specific CSA showed evolution in answer to training, but not in all studied muscles. No change in fiber type specific CSA of type I, IIA, IIAX, and IIX could be detected in the m. pectoralis in answer to training. The m. vastus lateralis, however, showed much more plasticity: there was a significant decrease in fiber specific CSA of all fiber types in answer to 8 weeks of harness training (type I: p = 0.0046; type IIA: p = 0.0061; type IIAX: p = 0.0134; and type IIX: p = 0.0081). Also, the m. semitendinosus showed plasticity, though to a much lesser extent. In that muscle, there was only a significant decrease of the type IIAX fiber type specific CSA (p = 0.0342) (Table 2).
Table 2. Comparison of fiber type specific cross-sectional area (CSA) of the different muscle fiber types in the m. pectoralis, the m. vastus lateralis, and m. semitendinosus before training (untrained) and after 8 weeks of harness training (trained).
When comparing mean mathematical and mean frame CSA before and after training per studied muscle group, the m. vastus lateralis showed significant training-induced shifts. There was a significant decrease in both the mean frame CSA [from 5,740 (3,487–8,568) μm2 to 3,818 (3,281–4,640) μm2] (p = 0.0007) and the mean mathematical CSA [from 5,351 (3,562–7,929) μm2 to 3,985 (3,025–4,922) μm2] (p = 0.0046) in answer to training. No significant changes were recorded for the m. pectoralis and the m. semitendinosus (Figure 1).
Figure 1. Evolution of mean mathematical CSA (A) and mean frame CSA (B) in the m. pectoralis, m. vastus lateralis, and m. semitendinosus before (white plots) and after (gray plots) 8 weeks of harness training. *Significantly different from untrained, p < 0.01.
Training-Induced Shifts in the Fiber Type Specific Minimum Feret's Diameter and Mean Minimum Feret's Diameter in the m. semitendinosus
There was a significant effect of harness training on the fiber type specific minimum Feret's diameter in the m. semitendinosus of fiber type I (p = 0.0214), which shows a decrease in diameter after 8 weeks of training. Fiber type IIA, IIAX, and IIX specific minimum Feret's diameter did not show any significant change in answer to training (Table 3).
Table 3. Comparison of fiber type specific and mean minimum Feret's diameter in the m. semitendinosus before training (untrained) and after 8 weeks of harness training (trained).
Across all fiber types, there was a significant training-induced decrease in the mean minimum Feret's diameter in the m. semitendinosus after 8 weeks of training (p = 0.0398).
Comparison of Mean Mathematical CSA and Mean Frame CSA Across All Three Studied Muscles After 8 Weeks of Harness Training
The mean frame CSA was significantly different between the three muscle groups after 8 weeks of harness training (p = 0.0029). The mean frame CSA was smaller in the m. vastus lateralis [3,818 (3,281–4,620) μm2] when comparing to the m. pectoralis (p = 0.0105) [4,783 (3,209–6,238) μm2] and to the m. semitendinosus (p = 0.0017) [4,680 (3,642–8,794) μm2].
When looking at mean mathematical CSA, no significant difference was observed between the three muscle groups.
Training-Induced Shifts in Capillary Supply in the m. vastus lateralis
Since no shifts were seen in muscle fiber type composition, or in fiber specific and mean fiber CSA of the m. pectoralis in answer to 8 weeks of harness training, and only few changes were observed in the m. semitendinosus, it was decided to focus on the m. vastus lateralis to further map out changes in capillary supply and mitochondrial density.
Capillary density (1.58-fold change; p = 0.0015), the number of capillaries surrounding type I (1.17-fold change; p = 0.0068) and type IIA fibers (1.21-fold change; p < 0.001), and the capillary-to-fiber ratio (C/F) (1.24-fold change; p < 0.001) were all significantly increased in answer to 8 weeks of harness training (Table 4).
Table 4. Comparison of capillary density, number of capillaries surrounding type I and type IIA fibers, and capillary-to-fiber ratio before training (untrained) and after 8 weeks of harness training (trained).
Training-Induced Shifts in Mitochondrial Density in the m. vastus lateralis
Significant differences in baseline fiber type specific mitochondrial subpopulations could be detected in the m. vastus lateralis (p < 0.001). However, no significant training-induced shifts could be detected with respect to general and subpopulation specific mitochondrial density in the m. vastus lateralis.
No significant changes in mitochondrial density could be detected in answer to training (Table 5). The ratio between SS and IMF mitochondria in type I fibers and type IIA fibers did not change in answer to harness training. When looking at the effect of harness training on mitochondrial density between type I and type IIA fibers, no significant effect could be detected for either SS or IMF mitochondria (Table 5).
Table 5. Ratio between subsarcolemmal and intermyofibrillar mitochondrial density in type I and type IIA fibers and ratio between type I and type IIA fibers of subsarcolemmal and intermyofibrillar mitochondria before training (untrained) and after 8 weeks of harness training (trained).
Before training, both SS (p < 0.001) and IMF (p < 0.001) mitochondria were significantly more represented in type I fibers when compared to type IIA fibers (1.44- and 1.75-fold, respectively), and this was still the case after training (1.63- and 1.98-fold, respectively) (p < 0.001). Note that type I fibers contain twice as many IMF mitochondria as type IIA fibers. Furthermore, the mitochondrial density of the SS region was higher than the density at the IMF region, and this was the case in both fiber types before (1.52-fold for type I and 1.85-fold for type IIA; both p < 0.001) and after 8 weeks of harness training [1.49-fold for type I (p < 0.001) and 1.80-fold for type IIA (p < 0.001)] (Table 6).
Table 6. Comparison of subsarcolemmal and intermyofibrillar mitochondrial density, expressed as intensity of fluorescence, between type I and type IIA fibers before training (untrained) and after 8 weeks of harness training (trained).
Discussion
To the best of our knowledge, this is the first training study to monitor four different muscle plasticity parameters (muscle fiber type composition, fiber type specific and mean fiber CSA, capillary supply, and mitochondrial density) in three muscles (posture vs. locomotor predominant). This study approach has enabled a better understanding of which muscle plasticity parameters change most rapidly in response to harness training. By combining locomotion dominant vs. posture dominant muscles, the observed changes can also be more reliably translated to the metabolic metamorphosis undergone by the muscle group involved. Results of the study show that of all four studied muscle plasticity parameters, the mean CSA, and the capillary supply show pronounced shifts in answer to harness training and thus can function as early indicators of training adaptation. With that respect, the m. vastus lateralis showed most plasticity, showing a decrease in CSA of all studied muscle fiber types, followed by the m. semitendinosus, showing only a decrease in CSA of the type IIAX fibers. Both, the m. vastus lateralis and m. semitendinosus are suitable for repeated harvesting of muscle biopsies without the occurrence of complications.
When looking into scientific literature, a large body of research has primarily focused on the m. gluteus medius (4, 7, 58, 59), and only a few studies have included more than one muscle (60–63); in some cases, even a combination of postural and locomotor muscles have been included (60, 64), but never focusing on all four muscle plasticity parameters at the same time.
With respect to sample size, the involvement of 12 horses seems limited compared to rodent training studies; however, because of the intensive management and costs associated with equine training studies, it is practically not realizable to include many more horses in a highly standardized training study. When looking in literature, most existing training studies performed on horses include five to maximum 10 horses (31, 65).
Hodgson and coworkers analyzed muscle fiber type composition and capillarization, together with mitochondrial enzymatic activity in the m. gluteus medius of seven horses before and after 9 months of endurance training (5). They found an increase in the percentage of type IIX fibers showing an increased concentration of mitochondria along with an increase in the number of capillaries surrounding these fibers (5). An increase in mitochondrial enzymatic activity was already observed after 3 months of training (5). Karlström et al. studied the training-induced shift in fiber type composition and CSA in the m. gluteus medius of 18 racehorses in order to compare between young (2 years old) and more aged trained horses (4–8 years old) without subjecting these horses to a standardized training protocol (6). They found in the more aged and trained horses a larger CSA for fiber type I and a smaller CSA for fiber type IIX, as well as a more numerous representation of fiber type IIA and fewer fiber type IIX compared to younger horses. After performing NADH tetrazolium reductase staining of muscle sections to measure dehydrogenase activity, an age- and training-dependent increase of the oxidative capacity of type IIX fibers was also observed (6). Revold et al. analyzed training-induced shifts in muscle fiber type composition and CSA together with metabolic changes in the m. gluteus medius of nine young Coldblooded trotter horses (7). Despite the lack of standardized training, these researchers found a training-dependent decrease in fibers of type IIX and an increase in fibers of type IIAX, together with an increase in cytoplasmic MCT 1 and basigin expression in fibers of respectively type IIAX and IIA, and an increase in CS and HAD of the m. gluteus medius (7).
Studies that have involved multiple muscle groups predominantly focus on comparing muscle plasticity parameters between muscle groups under baseline conditions (60, 61). For instance, Karlström et al. explored the baseline differences in fiber type composition, capillarization, and enzymatic profile between 10 locomotor muscles and 3 non-locomotor muscles in three Swedish Standardbred untrained horses (60). In that study, it has been shown that non-locomotor muscles such as the m. masseter and diaphragm have a high percentage of type I fibers (70–100%) and almost no type IIX fibers (<1%), and this was accompanied by a high oxidative and low glycolytic capacity, in opposite of what was seen in locomotor muscles such as m. vastus lateralis and semitendinosus. Regarding locomotor muscles, the m. vastus lateralis showed a higher percentage of type I fibers and a lower percentage of type IIX fibers compared to the m. semitendinosus (60). Essén and partners described the effect of age associated with unspecified training in 55 Standardbred horses divided into five different age groups, chosen to represent ages of interest according to muscle growth and activity (61). They analyzed the baseline fiber type composition of four different locomotor muscles: the m. gluteus medius, the m. vastus lateralis, the m. semitendinosus, and the m. triceps brachii. They found a predominant presence of fibers of type II in all horses, with the m. semitendinosus having a higher proportion of this fiber type compared to the other muscles. The m. gluteus medius showed a higher type I/type II ratio in older horses compared to foals and the m. vastus lateralis had a higher content in type IIA compared to other muscles, whereas the m. semitendinosus presented a higher content of type IIX fibers. The glycolytic capacity was found to be higher in the m. gluteus medius and the m. semitendinosus, but with age, a higher oxidative capacity was seen to develop in these two muscles. This study used the ATPase staining method instead of the more widely used myosin heavy chain (MHC) method (61). In a different study, Grotmol and his team also examined the baseline fiber type composition and also the spatial distribution of muscle fibers in different muscle layers of the m. gluteus medius and the m. semitendinosus in four American Standardbred horses that underwent unspecified racetrack training (63). The results showed a predominance of type IIX fibers at the surface of the fascicles, and of types I and IIA in the more central parts. This type II/type I difference between the superficial and the deeper sites was especially more pronounced in the m. gluteus medius (63). In a different study, Kawai et al. described the baseline muscle fiber type composition and metabolic properties of 46 different muscles from the whole body in six Thoroughbred horses, including the m. pectoralis profundus, the m. vastus lateralis, and the m. semitendinosus. In the hindlimb muscles, the fiber type IIA was predominant, and the m. vastus lateralis was one with the highest proportion of this fiber type. Additionally, the mean percentage of type IIX fibers in the m. gluteus medius was the highest compared to the rest of the muscles in this group. The hindlimb group had overall the lowest percentage of type I and IIA fibers (64). Additionally, only a few studies have involved the comparison of two locomotor muscles before and after training (61–63). For example, Dingboom et al. studied in 38 young Dutch Warmblood foals divided into three exercise groups (box-rest, training, and pasture) the difference in muscle fiber type composition of the deep m. gluteus medius and the m. semitendinosus. In that study, no effect of exercise could be seen; however, changes in composition could be seen with age. From 22 to 48 weeks, a decrease in type IIX fibers was observed in the m. semitendinosus, whereas an age-dependent increase of type I and IIA fibers, as well as an increase in type IIAX fibers, was found in the m. gluteus medius (62).
In the current study, the multimodal approach allowed for obtaining a view on and comparison of baseline muscle fiber type composition and fiber type specific CSA, mean mathematical CSA, and mean frame CSA between all three studied muscle groups. With that respect, no baseline CSA differences could be detected. On the other hand, significant differences were detected in the baseline representation of type I and type IIX fibers, while no differences were found when looking at type IIA fibers, the predominant fiber type in horses.
Both the m. pectoralis and the m. vastus lateralis harbored significantly more type I fibers, when compared to the m. semitendinosus (66). The m. semitendinosus on its turn, showed significantly more type IIX fibers when compared to both the m. pectoralis and the m. vastus lateralis. This is in accordance with previous studies (60, 61, 63, 64). Based on the baseline fiber type composition, the most pronounced aerobic metabolic profile could be assigned to the m. pectoralis, the least aerobic to the m. semitendinosus. This is also the first study to show that both the mean mathematical and mean frame CSA is quite comparable between these muscle groups.
It is important to notice that the baseline fiber type composition of the m. vastus lateralis in the current study is comparable with that reported for the m. gluteus medius in Standardbred trotter horses in other studies (4, 67, 68). The ease with which biopsies can be harvested in the m. vastus lateralis is an added value for studies in which a longitudinal follow-up of several weeks to months is performed and in which for example also harvest of acute biopsies (taken immediately after exercise) is involved. Besides that, the m. vastus lateralis is often involved in human muscular studies, which could open up new avenues from a comparative point of view (69–71).
When looking at training adaptation, the m. vastus lateralis showed the most plasticity, which was manifested by a significant decrease in the fiber type specific CSA of all studied muscle fiber types and a significant decrease of the mean CSA, and this was applicable for both the frame and mathematical methods. Training adaptation also occurred in the m. semitendinosus, but to a much lesser extent. In that muscle, there was a significant increase in the percentage of type I fibers and a decrease in the fiber type specific CSA of type IIAX fibers. No significant change in the mean CSA of all fibers was detected with either method. The m. pectoralis showed no training adaptation with the applied histological methods. So, despite the fact that baseline fiber type specific CSA and mean CSA did not differ between all three studied muscle groups at the start, a clear training-induced response was visible, and this was the muscle plasticity parameter that showed most pronounced changes by far. It shows that monitoring of evolution of muscle fiber CSA (either fiber type specific or mean fiber CSA) can be very valuable to effectively underpin results obtained with, for example, “omics” studies. Of course, for that purpose, it is important to know which physiological meaning can be assigned to either an increase or a decrease of fiber CSA.
Essén-Gustavsson and Lindholm reported that the training-induced decrease in mean CSA in the m. gluteus medius was associated with increased citrate synthase activity (72). In a follow-up study, they reported a decrease in mean CSA in the m. vastus lateralis, to be associated with an increased oxidative and decreased glycolytic capacity (61). A smaller mean CSA in response to training is thought to be associated with a reduced glycogen storage capacity (18), and occurrence of muscle hyperplasia (muscle growth), as previously shown in Friesian horses subjected to aquatraining and dry treadmill training (73). In that study, it was shown that these morphological changes also coincide with an increased capacity to process alternative fuels such as branched-chain amino acids, aromatic amino acids, and gut microbiome-derived xenobiotics, which can more easily reach the muscular metabolic machinery, because of a smaller fiber type specific CSA. More research is needed with that respect.
The fact that no significant effect of harness training was seen on the m. pectoralis is not that surprising, since muscle kinematic studies clearly show that during trotting, the hind limb muscles show the greatest activity, creating propulsion (74). Eight weeks of dry treadmill training, another type of training, induces important changes in the m. pectoralis of Friesian horses, both from a histological and from an untargeted metabolomics point of view (22, 75–77). Dry treadmill training increases the muscle mass of the m. pectoralis, assessed by means of ultrasound, and induces a shift toward a more aerobic profile, which was further confirmed by the increase in type I and a decrease in type IIX fibers, together with a decrease in mean CSA. It seems that dry treadmill training challenges the m. pectoralis, whereas trotter harness training has a major impact on the m. vastus lateralis. This knowledge needs to be taken into account when designing muscle-specific training programs and optimizing the breed- and discipline-specific management of equine athletes. On another note, it is important to keep in mind that despite the absence of a detectable change in fiber type composition as seen in the m. pectoralis and the m. vastus lateralis in response to training in the current study, alterations at the metabolic level may have already occurred. In line with this, Hodgson et al. did not observe changes in fiber type composition of the m. gluteus medius of Thoroughbred horses after up to 12 weeks of race training, but they did find a clear upregulation of the oxidative metabolism after 6 weeks of training, as indicated by an increased CS enzyme activity (78).
In the current study, we assessed the minimal Feret's diameter that has been applied in various human and rodent muscular studies (79–86) and few equine studies (87–89). It is interesting to notice that on both occasions (fiber type specific and overall minimum Feret's diameter), different results were obtained when compared to the fiber specific and mean CSA (frame and mathematical approach). In the m. semitendinosus, a decreased CSA of the type IIAX fibers vs. a decrease of the type I minimum Feret's diameter was obtained. Furthermore, both mean frame and mean mathematical CSA showed no changes in answer to training, whereas the minimum Feret's diameter approach showed a training-induced decrease across all muscle fiber types. It would be interesting to further examine the added value of involving the minimum Feret's diameter approach in future training studies.
In the m. vastus lateralis, which showed the most overall pronounced harness training-induced adaptation in comparison to the other two muscle groups, besides a decrease in fiber type specific CSA of all muscle fiber types and a decrease in overall mean CSA for both approaches (mathematical and frame CSA), this was accompanied by a significant increase in capillary supply. No changes in mitochondrial density were detected.
In the past, it was believed that only aerobic exercise could improve capillarization in order to bring oxygen to the muscle fibers and remove waste products. It is clear that also anaerobic exercise improves capillarization through angiogenesis (90, 91). In the current study, capillary density was significantly increased in the m. vastus lateralis after 8 weeks of harness training and the number of capillaries surrounding both type I and type IIA fibers was significantly increased (Table 4). This is in accordance with a study of Chanda and coworkers who showed an increased capillary density in polo ponies subjected to 22 weeks of low-to-moderate intensity exercise and polo competitions (29). However, looking at capillary density per se is not enough because it does not confirm anything with respect to angiogenesis if C/F remains unchanged. Therefore, using the C/F ratio is advisable since it takes the muscle fiber size into account and thus reflects true changes in the capillary bed. In our study, C/F increased with training and, combined with a decreased mean fiber CSA and an increased capillary density, this indicates that 8 weeks of harness training induced angiogenesis. This is in agreement with a study of Serrano et al., in which Andalusian horses were subjected to 8 months of endurance training followed by 3 months of detraining. In that study, capillary density increased in the m. gluteus medius after 8 months of training and returned to baseline after 3 months of detraining (92). The same conclusions were drawn by Rivero et al., who looked into the effect of different training protocols on muscle capillarization and concluded that all training regimens (respectively 5, 15, and 25 min at speeds at blood lactate concentrations of 2.5 and 4 mmol/L for 22 days) significantly increased mean capillary number in the m. gluteus medius, whereas C/F only increased in training protocols of 15 and 25 min, respectively, but not 5 min, which means capillary changes occur only after a minimal duration and load of training (31). When comparing the increase in capillary density of type I and type IIA fibers, respectively, the increase in capillaries surrounding type IIA fibers (1.21-fold) showed a more pronounced trend when compared with type I fibers (1.17-fold) (p = 0.10), suggesting that capillarization targets more type IIA fibers with the applied training protocol.
In most species, baseline mitochondrial density is higher in type I fibers when compared to type IIA fibers and the subsarcolemmal region is richer in mitochondria when compared to the intermyofibrillar region (44, 93–95). This was also seen in the current study, where mitochondrial density was higher in type I fibers. Besides that, baseline mitochondrial density was significantly higher for the SS population when compared to the IMF population, in both fiber types of the m. vastus lateralis (Table 6).
In the current study, no training-induced shifts in mitochondrial density were observed in the m. vastus lateralis (Table 5). Also, other studies have reported such results, though those studies involved additional parameters to assess mitochondrial functional activity and they did, for example, find increased mitochondrial respiration despite lack of changes in histological mitochondrial density (46, 52). It is well known that early adaptations occur at the mitochondrial gene level in answer to exercise (96, 97), whereas phenotypic mitochondrial changes can occur at a later stage (36, 51). Exercise upregulates mitochondrial-related genes, such as genes involved in mitochondrial biogenesis, as has been shown in a recent study (98). In that study, it was also shown that gut microbiome composition was associated with mitochondrial gene expression in race horses, a finding that underlines the importance of the gut–muscle axis (96, 98–100). A possible explanation for the lack of change in mitochondrial density seen in the current training study could be that, in horses, the first steps of physiological training adaptation are reflected by early metabolic changes and gene expression, followed by changes in capillarization and shifts in muscle fiber CSA, prior to increases in mitochondrial density. This would be the preferential scenario for an accelerated supply of, for example, gut-derived xenobiotics to the muscular oxidative machinery.
The observed training-induced changes in the current study suggest that the muscular machinery models itself toward an optimal smaller individual muscle fiber structure to receive and process fuels that can be swiftly delivered by the circulatory system and do not require a major transformation of the mitochondrial machinery. In a previous training study, we have shown that branched-chain amino acids, aromatic amino acids, and gut microbiome-related xenobiotics function as important fuels (73). These metabolites also showed remarkable changes in other previous training studies, however, without receiving further attention (40, 101–103). Most probably, these fuels are processed in both type I and type IIA fibers. Unlike glycogen and fatty acids, these fuels are not stored within muscle fibers and need a rapid gut–muscle gateway to reach these fibers. It is more than feasible that the circulatory system may act as a conduit to provide a continuous flow of these gut microbiome-derived metabolites toward the trained muscle fibers. As a consequence, the storage of fuels is unnecessary because the resource is inexhaustible and transport is ensured, and the stress on the mitochondrial system is modest since these fuels can feed the TCA cycle at steps further downstream from acetyl-CoA. More research is needed with that respect since this knowledge has an important impact on existing paradigms in training and feeding management.
Conclusions
This is the first training study to involve four muscle plasticity parameters in three muscles with a different function (posture vs. locomotion). Though no baseline (pre-training) fiber type specific CSA and overall mean CSA differences could be detected between all muscle groups, there was a clear decrease in fiber type specific CSA, and this was most pronounced for the m. vastus lateralis and was accompanied by a clear increase in capillary supply. The two other muscle plasticity parameters, shifts in fiber type composition and mitochondrial density, showed little to no training-induced changes. Therefore, it will be interesting to involve both muscle fiber CSA and capillary supply in follow-up harness training studies that involve “omics” analysis techniques, which are deemed to be very sensitive, but difficult to link to solid physiological deductions, without involvement of conventional techniques.
Data Availability Statement
The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author/s.
Ethics Statement
The animal study was reviewed and approved by Animal Ethics Committee of the Ghent University (EC 2016/40).
Author Contributions
CM, BB, LM, and CD were responsible for the conception and design of the research. CM, BB, CV, LM, MO, and CD performed the experiments. CM, BB, LM, CV, CD, and KG performed the analyses. CM, BB, LM, CV, WD, LL, and CD interpreted the results. CM, BB, CV, and CD drafted the manuscript. CM, BB, MO, LL, WD, KG, and CD edited and revised the manuscript. All authors read and approved the final version of the manuscript and accountable for all aspects of the work.
Funding
This study was partly funded by the Fonds Wetenschappelijk Onderzoek (FWO) Flanders, Belgium (grant number 1S57617N). The experiments were performed at the Department of Translational Physiology, Infectiology and Public Health. Research Group of Comparative Physiology, Faculty of Veterinary Medicine, Ghent University, Belgium and at a training facility (EC 2016/40).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher's Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1. Flück M. Functional, structural and molecular plasticity of mammalian skeletal muscle in response to exercise stimuli. J Exp Biol. (2006) 209:2239–48. doi: 10.1242/jeb.02149
2. Flück M, Hoppeler H. Molecular basis of skeletal muscle plasticity–from gene to form and function. In: Reviews of Physiology, Biochemistry and Pharmacology., vol 146. Berlin, Heidelberg: Springer (2003). p. 159–216. doi: 10.1007/s10254-002-0004-7
3. Brook MS, Wilkinson DJ, Phillips BE, Perez-Schindler J, Philp A, Smith K, et al. Skeletal muscle homeostasis and plasticity in youth and ageing: Impact of nutrition and exercise. Review. (2016) 216:15–41. doi: 10.1111/apha.12532
4. Tyler CM, Golland LC, Evans DL, Hodgson DR, Rose RJ. Skeletal muscle adaptations to prolonged training, overtraining and detraining in horses. Pflügers Arch. (1998) 436:391–7. doi: 10.1007/s004240050648
5. Hodgson DR, Rose RJ. Effects of a nine-month endurance training programme on muscle composition in the horse. Vet Rec. (1987) 121:271–4. doi: 10.1136/vr.121.12.271
6. Karlström K, Lindholm A, Collinder E, Essén-Gustavsson B. Muscle fibre type composition in young and racing Swedish cold-blooded trotters. Comp Exerc Physiol. (2009) 6:27–32. doi: 10.1017/S1478061509356145
7. Revold T, Mykkänen AK, Karlström K, Ihler CF, Pösö AR, Essén-Gustavsson B. Effects of training on equine muscle fibres and monocarboxylate transporters in young Coldblooded Trotters. Equine Vet J. (2010) 42:289–95. doi: 10.1111/j.2042-3306.2010.00274.x
8. Medler S. Mixing it up: The biological significance of hybrid skeletal muscle fibers. J Exp Biol. (2019) 222:jeb200832. doi: 10.1242/jeb.200832
9. Rivero J-LL, Hill EW. Skeletal muscle adaptations and muscle genomics of performance horses. Vet J. (2016) 209:5–13. doi: 10.1016/j.tvjl.2015.11.019
10. Barlow DA, Lloyd JM, Hellhake P, Seder JA. Equine muscle fiber types: A histological and histochemical analysis of select thoroughbred yearlings. J Equine Vet Sci. (1984) 4:60–6. doi: 10.1016/S0737-0806(84)80083-0
11. Simoneau JA, Bouchard C. Human variation in skeletal muscle fiber-type proportion and enzyme activities. Am J Physiol - Endocrinol Metab. (1989) 257:567–72. doi: 10.1152/ajpendo.1989.257.4.e567
12. Zierath JR, Hawley JA. Skeletal Muscle Fiber Type: Influence on Contractile and Metabolic Properties. PLoS Biol. (2004) 2:e348. doi: 10.1371/journal.pbio.0020348
13. Baldwin KM, Haddad F. Invited Review: Effects of different activity and inactivity paradigms on myosin heavy chain gene expression in striated muscle. J Appl Physiol. (2001) 90:345–57. doi: 10.1152/jappl.2001.90.1.345
14. Acevedo LM, Rivero J-LL. New insights into skeletal muscle fibre types in the dog with particular focus towards hybrid myosin phenotypes. Cell Tissue Res. (2006) 323:283–303. doi: 10.1007/s00441-005-0057-4
15. Methenitis S, Karandreas N, Spengos K, Zaras N, Stasinaki AN, Terzis G. Muscle fiber conduction velocity, muscle fiber composition, and power performance. Med Sci Sports Exerc. (2016) 48:1761–71. doi: 10.1249/MSS.0000000000000954
16. Rivero JL, Serrano AL, Henckel P, Aguera E. Muscle fiber type composition and fiber size in successfully and unsuccessfully endurance-raced horses. J Appl Physiol. (1993) 75:1758–66. doi: 10.1152/jappl.1993.75.4.1758
17. Herman JR, Rana SR, Chleboun GS, Gilders RM, Hageman FC, Hikida RS, et al. Correlation between muscle fiber cross-sectional area and strength gain using three different resistance-training programs in college-aged women. J Strength Cond Res. (2010) 24. doi: 10.1097/01.JSC.0000367128.04768.0a
18. Lefaucheur L, Lebret B, Ecolan P, Louveau I, Damon M, Prunier A, et al. Muscle characteristics and meat quality traits are affected by divergent selection on residual feed intake in pigs. J Anim Sci. (2010) 89:996–1010. doi: 10.2527/jas.2010-3493
19. Votion D-M, Navet R, Lacombe VA, Sluse F, Essén-Gustavsson B, Hinchcliff KW, et al. Muscle energetics in exercising horses. Equine Comp Exerc Physiol. (2007) 4:105–118. doi: 10.1017/S1478061507853667
20. Palencia P, Quiroz-Rothe E, Rivero JLL. New insights into the skeletal muscle phenotype of equine motor neuron disease: A quantitative approach. Acta Neuropathol. (2005) 109:272–84. doi: 10.1007/s00401-004-0940-1
21. Kinsey ST, Locke BR, Dillaman RM. Molecules in motion: influences of diffusion on metabolic structure and function in skeletal muscle. J Exp Biol. (2011) 214:263–74. doi: 10.1242/jeb.047985
22. Boshuizen B, De meeûs C, van de Winkel D, van Hauwe L, de Bruijn M, Touwen N, et al. The effect of treadmill training on equine muscle morphometrics and muscle metabolomics. Comp Exerc Physiol. (2018) 14:4.
23. Lindner A, Dag Erginsoy S, Kissenbeck S, Mosen H, Hetzel U, Drommer W, et al. Effect of different blood-guided conditioning programmes on skeletal muscle ultrastructure and histochemistry of sport horses. J Anim Physiol Anim Nutr (Berl). (2013) 97:374–86. doi: 10.1111/j.1439-0396.2012.01283.x
24. Karlström K, Essén-Gustavsson B, Lindholm A, Persson SGB. Capillary supply in relation to muscle metabolic profile and cardiocirculatory parameters. Equine Exerc Physiol. (1991) 3:239–244.
25. McCall GE, Byrnes WC, Dickinson A, Pattany PM, Fleck SJ. Muscle fiber hypertrophy, hyperplasia, and capillary density in college men after resistance training. J Appl Physiol. (1996) 81:2004–12. doi: 10.1152/jappl.1996.81.5.2004
26. Andersen P, Henriksson J. Capillary supply of the quadriceps femoris muscle of man: adaptive response to exercise. J Physiol. (1977) 270:677–90. doi: 10.1113/jphysiol.1977.sp011975
27. Ingjer F. Capillary supply and mitochondrial content of different skeletal muscle fiber types in untrained and endurance-trained men. A histochemical and ultrastructural study. Eur J Appl Physiol Occup Physiol. (1979) 40:197–209. doi: 10.1007/BF00426942
28. Hoier B, Hellsten Y. Exercise-induced capillary growth in human skeletal muscle and the dynamics of VEGF. Microcirculation. (2014) 21:301–14. doi: 10.1111/micc.12117
29. Chanda M, Srikuea R, Cherdchutam W, Chairoungdua A, Piyachaturawat P. Modulating effects of exercise training regimen on skeletal muscle properties in female polo ponies. BMC Vet Res. (2016) 12:245. doi: 10.1186/s12917-016-0874-6
30. Rietbroek N, Dingboom E, Schuurman S, Wiel E, Eizema K, Everts M. Effect of exercise on development of capillary supply and oxidative capacity in skeletal muscle of horses. Am J Vet Res. (2007) 68:1226–31. doi: 10.2460/ajvr.68.11.1226
31. Rivero J-LL, Ruz A, Martí-Korff S, Estepa J-C, Aguilera-Tejero E, Werkman J, et al. Effects of intensity and duration of exercise on muscular responses to training of thoroughbred racehorses. J Appl Physiol. (2007) 102:1871–82. doi: 10.1152/japplphysiol.01093.2006
32. Egan B, Zierath JR. Exercise metabolism and the molecular regulation of skeletal muscle adaptation. Cell Metab. (2013) 17:162–84. doi: 10.1016/j.cmet.2012.12.012
33. Tonkonogi M, Sahlin K. Physical exercise and mitochondrial function in human skeletal muscle. Exerc Sport Sci Rev. (2002) 30:129–137. doi: 10.1097/00003677-200207000-00007
34. Granata C, Jamnick NA, Bishop DJ. Training-induced changes in mitochondrial content and respiratory function in human skeletal muscle. Sport Med. (2018) 48:1809–28. doi: 10.1007/s40279-018-0936-y
35. Jacobs RA, Flück D, Bonne TC, Bürgi S, Christensen PM, Toigo M, et al. Improvements in exercise performance with high-intensity interval training coincide with an increase in skeletal muscle mitochondrial content and function. J Appl Physiol. (2013) 115:785–93. doi: 10.1152/japplphysiol.00445.2013
36. Menshikova E V, Ritov VB, Fairfull L, Ferrell RE, Kelley DE, Goodpaster BH. Effects of exercise on mitochondrial content and function in aging human skeletal muscle. Journals Gerontol - Ser A Biol Sci Med Sci. (2006) 61:534–40. doi: 10.1093/gerona/61.6.534
37. Phielix E, Meex R, Moonen-Kornips E, Hesselink MKC, Schrauwen P. Exercise training increases mitochondrial content and ex vivo mitochondrial function similarly in patients with type 2 diabetes and in control individuals. Diabetologia. (2010) 53:1714–21. doi: 10.1007/s00125-010-1764-2
38. Hoppeler H, Howald H, Conley K, Lindstedt SL, Claassen H, Vock P, et al. Endurance training in humans: Aerobic capacity and structure of skeletal muscle. J Appl Physiol. (1985) 59:320–7. doi: 10.1152/jappl.1985.59.2.320
39. Hey-Mogensen M, Højlund K, Vind BF, Wang L, Dela F, Beck-Nielsen H, et al. Effect of physical training on mitochondrial respiration and reactive oxygen species release in skeletal muscle in patients with obesity and type 2 diabetes. Diabetologia. (2010) 53:1976–85. doi: 10.1007/s00125-010-1813-x
40. Klein DJ, McKeever KH, Mirek ET, Anthony TG. Metabolomic response of equine skeletal muscle to acute fatiguing exercise and training. Front Physiol. (2020) 11:1–15. doi: 10.3389/fphys.2020.00110
41. Kühlbrandt W, Fox, Thomas D, Kenyon C, Larsson N-G, Nunnari J, O'Farrel P, et al. Chapter 14 energy conversion: mitochondria and chloroplasts. 6th ed. Molecular Biology of the Cell. 6th ed. Abingdom: Garland Science, Taylor & Francis Group, LLC (2015). p. 753–811.
42. Harwig MC, Viana MP, Egner JM, Harwig JJ, Widlansky ME, Rafelski SM, et al. Methods for imaging mammalian mitochondrial morphology: A prospective on MitoGraph. Anal Biochem. (2018) 552:81–99. doi: 10.1016/j.ab.2018.02.022
43. Skulachev VP. Mitochondrial filaments and clusters as intracellular power-transmitting cables. Trends Biochemi Sci. (2001) 23–29. doi: 10.1016/S0968-0004(00)01735-7
44. Chilibeck PD, Syrotuik DG, Bell GJ. The effect of strength training on estimates of mitochondrial density and distribution throughout muscle fibres. Eur J Appl Physiol Occup Physiol. (1999) 80:604–9. doi: 10.1007/s004210050641
45. Hawley JA. Adaptations of skeletal muscle to prolonged, intense endurance training. Clin Exp Pharmacol Physiol. (2002) 29:218–22. doi: 10.1046/j.1440-1681.2002.03623.x
46. White SH, Warren LK Li C, Wohlgemuth SE. Submaximal exercise training improves mitochondrial efficiency in the gluteus medius but not in the triceps brachii of young equine athletes. Sci Rep. (2017) 7:14389. doi: 10.1038/s41598-017-14691-4
47. Geof RJ, McCutcheon LJ, Shen HUA. Muscular and metabolic responses to moderate-intensity short-term training. Equine Vet J. (1999) 31:311–7. doi: 10.1111/j.2042-3306.1999.tb05240.x
48. McCutcheon LJ, Byrd SK, Hodgson DR. Ultrastructural changes in skeletal muscle after fatiguing exercise. J Appl Physiol. (1992) 72:1111–7. doi: 10.1152/jappl.1992.72.3.1111
49. Kayar SR, Hoppeler H, Mermod L, Weibel ER. Mitochondrial size and shape in equine skeletal muscle: A three-dimensional reconstruction study. Anat Rec. (1988) 222:333–9. doi: 10.1002/ar.1092220405
50. Prince FP, Hikida RS, Hagerman FC, Staron RS, Allen WH. A. morphometric analysis of human muscle fibers with relation to fiber types and adaptations to exercise. J Neurol Sci. (1981) 49:165–79. doi: 10.1016/0022-510X(81)90076-9
51. Axelrod CL, Fealy CE, Mulya A, Kirwan JP. Exercise training remodels human skeletal muscle mitochondrial fission and fusion machinery towards a pro-elongation phenotype. Acta Physiol. (2019) 225:13216. doi: 10.1111/apha.13216
52. White SH, Wohlgemuth S, Li C, Warren LK. Rapid Communication: Dietary selenium improves skeletal muscle mitochondrial biogenesis in young equine athletes. J Anim Sci. (2017) 95:4078–84. doi: 10.2527/jas2017.1919
53. Desplanches D, Amami M, Dupré-Aucouturier S, Valdivieso P, Schmutz S, Mueller M, et al. Hypoxia refines plasticity of mitochondrial respiration to repeated muscle work. Eur J Appl Physiol. (2014) 114:405–17. doi: 10.1007/s00421-013-2783-8
54. Votion DM, Fraipont A, Goachet AG, Robert C, Van Erck E, Amory H, et al. Alterations in mitochondrial respiratory function in response to endurance training and endurance racing. Equine Vet J. (2010) 42:268–74. doi: 10.1111/j.2042-3306.2010.00271.x
55. Sewell DA, Harris RC, Marlin DJ. Skeletal muscle characteristics in 2 year-old race-trained thoroughbred horses. Comp Biochem Physiol Part A Physiol. (1994) 108:87–96. doi: 10.1016/0300-9629(94)90059-0
56. Zsoldos RR, Voegele A, Krueger B, Schroeder U, Weber A, Licka TF. Long term consistency and location specificity of equine gluteus medius muscle activity during locomotion on the treadmill. BMC Vet Res. (2018) 14. doi: 10.1186/s12917-018-1443-y
57. Latham CM, White SH. Validation of primary antibodies for multiple immunofluorescent labeling of horse skeletal muscle fiber type. J Anim Sci. (2017) 95:53. doi: 10.2527/asasann.2017.107
58. Hodgson DR, Rose RJ, DiMauro J, Allen JR. Effects of a submaximal treadmill training programme on histochemical properties, enzyme activities and glycogen utilisation of skeletal muscle in the horse. Equine Vet J. (1985) 17:300–5. doi: 10.1111/j.2042-3306.1985.tb02504.x
59. Henckel P. A histochemical assessment of the capillary blood supply of the middle gluteal muscle of Thoroughbred horses. Equine Exerc Physiol. (1983) 225–8.
60. Karlstrom K, Essen-Gustavsson B, Lindholm A. Fibre type distribution, capillarization and enzymatic profile of locomotor and nonlocomotor muscles of horses and steers. Acta Anat (Basel). (1994) 151:97–106. doi: 10.1159/000147649
61. Essén B, Lindholm A, Thornton J. Histochemical properties of muscle fibre types and enzyme activities in skeletal muscles of Standardbred trotters of different ages. Equine Vet J. (1980) 12:175–80. doi: 10.1111/j.2042-3306.1980.tb03420.x
62. Dingboom EG, Van Oudheusden H, Eizema K, Weijs WA. Changes in fibre type composition of gluteus medius and semitendinosus muscles of Dutch Warmblood foals and the effect of exercise during the first year postpartum. Equine Vet J. (2002) 34:177–83. doi: 10.2746/042516402776767312
63. Grotmol S, Totland GK, Kryvi H, Breistøl A, Essén-Gustavsson B, Lindholm A. Spatial distribution of fiber types within skeletal muscle fascicles from Standardbred horses. Anat Rec. (2002) 268:131–6. doi: 10.1002/ar.10140
64. Kawai M, Minami Y, Sayama Y, Kuwano A, Hiraga A, Miyata H. Muscle fiber population and biochemical properties of whole body muscles in thoroughbred horses. Anat Rec. (2009) 292:1663–9. doi: 10.1002/ar.20961
65. Yamano S, Kawai M, Minami Y, Hiraga A, Miyata H. Differences in muscle fiber recruitment patterns between continuous and interval exercises. J Equine Sci. (2010) 21:59–65. doi: 10.1294/jes.21.59
66. Payne RC, Hutchinson JR, Robilliard JJ, Smith NC, Wilson AM. Functional specialisation of pelvic limb anatomy in horses (Equus caballus). J Anat. (2005) 206:557–74. doi: 10.1111/j.1469-7580.2005.00420.x
67. Bechtel PJ, Kline KH. Muscle fiber type changes in the middle gluteal of Quarter and Standardbred horses from birth through one year of age. In: Equine Exercise Physiology Proceedings of the First International Conference. (1983). p. 265–70. Available online at: http://www.iceep.org/pdf/iceep2/_1129101133_001.pdf
68. Ronéus M, Essén-Gustavsson B, Arnason T. Racing performance and longitudinal changes in muscle characteristics in standardbred trotters. J Equine Vet Sci. (1993) 13:355–61. doi: 10.1016/S0737-0806(06)81124-X
69. Phung LA, Foster AD, Miller MS, Lowe DA, Thomas DD. Super-relaxed state of myosin in human skeletal muscle is fiber-type dependent. Am J Physiol - Cell Physiol. (2020) 319:C1158–62. doi: 10.1152/ajpcell.00396.2020
70. Nederveen JP, Ibrahim G, Fortino SA, Snijders T, Kumbhare D, Parise G. Variability in skeletal muscle fibre characteristics during repeated muscle biopsy sampling in human vastus lateralis. Appl Physiol Nutr Metab. (2020) 45:368–75. doi: 10.1139/apnm-2019-0263
71. Jeon Y, Choi J, Kim HJ, Lee H, Lim JY, Choi SJ. Sex- and fiber-type-related contractile properties in human single muscle fiber. J Exerc Rehabil. (2019) 15:537–45. doi: 10.12965/jer.1938336.168
72. Essén-Gustavsson B, Lindholm A. Muscle fibre characteristics of active and inactive Standardbred horses. Equine Vet J. (1985) 17:434–8. doi: 10.1111/j.2042-3306.1985.tb02549.x
73. de Meeûs D'Argenteuil C, Boshuizen B, Oosterlinck M, van de Winkel D, De Spiegelaere W, et al. Flexibility of equine bioenergetics and muscle plasticity in response to different types of training : An integrative approach, questioning existing paradigms. PLoS ONE. (2021) 16:e0249922. doi: 10.1371/journal.pone.0249922
74. Robert C, Valette JP, Pourcelot P, Audigié F, Denoix JM. Effects of trotting speed on muscle activity and kinematics in saddlehorses. Equine Vet J Suppl. (2002) 295–301. doi: 10.1111/j.2042-3306.2002.tb05436.x
75. De Meeus d'Argenteuil C, Boshuizen B, Van de Winkel D, Van Hauwe L, de Bruijn M, Touwen N, et al. The effect of aquatraining and dry treadmill training on muscle morphometric and muscle metabolomics in horses. Equine Vet J Suppl. (2018) 50:5–35. doi: 10.1111/evj.26_13008
76. De meeûs C, Van de Winkel D, Boshuizen B, Loose D, de Bruijn M, Touwen N, et al. Longitudinal follow-up of equine muscle morphometrics and associated metabolic properties induced by 8 weeks of treadmill training. In: Physiological Bioenergetics : Mitochondria From Bench to Bedside, 2nd American Physiological Society Conference. San Diego, CA, USA: American Physiologycal Society (2017). Available online at: http://hdl.handle.net/1854/LU-8546474
77. De meeûs C, Boshuizen B, Van de Winkel D, Van Hauwe L, de Bruijn M, Touwen N, et al. Longitudinal follow-up of equine muscle morphometrics and associated metabolic properties induced by 8 weeks of treadmill training. In: Abstracts European Veterinary Conference Voorjaarsdagen 2018. The Haghe, The Netherlands: European Veterinary Conference (2018). Available online at: http://hdl.handle.net/1854/LU-8560946
78. Hodgson DR, Rose RJ, Dimauro J, Allen JR. Effects of training on muscle composition in horses. Am J Vet Res. (1986) 47:12–15.
79. Pertl C, Eblenkamp M, Pertl A, Pfeifer S, Wintermantel E, Lochmüller H, et al. A new web-based method for automated analysis of muscle histology. BMC Musculoskelet Disord. (2013) 14:26. doi: 10.1186/1471-2474-14-26
80. Briguet A, Courdier-Fruh I, Foster M, Meier T, Magyar JP. Histological parameters for the quantitative assessment of muscular dystrophy in the mdx-mouse. Neuromuscul Disord. (2004) 14:675–82. doi: 10.1016/j.nmd.2004.06.008
81. Birkbeck MG, Heskamp L, Schofield IS, Blamire AM, Whittaker RG. Non-invasive imaging of single human motor units. Clin Neurophysiol. (2020) 131:1399–406. doi: 10.1016/j.clinph.2020.02.004
82. Sayed RKA, de Leonardis EC, Guerrero-Martínez JA, Rahim I, Mokhtar DM, Saleh AM, et al. Identification of morphological markers of sarcopenia at early stage of aging in skeletal muscle of mice. Exp Gerontol. (2016) 83:22–30. doi: 10.1016/j.exger.2016.07.007
83. Mayeuf-Louchart A, Hardy D, Thorel Q, Roux P, Gueniot L, Briand D, et al. A high-content analysis method to study skeletal muscle with a new Fiji tool. Skelet Muscle. (2018) 8:25. doi: 10.1186/s13395-018-0171-0
84. Lin C-H, Hsu H-C, Hou Y-J, Chen K-H, Lai S-H, Chang W-M. Relationship between sonography of sternocleidomastoid muscle and cervical passive range of motion in infants with congenital muscular torticollis. Biomed J. (2018) 41:369–75. doi: 10.1016/j.bj.2018.10.001
85. Lau YS, Xu L, Gao Y, Han R. Automated muscle histopathology analysis using CellProfiler. Skelet Muscle. (2018) 8:32. doi: 10.1186/s13395-018-0178-6
86. De Meulemeester K, Calders P, Van Dorpe J, De Pauw R, Petrovic M, Cagnie B. Morphological differences in the upper trapezius muscle between female office workers with and without trapezius myalgia: facts or fiction?: A cross-sectional study. Am J Phys Med Rehabil. (2019) 98:117–124. doi: 10.1097/PHM.0000000000001029
87. Reesink HL, Hermanson JW, Cheetham J, Mu L, Mitchell LM, Soderholm LV, et al. Anatomic and neuromuscular characterisation of the equine cricothyroid muscle. Equine Vet J. (2013) 45:630–6. doi: 10.1111/evj.12023
88. Chalmers HJ, Caswell J, Perkins J, Goodwin D, Viel L, Ducharme NG, et al. Ultrasonography detects early laryngeal muscle atrophy in an equine neurectomy model. Muscle Nerve. (2016) 53:583–92. doi: 10.1002/mus.24785
89. Steel CM, Walmsley EA, Anderson GA, Coles CA, Ahern B, White JD. Immunohistochemical analysis of laryngeal muscle of horses clinically affected with recurrent laryngeal neuropathy. Equine Vet J. (2020) 53:710–17. doi: 10.1111/evj.13362
90. Wells GD, Selvadurai H, Tein I. Bioenergetic provision of energy for muscular activity. Paediatr Respir Rev. (2009) 10:83–90. doi: 10.1016/j.prrv.2009.04.005
91. Gavin TP. Basal and exercise-induced regulation of skeletal muscle capillarization. Exerc Sport Sci Rev. (2009) 37:86–92. doi: 10.1097/JES.0b013e31819c2e9b
92. Serrano AL, Quiroz-Rothe E, Rivero JLL. Early and long-term changes of equine skeletal muscle in response to endurance training and detraining. Pflugers Arch Eur J Physiol. (2000) 441:263–74. doi: 10.1007/s004240000408
93. Sullivan SM, Pittman RN. Relationship between mitochondrial volume density and capillarity in hamster muscles. Am J Physiol Circ Physiol. (1987) 252:149–55. doi: 10.1152/ajpheart.1987.252.1.H149
94. Shaw CS, Jones DA, Wagenmakers AJM. Network distribution of mitochondria and lipid droplets in human muscle fibres. Histochem Cell Biol. (2008) 129:65–72. doi: 10.1007/s00418-007-0349-8
95. van Wessel T, de Haan A, van der Laarse WJ, Jaspers RT. The muscle fiber type-fiber size paradox: hypertrophy or oxidative metabolism? Eur J Appl Physiol. (2010) 110:665–694. doi: 10.1007/s00421-010-1545-0
96. Perry CGR, Lally J, Holloway GP, Heigenhauser GJF, Bonen A, Spriet LL. Repeated transient mRNA bursts precede increases in transcriptional and mitochondrial proteins during training in human skeletal muscle. J Physiol. (2010) 588:4795–810. doi: 10.1113/jphysiol.2010.199448
97. Sriwijitkamol A, Ivy JL, Christ-Roberts C, DeFronzo RA, Mandarino LJ, Musi N. LKB1-AMPK signaling in muscle from obese insulin-resistant Zucker rats and effects of training. Am J Physiol Metab. (2006) 290:E925–E932. doi: 10.1152/ajpendo.00429.2005
98. Mach N, Moroldo M, Rau A, Lecardonnel J, Le Moyec L, Robert C, et al. Understanding the holobiont: crosstalk between gut microbiota and mitochondria during long exercise in horse. Front Molecular Biosci. (2021) 126. doi: 10.3389/fmolb.2021.656204
99. Lahiri S, Kim H, Garcia-Perez I, Reza MM, Martin KA, Kundu P, et al. The gut microbiota influences skeletal muscle mass and function in mice. Sci Transl Med. (2019) 11:eaan5662. doi: 10.1126/scitranslmed.aan5662
100. Lustgarten MS. The Role of the Gut Microbiome on Skeletal Muscle Mass and Physical Function: 2019 Update. Front Physiol. (2019) 1435. doi: 10.3389/fphys.2019.01435
101. Mach N, Ramayo-Caldas Y, Clark A, Moroldo M, Robert C, Barrey E, et al. Understanding the response to endurance exercise using a systems biology approach: combining blood metabolomics, transcriptomics and miRNomics in horses. BMC Genomics. (2017) 18:187. doi: 10.1186/s12864-017-3571-3
102. Plancade S, Clark A, Philippe C, Helbling J-C, Moisan M-P, Esquerré D, et al. Unraveling the effects of the gut microbiota composition and function on horse endurance physiology. Sci Rep. (2019) 9:9620. doi: 10.1038/s41598-019-46118-7
103. Le Moyec L, Robert C, Triba MN, Bouchemal N, Mach N, Rivière J, et al. A First Step Toward Unraveling the Energy Metabolism in Endurance Horses: Comparison of Plasma Nuclear Magnetic Resonance Metabolomic Profiles Before and After Different Endurance Race Distances. Front Mol Biosci. (2019) 6:45. doi: 10.3389/fmolb.2019.00045
104. R Development Core Team (2019). R: A Language and Environment for Statistical Computing. Vienna, Austria: R Foundation for Statistical Computing. Available online at: https://www.R-project.org/
Keywords: exercise, metabolism, muscle fiber typing, capillarization, horse (Equus caballus), mitochondria, muscle plasticity, fuel
Citation: de Meeûs d'Argenteuil C, Boshuizen B, Vidal Moreno de Vega C, Leybaert L, de Maré L, Goethals K, De Spiegelaere W, Oosterlinck M and Delesalle C (2021) Comparison of Shifts in Skeletal Muscle Plasticity Parameters in Horses in Three Different Muscles, in Answer to 8 Weeks of Harness Training. Front. Vet. Sci. 8:718866. doi: 10.3389/fvets.2021.718866
Received: 01 June 2021; Accepted: 07 September 2021;
Published: 18 October 2021.
Edited by:
Barbara Padalino, University of Bologna, ItalyReviewed by:
F. Capela e Silva, University of Evora, PortugalGeorge M. Strain, Louisiana State University, United States
Copyright © 2021 de Meeûs d'Argenteuil, Boshuizen, Vidal Moreno de Vega, Leybaert, de Maré, Goethals, De Spiegelaere, Oosterlinck and Delesalle. 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: Cathérine Delesalle, Catherine.Delesalle@UGent.be
†These authors share first authorship