Skip to main content

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



Heat stress (HS) affects the ruminal microbiota and decreases the lactation performance of dairy cows. Because HS decreases feed intake, the results of previous studies were confounded by the effect of HS on feed intake. This study examined the direct effect of HS on the ruminal microbiota using lactating Holstein cows that were pair-fed and housed in environmental chambers in a 2 × 2 crossover design. The cows were pair-fed the same amount of identical total mixed ration to eliminate the effect of feed or feed intake. The composition and structure of the microbiota of prokaryotes, fungi, and protozoa were analyzed using metataxonomics and compared between two thermal conditions: pair-fed thermoneutrality (PFTN, thermal humidity index: 65.5) and HS (87.2 for daytime and 81.8 for nighttime).


The HS conditions altered the structure of the prokaryotic microbiota and the protozoal microbiota, but not the fungal microbiota. Heat stress significantly increased the relative abundance of Bacteroidetes (primarily Gram-negative bacteria) while decreasing that of Firmicutes (primarily Gram-positive bacteria) and the Firmicutes-to-Bacteroidetes ratio. Some genera were exclusively found in the heat-stressed cows and thermal control cows. Some co-occurrence and mutual exclusion between some genera were also found exclusively for each thermal condition. Heat stress did not significantly affect the overall functional features predicted using the 16S rRNA gene sequences and ITS1 sequences, but some enzyme-coding genes altered their relative abundance in response to HS.


Overall, HS affected the prokaryotes, fungi, and protozoa of the ruminal microbiota in lactating Holstein cows to a different extent, but the effect on the structure of ruminal microbiota and functional profiles was limited when not confounded by the effect on feed intake. However, some genera and co-occurrence were exclusively found in the rumen of heat-stressed cows. These effects should be attributed to the direct effect of heat stress on the host metabolism, physiology, and behavior. Some of the “heat-stress resistant” microbes may be useful as potential probiotics for cows under heat stress.


Global warming as a result of the increased release of greenhouse gases (GHG) from different sources, including agriculture and livestock, is of great concern worldwide. Animal production, including dairy production, is impaired considerably by rising environmental temperature. The elevated ambient temperature accompanying global warming has been increasing the frequency and duration of heat stress (HS) episodes in dairy cows, especially in tropical, subtropical, and Mediterranean regions [1]. Dairy cows are more susceptible to HS than other farm animals, and they suffer from HS when the average temperature-humidity index (THI) exceeds 68 (approximately 22 °C at 50% relative humidity) [2]. Heat stress can have profound adverse effects on many aspects of dairy cows, including decreasing feed intake, feed efficiency, nutrient digestibility, milk production, milk quality; creating negative energy balance; and impairing reproductive performance among others [3,4,5,6]. Any of these adverse effects can lead to significant economic loss to the producers [7]. Moreover, HS also leads to poor animal welfare [8]. Heat stress in dairy cows is worsening as global warming continues, and HS is now well recognized as an environmental stressor that undermines the sustainable development of the dairy industry [9]. Therefore, HS in dairy cows has attracted tremendous research interest in the past decade. It has been shown that HS can directly impair the normal metabolism, physiology, and immune system in lactating dairy cows, including the metabolism of lipid in bovine primary adipocytes [10], the metabolism of carbohydrate and lipid, and milk protein synthesis in mammary tissues [11, 12], the metabolism of glucose and energy and energy balance [13], and immune system function [14]. Reduced voluntary feed intake is inherently associated with HS, and it aggravates these adverse effects indirectly [15,16,17].

The complex and diverse ruminal microbiota [18,19,20] plays a pivotal role in supplying the majority of the energy and nutrients required by dairy cows [21]. This microbiota is also dynamic responding to many internal and external factors [22,23,24], including HS [25, 26]. Heat stress can affect the ruminal microbiota because it reduces feed intake and impacts the metabolism, physiology, immune system, and behavior of the host. Using 16S rRNA or its gene-based analysis, the earliest studies showed that HS significantly affected the species richness and composition of the ruminal microbiota in dairy heifers with a BW ≥ 430 kg [27] and increased the genus Streptococcus and Clostridium coccoidesEubacterium rectale group but decreased Fibrobacter in the rumen of dairy heifers [28]. In recent years, metataxonomics and metagenomics have been used to comprehensively examine how the ruminal microbiota responds to HS. Zhao et al. [25] found that HS did not significantly alter the bacterial community structure but affected some bacteria, including lactate-producing bacteria (increase) and acetate-producing bacteria (decrease) in the rumen of lactating dairy cows. Another recent study showed that increasing THI did not significantly affect the α-diversity of the ruminal microbiota in goats but shifted the population of some rumen bacteria, including an increase of Bacteroidetes at the expense of Firmicutes [29]. This study also showed that elevated THI could affect the function of the ruminal microbiota as predicted using PICRUSt. Baek et al. [30] evaluated the effect of HS on the ruminal microbiota of Hanwoo steers, and they showed delayed ruminal microbiota response: significant effects after six days but not after three days. In a comparative study [26], HS was shown to affect the ruminal microbiota and its predicted function in lactating Holstein cows and Jersey cows to a different extent. The reported studies provided useful insights into the effects of HS on the ruminal microbiota, but the responses of the ruminal microbiota to HS varied considerably.

Many factors, such as THI, feed and feeding, species or breed of animals, can affect the responses of the ruminal microbiota to HS. In particular, HS decreases feed intake substantially [16, 17], and feed intake can profoundly affect the ruminal microbiota [31,32,33]. To the best of our knowledge, the animals were fed ad libitum in all the published studies that examined the effect of HS on the ruminal microbiota. Therefore, the decreased feed intake associated with heat-stressed animals but not control animals is a significant confounding factor that prevents the appraisal of the direct HS effect on the ruminal microbiota. We hypothesized that the actual effect of HS on the ruminal microbiota could only be determined when cows were fed the same diet and at the same feed intake. We also hypothesized that HS could also affect rumen methanogens, protozoa, and fungi differently, none of which has been targetedly analyzed in heat-stressed ruminant animals. In the present study, we tested the above hypotheses using lactating cows that were pair-fed the same diet using metataxonomics targeting all the kingdoms of the ruminal microbiota. Our results provided the first insights into how HS could affect the ruminal microbiota independent of the indirect effect on feed intake.


Animals and heat stress treatments

