- Agricultural Research Station, Fort Valley State University, Fort Valley, GA, United States
Black cumin or black seed (Nigella sativa) has many beneficial biological properties, and its processing for oil extraction produces a byproduct known as black seed meal (BSM), which is utilized as an animal feed supplement. An experiment was conducted on a commercial farm to determine the effects of BSM supplementation and long-duration transportation on stress and metabolomic responses and antioxidant and immune capacities in goats. Ninety-six uncastrated male Spanish goats (4–5 months old) were randomly divided into two treatment (TRT) groups. Forty-eight goats were fed a concentrate diet containing 15% BSM, and 48 goats were fed the same diet with no BSM (control, C) in separate corrals for 3 weeks with ad libitum water. On the day of the experiment, goats were loaded onto two identical trailers (5 × 2.3 m), with 40 goats/trailer (20 goats/TRT), and were transported for 16 h to simulate a commercial situation. Blood samples were collected at 0 h (15 min after loading), 2 h, 4 h, 10 h, and 16 h of transportation (Time; n = 8 goats/Time/TRT) by jugular venipuncture. The dietary BSM supplementation in goats did not affect stress responses, except for tyramine (p < 0.05), but Time significantly affected (p < 0.05) plasma epinephrine, metanephrine, and normetanephrine. The BSM supplement did not significantly affect the antioxidant and immune status variables. At the metabolome level, 15 amino acids, 4 acylcarnitines, 24 phosphatidylcholines and sphingomyelins, and 13 other metabolites were significantly affected (p < 0.05) by TRT. Acylcarnitine (C2), hexadecenoylcarnitine (C16:1), hydroxybutyrylcarnitine (C4OH), β-hydroxybutyric acid, and iso-butyric acid concentrations were higher (p < 0.05) in the BSM goats, indicating energy supply was mainly through lipid metabolism. The BSM group had lower (p < 0.05) concentrations of glucose, 11 of the amino acids, and TCA cycle metabolites compared to the C group. Supplementation of BSM in the meat goat diet prior to extended road transportation may help them use fat as an energy source instead of breaking down protein. However, at a 15% level, there were no significant effects on antioxidant and immune status indicators determined.
Introduction
Plant by-product-based feeds are a cost-effective alternative dietary supplement for livestock and have been shown to have several beneficial effects, including improving animal health, welfare, and productivity. Black cumin (Nigella sativa) is a small annual herb that is part of the Ranunculaceae family and is used as a source of industrial seed oil and medicine. The presence of alkaloids, coumarins, saponins, flavonoids, terpenes (thymoquinone), fixed oils, and phenolics is responsible for the medicinal activity of N. sativa. The active compounds of black cumin seed have been reported to have antioxidant, neuroprotective, anti-inflammatory, immunomodulatory, and analgesic activities, among others (1, 2).
Black cumin seed processing produces a significant amount of byproduct known as black cumin meal (black seed meal, BSM), which may be utilized as an animal feed supplement (3). Due to its high crude protein (330 g/kg), essential amino acids, fat (127 g/kg), and energy contents, BSM is an excellent alternative ruminant feed supplement of high nutritional value (4). Black cumin addition to the diet of goats has been reported to significantly decrease the levels of cortisol, glucose, and total cholesterol during the months of the hot summer season (5). A 10–15% BSM inclusion is generally recommended for small ruminants (6). These authors also noted that BSM does not influence palatability, as an increase in average daily gain was associated with increased dry matter intake. Among other compounds, black cumin contains thymoquinone, a major phytochemical bioactive ingredient (7), which makes up about 30–48% of the volatile oil fraction from black seed (8), although its content in the extract varies depending on the geographical region where the black cumin seeds are produced (9, 10).
Thymoquinone has been demonstrated to modulate the hypothalamic–pituitary–adrenal (HPA) axis and reduce hyperactivity, leading to decreased stress levels in animal models (11). N. sativa enhances the scavenging of free radicals, increases activities of superoxide dismutase, glutathione peroxidase, and catalase, and inhibits lipid peroxidation and NF-kB activity, thus imparting a protective effect against oxidative stress (12). This compound inhibits NF-kB by activating the Nrf2/ARE signaling pathway (13). Treatment with black seed and thymoquinone suppresses pro-inflammatory cytokines, such as IL-6 and TNF-alpha, in rats, indicating reduced inflammation and enhanced cell-mediated immunity (8, 14). Thus, thymoquinone can act as an antioxidant agent, inhibiting non-enzymatic peroxidation, increasing immunity, and helping animals tolerate stress (15).
Goats can perform well in places not ideally suited for other livestock and thus are hardy in nature (16). Spanish goats, commonly raised for meat production in the southeastern US, are known to browse a variety of plants and tolerate bitter tastes reasonably well, like most other breeds. Meat goats are frequently moved across vast distances in commercial settings that invariably cause stress and protracted negative aftereffects. During stressful situations, animals experience physiological changes that help them maintain bodily homeostasis; however, when stress becomes too intense and prolonged, these physiological systems are unable to cope, and the negative impacts of stress become obvious. The well-being of animals during transit is becoming an increasingly important societal concern. Antioxidant status, immune function, and physiological and metabolomic responses are all affected by intense stress, such as transportation (17, 18).
Researchers have used numerous physiological markers to assess stress in food animals, and studies have also focused on the accuracy of these animal welfare indices (19). Catecholamines (epinephrine and norepinephrine) and glucocorticoids are thought to form the central components of the endocrine response during stress. Together, these hormones help to orchestrate the body’s response to stress. The plasma levels of catecholamines increase during transportation stress, which may indicate simultaneous activity of both the adrenal cortex and adrenal medulla (20).
Metabolomics, a large-scale study of small molecules, commonly known as metabolites, within cells, biofluids, tissues, or organisms, is emerging as an effective tool to assess animal stress (17, 21, 22). These metabolites can most directly reflect an animal’s physiological status and can be influenced by genetic and environmental factors (22). For example, the plasma metabolome in goats is significantly impacted by stress, with amino acid levels decreasing and medium- and long-chain acylcarnitine concentrations increasing with increasing stress duration (17).
To our knowledge, the information available on the effects of dietary BSM supplementation on the physiological responses to intense and prolonged stress and distress in goats is scant. This study was conducted under commercial conditions to determine the effects of feeding a BSM supplement for 3 weeks before long-duration transportation on stress responses, including antioxidant and immune statuses and metabolomic profiles in goats.
Materials and methods
Study design
Ninety-six uncastrated male Spanish goats (4–5 months old) were used in this experiment conducted at a commercial production facility. Goats were randomly divided into two treatments (TRT) groups and fed either a concentrate diet containing BSM (BSM group; n = 48 goats) or a regular concentrate feed mixture with no BSM [control (C) group; n = 48 goats] for 3 weeks. The ingredients for both dietary treatments are shown in Table 1. The ingredients were the same for both treatments, but one had BSM included at 15% of the diet, since anecdotal evidence shows that it is a common practice among meat goat producers who have access to this product to include it in the goat feed at 15%. Both TRT groups were kept in separate corrals for 3 weeks with ad libitum water and feed (large feeder with adequate feeder space) and shelter. Animals in both treatment groups were deprived of feed overnight prior to transportation. One day prior to transportation, 8 goats from each dietary group were blood sampled to establish baseline values before stress-related changes occurred in the response variables; however, these 16 goats were not transported. On the day of the experiment, goats were loaded onto two identical trailers (5.0 × 2.3 meters) pulled by identical trucks, with 40 goats in each trailer (20 from each TRT), and were transported for 16 h. The sides of the trailers were provided with windows to allow adequate ventilation. The average temperature and humidity measured every 2 h were 26.1 °C and 29.3%, respectively, in trailer 1, and 26.6 °C and 25.6%, respectively, in trailer 2. The goats from the two treatments were kept separate in the two compartments of the trailer (Figure. 1). Blood samples were collected from randomly selected individual animals at 15 min after loading (0 h), 2, 4, 10, and 16 h of transportation (Time; n = 8 goats/Time/TRT). Representative samples of the two concentrate feeds were collected and analyzed for nutrient content at a commercial laboratory (Dairy One Forage Laboratory, Ithaca, New York, Table 1).
Table 1. Ingredients and chemical composition of experimental diets (Black seed meal, BSM, diet; Control diet, C).
Figure 1. Diagram showing the assignment of goats to treatments (TRT), trailer compartments, and time sampling (time).
Blood sampling
For 2, 4, and 10 h sampling, the trucks were stopped for 10 min each time, one experienced animal handler and one trained blood sampler entered the trailer, and a blood sample was collected from each goat by jugular venipuncture into K2EDTA-coated vacutainer tubes. The two individuals were deliberate in their approach to avoid any loud noise or rough handling that could agitate the animals. A different set of goats was sampled at each time, and each goat was marked on the horns with a colored marker after sampling to avoid being sampled again. Blood samples were kept on ice until the separation of plasma. The samples were centrifuged (LW Scientific E8 Portafuge, Lawrenceville, GA) immediately at 2,650 × g for 10 min using a portable centrifuge to separate plasma, and the samples were kept frozen on dry ice and transported to the lab. Plasma aliquots were stored at −80 °C for analysis.
All 96 plasma samples from baseline (8/TRT) and experimental animals (40/TRT) were shipped to The Metabolomics Innovation Center (TMIC) at the University of Alberta (Edmonton, Alberta, Canada) for targeted metabolomics and catecholamine analyses. Aliquots of 400 μL of plasma samples were transferred into labeled Eppendorf tubes and then frozen. The Eppendorf tubes were placed in the appropriately labeled cardboard cryo/freezer box, sealed with tape, and then placed in the ziploc bags. The cardboard cryo/freezer box and/or ziploc bag(s) were placed into the styrofoam container. Enough dry ice (at least 2 pounds for every shipping day) was placed in the styrofoam box to keep the samples cold for approximately 48 h. The styrofoam box was sealed, but not taped, to prevent any explosion due to the buildup of CO2 gas. The samples and dry ice were prepared and packaged in accordance with the International Air Transport Association or the US equivalent rules.
Glutathione peroxidase
Plasma glutathione peroxidase (GPx) concentrations were determined using a commercially available Goat GPx ELISA kit (Cat. No. MBS1601455, MyBioSource, Inc., San Diego, CA, United States). Briefly, the plate was precoated with goat GPx antibody, and the sample containing GPx was added, which resulted in the binding of GPX to antibodies coated on the wells. Then, a biotinylated goat GPx antibody was added that bound to GPx in the sample, followed by the addition of Streptavidin-Horseradish Peroxidase (HRP) that bound to the biotinylated GPx antibody. After incubation, the unbound Streptavidin-HRP was washed away during the washing step. A substrate solution was then added, and the color that developed was proportional to the amount of goat GPx. The reaction was terminated by the addition of the acidic stop solution, and the absorbance was measured at 450 nm. The concentrations were determined using a standard curve (1.0–300.0 ng/mL). The sensitivity of the assay was 0.58 ng/mL, and the intra- and inter-assay coefficients of variation were < 8 and <10%, respectively.
Superoxide dismutase
Plasma superoxide dismutase (SOD) concentrations were determined using a commercially available Goat SOD ELISA kit (Cat. No. MBS268466, MyBioSource, Inc., San Diego, CA, United States). The kit employs the double antibody sandwich technique, which is based on the characteristics of a target analyte with more than two possible epitopes identified by both the pre-coated capture antibody and the detection antibody simultaneously. In this assay, the pre-coated antibody is an anti-goat SOD monoclonal antibody, while the detection antibody is a biotinylated polyclonal antibody. Briefly, the samples and biotinylated antibodies were added to the ELISA plate wells and washed out with buffer solution after respective addition to the wells. Avidin-peroxidase conjugates were then added to the wells, and after thorough washing, TMB substrate was used for color development. The blue coloration, because of peroxidase activity, finally turned yellow after the addition of the stop solution. The color intensity and the quantity of SOD in the sample were positively correlated. The sensitivity of the assay is up to 0.06 ng/mL according to the kit manufacturer.
Thiobarbituric acid reactive substances
Thiobarbituric acid reactive substances (TBARS) in plasma samples were determined using the Goat TBARS ELISA kit (Cat. No. MBS2614837, MyBioSource, Inc., San Diego, CA, United States). The principle involved in this double antibody sandwich technique is similar to that explained above for the SOD assay.
Immunoglobulin A
Plasma immunoglobulin A (IgA) was determined using the Goat IgA ELISA kit (Cat. No. MBS734771, MyBioSource, Inc., San Diego, CA, United States), which applies the competitive enzyme immunoassay technique utilizing a polyclonal anti-IgA antibody and an IgA-HRP conjugate. The assay sample and buffer were incubated together with OgA-HRP conjugate in pre-coated plate for 1 h. After the incubation, the wells were decanted and washed five times. The wells were then incubated with a substrate HRP enzyme. The product of the enzyme-substrate reaction formed a blue colored complex. A stop solution was added to stop the reaction, which turned the solution yellow. The intensity of the color was measured spectrophotometrically at 450 nm in a microplate reader (Synergy HTX Microplate Reader, Bio-Tek, Winooski, VT). A standard curve was created relating the intensity of the color to the concentration of the standards, and the IgA concentration in each sample was interpolated from this standard curve. Sensitivity in this assay is 0.1 mg/mL.
Immunoglobulin G
A commercially available Goat Immunoglobulin G (IgG) ELISA kit (Cat. No. MBS734771, MyBioSource, Inc., San Diego, CA, United States) was used to determine the IgG concentrations in plasma samples. The capture antibody was pre-coated onto 96-well plates, and the biotin-conjugated antibody was used as a detection antibody. The standards, test samples, and biotin-conjugated detection antibody were added to the wells subsequently and washed with wash buffer. The TMB substrates were used to visualize the HRP enzymatic reaction. The TMB was catalyzed by HRP to produce a blue color product that changed to yellow after the addition of the acidic stop solution. The absorbance values were read using a microplate reader, and the unknown concentrations in samples were calculated. The sensitivity of this assay is 0.938 ng/mL.
Interleukin 6
A Goat Interleukin 6 (IL6) ELISA commercial kit (Cat. No. MBS025544, MyBioSource, Inc., San Diego, CA, United States) was also used for determining plasma IL6 concentrations. The sensitivity of this kit is 2.0 pg./mL, the detection range is 12.5–400.0 pg./mL, and both intra- and inter-assay coefficients of variation are less than 15%.
Catecholamines and derivatives
Catecholamine profiling was done on plasma samples at TMIC using the method described by Zheng et al. (23). A combined direct injection mass spectrometry (DI-MS) with a reverse-phase mass spectrometry combined with liquid chromatography (LC–MS/MS) custom assay using an ABSciex 4,000 Qtrap (Applied Biosystems/MDS Sciex) tandem mass spectrometry instrument (Applied Biosystems/MDS Analytical Technologies, Foster City, CA, United States) with an Agilent 1,260 series UHPLC system (Agilent Technologies, Palo Alto, CA, United States) was used for analysis of catecholamines and other biogenic amines. The data were analyzed using Analyst 1.6.2.
Metabolomics
At TMIC, a targeted quantitative metabolomics approach was employed to analyze the samples using a combination of DI-MS with a reverse-phase LC–MS/MS assay. This custom assay, in combination with mass spectrometry, 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 (24). The method combines the derivatization and extraction of analytes and the selective mass-spectrometric detection using multiple reaction monitoring (MRM) pairs. The procedure was followed as described in an earlier publication from our lab (21).
Statistical analysis
In this study, there were only two corrals, one assigned to each of the dietary treatments. Since the corral constitutes the experimental unit for the diet effect, and any effect of TRT cannot be statistically separated from other factors associated with the corral, the results of the diet effect from this exploratory study should be interpreted with caution. For the Time effect, each set of 4 goats that was time-sampled/trailer is considered the experimental unit. The baseline values were not used in the experimental model, but the means and SEM were presented where appropriate. Statistical analysis of antioxidants and immune function data was conducted using mixed procedures in SAS (release 9.4, SAS Institute, Cary, NC, United States) with TRT, Time, and TRT × Time as fixed effects. The data were first examined for normality and homogeneity of variance using the Shapiro–Wilk’s test and Levene’s test, respectively. Log transformation was used when the data did not meet the assumptions of ANOVA; however, they were backtransformed to the original scale before being presented. When significant by ANOVA, the means were separated using the pdiff procedure.
Data from all 80 samples were used for metabolomics and catecholamine statistical analyses. Preprocessing of data and general approach to statistical analysis are similar to those described in Batchu et al. (21). Briefly, since the data for all groups were not normally distributed, Mann–Whitney U rank method was used for univariate analysis, and the effect size (Cliff’s Delta method) and fold change (ratio between group medians) were determined. Kruskal–Wallis test was used for one-way ANOVA, and Dunn’s test with Benjamini-Hochberg False Discovery Rate (FDR) correction was used for multiple comparisons. Two-way ANOVA and post hoc tests were conducted on log-transformed data, using the Benjamini-Hochberg FDR method to correct p-values for multiple comparisons. For all types of comparisons, principal component analysis (PCA) and partial least squares discriminant analysis (PLS-DA) were performed using Metaboanalyst R. Cross-validation and permutation testing were performed for the PLS-DA model. In PLS-DA, Metaboanalyst reports R2 as R2Y (variance in Y explained), Q2 (cross-validated predictive ability). The accuracy, R2, and Q2 of the PLS-DA model were 0.98, 0.92, and 0.80, respectively. The 0.98 value indicates very high prediction accuracy. The high R2 of 0.92 indicates that the model explains the classes well. An R2Y of 0.92 means the model explains 92% of the variation in the class labels (BSM vs. C). The Q2 of 0.80 suggests the model is robust and can reliably predict new samples. The Q2 reflects out-of-sample performance and, together with the permutation test, serves as a key guard against overfitting. The permutation test p-value was 0.5e-4, indicating the PLS-DA results were statistically significant. Variable importance in projection (VIP) score plots were then created, with a score of >1.0 indicating that the metabolite is significantly involved in the separation of the classes.
A pathway enrichment analysis was performed using the Metaboanalyst 6.0 platform. The input consisted of a combined list of all metabolites identified as significantly different (p < 0.05) between the BSM and C groups. This type of analysis identifies which metabolic pathways are the most over-represented by this list of significant metabolites, highlighting the biological processes most impacted by dietary treatment.
Results
Antioxidant status
Glutathione peroxidase, SOD, and TBARS levels were not affected by TRT or TRT × Time interaction. However, Time had a significant effect on GPx (p < 0.01; Figure 2A), SOD (p < 0.05; Figure 2B), and TBARS (p < 0.01; Figure 2C) concentrations in goats. The mean GPx values increased at 2 and 4 h, reached low levels at 10 h, and remained at that level at 16 h. Plasma SOD concentrations gradually decreased over transportation time and reached low levels at 10 and 16 h. The TBARS values were higher at 10 and 16 h than at 0, 2, and 4 h.
Figure 2. Effects of treatment (BSM = Black seed meal; C = Control) and transportation time on plasma (A) glutathione peroxidase, (B) superoxide dismutase, and (C) thiobarbituric acid reactive substances (TBARS) concentrations in goats. (A–C) Time main effect means (n = 16 goats/time) with different letters differ significantly by pdiff procedure at p < 0.05. Baseline values were not included in the statistical analysis.
Immune status
The overall IgA values were noted to be higher (p < 0.01) in the C group than in the BSM group (Figure 3A), primarily due to the differences in values between TRT at 0 and 2 h. Time also had a significant effect (p < 0.01; Figure 3B) on plasma IgA; however, there was no clear pattern with the values moving up and down over transportation time. The IgG status was not affected by any of the factors studied (Figure 3C). Interleukin 6 values in both TRT groups decreased over time and reached their lowest values at 10 and 16 h (p < 0.01; Figure 3D), while TRT or TRT × Time did not have any significant effect.
Figure 3. Effects of treatment (BSM = black seed meal; C = control) and transportation time on plasma (A,B) immunoglobulin A, (C) immunoglobulin G, and (D) interleukin 6 concentrations in goats. a,bBars representing treatment main effect means (n = 40 goats/TRT) with different letters are significantly different by PDIFF procedure at p < 0.05. a,b,c,dTime main effect means (n = 16 goats/Time) with different letters differ significantly by procedure at p < 0.05. Baseline values were not included in the statistical analysis.
Catecholamines and derivatives
The effects of TRT on the concentrations of catecholamines and their derivatives are depicted in a heatmap (Supplementary Figure 1) created using normalized concentration ranges (from 0 to 1). High, medium, and low points on the color bar correspond to 0.95, 0.5, and 0.05, respectively. Tyramine concentration was observed to be higher (p < 0.05) in the C group compared to the BSM group of goats (Supplementary Figure 2). The BSM treatment did not have any significant effect on dopamine, phenylethylamine, 5-methoxytryptamine, epinephrine, norepinephrine, metanephrine, and normetanephrine concentrations.
Time had a significant effect (p < 0.05) on metanephrine, normetanephrine, and epinephrine concentrations (Supplementary Figure 3). The concentration of normetanephrine decreased with increasing transportation time, while the concentrations of epinephrine and metanephrine increased during the first 4 h of transportation and then decreased. The effects of transportation time on the concentrations of catecholamines and their derivatives are shown using a heatmap (Figure 4). The interaction effects were not significant.
Metabolomics
Involvement in biological processes related to antioxidant capacity, immune function, and/or energy metabolism was used as the primary metabolite selection criterion in this study. Intense stress in goats has been shown to negatively affect these biological processes, while thymoquinone, present in black seed, has been reported to help animals combat these negative effects.
At the metabolome level, 15 amino acids, 4 acylcarnitines, 24 phosphatidylcholines and sphingomyelins, and 13 other metabolites were affected (p < 0.05) by TRT. Of the significantly affected amino acids, 11 (phenylalanine, tyrosine, tryptophan, alanine, glutamine, methionine, proline, glutamic acid, taurine, threonine, and asparagine) were higher (p < 0.05) in the C group compared to BSM group, while isoleucine, valine, ornithine, and methyl histidine were lower (p < 0.05) in the C group compared to the BSM group (Table 2 and Supplementary Figure 4). The concentrations of three of the acylcarnitines [acetylcarnitine (C2), hexadecenoylcarnitine (C16:1), and hydroxybutyrylcarnitine (C4OH)] were lower (p < 0.05) in the C group than the BSM group, while the propionylcarnitine (C3) concentration was higher (p < 0.05) in the C group compared to the BSM group (Table 3 and Supplementary Figure 5). In this study, treatment appeared to have effects (p < 0.05) on 15 phosphatidylcholines and 9 sphingomyelins (Table 4 and Supplementary Figure 6), in addition to several other metabolites (Table 5 and Supplementary Figure 7).
Table 2. Amino acids significantly (P < 0.05) affected by treatment (BSM = Black seed meal; C = Control; n = 40 goats/TRT) in goats.
Table 3. Acylcarnitines significantly (P < 0.05) affected by treatment (BSM = Black seed meal; C = Control; n = 40 goats/TRT) in goats.
Table 4. Phosphatidylcholines and sphingomyelins significantly (p < 0.05) affected by treatment (BSM = Black seed meal; C = Control; n = 40 goats/TRT) in goats.
Table 5. Other metabolites significantly (P < 0.05) affected by treatment (BSM = Black seed meal; C = Control; n = 40 goats/TRT) in goats.
A volcano plot identified that 19 metabolites were affected by TRT (Figure 5). The relative abundance of metabolites affected by TRT was visualized by means of a heat map (Figure 6). The PCA plot created to separate metabolites by TRT (by principal components 1 and 2) revealed partial clustering of the BSM and C groups (Supplementary Figure 8). The total variance of the principal components for the two groups contributed to 33.3% of the PCA model (PC 1 = 21.9%, PC 2 = 11.4%). To maximize the separation of the groups observed by PCA, a PLS-DA was applied. In the PLS-DA model for the two groups, the two components contributed to a total variance of 31.3% (Component 1 = 16.9%, Component 2 = 14.4%). The two groups showed minimal overlapping, with the cluster larger in C compared to the BSM group (Figure 7). When averaged across all time points, lysophosphatidylcholine C17:0, acyl-alkyl-phosphatidylcholine C36:0, lysophosphatidylcholine C18:2, acetylcarnitine, β-hydroxybutyric acid, hydroxysphingomyelin C22:2, tyrosine, sphingomyelin C20:2, carnosine, methylmalonic acid, diacyl-phosphatidylcholine C32:2, propionic acid, alanine, hexadecenoylcarnitine, and phenylalanine were the top 15 metabolites identified by the PLS-DA multivariate model (p < 0.05) and VIP values as having the greatest influence (VIP scores >1.5) in separating the BSM and C group (Figure 8).
Figure 5. Volcano plot showing metabolites significantly affected by dietary treatment. The names of metabolites with p-values below 0.05 and fold change < 0.77 or > 1.3 are shown within the plot, and the absence of metabolite names indicates that none of the metabolites satisfied these criteria.
Figure 7. PLS-DA plot of principal components 1 and 2 for treatment classes (BSM = Black seed meal, C = Control; p < 0.05) in goats.
Figure 8. PLS-DA VIP plot showing the metabolites (VIP scores > 1.5) that significantly contribute to the differences between treatment (BSM = Black seed meal, C = Control) groups in goats. The metabolite concentrations averaged across all time points were used in PLS-DA model (p < 0.05).
Time had a significant effect (p < 0.05) on 18 amino acids (Supplementary Figure 9), 5 acylcarnitine (Supplementary Figure 10), and 7 phosphatidylcholines and sphingomyelins (Supplementary Figure 11). In addition, several other metabolites were also significantly affected by Time (p < 0.05; Supplementary Figure 12). The relative abundance of metabolites significantly affected by Time was visualized by means of a heat map (Supplementary Figure 13). The PLS-DA plot showed a total variance of 28.8% (Component 1 = 13.2%, Component 2 = 15.6%), as shown in Supplementary Figure 14.
Interaction effects were significant (p < 0.05) for pyruvic acid and eight amino acids (phenylalanine, tyrosine, tryptophan, glutamine, methionine, threonine, asparagine, and glutamic acid). Examples of amino acids with significant TRT × Time effects, along with pyruvic acid, are shown in Figure 9. The interaction effect was also significant (p < 0.05; Supplementary Figure 15) for one acylcarnitine [hexadecanoylcarnitine (C16:1)] and one phosphatidylcholine (lysophosphatidylcholine C17:0). The differences in energy-sourcing strategies in BSM and C groups are presented in the form of a mechanistic diagram in Supplementary Figure 16.
Figure 9. Examples of amino acids (A phenylalanine, B asparagine, C glutamine, D glutamic acid) and pyruvic acid (E) showing significant TRT × time effect (p < 0.05; n = 8 goats/time/TRT). Treatment means are different at time points with asterisk(s) (*p < 0.05, **p < 0.01). Baseline values were not included in the statistical analysis.
The pathway analysis identified that the most significantly impacted pathways were overwhelmingly related to amino acid metabolism and central energy production. The most prominent pathways included (i) the alanine, aspartate, and glutamate metabolism, (ii) the arginine and proline metabolism, (iii) the valine, leucine, and isoleucine degradation, and (iv) the phenylalanine metabolism. The tricarboxylic acid cycle (TCA cycle) was also identified as a highly significant pathway (Figure 10).
Figure 10. Results of the pathway analysis conducted using a list of all metabolites identified as significantly different (p < 0.05) between the BSM and C groups.
Discussion
In the current study, GPx and SOD concentrations and TBARS values were not differentially influenced by TRT, and the time patterns in both TRT groups were similar. These patterns may suggest that transportation caused extreme stress, which could have overridden the differential effects of dietary treatment. Among other compounds, black cumin contains thymoquinone, which has been reported to act as an antioxidant agent, as it inhibits the non-enzymatic peroxidation, increases immunity, and helps animals in tolerating heat stress (15). The effect noted in our study was the reverse of what has been reported in the literature. Furthermore, IgA levels were higher in the C group compared to the BSM group in the present study. Nigella sativa extracts at high concentrations have been reported to suppress lymphocyte response to mitogens and phagocytic activity of polymorphonuclear leukocytes in vitro (25). Thymoquinone has been reported to have higher antioxidant activity at lower levels and has a prooxidant effect at higher levels (26), which may suggest that the inclusion of 15% BSM in goat feed, as practiced by meat goat producers, could be higher than optimum. Similar products have been tested in other countries in livestock and poultry. For example, Ramdani et al. (6) studied the effect of 5, 10, and 15% inclusion levels in lambs and observed that approximately 10% BSM in the diet is optimum for increasing average daily gain, although the optimal level of BSM in small ruminants that results in the most favorable effect on antioxidant and immune capacities is not known. The effects of BSM on antioxidant and immune status indicators need further investigation, particularly to determine an ideal level for inclusion in the diet that could help animals combat intense stress conditions.
The elevated epinephrine and metanephrine concentrations at 2 and 4 h of transportation and their gradual decline thereafter were consistent with previous findings in goats, for example, by Nwe et al. (27) and Kannan et al. (28). This reduction in blood concentrations of these catecholamines is due to the decrease in fear and emotional stress caused by the initial novelty of the environment and experience, which also reflects variations in the activity of the sympatho-adrenal medullary axis over time during transportation in goats (27). Metanephrine and normetanephrine are O-methylated metabolites of catecholamines, epinephrine and norepinephrine, respectively, and are reliable indicators of stress under practical circumstances (28). Among the 8 catecholamines analyzed, tyramine concentrations were significantly affected by TRT. The BSM group had lower concentrations than the C group. Tyramine is a naturally occurring trace amine that acts as a catecholamine-releasing agent. As an indirect sympathomimetic agent, it releases norepinephrine from sympathetic nerve terminals (29). Tyramine is part of a stress response, and its concentration increases when an animal encounters stress. In goats, a previous experiment showed that tyramine concentrations increased with increasing levels of stress (30). The lower mean tyramine concentrations in the BSM goats compared to the C goats may indicate that BSM supplementation likely reduced plasma tyramine concentrations, thus indicating a lower stress response.
Phenylalanine or tyrosine serves as the precursor to catecholamines dopamine, norepinephrine, and epinephrine (31, 32). Tyrosine is also converted to tyramine via decarboxylation catalyzed by the enzyme tyrosine decarboxylase (33). In our study, the levels of amino acids phenylalanine and tyrosine, along with 9 other amino acids, were lower in the BSM group compared to the C group of goats, which may also explain the lower tyramine concentrations in the BSM group. Overall, BSM had no significant effect on catecholamine responses, except for tyramine, while the responses over time during transportation confirmed the effects noted in previously reported studies.
The primary effect that became apparent from metabolomics analysis was the change in energy-sourcing strategies used by the two groups of goats when exposed to transportation stress. Lipid metabolism predominated in the BSM goats, and protein catabolism was pronounced in the C goats, resulting in elevated amino acid concentrations in the latter group (Supplementary Figure 16). The higher concentrations of phenylalanine, tyrosine, tryptophan, glutamine, methionine, threonine, asparagine, and glutamic acid in the C group than the BSM group during the first 4 h of transportation, and the absence of differences between the two groups as the journey prolonged, resulted in significant TRT × Time interaction effects. Not only does this time pattern indicate that the BSM diet reduced stress during the first 4 h of transportation, but it also suggests that in the C goats, there was increased mobilization of amino acids from muscles for gluconeogenesis. The higher pyruvic acid concentrations in the C group compared to the BSM group at 10 h and 16 h of transportation, noticed in our study, further confirm that the amino acids are being used for glucose production, as these are both glucogenic and ketogenic (phenylalanine, isoleucine, threonine, tryptophan, and tyrosine) in function (34). The increase in some of these amino acids noticed in the C group could also be due to a higher need for neurotransmitter production during stress. The overall concentrations of 11 amino acids in the BSM group were lower compared to the C group, despite the fact that black cumin seeds contain high amounts of essential amino acids (9). This may suggest that short-term dietary BSM supplementation to goats may have the beneficial effect of lowering the stress response and resultant protein catabolism during long-duration transportation, as shown by lower plasma glucose concentration in the BSM group.
Transportation and treatment had a significant interaction effect for hexadecanoylcarnitine (C16:1), which has been associated with fasting and ketosis, among other conditions (35). Acetylcarnitine (C2), hexadecenoylcarnitine (C16:1), and hydroxybutyrylcarnitine (C4OH) were higher in the BSM group compared to the C group. The surge in the end products of short-chain acylcarnitine (C4OH) correlates with the more commonly measured ketone species β-hydroxybutyrate (36). The β-hydroxybutyrate levels were higher in the BSM group in our study. The rise in ketone levels supports the idea that there has been a switch from carbohydrate to lipid metabolism (37). Ruminants utilize ketone bodies as a source of energy in their peripheral tissues and small intestines, which helps them avoid being deprived of glucose (38). A strong positive correlation has been reported between β-hydroxybutyric acid and C4OH levels, and a moderate correlation between β-hydroxybutyric acid and C2 in the blood (36). The authors also observed that, in accordance with the changes in carnitine, there was an increase in C2 and C4OH with increasing ketosis. These results suggest that short-term dietary BSM supplementation in goats before long-duration transportation may help the animals utilize fat reserves as a source of energy, rather than relying on muscle breakdown, for energy production.
Hexadecenoylcarnitine (C16:1), a long-chain acylcarnitine with unsaturated fatty acid moiety, and hexenoylcarnitine (C6:1), a medium-chain acylcarnitine with unsaturated fatty acid moiety, increased with increasing transportation time in our study. During intense physical stress, acylcarnitines are exported from muscle cells as the carnitine acylation state is higher in cytoplasm than in mitochondria (39). In a moving vehicle, muscles are constantly active in animals to maintain posture and balance. Muscle contraction, which increases glucose and fatty acid oxidation, is the main cause of the increase in long- and medium-chain acylcarnitines in the blood (40). Increased lipolysis and the resulting faster β-oxidation rate than the TCA cycle could be the cause of the increase in long-chain acylcarnitines (41, 42). The ensuing accumulation of fatty acids in the matrix could result in mitochondrial stress and incomplete fatty acid oxidation, leading to acylcarnitines entering circulation (43), as evidenced by the increasing concentrations of long- and medium-chain acylcarnitines with increasing transportation duration.
Sphingolipids are involved in various cellular, developmental, and stress-response processes. Sphingomyelin is produced when a phosphocholine headgroup is transferred from phosphatidylcholine to ceramide, resulting in the production of diacylglycerol (DAG) and sphingomyelin. The most prevalent complex sphingolipids are the sphingomyelin species, which are required for cell survival (44). In our study, TRT had an impact on nine sphingomyelins, with five of these being higher in the BSM group. However, the hydroxylated forms were higher in the C group, and they are known to play a role in cell signaling, including responses to stress and inflammation. A previous study in goats suggested that increased fatty acid metabolism may contribute to higher concentrations of sphingomyelins (21). Ceramide, a substance that can lead to depression, may also increase as plasma sphingosine and sphinganine levels rise (45). In the present study, diet had an impact on 15 phosphatidylcholines, and long-duration transportation stress affected the concentrations of six phosphatidylcholines in goats. The decrease in phosphatidylcholines and their derivatives can be attributed to abnormal metabolism of phospholipids and cell membrane injury, and the increased levels of phosphatidylcholines, sphingomyelins, and their derivatives can be associated with increased fatty acid metabolism due to transportation stress. However, because there were changes in the concentrations of phosphatidylcholines and their derivatives in both TRT groups, the effect of the BSM diet on goats is not clear under the conditions of the experiment and requires future investigation.
Carnosine levels were higher in the C group compared to the BSM group. Carnosine is useful in preventing the damage caused by stress in animals. It is formed from the binding of the amino acids, alanine and histidine, which present a binding site for glucose. Carnosine has been reported to improve glucose metabolism in stressed animals. The increased carnosine concentration in the C group was likely because these goats experienced higher stress levels and thus required an elevated carnosine-mediated coping mechanism. In our study, creatinine decreased during the initial 4 h of transportation and then gradually increased with increasing transportation time. Creatinine is a non-protein nitrogenous compound that is produced by the breakdown of creatine in muscle. Creatine conversion to phosphocreatine is catalyzed by creatine kinase, with spontaneous formation of creatinine during the reaction (46). Muscle activity, heat stress, dehydration, and glomerular filtration rate have all been reported to increase blood creatinine concentrations. In the present study, creatinine levels began to increase after 4 h of transportation, while stress levels decreased based on epinephrine concentrations. Dehydration, rather than stress response, likely caused the increase in creatinine after 4 h, as the goats had access to ad libitum water until they were loaded onto the trailers, and dehydration did not occur until after the initial few hours. Previous studies in small ruminants have shown increases in creatinine concentrations after 20 h of transportation (47), severe water restriction (48), or heat stress (49).
Black cumin has been reported to reduce blood glucose by lowering hepatic gluconeogenesis (50). In our study, gluconeogenic metabolites (propionic acid and methylmalonic acid) and glucose were found to be lower in the BSM group than in the C group. The lower glucose and amino acid concentrations in the BSM group compared to the C group may be due to lower stress, proper glucose uptake by the muscle cells, lower hepatic gluconeogenesis, and lower protein degradation. Methylmalonic acid is produced in the body when it is necessary for energy production by breaking down proteins (51). The higher unsaturated acylcarnitines and low amino acid concentrations in the BSM group may also indicate the use of fat reserves for energy, instead of muscle breakdown or protein catabolism. As BSM lowers hepatic gluconeogenesis, the amino acids are likely not catabolized into TCA cycle metabolites (α-ketoglutaric acid, succinic acid, and fumaric acid), as noticed in the present study. The higher isobutyric acid concentrations in the BSM group further confirm the predominance of lipid metabolism in the BSM group.
Pathway analysis results confirm the paper’s central finding: the two groups of goats used different energy-sourcing strategies during transportation stress. The C group had significantly higher concentrations of numerous amino acids and TCA cycle intermediates. This is interpreted as the C group breaking down muscle protein for energy. The pathway analysis strongly corroborates this by flagging these exact amino acid and energy pathways as the primary areas of metabolic difference between the BSM and C groups. Since the metabolomics analysis revealed prominent effects mainly on energy-sourcing strategies in the two groups of goats, the metabolomics results could not be related to the antioxidant and immune status variables determined in this study. The catecholamine analysis to assess stress levels during transportation revealed significant effects on tyramine concentrations, which were attributed to dietary treatment. While the effects of stress on antioxidant and immune capacities in animals are well established, the differences in antioxidant and immune capacity indicators between the two groups of goats after exposure to intense stress were not immediately evident in the present experiment. It can be speculated that changes in the markers related to these variables may manifest 24–48 h after transportation under the conditions of this study.
Limitations of this study were that it was conducted on a commercial farm under routine practices typically followed on that farm. Therefore, replications of corrals could not be established, nor could the feed intake by individual animals during the three-week treatment period be recorded. Therefore, dietary effects should be interpreted with caution. Additionally, ether extract, crude protein, and thymoquinone contents in the BSM batch used in this study were not available. Since the primary clients of this farm were other meat goat farmers nationwide and small ruminant researchers at universities, the goats were fed a concentrate diet exclusively, particularly during the summer. This is not representative of meat goat farms in the southeastern US, where goats are typically raised on pastures with concentrate and hay supplements and sent to processing plants either directly or through livestock auctions.
Conclusion
Black seed meal supplementation did not differentially affect physiological responses, except tyramine concentrations. The lower glucose and TCA cycle metabolites in the BSM group suggest an ability of BSM to lower stress, enhance proper glucose uptake by the muscle cells, lower hepatic gluconeogenesis, and lower protein degradation. Metabolomic profiles suggest a shift toward lipid utilization when goats were fed BSM for a brief period before long-duration transportation. However, physiological and immune status indicators could not be related to this conclusion. Despite higher concentrations of amino acids that can help enhance immune function, energy balance, and anti-inflammatory activity, the antioxidant and immune status indicators determined did not support this conclusion. Further studies are required under controlled conditions that enable monitoring of feed intake and weight gain at different levels of BSM in the diet. The data will help us better understand the optimal level of inclusion of this product in the diet, which has a positive impact on animal welfare, stress resistance, and productivity.
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 approved by FVSU Agricultural and Laboratory Animal Care and Use Committee. The study was conducted in accordance with the local legislation and institutional requirements.
Author contributions
PG: Formal analysis, Writing – original draft, Visualization, Data curation, Methodology, Conceptualization, Investigation, Writing – review & editing. PB: Validation, Data curation, Methodology, Supervision, Investigation, Writing – review & editing. AS: Methodology, Formal analysis, Investigation, Writing – review & editing. TT: Project administration, Conceptualization, Methodology, Investigation, Writing – review & editing, Resources, Funding acquisition. GK: Funding acquisition, Supervision, Formal analysis, Writing – review & editing, Investigation, Visualization, Resources, Project administration, Methodology, Conceptualization.
Funding
The author(s) declared that financial support was received for this work and/or its publication. This research was partially funded by the USDA-NIFA Evans-Allen Program.
Acknowledgments
We thank Gregory Dykes, Carlton Green, and Faythe Robinson for their technical assistance. We also thank Mark Berjanskii for assistance with statistical analysis of metabolomics data.
Conflict of interest
The author(s) declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The author(s) declared that Generative AI was not used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
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.2025.1721007/full#supplementary-material
References
1. Ahmad, MF, Ahmad, FA, Ashraf, SA, Saad, HH, Wahab, S, Khan, MI, et al. An updated knowledge of black seed (Nigella sativa Linn.): review of phytochemical constituents and pharmacological properties. J Herb Med. (2021) 25:100404. doi: 10.1016/j.hermed.2020.100404,
2. Cherif, M, Ben Salem, H, and Abidi, S. Effect of the addition of Nigella sativa seeds to low or high concentrate diets on intake, digestion, blood metabolites, growth, and carcass traits of Barbarine lamb. Small Rumin Res. (2018) 158:1–8. doi: 10.1016/j.smallrumres.2017.11.008,
3. Paarakh, PM. Nigella sativa Linn. –A comprehensive review. Indian J Nat Prod Resour. (2010) 1:409–29.
4. Mahmoud, AEM, and Bendary, MM. Effect of whole substitution of protein source by Nigella sativa meal and sesame seed meal in ration on the performance of growing lambs and calves. Glob Vet. (2014) 13:391–6.
5. Habeeb, AAM, and El Tarabany, AA. Effect of Nigella sativa or curcumin on daily body weight gain, feed intake and some physiological functions in growing Zaraibi goats during hot summer season. J Radiat Res Appl Sci. (2012) 5:60–78.
6. Ramdani, D, Juandita, KN, Hernaman, I, and Alhuur, KRG. Effects of dietary black cumin (Nigella sativa L.) meal on performance, blood metabolites, and digestibility in a rice straw-based diet of fattening Garut lambs. Vet World. (2024) 17:2152–8. doi: 10.14202/vetworld.2024.2152-2158,
7. Schneider-Stock, R, Fakhoury, IH, Zaki, AM, El-Baba, CO, and Gali-Muhtasib, HU. Thymoquinone: fifty years of success in the battle against cancer models. Drug Discov Today. (2014) 19:18–30. doi: 10.1016/j.drudis.2013.08.021,
8. Hamdan, AM, Al-Gayyar, MM, Shams, MEE, Alshaman, US, Prabahar, K, Bagalagel, A, et al. Thymoquinone therapy remediates elevated brain tissue inflammatory mediators induced by chronic administration of food preservatives. Sci Rep. (2019) 9:7026. doi: 10.1038/s41598-019-43568-x,
9. Albakry, Z, Karrar, E, Ahmed, IAM, Oz, E, Proestos, C, El Sheikha, AF, et al. Nutritional composition and volatile compounds of black cumin (Nigella sativa L.) seed, fatty acid composition and tocopherols, polyphenols, and antioxidant activity of its essential oil. Horticulturae. (2022) 8:575. doi: 10.3390/horticulturae8070575
10. Kabir, Y, Akasaka-Hashimoto, Y, and Komai, M. Volatile compounds of black cumin (Nigella sativa L.) seeds cultivated in Bangladesh and India. Heliyon. (2020) 6:e05343. doi: 10.1016/j.heliyon.2020.e05343,
11. Saneha, OR, Krishna, P, Das, S, Krishnakumar, IM, and Joseph, L. Exploring the possibility of a proprietary black cumin oil extract as a dual orexin receptor antagonist in restoring stress-sleep balance on stress-induced sleep deprived animals. Pharma Nutr. (2023) 26:100357
12. Forouzanfar, F., and Hosseinzadeh, H. (2020). “Protective role of Nigella sativa and thymoquinone in oxidative stress: a review,” in Nuts and seeds in health and disease prevention, ed. Preedy, V. R., and Watson, R. R. (Cambridge, MA: Academic Press), 127–146, doi: 10.1016/B978-0-12-818553-7.00011-5.
13. Velagapudi, R, Kumar, A, Bhatia, HS, El-Bakoush, A, Lepiarz, I, Fiebich, BI, et al. Inhibition of neuroinflammation by thymoquinone requires activation of Nrf2/ARE signalling. Int Immunopharmacol. (2017) 48:17–29. doi: 10.1016/j.intimp.2017.04.018,
14. Badr, G, Alwasel, S, Ebaid, H, Mohany, M, and Alhazza, I. Perinatal supplementation with thymoquinone improves diabetic complications and T cell immune responses in rat offspring. Cell Immunol. (2011) 267:133–40. doi: 10.1016/j.cellimm.2011.01.002,
15. Azab Awad-Allah, M. Effect of supplementation with niacin and Nigella sativa seeds on Friesian calves under heat stress conditions. J Anim Poult Prod. (2002) 27:791–801.
16. Kouakou, B, Gelaye, S, Kannan, G, Pringle, TD, and Amoah, EA. Blood metabolites, meat quality and muscle calpain-calpastatin activities in goats treated with low doses of recombinant bovine somatotropin. Small Rumin Res. (2005) 57:203–12. doi: 10.1016/j.smallrumres.2004.08.001
17. Batchu, P, Terrill, TH, Kouakou, B, Estrada-Reyes, ZM, and Kannan, G. Plasma metabolomic profiles as affected by diet and stress in Spanish goats. Sci Rep. (2021) 11:12607. doi: 10.1038/s41598-021-91893-x,
18. Kannan, G, Terrill, TH, Kouakou, B, Gazal, OS, Gelaye, S, Amoah, EA, et al. Transportation of goats: effects on physiological stress responses and live weight loss. J Anim Sci. (2000) 78:1450–7. doi: 10.2527/2000.7861450x,
19. Verbeke, WAJ, and Viaene, J. Ethical challenges for livestock production: meeting consumer concerns about meat safety and animal welfare. JAGE. (2000) 12:141–51.
20. Dalin, A-M, Magnusson, U, Häggendal, J, and Nyberg, L. The effect of transport stress on plasma levels of catecholamines, cortisol, corticosteroid-binding globulin, blood cell count, and lymphocyte proliferation in pigs. Acta Vet Scand. (1993) 34:59–68. doi: 10.1186/BF03548224,
21. Batchu, P, Naldurtiker, A, Kouakou, B, Terrill, TH, McCommon, GW, and Kannan, G. Metabolomic exploration of the effects of habituation to livestock trailer and extended transportation in goats. Front Mol Biosci. (2022) 9:1027069. doi: 10.3389/fmolb.2022.1027069,
22. Kannan, G, and Batchu, P. Omics approach for assessing welfare in sheep and goats: a focus on metabolomics In: G Kannan, editor. Small ruminant welfare, production and sustainability. Cambridge, MA: Academic Press (2025). 149–90. doi: 10.1016/B978-0-443-22201-6.00006-2
23. Zheng, J, Mandal, R, and Wishart, DS. A sensitive, high-throughput LC-MS/MS method for measuring catecholamines in low volume serum. Anal Chim Acta. (2018) 1037:159–16. doi: 10.1016/j.aca.2018.01.021,
24. Foroutan, A, Guo, AC, Vazquez-Fresno, R, Lipfert, M, Zhang, L, Zheng, J, et al. Chemical composition of commercial cow's milk. J Agric Food Chem. (2019) 67:4897–914. doi: 10.1021/acs.jafc.9b00204,
25. Haq, A, Abdullatif, M, Lobo, PI, Khabar, KS, Sheth, KV, and Al-Sedairy, ST. Nigella sativa: effect on human lymphocytes and polymorphonuclear leukocyte phagocytic activity. Immunopharmacology. (1995) 30:147–55. doi: 10.1016/0162-3109(95)00016-M,
26. Asaduzzaman Khan, M, Tania, M, Fu, S, and Fu, J. Thymoquinone, as an anticancer molecule: from basic research to clinical investigation. Oncotarget. (2017) 8:51907–19. doi: 10.18632/oncotarget.17206,
27. Nwe, TM, Hori, E, Manda, M, and Watanabe, S. Significance of catecholamines and cortisol levels in blood during transportation stress in goats. Small Rumin Res. (1996) 20:129–35. doi: 10.1016/0921-4488(95)00781-4
28. Kannan, G, Batchu, P, Naldurtiker, A, Dykes, GS, Gurrapu, P, Kouakou, B, et al. Habituation to livestock trailer and its influence on stress responses during transportation in goats. Animals. (2023) 13:1191. doi: 10.3390/ani13071191,
29. Gillam, LK, Palmer, JP, and Taborsky, GJ Jr. Tyramine-mediated activation of sympathetic nerves inhibits insulin secretion in humans. J Clin Endocrinol Metab. (2007) 92:4035–8. doi: 10.1210/jc.2007-0536,
30. Batchu, P, Hazard, T, Lee, JH, Terrill, TH, Kouakou, B, and Kannan, G. High-condensed tannin diet and transportation stress in goats: effects on physiological responses, gut microbial counts and meat quality. Animals. (2021) 11:2857. doi: 10.3390/ani11102857,
32. Purves, D, Brannon, EM, Cabeza, R, Huettel, SA, LaBar, KS, Platt, ML, et al. Principles of cognitive neuroscience. 2nd ed. Sunderland, MA: Sinauer Associates (2013).
33. David, J-C, Dairman, W, and Udenfriend, S. Decarboxylation to tyramine: a major route of tyrosine metabolism in mammals. Proc Natl Acad Sci USA. (1974) 71:1771–5. doi: 10.1073/pnas.71.5.1771,
34. Brosnan, JT. Interorgan amino acid transport and its regulation. J Nutr. (2003) 133:2068S–72S. doi: 10.1093/jn/133.6.2068S,
35. Soeters, MR, Serlie, MJ, Sauerwein, HP, Duran, M, Ruiter, JP, Kulik, W, et al. Characterization of D-3-hydroxybutyrylcarnitine (ketocarnitine): an identified ketosis-induced metabolite. Metabolism. (2012) 61:966–73. doi: 10.1016/j.metabol.2011.11.009,
36. Hack, A, Busch, V, Pascher, B, Busch, R, Bieger, I, Gempel, K, et al. Monitoring of ketogenic diet for carnitine metabolites by subcutaneous microdialysis. Pediatr Res. (2006) 60:93–6. doi: 10.1203/01.pdr.0000219479.95410.79,
37. Steinhauser, ML, Olenchock, BA, O’Keefe, J, Lun, M, Pierce, KA, Lee, H, et al. The circulating metabolome of human starvation. JCI Insight. (2018) 3:e121434. doi: 10.1172/jci.insight.121434,
38. Penner, GB, Steele, MA, Aschenbach, JR, and McBride, BW. Ruminant nutrition symposium: molecular adaptation of ruminal epithelia to highly fermentable diets. J Anim Sci. (2011) 89:1108–19. doi: 10.2527/jas.2010-3378,
39. Ramsay, RR, and Arduini, A. Carnitine acyltransferases and their role in modulating acyl-CoA pools. Arch Biochem Biophys. (1993) 302:307–14. doi: 10.1006/abbi.1993.1216,
40. Hiatt, WR, Regensteiner, JG, Wolfel, EE, Ruff, L, and Brass, EP. Carnitine and acylcarnitine metabolism during exercise in humans. Dependence on skeletal muscle metabolic state. J Clin Invest. (1989) 84:1167–73. doi: 10.1172/JCI114281,
41. Ghaffari, MH, Sadri, H, Schuh, K, Dusel, G, Prehn, C, Adamski, J, et al. Alterations of the acylcarnitine profiles in blood serum and in muscle from periparturient cows with normal or elevated body condition. J Dairy Sci. (2020) 103:4777–94. doi: 10.3168/jds.2019-17713,
42. Yang, Y, Sadri, H, Prehn, C, Adamski, J, Rehage, J, Danicke, S, et al. Acylcarnitine profiles in serum and muscle of dairy cows receiving conjugated linoleic acids or a control fat supplement during early lactation. J Dairy Sci. (2019) 102:754–67. doi: 10.3168/jds.2018-14685,
43. Koves, TR, Ussher, JR, Noland, RC, Slentz, D, Mosedale, M, Ilkayeva, O, et al. Mitochondrial overload and incomplete fatty acid oxidation contribute to skeletal muscle insulin resistance. Cell Metab. (2008) 7:45–56. doi: 10.1016/j.cmet.2007.10.013,
44. Gault, CR, Obeid, LM, and Hannun, YA. An overview of sphingolipid metabolism: from synthesis to breakdown. Adv Exp Med Biol. (2010) 688:1–23. doi: 10.1007/978-1-4419-6741-1_1,
45. Gulbins, E, Palmada, M, Reichel, M, Lüth, A, Böhmer, C, Amato, D, et al. Acid sphingomyelinase–ceramide system mediates effects of antidepressant drugs. Nat Med. (2013) 19:934–8. doi: 10.1038/nm.3214,
46. Allen, PJ. Creatine metabolism and psychiatric disorders: does creatine supplementation have therapeutic value? Neurosci Biobeh Rev. (2012) 36:1442–62. doi: 10.1016/j.neubiorev.2012.03.005,
47. Li, R, Wang, L, Chen, B, and Qi, P. Effects of transportation on blood indices, oxidative stress, rumen fermentation parameters and rumen microbiota in goats. Animals. (2024) 14:1616. doi: 10.3390/ani14111616,
48. Kaliber, M, Koluman, N, and Silanikove, N. Physiological and behavioral basis for the successful adaptation of goats to severe water restriction under hot environmental conditions. Animal. (2016) 10:82–8. doi: 10.1017/S1751731115001652,
49. Mehaba, N, Coloma-Garcia, W, Such, X, Caja, G, and Salama, AAK. Heat stress affects some physiological and productive variables and alters metabolism in dairy ewes. J Dairy Sci. (2020) 104:1099–110. doi: 10.3168/jds.2020-18943,
50. Al-Awadi, F, Fatania, H, and Shamte, U. The effect of a plant mixture extract on liver gluconeogenesis in streptozotocin induced diabetic rats. Diabetes Res. (1991) 18:163–8.
Keywords: antioxidant activity, black seed meal, catecholamines, metabolomics, stress
Citation: Gurrapu P, Batchu P, Shaik A, Terrill TH and Kannan G (2026) On-farm dietary supplementation of black seed (Nigella sativa) meal in goats: effects on physiological and metabolomic responses during transportation. Front. Vet. Sci. 12:1721007. doi: 10.3389/fvets.2025.1721007
Edited by:
Maria Giovanna Ciliberti, University of Foggia, ItalyReviewed by:
Chanadol Supapong, Rajamangala University of Technology Isan, Nakhon Ratchasima, ThailandChirasak Phoemchalard, Mahidol University, Amnat Charoen Campus, Thailand
Copyright © 2026 Gurrapu, Batchu, Shaik, Terrill and Kannan. 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: Govind Kannan, Z2tAYXVidXJuLmVkdQ==
Present address: Govind Kannan, Department of Poultry Science, Auburn University, Auburn, AL, United States
Priyanka Gurrapu