Background: Recently, we reported that some dairy cows could produce high amounts of milk with high amounts of protein (defined as milk protein yield [MPY]) when a population was raised under the same nutritional and management condition, a potential new trait that can be used to increase high-quality milk production. It is unknown to what extent the rumen microbiome and its metabolites, as well as the host metabolism, contribute to MPY. Here, analysis of rumen metagenomics and metabolomics, together with serum metabolomics was performed to identify potential regulatory mechanisms of MPY at both the rumen microbiome and host levels. Results: Metagenomics analysis revealed that several Prevotella species were significantly more abundant in the rumen of high-MPY cows, contributing to improved functions related to branched-chain amino acid biosynthesis. In addition, the rumen microbiome of high-MPY cows had lower relative abundances of organisms with methanogen and methanogenesis functions, suggesting that these cows may produce less methane. Metabolomics analysis revealed that the relative concentrations of rumen microbial metabolites (mainly amino acids, carboxylic acids, and fatty acids) and the absolute concentrations of volatile fatty acids were higher in the high-MPY cows. By associating the rumen microbiome with the rumen metabolome, we found that specific microbial taxa (mainly Prevotella species) were positively correlated with ruminal microbial metabolites, including the amino acids and carbohydrates involved in glutathione, phenylalanine, starch, sucrose, and galactose metabolism. To detect the interactions between the rumen microbiome and host metabolism, we associated the rumen microbiome with the host serum metabolome and found that Prevotella species may affect the host’s metabolism of amino acids (including glycine, serine, threonine, alanine, aspartate, glutamate, cysteine, and methionine). Further analysis using the linear mixed effect model estimated contributions to the variation in MPY based on different omics and revealed that the rumen microbial composition, functions, and metabolites, and the serum metabolites contributed 17.81, 21.56, 29.76, and 26.78%, respectively, to the host MPY. Conclusions: These findings provide a fundamental understanding of how the microbiome-dependent and host-dependent mechanisms contribute to varied individualized performance in the milk production quality of dairy cows under the same management condition. This fundamental information is vital for the development of potential manipulation strategies to improve milk quality and production through precision feeding.

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. 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 (h2) ≥ 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 h2 ≥ 0.15) rumen microbial taxa, and five of these SNPs were known quantitative trait loci for feed efficiency in cattle. 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.

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.

Cattle with high feed efficiencies (designated “efficient”) produce less methane gas than those with low feed efficiencies (designated “inefficient”); however, the role of the methane producers in such difference is unknown. This study investigated whether the structures and populations of methanogens in the rumen were associated with differences in cattle feed efficiencies by using culture-independent methods. Two 16S rRNA libraries were constructed using ∼800-bp amplicons generated from pooled total DNA isolated from efficient ( n = 29) and inefficient ( n = 29) animals. Sequence analysis of up to 490 randomly selected clones from each library showed that the methanogenic composition was variable: less species variation (22 operational taxonomic units [OTUs]) was detected in the rumens of efficient animals, compared to 27 OTUs in inefficient animals. The methanogenic communities in inefficient animals were more diverse than those in efficient ones, as revealed by the diversity indices of 0.84 and 0.42, respectively. Differences at the strain and genotype levels were also observed and found to be associated with feed efficiency in the host. No difference was detected in the total population of methanogens, but the prevalences of Methanosphaera stadtmanae and Methanobrevibacter sp. strain AbM4 were 1.92 ( P < 0.05) and 2.26 ( P < 0.05) times higher in inefficient animals, while Methanobrevibacter sp. strain AbM4 was reported for the first time to occur in the bovine rumen. Our data indicate that the methanogenic ecology at the species, strain, and/or genotype level in the rumen may play important roles in contributing to the difference in methane gas production between cattle with different feed efficiencies.

Microorganisms 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. Rumen 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. Metatranscriptomes 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.

Linkage of rumen microbial structure to host phenotypical traits may enhance the understanding of host–microbial interactions in livestock species. T