- 1Laboratory of Animal (Poultry) Genetics Breeding and Reproduction, Ministry of Agriculture, Institute of Animal Science, Chinese Academy of Agricultural Sciences (CAAS), Beijing, China
- 2College of Animal Science, Xinjiang Agricultural University, Urumqi, China
- 3CAAS-ILRI Joint Laboratory on Livestock and Forage Genetic Resources, Institute of Animal Science, Chinese Academy of Agricultural Sciences (CAAS), Beijing, China
- 4Centre d'Anthropobiologie et de Génomique de Toulouse, Université Paul Sabatier, Toulouse, France
With their enormous muscle mass and athletic ability, horses are well-positioned as model organisms for understanding muscle metabolism. There are two different types of horse breeds—Guanzhong (GZ) horses, an athletic breed with a larger body height (~148.7 cm), and the Ningqiang pony (NQ) horses, a lower height breed generally used for ornamental purposes—both inhabited in the same region of China with obvious differences in muscle content. The main objective of this study was to evaluate the breed-specific mechanisms controlling muscle metabolism. In this study, we observed muscle glycogen, enzyme activities, and LC–MS/MS untargeted metabolomics in the gluteus medius muscle of six, each of GZ and NQ horses, to explore differentiated metabolites that are related to the development of two muscles. As expected, the glycogen content, citrate synthase, and hexokinase activity of muscle were significantly higher in GZ horses. To alleviate the false positive rate, we used both MS1 and MS2 ions for metabolite classification and differential analysis. As a result, a total of 51,535 MS1 and 541 MS2 metabolites were identified, and these metabolites can separate these two groups from each other. Notably, 40% of these metabolites were clustered into lipids and lipid-like molecules. Furthermore, 13 significant metabolites were differentially detected between GZ and NQ horses (fold change [FC] value ≥ 2, variable important in projection value ≥1, and Q value ≤ 0.05). They are primarily clustered into glutathione metabolism (GSH, p = 0.01), taurine, and hypotaurine metabolism (p < 0.05) pathways. Seven of the 13 metabolites were also found in thoroughbred racing horses, suggesting that metabolites related to antioxidants, amino acids, and lipids played a key role in the development of skeleton muscle in horses. Those metabolites related to muscle development shed a light on racing horses' routine maintenance and improvement of athletic performance.
1. Introduction
The excellent locomotion and racing performance of horses depend on their skeletal muscle metabolism. Therefore, the horse is a supremely effective animal model for understanding the effects of metabolism on muscle development (1). However, metabolite regulation mechanism associated with the training or development of horses remains unclear (2, 3). Previous studies have focused on diseases-related and exercise-related metabolites, for instance, the hematometabolic characteristics of blood metabolism in horses with atypical myopathy (AM) (4) and the identification of the biomarkers for acute laminitis by serum (5) metabolomics, and provide indexes for clinical diagnosis. Furthermore, the metabolomics of blood (6, 7) and muscle (8) for horses during endurance exercise reasons out the energetic demands for locomotion. Undoubtedly, metabolomics analysis has practical use for tracking stimulants from equine urine (9–11), and the analysis can also reflect exercise fatigue (12, 13) during exercise and training. Nevertheless, its utility is being increasingly acknowledged between different species [e.g., dogs (14, 15), chicken (16, 17), and cattle (18)]. Moreover, metabolomics has been applied to other livestock studies (19–21). Thus, metabolomics can afford opportunities to understand metabolic mechanisms between different breeds.
Many factors will affect muscle metabolism, such as exercise (8, 22), climate (23), and feed ingredients (24). To understand the muscle development related to breeds, we focused on two different horse breeds living in the Shaanxi Province of China that share similar herding management but exhibit different phenotypes (25). One is Guanzhong (GZ) horses, which have the characteristics of strong adaptability, tall body, easy feeding, fast growth and development, strong resistance to stress, and high reproduction rate (26), and the other one is the Ningqiang pony (NQ) horses, for which the average adult wither height is <106 cm as documented in a previous study (27). It adapts to the mountain areas with agile and small body size to carry heavy goods (28, 29). The variation in location between these two breeds requires muscle to provide power. Since the larger the muscle, the greater is its potential force output, skeletal muscle mass is important for a horse's performance and can vary greatly from breed to breed (30, 31). The huge difference in body size leads to different physiological changes in muscle composition and metabolism. There has been evidence that different breeds significantly affect muscle fiber structures, and regarding physiological changes, metabolism probably plays a key role during muscle formation (31).
Untargeted metabolomics is a comprehensive study that describes the entire endogenous metabolites in various organisms or tissues through the identification of metabolites in search of differential metabolites (32). The untargeted metabolomics can reflect the overall metabolite level and has been widely used to reveal the dynamic characteristics of metabolites under different conditions (33). As an approach to untargeted metabolomics, liquid chromatography–mass spectrometry (LC-MS) is considered to be the most stable and suitable technique for studying metabolism due to its wide coverage of substances (34). In the process of identifying metabolites, untargeted metabolomics may produce many false positive signals, when using the fragment ion of MS1 (35), whereas when combined with MS2, this bias can be alleviated. In this study, both MS1 and MS2 were used for the classification of metabolites and in the analysis; it was further enriched by the pathway, as reported in a previous study (36).
In the present study, we aimed to evaluate the breed-specific mechanisms related to horse muscle development. Thus, we performed biochemistry index detection and untargeted metabolomics analysis of muscle tissue from GZ and NQ horses. After the identification of metabolites, the significantly different metabolites between Ningqiang pony and Guanzhong pony muscles were analyzed. Our study would help understand the key components of metabolic networks and interactions during muscle formation. Elucidation of the underlying mechanisms contributes to the processes of feeding and exercise.
2. Materials and methods
2.1. Ethics statement
All animal experiments were conducted according to the regulations and guidelines established by the Animal Care Committee of the Institute of Animal Sciences, Chinese Academy of Agricultural Sciences (approval number: IAS2019-24). All incisions for muscle collection were performed under local anesthesia using 1% pentobarbital sodium, which is mainly used as an anesthetic agent for animals. Every effort was made to reduce animal suffering.
2.2. Sample collection by muscle biopsies
A total of 12 muscle samples were collected from Guanzhong horse (n = 6; male; age, 5.50 ± 2.42 years) and Ningqiang pony breeding center (n = 6; male; age, 4.89 ± 2.81 years), which are located in Shaanxi province, China (Supplementary Figure 1A). These two breeding centers share similar feeding and racing management, since they are near to each other (29). The feed ingredients are listed in Supplementary Table 1. Horses were weighed and examined by a clinician in a calm and quiet state prior to the muscle biopsy.
After preoperative disinfection with alcohol and iodophor, horses were subjected to local anesthesia (1 ml 1% pentobarbital sodium) in the gluteus medius (Supplementary Figure 1B). An appropriate amount of gluteus medius muscle was taken using an 8-mm biopsy punch (ACU-Punch®, Acuderm, Inc., Ft. Lauderdale, USA). The sample was placed into normal saline (0.9%NaCl solution) (Zhonghai Biotech, Beijing, China) to wash the blood, and tissues such as fat and fascia were removed. The final specimen was first stored in liquid nitrogen and then stored at −80°C until the metabolomic analysis. After sampling, we disinfected the wound with iodophor. For better postoperative recovery, Yunnan Baiyao powder was sprinkled over the incision until the dosage was sufficient to cover the wound (Yunnan Baiyao Group Co., Ltd, Kunming, China). Yunnan Baiyao [YNBY] is a traditional Chinese herbal medicinal formula widely used for its hemostatic effects and has been approved by the China Food and Drug Administration [CFDA] (37).
2.3. Detecting enzyme activity and muscle glycogen in the horse muscle
The glycogen content (nmol/L), hexokinase (HK) (mU/L) activity, and citrate synthase (CS) (U/L) activity of muscle were measured using the ELISA Kit (MLBio®, Shanghai, China) following their guidelines (38, 39). The substance extracted by the kit was determined under the condition of a 450-nm microplate reader.
2.4. Liquid chromatography mass spectrometry (LC–MS) measurement of metabolites
Metabolome measurements were performed by Lianchuan Biotechnology Co., Ltd (Hangzhou, China). First, a muscle biopsy was carried out to retrieve approximately 100-mg specimens, which were then ground and placed in liquid nitrogen. Then, the metabolites in the sample were extracted by shaking after adding 120 μl of 50% methanol (LC–MS grade, Honeywell®, Charlotte, USA) and treated at −20°C overnight to allow the proteins precipitated in the sample. After centrifugation (Heraeus™ Fresco17, Thermo Fisher Scientific®, Waltham, USA), the supernatants were transferred into a 96-well plate. To confirm its stability and repeatability, 10 μl was extracted from each sample and prepared as a quality control (QC) sample.
Compounds were analyzed in positive and negative ion modes by the TripleTOF 5600+ high-resolution tandem mass spectrometer (SCIEX®, Warrington, UK). Chromatographic separation was performed using an ultraperformance liquid chromatography (UPLC) system (SCIEX®, UK) with default parameters. The ACQUITY UPLC T3 column (100 mm × 2.1 mm, 1.8 μm, Waters, UK) was used for reverse-phase separation. Solvent A (water, 0.1% formic acid) and solvent B (acetonitrile, 0.1% formic acid) were used as the mobile phase to separate metabolites. A high-resolution mass spectrometer, TripleTOF 5600+ (SCIEX, UK) time-of-flight mass spectrometry was selected for the detection of eluted metabolites in the column. The ion spray floating voltage is set to 5 kV in positive ion mode and −4.5 kV in negative ion mode. The mass spectroscopy (MS) data were collected using Information Dependent Acquisition (IDA) mode. The mass range of TOF is 60–1,200 Da and survey scans are acquired every 150 ms. The first 12 signal ions with a positive charge (or a negative charge in the negative ion mode) and signal accumulation intensity of more than 100 times/s were selected from the MS1 map for MS2 fragmentation scanning. At the same time, every 10 QC samples were scanned. The quality difference between QC was used to correct the systematic error of the entire batch experiments.
2.5. Metabolomics data proceeding
The raw data were first processed using XCMS (40), which took time–m/z pairs, retention time, and peak area into consideration. Then, they were extracted and identified by msConvert (41) of the Proteowizard (42), XCMS (40), CAMERA (43), and metaX (44) packages that were included in R (v.4.1.2). The Total Ion Chromatogram (TIC) diagram utilizes retention time as the abscissa and the total intensity of each ion at each time point in the mass spectrum as the ordinate. The TIC peak retention time and peak area of the QC samples in this study overlapped well, indicating that the instrument was stable and reliable (Supplementary Figure 2). The m/z values of metabolites mostly range from 200 to 800 and the retention time mostly ranges from 0 to 4 min. The main workflow design of the study is shown in Figure 1.
For the identification of metabolites, we relied on three criteria: retention time (RT), accurate mass match to the library ± 10 ppm (MS1), MS/MS spectra (MS2) with high forward and reverse scores, as reported in previous studies (45–47). MS1 compounds were identified using the database of The Human Metabolome Database (HMDB) (www.hmdb.ca) (48) and Kyoto Encyclopedia of Genes and Genomes (KEGG) (https://www.kegg.jp/) (49), while MS2 metabolites were further identified using multiple databases (HMDB, KEGG; Massbank (http://www.massbank.jp/) (50); LipidBlast (51); and in-house database) as shown in Figure 1. Among them, the in-house database (self-built database) contains MS2 spectra of about 20,000 kinds of standard products (positive ions: 66,124; negative ions: 20,088). Due to certain differences in the species and quantities of metabolites between the two modes (Supplementary Figure 3), taken MS1 and MS2 together can alleviate the false positive rate to a large extent. We annotated these metabolites using CAMERA (43) and metaX (44). Specific parameters of the software are provided in Supplementary Tables 2–4.
2.6. Quality control analysis (QC)
Features with missing rate—more than 50% in QC samples or 80% in test samples—were removed, and the values for missing peaks were extrapolated with the k-nearest neighbor algorithm to further improve the data quality. In addition, the relative standard deviations of the metabolic features were calculated across all QC samples, and those with standard deviations >30% were removed. Data normalization was performed on all samples using the probabilistic quotient normalization algorithm. Then, QC-robust spline batch correction was performed using QC samples. Principal component analysis (PCA) was performed to detect outliers and batch effects using the pre-processed dataset.
2.7. Statistical analysis
Student's t-tests were used to identify metabolites that differed significantly between GZ and NQ horses. Two estimates of the false discovery rate (FDR, q-value) and the Benjamin–Hochberg procedure were calculated to correct the p-values. We also conducted the supervised Partial least squares discriminant analysis (PLS-DA) using metaX (44) for variables that discriminated the profiling statistical method. To verify the PLA-DA results, the displacement tests of R2 and Q2 were carried out for 200 times. R2 is the goodness-of-fit parameter and Q2 is the predictive ability parameter, which are used to verify the accuracy and predictability of the PLS-DA model. Multivariate analysis and PLS-DA discriminant analysis were used to obtain the variable important in projection (VIP) value of each metabolite variable. Differentially expressed metabolites were screened using VIP values of the first two principal components of the multivariate PLS-DA model, combined with p-values and fold changes, following the criteria, namely, fold change of ≥ 2, VIP of ≥ 1, and q-value of ≤ 0.05. The result was showed by a volcano plot using the R (v.4.1.2) packages of “ggplot2”.
A mental test was performed to establish connections between phenotype indexes and metabolites using the R (v.4.1.2) packages of “ggcor.”. The Spearman correlation coefficient was used to analyze the correlations between metabolome and functional enzymes with a threshold of p of < 0.05. Heatmaps were used to visualize different patterns of those differential metabolites between GZ and NQ horses. The distance was calculated as the Euclidean distance, with the clustering method of “ward D2”. Similarly, a scatter bubble diagram was introduced to show the metabolic pathways from the result of MetPA and pathway analysis module of the MetaboAnalyst 5.0 online platform (52).
3. Results
3.1. The variability muscle glycogen concentration in GZ and NQ horses
We found that the glycogen concentration was significantly high in the muscles of GZ horses (p = 4.8e–6) (Figure 2C) and that the citrate synthase (p = 0.02) and hexokinase (p = 0.006) activity were also high in GZ horses (Figures 2A, B). The citrate synthase and hexokinase activity are the core enzymes of the TCA cycle. Mental analysis showed a significant correlation between the identified metabolites and the muscle glycogen content (Supplementary Figure 4).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1162953/fvets-10-1162953-HTML/image_m/fvets-10-1162953-g002.jpg)
Figure 2. Muscle glycogen content and TCA oxidase activity of GZ and NQ horses. (A) The variability of citrate synthase activity in GZ and NQ horses. (B) The variability hexokinase activity in GZ and NQ horses. (C) The variability muscle glycogen concentration in GZ and NQ horses. *: p ≤ 0.05; **: p ≤ 0.01; ***: p ≤ 0.001; ****: p ≤ 0.0001.
3.2. Metabolite detection and classification using MS1 and MS2
A total of 51,535 metabolites were annotated by the HMDB database, which included 30,539 MS1 metabolites in positive mode, and included 20,996 in negative mode (Supplementary Table 5). All these metabolites can be classified into 24 superclass categories, and the first six groups accounted for more than 95% of all these metabolites; they were listed as lipids and lipid-like molecules (55.79%); phenylpropanoids and polyketides (14.28%); organic acids and derivatives (7.08%); organoheterocyclic compounds (6.58%); organic oxygenated compounds (5.86%); and benzenoids (5.72%). It is worth mentioning that these six categories accounted for more than 95% of all these metabolites (Figure 3A and Supplementary Figure 5). A total of 15,978 out of 51,535 MS1 metabolites were annotated into 40 pathways and 6 categories according to the KEGG database. Similarly, metabolites were most abundant in the amino acids and lipid metabolism pathways (Supplementary Figure 6 and Supplementary Table 5).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1162953/fvets-10-1162953-HTML/image_m/fvets-10-1162953-g003.jpg)
Figure 3. (A) MS1 metabolite classification bar chart. (B) MS2 metabolite classification bar chart. The abscissa represents the species of MS2 metabolites annotated to the class. Different colors in the bar chart represent different superclass categories. (C) A bubble map of MS2 metabolite pathway. The bubble color represents the size of p-value, and the smaller the p-value, the redder the color, the larger the p-value, and the bluer the color. Bubble size represents the amount of MS2 metabolites enriched. (D) PCA plot. The blue circle represents the Guanzhong (GZ) horse muscle individual, the yellow triangle represents the Ningqiang pony (NQ) horse muscle individual, and the red square represents the quality control (QC) individual. Each scatter represents a sample.
Additionally, these different substances were fragmented and were used to generate the MS2 mass spectrogram, according to the standard substance from the database. Then, they were scored to yield the matching result that ultimately provides the MS2 identification results of the metabolites (Supplementary Table 4). As a result, a total of 374 metabolites defined by MS2 were detected in positive mode and 167 in negative mode (Supplementary Table 5). Interestingly, nearly 44% of metabolites were clustered into lipids and lipid-like molecules in the superclass classification of these MS2 metabolites. Further annotation showed that they were clustered into glycerophospholipids (~20%), organic acids and derivatives (~18%), fatty acyls (~17%), and carboxylic acids and derivatives (~14%), as well as unknown (31%), and the detailed list is shown in Figure 3B and Supplementary Tables 6, 7. When these MS2 metabolites were enriched by the KEGG database, they were annotated into 47 pathways and 6 major groups (Figure 3C). Most of them were enriched in the metabolism pathway, which is similar to the result identified by MS1.
Both metabolites identified by MS1 and MS2 models were used for PCA. The results showed that GZ and NQ horses could be differentiated clearly and that PCA1 and PCA2 could explain 37.9% of metabolite variation between the two populations (Figure 3D). Reproducible aggregation of QC samples indicates good stability of the instrument and high quality of the measurement data.
3.3. Different muscle metabolites identified between GZ and NQ horses
The PLS-DA score plot showed that GZ and NQ horses could be separated at PC1 with a variation of ~20%, which mirrors the result of PCA (Figure 4A). In the permutation test plot of the PLS-DA model, the results showed that R2 and Q2 both decreased gradually with the gradual decrease in replacement retention, proving appropriate robustness of the model, and most of the Q2 values <0 indicated that the analysis of differential metabolites was accurate (Figure 4B). Of the MS1 metabolites identified in the muscle biopsy samples, only 149 were significantly altered between GZ and NQ horses (60 upregulation; 89 downregulation; Q < 0.05; FC > 2, VIP ≥ 1) (Figure 4C).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1162953/fvets-10-1162953-HTML/image_m/fvets-10-1162953-g004.jpg)
Figure 4. (A) PLS-DA score plot. Red represents the Guanzhong horse muscle group (GZ), blue represents the Ningqiang pony horse muscle group (NQ). (B) PLS-DA permutation test diagram showing the plots of R2, Q2, and accuracy. (C) A volcano plot of MS1 differential metabolites. Red represents upregulated metabolites and blue represents downregulated metabolites; Q = 0.05 was used as the threshold of whether the metabolites were significantly different. (D) A volcano plot of MS2 differential metabolites. (E) Cluster analysis of 13 differential. Six samples were Guanzhong (GZ) in blue group, and six samples were Ningqiang (NQ) in red group. (F) MetPA analysis of 13 differential MS2 metabolites pathway. The larger the bubble and abscissa values, the more important the path. The red color indicates the high impact of the respective metabolic pathway. p < 0.05 (significant).
Furthermore, by using the same strategy, 13 significant metabolites were identified, including 3 upregulated and 10 downregulated metabolites (Figure 4D). We found that there are 13 common differential metabolites from both MS1 and MS2 (Supplementary Table 8). Moreover, these 13 MS2 metabolites were classified into three superclass and five class categories (Supplementary Table 10).
3.4. Glutathione metabolism and taurine metabolism alter between GZ and NQ horses
With these abovementioned metabolites, we then made a cluster analysis to explore functional components indicative of muscle development (Figure 4E). In this way, a total of five metabolic pathways were detected, including glutathione metabolism (p = 0.01) and taurine and hypotaurine metabolism pathways (p = 0.04) (Figure 4F and Supplementary Table 9), suggesting that antioxidant supplementation responds to oxidative stress during muscle development.
3.5. Muscle metabolic network inferred for GZ and NQ horses
For a clear understanding of the interactive role of metabolites that contribute to muscle development, we demonstrate a metabolic network diagram to visualize the result in Supplementary Figure 7. As shown in Figures 5A, B, these metabolites are not independent but act as hub roles in the metabolite networks. For instance, the glutathione (GSH) metabolism pathway regulates taurine (TU), the hypotaurine (HTU) metabolism pathways regulates antioxidative mechanisms, and bile acid biosynthesis regulates GSH generation to modulate the oxidative response. The valine, leucine, and isoleucine degradation pathways are involved in the tricarboxylic acid cycle (TCA). Among them, 3-methyl-2-oxopentanoic acid can regulate the TCA cycle by producing tiglyl-CoA (Figure 5A). As intermediate metabolites, L-hydroxycarnitine and 3-hydroxycarnitine metabolites serve important roles in carnitine synthesis and oxidation of branched-chain fatty acid pathways (Figure 5B). These key metabolites were significantly more abundant in NQ horses than in GZ horses (Figures 5A, B).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1162953/fvets-10-1162953-HTML/image_m/fvets-10-1162953-g005.jpg)
Figure 5. (A, B) Thirteen significantly different metabolites pathway network graphs. Red color indicates the differential metabolites analyzed. The ordinate of the bar represents the relative abundance equivalent to the metabolite content. GZ, Guanzhong horse; NQ, Ningqiang pony horse. *: p ≤ 0.05; **: p ≤ 0.01; ***: p ≤ 0.001; ****: p ≤ 0.0001.
4. Discussion
With their large skeletal muscle mass and innate ability to exercise, horses are positioned as a unique model organism, and they can be used to understand the effects of breed differences, specifically on muscle metabolism. As one of the largest skeletal muscle of horses, the gluteal muscle has a function in stabilizing the limb during stance and locomotion (53). This study shows the presence of higher glycogen concentration in GZ horses than in NQ horses, which is the major source of energy during muscle development (54); the higher level of glycogen reflects stronger athletic performance in GZ horses. Besides, citrate synthase and hexokinase activity that played key roles in the TCA cycle (Figures 2A, B) were also associated with glycogen utilization to produce more energy for muscle growth and racing activities. These measurements may provide meaningful markers for the prediction of the strength of horse muscle locomotion ability. We then found metabolites in GZ and NQ horses, which were clustered into lipids and lipid-like molecules; these components are related to energy storage and production. These findings demonstrate that glycolytic and energy-related metabolisms have better performance in the racing breed horses to store and utilize more energy than short-stature horses (55).
There were 13 MS2 significantly differential metabolites detected in this study between the NQ and GZ horse breeds (Figure 4D); among them, L-glutathione (oxidized form), pyroglutamic acid, taurine, 3-dehydroxycarnitine, and L-propionylcarnitine metabolites were significantly associated with muscle movement, as reported in previous studies (56–60). Seven of these metabolites can be found in exercise training horses: taurine, L-propionylcarnitine, 3-dehydroxycarnitine (deoxycarnitine), 3-methyl-2-oxopentanoic acid (3-methyl-2-oxovalerate), Leu-Gly (L-leucylglycine), L-glutathione (oxidized form), and pyroglutamic acid (5-oxoproline) (58, 59, 61, 62); this indicates that some “housekeeping” ingredients functioned well during this process. Notably, three of them are related to the lipids and lipid-like molecules. Carnitine is essential for the transport of activated fatty acids across the mitochondrial membrane during mitochondrial β-oxidation (63). Muscles take up carnitine from blood in an exchange–diffusion process with endogenous deoxycarnitine, the immediate precursor of carnitine (64), and L-propionylcarnitine is an acylcarnitine. The general role of acylcarnitines is to transport acyl groups (organic acids and fatty acids) from the cytoplasm into the mitochondria so that they can be broken down to produce energy. This process is known as β-oxidation. L-Propionylcarnitine is a carnitine derivative that has a high affinity for muscular carnitine transferase: it increases cellular carnitine content, thereby allowing free fatty acid transport into the mitochondria. Long-term oral treatment with propionyl-L-carnitine improves maximum exercise duration and maximum oxygen consumption over placebo and indicates a specific propionyl-L-carnitine effect on peripheral muscle metabolism (65, 66). All of these indications suggest that lipid metabolism play an important role in muscle development and provides energy during this process (12). This may contribute to the differentiation of muscle development between GZ and NQ horses.
Further analysis demonstrated that glutathione metabolism (p = 0.01; Figure 4F) and the taurine and hypotaurine metabolism pathways (p = 0.04; Figure 4F) were significantly associated with muscle metabolic mechanisms. Moreover, these metabolites interact with each other and treated as rate-limiting metabolites in the TCA cycle, which generate energy for biological process. The glutathione metabolism pathway is of important metabolite for antioxygen. The L-glutathione (oxidized form) is a ubiquitous and powerful cellular antioxidant essential for protecting cells from the harmful effects of oxidative damage and free radical damage. It is a highly abundant metabolite in cells (67). Furthermore, the taurine pathway, which is widely distributed in animal tissues and has many biological functions, such as an important role in osmotic pressure regulation, act as a stabilizer for the cell membrane and protein, having antioxidant and anti-inflammatory functions, regulating mitochondrial tRNA activity, and participating in calcium homeostasis (68–70). Its deficiency leads to muscle dysfunction and affects athletic performance (70–72), and is related to skeletal muscle aging (73, 74). This is reasonable along with the notion that the muscle glycogen content is significantly the Pearson association with those metabolites (Supplementary Figure 4) because the more the energy resource, the more is the metabolites produced. These significant metabolites reveal differences in glycolytic metabolism and exercise capacity in horse muscles of different breeds. It provided a theoretical basis for the muscle metabolism analysis of different species and for the direction of identification and breeding.
5. Conclusion
In this study, the differential metabolites in muscle were detected between two different horse breeds (GZ and NQ) living in the same area. We found that the GZ athletic horses exhibit higher muscle glycogen contents and higher enzyme activities than NQ horses. While, notably, the lipid and lipid-like taurine and taurine hypotaurine metabolism and glutathione metabolism, which were interactive in the TCA cycle, were significantly decreased in GZ horses. This indicates that GZ horses possess higher energy storage and utilization ability to meet its requirement for exercise and racing. In summary, we conclude that these metabolites have the potential to be used as an indicator to understand the metabolic aspects of exercise capacity in different breeds of horses. In future, in vivo and in vitro validations will be needed to understand the biological mechanisms.
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.
Ethics statement
The animal study was reviewed and approved by Animal Care Committee of the Institute of Animal Sciences, Chinese Academy of Agricultural Sciences (approval number: IAS2019-24).
Author contributions
SM prepared samples, performed statistical analyses, and wrote the manuscript. YZ collected samples, prepared samples, and reviewed the manuscript. SL and ZZ collected samples. XL idealized the project and reviewed the manuscript. LJ and XL supervised this work. All authors contributed to and approved the final manuscript.
Funding
The project was supported by the National Natural Science Foundation of China (32002144, 31560620, 31972530, and 31772553). XL received the support of the International Postdoctoral Exchange Fellowship Program (20190102) and the European Union's Horizon 2020 research and innovation program under the Marie Sokolowskis–Curie grant agreement 101027750.
Acknowledgments
The authors acknowledge the financial support from the National Natural Science Foundation of China (31560620, 32002144, 31972530, and 31772553). The authors thank all owners of affected horses for their collaboration. The authors thank Lianchuan Biotechnology Co., Ltd. for its support in liquid chromatography–mass spectrometry analysis, peak extraction, and compound identification.
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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fvets.2023.1162953/full#supplementary-material
References
1. Klein DJ, Anthony TG, McKeever KH. Metabolomics in equine sport and exercise. J Anim Physiol Anim Nutr. (2021) 105:140–8. doi: 10.1111/jpn.13384
2. Goldansaz SA, Guo AC, Sajed T, Steele MA, Plastow GS, Wishart DS, et al. Livestock metabolomics and the livestock metabolome: a systematic review. PLoS One. (2017) 12:e177675. doi: 10.1371/journal.pone.0177675
3. Liu X, Zhang Y, Li Y, Pan J, Wang D, Chen W, et al. EPAS1 gain-of-function mutation contributes to high-altitude adaptation in Tibetan horses. Mol Biol Evol. (2019) 36:2591–03. doi: 10.1093/molbev/msz158
4. Wouters CP, Toquet MP, Renaud B, Francois AC, Fortier-Guillaume J, Marcillaud-Pitel C, et al. Metabolomic signatures discriminate horses with clinical signs of atypical myopathy from healthy co-grazing horses. J Proteome Res. (2021) 20:4681–92. doi: 10.1021/acs.jproteome.1c00225
5. Steelman SM, Johnson P, Jackson A, Schulze J, Chowdhary BP. Serum metabolomics identifies citrulline as a predictor of adverse outcomes in an equine model of gut-derived sepsis. Phys Genom. (2014) 46:339–47. doi: 10.1152/physiolgenomics.00180.2013
6. Luck MM, Le Moyec L, Barrey E, Triba MN, Bouchemal N, Savarin P, et al. Energetics of endurance exercise in young horses determined by nuclear magnetic resonance metabolomics. Front Physiol. (2015) 6:198. doi: 10.3389/fphys.2015.00198
7. 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 Genom. (2017) 18:187. doi: 10.1186/s12864-017-3571-3
8. de Meeûs d'Argenteuil C, Boshuizen B, Vidal Moreno de Vega C, Leybaert L, de Maré L, Goethals K, et al. (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
9. Cloteau C, Dervilly G, Kaabia Z, Bagilet F, Delcourt V, Loup B, et al. From a non-targeted metabolomics approach to a targeted biomarkers strategy to highlight testosterone abuse in equine. Illustration of a methodological transfer between platforms and laboratories. Drug Test Anal. (2022) 14:864–78. doi: 10.1002/dta.3221
10. Keen B, Cawley A, Reedy B, Fu S. Metabolomics in clinical and forensic toxicology, sports anti-doping and veterinary residues. Drug Test Anal. (2022) 14:794–807. doi: 10.1002/dta.3245
11. Stojiljkovic N, Leroux F, Bubanj S, Popot MA, Paris A, Tabet JC, et al. Tracking main environmental factors masking a minor steroidal doping effect using metabolomic analysis of horse urine by liquid chromatography-high-resolution mass spectrometry. Eur J Mass Spectrom. (2019) 25:339–53. doi: 10.1177/1469066719839034
12. Klein DJ, McKeever KH, Mirek ET, Anthony TG. Metabolomic response of equine skeletal muscle to acute fatiguing exercise and training. Front Physiol. (2020) 11:110. doi: 10.3389/fphys.2020.00110
13. Ohmura H, Mukai K, Takahashi Y, Takahashi T. Metabolomic analysis of skeletal muscle before and after strenuous exercise to fatigue. Sci Rep. (2021) 11:11261. doi: 10.1038/s41598-021-90834-y
14. Middleton RP, Lacroix S, Scott-Boyer MP, Dordevic N, Kennedy AD, Slusky AR, et al. Metabolic differences between dogs of different body sizes. J Nutr Metab. (2017) 2017:4535710. doi: 10.1155/2017/4535710
15. Brookes PS, Jimenez AG. Metabolomics of aging in primary fibroblasts from small and large breed dogs. GeroScience. (2021) 43:1683–96. doi: 10.1007/s11357-021-00388-0
16. Shi K, Zhao Q, Shao M, Duan Y, Li D, Lu Y, et al. Untargeted metabolomics reveals the effect of selective breeding on the quality of chicken meat. Metabolites. (2022) 12. doi: 10.3390/metabo12050367
17. Liao X, Shi X, Hu H, Han X, Jiang K, Liu Y, et al. Comparative metabolomics analysis reveals the unique nutritional characteristics of breed and feed on muscles in chinese taihe black-bone silky fowl. Metabolites. (2022) 12. doi: 10.3390/metabo12100914
18. Ritota M, Casciani L, Failla S, Valentini M. HRMAS-NMR spectroscopy and multivariate analysis meat characterisation. Meat Sci. (2012) 92:754–61. doi: 10.1016/j.meatsci.2012.06.034
19. Zhu C, Petracci M, Li C, Fiore E, Laghi L. An untargeted metabolomics investigation of Jiulong Yak (Bos grunniens) Meat by (1)H-NMR. Foods. (2020) 9:481. doi: 10.3390/foods9040481
20. Lisuzzo A, Bonelli F, Sgorbini M, Nocera I, Cento G, Mazzotta E, et al. Differences of the plasma total lipid fraction from pre-foaling to post-foaling period in donkeys. Animals. (2022) 12. doi: 10.3390/ani12030304
21. Lisuzzo A, Laghi L, Fiore F, Harvatine K, Mazzotta E, Faillace V, et al. Evaluation of the metabolomic profile through (1)H-NMR spectroscopy in ewes affected by postpartum hyperketonemia. Sci Rep. (2022) 12:16463. doi: 10.1038/s41598-022-20371-9
22. Mukai K, Kitaoka Y, Takahashi Y, Takahashi T, Takahashi K, Ohmura H, et al. Moderate-intensity training in hypoxia improves exercise performance and glycolytic capacity of skeletal muscle in horses. Physiol Rep. (2021) 9:e15145. doi: 10.14814/phy2.15145
23. He Y, Nadeau J, Reed S, Hoagland T, Bushmich S, Aborn S, et al. The effect of season on muscle growth, fat deposition, travel patterns, and hoof growth of domestic young horses. J Equine Vet Sci. (2020) 85:102817. doi: 10.1016/j.jevs.2019.102817
24. Raspa F, Dinardo FR, Vervuert I, Bergero D, Bottero MT, Pattono D, et al. A Fibre- vs. cereal grain-based diet: Which is better for horse welfare? Effects on intestinal permeability, muscle characteristics and oxidative status in horses reared for meat production. J Anim Physiol Anim Nutr. (2022) 106:313–26. doi: 10.1111/jpn.13643
25. Liu X, Zhang Y, Liu W, Li Y, Pan J, Pu Y, et al. A single-nucleotide mutation within the TBX3 enhancer increased body size in Chinese horses. Curr Biol. (2022) 32:480–7. doi: 10.1016/j.cub.2021.11.052
26. Zeng L, Chen N, Yao Y, Dang R, Chen H, Lei C, et al. Analysis of genetic diversity and structure of guanzhong horse using microsatellite markers. Anim Biotechnol. (2019) 30:95–8. doi: 10.1080/10495398.2017.1416392
27. Pu Y, Zhang Y, Zhang T, Han J, Ma Y, Liu X, et al. Identification of novel lncRNAs differentially expressed in placentas of Chinese Ningqiang Pony and Yili horse breeds. Animals. (2020) 10. doi: 10.3390/ani10010119
28. Zhang T, Lu H, Cao L, Zeng Z. Genetic diversity of microsatellite on ningqiang pony. Hubei Agricultural Sciences. (2008) 868–70.
29. National Livestock and Poultry Genetic Resources. (2011). Animal genetic resources in China-Horses Donkeys. Camels: China Agricultural Press.
30. Kearns CF, McKeever KH, Abe T. Overview of horse body composition and muscle architecture: implications for performance. Vet J. (2002) 164:224–34. doi: 10.1053/tvjl.2001.0702
31. Rivero JL, 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
32. Gika HG, Theodoridis GA, Plumb RS, Wilson ID. Current practice of liquid chromatography-mass spectrometry in metabolomics and metabonomics. J Pharm Biomed Anal. (2014) 87:12–25. doi: 10.1016/j.jpba.2013.06.032
33. Ribbenstedt A, Ziarrusta H, Benskin JP. Development, characterization and comparisons of targeted and non-targeted metabolomics methods. PLoS One. (2018) 13:e207082. doi: 10.1371/journal.pone.0207082
34. Otter D, Cao M, Lin HM, Fraser K, Edmunds S, Lane G, et al. Identification of urinary biomarkers of colon inflammation in IL10-/- mice using Short-Column LCMS metabolomics. J Biomed Biotechnol. (2011) 2011:974701. doi: 10.1155/2011/974701
35. Myers OD, Sumner SJ, Li S, Barnes S, Du X. Detailed investigation and comparison of the XCMS and MZmine 2 chromatogram construction and chromatographic peak detection methods for preprocessing mass spectrometry metabolomics data. Anal Chem. (2017) 89:8689–95. doi: 10.1021/acs.analchem.7b01069
36. Wu IL, Turnipseed SB, Storey JM, Andersen WC, Madson MR. Comparison of data acquisition modes with Orbitrap high-resolution mass spectrometry for targeted and non-targeted residue screening in aquacultured eel. Rapid Commun Mass Spectrom. (2020) 34:e8642. doi: 10.1002/rcm.8642
37. Thies C. Insight into traditional Chinese medicine obtained from an October 2004 visit to china organized by the people to people ambassador program. Pharm Dev Technol. (2005) 10:447–50. doi: 10.1080/10837450500411892
38. Bergh A, Nordlof H, Essen-Gustavsson B. Evaluation of neuromuscular electrical stimulation on fibre characteristics and oxidative capacity in equine skeletal muscles. Equine Vet J Suppl. (2010) 671–5. doi: 10.1111/j.2042-3306.2010.00180.x
39. Kitaoka Y, Mukai K, Aida H, Hiraga A, Masuda H, Takemasa T, et al. Effects of high-intensity training on lipid metabolism in Thoroughbreds. Am J Vet Res. (2012) 73:1813–8. doi: 10.2460/ajvr.73.11.1813
40. Smith CA, Want EJ, O'Maille G, Abagyan R, Siuzdak G. XCMS: processing mass spectrometry data for metabolite profiling using nonlinear peak alignment, matching, and identification. Anal Chem. (2006) 78:779–87. doi: 10.1021/ac051437y
41. Adusumilli R, Mallick P. Data conversion with proteowizard msconvert. Methods Mol Biol. (2017) 1550:339–68. doi: 10.1007/978-1-4939-6747-6_23
42. Kessner D, Chambers M, Burke R, Agus D, Mallick P. ProteoWizard: open source software for rapid proteomics tools development. Bioinformatics. (2008) 24:2534–6. doi: 10.1093/bioinformatics/btn323
43. Kuhl C, Tautenhahn R, Bottcher C, Larson TR, Neumann S. CAMERA: an integrated strategy for compound spectra extraction and annotation of liquid chromatography/mass spectrometry data sets. Anal Chem. (2012) 84:283–9. doi: 10.1021/ac202450g
44. Wen B, Mei Z, Zeng C, Liu S. metaX: a flexible and comprehensive software for processing metabolomics data. BMC Bioinformatics. (2017) 18:183. doi: 10.1186/s12859-017-1579-y
45. Shen B, Yi X, Sun Y, Bi X, Du J, Zhang C, et al. Proteomic and metabolomic characterization of COVID-19 Patient Sera. Cell. (2020) 182:59–72. doi: 10.1016/j.cell.2020.05.032
46. Coletto E, Latousakis D, Pontifex MG, Crost EH, Vaux L, Perez SE, et al. The role of the mucin-glycan foraging Ruminococcus gnavus in the communication between the gut and the brain. Gut Microbes. (2022) 14:2073784. doi: 10.1080/19490976.2022.2073784
47. Liang L, Rasmussen MH, Piening B, Shen X, Chen S, Rost H, et al. Metabolic dynamics and prediction of gestational age and time to delivery in pregnant women. Cell. (2020) 181:1680–92. doi: 10.1016/j.cell.2020.05.002
48. Djoumbou FY, Eisner R, Knox C, Chepelev L, Hastings J, Owen G, et al. ClassyFire: automated chemical classification with a comprehensive, computable taxonomy. J Cheminformatics. (2016) 8:61. doi: 10.1186/s13321-016-0174-y
49. Aoki-Kinoshita KF, Kanehisa M. Gene annotation and pathway mapping in KEGG. Methods Mol Biol. (2007) 396:71–91. doi: 10.1007/978-1-59745-515-2_6
50. Horai H, Arita M, Kanaya S, Nihei Y, Ikeda T, Suwa K, et al. MassBank: a public repository for sharing mass spectral data for life sciences. J Mass Spectrom. (2010) 45:703–14. doi: 10.1002/jms.1777
51. Kind T, Liu KH, Lee DY, DeFelice B, Meissen JK, Fiehn O, et al. LipidBlast in silico tandem mass spectrometry database for lipid identification. Nat Methods. (2013) 10:755–8. doi: 10.1038/nmeth.2551
52. Pang, Z., Chong J., Zhou G., de Lima Morais D. A., Chang L., Barrette M., et al. (2021). MetaboAnalyst 5.0: narrowing the gap between raw spectra and functional insights. Nucleic Acids Res. 49, W388–W396. doi: 10.1093/nar/gkab382
53. Zsoldos RR, Voegele A, Krueger B, Schroeder U, Weber A, Licka TF, et al. Long term consistency and location specificity of equine gluteus medius muscle activity during locomotion on the treadmill. BMC Vet Res. (2018) 14:126. doi: 10.1186/s12917-018-1443-y
54. Vigh-Larsen JF, Ortenblad N, Spriet LL, Overgaard K, Mohr M. Muscle glycogen metabolism and high-intensity exercise performance: a narrative review. Sports Med. (2021) 51:1855–74. doi: 10.1007/s40279-021-01475-0
55. Valberg SJ, Macleay JM, Billstrom JA, Hower-Moritz MA, Mickelson JR. Skeletal muscle metabolic response to exercise in horses with 'tying-up' due to polysaccharide storage myopathy. Equine Vet J. (1999) 31:43–7. doi: 10.1111/j.2042-3306.1999.tb03789.x
56. Vaz FM, Wanders RJ. Carnitine biosynthesis in mammals. Biochem J. (2002) 361:417–29. doi: 10.1042/bj3610417
57. Brevetti G, Fanin M, De Amicis V, Carrozzo R, Di Lello F, Martone VD, et al. Changes in skeletal muscle histology and metabolism in patients undergoing exercise deconditioning: effect of propionyl-L-carnitine. Muscle Nerve. (1997) 20:1115–20.
58. Ma Y, Maruta H, Sun B, Wang C, Isono C, Yamashita H, et al. Effects of long-term taurine supplementation on age-related changes in skeletal muscle function of Sprague-Dawley rats. Amino Acids. (2021) 53:159–70. doi: 10.1007/s00726-020-02934-0
59. Seidel U, Huebbe P, Rimbach G. Taurine: a regulator of cellular redox homeostasis and skeletal muscle function. Mol Nutr Food Res. (2019) 63:e1800569. doi: 10.1002/mnfr.201800569
60. Vargiu R, Licheri D, Carcassi AM, Naimi S, Collu M, Littarru GP, et al. Enhancement of muscular performance by a coformulation of propionyl-L-carnitine, coenzyme Q10, nicotinamide, riboflavin and pantothenic acid in the rat. Physiol Behav. (2002) 76:257–63. doi: 10.1016/S0031-9384(02)00717-5
61. Carbonin PU, Ramacci MT, Pahor M, Di Gennaro M, Gambassi Jr G, Giudice L, et al. Antiarrhythmic effect of L-propionylcarnitine in isolated cardiac preparations. Cardioscience. (1991) 2:109–14.
62. Fischer M, Keller J, Hirche F, Kluge H, Ringseis R, Eder K, et al. Activities of gamma-butyrobetaine dioxygenase and concentrations of carnitine in tissues of pigs. Comp Biochem Physiol A-Mol Integr Physiol. (2009) 153:324–31. doi: 10.1016/j.cbpa.2009.03.005
63. Hoppel C. The role of carnitine in normal and altered fatty acid metabolism. Am J Kidney Dis. (2003) 41:S4–S12. doi: 10.1016/S0272-6386(03)00112-4
64. Siliprandi N, Di Lisa F, Menabo R, Ciman M, Sartorelli L. Transport and functions of carnitine in muscles. J Clin Chem Clin Biochem. (1990) 28:303–6.
65. Ferrari R, Merli E, Cicchitelli G, Mele D, Fucili A, Ceconi C, et al. Therapeutic effects of L-carnitine and propionyl-L-carnitine on cardiovascular diseases: a review. Ann NY Acad Sci. (2004) 1033:79–91. doi: 10.1196/annals.1320.007
66. Brevetti G, Perna S, Sabba C, Martone VD, Condorelli M. Propionyl-L-carnitine in intermittent claudication: double-blind, placebo-controlled, dose titration, multicenter study. J Am Coll Cardiol. (1995) 26:1411–6. doi: 10.1016/0735-1097(95)00344-4
67. Tomin T, Schittmayer M, Birner-Gruenberger R. Addressing glutathione redox status in clinical samples by two-step alkylation with n-ethylmaleimide isotopologues. Metabolites. (2020) 10. doi: 10.3390/metabo10020071
68. Wen C, Li F, Zhang L, Duan Y, Guo Q, Wang W, et al. Taurine is involved in energy metabolism in muscles, adipose tissue, and the liver. Mol Nutr Food Res. (2019) 63:e1800536. doi: 10.1002/mnfr.201800536
69. Veeravalli S, Phillips IR, Freire RT, Varshavi D, Everett JR, Shephard EA, et al. Flavin-containing monooxygenase 1 catalyzes the production of taurine from hypotaurine. Drug Metab Dispos. (2020) 48:378–85. doi: 10.1124/dmd.119.089995
70. Warskulat U, Flogel U, Jacoby C, Hartwig HG, Thewissen M, Merx MW, et al. Taurine transporter knockout depletes muscle taurine levels and results in severe skeletal muscle impairment but leaves cardiac function uncompromised. FASEB J. (2004) 18:577–9. doi: 10.1096/fj.03-0496fje
71. Hamilton EJ, Berg HM, Easton CJ, Bakker AJ. The effect of taurine depletion on the contractile properties and fatigue in fast-twitch skeletal muscle of the mouse. Amino Acids. (2006) 31:273–8. doi: 10.1007/s00726-006-0291-4
72. Ito T, Oishi S, Takai M, Kimura Y, Uozumi Y, Fujio Y, et al. Cardiac and skeletal muscle abnormality in taurine transporter-knockout mice. J Biomed Sci. (2010) 17 Suppl 1:S20. doi: 10.1186/1423-0127-17-S1-S20
73. Ito T, Yoshikawa N, Inui T, Miyazaki N, Schaffer SW, Azuma J, et al. Tissue depletion of taurine accelerates skeletal muscle senescence and leads to early death in mice. PLoS One. (2014) 9:e107409. doi: 10.1371/journal.pone.0107409
Keywords: horse (Equus caballus), untargeted metabolomics, muscle metabolism, lipid, glutathione metabolism
Citation: Meng S, Zhang Y, Lv S, Zhang Z, Liu X and Jiang L (2023) Comparison of muscle metabolomics between two Chinese horse breeds. Front. Vet. Sci. 10:1162953. doi: 10.3389/fvets.2023.1162953
Received: 10 February 2023; Accepted: 03 April 2023;
Published: 05 May 2023.
Edited by:
Huansheng Yang, Hunan Normal University, ChinaCopyright © 2023 Meng, Zhang, Lv, Zhang, Liu and Jiang. 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: Xuexue Liu, lxx_caas@163.com; Lin Jiang, jianglin@caas.cn
†These authors have contributed equally to this work and share first authorship