ORIGINAL RESEARCH article
Identification of Rumen Microbial Genes Involved in Pathways Linked to Appetite, Growth, and Feed Conversion Efficiency in Cattle
- 1Beef and Sheep Research Centre, Future Farming Systems Group, Scotland’s Rural College, Edinburgh, United Kingdom
- 2Division of Genetics and Genomics, The Roslin Institute and R(D)SVS, University of Edinburgh, Edinburgh, United Kingdom
- 3The Rowett Institute, University of Aberdeen, Aberdeen, United Kingdom
The rumen microbiome is essential for the biological processes involved in the conversion of feed into nutrients that can be utilized by the host animal. In the present research, the influence of the rumen microbiome on feed conversion efficiency, growth rate, and appetite of beef cattle was investigated using metagenomic data. Our aim was to explore the associations between microbial genes and functional pathways, to shed light on the influence of bacterial enzyme expression on host phenotypes. Two groups of cattle were selected on the basis of their high and low feed conversion ratio. Microbial DNA was extracted from rumen samples, and the relative abundances of microbial genes were determined via shotgun metagenomic sequencing. Using partial least squares analyses, we identified sets of 20, 14, 17, and 18 microbial genes whose relative abundances explained 63, 65, 66, and 73% of the variation of feed conversion efficiency, average daily weight gain, residual feed intake, and daily feed intake, respectively. The microbial genes associated with each of these traits were mostly different, but highly correlated traits such as feed conversion ratio and growth rate showed some overlapping genes. Consistent with this result, distinct clusters of a coabundance network were enriched with microbial genes identified to be related with feed conversion ratio and growth rate or daily feed intake and residual feed intake. Microbial genes encoding for proteins related to cell wall biosynthesis, hemicellulose, and cellulose degradation and host–microbiome crosstalk (e.g., aguA, ptb, K01188, and murD) were associated with feed conversion ratio and/or average daily gain. Genes related to vitamin B12 biosynthesis, environmental information processing, and bacterial mobility (e.g., cobD, tolC, and fliN) were associated with residual feed intake and/or daily feed intake. This research highlights the association of the microbiome with feed conversion processes, influencing growth rate and appetite, and it emphasizes the opportunity to use relative abundances of microbial genes in the prediction of these performance traits, with potential implementation in animal breeding programs and dietary interventions.
The global population is expected to reach 9.8 billion by 2050 (United Nations–Department of Economic and Social Affairs/Population Division, 2017), resulting in an escalation of the global demand for food and of the need for economically and environmentally sustainable livestock production systems (Godfray et al., 2010; Gerber et al., 2013). A large portion of livestock production is based on ruminants. In 2017, the EU-28 had a population of 88 million bovine animals, including cattle and water buffalo (Eurostat, 2018). Ruminants are particularly interesting due to their ability to convert human-indigestible plant biomass into high-quality products for human consumption such as meat and milk. Ruminants live in a symbiotic relationship with their rumen microbiota (comprising bacteria, protozoa, fungi, and archaea), which produce enzymes able to digest their food by breaking down complex polysaccharides of the plant biomass into volatile fatty acids (VFA), microbial proteins, and vitamins (Russell and Hespell, 1981; Bergman, 1990; Van Soest, 1994). Thus, the rumen microbiota fermentation profile has a significant influence on the feed conversion efficiency of the host (Russell, 2001; Li et al., 2009; Hernandez-Sanabria et al., 2011; Jami et al., 2014; Sasson et al., 2017; Meale et al., 2018) and is accountable for up to 70% of the host’s daily energy requirements (Bergman, 1990).
In beef cattle production systems, expenses associated with feed account for up to 75% of the total production costs (Moran, 2005a; Nielsen et al., 2013), which makes the improvement of feed conversion efficiency very economically compelling. There is consequently great interest in understanding the host–microbial symbiotic relationships responsible for the conversion of feed into energy, protein, and vitamins usable by the host animal, but the mechanisms and degree to which the rumen microbiome impacts on animal production, health, and efficiency remain undercharacterized (Brulc et al., 2009; Creevey et al., 2014). Although the rumen harbors a core microbiome (Jami and Mizrahi, 2012; Henderson et al., 2015), in agreement with studies performed in the human gastrointestinal tract (Tap et al., 2009; Qin et al., 2010), the structure, and composition of the rumen microbiome varies within and between animals with differing performance traits. For example, in lactating dairy cattle, the increased methane yield during late lactation in comparison to early lactation within the same individual was found to be associated with significant changes in the ruminal microbial community structure (Lyons et al., 2018); Myer et al. (2015) showed different relative abundances of some microbial taxa and operational taxonomic units in animals with different average daily gain (ADG); Shabat et al. (2016) focused on residual feed intake (RFI) to demonstrate that highly efficient animals had a less diverse microbiota, being dominated by specific taxa and microbial genes which were involved in simpler metabolic pathway networks when compared to their less efficient counterparts. Other authors have reported that the rumen microbiome varies more between animals than within animals, proposing that the host itself and its physiological parameters have a significant influence on its own rumen microbiome (Li et al., 2009) and, therefore, on the efficiency of feed conversion into energy. In a mouse study, Benson et al. (2010) found that there is a well-defined portion of the gut microbiota that is subject to host genetic control, proposing it to be regarded as a host trait, rather than an environmental trait affecting the host. In agreement, in a beef cattle study, Roehe et al. (2016) confirmed the host genetic influence on the rumen bacterial composition using a genetic model based on sire progeny groups. The differences between sire progeny groups in methane emissions were in some cases larger than the differences found between diets differing largely in plant fiber content, suggesting a substantial host genetic influence on the microbial communities.
Selecting animals for breeding based on their ability to harvest energy from feed, together with nutritional interventions, could be the basis for an effective strategy to produce faster growing and more efficient animals (Gerber et al., 2013; Scollan et al., 2018). Given that the host has influence over the ruminal microbiome, which impacts the animals’ feed conversion efficiency, this selection may be further improved by the inclusion of rumen metagenomic information into predictive models, as previously suggested by Ross et al. (2013). Feed conversion efficiency is very often estimated by either feed conversion ratio (FCR) or RFI; the latter is independent of growth and maturity patterns and is expected to be more sensitive and precise in measurements of feed utilization (Arthur and Herd, 2008). The use of microbial genes as proxies for feed conversion efficiency traits may be much more cost effective, rapid, and less labor intensive than their recording (Ross et al., 2013; Roehe et al., 2016). Our earlier research was the first proposing that the inclusion of relative abundance of microbial genes as proxies for FCR may be favorable, allowing their use as selection criteria for breeding animals, by identifying 49 microbial genes that explained 88.3% of the variation observed in FCR (Roehe et al., 2016). To our knowledge, no other studies have focused on the relationship between microbial gene abundances and RFI, daily feed intake (DFI), and ADG, which highlights the importance and novelty of the present work.
This study aimed at validating whether rumen microbial gene abundances are suitable proxies for feed conversion efficiency traits such as FCR; the analysis was further extended by focusing on RFI. Based on the previous evidence of strong interactions between the rumen microbiome and the host animal with consequences for feed conversion efficiency (Guan et al., 2008; Roehe et al., 2016; Shabat et al., 2016), we hypothesized that microbial gene abundances are linked to the animals’ appetite and, consequently, to feed intake. A further aim of this research was to gain insight into the association of growth rate with the microbial gene abundances. Building on this, we aimed at better understanding the rumen microbial functional network associated with feed conversion efficiency and its component traits. This research will improve on the current knowledge about the impact of the rumen microbiome on appetite, growth, and efficiency of feed conversion processes.
Materials and Methods
This study was conducted at the Beef and Sheep Research Centre, SRUC, UK. The study was carried out in accordance with the requirements of the UK Animals (Scientific Procedures) Act 1986. The protocol was approved by the Animal Experiment Committee of SRUC. All standard biosecurity and institutional safety procedures were applied during the animal experiment and the laboratory analysis.
Animals, Adaptation Period, and Measurement of Traits
Two experiments were carried out to determine the effect of nitrate or lipid additives within different basal diets on methane emissions from beef cattle. The first experiment was conducted in 2013, and it consisted of a 2 × 2 × 3 factorial design including 84 steers of two breed types (crossbreed Charolais, CHx and Luing); two basal diets, forage (FOR) and concentrate (CONC), which consisted respectively of ratios of 520:480 and 84:916 forage to concentrate (g/kg dry matter); and three treatments, nitrate and lipid feed additives, as well as the control. From these animals, 24 animals were selected with extreme high and low FCR values within breed type and basal diet (two animals per feed additive and control). More details related to this experiment can be found in Duthie et al. (2015) and Troy et al. (2015). The second experiment was a 2 × 4 factorial design experiment, conducted in 2014, involving 80 animals. There were two breed types—40 crossbred Limousin (LIMx) and 40 crossbred Aberdeen Angus (AAx)—which were subject to a balanced design consisting of four dietary treatments using one basal diet (550:450 forage to concentrate ratio g/kg dry matter, FOR) and testing the effects of feed additives nitrate, lipid, or their combination in comparison to the control on methane output. Full details of the experiment are presented in Duthie et al. (2017). From this experiment, 18 animals were selected within each combination of breed type and diet: nine for the high FCR group and nine for the low FCR group. DFI was assessed by measuring dry matter intake (DMI, kg/day), which was recorded in both experiments using electronic feeding equipment (HOKO, Insentec, Marknesse, The Netherlands). Body weight (BW) was measured weekly using a calibrated weight scale (before fresh feed was offered). Growth was modeled by linear regression of BW against test date to obtain ADG, mid-test BW, and mid-test metabolic BW (MBW = BW0.75). FCR was calculated as average DMI (kg/day) divided by ADG. RFI was estimated as deviation of actual DMI (kg/day) from DMI predicted based on linear regression of actual DMI on ADG, mid-MBW, and fat depth at 12th/13th rib at the end of the 56-day test (Duthie et al., 2015; Troy et al., 2015; Duthie et al., 2017).
A flowchart summarizing the methods for generation of data and subsequent statistical analyses is presented in Figure 1.
Figure 1 Flowchart summarizing methods for generation of data and their statistical analyses: This flowchart summarizes how the data were generated and which statistical analyses were used to identify the associations between gene abundances and performance traits of animals to understand the rumen microbial functional pathways associated with these traits. KEGG, Kyoto Encyclopedia of Genes and Genomes; FCR, feed conversion ratio; ADG, average daily gain; RFI, residual feed intake; DFI, daily feed intake; PLS, partial least squares; LDA, linear discriminant analysis.
Sampling of Rumen Digesta and Whole Metagenomic Sequencing
As described in Duthie et al. (2015) and Auffret et al. (2017), animals from both experiments were slaughtered in a commercial abattoir where two samples of rumen digesta (∼50 ml) were collected immediately after the rumen was opened to be drained. The slaughter house sample collection process results in well-mixed samples of rumen contents. DNA was extracted from the rumen samples of 42 animals following the methodology described in Rooke et al. (2014). Illumina TruSeq libraries were prepared from genomic DNA and sequenced on Illumina HiSeq systems 4000 by Edinburgh Genomics (Edinburgh, UK). Paired-end reads (2 × 150 bp) were generated, resulting in between 8 and 15 GB per sample (between 40 and 73 million paired reads). The raw data can be downloaded from the European Nucleotide Archive under accession PRJEB21624.
Identification of the Rumen Microbial Gene Abundances
Bioinformatics analysis for identification of rumen microbial genes was carried out as previously described by Wallace et al. (2015). Briefly, to measure the abundance of known functional microbial genes in the rumen samples, reads from whole metagenome sequencing were aligned to the Kyoto Encyclopedia of Genes and Genomes (KEGG) database (Kanehisa and Goto, 2000) using Novoalign (www.novocraft.com). Parameters were adjusted such that all hits were reported that were equal in quality to the best hit for each read and allowing up to a 10% mismatch across the fragment. The KEGG Orthologue groups (KO) of all hits that were equal to the best hit were examined. If we were unable to resolve the read to a single KO, the read was ignored; otherwise, the read was assigned to the unique KO. Read counts were summed and normalized to the total number of hits. This mapping of the whole metagenomic data to the KEGG database resulted in a dataset comprising of 4,966 KEGG genes. Microbial genes were removed from the dataset when they were absent from three or more animals and when the mean relative abundance was lower than 0.001%, leaving 1,692 microbial genes for further analyses.
For each of the 1,692 microbial genes, a linear model was fitted, including as fixed effects a combined class variable of breed, diet, and year of experiment (six levels) and the FCR groups (high FCR, FCR-H and low FCR, FCR-L) using the lm() function in R version 3.4.2. The microbial genes which resulted in P ≥ 0.1 for the differences in FCR groups were not considered in the partial least squares analyses (PLS, SAS version 9.3 for Windows, SAS Institute Inc., Cary, NC, USA) to avoid excessive noise of microbial genes uncorrelated to the traits of interest. In the linear model, FCR groups were replaced successively by ADG, RFI, and DFI as covariables to identify only potentially relevant microbial genes of these traits for further PLS analyses. In addition, genes with unknown function were removed from these datasets.
Microbial genes whose relative abundances were significantly associated to each trait in the linear models were analyzed using a sequential PLS-based methodology. First, PLS models were calculated in which the number of latent variables was determined by “leave-one-out” cross-validation, and genes with lower variable importance in projection (VIP) were removed. Second, the sets of genes created in the first step were evaluated by PLS models using three latent variables to determine the smaller set of genes leading to higher explained variation of both independent and dependent variables.
Each set of microbial genes identified in the PLS analyses as best predicting the trait was then used in a linear discriminant analysis (LDA), performed in R version 3.5.1 (2018-07-02) package MASS_7.3-51.4. In these analyses, the categories were for FCR those described previously as FCR-H and FCR-L; for all other traits, animals were classified as high or low, depending on their observations being higher or lower than the median (balanced for trial, breed, and diet).
The microbial genes identified to be significantly associated with each trait were submitted to an extensive review about their functionality based on databases such as KEGG (Kanehisa Laboratories, 2018), BioCyc (Karp et al., 2017), and UniProt (Bateman, 2019) and information from the literature.
The coabundances between microbial genes were investigated in a stepwise network analysis using the Graphia Professional software (Kajeka Ltd, Edinburgh; Freeman et al., 2007), in which nodes represent microbial genes and edges represent a correlation value above a defined value of r. In the first step, the correlation threshold of r = 0.45 was selected such that all microbial genes (n = 1,692) were included in the network. The microbial genes identified by PLS to be associated with a trait of interest were then located in the network. Clustering was performed using the Markov clustering method (MCL) available in Graphia Professional using the default settings (inflation, preinflation, and scheme values of 6). All clusters that held at least one microbial gene previously identified in the PLS analysis to be associated with a trait of interest were identified. These were incorporated into a new network generated at correlation threshold of r = 0.80 containing 1,135 microbial genes. MCL was then performed on this network, with inflation and preinflation values of 2 and scheme value of 6, reflecting the clustering structure suggested in the network itself. Analyses of enrichment of genes identified in the PLS as associated to each trait were performed on the clusters, and significance was assessed at P < 0.05.
Performance Traits Related With Feed Conversion Efficiency
The average FCR values observed for animals selected into FCR-H (inefficient) and FCR-L (efficient) groups differed significantly by 2.3 kg DFI/kg ADG (Figure 2). When comparing these two groups for other traits, the FCR-H group had significantly higher values of RFI (0.8 kg) and significantly lower ADG (0.39 kg); in the case of DFI, no significant difference was observed between the FCR groups.
Figure 2 Distribution of variation and range of performance traits: (A) feed conversion ratio, (B) average daily gain, (C) residual feed intake, and (D) daily feed intake within feed conversion ratio groups (high and low). The boxplots show the variation and range of each trait within each feed conversion ratio group. FCR, feed conversion ratio; AAx, crossbred Aberdeen Angus; CHx, crossbred Charolais; LIMx, crossbred Limousin.
ADG and FCR had a strong significant negative correlation of 0.80, suggesting that high growth rate is associated with efficient animals, using less feed per kilogram of weight gain. FCR and RFI were significantly positively correlated, but at a low level of 0.32. DFI was significantly correlated with RFI and ADG at high and moderate levels of 0.77 and 0.53, respectively.
Rumen Microbial Genes Associated With Feed Conversion Efficiency Traits
The PLS analyses identified sets of 20 and 14 microbial genes whose relative abundances explained 63.4 and 65.4% of the variation in FCR and ADG, respectively, and sets of 17 and 18 microbial genes whose relative abundances explained 65.6 and 72.9% of the variation in RFI and DFI, respectively, including the combined fixed effect of diet, breed, and year of experiment (Table 1). Without this combined fixed effect, the variances explained by microbial genes in FCR and ADG decreased to 54.2 and 61.4%, while in RFI and DFI, they decreased to 50.8 and 67.7%, respectively. A discriminant analysis between groups of high- and low-performing animals, using the set of microbial genes identified in the PLS analysis to best predict each trait, resulted in prediction accuracies of 90, 79, 86, and 86% for FCR, ADG, RFI, and DFI (Figure 3).
Table 1 Percentage of variation in each trait explained by the microbial genes identified in the partial least squares (PLS).
Figure 3 Linear discriminant analysis density plots: Microbial genes identified in the PLS analyses to be significantly associated with the trait were used in a linear discriminant analysis of high- and low-performing animals. The density plots represent the predicted categories for each trait. The accuracy value represents the percentage of animals that were correctly assigned to their category. FCR, feed conversion ratio; ADG, average daily gain; RFI, residual feed intake; DFI, daily feed intake; LD1, linear discriminant 1.
The Venn diagram presented in Figure 4 illustrates the overlap between the sets of genes identified for the prediction of each of the four traits. For the prediction of FCR and ADG, six microbial genes were simultaneously selected: UDP-N-acetylmuramoylalanine-D-glutamate ligase, glycine cleavage system H protein, translation initiation factor IF-1, N utilization substance protein A, DNA-binding protein HU-beta, and diphthamide synthase subunit dph2 (murD, gcvH, infA, nusA, hupB, and dph2, respectively). Three microbial genes were simultaneously selected for the prediction of traits RFI and DFI: glucose-1-phosphate cytidylyltransferase, CDP-glucose 4,6-dehydratase, and energy-converting hydrogenase B subunit D (rfbF, rfbG, and ehbD, respectively). The microbial genes identified for the prediction of more than one trait are highlighted in the shaded rows in Tables 2–5, in which a more detailed information about their function and importance for prediction is provided.
Figure 4 Overlap analysis of identified microbial genes: The image illustrates the number of microbial genes identified in the partial least squares analysis as fitted for prediction of each animal performance trait exclusively, and the number of microbial genes simultaneously predicted for multiple traits: six microbial genes were simultaneously identified for prediction of FCR (feed conversion ratio) and ADG (average daily gain), and three for both RFI (residual feed intake) and DFI (daily feed intake).
Based on the relative abundance of 1,135 microbial genes across rumen samples, a coabundance network was developed (Figure 5), and clusters were identified. The clustering pattern evidences the microbial genes that are more closely connected to microbial genes previously identified in the PLS analyses. The network cluster to which each microbial gene belongs to is presented in Tables 2–5. Cluster 2 was significantly enriched for microbial genes predicting DFI and RFI&DFI (RFI and/or DFI). Cluster 4 was enriched for microbial genes predicting RFI and RFI&DFI. Microbial genes simultaneously predicting FCR and ADG were enriched in clusters 20 and 21, while those predicting FCR&ADG (FCR and/or ADG) were enriched in clusters 21 and 25. ADG-predicting microbial genes were enriched in clusters 21 and 25, whereas FCR-predicting genes were only enriched in cluster 25. Other genes previously identified in the PLS analysis were scattered across the graph.
Figure 5 Correlation network analysis of metagenomic data: Each node represents a vector of relative abundances of each microbial gene in all 42 animals, and the edges represent a correlation between the microbial genes. A minimum correlation threshold of 0.80 was applied to the network. Different colors illustrate different clusters, which were calculated using MCL method (inflation: 2; preinflation: 2; scheme: 6). Clusters identified by numbers were found to be significantly (P < 0.05) enriched for microbial genes identified for the traits whose abbreviations are between brackets (FCR, feed conversion ratio; ADG, average daily gain; RFI, residual feed intake; DFI, daily feed intake; FCR&ADG, set including microbial genes identified for prediction of either FCR and/or ADG; RFI&DFI, set including microbial genes identified for prediction of RFI and/or DFI; FCR+ADG, set including microbial genes simultaneously identified for prediction of both traits FCR and ADG).
Most microbial genes identified exclusively for the prediction of FCR are related to carbohydrate metabolism and transport: fructuronate reductase, galactokinase, alpha-glucuronidase, beta-glucuronidase, beta-glucosidase, phosphate butyryltransferase P, UDP-N-acetylglucosamine acyltransferase, gluconate 5-dehydrogenase, and lactate permease (respectively uxuB, galK, aguA, uidA, K01188, ptb, lpxA, idnO, and lctP) were proportionally more abundant in efficient animals (lower FCR, Supplementary Figure S1A). The microbial gene lactoylglutathione lyase (glo1) is also associated with carbohydrate metabolism and identified for predicting FCR, but it had higher relative abundance in less efficient animals (higher FCR). Microbial genes galK and xylE (i.e., MFS transporter, SP family, xylose:H+ symporter) were both located in cluster 5, but this cluster was not significantly enriched for microbial genes associated to FCR. On the other hand, cluster 25 was enriched due to the presence of microbial genes uxuB and lpxA.
Microbial genes associated with amino acid metabolism and transport pathways were identified for the prediction of ADG and found to be relatively more abundant in animals with higher ADG (see Supplementary Figure S1B), e.g., aspartate-semialdehyde dehydrogenase and phenylacetate-CoA ligase (asd and paak, respectively). Some housekeeping genes were also identified for this set, including large subunit ribosomal protein L17 and L36, F-type H+-transporting ATPase subunit delta and FKBP-type peptidyl-prolyl cis-trans isomerase slyD (rplQ and rpmJ, atpH, and slyD). Genes rplQ, atpH, and slyD were relatively more abundant in animals with higher ADG, and rpmJ was relatively more abundant in animals with lower ADG. The microbial gene N-acetylmuramoyl-L-alanine amidase (amiABC) was identified for prediction of ADG, being negatively correlated with the trait.
All microbial genes simultaneously identified for predicting FCR and ADG showed a negative correlation to FCR and a positive correlation to ADG. These included housekeeping genes (infA, hupB, and dph2), a gene related to carbohydrate metabolism (gcvH), murD, which was associated with peptidoglycan metabolism and D-glutamine and D-glutamate metabolism, and nusA, associated with transcription regulation. Cluster 21 was enriched in ADG- and FCR&ADG-predicting microbial genes due to the presence of atpH, rplQ (ADG), and infA (FCR&ADG).
Five microbial genes identified for the prediction of RFI were associated with environmental sensing, bacterial chemotaxis, and motility: sensor kinase cheA, response regulator cheY, methyl accepting chemotaxis protein, flagellar motor switch protein fliN/fliY, and flagellar hook protein flgE (cheA, cheY, mcp, fliN, and flgE, respectively) were found to be relatively more abundant in more efficient animals, i.e., lower RFI. Other microbial genes associated with RFI are involved in the biosynthesis of cofactors and vitamins, particularly vitamin B12 production, for example, cobalt transport protein, threonine-phosphate decarboxylase, and precorrin-6Y C5,15-methyltransferase (decarboxylating), which correspond respectively to cbiN, cobD, and cobL (Supplementary Figure S1C). Finally, three genes that encode proteins related to carbohydrate transport and metabolism were relatively more abundant in more efficient animals (i.e., lower RFI): the simple sugar transport system permease protein, oxaloacetate decarboxylase, alpha subunit, and aldehyde:ferredoxin oxidoreductase (respectively ABC.SS.P, oadA, and aor). Cluster 4 was significantly enriched in microbial genes associated with RFI due to the presence of microbial genes cobD, cobL, mcp, and oadA, and serine-type D-Ala-D-Ala carboxypeptidase (penicillin-binding protein 5/6), inner membrane protein, and Cd2+/Zn2+-exporting ATPase (respectively, dacC, ybrG, and zntA).
The set of microbial genes identified for prediction of DFI included four microbial genes, proportionally more abundant in animals with higher DFI, which encoded proteins associated with environmental sensing, i.e., nitrogen regulatory protein P-II 1, outer membrane channel protein TolC, and preprotein translocase subunit YajC (glnB, tolC, and yajC, respectively). Nitrate reductase 1, alpha subunit (narG) was related to denitrification, releasing nitrite, and it was found to be relatively more abundant in animals with lower DFI (Supplementary Figure S1D). DNA-directed RNA polymerase subunit beta (rpoB, proportional higher abundance in animals with lower DFI), ribosomal large subunit pseudouridine synthase B, exodeoxyribonuclease VII small subunit, ribonuclease III, N utilization substance protein B, and integration host factor subunit alpha (respectively rluB, xseB, rnc, nusB, and ihfA, proportionally more abundant in animals with higher DFI) are housekeeping genes identified in this work for the prediction of DFI. Cluster 2 was significantly enriched with microbial genes associated with DFI due to the presence of glnB, infA, mrdA, nusB, rdgB, rluB, tolC, and xseB.
RFI- and DFI-predicting genes include glucose-1-phosphate cytidylyltransferase, CDP-glucose 4,6-dehydratase (respectively rfbF and rfbG, related to amino sugar and nucleotide sugar metabolism), and energy-converting hydrogenase B subunit D (ehbD, housekeeping). These three genes were proportionally more abundant in less efficient animals (higher RFI associated with increased DFI).
Rumen Microbial Gene Abundances Associated With Efficiency Traits
Our research indicates that there is a substantial link between rumen microbial gene abundances and appetite (measured as feed intake), growth rate, and feed conversion efficiency (Figure 6). The relative abundances of 20 and 17 microbial genes accounted for substantial variation (>60%) in FCR and RFI, respectively. The discriminant analyses of high- and low-performing animals indicated that accurate classification (>85% correct assignment of FCR and RFI categories) could be achieved using the microbial genes identified in the PLS for the prediction of the traits. Roehe et al. (2016) also found an association of microbial gene abundances with FCR, but their results were based on a smaller number of animals selected for their extreme values in methane emissions. In the present study, animals were selected based on their extreme FCR values, yielding a statistically more powerful estimate of this trait. Whereas FCR is calculated as a ratio between DFI and ADG and is therefore highly affected by growth rate and body composition, RFI is independent of these traits (Berry and Crowley, 2013). The low phenotypic correlation (r = 0.32) between FCR and RFI suggests that these traits capture substantially distinct characteristics.
Figure 6 Summary of microbial genes identified for the prediction of each trait: Traits are located in the four central boxes: FCR, feed conversion ratio; ADG, average daily gain; RFI, residual feed intake; DFI, daily feed intake. Solid lines represent positive correlations, and dotted lines represent negative correlations. Microbial genes are listed in the outside boxes, organized by general function, and each general function is represented by a different color.
For ADG and DFI, the relative abundances of 14 and 18 microbial genes, respectively, also explained substantial variation (>65%), and the discriminant analyses of high- and low-performing animals resulted in high prediction accuracies of 79 and 86%, respectively. These component traits were moderately correlated, agreeing with the report by Berry and Crowley (2013) of a large independent variation of feed intake and weight gain.
The animals’ appetite, feeding behaviour, and gastrointestinal motility (among other traits) are thought to be regulated by several mechanisms, including a communication between the rumen microbiome and the brain, through the gut–liver–brain axis (vagus nerve). This communication has been proposed to be mediated by multiple mechanisms, such as insulin/glucagon homeostasis, oxidation of acetyl coenzyme A, and release of VFA by the rumen microbiota (like propionate, associated with hypophagic behavior in ruminants, or butyrate and acetate, associated with motility of the gastrointestinal tract in monogastric animals; Sakata and Tamate, 1979; Cherbut, 2003; Oba and Allen, 2003; Arora et al., 2011; Maldini and Allen, 2018). Given the predictability of performance traits using relative abundances of rumen microbial genes observed in the present research (particularly that of DFI) and the high impact of the rumen microbiome on feed intake regulation (as discussed in the literature), we hypothesize that rumen microbial genes are closely involved in the metabolic pathways that regulate feed intake.
Differential Microbial Gene Sets Predicting Distinct Trait Complexes
The coabundance microbial gene network (Figure 5) identified two separate trait complexes. While microbial genes identified for the prediction of FCR were grouped with ADG-predicting genes, microbial genes identified for the prediction of RFI were grouped with DFI-predicting genes, as revealed by differential enrichment in separate clusters (Supplementary Figure S2). For example, beta-glucosidase is encoded by microbial genes bglX and K01188, which were associated to different traits (DFI and FCR, respectively). This type of differential clustering was previously observed for microbial genes associated with methane emissions and FCR by Roehe et al. (2016). The trait complexes associated with feed conversion efficiency were further evidenced when analyzing the overlapping genes identified for the prediction of each trait (Figure 4 and shaded rows in Tables 2–5), i.e., six microbial genes were identified for the prediction of both FCR and ADG and three genes for the prediction of both RFI and DFI. In agreement, strong correlations were observed for each pair of traits, as shown previously in the literature with the literature (Arthur and Herd, 2008; Herd et al., 2014). These results suggest that different microbial genes can be used to predict each trait. Furthermore, microbial genes overlapping for the prediction of more than one trait might be useful for the interpretation of biological processes explaining the correlation between phenotypes.
Metabolic Pathways of Microbial Genes Associated With Efficiency Traits
Our results indicate that most proteins encoded by microbial genes identified for the prediction of FCR were generally involved in carbohydrates metabolism and transport. For example, aguA and K01188 are involved in biomass conversion, through the degradation of hemicelluloses and lignocelluloses and lactate biosynthesis (Cairns and Esen, 2010; Lee et al., 2012; Michlmayr and Kneifel, 2014; Li, 2015). Microbial genes xylE, aguA, and uidA are involved in xylan degradation, the main component of hemicellulose (Lee et al., 2012; Fliegerova et al., 2015). Xylose needs to be taken up by a transporter (putatively associated with xylE) before it is metabolized, and it has been recognized as a rate-controlling step in bacterial metabolism (Chaillou and Pouwels, 1999). Furthermore, microbial genes such as uidA [previously identified by Roehe et al. (2016)], directly involved in carbohydrate metabolism pathways like pentose and glucuronate interconversions and galactose metabolism, are coupled with NAD or NADP oxidoreduction, important for regulating the flux of carbon and energy sources in microorganisms (Spaans et al., 2015). In addition, punA (i.e., purine-nucleoside phosphorylase) is involved in the metabolism of nucleotides, nicotinate and nicotinamide (vitamin B3), which also contain NAD and NADP, and is therefore important in carbohydrate, protein, and lipid metabolism reactions. Positive effects of vitamin B3 have been previously observed in healthy rumen microbiomes in beef and dairy cattle (Aschemann et al., 2012; Luo et al., 2017). Microbial genes uidA and punA were more abundant in efficient animals.
Proteins encoded by lctP, K01188, and ptb, involved in lactate transport and cellulose and butyrate metabolism, respectively, could be involved in host–microbiome crosstalk mechanisms in cattle due to their participation in metabolic pathways that involve the release of H+, such as lactate metabolism, potentially reducing microbial fiber-degrading activity and consequently slowing digestion and rumen emptying rate, causing a decrease in appetite (Moran, 2005b). Furthermore, beta-glucosidase is widely present in lactic acid bacteria and is thought to interact with the human host (Michlmayr and Kneifel, 2014). Butyrate has been shown in rats to directly activate the intestinal gluconeogenesis genes in enterocytes via an increase in cationic antimicrobial peptides (cAMP, De Vadder et al., 2014). In contrast, glo1 (more abundant in FCR-H) is involved in methylglyoxal degradation, which is a highly toxic substance that decreases bacterial cell viability, and is produced by bacteria when there is carbohydrate excess and nitrogen limitation (Russell, 1993). Therefore, glo1 is a strong candidate biomarker of rumen microbiome difference in less efficient animals (i.e., FCR-H).
The microbial gene with highest impact in prediction of ADG was amiABC, which is mainly involved in the peptidoglycan turnover through cleavage of glyosidic bonds and release of amino acids and cAMP resistance (Uehara and Park, 2008; Uehara et al., 2010). Some bacteria (mostly pathogenic) have evolved mechanisms of resistance, such as decreased affinity to cAMPs (Anaya-López et al., 2013), and the higher abundance of amiABC in animals with lower ADG may be indicative of higher abundance of pathogens, which can cause inflammatory response in the rumen potentially reducing nutrient use and absorption (Reynolds et al., 2017). Brown et al. (2003) demonstrated that acetate and propionate are agonists of the human receptors GPR43 and GPR41, and Hong et al. (2005) proposed that acetate and propionate induce lipid accumulation and inhibition of lipolysis through the GPR43 receptor in mice. These genes are also part of the bovine genome, where they mediate an inhibitory effect of acetate, propionate, and butyrate on cAMP signaling (Wang et al., 2009). This could indicate that, in less efficient animals (lower ADG), the lower amount of acetate, propionate, and butyrate may lead to decreased inhibition of lipolysis by the host, which potentially results in lower ADG. Alternatively, the lower amount of VFAs in these animals may lead to decreased inhibition of cAMP signaling and increased release of cAMPs by the host to the rumen. The cAMPs act primarily on organisms without effective resistance mechanisms, consequently increasing the relative abundance of cAMP-resisting organisms and of the microbial genes encoding for the resistance. Two other microbial genes identified in the present research are part of the cAMP resistance pathway—lpxA and tolC (associated with FCR and DFI, respectively). Although all three genes (amiABC, lpxA, and tolC) are part of the same pathway, they present opposite tendencies—while lpxA and tolC are proportionally highly abundant in animals with higher ADG and lower FCR, amiABC is relatively highly abundant in animals with lower ADG and higher FCR. The gene lpxA is related to lipid A integration in the cell wall, as a preventive measure against the hosts’ immune system, and tolC is involved in the efflux of antibiotics (Raetz et al., 2007; Zgurskaya et al., 2011). This could be indicative of the different cAMP resistance mechanisms evolved by bacterial organisms, which include modification of the cell external surface, efflux pumps, and biosynthesis and crosslinking of cell envelope components (Nizet, 2006).
The set of microbial genes associated with ADG included mostly housekeeping genes and genes related to amino acid metabolism and transport. Artegoitia et al. (2017) found a link between ruminal aromatic amino acids synthesis such as phenylalanine and high ADG in beef steers. For example, paak [previously mentioned by Kamke et al. (2016) related to sheep with high production of methane] and asd encode proteins that respectively catalyze phenylalanine and phenylacetate (related to aspartate degradation and biosynthesis of amino acids including threonine), with release of H+. In the current research, both of these genes were positively correlated to ADG, which is supported by the positive correlations between ADG and dry matter intake (DMI), between DMI and methane emissions, and between methane emissions and body weight measurements (weaning weight, yearling weight, and final weight), previously observed in cattle (Koots et al., 1994; Arthur et al., 2001; Herd et al., 2014).
Some housekeeping genes were simultaneously identified for the prediction of FCR and ADG, such as protein translation from diphthamide (dph2) or peptidoglycan biosynthesis (murD), both more abundant in efficient animals (higher ADG and lower FCR). The importance of diphthamide biosynthesis in archaea is not yet fully known (Narrowe et al., 2018). Microbial gene murD is related to the glutamate–glutamine cycle, an important appetite regulator in humans (Delgado, 2013), but in the present research, it was not associated to DFI.
Proteins encoded by microbial genes associated with RFI are mostly related to chemotaxis (cheA and cheY), detoxification (Cd2+/Zn2+-exporting ATPase, zntA), and vitamin B12 production (cbiN, cobD, and cobL). The negative correlation of microbial genes involved in chemotaxis and motility with RFI may suggest an increased microbial metabolism in efficient animals, derived from their ability to sense chemical gradients in their surrounding environment and to react accordingly, i.e., moving closer to nutrients (Rajagopala et al., 2007). Microbial gene zntA was also more abundant in efficient animals and plays a role in the homeostasis of transition metals (Cd2+, Zn2+), participating in functional pathways ranging from cellular respiration to gene expression (Fraústro da Silva and Williams, 2001). Finally, higher relative abundance of microbial genes involved in vitamin B12 production (cbiN, cobD, and cobL) was observed in more efficient animals. This essential cofactor needs to be taken up directly from the diet or to be made available for animal absorption by the rumen microbial organisms because it is not produced by eukaryotes (Warren et al., 2002). Furthermore, vitamin B12 has been previously associated with increased cobalt content on high-fiber diets and increased VFA, such as acetate (Beaudet et al., 2017), which may affect the animals’ appetite (Frost et al., 2014), in line with our observation of higher relative abundance of these genes in more efficient animals, i.e., animals with lower feed intake than expected.
The four most important microbial genes identified for the prediction of DFI included the three microbial genes also identified for prediction of RFI (rfbG, rfbF, and ehbD) and narG. Microbial genes rfbG and rfbF (VIP > 1.4) are part of the rfc region (Morona et al., 1994) and are related to nucleotide sugar metabolism, which is necessary for the production of microbial lipopolysaccharide (LPS). LPS is a major virulence factor of Gram-negative bacteria, particularly due to the O-antigen, paramount for host colonization and niche adaptation by bacterial organisms, due to its part in the protection from host immune response (Reeves, 1995; Samuel and Reeves, 2003; Geue et al., 2017). Both genes rfbG and rfbF showed a positive correlation to RFI and DFI, supporting our hypothesis that the use of energy to stimulate the innate immune system against pathogens increases DFI and reduces feed conversion efficiency as determined by RFI (Neal et al., 1991; Jing et al., 2014; Vigors et al., 2016). Other microbial genes positively correlated to DFI were found to be involved in resistance mechanisms, such as the penicillin-binding protein 2-encoding gene (mrdA), which belongs to the peptidoglycan and beta-lactam resistance metabolic pathways. These proteins are transpeptidases or carbopeptidases involved in peptidoglycan metabolism and have an important role against beta-lactam resistance (Zapun et al., 2008). The microbial gene myo-inositol-1-phosphate synthase (INO1) is related to antibiotic biosynthesis, including streptomycin. Microbial gene ehbD is a subunit of the energy-converting hydrogenase B, found in methanogens such as Methanococcus maripaludis. This microbial gene is important due to its role in autotrophic CO2 assimilation (Porat et al., 2006), having implications for microbial growth. Furthermore, narG, part of the narGHIJ operon, essential for some microorganisms to gather energy under anaerobic conditions by the reduction in nitrate to nitrite in a denitrification process (Blasco et al., 1990; Latham et al., 2016), was proportionally more abundant in animals with low DFI.
The microbial gene nusB (associated with DFI) is part of a set of nus genes, which also includes nusA (identified for prediction of FCR and ADG). Genes in the nus complex are involved in transcription termination and antitermination processes, such as Rho-dependent transcriptional termination (Torres et al., 2004), which is the regulatory mechanism involved in the efficient transcription of the tryptophan operon (Farnham et al., 1982; Kuroki et al., 1982; Prasch et al., 2009). The nus-complex microbial genes were found to be relatively more abundant in efficient animals. This association may be due to the influence of the nus genes, which extends from the ribosomal operons to the tryptophan operon and constitutes a good example of how termination and antitermination processes can control gene expression, occurring during RNA transcription, and potentially positively impacting bacterial growth and rumen fermentation processes.
Although microbial genes amiABC, tolC, glo1, rfbF, rfbG, and lpxA were identified in the present research for the prediction of different traits, all are associated with bacterial defense mechanisms either from other bacteria or from the host. The majority of these genes had higher abundance in less efficient animals. This suggests that the presence of either bacterial pathogens in the rumen or antibiotics produced as host immune responses might represent a significant energy sink, impairing feed conversion efficiency.
Further improvement of prediction of feed conversion traits using metagenomic information may be achieved through the integration of protein, enzyme, and pathway data from the Hungate collection (Seshadri et al., 2018) and the large rumen metagenomic reference dataset (Stewart et al., 2018).
The results presented here suggest that relative abundances of rumen microbial genes may be highly informative predictors of feed conversion efficiency, growth rate, and feed intake, which are labor intensive, time consuming, and expensive traits to record. Most microbial genes identified for the prediction of traits in this research were trait specific. Microbial genes related to cellulose and hemicellulose degradation, vitamin B12 synthesis, and amino acids metabolism were associated to enhanced feed conversion efficiency (FCR or RFI), while those involved in nucleotide sugars metabolism, pathogen LPS synthesis, cAMP resistance, and degradation of toxic compounds were associated with inefficient feed conversion. Furthermore, we identified specific microbial genes encoding proteins related to the crosstalk between the microbiome and the host cells, such as murD and amiABC, and associated to gene expression regulatory mechanisms, such as nusA and nusB. Thus, our results provide a deeper understanding of the potential influence of the rumen microbiome on the feed conversion efficiency of its host, highlighting specific enzymes involved in metabolic pathways that reflect the complex functional networks impacting the conversion of feed into animal products such as meat.
RR and MW conceived and designed the overall study, and JL, MA, TF, and RR conceived, designed, and executed the bioinformatics analysis. RS and MW carried out the bioinformatics to obtain the rumen microbial gene abundances. C-AD, TS, RD, and AW provided essential insight into feed conversion efficiency, rumen metabolism, nutrition, and microbiology. JL, MA, and RR wrote the initial draft, and subsequently, all authors contributed intellectually to the interpretation and presentation of the results.
The project was supported by grants from the Biotechnology and Biological Sciences Research Council (BBSRC BB/N01720X/1 and BB/N016742/1) and by the Scottish Government (RESAS Division) as part of the 2016–2021 commission. The research is based on data from experiments funded by the Scottish Government as part of the 2011–2016 commission, Agriculture and Horticulture Development Board (AHDB) Beef & Lamb, Quality Meat Scotland (QMS), and Department for Environment Food & Rural Affairs (Defra).
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
We thank Dr Irene Cabeza Luna, Laura Nicoll, Lesley Deans, and Claire Broadbent for the excellent technical support.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2019.00701/full#supplementary-material
Artegoitia, V. M., Foote, A. P., Lewis, R. M., Freetly, H. C. (2017). Rumen fluid metabolomics analysis associated with feed efficiency on crossbred steers. Sci. Rep. 7, 1–14. doi: 10.1038/s41598-017-02856-0
Arthur, P. F., Renand, G., Krauss, D. (2001). Genetic and phenotypic relationships among different measures of growth and feed efficiency in young Charolais bulls. Livest. Prod. Sci. 68, 131–139. doi: 10.1016/S0301-6226(00)00243-8
Aschemann, M., Lebzien, P., Hüther, L., Südekum, K.-H., Dänicke, S. (2012). Effect of niacin supplementation on rumen fermentation characteristics and nutrient flow at the duodenum in lactating dairy cows fed a diet with a negative rumen nitrogen balance. Arch. xAnim. Nutr. 66, 303–318. doi: 10.1080/1745039x.2012.697353
Auffret, M. D., Dewhurst, R. J., Duthie, C. A., Rooke, J. A., John Wallace, R., Freeman, T. C., et al. (2017). The rumen microbiome as a reservoir of antimicrobial resistance and pathogenicity genes is directly affected by diet in beef cattle. Microbiome 5, 159–169. doi: 10.1186/s40168-017-0378-z
Beaudet, V., Gervais, R., Graulet, B., Nozière, P., Doreau, M., Fanchone, A., et al. (2017). Effects of dietary nitrogen levels and carbohydrate sources on apparent ruminal synthesis of some B vitamins in dairy cows. J. Dairy Sci. 99, 2730–2739. doi: 10.3168/jds.2017-100-1-0849
Benson, A. K., Kelly, S. A., Legge, R., Ma, F., Low, S. J., Kim, J., et al. (2010). Individuality in gut microbiota composition is a complex polygenic trait shaped by multiple environmental and host genetic factors. Proc. Natl. Acad. Sci. 107, 18933–18938. doi: 10.1073/pnas.1007028107
Blasco, F., Iobbi, C., Ratouchniak, J., Bonnefoy, V., Chippaux, M. (1990). Nitrate reductases of Escherichia coli: sequence of the second nitrate reductase and comparison with that encoded by the narGHJI operon. MGG Mol. Gen. Genet. 222, 104–111.
Brown, A. J., Goldsworthy, S. M., Barnes, A. A., Eilert, M. M., Tcheang, L., Daniels, D., et al. (2003). The orphan G protein-coupled receptors GPR41 and GPR43 are activated by propionate and other short chain carboxylic acids. J. Biol. Chem. 278, 11312–11319. doi: 10.1074/jbc.m211609200
Brulc, J. M., Antonopoulos, D. A., Miller, M. E. B., Wilson, M. K., Yannarell, A. C., Dinsdale, E. A., et al. (2009). Gene-centric metagenomics of the fiber-adherent bovine rumen microbiome reveals forage specific glycoside hydrolases. Proc. Natl. Acad. Sci. 106, 1948–1953. doi: 10.1073/pnas.0806191105
Chaillou, P., Pouwels, P. H. (1999). Transport of D-xylose in Lactobacillus pentosus, Lactobacillus casei, and Lactobacillus plantarum: evidence for a mechanism of facilitated diffusion via the phosphoenolpyruvate: mannose phosphotransferase system. Society 181, 4768–4773.
De Vadder, F., Kovatcheva-Datchary, P., Goncalves, D., Vinera, J., Zitoun, C., Duchampt, A., et al. (2014). Microbiota-generated metabolites promote metabolic benefits via gut–brain neural circuits. Cell 156, 84–96. doi: 10.1016/j.cell.2013.12.016
Duthie, C.-A., Rooke, J. A., Troy, S., Hyslop, J. J., Ross, D. W., Waterhouse, A., et al. (2015). Impact of adding nitrate or increasing the lipid content of two contrasting diets on blood methaemoglobin and performance of two breeds of finishing beef steers. Animal 10, 786–795. doi: 10.1017/S1751731115002657
Duthie, C. A., Troy, S. M., Hyslop, J. J., Ross, D. W., Roehe, R., Rooke, J. A. (2017). The effect of dietary addition of nitrate or increase in lipid concentrations, alone or in combination, on performance and methane emissions of beef cattle. Animal 12, 280–287. doi: 10.1017/S175173111700146X
Eurostat. (2018). Eurostat online table Bovine, alive, EU-28. Luxembourg: Eurostat. Available at: http://appsso.eurostat.ec.europa.eu/nui/show.do?dataset=apro_mt_lscatl&lang=en (Accessed October 16, 2018).
Farnham, P. J., Greenblatt, J., Platt, T. (1982). Effects of NusA protein on transcription termination in the tryptophan operon of Escherichia coli. Cell 29, 945–951. doi: 10.1016/0092-8674(82)90457-3
Freeman, T. C., Goldovsky, L., Brosch, M., van Dongen, S., Maziere, P., Grocock, R. J., et al. (2007). Construction, visualisation, and clustering of transcription networks from microarray expression data. PLoS Comput. Biol. 3, e206 (11 pages). doi: 10.1371/journal.pcbi.0030206
Frost, G., Sleeth, M. L., Sahuri-Arisoylu, M., Lizarbe, B., Cerdan, S., Brody, L., et al. (2014). The short-chain fatty acid acetate reduces appetite via a central homeostatic mechanism. Nat. Commun. 5, 3611. doi: 10.1038/ncomms4611
Gerber, P. J., Steinfeld, H., Henderson, B., Mottet, A., Opio, C., Dijkman, J., et al. (2013). Tackling climate change through livestock - A global assessment of emissions and mitigation opportunities. Rome: Food and Agriculture Organization of the United Nations (FAO).
Geue, L., Menge, C., Eichhorn, I., Semmler, T., Wieler, L. H., Pickard, D., et al. (2017). Evidence for contemporary switching of the O-antigen gene cluster between shiga toxin-producing escherichia coli strains colonizing cattle. Front. Microbiol. 8, 424. doi: 10.3389/fmicb.2017.00424
Godfray, H. C. J., Beddington, J. R., Crute, I. R., Haddad, L., Lawrence, D., Muir, J. F., et al. (2010). Food security: the challenge of feeding 9 billion people. Science 327, 812–818. doi: 10.1126/science.1185383
Guan, L. L., Nkrumah, J. D., Basarab, J. A., Moore, S. S. (2008). Linkage of microbial ecology to phenotype: correlation of rumen microbial ecology to cattle’s feed efficiency. FEMS Microbiol. Lett. 288, 85–91. doi: 10.1111/j.1574-6968.2008.01343.x
Henderson, G., Cox, F., Ganesh, S., Jonker, A., Young, W., Global Rumen Census Collaborators G., et al. (2015). Rumen microbial community composition varies with diet and host, but a core microbiome is found across a wide geographical range. Sci. Rep. 5, 14567. doi: 10.1038/srep14567
Herd, R. M., Arthur, P. F., Donoghue, K. A., Bird, S. H., Bird-Gardiner, T., Hegarty, R. S. (2014). Measures of methane production and their phenotypic relationships with dry matter intake, growth, and body composition traits in beef cattle. J. Anim. Sci. 92, 5267–5274. doi: 10.2527/jas.2014-8273
Hernandez-Sanabria, E., Goonewardene, L. A., Wang, Z., Durunna, O. N. (2011). Impact of feed efficiency and diet on the adaptive variations in the bacterial community in the rumen fluid of cattle. Appl. Environ. Microbiol. 1, 1203–1214. doi: 10.1128/AEM.05114-11
Hong, Y.-H., Nishimura, Y., Hishikawa, D., Tsuzuki, H., Miyahara, H., Gotoh, C., et al. (2005). Acetate and propionate short chain fatty acids stimulate adipogenesis via GPCR43. Endocrinology 146, 5092–5099. doi: 10.1210/en.2005-0545
Jiang, X., Neal, B., Santiago, F., Lee, S. J., Reeves, P. R. (2006). Structure and sequence of the rfb (O antigen) gene cluster of Salmonella serovar typhimuriunfi (strain LT2). Mol. Microbiol. 5, 695–713. doi: 10.1111/j.1365-2958.1991.tb00741.x
Jing, L., Zhang, R., Liu, Y., Zhu, W., Mao, S. (2014). Intravenous lipopolysaccharide challenge alters ruminal bacterial microbiota and disrupts ruminal metabolism in dairy cattle. Br. J. Nutr. 112, 170–182. doi: 10.1017/S000711451400066X
Kamke, J., Kittelmann, S., Soni, P., Li, Y., Tavendale, M., Ganesh, S., et al. (2016). Rumen metagenome and metatranscriptome analyses of low methane yield sheep reveals a Sharpea-enriched microbiome characterised by lactic acid formation and utilisation. Microbiome 4, 56. doi: 10.1186/s40168-016-0201-2
Kanehisa Laboratories (2018). KEGG MAPPER. Available at: https://www.genome.jp/kegg/tool/map_pathway1.html [Accessed October 1, 2018]
Karp, P. D., Billington, R., Caspi, R., Fulcher, C. A., Latendresse, M., Kothari, A., et al. (2017). The BioCyc collection of microbial genomes and metabolic pathways. Brief. Bioinform. 201, 1–9. doi: 10.1093/bib/bbx085
Koots, K. R., Gibson, J. P., Wilton, J. W. (1994). Analyses of published genetic parameter estimates for beef production traits. 2. Phenotypic and genetic correlations. Anim. Breed. Abstr. 62, 825–853. Available at: http://agris.fao.org/agris-search/search.do?recordID=GB9508564
Kuroki, K., Ishii, S., Kano, Y., Miyashita, T., Nishi, K., Imamoto, F. (1982). Involvement of the nusA and nusB gene products in transcription of Escherichia coli tryptophan operon in vitro. MGG Mol. Gen. Genet. 185, 369–371. doi: 10.1007/BF00330816
Latham, E. A., Anderson, R. C., Pinchak, W. E., Nisbet, D. J. (2016). Insights on alterations to the rumen ecosystem by nitrate and nitrocompounds. Front. Microbiol. 7, 228. doi: 10.3389/fmicb.2016.00228
Li, M., Penner, G. B., Hernandez-Sanabria, E., Oba, M., Guan, L. L. (2009). Effects of sampling location and time, and host animal on assessment of bacterial diversity and fermentation parameters in the bovine rumen. J. Appl. Microbiol. 107, 1924–1934. doi: 10.1111/j.1365-2672.2009.04376.x
Luo, D., Gao, Y., Lu, Y., Qu, M., Xiong, X., Xu, L., et al. (2017). Niacin alters the ruminal microbial composition of cattle under high-concentrate condition. Anim. Nutr. 3, 180–185. doi: 10.1016/j.aninu.2017.04.005
Lyons, T., Bielak, A., Doyle, E., Kuhla, B. (2018). Variations in methane yield and microbial community profiles in the rumen of dairy cows as they pass through stages of first lactation. J. Dairy Sci. 101, 5102–5114. doi: 10.3168/jds.2017-14200
Maldini, G., Allen, M. S. (2018). Temporal effects of ruminal propionic acid infusion on feeding behavior of Holstein cows in the postpartum period. J. Dairy Sci. 101, 3077–3084 doi: 10.3168/jds.2017-13857
Meale, S. J., Auffret, M. D., Watson, M., Morgavi, D. P., Cantalapiedra-Hijar, G., Duthie, C. A., et al. (2018). Fat accretion measurements strengthen the relationship between feed conversion efficiency and nitrogen isotopic discrimination while rumen microbial genes contribute little. Sci. Rep. 8, 3854. doi: 10.1038/s41598-018-22103-4
Michlmayr, H., Kneifel, W. (2014). β-Glucosidase activities of lactic acid bacteria: mechanisms, impact on fermented food and human health. FEMS Microbiol. Lett. 352, 1–10. doi: 10.1111/1574-6968.12348
Moran, J. (2005a). “Chapter 5: How the rumen works,” in Tropical dairy farming: feeding management for small holder dairy farmers in the humid tropics. (Collingwood: Csiro Publishing), 41–49. doi: 10.1016/0022-1694(84)90033-7
Narrowe, A. B., Spang, A., Stairs, C. W., Caceres, E. F., Baker, B. J., Miller, C. S., et al. (2018). Complex evolutionary history of translation elongation factor 2 and diphthamide biosynthesis in archaea and parabasalids. Genome Biol. Evol. 10, 2380–2393. doi: 10.1093/gbe/evy154
Neal, B., Jiang, X.-M., Lee, S. J., Santiago, F., Reeves, P. R., Romana, L. K. (1991). Structure and sequence of the rfb (O antigen) gene cluster of Salmonella serovar typhimurium (strain LT2). Mol. Microbiol. 5, 695–713. doi: 10.1111/j.1365-2958.1991.tb00741.x
Nielsen, M. K., MacNeil, M. D., Dekkers, J. C. M., Crews, D. H., Jr., Rathje, T. A., Enns, R. M., et al. (2013). Review: Life-cycle, total-industry genetic improvement of feed efficiency in beef cattle: Blueprint for the Beef Improvement Federation. Prof. Anim. Sci. 29, 559–565. doi: 10.15232/s1080-7446(15)30285-0
Nizet, V. (2006). Antimicrobial peptide resistance mechanisms of human bacterial pathogens. Curr Issues Mol Biol 8, 11–26 Available at: https://www.ncbi.nlm.nih.gov/pubmed/16450883.
Oba, M., Allen, M. S. (2003). Intraruminal infusion of propionate alters feeding behavior and decreases energy intake of lactating dairy cows. J. Nutr 133, 1094–1099. doi: 10.3168/jds.S0022-0302(03)73889-2
Porat, I., Kim, W., Hendrickson, E. L., Xia, Q., Zhang, Y., Wang, T., et al. (2006). Disruption of the operon encoding Ehb hydrogenase limits anabolic CO2 assimilation in the archaeon Methanococcus maripaludis. J. Bacteriol. 188, 1373–1380. doi: 10.1128/JB.188.4.1373-1380.2006
Qin, J., Li, R., Raes, J., Arumugam, M., Burgdorf, S., Manichanh, C., et al., (2010). A human gut microbial gene catalog established by metagenomic sequencing. Nature 464, 59–65. doi: 10.1038/nature08821
Raetz, C. R. H., Reynolds, C. M., Trent, M. S., Bishop, R. E. (2007). Lipid A modification systems in Gram-negative bacteria. NIH Public Access 76, 295–329. doi: 10.1146/annurev.biochem.76.010307.145803
Reynolds, J. G., Foote, A. P., Freetly, H. C., Oliver, W. T., Lindholm-Perry, A. K. (2017). Relationships between inflammation- and immunity-related transcript abundance in the rumen and jejunum of beef steers with divergent average daily gain. Anim. Genet. 48, 447–449. doi: 10.1111/age.12546
Roehe, R., Dewhurst, R. J., Duthie, C. A., Rooke, J. A., McKain, N., Ross, D. W., et al. (2016). Bovine host genetic variation influences rumen microbial methane production with best selection criterion for low methane emitting and efficiently feed converting hosts based on metagenomic gene abundance. PLoS Genet. 12, e1005846. doi: 10.1371/journal.pgen.1005846
Rooke, J. A., Wallace, R. J., Duthie, C.-A., McKain, N., Souza, S. M., Hyslop, J. J., et al. (2014). Hydrogen and methane emissions from beef cattle and their rumen microbial community vary with diet, time after feeding and genotype. Br. J. Nutr. 112, 398–407. doi: 10.1017/s0007114514000932
Ross, E. M., Moate, P. J., Marett, L. C., Cocks, B. G., Hayes, B. J. (2013). Metagenomic predictions: from microbiome to complex health and environmental phenotypes in humans and cattle. PLoS One 8, e73056. doi: 10.1371/journal.pone.0073056
Russell, J. B. (1993). Glucose toxicity in Prevotella ruminicola: methylglyoxal accumulation and its effect on membrane physiology. Appl. Environ. Microbiol. 59, 2844–2850. Available at: https://www.ncbi.nlm.nih.gov/pubmed/8215358
Samuel, G., Reeves, P. (2003). Biosynthesis of O-antigens: genes and pathways involved in nucleotide sugar precursor synthesis and O-antigen assembly. Carbohydr. Res. 338, 2503–2519. doi: 10.1016/j.carres.2003.07.009
Sasson, G., Ben-Shabat, S., Seroussi, E., Doron-Faigenboim, A., Shterzer, N., Yaacoby, S., et al. (2017). Heritable bovine rumen bacteria are phylogenetically related and correlated with the cow’s capacity to harvest. Am. Soc. Microbiol. 8, e00703–17. doi: 10.1128/mBio.00703-17
Scollan, N. D., Morrisson, S., Popova, M., Attwood, G. T., Creevey, C. J., Tapio, I., et al. (2018). Addressing global ruminant agricultural challenges through understanding the rumen microbiome: past, present, and future. Front. Microbiol. 9, 2161. doi: 10.3389/fmicb.2018.02161
Seshadri, R., Leahy, S. C., Attwood, G. T., Teh, K. H., Lambie, S. C., Cookson, A. L., et al. (2018). Cultivation and sequencing of rumen microbiome members from the Hungate1000 Collection. Nat. Biotechnol. 36, 359–367. doi: 10.1038/nbt.4110
Shabat, S., Sasson, G., Doron-Faigenboim, A., Durman, T., Yaacoby, S., Berg Miller, M. E., et al. (2016). Specific microbiome-dependent mechanisms underlie the energy harvest efficiency of ruminants. ISME J. 10, 2958–2972. doi: 10.1038/ismej.2016.62
Stewart, R. D., Auffret, M. D., Warr, A., Wiser, A. H., Press, M. O., Langford, K. W., et al. (2018). Assembly of 913 microbial genomes from metagenomic sequencing of the cow rumen. Nat. Commun. 9, 870. doi:10.1038/s41467-018-03317-6
Tap, J., Mondot, S., Levenez, F., Pelletier, E., Caron, C., Furet, J. P., et al. (2009). Towards the human intestinal microbiota phylogenetic core. Environ. Microbiol. 11, 2574–2584. doi: 10.1111/j.1462-2920.2009.01982.x
Torres, M., Balada, J. M., Zellars, M., Squires, C., Squires, C. L. (2004). In vivo effect of NusB and NusG on rRNA transcription antitermination. J. Bacteriol. 186, 1304–1310. doi: 10.1128/JB.186.5.1304-1310.2004
Troy, S. M., Duthie, C. A., Hyslop, J. J., Roehe, R., Ross, D. W., Wallace, R. J., et al. (2015). Effectiveness of nitrate addition and increased oil content as methane mitigation strategies for beef cattle fed two contrasting basal diets. J. Anim. Sci. 93, 1815–1823. doi: 10.2527/jas.2014-8688
Uehara, T., Parzych, K. R., Dinh, T., Bernhardt, T. G. (2010). Daughter cell separation is controlled by cytokinetic ring-activated cell wall hydrolysis. EMBO J. 29, 1412–1422. doi: 10.1038/emboj.2010.36
United Nations–Department of Economic and Social Affairs/Population Division (2017). World Population Prospects: the 2017 revision. Key findings and advance tables, ed. United Nations New York: United Nations.
Vigors, S., O’Doherty, J. V., Kelly, A. K., O’Shea, C. J., Sweeney, T. (2016). The Effect of Divergence in feed efficiency on the intestinal microbiota and the intestinal immune response in both unchallenged and lipopolysaccharide challenged ileal and colonic explants. PLoS One 11, e0148145. doi: 10.1371/journal.pone.0148145
Wallace, R. J., Rooke, J. A., McKain, N., Duthie, C. A., Hyslop, J. J., Ross, D. W., et al. (2015). The rumen microbial metagenome associated with high methane production in cattle. BMC Genomics 16, 839. doi: 10.1186/s12864-015-2032-0
Wang, A., Gu, Z., Heid, B., Akers, R. M., Jiang, H. (2009). Identification and characterization of the bovine G protein-coupled receptor GPR41 and GPR43 genes. J. Dairy Sci. 92, 2696–2705. doi: 10.3168/jds.2009-2037
Keywords: feed conversion efficiency, appetite, metagenomics, rumen microbiome, microbial gene networks
Citation: Lima J, Auffret MD, Stewart RD, Dewhurst RJ, Duthie C-A, Snelling TJ, Walker AW, Freeman TC, Watson M and Roehe R (2019) Identification of Rumen Microbial Genes Involved in Pathways Linked to Appetite, Growth, and Feed Conversion Efficiency in Cattle. Front. Genet. 10:701. doi: 10.3389/fgene.2019.00701
Received: 07 February 2019; Accepted: 03 July 2019;
Published: 08 August 2019.
Edited by:Robert J. Schaefer, University of Minnesota Twin Cities, United States
Reviewed by:Robert W. Li, Agricultural Research Service (USDA), United States
Prakash G. Koringa, Anand Agricultural University, India
Copyright © 2019 Lima, Auffret, Stewart, Dewhurst, Duthie, Snelling, Walker, Freeman, Watson and Roehe. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
†ORCID: Tom Freeman, orcid.org/0000-0001-5235-8483