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

ORIGINAL RESEARCH article

Front. Microbiol., 29 January 2026

Sec. Food Microbiology

Volume 16 - 2025 | https://doi.org/10.3389/fmicb.2025.1712904

Effects of different doses of oregano essential oil on growth performance, health parameters, and the gut microbiome-metabolome profile in heat-stressed Pinan cattle (Piedmontese × Nanyang)

Jiangge WangJiangge Wang1Lei LiuLei Liu1Jiashun ShengJiashun Sheng1Liyang ZhangLiyang Zhang1Qiaozhen LiQiaozhen Li2Tengyun Gao
Tengyun Gao1*Xian Liu
Xian Liu3*
  • 1Henan International Joint Laboratory of Nutrition Regulation and Ecological Raising of Domestic Animal, College of Animal Science and Technology, Henan Agricultural University, Zhengzhou, China
  • 2Animal Husbandry Bureau of Xinye County, Xinye, China
  • 3Henan Animal Husbandry Technology Extension Station, Zhengzhou, China

Introduction: This study investigated the dose–response efficacy of dietary oregano essential oil (OEO) in mitigating severe heat stress (THI ≈ 86) in beef cattle.

Methods: Thirty-six Pinan bulls were fed a basal diet alone (control) or supplemented with 7 (L-OEO) or 14 g/d (H-OEO) of OEO for 60 days.

Results: The low-dose OEO (7 g/d) significantly improved hepatic function (reduced ALT, AST; increased ALB, TP), enhanced immune (increased IgA, IgM) and antioxidant status (decreased MDA, increased SOD, CAT), and increased the ruminal abundance of Bacteroidota and Prevotella, which correlated negatively with acyl-glycine metabolites. In contrast, the high dose (14 g/d) only increased GSH-Px and T4, resulted in higher MDA than L-OEO, and did not significantly affect the rumen microbiota.

Conclusion: Supplementation with 7 g/d OEO optimally improved overall health and metabolic function in heat-stressed bulls, whereas a 14 g/d dose offered no additional benefits.

1 Introduction

Summer heat stress poses a significant challenge to global beef production. Beef cattle are particularly susceptible due to the metabolic heat generated from ruminal fermentation and their insulating hair coats (Nakamura and Morrison, 2010). When the temperature-humidity index (THI) exceeds a critical threshold, cattle exhibit characteristic signs of heat stress, including reduced feed intake and impaired ruminal function (Mader et al., 2006). To maintain thermal homeostasis, physiological responses are activated, such as the suppression of the hypothalamic feeding center to reduce appetite and minimize endogenous heat production from feeding and rumination (Morrison et al., 2014; Bernabucci et al., 2010).

This systemic stressor profoundly disrupts the gastrointestinal ecosystem. Heat stress induces ruminal and intestinal dysbiosis through a combination of elevated core body temperature, reduced feed intake, and altered water consumption (Bernabucci et al., 2010; Zhao et al., 2019; Park et al., 2022). This shift often depletes beneficial cellulolytic bacteria while enriching potential pathogens and lactate producers, compromising volatile fatty acid profiles, fiber digestibility, and gut barrier integrity (Bernabucci et al., 2010; He et al., 2021). The resulting endotoxin translocation and inflammatory response exacerbate systemic metabolic dysregulation.

As the central metabolic organ, the liver is severely impaired under heat stress. It faces the dual challenge of maintaining gluconeogenesis while counteracting oxidative damage (Bernabucci et al., 2010; Abbas et al., 2020). Persistent oxidative stress inhibits fatty acid oxidation and export, leading to hepatic lipidosis (Skibiel et al., 2018). Concurrently, mitochondrial dysfunction increases reactive oxygen species (ROS) production, while key antioxidant enzymes (e.g., SOD, GSH-Px, CAT) are depleted, culminating in hepatocellular injury (Wang et al., 2020; Eslamizad et al., 2024).

Given this multifaceted pathophysiology, natural interventions with multi-target effects are needed. Oregano essential oil (OEO), rich in bioactive phenolics like thymol and carvacrol, shows promise in mitigating heat stress in livestock (Garcia-Galicia et al., 2020; Nazar et al., 2019; Yousefi et al., 2024). Its potential mechanisms include modulating gut microbiota, alleviating oxidative stress, and regulating immune responses (Cui et al., 2024; Habibi et al., 2015; Terenina et al., 2011; Zeng et al., 2015). However, its efficacy and mechanisms in heat-stressed beef cattle require systematic elucidation.

This knowledge gap is particularly relevant for the Pinan cattle (Piedmontese × Nanyang), a breed renowned in central China for its rapid growth and high meat quality (Wei et al., 2022; Bo et al., 2024). Despite its economic importance, the physiological responses of this newly developed breed to heat stress and the efficacy of mitigants like OEO remain unexplored.

Therefore, this study provides the first comprehensive evaluation of OEO supplementation in Pinan cattle. We employed an integrated multi-omics approach, moving beyond conventional growth and antioxidant assessments to unravel the complex host-microbiota-metabolite interactions. Our goal was to determine the dose-dependent effects of OEO on heat-stress-impaired pathways and provide evidence-based recommendations for enhancing resilience in this unique breed.

2 Materials and methods

2.1 Experimental design and animal management

The experiment was conducted from June to August 2024. Thirty-six healthy 12-month-old Pinan (Piedmontese × Nanyang) bulls with uniform initial body weight (380 ± 25 kg) were obtained from the Xinye County Hengyong Beef Cattle Breeding Cooperative, Henan Province (32°17′N, 110°58′E). In a completely randomized design, the cattle were assigned to three dietary treatments (n = 12 cattle/treatment). The bulls were housed in a well-ventilated barn equipped with individual stalls, allowing for strict individual management, feeding, and supplementation. This design ensured that each animal received its precise dietary treatment without cross-contamination and that all measurements could be unequivocally attributed to the individual animal. Treatments were: (Nakamura and Morrison, 2010) CON, basal diet without OEO; (Mader et al., 2006) LOEO, basal diet + 7 g OEO/animal/day; (Morrison et al., 2014) HOEO, basal diet + 14 g OEO/animal/day. Following a 7-day dietary adaptation, all bulls entered a 60-day experimental period during which OEO supplementation was initiated and maintained as per treatment allocation.

Cattle were offered a total mixed ration (TMR) twice daily at 06:00 and 18:00 h. The daily OEO dose for each bull was first evenly mixed with a small quantity of the TMR (carrier), and then this premix was thoroughly incorporated into the individual’s entire 18:00 TMR portion using a mechanical mixer to ensure homogeneous distribution. This fixed feeding regimen (10 kg per feeding, 20 kg daily) was implemented to standardize nutrient and OEO intake. Feed refusals were recorded daily and were negligible (<5%), confirming that the actual OEO intake matched the intended dosage.

The oregano essential oil (OEO) used in this study was procured from Feidi Technology (Product Standard Code: Q/HFD0103-2023). Chemical characterization by gas chromatography–mass spectrometry (GC-MS) established that the oil was overwhelmingly dominated by phenolic compounds. The two primary bioactive phenols, thymol, and carvacrol, were identified as the major constituents, constituting 24.97% and 64.96% of the total composition, respectively. Their combined relative abundance of 89.93% definitively classifies this OEO as the carvacrol-thymol chemotype, which is renowned for its potent biological activities. Other notable components included D-limonene (6.49%) and 2,4-di-tert-butylphenol (1.79%), while the remaining identified compounds each accounted for less than 0.5% of the total profile.

The feed composition and nutritional components are shown in Table 1.

Table 1
www.frontiersin.org

Table 1. Feed ingredient composition and nutritional content.

Given that the animals were individually stall-housed, fed a fixed amount, and supplemented precisely on an individual basis, the experimental unit for all subsequent statistical analyses was unequivocally defined as the individual animal.

2.2 Field Indicator measurement

2.2.1 Environmental parameter measurements

At the daily heat-peak (14:00 h), ambient temperature (T, °C) and relative humidity (RH, %) were recorded at 1.5 m above the barn floor using a handheld digital thermo-hygrometer (TASI TA620A; accuracy ±0.5 °C and ±3% RH). The temperature–humidity index (THI) was calculated according to ASABE (citation number):

Thi=0.8T+RH(T14.4)+46.4

Heat-stress severity was classified as THI ≤ 74 (no stress), 75–78 (mild), 79–83 (moderate), and ≥ 84 (severe) (Gaughan et al., 2008).

2.2.2 Growth performance

On days 1 and 60, body weight was recorded at 06:00 following a 12-h fast (water allowed). Average daily gain (ADG) and feed conversion ratio (FCR) were calculated as:

ADG(kgd1)=(BWfinalBWinitial)/60
FCR=daily feed intake/ADG

2.3 Sample collection and processing

Blood samples (5 mL) were collected from the tail vein into heparinized vacuum tubes. After 30 min clotting at 25 °C, samples were centrifuged (3,000 × g, 15 min, 131 4 °C). Serum was aliquoted (1 mL) into sterile cryovials, snap-frozen on dry ice (−78 °C), and stored at −40 °C pending untargeted metabolomics and biochemical analyses.

On the final day before morning feeding, ~100 mL rumen fluid was collected via an orogastric tube from each animal. Fluid was immediately filtered through four layers of sterile cheesecloth, and pH was recorded with a calibrated portable pH meter. Aliquots (~2 mL) of filtrate were snap-frozen in liquid nitrogen and stored at −80 °C pending 16S rRNA sequencing of the rumen microbiota.

Fresh fecal samples were obtained by digital rectal stimulation using sterile long-arm. In review gloves, snap-frozen in liquid nitrogen, and stored at −80 °C pending 16S rRNA sequencing of the gut microbiota.

2.4 Sample analysis

2.4.1 Serum biochemical parameter assays

Serum concentrations of immunoglobulin A (IgA), immunoglobulin M (IgM), immunoglobulin G (IgG), total antioxidant capacity (T-AOC), total superoxide dismutase (T-SOD), glutathione peroxidase (GSH-Px), malondialdehyde (MDA), catalase (CAT), heat shock protein 70 (HSP70), cortisol (COR), triiodothyronine (T3), and thyroxine (T4) were determined using corresponding commercial assay kits (Nanjing Jiancheng Bioengineering Institute, China) in accordance with the manufacturer’s instructions.

The concentrations of alanine aminotransferase (ALT), aspartate aminotransferase (AST), alkaline phosphatase (ALP), urea (UREA), albumin (ALB), total protein (TP), glucose (GLU), triglycerides (TG), and free hemoglobin (Hb) were measured using a fully automated biochemical analyzer (Mindray BS-330E, Shenzhen Mindray Bio-Medical Electronics Co., Ltd., China).

2.4.2 Serum untargeted metabolomics analysis

