
==== Front
Anim Microbiome
Anim Microbiome
Animal Microbiome
2524-4671
BioMed Central London

337
10.1186/s42523-024-00337-0
Research
Relationship between the rumen microbiome and liver transcriptome in beef cattle divergent for feed efficiency
Keogh Kate kate.a.keogh@teagasc.ie

12
Kenny David A. 1
Alexandre Pamela A. 2
Waters Sinead M. 3
McGovern Emily 1
McGee Mark 4
Reverter Antonio 2
1 https://ror.org/03sx84n71 grid.6435.4 0000 0001 1512 9569 Animal and Bioscience Research Department, Teagasc, Animal & Grassland Research and Innovation Centre, Grange, Dunsany, Co. Meath, Ireland
2 CSIRO Agriculture & Food, Queensland Bioscience Precinct, 306 Carmody Rd., St. Lucia, Brisbane, QLD 4067 Australia
3 https://ror.org/03bea9k73 grid.6142.1 0000 0004 0488 0789 School of Biological and Chemical Sciences, Ryan Institute, University of Galway, Galway, Ireland
4 https://ror.org/03sx84n71 grid.6435.4 0000 0001 1512 9569 Livestock Systems Research Department, Teagasc, Animal & Grassland Research and Innovation Centre, Grange, Dunsany, Co. Meath, Ireland
20 9 2024
20 9 2024
2024
6 5216 7 2024
24 8 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by/4.0/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Background

Feed costs account for a high proportion of the variable cost of beef production, ultimately impacting overall profitability. Thus, improving feed efficiency of beef cattle, by way of determining the underlying genomic control and selecting for feed efficient cattle provides a method through which feed input costs may be reduced whilst also contributing to the environmental sustainability of beef production. The rumen microbiome dictates the feed degradation capacity and consequent nutrient supply in ruminants, thus potentially impacted by feed efficiency phenotype. Equally, liver tissue has been shown to be responsive to feed efficiency phenotype as well as dietary intake. However, although both the rumen microbiome and liver transcriptome have been shown to be impacted by host feed efficiency phenotype, knowledge of the interaction between the rumen microbiome and other peripheral tissues within the body, including the liver is lacking. Thus, the objective of this study was to compare two contrasting breed types (Charolais and Holstein-Friesian) divergent for residual feed intake (RFI) over contrasting dietary phases (zero-grazed grass and high-concentrate), based on gene co-expression network analysis of liver transcriptome data and microbe co-abundance network of rumen microbiome data. Traits including RFI, dry matter intake (DMI) and growth rate (ADG), as well as rumen concentrations of volatile fatty acids were also included within the network analysis.

Results

Overall, DMI had the greatest number of connections followed by RFI, with ADG displaying the fewest number of significant connections. Hepatic genes related to lipid metabolism were correlated to both RFI and DMI phenotypes, whilst genes related to immune response were correlated to DMI. Despite the known relationship between RFI and DMI, the same microbes were not directly connected to these phenotypes, the Succiniclasticum genus was however, negatively connected to both RFI and ADG. Additionally, a stepwise regression analysis revealed significant roles for both Succiniclasticum genus and Roseburia.faecis sp. in predicting RFI, DMI and ADG.

Conclusions

Results from this study highlight the interactive relationships between rumen microbiome and hepatic transcriptome data of cattle divergent for RFI, whilst also increasing our understanding of the underlying biology of both DMI and ADG in beef cattle.

Supplementary Information

The online version contains supplementary material available at 10.1186/s42523-024-00337-0.

Keywords

Feed efficiency
Beef cattle
Gene co-expression network analysis
http://dx.doi.org/10.13039/100010665 H2020 Marie Skłodowska-Curie Actions 754380 issue-copyright-statement© BioMed Central Ltd. 2024
==== Body
pmcBackground

Within beef production systems, feed costs alone account for up to 75% of the total variable cost of production [1], thus impacting overall profitability. Therefore, selecting and breeding cattle that are more feed efficient, for example through residual feed intake (RFI) phenotype, may increase production profitability by reducing feed input costs, whilst also contributing to the environmental sustainability of beef production. Moreover, feed intake (measured as average daily dry matter intake, DMI) and growth rate (measured as average daily gain, ADG) are important traits within the beef industry due to their direct and indirect effects on overall productivity and consequently profitability and sustainability. Thus, it is of interest for beef producers to rear faster growing animals with reduced dietary intake to optimise productivity of production systems [2]. Indeed, the aforementioned traits (DMI, ADG and RFI) have been shown to be moderately heritable [3–6], thus offering an opportunity to reduce feed costs through breeding cattle that are more feed efficient through genomic selection breeding programs. However, despite the clear benefit of breeding cattle that are feed efficient, evaluations of the molecular control regulating feed efficiency in beef cattle are not conclusive, such that key genes or genomic regions contributing to the trait are yet to be identified [1]. The contrasting results from molecular and genomic based studies across the literature are undoubtedly due to the multifaceted nature of the feed efficiency trait as well as the various confounding experimental parameters employed across studies such as breed, dietary management system and stage of development during which feed efficiency was evaluated. If reliable genomic selection processes are to be implemented for feed efficiency in beef cattle it is essential that genomic regions contributing to the trait are reliable across these various confounding factors.

The availability of nutrients for both growth and maintenance purposes in ruminants is dependent on the functionality of the rumen microbiome. Ruminants rely on the complement of bacteria, archaea and protozoa, amongst others, in the rumen microbiota for the degradation of feed and the production of microbial protein, vitamins and volatile fatty acids (VFAs), the latter of which provide up to 70% of the host’s energy requirements [7]. Research from our own group as well as others has identified links between the rumen microbiome and feed efficiency phenotype [8–13]. Moreover, Shabat et al. [14] reported that the rumen microbiome could predict variation in an animal’s feed efficiency phenotype and concluded that reduced rumen microbial community diversity may support a more feed efficient animal. However, similar to the identification of genomic regions associated with feed efficiency in beef cattle, results related to the contribution of the rumen microbiome to the feed efficiency phenotype are inconsistent, again, most likely due to the aforementioned confounding experimental parameters employed. Thus, it is important that such factors are considered and examined together [13]. Determination of the relationship between the rumen microbiome and host feed efficiency phenotype has the potential to not only facilitate the selection of cattle with enhanced nutrient utilisation, but also to enable the manipulation of the rumen microbiome to enhance its energy harvesting capacity [11, 14, 15].

The VFAs produced within the rumen following microbial degradation of feed are responsible for a large proportion of the host’s energy requirements. Furthermore, the liver, a highly metabolically active organ, has been shown to be responsive to dietary intake [16] as well as being affected by feed efficiency phenotype [2]. Thus, the objective of this study was to evaluate the interactions between the rumen microbiome, VFA concentrations and hepatic gene expression profiles with three production and efficiency traits of interest, namely DMI, RFI and ADG, through a network based systems biology analysis. In order to account for the confounding effects of both breed and dietary source, data utilised in this study (rumen microbiome, liver transcriptome, VFA and phenotype data) were derived from two contrasting breed types (Charolais and Holstein-Friesian) divergent for RFI across contrasting dietary phases (high concentrate and zero-grazed grass).

Results

Animal performance

Descriptive results pertaining to growth, dietary intake and RFI values are outlined in full in Higgins et al. [17]. Groups selected as divergent for RFI were significantly different from one another (P < 0.05). Across all dietary phases for each breed, High-RFI steers consumed more feed on average than their Low-RFI counterparts (P < 0.001), whilst there was no difference (P > 0.05) in ADG across groups for each breed and each dietary phase.

Microbial sequencing analysis

Complete results related to microbial sequencing are presented in full in McGovern et al. [13]. Briefly an average of 272,460 (± 69,596) reads were generated from the rumen fluid samples. Following merging of sequences and quality filtering, an average of 217,817 (± 55,519) reads were retained. The average number of counts per sample that were assigned to an open taxonomic unit (post filtering) was 175,304 ± 74,272. Microbial taxa identified within the rumen fluid samples and utilised within the current study are presented in Additional Table 1.

Liver gene expression

For the RFI, breed and diet contrasts, 12,161, 12,114 and 12,581 genes, respectively, were identified as expressed of which 608, 605 and 629 were classified as differentially expressed and used for subsequent co-expression network analysis. Of these differentially expressed genes, only 2 were common across the three contrasts, namely SPP1, which encodes a cytokine and ABHD2 which encodes an acylglycerol lipase protein. Genes identified as differentially expressed and included in the co-expression network analysis are presented in Additional Table 2.

Co-expression network analysis

Of the 2,017 nodes (1,842 differentially expressed genes, 159 microbial taxa, 13 VFA and 3 phenotypes) used for network analysis, 1,534  displayed significant correlations, resulting in a total of 35,549 significant (P < 0.05) connections between nodes (Fig. 1). Of the three phenotypes examined, DMI had the highest number of first neighbour connections, followed by RFI, with ADG displaying the fewest number of first neighbour connections (Table 1). Connections between hepatic genes accounted for the highest proportion of first neighbour connections for each of the three phenotypes, followed by microbial interactions, and VFA connections (Table 2). The co-expression network pertaining specifically to direct connections of the DMI, RFI and ADG phenotypes is presented in Fig. 2 and a comparison of nodes directly connected to DMI, RFI and ADG is provided in Fig. 3 and Additional Table 3. Direct connections between the phenotypes and microbial taxa revealed a lack of commonality between DMI and RFI, as well as DMI and ADG, with only the Succiniclasticum genus negatively connected to both RFI and ADG (Fig. 4). Details related to interactions of microbial taxa directly connected to DMI, RFI and ADG are presented in Table 3, along with enriched (P < 0.05) gene ontology and pathway analysis results. Full details related to interactions of microbial taxa connected to DMI, RFI and ADG are presented in full in Additional Table 4.

