Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Vet. Sci., 22 January 2026

Sec. Animal Behavior and Welfare

Volume 12 - 2025 | https://doi.org/10.3389/fvets.2025.1721007

This article is part of the Research TopicAdvancing Animal Welfare Assessment: From Biomarkers to Smart Monitoring TechnologiesView all 10 articles

On-farm dietary supplementation of black seed (Nigella sativa) meal in goats: effects on physiological and metabolomic responses during transportation

Priyanka GurrapuPriyanka GurrapuPhaneendra BatchuPhaneendra BatchuArshad ShaikArshad ShaikThomas H. TerrillThomas H. TerrillGovind Kannan
&#x;Govind Kannan*
  • 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
www.frontiersin.org

Table 1. Ingredients and chemical composition of experimental diets (Black seed meal, BSM, diet; Control diet, C).

Figure 1
Flowchart illustrating the allocation of 80 goats. The goats are divided into two corrals of 40 each, labeled BSM group. Each corral's goats are then randomly allocated into two trailers, with compartments noted as front and rear. An example of time sampling is shown for both BSM group and control group.

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
Three line graphs showing the effects of transportation time on different antioxidant levels. A: Glutathione peroxidase levels peak at 4 hours, higher for C (orange) than BSM (blue), then decline. B: Superoxide dismutase levels decrease over time, with slight fluctuations, BSM generally higher. C: TBARS levels stay low until 10 hours, then rise, showing a higher increase for C. Error bars represent SEM. Baseline values are provided.

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
Four graphs labeled A to D show immunological responses over transportation time. Graph A is a bar chart comparing Immunoglobulin A in BSM and C, with C higher. Graph B is a line chart showing Immunoglobulin A changes over time, with notable fluctuations at 2, 4, 10, and 16 hours. Graph C depicts Immunoglobulin G variations, peaking at 10 hours. Graph D shows a decrease in Interleukin 6 levels over time. BSM and C data are represented by lines with error bars indicating SEM.

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.

Figure 4
Heatmap displaying concentration levels of various compounds over time intervals (0, 2, 4, 10, 16 hours). Compounds include S-Methoxytryptamine, Epinephrine, Metanephrine, Dopamine, and others. Color scale ranges from low (yellow) to high (red) concentration.

Figure 4. Heatmap of plasma catecholamines and derivatives clustered by transportation time.

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
www.frontiersin.org

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
www.frontiersin.org

Table 3. Acylcarnitines significantly (P < 0.05) affected by treatment (BSM = Black seed meal; C = Control; n = 40 goats/TRT) in goats.

Table 4
www.frontiersin.org

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
www.frontiersin.org

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 metabolite significance versus fold change, with red and blue dots indicating significant increases or decreases, respectively. Metabolites are labeled, including LysoPC, SM, alanine, tryptophan, and methylmalonic acid.

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 6
Heatmap showing metabolite abundance across samples labeled BSM and C. Rows represent metabolites; columns are samples. Colors range from yellow (low) to red (high) abundance. A legend indicates intensity levels.

Figure 6. Heatmap of significantly affected (p < 0.05) plasma metabolites clustered by treatment.

Figure 7
Scatter plot showing two overlapping clusters represented by ellipses. One cluster, labeled BSM, is in red, centered around component values -0.5, 0. The other, labeled C, is in green, centered around axis 0.5, 0. Each axis represents different components, with percentages indicating explained variance: Component 1 (16.9%) and Component 2 (14.4%). Points are distributed within each ellipse.

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
Scatter plot displaying VIP scores for various metabolites on the vertical axis, ranging from Phenylalanine to lysophosphatidylcholine C17:0. Scores vary from 1.4 to 2.4, represented by black dots. A color-coded sidebar indicates levels from low (blue) to high (red) with labels BSM and C.

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
Line graphs show the concentrations of five substances over time: phenylalanine (A), asparagine (B), glutamine (C), glutamic acid (D), and pyruvic acid (E). Each graph compares two conditions, BSM and C, across transportation times of 0, 2, 4, 10, and 16 hours. The graphs include errors bars for standard error of the mean (SEM). Statistical significance at various points is indicated by asterisks. The initial baselines for each graph are provided: phenylalanine BSM 72.1 ± 6.35, C 75.6 ± 5.75; asparagine BSM 52.9 ± 5.69, C 47.5 ± 5.94; glutamine BSM 331.9 ± 27.50, C 288.9 ± 30.84; glutamic acid BSM 105.9 ± 9.27, C 108.4 ± 9.95; pyruvic acid BSM 55.3 ± 2.37, C 55.6 ± 1.73.

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
Network diagram of metabolic pathways represented as labeled nodes. Each node is connected by lines indicating relationships among pathways, with varying node sizes and colors highlighting the importance and involvement of each metabolism. Key pathways include alanine, aspartate, and glutamate metabolism, and arginine biosynthesis.

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,