Serum untargeted metabolomics was performed by Personalbio Biotechnology Co., Ltd. (Shanghai, China). Ultra-high performance liquid chromatography coupled with high-resolution mass spectrometry (UHPLC-HRMS) analysis was performed using a Thermo Vanquish Flex UHPLC system (Thermo Fisher Scientific, USA) coupled to a Thermo Orbitrap Exploris 120 mass spectrometer (Thermo Fisher Scientific, USA).

Serum metabolites were extracted using a standardized protocol. Briefly, a 50 μL aliquot of each serum sample was mixed with 200 μL of pre-cooled methanol:acetonitrile (1:1, v/v). The mixture was vortexed for 30 s and then incubated at −20 °C for 30 min to precipitate proteins. After centrifugation at 12,000 × g and 4 °C for 10 min, 200 μL of the supernatant was collected and completely dried under a vacuum concentrator. The dried metabolite extract was reconstituted in 150 μL of 50% methanol containing 5 ppm L-2-chlorophenylalanine as an internal standard. The solution was vortexed for 30 s, centrifuged again under the same conditions, and the final supernatant was filtered through a 0.22 μm membrane prior to UHPLC-MS analysis.

A pooled quality control (QC) sample was generated by combining equal volumes of the filtrate from all individual samples. This QC sample was processed identically to the analytical samples and was injected repeatedly throughout the analytical sequence to monitor instrument stability and data quality.

Chromatographic separation was performed on an ACQUITY UPLC HSS T3 column (100Å, 1.8 μm, 2.1 mm × 100 mm; Waters, USA) maintained at 40 °C. The mobile phase flow rate was 0.4 mL/min, and the injection volume was 2 μL. The metabolites were detected using a Thermo Orbitrap Exploris 120 mass spectrometer (Thermo Fisher Scientific, USA) controlled by Xcalibur software (v4.7). The mass spectrometer was equipped with a HESI source and operated in both positive and negative ionization modes with a spray voltage of 3.5 kV and −3.0 kV, respectively. The sheath gas and auxiliary gas were set to 40 arb and 10 arb, respectively. The capillary temperature was 320 °C, and the auxiliary gas heater temperature was 300 °C. Full MS scans were acquired from 70 to 1,000 m/z at a resolution of 60,000. Data-dependent acquisition (DDA) mode was used for MS/MS fragmentation of the top 4 most intense ions at a resolution of 15,000.

Subsequent statistical analyses—including univariate analysis, multivariate analysis, and correlation analysis of differential metabolites—were conducted using the GenesCloud platform.1

2.4.3 Rumen and gut microbiota analysis

Rumen fluid and fecal samples were shipped on dry ice to Personalbio Biotechnology In review Co., Ltd. (Shanghai, China) for 16S rRNA gene sequencing. After thawing and homogenization, total genomic DNA was extracted with the OMEGA Soil DNA Kit (M5635-02, Omega Bio-Tek, USA) following the manufacturer’s instructions. Extracted DNA was checked for integrity and fragment size on a 0.8% agarose gel, and concentration and purity were determined with a NanoDrop 2000/2000c spectrophotometer (Thermo Fisher Scientific, USA). High-quality DNA was used as template for PCR amplification of the V3–V4 hypervariable region of the bacterial 16S rRNA gene with primers 338F (5′-ACTCCTACGGGAGGCAGCA-3′) and 806R (5′-GGACTACHVGGGTWTCTAAT-3′). Amplicons were purified with AMPure XP beads and used to construct sequencing libraries. Libraries were sequenced on an Illumina NovaSeq 6000 platform (Illumina, USA) generating 2 × 150 bp paired-end reads. Raw sequencing reads were processed and analyzed as described below.

2.4.4 Correlation analysis

Pearson correlation analysis was conducted to explore the relationships between differentially abundant gut/rumen microbes and differential serum metabolites. For each comparison, both the correlation coefficient (r) and the corresponding p-value are reported. To account for the multiple comparisons inherent in this high-dimensional data, the False Discovery Rate was controlled using the Benjamini–Hochberg procedure, and the resulting FDR-adjusted p-values (q-values) are reported. A correlation was considered statistically significant if the q-value was <0.05.

Furthermore, to aid in the biological interpretation of the magnitude of the correlations, the absolute value of the correlation coefficient r was used to define the strength of the relationship as follows: r ≥ 0.60: Strong correlation, 0.30 ≤ r < 0.60: Moderate correlation, r < 0.30: Weak correlation. Only correlations that were both statistically significant (q < 0.05) and at least of moderate strength (r ≥ 0.30) are discussed in the main text.

2.5 Data analysis and statistics

All data obtained from individual animals were analyzed using one-way analysis of variance (ANOVA) in SPSS software (Version 26.0). The individual animal was considered the experimental unit (n = 12 per treatment). The statistical model included the fixed effect of the dietary treatment (CON, LOEO, HOEO). When the ANOVA indicated a significant overall treatment effect (p < 0.05), differences between individual treatment means were further evaluated using the Least Significant Difference (LSD) post hoc test. Data are expressed as means ± standard deviation (SD) of the individual measurements. A probability value of p < 0.05 was considered statistically significant.

The individual animal was unequivocally considered the experimental unit for all statistical analyses (n = 12 per treatment), as treatments were applied, and all data (feed intake, blood parameters, microbial communities, etc.) were collected and measured on an individual basis.

3 Results

3.1 Effects of essential oil supplementation on growth performance of heat-stressed Pinan cattle

The Temperature–Humidity Index (THI), respiratory rate, and growth performance are summarized in Tables 2, 3. Throughout the trial, the mean THI exceeded 84 in all groups, indicating persistent severe heat stress. Respiratory rate did not differ among groups (p > 0.05). Likewise, initial body weight, final body weight, and average daily gain (ADG) were not significantly different (p > 0.05).

Table 2
www.frontiersin.org

Table 2. Temperature, relative humidity, THI, and respiratory rate of each group of cattle house.

Table 3
www.frontiersin.org

Table 3. Effect of OEO on the ADG of Pinan cattle under heat stress.

3.2 Effects of OEO on serum biochemical parameters in heat-stressed Pinan cattle

As summarized in Table 4, 60 days of OEO supplementation significantly improved markers of liver injury, protein metabolism, and synthetic function. Compared with CON, both OEO-treated groups exhibited significantly lower serum alanine aminotransferase (ALT, p < 0.05) and aspartate aminotransferase (AST, p < 0.01) concentrations. Albumin (ALB) concentration was significantly higher in the low-dose group than in both the high-dose group (p < 0.01) and the control (p < 0.01). The high-dose group also exhibited significantly higher ALB than the control (p < 0.05). Total protein (TP) was also markedly elevated in L-OEO and H-OEO (p < 0.01). Furthermore, urea (UREA) concentration was significantly reduced in L-OEO (p < 0.05) and H-OEO (p < 0.01) relative to CON.

Table 4
www.frontiersin.org

Table 4. Effect of OEO on the serum biochemical indicators of Pinan cattle under heat stress.

Regarding immune function, OEO supplementation significantly enhanced early immune responses. Serum IgA concentrations were significantly higher in both L-OEO and H-OEO compared to CON (p < 0.01), and IgM concentrations were significantly higher than those in CON (p < 0.05).

In terms of antioxidant capacity, H-OEO showed a significant increase in GSH-Px activity relative to both L-OEO and CON (p < 0.01). SOD activity was significantly higher in both OEO-supplemented groups than in the control (p < 0.05). CAT activity in L-OEO was significantly higher than in CON (p < 0.05), whereas no significant difference was observed between H-OEO and either L-OEO or CON (p > 0.05).

L-OEO showed the lowest MDA, differing from both CON (p < 0.01) and H-OEO (p < 0.05), while H-OEO was also lower than CON (p < 0.05). For other serum biochemical indicators, the hemoglobin (Hb) concentration in H-OEO was significantly lower than that in CON (p < 0.05). In contrast, the thyroxine (T4) concentration in H-OEO was significantly higher than in CON (p < 0.05). For both Hb and T4, no significant difference was detected between H-OEO and L-OEO (p > 0.05).

3.3 Serum untargeted metabolomics

The untargeted metabolomics analysis detected a total of 31,747 metabolite features in the serum samples, with 18,816 features identified in positive ion mode and 12,931 features in negative ion mode. While a significant proportion of these features could not be unequivocally assigned to known metabolic pathways in current databases, they represent a comprehensive snapshot of the systemic metabolic state. The identified metabolites were dominated by features consistent with lipid and phospholipid species (e.g., various phosphatidylcholines detected primarily in positive ion mode), amino acid derivatives, and other organic acids. This extensive coverage of the serum metabolome provides a robust basis for detecting relative changes in metabolite abundance between experimental groups.

3.3.1 Principal component analysis (PCA)

As shown in Figure 1, under positive ion mode, a clear separation in metabolic profiles was observed between both OEO-supplemented groups (L-OEO and H-OEO) and the control group (CON), whereas the metabolic profiles of L-OEO and H-OEO were not clearly separated from each other. In contrast, no obvious separation was detected among any of the groups under negative ion mode.

Figure 1
Two principal component analysis (PCA) plots labeled

Figure 1. PCA analysis of serum metabolites.

3.3.2 Differential metabolite screening

Differential metabolites were identified based on the criteria of p < 0.05 and fold change (FC) > 1.5 or FC < 0.58, and were visualized using volcano plots (Figure 2). The results indicated that the differential metabolites between L-OEO and H-OEO consisted primarily of various phosphatidylcholines (PCs). Comparisons between L-OEO and CON revealed significant alterations in metabolites including multiple PCs, isobutyrylglycine, N-acetylvaline, hydroxyphenylacetylglycine, 2-methylbutyrylglycine (2-MBG), thymol sulfate, and docosahexaenoic acid ethyl ester. Meanwhile, differential metabolites between H-OEO and CON included isobutyrylglycine, N-acetylvaline, hydroxyphenylacetylglycine, thymol sulfate, glycodeoxycholic acid, taurochenodeoxycholic acid, and cholylalanine.

Figure 2
Six scatter plots labelled (a) to (f) compare log2 fold change (FC) against negative log10 P-values. Each plot illustrates differential expression with dots: red for upregulated, blue for no difference, yellow for downregulated. Vertical and horizontal dashed lines denote significance thresholds. Plots compare conditions such as L-OEO vs H-OEO, L-OEO vs CON, H-OEO vs CON, with both positive and negative expressions analyzed.

Figure 2. Volcano plot for screening differential metabolites in the serum of heat stressed Pinang cattle with different doses of OEO. (a) Positive ion mode: L-OEO vs. H-OEO. (b) Positive ion mode: L-OEO vs. CON. (c) Positive ion mode: H-OEO vs. CON. (d) Negative ion mode: L-OEO vs. H-OEO. (e) Negative ion mode: L-OEO vs. CON. (f) Negative ion mode: H-OEO vs. CON.

