Rumen Microbial Metabolic Responses of Dairy Cows to the Honeycomb Flavonoids Supplement Under Heat-Stress Conditions

Introduction

Heat stress traditionally refers to the combination of the excessive heat accumulation and the reductive heat dissipation, which may seriously induce fervescence (1, 2), the reduction of feed intake (3), the turbulence of rumen microbiota, and the reduction of milk production of dairy cows (4). During the last few decades, researchers dedicated to find the proper methods in ameliorating heat stress during dairy cow production and proposed certain solutions including spraying, fans, shading, and nutritional methods (3).

The most commonly applied nutritional methods to ameliorate heat stress is to enhance rumen digestibility of dairy cows through the addition of minerals, vitamins, prebiotics, and probiotics with balanced nutrition (5). Compared with the initial supplement of nutritional elements, microbial secondary metabolites, such as flavonoids, played an essential role in shaping microbial communities, maintaining homoeostasis, and promoting host growth (6) by their splendid antioxidation and free radical scavenging capabilities (7–9). Moreover, flavonoids also regulated nutrient absorptivity by enhancing ruminal epithelial development (10). Li et al. (11) proved the ameliorative effects of flavonoids on dairy cows under heat stress conditions. In contrast, how flavonoids worked on ruminal microbial communities and metabolism is still not clear.

Rumen microbiota is the stabilized and powerful self-regulated digestive ecosystem that converts feed nutrients into absorbable compounds for the host (12). Ruminal microbiota communities also showed the accommodated responses to heat stress, such as the decrease of higher fiber degradative bacteria, Fibrobacter, and the increase of starch-utilizing bacteria, Clostridium and Streptococcus, which further reduced heat generation and utilized more easily fermented carbohydrates to maintain the energy provision (13). Microbial metabolites are subsequently altered to adapt to the change of microbial communities, and further regulated the growth and production of the host (14). Disorders of ruminal microbial ecosystem and the following declines of nutritional digestibility when subjected to heat stress mainly caused the reduction of dairy production (15, 16). Therefore, the most critical point in attenuating heat stress is the re-establishment of the physiological thermal energy balance between heat acquisition and heat dissipation through the regulation of ruminal microbial community and nutrient metabolism.

In the present article, flavonoids extracted from honeycomb were supplemented to investigate the productive, rumen microbial, and rumen metabolism responses of cows under heat stress conditions, which aimed to provide a selectable suggestion for dairy production. We hypothesized that flavonoids help promote the proliferation of fiber-degrading bacteria, following the regulation of energy metabolism, and ultimately enhanced the production performances of dairy cows.

Materials and Methods

Animal care and procedures followed the Chinese Guidelines for Animal Welfare, which was approved by the Animal Care and Use Committee of Jiangxi Agricultural University, with the approval number JXAULL-20210619.

Experimental Design

The experiment was conducted in the Bengbu dairy farm, Modern Farming (Wuhe) Co. Ltd., Anhui Province, China (32.92N, 117.38E) from July 7, 2021 to August 7, 2021. A total of 600 multiparous dairy herds with an average live weight of 643.6 ± 37.3 kg, average lactation of 159.3 ± 18.8 days, and average calving parities of 2.58 ± 0.49 were randomly allocated into the control treatment (CON), the heat stress treatment (HS), and honeycomb flavonoids supplement under heat stress conditions treatment (HF). Each treatment was distributed to an individual stall, and with the use of a transverse ventilation system, the wind intensity could be regulated by automatic ventilators. Each stall contains four pens and 50 cows in each pen. Each pen was considered as a replicate and arranged along the wind direction to ensure that each pen receives similar wind intensity. During the trial period, cows in control treatment received the spraying treatment to ensure non-occurrence of heat stress. No cooling methods were provided for the HS and HF treatments.

All cows were provided the same diets, which were formulated according to NRC (2001) to meet or exceed the energy requirements of Holstein dairy cows yielding 30 kg of milk/day with 3.5% milk fat and 3.0% true protein. Diets were fed three times per day at 06:00, 13:00, and 21:00. Details of ingredient analysis and chemical composition of diets are shown in Table 1.

TABLE 1

Table 1. Ingredients and chemical composition of the TMR (dry-matter basis).

Temperature and humidity of each feeding lair were recorded every day, and the temperature humidity indexes (THI) were calculated through the following equation:

$\begin{array}{l}\text{THI}=\left(1.8\times \text{T}+32\right)-\left[\left(0.55-0.0055\times \text{RH}\right)\right]\times \\ \left(1.8\times \text{T}-26\right)\end{array}$

where T = temperature (°C), RH = relative humidity (%).

THI ≥70 was considered as heat stress occurring in high-yielding dairy cows (17). Body temperature was measured through rectal thermometry once a week through a 30-size sampling inspection of each treatment.