Fig. 1 Gene co-expression network constructed using PCIT algorithm on 2,017 input nodes related to hepatic genes differentially expressed, rumen microbial taxa, volatile fatty acid concentrations and DMI, RFI and ADG phenotypes

Table 1 Proportion of first neighbour nodes and direct connections pertaining to the DMI, RFI and ADG phenotypes

Phenotype	Number of first neighbour nodes	Number of direct connections	Percentage of first neighbour nodes from main network	Percentage of direct connections from main network	
DMI	180	811	11.7%	2.3%	
RFI	49	131	3.2%	0.37%	
ADG	37	66	2.4%	0.19%	

Table 2 Proportion of first neighbour phenotype nodes that are microbes, genes, VFAs or phenotypes

Node type	ADG	DMI	RFI	
Number of interactions	37	180	49	
Microbe	4 (10.8%)	10 (5.55%)	9 (18.4%)	
VFA	2 (5.4%)	5 (2.8%)	0 (0%)	
Phenotype	0 (0%)	1 (0.55%)	1 (2%)	
Gene	31 (83.8%)	164 (91.1%)	39 (79.6%)	

Fig. 2 Genes, microbial taxa and volatile fatty acids directly connected to the DMI, RFI and ADG phenotypes

Fig. 3 Venn diagram depicting commonality of first neighbour connections of the DMI, RFI and ADG phenotypes

Fig. 4 Microbial taxa directly interacting with (a) RFI, (b) DMI and (c) ADG. Orange lines depict positive connections, whist green lines indicate negative connections

Table 3 Connections of microbial taxa directly interacting with DMI, RFI and ADG phenotypes

Microbe	Connections	Microbe	VFA	Gene	Enriched term	
DMI						
 -C.aminophilum sp.	60	6	5	48	Steroid biosynthesis

Proteasomal processes

	
 -Butyrivibrio	69	24	4	40	Nucleic acid binding	
 -S.copri sp.	57	12	4	40	Cholesterol biosynthesis	
 -R.faecis sp.	186	83	9	93	Cholesterol homeostasis	
 -Lachnobacterium	53	6	2	44	Transcription	
 -Prevotella	73	6		66	Carbohydrate metabolism	
 -Roseburia	187	86	9	91	Steroid biosynthesis

Metabolic pathways

	
 -L.ruminis sp.	65	19	4	41	Retinoid binding	
 -L.salivarius sp.	160	74	4	81	Steroid biosynthesis

Metabolic pathways

	
 -Succinivibrio	42	20	1	20	Plasma membrane	
RFI						
 -Methanobrevibacter	426	78	2	346	Metabolic pathways

MAPK signaling

Lipid transporter activity

Catabolic processes

	
 -Pirellulales	331	72	5	253	Ras signaling pathway

Interferon signaling

Steroid catabolic process

Inflammatory response

	
 -Planctomycetes	323	72	5	245	Interferon signaling	
 -Planctomycetia	331	72	5	253	Ras signaling pathway

Interferon signaling

Steroid catabolic process

Inflammatory response

	
 -Pyramidobacter	89	21		67	B-cell receptor signaling

Chemokine signaling

T-cell receptor signaling

Immunodeficiency

	
 -Succiniclasticum	43	16		25	Plasma membrane	
 -Synergistales	185	33		151	Hippo signaling

Immunodeficiency

	
 -Synergistetes	185	33		151	Hippo signaling

Immunodeficiency

	
 -Synergistia	185	33		151	Hippo signaling

Immunodeficiency

	
ADG						
 -CW040	217	72	6	138	Steroid biosynthesis

Metabolic pathways

Citrate cycle

	
 -D.D168 sp.	26	5		20	Ras signaling pathway	
 -Desulfovibrio	27	4		22	Endoplasmic reticulum	
 -Succiniclasticum	43	16		25	Plasma membrane	
Taxonomic levels for all microbes listed are presented in Additional Table 1

Stepwise regression and correlation analyses

When each phenotype and their associated microbial connections were analysed separately through the regression analysis, 83% of the variation was accounted for through the Pyramidobacter and Succiniclasticum genera for the RFI phenotype. Similarly, the Succiniclasticum genus, together with the CW040 order explained 73% within the ADG trait, whilst the Butyrivibrio and Prevotella genera, together with S.copri sp. and R.faecis sp. explained 89.9% of the variation within the DMI trait. When all microbes directly connected to the phenotypes examined were analysed together, 90% of the variation for RFI, DMI and ADG was explained by the following microbial taxa: the Planctomycetes phylum; CW040 order; the Prevotella, Roseburia, Succiniclasticum, Methanobrevibacter, Pyramidobacter, Butyrivibrio genera and; S.copri sp., R.faecis sp., L.ruminis sp., C.aminophilum sp. and D.D168 sp. (Table 4). Of the significantly correlated microbes, only the Succiniclasticum genus and R.faecis sp. were commonly contributing to variation across the three traits examined.

Table 4 Stepwise regression analysis results of RFI, DMI and ADG phenotypes with significantly correlated microbial taxa. For each phenotype the microbes listed account for 90% of the variation

RFI	DMI	ADG	
Prevotella	Butyrivibrio	Succiniclasticum	
Roseburia	Prevotella	Methanobrevibacter	
Succiniclasticum	Succiniclasticum	C.aminophilum	
Methanobrevibacter	Planctomycetes	R.faecis sp.	
Pyramidobacter	Aminophilum	D.D168 sp.	
S.copri sp.	S. copri sp.	CW040	
R.faecis sp.	R.faecis sp.		
Ruminis			
Taxonomic levels for all microbes listed are presented in Additional Table 1

Results of the correlation analysis undertaken in SAS are presented in Fig. 5. Correlation results followed the same pattern as per the co-expression network results for RFI (Fig. 4), with negative correlations (P < 0.05) apparent between RFI and the Succiniclasticum and Pyramidobacter genera, the Synergistetes and Planctomycetes phyla, Planctomycetia and Synergistia classes and Synergistales and Pirellulales orders. In addition to the positive connections evident through the co-expression network analysis, through the CORR procedure in SAS, negative associations (P < 0.05) were also apparent between DMI and the Succiniclasticum and Methanobrevibacter genera and the CW040 order. Similar to both RFI and DMI, a negative association (P < 0.05) was also apparent between the Succiniclasticum genus and ADG, representing the only significant association between ADG and the microbial taxa from the correlation analysis.

Fig. 5 Results from correlation analysis between the RFI, DMI and ADG phenotypes and microbial taxa directly connected to the phenotypes examined. Green cells represent negative correlations, with positive correlations in red. Cells in bold font represent significant associations (P < 0.05)

Hierarchical clustering of the phenotypes and first neighbour microbial connections is presented in Fig. 6. From Fig. 6, it is evident that the RFI and DMI phenotypes were clustered with the Roseburia and Butyrivibrio genera, as well as with L.Salivarius sp. and R.Faecis sp. The ADG phenotype was separately clustered together with the Lachnobacterium and Prevotella genera, whilst also clustered with S.Copri sp., C.Aminophilium sp. and L.Ruminis sp. Hierarchical clustering also revealed a relationship between members of the same clade, for example the Pyramidobacter genus, Synergistetes phylum, Synergistia class and Synergistales order were all clustered together, with the Planctomycetia class, Planctomycetes phylum and Pirellulales order also clustered together separately.

Fig. 6 Heatmap of the hierarchical cluster analysis conducted in PermutMatrix between the ADG, RFI and DMI phenotypes with their first neighbour microbial connections derived from the PCIT results

Discussion

Within the published literature, biological processes related to immune function and lipid metabolism have been attributed to the underlying biology governing the RFI phenotype [17–21]. Indeed, results from the current study further corroborate that finding, where lipid metabolism related genes pertaining to the RFI differential expression contrast were connected to the RFI phenotype, whilst also connected to both DMI and ADG too. For example, RFI was positively connected to CYP26B1, whilst also negatively related to both PDK4 and OSBPL11. Indeed, all three of these differentially expressed genes were pertaining to the RFI contrast, whilst CYP26B1 was related to both RFI and diet contrasts, indicating a role for this gene in mediating the intersection of RFI with diet. CYP26B1 encodes a member of the cytochrome P450 superfamily which are responsible for catalysing reactions involved in the synthesis of cholesterol, steroids and other lipids. Additionally, CYP26A1, which functions similarly to CYP26B1 was also positively connected to the DMI phenotype. OSBPL11 encodes an intracellular lipid receptor, which plays a role in regulating ADIPOQ and FABP4 levels in differentiating adipocytes and is also involved in regulating adipocyte triglyceride storage. Whilst PDK4 encodes a mitochondrial protein which functions in the regulation of both glucose and fatty acid metabolism. The relevance of the PDK4 gene towards the RFI phenotype is further apparent through previous reports of this gene in relation to RFI, specifically PDK4 was down-regulated in skeletal muscle tissue of efficient Holstein-Friesian bulls following a high-concentrate finishing diet [21], as well as in the liver tissue of efficient Charolais steers [17]. Of the lipid related genes connected to the DMI phenotype, all were pertaining specifically to the RFI differential expression contrast, whilst also being directly connected to the DMI phenotype, suggesting a role for these genes towards both RFI and DMI phenotypes. Moreover, a number of these genes have previously been implicated with variation in RFI, including ASCL1 [20, 22]; APOA1 [19, 23], ELOVL2 [19] and FADS1 [20]. The FADS1 gene, which encodes a desaturase enzyme and regulates the unsaturation of fatty acids was also connected to ADG. Interestingly of the lipid metabolism genes connected to ADG (HMGCS1, FDFT1, SQLE and FADS1), all were pertaining to the diet differential expression contrast, with the exception of FADS1, highlighting the influence of contrasting diets with growth rates in beef cattle.

