- Research article
- Open Access
Changes in ruminal and reticular pH and bacterial communities in Holstein cattle fed a high-grain diet
BMC Veterinary Research volume 14, Article number: 310 (2018)
Subacute ruminal acidosis (SARA) is characterized by a ruminal pH depression, and microbiota can also be affected by a higher acidity and/or dietary changes. Previous studies have revealed similar patterns in pH reduction in the rumen and reticulum, whereas changes in reticular pH and bacterial community following a high-grain diet are not fully understood. The aim of this study was to investigate the changes in reticular pH and bacterial community structure following a high-grain diet simultaneously with those in the rumen.
SARA was diagnosed when ruminal and reticular pH remained under 5.6 for 350 ± 14 and 312 ± 24 min/d, respectively, on the last day of the CON period. During the CON period, lower proportion of acetic acid and higher proportion of butyric acid were observed compared with the HAY period. The proportions of acetic acid and propionic acid were lower and higher, respectively, in the rumen compared with the reticulum. From 454 pyrosequencing analysis, the relative abundance of several genera differed significantly between the two periods and the two locations. During the HAY period, higher relative abundances of Prevotella, Eubacterium, Oscillibacter, and Succiniclasticum and lower relative abundances of Ruminococcus, Clostridium, and Olsenella were identified compared with the CON period. Furthermore, the relative abundance of Eubacterium was lower in the rumen compared with the reticulum. Bacterial diversity indices were significantly different between the HAY and CON periods, being higher in the HAY period. The quantitative real-time PCR showed that the copy numbers of several cellulolytic bacteria (Fibrobacter succinogenes and Ruminococcus albus) were higher during the HAY period.
A high-grain diet showed similar impacts on the pH, fermentation, and bacterial community structure in the rumen and reticulum. During the CON period, ruminal and reticular pH decreased following the high-grain challenge, and lower bacterial diversity and changes in the bacterial composition, similarity, and bacterial copy numbers were observed due to a higher acidity and dietary changes compared with the HAY period. These changes may influence the fermentative ability of the rumen and reticulum.
Subacute ruminal acidosis (SARA) is a common health and production problem in dairy herds, and is defined as a condition characterized by ruminal pH < 5.6 for more than 3 h per day . Monitoring of forestomach pH has been suggested as a potentially valuable tool for diagnosing SARA [2, 3]. Sato et al.  measured ruminal and reticular pH simultaneously using a radio transmission pH measurement system and observed similar patterns in pH reduction in the rumen and reticulum after feeding. More recently, Falk et al.  compared the ruminal and reticular pH of cows in early stages of lactation. They found that reticular pH was higher than ruminal pH in averaged pH profiles at all weeks of lactation. Collectively, these two studies demonstrated that reticular pH is closely correlated with ruminal pH.
Pyrosequencing technology is an alternative to traditional gel-based sequencing techniques for de novo DNA sequencing, and has been used to analyze ruminant gastrointestinal microbiota, especially in the rumen [4,5,6]. A decrease in ruminal pH following a high-grain diet is quite common; microbiota and epithelial bacterial communities in the rumen can also change due to higher acidity and/or dietary changes [5, 7]. Sato  reported that a high-grain diet caused not only decreased ruminal pH but also a reduction in the relative abundance of the ruminal phyla Bacteroidetes and Proteobacteria in cattle with repeatedly induced SARA. Furthermore, the addition of starch source reduces bacterial diversity indices [8, 9], and the relative abundances of ruminal phyla are affected by starch and oil additions .
Although rumen microbiota are influenced by rumen conditions such as pH, feeding materials, and nutritional compositions [5, 6], changes in reticular pH and bacterial communities following a high-grain diet are not fully understood. Therefore, the objective of this study was to identify the continuous changes in reticular pH and its effect on bacterial community structure following a high-grain diet simultaneously with those in the rumen.
Ruminal and reticular pH and volatile fatty acids (VFAs)
The 24-h mean ruminal and reticular pH decreased following the high-grain challenge (Fig. 1), and reticular pH was higher than ruminal pH during the HAY and CON periods, except on day 12. No significant differences were observed in 1-h mean pH and duration of pH < 5.6 in the rumen and reticulum. A significant difference (P < 0.05) between the ruminal and reticular pH was identified only on day 1, and ruminal pH was significantly (P < 0.05) decreased on day 9 compared to that on day 7 (Fig. 1). Linear regression analysis revealed a significant relationship between the 24-h mean pH and days in both rumen and reticulum (R2 = 0.79 and R2 = 0.81, respectively; P < 0.05). Diurnal changes in the 1-h mean ruminal and reticular pH observed on the last days of the HAY (day 7) and CON (day 14) periods are shown in Fig. 2. In day 14, linear regression analysis revealed that the ruminal and reticular pH decreased up to 4 h after the morning feeding (R2 = 0.92 and R2 = 0.89, respectively; P < 0.05 and P = 0.057, respectively), and increased from 4 h after the evening feeding (R2 = 0.99 and R2 = 0.97, respectively; P < 0.05). The duration of time where ruminal and reticular pH < 5.6 were 350 ± 14 and 312 ± 24 min/d, respectively, on the last day of the CON period (Fig. 3).
No significant differences were observed in total VFA and NH3-N concentrations in the rumen and reticulum (Table 1). However, a lower proportion of acetic acid and higher proportions of butyric acid and other VFAs (95% confidence interval: 2.14 to 3.46), excluding acetic acid, propionic acid, and butyric acid from the total VFAs, were observed during the CON period compared with the HAY period. The proportions of ruminal acetic acid and propionic acid were significantly lower and higher compared with the reticulum values, respectively, in both periods (P < 0.05). The proportions of ruminal butyric acid and other VFAs were significantly higher in the CON than in the HAY period (P < 0.05). The ratios of ruminal acetic acid to propionic acid in both periods were significantly lower than those in the reticulum (P < 0.05). The concentration of ruminal lactic acid was significantly higher in the CON than in the HAY period (P < 0.05).
Of the major phyla, Firmicutes, Bacteroidetes, and Actinobacteria were the most abundant, accounting for 74.3% of the total ruminal sequences and 74.0% of the total reticular sequences (Fig. 4). The remaining phyla had low relative abundances of < 1%. A total of 337 bacterial genera were identified; the relative abundances of 323 genera comprised < 1% of the total sequences. Of the major genera, Prevotella was the most abundant in both the rumen and reticulum, followed by Ruminococcus and Clostridium. The relative abundance of Prevotella (18.4%) in the rumen was lower than that in the reticulum (21.2%), whereas the relative abundances of Ruminococcus (10.5%) and Clostridium (7.7%) in the rumen were higher than those of the reticulum (7.7% and 6.8%, respectively) regardless of diet.
At the phylum level, the relative abundances of Firmicutes and Bacteroidetes differed significantly between the HAY and CON periods (P < 0.05; Table 2). At the genus level, the relative abundances of Prevotella, Eubacterium (95% confidence interval: 2.02 to 3.18), Oscillibacter, and Succiniclasticum were lower during the CON period, whereas the relative abundances of Ruminococcus, Clostridium, and Olsenella (95% confidence interval: 0.78 to 1.80) were higher during this period; these genera differed significantly between the HAY and CON periods (P < 0.05; Table 2). The relative abundance of Barnesiella exhibited a trend for difference (P = 0.077) between the two periods, being higher during the CON period compared with the HAY period. Considering the relative abundance of Eubacterium spp., a lower relative abundance (P < 0.05) was detected for rumen samples than reticulum samples. The period × location interaction was significant (P < 0.05) for the relative abundance of Olsenella. The relative abundances of Olsenella in the rumen and reticulum during the CON period were significantly (P < 0.05) lower compared with those during the HAY period. During the CON period, the relative abundance of Olsenella was significantly (P < 0.05) higher in the rumen compared with the reticulum. However, there was no statistical significance in the relative abundances of Butyrivibrio (95% confidence interval: 3.20 to 12.0), unclassified Clostridiales (95% confidence interval: 0.73 to 2.93), Acetivibrio (95% confidence interval: 0.63 to 2.27), unclassified Lachnospiraceae (95% confidence interval: 0.33 to 1.95), unclassified Clostridiaceae (95% confidence interval: 0.81 to 1.45), and Bacteroides (95% confidence interval: 0.43 to 1.55) between the two periods and two locations.
Bacterial diversity analysis
Rarefaction curves calculated at a 97% similarity level indicated that the reticulum had higher bacterial diversity than did the rumen (Additional file 1: Figure S1), and the HAY period had higher bacterial diversity than did the CON period, regardless of the location (Table 3). The PCoA results indicated that the HAY period was separate from the CON period in the plot regardless of location, whereas the individual plots for cattle showed close similarity between the rumen and reticulum (PC1 + PC2 = 27.6%; Fig. 5). Moreover, the operational taxonomic units (OTUs), abundance-based coverage estimator (ACE), Chao1, Shannon (95% confidence interval: 4.70 to 5.52), Simpson indices were higher in the HAY period, while Simpson index was higher in the CON period. The effect of location and the interaction effect of period × location were not significant for indices.
Copy numbers of bacterial 16S genes
The copy numbers of total methanogens, Fibrobacter succinogenes, and Selenomonas ruminantium differed significantly (P < 0.05) between the HAY and CON periods, being higher during the HAY period (Table 4). The copy number of Ruminococcus albus was tended (P = 0.068) to be higher during the HAY period compared with the CON period. No statistical difference was identified in the copy numbers of Megasphaera elsdenii (95% confidence interval: 2.93 to 4.87), Streptococcus bovis (95% confidence interval: 5.65 to 7.47), and Ruminococcus flavefaciens (95% confidence interval: 5.92 to 7.042). The copy number of total methanogens was significantly different between the rumen and reticulum (P < 0.05), being higher in the rumen. The period × location interaction was significant for Ruminococcus flavefaciens (P < 0.05). The copy number of Ruminococcus flavefaciens in the rumen during the HAY period was significantly (P < 0.05) higher compared with that during the CON period.
This study aimed to identify the continuous changes in reticular pH following a high-grain diet and to investigate its effect on the reticular bacterial composition, diversity, and similarity simultaneously with those in the rumen. Sato et al.  previously identified significant positive correlations between reticular and ruminal pH in both healthy and SARA-induced cows. SARA was diagnosed based on a condition characterized by ruminal pH < 5.6 for more than 3 h per day  in our study, and the ruminal and reticular pH remained under 5.6 over 3 h from 2 days after the high-grain diet feeding to the end of the experiment. Therefore, SARA was successfully induced in the present study. Furthermore, the lower proportion of acetic acid and higher proportion of butyric acid following the high-grain diet were consistent with a previous study . Although a temporarily higher concentration of lactic acid at 20:00 was identified in the CON period, this increase may not accurately reflect the diurnal changes in the 1-h mean ruminal and reticular pH throughout a day. Therefore, VFA production by microbes and its removal by absorption through ruminal epithelial cells, neutralization with salivary buffer, and passage to the lower digestive tracts  might be a more plausible determinant of the ruminal and reticular pH. However, no statistically significant differences in total VFA concentration were identified, so we can only speculate that fatty acid removal was encouraged by the marked decrease in ruminal and reticular pH in this study.
In our study, significant change in the relative abundance at the phylum level was only identified between the HAY and CON periods. Firmicutes and Bacteroidetes were the most abundant bacterial phyla in both the rumen and reticulum, and they were affected by a high-grain diet in our study. Gram-negative Bacteroidetes bacteria have been shown to decrease in the rumen of cattle with SARA induced by either grain or alfalfa pellets , and their relative abundance was reduced in cattle with induced SARA . Because a low rumen pH can lead to death and lysis of gram-negative bacteria , higher acidity during the CON period may cause a decrease in the relative abundance of Bacteroidetes. The gram-positive bacteria phylum Firmicutes is composed of species metabolically capable of consuming newly available fermentable carbohydrates ; therefore, an increase in their relative abundance may be related to the starch-rich high-grain diet and the decrease in the relative abundance of Bacteroidetes during the CON period.
At the genus level, changes in relative abundance were primarily identified between the high-forage and high-grain diets. Prevotella bacteria can use a variety of substrates in the rumen [12, 13], and was the predominant genus in our study. The lower relative abundance of Prevotella observed during the CON period is consistent with previous studies demonstrating that starch addition lowers the relative abundance of Prevotella . Although the relationship between substrate and growth of Prevotella is not well studied, a higher acidity during the CON period may affect the relative abundance of Prevotella. By contrast, the relative abundance of Ruminococcus was higher during the CON period. Although Ruminococcus is commonly known for cellulolytic species, several others (e.g., Ruminococcus bromii) can exploit starch . Therefore, the increase in the relative abundance of Ruminococcus during the CON period may be related to an increase in starch-fermenting Ruminococcus species. In addition, the genus Clostridium involved in the digestion of various substrates in the rumen  was higher during the CON period than during the HAY period, and it was consistent with a previous study .
Previously, Taguchi et al.  reported that Eubacterium ruminantium is associated with hemicellulose (mainly xylan) degradation in the rumen, and lower relative amount of Eubacterium ruminantium in the concentrate diet than the hay diet Tajima et al.  was consistent with the present study. In addition, Eubacterium was the only genus that differed significantly between the rumen and reticulum in this study. However, the only difference between the rumen and reticulum conditions that we suggest in this study was a slightly higher reticular pH than ruminal pH throughout the experiment. Although, the xylanase activity of Eubacterium ruminantium decreased in accordance with a decrease in pH (from 7.0 to 5.0; ), this could not be linked directly to the relative abundance of Eubacterium genus. Therefore, further studies are required to clarify the effects of acidic conditions on the changes in the relative abundance of Eubacterium genus. Collectively, these results suggest that changes in the relative abundances of bacterial genera may be mainly influenced by changes in useable substrates, as well as by the acidic environment.
The passage rate of ruminal contents from the rumen to the duodenum is increased in cattle by feeding a concentrated diet, resulting in less time for microbial fermentation . The increased passage rate could cause methanogenesis to shift to the hind gut and manure , and it was consistent with the lower copy number of total methanogens in the CON period. In this study, several bacteria well known for their cellulolytic activity (Fibrobacter succinogenes, Ruminococcus albus, and Ruminococcus flavefaciens; ) had higher copy numbers during the HAY period, which may be related to the higher fiber contents during this period. Meanwhile, Megasphaera elsdenii, a lactate-metabolizing species, increases in abundance as the bacterial community adapts to more readily fermentable carbohydrates , and the higher copy number of Megasphaera elsdenii during the CON period may be associated with higher relative abundance of lactate-producing bacterium, such as Olsenella, during the same period. Moreover, dilution of the reticulum contents occurs with fresh and less-fermented feed [2, 3], which might cause the lower bacterial copy numbers in the reticulum compared with the rumen.
Feeding a high-grain diet to cattle decreases the ruminal pH and bacterial diversity of the rumen epithelial community , and calves fed a high-grain diet show significantly lower ruminal pH and bacterial diversity compared with calves fed a high-grain diet with forage . In the present study, ruminal and reticular pH decreased after cattle transitioned from a high-forage to a high-grain diet, and bacterial diversity indices, such as OTUs, ACE, Chao1, and Shannon index, were lower during the CON period than during the HAY period. Furthermore, PCoA results showed that the HAY period was separate from the CON period in the plot regardless of location, and the rarefaction curve was consistent with the bacterial diversity analysis. Therefore, we assumed that the lower ruminal and reticular pH during the CON period resulted in the lower bacterial diversity. Although no statistically significant difference was identified in bacterial diversity between the rumen and reticulum, the higher diversity indices identified in the reticulum were consistent with the higher pH there compared to that in the rumen, which supports our hypothesis.
To the best of our knowledge, this is the first study to examine changes in reticular pH continuously and bacterial community structure following a high-grain diet feeding. Ruminal and reticular pH decreased after feeding a high-grain diet, and SARA was successfully induced during the CON period. Changes in the bacterial community structure and copy number were mainly identified between the two periods, which might be related to the changes in pH and diet. Furthermore, the decrease in pH identified during the CON period was consistent with the low bacterial diversity during the same period. These findings suggest that bacterial composition and diversity in both the rumen and reticulum were affected mainly by the acidic conditions and substrates useable for their growth, all of which may influence the fermentative ability of the rumen and reticulum.
Animals and experimental design
All animals were cared for according to protocols approved by Iwate University Laboratory Animal Care and Use Committee (A201401). Four rumen-cannulated Holstein bull cattle (192 ± 12 kg; 9.0 ± 1.4 months of age) were used in this study. Following the 7-day high-forage adaptation period, all cattle were fed a high-forage diet for 7 days (HAY period), and then a high-grain diet for 7 days (CON period). A SARA challenge was defined a HAY period followed by a CON period. During the experimental periods, cattle were fed mixed hay (orchard and timothy hay) in the HAY period and a high-grain diet with a forage-to-concentrate ratio of 19:81 (dry matter basis) during the CON period. The diet was supplied daily in two equal portions at 08:00 and 17:00. All cattle had free access to water and mineral salt block (E100TZ; ZENOAQ, Koriyama, Fukushima, Japan) throughout the study period. Daily total DMI was recorded for individual cattle throughout the experimental period, and all feed offered to the cattle was consumed. The amount and content of a high-forage (control) and a high-grain (SARA-inducing) diet were based on our previous study . The chemical compositions of the mixed hay and high-grain diet fed to the cattle are shown in Additional file 1: Table S1.
Sampling, measurements, and diagnosis of SARA
Ruminal and reticular pH were simultaneously measured every 10 min using a radio transmission system (YCOW-S; DKK-TOA Yamagata, Yamagata, Japan), as reported previously . The pH sensors were placed in the ventral sac of the rumen and the reticulum through the rumen fistula. Ruminal and reticular pH were continuously measured for 7 days during the HAY period (days 1–7) and consecutively for 7 days during the CON period (days 8–14). Using the rumen and reticulum accessible tubes, the fluid samples were collected from the rumen and reticulum adjacent to the pH sensor at 8:00, 14:00, and 20:00 on the last days of the HAY (day 7) and CON (day 14) periods, respectively.
The fluid samples were filtered immediately through two layers of cheesecloth after sampling. For the VFA analysis, 1 mL 25% HO3P in 3 N H2SO4 was added to 5 mL rumen and reticulum fluid. Total VFA and VFA components (i.e., acetic acid, propionic acid, and butyric acid) were separated and quantified by gas chromatography (Model 135, Hitachi, Tokyo, Japan) using a packed glass column (Thermon-3000, 3%) on a Shimalite TPA 60–80-mesh support (Shinwa Chemical Industries, Ltd., Kyoto, Japan). For the lactic acid analysis, the fluid samples were centrifuged at 2000×g for 15 min, and concentration in the supernatant were determined using a commercially available kit (F-kit; D-lactate/L-lactate, J. K. International Co., Tokyo, Japan). To measure NH3-N concentration, fluid samples were analyzed using the steam distillation method with an NH3-N analyzer (Kjeltec Auto Sampler System 1035 Analyzer, Tecator, Sweden).
Total bacterial DNA was extracted as described previously . Briefly, the fluid samples were incubated with 750 μg/mL lysozyme (Sigma-Aldrich, St. Louis, MO, USA) at 37 °C for 90 min. Then, 10 μL purified achromopeptidase (Wako Pure Chemical Industries, Ltd., Osaka, Japan) was added at a concentration of 10,000 U/mL and incubated at 37 °C for 30 min. The suspension was treated with 60 μL 1% sodium dodecyl sulfate and 1 mg/mL proteinase K (Merck Japan, Tokyo, Japan), and incubated at 55 °C for 5 min. The lysate was treated with phenol/chloroform/isoamyl alcohol (Wako Pure Chemical Industries, Ltd.) and chloroform (Life Technologies Japan, Ltd., Tokyo, Japan). DNA was precipitated by adding 5 M NaCl and 100% ethanol and centrifuged at 21,900×g for 15 min. The DNA pellet was rinsed with 70% ethanol, dried, and dissolved in TE buffer. The purified DNA was quantified using a Biospec-nano (Shimadzu, Kyoto, Japan) and stored at − 80 °C until further analysis.
The V1/V2 region of the 16S rRNA gene was amplified using a forward primer (5′-CCATCTCATCCCTGCGTGTCTCCGACTCAGNNNNNNNNNNAGRGTTTGATYMTGGCTCAG-3′) containing 454 primer A, a unique 10-bp barcode sequence for each sample (indicated as N), and 27Fmod (5′-AGRGTTTGATYMTGGCTCAG) in which the third base, A, in the original primer 27F was changed to R, as well as the reverse primer (5′-CCTATCCCCTGTGTGCCTTGGCAGTCTCAGTGCTGCCTCCCGTAGGAGT-3′) containing 454 primer B and reverse primer 338R (5′-TGCTGCCTCCCGTAGGAGT). Amplified products of ~ 370 bp were confirmed using agarose gel electrophoresis, purified using AMPure XP magnetic purification beads (Beckman Coulter, Inc., Brea, CA, USA), and quantified using the Quant-iT PicoGreen dsDNA Assay Kit (Life Technologies Japan). Mixed samples were prepared by pooling approximately equal amounts of PCR amplicons from each sample, and then subjected to 454 GS Junior (Roche Applied Science, Indianapolis, IN, USA) sequencing following the manufacturer’s instructions.
Pyrosequencing data analysis
All pyrosequencing reads were filtered according to the procedure of Kim et al. , who developed an analysis pipeline for barcoded 454 pyrosequencing of PCR amplicons in V1/V2, the region amplified by the 27Fmod/338R primers. A total of 151,797 filter-passed reads were processed using MOTHUR (ver. 1.35, University of Michigan; http://www.mothur.org/wiki/; ), and all samples were standardized by random subsampling to 4522 sequences per sample using the “sub.sample” command to generate rarefaction curves and calculate ACE, Chao1 richness estimator, and Shannon diversity index, according to the Illumina MiSeq protocol described previously . Unique sequences were determined and used to align against the SILVA reference alignment database ; chimera were removed using chimera.uchime (http://drive5.com/uchime); sequences identified as being of eukaryotic origin were removed; the candidate sequences were screened and preclustered to eliminate outliers; and a distance matrix was generated from the resulting sequences. Sequences were clustered into OTU, with a cutoff of 97% similarity. A rarefaction curve was generated at the 97% similarity level, which was calculated by the distance-based OTU . To calculate the nonparametric species richness estimators Chao 1 and ACE and the Shannon diversity index, the “summary.single” command was used. The unweighted UniFrac distance method  was used to perform a principal coordinates analysis (PCoA) with all OTU.
Real-time quantitative PCR
Quantitative real-time PCR (qRT-PCR) was performed to evaluate the copy number of methyl-coenzyme M reductase α-subunit (mcrA) from total methanogens, and 16S rRNA genes from Fibrobacter succinogenes, Megasphaera elsdenii, Ruminococcus albus, Ruminococcus flavefaciens, Streptococcus bovis, and Selenomonas ruminantium using SYBR green (iQ SYBR Green Supermix, Bio-Rad, Hercules, CA, USA) with the MiniOpticon Real-Time PCR system (Bio-Rad). Primer pairs (Additional file 1: Table S2) were selected to detect bacterial species closely associated with dietary changes and other bacterial species. Each sample contained 10 ng DNA, 2× SYBR green, and 0.6 μM each primer in a final volume of 20 μL. Amplification conditions were as follows: 95 °C for 3 min; 40 cycles of 10 s at 95 °C; 20 s at 63 °C (for total methanogens), 60 °C (for Fibrobacter succinogenes), 58 °C (for Megasphaera elsdenii), 57 °C (for Streptococcus bovis and Selenomonas ruminantium), or 55 °C (for Ruminococcus albus and Ruminococcus flavefaciens); and 30 s at 72 °C. The fluorescence signal was collected at the end of each cycle. To obtain melting curve data, the temperature was increased in 0.5 °C increments from 65 to 94 °C. A standard curve for each primer pair was constructed from recombinant plasmid DNA containing 16S rRNA inserts of DNA purified from a pure culture of the target species. The strains used for plasmid preparation were as follows: Methanobrevibacter ruminantium JCM13430 (DSM1093), Fibrobacter succinogenes ATCC19169, Megasphaera elsdenii ATCC25940, Ruminococcus albus ATCC27210, Ruminococcus flavefaciens ATCC19208, Streptococcus bovis ATCC33317, and Selenomonas ruminantium ATCC12561. Plasmid DNA was quantified and subjected to seven sequential 10-fold dilutions. Data were collected and processed using CFX Manager software ver. 1.5 (Bio-Rad).
Ruminal and reticular pH data and the duration of pH < 5.6 were summarized as 24-h means during the HAY and CON periods. The normality of the distribution of variables was tested using Shapiro-Wilk test, and non-normal data (the proportion of other VFAs, Shannon index, the relative abundances of Butyrivibrio, Eubacterium, unclassified Clostridiales, Acetivibrio, unclassified Lachnospiraceae, Olsenella, unclassified Clostridiaceae, Bacteroides, and the copy numbers of Megasphaera elsdenii, Streptococcus bovis, and Ruminococcus flavefacience) were root-square transformed before analysis. One-way repeated measures ANOVA followed by Tukey multiple comparison method was used to determine the significance of the difference in 24-h and 1-h means pH, the duration of pH < 5.6, total VFA concentration, the proportions of individual VFA, and the concentrations of NH3-N and lactic acid to compare rumen × reticulum in a same period and HAY × CON periods in a same location. The relative abundances of bacteria phyla and genera, bacterial diversity indices, and bacterial species copy numbers were analyzed using two-way repeated-measures ANOVA and Tukey’s multiple range tests. The statistical model included the fixed effects of location (rumen and reticulum), period (HAY and CON periods), and their interactions, plus the random effect of animal. Linear regression analysis of 24-h and 1-h mean ruminal and reticular pH against days and hours was performed to determine whether there are significant changes in pH during the experimental period. All numerical data were analyzed using Prism ver. 7.01 (GraphPad Software, Inc., La Jolla, CA, USA), and are expressed as means ± standard error of mean (SEM). Differences were considered significant at P < 0.05, and trends suggesting possible significance were determined at 0.05 ≤ P < 0.10.
Abundance-based coverage estimator
Analysis of variance
Operational taxonomic unit
Principal coordinate analysis
Quantitative real time polymerase chain reaction
Subacute ruminal acidosis
Volatile fatty acid
Gozho GN, Plaizier JC, Krause DO, Kennedy AD, Wittenberg KM. Subacute ruminal acidosis induces ruminal lipopolysaccharide endotoxin release and triggers an inflammatory response. J Dairy Sci. 2005;88:1399–403.
Sato S, Ikeda A, Tsuchiya Y, Ikuta K, Murayama I, Kanehira M, Okada K, Mizuguchi H. Diagnosis of subacute ruminal acidosis (SARA) by continuous reticular pH measurements in cows. Vet Res Commun. 2012a;36:201–5. https://doi.org/10.1007/s11259-012-9528-8.
Falk M, Münger A, Dohme-Meier F. Technical note: a comparison of reticular and ruminal pH monitored continuously with 2 measurement systems at different weeks of early lactation. J Dairy Sci. 2016;99:1951–5. https://doi.org/10.3168/jds.2015-9725.
Zened A, Combes S, Cauquil L, Mariette J, Klopp C, Bouchez O, Troegeler-Meynadier A, Enjalbert F. Microbial ecology of the rumen evaluated by 454 GS FLX pyrosequencing is affected by starch and oil supplementation of diets. FEMS Microbiol Ecol. 2013;83:504–14. https://doi.org/10.1111/1574-6941.12011.
Sato S. Pathophysiological evaluation of subacute ruminal acidosis (SARA) by continuous ruminal pH monitoring. Anim Sci J. 2016;87:168–77. https://doi.org/10.1111/asj.12415.
Kim Y-H, Nagata R, Ohtani N, Ichijo T, Ikuta K, Sato S. Effects of dietary forage and calf starter diet on ruminal pH and bacteria in Holstein calves during weaning transition. Front Microbiol. 2016;7:1575. https://doi.org/10.3389/fmicb.2016.01575.
Liu JH, Bian GR, Zhu WY, Mao SY. High-grain feeding causes strong shifts in ruminal epithelial bacterial community and expression of toll-like receptor genes in goats. Front Microbiol. 2015;6:167. https://doi.org/10.3389/fmicb.2015.00167.
Kimura A, Sato S, Kato T, Ikuta K, Yamagishi N, Okada K, Mizuguchi H, Ito K. Relationship between pH and temperature in the ruminal fluid of cows, based on a radio-transmission pH-measurement system. J Vet Med Sci. 2012;74:1023–8. https://doi.org/10.1292/jvms.12-0084.
Mao SY, Zhang RY, Wang DS, Zhu WY. Impact of subacute ruminal acidosis (SARA) adaptation on rumen microbiota in dairy cattle using pyrosequencing. Anaerobe. 2013;24:12–9. https://doi.org/10.1016/j.anaerobe.2013.08.003.
Allen MS. Relationship between fermentation acid production in the rumen and the requirement for physically effective fiber. J Dairy Sci. 1997;80:1447–62. https://doi.org/10.3168/jds.S0022-0302(97)76074-0.
Khafipour E, Li S, Plaizier JC, Krause DO. Rumen microbiome composition determined using two nutritional models of subacute ruminal acidosis. Appl Environ Microbiol. 2009;75:7115–24. https://doi.org/10.1128/AEM.00739-09.
Stevenson DM, Weimer PJ. Dominance of Prevotella and low abundance of classical ruminal bacterial species in the bovine rumen revealed by relative quantification real-time PCR. Appl Microbiol Biotechnol. 2007;75:165–74. https://doi.org/10.1007/s00253-006-0802-y.
Bekele AZ, Koike S, Kobayashi Y. Genetic diversity and diet specificity of ruminal Prevotella revealed by 16S rRNA gene-based analysis. FEMS Microbiol Lett. 2010;305:49–57. https://doi.org/10.1111/j.1574-6968.2010.01911.x.
Klieve AV, O’Leary MN, McMillen L, Ouwerkerk D. Ruminococcus bromii, identification and isolation as a dominant community member in the rumen of cattle fed a barley diet. J Appl Microbiol. 2007;103:2065–73. https://doi.org/10.1111/j.1365-2672.2007.03492.x.
Pitta DW, Pinchak WE, Dowd S, Dorton K, Yoon I, Min BR, Fulford JD, Wickersham TA, Malinowski DP. Longitudinal shifts in bacterial diversity and fermentation pattern in the rumen of steers grazing wheat pasture. Anaerobe. 2014;30:11–7. https://doi.org/10.1016/j.anaerobe.2014.07.008.
Taguchi H, Koike S, Kobayashi Y, Cann IK, Karita S. Partial characterization of structure and function of a xylanase gene from the rumen hemicellulolytic bacterium Eubacterium ruminantium. Anim Sci J. 2004;75:325–32. https://doi.org/10.1111/j.1740-0929.2004.00193.x.
Tajima K, Aminov RI, Nagamine T, Matsui H, Nakamura M, Benno Y. Diet-dependent shifts in the bacterial population of the rumen revealed with real-time PCR. Appl Environ Microbiol. 2001;67:2766–74. https://doi.org/10.1128/AEM.67.6.2766-2774.2001.
Huhtanen P, Kukkonen U. Comparison of methods, markers, sampling sites and models for estimating digesta passage kinetics in cattle fed at two levels of intake. Anim Feed Sci Technol. 1995;52:141–58. https://doi.org/10.1016/0377-8401(94)00699-A.
Hindrichsen IK, Wettstein HR, Machmüller A, Kreuzer M. Methane emission, nutrient degradation and nitrogen turnover in dairy cows and their slurry at different milk production scenarios with and without concentrate supplementation. Agric Ecosys Environ. 2006;113:150–61. https://doi.org/10.1016/j.agee.2005.09.004.
Koike S, Kobayashi Y. Development and use of competitive PCR assays for the rumen cellulolytic bacteria: Fibrobacter succinogenes, Ruminococcus albus, and Ruminococcus flavefaciens. FEMS Microbiol Lett. 2001;204:361–6. https://doi.org/10.1111/j.1574-6968.2001.tb10911.x.
Huber TL. Physiological effects of acidosis on feedlot cattle. J Anim Sci. 1976;43:902–9.
Sato S, Kimura A, Anan T, Yamagishi N, Okada K, Mizuguchi H, Ito K. A radio transmission pH measurement system for continuous evaluation of fluid pH in the rumen of cows. Vet Res Commun. 2012b;36:85–9. https://doi.org/10.1007/s11259-012-9518-x.
Kim SW, Suda W, Kim S, Oshima K, Fukuda S, Ohno H, Morita H, Hattori M. Robustness of gut microbiota of healthy adults in response to probiotic intervention revealed by high-throughput pyrosequencing. DNA Res. 2013;20:241–53. https://doi.org/10.1093/dnares/dst006.
Schloss PD, Westcott SL, Ryabin T, Hall JR, Hartmann M, Hollister EB, Lesniewski RA, Oakley BB, Parks DH, Robinson CJ, Sahl JW, Stres B, Thallinger GG, Van Horn DJ, Weber CF. Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl Environ Microbiol. 2009;75:7537–41. https://doi.org/10.1128/AEM.01541-09.
Kozich JJ, Westcott SL, Baxter NT, Highlander SK, Schloss PD. Development of a dual-index sequencing strategy and curation pipeline for analyzing amplicon sequence data on the MiSeq Illumina sequencing platform. Appl Environ Microbiol. 2013;79:5112–20. https://doi.org/10.1128/AEM.01043-13.
Pruesse E, Quast C, Knittel K, Fuchs BM, Ludwig W, Peplies J, Glöckner FO. SILVA: a comprehensive online resource for quality checked and aligned ribosomal RNA sequence data compatible with ARB. Nucl Acids Res. 2007;35:7188–96. https://doi.org/10.1093/nar/gkm864.
Schloss PD, Gevers D, Westcott SL. Reducing the effects of PCR amplification and sequencing artifacts on 16S rRNA-based studies. PLoS One. 2011;6:e27310. https://doi.org/10.1371/journal.pone.0027310.
Lozupone C, Knight R. UniFrac: a new phylogenetic method for comparing microbial communities. Appl Environ Microbiol. 2005;71:8228–35. https://doi.org/10.1128/AEM.71.12.8228-8235.2005.
The authors thank Dr. Yasuo Kobayashi, Research Faculty of Agriculture, Hokkaido University, who kindly provided the recombinant DNA for qRT-PCR analysis.
This work was supported by the Cabinet Office, Government of Japan, Cross-ministerial Strategic Innovation Promotion Program (SIP), “Technologies for creating next-generation agriculture, forestry and fisheries.”
Availability of data and materials
The sequencing data were deposited into the Sequence Read Archive (SRA) of NCBI and can be accessed via accession number SRP093895.
All animals were cared for according to protocols approved by Iwate University Laboratory Animal Care and Use Committee (A201401).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. Compositions of the high-forage and high-grain diet on percentage (%) and dry matter bases. Table S2. Primers sequences used for qRT-PCR. Figure S1. Rarefaction curves calculated from the 454 pyrosequencing data at a 97% similarity level in the rumen and reticulum. Cattle were fed a high-forage diet until day 7 (HAY period) and a high-grin diet until day 14 (CON period). Solid and dash dot lines represent the rumen of the HAY and CON periods, respectively. Dash and dot lines represent the reticulum of the HAY and CON periods, respectively. (DOCX 441 kb)