Background
The global temperature continues to rise; the hot and humid summer climate of Korea can cause heat stress (HS) in cattle. Beef cattle are particularly vulnerable to HS because they are unable to dissipate heat [1, 2]. Cattle exposed to HS consume less feed to prevent further increases in body temperature from metabolic heat generation [3, 4]. The loss of energy to maintain body temperature results in a lack of nutrients for growth and production [5], thereby maintaining cattle homeothermy at the expense of productivity and profitability [2]. Additionally, external indicators such as respiration rate and core body temperature are commonly used to determine HS damage [6, 7]. HS alters the gut microbiome, resulting in dysbiosis, impaired barrier and transport functions, and intestinal microstructures [8]. Moreover, at high temperatures, numerous reactive oxygen species (ROS), such as superoxide, hydrogen peroxide, and hydroxyl radicals, are generated in the body as by-products of the respiratory metabolic process [9]. ROS are capable to damage the antioxidant system in cattle and induce oxidative stress [10, 11]. At the cellular level, oxidative stress is an imbalance between oxidant species and antioxidants; cellular proteins, lipids, and DNA can be damaged by reactive oxygen, nitrogen, and other species formed from a compromised oxidative balance [12]. Various substances, such as minerals and vitamins, are known antioxidants that prevent oxidative cell damage; selenium, an antioxidant, affects the production of glutathione peroxidase, an antioxidant enzyme that prevents damage to cell membranes and peroxidation of fats and prevents cell function damage [13]. Zinc is also a component of superoxide dismutase, which scavenges superoxide, a component of ROS, in immune cells and affects protein and gene expression, animal reproduction, growth, and the immune system [14, 15]. Iron (Fe), a component of heme, is found primarily in hemoglobin and myoglobin; these proteins require iron as a cofactor for enzymes in the electron transport chain, specifically for myeloperoxidase, catalase, and cytochrome P-450. Iron deficiency results in hypochromic microcytic anemia due to failure in hemoglobin production, and the iron requirement is higher in young cattle than in mature cattle [16]. Iron also plays a beneficial role against oxidative stress in ruminants [17]. Traditionally, trace minerals have been supplied in the form of inorganic salts such as chlorides, sulphates, carbonates and oxides in ruminant diets, but their absorption rates are low and their excretion rates are high, which can cause environmental concerns [18, 19]. Supplying organic minerals can produce stable soluble molecules with high bioavailability and are better intestinal absorbed than inorganic minerals, due to the fact that organic form minerals are metal ions bound to organic substances such as peptides, amino acids, or polysaccharides and are then absorbed [20, 21]. Although the National Research Council (NRC) has advised the dose of the aforementioned minerals under normal circumstances [22], the ideal dose needed to counteract the negative effects of HS has not yet been established, which is a significant barrier for calf-raising. We hypothesized that high-concentration organic minerals (HOM) with antioxidant potential—as long as they don’t exceed the maximum permissible concentration advised by the NRC—can effectively and sustainably mitigate the negative effects of HS in calves during HS and the time afterwards. Our earlier study found that dairy steers with dietary higher mineral concentrations had lower levels of metabolic changes and oxidative damage related to HS [23]. Therefore, the objective of the present study was to evaluate the effect of inorganic and organic mineral supply on the rumen microbial community, blood oxidation state, and heat shock protein levels under stressful and non-stressful temperature conditions.
Materials and methods
Animal care
The study was conducted at the animal farm of the Sunchon National University (SCNU) and in the Ruminant Nutrition and Anaerobe Laboratory of the Department of Animal Science and Technology at SCNU, located in Jeonnam, South Korea. The handling of all experimental animals and relevant protocols were performed in accordance with the guidelines of the SCNU Institutional Animal Care and Use Committee (IACUC approval number: SCNU-IACUC-2020-06).
Animals, experimental design, and diet
The experimental design was a 4-period crossover with a 4 × 2 factorial arrangement involving 4 treatment groups (Con, IM, OM, and HOM) and two thermal environments (HS and recovery). The experiment trial was conducted for 140 d (July 2021 to November 2021) with each period lasting 35 d, including 14 d of HS, 14 d of recovery (thermal neutral) condition and a 7-d washing period (Fig. 1). Eight Holstein bull calves (123.81 ± 9.76 kg; 5 months old) were used in the study. The calves were randomly distributed into four groups of 2 calves each to evaluate the four treatments. Calves were kept in individual pens in a temperature-controlled barn. Calves were fed once daily (0900 h) with Timothy hay and concentrate feed (CJ Feed & Care, Korea) separately at a ratio of 6:4 and had free access to water. The treatments included a control (without mineral supplementation), inorganic minerals (IM), organic minerals (OM), and high-concentration organic minerals (HOM), and the top dressing was performed at 0.4% of the basal diet (Table 1). Inorganic minerals used were ZnSO4 and MnSO4 (Seoan Chemtec Co., Ltd., Korea), FeSO4, CuSO4 and Na2SeO3 (TMC Co., Ltd., Korea), Ca(IO3)2 and CoSO4 (Innobatech Co., Ltd., Korea), while the organic minerals used were amino acid complex Zn, Mn, Fe, Cu, I, Co, and Se (Prochem Co., Ltd., Korea). The ingredients and chemical composition of the basal diet are listed in Table 2 which were formulate according to NRC [22]. Feed intake was determined daily by weekly collection of feed residues and ingredients. The chemical composition of the basal diet was analyzed using standard methods [24]. Neutral detergent fiber (NDF) and acid detergent fiber (ADF) contents were determined according to the protocols described by Van Soest et al. [25] and Van Soest [26], respectively.