A role for altered immune response towards the RFI phenotype across various tissues throughout the body including the liver, has already been established [17–21]. Interestingly within the current study we did not observe such a relationship between the RFI phenotype and immune response genes, with the exception of, C1QTNF7 which was negatively connected to RFI. However, although not strongly apparent within the RFI phenotype, a more pronounced relationship of immune response genes was evident with the DMI phenotype. This was apparent through connections between DMI and the following genes: A2M, C1QBP, C7, CCL21, CD19, CD22, CRP, FCER2, FCRL5, FNDC4, HMCN1, IFNLR1 and IRF2BP2. Indeed, A2M, CCL21 and CRP were previously reported as differentially expressed between cattle divergent for RFI [22, 23], suggesting the importance of these genes to both the RFI and DMI traits. The relationship between the DMI phenotype and immune response is further established through the connection between the Prevotella genus and S.copri sp. (Segatella copri sp. formerly Prevotella copri sp.) of bacteria within the rumen microbiome, both of which have been implicated in immune function. For example, high abundance of P.copri sp. was correlated with increased concentrations of serum metabolites related to chronic inflammation in the gut of pigs in the data of Chen et al. [24]. Whilst the Prevotella genus was reported to be associated with gut mucosal inflammation in mice [25]. Moreover, a positive relationship was apparent in the current study between the Prevotella genus and DMI, whilst Jiang et al. [26] and Carberry et al. [8] both conversely reported negative correlations between relative abundance of ruminal Prevotella genus. Moreover, both Prevotella genus and S.copri sp. have been related to RFI within various species too. For example, Jiang et al. [27] reported that P.copri sp. was one of the most abundant microbes in low-feed efficiency pigs and also noted the importance of that particular species to feed efficiency during daily-phase feeding strategies in pigs [28]. Furthermore, in cattle, Brooke et al. [29] suggested that P.copri sp. may be a potential microbial marker for the identification of cattle with improved feed efficiency in their life-span and in the production cycle. Additionally, Carberry et al. [8] reported an effect of RFI phenotype on ruminal abundance of Prevotella genus, irrespective of the varied diets offered in that study, but did note greater abundance of Prevotella genus when cattle were offered a low forage diet compared to a high forage diet. Conversely though, Lopes et al. [30] reported that the correlation between the Prevotella genus and feed efficiency was dependent on the diet offered as well as the specific Prevotella spp. Indeed, different Prevotella spp. have been associated with both higher and lower feed efficiency in cattle and sheep [29, 31–35]. Furthermore, Zhou et al. [36] reported that the Prevotella genus was the most abundant genus in both rumen content-associated and epithelial tissue-attached bacterial communities suggesting a role for this genus in serving as a marker for host RFI classification. Additionally, Jewell et al. [9] reported that the abundance of specific Prevotella spp. is host specific, with Yang et al. [37] suggested that Prevotella may be a key microbe increasing host feed intake, suggesting that Prevotella spp. could promote the host’s appetite and decrease feed efficiency. Overall, results indicate that the role of varied immune response in relation to RFI may also be attributed to variation in DMI. 

Despite the known relationship between RFI and DMI, namely RFI being computed after accounting for DMI, results from the current study highlight differential connections between these two phenotypes and ruminal microbiome taxa, whereby microbial taxa directly connected to RFI were different to those directly connected to DMI. For the DMI phenotype all direct connections between DMI and microbial taxa were positive, with interactions amongst these specific taxa also positive. Conversely though, for the RFI trait, all connections between RFI and microbial taxa were negative, whilst connections between the RFI taxa were positive. Overall suggesting differential microbial responses, despite the relationship between the two traits.

The Butyrivibrio genus, which was positively related to DMI in the current study, is involved in a number of ruminal functions in addition to butyrate production including fibre degradation, protein breakdown, biohydrogenation of lipids and the production of microbial inhibitors. Of particular importance to ruminant digestion, and therefore productivity is the contribution of this bacteria to the degradation of plant structural carbohydrates, principally hemicellulose. Together with the Prevotella genus, both Butyrivibrio and Prevotella genera are among the most abundant bacteria found within the rumen and have important functions in the metabolism of proteins and peptides [34]. Both of these genera breakdown protein and carbohydrates in feed, undertake de novo peptide synthesis and use products of cellulose degradation from other cellulotyic bacteria as an energy source. In addition to being related to DMI, hierarchical clustering analysis revealed a relationship between the Butyrivibrio genus with RFI. Indeed, within the context of feed efficiency, Myer et al. [10] and McGovern et al. [12] reported greater abundance of the Butyrivibrio genus in efficient steers and a negative correlation with RFI, respectively. Conversely, Jewell et al. [9] reported greater abundance in the rumen of High-RFI dairy cows, indicating a positive relationship between RFI and rumen Butyrivibrio genus. The relative importance of the Butyrivibrio genus towards feed efficiency is also apparent in monogastrics, specifically the Butyrivibrio genus was strongly correlated with feed efficiency in pigs during the weaning phase, which the authors attributed to potentially be due to an enhanced ability to ferment complex carbohydrates [38]. Similarly, Kubasova et al. [39] also reported greater abundance of Butyrivibrio genus in fecal samples of pigs. Although the Butyrivibrio genus is a primarily butyrate producing bacteria, results from the current interaction study reported direct connections with only the lactic acid VFAs, which in turn were connected to hepatic genes involved in functions related to gene expression, lipid metabolism, growth, as well as protease activity, which given the function of the Butyrivibrio genus towards protein breakdown is of interest. Specifically, genes involved in the regulation of peptidase activity (PI16) as well as those with an associated protease function (CAPN6, SERPINA1, SERPINA3-1) were connected to lactic acid VFAs within the rumen. Moreover, of these, PI16, SERPINA1 and SERPINA3 were previously reported as differentially expressed in the RFI based datasets of Keogh et al. [21], Weber et al. [22], and Alexandre et al. [23]. The lactic acid VFAs were also connected to the leptin receptor gene (LEPR) highlighting the role of leptin towards mediating satiety status within the body. Moreover, the LEPR gene was differentially expressed through the breed contrast, highlighting differential hepatic expression of this gene between Charolais and Holstein-Friesian steers. Additionally, genes involved in TGF-beta growth signaling (BMP2, CRM1A, FNDC4 and HMCN1) were also connected to ruminal lactic acid concentrations, which was of interest as Alexandre et al. [23] previously reported TGFB1 as a key regulator for feed efficiency in skeletal muscle of Nellore cattle.

Similar to the Butyrivibrio genus, the R.faecis sp. is also a primarily butyrate producing bacteria. Indeed, although correlated with DMI in the current study, as well as clustered with both DMI and RFI through the hierarchical clustering analysis, R.faecis sp. was previously associated with growth rate in pigs, with a greater abundance of this microbe in pigs supplemented with a carbohydrate complex diet, which the authors of that study suggested could have accounted for the improved feed efficiency observed in those animals [40]. However, through the regression analysis in the current study, R.faecis sp. was identified, together with the Succiniclasticum genus, as playing a role in determining variation in the RFI, DMI and ADG traits examined in this study. Volatile fatty acids connected to R.faecis sp. (lactic acid, ammonia, valeric acid, acetic:propionic, isobutyric acid and total-SCFA), were also connected to hepatic genes involved in processes related to gene expression, growth, immune and lipid metabolism. In the current study, the Roseburia genus was positively associated with DMI, this microbe utilises carbohydrates for growth and its abundance is known to increase with greater proportions of concentrates within the diet [32, 41]. Similarly, Li et al. [42] and Ellison et al. [33] reported greater abundance of the Roseburia genus in cattle fed a high-energy diet and lambs fed a concentrate diet, respectively, indicating a role for this bacterium depending on diet composition. Additionally through its connections with isobutyric acid, the R.faecis sp. was connected to CMKLR2 and ND6 which encode proteins involved in adipokinetic hormone activity and glucose homeostasis as well as a mitochondrial gene, respectively, potentially suggesting a role for R.faecis sp. in mediating hepatic glucose homeostasis and mitochondrial function as a consequence of dietary intake. Moreover, CMKLR2 was specifically pertaining to the diet contrast, whilst ND6 was differentially expressed in both diet and RFI contrasts, overall highlighting the effect of diet on these genes related to R.faecis sp. abundance. Whilst a role for these genes related to RFI is already established through their differential expression within the published literature [23, 43].