The animal trial of heat stress has been reported previously [11]. Briefly, animals were cared for per the guidelines of the institutional animal care and use committee of the Institute of Animal Science, Chinese Academy of Agricultural Sciences. Four multiparous, lactating Holstein cows (101 ± 10 d in milk, 574 ± 36 kg of body weight, 38 ± 2 kg of milk/d, second parity, 1–2 month pregnant) were individually housed in four environmental chambers (temperature fluctuation: ± 0.5 °C, relative humidity fluctuation: ± 5%). The chambers had 12 h light:12 h dark cycles. The cows were initially fed twice daily (0500 and 1700 h) a total mixed ration (TMR: hay, 24.0%; whole corn silage, 26.0%; steam-flaked corn, 22.0%; extruded soybean, 2.1%; bean pulp, 11.4%; rapeseed and corn meal, 12.4%; limestone and salt, 1.2%; salt, 0.38%; and supplement 0.6%) ad libitum and allowed to adapt to the chamber (maintained at thermal neutral conditions: 20 °C, 40% humidity) for 9 d. Then the four cows were randomly allocated to 1 of 2 treatments in a 2 × 2 crossover design: PFTN (20 °C, 55% humidity; THI = 65.5) and cyclic HS conditions (55% humidity; 06:00–18:00 h at 36 °C while 18:00–06:00 h at 32 °C; THI = 87.2 and 81.8, respectively). THI was calculated as THI = (0.8 × air temperature in °C) + relative humidity × (air temperature °C – 14.3) + 46.4 [34, 35]. All the cows were fed the same TMR formulated to meet or exceed the predicted requirements [36] of energy, protein, minerals, and vitamins. The animals in the HS group were fed ad libitum, while those in the PFTN group were fed the amount based on the intake of the HS group of the previous day. During the experiment, all the animals had free access to drinking water, and they were milked twice daily (05:00 and 17:00 h). Each of the two experimental periods lasted for 18 d, and a 30-d wash-out period (thermoneutral conditions at 20 °C and 55% humidity, THI = 65.5, with ad libitum feeding) separated the two experimental periods. Rumen fluid samples were collected from individual cows three hours after the morning feeding on d 9 of each treatment period using stomach tubing (the first 50 mL of rumen fluid was discarded to avoid saliva contamination) [37]. The rumen samples were immediately stored at − 80 °C until metagenomic DNA extraction.

Metagenomic DNA extraction and metataxonomic analysis of the ruminal microbiota

Metagenomic DNA was extracted from each rumen sample using the RBB + C method [38]. The quality and quantity of the extracted DNA were determined using a NanoDrop ND-2000 spectrophotometer (Thermo Scientific, NanoDrop Technologies, Wilmington, DE, USA) followed by agarose gel (1%, w/v) electrophoresis. The rumen prokaryotic, protozoal, and fungal microbiotas were analyzed using metataxonomics and separate amplicon libraries prepared using PCR primers specific for the V4 region of 16S rRNA genes of bacteria and archaea, ITS1 of fungi, and the V3-V4 region of the 18S rRNA gene of protozoa, respectively [22]. The primer information and respective PCR conditions were presented in our previous study [22]. The amplicon libraries containing unique barcodes for multiplexing were pooled and sequenced using the 2 × 300 paired-end protocol on an Illumina MiSeq platform at the Molecular and Cellular Imaging Center of The Ohio State University following the Illumina protocol [39].

The sequence reads were analyzed using QIIME2 (version 2019.4) [40] as done in a previous study [22]. Briefly, after trimming off the adapter and primer sequences using Cutadapt [41], DADA2 was used to perform quality filtering (Q-score > 25), merging of paired-end reads, and removal of chimeric sequences [42]. For the fungal ITS1 sequence reads, ITSxpress [43] was used to trim the ITS1 region from paired-end reads resulting in a single merged FASTQ file followed by using the q2-dada2 denoise-single method [42] for further dereplication and chimeric sequence removal. For prokaryotic sequences, further taxonomic filtration was done to remove the amplicon sequence variants (ASVs) labeled as ‘unassigned’, ‘chloroplast’, or ‘mitochondria’. The resulting prokaryotic and fungal ASVs from each denoised multi-kingdom library were taxonomically assigned using respective taxonomic classifiers that had been manually trained using the Naïve Bayes classifier [44] with the Greengenes 16S reference database (13_8 version) [45] for bacteria and archaea, UNITE ITS reference sequences (v.8.0 11.18.2018 database) [46] for fungi. Protozoal ASVs were classified using BLASTn against the NCBI nucleotide collections excluding uncultured/environmental sample sequences (accessed March 6, 2021), and the best BLASTn hits were selected for assigning the taxonomy of each protozoal ASV. Prevalent phyla and genera found in at least 50% of the rumen samples were mainly discussed in this study.

