Background The symbiotic rumen microbiota is essential for the digestion of plant fibers and contributes to the variation of production and health traits in ruminants. However, to date, the heritability of rumen microbial features and host genetic components associated with the rumen microbiota, as well as whether such genetic components are animal performance relevant, are largely unknown. Results In the present study, we assessed rumen microbiota from a cohort of 709 beef cattle and showed that multiple factors including breed, sex, and diet drove the variation of rumen microbiota among animals. The diversity indices, the relative abundance of ~ 34% of microbial taxa (59 out of 174), and the copy number of total bacteria had a heritability estimate ( h 2 ) ≥ 0.15, suggesting that they are heritable elements affected by host additive genetics. These moderately heritable rumen microbial features were also found to be associated with host feed efficiency traits and rumen metabolic measures (volatile fatty acids). Moreover, 19 single nucleotide polymorphisms (SNPs) located on 12 bovine chromosomes were found to be associated with 14 (12 of them had h 2 ≥ 0.15) rumen microbial taxa, and five of these SNPs were known quantitative trait loci for feed efficiency in cattle. Conclusions These findings suggest that some rumen microbial features are heritable and could be influenced by host genetics, highlighting a potential to manipulate and obtain a desirable and efficient rumen microbiota using genetic selection and breeding. It could be a useful strategy to further improve feed efficiency and optimize rumen fermentation through targeting both cattle and their rumen microbiota. Electronic supplementary material The online version of this article (10.1186/s40168-019-0699-1) contains supplementary material, which is available to authorized users.

Exploring compositional and functional characteristics of the rumen microbiome can improve the understanding of its role in rumen function and cattle feed efficiency. In this study, we applied metatranscriptomics to characterize the active rumen microbiomes of beef cattle with different feed efficiencies (efficient, n ϭ 10; inefficient, n ϭ 10) using total RNA sequencing. Active bacterial and archaeal compositions were estimated based on 16S rRNAs, and active microbial metabolic functions including carbohydrate-active enzymes (CAZymes) were assessed based on mRNAs from the same metatranscriptomic data sets. In total, six bacterial phyla (Proteobacteria, Firmicutes, Bacteroidetes, Spirochaetes, Cyanobacteria, and Synergistetes), eight bacterial families (Succinivibrionaceae, Prevotellaceae, Ruminococcaceae, Lachnospiraceae, Veillonellaceae, Spirochaetaceae, Dethiosulfovibrionaceae, and Mogibacteriaceae), four archaeal clades (Methanomassiliicoccales, Methanobrevibacter ruminantium, Methanobrevibacter gottschalkii, and Methanosphaera), 112 metabolic pathways, and 126 CAZymes were identified as core components of the active rumen microbiome. As determined by comparative analysis, three bacterial families (Lachnospiraceae, Lactobacillaceae, and Veillonellaceae) tended to be more abundant in low-feed-efficiency (inefficient) animals (P Ͻ 0.10), and one archaeal taxon (Methanomassiliicoccales) tended to be more abundant in high-feed-efficiency (efficient) cattle (P Ͻ 0.10). Meanwhile, 32 microbial metabolic pathways and 12 CAZymes were differentially abundant (linear discriminant analysis score of Ͼ2 with a P value of Ͻ0.05) between two groups. Among them, 30 metabolic pathways and 11 CAZymes were more abundant in the rumen of inefficient cattle, while 2 metabolic pathways and 1 CAZyme were more abundant in efficient animals. These findings suggest that the rumen microbiomes of inefficient cattle have more diverse activities than those of efficient cattle, which may be related to the host feed efficiency variation.IMPORTANCE This study applied total RNA-based metatranscriptomics and showed the linkage between the active rumen microbiome and feed efficiency (residual feed intake) in beef cattle. The data generated from the current study provide fundamental information on active rumen microbiome at both compositional and functional levels, which serve as a foundation to study rumen function and its role in cattle feed efficiency. The findings that the active rumen microbiome may contribute to variations in feed efficiency of beef cattle highlight the possibility of enhancing nutrient utilization and improve cattle feed efficiency through modification of rumen microbial functions.KEYWORDS metatranscriptome, RNA-seq, rumen microbiome, beef cattle, feed efficiency R uminants can convert fibrous plant materials into edible meat and milk products for human consumption through rumen microbial fermentation. The symbiotic rumen microbiota mainly consists of bacteria, archaea, fungi, and ciliated protozoa (1)

