Plasma metabolomic profiles as affected by diet and stress in Spanish goats

The effects of high-condensed tannin (CT) diet combined with preslaughter stress have not been studied at the metabolome level in goats. This study was conducted to determine the effects of feeding sericea lespedeza (SL; Lespedeza cuneata), a high-CT legume, and transportation stress on plasma metabolome in goats. Uncastrated male Spanish goats (age = 8 months; BW = 26.0 ± 0.48 kg) were either fed ground ‘Serala’ SL hay (SER), bermudagrass (Cynodon dactylon) hay (BG), or bermudagrass hay-dewormed goats (BG-DW; Control) at 75% of intake, with a corn-based supplementation (25%) for 8 weeks (n = 12/Diet). At the end of the trial, goats were subjected to one of two stress treatments (ST): transported for 90 min to impose stress (TS) or held in pens (NTS) before slaughtering, in two replicates. Live and carcass weights, and blood samples were collected at 0, 30, 60 and 90 min of transportation or holding time (Time). The data were analyzed using MIXED procedures in SAS and metabolomics data were analyzed using the R software. When measured after ST, SER group had the lowest body weight (P < 0.05) among the three diet groups. Carcass weights were high in the BG-DW, low in SER, and intermediate in BG group. Plasma creatine concentrations decreased over Time (P < 0.01) in the TS goats in all diet groups. Meat crude protein percentages were higher (P < 0.05) in SER (22.5 ± 0.22) and BG-DW (22.3 ± 0.22) groups compared to the BG group (21.6 ± 0.22). At the metabolome level, SER group had the lowest (P < 0.05) glycine, alanine, threonine, taurine, trans-hydroxyproline, methionine, and histidine concentrations and highest (P < 0.01) lysine and citrulline concentrations among the Diet groups. Butyric acid, concentration was higher (P < 0.05) in the SER group compared to BG group. Eight medium- and long-chained acylcarnitines were higher (P < 0.05) in the BG-DW group than SER or BG groups. In general, amino acid levels decreased and acylcarnitine increased with Time (P < 0.05) in all groups. Sericea diet can be beneficial in enhancing stress coping abilities in goats due to elevated butyrate, lysine, and citrulline levels; however, SER resulted in lower energy level in goats compared to BG or BG-DW groups. Fatty acid metabolism is the main energy pathway in all groups during prolonged stress. Inclusion of certain varieties of SL in the diet must be carefully controlled to prevent possible negative effect.


Scientific Reports
| (2021) 11:12607 | https://doi.org/10.1038/s41598-021-91893-x www.nature.com/scientificreports/ Metabolomics has emerged as a robust high-output analytical method for identifying metabolites in ruminants in response to stress and diet 11,12 . It is the large-scale study of low molecular weight molecules within cells, biofluids, tissues, or organisms, and these metabolites can be influenced by genetic and environmental factors 13 . Metabolomics can be used to identify the metabolite biomarkers and biological pathways of importance in both basic and applied research 14 . Typically, compounds are identified by comparing the measured masses to those of known metabolites stored in databases such as HMDB 15 , LipidMaps 16 and Metlin 17 .
Studies on the effects of diet and stress on plasma metabolomics in goats are very limited. Although it has been reported that SL diet can have several positive effects in goats, to the best of our knowledge, there has been no study on the combined effects of SL and preslaughter stress at the metabolome level in goats. The overall aim of the present study was to determine the separate and combined effects of feeding a diet rich in CT and shortterm simulated preslaughter stress on live and carcass weights, meat composition, and plasma metabolome in Spanish goats.