The Synergistetes phylum was negatively connected with RFI in the current study. Similarly, McLoughlin et al. [44] also reported a negative correlation between this phylum and feed efficiency in the solid rumen fraction of sheep. Additionally, McCormack et al. [45] also reported a role for this microbe towards variation in RFI in pigs. In addition to the negative relationship between the Synergistetes phylum and RFI; class, order and genus members of this clade were also negatively connected to RFI. Moreover, the members of this specific phylum were all positively connected to each other. Specifically other members of this clade negatively connected to RFI included the Synergistia class, Synergistales order and Pyramidobacter genus. Indeed, the relationship between the members of this specific phylum was apparent through the PCIT; correlation and; hierarchical clustering analyses undertaken. Of these microbial taxa, the Pyramidobacter genus has previously been implicated towards variation in feed efficiency. For example, McLoughlin et al. [44] identified positive associations between the relative abundance of the Pyramidobacter genus and feed conversion ratio in sheep, whilst also reporting a negative association with ADG. Similarly in Simmental bulls, McGovern et al. [12] identified an association between RFI and abundance of the Pyramidobacter genus. Whilst in pigs McCormack et al. [45] and Kubasova et al. [39] identified an association between the Pyramidobacter genus and RFI in the cecal digesta of pigs and reported greater abundance of the Pyramidobacter genus in the fecal microbiota of Low-RFI pigs compared to High-RFI pigs, respectively. Similar to the Synergistetes phylum, the Planctomycetes phylum as well as class sub-member (Planctomycetia), and order member (Pirellulales) were also negatively correlated with RFI, whilst positively correlated with one another. Indeed, the Planctomycetes phylum was shown to be altered in both pigs and hens divergent in feed efficiency potential [45, 46]. In cattle, Freetly et al. [47] reported an effect of ADG on abundance of the Planctomyetes phylum, Planctomycetia class and Pirellulales order in the rumen of beef cattle, with results from this current study highlighting a role for these microbes towards RFI in beef cattle.

Of all the microbes directly connected to the phenotypes examined in this study only one was commonly significantly correlated across more than one phenotype. Specifically, the Succiniclasticum genus was negatively correlated with both RFI and ADG. Moreover, through the regression analysis, the Succiniclasticum genus was observed to contribute to variation in all three traits examined in this study. A role for the Succiniclasticum genus towards variation in feed efficiency has been established previously within the published literature, for example, in a study by Myer et al. [10] more efficient steers were observed to have greater abundance of succinate producing bacteria including the Succiniclasticum genus. Auffret et al. [48] also reported significantly greater abundance of the Succiniclasticum genus in high feed efficiency beef cattle. Conversely though, Manzanares-Miranda et al. [49] reported lower abundance of the Succiniclasticum genus in Low-RFI bulls. Interestingly, an effect of diet on the abundance of Succiniclasticum genus is evident within the literature. For example, McCann [50] reported greater proportions of the Succiniclasticum genus in steers consuming low quality forage, however the same authors observed abundance to be mostly undetected in forage diets but more abundant in a high grain diet [32]. Additionally, McCabe et al. [51] reported alterations to the abundance of the Succiniclasticum genus dependent on the amount of feed consumed, whilst Luo et al. [52] showed that a high concentrate diet increased the abundance of the Succiniclasticum genus. Moreover, there is also evidence for an effect of breed on ruminal abundance of Succiniclasticum genus. For example, the Succiniclasticum genus was higher in the High-RFI beef cattle in Li et al. [42], however this difference was only apparent in Charolais and not in the other breeds examined in that study. Similarly, in sheep, McLoughlin et al. [53] identified higher abundance of the Succiniclasticum genus in the Connemara breed of sheep compared to other breeds. Together these results indicate towards a role for the Succiniclasticum genus towards the phenotypes examined in this study, however the contribution may be dependent on both individual genotype as well as the dietary management system in place.

Conclusions

The rumen microbiome influences the availability of nutrients for subsequent growth purposes through the degradation of ingested feed. Equally the liver is a highly metabolically active organ and both the rumen microbiome and liver have been shown to be affected by RFI phenotype as well as dietary intake. Results from this study highlight the interaction amongst the rumen microbiome, VFA concentrations and hepatic gene expression profiles with three production and efficiency traits of interest, namely DMI, ADG and RFI. Results from this interactive study show a clear relationship between hepatic genes related to lipid metabolism towards RFI, whilst genes with an associated immune function were reported as primarily related to DMI. This study also indicated towards differential microbiome interactions between RFI and DMI despite the known correlation between these two traits, whilst the Succiniclasticum genus was identified as the only microbe connected to more than one trait, namely ADG and RFI. However, results highlight a potential role for both the Succiniclasticum genus and R.faecis sp. towards RFI, DMI and ADG phenotypes in beef cattle. Taken together, this study provides insights into the interaction amongst rumen microbiome and hepatic gene expression, which may be contributing to the underlying biology of DMI, ADG and RFI in beef cattle.

Methods

Animal management and phenotype collection

The animal model utilised in this study was conducted as part of a larger research programme designed to investigate the within-animal repeatability of feed intake, growth and feed efficiency in two contrasting breeds (Charolais and Holstein-Friesian) of beef steers, which were offered contrasting diets over separate dietary intake test periods [54, 55]. Details related to the original animal model are described previously in Higgins et al. [17] and are only briefly outlined here. Charolais (n = 90) and Holstein-Friesian (n = 77) steers were offered contrasting diets over different stages of development as follows: (i) a high-concentrate diet during the growing phase; (ii) zero-grazed grass diet during the growing phase and; (iii) high-concentrate diet during the finishing phase. Charolais steers were on average 373 (± 18) days of age and weighed 485 (± 38) kg, whilst Holstein-Friesian steers were 399 (± 7.6) days of age and weighed 401 (± 43.3) kg at the start of the trial. Upon completion of a dietary adaptation period, lasting 14 days, individual animal intakes were recorded (using an electronic Calan gate system; American Calan Inc., Northwood, NH, USA) over the three feeding phases, which each lasted for 70 days. Steers were weighed at the beginning and end of each dietary phase as well as on a fortnightly basis throughout. All steers were offered the same concentrate diet ad libitum during each of the two high-concentrate phases, with a restricted allowance of grass silage also provided. For the interim zero-grazed grass phase, steers were individually offered fresh herbage, harvested twice daily from Lolium perenne dominant swards, ad libitum. All steers had unrestricted access to fresh, clean drinking water. Upon completion of each dietary phase, individual RFI values were determined within breed for all steers as previously described in Higgins et al. [17], and animals were ranked as either High-RFI or Low-RFI, selecting the highest High-RFI and lowest Low-RFI for subsequent analyses.

Rumen digesta sampling and sequencing

Full details related to rumen digesta sampling and subsequent sequencing analysis are described in full in McGovern et al. [13]. At the midpoint of each dietary phase, following a dietary adaptation period, a single rumen fluid sample was collected from all steers via stomach intubation (Flora Rumen Scoop, Profs-Products, Guelph, Canada). All samples were harvested approximately 2–4 h post-feeding, and were immediately snap frozen in liquid nitrogen and subsequently stored at -80 °C, pending further analysis. Ten rumen digesta samples per High- and Low-RFI groups from each breed and diet were used for subsequent microbial DNA isolation, with the exception of the Low-RFI Charolais and High-RFI Holstein-Friesian steers during the zero-grazed grass diet and the Low-RFI Holstein Friesian steers during the second high-concentrate diet, where only 9 samples were available for each. Frozen rumen liquid samples (20 g) was homogenised to a fine powder under liquid nitrogen using a pestle and mortar and stored at -80 °C. Approximately 250 mg of the homogenised frozen powder was then used for DNA isolation, which was undertaken using the repeated bead beating and column purification method [56]. The quality of the resultant DNA samples was assessed on an agarose gel, with DNA yield and purity also assessed on a Nanodrop 1000 spectrophotometer. Amplicon libraries were prepared through PCR amplification, targeting the V4 region of the 16 S rRNA gene in both bacteria and archaea. Full details of library preparation are previously outlined in McGovern et al. [13]. Amplicon generation was validated through visualisation on an agarose gel. Amplicons were pooled in equal concentrations and gel purified to remove unwanted products using the Qiagen Gel Extraction Kit (Qiagen, Manchester, UK). The pooled purified libraries were measured for purity and quantity on the Nanodrop 1000 spectrophotometer and further quantified using the KAPA SYBR FAST universal kit with Illumina Primer Premix (Roche Diagnostics, West Sussex, UK). The library pool was then diluted and denatured according to the Illumina MiSeq library preparation guide. The sequencing was conducted using 500 cycle MiSeq reagent kits (Illumins, San Diego, CA, USA). The sequencing reads generated were imported into Qiime2 [57], where the DADA2 pipeline [58] was used for the detection of operational taxonomic units [59]. Taxonomy was assigned using a naïve Bayes classifier trained on the RefSeq database [60]. Sequence read files associated with this analysis are available through the NCBI Sequence Read Archive (Accession no. PRJNA483745).

Rumen fermentation profiling

The concentration of ruminal VFA composition was measured using a gas chromatograph (model 3800 Varian gas chromatograph) as per McGovern et al. [13]. The concentration of the following acids was determined: acetic, propionic, isobutyric, butyric, isovaleric, valeric, total short chain fatty acids, acetic:propionic acid ratio; D-lactic acid; L-lactic acid; DL-lactic acid; g-lactic acid and ammonia (NH3).

Liver tissue sample collection and RNA-sequencing

Liver tissue sample collection and RNA-sequencing methodology is described in full in Higgins et al. [17]. Briefly, at the end of each dietary phase and within breed, steers were ranked as either High-RFI (feed-inefficient; n = 12) or Low-RFI (feed-efficient; n = 12) and were subsequently used for the collection of liver tissue biopsies. All steers selected for biopsy collection were administered a local anaesthetic (5 ml; Adrenacaine, Norbrook Laboratories, Ireland Ltd.) to the biopsy site location. Following anaesthetisation, liver tissue was harvested through percutaneous punch between the 11th and 12th ribs as previously described by McCarthy et al. [61]. Care was taken to ensure that all samples were consistently harvested from the same location from each animal. All instruments used for biopsy collection were sterilised, washed with 70% ethanol and treated with RNaseZap (Ambion, Applera Ireland, Dublin, Ireland), prior to use. Following collection, all tissue samples were washed with sterile DPBS and immediately snap frozen in liquid nitrogen before subsequent storage at -80 °C pending further processing.