Raw honeycomb was provided by the Institute of Bee Research, Jiangxi Agricultural University, Jiangxi, China. HF was further extracted in the laboratory using an ultrasonic (Elmasonic X-tra Flexl; Elma, Konstanz, Germany) extracting procedure combined with the solid-to-liquid ratio of 1:29.76, 56.33% of ethanol concentration, and the reaction time of 28.88 min. Extractions were further centrifugated at 3,000 rpm for 10 min. The postextracting liquor was filtered after cooling and collected into an Erlenmeyer flask.

Flavonoids content measurement were subsequently applied based on the Chinese National standard GB/T 20574-2006. To be simply stated, flavonoids standard samples were first precisely weighed at 10 mg and dissolved into 50 ml ethanol. Gradients of 1, 2, 3, 4, 5, and 6 ml of the standard solution were collected into a 50-ml volumetric flask, and isochorically added to 15 ml of 95% ethanol. Then, 1 ml of 100 g/L aluminum nitrate solution and 1 ml of 9.8 g/L potassium acetate solution were added into the mixture, respectively. Absorption of each standard solution was measured at 415 nm wavelength after standing for 60 min. The standard curve was then made as follows:

$\begin{array}{l}\text{y}=0.4145\text{x}+0.001({\text{R}}^2=0.9997)\end{array}$

where y = flavonoids content and x = absorption.

Honeycomb flavonoids content was measured about 45% of the standard flavonoids and supplemented with 10 g/day for each cow.

Production Performance Measurement

Average daily intake was shown as the dry matter intake (DMI) of each treatment, which were calculated through the deviation between the supply and residue. DMI were calculated and displayed as the average of four replicates. Cows were milked three times per day (08:00, 14:00, and 20:00, respectively) and the daily milk production was automatically recorded through the rotary milking facilities (9JRP-50P2100; Delaval, Israel). Milk samples were collected in 100-ml vials from each replicate at the last day and subsequently stored at 4°C with the addition of 2-bromo-2-nitropropan-1,3-diol for further milk quality analysis including milk protein, milk fat, and somatic cell count.

Ruminal Fermentation Parameter Measurement

Rumen fluid samples of 24 cows (two cows in each replicate, eight cows per treatment, second lactation, similar body weight) were collected at the last day through esophageal tubing at 3 h after morning feeding. All samples were divided into two portions. One was conducted to analyze ruminal pH, rumen volatile fatty acids (VFAs), and ammonia-N (NH₃-N). The other portion was frozen in the liquid nitrogen immediately and then stored at −80°C for further microbiota and metabolite measurement. Portable type pH meter (Testo 205; Testo AG, Lenzkirch, Germany) was applied for the measurement of ruminal pH immediately after rumen fluid sample was collected. Individual and total VFAs in the aliquots were measured using a gas chromatograph (GC-2010; Shimadzu, Kyoto, Japan). Concentration of NH₃-N was determined by indophenol method and the absorbance value was measured through UV-2600 ultraviolet spectrophotometer (Tianmei Ltd., China) at the 700 nm wavelength (12).

Digestibility Measurement

Digestibility measurement was conducted during the last 3 days. To be simply stated, total fecal samples were collected at intervals of 6 h (18 samples in total) after 12-h-long fasting treatment. Digestibility of dry matter (DM), crude protein (CP), neutral detergent fiber (NDF), and fatty acids (EE) were subsequently measured based on the methods we stated previously (18).

Rumen Microbial Community Measurement

Rumen fluid DNA was extracted using Bacterial Genome DNA Extraction Kit (DP302, TIANGEN; TIANGEN BIOTECH (BEIJING) Co., Ltd). The V4 region of the 16S rRNA gene was amplified using the universal primers 520F and 802R (F: GTGCCAGCMGCCGCGGTAA and R: GGACTACHVGGGTWTCTAAT) (18). All PCR reactions were carried out with Phusion High-Fidelity PCR Master Mix (New England Biolabs). The mixture of PCR products was purified with Qiagen Gel Extraction Kit (Qiagen, Hilden, Germany), followed by the generation of sequencing libraries using TruSeq DNA PCR-Free Sample Preparation Kit (Illumina, USA). The library quality was assessed on the Qubit@ 2.0 Fluorometer (Thermo Scientific) and Agilent Bioanalyzer 2100 system. Finally, the library was sequenced on Illumina HiSeq 4000 platform (Illumina Inc., San Diego, USA).

Quality filtering of raw tags were performed under specific filtering conditions to obtain the high-quality clean tags according to the Quantitative Insights Into Microbial Ecology (QIIME, V1.7.0) quality controlling process. Sequences with >97% similarity were assigned into the same operational taxonomic unit (OTU).

Rumen Metabolite Measurement