To gain a systemic overview of the metabolic perturbations, KEGG pathway enrichment analysis was performed on the differential metabolites identified in all pairwise comparisons (L-OEO vs. CON, H-OEO vs. CON, and L-OEO vs. H-OEO). Notably, no pathways met the threshold for statistical significance (FDR-corrected p-value < 0.05) in any comparison. This indicates that the metabolic response to OEO supplementation was not characterized by the strong enrichment or depletion of any single canonical pathway. Instead, the effects appear to be distributed across a broad network of interconnected metabolic processes, as reflected in the diverse classes of metabolites altered (e.g., acylglycines, bile acids, and phospholipids).

3.4 16S microbiomics

3.4.1 Microbial composition

The effects of OEO supplementation on the ruminal and gut microbiota are summarized in Tables 58. OEO administration exerted distinct modulatory effects on the microbial communities in the rumen and intestine.

Table 5
www.frontiersin.org

Table 5. Effects of OEO on the abundance of phylum-level gut microbiota in heat-stressed Pinan cattle.

Table 6
www.frontiersin.org

Table 6. Effects of OEO on the abundance of phylum-level rumen microbiota in heat-stressed Pinan cattle.

Table 7
www.frontiersin.org

Table 7. Effects of OEO on the abundance of genus-level gut microbiota in heat-stressed Pinan cattle.

Table 8
www.frontiersin.org

Table 8. Effects of OEO on the abundance of genus-level rumen microbiota in heat-stressed Pinan cattle.

At the phylum level, the dominant bacterial phyla in the gut were Firmicutes_A and Bacteroidota, which collectively accounted for over 90% of the total bacterial abundance. The relative abundance of Verrucomicrobiota was significantly lower in both OEO-supplemented groups (L-OEO and H-OEO) than in the CON (p < 0.05). No significant differences were observed among the groups for the remaining phyla (p > 0.05).

In the rumen, the predominant bacterial phyla were similarly Bacteroidota and Firmicutes_A. Compared to CON, L-OEO exhibited a significant increase in the abundance of Bacteroidota (p < 0.05) and a highly significant reduction in Firmicutes_A (p < 0.01). In contrast, no significant differences in microbial composition were observed between H-OEO and the CON (p > 0.05). Moreover, the abundance of Fibrobacterota was significantly higher in L-OEO than in both H-OEO and CON.

At the genus level, Faecalicoccus and Cryptobacteroides were the dominant genera in the gut across all groups. The abundance of Cryptobacteroides was significantly higher in L-OEO than in H-OEO and C (p < 0.05).

In the rumen, the predominant genera were Prevotella and Cryptobacteroides. Compared with H-OEO and CON, L-OEO showed a significantly higher abundance of Prevotella (p < 0.01). Conversely, the abundance of Cryptobacteroides was markedly lower in L-OEO than in both H-OEO and CON (p < 0.01). Additionally, the abundance of Sodaliphilus was significantly lower in L-OEO than in CON (p < 0.05), while no significant difference was observed between L-OEO and H-OEO (p > 0.05).

3.4.2 Microbial diversity

Analysis of α-diversity indices (Chao1 and Shannon) revealed no significant differences among the three treatment groups in either the ruminal or gut microbiota (all p > 0.05; Figure 3), indicating that OEO supplementation did not alter the overall richness or evenness of the microbial communities.

Figure 3
Two panels display box plots comparing microbial diversity indexes across different groups. Panel (a) shows gut groups with Chao1, Simpson, Shannon, and Observed Species indexes. Panel (b) displays rumen groups with the same indexes. Each plot includes p-values and color-coded groups labeled as L-OEO, H-OEO, and CON.

Figure 3. Alpha diversity analysis of gut and ruminal microbiota. (a) Alpha diversity analysis of gut microbiota. (b) Alpha diversity analysis of ruminal microbiota.

Beta diversity of the gut and rumen microbiota was assessed using non-metric multidimensional scaling (NMDS). As illustrated in Figure 4, no distinct clustering pattern was observed among the groups in the gut microbial communities. In contrast, clear separation was evident in the rumen microbiota, with samples from L-OEO and H-OEO forming clusters distinct from those of CON.

Figure 4
Panel a shows a non-metric multidimensional scaling (NMDS) plot with three groups: L-OEO (Gut), H-OEO (Gut), and CON (Gut), represented by distinct ellipses. Panel b displays a LefSe analysis bar chart of different taxa with an LDA threshold of two, indicating group-specific bacterial abundance. Panel c illustrates another NMDS plot for L-OEO (Rumen), H-OEO (Rumen), and CON (Rumen), with overlapping ellipses. Panel d shows a corresponding LefSe analysis bar chart for different rumen-associated taxa with varying LDA scores by group.

Figure 4. NMDS analysis and LDA score plots of gut and rumen microbiota. (a) NMDS of gut microbiota. (b) LEfSe results for differential gut microbial taxa. Bar length represents log10-transformed LDA score. (c) NMDS of rumen microbiota. (d) LEfSe results for differential rumen bacterial taxa. Bar length represents log10-transformed LDA score.

Linear discriminant analysis (LDA) revealed that in the gut, the genus Cryptobacteroides was significantly enriched in the lL-OEO. No characteristic microbial taxa were identified as significantly enriched in the H-OEO, while the phylum Verrucomicrobiota was significantly enriched in the cCON.

In the rumen, low-dose OEO supplementation resulted in significant enrichment of Bacteroidaceae and Prevotella. In contrast, high-dose OEO supplementation was associated with significant enrichment of Lachnospirales and Butyrivibrio.

To infer whether the observed compositional shifts in the microbiota translated to changes in their metabolic potential, we performed phylogenetic investigation of communities by reconstruction of unobserved states (PICRUSt2) analysis. The predicted abundances of Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were compared across the three dietary groups. No KEGG pathways at Level 3 were found to be significantly different between any of the group comparisons (CON vs. L-OEO, CON vs. H-OEO, or L-OEO vs. H-OEO) after false discovery rate (FDR) correction for multiple comparisons (all q > 0.05). This indicates that while OEO supplementation altered the relative abundance of specific bacterial taxa, it did not induce a major restructuring of the core predicted functional potential of the microbial community at the resolution of this analysis.

3.5 Integrated analysis

Pearson correlation analysis was conducted to explore the relationships between differentially abundant gut/rumen microbes and serum metabolites across the experimental groups. As illustrated in Figure 5, within the gut microbiota, serum concentrations of phosphatidylcholine (PC), hydroxyphenylacetylglycine, and isobutyrylglycine showed significant negative correlations (p < 0.05) with the abundance of Cryptobacteroides (gut) when comparing L-OEO and CON In the rumen, a significant negative correlation (p < 0.05) was identified between serum PC levels and the abundance of Prevotella (rumen) in both L-OEO and H-OEO. Furthermore, in the comparison between L-OEO and CON, serum levels of PC, N-acetylvaline, hydroxyphenylacetylglycine, and isobutyrylglycine were significantly negatively correlated (p < 0.05) with Prevotella (rumen) abundance, but significantly positively correlated (p < 0.05) with the abundance of Cryptobacteroides (rumen) and Saccharofermentans (rumen). Conversely, serum docosahexaenoic acid ethyl ester concentration was positively correlated (p < 0.05) with Prevotella (rumen) and negatively correlated (p < 0.05) with Saccharofermentans (rumen). Between H-OEO and CON, serum thymol sulfate exhibited a significant positive correlation (p < 0.05) with the abundance of Sodaliphilus (rumen), while serum N-acetylvaline was significantly negatively correlated (p < 0.05) with Prevotella (rumen).

Figure 5
Heatmaps labeled a to f display clustered data of various compounds along the vertical axis, with conditions or samples along the horizontal axis. The color scale ranges from blue to red, representing values from negative zero point five to positive zero point five. Each heatmap is slightly distinct in layout and data distribution, indicating different experimental or analytical conditions.

Figure 5. Correlation analysis of differential gut and rumen microbes with differential serum metabolites. (a) Correlation analysis of differential gut microbes and serum metabolites between L-OEO and H-OEO. (b) Correlation analysis of differential rumen microbes and serum metabolites between L-OEO and L-OEO. (c) Correlation analysis of differential gut microbes and serum metabolites between L-OEO and CON. (d) Correlation analysis of differential rumen microbes and serum metabolites between L-OEO and CON. (e) Correlation analysis of differential gut microbes and serum metabolites between H-OEO and CON. (f) Correlation analysis of differential rumen microbes and serum metabolites between H-OEO and CON.

To mitigate the risk of false positives from multiple comparisons, we performed FDR correction (Benjamini–Hochberg method) on all Pearson correlation analyses between microbes and metabolites. After rigorous correction, robust microbe-metabolite associations were identified in both the intestine and rumen. The most prominent association network emerged from comparing the rumen between the L-OEO and CON groups, where key bacterial genera (e.g., Prevotella, Cryptobacteroides, UBA4334) showed extremely strong correlations with multiple phospholipid and amino acid metabolites (all r > 0.75, q < 0.01). All significant associations are detailed in Table 9.

Table 9
www.frontiersin.org

Table 9. Significant microbiota-metabolite correlations after FDR correction.

4 Discussion

Heat stress reduces feed intake, increases respiratory rate, and disrupts ruminal In review fermentation in Pinan cattle, collectively leading to a negative energy balance. Concurrently, high ambient temperatures promote lipid peroxidation of the inner mitochondrial membrane and enhance electron leakage from Complexes I and III, thereby increasing superoxide anion (O₂•) production. In addition, suppressed ATP synthase activity lowers energy production efficiency. These metabolic alterations collectively trigger excessive reactive oxygen species (ROS) generation (Khan et al., 2021), resulting in widespread oxidative damage.

Under heat stress, reduced feed intake is normally the primary constraint on weight gain. By fixing daily intake at 20 kg TMR per head we removed this constraint; consequently, the numerical rise in ADG (L-OEO 0.29 H-OEO 0.34 CON 0.51 kg/d) reflects a true improvement in the efficiency with which dietary nutrients were converted to live-weight. It should be noted that the L-OEO began the trial significantly heavier than both the H-OEO and CON. Because initial body weight is a key determinant of subsequent growth potential, heavier animals typically exhibit a flatter growth curve and devote a larger fraction of energy to maintenance. Thus, the fact that L-OEO still achieved a higher ADG than CON, despite its higher starting weight and associated maintenance burden, underscores the positive biological impact of the low-dose essential-oil supplement.

