Differences in the rumen methanogen populations of lactating Jersey and Holstein dairy cows under the same diet regimen

Research progress in isolation and identification of rumen probiotics

With the increasing research on the exploitation of rumen microbial resources, rumen probiotics have attracted much attention for their positive contributions in promoting nutrient digestion, inhibiting pathogenic bacteria, and improving production performance. In the past two decades, macrogenomics has provided a rich source of new-generation probiotic candidates, but most of these "dark substances" have not been successfully cultured due to the restrictive growth conditions. However, fueled by high-throughput culture and sorting technologies, it is expected that the potential probiotics in the rumen can be exploited on a large scale, and their potential applications in medicine and agriculture can be explored. In this paper, we review and summarize the classical techniques for isolation and identification of rumen probiotics, introduce the development of droplet-based high-throughput cell culture and single-cell sequencing for microbial culture and identification, and finally introduce promising cultureomics techniques. The aim is to provide technical references for the development of related technologies and microbiological research to promote the further development of the field of rumen microbiology research.

Rumen Development of Tianhua Mutton Sheep Was Better than That of Gansu Alpine Fine Wool Sheep under Grazing Conditions

The purpose of this experiment was to investigate the differences in rumen tissue morphology, volatile fatty acid content, and rumen microflora between Tianhua mutton sheep and Gansu alpine fine wool sheep under the same grazing conditions. Twelve 30-day-old lambs were randomly selected from two different flocks in Duolong Village and grazed together for a period of 150 days. The rumen tissue was fixed with 4% paraformaldehyde and brought back to the laboratory for H&E staining, the volatile fatty acid content of the rumen contents was detected by gas chromatography, and the rumen flora structure was sequenced by full-length sequencing of the bacterial 16S rRNA gene using the PacBio sequencing platform. The acetic acid and total acid contents of the rumen contents of Tianhua mutton sheep were significantly higher than those of Gansu alpine fine wool sheep (p < 0.05). The rumen papillae height of Tianhua mutton sheep was significantly higher than that of Gansu alpine fine wool sheep (p < 0.05). The diversity and richness of the rumen flora of Tianhua mutton sheep were higher than those of Gansu alpine fine wool sheep, and Beta analysis showed that the microflora structure of the two fine wool sheep was significantly different. At the phylum level, Firmicutes and Bacteroidetes dominated the rumen flora of Tianhua mutton sheep and Gansu alpine fine wool sheep. At the genus level, the dominant strains were Christensenellaceae_R_7_group and Rikenellaceae_RC9_gut_group. LEfSe analysis showed that Prevotella was a highly abundant differential species in Tianhua mutton sheep and lachnospiraccac was a highly abundant differential species in Gansu alpine fine wool sheep. Finally, both the KEGG and COG databases showed that the enrichment of biometabolic pathways, such as replication and repair and translation, were significantly higher in Tianhua mutton sheep than in Gansu alpine fine wool sheep (p < 0.05). In general, there were some similarities between Tianhua mutton sheep and Gansu alpine fine wool sheep in the rumen tissue morphology, rumen fermentation ability, and rumen flora structure. However, Tianhua mutton sheep had a better performance in the rumen acetic acid content, rumen papillae height, and beneficial bacteria content. These differences may be one of the reasons why Tianhua mutton sheep are more suitable for growing in alpine pastoral areas than Gansu alpine fine wool sheep.

Unraveling the phylogenomic diversity of Methanomassiliicoccales and implications for mitigating ruminant methane emissions

BACKGROUND

Methanomassiliicoccales are a recently identified order of methanogens that are diverse across global environments particularly the gastrointestinal tracts of animals; however, their metabolic capacities are defined via a limited number of cultured strains.

RESULTS

Here, we profile and analyze 243 Methanomassiliicoccales genomes assembled from cultured representatives and uncultured metagenomes recovered from various biomes, including the gastrointestinal tracts of different animal species. Our analyses reveal the presence of numerous undefined genera and genetic variability in metabolic capabilities within Methanomassiliicoccales lineages, which is essential for adaptation to their ecological niches. In particular, gastrointestinal tract Methanomassiliicoccales demonstrate the presence of co-diversified members with their hosts over evolutionary timescales and likely originated in the natural environment. We highlight the presence of diverse clades of vitamin transporter BtuC proteins that distinguish Methanomassiliicoccales from other archaeal orders and likely provide a competitive advantage in efficiently handling B12. Furthermore, genome-centric metatranscriptomic analysis of ruminants with varying methane yields reveal elevated expression of select Methanomassiliicoccales genera in low methane animals and suggest that B12 exchanges could enable them to occupy ecological niches that possibly alter the direction of H2 utilization.

CONCLUSIONS

We provide a comprehensive and updated account of divergent Methanomassiliicoccales lineages, drawing from numerous uncultured genomes obtained from various habitats. We also highlight their unique metabolic capabilities involving B12, which could serve as promising targets for mitigating ruminant methane emissions by altering H2 flow.

Evolving understanding of rumen methanogen ecophysiology

Production of methane by methanogenic archaea, or methanogens, in the rumen of ruminants is a thermodynamic necessity for microbial conversion of feed to volatile fatty acids, which are essential nutrients for the animals. On the other hand, methane is a greenhouse gas and its production causes energy loss for the animal. Accordingly, there are ongoing efforts toward developing effective strategies for mitigating methane emissions from ruminant livestock that require a detailed understanding of the diversity and ecophysiology of rumen methanogens. Rumen methanogens evolved from free-living autotrophic ancestors through genome streamlining involving gene loss and acquisition. The process yielded an oligotrophic lifestyle, and metabolically efficient and ecologically adapted descendants. This specialization poses serious challenges to the efforts of obtaining axenic cultures of rumen methanogens, and consequently, the information on their physiological properties remains in most part inferred from those of their non-rumen representatives. This review presents the current knowledge of rumen methanogens and their metabolic contributions to enteric methane production. It also identifies the respective critical gaps that need to be filled for aiding the efforts to mitigate methane emission from livestock operations and at the same time increasing the productivity in this critical agriculture sector.

Effects of wheat bran replacement with pomegranate seed pulp on rumen fermentation, gas production, methanogen and protozoa populations of camel and goat rumen using competitive PCR technique: An in vitro study

BACKGROUND

Microbial populations in the rumen play an essential role in the degradation of Cellulosic dietary components and in providing nutrients to the host animal.

OBJECTIVE

This study aims to detect the effect of pomegranate seed pulp (PSP) on rumen fermentation, digestibility and methanogens and the protozoa population (by competitive polymerase chain reaction [PCR]) of the camel and goat rumen fluid.

MATERIALS AND METHODS

PSP was added to the experimental treatments and replaced by wheat bran (0%, 5% and 10%). Rumen fluid was collected from three goats and two camels according to the similarity of sex, breed, origin and time and used for three gas production studies. DNA extraction was performed by the RBB + c method, the ImageJ programme calculated band intensities (target and competing DNA), and line gradients were plotted based on the number of copies and intensity.

RESULTS

Our result showed that diets did not significantly affect the methanogen and protozoa population. Animal species affected microbial populations so that both populations in camels were less than goats. The production of gas and volatile fatty acids was not affected by diets. These two parameters and NH3 concentration and methane production in goats were higher than in camel. The pH of digested dry matter and microbial protein in camels was higher than in goats.

CONCLUSIONS

Therefore, the competitive PCR technique is an effective method for enumerating rumen microbiota. This supplementation can be considered a strategy to achieve performance and environmental benefits.

A Longitudinal Study with a Laser Methane Detector (LMD) Highlighting Lactation Cycle-Related Differences in Methane Emissions from Dairy Cows

Reversing climate change requires broad, cohesive, and strategic plans for the mitigation of greenhouse gas emissions from animal farming. The implementation and evaluation of such plans demand accurate and accessible methods for monitoring on-field CH4 concentration in eructating breath. Therefore, this paper describes a longitudinal study over six months, aiming to test a protocol using a laser methane detector (LMD) to monitor CH4 emissions in semi-extensive dairy farm systems. Over 10 time points, CH4 measurements were performed in dry (late gestation) and lactating cows at an Azorean dairy farm. Methane traits including CH4 concentration related to eructation (E_CH4) and respiration (R_CH4), and eructation events, were automatically computed from CH4 measured values using algorithms created for peak detection and analysis. Daily CH4 emission was estimated from each profile’s mean CH4 concentration (MEAN_CH4). Data were analyzed using a linear mixed model, including breed, lactation stage, and parity as fixed effects, and cow (subject) and time point as random effects. The results showed that Holsteins had higher E_CH4 than Jersey cows (p < 0.001). Although a breed-related trend was found in daily CH4 emission (p = 0.060), it was not significant when normalized to daily milk yield (p > 0.05). Methane emissions were lower in dry than in lactation cows (p < 0.05) and increased with the advancement of the lactation, even when normalizing it to daily milk yield (p < 0.05). Primiparous cows had lower daily CH4 emissions related to R_ CH4 compared to multiparous (p < 0.001). This allowed the identification of periods of higher CH4 emissions within the milk production cycle of dairy cows, and thus, the opportunity to tailor mitigation strategies accordingly.