Rumen metabolomics in the present study were measured using GC/MS analysis method through Agilent 7890 gas chromatograph system coupled with a Pegasus HT time-of-flight mass spectrometer (LECO, St. Joseph, MI). The main process was detailed as stated in our previous study (19). To simply state, in all samples were added 55 μl of methoxy amination reagent (20 mg/ml dissolved in pyridine) followed by a 20-min-long incubation at 80°C. Subsequently, all samples were incubated for 1 h at 70–80°C after adding 75 μl BSTFA reagent (1% TMCS, v/v). Thereafter, 10 μl fatty acid methyl ester (FAMEs) (standard mixture of fatty acid methyl esters, C8–C16:1 mg/ml; C18–C24:0.5 mg/ml in chloroform) was added into each sample after all samples were cooled to room temperature. Finally, all samples were prepared for GC–MS analysis after mixing them well.

Statistical Analysis

Body temperature; production performances including DMI and milk yield; milk quality parameters including milk fat, milk protein, and SCC; and ruminal fermentation variables including ruminal pH, VFAs, and ammonia-N were first determined by conducting a normal distribution test using SAS procedure “proc univariate data=test normal” and subsequently received one-way ANOVA S-N-K test by SAS (SAS Institute, Inc., Cary, NC, USA). Significance would be considered when p-value <0.05 while a tendency was considered when 0.05 ≤ p < 0.10. OTU abundances of each rumen bacteria were first determined by conducting a percentage transformation, and then one-way ANOVA S-N-K test of SAS 9.2 was applied for the differential analysis. Alpha diversity and beta diversity in our samples were calculated with QIIME 2 (20) and displayed with R software (version 3.3.1; R Core Team, Vienna, Austria). Principal coordinate analysis (PCoA) for different rumen methanogens were conducted using R “vegan package”.