Full details related to RNA isolation and subsequent RNA sequencing and bioinformatic analysis are provided in Higgins et al. [17] and are only summarised here. Briefly, 50 mg of liver tissue from each biopsy sample was used for the isolation of total RNA. RNA was isolated from tissue samples in 3 ml of QIAzol reagent using a rotor-stator tissue lyser (Qiagen, UK). RNA was subsequently precipitated and purified using the RNeasy plus Universal kit (Qiagen, UK) according to the manufacturer’s instructions. Quality and quantity of RNA isolated were determined using the RNA 6000 Nano Lab Chip kit (Aglient Technologies Ireland Ltd., Dublin, Ireland) on an Aglient Bioanalyser 2100 and using a Nanodrop spectrophotometer (Nanodrop Technologies, Wilmington, DE, USA), respectively. All RNA samples displayed RNA integrity numbers (RIN) greater than 8 and thus were deemed to be of suitable quality for subsequent RNA-sequencing. Individual cDNA libraries were prepared from each separate liver RNA sample for cattle divergent for RFI across each breed and dietary phase, using the Illumina TruSeq stranded mRNA sample prep kit (Illumina, San Diego, CA, USA) according to the manufacturer’s instructions. Resultant cDNA libraries were validated using the DNA 1000 Nano Lab Chip kit on the Aglient Bioanalyser 2100. Sequencing was subsequently undertaken on an Illumina HiSeq 2500 sequencer. All sequencing data used in this study are publicly available in NCBI’s Gene Expression Omnibus and can be accessed through GEO ID GSE111464.

Quality control of sequencing reads was undertaken using FastQC (v 0.11.5; [62]), followed by removal of sequencing adapters and any low quality reads using Cutadapt (v 1.13; [63]). Trimmed sequencing reads were mapped to the bovine reference genome (ARS-UCD1.2; [64]) and also quantified using STAR (v.2.5.1; [65]). Differential expression was undertaken using the edgeR package within the R environment [66]. Within edgeR, gene expression reads were estimated as Counts Per Million (CPM) and genes which presented with at least 1 CPM in at least half of the samples were retained for differential expression analysis. Differentially expressed genes were identified for each of the main contrasts of RFI phenotype (Low-RFI versus High-RFI), breed (Charolais versus Holstein-Friesian) and dietary source (high concentrate versus zero-grazed grass). The model for differential expression due to a given contrast (e.g., RFI) contained the other two contrasts as main effects (e.g., breed and diet).

Co-expression network analysis

The following datasets were utilised for the gene co-expression network analysis: liver transcriptomics; 16 S rumen microbial abundance taxonomy; VFA concentrations; and phenotype data (DMI, RFI and ADG; Fig. 7). From the differential expression analysis conducted in edgeR, the top 5% differentially expressed genes, for each of the three contrasts (RFI, breed, and diet) based on corrected p-value were selected for subsequent inclusion in the co-expression network analysis. Thus, nodes selected for subsequent co-expression analysis included: (i) genes differentially expressed based on the RFI contrast; (ii) genes differentially expressed based on the breed contrast; (iii) genes differentially expressed based on the diet contrast; (iv) ruminal microbial abundance; (v) VFA concentrations and; (vi) three phenotypes of interest: RFI, DMI and ADG. Significant connections (edges) between nodes were identified using the Partial Correlation and Information Theory (PCIT) algorithm [67]. The PCIT algorithm determines the significance of the correlation between a pair of nodes after accounting for all other nodes within the network [67]. The resultant network of co-expressed genes was imported into Cytoscape software [68] for visualisation. In order to assign biological annotation of the generated network, co-expression network analysis results were further evaluated for functional enrichment using Gprofiler and David gene ontology.

Fig. 7 Overview of data used for co-expression network analysis. ([VFA]: volatile fatty acid concentrations; H-RFI: High-RFI; L-RFI: Low-RFI; ZG: zero-grazed grass diet; HC: high concentrate diet; CH: Charolais; HF: Holstein-Friesian; DMI: dry matter intake; RFI: residual feed intake; ADG: average daily gain)

Stepwise regression and correlation analyses