Effect of Metabolizable Protein Supply on Milk Performance, Ruminal Fermentation, Apparent Total-Tract Digestibility, Energy and Nitrogen Utilization, and Enteric Methane Production of Ayrshire and Holstein Cows

In North America, the nutrient requirements of dairy cattle are predicted using the Cornell Net Carbohydrate and Protein System (CNCPS) or the National Research Council (NRC). As Holstein is the most predominant dairy cattle breed, these models were developed based on the phenotypic, physiological, and genetic characteristics of this breed. However, these models may not be appropriate to predict the nutrient requirements of other breeds, such as Ayrshire, that are phenotypically and genetically different from Holstein. The objective of this study was to evaluate the effects of increasing the metabolizable protein (MP) supply using CNCPS on milk performance, ruminal fermentation, apparent total-tract digestibility, energy and N utilization, and enteric methane production in Ayrshire vs. Holstein lactating dairy cows. Eighteen (nine Ayrshire; nine Holstein) lactating cows were used in a replicated 3 × 3 Latin square design (35-d periods) and fed diets formulated to meet 85%, 100%, or 115% of MP daily requirement. Except for milk production, no breed × MP supply interaction was observed for the response variables. Dry matter intake (DMI) and the yields of energy-corrected milk (ECM), fat, and protein were less (p < 0.01) in Ayrshire vs. Holstein cows. However, feed efficiency and N use efficiency for milk production did not differ between the two breeds, averaging 1.75 kg ECM/kg DMI and 33.7 g milk N/100 g N intake, respectively. Methane yield and intensity and urinary N also did not differ between the two breeds, averaging 18.8 g CH₄ /kg DMI, 10.8 g CH₄ /kg ECM, and 27.6 g N/100 g N intake, respectively. Yields of ECM and milk protein increased (p ≤ 0.01) with increasing MP supply from 85% to 100% but no or small increases occurred when MP supply increased from 100 to 115%. Feed efficiency increased linearly with an increasing MP supply. Nitrogen use efficiency (g N milk/100g N intake) decreased linearly (by up to 5.4 percentage units, (p < 0.01) whereas urinary N excretion (g/d or g/100 g N intake) increased linearly (p < 0.01) with an increasing MP supply. Methane yield and emission intensity were not affected by MP supply. This study shows that feed efficiency, N use efficiency, CH₄ (yield and intensity), and urinary N losses did not differ between Ayrshire and Holstein cows. Energy-corrected milk yield and feed efficiency increased, but N use efficiency decreased and urinary N losses increased with increasing dietary MP supply regardless of breed. Ayrshire and Holstein breeds responded similarly to increasing MP levels in the diet.

Characterization and Comparison of the Rumen Luminal and Epithelial Microbiome Profiles Using Metagenomic Sequencing Technique

Bacterial dysbiosis as a result of nutritional, bacterial, viral, and parasitic gastrointestinal infections can adversely affect the metabolism, productivity, and overall health of cattle. The purpose of this project was to characterize the commensal microbiota present in two locations of the rumen concomitantly in vivo with the animals undergoing habitual husbandry, as it was hypothesized that there are major differences in the commensal microbiota present in the two locations of the adult bovine major forestomach. A surgically fitted rumen cannula was used to allow ruminal lumen contents and mucosal biopsies to be collected from six crossbred yearling steers. In order to assess as much environmental and individual steer microbiota variation as possible, each animal was randomly sampled three times over a 3 week period. 16S rRNA sequencing was performed to provide a detailed descriptive analysis from phylum to genus taxonomic level. Significant differences were observed between luminal and epimural bacterial populations in the bovine rumen. As expected, a core microbiome composed by Firmicutes and Bacteroidetes represented over 90% of the microbiome, however, further analysis showed distinct diversity and distribution of the microbiome between the two locations. Characterizing the gastrointestinal microbiome in vivo is imperative. The novelty and the contribution of this study to the literature is the use of live cattle which allowed real-time sample collections and analysis of the rumen microbiome providing an understanding of what is normal in the live animal.

Vertebrate host phylogeny influences gut archaeal diversity

Commonly used 16S rRNA gene primers do not detect the full range of archaeal diversity present in the vertebrate gut. As a result, several questions regarding the archaeal component of the gut microbiota remain, including which Archaea are host-associated, the specificities of such associations and the major factors influencing archaeal diversity. Using 16S rRNA gene amplicon sequencing with primers that specifically target Archaea, we obtained sufficient sequence data from 185 gastrointestinal samples collected from 110 vertebrate species that span five taxonomic classes (Mammalia, Aves, Reptilia, Amphibia and Actinopterygii), of which the majority were wild. We provide evidence for previously undescribed Archaea-host associations, including Bathyarchaeia and Methanothermobacter, the latter of which was prevalent among Aves and relatively abundant in species with higher body temperatures, although this association could not be decoupled from host phylogeny. Host phylogeny explained archaeal diversity more strongly than diet, while specific taxa were associated with both factors, and cophylogeny was significant and strongest for mammalian herbivores. Methanobacteria was the only class predicted to be present in the last common ancestors of mammals and all host species. Further analysis indicated that Archaea-Bacteria interactions have a limited effect on archaeal diversity. These findings expand our current understanding of Archaea-vertebrate associations.

Holstein and Jersey Steers Differ in Rumen Microbiota and Enteric Methane Emissions Even Fed the Same Total Mixed Ration

Previous studies have focused on the rumen microbiome and enteric methane (CH₄) emissions in dairy cows, yet little is known about steers, especially steers of dairy breeds. In the present study, we comparatively examined the rumen microbiota, fermentation characteristics, and CH₄ emissions from six non-cannulated Holstein (710.33 ± 43.02 kg) and six Jersey (559.67 ± 32.72 kg) steers. The steers were fed the same total mixed ration (TMR) for 30 days. After 25 days of adaptation to the diet, CH₄ emissions were measured using GreenFeed for three consecutive days, and rumen fluid samples were collected on last day using stomach tubing before feeding (0 h) and 6 h after feeding. CH₄ production (g/d/animal), CH₄ yield (g/kg DMI), and CH₄ intensity (g/kg BW^0.75) were higher in the Jersey steers than in the Holstein steers. The lowest pH value was recorded at 6 h after feeding. The Jersey steers had lower rumen pH and a higher concentration of ammonia-nitrogen (NH₃-N). The Jersey steers had a numerically higher molar proportion of acetate than the Holstein steers, but the opposite was true for that of propionate. Metataxonomic analysis of the rumen microbiota showed that the two breeds had similar species richness, Shannon, and inverse Simpson diversity indexes. Principal coordinates analysis showed that the overall rumen microbiota was different between the two breeds. Both breeds were dominated by Prevotella ruminicola, and its highest relative abundance was observed 6 h after feeding. The genera Ethanoligenens, Succinivibrio, and the species Ethanoligenens harbinense, Succinivibrio dextrinosolvens, Prevotella micans, Prevotella copri, Prevotella oris, Prevotella baroniae, and Treponema succinifaciens were more abundant in Holstein steers while the genera Capnocytophaga, Lachnoclostridium, Barnesiella, Oscillibacter, Galbibacter, and the species Capnocytophaga cynodegmi, Galbibacter mesophilus, Barnesiella intestinihominis, Prevotella shahii, and Oscillibacter ruminantium in the Jersey steers. The Jersey steers were dominated by Methanobrevibacter millerae while the Holstein steers by Methanobrevibacter olleyae. The overall results suggest that sampling hour has little influence on the rumen microbiota; however, breeds of steers can affect the assemblage of the rumen microbiota and different mitigation strategies may be needed to effectively manipulate the rumen microbiota and mitigate enteric CH₄ emissions from these steers.

Influence of industry standard feeding frequencies on behavioral patterns and rumen and fecal bacterial communities in Holstein and Jersey cows

This study aimed to evaluate the effects of feeding frequency on behavioral patterns and on diurnal fermentation and bacteriome profiles of the rumen and feces in Holstein and Jersey cows. Ten Holstein and 10 Jersey cows were offered a TMR (53:47 forage-to-concentrate ratio dry matter basis) for ad libitum consumption and were randomly allocated within breed to one of the following feeding frequencies: (1) TMR delivered 1×/d (at 0600 h) or (2) TMR delivered 2×/d (at 0600 and 1800 h). The experiment lasted for 28 d with the first 14 d for cow adaptation to the Calan gates and the next 14 d for data collection. On d 23 and 24, an observer manually recorded the time budget (time spent lying, eating, drinking, standing, and milking), rumination activity, and number of visits to the feeding gate from each animal. On d 28, 5 concomitant collections of rumen and fecal samples were performed at intervals of 6 h via esophageal tubing and fecal grab, respectively. The bacteriome composition from these samples was determined through sequencing of the V4 region of the 16S rRNA gene. Feeding frequency did not affect behavioral patterns; however, Holstein cows spend more time lying (15.4 vs. 13.5 ± 0.8 h) and ruminating (401 vs. 331 ± 17.5 min) than Jersey cows. Fermentation profiles were similar by feeding frequency in both breeds. While no major diurnal fluctuations were observed in the fecal bacterial community from both breeds, diurnal fluctuations were identified in the rumen bacterial community from Holstein cows which appeared to follow pH responses. Overall, the bacterial community composition was not differentiated by industry standard feeding frequencies but was differentiated by breed and sample type.

Interactions between Anaerobic Fungi and Methanogens in the Rumen and Their Biotechnological Potential in Biogas Production from Lignocellulosic Materials

Anaerobic fungi in the digestive tract of herbivores are one of the critical types of fiber-degrading microorganisms present in the rumen. They degrade lignocellulosic materials using unique rhizoid structures and a diverse range of fiber-degrading enzymes, producing metabolic products such as H2/CO2, formate, lactate, acetate, and ethanol. Methanogens in the rumen utilize some of these products (e.g., H2 and formate) to produce methane. An investigation of the interactions between anaerobic fungi and methanogens is helpful as it provides valuable insight into the microbial interactions within the rumen. During the last few decades, research has demonstrated that anaerobic fungi stimulate the growth of methanogens and maintain methanogenic diversity. Meanwhile, methanogens increase the fiber-degrading capability of anaerobic fungi and stimulate metabolic pathways in the fungal hydrogenosome. The ability of co-cultures of anaerobic fungi and methanogens to degrade fiber and produce methane could potentially be a valuable method for the degradation of lignocellulosic materials and methane production.

The Rumen Specific Bacteriome in Dry Dairy Cows and Its Possible Relationship with Phenotypes

Most microbiome studies of dairy cows have investigated the compositions and functions of rumen microbial communities in lactating dairy cows. The importance of the relationships among hosts, microbiota, diet composition, and milk production remains unknown in dry dairy cows. Thus, in the present study, the composition of the rumen microbiome in cows from three dairy farms was investigated to identify core bacteria contributing to various physiological roles during rumen fermentation in dry dairy cows. The results indicated that ruminal fluid in dry dairy cows from different regional farms had core rumen microbiota that could be clearly distinguished from that of cows of the other farms. Further identification of key microorganisms associated with each farm revealed that Prevotella, Methanobrevibacter, Pseudobutyrivibrio, Ruminococcus, Bacteroides, and Streptococcus were major contributors. Spearman's correlation indicated that the abundance of genera such as Prevotella and Ruminococcus in dry dairy cows could indicate milk yield in the previous lactating period. Functional pathway analysis of the rumen bacterial communities demonstrated that amino acid metabolism and carbohydrate metabolism were the major pathways. Our findings provide knowledge of the composition and predicted functions of rumen microbiota in dry dairy cows from regional farms, which underscore the importance of the relationships among hosts, microbiota, diet composition, and milk production.

Impact of inulin and yeast containing synbiotic on calves' productivity and greenhouse gas production

Aim:

The research aimed to determine the impact of synbiotic: 6 g of prebiotic inulin and 5 g of probiotic Saccharomyces cerevisiae strain 1026 on calves’ productivity and greenhouse gas (GHG) production.

Materials and Methods:

The research was conducted with 10 Holstein Friesian and Red Holstein (Bos taurus L.) crossbreed calves of mean age 33±6 days and initial body weight 73.4±12.75 kg. We added the synbiotic into the diet of five dairy crossbreed calves (SynG) and five calves in control group (CoG) received non-supplemented diet. The duration of the experiment was 56 days. The weight of calves and amount of methane (CH₄) and carbon dioxide (CO₂) in the rumen were determined on day 1, 28, and 56. On day 56, three calves from each group were slaughtered. Meat samples were assessed for some indicators of meat quality. The main methanogens were detected in the rumen fluid and feces.

Results:

The weight gain during the whole experiment period of 56 days was higher in the SynG (62.6±13.75 kg) compared to CoG (36.8±7.98 kg) calves (p<0.01). There were no significant differences in the levels of protein (%), fat (unsaturated and saturated – %), and cholesterol (mg/100 g) in meat samples from both groups. At the end of the experiment, the amount of CH₄ in calves’ rumen in CoG was higher (Me=792.06 mg/m³, interquartile range [IQR] 755.06-873.59) compared to SynG (Me=675.41 mg/m³, IQR 653.46-700.50) group (p<0.01). The values for CO₂ were also increased in CoG (Me=4251.28 mg/m³, IQR 4045.58-4426.25) compared to SynG (Me=3266.06 mg/m³, IQR 1358.98-4584.91) group (p=0.001). There were no significant differences in the calves’ weight and certain methanogen species in rumen liquid and feces on the 56^th day of the experiment. Significantly higher results in the parameter total prokaryotes (V3) (bacteria+archaea) in rumen fluid were in SynG, whereas significantly higher results in the parameter total methanogens Met630/803 in rumen fluid were in CoG, p<0.05.

Conclusion:

The main results showed that the synbiotic can increase the daily weight gain in calves and decrease the amount of GHG in rumen but does not impact different methanogen species in rumen liquid and feces and meat protein, fat, and cholesterol levels.

Comparative diversity of microbiomes and Resistomes in beef feedlots, downstream environments and urban sewage influent

Background

Comparative knowledge of microbiomes and resistomes across environmental interfaces between animal production systems and urban settings is lacking. In this study, we executed a comparative analysis of the microbiota and resistomes of metagenomes from cattle feces, catch basin water, manured agricultural soil and urban sewage.

Results

Metagenomic DNA from composite fecal samples (FC; n = 12) collected from penned cattle at four feedlots in Alberta, Canada, along with water from adjacent catchment basins (CB; n = 13), soil (n = 4) from fields in the vicinity of one of the feedlots and urban sewage influent (SI; n = 6) from two municipalities were subjected to Illumina HiSeq2000 sequencing. Firmicutes exhibited the highest prevalence (40%) in FC, whereas Proteobacteria were most abundant in CB (64%), soil (60%) and SI (83%). Among sample types, SI had the highest diversity of antimicrobial resistance (AMR), and metal and biocide resistance (MBR) classes (13 & 15) followed by FC (10 & 8), CB (8 & 4), and soil (6 & 1). The highest antimicrobial resistant (AMR) gene (ARG) abundance was harboured by FC, whereas soil samples had a very small, but unique resistome which did not overlap with FC & CB resistomes. In the beef production system, tetracycline resistance predominated followed by macrolide resistance. The SI resistome harboured β-lactam, macrolide, tetracycline, aminoglycoside, fluoroquinolone and fosfomycin resistance determinants. Metal and biocide resistance accounted for 26% of the SI resistome with a predominance of mercury resistance.

Conclusions

This study demonstrates an increasing divergence in the nature of the microbiome and resistome as the distance from the feedlot increases. Consistent with antimicrobial use, tetracycline and macrolide resistance genes were predominant in the beef production system. One of the feedlots contributed both conventional (raised with antibiotics) and natural (raised without antibiotics) pens samples. Although natural pen samples exhibited a microbiota composition that was similar to samples from conventional pens, their resistome was less complex. Similarly, the SI resistome was indicative of drug classes used in humans and the greater abundance of mercury resistance may be associated with contamination of municipal water with household and industrial products.

Electronic supplementary material

The online version of this article (10.1186/s12866-019-1548-x) contains supplementary material, which is available to authorized users.

Rumen methanogen and protozoal communities of Tibetan sheep and Gansu Alpine Finewool sheep grazing on the Qinghai-Tibetan Plateau, China

Background

Tibetan sheep (TS) and Gansu Alpine Finewool sheep (GS) are both important plateau sheep raised and fed on the harsh Qinghai–Tibetan Plateau, China. Rumen methanogen and protozoal communities of plateau sheep are affected by their hosts and living environments, and play important roles in ruminant nutrition and greenhouse gas production. However, the characteristics, differences, and associations of these communities remain largely uncharacterized.

Results

The rumen methanogen and protozoal communities of plateau sheep were investigated by 16S/18S rRNA gene clone libraries. The predominant methanogen order in both sheep species was Methanobacteriales followed by Methanomassiliicoccales, which is consistent with those seen in global ruminants. However, the most dominant species was Methanobrevibacter millerae rather than Methanobrevibacter gottschalkii seen in most ruminants. Compared with GS and other ruminants, TS have more exclusive operational taxonomic units and a lower proportion (64.5%) of Methanobrevibacter. The protozoa were divided into Entodiniomorphida and Vestibuliferida, including nine genera and 15 species. The proportion of holotrich protozoa was much lower (1.1%) in TS than ordinary sheep. The most predominant genus was Entodinium (70.0%) in TS and Enoploplastron (48.8%) in GS, while the most common species was Entodinium furca monolobum (43.9%) and Enoploplastron triloricatum (45.0%) in TS and GS, respectively; Entodinium longinucleatum (22.8%) was only observed in TS. LIBSHUFF analysis indicated that the methanogen communities of TS were significantly different from those of GS, but no significant differences were found in protozoal communities.

Conclusion

Plateau sheep have coevolved with unique rumen methanogen and protozoal communities to adapt to harsh plateau environments. Moreover, the host appears to have a greater influence on rumen methanogen communities than on rumen protozoal communities. The observed associations of methanogens and protozoa, together with the findings of previous studies on methane emissions from ruminant livestock, revealed that the lower proportion of Methanobrevibacter and holotrich protozoa may be responsible for the lower methane emission of TS. These findings facilitate our understanding of the rumen microbial ecosystem in plateau sheep, and could help the development of new strategies to manipulate rumen microbes to improve productivity and reduce the emission of greenhouse gases.

Electronic supplementary material

The online version of this article (10.1186/s12866-018-1351-0) contains supplementary material, which is available to authorized users.

Lower dietary concentrate level increases bacterial diversity in the rumen of Cervus elaphus yarkandensis

The ruminal microbiota plays major roles in feed digestion. The composition and fermentation of the bacterial communities in 3 important ruminant species have been studied previously. Here, we extended this research to the effect of concentrate-to-forage ratios on ruminal bacterial communities in Tarim red deer (Cervus elaphus yarkandensis). Different concentrate-to-forage ratios (2:8, 3:7, 4:6, and 5:5) were fed to Tarim red deer for 20 days. Ruminal bacterial communities were elucidated by 16S ribosomal RNA gene sequencing on an Illumina HiSeq 2500 platform. The microbial composition and biodiversity at the different concentrate-to-forage ratio levels were analyzed using clustering of operational taxonomic units based on 97% sequence identity, taxonomic classification at the phylum and genus levels, α diversity, and β diversity. Rumen microorganisms of deer fed a diet with a concentrate-to-forage ratio of 2:8 had the highest species diversity, followed by ratios of 3:7, 4:6, and 5:5. The community structure of the A1 and A2 samples and the A3 and A4 samples was similar. The bacterial composition appeared to be affected by diet, with a lower dietary concentrate level tending to increase the richness and diversity of ruminal bacteria in the rumen of Tarim red deer.

Genetic Markers Are Associated with the Ruminal Microbiome and Metabolome in Grain and Sugar Challenged Dairy Heifers

Dairy heifers were subjected to a non-life-threatening challenge designed to induce ruminal acidosis by feeding grain and sugar. Large among animal variation in clinical signs of acidosis, rumen metabolite concentrations, and the rumen microbiome occurred. This exploratory study investigates sources of the variation by examining associations between the genome, metabolome, and microbiome, albeit with a limited population. The broader objective is to provide a rationale for a larger field study to identify markers for susceptibility to ruminal acidosis. Initially, heifers (n = 40) allocated to five feed additive groups were fed 20-days pre-challenge with a total mixed ration and additives. Fructose (0.1% of bodyweight/day) was added for the last 10 days pre-challenge. On day 21 heifers were challenged with 1.0% of bodyweight dry matter wheat + 0.2% of bodyweight fructose + additives. Rumen samples were collected via stomach tube weekly (day 0, 7, and 14) and at five times over 3.6 h after challenge and analyzed for pH and volatile fatty acid, ammonia, D-, and L-lactate concentrations. Relative abundance of bacteria and archaea were determined using Illumina MiSeq. Genotyping was undertaken using a 150K Illumina SNPchip. Genome-wide association was performed for metabolite and microbiome measures (n = 33). Few genome associations occurred with rumen pH, concentration of acetate, propionate, total volatile fatty acids, or ammonia, or the relative abundance of the Firmicutes, Bacteroidetes, and Spirochaetes phyla. Metabolites and microbial phyla that had markers associated and quantitative trait loci (QTL) were: acetate to propionate ratio (A:P), D-, L-, and total lactate, butyrate, acidosis eigenvalue, Actinobacteria, Chloroflexi, Euryarchaeota, Fibrobacteres, Planctomycetes, Proteobacteria, and Tenericutes. A putative genomic region overlapped for Actinobacteria, Euryarchaeota, and Fibrobacteres and covered the region that codes for matrix extracellular phosphoglycoprotein (MEPE). Other overlapping regions were: (1) Chloroflexi, Tenericutes, and A:P, (2) L- and total lactate and Actinobacteria, and (3) Actinobacteria, Euryarchaeota, Fibrobacteres, and A:P. Genome-wide associations with the metabolome and microbiome occurred despite the small population, suggesting that markers for ruminal acidosis susceptibility exist. The findings may explain some of the variation in metabolomic and microbial data and provide a rationale for a larger study with a population that has variation in acidosis.

Comparison of Rumen and Manure Microbiomes and Implications for the Inoculation of Anaerobic Digesters

Cattle manure is frequently used as an inoculum for the start-up of agricultural biogas plants or as a co-substrate in the anaerobic digestion of lignocellulosic feedstock. Ruminal microbiota are considered to be effective plant fiber degraders, but the microbes contained in manure do not necessarily reflect the rumen microbiome. The aim of this study was to compare the microbial community composition of cow rumen and manure with respect to plant fiber-digesting microbes. Bacterial and methanogenic communities of rumen and manure samples were examined by 454 amplicon sequencing of bacterial 16S rRNA genes and mcrA genes, respectively. Rumen fluid samples were dominated by Prevotellaceae (29%), whereas Ruminococcaceae was the most abundant family in the manure samples (31%). Fibrobacteraceae (12%) and Bacteroidaceae (13%) were the second most abundant families in rumen fluid and manure, respectively. The high abundances of fiber-degrading bacteria belonging to Prevotellaceae and Fibrobacteraceae might explain the better performance of anaerobic digesters inoculated with rumen fluid. Members of the genus Methanobrevibacter were the predominant methanogens in the rumen fluid, whereas methanogenic communities of the manure samples were dominated by the candidate genus Methanoplasma. Our results suggest that inoculation or bioaugmentation with fiber-digesting rumen microbiota can enhance the anaerobic digestion of lignocellulosic biomass.

Compositional mixed modeling of methane emissions and ruminal volatile fatty acids from individual cattle and multiple experiments

The aim of the study was to investigate the association of methane (CH) yields (g/kg DMI) with rumen VFA molar proportions and animal and diet-related covariates from individual animals and multiple experiments. The dataset available consisted of 284 measurements of CH yields for beef cattle from 6 experiments measured in indirect respiration chambers. A compositional modeling approach was employed where VFA measurements were considered as a whole, instead of in isolation, emphasizing their multivariate relative scale. The analysis revealed expected close groupings of acetate and butyrate; propionate and valerate; iso-butyrate and iso-valerate. Linear mixed models were then fitted to examine relationships between CH yield and VFA, represented by meaningful log-contrasts of components called compositional balances, while accounting for other animal and diet-related covariates and random variability between experiments. A compositional balance representing (acetate × butyrate)/propionate best explained the contribution of VFA to variation in CH yield. The covariates DMI, forage:concentrate proportion (expressed as a categorical variable diet type: high concentrate, mixed forage:concentrate or high forage), and diet ME were also statistically significant. These results provided new insights into the relative inter-relationships among VFA measurements and also between VFA and CH yield. In conclusion, VFA molar proportions as represented by compositional balances were a significant contributor to explaining variation in CH yields from individual cattle.

Taxon abundance, diversity, co-occurrence and network analysis of the ruminal microbiota in response to dietary changes in dairy cows

The ruminal microbiome, comprising large numbers of bacteria, ciliate protozoa, archaea and fungi, responds to diet and dietary additives in a complex way. The aim of this study was to investigate the benefits of increasing the depth of the community analysis in describing and explaining responses to dietary changes. Quantitative PCR, ssu rRNA amplicon based taxa composition, diversity and co-occurrence network analyses were applied to ruminal digesta samples obtained from four multiparous Nordic Red dairy cows fitted with rumen cannulae. The cows received diets with forage:concentrate ratio either 35:65 (diet H) or 65:35 (L), supplemented or not with sunflower oil (SO) (0 or 50 g/kg diet dry matter), supplied in a 4 × 4 Latin square design with a 2 × 2 factorial arrangement of treatments and four 35-day periods. Digesta samples were collected on days 22 and 24 and combined. QPCR provided a broad picture in which a large fall in the abundance of fungi was seen with SO in the H but not the L diet. Amplicon sequencing showed higher community diversity indices in L as compared to H diets and revealed diet specific taxa abundance changes, highlighting large differences in protozoal and fungal composition. Methanobrevibacter ruminantium and Mbb. gottschalkii dominated archaeal communities, and their abundance correlated negatively with each other. Co-occurrence network analysis provided evidence that no microbial domain played a more central role in network formation, that some minor-abundance taxa were at nodes of highest centrality, and that microbial interactions were diet specific. Networks added new dimensions to our understanding of the diet effect on rumen microbial community interactions.

Methane Production in Dairy Cows Correlates with Rumen Methanogenic and Bacterial Community Structure

Methane (CH₄) is produced as an end product from feed fermentation in the rumen. Yield of CH₄ varies between individuals despite identical feeding conditions. To get a better understanding of factors behind the individual variation, 73 dairy cows given the same feed but differing in CH₄ emissions were investigated with focus on fiber digestion, fermentation end products and bacterial and archaeal composition. In total 21 cows (12 Holstein, 9 Swedish Red) identified as persistent low, medium or high CH₄ emitters over a 3 month period were furthermore chosen for analysis of microbial community structure in rumen fluid. This was assessed by sequencing the V4 region of 16S rRNA gene and by quantitative qPCR of targeted Methanobrevibacter groups. The results showed a positive correlation between low CH₄ emitters and higher abundance of Methanobrevibacter ruminantium clade. Principal coordinate analysis (PCoA) on operational taxonomic unit (OTU) level of bacteria showed two distinct clusters (P < 0.01) that were related to CH₄ production. One cluster was associated with low CH₄ production (referred to as cluster L) whereas the other cluster was associated with high CH₄ production (cluster H) and the medium emitters occurred in both clusters. The differences between clusters were primarily linked to differential abundances of certain OTUs belonging to Prevotella. Moreover, several OTUs belonging to the family Succinivibrionaceae were dominant in samples belonging to cluster L. Fermentation pattern of volatile fatty acids showed that proportion of propionate was higher in cluster L, while proportion of butyrate was higher in cluster H. No difference was found in milk production or organic matter digestibility between cows. Cows in cluster L had lower CH₄/kg energy corrected milk (ECM) compared to cows in cluster H, 8.3 compared to 9.7 g CH₄/kg ECM, showing that low CH₄ cows utilized the feed more efficient for milk production which might indicate a more efficient microbial population or host genetic differences that is reflected in bacterial and archaeal (or methanogens) populations.

The ruminal microbiome associated with methane emissions from ruminant livestock

Methane emissions from ruminant livestock contribute significantly to the large environmental footprint of agriculture. The rumen is the principal source of methane, and certain features of the microbiome are associated with low/high methane phenotypes. Despite their primary role in methanogenesis, the abundance of archaea has only a weak correlation with methane emissions from individual animals. The composition of the archaeal community appears to have a stronger effect, with animals harbouring the Methanobrevibacter gottschalkii clade tending to be associated with greater methane emissions. Ciliate protozoa produce abundant H₂, the main substrate for methanogenesis in the rumen, and their removal (defaunation) results in an average 11% lower methane emissions in vivo, but the results are not consistent. Different protozoal genera seem to result in greater methane emissions, though community types (A, AB, B and O) did not differ. Within the bacteria, three different ‘ruminotypes’ have been identified, two of which predispose animals to have lower methane emissions. The two low-methane ruminotypes are generally characterized by less abundant H₂-producing bacteria. A lower abundance of Proteobacteria and differences in certain Bacteroidetes and anaerobic fungi seem to be associated with high methane emissions. Rumen anaerobic fungi produce abundant H₂ and formate, and their abundance generally corresponds to the level of methane emissions. Thus, microbiome analysis is consistent with known pathways for H₂ production and methanogenesis, but not yet in a predictive manner. The production and utilisation of formate by the ruminal microbiota is poorly understood and may be a source of variability between animals.

Dietary pea fiber increases diversity of colonic methanogens of pigs with a shift from Methanobrevibacter to Methanomassiliicoccus-like genus and change in numbers of three hydrogenotrophs

BACKGROUND

Pea fiber (PF) is a potential fibrous supplement in swine production. The influence of dietary PF on microbial community in the colon of pigs remains largely unexplored. Methanogens in the hindgut of monogastric animals play important roles in degradation of dietary fibers and efficient removal of microbial metabolic end product H2. Understanding the impact of dietary PF on the structure of colonic methanogens may help understand the mechanisms of microbe-mediated physiological functions of PF. This study investigated the influence of PF on the diversity and quantity and/or activity of colonic methanongens of piglets and finishing pigs. Four archaeal 16S rRNA clone libraries were constructed for piglets and finishers fed with control (Piglet-C and Finisher-C) or PF diet (Piglet-P and Finisher-P).

RESULTS

There were 195, 190, 194 and 196 clones obtained from the library Piglet-C, Piglet-P, Finisher-C and Finisher-P, respectively, with corresponding 12, 11, 11 and 16 OTUs (operational taxonomic units). Significant differences of Shannon Index among the four libraries were found (P < 0.05). Libshuff analysis showed that the archaeal community structure among the four libraries were significantly different (P < 0.0001). The predominant methanogens shifted from Methanobrevibacter to Methanobrevibacter and Methanomassiliicoccus-like genus as a result of dietary PF. Supplementation of PF significantly increased the copy numbers of mcrA and dsrA genes (P < 0.05).

CONCLUSIONS

Alteration of methanogenic community structure may lead to functional transition from utilization of H2/CO2 to employment of both H2/CO2 and methanol/CO2. Quantification of three functional genes (mcrA, dsrA and fhs) of methanogens, sulfate-reducing bacteria (SRB) and acetogens revealed that dietary PF also increased the activity of methanogens and SRB,probably associated with increased proportion of Methanomassiliicoccus luminyensis-species. Further study is required to examine the interaction between specific methanogens and SRB during fermentation of dietary PF.

Rumen and Cecum Microbiomes in Reindeer (Rangifer tarandus tarandus) Are Changed in Response to a Lichen Diet and May Affect Enteric Methane Emissions

Reindeer (Rangifer tarandus tarandus) are large Holarctic herbivores whose heterogeneous diet has led to the development of a unique gastrointestinal microbiota, essential for the digestion of arctic flora, which may include a large proportion of lichens during winter. Lichens are rich in plant secondary metabolites, which may affect members of the gut microbial consortium, such as the methane-producing methanogenic archaea. Little is known about the effect of lichen consumption on the rumen and cecum microbiotas and how this may affect methanogenesis in reindeer. Here, we examined the effects of dietary lichens on the reindeer gut microbiota, especially methanogens. Samples from the rumen and cecum were collected from two groups of reindeer, fed either lichens (Ld: n = 4), or a standard pelleted feed (Pd: n = 3). Microbial densities (methanogens, bacteria and protozoa) were quantified using quantitative real-time PCR and methanogen and bacterial diversities were determined by 454 pyrosequencing of the 16S rRNA genes.

In general, the density of methanogens were not significantly affected (p>0.05) by the intake of lichens. Methanobrevibacter constituted the main archaeal genus (>95% of reads), with Mbr. thaueri CW as the dominant species in both groups of reindeer. Bacteria belonging to the uncharacterized Ruminococcaceae and the genus Prevotella were the dominant phylotypes in the rumen and cecum, in both diets (ranging between 16–38% total sequences). Bacteria belonging to the genus Ruminococcus (3.5% to 0.6%; p = 0.001) and uncharacterized phylotypes within the order Bacteroidales (8.4% to 1.3%; p = 0.027), were significantly decreased in the rumen of lichen-fed reindeer, but not in the cecum (p = 0.2 and p = 0.087, respectively). UniFrac-based analyses showed archaeal and bacterial libraries were significantly different between diets, in both the cecum and the rumen (vegan::Adonis: pseudo-F<0.05). Based upon previous literature, we suggest that the altered methanogen and bacterial profiles may account for expected lower methane emissions from lichen-fed reindeer.

Influence of periparturient and postpartum diets on rumen methanogen communities in three breeds of primiparous dairy cows

BACKGROUND

Enteric methane from rumen methanogens is responsible for 25.9 % of total methane emissions in the United States. Rumen methanogens also contribute to decreased animal feed efficiency. For methane mitigation strategies to be successful, it is important to establish which factors influence the rumen methanogen community and rumen volatile fatty acids (VFA). In the present study, we used next-generation sequencing to determine if dairy breed and/or days in milk (DIM) (high-fiber periparturient versus high-starch postpartum diets) affect the rumen environment and methanogen community of primiparous Holstein, Jersey, and Holstein-Jersey crossbreeds.

RESULTS

When the 16S rRNA gene sequences were processed and assigned to operational taxonomic units (OTU), a core methanogen community was identified, consisting of Methanobrevibacter (Mbr.) smithii, Mbr. thaueri, Mbr. ruminantium, and Mbr. millerae. The 16S rRNA gene sequence reads clustered at 3 DIM, but not by breed. At 3 DIM, the mean % abundance of Mbr. thaueri was lower in Jerseys (26.9 %) and higher in Holsteins (30.7 %) and Holstein-Jersey crossbreeds (30.3 %) (P < 0.001). The molar concentrations of total VFA were higher at 3 DIM than at 93, 183, and 273 DIM, whereas the molar proportions of propionate were increased at 3 and 93 DIM, relative to 183 and 273 DIM. Rumen methanogen densities, distributions of the Mbr. species, and VFA molar proportions did not differ by breed.

CONCLUSIONS

The data from the present study suggest that a core methanogen community is present among dairy breeds, through out a lactation. Furthermore, the methanogen communities were more influenced by DIM and the breed by DIM interactions than breed differences.

Breed and Lactation Stage Alter the Rumen Protozoal Fatty Acid Profiles and Community Structures in Primiparous Dairy Cattle

The protozoal fatty acid (FA) composition and community structure are important to dairy cattle nutrition and their products. The purpose of the study was to observe if the rumen protozoal FA profiles and protozoal community structure differed by breed and lactation stage. At 93, 183, and 273 days in milk (DIM), whole rumen digesta samples were collected from seven co-housed Holstein (H), eight Jersey (J), and seven Holstein-Jersey crossbreed (C) cows. Rumen protozoal linoleic acid was higher at 183 DIM (8.1%) and 273 DIM (8.3%) than at 93 DIM (5.7%). Oleic acid was the most abundant protozoal unsaturated FA (10.1%). Protozoal rumenic acid and protozoa of the genus Metadinium were higher in J (9.9%) than in H (0.52%) and C (0.96%). Protozoa belonging to the genus Entodinium were more abundant in H (45.2%) than in J (23.4%) and C (30.2%). In conclusion, breed and DIM affected several protozoal FAs and genera.

Nitrate and Inhibition of Ruminal Methanogenesis: Microbial Ecology, Obstacles, and Opportunities for Lowering Methane Emissions from Ruminant Livestock

Ruminal methane production is among the main targets for greenhouse gas (GHG) mitigation for the animal agriculture industry. Many compounds have been evaluated for their efficacy to suppress enteric methane production by ruminal microorganisms. Of these, nitrate as an alternative hydrogen sink has been among the most promising, but it suffers from variability in efficacy for reasons that are not understood. The accumulation of nitrite, which is poisonous when absorbed into the animal’s circulation, is also variable and poorly understood. This review identifies large gaps in our knowledge of rumen microbial ecology that handicap the further development and safety of nitrate as a dietary additive. Three main bacterial species have been associated historically with ruminal nitrate reduction, namely Wolinella succinogenes, Veillonella parvula, and Selenomonas ruminantium, but others almost certainly exist in the largely uncultivated ruminal microbiota. Indications are strong that ciliate protozoa can reduce nitrate, but the significance of their role relative to bacteria is not known. The metabolic fate of the reduced nitrate has not been studied in detail. It is important to be sure that nitrate metabolism and efforts to enhance rates of nitrite reduction do not lead to the evolution of the much more potent GHG, nitrous oxide. The relative importance of direct inhibition of archaeal methanogenic enzymes by nitrite or the efficiency of capture of hydrogen by nitrate reduction in lowering methane production is also not known, nor are nitrite effects on other members of the microbiota. How effective would combining mitigation methods be, based on our understanding of the effects of nitrate and nitrite on the microbiome? Answering these fundamental microbiological questions is essential in assessing the potential of dietary nitrate to limit methane emissions from ruminant livestock.

Potential of tannin-rich plants for modulating ruminal microbes and ruminal fermentation in sheep

The objective of this work was to study nutritional strategies for decreasing methane production by ruminants fed tropical diets, combining in vitro and in vivo methods. The in vitro approach was used to evaluate the dose effect of condensed tannins (CT) contained in leaves of Gliricidia sepium, Leucaena leucocephala, and Manihot esculenta (39, 75, and 92 g CT/kg DM, respectively) on methane production and ruminal fermentation characteristics. Tannin-rich plants (TRP) were incubated for 24 h alone or mixed with a natural grassland hay based on Dichanthium spp. (control plant), so that proportions of TRP were 0, 0.25, 0.5, 0.75, and 1.0. Methane production, VFA concentration, and fermented OM decreased with increased proportions of TRP. Numerical differences on methane production and VFA concentration among TRP sources may be due to differences in their CT content, with greater effects for L. leucocephala and M. esculenta than for G. sepium. Independently of TRP, the response to increasing doses of CT was linear for methane production but quadratic for VFA concentration. As a result, at moderate tannin dose, methane decreased more than VFA. The in vivo trial was conducted to investigate the effect of TRP on different ruminal microbial populations. To this end, 8 rumen-cannulated sheep from 2 breeds (Texel and Blackbelly) were used in two 4 × 4 Latin square designs. Diets were fed ad libitum and were composed of the same feeds used for the in vitro trial: control plant alone or combined with pellets made from TRP leaves at 44% of the diet DM. Compared to TRP, concentration of Ruminococcus flavefaciens was greater for the control diet and concentration of Ruminococcus albus was least for the control diet. The methanogen population was greater for Texel than for Blackbelly. By contrast, TRP-containing diets did not affect protozoa or Fibrobacter succinogenes numbers. Hence, TRP showed potential for mitigating methane production by ruminants. These findings suggest that TRP fed as pellets could be used to decrease methane production.

Associative patterns among anaerobic fungi, methanogenic archaea, and bacterial communities in response to changes in diet and age in the rumen of dairy cows

The rumen microbiome represents a complex microbial genetic web where bacteria, anaerobic rumen fungi (ARF), protozoa and archaea work in harmony contributing to the health and productivity of ruminants. We hypothesized that the rumen microbiome shifts as the dairy cow advances in lactations and these microbial changes may contribute to differences in productivity between primiparous (first lactation) and multiparous (≥second lactation) cows. To this end, we investigated shifts in the ruminal ARF and methanogenic communities in both primiparous (n = 5) and multiparous (n = 5) cows as they transitioned from a high forage to a high grain diet upon initiation of lactation. A total of 20 rumen samples were extracted for genomic DNA, amplified using archaeal and fungal specific primers, sequenced on a 454 platform and analyzed using QIIME. Community comparisons (Bray-Curtis index) revealed the effect of diet (P < 0.01) on ARF composition, while archaeal communities differed between primiparous and multiparous cows (P < 0.05). Among ARF, several lineages were unclassified, however, phylum Neocallimastigomycota showed the presence of three known genera. Abundance of Cyllamyces and Caecomyces shifted with diet, whereas Orpinomyces was influenced by both diet and age. Methanobrevibacter constituted the most dominant archaeal genus across all samples. Co-occurrence analysis incorporating taxa from bacteria, ARF and archaea revealed syntrophic interactions both within and between microbial domains in response to change in diet as well as age of dairy cows. Notably, these interactions were numerous and complex in multiparous cows, supporting our hypothesis that the rumen microbiome also matures with age to sustain the growing metabolic needs of the host. This study provides a broader picture of the ARF and methanogenic populations in the rumen of dairy cows and their co-occurrence implicates specific relationships between different microbial domains in response to diet and age.

Toward the identification of methanogenic archaeal groups as targets of methane mitigation in livestock animalsr

In herbivores, enteric methane is a by-product from the digestion of plant biomass by mutualistic gastrointestinal tract (GIT) microbial communities. Methane is a potent greenhouse gas that is not assimilated by the host and is released into the environment where it contributes to climate change. Since enteric methane is exclusively produced by methanogenic archaea, the investigation of mutualistic methanogen communities in the GIT of herbivores has been the subject of ongoing research by a number of research groups. In an effort to uncover trends that would facilitate the development of efficient methane mitigation strategies for livestock species, we have in this review summarized and compared currently available results from published studies on this subject. We also offer our perspectives on the importance of pursuing current research efforts on the sequencing of gut methanogen genomes, as well as investigating their cellular physiology and interactions with other GIT microorganisms.

Response of the Rumen Microbiota of Sika Deer (Cervus nippon) Fed Different Concentrations of Tannin Rich Plants

High throughput sequencing was used to examine the rumen microbiota of sika deer fed high (OLH) and low concentration (OLL) of tannin rich oak leaves. The results showed that Prevotella spp. were the most dominant bacteria. The most predominant methanogens were the members of the order Methanoplasmatales. The dominant rumen protozoa were Entodinium longinucleatum, Eudiplodinium maggii, and Epidinium caudatum, and the fungal communities were mostly represented by Piromyces spp. Moreover, the relative abundance of Pseudobutyrivibrio spp. (P=0.026), unidentified bacteria (P=0.028), and Prevotella spp. (P=0.022) was lower in the OLH group than in the OLL group. The concentration of propionate in the OLH group was greater than in the OLL group (P=0.006). Patterns of relationships showed that methanogens belonging to the order Methanoplasmatales were negatively correlated with Treponema spp., Ent. Longinucleatum, and acetate. Methanosphaera stadtmanae was positively correlated to propionate, while Methanobrevibacter ruminantium was negatively associated with Methanobrevibacter thaueri and Methanobrevibacter millerae. Tannins altered the rumen microbes and fermentation patterns. However, the response of the entire rumen microbiota and the relationship between rumen microorganisms and the fermentation parameters were not fully understood.

Methyl-coenzyme M reductase A as an indicator to estimate methane production from dairy cows

The evaluation of greenhouse gas mitigation strategies requires the quantitative assessment of individual methane production. Because methane measurement in respiration chambers is highly accurate, but also comprises various disadvantages such as limited capacity and high costs, the establishment of an indicator for estimating methane production of individual ruminants would provide an alternative to direct methane measurement. Methyl-coenzyme M reductase is involved in methanogenesis and the subunit α of methyl-coenzyme M reductase is encoded by the mcrA gene of rumen archaea. We therefore examined the relationship between methane emissions of Holstein dairy cows measured in respiration chambers with 2 different diets (high- and medium-concentrate diet) and the mcrA DNA and mcrA cDNA abundance determined from corresponding rumen fluid samples. Whole-body methane production per kilogram of dry matter intake and mcrA DNA normalized to the abundance of the rrs gene coding for 16S rRNA correlated significantly when using qmcrA primers. Use of qmcrA primers also revealed linear correlation between mcrA DNA copy number and methane yield. Regression analyses based on normalized mcrA cDNA abundances revealed no significant linear correlation with methane production per kilogram of dry matter intake. Furthermore, the correlations between normalized mcrA DNA abundance and the rumen fluid concentration of acetic and isobutyric acid were positive, whereas the correlations with propionic and lactic acid were negative. These data suggest that the mcrA DNA approach based on qmcrA primers could potentially be a molecular proxy for methane yield after further refinement.

Examination of the rumen bacteria and methanogenic archaea of wild impalas (Aepyceros melampus melampus) from Pongola, South Africa

Although the rumen microbiome of domesticated ruminants has been evaluated, few studies have explored the rumen microbiome of wild ruminants, and no studies have identified the rumen microbiome in the impala (Aepyceros melampus melampus). In the present study, next-generation sequencing and real-time polymerase chain reaction were used to investigate the diversity and density of the bacteria and methanogenic archaea residing in the rumen of five adult male impalas, culled during the winter dry season in Pongola, South Africa. A total of 15,323 bacterial 16S rRNA gene sequences (from five impala), representing 3,892 different phylotypes, were assigned to 1,902 operational taxonomic units (OTUs). A total of 20,124 methanogen 16S rRNA gene sequence reads (from four impala), of which 5,028 were unique, were assigned to 344 OTUs. From the total sequence reads, Bacteroidetes, Proteobacteria, and Firmicutes were the most abundant bacterial phyla. While the majority of the bacterial genera found were unclassified, Prevotella and Cupriavidus were the most abundant classified genera. For methanogens, the genera Methanobrevibacter and Methanosphaera represented 94.3% and 4.0% of the classified sequences, respectively. Most notable was the identification of Methanobrevibacter thaueri-like 16S rRNA gene sequence reads in all four impala samples, representing greater than 30% of each individual's total sequences. Both data sets are accessible through NCBI's Sequence Read Archive (SRA), under study accession number SRP [048619]. The densities of bacteria (1.26 × 10(10)-3.82 × 10(10) cells/ml whole rumen contents) and methanogens (4.48 × 10(8)-7.2 × 10(9) cells/ml of whole rumen contents) from five individual impala were similar to those typically observed in domesticated ruminants.

Bacteria and methanogens differ along the gastrointestinal tract of Chinese roe deer (Capreolus pygargus)

The current study provides the insight into the bacteria in the gastrointestinal tract (GIT) and methanogens presented in the rumen and cecum of the Chinese roe deer (Capreolus pygargus). The ruminal, ileal, cecal, and colonic contents, as well as feces, were obtained from each of the three, free-range, roe deer ingesting natural pasture after euthanasia. For the bacterial community, a total of 697,031 high-quality 16S rRNA gene sequences were generated using high-throughput sequencing, and assigned to 2,223 core operational taxonomic units (OTUs) (12 bacterial phyla and 87 genera). The phyla Firmicutes (51.2%) and Bacteroidetes (39.4%) were the dominant bacteria in the GIT of roe deer. However, the bacterial community in the rumen was significantly (P<0.01) different from the other sampled regions along the GIT. Secondly, Prevotella spp., Anaerovibrio spp., and unidentified bacteria within the families Veillonellaceae and Paraprevotellaceae were more abundant in the rumen than in the other regions. Unidentified bacteria within the family Enterobacteriaceae, Succinivibrio spp., and Desulfovibrio spp. were more predominant in the colon than in other regions. Unidentified bacteria within the family Ruminococcaceae, and Bacteroides spp. were more prevalent in the ileum, cecum and fecal pellets. For methanogens in the rumen and cecum, a total of 375,647 high quality 16S rRNA gene sequences were obtained and assigned to 113 core OTUs. Methanobrevibacter millerae was the dominant species accounting for 77.3±7.4 (S.E) % and 68.9±4.4 (S.E) % of total sequences in the rumen and cecum of roe deer, respectively. However, the abundance of Methanobrevibacter smithii was higher in the rumen than in the cecum (P = 0.004). These results revealed that there was intra variation in the bacterial community composition across the GIT of roe deer, and also showed that the methanogen community in the rumen differed from that in the cecum.

Few highly abundant operational taxonomic units dominate within rumen methanogenic archaeal species in New Zealand sheep and cattle

Sequencing and analyses of 16S rRNA gene amplicons were performed to estimate the composition of the rumen methanogen community in 252 samples from eight cohorts of sheep and cattle, separated into 16 different sample groups by diet, and to determine which methanogens are most prominent in the rumens of farmed New Zealand ruminants. Methanobacteriales (relative abundance ± standard deviation, 89.6% ± 9.8%) and Methanomassiliicoccales (10.4% ± 9.8%) were the two major orders and contributed 99.98% (±0.1%) to the rumen methanogen communities in the samples. Sequences from Methanobacteriales were almost entirely from only four different species (or clades of very closely related species). Each was detectable in at least 89% of the samples. These four species or clades were the Methanobrevibacter gottschalkii clade and Methanobrevibacter ruminantium clade with a mean abundance of 42.4% (±19.5% standard deviation) and 32.9% (±18.8%), respectively, and Methanosphaera sp. ISO3-F5 (8.2% ± 6.7%) and Methanosphaera sp. group5 (5.6% ± 5.7%). These four species or clades appeared to be primarily represented by only one or, in one case, two dominant sequence types per species or clade when the sequences were grouped into operational taxonomic units (OTUs) at 99% sequence identity. The mean relative abundance of Methanomassiliicoccales in the samples was relatively low but exceeded 40% in some of the treatment groups. Animal feed affected the apparent methanogen community structure of both orders, as evident from differences in relative abundances of the major OTUs in animals under different feeding regimens.

Diversity and community composition of methanogenic archaea in the rumen of Scottish upland sheep assessed by different methods

Ruminal archaeomes of two mature sheep grazing in the Scottish uplands were analysed by different sequencing and analysis methods in order to compare the apparent archaeal communities. All methods revealed that the majority of methanogens belonged to the Methanobacteriales order containing the Methanobrevibacter, Methanosphaera and Methanobacteria genera. Sanger sequenced 1.3 kb 16S rRNA gene amplicons identified the main species of Methanobrevibacter present to be a SGMT Clade member Mbb. millerae (≥91% of OTUs); Methanosphaera comprised the remainder of the OTUs. The primers did not amplify ruminal Thermoplasmatales-related 16S rRNA genes. Illumina sequenced V6–V8 16S rRNA gene amplicons identified similar Methanobrevibacter spp. and Methanosphaera clades and also identified the Thermoplasmatales-related order as 13% of total archaea. Unusually, both methods concluded that Mbb. ruminantium and relatives from the same clade (RO) were almost absent. Sequences mapping to rumen 16S rRNA and mcrA gene references were extracted from Illumina metagenome data. Mapping of the metagenome data to16S rRNA gene references produced taxonomic identification to Order level including 2–3% Thermoplasmatales, but was unable to discriminate to species level. Mapping of the metagenome data to mcrA gene references resolved 69% to unclassified Methanobacteriales. Only 30% of sequences were assigned to species level clades: of the sequences assigned to Methanobrevibacter, most mapped to SGMT (16%) and RO (10%) clades. The Sanger 16S amplicon and Illumina metagenome mcrA analyses showed similar species richness (Chao1 Index 19–35), while Illumina metagenome and amplicon 16S rRNA analysis gave lower richness estimates (10–18). The values of the Shannon Index were low in all methods, indicating low richness and uneven species distribution. Thus, although much information may be extracted from the other methods, Illumina amplicon sequencing of the V6–V8 16S rRNA gene would be the method of choice for studying rumen archaeal communities.

Quantitative analysis of ruminal methanogenic microbial populations in beef cattle divergent in phenotypic residual feed intake (RFI) offered contrasting diets

Background

Methane (CH₄) emissions in cattle are an undesirable end product of rumen methanogenic fermentative activity as they are associated not only with negative environmental impacts but also with reduced host feed efficiency. The aim of this study was to quantify total and specific rumen microbial methanogenic populations in beef cattle divergently selected for residual feed intake (RFI) while offered (i) a low energy high forage (HF) diet followed by (ii) a high energy low forage (LF) diet. Ruminal fluid was collected from 14 high (H) and 14 low (L) RFI animals across both dietary periods. Quantitative real time PCR (qRT-PCR) analysis was conducted to quantify the abundance of total and specific rumen methanogenic microbes. Spearman correlation analysis was used to investigate the association between the relative abundance of methanogens and animal performance, rumen fermentation variables and diet digestibility.

Results

Abundance of methanogens, did not differ between RFI phenotypes. However, relative abundance of total and specific methanogen species was affected (P < 0.05) by diet type, with greater abundance observed while animals were offered the LF compared to the HF diet.

Conclusions

These findings suggest that differences in abundance of specific rumen methanogen species may not contribute to variation in CH₄ emissions between efficient and inefficient animals, however dietary manipulation can influence the abundance of total and specific methanogen species.

RIM-DB: a taxonomic framework for community structure analysis of methanogenic archaea from the rumen and other intestinal environments

Methane is formed by methanogenic archaea in the rumen as one of the end products of feed fermentation in the ruminant digestive tract. To develop strategies to mitigate anthropogenic methane emissions due to ruminant farming, and to understand rumen microbial differences in animal feed conversion efficiency, it is essential that methanogens can be identified and taxonomically classified with high accuracy. Currently available taxonomic frameworks offer only limited resolution beyond the genus level for taxonomic assignments of sequence data stemming from high throughput sequencing technologies. Therefore, we have developed a QIIME-compatible database (DB) designed for species-level taxonomic assignment of 16S rRNA gene amplicon data targeting methanogenic archaea from the rumen, and from animal and human intestinal tracts. Called RIM-DB (Rumen and Intestinal Methanogen-DB), it contains a set of 2,379 almost full-length chimera-checked 16S rRNA gene sequences, including 20 previously unpublished sequences from isolates from three different orders. The taxonomy encompasses the recently-proposed seventh order of methanogens, the Methanomassiliicoccales, and allows differentiation between defined groups within this order. Sequence reads from rumen contents from a range of ruminant-diet combinations were taxonomically assigned using RIM-DB, Greengenes and SILVA. This comparison clearly showed that taxonomic assignments with RIM-DB resulted in the most detailed assignment, and only RIM-DB taxonomic assignments allowed methanogens to be distinguished taxonomically at the species level. RIM-DB complements the use of comprehensive databases such as Greengenes and SILVA for community structure analysis of methanogens from the rumen and other intestinal environments, and allows identification of target species for methane mitigation strategies.

Effects on enteric methane production and bacterial and archaeal communities by the addition of cashew nut shell extract or glycerol-an in vitro evaluation

The objective of the study was to evaluate the effect of cashew nut shell extract (CNSE) and glycerol (purity >99%) on enteric methane (CH4) production and microbial communities in an automated gas in vitro system. Microbial communities from the in vitro system were compared with samples from the donor cows, in vivo. Inoculated rumen fluid was mixed with a diet with a 60:40 forage:concentrate ratio and, in total, 5 different treatments were set up: 5mg of CNSE (CNSE-L), 10mg of CNSE (CNSE-H), 15mmol of glycerol/L (glycerol-L), and 30mmol of glycerol/L (glycerol-H), and a control without feed additive. Gas samples were taken at 2, 4, 8, 24, 32, and 48h of incubation, and the CH4 concentration was measured. Samples of rumen fluid were taken for volatile fatty acid analysis and for microbial sequence analyses after 8, 24, and 48h of incubation. In vivo rumen samples from the cows were taken 2h after the morning feeding at 3 consecutive days to compare the in vitro system with in vivo conditions. The gas data and data from microbial sequence analysis (454 sequencing) were analyzed using a mixed model and principal components analysis. These analyses illustrated that CH4 production was reduced with the CNSE treatment, by 8 and 18%, respectively, for the L and H concentration. Glycerol instead increased CH4 production by 8 and 12%, respectively, for the L and H concentration. The inhibition with CNSE could be due to the observed shift in bacterial population, possibly resulting in decreased production of hydrogen or formate, the methanogenic substrates. Alternatively the response could be explained by a shift in the methanogenic community. In the glycerol treatments, no main differences in bacterial or archaeal population were detected compared with the in vivo control. Thus, the increase in CH4 production may be explained by the increase in substrate in the in vitro system. The reduced CH4 production in vitro with CNSE suggests that CNSE can be a promising inhibitor of CH4 formation in the rumen of dairy cows.

Hydrogen and methane emissions from beef cattle and their rumen microbial community vary with diet, time after feeding and genotype

The aims of the present study were to quantify hydrogen (H2) and methane (CH4) emissions from beef cattle under different dietary conditions and to assess how cattle genotype and rumen microbial community affected these emissions. A total of thirty-six Aberdeen Angus-sired (AAx) and thirty-six Limousin-sired (LIMx) steers were fed two diets with forage:concentrate ratios (DM basis) of either 8:92 (concentrate) or 52:48 (mixed). Each diet was fed to eighteen animals of each genotype. Methane (CH4) and H2 emissions were measured individually in indirect respiration chambers. H2 emissions (mmol/min) varied greatly throughout the day, being highest after feed consumption, and averaged about 0·10 mol H2/mol CH4. Higher H2 emissions (mol/kg DM intake) were recorded in steers fed the mixed diet. Higher CH4 emissions (mol/d and mol/kg DM intake) were recorded in steers fed the mixed diet (P< 0·001); the AAx steers produced more CH4 on a daily basis (mol/d, P< 0·05) but not on a DM intake basis (mol/kg DM intake). Archaea (P= 0·002) and protozoa (P< 0·001) were found to be more abundant and total bacteria (P< 0·001) less abundant (P< 0·001) on feeding the mixed diet. The relative abundance of Clostridium cluster IV was found to be greater (P< 0·001) and that of cluster XIVa (P= 0·025) lower on feeding the mixed diet. The relative abundance of Bacteroides plus Prevotella was greater (P= 0·018) and that of Clostridium cluster IV lower (P= 0·031) in the LIMx steers. There were no significant relationships between H2 emissions and microbial abundance. In conclusion, the rate of H2 production immediately after feeding may lead to transient overloading of methanogenic archaea capacity to use H2, resulting in peaks in H2 emissions from beef cattle.

Discovery of a novel rumen methanogen in the anaerobic fungal culture and its distribution in the rumen as revealed by real-time PCR

Background

The novel archaea belonging to Rumen Cluster C (RCC), which may play an important role in methane production in the rumen have received increased attention. However, the present information on RCC in the rumen is limited by the unsuccessful isolation of axenic pure RCC from the rumen. In the present study, RCC grown in anaerobic fungal subcultures was identified by the molecular and culture methods.

Results

A novel RCC species existing in the fungal subcultures was identified and demonstrated by the 16S rRNA gene clone library. Interestingly, the novel RCC species survived in the fungal cultures over all the subculture transferring, even in the 62^nd subculture, in contrast to the other methanogens, which disappeared during subcultures. Further work showed that subculture transfer frequency significantly affected the relative abundance of the novel RCC species in the fungal subcultures. The five-day and seven-day transfer frequencies increased the relative abundance of the RCC species (P<0.05). In addition, quantitative real-time PCR revealed that high concentrate diets did not affect the abundance of archaea, but numerically reduced the abundance of the novel RCC species in the rumen. In addition, the relative abundance of the RCC species was numerically higher in the rumen liquid fraction than in the rumen epithelium and solid fractions. Finally, a purified fungal culture containing the RCC species was successfully obtained. PCR and sequencing analysis showed that the novel RCC species contained a mcrA gene, which is known to play a crucial role in methanogenesis, and thus could be identified as a methanogen.

Conclusion

In this study, a novel RCC species was identified as a methanogen and closely associated with anaerobic fungi. This novel approach by using co-culture with anaerobic fungi may provide a feasible way to culture and investigate not yet identified methanogens.

An antimethanogenic nutritional intervention in early life of ruminants modifies ruminal colonization by Archaea

The aim of this work was to study whether feeding a methanogen inhibitor from birth of goat kids and their does has an impact on the archaeal population colonizing the rumen and to what extent the impact persists later in life. Sixteen goats giving birth to two kids were used. Eight does were treated (D+) with bromochloromethane after giving birth and over 2 months. The other 8 goats were not treated (D−). One kid per doe in both groups was treated with bromochloromethane (k+) for 3 months while the other was untreated (k−), resulting in four experimental groups: D+/k+, D+/k−, D−/k+, and D−/k−. Rumen samples were collected from kids at weaning and 1 and 4 months after (3 and 6 months after birth) and from does at the end of the treating period (2 months). Pyrosequencing analyses showed a modified archaeal community composition colonizing the rumen of kids, although such effect did not persist entirely 4 months after; however, some less abundant groups remained different in treated and control animals. The different response on the archaeal community composition observed between offspring and adult goats suggests that the competition occurring in the developing rumen to occupy different niches offer potential for intervention.