For metabolomics data analysis, Chroma TOF 4.3X software of LECO Corporation and LECO-Fiehn Rtx5 database were used for peak identification and integration of the peak area. Multivariate analysis including principal component analysis (PCA) were conducted using SIMCA-P software (V.14.0; Umetrics, Umea, Sweden). Differentially expressed metabolites between two treatments were identified based on variable importance in projection (VIP) and statistical analysis (VIP > 1 and p < 0.05). Functional enrichment analysis was conducted using Metaboanalyst 5.0 (https://www.metaboanalyst.ca/). Significant function enrichment would be considered when p-value <0.05.

Results

THI Recording

In our study, temperature and relative humidity of each treatment were recorded in five different locations, and the THI was calculated based on the average temperature and humidity. THI was recorded throughout the whole experimental period and displayed in Figure 1. During the entire experimental period, environmental THI exceeded 70 for HS treatments, while during 27 days of CON treatment, it was below 70. HF treatment showed the same THI with HS for the cows that received HF treatment and reared in the same stall with HS. THI provided a confirmation of the occurrence of heat stress.

FIGURE 1

Figure 1. THI record of both heat stress and control treatments. HS, heat stress; CON, control treatment.

Production Performances and Milk Quality

The production parameters for dairy cows were first recorded, and the results are shown in Table 2. Heat stress significantly restrained DMI and milk yield compared with CON treatment, while both parameters increased after honeycomb flavonoids supplement, however, not significantly. Besides, honeycomb flavonoids significantly reduced the milk somatic cell counts, which might be the indicator of mastitis, and remarkably increased under heat stress conditions. Body temperature in HS treatment significantly increased compared with CON, while it decreased after honeycomb flavonoids supplement. There is an increasing trend of milk fat after HF treatment, while no significant changes were found for milk protein content.

TABLE 2

Table 2. Effects of honeycomb flavonoids supplement on production performances and milk quality during heat stress conditions.

Rumen Fermentability and Nutrient Digestibility

Rumen fermentability, including ruminal pH, VFAs, and ammonia-N, were measured, and the results are shown in Table 3. Heat stress reduced ruminal acetate, propionate, and butyrate content while the acetate and ratio of acetate/propionate significantly decreased (p < 0.05). Acetate content increased significantly after HF supplement (p < 0.05), although still inferior to CON treatment. In addition, the ratio of acetate/propionate in HF were significantly higher than that in HS. No significant differences were detected for other parameters among CON, HS, and HF treatments.

TABLE 3

Table 3. Effects of honeycomb flavonoids supplement on rumen fermentable parameters during heat stress conditions.

Furthermore, nutrient digestibility was measured through total feces collection method, and the results are shown in Table 4. Digestion of dry matter, NDF, and fatty acids was enhanced after HF supplement while the digestibility of NDF increased significantly compared with HS. No significant changes were found for the CP digestion.

TABLE 4

Table 4. Effects of honeycomb flavonoids supplement on cow digestibility under heat stress conditions.

Rumen and Microbial Communities

A total of 4,700 OTUs, 15 phyla, and 2,600 genera were identified after quality control, and all the taxonomic information is displayed in Supplementary Table 1. All identified bacteria were chosen for further investigation of α-diversity, β-diversity, and differential communities.

α-Diversity

Alpha diversity was first investigated to determine the complexity of rumen microbial diversity through Chao1, Shannon, Simpson, and ACE indexes; all results are shown in Table 5. Alpha-diversity indexes including ACE and Chao1 significantly decreased under heat stress conditions (p < 0.05). Meanwhile, the aforementioned indexes increased after HF supplement, however, not significantly. No changes were detected for other indexes among all the three treatments.

TABLE 5

Table 5. Effects of honeycomb flavonoids supplement on α-diversity of ruminal microbiota during heat stress conditions.

β-Diversity

Subsequently, differential analyses on rumen microbial communities among the three treatments were applied. Principal coordinates analysis (PCoA) analysis was first proceeded. As shown in Figure 2, PCoA axes 1 and 2 accounted for 52.21 and 16.83%, respectively. Bacterial communities in HS treatment could be significantly separated from those in CON through PCoA axes 1 and 2, except HS6. Meanwhile, bacterial communities were significantly altered after HF supplement compared with heat stress conditions.

FIGURE 2

Figure 2. Principal coordinate analysis (PCoA) on ruminal bacterial community structures of CON, HS, and HF treatments. CON, control treatment; HS, heat stress treatment; HF, honeycomb flavonoids supplement treatment.

Thereafter, differential analyses on relative abundances of ruminal bacterial communities in the level of phyla and genera were determined. As the results have shown in Table 6, Firmicutes and Bacteroidetes contributed the top 2 bacterial biomass. Heat stress treatment significantly proliferated the relative abundance of Firmicutes, while it suppressed that of Bacteroidetes. Honeycomb flavonoids slightly increased the relative abundances of Bacteroidetes while it reduced Firmicutes, however, not significantly. It was especially noticed that HF supplement significantly increases the abundance of Proteobacteria. No changes were found for other phyla among all the three treatments.

TABLE 6

Table 6. Effects of honeycomb flavonoids supplement on relative abundances of ruminal microbiota communities during heat stress conditions (level of phyla, %).

As referred to the genera level, shown in Table 7, Prevotella, Ruminococcaceae, Succiniclasticum, Lachnospiraceae, and Eubacterium located the five most abundant genera and accounted nearly 90% of all microbiota profiles. Microbial communities including Succiniclasticum, Eubacterium, Acetitomaculum, Bifidobacterium, and Succinivibrio significantly proliferated while the growth of Ruminococcaceae, Ruminococcus, Pseudobutyrivibrio, and Streptococcus significantly suppressed under heat stress conditions. However, HF supplement effectively attenuated these changes mainly manifested in the increasing abundance of Succiniclasticum, Acetitomaculum, Bifidobacterium, and Succinivibrio, and the decrease of Ruminococcaceae and Ruminococcus. Noticeably, the main starch-degrading bacteria such as Pseudobutyrivibrio and Streptococcus increased after HF supplement compared with HS treatment.

TABLE 7

Table 7. Effects of honeycomb flavonoids supplement on relative abundances of ruminal microbiota communities during heat stress conditions (level of genera, %).

Effects of Honeycomb Flavonoids Supplement on Rumen Metabolite Alterations Under Heat Stress Conditions

A total of 1,169 peaks were detected with GC–TOF–MS, and 1,150 of metabolites were further identified across all samples by retaining the treatments with null value ≤ 50%. All identified metabolites are provided in Supplementary Table 2. The metabolites mainly clustered into organic acids, fatty acids, carbohydrates, amino acids, purines, and biogenic amines.

Differential analysis on the integrity of metabolite alteration of honeycomb flavonoids supplement under heat stress treatment was first conducted through PCA. As shown in Figure 3A, PC1 and PC2 accounted for 55.9 and 20.8% of the total variation, respectively. Integrally speaking, samples of HS treatment could clearly separate from that of CON, except HS 6, which indicated that rumen metabolites showed a significant alteration under heat stress conditions compared with CON. In addition, metabolites also showed a significant difference after HF supplement compared with those under heat stress treatment as shown in Figure 3B. PC1 and PC2 accounted for 70.7 and 12.6% of the total variation, respectively. Samples of HF treatment could clearly separate from HS based on PC1 and PC2, except HF6.

FIGURE 3

Figure 3. Principal components analysis (PCA) on ruminal metabolite content between different treatments. (A) Principal components analysis (PCA) on ruminal metabolites content between CON and HS treatments. CO, control treatment; HS, heat stress treatment. (B) Principal components analysis (PCA) on ruminal metabolite structures between HS and HF treatments. HS, heat stress treatment; HF, honeycomb flavonoids supplement treatment.

Thereafter, significantly altered metabolites under heat-stressing treatment compared with control, and those after HF supplement compared with HS, were analyzed based on the statistical standard of VIP >1 and p-value <0.05. All the results are shown in Tables 8, 9, respectively. A total of 35 significant differential altered metabolites including 16 downregulated and 19 upregulated metabolites were detected based on the filtering standard. The downregulated metabolites mainly related with the carbohydrate metabolism such as FMNH and cis-acetylacrylate; amino acid metabolism including isoleucyl-asparagine, dimethyl-L-arginine, and serinyl-threonine; and certain organic acids. Besides, most upregulated metabolites related with nucleotide metabolism, such as hypoxanthine, pyrimidin-2-one-ribonucleoside, 2-phenylacetamide, xanthine, thymidine, and a number of bioamines.

TABLE 8

Table 8. Ruminal metabolite alterations under heat stress conditions.

TABLE 9

Table 9. Ruminal metabolite alterations to the honeycomb flavonoids supplement under heat stress conditions.

The investigation on rumen metabolites changes after HF supplements were subsequently provided. As shown in Table 8, a total of 20 significant differential altered metabolites, including 11 upregulated metabolites and nine downregulated metabolites, were identified. The mainly upregulated metabolites are clustered into lipid metabolism, such as phosphatidylglycerol (PG), phosphatidylinositol (PI), phosphatidylserine (PS), and phosphatidylethanolamine (PE). Besides, the downregulated metabolites were mostly bioamines, which formed from protein degradation, and might insult ruminal epithelium when excessively accumulated.

Finally, the functional enrichment analysis was applied based on all differential metabolites to determine the mainly altered functions under heat stress conditions and after honeycomb flavonoids supplement. All enrichment results are shown in Figure 4.

FIGURE 4

Figure 4. Functional enrichment analysis on differential ruminal metabolites. (A) Functional enrichment analysis on differential ruminal metabolites between CON and HS treatments. CON, control treatment; HS, heat stress treatment; p < 0.05 was considered as significant difference. (B) Functional enrichment analysis on differential ruminal metabolites after honeycomb flavonoids supplement under heat stress conditions. CON, control treatment; HS, heat stress treatment; p < 0.05 was considered as significant difference.

Figure 4A shows the critical altered functions enriched between heat stress and control treatments. Vitamin B₆ metabolism, phosphatidylinositol phosphate metabolism, and spermidine biosynthesis accounted for the three most enriched functions, which might indicate the suppression of both carbohydrate metabolism and regular protein metabolism under heat stress conditions. Furthermore, functional investigation of honeycomb flavonoids supplement on ruminal metabolism was applied and shown in Figure 4B. The energy metabolism–related functions such as beta-oxidation of fatty acids, and the bacterial proliferation–related functions such as purine metabolism presented a major difference after honeycomb flavonoids supplementation. Besides, amino acids and their derivate-related function also showed a great promotion after HF treatment.

Discussion

Heat stress is seriously detrimental to dairy production and economic benefits, especially in South China. When ambient temperature and humidity do not decrease enough, dairy cows are unable to lose heat and further cause them to be in a constant state of heat stress. Cows supplemented with Chinese medicines, which contain abundant flavonoids, showed better immune condition and alleviation effects under the heat stress environment (21). Therefore, in our present study, we purified flavonoids from honeycomb and supplemented it into dairy cows to determine both the nutrient digestibility and milk productivity responses under heat stress environments.

Effects of Honeycomb Flavonoids Supplement on Rumen Microbial Homeostasis and Nutrient Degradability

Energy, which is mainly derived from nutrient ingestion, restricted neither the growth nor the production performances of dairy cows. Absence of energy provision under heat stress conditions directly suppressed rumen microbial proliferation and degradability, and thereafter disrupted ruminal homeostasis.

Indeed, ruminal homeostasis provided an ideal environment for nutrient digestion, transportation, energy generation, and the microbial proliferation (22, 23). Hence, ruminal microbiota generally showed adaptive responses to heat stress. Relative abundances of fiber-degrading bacteria including Fibrobacteres and Bacteroidetes considerably decreased because of the high energy dissipation and exothermic effects during the fiber degradation process, while starch-utilizing bacteria such as Firmicutes and Streptococcus showed more degradable activities (24, 25). Besides, heat-stressed cows also had a significantly higher relative abundance of lactate-producing bacteria (e.g., Streptococcus), which utilize soluble carbohydrate as an energy source (26). All these adaptive changes were also observed in the HS treatment compared with CON, which affirmed the occurrence of heat stress. Microbial changes consequentially resulted in lower VFA production and energy provision. The relative abundance of the acetate-producing bacterium Acetobacter decreased with heat stress treatment (16). Such a finding was also observed in the present study and may give proof of the decrease of acetate content.

All the aforementioned drawbacks may be conversed after flavonoids supplement for the regulatory effects on microbial proliferation. Generally, flavonoids are presented in their glycoside forms, which were further deglycosylated by β-glucosidases secreted by microbiota and resulted in the increase of aglycons (27). The aglycons were further utilized by ruminal bacteria to enhance proliferation of bacteria. Interestingly, HF supplement in the present study significantly increases both cellulolytic and amylolytic bacteria. These alterations help recover microbial construction and further energy metabolism systems, which were reflected in the more abundant content of total VFAs and acetate. These changes might be attributed to the increase of α-diversity that utilized more carbohydrates into VFAs and the increase of rumen cellulolytic bacteria. For the ruminal acetate, it was primarily derived from cellulose degradation (28).

Besides, as shown by metabolomics result, beta oxidation of fatty acids, which could further generate acetyl-CoA for energy provision (29), was significantly enhanced after HF supplement. Fatty acids possessed superior energy to carbohydrates and might supply more energy for rumen microbial proliferation; thus, the re-establishment of micro-ecosystem could be recovered.

In addition, proliferated microbiota enhanced ruminal protein utilization. Rumen cellulolytic bacteria have been proven to use ammonia as their sole nitrogen source because the incorporation of preformed amino acids is minimal (30). The proliferated cellulolytic bacteria consumed more ammonia, which causatively induced the reduction of bioamines after HF supplement, as aforementioned in the metabolomics results. Moreover, the promoted ammonia utilization catalyzed protein degradation; therefore, nitrogen utilization was enhanced.

Effects of Honeycomb Flavonoids Supplement on Milk Production and Milk Quality

Milk production is superior to any other parameters for the dairy cattle industry because of both economic and effective factors. DMI and gastrointestinal absorptivity contributed the most on milk production. In our research, honeycomb flavonoids promoted milk production; reasons could be attributed to enhancing nutrient absorptivity.

Nutrient absorptivity is mainly affected by epithelial transporting capacity. The epithelial barrier and transporting functions were disrupted under heat-stressed conditions by severe inflammatory responses because of the lipopolysaccharide release caused by the death of bacteria (31). Flavonoids may benefit from the epithelial integrity and transporting efficiency because of their mighty anti-inflammatory capacities, which therefore enhanced nutrient absorptivity (32). In addition, bile acid (BA) metabolism might be another causative factor. In recent research, the BA signaling pathway was proved to be critical in both lipid and carbohydrate metabolism, also in maintaining esoteric glucose and cholesterol homeostasis as well as immune states (33). BA is also a crucial secondary metabolite that is utilized by gut microbiota and communicated with epithelial cells. Intriguingly, flavonoids could be responsible for increased BA production (6), and further help promote lipid and carbohydrate metabolism to provide more energy for milk production.

Milk fat content increased after HF supplement. Milk fat is a necessary parameter of milk quality and beneficial to human health, and is mostly formed from acetate generated by ruminal fermentation (28). Based on the previous finding, acetate is mainly generated by cellulose degradation (16), and the increasing cellulose-degrading bacteria, especially the mighty acetate-generating bacteria Acetitomaculum, significantly increased after HF supplement, which finally resulted in the increase in acetate content and further milk fat content.

Moreover, improved ruminal lipid metabolism, as showed in our metabolomics result, may be another causative factor for the increase of milk fat content. When dietary lipids enter the rumen, the initial step in lipid metabolism is the hydrolysis of the ester linkages found in TG, phospholipids, and glycolipids, and this is primarily carried out by hydrolases produced by rumen bacteria (34). The hydrolyzed lipid was further transported into mammary gland and synthesized milk fat as substrates. Therefore, milk fat content greatly increased after honeycomb flavonoids supplement.

Conclusion

In summary, production performances were seriously suppressed under heat stress conditions and caused huge economic loss. Honeycomb flavonoids supplement help proliferate microbial abundances and fiber-degrading bacteria, which further promoted nutrient degradability and energy metabolism, and ultimately enhanced the production performances of dairy cows. These findings in the present study might provide a proper selection for dairy production under heat stress conditions.

Data Availability Statement

The data presented in the study are deposited in the NCBI Sequence Read Archive (SRA, http://www.ncbi.nlm.nih.gov/Traces/sra/), accession number PRJNA753017.

Ethics Statement

The animal study was reviewed and approved by Animal care and procedures followed the Chinese Guidelines for Animal Welfare, which was approved by the Animal Care and Use Committee of Jiangxi Agricultural University, with the approval number JXAULL-20210619.

Author Contributions

EL and HW designed the study. EL, MS, CH, KM, and QL conducted the experiment. JZ, DW, SW, CZ, WL, and SG participated in parameter measurement and data analysis. FX wrote the article and contributed to English editing. All authors contributed to the article and approved the submitted version.

Funding

The study was funded by the Science and Technology Planning Project of Jiangxi Educational Department (GJJ200414).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher's Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Acknowledgments

We thank the Nanchang Key Laboratory of Animal Health and Safety Production and Jiangxi Province Key Laboratory of Animal Nutrition/Engineering Research Center of Feed Development.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fvets.2022.845911/full#supplementary-material

Supplementary Table 1. Taxonomy results of ruminal bacteria community.

Supplementary Table 2. Identification results of ruminal metabolite responses to honeycomb flavonoids supplement.

References

1. Bagath M, Krishnan G, Devaraj C, Rashamol VP, Pragna P, Lees AM, et al. The impact of heat stress on the immune system in dairy cattle: a review. Res Vet Sci. (2019) 126:94–102. doi: 10.1016/j.rvsc.2019.08.011

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Gonzalez-Rivas PA, Chauhan SS, Ha M, Fegan N, Dunshea FR, Warner RD. Effects of heat stress on animal physiology, metabolism, and meat quality: a review. Meat Sci. (2020) 162:108025. doi: 10.1016/j.meatsci.2019.108025

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Becker CA, Collier RJ, Stone AE. Invited review: physiological and behavioral effects of heat stress in dairy cows. J Dairy Sci. (2020) 103:6751–70. doi: 10.3168/jds.2019-17929

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Takahashi M. Heat stress on reproductive function and fertility in mammals. Reprod Med Biol. (2012) 11:37–47. doi: 10.1007/s12522-011-0105-6

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Guo J, Gao S, Quan S, Zhang Y, Bu D, Wang J. Blood amino acids profile responding to heat stress in dairy cows. Asian Aust J Anim Sci. (2018) 31:47–53. doi: 10.5713/ajas.16.0428

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Visconti A, Le Roy CI, Rosa F, Rossi N, Martin TC, Mohney RP, et al. Interplay between the human gut microbiome and host metabolism. Nat Commun. (2019) 10:4505. doi: 10.1038/s41467-019-12476-z

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Claudie B, Fatoumata D, Roland B, Luke M, Edward T, Diarra MS. Pathotype and antibiotic resistance gene distributions of Escherichia coli isolates from broiler chickens raised on antimicrobial-supplemented diets. Appl Environ Microbiol. (2009) 75:6955–62. doi: 10.1128/AEM.00375-09

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Espley R, Butts C, Laing W, Martell S, Smith H, McGhie T, et al. Dietary flavonoids from modified apple reduce inflammation markers and modulate gut microbiota in mice. J Nutr. (2013) 144:146–54. doi: 10.3945/jn.113.182659

PubMed Abstract | CrossRef Full Text | Google Scholar

9. De Wasch K, Okerman L, De Brabander H, Van Hoof J, Croubels S, De Backer P. Detection of residues of tetracycline antibiotics in pork and chicken meat: correlation between results of screening and confirmatory tests. Analyst. (1998) 123:2737–41. doi: 10.1039/a804909b

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Pei R, Liu X, Bolling B. Flavonoids and gut health. Curr Opin Biotechnol. (2020) 61:153–9. doi: 10.1016/j.copbio.2019.12.018

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Li Y, Fang L, Xue F, Mao S, Jiang L. Effects of bamboo leaf extract on the production performance, rumen fermentation parameters, and rumen bacterial communities of heat-stressed dairy cows. Anim Biosci. (2021) 34:1784–93. doi: 10.5713/ab.20.0527

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Xue F, Sun F, Jiang L, Hua D, Wang Y, Nan X, et al. Effects of partial replacment of dietary forage using kelp powder (Thallus laminariae) on ruminal fermentation and lactation performances of dairy cows. Animals (Basel). (2019) 9:852. doi: 10.3390/ani9100852

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Wang Z, Yang DS, Li XY, Yu YN, Yong LY, Zhang PH, et al. Modulation of rumen fermentation and microbial community through increasing dietary cation-anion difference in Chinese Holstein dairy cows under heat stress conditions. J Appl Microbiol. (2021) 130:722–35. doi: 10.1111/jam.14812

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Emanuelson M, Murphy M, Lindberg JE. Effects of heat-treated and untreated full-fat rapeseed and tallow on rumen metabolism, digestibility, milk composition and milk yield in lactating cows. Anim Feed Sci Technol. (1991) 34:291–309.

Google Scholar

15. Baek YC, Choi H, Jeong JY, Lee SD, Kim MJ, Lee S, et al. The impact of short-term acute heat stress on the rumen microbiome of Hanwoo steers. J Anim Sci Technol. (2020) 62:208–17. doi: 10.5187/jast.2020.62.2.208

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Pinto S, Hoffmann G, Ammon C, Amon T. Critical THI thresholds based on the physiological parameters of lactating dairy cows. J Therm Biol. (2020) 88:102523. doi: 10.1016/j.jtherbio.2020.102523

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Xu L, Wen L, Ge Y, Wan G, Qu M, Xue F. Metagenomic insights into the effects of rare-earth elements supplementation on rumen digestibility and meat quality of beef cattle. Front Microbiol. (2020) 11:1933. doi: 10.3389/fmicb.2020.01933

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Xue F, Pan X, Jiang L, Guo Y, Xiong B. GC-MS analysis of the ruminal metabolome response to thiamine supplementation during high grain feeding in dairy cows. Metabolomics. (2018) 14:67. doi: 10.1007/s11306-018-1362-8

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Bolyen E, Rideout JR, Dillon MR, Bokulich NA, Abnet CC, Al-Ghalith GA, et al. Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. Nat Biotechnol. (2019) 37:852–7.

PubMed Abstract | Google Scholar

21. Ran M, Cha C, Xu Y, Zhang H, Yang Z, Li Z, et al. Traditional Chinese herbal medicine complex supplementation improves reproductive performance, serum biochemical parameters, and anti-oxidative capacity in periparturient dairy cows. Anim Biotechnol. (2020) 15:1–10. doi: 10.1080/10495398.2020.1819823

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Martens H, Leonhard-Marek S, Röntgen M, Stumpff F. Magnesium homeostasis in cattle: absorption and excretion. Nutr Res Rev. (2018) 31:114–30. doi: 10.1017/S0954422417000257

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Xue Y, Lin L, Hu F, Zhu W, Mao S. Disruption of ruminal homeostasis by malnutrition involved in systemic ruminal microbiota-host interactions in a pregnant sheep model. Microbiome. (2020) 8:138. doi: 10.1186/s40168-020-00916-8

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Correia Sales GF, Carvalho BF, Schwan RF, de Figueiredo Vilela L, Moreno Meneses JA, Gionbelli MP, et al. Heat stress influence the microbiota and organic acids concentration in beef cattle rumen. J Therm Biol. (2021) 97:102897. doi: 10.1016/j.jtherbio.2021.102897

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Zhong S, Ding Y, Wang Y, Zhou G, Guo H, Chen Y, et al. Temperature and humidity index (THI)-induced rumen bacterial community changes in goats. Appl Microbiol Biotechnol. (2019) 103:3193–203. doi: 10.1007/s00253-019-09673-7

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Patra AK, Kar I. Heat stress on microbiota composition, barrier integrity, and nutrient transport in gut, production performance, and its amelioration in farm animals. J Anim Sci Technol. (2021) 63:211–47. doi: 10.5187/jast.2021.e48

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Murota K, Nakamura Y, Uehara M. Flavonoid metabolism: the interaction of metabolites and gut microbiota. Biosci Biotechnol Biochem. (2018) 82:600–10. doi: 10.1080/09168451.2018.1444467

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Lock AL, Bauman DE. Modifying milk fat composition of dairy cows to enhance fatty acids beneficial to human health. Lipids. (2004) 39:1197–206. doi: 10.1007/s11745-004-1348-6

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Mackie RI, White BA, Bryant MP. Lipid metabolism in anaerobic ecosystems. Crit Rev Microbiol. (1991) 17:449–79. doi: 10.3109/10408419109115208

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Pengpeng W, Tan Z. Ammonia assimilation in rumen bacteria: a review. Anim Biotechnol. (2013) 24:107–28. doi: 10.1080/10495398.2012.756402

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Koch F, Thom U, Albrecht E, Weikard R, Nolte W, Kuhla B, et al. Heat stress directly impairs gut integrity and recruits distinct immune cell populations into the bovine intestine. Proc Natl Acad Sci USA. (2019)116:10333–8. doi: 10.1073/pnas.1820130116

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Maleki SJ, Crespo JF, Cabanillas B. Anti-inflammatory effects of flavonoids. Food Chem. (2019) 299:125124. doi: 10.1016/j.foodchem.2019.125124

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Chávez-Talavera O, Tailleux A, Lefebvre P, Staels B. Bile acid control of metabolism and inflammation in obesity, type 2 diabetes, dyslipidemia, and nonalcoholic fatty liver disease. Gastroenterology. (2017) 152:1679–94.e1673. doi: 10.1053/j.gastro.2017.01.055

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Conte G, Mele M, Chessa S, Castiglioni B, Serra A, Pagnacco G, et al. Diacylglycerol acyltransferase 1, stearoyl-CoA desaturase 1, and sterol regulatory element binding protein 1 gene polymorphisms and milk fatty acid composition in Italian Brown cattle. J Dairy Sci. (2010) 93:753–63. doi: 10.3168/jds.2009-2581

PubMed Abstract | CrossRef Full Text | Google Scholar