Analysis of serum biochemical parameters indicated that OEO supplementation significantly reduced markers of liver injury (ALT and AST) and protein catabolism (UREA), while enhancing indicators of synthetic function (ALB and TP). These results suggest that OEO alleviates heat stress-induced hepatic damage by protecting hepatocytes and improving protein metabolism (Cui et al., 2024; Song et al., 2014; Zhou et al., 2020). Notably, the L-OEO showed greater efficacy in promoting protein synthesis, whereas the H-OEO resulted in a more pronounced reduction in blood urea nitrogen.

Consistent with the results of the present study, in which OEO improved liver function in heat-stressed cattle, Zhang et al. (2024) reported that dietary OEO supplementation in sows also exerted significant hepatoprotective effects by enhancing serum antioxidant capacity (significantly increasing CAT and GSH), reducing the level of the liver-specific enzyme AST, and optimizing the structure of the gut microbiota.

In terms of immune response, OEO supplementation significantly elevated early immune mediators (IgM) and mucosal immunity-related antibodies (IgA), but did not significantly affect longer-term immunoglobulins (IgG). These observations are consistent with the established immune shift under heat stress, which suppresses Th1-mediated cellular immunity while promoting Th2-driven humoral responses. This suggests that OEO may counter heat-induced immunosuppression primarily through enhancement of the humoral immune pathway (Bernabucci et al., 2010; Gupta et al., 2022).

Consistent with the findings of Luo et al. (2024), dietary supplementation with oregano essential oil under summer feeding conditions significantly increased serum immunoglobulin (IgA and IgM) levels in calves by modulating the ruminal microbiota, suggesting that oregano essential oil can effectively enhance humoral immunity in ruminants under heat stress.

Regarding antioxidant capacity, the L-OEO exhibited a significant reduction in the lipid peroxidation product MDA and an increase in CAT activity, whereas the H-OEO demonstrated a pronounced enhancement in GSH-Px activity. This dose-dependent variation may reflect the concentration-related bioactivity of carvacrol and thymol, the primary active components of OEO (Nehme et al., 2021). The simultaneous increase in SOD activity further supports the synergistic mitigation of oxidative stress, although overall antioxidant capacity did not improve significantly. In review the absence of a significant change in COR levels suggests that OEO may help attenuate overactivation of the hypothalamic–pituitary–adrenal (HPA) axis.

The significant increase in T4 observed in the H-OEO suggests the involvement of thyroid hormone in the metabolic regulation mediated by OEO supplementation. Meanwhile, the decrease in free Hb indicates enhanced erythrocyte membrane stability (Das et al., 2016), which is consistent with the overall improvement in antioxidant status.

These findings indicate that OEO effectively mitigates liver injury, improves metabolic function, enhances immune responses, and bolsters antioxidant capacity in heat-stressed Pinan cattle. The benefits are likely mediated by thymol and carvacrol, the principal bioactive of OEO, which donate electrons to neutralize superoxide anions, providing direct antioxidant protection and upregulating key enzymes (GSH-Px, SOD, CAT). Consequently, the endogenous antioxidant defense system is comprehensively reinforced. The divergent responses between low- and high-dose supplementation suggest that lower dosage may be optimal for specific physiological endpoints.

Notably, the improvements in antioxidant and immune parameters occurred in the absence of significant changes in respiratory rate, a conventional clinical indicator of heat stress. This apparent discrepancy can be explained by the distinct levels of biological organization at which these measurements are taken and the primary mechanisms of action of OEO. Respiratory rate is a whole-animal, integrative physiological response coordinated by the hypothalamic thermoregulatory center to maximize evaporative heat loss (West, 2003). Its modulation is often a refractory endpoint that may not respond until core body temperature is significantly altered. In contrast, the antioxidant and immune parameters measured (e.g., SOD, CAT, IgA) reflect the cellular and biochemical milieu, which are more directly influenced by dietary bioactive compounds like thymol and carvacrol.

The observed divergence suggests that the beneficial effects of OEO supplementation in heat-stressed cattle may operate primarily at the sub-physiological level. By directly scavenging ROS and potentially upregulating the expression of endogenous antioxidant enzymes, OEO mitigates oxidative damage at its source within tissues. Similarly, its immunomodulatory effects are likely mediated through direct interactions with immune cells and the gut-associated lymphoid tissue. This intracellular and molecular improvement can occur independently of, and potentially prior to, changes in the animal’s effort to dissipate heat. Therefore, the enhancement of the antioxidant and immune systems represents a more sensitive and early indicator of OEO’s efficacy, underscoring its role in fortifying the host’s internal defense mechanisms against the cellular assaults of heat stress, even when the external stressor (high THI) persists and the overall physiological heat-dissipation response remains fully engaged.

Previous studies have shown that heat stress increases insulin, which in turn suppresses lipolysis and curtails lipid mobilization (Abbas et al., 2020; Ma et al., 2021; Delgado-Martin et al., 2024). To compensate for the ensuing energy deficit, muscle protein is catabolized to release branched-chain amino acids (BCAAs) such as valine and isoleucine. These BCAAs are subsequently metabolized in the liver by branched-chain amino acid transaminase (BCAT) and branched-chain α-keto acid dehydrogenase (BCKDH). This process drives deamination pathways—leucine → isovaleryl-CoA → 3-methylcrotonyl-CoA → acetyl-CoA, and valine → succinyl-CoA—thereby furnishing substrates for the tricarboxylic acid (TCA) cycle and sustaining energy production.

In the present study, BCAA-derived acylglycines—N-acetylvaline, isobutyrylglycine, and 2-methylbutyrylglycine (2-MBG)—together with hydroxyphenylacetylglycine, 3-methyl-thymine (3-MT), and thymol sulfate, were all significantly elevated. The accumulation of acylglycines points to impaired hepatocellular detoxification, most likely owing to oxidative inhibition of glycine N-acyltransferase (GLYAT); this inhibition mirrors both liver dysfunction and accelerated muscle proteolysis (Guo et al., 2020; Alves et al., 2019). Elevated 2-MBG further signals disrupted isoleucine catabolism (Newgard et al., 2009), which may curtail TCA-cycle flux and ATP synthesis. Detection of thymol sulfate, a phase-II conjugate of thymol, provides an indirect indicator of the parent compound’s in-vivo antioxidant and anti-inflammatory activity (Zhou et al., 2020; Yu et al., 2016; Placha et al., 2022).

These alterations underscore a metabolic vicious cycle under heat stress: BCKDH inhibition blocks BCAA oxidation to succinyl-CoA, curtailing TCA-cycle entry and causing intermediate accumulation (Sun et al., 2016). Simultaneously, oxidative stress depletes cysteine, the rate-limiting precursor of glutathione (GSH) (Ulrich and Jakob, 2019; Nguyen et al., 2023), perpetuating the In review cycle: energy deficit drives BCAA catabolism; BCAA depletion impairs GSH synthesis; weakened antioxidant defense exacerbates mitochondrial damage and ROS generation; and further compromise of energy metabolism ensues (Sharma et al., 2024; Circu and Aw, 2010; Kim et al., 2015; Bifari and Nisoli, 2017).

In the H-OEO, bile acid-related compounds—glycodeoxycholic acid, taurochenodeoxycholic acid and cholylalanine—were also markedly up-regulated alongside the BCAA-derived metabolites. Glycodeoxycholic acid, a secondary bile acid, controls lipid metabolism and bile acid homeostasis via the FXR–FGF19 pathway (Meessen et al., 2025), but can exert pro-inflammatory effects in settings such as acute pancreatitis (Fauzi et al., 2022). Taurochenodeoxycholic acid displays anti-inflammatory and immunomodulatory activity through up-regulation of SRSF9 and potentiation of glucocorticoid receptor (GR) signaling (Bao et al., 2021). Cholylalanine, a conjugated bile acid, may modulate nuclear receptor activity and reshape the gut microbiota (Zhao et al., 2025).

The KEGG enrichment analysis, while not identifying any single significantly altered pathway, provides an important systemic perspective. The absence of a starkly dysregulated pathway suggests that OEO does not act like a targeted pharmaceutical agent that blocks a specific enzymatic reaction. Instead, its supplementation in heat-stressed cattle facilitates a global rebalancing of metabolism that cuts across multiple pathway boundaries.

This finding is entirely consistent with the known pleiotropic nature of plant essential oils. Rather than overwhelming a single pathway, the components of OEO, particularly thymol and carvacrol, likely exert subtle modulatory effects on multiple cellular targets—including membrane integrity, enzyme activity, and signal transduction—which collectively manifest as the observed improvements in liver health, antioxidant status, and immune function (Walasek-Janusz et al., 2024; Cui et al., 2024). The most sensitive biomarkers of this system-wide rebalancing were the specific metabolite changes we detected, such as the accumulation of acylglycines signaling a relief of hepatic detoxification pressure and the presence of thymol sulfate confirming the direct bioactivity of OEO constituents. Thus, the mechanism of action is best described as a “systems-level stabilization” rather than a “pathway-specific intervention”.

Collectively, our data show that oregano essential oil (OEO) markedly reshapes metabolism in heat-stressed Pinan cattle. By accelerating branched-chain amino acid (BCAA) catabolism, the supplement partially offsets the energy deficit; however, concomitant inhibition of branched-chain α-keto acid dehydrogenase (BCKDH) causes intermediate metabolites to accumulate and exacerbates oxidative injury. Moreover, high-dose OEO activates bile-acid metabolism, potentially providing an alternative energy source during thermal stress (Van Passel et al., 2011).

The rumen and gut microbial communities in ruminants form a complex and highly integrated ecosystem that plays a vital role in nutrient metabolism, immune defense, developmental regulation, and environmental adaptation. In the present study, at the phylum level, only Verrucomicrobiota showed a significant decrease in abundance in the gut of treated animals. Akkermansia muciniphila, a representative species of this phylum, utilizes mucin as its sole carbon and nitrogen source by secreting glycoside hydrolases, sulfatases, and carbohydrate esterases to degrade mucin (Van Passel et al., 2011). The decline in its abundance is likely attributable to oregano-essential-oil-induced disintegration of the cell envelope, as thymol and carvacrol insert into lipid bilayers and destabilise membrane integrity, consequently reducing Verrucomicrobiota populations.

In the rumen, supplementation enriched Bacteroidota and Fibrobacterota while depressing Firmicutes_A. Bacteroidota drive polysaccharide degradation (cellulose, starch) (Hao et al., 2021), and their expansion likely improves dietary energy harvest. Fibrobacterota, the primary cellulolytic rumen taxon (Yan et al., 2024), directly boosts forage conversion. Firmicutes_A harbours numerous glycoside hydrolases and sugar transporters (e.g., Lachnospiraceae) that cleave complex polysaccharides and import monosaccharides; consequently, its decline may lower intestinal carbohydrate uptake efficiency (Mohamed, 2025).

