Metagenomics reveals differences in microbial composition and metabolic functions in the rumen of dairy cows with different residual feed intake
Animal Microbiome volume 4, Article number: 19 (2022)
Rumen microbial composition and functions have vital roles in feed digestion and fermentation and are linked to feed efficiency in cattle. This study selected Holstein cows, which are high in both milk protein content and milk yield, to analyse the relationship between the rumen microbiota and residual feed intake (RFI). Eighteen multiparous lactating cows were divided into low RFI (LRFI, high efficiency, n = 9) and high RFI (HRFI, low efficiency, n = 9) groups to investigate the differences in microbial composition and functions.
The relative abundances of butyrate producers, including the Clostridium, Butyrivibrio, Eubacterium and Blautia genera, were higher in HRFI cows than in LRFI cows (P < 0.05). Four carbohydrate metabolic pathways (glycolysis/gluconeogenesis, pentose phosphate pathway, fructose and mannose metabolism, and butanoate metabolism) and one energy metabolism pathway (methane metabolism), were more abundant in HRFI animals (P < 0.05). Quorum sensing and DNA replication pathways were more abundant in HRFI cows. For CAZyme profiles, 14 out of 19 genes encoding carbohydrates-deconstructing enzymes were more abundant in HRFI cows (P < 0.05). Seven Lachnospiraceae species associated with carbohydrate metabolism and quorum sensing may contribute to the difference in feed efficiency. Moreover, the LRFI cows had lower abundances of Methanosphaera (P < 0.01), Methanobrevibacter ruminantium (P = 0.09) and methanogenesis functions (P = 0.04).
The rumen microbiota of low-efficiency cows has stronger abilities to degrade carbohydrates and produce methane, and quorum sensing pathways could also be associated with differences in feed efficiency. This study provides a deeper understanding of the microbial ecology of dairy cows with different feed efficiencies and highlights the possibility of modulating the rumen microbiome or microbial functions to improve the feed efficiency of dairy cows.
Feed efficiency is an important trait in the dairy industry since feed inputs account for 40–60% of the total variable costs in dairy farms . Improving the feed efficiency of dairy cows has profound significance, as it is expected to decrease the total feed cost and reduce greenhouse gas emissions and urinary nitrogen excretion. Additionally, there are various methods to evaluate feed efficiency, among which residual feed intake (RFI) is one of the most common measurements in cattle [2, 3]. The RFI is defined as the difference between the actual and predicted dry matter intake (DMI) for an animal. A low or negative RFI value represents high efficiency, while a high RFI value represents low efficiency. The RFI is independent of the level of production, and animals with lower values are deemed more efficient .
Physiological processes have been identified as important factors contributing to the differences in RFI among dairy cows, including feeding behaviour, activity, digestibility, and rumen-temperature variables [5, 6]. Several studies have reported that rumen fermentation variables and nutrient absorption affect the RFI of cattle [7, 8]. Rumen fermentation is largely mediated by diverse rumen microbes to produce the major products volatile fatty acids (VFAs), which can provide greater than 70% of the energy requirement for ruminants . Therefore, the rumen microbiome plays a vital role in variations in RFI.
Previous studies have reported differences in the rumen microbial composition between efficient and inefficient beef cattle [10, 11]. McGovern et al.  found that the abundances of Lactobacillales and Ruminobacter genera were associated with RFI. Lower abundances of Methanosphaera stadtmanae and Methanobrevibacter sp. has been identified in high RFI animals . Most of these studies used beef cattle or growing heifers, but there are distinct differences between dairy and beef cattle. Limited studies have shown a correlation between rumen microbial ecology and RFI in lactating dairy cattle . However, a high concentrate diet (70% concentrate) may lead to abnormal ruminal microbial metabolism and further affect the rumen microbiota. Thus, it is important to examine the relationship between RFI and the microbiome in cows fed an optimal forage-to-concentrate ratio. Furthermore, a previous study did not address the impact of differences in individual performance (e.g., milk protein yield) on microorganisms . Xue et al.  revealed the differences in the rumen microbiome and its metabolites between individual cows with different milk protein yields and indicated a correlation between rumen ecology and host performance. In the present study, we profiled the rumen microbial composition and functions in dairy cows with relatively high milk yield and milk protein content that had divergent RFI to uncover the relationship between RFI and microbiome components. The results will provide a better understanding of the rumen microbiome and its impact on the RFI of lactating dairy cattle.
Characterization of phenotypes
The low RFI (LRFI) cows consumed 2.43 kg dry matter/d less than their high RFI (HRFI) counterparts. The ratios of milk yield to DMI (P < 0.01) and of energy-corrected milk to DMI (P = 0.03) were greater in LRFI group than in HRFI group (Additional file 1: Table S1). Although the percentage of milk protein was significantly lower (P = 0.04) in the LRFI group, no significant differences (P > 0.10) were detected in the yields of milk and milk protein between the two groups.
The concentrations of total VFAs (P = 0.08), acetate (P = 0.10), butyrate (P = 0.10), valerate (P = 0.08) and isovalerate (P = 0.09) tended to be lower in the LRFI group (P < 0.10, Table 1). However, the molar proportion of propionate tended to be higher in LRFI cows (24.0 vs. 22.7, P = 0.09), and the molar proportions of isobutyrate (0.90 vs. 0.74, P = 0.05) and isovalerate (1.50 vs. 1.23, P < 0.01) were significantly greater in LRFI cows than HRFI cows. Methane (CH4) production was higher in HRFI cows than in LRFI cows (P = 0.04).
The rumen metagenome profile
A total of 909,094,550 raw reads were obtained from the 18 rumen digesta samples, with an average of 50,505,253 ± 928,365 (SEM) raw reads per sample. After quality control and host gene removal, 901,313,532 clean reads, with 50,072,974 ± 921,079 raw reads per sample, were retained (Additional file 1: Table S2). Multiple megahit assembly of the raw sequencing reads resulted in 11,670,011 contigs, with an average of 648,334 ± 44,632 contigs per sample. At the domain level, the rumen metagenome consisted of 94.84% bacteria, 3.10% eukaryotes, 1.33% archaea, and 0.51% viruses (Additional file 1: Fig. S1).
Differences in microbial community compositions between LRFI and HRFI cows
The permutational multivariate analysis of variance (PERMANOVA) showed no differences between HRFI and LRFI groups in bacteria (P = 0.17), eukaryotes (P = 0.25) archaea (P = 0.15) or viruses (P = 0.13, Additional file 1: Table S3). The principal coordinate analysis (PCoA) did not show clear separation between the HRFI and LRFI groups based on archaeal species (Additional file 1: Fig. S2a). Meanwhile, no significant difference was found between the HRFI and LRFI groups in archaeal phyla. However, the abundance of Methanosphaera was higher in HRFI group (4.66 vs. 3.28, P < 0.01). The relative abundance of Methanobrevibacter ruminantium, the most abundant archaeal species, tended to be higher in HRFI cows than in LRFI cows (22.45 vs. 20.14, P = 0.09, Fig. 1).
The PCoA of bacterial species showed no clear separation between the HRFI and LRFI groups (Additional file 1: Fig. S2b). At the phylum level, compared with LRFI group, the relative abundance of Firmicutes was significantly greater in HRFI group (36.57 vs. 31.91, P = 0.04), and Proteobacteria tended to be lower in the HRFI cows than in LRFI cows (4.45 vs. 7.61., P = 0.09, Fig. 2A). At the genus level, the abundances of Clostridium, Butyrivibrio, Eubacterium and Blautia were significantly higher in HRFI cows (P < 0.05, Fig. 2B). At the species level, 28 species, including 24 Firmicutes sp., showed significantly higher abundances in HRFI group (LDA > 2, P < 0.05), while 10 species were more abundant in the LRFI group (LDA > 2, P < 0.05, Fig. 2C).
Differences in functions of the rumen microbiome between LRFI and HRFI cows
For Kyoto Encyclopedia of Genes and Genomes (KEGG) profiles, a total of 206 endogenous level-3 pathways were observed in our study. These pathways belonged to 4 level-1 KEGG functional categories, including “metabolism” (69.67 ± 0.47%), “genetic information processing” (15.77 ± 0.07%), “cellular processes” (7.38 ± 0.02%) and “environment information processing” (7.18 ± 0.02%). At the second level of KEGG pathways, a total of 24 functional categories were identified. The categories “carbohydrate metabolism” (16.54 ± 0.14%), “global and overview maps” (11.62 ± 0.13%), “amino acid metabolism” (10.25 ± 0.10%), “replication and repair” (8.07 ± 0.10%), and “nucleotide metabolism” (10.25 ± 0.10%) were the most abundant functions.
The PCoA indicated no separation in functional potential between the two feed efficiency cohorts (Additional file 1: Fig. S3a). At KEGG level 2, the category “carbohydrate metabolism” tended to be higher in the HRFI animals than in the LRFI animals (P = 0.06, Fig. 3A). At KEGG level 3, 13 pathways, including 10 “metabolism” pathways, 2 “genetic information processing” pathways, and 1 “cellular processes” pathway, were upregulated in HRFI cows (P < 0.05; Fig. 3B), and two KEGG pathways were significantly enriched in LRFI cows (P < 0.05; Fig. 3B). Within the category “carbohydrate metabolism”, “ko00010: Glycolysis/gluconeogenesis”, “ko00030: pentose phosphate pathway”, “ko00051: fructose and mannose metabolism”, and “ko00650: butanoate metabolism” were significantly enriched in the HRFI cows (P < 0.05; Fig. 3B). For the category “energy metabolism”, only “ko00680: methane metabolism” was enriched in HRFI cows (P = 0.04, Fig. 3B). For the other highly abundant biological processes, “DNA replication” and “quorum sensing” were higher in the HRFI group (P < 0.05; Fig. 3). However, “ko00100: steroid biosynthesis” and “ko00909: sesquiterpenoid and triterpenoid biosynthesis” were significantly enriched in LRFI cows (P < 0.05; Fig. 3B).
CAZyme functions of the rumen microbiome between LRFI and HRFI cows
In terms of CAZyme profiles, 199 genes were detected, including 97 glycoside hydrolases (GHs), 42 glycosyltransferases (GTs), 31 carbohydrate-binding modules (CBMs), 14 carbohydrate esterases (CEs), 10 polysaccharide lyases (PLs), and 5 auxiliary activities (AAs). Among them, GT2 (6.99% ± 0.15%), GH2 (6.33% ± 0.08%), and GT4 (3.85% ± 0.05%) were the most abundant CAZymes in the rumen of these cows.
No obvious separation was observed for the CAZyme profiles between feed efficiency groups based on PCoA (Additional file 1: Fig. S3b). We then compared the relative abundance of CAZymes, and 19 differentially abundant CAZymes related to deconstructing carbohydrates (GH, CE, PL, AA, and CBM) were identified between the two RFI groups. Only five CAZymes (4 GHs, 1 CE) were more abundant in LRFI animals, while 14 CAZymes (9 GHs, 3 CEs, and 2 AAs) were enriched in HRFI animals. Among the GTs, five CAZymes were more abundant in the LRFI group, while four CAZymes were enriched in HRFI cows. Regarding the CBMs, three CAZymes were more abundant in the LRFI group, while only one was enriched in the HRFI group (P < 0.05; Fig. 4).
Relationships between microbial species and functions
Spearman’s rank correlation was constructed between differential microbial species and metabolic pathways. A total of 38 species showed significant relationships with 15 differentially abundant pathways (R > 0.50 and P < 0.05, Fig. 5), and 22 species had positive relationships with “quorum sensing” pathway (R > 0.50 and P < 0.05. Figure 5). Among these 22 species, 15 also had a positive relationship with “carbohydrate metabolism” pathways (ko00010, ko00030, ko00051, and ko00650, Fig. 5). Among these 19 species, 7 belonged to the family Lachnospiraceae, including Ruminococcus torques, Clostridium clostridioforme, Blautia obeum, B. producta, L. bacterium FE2018, Marvinbryantia formatexigens and an unclassified Butyrivibrio sp. (Fig. 5).
Microbial functions and species related to carbohydrate metabolism, energy metabolism and other processes in the rumen of the animals are summarized in Fig. 6. Among the genes encoding enzymes involved in the fructose and mannose metabolism pathway, 8 genes, EC 126.96.36.199, EC 188.8.131.52, EC 184.108.40.206, EC 220.127.116.11, EC 18.104.22.168, EC 22.214.171.124, EC 126.96.36.199 and EC 188.8.131.52, tended to be more abundant in the HRFI cows (P < 0.10). In the glycolysis pathway, the abundance of EC 184.108.40.206, EC 220.127.116.11 and EC 18.104.22.168 tended to be higher in the HRFI cows (P < 0.10), and EC 22.214.171.124 was more abundant in the HRFI cows (P = 0.03). In the pentose phosphate pathway, two genes (EC 126.96.36.199 and EC 188.8.131.52) tended to be more abundant in the HRFI group (P < 0.10), and 5 genes (EC 184.108.40.206, EC 220.127.116.113, EC 18.104.22.168, EC 22.214.171.124 and EC 126.96.36.199) were enriched in the HRFI group (P < 0.05). Three genes (EC 188.8.131.52, EC 184.108.40.206 and EC 220.127.116.11) involved in quorum sensing pathways, and 6 genes (EC 18.104.22.168, EC 22.214.171.124, EC 126.96.36.199, EC 188.8.131.52, EC 184.108.40.206 and EC 220.127.116.11) involved in methane metabolism were significantly enriched in the HRFI group (P < 0.05). In the butyrate metabolism pathway, four genes (EC 18.104.22.168, EC 22.214.171.124, EC 126.96.36.199 and EC 188.8.131.52) tended to be more abundant in HRFI cows (P < 0.10), and 3 genes (EC 184.108.40.206, EC 220.127.116.11, EC 18.104.22.168) related to butyrate production were significantly enriched in HRFI cows (P < 0.05). However, in the succinate pathway for propionate production, three genes encoding enzymes (EC 22.214.171.124, EC 126.96.36.199, EC 188.8.131.52) were more abundant in the LRFI cows.
Improvement in feed efficiency could greatly reduce feed costs and increase profits. In our study, LRFI cows with similar milk production consumed 2.50 kg less dry matter per day than the HRFI cows did, which is consistent with previous studies [15, 16]. The conversion of feed to VFAs is dependent on rumen microbes, and the rumen microbiome plays an important role in determining an animal’s RFI [17, 18]. In our study, the microbial differences at the taxonomic and functional levels between the two RFI groups explain the underlying mechanisms.
A greater amount of bacterial species (28 vs. 10) were more abundant in HRFI cows than in LRFI cows, which is in line with previous findings showing that HRFI animals tended to have higher microbiome richness . Our results found higher abundances of Clostridium, Butyrivibrio, Eubacterium and Blautia in HRFI cows. Clostridium species ferment carbohydrates to form hydrogen and butyrate . The genus Butyrivibrio plays an important role in hemicellulose and pectin breakdown in the rumen, and the end product of this fermentation is butyrate [20, 21]. The genera Eubacterium and blautia possess the capability to produce butyrate . Moreover, most of the species that abundances were higher in HRFI cows were butyrate producers and showed positive associations with the “ko00650: butanoate metabolism” pathway, suggesting their important roles in butyrate biosynthesis. Indeed, we detected a trend towards higher butyrate concentrations in HRFI cows.
Methanosphaera is a methanogenic archaeon that utilizes methyl group compounds to generate methane . In the rumen, methanol is a product of pectin hydrolysis by protozoa and esterase activity of bacteria . The higher DMI in HRFI cows means more substrates available for fermentation in the rumen. The greater potential for increased availability of methanol is associated with the higher abundance of Methanosphaera in HRFI group in the current study, suggesting the importance of Methanosphaera in affecting host feed efficiency. The higher abundance of Methanobrevibacter ruminantium in HRFI cows was consistent with the study of Shabat et al. . These results suggest that HRFI cows may produce more methane, leading to energy waste and lower efficiency. Moreover, methane metabolism and KO genes related to methanogenesis were more abundant in HRFI cows. In line with our results, previous studies also reported a higher abundance of methanogenesis pathways in HRFI cows . Methane emissions cause a significant loss of dietary energy and reduce the efficiency of animals . A higher abundance of KO genes related to propionate production and a higher tendency of the molar proportion of propionate in LRFI cows indicates that the rumen microbiomes of LRFI cows are more likely to generate substances with greater usable energy values. As mentioned above, we found a lower tendency of butyrate concentrations and higher tendency of the molar proportion of propionate in LRFI cows. Moss et al.  reported that formation of butyrate results in production of additional methanogenic substrates (formate and H2), while propionate formation is involved in hydrogen utilisation, which can be considered as a competitive pathway for hydrogen use in the rumen. Thus, the lower butyrate concentrations and higher molar proportion of propionate in LRFI cows may decrease methanogenesis by competing for hydrogen.
The functional genes of the rumen microbiome identified by metagenomics provide a way to evaluate the functions of rumen microorganisms [26, 27]. In a recent study, Xue et al.  found that rumen microbial functions (21.56%) make a greater contribution to milk protein yield than rumen microbial composition (17.81%). These findings suggest that the rumen microbial function may reflect the mechanism by which microorganisms influence phenotypes. Previous studies used metagenomics  or metatranscriptomics  to investigate the functional profiles of rumen microbiomes in cows with different RFI. A metagenomic study using rumen samples from beef cattle reported a higher number of KEGG pathways enriched in HRFI cattle , consistent with the results in our study. Shabat et al.  investigated the linkages between the rumen microbiome and host efficiency in dairy cows with different RFI and revealed less diverse metabolic pathways in efficient cows, suggesting that efficient microbiomes are more specialized to meet the energetic need of ruminants. Additionally, we found that the KEGG functions of carbohydrate degradation, including “glycolysis”, “pentose phosphate pathway” and “fructose and mannose metabolism”, were enriched in the rumen of HRFI animals. These results indicate that HRFI microbiome may generate more diverse products such as pyruvate, acetyl-CoA and hydrogen. In the current study, HRFI cows had a higher number of CAZyme-encoding genes. Similar results were reported at the RNA level in beef cattle in a metatranscriptomics study . The higher abundance of KO genes related to carbohydrate degradation in HRFI cows indicates higher activities of degrading complex substrates in inefficient cows. Most species showing positive relationships with carbohydrate metabolism belonged to the Lachnospiraceae family. Many members of the family Lachnospiraceae have cellulolytic activity and are related to butyrate production [28, 29]. In the current study, we found increased butyrate metabolism in HRFI cows. Therefore, the Lachnospiraceae species may play an essential role in carbohydrate metabolism and further affect feed efficiency in dairy cows. However, it is unclear whether the rumen microbiota is either the causes or results of cattle feed efficiency. Thus, our results warrant further research to fully elucidate the relationship between the alteration in microbiota and host functional changes.
The higher activity of rumen microbiomes in HRFI animals could be attributed to their higher feed intake. Pathak  reported a strong positive correlation of microbial activities with feed intake because more available substrates and nutrients from animals with higher feed intake were provided for microbial growth. Indeed, the higher abundance of the DNA replication pathway in our current study provides further evidence for the greater microbial growth rate in HRFI cows. The higher microbial growth rate may result in higher population density in bacteria. Enrichment of the quorum sensing pathway in HRFI cows has not been reported previously. In a recent study, some of the bacterial species in rumen were identified the existence of luxS genes, and these species showed positive relationships with quorum sensing pathway in our current study and others . Moreover, LuxS proteins, which are encoded by the luxS gene, were proven to be critical enzymes for activating the bacterial methyl cycle and for producing AI-2 , a member of a family of signalling molecules used in quorum sensing. These results indicate that quorum sensing could be associated with the difference in feed efficiency between dairy cows with divergent RFI values. Additionally, the bacterial species showing positive associations with carbohydrate metabolism were all positively correlated with quorum sensing pathways. In short, several Lachnospiraceae species, including Ruminococcus torques, Blautia obeum, Blautia producta, Blautia schinkii, Blautia wexlerae, Dorea longicatena, Clostridium clostridioforme, Clostridium symbiosum, Marvinbryantia formatexigens and Lachnospiraceae bacterium FE2018, may play a vital role in both carbohydrate metabolism and quorum sensing. Further studies are required to fully elucidate the relationship between the key microbes and pathways described above to validate our current findings.
The current study revealed the taxonomic features and functions of rumen microbiomes. Although carbohydrate degradation functions were more abundant in HRFI cows (low efficiency), the higher abundances of butyrate-producing species, functions in butyrate metabolism, methanogenic archaea and methanogenesis pathways indicate that HRFI cows may promote methane production because of their tendency towards butyrate production. Compared with HRFI group, the higher abundances of KO genes related to propionate production and higher tendency of the molar proportion of propionate in LRFI cows may partly explain the higher efficiency in these cows. Several Lachnospiraceae species associated with carbohydrate metabolism and quorum sensing may contribute to the differences in RFI. This study provides a in-depth understanding of the microbial ecology of dairy cows with different RFI, and highlights the possibility of modulating the rumen microbiome or microbial functions to improve the feed efficiency of dairy cows. Future studies are required to assess information on the metabolic intermediates related to carbohydrate metabolism, which will provide evidence for the effect of microbial carbohydrate metabolism on feed efficiency.
Materials and methods
Experimental design and sample collection
All the procedures involving animals in this study were approved by the Animal Use and Care Committee of Zhejiang University (Hangzhou, China, No. 12410). The RFI value of each animal was calculated as described previously . Briefly, the RFI was estimated as the difference between expected feed intake and actual feed intake, where the expected feed intake was computed through a multiple linear regression model using the regression of actual feed intake on energy-corrected milk yield, metabolic body weight (BW0.75), and average daily gain over the measurement period. Based on previous RFI value calculations , the 9 dairy cows with the lowest RFI (LRFI) and 9 with the highest RFI (HRFI) were selected from the cohort of 53 multiparous mid-lactation Holstein dairy cows. Power calculations revealed that our sample size will produce 99.0% power while controlling type I error under 5%, based on a t test of RFI. The RFI coefficients distribution of the LRFI and HRFI groups are depicted in Additional file 1: Fig. S4. Cows were housed in a free-stall barn equipped with electronic recognition feeding system for each cow (Zhenghong Co., Shanghai, China). Cows were allowed free access to a total mixed ration consisting of 41% roughage and 59% concentrate, and diet composition has been described previously . Rumen fluid was collected using oral stomach tubes before morning feeding as described by Shen et al. . The pH was immediately measured using a calibrated pH meter (Starter 300; Ohaus Instruments Co. Ltd., Nanjing, China). Samples were stored at − 80 °C until further processing. Two mL of rumen fluid was mixed with 20 μL of 25% orthophosphate acid and then centrifuged at 20,000 × g at 4 °C for 10 min. The supernatant was harvested to use for analysis of VFA by gas chromatograph (GC-2010, Shimadzu, Kyoto, Japan) . Mathematical model was employed to predict methane production by the equation of Ellis et al. : CH4 (MJ/d) = 3.23 (± 1.12) + 0.81 (± 0.086) × DMI (kg/d).
DNA extraction and metagenome sequencing
Microbial DNA was extracted from rumen fluid samples by the repeat bead-beating plus column method . The DNA concentration and purity were determined by using NanoDrop spectrophotometer (Thermo Fisher Scientific, Wilmington, DE, USA). The paired-end library was constructed using TruSeq™ DNA Sample Prep Kit (Illumina, San Diego, CA, USA) and paired-end sequencing was performed on the Illumina HiSeq 4000 platform (150 bp paired-end).
The 3’ and 5’ends were then stripped with SeqPrep, and low-quality reads (length < 50 bp, quality value < 20, or N bases) were removed with Sickle (version 1.33). The reads were aligned to the bovine genome (bosTau8 3.7) with BWA (Version 0.7.9a), and any hit that associated with a read with a corresponding reads was removed . Reads from metagenomic trimmed data were assembled by multiple megahit using Megahit (Version 1.1.2) . Open reading frames (ORFs) in each metagenomic sample were predicted with MetaGene, and the predicted ORFs with lengths > 100 bp were extracted and translated to amino acid sequences using the NCBI translation table . The assembled contigs were clustered using CD-HIT, and non-redundant gene catalogue were constructed based on the sequences from gene sets with a 95% identity (90% coverage) . After quality control, SOAPaligner (Version 2.2.1) was used to map the reads to predicted genes with 95% identity, and the gene abundance in each sample was evaluated .
Taxonomic and functional annotation of rumen microbiota
We used BLASTP (version 2.2.28 +) for taxonomic annotation , by aligning the nonredundant gene catalogues against the NCBI NR database . The PCoA based on Bray–Curtis dissimilarity matrices was applied to visualize the taxonomic composition at species level. Microbial taxa with a relative abundance greater than 0.1% in at least 50% of animals within each efficiency group were used for downstream analysis. The KEGG pathway annotation was performed with Diamond against the KEGG database, with an e-value of 1e-5 . Carbohydrate-active enzyme annotation was made with hmmscan (Version 3.1b2). The abundances of pathways, KEGG Orthology (KO), KEGG enzymes, and CAZymes were normalized into counts per million reads (cpm) for downstream analysis.
Correlation analyses between rumen microbiota and pathways
Correlations between significantly different species and significantly different metabolic pathways were analysed using Spearman’s rank correlation, and significant correlation coefficients with R > 0.5 or < -0.5 (P < 0.05) were determined to generate the correlation network. The correlation networks were visualized by Gephi software (version 0.9.2).
Lactation performance and rumen VFAs concentrations were analysed by Student’s t test. Taxonomic and functional data were analysed on the online platform of the Majorbio Cloud Platform . Differential abundances of phylum, family, genus, and CAZymes were tested by Wilcoxon test using stats R package in R software (version 3.3.1). Significance was considered at P ≤ 0.05, and a tendency was defined as 0.05 < P ≤ 0.10.
Availability of data and materials
The Illumina sequencing raw data for our samples have been deposited in the NCBI Sequence Read Archive (SRA) under accession numbers: PRJNA597489.
Shalloo L, Kennedy J, Wallace M, Rath M, Dillon P. The economic impact of cow genetic potential for milk production and concentrate supplementation level on the profitability of pasture based systems under different EU milk quota scenarios. J Agric Sci. 2004;142(3):357–69.
Li B, VanRaden PM, Guduk E, O’Connell JR, Null DJ, Connor EE, et al. Genomic prediction of residual feed intake in US Holstein dairy cattle. J Dairy Sci. 2020;103(3):2477–86.
Connor EE, Hutchison JL, Norman HD, Olson KM, Van Tassell CP, Leith JM, et al. Use of residual feed intake in Holsteins during early lactation shows potential to improve feed efficiency through genetic selection1. J Anim Sci. 2013;91(8):3978–88.
Basarab JA, Price MA, Aalhus JL, Okine EK, Snelling WM, Lyle KL. Residual feed intake and body composition in young growing cattle. Can J Anim Sci. 2003;83(2):189–204.
Fischer A, Delagarde R, Faverdin P. Identification of biological traits associated with differences in residual energy intake among lactating Holstein cows. J Dairy Sci. 2018;101(5):4193–211.
Connor EE. Invited review: improving feed efficiency in dairy production: challenges and possibilities. Animal. 2015;9(3):395–408.
Kong RS, Liang G, Chen Y, Stothard P, Guan LL. Transcriptome profiling of the rumen epithelium of beef cattle differing in residual feed intake. BMC Genomics. 2016;17:592.
McDonnell RP, Hart KJ, Boland TM, Kelly AK, McGee M, Kenny DA. Effect of divergence in phenotypic residual feed intake on methane emissions, ruminal fermentation, and apparent whole-tract digestibility of beef heifers across three contrasting diets. J Anim Sci. 2016;94(3):1179–93.
Bergman E. Energy contributions of volatile fatty acids from the gastrointestinal tract in various species. Physiol Rev. 1990;70(2):567–90.
Guan LL, Nkrumah JD, Basarab JA, Moore SS. Linkage of microbial ecology to phenotype: correlation of rumen microbial ecology to cattle’s feed efficiency. FEMS Microbiol Lett. 2008;288(1):85–91.
Jami E, White BA, Mizrahi I. Potential role of the bovine rumen microbiome in modulating milk composition and feed efficiency. PLoS ONE. 2014;9(1):e85423.
McGover E, McGee M, Byrne CJ, Kenny DA, Kelly AK, Waters SM. Investigation into the effect of divergent feed efficiency phenotype on the bovine rumen microbiota across diet and breed. Sci Rep. 2020;10(1):1–10.
Shabat SK, Sasson G, Doron-Faigenboim A, Durman T, Yaacoby S, Berg Miller ME, et al. Specific microbiome-dependent mechanisms underlie the energy harvest efficiency of ruminants. ISME J. 2016;10(12):2958–72.
Xue MY, Sun HZ, Wu XH, Liu JX, Guan LL. Multi-omics reveals that the rumen microbiome and its metabolome together with the host metabolome contribute to individualized dairy cow performance. Microbiome. 2020;8:64.
Liu E, VandeHaar MJ. Relationship of residual feed intake and protein efficiency in lactating cows fed high- or low-protein diets. J Dairy Sci. 2020;103(4):3177–90.
Potts SB, Boerman JP, Lock AL, Allen MS, VandeHaar MJ. Residual feed intake is repeatable for lactating Holstein dairy cows fed high and low starch diets. J Dairy Sci. 2015;98(7):4735–47.
Li F, Hitch TCA, Chen Y, Creevey CJ, Guan LL. Comparative metagenomic and metatranscriptomic analyses reveal the breed effect on the rumen microbiome and its associations with feed efficiency in beef cattle. Microbiome. 2019;7:6.
Delgado B, Bach A, Guasch I, González C, Elcoso G, Pryce JE, et al. Whole rumen metagenome sequencing allows classifying and predicting feed efficiency and intake levels in cattle. Sci Rep. 2019;9:11.
Lin P-Y, Whang L-M, Wu Y-R, Ren W-J, Hsiao C-J, Li S-L, et al. Biological hydrogen production of the genus Clostridium: metabolic study and mathematical model simulation. Int J Hydrogen Energ. 2007;32(12):1728–35.
Moon CD, Pacheco DM, Kelly WJ, Leahy SC, Li D, Kopecny J, et al. Reclassification of Clostridium proteoclasticum as Butyrivibrio proteoclasticus comb. nov., a butyrate-producing ruminal bacterium. Int J Syst Evol Microbiol. 2008;58(Pt 9):2041–5.
Palevich N, Kelly WJ, Leahy SC, Denman S, Altermann E, Rakonjac J, et al. Comparative genomics of rumen butyrivibrio spp. uncovers a continuum of polysaccharide-degrading capabilities. Appl Environ Microbiol. 2019;86(1):e01993-e2019.
Berni Canani R, Sangwan N, Stefka AT, Nocerino R, Paparo L, Aitoro R, et al. Lactobacillus rhamnosus GG-supplemented formula expands butyrate-producing bacterial strains in food allergic infants. ISME J. 2016;10(3):742–50.
Zhou M, Hernandez-Sanabria E, Guan LL. Assessment of the microbial ecology of ruminal methanogens in cattle with different feed efficiencies. Appl Environ Microbiol. 2009;75(20):6524–33.
Jeyanathan J, Martin C, Morgavi D. The use of direct-fed microbials for mitigation of ruminant methane emissions: a review. Animal. 2014;8:250–61.
Moss AR, Jouany J-P, Newbold J. Methane production by ruminants: its contribution to global warming. Ann Zootech. 2000;49:231–53.
Li F, Guan LL. Metatranscriptomic profiling reveals linkages between the active rumen microbiome and feed efficiency in beef cattle. Appl Environ Microbiol. 2017;83(9):e00061-e117.
Shen J, Zheng L, Chen X, Han X, Cao Y, Yao J. Metagenomic analyses of microbial and carbohydrate-active enzymes in the rumen of dairy goats fed different rumen degradable starch. Front Microbiol. 2020;11:1003.
Meehan CJ, Beiko RG. A phylogenomic view of ecological specialization in the lachnospiraceae, a family of digestive tract-associated bacteria. Genome Biol Evol. 2014;6:703–13.
Vacca M, Celano G, Calabrese FM, Portincasa P, Gobbetti M, De Angelis M. The controversial role of human gut lachnospiraceae. Microorganisms. 2020;8:573.
Pathak A. Various factors affecting microbial protein synthesis in the rumen. Vet World. 2008;1(6):186.
Won MY, Oyama LB, Courtney SJ, Creevey CJ, Huws SA. Can rumen bacteria communicate to each other? Microbiome. 2020;8(1):23.
De Keersmaecker SC, Sonck K, Vanderleyden J. Let LuxS speak up in AI-2 signaling. Trends Microbiol. 2006;14(3):114–9.
Xie YY, Miao C, Lu Y, Sun HZ, Liu JX. Nitrogen metabolism and mammary gland amino acid utilization in lactating dairy cows with different residual feed intake. Anim Biosci. 2021;34(10):1600–6.
Shen JS, Chai Z, Song LJ, Liu JX, Wu YM. Insertion depth of oral stomach tubes may affect the fermentation parameters of ruminal fluid collected in dairy cows. J Dairy Sci. 2012;95(10):5978–84.
Xue MY, Sun HZ, Wu XH, Guan LL, Liu JX. Assessment of rumen microbiota from a large dairy cattle cohort reveals the pan and core bacteriomes contributing to varied phenotypes. Appl Environ Microbiol. 2018;84(19):e00970-e1018.
Ellis JL, Kebreab E, Odongo NE, McBride BW, Okine EK, France J. Prediction of methane production from dairy and beef cattle. J Dairy Sci. 2007;90:3456–66.
Yu Z, Morrison M. Improved extraction of PCR-quality community DNA from digesta and fecal samples. Biotech. 2004;36(5):808–12.
Li H, Durbin R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009;25(14):1754–60.
Li D, Liu C-M, Luo R, Sadakane K, Lam T-W. MEGAHIT: an ultra-fast single-node solution for large and complex metagenomics assembly via succinct de Bruijn graph. Bioinformatics. 2015;31(10):1674–6.
Noguchi H, Park J, Takagi T. MetaGene: prokaryotic gene finding from environmental genome shotgun sequences. Nucleic Acids Res. 2006;34(19):5623–30.
Fu L, Niu B, Zhu Z, Wu S, Li W. CD-HIT accelerated for clustering the next-generation sequencing data. Bioinformatics. 2012;28(23):3150–2.
Li R, Yu C, Li Y, Lam T-W, Yiu S-M, Kristiansen K, et al. SOAP2: an improved ultrafast tool for short read alignment. Bioinformatics. 2009;25(15):1966–2017.
Altschul SF, Madden TL, Schafer AA, Zhang J, Zhang Z, Miller W, et al. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997;25(17):3389–402.
Pruitt KD, Tatusova T, Maglott DR. NCBI reference sequences (RefSeq): a curated non-redundant sequence database of genomes, transcripts and proteins. Nucleic Acids Res. 2007;35(Suppl 1):61–5.
Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28(1):27–30.
Berlemont R, Martiny AC. Genomic potential for polysaccharide deconstruction in bacteria. Appl Environ Microbiol. 2015;81(4):1513–9.
The authors thank the staff of the Hangzhou Zhengxing Animal Industry Company (Hangzhou, China) for their animal raising and sampling assistance. We also acknowledge the members of the Institute of Dairy Science of Zhejiang University (Hangzhou, China) for their help in the field sampling and data analysis.
This research was supported by grants from the National Natural Science Foundation of China (31872380) and the China Agricultural (Dairy) Research System (CARS-36, Beijing).
Ethics approval and consent to participate
Animal care and experimental procedures were approved by the Animal Care Committee of Zhejiang University (Hangzhou, China, No. 12410) and were under the university’s guidelines for animal research.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
. Comparison of lactation performance between the cows with high (HRFI) and low residual feed intake (LRFI). Table S2. Summary of sequence data generated from rumen samples of 9 high (H) and 9 low (L) residual feed intake cows. Table S3. Permutational multivariate analysis of variance (PERMANOVA) for the Bray–Curtis dissimilarity matrices for microbial taxonomy between the cows with high (HRFI) and low residual feed intake (LRFI). Figure S1. Rumen microbial composition based on the domain level taxonomy. Figure S2. Microbial profiles of the cows with high (HRFI) and low residual feed intake (LRFI). Figure S3. Microbial metabolic pathways and carbohydrate-active enzymes (CAZymes) of the cows with high (HRFI) and low residual feed intake (LRFI). Figure S4. Distribution of residual feed intake in high (HRFI) and low residual feed intake (LRFI) cows.
About this article
Cite this article
Xie, Y., Sun, H., Xue, M. et al. Metagenomics reveals differences in microbial composition and metabolic functions in the rumen of dairy cows with different residual feed intake. anim microbiome 4, 19 (2022). https://doi.org/10.1186/s42523-022-00170-3
- Dairy cattle
- Residual feed intake
- Rumen microbiota