Methods
Animals. The animal care protocol for this research was reviewed and approved by the Fort Valley State University's Agricultural and Laboratory Animal Care and Use Committee (ALACUC) prior to beginning of the experiment. In a Completely Randomized Design, 36 uncastrated male Spanish goats (body weight 26.1 ± 0.16 kg) were assigned to 3 treatments: ground 'Serela' SL hay (SER), bermudagrass hay (BG), or bermudagrass hay-dewormed goats (BG-DW; Control) at 75% of intake, with a corn-based supplementation (25%) for 8 weeks (n = 12 goats/group). Goats were kept in individual pens during the feeding trial and feed intake and average daily gain were monitored. At the end of the feeding period, all animals were weighed and transported to the processing facility and held overnight in holding pens prior to slaughter. On the day of slaughter, half the goats were transported (TS) for 90 min to impose stress (stress treatment, ST) and the remaining goats were held in pens (non-transported, NTS) before slaughtering. The trial was conducted on two different days (replicate) using 18 goats on each day. All goats were weighed prior to beginning of transportation and after transportation. Animals were slaughtered using humane procedures.
Feed sampling. Samples of hay and concentrates collected at the beginning, middle, and end of the experiment were pooled and analyzed for nutrient content at a commercial laboratory (Dairy One Forage Laboratory, Ithaca, New York) and estimates of diet nutrient analyses were calculated. These values are presented along with the ingredients used to formulate concentrate portions in Table 1. The 'Serala' sericea lespedeza hay samples analyzed had a total CT content of 5.1% and the bermudagrass hay samples had virtually no CT content.
Blood and meat sampling. Blood samples were collected every 30 min from individual animals by jugular venipuncture into the vacuum tubes containing K 2 EDTA and kept on ice until separation of plasma. Blood tubes were centrifuged at 3000 rpm for 10 min to separate the plasma. Plasma aliquots were stored at − 20 °C until analysis. Goats were slaughtered and the carcasses were weighed and stored at 2 °C. After 24 h of cooler storage, each carcass was then weighed again, and loin/rib chops were collected, vacuum packed, and stored at − 20 °C for proximate analysis.
Proximate analysis. Percent dry matter, ash, fat, and protein were analyzed using AOAC 18 procedures.
Samples for proximate analysis were thawed at 4 °C. Total protein was determined using a carbon/nitrogen analyzer (Vario MAX cube, Elementar, NJ), standardized with l-glutamic acid. Approximately 0.5 g of meat samples were taken in metal crucibles and placed in the auto sampler. Percent nitrogen was determined by combustion (900 °C), and these values were multiplied by 6.25 to obtain protein percentages. For dry matter analysis, approximately 2 g samples were placed in crucibles and kept in an oven at 100 °C overnight. The crucibles were taken out, placed in a desiccator to cool, and then weighed to calculate the percent dry matter. The crucible with dried sample from the dry matter determination was placed in a furnace (Thermolyne, Type 48000 Furnace, Augusta, Georgia) at 600 °C for at least 2 h. The crucibles with ashed samples were taken out, placed in the desiccator to cool, and then weighed to calculate percent ash. Ether extract was determined using a Soxtec 8000 (Foss North America, Eden Prairie, MN) Extraction Unit.
Metabolomics. An aliquot (200 µL) of plasma samples from SER, BG, and BG-DW groups collected at 0, 30, 60 and 90 min of transportation from TS and NTS animals were shipped on dry ice to the University of Alberta (Edmonton, AB Canada) for targeted metabolomics analysis. A total of 136 metabolites were analyzed using DI/LC-MS/MS TMIC PRIME Assay.
A targeted quantitative metabolomics approach was employed to analyze the samples using a combination of direct injection mass spectrometry with a reverse-phase LC-MS/MS custom assay. This custom assay, in combination with an ABSciex 4000 Qtrap (Applied Biosystems/MDS Sciex) mass spectrometer, can be used for the targeted identification and quantification of up to 150 different endogenous metabolites, including amino acids, acylcarnitines, biogenic amines and derivatives, uremic toxins, glycerophospholipids, sphingolipids, and sugars 19 . The method combines the derivatization and extraction of analytes, and the selective mass-spectrometric detection using multiple reaction monitoring (MRM) pairs. Isotope-labeled internal standards and other internal standards were used for metabolite quantification. The custom assay contains a 96 deep-well plate with a filter plate attached with sealing tape, and reagents and solvents used to prepare the plate assay. The first 14 wells were used for one blank, three zero samples, seven standards and three quality control samples. For all metabolites except organic acids, samples were thawed on ice and were vortexed and centrifuged at 13,000×g. Each sample (10 µL) was loaded onto the center of the filter on the upper 96-well plate and dried in a stream of nitrogen. www.nature.com/scientificreports/ Subsequently, phenyl-isothiocyanate was added for derivatization. After incubation, the filter spots were dried again using an evaporator. Extraction of the metabolites was then achieved by adding 300 µL of extraction solvent. The extracts were obtained by centrifugation into the lower 96-deep well plate, followed by a dilution step with MS running solvent. For organic acid analysis, 150 µL of ice-cold methanol and 10 µL of isotope-labeled internal standard mixture was added to 50 µL of sample for overnight protein precipitation, followed by centrifugation at 13,000×g for 20 min. The supernatant (50 µL) was loaded into the center of wells of a 96-deep well plate, followed by the addition of 3-nitrophenylhydrazine (NPH) reagent. After incubation for 2 h, BHT stabilizer and water were added before LC-MS injection.
Mass spectrometric analysis was performed on an ABSciex 4000 Qtrap ® tandem mass spectrometry instrument (Applied Biosystems/MDS Analytical Technologies, Foster City, CA) equipped with an Agilent 1260 series UHPLC system (Agilent Technologies, Palo Alto, CA). The samples were delivered to the mass spectrometer by an LC method followed by a direct injection method. Data analysis was done using Analyst 1.6.2.
Statistical analysis. The body weight, carcass weight, and meat proximate concentration data were analyzed using Mixed Procedures in SAS (release 9.1, SAS Institute, Cary, NC, USA) as a Completely Randomized  20 to each metabolite in order to create intuitive figures for visualizing the data. Principal Component Analysis (PCA), an unsupervised clustering technique, was used to reduce dimensionality of data and examine the intrinsic variation in the data set. To analyze selected metabolite groups of interest (amino acids/organic acids, and acyl carnitines), a classical analytical approach to assess group-wise differences in a univariate parameter-by-parameter fashion (t test) was used. The Student's t test (α = 0.05) was utilized to analyze each metabolite between two experimental groups: BG vs SER and BG-DW vs SER. Univariate volcano plots were generated by plotting the fold change against false discovery rate adjusted P-value. Volcano plots show data based on their t test P-values vs. their fold change and are generally used to rapidly identify and visually display significant differences in metabolites.
To perform a more robust analysis, metabolomics data within TS and NTS groups were analyzed using the R software using a Repeated Measures Model. Fixed effects included diet and time, and animal was added as random effect. Interaction effects were also included in the model. To evaluate changes in plasma metabolome caused by Diet (SER vs. BG) effect with an ST group, a hierarchical clustering analysis using the degree of abundance was conducted in which metabolites with identical abundance patterns were clustered together. Clustering was visualized by means of heatmap dendograms. Proc Mixed in SAS provides a very flexible environment to model many types of repeated measures data, including repeated in time. The metabolites that were found significant were further analyzed using Mixed procedure in SAS to compare the main effects of Diet (SER, BG, BG-DW), Time (0, 30, 60, and 90 min), and their interaction, and the means were separated using the LSD test when significant by ANOVA at P < 0.05.
Ethics approval and consent to participate. All methods in this study were carried out in accordance with relevant guidelines and regulations. The animal care protocols were approved by the FVSU Agricultural and Laboratory Animal Care and Use Committee prior to beginning of the study following the ADSA-ASAS-PSA Guide for Care and Use of Agricultural Animals in Research and Teaching.