At the genus level, only low-dose OEO increased gut Cryptobacteroides. Members of this genus encode a broad suite of carbohydrate-active enzymes (CAZymes) that target complex hemicelluloses (Lin et al., 2024), implying improved energy harvest from the diet. In review in the rumen, the same dose elevated Prevotella while depressing Sodaliphilus and Cryptobacteroides. Although Prevotella hydrolyses cellulose and protein, its overgrowth can raise ruminal ammonia nitrogen (Chiquette et al., 2008). These dose-dependent shifts underscore the fine-tuning of microbial function and host metabolic homeostasis by OEO.

Both Cryptobacteroides and Sodaliphilus extensive polysaccharide-utilization loci (PULs) that cleave complex polysaccharides into simple sugars, boosting fiber digestion (Lin et al., 2024). Their decline likely reflects direct membrane disruption by OEO bioactive—thymol and carvacrol—which compromise cell-wall integrity.

A key finding of our microbiome analysis is the dissociation between the significant compositional shifts (β-diversity and specific taxa) and the stability observed in both α-diversity and predicted metagenomic functions (PICRUSt2). The absence of α-diversity changes demonstrates that OEO, at the doses tested, acts as a precision modulator of the gut ecosystem rather than a broad-spectrum disruptor. It selectively enriches beneficial taxa (e.g., Prevotella) without compromising the overall ecological stability.

The specific enrichment of Bacteroidaceae and Prevotella under low-dose OEO supplementation can be attributed to a combination of direct and indirect mechanisms. Firstly, certain members of the Bacteroidetes phylum, particularly Prevotella, possess efflux pump systems and membrane adaptations that may confer a degree of intrinsic tolerance to low concentrations of plant secondary compounds, allowing them to thrive when more sensitive taxa are mildly inhibited (Song et al., 2014; Zhang et al., 2021).

Secondly, and perhaps more importantly, the primary bioactive components of OEO, thymol and carvacrol, are known to disrupt the membrane of Gram-negative bacteria. At sub-inhibitory (low) concentrations, this membrane perturbation may not be lethal but can instead act as an environmental cue or stressor. This can trigger a shift in microbial metabolism towards energy storage and competitive growth. Prevotella and other Bacteroidetes are adept at rapidly uptaking and fermenting soluble sugars and starches. The mild stress from low-dose OEO might channel microbial metabolism towards the efficient degradation of readily available polysaccharides, a niche where these bacteria excel.

Indirectly, by creating a mildly selective environment, low-dose OEO likely reduces competition from OEO-sensitive bacteria (e.g., some Firmicutes), freeing up ecological space and resources. This allows resilient and metabolically versatile taxa like Prevotella to dominate. This enriched community is highly proficient at degrading dietary polysaccharides, leading to increased propionate production, which aligns with our observed improvements in host energy metabolism.

In contrast, the lack of additional benefits and the onset of potential adverse effects at the high OEO dose (14 g/d) can be mechanistically explained by a shift from a selective, stimulatory effect to a broad-spectrum antimicrobial one. The ‘hormetic effect’—where a low dose is beneficial while a high dose is detrimental—is well-documented for plant essential oils (Zhou et al., 2019; Hristov et al., 2013).

At high concentrations, the phenolic compounds thymol and carvacrol exert their full antimicrobial power by: 1. Integrating into and severely disrupting microbial membranes, causing irreversible damage, leakage of cellular contents, and cell death. This non-selective action does not spare beneficial, functionally important bacteria. 2. Quenching microbial quorum-sensing signals, thereby disrupting the microbial cross-talk necessary for a stable and cooperative ecosystem.

This broad-spectrum suppression is evidenced in our study by the absence of characteristic beneficial taxa enrichment in the H-OEO and the minimal impact on ruminal community structure. The collapse of the delicate symbiotic network, particularly the suppression of key fibrolytic bacteria, could impair the digestion of structural carbohydrates (e.g., cellulose and hemicellulose), potentially explaining the numerical (though not significant) reduction in energy-producing bacteria.

Furthermore, the high dose may have direct negative effects on the host. The rumen epithelium itself can be irritated by high concentrations of potent phenols, potentially compromising its absorptive and barrier functions. Additionally, an excessive dose may overwhelm the host’s liver detoxification pathways, increasing the metabolic burden rather than alleviating it. This shift from a precision-guided modulation at low dose to an indiscriminate suppression at high dose provides a clear mechanistic rationale for the observed “threshold effect,” where doubling the dosage failed to yield additional benefits and instead showed a tendency to suppress some positive physiological responses.

Microbial-diversity analysis revealed high resilience of the gut microbiota to OEO, whereas the rumen community responded in a dose-dependent manner: low-dose OEO significantly improved ruminal structure, but high-dose had minimal impact. Linear discriminant analysis (LDA) identified gut Cryptobacteroides as significantly enriched under low-dose OEO, likely enhancing energy harvest through degradation of complex polysaccharides (Sidney et al., 2024). No characteristic taxa were detected in the H-OEO, consistent with a broad-spectrum inhibitory effect or microbial adaptive resistance at elevated OEO concentrations (Zhou et al., 2020). In controls, Verrucomicrobiota was markedly enriched; members of this phylum encode GH29 (fucosidase), GH141 (rhamnosidase) and sulfotransferases that specialize in degrading sulfated polysaccharides such as fucoidan (Orellana et al., 2022).

Low-dose OEO supplementation significantly enriched Bacteroidaceae and Prevotella in the rumen, suggesting enhanced polysaccharide degradation and propionate production. This metabolic shift may compete for H₂, thereby reducing the substrate availability for methanogens and potentially mitigating methane emissions. In contrast, high-dose OEO enriched Lachnospirales and Butyrivibrio, which are core taxa involved in hemicellulose-pectin degradation and butyrate production via PUL-enzyme systems and the butyrate kinase pathway. These metabolic activities not only provide energy to the host but also contribute to milk fat synthesis (Wirth et al., 2018; Palevich et al., 2019). The ruminal microbiota in the CON was predominantly composed of Firmicutes and Ruminococcaceae, which specialize in degrading lignin and cellulose. This process produces H₂ and butyrate, thereby supplying substrates to methanogens (Malik et al., 2023).

The lack of significant PICRUSt2 results, despite clear taxonomic changes, suggests that the primary mechanism behind the improved host physiology may not be a major rewiring of microbial catabolic potential. Instead, the benefits are likely mediated through alternative pathways: firstly, via direct host effects of the absorbed bioactive OEO compounds (thymol and carvacrol), which our serum biochemistry and metabolomics strongly support; and secondly, through potentially fine-tuned changes in microbial gene expression or metabolic output that are not captured by a static genetic potential prediction. This holistic view underscores that the health benefits of phytogenic feed additives can be achieved without destabilizing the functional core of the gut microbiota.

In summary, OEO induced only minor gut-microbiota shifts: low-dose enrichment of Cryptobacteroides improved energy harvest, whereas the high dose lacked discernible benefits, plausibly owing to broad-spectrum inhibition or microbial resistance. In the rumen, a clear dose–response prevailed: the low dose selectively expanded polysaccharide-degrading and propionate-producing taxa, boosting energy yield and curtailing methanogenesis; the high dose favored fiber degraders and elevated butyrate. These results underscore the need for precise dosing when using OEO to refine rumen fermentation and mitigate methane emissions in cattle.

In addition to altering the composition of intestinal and rumen microbial communities, our analysis of microbiota–metabolite correlations, following strict FDR correction, revealed the mechanisms by which these microbial changes influence host metabolic status.

In the intestinal context, a significant positive correlation was observed between Paraprevotella and PC [16:0/18:1(9Z)] in both the L-OEO and H-OEO. This suggests that Paraprevotella may participate in or promote this phospholipid metabolic pathway (Yao et al., 2022), indicating that L-OEO and H-OEO treatments may influence host phospholipid metabolism by modulating Paraprevotella, with a more pronounced effect at higher concentrations.

In the L-OEO and CON, Cryptobacteroides showed a significant negative correlation with several phospholipids (PCs). This implies that Cryptobacteroides may inhibit the synthesis or promote the degradation of these phospholipids, suggesting that L-OEO treatment may influence phospholipid metabolism by downregulating Cryptobacteroides. These phospholipids may be associated with cell membrane structure or signal transduction (O’Donnell, 2022).

Moreover, UBA737 was significantly negatively correlated with Isobutyrylglycine, suggesting that UBA737 may inhibit the metabolism of Isobutyrylglycine. RF16 exhibited significant negative correlations with 3-Methyl-thymine and PC [20:1(11Z)/16:0], indicating that RF16 may be involved in the negative regulation of pyrimidine metabolism and phospholipid metabolism. Thus, L-OEO treatment may affect DNA synthesis and cell membrane function by modulating RF16 (Liu et al., 2019; Li et al., 2024).

In the H-OEO and CON, CAG-41 was significantly negatively correlated with Glycodeoxycholic acid, suggesting that CAG-41 may inhibit the metabolism of glycodeoxycholic acid. CAG-41belongs to Firmicutes_A, many members of which are known to possess bile salt hydrolase (BSH) activity and are capable of deconjugating bile acids, thereby altering bile acid composition and concentration (Ridlon et al., 2006). This indicates that H-OEO treatment may influence host bile acid metabolism by modulating CAG-41, subsequently affecting lipid absorption or intestinal barrier function.

In the rumen, within the L-OEO and H-OEO, Prevotella showed a significant negative correlation with 1-Oleoyl-2-Stearoyl-sn-Glycero-3-Phosphocholine. Prevotella bivia has been experimentally confirmed to secrete PLA₂ or other phospholipases capable of hydrolyzing PC into lysophosphatidylcholine and free fatty acids (Mikamo et al., 1998). It is inferred that Prevotella may inhibit the synthesis or promote the degradation of this phospholipid, suggesting that both L-OEO and H-OEO treatments might influence ruminal phospholipid metabolism by modulating Prevotella, with the higher dose exhibiting a more pronounced inhibitory effect.

In the L-OEO and CON, Prevotella was significantly negatively correlated with 1-stearoyl-2-oleoyl-sn-glycero-3-phosphoserine (r = −0.922, q = 0) and Thymol Sulfate (r = −0.889, q = 0). A mucosal desulfating sulfatase (mdsA) gene was cloned and expressed from Prevotella strain RS2 isolated from pig colon, and its product can remove sulfate from sugars or mucins (Wright et al., 2000). It is hypothesized that Prevotella may inhibit the metabolism of Thymol Sulfate, and L-OEO treatment might affect rumen membrane lipid structure and sulfate detoxification processes by downregulating Prevotella.

