Microbial composition, rumen fermentation parameters, enteric methane emissions, and lactational performance of phenotypically high and low methane-emitting dairy cows.

This experiment was designed to investigate the relation of high and low methane-yield phenotypes with body weight (BW), dry matter intake (DMI), lactation performance, enteric CH4 emissions, and rumen fermentation parameters in lactating dairy cows. A total of 130 multi- and primiparous Holstein cows were screened for enteric CH4 emissions using the GreenFeed system (C-Lock Inc.). Out of these 130 cows, 5 were identified as phenotypically high (HM) and 5 as phenotypically low (LM) CH4 emitters. Cows in the LM group had lower daily enteric CH4 emissions than cows in the HM group (on average 346 vs. 439 g/d, respectively), lower CH4 yield (15.5 vs. 20.4 g of CH4/kg of DMI), and CH4 intensity (13.2 vs. 17.0 g of CH4/ kg of energy-corrected milk yield). Enteric emissions of CO2 and H2 did not differ between HM and LM cows. These 10 cows were blocked by parity, days in milk, and milk production, and were used in a 5-wk randomized complete block design experiment. Milk composition, production, and BW were also not different between LM and HM cows. The concentration of total volatile fatty acids in ruminal contents did not differ between CH4 phenotypes, but LM cows had a lower molar proportion of acetate (57 vs. 62.1%), a higher proportion of propionate (27.5 vs. 21.6%, respectively), and therefore a lower acetate-to-propionate ratio than HM cows. Consistently, the 16S cDNA analysis revealed the abundance of Succinivibrionaceae and unclassified Veillonellaceae to be higher in LM cows compared with HM cows, bacteria that were positively correlated with ruminal propionate concentration. Notably, Succinivibrionaceae trigger the formation of propionate via oxaloacetate pathway from phosphoenolpyruvate via Enzyme Commission: 4.1.1.49, which showed a trend to be higher in LM cows compared with HM cows. Additionally, LM cows possessed fewer transcripts of a gene encoding for methyl-CoM reductase enzyme compared with HM. In this study, low and high CH4-yield cows have similar production performance and milk composition, but total-tract apparent digestibility of organic matter and fiber fractions was lower in the former group of animals.

Temporal changes in total and metabolically active ruminal methanogens in dairy cows supplemented with 3-nitrooxypropanol.

The purpose of this study was to investigate the effect of 3-nitrooxypropanol (3-NOP), a potent methane inhibitor, on total and metabolically active methanogens in the rumen of dairy cows over the course of the day and over a 12-wk period. Rumen contents of 8 ruminally cannulated early-lactation dairy cows were sampled at 2, 6, and 10 h after feeding during wk 4, 8, and 12 of a randomized complete block design experiment in which 3-NOP was fed at 60 mg/kg of feed dry matter. Cows (4 fed the control and 4 fed the 3-NOP diet) were blocked based on their previous lactation milk yield or predicted milk yield. Rumen samples were extracted for microbial DNA (total) and microbial RNA (metabolically active), PCR amplified for the 16S rRNA gene of archaea, sequenced on an Illumina platform, and analyzed for archaea diversity. In addition, the 16S copy number and 3 ruminal methanogenic species were quantified using the real-time quantitative PCR assay. We detected a difference between DNA and RNA (cDNA)-based archaea communities, revealing that ruminal methanogens differ in their metabolic activities. Within DNA and cDNA components, methanogenic communities differed by sampling hour, week, and treatment. Overall, Methanobrevibacter was the dominant genus (94.3%) followed by Methanosphaera, with the latter genus having greater abundance in the cDNA component (14.5%) compared with total populations (5.5%). Methanosphaera was higher at 2 h after feeding, whereas Methanobrevibacter increased at 6 and 10 h in both groups, showing diurnal patterns among individual methanogenic lineages. Methanobrevibacter was reduced at wk 4, whereas Methanosphaera was reduced at wk 8 and 12 in cows supplemented with 3-NOP compared with control cows, suggesting differential responses among methanogens to 3-NOP. A reduction in Methanobrevibacter ruminantium in all 3-NOP samples from wk 8 was confirmed using real-time quantitative PCR. The relative abundance of individual methanogens was driven by a combination of dietary composition, dry matter intake, and hydrogen concentrations in the rumen. This study provides novel information on the effects of 3-NOP on individual methanogenic lineages, but further studies are needed to understand temporal dynamics and to validate the effects of 3-NOP on individual lineages of ruminal methanogens.

Comparing noninvasive sampling techniques with standard cannula sampling method for ruminal microbial analysis.