Results
Feed intake and average daily gain. Diet did not influence feed intake; however, average daily gain was affected by Diet (P < 0.01). The average daily gains were 57.0 ± 4.72, 94.1 ± 5.12, and 127.3 ± 4.94 g, respectively, in SER, BG, and BG-DW groups.
Live and carcass weights and carcass yield. Diet main effect was significant for the goats' live weight when measured after ST (P < 0.05; Table 2) with SER group having the lowest body weight among the three dietary treatment groups. Interaction means revealed that the goats subjected to BG-DW plus NTS had the highest body weight and those subjected to SER plus TS had the lowest body weight. Hot and cold carcass weights were influenced by Diet (P < 0.01; Table 2) but not by ST or Diet × ST interaction. Both hot and cold carcass weights  (Table 3).
Plasma glucose, creatine, and creatinine concentrations. Plasma glucose concentrations were lower (P < 0.05) in SER diet group than in BG-DW or BG groups. There were significant Time (P < 0.01) and ST × Time interaction (P < 0.01) effects on glucose concentrations, since the levels in the TS group increased steadily over transportation time, while such a pattern was absent in the NTS group (Fig. 1A). Plasma creatine levels were not influenced by any interaction effects studied other than ST × Time. Overall mean creatine concentrations decreased over transportation time (Time, P < 0.01) in the TS group but not in the NTS group (Fig. 1B). Plasma creatinine levels were not influenced by any of the factors studied other than Diet effect. Plasma creatinine concentrations were high in SER diet group, low in BG, and intermediate in the BG-DW group (P < 0.01; Fig. 1C).