Cryptobacteroides was significantly positively correlated with PC [20:1(11Z)/16:0] (r = 0.756, q = 0.001) and 1-stearoyl-2-oleoyl-sn-glycero-3-phosphoserine (r = 0.846, q = 0.001). Cryptobacteroides may synthesize choline, a precursor for PC synthesis (Broad and Dawson, 1976; Sommese et al., 2012), suggesting it could promote the synthesis or accumulation of these phospholipids. This indicates that L-OEO treatment might influence ruminal phospholipid metabolism by upregulating Cryptobacteroides – a finding contrary to the intestinal results, highlighting site-specific effects.

Succiniclasticum was significantly positively correlated with Docosahexaenoic acid ethyl ester (DHA-EE) (r = 0.811, q = 0.003). Succiniclasticum may be involved in the synthesis or metabolism of DHA-EE, suggesting that L-OEO treatment could affect rumen unsaturated fatty acid metabolism (Carreño et al., 2019) by modulating Succiniclasticum, thereby influencing host lipid nutrition.

UBA4334 showed significant negative correlations with Hydroxyphenylacetylglycine, 1-stearoyl-2-oleoyl-sn-glycero-3-phosphoserine, Isobutyrylglycine, N-Acetylvaline, and 2-Methylbutyrylglycine. UBA4334 might inhibit the metabolism of these amino acids and phospholipids. L-OEO treatment could potentially impact ruminal amino acid metabolism and membrane lipid synthesis by downregulating UBA4334. Further isolation and genomics-centered characterization of the UBA4334 strain are necessary to validate this causal relationship.

WRMH01 was significantly positively correlated with 3-Methyl-thymine (r = 0.686, q = 0.028), 1-stearoyl-2-oleoyl-sn-glycero-3-phosphoserine (r = 0.931, q = 0), Hydroxyphenylacetylglycine (r = 0.808, q = 0), and Thymol Sulfate (r = 0.893, q = 0). WRMH01 may promote pyrimidine, phospholipid, and sulfate metabolism. L-OEO treatment might influence ruminal DNA synthesis, membrane lipid structure, and detoxification processes by upregulating WRMH01. Further isolation and identification of the WRMH01 strain are required to confirm the causal link.

In the H-OEO and CON, Sodaliphilus was significantly positively correlated with Cholylalanine (r = 0.754, q = 0.018). However, no publicly available studies have directly reported that Sodaliphilus can promote or participate in Cholylalanine metabolism, indicating that further research is still needed.

In summary, both L-OEO and H-OEO influence phospholipid, bile acid, amino acid, and fatty acid metabolism in both the intestine and rumen by modulating specific microorganisms (such as Paraprevotella, Cryptobacteroides, and Prevotella). H-OEO demonstrates superior efficacy in regulating phospholipid metabolism in the intestine and exhibits stronger inhibition of phospholipid degradation in the rumen. In contrast, L-OEO exerts opposite effects in the rumen compared to the intestine (e.g., in the regulation of Cryptobacteroides), highlighting site-specific effects. Further studies are needed to validate the causal relationships for certain microbiota–metabolite associations (e.g., those involving UBA4334 and WRMH01).

While this study provides clear evidence for the dose-dependent effects of OEO in heat-stressed cattle, certain limitations should be considered. Firstly, although the individual stall-feeding design was ideal for precisely controlling OEO intake, the number of biological replicates (n = 12 per treatment), while standard for such studies, inherently limits the statistical power to detect more subtle effects. Future studies with a larger sample size would be valuable to confirm and extend our findings. Secondly, the fixed feeding regimen (20 kg/head/day), while crucial for standardizing nutrient and OEO intake, means that the potential interplay between OEO supplementation and voluntary feed intake under heat stress could not be assessed. Investigations involving ad libitum feeding are needed to explore this interaction and its implications for practical production settings.

5 Conclusion

This study provides the first systematic evaluation of the dose-dependent effects of oregano essential oil (OEO) on Pinan cattle under sustained severe heat stress. The results demonstrate that a low dose of 7 g/animal/day optimally enhances liver function, antioxidant capacity, and immune response, while enriching beneficial ruminal taxa such as Prevotella and Bacteroidota to support energy harvest. Notably, low-dose OEO also differentially regulates phospholipid and amino acid metabolism across intestinal and ruminal environments, reflecting site-specific microbial modulation. In contrast, doubling the dosage to 14 g/animal/day not only failed to yield additional benefits but also induced broad-spectrum microbial inhibition and potential disruption of host metabolic homeostasis, revealing a clear “threshold effect.” Therefore, 7 g/animal/day is recommended as a precise and effective supplementation level of OEO for alleviating heat stress in Pinan cattle. Future studies should validate its economic and environmental potential in production settings under ad libitum feeding and explore synergistic models combining OEO with fermentable energy sources.

Data availability statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found at: Z2FvdGVuZ3l1bkBoZW5hdS5lZHUuY24=.

Ethics statement

The animal study was approved by Henan International Joint Laboratory of Nutrition Regulation and Ecological Raising of Domestic Animal. The study was conducted in accordance with the local legislation and institutional requirements.

Author contributions

JW: Data curation, Formal analysis, Visualization, Writing – original draft. LL: Investigation, Writing – review & editing. JS: Investigation, Writing – review & editing. LZ: Conceptualization, Methodology, Supervision, Writing – review & editing. QL: Conceptualization, Funding acquisition, Project administration, Supervision, Visualization, Writing – review & editing. TG: Conceptualization, Funding acquisition, Methodology, Supervision, Writing – review & editing. XL: Funding acquisition, Project administration, Resources, Supervision, Writing – review & editing.

Funding

The author(s) declared that financial support was received for this work and/or its publication. This study was funded by the Henan Beef Cattle Industrial Technology System (HARS-22-13-T) and the National Key Research and Development Program of China - Project 4 (2018YFD0501704).

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.

Footnotes

References

Abbas, Z., Sammad, A., Hu, L., Fang, H., Xu, Q., and Wang, Y. (2020). Glucose metabolism and dynamics of facilitative glucose transporters (GLUTs) under the influence of heat stress in dairy cattle. Meta 10:312. doi: 10.3390/metabo10080312,

PubMed Abstract | Crossref Full Text | Google Scholar

Alves, A., Bassot, A., Bulteau, A.-L., Pirola, L., and Morio, B. (2019). Glycine metabolism and its alterations in obesity and metabolic diseases. Nutrients 11:1356. doi: 10.3390/nu11061356,

PubMed Abstract | Crossref Full Text | Google Scholar

Bao, L., Hao, D., Wang, X., He, X., Mao, W., and Li, P. (2021). Transcriptome investigation of anti-inflammation and immuno-regulation mechanism of taurochenodeoxycholic acid. BMC Pharmacol. Toxicol. 22:23. doi: 10.1186/s40360-021-00491-0,

PubMed Abstract | Crossref Full Text | Google Scholar

Bernabucci, U., Lacetera, N., Baumgard, L. H., Rhoads, R. P., Ronchi, B., and Nardone, A. (2010). Metabolic and hormonal acclimation to heat stress in domesticated ruminants. Animal 4, 1167–1183. doi: 10.1017/S175173111000090X,

PubMed Abstract | Crossref Full Text | Google Scholar

Bifari, F., and Nisoli, E. (2017). Branched-chain amino acids differently modulate catabolic and anabolic states in mammals: a pharmacological point of view. Br. J. Pharmacol. 174, 1366–1377. doi: 10.1111/bph.13624,

PubMed Abstract | Crossref Full Text | Google Scholar

Bo, D., Feng, Y., Bai, Y., Li, J., Wang, Y., You, Z., et al. (2024). Whole-genome resequencing reveals genetic diversity and growth trait-related genes in Pinan cattle. Animals 14:2163. doi: 10.3390/ani14152163,

PubMed Abstract | Crossref Full Text | Google Scholar

Broad, T. E., and Dawson, R. M. C. (1976). Role of choline in the nutrition of the rumen protozoon Entodinium caudatum. J. Gen. Microbiol. 92, 391–397. doi: 10.1099/00221287-92-2-391,

PubMed Abstract | Crossref Full Text | Google Scholar

Carreño, D., Toral, P. G., Pinloche, E., Belenguer, A., Yáñez-Ruiz, D. R., Hervás, G., et al. (2019). Rumen bacterial community responses to DPA, EPA and DHA in cattle and sheep: a comparative in vitro study. Sci. Rep. 9:11857. doi: 10.1038/s41598-019-48294-y,

PubMed Abstract | Crossref Full Text | Google Scholar

Chiquette, J., Allison, M. J., and Rasmussen, M. A. (2008). Prevotella bryantii 25A used as a probiotic in early-lactation dairy cows: effect on ruminal fermentation characteristics, Milk production, and Milk composition. J. Dairy Sci. 91, 3536–3543. doi: 10.3168/jds.2007-0849,

PubMed Abstract | Crossref Full Text | Google Scholar

Circu, M. L., and Aw, T. Y. (2010). Reactive oxygen species, cellular redox systems, and apoptosis. Free Radic. Biol. Med. 48, 749–762. doi: 10.1016/j.freeradbiomed.2009.12.022,

PubMed Abstract | Crossref Full Text | Google Scholar

Cui, H., Zhang, C., Su, K., Fan, T., Chen, L., Yang, Z., et al. (2024). Oregano essential oil in livestock and veterinary medicine. Animals 14:1532. doi: 10.3390/ani14111532,

PubMed Abstract | Crossref Full Text | Google Scholar

Das, R., Sailo, L., Verma, N., Bharti, P., Saikia, J., Imtiwati,, et al. (2016). Impact of heat stress on health and performance of dairy animals: a review. Vet World 9, 260–268. doi: 10.14202/vetworld.2016.260-268

Crossref Full Text | Google Scholar

Delgado-Martin, M., Zhang, Q., and Kazak, L. (2024). BAT-tling oxidative stress through BCAA catabolism. Life Metab. 3:loae028. doi: 10.1093/lifemeta/loae028,

PubMed Abstract | Crossref Full Text | Google Scholar

Eslamizad, M., Albrecht, D., Kuhla, B., and Koch, F. (2024). Cellular and mitochondrial adaptation mechanisms in the colon of lactating dairy cows during hyperthermia. J. Dairy Sci. 107, 3292–3305. doi: 10.3168/jds.2023-24004,

PubMed Abstract | Crossref Full Text | Google Scholar

Fauzi, A., Suhendro, S. M., Rani, A. A., Sajuthi, D., Permanawati, A. R., and Makmun, D. (2022). Role of glycodeoxycholic acid to induce acute pancreatitis in Macaca nemestrina. J. Med. Primatol. 51, 134–142. doi: 10.1111/jmp.12577,

PubMed Abstract | Crossref Full Text | Google Scholar