To determine the independent rumen microbial predictors of RFI, DMI and ADG, a stepwise regression analysis was conducted. This analysis was undertaken in SAS (version 9.4) using the REG procedure, incorporating the three phenotypes examined in this study as well as any microbe directly connected to the phenotypes, based on the results from the co-expression network analysis within the model and utilising R2 as the selection criteria. A correlation analysis was undertaken on the aforementioned phenotype and microbe variables using the CORR procedure of SAS. Additionally, a hierarchical clustering analysis was performed between the phenotypes and first neighbour microbial connections using PermutMatrix (version 1.9.4; http://www.atgc-montpellier.fr/permutmatrix/, [69]).

Electronic supplementary material

Below is the link to the electronic supplementary material.

Additional Table 1: Microbial taxa identified within the rumen fluid and included within the co-expression network analysis. Additional Table 2: Top 5% differentially expressed genes for each RFI, breed and diet contrasts in liver tissue. Additional Table 3: First neighbours for each of the three phenotypes examined (DMI, RFI and ADG). Additional Table 4: First neighbours for microbial taxa directly connected to each of the three phenotypes examined (DMI, RFI and ADG)

Acknowledgements

The authors would like to acknowledge receipt of funding for the animal model used in this study from the Irish Department of Agriculture, Food and the Marine (DAFM) via the IdentiFEED project (13/S/519). Kate Keogh received funding for the work conducted in this study from the Research Leaders 2025 programme co-funded by Teagasc and the European Union’s Horizon 2020 research and innovation programme under Marie Skłodowska-Curie grant agreement number 754380.

Author contributions

Conceived and designed the experiments: KK and AR. Management and design of animal model: DAK, MM. Data analysis: KK, AR, PA, EM, SM. Prepared manuscript: KK. Edited manuscript: KK, DAK, MM, PA, AR, EM, SM. All authors have read and approved the final manuscript.

Funding

The animal model utilised in this study was funding by the Irish Department of Agriculture, Food and the Marine (DAFM) via the IdentiFEED project (13/S/519). Kate Keogh received funding for the work conducted in this study from the Research Leaders 2025 programme co-funded by Teagasc and the European Union’s Horizon 2020 research and innovation programme under Marie Skłodowska-Curie grant agreement number 754380.

Data availability

The transcriptomic datasets utilised for this study can be found in the NCBI’s Gene Expression Omnibus (GEO) database [https://www.ncbi.nlm.nih.gov/geo/] (GEO accession ID: GSE111464). Sequence files related to the 16s sequencing are available in NCBI Sequence Read Archive (Accession no. PRJNA483745).

Declarations

Ethics approval and consent to participate

This study was conducted at the Teagasc Animal & Grassland Research and Innovation Centre Grange, Co. Meath in Ireland. All procedures involving animals were reviewed and approved by the Teagasc Animal Ethics Committee and all procedures involving animals were conducted under an experimental license (AE19132/P029) issued by the Health Products Regulatory Authority in Ireland in accordance with the cruelty to Animals Act 1876 and the European Communities (Amendment of Cruelty to Animals Act 1876) Regulations 2002 and 2005. All experiments were performed in accordance with relevant regulations and the ARRIVE (Animal Research: Reporting on In Vivo Experiments) guidelines.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Abbreviations

ADG Average daily gain

CPM Counts per millions

DMI Dry matter intake

DPBS Dulbecco’s phosphate buffered saline

GEO Gene expression omnibus

H1 High concentrate diet during the growing phase

H2 High concentrate diet during the finishing phase

HC High concentrate

RFI Residual feed intake

RIN RNA integrity number

VFA Volatile fatty acid

ZG Zero-grazed grass diet during the growing phase

Publisher’s note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

1. Kenny DA Fitzsimons C Waters SM McGee M Improving feed efficiency of beef cattle; current state of the art and future challenges Animal 2018 12 1815 26 10.1017/S1751731118000976 29779496
Kenny DA, Fitzsimons C, Waters SM, McGee M. Improving feed efficiency of beef cattle; current state of the art and future challenges. Animal. 2018;12:1815–26.29779496
2. Hill RA Ahola JK Hill RA Feed efficiency interactions with other traits: growth and product quality Feed efficiency in the beef industry 2012 Ames, IA Wiley-Blackwell 148
Hill RA, Ahola JK. Feed efficiency interactions with other traits: growth and product quality. In: Hill RA, editor. Feed efficiency in the beef industry. Ames, IA: Wiley-Blackwell; 2012. p. 148.
3. Schenkel F Miller S Wilton J Genetic parameters and breed differences for feed efficiency, growth, and body composition traits of young beef bulls Can J Anim Sci 2004 84 177 86 10.4141/A03-085
Schenkel F, Miller S, Wilton J. Genetic parameters and breed differences for feed efficiency, growth, and body composition traits of young beef bulls. Can J Anim Sci. 2004;84:177–86.
4. Nkrumah J Basarab J Wang Z Li C Price M Okine E Genetic and phenotypic relationships of feed intake and measures of efficiency with growth and carcass merit of beef cattle J Anim Sci 2007 85 2711 20 10.2527/jas.2006-767 17526662
Nkrumah J, Basarab J, Wang Z, Li C, Price M, Okine E, et al. Genetic and phenotypic relationships of feed intake and measures of efficiency with growth and carcass merit of beef cattle. J Anim Sci. 2007;85:2711–20.17526662
5. Mao F Chen L Vinsky M Okine E Wang Z Basarab J Phenotypic and genetic relationships of feed efficiency with growth performance, ultrasound, and carcass merit traits in Angus and Charolais steers J Anim Sci 2013 91 2067 76 10.2527/jas.2012-5470 23463551
Mao F, Chen L, Vinsky M, Okine E, Wang Z, Basarab J, et al. Phenotypic and genetic relationships of feed efficiency with growth performance, ultrasound, and carcass merit traits in Angus and Charolais steers. J Anim Sci. 2013;91:2067–76.23463551
6. Berry DP, Crowley JJ. Cell biology symposium: genetics of feed efficiency in dairy and beef cattle. J Anim Sci. 2013;91:1594–1613.
7. Flint HJ Bayer EA Rincon MT Lamed R White BA Polysaccharide utilization by gut bacteria: potential for new insights from genomic analysis Nat Rev Microbiol 2008 6 121 31 10.1038/nrmicro1817 18180751
Flint HJ, Bayer EA, Rincon MT, Lamed R, White BA. Polysaccharide utilization by gut bacteria: potential for new insights from genomic analysis. Nat Rev Microbiol. 2008;6:121–31.18180751
8. Carberry CA Kenny DA Han S McCabe MS Waters SM Effect of phenotypic residual feed intake and dietary forage content on the rumen microbial community of beef cattle Appl Environ Microbiol 2012 78 4949 58 10.1128/AEM.07759-11 22562991
Carberry CA, Kenny DA, Han S, McCabe MS, Waters SM. Effect of phenotypic residual feed intake and dietary forage content on the rumen microbial community of beef cattle. Appl Environ Microbiol. 2012;78:4949–58.22562991
9. Jewell KA McCormick CA Odt CL Weimer PJ Suen G Ruminal Bacterial Community composition in dairy cows is dynamic over the course of two lactations and correlates with feed efficiency Appl Environ Microbiol 2015 81 4697 710 10.1128/AEM.00720-15 25934629
Jewell KA, McCormick CA, Odt CL, Weimer PJ, Suen G. Ruminal Bacterial Community composition in dairy cows is dynamic over the course of two lactations and correlates with feed efficiency. Appl Environ Microbiol. 2015;81:4697–710.25934629
10. Myer PR Smith TP Wells JE Kuehn LA Freetly HC Rumen microbiome from steers differing in feed efficiency PLoS ONE 2015 10 e0129174 10.1371/journal.pone.0129174 26030887
Myer PR, Smith TP, Wells JE, Kuehn LA, Freetly HC. Rumen microbiome from steers differing in feed efficiency. PLoS ONE. 2015;10:e0129174.26030887
11. Roehe R Dewhurst RJ Duthie CA Rooke JA McKain N Ross DW 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 2016 12 e1005846 10.1371/journal.pgen.1005846 26891056
Roehe R, Dewhurst RJ, Duthie CA, Rooke JA, McKain N, Ross DW, et al. 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. 2016;12:e1005846.26891056
12. McGovern E Kenny DA McCabe MS Fitzsimons C McGee M Kelly AK 16S rRNA sequencing reveals relationship between Potent Cellulolytic Genera and feed efficiency in the Rumen of bulls Front Microbiol 2018 9 1842 10.3389/fmicb.2018.01842 30147683
McGovern E, Kenny DA, McCabe MS, Fitzsimons C, McGee M, Kelly AK, et al. 16S rRNA sequencing reveals relationship between Potent Cellulolytic Genera and feed efficiency in the Rumen of bulls. Front Microbiol. 2018;9:1842.30147683
13. McGovern 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 15317 10.1038/s41598-020-71458-0 32948787
McGovern 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:15317.32948787
14. Shabat SK Sasson G Doron-Faigenboim A Durman T Yaacoby S Berg Miller ME Specific microbiome-dependent mechanisms underlie the energy harvest efficiency of ruminants ISME J 2016 10 2958 72 10.1038/ismej.2016.62 27152936
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:2958–72.27152936
15. Li F Metatranscriptomic profiling reveals linkages between the active rumen microbiome and feed efficiency in beef cattle Appl Environ Microbiol 2017 83 e00061 17 10.1128/AEM.00061-17 28235871
Li F. Metatranscriptomic profiling reveals linkages between the active rumen microbiome and feed efficiency in beef cattle. Appl Environ Microbiol. 2017;83:e00061–17.28235871
16. Johnson DE Johnson KA Baldwin RL Changes in liver and gastrointestinal tract energy demands in response to physiological workload in ruminants J Nutr 1990 120 649 55 10.1093/jn/120.6.649 2191096
Johnson DE, Johnson KA, Baldwin RL. Changes in liver and gastrointestinal tract energy demands in response to physiological workload in ruminants. J Nutr. 1990;120:649–55.2191096
17. Higgins MG Kenny DA Fitzsimons C Blackshields G Coyle S McKenna C The effect of breed and diet type on the global transcriptome of hepatic tissue in beef cattle divergent for feed efficiency BMC Genomics 2019 20 525 10.1186/s12864-019-5906-8 31242854
Higgins MG, Kenny DA, Fitzsimons C, Blackshields G, Coyle S, McKenna C, et al. The effect of breed and diet type on the global transcriptome of hepatic tissue in beef cattle divergent for feed efficiency. BMC Genomics. 2019;20:525.31242854
18. Fitzsimons C, McGee M, Keogh K, Waters SM, Kenny DA. Molecular physiology of feed efficiency in beef cattle in Biology of Domestic Animals (ed. Hill, R.) CRC Press. 2017;180–231.
19. Mukiibi R Vinsky M Keogh K Fitzsimmons C Stothard P Waters SM Transcriptome analyses reveal reduced hepatic lipid synthesis and accumulation in more feed efficient beef cattle Sci Rep 2018 8 7303 10.1038/s41598-018-25605-3 29740082
Mukiibi R, Vinsky M, Keogh K, Fitzsimmons C, Stothard P, Waters SM, et al. Transcriptome analyses reveal reduced hepatic lipid synthesis and accumulation in more feed efficient beef cattle. Sci Rep. 2018;8:7303.29740082
20. Taiwo G Idowu MD Wilson M Pech-Cervantes A Estrada-Reyes ZM Ogunade IM Residual feed intake in beef cattle is associated with differences in hepatic mRNA expression of fatty acid, amino acid and mitochondrial energy metabolism genes Front Anim Sci 2022 3 838591 10.3389/fanim.2022.828591
Taiwo G, Idowu MD, Wilson M, Pech-Cervantes A, Estrada-Reyes ZM, Ogunade IM. Residual feed intake in beef cattle is associated with differences in hepatic mRNA expression of fatty acid, amino acid and mitochondrial energy metabolism genes. Front Anim Sci. 2022;3:838591.
21. Keogh K McKenna C Waters SM Porter RK Fitzsimons C McGee M Effect of breed and diet on the M. Longissimus Thoracis et lumborum transcriptome of steers divergent for residual feed intake Sci Rep Sci Rep 2023 13 9034 10.1038/s41598-023-35661-z 37270611
Keogh K, McKenna C, Waters SM, Porter RK, Fitzsimons C, McGee M, et al. Effect of breed and diet on the M. Longissimus Thoracis et lumborum transcriptome of steers divergent for residual feed intake. Sci Rep Sci Rep. 2023;13:9034.37270611
22. Weber KL Welly BT Van Eenennaam AL Young AE Porto-Neto LR Reverter A Identification of gene networks for residual feed intake in Angus cattle using genomic prediction and RNA-seq PLoS ONE 2016 11 e0152274 10.1371/journal.pone.0152274 27019286
Weber KL, Welly BT, Van Eenennaam AL, Young AE, Porto-Neto LR, Reverter A, et al. Identification of gene networks for residual feed intake in Angus cattle using genomic prediction and RNA-seq. PLoS ONE. 2016;11:e0152274.27019286
23. Alexandre PA Naval-Sanchez M Porto-Neto LR Ferraz JBS Reverter A Fukumasu H Systems Biology reveals NR2F6 and TGFB1 as key regulators of feed efficiency in beef cattle Front Genet 2019 10 230 10.3389/fgene.2019.00230 30967894
Alexandre PA, Naval-Sanchez M, Porto-Neto LR, Ferraz JBS, Reverter A, Fukumasu H. Systems Biology reveals NR2F6 and TGFB1 as key regulators of feed efficiency in beef cattle. Front Genet. 2019;10:230.30967894
24. Chen C Fang S Wei H He M Fu H Xiong X Prevotella copri increases fat accumulation in pigs fed with formula diets Microbiome 2021 9 175 10.1186/s40168-021-01110-0 34419147
Chen C, Fang S, Wei H, He M, Fu H, Xiong X, et al. Prevotella copri increases fat accumulation in pigs fed with formula diets. Microbiome. 2021;9:175.34419147
25. Iljazovic A Roy U Gálvez EJC Lesker TR Zhao B Gronow A Perturbation of the gut microbiome by Prevotella spp. enhances host susceptibility to mucosal inflammation Mucosal Immunol 2021 14 113 24 10.1038/s41385-020-0296-4 32433514
Iljazovic A, Roy U, Gálvez EJC, Lesker TR, Zhao B, Gronow A, et al. Perturbation of the gut microbiome by Prevotella spp. enhances host susceptibility to mucosal inflammation. Mucosal Immunol. 2021;14:113–24.32433514
26. Jiang Y Ogunade IM Arriola KG Qi M Vyas D Staples CR Effects of the dose and viability of Saccharomyces cerevisiae. 2. Ruminal fermentation, performance of lactating dairy cows, and correlations between ruminal bacteria abundance and performance measures J Dairy Sci 2017 100 8102 18 10.3168/jds.2016-12371 28822545
Jiang Y, Ogunade IM, Arriola KG, Qi M, Vyas D, Staples CR, et al. Effects of the dose and viability of Saccharomyces cerevisiae. 2. Ruminal fermentation, performance of lactating dairy cows, and correlations between ruminal bacteria abundance and performance measures. J Dairy Sci. 2017;100:8102–18.28822545
27. Jiang H Fang S Yang H Chen C Identification of the relationship between the gut microbiome and feed efficiency in a commercial pig cohort J Anim Sci 2021 99 skab045 10.1093/jas/skab045 33570553
Jiang H, Fang S, Yang H, Chen C. Identification of the relationship between the gut microbiome and feed efficiency in a commercial pig cohort. J Anim Sci. 2021;99:skab045.33570553
28. Jiang Q Xie C Chen L Xiao H Xie Z Zhu X Identification of gut microbes associated with feed efficiency by daily-phase feeding strategy in growing-finishing pigs Anim Nutr 2023 12 42 53 10.1016/j.aninu.2022.09.005 36381065
Jiang Q, Xie C, Chen L, Xiao H, Xie Z, Zhu X, et al. Identification of gut microbes associated with feed efficiency by daily-phase feeding strategy in growing-finishing pigs. Anim Nutr. 2023;12:42–53.36381065
29. Brooke CG Najafi N Dykier KC Hess M Prevotella copri, a potential indicator for high feed efficiency in western steers Anim Sci J 2019 90 696 701 10.1111/asj.13197 30848016
Brooke CG, Najafi N, Dykier KC, Hess M. Prevotella copri, a potential indicator for high feed efficiency in western steers. Anim Sci J. 2019;90:696–701.30848016
30. Lopes DRG de Souza Duarte M La Reau AJ Chaves IZ de Oliveira Mendes TA Detmann E Assessing the relationship between the rumen microbiota and feed efficiency in Nellore steers J Anim Sci Biotechnol 2021 12 79 10.1186/s40104-021-00599-7 34261531
Lopes DRG, de Souza Duarte M, La Reau AJ, Chaves IZ, de Oliveira Mendes TA, Detmann E, et al. Assessing the relationship between the rumen microbiota and feed efficiency in Nellore steers. J Anim Sci Biotechnol. 2021;12:79.34261531
31. Hernandez-Sanabria E Goonewardene LA Wang Z Durunna ON Moore SS Guan LL Impact of feed efficiency and diet on adaptive variations in the bacterial community in the rumen fluid of cattle Appl Environ Microbiol 2012 78 1203 14 10.1128/AEM.05114-11 22156428
Hernandez-Sanabria E, Goonewardene LA, Wang Z, Durunna ON, Moore SS, Guan LL. Impact of feed efficiency and diet on adaptive variations in the bacterial community in the rumen fluid of cattle. Appl Environ Microbiol. 2012;78:1203–14.22156428
32. McCann JC Wiley LM Forbes TD Rouquette FM Jr Tedeschi LO Relationship between the rumen microbiome and residual feed intake-efficiency of Brahman bulls stocked on bermudagrass pastures PLoS ONE 2014 9 e91864 10.1371/journal.pone.0091864 24642871
McCann JC, Wiley LM, Forbes TD, Rouquette FM Jr, Tedeschi LO. Relationship between the rumen microbiome and residual feed intake-efficiency of Brahman bulls stocked on bermudagrass pastures. PLoS ONE. 2014;9:e91864.24642871
33. Ellison MJ Conant GC Lamberson WR Cockrum RR Austin KJ Rule DC Diet and feed efficiency status affect rumen microbial profiles of sheep Small Rumin Res 2017 156 12 9 10.1016/j.smallrumres.2017.08.009
Ellison MJ, Conant GC, Lamberson WR, Cockrum RR, Austin KJ, Rule DC, et al. Diet and feed efficiency status affect rumen microbial profiles of sheep. Small Rumin Res. 2017;156:12–9.
34. Delgado B Bach A Guasch I González C Elcoso G Pryce JE Whole rumen metagenome sequencing allows classifying and predicting feed efficiency and intake levels in cattle Sci Rep 2019 9 11 10.1038/s41598-018-36673-w 30626904
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.30626904
35. Parra MC Costa DF Meale SJ Silva LFP Rumen bacteria and feed efficiency of beef cattle fed diets with different protein content Anim Prod Sci 2022 62 1029 39 10.1071/AN21508
Parra MC, Costa DF, Meale SJ, Silva LFP. Rumen bacteria and feed efficiency of beef cattle fed diets with different protein content. Anim Prod Sci. 2022;62:1029–39.
36. Zhou M Ghoshal B Stothard P Guan LL Distinctive roles between rumen epimural and content bacterial communities on beef cattle feed efficiency: a combined analysis Curr Res Microb Sci 2021 2 100085 34934993
Zhou M, Ghoshal B, Stothard P, Guan LL. Distinctive roles between rumen epimural and content bacterial communities on beef cattle feed efficiency: a combined analysis. Curr Res Microb Sci. 2021;2:100085.34934993
37. Yang H Yang M Fang S Huang X He M Ke S Evaluating the profound effect of gut microbiome on host appetite in pigs BMC Microbiol 2018 18 215 10.1186/s12866-018-1364-8 30547751
Yang H, Yang M, Fang S, Huang X, He M, Ke S, et al. Evaluating the profound effect of gut microbiome on host appetite in pigs. BMC Microbiol. 2018;18:215.30547751
38. McCormack UM Curião T Buzoianu SG Prieto ML Ryan T Varley P Exploring a possible link between the intestinal microbiota and feed Efficiency in pigs Appl Environ Microbiol 2017 83 e00380 17 10.1128/AEM.00380-17 28526795
McCormack UM, Curião T, Buzoianu SG, Prieto ML, Ryan T, Varley P, et al. Exploring a possible link between the intestinal microbiota and feed Efficiency in pigs. Appl Environ Microbiol. 2017;83:e00380–17.28526795
39. Kubasova T, Davidova-Gerzova L, Babak V, Cejkova D, Montagne L et al. Le-Floc’h N,. Effects of host genetics and environmental conditions on fecal microbiota composition of pigs. PLoS One. 2018;13:e0201901.
40. Torres-Pitarch A Gardiner GE Cormican P Rea M Crispie F O’Doherty JV Effect of cereal fermentation and carbohydrase supplementation on growth, nutrient digestibility and intestinal microbiota in liquid-fed grow-finishing pigs Sci Rep 2020 10 13716 10.1038/s41598-020-70443-x 32792575
Torres-Pitarch A, Gardiner GE, Cormican P, Rea M, Crispie F, O’Doherty JV, et al. Effect of cereal fermentation and carbohydrase supplementation on growth, nutrient digestibility and intestinal microbiota in liquid-fed grow-finishing pigs. Sci Rep. 2020;10:13716.32792575
41. Zhang J Shi H Wang Y Cao Z Yang H Li S Effect of limit-fed diets with different forage to concentrate ratios on fecal bacterial and archaeal community composition in Holstein heifers Front Microbiol 2018 9 976 10.3389/fmicb.2018.00976 29867879
Zhang J, Shi H, Wang Y, Cao Z, Yang H, Li S. Effect of limit-fed diets with different forage to concentrate ratios on fecal bacterial and archaeal community composition in Holstein heifers. Front Microbiol. 2018;9:976.29867879
42. 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 10.1186/s40168-019-0618-5 30642389
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.30642389
43. McKenna C Keogh K Porter RK Waters SM Cormican P Kenny DA An examination of skeletal muscle and hepatic tissue transcriptomes from beef cattle divergent for residual feed intake Sci Rep 2021 11 8942 10.1038/s41598-021-87842-3 33903612
McKenna C, Keogh K, Porter RK, Waters SM, Cormican P, Kenny DA. An examination of skeletal muscle and hepatic tissue transcriptomes from beef cattle divergent for residual feed intake. Sci Rep. 2021;11:8942.33903612
44. McLoughlin S Spillane C Claffey N Smith PE O’Rourke T Diskin MG Rumen Microbiome composition is altered in Sheep Divergent in feed efficiency Front Microbiol 2020 11 1981 10.3389/fmicb.2020.01981 32983009
McLoughlin S, Spillane C, Claffey N, Smith PE, O’Rourke T, Diskin MG, et al. Rumen Microbiome composition is altered in Sheep Divergent in feed efficiency. Front Microbiol. 2020;11:1981.32983009
45. McCormack UM Curião T Metzler-Zebeli BU Magowan E Berry DP Reyer H Porcine feed efficiency-Associated Intestinal Microbiota and physiological traits: finding consistent cross-locational biomarkers for residual feed intake mSystems 2019 4 e00324 18 10.1128/mSystems.00324-18 31213524
McCormack UM, Curião T, Metzler-Zebeli BU, Magowan E, Berry DP, Reyer H, et al. Porcine feed efficiency-Associated Intestinal Microbiota and physiological traits: finding consistent cross-locational biomarkers for residual feed intake. mSystems. 2019;4:e00324–18.31213524
46. Yan W Sun C Yuan J Yang N Gut metagenomic analysis reveals prominent roles of Lactobacillus and cecal microbiota in chicken feed efficiency Sci Rep 2017 7 45308 10.1038/srep45308 28349946
Yan W, Sun C, Yuan J, Yang N. Gut metagenomic analysis reveals prominent roles of Lactobacillus and cecal microbiota in chicken feed efficiency. Sci Rep. 2017;7:45308.28349946
47. Freetly HC Dickey A Lindholm-Perry AK Thallman RM Keele JW Foote AP Digestive tract microbiota of beef cattle that differed in feed efficiency J Anim Sci 2020 98 2 skaa008 10.1093/jas/skaa008 31930312
Freetly HC, Dickey A, Lindholm-Perry AK, Thallman RM, Keele JW, Foote AP, et al. Digestive tract microbiota of beef cattle that differed in feed efficiency. J Anim Sci. 2020;98(2):skaa008.31930312
48. Auffret MD Stewart RD Dewhurst RJ Duthie CA Watson M Roehe R Identification of microbial genetic capacities and potential mechanisms within the rumen microbiome explaining differences in beef cattle feed efficiency Front Microbiol 2020 11 1229 10.3389/fmicb.2020.01229 32582125
Auffret MD, Stewart RD, Dewhurst RJ, Duthie CA, Watson M, Roehe R, et al. Identification of microbial genetic capacities and potential mechanisms within the rumen microbiome explaining differences in beef cattle feed efficiency. Front Microbiol. 2020;11:1229.32582125
49. Manzanares-Miranda N Garcia-Mazcorro JF Pérez-Medina EB Vaquera-Vázquez A Martínez-Ruiz A Ramos-Zayas Y Microbial populations in ruminal liquid samples from young beefmaster bulls at both extremes of RFI values Microorganisms 2023 11 663 10.3390/microorganisms11030663 36985235
Manzanares-Miranda N, Garcia-Mazcorro JF, Pérez-Medina EB, Vaquera-Vázquez A, Martínez-Ruiz A, Ramos-Zayas Y, et al. Microbial populations in ruminal liquid samples from young beefmaster bulls at both extremes of RFI values. Microorganisms. 2023;11:663.36985235
50. McCann JC. 2013 Effect of post-extraction algla residue supplementation on the rumen microbiome of steers consuming low-quality forage [M.S. thesis]. College Station, TX: Texas A&M University.
51. McCabe MS Cormican P Keogh K O’Connor A O’Hara E Palladino RA Illumina MiSeq Phylogenetic Amplicon sequencing shows a large reduction of an uncharacterised Succinivibrionaceae and an increase of the Methanobrevibacter gottschalkii clade in feed restricted cattle PLoS ONE 2015 10 e0133234 10.1371/journal.pone.0133234 26226343
McCabe MS, Cormican P, Keogh K, O’Connor A, O’Hara E, Palladino RA, et al. Illumina MiSeq Phylogenetic Amplicon sequencing shows a large reduction of an uncharacterised Succinivibrionaceae and an increase of the Methanobrevibacter gottschalkii clade in feed restricted cattle. PLoS ONE. 2015;10:e0133234.26226343
52. Luo D Gao Y Lu Y Qu M Xiong X Xu L Niacin alters the ruminal microbial composition of cattle under high-concentrate condition Anim Nutr 2017 3 180 5 10.1016/j.aninu.2017.04.005 29767159
Luo D, Gao Y, Lu Y, Qu M, Xiong X, Xu L, et al. Niacin alters the ruminal microbial composition of cattle under high-concentrate condition. Anim Nutr. 2017;3:180–5.29767159
53. McLoughlin S Spillane C Campion FP Claffey N Sosa CC McNicholas Y Smith PE Diskin MG Waters SM Breed and ruminal fraction effects on bacterial and archaeal community composition in sheep Sci Rep 2023 13 3336 10.1038/s41598-023-28909-1 36849493
McLoughlin S, Spillane C, Campion FP, Claffey N, Sosa CC, McNicholas Y, Smith PE, Diskin MG, Waters SM. Breed and ruminal fraction effects on bacterial and archaeal community composition in sheep. Sci Rep. 2023;13:3336.36849493
54. Coyle S Fitzsimons C Kenny DA Kelly AK McGee M Feed efficiency correlations in beef cattle offered a zero-grazed grass and a high concentrate diet Adv Anim Biosci 2017 8 121
Coyle S, Fitzsimons C, Kenny DA, Kelly AK, McGee M. Feed efficiency correlations in beef cattle offered a zero-grazed grass and a high concentrate diet. Adv Anim Biosci. 2017;8:121.
55. Coyle S Fitzsimons C Kenny DA Kelly AK McGee M Repeatability of feed efficiency in steers offered a high-concentrate diet J Anim Sci 2016 94 supplement5 719
Coyle S, Fitzsimons C, Kenny DA, Kelly AK, McGee M. Repeatability of feed efficiency in steers offered a high-concentrate diet. J Anim Sci. 2016;94(supplement5):719.
56. Yu Z Morrison M Improved extraction of PCR-quality community DNA from digesta and fecal samples Biotechniques 2004 36 808 12 10.2144/04365ST04 15152600
Yu Z, Morrison M. Improved extraction of PCR-quality community DNA from digesta and fecal samples. Biotechniques. 2004;36:808–12.15152600
57. Bolyen E, Ram Rideout J, Dillon MR, Bokulich NA, Abnet C, Al-Ghalith GA, et al. QIIME 2: reproducible, interactive, scalable, and extensible microbiome data science. Peer Preprints; 2018.
58. Callahan BJ McMurdie PJ Rosen MJ Han AW Johnson AJ Holmes SP DADA2: high-resolution sample inference from Illumina amplicon data Nat Methods 2016 13 581 3 10.1038/nmeth.3869 27214047
Callahan BJ, McMurdie PJ, Rosen MJ, Han AW, Johnson AJ, Holmes SP. DADA2: high-resolution sample inference from Illumina amplicon data. Nat Methods. 2016;13:581–3.27214047
59. Callahan BJ McMurdie PJ Holmes SP Exact sequence variants should replace operational taxonomic units in marker-gene data analysis ISME J 2017 11 2639 10.1038/ismej.2017.119 28731476
Callahan BJ, McMurdie PJ, Holmes SP. Exact sequence variants should replace operational taxonomic units in marker-gene data analysis. ISME J. 2017;11:2639.28731476
60. O’Leary NA Wright MW Brister JR Ciufo S Haddad D McVeigh R Reference sequence (RefSeq) database at NCBI: current status, taxonomic expansion, and functional annotation Nucleic Acids Res 2016 44 D733 45 10.1093/nar/gkv1189 26553804
O’Leary NA, Wright MW, Brister JR, Ciufo S, Haddad D, McVeigh R, et al. Reference sequence (RefSeq) database at NCBI: current status, taxonomic expansion, and functional annotation. Nucleic Acids Res. 2016;44:D733–45.26553804
61. McCarthy SD Butler ST Patton J Daly M Morris DG Kenny DA Differences in the expression of genes involved in the somatotropic axis in divergent strains of holstein-friesian dairy cows during early and mid lactation J Dairy Sci 2009 92 5229 38 10.3168/jds.2008-1856 19762841
McCarthy SD, Butler ST, Patton J, Daly M, Morris DG, Kenny DA, et al. Differences in the expression of genes involved in the somatotropic axis in divergent strains of holstein-friesian dairy cows during early and mid lactation. J Dairy Sci. 2009;92:5229–38.19762841
62. Andrews S. FastQC: a quality control tool for high throughput sequence data. 2010. [ https://www.bioinformatics.babraham.ac.uk/projects/fastqc/].
63. Martin M Cutadapt removes adapter sequences from high-throughput sequencing reads EMBnet J 2011 17 1 10.14806/ej.17.1.200
Martin M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet J. 2011;17:1.
64. Rosen BD Bickhart DM Schnabel RD Koren S Elsik CG Tseng E De novo assembly of the cattle reference genome with single-molecule sequencing Gigascience 2020 9 giaa021 10.1093/gigascience/giaa021 32191811
Rosen BD, Bickhart DM, Schnabel RD, Koren S, Elsik CG, Tseng E, et al. De novo assembly of the cattle reference genome with single-molecule sequencing. Gigascience. 2020;9:giaa021.32191811
65. Dobin A Davis CA Schlesinger F Drenkow J Zaleski C Jha S STAR: ultrafast universal RNA-seq aligner Bioinformatics 2013 29 15 21 10.1093/bioinformatics/bts635 23104886
Dobin A, Davis CA, Schlesinger F, Drenkow J, Zaleski C, Jha S, et al. STAR: ultrafast universal RNA-seq aligner. Bioinformatics. 2013;29:15–21.23104886
66. Robinson MD McCarthy DJ Smyth GK edgeR: a Bioconductor package for differential expression analysis of digital gene expression data Bioinformatics 2010 26 139 40 10.1093/bioinformatics/btp616 19910308
Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26:139–40.19910308
67. Reverter A Chan EKF Combining partial correlation and an information theory approach to the reversed engineering of gene co-expression networks Bioinformatics 2008 24 2491 7 10.1093/bioinformatics/btn482 18784117
Reverter A, Chan EKF. Combining partial correlation and an information theory approach to the reversed engineering of gene co-expression networks. Bioinformatics. 2008;24:2491–7.18784117
68. Shannon P Markiel A Ozier O Baliga NS Wang JT Ramage D Cytoscape: a software environment for integrated models of biomolecular interaction networks Genome Res 2003 13 2498 504 10.1101/gr.1239303 14597658
Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13:2498–504.14597658
69. Caraux G Pinloche S PermutMatrix: a graphical environment to arrange gene expression profiles in optimal linear order Bioinformatics 2005 21 1280 1 10.1093/bioinformatics/bti141 15546938
Caraux G, Pinloche S. PermutMatrix: a graphical environment to arrange gene expression profiles in optimal linear order. Bioinformatics. 2005;21:1280–1.15546938