BackgroundMicroorganisms are responsible for fermentation within the rumen and have been reported to contribute to the variation in feed efficiency of cattle. However, to what extent the breed affects the rumen microbiome and its association with host feed efficiency is unknown. Here, rumen microbiomes of beef cattle (n = 48) from three breeds (Angus, Charolais, Kinsella composite hybrid) with high and low feed efficiency were explored using metagenomics and metatranscriptomics, aiming to identify differences between functional potentials and activities of same rumen microbiomes and to evaluate the effects of host breed and feed efficiency on the rumen microbiome.ResultsRumen metagenomes were more closely clustered together and thus more conserved among individuals than metatranscriptomes, suggesting that inter-individual functional variations at the RNA level were higher than those at the DNA level. However, while mRNA enrichment significantly increased the sequencing depth of mRNA and generated similar functional profiles to total RNA-based metatranscriptomics, it led to biased abundance estimation of several transcripts. We observed divergent rumen microbial composition (metatranscriptomic level) and functional potentials (metagenomic level) among three breeds, but differences in functional activity (metatranscriptomic level) were less apparent. Differential rumen microbial features (e.g., taxa, diversity indices, functional categories, and genes) were detected between cattle with high and low feed efficiency, and most of them were breed-specific.ConclusionsMetatranscriptomes represent real-time functional activities of microbiomes and have the potential to better associate rumen microorganisms with host performances compared to metagenomics. As total RNA-based metatranscriptomics seem to avoid potential biases caused by mRNA enrichment and allow simultaneous use of rRNA for generation of compositional profiles, we suggest their use for linking the rumen microbiome with host phenotypes in future studies. However, if exploration of specific lowly expressed genes is desired, mRNA enrichment is recommended as it will enhance the resolution of mRNA. Finally, the differential microbial features observed between efficient and inefficient steers tended to be specific to breeds, suggesting that interactions between host breed genotype and the rumen microbiome contribute to the variations in feed efficiency observed. These breed-associated differences represent an opportunity to engineer specific rumen microbiomes through selective breeding of the hosts.Electronic supplementary materialThe online version of this article (10.1186/s40168-019-0618-5) contains supplementary material, which is available to authorized users.

Taxonomic characterization of active gastrointestinal microbiota is essential to detect shifts in microbial communities and functions under various conditions. This study aimed to identify and quantify potentially active rumen microbiota using total RNA sequencing and to compare the outcomes of this approach with the widely used targeted RNA/DNA amplicon sequencing technique. Total RNA isolated from rumen digesta samples from five beef steers was subjected to Illumina paired-end sequencing (RNA-seq), and bacterial and archaeal amplicons of partial 16S rRNA/rDNA were subjected to 454 pyrosequencing (RNA/DNA Amplicon-seq). Taxonomic assessments of the RNA-seq, RNA Amplicon-seq, and DNA Amplicon-seq datasets were performed using a pipeline developed in house. The detected major microbial phylotypes were common among the three datasets, with seven bacterial phyla, fifteen bacterial families, and five archaeal taxa commonly identified across all datasets. There were also unique microbial taxa detected in each dataset. Elusimicrobia and Verrucomicrobia phyla; Desulfovibrionaceae, Elusimicrobiaceae, and Sphaerochaetaceae families; and Methanobrevibacter woesei were only detected in the RNA-Seq and RNA Amplicon-seq datasets, whereas Streptococcaceae was only detected in the DNA Amplicon-seq dataset. In addition, the relative abundances of four bacterial phyla, eight bacterial families and one archaeal taxon were different among the three datasets. This is the first study to compare the outcomes of rumen microbiota profiling between RNA-seq and RNA/DNA Amplicon-seq datasets. Our results illustrate the differences between these methods in characterizing microbiota both qualitatively and quantitatively for the same sample, and so caution must be exercised when comparing data.

The objective of the current study was to employ a DNA-based sequencing technology to study the effect of active dry yeast (ADY) supplementation, diet type, and sample location within the rumen on rumen bacterial community diversity and composition, and to use an RNA-based method to study the effect of ADY supplementation on rumen microbial metabolism during high-grain feeding (HG). Our previous report demonstrated that the supplementation of lactating dairy cows with ADY attenuated the effect of subacute ruminal acidosis. Therefore, we used samples from that study, where 16 multiparous, rumen-cannulated lactating Holstein cows were randomly assigned to 1 of 2 dietary treatments: ADY (Saccharomyces cerevisiae strain Y1242, 80 billion cfu/animal per day) or control (carrier only). Cows received a high-forage diet (77:23, forage:concentrate), then were abruptly switched to HG (49:51, forage:concentrate). Rumen bacterial community diversity and structure were highly influenced by diet and sampling location (fluid, solids, epimural). The transition to HG reduced bacterial diversity, but epimural bacteria maintained a greater diversity than fluid and solids. Analysis of molecular variance indicated a significant separation due to diet × sampling location, but not due to treatment. Across all samples, the analysis yielded 6,254 nonsingleton operational taxonomic units (OTU), which were classified into several phyla: mainly Firmicutes, Bacteroidetes, Fibrobacteres, Tenericutes, and Proteobacteria. High forage and solids were dominated by OTU from Fibrobacter, whereas HG and fluid were dominated by OTU from Prevotella. Epimural samples, however, were dominated in part by Campylobacter. Active dry yeast had no effect on bacterial community diversity or structure. The phylum SR1 was more abundant in all ADY samples regardless of diet or sampling location. Furthermore, on HG, OTU2 and OTU3 (both classified into Fibrobacter succinogenes) were more abundant with ADY in fluid and solids than control samples. This increase with ADY was paralleled by a reduction in prominent Prevotella OTU. Metatranscriptomic profiling of rumen microbiome conducted on random samples from the HG phase showed that ADY increased the abundance of the cellulase endo-β-1,4-glucanase and had a tendency to increase the hemicellulase α-glucuronidase. In conclusion, the shift from high forage to HG and sampling location had a more significant influence on ruminal bacterial community abundance and structure compared with ADY. However, evidence suggested that ADY can increase the abundance of some dominant anaerobic OTU belonging to F. succinogenes and phylum SR1. Further, microbial mRNA-based evidence suggested that ADY can increase the abundance of a specific microbial fibrolytic enzymes.