Garcia-Galicia, I. A., Arras-Acosta, J. A., Huerta-Jimenez, M., Rentería-Monterrubio, A. L., Loya-Olguin, J. L., Carrillo-Lopez, L. M., et al. (2020). Natural oregano essential oil may replace antibiotics in lamb diets: effects on meat quality. Antibiotics 9:248. doi: 10.3390/antibiotics9050248,

PubMed Abstract | Crossref Full Text | Google Scholar

Gaughan, J. B., Mader, T. L., Holt, S. M., and Lisle, (2008). A new heat load index for feedlot cattle. Journal of Animal Science, 86, 226–234. doi: 10.2527/jas.2007-0305,

PubMed Abstract | Crossref Full Text | Google Scholar

Guo, Y. S., Tao, J. Z., Xu, L. H., Wei, F. H., and He, S. H. (2020). Identification of disordered metabolic networks in postpartum dairy cows with left displacement of the abomasum through integrated metabolomics and pathway analyses. J. Vet. Med. Sci. 82, 115–124. doi: 10.1292/jvms.19-0378,

PubMed Abstract | Crossref Full Text | Google Scholar

Gupta, S., Sharma, A., Joy, A., Dunshea, F. R., and Chauhan, S. S. (2022). The impact of heat stress on immune status of dairy cattle and strategies to ameliorate the negative effects. Animals 13:107. doi: 10.3390/ani13010107,

PubMed Abstract | Crossref Full Text | Google Scholar

Habibi, E., Shokrzadeh, M., Chabra, A., Naghshvar, F., Keshavarz-Maleki, R., and Ahmadi, A. (2015). Protective effects of Origanum vulgare ethanol extract against cyclophosphamide-induced liver toxicity in mice. Pharm. Biol. 53, 10–15. doi: 10.3109/13880209.2014.908399,

PubMed Abstract | Crossref Full Text | Google Scholar

Hao, Y., Gong, Y., Huang, S., Ji, S., Wang, W., Wang, Y., et al. (2021). Effects of age, diet CP, NDF, EE, and starch on the rumen Bacteria community and function in dairy cattle. Microorganisms 9:1788. doi: 10.3390/microorganisms9081788,

PubMed Abstract | Crossref Full Text | Google Scholar

He, Y., Maltecca, C., and Tiezzi, F. (2021). Potential use of gut microbiota composition as a biomarker of heat stress in monogastric species: a review. Animals 11:1833. doi: 10.3390/ani11061833,

PubMed Abstract | Crossref Full Text | Google Scholar

Hristov, A. N., Lee, C., Cassidy, T., Heyler, K., Tekippe, J. A., Varga, G. A., et al. (2013). Effect of Origanum vulgare L. leaves on rumen fermentation, production, and milk fatty acid composition in lactating dairy cows. J. Dairy Sci. 96, 1189–1202. doi: 10.3168/jds.2012-5975,

PubMed Abstract | Crossref Full Text | Google Scholar

Khan, R. I. N., Sahu, A. R., Malla, W. A., Praharaj, M. R., Hosamani, N., Kumar, S., et al. (2021). Systems biology under heat stress in Indian cattle. Gene 805:145908. doi: 10.1016/j.gene.2021.145908,

PubMed Abstract | Crossref Full Text | Google Scholar

Kim, G. H., Kim, J. E., Rhie, S. J., and Yoon, S. (2015). The role of oxidative stress in neurodegenerative diseases. Exp. Neurobiol. 24, 325–340. doi: 10.5607/en.2015.24.4.325,

PubMed Abstract | Crossref Full Text | Google Scholar

Li, L., Qu, J., Zhu, H., Liu, Y., Wu, J., Shao, G., et al. (2024). Effects of feeding different levels of dietary corn silage on growth performance, rumen fermentation and bacterial community of post-weaning dairy calves. Anim Biosci 37, 261–273. doi: 10.5713/ab.23.0174,

PubMed Abstract | Crossref Full Text | Google Scholar

Lin, L., Ma, H., Zhang, J., Yang, H., Zhang, J., Lai, Z., et al. (2024). Lignocellulolytic microbiomes orchestrating degradation cascades in the rumen of dairy cattle and their diet-influenced key degradation phases. Animal Adv. 1:0–0. doi: 10.48130/animadv-0024-0002

Crossref Full Text | Google Scholar

Liu, C., Wu, H., Liu, S., Chai, S., Meng, Q., and Zhou, Z. (2019). Dynamic alterations in yak rumen Bacteria community and metabolome characteristics in response to feed type. Front. Microbiol. 10:1116. doi: 10.3389/fmicb.2019.01116,

PubMed Abstract | Crossref Full Text | Google Scholar

Luo, Z., Liu, T., Cairang, D., Cheng, S., Hu, J., Shi, B., et al. (2024). Oregano essential oil as a natural plant additive affects growth performance and serum antibody levels by regulating the rumen microbiota of calves. Animals 14:820. doi: 10.3390/ani14060820,

PubMed Abstract | Crossref Full Text | Google Scholar

Ma, B., Zhang, L., Li, J., Xing, T., Jiang, Y., and Gao, F. (2021). Heat stress alters muscle protein and amino acid metabolism and accelerates liver gluconeogenesis for energy supply in broilers. Poult. Sci. 100, 215–223. doi: 10.1016/j.psj.2020.09.090,

PubMed Abstract | Crossref Full Text | Google Scholar

Mader, T. L., Davis, M. S., and Brown-Brandl, T. (2006). Environmental factors influencing heat stress in feedlot cattle1,2. J. Anim. Sci. 84, 712–719. doi: 10.2527/2006.843712x,

PubMed Abstract | Crossref Full Text | Google Scholar

Malik, P. K., Trivedi, S., Kolte, A. P., Mohapatra, A., Biswas, S., Bhattar, A. V. K., et al. (2023). Comparative analysis of rumen metagenome, metatranscriptome, fermentation and methane yield in cattle and buffaloes fed on the same diet. Front. Microbiol. 14:1266025. doi: 10.3389/fmicb.2023.1266025,

PubMed Abstract | Crossref Full Text | Google Scholar

Meessen, E. C. E., Majait, S., Ay, Ü., Olde Damink, S. W., Romijn, J. A., Holst, J. J., et al. (2025). Glycodeoxycholic acid inhibits primary bile acid synthesis with minor effects on glucose and lipid homeostasis in humans. J. Clin. Endocrinol. Metabol. 110, 1468–1477. doi: 10.1210/clinem/dgae399,

PubMed Abstract | Crossref Full Text | Google Scholar

Mikamo, H., Kawazoe, K., Sato, Y., Imai, A., and Tamaya, T. (1998). Preterm labor and bacterial intra-amniotic infection: arachidonic acid liberation by Phospholiphase A2 of prevotella bivia. Anaerobe 4, 209–212. doi: 10.1006/anae.1998.0165,

PubMed Abstract | Crossref Full Text | Google Scholar

Mohamed, M. (2025). Disruption of the spinal cord-gut axis alters gut microbial dynamics and carbohydrate cross-feeding. bioRxiv. doi: 10.1101/2025.05.27.656368

Crossref Full Text | Google Scholar

Morrison, S. F., Madden, C. J., and Tupone, D. (2014). Central neural regulation of Brown adipose tissue thermogenesis and energy expenditure. Cell Metab. 19, 741–756. doi: 10.1016/j.cmet.2014.02.007,

PubMed Abstract | Crossref Full Text | Google Scholar

Nakamura, K., and Morrison, S. F. (2010). A thermosensory pathway mediating heat-defense responses. Proc. Natl. Acad. Sci. USA 107, 8848–8853. doi: 10.1073/pnas.0913358107,

PubMed Abstract | Crossref Full Text | Google Scholar

Nazar, F. N., Videla, E. A., and Marin, R. H. (2019). Thymol supplementation effects on adrenocortical, immune and biochemical variables recovery in Japanese quail after exposure to chronic heat stress. Animal 13, 318–325. doi: 10.1017/S175173111800157X,

PubMed Abstract | Crossref Full Text | Google Scholar

Nehme, R., Andrés, S., Pereira, R. B., Ben Jemaa, M., Bouhallab, S., Ceciliani, F., et al. (2021). Essential oils in livestock: from health to food quality. Antioxidants 10:330. doi: 10.3390/antiox10020330,

PubMed Abstract | Crossref Full Text | Google Scholar

Newgard, C. B., An, J., Bain, J. R., Muehlbauer, M. J., Stevens, R. D., Lien, L. F., et al. (2009). A branched-Chain amino acid-related metabolic signature that differentiates obese and lean humans and contributes to insulin resistance. Cell Metab. 9, 311–326. doi: 10.1016/j.cmet.2009.02.002,

PubMed Abstract | Crossref Full Text | Google Scholar

Nguyen, M., Le Mignon, M., Schnellbächer, A., Wehsling, M., Braun, J., Baumgaertner, J., et al. (2023). Mechanistic insights into the biological activity of S-Sulfocysteine in CHO cells using a multi-omics approach. Front. Bioeng. Biotechnol. 11:1230422. doi: 10.3389/fbioe.2023.1230422,

PubMed Abstract | Crossref Full Text | Google Scholar

O’Donnell, V. B. (2022). New appreciation for an old pathway: the lands cycle moves into new arenas in health and disease. Biochem. Soc. Trans. 50, 1–11. doi: 10.1042/BST20210579,

PubMed Abstract | Crossref Full Text | Google Scholar

Orellana, L. H., Francis, T. B., Ferraro, M., Hehemann, J.-H., Fuchs, B. M., and Amann, R. I. (2022). Verrucomicrobiota are specialist consumers of sulfated methyl pentoses during diatom blooms. ISME J. 16, 630–641. doi: 10.1038/s41396-021-01105-7,

PubMed Abstract | Crossref Full Text | Google Scholar

Palevich, N., Kelly, W. J., Leahy, S. C., Denman, S., Altermann, E., Rakonjac, J., et al. (2019). Comparative genomics of rumen Butyrivibrio spp. uncovers a continuum of polysaccharide-degrading capabilities. Appl. Environ. Microbiol. 86:e01993-19. doi: 10.1128/AEM.01993-19,

PubMed Abstract | Crossref Full Text | Google Scholar

Park, T., Ma, L., Gao, S., Bu, D., and Yu, Z. (2022). Heat stress impacts the multi-domain ruminal microbiota and some of the functional features independent of its effect on feed intake in lactating dairy cows. J Animal Sci Biotechnol 13:71. doi: 10.1186/s40104-022-00717-z,

PubMed Abstract | Crossref Full Text | Google Scholar

Placha, I., Bacova, K., and Plachy, L. (2022). Current knowledge on the bioavailability of thymol as a feed additive in humans and animals with a focus on rabbit metabolic processes. Animals 12:1131. doi: 10.3390/ani12091131,