Timeline describing overall study design, data, and sample collection for each period. IBW, Initial body weight; RR, Respiration rate; RT, Rectal temperature; BT, Body temperature; ADG, Average daily gain
The ambient temperature (°C) and relative humidity (%) of the experimental shed were recorded throughout the experimental period using a Testo 174H Mini data logger (West Chester, PA, USA) and are shown in Fig. 2. The temperature humidity index (THI) was calculated as THI = (0.8 × ambient temperature) + [(% relative humidity/100) × (ambient temperature − 14.4)] + 46.4 [27].

Recorded ambient temperature, relative humidity, and temperature humidity index during study periods. RH, relative humidity; THI, temperature humidity index
Sample collection and measurements
Body weight before feeding was measured on d 14 and 28 of each period, and 15 mL of blood was collected from the jugular vein of each calf 2 h after feeding. The collected blood was immediately transferred to 3 vacutainers (SSTTM II Advance, Trace Element serum, K2 EDTA; BD vacutainer®) of 5 mL each and centrifuged at 3,000 r/min for 15 min, and the serum and plasma samples were transferred into 2-mL microtubes and stored at −80 °C until analysis. After blood collection, rumen fluid samples were collected using a stomach tube, and pH was immediately measured using a pH meter (Seven CompactTM pH/Ion meter S220, Mettler Toledo, Switzerland). Simultaneously, the rumen fluid was transported to the laboratory using dry ice and stored at −80 °C until ammonia nitrogen (NH3-N), volatile fatty acids (VFAs), and rumen microbiota could be analyzed. In addition, the respiration rate, rectal temperature, and body temperature were measured twice a day (09:00 and 15:00) for three consecutive days starting from d 11 and 25 in the HS and recovery periods, respectively. The respiration rate was converted to the respiration rate per minute by multiplying by 2 after observing the movement of the flank for 30 s. The rectal temperature was recorded 30 s after the digital thermometer (WPT-1; CAS, Korea; accuracy ± 1 ºC) was inserted into the rectum, and the body temperature was measured using an infrared thermometer (SO-20-01; Apollo, Korea; accuracy ± 0.2 ºC).
Sample analysis
Analysis of serum biochemistry, mineral concentrations, cortisol, total oxidant status, total antioxidant status, and oxidative stress index
Serum calcium (Ca), magnesium (Mg), phosphorus (P), blood urea nitrogen (BUN), total protein (TP), aspartate aminotransferase (AST), total bilirubin (TBIL), and cholesterol (CHOL) levels were analyzed using a biochemical analyzer (IDEXX Catalyst One, USA). Serum glucose levels were analyzed using a portable glucose test meter (FreeStyle Optium Neo H Blood Glucose and Ketone System; Abbott Diabetes Care Ltd., UK). Serum concentrations of Zn and Se were analyzed using an inductively coupled plasma mass spectrometry (ICP-MS) system (ELAN DRCe, PerkinElmer, Germany), and Fe was analyzed using a Cobas 8000 (Roche, Germany). Serum cortisol levels were measured using a bovine cortisol enzyme-linked immunosorbent assay (ELISA) kit (cat. MBS2608983, Mybiosource), and total oxidant status (TOS) and total antioxidant status (TAS) were measured using the Rel Assay kit (Rel Assay Diagnostics, Turkey) following the manufacturer’s instructions. The oxidative stress index (OSI) was calculated according to the formula OSI = [(TOS, mol/L)/(TAS, mmol Trolox equivalent/L) × 100] [28].
Analysis of plasma superoxide dismutase, glutathione peroxidase, malondialdehyde, and heat shock proteins
Plasma superoxide dismutase (SOD, Item # 706002) and glutathione peroxidase (GPX, Item # 703102) levels were measured using an ELISA kit (Cayman Chemical Company, Ann Arbor, USA) following the manufacturer’s instructions, and malondialdehyde (MDA) was measured using an MDA Assay kit (ab238537, Abcam) according to the manufacturer’s instructions. HSP levels (HSP27, HSP70, and HSP90) were measured using MyBioSource (San Diego, CA, USA) bovine ELISA kits (MBS011935, MBS2882245, and MBS094979) following the manufacturer’s instructions.
Analysis of NH3-N and VFA concentrations
The NH3-N concentration was measured colorimetrically using a Libra S22 spectrophotometer (CB40FJ; Biochrom Ltd., Cambourne, UK), following the protocol described by Chaney and Marbach [29]. VFA concentration was measured according to the method described by Han et al. [30]. and Tabaru et al. [31] using high-performance liquid chromatography (HPLC; Agilent Technologies 1200 series, Waldbronn, Germany). A UV detector (set at 210 nm and 220 nm), METACARB87H column (Varian, Palo Alto, CA, USA), and buffered solvent (0.00425 mol/L H2SO4) at a flow rate of 0.6 mL/min were used to perform HPLC.
DNA extraction and PCR amplification and 16S rRNA amplicon sequencing
All rumen fluid samples were sent to Macrogen, Inc. (Seoul, Korea) for DNA extraction and metataxonomic analysis of the rumen microbiota. Briefly, DNA was extracted using a DNeasy PowerSoil kit (Qiagen, Hilden, Germany), following the manufacturer’s protocol [32]. The quality and quantity of the DNA were assessed using PicoGreen and NanoDrop, respectively. Sequencing libraries were prepared according to the Illumina 16S Metagenomic Sequencing Library protocols to amplify the V3 and V4 regions. The PCR assay was performed with 2 ng of gDNA, 1× reaction buffer, 1 nmol/L of dNTP mix, 500 nmol/L of each of the universal F/R PCR primers, and 2.5 U of Herculase II fusion DNA polymerase (Agilent Technologies, Santa Clara, CA, USA). The cycle condition for 1st PCR was 3 min at 95 °C for heat activation, and 25 cycles of 30 s at 95 °C, 30 s at 55 °C, and 30 s at 72 °C, followed by a 5 min final extension at 72 °C. The universal primer pair with Illumina adapter overhang sequences used for the first amplification were as follows: V3-F: 5′-TCG TCG GCA GCG TCA GAT GTG TAT AAG AGA CAG CCT ACG GGN GGC WGC AG-3′, V4-R: 5′-GTC TCG TGG GCT CGG AGA TGT GTA TAA GAG ACA GGA CTA CHV GGG TAT CTA ATC C-3′. AMPure beads (Agencourt Bioscience, Beverly, MA, USA) were used to purify the products of the first and second PCR. Individual amplicon libraries were normalized after quantification using PicoGreen. They were then size-verified using a TapeStation DNA ScreenTape D1000 (Agilent Technologies), pooled at an equimolar ratio, and sequenced on a MiSeq system (Illumina, San Diego, CA, USA) using a 2 × 300 bp kit.
Sequence data processing and metataxonomic analysis
After sequencing, the raw Illumina MiSeq data were classified by sample using an index sequence, and a paired-end FASTQ file was generated for each sample. Sequencing adapter and F/R primer sequences of the target gene region were removed using Cutadapt (v3.2) [33].
The DADA2 (v1.18.0) [34] package of R (v4.0.3) was used for error correction in the amplicon sequencing process. For paired-end reads, the forward sequence (Read1) and reverse sequence (Read2) were cut to 250 bp and 200 bp, respectively, and sequences with expected errors of two or more were excluded. An error model was established for each batch to remove noise from each sample. After assembling the paired-end sequence corrected for sequencing errors into one sequence, the chimeric sequence was removed using the DADA2 Consensus method to form amplicon sequence variants (ASVs). In addition, for a comparative analysis of the microbial community, the QIIME (v1.9) [35] program was used for normalization by applying subsampling based on the number of reads of the sample with the minimum number of reads among all samples.
For each ASV sequence, BLAST + (v2.9.0) [36] was performed using the reference database (DB; NCBI 16S Microbial DB), and the taxonomic information for the organism of the subject with the highest similarity was assigned. At this time, if the query coverage of the best-hit matching in the DB was less than 85% or the identity of the matched area was less than 85%, taxonomy information was not allocated.
Various microbial communities were comparatively analyzed using QIIME with the aforementioned ASVs abundances and taxonomic information. To analyze beta diversity indices, we used principal coordinate analysis (PCoA) based on both unweighted and weighted UniFrac distances for visualization [37]. Before analysis, we assessed the homogeneity of multivariate dispersion among treatments and groups using permutational multivariate analysis of dispersion (PERMDISP) [38]. Subsequently, a permutational analysis of variance (PERMANOVA) was performed to detect any potential differences in the microbial community structures between the minerals and temperatures. All procedures, including PCoA, PERMDISP, and PERMANOVA, were conducted using the vegan package in R software (v4.0.3).
Statistical analysis
Before analysis, the homoscedasticity and normality of the distribution assumptions were tested. Growth performance, respiration rate, body temperature, blood parameters, and rumen fermentation parameters were analyzed by repeated measures using MIXED procedure of the SAS statistical package (version 9.4; SAS Institute Inc., Cary, NC, USA) with an autoregressive covariance structure. The model included the main fixed effects of temperature, minerals, and interactions between the two variables, whereas the effect of each calf was considered random [39]. Statistical significance was set at P < 0.05 unless otherwise stated. Following a significant effect, Duncan’s post-hoc test was performed to assess the differences between the means.
Results
Growth performance, body temperature, and respiration rate
The growth performance of Holstein calves supplemented with inorganic and organic minerals at different temperatures is shown in Table 3. The average daily gain (ADG; kg), dry matter intake (DMI; kg/d), and FE were not significantly different between the treatment groups. The body temperatures and respiration rates of Holstein calves supplemented with inorganic and organic minerals at different temperatures are presented in Table 4 and Fig. 3. Body temperature was higher in the HS condition than in the recovery condition for all parameters (rectal, flank, rump, perineum, belly, back, and thigh) at 15:00 h (P < 0.05). The respiration rate was also higher in the HS condition than during the recovery period at both time points (P < 0.05).

Respiration rate of Holstein calves supplemented with inorganic and organic minerals at different temperatures. Con: Basal diet (without mineral supplementation), IM: Basal diet + inorganic minerals supplementation, OM: Basal diet + organic minerals supplementation, HOM: Basal diet + higher organic minerals supplementation
Serum biochemistry, mineral concentrations, cortisol, TOS, TAS, and OSI
The serum biochemistry and mineral concentrations of Holstein calves supplemented with inorganic and organic minerals at different temperatures are presented in Table 5. None of the tested serum biochemical parameters were significantly different between the different temperature and treatment groups. Serum Se concentration was higher in the recovery condition than in the HS condition (P < 0.05), and in the HOM group, it was the highest at 90.38 μg/dL and 102.00 μg/dL in the heat and recovery conditions, respectively (P < 0.05). The results for serum cortisol, TOS, TAS, and OSI in Holstein calves supplemented with inorganic and organic minerals at different temperatures are presented in Fig. 4. The serum cortisol concentration was 20.26 ng/mL in the HOM group in the HS condition, 5.60 ng/mL lower than that in the control group (P < 0.05). The TAS was the highest with 2.71 mmol Trolox equivalent/L in the OM group, followed by the HOM group under HS conditions and was the highest at 2.58 mmol Trolox equivalent/L in the HOM group under the recovery conditions (P < 0.05). The OSI was lower in the mineral-supplemented groups than in the Con in both periods, and the lowest OSI was observed in the HOM group (149.61) under recovery conditions (P < 0.05).

Serum cortisol (A) and total oxidant status (TOS), total antioxidant status (TAS), and oxidative stress index (OSI) (B) of Holstein calves supplemented with inorganic and organic minerals at different temperatures. Con: Basal diet (without mineral supplementation), IM: Basal diet + inorganic minerals supplementation, OM: Basal diet + organic minerals supplementation, HOM: Basal diet + higher organic minerals supplementation
Plasma SOD, GPX, MDA, and HSPs
The results of plasma SOD, GPX, MDA, and HSPs in Holstein calves supplemented with inorganic and organic minerals at different temperatures are presented in Table 6. The plasma SOD and GPX concentrations were numerically lower, whereas HSPs were numerically higher in the mineral-supplemented groups than in the Con group under both heat and recovery conditions (P > 0.05). The HOM group under heat and recovery conditions showed lower plasma MDA and HSP70 levels and higher SOD and GPX activities (P > 0.05).
Rumen fermentation characteristics and composition of rumen microbiota
The rumen fermentation characteristics of Holstein calves supplemented with inorganic and organic minerals at different temperatures are presented in Table 7. No significant differences were observed in pH, NH3-N, or total VFA among the different temperatures and treatment groups. Figure 5A and B show the beta diversity results. Individuals were projected onto the first two dimensions of PCoA, which accounted for 9.1% and 6.9% of the variability in unweighted UniFrac distances, respectively. Similarly, the first two dimensions of the PCoA accounted for 34.1% and 20.1% of the variability in the weighted UniFrac distances observed between the samples (Fig. 5C and D). Regardless of the factors investigated (temperature and minerals) or the distance matrix utilized (unweighted or weighted UniFrac), no significant differences in the individual dispersions were observed (PERMDISP: P > 0.05). Moreover, the PERMANOVA procedure failed to detect substantial differences (P > 0.05) between the minerals using either unweighted or weighted UniFrac distances. Conversely, PERMANOVA revealed a significant influence of temperature on the overall microbial structure, based on both unweighted (P = 0.023) and weighted (P = 0.038) UniFrac distances. The ruminal bacterial abundances in Holstein calves supplemented with inorganic and organic minerals at different temperatures are presented in Fig. 6. Upon evaluating the composition of the rumen microbiota, at the phylum level, Bacteroidetes (accounting for 40.86% to 57.87%) and Firmicutes (32.40% to 45.96%) were the two major bacterial taxa among all treatment groups. Firmicutes and Actinobacteria decreased, whereas Fibrobacteres, Spirochaetes, and Tenericutes increased under HS conditions (Fig. 6A). At the genus level, Treponema increased during HS, whereas Christensenella increased during the recovery period (Fig. 6B).

PCoA of rumen microbiota of Holstein calves supplemented with inorganic and organic minerals at different temperatures. A Unweighted UniFrac distances (minerals). B Unweighted UniFrac distances (temperature). C Weighted UniFrac distances (minerals). D Weighted UniFrac distances (temperature). Con (C): Basal diet (without mineral supplementation), IM: Basal diet + inorganic minerals supplementation, OM: Basal diet + organic minerals supplementation, HOM: Basal diet + higher organic minerals supplementation

Rumen bacterial abundance of Holstein calves supplemented with inorganic and organic minerals at different temperatures. A Phylum level. B Genus level. Con: Basal diet (without mineral supplementation), IM: Basal diet + inorganic minerals supplementation, OM: Basal diet + organic minerals supplementation, HOM: Basal diet + higher organic minerals supplementation
Discussion
Cattle exposed to HS show reduced feed intake, increased water intake, and changes in endocrine status, resulting in decreased body weight and average daily gain, and worsening body condition due to increased maintenance requirements [40]. However, non-significant variations in ADG (kg) and FE were observed among the treatment groups, which may be due to the short-term HS conditions used in this study. Body temperature and respiration rate (RR) were unaffected by mineral supply but increased under HS conditions according to the change in THI. This is similar to the study in which rectal temperature and RR were considered visible physiological indicators for determining the thermal regulation response of dairy cows to heat stress and that they are altered due to high environmental temperatures [41, 42].
Among these minerals, selenium, iron, and zinc play important roles in animal growth and health by being involved in the antioxidant defense system [17, 43, 44]. Increased lipid peroxidation under stress conditions stimulates the stress axis to increase cortisol concentrations [45, 46]. In this study, the supply of organic minerals (OM and HOM groups) was beneficial in reducing cortisol levels and oxidative stress index during short-term HS.
The MDA concentration is an indicator of lipid peroxidation [47], and antioxidant enzymes, such as SOD and GPX, can be synthesized in the body to remove ROS generated during HS [48, 49]. HSP70 expression has been used as an indicator of heat stress in animals [50]. In this study, the mineral supplementation and temperature did not significantly affect the MDA and HSP70, while a tendency for an increase in SOD levels was observed in HOM under heat and recovery conditions. This suggests that high concentrations of organic minerals are beneficial in alleviating the side effects of HS.
No significant differences were observed in pH and NH3-N according to temperature changes or mineral supply, which explains why the supply of high concentrations of organic minerals did not negatively affect rumen fermentation. HS decreases the proportion of acetate in the rumen, and increases butyrate proportion and propionate concentration [51, 52]. In the present study, non-significant differences in the molar proportions of acetate, propionate, and butyrate were recorded, which might be due to the short-term HS conditions.
In Holstein rumen microbiota under heat stress, the increase in Fibrobacteres and the soluble carbohydrate-digesting bacterium Treponema and the decrease in the acetate-producing bacteria Actinobacteria were consistent with previous studies [53, 54]. Additionally, the relative abundance of Christensenella increased as the mineral concentration increased under HS conditions. According to Hu et al. [55] and Correia Sales et al. [56], Christensenella mainly affects rumen energy metabolism, and becomes relatively abundant when high-energy feed is consumed. This suggests that the addition of organic minerals may alter the composition of rumen bacterial communities. Although the overall rumen microbiome was not influenced by mineral supplementation, a significant alteration was observed between the HS and recovery periods, which is consistent with earlier studies [57].
Conclusion
A high concentration of organic mineral supplementation during HS reduced the cortisol concentration and increased the total antioxidant status in Holstein bull calves without altering the overall rumen microbiota. Overall, these results suggest that the side effects of heat stress can be alleviated by adding HOM.