α-Diversity metrics including species richness (Chao1 estimates, observed ASVs, observed genera, and observed species), evenness, Faith’s phylogenetic diversity, Shannon’s index, and Simpson’s index were calculated based on the rarefied ASV BIOM tables. Comparison of the overall microbiotas between the PFTN control cows and the HS cows were done using principal coordinates analysis (PCoA) based on weighted UniFrac distances [47]. The number of shared genera between the two treatments and exclusively found in each treatment were visualized using a Venn diagram. Multi-kingdom co-occurrence and mutual-exclusion relationships were analyzed based on prevalent microbial genera (present in at least 50% of the rumen samples) using FastSpar by computing the correlation with the SparCC algorithm for compositional data ( [48]. Of the microbial interactions with tendency (P ≤ 0.1), the microbial networks specific to each treatment were defined and visualized using co-expression differential network analysis (CoDiNA) [49].

Prediction of functional profiles of the prokaryotic and fungal microbiotas

Microbial functional profiles were predicted from 16S- and ITS1-based ASVs using PICRUSt2 [50] using the reference genomes of the Integrated Microbial Genomes (IMG) database implemented in PICRUSt2. The normalized counts of predicted Enzyme Commission (EC) numbers for both prokaryotic and fungal microbiotas were used to define the overall functional dissimilarities between the PFTN control cows and the HS cows and then analyzed using principal components analysis (PCA) based on the Bray-Curtis dissimilarity metrics [51]. The PCA results were visualized as plots using the ggfortify package of R (3.5.3) [52].

Statistical analysis

α-Diversity metrics of each kingdom microbiota (prokaryotic, fungal, or protozoal), Firmicutes-to-Bacteroidetes ratio, and the number of PICRUSt2-predicted functional features were statistically analyzed with treatment (i.e., thermal condition), period, and treatment × period interaction as fixed effects and cow as random effect using the GLIMMIX procedure of SAS 9.4 (SAS Institute Inc., Cary, NC, USA). The comparative analysis results of PCoA between the two thermal conditions were further assessed using permutational multivariate analysis of variance (PERMANOVA) with 999 random permutations. Differentially abundant microbial phyla, genera, and EC numbers between the PFTN control cows and HS cows were defined using linear discriminant analysis effect size (LefSe) [53] with logarithmic LDA score 2 as the cutoff. Statistical significance was declared at P ≤ 0.05 and tendency at 0.05 < P ≤ 0.10.


As presented in our previous study [11], the HS conditions significantly increased respiratory rates (33.3 bpm for FFTN vs. 77.8 bpm for HS, P < 0.0001) and rectal temperature (38.5 °C for FFTN vs. 40.0 °C for HS, P = 0.0026). By using the PFTN group as control to eliminate the confounding difference in DMI, comparison between PFTN and HS allowed for the revelation of the true effect of heat stress on the multi-kingdom rumen microbiome.

After denoising and quality filtration, on average > 12,000 sequences were each obtained for prokaryotes and fungi, and > 21,000 sequences were retained for protozoa for each sample (Additional file 1: Table S1). Based on the fewest numbers of sequences for each kingdom among all the samples, 8941, 10,691, and 12,120 sequences per sample were rarefied and analyzed for prokaryotes, fungi, and protozoa, respectively. Good’s coverage for kingdoms and all the samples reached over 99.9% (data not shown).

Heat stress affects the diversity and co-occurrence interactions of the ruminal microbiota

Heat stress did not affect any of the α-diversity metrics for prokaryotic or fungal microbiota, but it tended to lower the evenness (P = 0.057) and Shannon’s index (P = 0.075) of the protozoal microbiota compared to the pair-fed thermal neutral (PFTN) control cows, while more protozoal genera (P = 0.004) were observed in HS (Table 1). Period had some minor effect on the evenness of the prokaryotic microbiota (P = 0.045), observed species richness (P = 0.085), and Faith’s phylogenetic diversity (P = 0.084) of fungi. The evenness (P = 0.001), Shannon’s index (P = 0.001), and Simpson’s index (P = 0.009) of protozoal microbiota were all affected by period (Table 1).

Table 1 α-Diversity metrics of ruminal multi-kingdom microbiota in heat-stressed (HS) cows and pair-fed thermoneutral (PFTN) control cows

Based on PCoA and PERMANOVA analyses, compared to PFTN, HS altered the overall prokaryotic microbiota with a trend toward statistical significance (P = 0.068) and the protozoal microbiota significantly (P = 0.045), but it did not affect (P = 0.207) the overall fungal microbiota (Fig. 1). Co-occurrence network analysis using CoDiNA identified 17 and 10 significant interactions among 150 overall edges exclusively found in the ruminal microbiotas of the PFTN control and the HS cows, respectively (Fig. 2A). Three genus-level nodes were found only in the PFTN control cows; they included two bacterial genera: Ruminococcus and Selenomonas, and one protozoal genus Isotricha. Only two genus-level nodes were exclusively found in the HS cows: one archaeal genus Methanobrevibacter and one bacterial genus Treponema. No fungal nodes were exclusively found in either the PFTN control cows or the HS cows. Of the three genus nodes specific for the PFTN control cows (Fig. 2B), Ruminococcus was mutually exclusive with three bacterial genera (i.e., Prevotella, YRC22, and RFN20) and one protozoal genus (i.e., Diploplastron), but it showed co-occurrence with Aspergillus. On the other hand, Selenomonas was mutually exclusive with Aspergillus, and it had co-occurrence with four genera: YRC22, Succiniclasticum, Neocallimastix, and Diploplastron. The only specific genus node, Isotricha, had seven significant interactions, including co-occurrence with three genera of bacteria (i.e., Succiniclasticum, Sphaerochaeta, and Treponema) and three genera of fungi (i.e., Neocallimastix, Piromyces, and Caecomyces), and mutual exclusion with Kazachstania (a genus of yeast). In the specific network of HS cows (Fig. 2B), the exclusive node Treponema showed co-occurrence with three rumen fungal genera (Caecomyces, Neocallimastix, and Piromyces) and three bacterial genera (Prevotella, CF231, and Succiniclasticum). Another exclusive node of the network of HS cows, Methanobrevibacter, was mutually exclusive with two Bacteroidetes genera (Prevotella and CF231) while co-occurring with Mogibacterium. Those two exclusive nodes (Treponema and Methanobrevibacter) were mutually exclusive.

Fig. 1
figure 1

Principal coordinates analysis (PCoA) plots (based on weighted UniFrac distance metrics) showing the comparison of the ruminal microbiota of bacteria and archaea (A), fungi (B), and protozoa (C) between the heat-stressed cows (red dot) and the pair-fed thermal neutral control cows (blue dot). The P-values were based on PERMANOVA analysis

Fig. 2
figure 2

Co-occurrence and mutual exclusion network based on significant correlations (P ≤ 0.1) of relative abundance of major classified genera (each detected in at least 50% of the sample) of ruminal microbiota of pair-fed thermal neutral (PFTN) control cows and heat-stressed (HS) cows (A), and the genus-level nodes exclusively found in the PFTN cows and the HS cows (B). The size of the nodes represents the sum of the weights of the interactions connected. The exclusive nodes were marked as red asterisks within each exclusive microbial network. Blue and red lines indicate co-occurrence and mutual exclusion, respectively

Heat stress affects some of the taxa of the ruminal microbiota

Comparison using Venn diagrams revealed the microbial genera (both taxonomically classified and unclassified genus-level equivalents) shared by or exclusively detected in the PFTN control cows and the HS cows (Fig. 3). In the prokaryotic microbiota, 88 of the 142 detected genera were shared between the two thermal conditions, and so were all the three identified archaeal genera. Thermoneutrality and HS corresponded to 30 and 24 bacterial genera that were exclusively found therefor. Of the 164 detected fungal genera (including taxonomically unclassified genus-level equivalents), 81 were shared, whereas 47 and 36 were exclusively found in the PFTN control cows and the HS cows, respectively. Only seven protozoal genera were detected, and all of them were shared between the two thermal conditions.

Fig. 3
figure 3

Venn diagrams showing the genera of rumen microbes shared between and unique to the heat-stressed (HS) cows and the pair-fed thermal neutral (PFTN) control cows

Among the taxonomically classified genera, 24 genera in 13 phyla of prokaryotes, 29 genera in three phyla of fungi, and seven genera in the same subclass (Trichostomatia) of protozoa were found in at least 50% of all the rumen samples (Additional file 2: Fig. S1), and they were considered as prevalent taxa. Analysis using LEfSe identified differentially abundant microbial phyla and genera between PFTN control and HS cows (Fig. 4). At the phylum level, Firmicutes and Neocallimastigomycota were enriched in PFTN control cows, while Bacteroidetes was enriched in the HS cows. This resulted in a decreased (P < 0.001) Firmicutes-to-Bacteroidetes ratio in the HS cows (Fig. 5). There was also a period (before and after the animals were crossed over to the other thermal condition) effect on the Firmicutes-to-Bacteroidetes ratio (P = 0.002) though there was no thermal condition by period interaction with respect to this ratio (P = 0.20). Two of the cows (cows 5125 and 10,046) also exhibited a significant difference (P < 0.05) in Firmicutes-to-Bacteroidetes ratio under either thermal condition (Fig. 5). Different from the other cows, cow 5125 did not decrease the Firmicutes-to-Bacteroidetes ratio when heat stressed. Among the taxonomically classified genera of bacteria, Ruminococcus and Desulfovibrio were enriched in the PFTN control cows, while Shuttleworthia and Anaeroplasma expanded their relative abundance in the HS cows. Another six unclassified genera (one each assigned to Clostridiales, WCHB1–25, Lachnospiraceae, Ruminococcaceae, Bacteroidetes, and Paraprevotellaceae) were also enriched in the PFTN control cows (data not shown). Although not identified as a biomarker by the LEfSe analysis, Prevotella tended to increase (P < 0.093) its relative abundance in the HS cows, by more than 10% (17.1% vs. 27.9%). Of the prevalent fungal genera, Filobasidium, Mycosphaerella, and Issatchenkia were enriched in the HS cows, while Piromyces was enriched in the PFTN control cows (Fig. 4). Of the protozoa, the genus Isotricha was enriched in the PFTN control cows, and the HS did not increase the relative abundance of any of the protozoal genera.

Fig. 4
figure 4

Differentially abundant rumen microbial taxa between heat-stressed (HS) cows and the pair-fed thermal neutral (PFTN) control cows, which were detected using LEfSe with an LDA effect size ≥ 2. Only classified prevalent taxa (each detected in at least 50% of the samples) were visualized with additional taxonomic lineages of prokaryotes (B) and fungi (C) embedded in a respective cladogram

Fig. 5
figure 5

Firmicutes-to-Bacteroidetes ratio (A) among the four pair-fed thermal neutral (PFTN) control cows or the heat-stressed (HS) cows and (B) between the two groups of cows. * indicates P < 0.05. P-value: treatment: 0.001; period: 0.002; T × P: 0.198

Heat stress affects some predicted functions of prokaryotes

The overall numbers of functional features predicted using the 16S rRNA gene and PICRUSt2 with six different reference databases were not changed by HS (Additional file 1: Table S2). Based on the relative abundances of predicted features (indicated by EC numbers), the overall distributions of the prokaryotic enzyme profiles were not affected by HS (P = 0.117), while the fungal enzyme profiles tended (P = 0.094) to be affected (Fig. 6).

Fig. 6
figure 6

PCA plots (based on Bray-Curtis dissimilarity) comparing the overall functional features of the rumen microbiota of prokaryotes (A) and fungi (B). Functional features were predicted using PICRUSt2 and the ExPASy ENZYME database. P-values were based on PERMANOVA (9,999 permutations)

Comparative analysis using LEfSe identified the biomarker enzymes (indicated by EC numbers) that were differentially abundant (P < 0.05, LDA score > 2) between the HS cows and PFTN control cows, particularly among the major enzymes each with a relative abundance > 0.1% for at least one of the thermal conditions (Table 2 and Table 3). Overall, 18 and 9 EC numbers of prokaryotes were enriched and reduced, respectively by HS. These included three oxidoreductases (EC 1.-.-.-), four transferases (EC 2.-.-.-) including three transferases transferring phosphorus-containing groups (EC 2.7.-.-), seven hydrolases (EC 3.-.-.-) including three hydrolases acting on ester bonds (EC 3.1.-.-) and three hydrolases acting on peptide bonds (peptidases, EC 3.4.-.-), and two isomerases (EC 5.-.-.-). Of the fungal enzymes, six and 14 were enriched and reduced, respectively by HS (Table 2 and Table 3). These included three oxidoreductases (EC 1.-.-.-) including two oxidoreductases acting on the CH-OH group of donors.

Table 2 Predicted prokaryotic and fungal enzymes* significantly (P < 0.05) enriched (as detected using LEfSe) in the ruminal microbiota of heat-stressed (HS) cows relative to pair-fed thermal neutral (PFTN) control cows
Table 3 Predicted prokaryotic and fungal enzymes* significantly (P < 0.05) reduced (as detected using LEfSe) in the ruminal microbiota of heat-stressed (HS) cows relative to pair-fed thermal neutral (PFTN) control cows


Although not affecting ruminal temperature [54], heat stress has been shown to alter the ruminal microbiota in dairy cows [25, 29, 30]. However, all the studies in the literature did not eliminate the negative effect on feed intake [17, 55], failing to separate the direct effect from the indirect effect caused by decreased feed intake. In this study, we eliminated that confounding factor using pair-feeding so that we could investigate how heat stress can affect the ruminal microbiota independent of its effect on feed intake. We used a small number of animals (n = 4 per thermal conditions) because it is prohibitive to use many animals in large environmental chambers, which are costly to purchase and operate. We used a 2 × 2 crossover design to minimize the animal effect. Overall, the heat stress conditions (THI = 87.2 during daytime and THI = 81.8 during nighttime) affected the three kingdoms (prokaryotes, fungi, and protozoa) differently. Except for the genus richness of protozoal microbiota, the heat stress treatment did not affect any of the α-diversity metrics of any of the three kingdoms. This is in general agreement with the results (only bacterial microbiota was analyzed) of another study on Holstein cows, except for Chao1 richness estimates [25] and Hanwoo steers [30]. These results suggest that a decrease in feed intake may not exacerbate the effect of heat stress on the α-diversity of the ruminal microbiota.

The heat stress treatment significantly altered the protozoal microbiota as shown by PCoA and PERMONAVA analyses and tended to decrease Shannon’s index and evenness. Analysis using LEfSe revealed a significant decrease in the relative abundance of Isotricha. The relative abundance of Dasytricha was 3.92-fold lower in HS cows compared to that of PFTN cows. Heat stress is accompanied by low rumen pH due to decreased saliva secretion [3], and indeed, the pH of the rumen was lower in the HS cows than in the PFTN control cows (pH 6.7 vs. 6.3, P = 0.01) [11]. The lower rumen pH in the HS cows than the PFTN control cows might be a major reason explaining the observed effect of heat stress on rumen protozoal microbiota, which are inhibited by low ruminal pH [56]. The decrease of Isotricha and Dasytricha also agrees with the report that these two genera of holotrichs are probably more sensitive to ruminal pH than other protozoal genera [57]. In the present study, we did not determine if the heat stress treatment decreased total rumen protozoal counts. No study in the literature has included protozoa in the evaluation of the effect of heat stress on ruminal microbiota. Given the sensitivity to rumen pH, rumen protozoa probably could lose abundance in heat stressed cows, but quantitative analyses using microscopic counting or qPCR are needed to verify this surmise in future studies. It should be noted that there was a significant period effect with respect to the impact of heat stress on the protozoal Shannon’s index, Simpson’s index, and evenness. It remains to be determined how the residual effect of period (thermal vs. heat stress in a crossover design) can affect the above α-diversity metrics. Such period effect should be taken into consideration in the design of future studies.

Rumen fungi are recognized for their fibrolytic activity [58, 59], and many of the gut fungi have not been classified [60]. Although all the taxonomically known genera of rumen fungi were detected in the PFTN control cows, in the present study, heat stress did not affect the diversity (α or β) of the fungal microbiota. However, some genus-level taxa disappeared while others appeared in the HS cows. The heat stress also considerably enriched the genera Mycosphaerella, Filobasidium, and Issatchenkia, while diminishing Neocallimastigomycota and the genus Piromyces. Except for Piromyces, all those genera were not commonly found in the rumen and had low relative abundance. The animal study [11] showed the digestibility of all the nutrients (NDF, DM, ADF, CP, EE, and OM) were significantly higher (P < 0.05, without significant period effect) in the HS cows than in the PFTN control cows even though heat stress impaired lactation performance [decreasing milk yield by 17.0% (25.9 vs. 21.5 kg/d), milk protein by 4.1% (2.68% vs. 2.57% of milk), milk protein yield by 19% (705 vs. 571 g/d), 4% fat-corrected milk by 23% (28.0 vs. 21.6 kg/d), and fat yield by 19% (1096 vs. 884 g/d)]. The authors attributed the improvements of feed digestion to decreased rumen motility and thus passage rate [11], both of which are common in HS cows [61]. Heat stress also tended to alter the predicted fungal features and decreased many of the fungal enzymes, especially hydrolases including phosphotransferase and esterases. This discrepancy is consistent with the overall minor role of rumen fungi in feed digestion. To the best of our knowledge, no study in the literature has examined how heat stress could affect the rumen fungal microbiota. Future research is needed to quantitatively determine to what extent heat stress affects the rumen fungal populations and their function.

Prokaryotes, which account for the majority of the ruminal microbiota, were analyzed in all studies on ruminal microbiota, but ruminal fungi and protozoa were rarely analyzed. In the present study, heat stress altered the overall prokaryotic microbiota (significantly if based on unweighted UniFrac analysis, data not shown). This finding is consistent with the report of Kim et al. [26] but does not corroborate the report by Zhao et al. [25]. However, the three studies differ in the actual THI imposed on the cows, making comparison and interpretation of the results difficult. The heat stress conditions of the present study significantly increased the relative abundance of Bacteroidetes (primarily Gram-negative bacteria), while decreasing that of Firmicutes (primarily Gram-positive bacteria). Such a phylum-level shift was reported in the ruminal microbiota in heat-stressed goats [29] and the fecal microbiota of heat-stressed lactating dairy cows [62]. Future research is needed to verify if this Firmicutes to Bacteroidetes shift has anything to do with their Gram staining features and if the Firmicutes-to-Bacteroidetes ratio of the rumen or fecal microbiota can serve as a microbial biomarker of heat stress. In the recent studies, a shift of other bacterial phyla was reported, including a decrease in Fibrobacteres but an increase in Fusobacteria, Tenericutes, and Cyanobacteria in heat-stressed lactating Holstein cows [26] and a decrease in Proteobacteria and Chloroflexi in heat-stressed Hanwoo steers [30]. In the present study, these phyla were not affected by heat stress. The discrepancy might be attributable to the decreased feed intake by the HS cows in the two cited studies.

Many bacterial genera (including some that have no cultured species) were also detected exclusively under each thermal condition. Similar to the case of fungi, nearly all these genera had a very low relative abundance. Given the very high Good’s coverage (≥ 99.9%), those exclusively detected genera were probably very minor ones selected by each thermal condition. Only three taxonomically known genera (unclassified genus-level taxa were not discussed or compared with other studies because it is not possible) were enriched (Prevotella, Shuttleworthia, and Anaeroplasma) or diminished (Ruminococcus and Desulfovibrio) by heat stress. The observed decrease in Ruminococcus in the HS cows contradicts the report for heat-stressed goats [29] but concurs with the reports for heat-stressed dairy cows [25, 26] or Hanwoo steers [30]. On the other hand, none of the genera (except for Prevotella) detected in the HS cows in other recent studies was detected in the present study (Additional file 1: Table S3). Among other factors, the elimination of the confounding effect on feed intake might be attributable to these discrepancies. The “heat-stress resistant” microbes that were enriched or exclusively found in the rumen of heat-stressed cows may be useful as potential probiotics for cows under heat stress.

The decrease in Ruminococcus, which is a highly fibrolytic genus [63], and Desulfovibrio, which is a sulfate-reducing VFA-utilizing genus [64], in the HS cows, seems contrary to the increased fiber digestibility (by 21.4%, dry matter) and VFA concentration (by 82.9%, total VFA) in those HS cows. The increase in Shuttleworthia, which is saccharolytic producing acetic, butyric, and lactic acids as major fermentation products [65], and Anaeroplasma¸ which ferments sugars to primarily acetic, formic, propionic, lactic, and succinic acids [66], is consistent with the increased rumen VFA concentration in the HS cows, but their absolute abundance needs to be quantified to determine if and to what extent they can contribute to the increased VFA in the HS cows. As the largest genus, Prevotella increased its relative abundance by > 10% (27.9% vs. 17.1%) in the HS cows, and the large increase of this versatile genus might explain some of the observed changes in nutrient digestibility and fermentation characteristics [11]. It should be noted that many taxonomically unknown genus-level taxa were exclusively found in either the HS cows or the PFTN control cows, suggesting that they may adapt or be susceptible to the rumen environment created by heat stress. Their role and contribution to the rumen function in HS cows remain to be determined.

Using two different network analyses (FastSpar and CoDiNA), we identified some co-occurrence and mutual-exclusion relationships between some taxonomically known genera exclusively detected for each thermal condition. Interestingly, more co-occurrence or mutual exclusion relationships were found in the ruminal microbiota of the PFTN control cows than the HS cows (17 vs. 10). The two thermal conditions also corresponded to co-occurrence and mutual exclusion of very different genera. In the literature, only one study examined co-occurrence and mutual exclusion of rumen microbial taxa while determining the effect of heat stress on the ruminal microbiota in goats [29]. Mutual exclusion between Ruminococcus and Prevotella was reported, but it was not stated if it was in the rumen of thermal control goats or heat-stressed goats. Because co-occurrence and mutual exclusion were based on correlation of relative abundance between genera, they may not necessarily be attributable to mutualism or antagonism, respectively. However, the co-occurrence and mutual exclusion relationships detected exclusively for each thermal condition suggest that heat stress can alter the response of those genera.

All the recent studies evaluated the impact of heat stress on the taxonomic diversity and composition of the ruminal microbiota [25, 26, 29, 30], but only two studies [26, 29] evaluated the effect on the functional features of the ruminal microbiota. In the present study, we predicted the functional features based on the 16S rRNA gene sequence data using PICRUSt2. The overall functional features of prokaryotes did not differ between the HS cows and the PFTN control cows. The discrepancy between the effect of heat stress on the rumen prokaryotic microbiota (PERMANOVA P = 0.068) and on the predicted prokaryotic functional profiles (PERMANOVA P = 0.117) is consistent with the functional redundancy of the ruminal microbiota [67]. The heat stress did affect some of the predicted enzymes, particularly striking the enrichment of seven hydrolases including those acting on ester bonds and peptide bonds. This corroborates the increased nutrient digestibility reported in the HS cows [11]. The increased retention time corresponding to the slowed passage rate might have increased some hydrolytic bacteria. However, the minimal impact of heat stress on the rumen function does not mirror the impaired lactation performance (by 17.0%) in the HS cows. The prediction of functions using phylogenetic markers certainly has limitations. In future studies, direct and quantitative functional analyses using other technologies, such as metatranscriptomics, metabolomics, and guild-specific tests, can better help determine how heat stress affects rumen functions. Nonetheless, these results suggest that heat stress may decrease lactation performance in dairy cows by largely impairing the physiology and metabolism, not the rumen function. Indeed, heat stress was found to decrease milk protein synthesis in the mammary glands [11] and other important biological functions such as metabolism, immunity, inflammation, changes in behavior, and antioxidant capacity [8, 13, 14, 68]. A recent study using Holstein heifers showed that fermented herbal tea residues improved some physiological indices related to heat stress, including respiratory frequency, rectal temperature, concentrations of immunoglobulins, antioxidant capacity, and serum concentrations of heat stress-related parameters [69]. In finishing male pigs, Saccharomyces cerevisiae boulardii CNCM I-1079 supplementation was shown to help cope with heat stress by modulating their gut microbiome and feeding behavior [70]. Future studies are warranted to evaluate if these approaches can be applied to lactating dairy cows.

As global warming continues to worsen, heat stress will affect animals and humans, particularly those living in tropical regions, to a greater extent than what has been reported. In the present study, we focused on the ruminal microbiome because of its critical role in feed digestion and supplying much of the nutrients to cows. Other studies have shown that heat stress could also affect the fecal microbiome [71,72,73], a proxy of the microbiome of large intestines. Although not reported, heat stress can also affect the gut microbiome of animals and humans. Pair-feeding should be used in future studies to examine the effect of heat stress on the gut microbiome thereof.


Heat stress affected the ruminal microbiota structure of the three kingdoms (prokaryotes, fungi, and protozoa) to a different extent in the absence of the confounding factor of feed intake. Heat stress selected some taxa of bacteria, fungi, and protozoa while decreasing others. Heat stress also altered the co-occurrence and mutual exclusion between some genera. Firmicutes-to-Bacteroidetes ratio decreased in heat-stressed cows, and several genera decreased (Ruminococcus, Desulfovibrio, Piromyces, Isotricha) or increased (Anaeroplasma, Shuttleworthia, Filobasidium). The heat stress conditions had little impact on the predicted functions of the rumen prokaryotic or fungal microbiota. The effect of heat stress on the ruminal microbiota observed in the present study should be primarily attributed to its adverse effect on the host metabolism, physiology, and behavior. Some of the “heat-stress resistant” bacteria and fungi may be useful as potential probiotics for cows under heat stress.

Availability of data and materials

The 16S rRNA gene amplicon sequences, ITS1 amplicon sequences, and 18S rRNA gene sequences generated for this study can be found in the NCBI Sequence Read Archive (SRA) under BioProject PRJNA562317.



Amplicon sequence variants


Enzyme Commission


Greenhouse gases


Heat stress; Heat-stressed


Integrated Microbial Genomes


Linear discriminant analysis effect size


Pair-fed thermoneutrality; Pair-fed thermal neutral


Principal components analysis


Principal coordinates analysis


Permutational multivariate analysis of variance


Temperature-humidity index


  1. Becker C, Stone A. Graduate student literature review: heat abatement strategies used to reduce negative effects of heat stress in dairy cows. J Dairy Sci. 2020;103(10):9667–75.

    CAS  Article  PubMed  Google Scholar 

  2. Bouraoui R, Lahmar M, Majdoub A, Belyea R. The relationship of temperature-humidity index with milk production of dairy cows in a mediterranean climate. Anim Res. 2002;51(6):479–91.

    Article  Google Scholar 

  3. Das R, Sailo L, Verma N, Bharti P, Saikia J. Impact of heat stress on health and performance of dairy animals: a review. Vet World. 2016;9(3):260–8.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  4. Sammad A, Wang YJ, Umer S, Lirong H, Khan I, Khan A, et al. Nutritional physiology and biochemistry of dairy cattle under the influence of heat stress: consequences and opportunities. Animals. 2020;10(5):793.

    Article  PubMed Central  Google Scholar 

  5. West JW. Effects of heat-stress on production in dairy cattle. J Dairy Sci. 2003;86:2131–44.

    CAS  Article  PubMed  Google Scholar 

  6. Ranjitkar S, Bu D, Van Wijk M, Ma Y, Ma L, Zhao L, et al. Will heat stress take its toll on milk production in China? Clim Chang. 2020;161(4):637–52.

    Article  Google Scholar 

  7. St-Pierre N, Cobanov B, Schnitkey G. Economic losses from heat stress by us livestock industries. J Dairy Sci. 2003;86:E52–77.

  8. Becker C, Collier R, Stone A. Invited review: physiological and behavioral effects of heat stress in dairy cows. J Dairy Sci. 2020;103(8):6751–70.

    CAS  Article  PubMed  Google Scholar 

  9. Herbut P, Angrecka S, Walczak J. Environmental parameters to assessing of heat stress in dairy cattle—a review. Int J Biometeorol. 2018;62(12):2089–97.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Faylon M, Baumgard L, Rhoads R, Spurlock D. Effects of acute heat stress on lipid metabolism of bovine primary adipocytes. J Dairy Sci. 2015;98(12):8732–40.

    CAS  Article  PubMed  Google Scholar 

  11. Gao S, Guo J, Quan S, Nan X, Fernandez MS, Baumgard L, et al. The effects of heat stress on protein metabolism in lactating Holstein cows. J Dairy Sci. 2017;100(6):5040–9.

    CAS  Article  PubMed  Google Scholar 

  12. Gao S, Ma L, Zhou Z, Zhou Z, Baumgard LH, Jiang D, et al. Heat stress negatively affects the transcriptome related to overall metabolism and milk protein synthesis in mammary tissue of midlactating dairy cows. Physiol Genomics. 2019;51(8):400–9.

    CAS  Article  PubMed  Google Scholar 

  13. Abbas Z, Sammad A, Hu L, Fang H, Xu Q, Wang Y. Glucose metabolism and dynamics of facilitative glucose transporters (gluts) under the influence of heat stress in dairy cattle. Metabolites. 2020;10(8):312.

    CAS  Article  PubMed Central  Google Scholar 

  14. Dahl GE, Tao S, Laporta J. Heat stress impacts immune status in cows across the life cycle. Front Vet Sci. 2020;7:116.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Collier R, Renquist B, Xiao Y. A 100-year review: stress physiology including heat stress. J Dairy Sci. 2017;100(12):10367–80.

    CAS  Article  PubMed  Google Scholar 

  16. Garner J, Douglas M, Williams S, Wales W, Marett L, DiGiacomo K, et al. Responses of dairy cows to short-term heat stress in controlled-climate chambers. Anim Prod Sci. 2017;57(7):1233–41.

    Article  Google Scholar 

  17. Tao S, Rivas RMO, Marins TN, Chen Y-C, Gao J, Bernard JK. Impact of heat stress on lactational performance of dairy cows. Theriogenology. 2020;150:437–44.

    Article  PubMed  Google Scholar 

  18. Kim M, Morrison M, Yu Z. Status of the phylogenetic diversity census of ruminal microbiomes. FEMS Microbiol Ecol. 2011;76(1):49–63.

    CAS  Article  PubMed  Google Scholar 

  19. Creevey CJ, Kelly WJ, Henderson G, Leahy SC. Determining the culturability of the rumen bacterial microbiome. Microb Biotechnol. 2014;7(5):467–79.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Henderson G, Cox F, Ganesh S, Jonker A, Young W, Janssen PH. Rumen microbial community composition varies with diet and host, but a core microbiome is found across a wide geographical range. Sci Rep. 2015;5(1):1–15.

    CAS  Article  Google Scholar 

  21. Firkins JL, Yu Z. Characterisation and quantification of the microbial populations in the rumen. In: Sejrsen K, Hvelplund T, Nielsen MO, editors. Ruminant physiology, digestion, metabolism and impact of nutrition on gene expression, immunology and stress. Wageningen: Academic Publishers; 2006. p. 19–54.

  22. Park T, Ma L, Ma Y, Zhou X, Bu D, Yu Z. Dietary energy sources and levels shift the multi-kingdom microbiota and functions in the rumen of lactating dairy cows. J Anim Sci Biotechnol. 2020;11:66.

  23. Newbold C, Ramos-Morales E. Ruminal microbiome and microbial metabolome: effects of diet and ruminant host. Animal. 2020;14(S1):s78–86.

    CAS  Article  PubMed  Google Scholar 

  24. Menajovsky S, Walpole C, DeVries T, Schwartzkopf-Genswein K, Walpole M, Penner G. The effect of the forage-to-concentrate ratio of the partial mixed ration and the quantity of concentrate in an automatic milking system for lactating Holstein cows. J Dairy Sci. 2018;101(11):9941–53.

    CAS  Article  PubMed  Google Scholar 

  25. 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. 2019;9(11):925.

    Article  PubMed Central  Google Scholar 

  26. Kim D-H, Kim M-H, Kim S-B, Son J-K, Lee J-H, Joo S-S, et al. Differential dynamics of the ruminal microbiome of Jersey cows in a heat stress environment. Animals. 2020;10(7):1127.

    Article  PubMed Central  Google Scholar 

  27. Tajima K, Nonaka I, Higuchi K, Takusari N, Kurihara M, Takenaka A, et al. Influence of high temperature and humidity on rumen bacterial diversity in Holstein heifers. Anaerobe. 2007;13(2):57–64.

    Article  PubMed  Google Scholar 

  28. Uyeno Y, Sekiguchi Y, Tajima K, Takenaka A, Kurihara M, Kamagata Y. An rrna-based analysis for evaluating the effect of heat stress on the rumen microbial composition of Holstein heifers. Anaerobe. 2010;16(1):27–33.

    CAS  Article  PubMed  Google Scholar 

  29. 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 Biotehcnol. 2019;103(7):3193–203.

    CAS  Article  Google Scholar 

  30. Baek YC, Choi H, Jeong J-Y, 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(2):208–17.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  31. Rodríguez C, González J, Alvir M, Redondo R, Cajarville C. Effects of feed intake on composition of sheep rumen contents and their microbial population size. Brit J Nutr. 2003;89(1):97–103.

    Article  PubMed  Google Scholar 

  32. Kala A, Kamra D, Kumar A, Agarwal N, Chaudhary L, Joshi C. Impact of levels of total digestible nutrients on microbiome, enzyme profile and degradation of feeds in buffalo rumen. PLoS One. 2017;12(2):e0172051.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  33. McCabe MS, Cormican P, Keogh K, O’Connor A, O’Hara E, Palladino RA, et al. Illumina miseq phylogenetic amplicon sequencing shows a large reduction of an uncharacterised Succinivibrionaceae and an increase of the Methanobrevibacter gottschalkii clade in feed restricted cattle. PLoS One. 2015;10(7):e0133234.

  34. Thom EC. The discomfort index. Weatherwise. 1959;12(2):57–61.

    Article  Google Scholar 

  35. Alhidary I, Shini S, Al Jassim R, Gaughan J. Physiological responses of australian merino wethers exposed to high heat load. J Anim Sci. 2012;90(1):212–20.

    CAS  Article  PubMed  Google Scholar 

  36. NRC. Nutrient requirements of dairy cattle. 7th rev. ed. Washington, DC: The national academies press; 2001.

    Google Scholar 

  37. Paz HA, Anderson CL, Muller MJ, Kononoff PJ, Fernando SC. Rumen bacterial community composition in Holstein and Jersey cows is different under same dietary condition and is not affected by sampling method. Front Microbiol. 2016;7:1206.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Yu Z, Morrison M. Improved extraction of pcr-quality community DNA from digesta and fecal samples. Biotechniques. 2004;36(5):808–12.

    CAS  Article  PubMed  Google Scholar 

  39. Fadrosh DW, Ma B, Gajer P, Sengamalay N, Ott S, Brotman RM, et al. An improved dual-indexing approach for multiplexed 16s rrna gene sequencing on the illumina miseq platform. Microbiome. 2014;2(1):1–7.

    Article  Google Scholar 

  40. 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(8):852–7.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  41. Martin M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet J. 2011;17(1):10–2.

    Article  Google Scholar 

  42. Callahan BJ, McMurdie PJ, Rosen MJ, Han AW, Johnson AJA, Holmes SP. Dada2: high-resolution sample inference from illumina amplicon data. Nat Methods. 2016;13(7):581–3.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  43. Rivers AR, Weber KC, Gardner TG, Liu S, Armstrong SD. Itsxpress: Software to rapidly trim internally transcribed spacer sequences with quality scores for marker gene analysis. F1000Res. 2018;7:1418.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  44. Bokulich NA, Kaehler BD, Rideout JR, Dillon M, Bolyen E, Knight R, et al. Optimizing taxonomic classification of marker-gene amplicon sequences with qiime 2’s q2-feature-classifier plugin. Microbiome. 2018;6(1):1–17.

    Article  Google Scholar 

  45. DeSantis TZ, Hugenholtz P, Larsen N, Rojas M, Brodie EL, Keller K, et al. Greengenes, a chimera-checked 16s rrna gene database and workbench compatible with arb. Appl Environ Microbiol. 2006;72(7):5069–72.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  46. UNITE Community. Unite qiime release for fungi: UNITE Community; 2019.

    Google Scholar 

  47. Lozupone C, Knight R. Unifrac: a new phylogenetic method for comparing microbial communities. Appl Environ Microbiol. 2005;71(12):8228–35.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  48. Watts SC, Ritchie SC, Inouye M, Holt KE. Fastspar: rapid and scalable correlation estimation for compositional data. Bioinformatics. 2019;35(6):1064–6.

    CAS  Article  PubMed  Google Scholar 

  49. Gysi DM, Fragoso TM, Buskamp V, Almaas E, Nowick K. Comparing multiple networks using the co-expression differential network analysis (CoDiNA). arXiv preprint arXiv. 2018:180200828.

  50. Douglas GM, Maffei VJ, Zaneveld JR, Yurgel SN, Brown JR, Taylor CM, et al. Picrust2 for prediction of metagenome functions. Nat Biotechnol. 2020;38(6):685–8.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  51. Beals EW. Bray-curtis ordination: An effective strategy for analysis of multivariate ecological data. Adv Ecol Res. 1984;14:1–55.

    Article  Google Scholar 

  52. Tang Y, Horikoshi M, Li W. Ggfortify: unified interface to visualize statistical results of popular r packages. R J. 2016;8(2):474.

    Article  Google Scholar 

  53. Segata N, Waldron L, Ballarini A, Narasimhan V, Jousson O, Huttenhower C. Metagenomic microbial community profiling using unique clade-specific marker genes. Nat Methods. 2012;9(8):811–4.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  54. Lees A, Lees J, Lisle A, Sullivan M, Gaughan J. Effect of heat stress on rumen temperature of three breeds of cattle. Int J Biometeorol. 2018;62(2):207–15.

    CAS  Article  PubMed  Google Scholar 

  55. Liu J, Li L, Chen X, Lu Y, Wang D. Effects of heat stress on body temperature, milk production, and reproduction in dairy cows: a novel idea for monitoring and evaluation of heat stress—a review. Asian Austral J Anim Sci. 2019;32(9):1332–9.

    Article  Google Scholar 

  56. Purser D, Moir R. Ruminal flora studies in the sheep. Ix. The effect of ph on the ciliate population of the rumen in vivo. Austral J Agric Res. 1959;10(4):555–64.

  57. Santra A, Chaturvedi O, Tripathi M, Kumar R, Karim S. Effect of dietary sodium bicarbonate supplementation on fermentation characteristics and ciliate protozoal population in rumen of lambs. Small Ruminant Res. 2003;47(3):203–12.

    Article  Google Scholar 

  58. Paul SS, Deb SM, Punia BS, Singh D, Kumar R. Fibrolytic potential of anaerobic fungi (Piromyces sp.) isolated from wild cattle and blue bulls in pure culture and effect of their addition on in vitro fermentation of wheat straw and methane emission by rumen fluid of buffaloes. J Sci Food Agric. 2010;90:1218–26.

  59. Wang P, Qi M, Barboza P, Leigh MB, Ungerfeld E, Selinger LB, et al. Isolation of high-quality total rna from rumen anaerobic bacteria and fungi, and subsequent detection of glycoside hydrolases. Can J Microbiol. 2011;57(7):590–8.

    CAS  Article  PubMed  Google Scholar 

  60. Paul SS, Bu D, Xu J, Hyde KD, Yu Z. A phylogenetic census of global diversity of gut anaerobic fungi and a new taxonomic framework. Fungal Div. 2018;89(1):253–66.

    Article  Google Scholar 

  61. Bernabucci U, Lacetera N, Danieli PP, Bani P, Nardone A, Ronchi B. Influence of different periods of exposure to hot environment on rumen function and diet digestibility in sheep. Int J Biometeorol. 2009;53(5):387–95.

    CAS  Article  PubMed  Google Scholar 

  62. Chen S, Wang J, Peng D, Li G, Chen J, Gu X. Exposure to heat-stress environment affects the physiology, circulation levels of cytokines, and microbiome in dairy cows. Sci Rep. 2018;8:14606.

  63. Ezaki T. Ruminococcus. In: Trujillo ME, Dedysh S, Devos P, Hedlund B, Kämpfer P, Rainey FA, et al., editors. Bergey's Manual of Systematics of Archaea and Bacteria; 2015. p. 1–5.

  64. Galushko A, Kuever J. Desulfovibrio. In: Trujillo ME, Dedysh S, DeVos P, Hedlund B, Kämpfer P, Rainey FA, et al., editors. Bergey's Manual of Systematics of Archaea and Bacteria; 2015. p. 1–10.

  65. Wade WG, Downes J. Shuttleworthia. In: Trujillo ME, Dedysh S, DeVos P, Hedlund B, Kämpfer P, Rainey FA, et al., editors. Bergey's Manual of Systematics of Archaea and Bacteria; 2015. p. 1–3.

  66. Brown DR, Bradbury JM, Johnson K-E. Anaeroplasma. In: Trujillo ME, Dedysh S, DeVos P, Hedlund B, Kämpfer P, Rainey FA, et al., editors. Bergey's Manual of Systematics of Archaea and Bacteria; 2015. p. 1–5.

  67. Weimer PJ. Redundancy, resilience, and host specificity of the ruminal microbiota: implications for engineering improved ruminal fermentations. Front Microbiol. 2015;6:296.

    Article  PubMed  PubMed Central  Google Scholar 

  68. Bagath M, Krishnan G, Devaraj C, Rashamol V, Pragna P, Lees A, et al. The impact of heat stress on the immune system in dairy cattle: a review. Res Vet Sci. 2019;126:94–102.

    CAS  Article  PubMed  Google Scholar 

  69. Xie Y, Chen Z, Wang D, Chen G, Sun X, He Q, et al. Effects of fermented herbal tea residues on the intestinal microbiota characteristics of Holstein heifers under heat stress. Front Microbiol. 2020;11:1014.

    Article  PubMed  PubMed Central  Google Scholar 

  70. Labussière E, Achard C, Dubois S, Combes S, Castex M, Renaudeau D. Saccharomyces cerevisiae boulardii cncm i-1079 supplementation in finishing male pigs helps to cope with heat stress through feeding behaviour and gut microbiota modulation. Brit J Nutr. 2022;127(3):353–68.

  71. He J, Guo H, Zheng W, Xue Y, Zhao R, Yao W. Heat stress affects fecal microbial and metabolic alterations of primiparous sows during late gestation. J Anim Sci Biotechnol. 2019;10:84.

  72. Zhu L, Liao R, Wu N, Zhu G, Yang C. Heat stress mediates changes in fecal microbiome and functional pathways of laying hens. Appl Microbiol Biotech. 2019;103(1):461–72.

    CAS  Article  Google Scholar 

  73. Li H, Li R, Chen H, Gao J, Wang Y, Zhang Y, et al. Effect of different seasons (spring vs summer) on the microbiota diversity in the feces of dairy cows. Int J Biometeorol. 2020;64(3):345–54.

Download references


Not applicable.


This research was partially funded by the National Natural Science Foundation of China (award number: 31872383), the Scientific Research Project for Major Achievements of The Agricultural Science and Technology Innovation Program (award number: ASTIP-IAS07–1), and the Beijing Dairy Industry Innovation Team (Award number: BAIC06–2021).

Author information

Authors and Affiliations



DB and ZY conceived and designed the study. LM and SG conducted the animal experiment and collected the samples. TP conducted the metataxonomic analysis. TP and ZY wrote the paper. All the authors read and edited the paper. The author(s) read and approved the final manuscript.

Corresponding authors

Correspondence to Dengpan Bu or Zhongtang Yu.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

All the authors read and agree to the content of this paper and its publication.

Competing interests

The authors declare no competing interest.

Supplementary Information

Additional file 1: Table S1.

DADA2 denoising statistics of the ruminal microbiota of the heat-stressed (HS) cows and the pair-fed thermal neutral (PFTN) control cows. Table S2. Counts of PICTRUSt2-predicted functional features in the ruminal microbiota of heat-stressed (HS) cows and pair-fed thermal neutral (PFTN) control cows. Table S3. Comparison of the taxonomically known genera of bacteria affected by heat stress.

Additional file 2: Fig. S1.

Relative abundance of prevalent (found in at least 50% of all the samples) classified prokaryotic phyla (A), prokaryotic genera (B), fungal phyla (C), fungal genera (D), and protozoal genera (E).

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Park, T., Ma, L., Gao, S. et al. 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 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Functional profiles
  • Heat stress
  • Microbiome
  • Multi-kingdom
  • Network analysis
  • Ruminal microbiota