Plasma metabolomics.
After conducting a preliminary PCA to study the clustering of all three dietary treatments for all plasma metabolites, specific comparisons of interest were carried out to identify the metabolites that were significantly affected by the factors studied. The relative abundance of all metabolites when SER and BG diets were compared was visualized by means of a heat map (see Supplementary Fig. S1). A volacano plot of acylcarnitines and peptides comparing SER and BG diets ( Fig. 2A) revealed that SER diet group had significantly lower (P < 0.05) trans-hydroxyproline, carnosine, and aspartic acid, and higher (P < 0.05) butyryl carnitine, valeryl carnitine, alpha aminoadipic acid, butyric acid, and ornithine concentrations. The concentrations of these metabolites in the two diet groups and the related biological pathways they are involved in are shown in Table 4. The relative abundance of all metabolites when SER and BG-DW diets were compared was visualized by means of a heat map (see Supplementary Fig. S1). A volcano plot of acyl carnitines and peptides/organic acids comparing SER and BG-DW diets (Fig. 2B) revealed that Diet had effects (P < 0.05) on 21 acyl carnitines and 5 peptides/organic acids. In addition to higher butyryl carnitine, valeryl carnitine, alpha aminoadipic acid, and ornithine concentrations, carnitine, citrulline, and hydroxyvaleryl carnitine concentrations were also higher (P < 0.01) in the SER group when compared with BG-DW group (Table 5).
To evaluate changes in plasma metabolome caused by Diet (SER vs. BG) effect within an ST group, a hierarchical clustering analysis using the degree of abundance was conducted in which metabolites with identical abundance patterns were clustered together. Clustering was visualized by means of heatmap dendograms shown in Figs. 3 and 4 for NTS and TS groups, respectively. Repeated measures analysis of data separated by ST, and PCA (Fig. 5) plot revealed that a total of only 10 compounds were influenced (P < 0.05) by Diet when SER and BG groups were compared in the NTS group (Fig. 5A). However, in the TS group, 25 amino acids, 14 acylcarnitines, and 7 other organic acids and compounds were significantly (P < 0.05) affected by Diet (Fig. 5B).
Notable among the three metabolites influenced by Diet in the NTS group was the amino acid lysine concentration, which was higher (P < 0.01) in the SER group compared to BG and BG-DW groups ( Table 6). In the TS group, however, 10 amino acids (Table 7) and 15 acylcarnitines and other organic compounds (Table 8) were affected (P < 0.05) by Diet. The SER group had the lowest glycine, alanine, threonine, taurine, trans-hydroxyproline, methionine, and histidine concentrations and highest citrulline concentrations among the three diet treatments (Table 7). Carnitine, butyryl carnitine, and valeryl carnitine concentrations were higher (P < 0.05) in the SER groups compared to BG or BG-DW groups; however, 8 other acylcarnitines were highest in the BG-DW group (Table 8). Transportation time had an effect (P < 0.05) on 21 amino acids and 18 acyl carnitines and other organic acids. For a majority of amino acids, the concentrations decreased over Time in all three dietary treatments and for a majority of the acyl carnitines, the concentrations increased over Time. An example of a    Fig. 6. Diet × Time interaction effect was not significant for any of the metabolites.