Rumen microbes play an important role in the conversion of indigestible plant material to energy and protein in dairy cows. Sampling for ruminal contents via cannula is considered the gold standard technique for microbial analysis, but the technique requires ruminally cannulated animals and specialized animal facilities. The purpose of this study was to determine whether other sampling methods and locations along the digestive tract may serve as noninvasive proxies to the cannula method for microbial analysis. Six ruminally cannulated lactating Holstein dairy cows were adapted to a standard total mixed ration for 2 wk and sampled during the third week. Sampling locations and methods included salivary content, rumination bolus (regurgitated digesta collected from the cow's mouth), feces, and rumen contents via stomach tube and cannula. Stomach tube and cannula samples differ in proportions of solid and liquid material and were therefore separated into whole (as collected), liquid, and solid fractions. Samples were collected at 0 (before feeding), 2, 4, 6, 8, and 12 h after feeding over 2 d. All samples were extracted for total genomic DNA and selected samples for metabolically active DNA (RNA), PCR-amplified for the V1-V2 region of the 16S rRNA bacterial gene, and analyzed for bacterial diversity using the QIIME2 pipeline followed by statistical analysis in R (https://www.R-project.org/). In DNA-based analysis, at the community level, saliva, rumination bolus, and fecal samples clustered in separate groups, whereas all fractions of stomach tube and cannula samples clustered together, indicating that microbial communities of stomach tube and cannula samples were homogeneous. Rumination bolus samples at 6, 8, and 12 h after feeding clustered with stomach tube and cannula samples, indicating that rumination bolus samples may be an alternative for cannula samples; however, time of sampling is critical for sampling of bolus digesta. Results of the RNA-based analysis of rumination bolus samples and solid samples from cannula and stomach tube at 0 and 6 h after feeding were similar. We concluded that the solid fraction of samples obtained via the stomach tube method may serve as a proxy for the solid fraction of whole ruminal contents obtained via cannula for DNA-based microbial investigations. Both rumination bolus and stomach tube solid samples may serve as proxies for cannula solid samples for RNA-based microbial analysis.

Effect of 2-hydroxy-4-(methylthio) butanoate (HMTBa) supplementation on rumen bacterial populations in dairy cows when exposed to diets with risk for milk fat depression.

Diet-induced milk fat depression (MFD) is a condition marked by a reduction in milk fat yield experimentally achieved by increasing dietary unsaturated fatty acids and fermentable carbohydrates. 2-Hydroxy-4-(methylthio) butanoate (HMTBa) is a methionine analog observed to reduce diet-induced MFD in dairy cows. We hypothesize that the reduction in diet-induced MFD by HMTBa is due to changes in the rumen microbiota. To test this, 22 high-producing cannulated Holstein dairy cows were placed into 2 groups using a randomized block design and assigned to either control or HMTBa supplementation (0.1% of diet dry matter). All cows were then exposed to 3 different diets with a low risk (32% neutral detergent fiber, no added oil; fed d 1 to 7), a moderate risk (29% neutral detergent fiber and 0.75% soybean oil; fed d 8 to 24), or a high risk (29% neutral detergent fiber and 1.5% soybean oil; fed d 25 to 28) for diet-induced MFD. Rumen samples were collected on d 0, 14, 24, and 28, extracted for DNA, PCR-amplified for the V1-V2 region of the 16S rRNA gene, sequenced on an Illumina MiSeq (Illumina, San Diego, CA), and subjected to bacterial diversity analysis using the QIIME pipeline. The α diversity estimates (species richness and Shannon diversity) were decreased in the control group compared with the HMTBa group. Bacterial community composition also differed between control and HMTBa groups based on both weighted UniFrac (relative abundance of commonly detected bacteria) and unweighted UniFrac (presence/absence) distances. Within the HMTBa group, no differences were observed in bacterial community composition between d 0 and d 14, 24, and 28; however, in the control group, d 0 samples were different from d 14, 24, and 28. Certain bacterial genera including Dialister, Megasphaera, Lachnospira, and Sharpea were increased in the control group compared with the HMTBa group. Interestingly, these genera were positively correlated with milk fat trans-10,cis-12 conjugated linoleic acid and trans-10 C18:1, fatty acid isomers associated with biohydrogenation-induced MFD. It can be concluded that diet-induced MFD is accompanied by significant alterations in the rumen bacterial community and that HMTBa supplementation reduces these microbial perturbations.

Alterations in ruminal bacterial populations at induction and recovery from diet-induced milk fat depression in dairy cows.

Ten ruminally cannulated Holstein cows were used in a crossover design that investigated changes in ruminal bacterial populations in response to induction and recovery from diet-induced milk fat depression (MFD). Further, the effect on the ruminal microbiota of the cows with diet-induced milk fat depression inoculated with rumen contents from non-milk fat-depressed donor cows was evaluated. Milk fat depression was induced during the first 10 d of each period by feeding a low-fiber, high-starch, and high-polyunsaturated fatty acid diet (26.1% neutral detergent fiber, 28.1% starch, 5.8% total fatty acids, and 1.9% C18:2), resulting in a 30% decrease in milk fat yield. Induction was followed by a recovery phase, where all cows were switched to a high-fiber, low-starch, and low-polyunsaturated fatty acid diet (31.8% neutral detergent fiber, 23% starch, 4.2% total fatty acids, and 1.2% C18:2) and were allocated to (1) control (no inoculation) or (2) ruminal inoculation with donor cow digesta (8 kg/d for 6 d). Ruminal samples were collected at the end of induction (d 10) and during recovery (d 13, 16, and 28), separated to solid and liquid fractions, extracted for DNA, PCR- amplified for the V1-V2 region of the 16S rRNA gene, and analyzed for bacterial diversity. Results indicated that bacterial communities were different between fractions. In each fraction, differences were significant between the induction (d 10) and recovery (d 13, 16, and 28) periods; however, differences were less apparent with time during the recovery period. The MFD (d 10) was typified by a reduction in the relative sequence abundance of Bacteroidetes and an increase in the relative sequence abundance of Firmicutes and Actinobacteria across both fractions. At the genus level, relative sequence abundance of unclassified Lachnospiraceae, Butyrivibrio, Bulleidia, and Coriobacteriaceae were higher on d 10 and were positively correlated with trans-10,cis-12 CLA and the trans-10 isomer, suggesting their potential role in altered biohydrogenation reactions. A switch to the recovery diet resulted in a sharp increase in the Bacteroidetes lineages and a decrease in Firmicutes members on d 13; however, this shift appears to stabilize by d 28, indicating the restoration process for ruminal bacteria from an altered state is gradual and complex. Inoculation of 10% of rumen contents from non-MFD donor cows to MFD cows revealed this procedure had transient effects on only a few bacterial populations, and such effects disappeared after d 16 following cessation of inoculation. It can be concluded that alterations in milk FA profiles at induction are preceded by microbial alterations in the rumen driven by dietary changes.

Metagenomic Analysis of the Rumen Microbiome of Steers with Wheat-Induced Frothy Bloat.

Frothy bloat is a serious metabolic disorder that affects stocker cattle grazing hard red winter wheat forage in the Southern Great Plains causing reduced performance, morbidity, and mortality. We hypothesize that a microbial dysbiosis develops in the rumen microbiome of stocker cattle when grazing on high quality winter wheat pasture that predisposes them to frothy bloat risk. In this study, rumen contents were harvested from six cannulated steers grazing hard red winter wheat (three with bloat score "2" and three with bloat score "0"), extracted for genomic DNA and subjected to 16S rDNA and shotgun sequencing on 454/Roche platform. Approximately 1.5 million reads were sequenced, assembled and assigned for phylogenetic and functional annotations. Bacteria predominated up to 84% of the sequences while archaea contributed to nearly 5% of the sequences. The abundance of archaea was higher in bloated animals (P < 0.05) and dominated by Methanobrevibacter. Predominant bacterial phyla were Firmicutes (65%), Actinobacteria (13%), Bacteroidetes (10%), and Proteobacteria (6%) across all samples. Genera from Firmicutes such as Clostridium, Eubacterium, and Butyrivibrio increased (P < 0.05) while Prevotella from Bacteroidetes decreased in bloated samples. Co-occurrence analysis revealed syntrophic associations between bacteria and archaea in non-bloated samples, however; such interactions faded in bloated samples. Functional annotations of assembled reads to Subsystems database revealed the abundance of several metabolic pathways, with carbohydrate and protein metabolism well represented. Assignment of contigs to CaZy database revealed a greater diversity of Glycosyl Hydrolases dominated by oligosaccharide breaking enzymes (>70%) in non-bloated samples. However, the abundance and diversity of CaZymes were greatly reduced in bloated samples indicating the disruption of carbohydrate metabolism. We conclude that mild to moderate frothy bloat results from tradeoffs both within and between microbial domains due to greater competition for substrates that are of limited availability as a result of biofilm formation.

Longitudinal shifts in bacterial diversity and fermentation pattern in the rumen of steers grazing wheat pasture.

Grazing steers on winter wheat forage is routinely practiced in the Southern Great Plains of the US. Here, we investigated the dynamics in bacterial populations of both solid and liquid ruminal fractions of steers grazing on maturing wheat forage of changing nutritive quality. The relationship between bacterial diversity and fermentation parameters in the liquid fraction was also investigated. During the first 28 days, the wheat was in a vegetative phase with a relatively high crude protein content (CP; 21%), which led to the incidence of mild cases of frothy bloat among steers. Rumen samples were collected on days 14, 28, 56 and 76, separated into solid and liquid fractions and analyzed for bacterial diversity using 16S pyrotag technology. The predominant phyla identified were Bacteroidetes (59-77%) and Firmicutes (20-33%) across both ruminal fractions. Very few differences were observed in the rumen bacterial communities within solid and liquid fractions on day 14. However, by day 28, the relatively high CP content complemented a distinct bacterial and chemical composition of the rumen fluid that was characterized by a higher ratio (4:1) of Bacteroidetes:Firmicutes and a corresponding lower acetate:propionate (3:1) ratio. Further, a greater accumulation of biofilm (mucopolysaccharide complex) on day 28 was strongly associated with the abundance of Firmicutes lineages such as Clostridium, Ruminococcus, Oscillospira and Moryella (P<0.05) in the fiber fraction. Such changes were diminished as the CP concentration declined over the course of the study. The abundance of Firmicutes was noticeable by 76 d in both fractions which signifies the development of a core microbiome associated with digestion of a more recalcitrant fiber in the mature wheat. This study demonstrates dynamics in the rumen microbiome and their association with fermentation activity in the rumen of steers during the vegetative (bloat-prone) and reproductive stages of wheat forage.

Temporal dynamics in the ruminal microbiome of dairy cows during the transition period.

The transition period in dairy cows refers to the period from 3 wk before calving to 3 wk post-calving and is a critical time for influencing milk production and cow health. We hypothesize that the ruminal microbiome shifts as dairy cows transition from a non-lactation period into lactation due to changes in dietary regimen. The purpose of this study was to identify differences in the ruminal microbiome of primiparous and multiparous (study group) cows during the transition period. Five primiparous and 5 multiparous cows were randomly selected from a herd, and ruminal contents were sampled, via stomach tube, 4 times (study day) at 3 wk before calving date (S1), 1 to 3 d post-calving (S2), and 4 (S3) and 8 wk (S4) into lactation and were evaluated for bacterial diversity using 16S pyrotags. Both groups received the same pre-fresh diet (14.6% CP, 44.0% NDF, 21.9% starch) and 3 different lactation diets (L1, L2, and L3) varying in forage base but not amount and formulated to have similar nutrient specifications (16.8% to 17.7% CP; 32.5% to 33.6% NDF; 26.2% to 29.1% starch) post-calving. Forty bacterial communities were analyzed on the basis of annotations of 100,000 reads, resulting in 15,861 operational taxonomic units grouped into 17 bacterial phyla. The UniFrac distance metric revealed that both study group and study day had an effect on the community compositions (P < 0.05; permutational multivariate ANOVA test). The most abundant phyla observed were Bacteroidetes and Firmicutes across all the communities. As the cows transitioned into lactation, the ratio of Bacteroidetes to Firmicutes increased from 6:1 to 12:1 (P < 0.05; Mann-Whitney U test), and this ratio was greater in primiparous cows than in multiparous cows (P < 0.05). This report is the first to explore the effect of parity on dynamics in the ruminal microbiome of cows during the transition period.

Bacterial diversity associated with feeding dry forage at different dietary concentrations in the rumen contents of Mehshana buffalo (Bubalus bubalis) using 16S pyrotags.

Pyrosequencing of 16S rRNA gene targeting bacteria was applied to identify diet-induced shifts in the microbiome of both solid and liquid ruminal fractions retrieved from water buffalo fed different diets. The depth of coverage of metabolically active bacteria in a community using different primer pairs was also investigated. To assess reproducibility, animal to animal variation was considered in all phylogenetic and community comparisons. The experiment included four non-lactating water buffaloes fed three different diets for six weeks each; diets were M1 (50% concentrate: 50% dry roughage), M2 (25% concentrate: 75% dry roughage) and M3 (100% dry roughage). A total of 333, 851 pyrotags were analyzed in this study. Phylogenetic analysis revealed significant differences in the rumen microbiome mediated by primer and diet (P < 0.05). Differences in community composition due to primer, diet, fraction and animal were compared using unweighted and weighted UniFrac analysis. Clustering of communities was largely explained by primer differences in both weighted and unweighted UniFrac analyses (P < 0.001). In the weighted analysis, communities clustered by diets (P < 0.05) and fractions (P < 0.08) while no inter-animal variation was observed. The identified repertoire of bacterial populations was dependent on the primer pair, as targeting the V4-V5 region resulted in greater diversity profiles of the microbiome. Within each primer pair, dietary changes altered the community composition with noticeable shifts at genus level. Genera such as Ruminococcus and Fibrobacter (P < 0.05) were higher in abundance on M3 diet while Prevotella dominated (P < 0.05) on M1 diet.