PubMed Abstract | Crossref Full Text | Google Scholar

Ridlon, J. M., Kang, D.-J., and Hylemon, P. B. (2006). Bile salt biotransformations by human intestinal bacteria. J. Lipid Res. 47, 241–259. doi: 10.1194/jlr.R500013-JLR200,

PubMed Abstract | Crossref Full Text | Google Scholar

Sharma, S., Zhang, X., Azhar, G., Patyal, P., Verma, A., Kc, G., et al. (2024). Valine improves mitochondrial function and protects against oxidative stress. Biosci. Biotechnol. Biochem. 88, 168–176. doi: 10.1093/bbb/zbad169,

PubMed Abstract | Crossref Full Text | Google Scholar

Sidney, T., Treon, E., Taiwo, G., Felton, E., Fan, P., and Ogunade, I. M. (2024). Exploring the rumen microbial function in Angus bulls with divergent residual feed intake. Front. Anim. Sci. 5:1485447. doi: 10.3389/fanim.2024.1485447

Crossref Full Text | Google Scholar

Skibiel, A. L., Zachut, M., Do Amaral, B. C., Levin, Y., and Dahl, G. E. (2018). Liver proteomic analysis of postpartum Holstein cows exposed to heat stress or cooling conditions during the dry period. J. Dairy Sci. 101, 705–716. doi: 10.3168/jds.2017-13258,

PubMed Abstract | Crossref Full Text | Google Scholar

Sommese, L., Pagliuca, C., Avallone, B., Ippolito, R., Casamassimi, A., Costa, V., et al. (2012). Evidence of Bacteroides fragilis protection from Bartonella henselae-induced damage. PLoS One 7:e49653. doi: 10.1371/journal.pone.0049653

Crossref Full Text | Google Scholar

Song, X., Luo, J., Fu, D., Zhao, X., Bunlue, K., Xu, Z., et al. (2014). Traditional Chinese medicine prescriptions enhance growth performance of heat stressed beef cattle by relieving heat stress responses and increasing apparent nutrient digestibility. Asian Australas. J. Anim. Sci. 27, 1513–1520. doi: 10.5713/ajas.2014.14058,

PubMed Abstract | Crossref Full Text | Google Scholar

Sun, H., Olson, K. C., Gao, C., Prosdocimo, D. A., Zhou, M., Wang, Z., et al. (2016). Catabolic defect of branched-chain amino acids promotes heart failure. Circulation 133, 2038–2049. doi: 10.1161/CIRCULATIONAHA.115.020226,

PubMed Abstract | Crossref Full Text | Google Scholar

Terenina, M. B., Misharina, T. A., Krikunova, N. I., Alinkina, E. S., Fatkulina, L. D., and Vorob’yova, A. K. (2011). Oregano essential oil as an inhibitor of higher fatty acid oxidation. Appl. Biochem. Microbiol. 47, 445–449. doi: 10.1134/S0003683811040181,

PubMed Abstract | Crossref Full Text | Google Scholar

Ulrich, K., and Jakob, U. (2019). The role of thiols in antioxidant systems. Free Radic. Biol. Med. 140, 14–27. doi: 10.1016/j.freeradbiomed.2019.05.035,

PubMed Abstract | Crossref Full Text | Google Scholar

Van Passel, M. W. J., Kant, R., Zoetendal, E. G., Plugge, C. M., Derrien, M., Malfatti, S. A., et al. (2011). The genome of Akkermansia muciniphila, a dedicated intestinal mucin degrader, and its use in exploring intestinal metagenomes. PLoS One 6:e16876. doi: 10.1371/journal.pone.0016876,

PubMed Abstract | Crossref Full Text | Google Scholar

Walasek-Janusz, M., Grzegorczyk, A., Malm, A., Nurzyńska-Wierdak, R., and Zalewski, D. (2024). Chemical composition, and antioxidant and antimicrobial activity of oregano essential oil. Molecules 29:435. doi: 10.3390/molecules29020435,

PubMed Abstract | Crossref Full Text | Google Scholar

Wang, J., Li, J., Wang, F., Xiao, J., Wang, Y., Yang, H., et al. (2020). Heat stress on calves and heifers: a review. J. Animal Sci. Biotechnol 11:79. doi: 10.1186/s40104-020-00485-8,

PubMed Abstract | Crossref Full Text | Google Scholar

Wei, X., Zhu, Y., Zhao, X., Zhao, Y., Jing, Y., Liu, G., et al. (2022). Transcriptome profiling of mRNAs in muscle tissue of Pinan cattle and Nanyang cattle. Gene 825:146435. doi: 10.1016/j.gene.2022.146435,

PubMed Abstract | Crossref Full Text | Google Scholar

West, J. W. (2003). Effects of heat-stress on production in dairy cattle. J. Dairy Sci. 86, 2131–2144. doi: 10.3168/jds.S0022-0302(03)73803-X,

PubMed Abstract | Crossref Full Text | Google Scholar

Wirth, R., Kádár, G., Kakuk, B., Maróti, G., Bagi, Z., Szilágyi, Á., et al. (2018). The planktonic Core microbiome and Core functions in the cattle rumen by next generation sequencing. Front. Microbiol. 9:2285. doi: 10.3389/fmicb.2018.02285,

PubMed Abstract | Crossref Full Text | Google Scholar

Wright, D. P., Knight, C. G., Parkar, S. G., Christie, D. L., and Roberton, A. M. (2000). Cloning of a mucin-Desulfating sulfatase Gene from Prevotella strain RS2 and its expression using a Bacteroides recombinant system. J. Bacteriol. 182, 3002–3007. doi: 10.1128/JB.182.11.3002-3007.2000,

PubMed Abstract | Crossref Full Text | Google Scholar

Yan, Q., Li, X., Zhou, X., Chen, W., Tian, X., Wittayakun, S., et al. (2024). Macleaya cordata extract exhibits some potential as a surrogate antibiotic by improving gastrointestinal epithelial status and humoral response in goats. Animal Nutr. 18, 356–366. doi: 10.1016/j.aninu.2024.04.014,

PubMed Abstract | Crossref Full Text | Google Scholar

Yao, Q., Gao, Y., Fan, L., Wang, J., and Zheng, N. (2022). 2′-Fucosyllactose remits colitis-induced liver oxygen stress through the gut–liver–metabolites Axis. Nutrients 14:4186. doi: 10.3390/nu14194186,

PubMed Abstract | Crossref Full Text | Google Scholar

Yousefi, M., Hoseini, S. M., Vatnikov, Y. A., Karamyan, A., and Kulikov, E. V. (2024). Dietary thymol supplementation promotes antioxidant responses and thermal stress resistance in rainbow trout, Oncorhynchus mykiss. Animals 14:2988. doi: 10.3390/ani14202988,

PubMed Abstract | Crossref Full Text | Google Scholar

Yu, Y.-M., Chao, T.-Y., Chang, W.-C., Chang, M. J., and Lee, M.-F. (2016). Thymol reduces oxidative stress, aortic intimal thickening, and inflammation-related gene expression in hyperlipidemic rabbits. J. Food Drug Anal. 24, 556–563. doi: 10.1016/j.jfda.2016.02.004,

PubMed Abstract | Crossref Full Text | Google Scholar

Zeng, Z., Xu, X., Zhang, Q., Li, P., Zhao, P., Li, Q., et al. (2015). Effects of essential oil supplementation of a low-energy diet on performance, intestinal morphology and microflora, immune properties and antioxidant activities in weaned pigs. Anim. Sci. J. 86, 279–285. doi: 10.1111/asj.12277,

PubMed Abstract | Crossref Full Text | Google Scholar

Zhang, Y., Deng, Y., Hao, Y., Fang, J., and Feng, J. (2024). Effects of supplementation with oregano essential oil during late gestation and lactation on serum metabolites, antioxidant capacity and Fecal microbiota of sows. Animals 14:753. doi: 10.3390/ani14050753,

PubMed Abstract | Crossref Full Text | Google Scholar

Zhang, R., Wu, J., Lei, Y., Bai, Y., Jia, L., Li, Z., et al. (2021). Oregano essential oils promote rumen digestive ability by modulating epithelial development and microbiota composition in beef cattle. Front. Nutr. 8:722557. doi: 10.3389/fnut.2021.722557,

PubMed Abstract | Crossref Full Text | Google Scholar

Zhao, Q., Duan, Z., Lai, R., and Ma, P. (2025). Novel microbially transformed bile acids: biosynthesis, structure, and function. Pharmacol. Res. 216:107775. doi: 10.1016/j.phrs.2025.107775,

PubMed Abstract | Crossref Full Text | Google Scholar

Zhao, S., Min, L., Zheng, N., and Wang, J. (2019). Effect of heat stress on bacterial composition and metabolism in the rumen of lactating dairy cows. Animals 9:925. doi: 10.3390/ani9110925,

PubMed Abstract | Crossref Full Text | Google Scholar

Zhou, R., Wu, J., Lang, X., Liu, L., Casper, D. P., Wang, C., et al. (2020). Effects of oregano essential oil on in vitro ruminal fermentation, methane production, and ruminal microbial community. J. Dairy Sci. 103, 2303–2314. doi: 10.3168/jds.2019-16611,

PubMed Abstract | Crossref Full Text | Google Scholar

Zhou, R., Wu, J., Zhang, L., Liu, L., Casper, D. P., Jiao, T., et al. (2019). Effects of oregano essential oil on the ruminal pH and microbial population of sheep. PLoS One 14:e0217054. doi: 10.1371/journal.pone.0217054,

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: heat stress, metabonomics, microbiota, oregano essential oil, Pinan cattle (Piedmontese × Nanyang)

Citation: Wang J, Liu L, Sheng J, Zhang L, Li Q, Gao T and Liu X (2026) Effects of different doses of oregano essential oil on growth performance, health parameters, and the gut microbiome-metabolome profile in heat-stressed Pinan cattle (Piedmontese × Nanyang). Front. Microbiol. 16:1712904. doi: 10.3389/fmicb.2025.1712904

Received: 25 September 2025; Revised: 28 November 2025; Accepted: 09 December 2025;
Published: 29 January 2026.

Edited by:

Xiangfeng Kong, Institute of Subtropical Agriculture (CAS), China

Reviewed by:

Yang Liu, Hunan Institute of Animal Husbandry and Veterinary Medicine, China
Xinjun Qiu, Chinese Agricultural University, China

Copyright © 2026 Wang, Liu, Sheng, Zhang, Li, Gao and Liu. 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: Tengyun Gao, Z2FvdGVuZ3l1bkBoZW5hdS5lZHUuY24=; Xian Liu, bGl1eGlhbjY0MUAxNjMuY29t

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.