Discussion
Meat goat management practices that increase animal stress may alter the normal course of conversion of muscle to meat, leading to inferior meat quality 3 . Stress in animals can negatively affect animal weight gain, immune function, and physiological status that in turn can lead to poorer health and lower carcass yield and economic returns 21 . Diets high in CT, such as SL in the form of hay or pellets, or SL feed as diet supplement have shown to ameliorate the negative effects of stress in small ruminants, possibly due to the high polyphenol content 22 . Animal management during the 24-h period prior to slaughter is crucial not only for animal welfare reasons, but also for profitability. Diet had a significant effect when body weights were measured after ST when the gut  www.nature.com/scientificreports/ Creatine levels have been used to assess muscle protein mass in small ruminants 26 and is considered to have protective action against oxidative stress 27 . The decrease in creatine concentration in blood over Time may indicate increased energy depletion, as well as susceptibility to oxidative stress due to transportation. Creatine metabolism results in the production of creatinine. Plasma creatinine concentrations were high in SER diet group, low in BG, and intermediate in the BG-DW group.
Blood metabolites, such as glucose and creatinine, have been used by researchers to monitor nutrient status and muscle mass 28 . Plasma glucose concentrations were lower and creatinine concentrations higher in the SER diet group than in BG-DW or BG groups in our study. These results are consistent with the findings of Turner et al. 28 , who also reported that circulating glucose concentrations were lower and creatinine concentrations higher in SL-fed bucklings of three different breed groups. It is likely that CT interfered with fiber digestion, leading to lower VFA production. Lower propionate production leads to low gluconeogenesis precursor and lower glucose www.nature.com/scientificreports/ production. Lower glucose combined with ST could have triggered higher cortisol release, leading to protein breakdown, for gluconeogenesis to overcome the deficit due to diet (SER). Lespedeza has been reported to be a forage of poorer quality compared to alfalfa 28 .  www.nature.com/scientificreports/ Crude protein percentages were higher in SER and BG-DW groups compared to the BG group. Lee et al. 29 studied the effects of three different dietary regimens involving different hay-concentrate combinations on proximate composition of loin chops in Boer × Spanish kids and observed no changes in protein percentages. In the present study, the three diet treatments were comprised of equal proportions of roughage and concentrate. However, goats tend to consume a higher quantity of SER hay compared to BG hay that could result in greater protein content of longissimus muscle, although the effect may be more significant in adult goats than kids 30 . Barry et al. 31 found that the crude protein percentage of longissimus muscle from tannin-fed lambs was higher compared with lambs fed the same diet supplemented with polyethylene glycol, a binding agent that negates the effects of CT.
Plasma metabolome can be influenced by genetic and environmental factors 13 . When SER group was compared individually with either BG or BG-DW for acylcarnitines, peptides and organic acids, there were significant effects of Diet on certain metabolites. Notable among the metabolites that were affected was butyric acid concentrations that were higher in the SER compared with the BG diet group. Butyric acid, a short-chain fatty acid resulting from rumen microbial fermentation of dietary fiber, has been shown to have anti-inflammatory properties 32 , and circulating levels of butyric acid have been reported to inversely correlate with inflammatory markers 33 . Butyrate has several beneficial effects in energy homeostasis, inflammation, and immune function in humans 34 . Higher butyric acid concentration in the SL-fed group may indicate better ability of goats to withstand oxidative stress. Koenig and Beauchemin 35 reported that dietary condensed tannin supplementation increased rumen butyric acid concentrations in steers. Since butyrate has a greater potential for diffusional uptake 36 , a higher concentration of butyric acid in the rumen is likely to increase its circulating concentration.
Carnitine, butyryl carnitine, and valeryl carnitine concentrations were higher in the SER groups compared to BG or BG-DW groups; however, 8 other acylcarnitine concentrations were comparable to those in BG, but lower than in the BG-DW group. Carnitine is an important compound that has functions in intermediary metabolism. Carnitine plays a major role in importing acetyl-CoA into mitochondria and acylcarnitines function as a cofactor for the transport of long-chain fatty acids through the mitochondrial membrane 37 . In mitochondria, the acylcarnitine is reconverted to free carnitine and the respective long-chain acyl-CoA. The acyl-CoA is then used to produce ATP via β-oxidation and the TCA cycle 38 . Acylcarnitine may indicate the extent of β-oxidation of fatty acids and the energy status of animals 39 . Both medium-and long-chained acylcarnitines along with free carnitine have been identified in several plant species 40 . The higher concentrations of carnitine and two acylcarnitines in the SER group compared to the BG group suggests that goats could have obtained these metabolites Table 7. Effects of dietary treatment (Diet; n = 12 goats/diet group) and transportation time (Time) on concentrations (µM) of plasma amino acids and intermediates in goats subjected to transportation stress for 90 min. abc Diet main effect means within a row with different superscripts differ significantly (P < 0.05) by LSD test. SER Serala sericea hay-fed, BG bermudagrass hay-fed, BG-DW bermudagrass hay-fed dewormed. Lower plasma β-hydroxybutyrate, a ketone body produced from lipolysis in the liver mitochondria, may indicate inhibition of β-oxidation of fatty acids 41 . It is likely that gastrointestinal parasite infection level influenced the extent of fatty acid metabolism in goats. Although carnitine and acylcarnitine concentrations were higher in the SER group compared to the BG group, the β-hydroxybutyrate concentration was lower in the SER group, indicating only a moderate level of fatty acid metabolism in this group. Among the three metabolites influenced by Diet in the NTS group was the amino acid lysine concentration, which was higher in the SER group compared to BG and BG-DW groups. A diet deficient in lysine has been reported to decrease the whole-brain content of lysine and affect norepinephrine activity in the hypothalamus 42,43 . Lysine deficiency increases release of serotonin in the amygdala, with consequent changes in psychobehavioral responses to stress, such as increased stress-induced anxiety 44 , and also results in immunodeficiency 45 in animal models. The higher plasma lysine concentration in the SER group indicates that CT-diet can confer better stress coping ability when they encounter typical preharvest livehaul situations. The metabolite α-aminoadipic acid, which was higher in the SER group compared to BG groups, is used by plants and animals to catabolize the amino acid lysine, while yeast and fungi use this derivative in the synthesis of lysine 46 . In a series of biochemical steps, lysine is catabolized to saccharopine, which is then converted to α-aminoadipic acid and then to acetyl CoA 47 . Lysine is also one of the precursors for the synthesis of carnitine 48 . Since lysine is one of the essential amino acids, its higher plasma concentration in SL-fed goats is likely from this leguminous plant.
In the TS group, the SER group had the lowest glycine, alanine, threonine, taurine, trans-hydroxyproline, methionine, and histidine concentrations and highest citrulline concentrations among the three diet treatments. The lower concentration of several amino acids in the SER group compared to the other groups studied does not mean that protein synthesis was impaired, since this correlates with intracellular amino acid pool rather than plasma amino acid concentrations 49 . Citrulline is an amino acid that plays an important role in the urea cycle. Arginine, which is synthesized from citrulline, has been reported to play an important role in modulating immune response during inflammatory stress in dairy cows 50 , and thus upregulated arginine-citrulline-nitric oxide pathway can help SER-fed goats cope up with oxidative stress and immune challenge 51 . Histidine and alanine are major amino acids involved in glucose metabolism in ruminant animals 52 . During times of feed deprivation or stress in dairy cows, gluconeogenesis is regulated by alanine to replenish glucose levels 52 . The lower alanine concentrations in the SER group compared to BG or BG-DW group indicates that carbohydrate metabolism was negatively affected by the high-CT diet. This point is further supported by lower concentrations of glucose and several other amino acids in the SER group, suggesting gluconeogenesis was negatively affected. www.nature.com/scientificreports/ Transportation time had a significant effect on 21 amino acids and 18 acylcarnitines and other organic acids. With the exception of 3 amino acids, concentrations of all other amino acids decreased with Time in the TS goats, suggesting gluconeogenesis became inadequate to maintain energy levels with increasing stress. Stress due to transportation increases with transportation time in goats 3 . Feed was withdrawn from the goats in this experiment the previous evening; therefore, feed deprivation duration also increased with transportation time. In addition to depletion of amino acids, animals can also be susceptible to inflammatory reaction due to decrease in the concentration of histidine. Low plasma histidine concentrations are associated with inflammation 53 . Histidine supplementation has been shown to influence the concentrations of other metabolites in the blood, resulting in increased concentrations of several amino acids, including glycine and aspartic acid 54 . In contrast to down regulation of amino acid-related gluconeogenesis with increasing transportation time, long-chained acylcarnitine concentrations increased with increasing time in the TS group. Concentrations of free carnitine decreased and those of 9 acylcarnitines increased with time in the TS group, whereas only the concentrations of 2 acylcarnitines increased in the NTS group, suggesting that fatty acid metabolism is upregulated with increasing levels of stress in goats and is the predominant energy generating pathway during prolonged stressful situations. This is further supported by the increasing levels of β-hydroxybutyrate concentration with increasing time, although there was only a minimal increase in the SER group.
In conclusion, feeding Serala SL diet increased plasma concentrations of butyric acid, lysine, and citrulline, all of which can have beneficial effects in goats. Elevated butyric acid concentration can help with energy homeostasis, inflammation, and immune function, the amino acid lysine can be beneficial in reducing stress www.nature.com/scientificreports/ induced psychobehavioral changes and boost immune function, and citrulline can modulate immune response during inflammatory stress. There is evidence that fatty acid metabolism is the predominant energy production pathway in the SER group, although not at the level of dewormed, BG-fed goats. Several amino acid concentrations were lower in the SER group, suggesting that gluconeogenesis either did not sustain or happened only to a limited extent compared to the BG-DW group. The BG-DW goats were able to sustain energy production via elevated levels of both gluconeogenesis and β-oxidation of fatty acids as evidenced by higher amino acid and acylcarnitine concentrations. All three diet groups switched energy generation from gluconeogenesis to fatty acid metabolism when transportation time increased, as seen by decreasing amino acid concentrations and increasing medium-and long-chained acylcarnitines with increasing transportation time. Further studies are required to establish an ideal proportion of SL hay in the diet of meat goats, such that the beneficial effects are maximized, and the negative effects are minimized. Overall, SER diet can be beneficial to goats to some extent in enhancing stress coping abilities.

Data availability
The datasets used during the current study are available from the corresponding author on reasonable request.