PubMed Abstract | Crossref Full Text | Google Scholar

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,

PubMed Abstract | Crossref Full Text | Google Scholar

3. Paarakh, PM. Nigella sativa Linn. –A comprehensive review. Indian J Nat Prod Resour. (2010) 1:409–29.

Google Scholar

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.

Google Scholar

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.

Google Scholar

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,

PubMed Abstract | Crossref Full Text | Google Scholar

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,

PubMed Abstract | Crossref Full Text | Google Scholar

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,

PubMed Abstract | Crossref Full Text | Google Scholar

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

Crossref Full Text | Google Scholar

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,

PubMed Abstract | Crossref Full Text | Google Scholar

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

Google Scholar

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.

Crossref Full Text | Google Scholar

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,

PubMed Abstract | Crossref Full Text | Google Scholar

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,

PubMed Abstract | Crossref Full Text | Google Scholar

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.

Google Scholar

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

Crossref Full Text | Google Scholar

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,

PubMed Abstract | Crossref Full Text | Google Scholar

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,

PubMed Abstract | Crossref Full Text | Google Scholar

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.

Google Scholar

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,

PubMed Abstract | Crossref Full Text | Google Scholar

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,

PubMed Abstract | Crossref Full Text | Google Scholar

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

Crossref Full Text | Google Scholar

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,

PubMed Abstract | Crossref Full Text | Google Scholar

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,

PubMed Abstract | Crossref Full Text | Google Scholar

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,

PubMed Abstract | Crossref Full Text | Google Scholar

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,

PubMed Abstract | Crossref Full Text | Google Scholar

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

Crossref Full Text | Google Scholar

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,

PubMed Abstract | Crossref Full Text | Google Scholar

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,

PubMed Abstract | Crossref Full Text | Google Scholar

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,

PubMed Abstract | Crossref Full Text | Google Scholar

31. Bolander, FF. Molecular Endocrinology. 3rd ed. Cambridge, MA: Academic Press (2004).

Google Scholar

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).

Google Scholar

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,

PubMed Abstract | Crossref Full Text | Google Scholar

34. Brosnan, JT. Interorgan amino acid transport and its regulation. J Nutr. (2003) 133:2068S–72S. doi: 10.1093/jn/133.6.2068S,

PubMed Abstract | Crossref Full Text | Google Scholar

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,

PubMed Abstract | Crossref Full Text | Google Scholar

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,

PubMed Abstract | Crossref Full Text | Google Scholar

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,

PubMed Abstract | Crossref Full Text | Google Scholar

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,

PubMed Abstract | Crossref Full Text | Google Scholar

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,

PubMed Abstract | Crossref Full Text | Google Scholar

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,

PubMed Abstract | Crossref Full Text | Google Scholar

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,

PubMed Abstract | Crossref Full Text | Google Scholar

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,

PubMed Abstract | Crossref Full Text | Google Scholar

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,

PubMed Abstract | Crossref Full Text | Google Scholar

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,

PubMed Abstract | Crossref Full Text | Google Scholar

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,

PubMed Abstract | Crossref Full Text | Google Scholar

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,

PubMed Abstract | Crossref Full Text | Google Scholar

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,

PubMed Abstract | Crossref Full Text | Google Scholar

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,

PubMed Abstract | Crossref Full Text | Google Scholar

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,

PubMed Abstract | Crossref Full Text | Google Scholar

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.

PubMed Abstract | Google Scholar

51. Kundrapu, S, and Noguez, J. Laboratory assessment of anemia. Adv Clin Chem. (2018) 83:197–225. doi: 10.1016/bs.acc.2017.10.006,

PubMed Abstract | Crossref Full Text | Google Scholar

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

Received: 08 October 2025; Revised: 13 December 2025; Accepted: 16 December 2025;
Published: 22 January 2026.

Edited by:

Maria Giovanna Ciliberti, University of Foggia, Italy

Reviewed by:

Chanadol Supapong, Rajamangala University of Technology Isan, Nakhon Ratchasima, Thailand
Chirasak 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

Disclaimer: 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.