Transcriptome analyses reveal reduced hepatic lipid synthesis and accumulation in more feed efficient beef cattle

Relationship between the rumen microbiome and liver transcriptome in beef cattle divergent for feed efficiency

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.

Effect of breed and dietary composition on the miRNA profile of beef steers divergent for feed efficiency

Identifying and breeding cattle that are more feed efficient is of great benefit to beef production. Additionally, it is crucial that genes contributing to feed efficiency are robust across varying management settings including dietary source as well as being relevant across contrasting breeds of cattle. The aim of this study was to determine miRNAs that are contributing to the expression of residual feed intake (RFI) across two breeds and dietary sources. miRNA profiling was undertaken in Longissimus dorsi tissue of Charolais and Holstein-Friesian steers divergent for RFI phenotype following two contrasting consecutive diets (high-forage and high-concentrate). Ten miRNA were identified as differentially expressed (adj. P < 0.1) across the breed and diet contrasts examined. Of particular interest was the differential expression of miR-2419-5p and miR-2415-3p, both of which were up-regulated in the Low-RFI Charolais steers across each dietary phase. Pathway analysis of target mRNA genes of differentially expressed miRNA revealed enrichment (P < 0.05) for pathways including metabolic related pathways, insulin receptor signalling, adipogenesis as well as pathways related to skeletal muscle growth. These results provide insight into the skeletal muscle miRNAome of beef cattle and their potential molecular regulatory mechanisms relating to feed efficiency in beef cattle.

Characterization of rumen microbiome and immune genes expression of crossbred beef steers with divergent residual feed intake phenotypes

We investigated whole blood and hepatic mRNA expressions of immune genes and rumen microbiome of crossbred beef steers with divergent residual feed intake phenotype to identify relevant biological processes underpinning feed efficiency in beef cattle. Low-RFI beef steers (n = 20; RFI = - 1.83 kg/d) and high-RFI beef steers (n = 20; RFI = + 2.12 kg/d) were identified from a group of 108 growing crossbred beef steers (average BW = 282 ± 30.4 kg) fed a high-forage total mixed ration after a 70-d performance testing period. At the end of the 70-d testing period, liver biopsies and blood samples were collected for total RNA extraction and cDNA synthesis. Rumen fluid samples were also collected for analysis of the rumen microbial community. The mRNA expression of 84 genes related to innate and adaptive immunity was analyzed using pathway-focused PCR-based arrays. Differentially expressed genes were determined using P-value ≤ 0.05 and fold change (FC) ≥ 1.5 (in whole blood) or ≥ 2.0 (in the liver). Gene ontology analysis of the differentially expressed genes revealed that pathways related to pattern recognition receptor activity, positive regulation of phagocytosis, positive regulation of vitamin metabolic process, vascular endothelial growth factor production, positive regulation of epithelial tube formation and T-helper cell differentiation were significantly enriched (FDR < 0.05) in low-RFI steers. In the rumen, the relative abundance of PeH15, Arthrobacter, Moryella, Weissella, and Muribaculaceae was enriched in low-RFI steers, while Methanobrevibacter, Bacteroidales_BS11_gut_group, Bacteroides and Clostridium_sensu_stricto_1 were reduced. In conclusion, our study found that low-RFI beef steers exhibit increased mRNA expression of genes related to immune cell functions in whole blood and liver tissues, specifically those involved in pathogen recognition and phagocytosis regulation. Additionally, these low-RFI steers showed differences in the relative abundance of some microbial taxa which may partially account for their improved feed efficiency compared to high-RFI steers.

An across breed, diet and tissue analysis reveals the transcription factor NR1H3 as a key mediator of residual feed intake in beef cattle

BACKGROUND

Provision of feed is a major determinant of overall profitability in beef production systems, accounting for up to 75% of the variable costs. Thus, improving cattle feed efficiency, by way of determining the underlying genomic control and subsequently selecting for feed efficient cattle, provides a method through which feed input costs may be reduced. The objective of this study was to undertake gene co-expression network analysis using RNA-Sequence data generated from Longissimus dorsi and liver tissue samples collected from steers of two contrasting breeds (Charolais and Holstein-Friesian) divergent for residual feed intake (RFI), across two consecutive distinct dietary phases (zero-grazed grass and high-concentrate). Categories including differentially expressed genes (DEGs) based on the contrasts of RFI phenotype, breed and dietary source, as well as key transcription factors and proteins secreted in plasma were utilised as nodes of the gene co-expression network.

RESULTS

Of the 2,929 DEGs within the network analysis, 1,604 were reported to have statistically significant correlations (≥ 0.80), resulting in a total of 43,876 significant connections between genes. Pathway analysis of clusters of co-expressed genes revealed enrichment of processes related to lipid metabolism (fatty acid biosynthesis, fatty acid β-oxidation, cholesterol biosynthesis), immune function, (complement cascade, coagulation system, acute phase response signalling), and energy production (oxidative phosphorylation, mitochondrial L-carnitine shuttle pathway) based on genes related to RFI, breed and dietary source contrasts.

CONCLUSIONS

Although similar biological processes were evident across the three factors examined, no one gene node was evident across RFI, breed and diet contrasts in both liver and muscle tissues. However within the liver tissue, the IRX4, NR1H3, HOXA13 and ZNF648 gene nodes, which all encode transcription factors displayed significant connections across the RFI, diet and breed comparisons, indicating a role for these transcription factors towards the RFI phenotype irrespective of diet and breed. Moreover, the NR1H3 gene encodes a protein secreted into plasma from the hepatocytes of the liver, highlighting the potential for this gene to be explored as a robust biomarker for the RFI trait in beef cattle.

Analyses of plasma metabolites using a high performance four-channel CIL LC-MS method and identification of metabolites associated with enteric methane emissions in beef cattle

Reducing enteric methane (one greenhouse gas) emissions from beef cattle not only can be beneficial in reducing global warming, but also improve efficiency of nutrient utilization in the production system. However, direct measurement of enteric methane emissions on individual cattle is difficult and expensive. The objective of this study was to detect plasma metabolites that are associated with enteric methane emissions in beef cattle. Average enteric methane emissions (CH4) per day (AVG_DAILYCH4) for each individual cattle were measured using the GreenFeed emission monitoring (GEM) unit system, and beef cattle with divergent AVG_DAILYCH4 from Angus (n = 10 for the low CH4 group and 9 for the high CH4 group), Charolais (n = 10 for low and 10 for = high), and Kinsella Composite (n = 10 for low and 10 for high) populations were used for plasma metabolite quantification and metabolite-CH4 association analyses. Blood samples of these cattle were collected near the end of the GEM system tests and a high performance four-channel chemical isotope labeling (CIL) liquid chromatography (LC) mass spectrometer (MS) method was applied to identify and quantify concentrations of metabolites. The four-channel CIL LC-MS method detected 4235 metabolites, of which 1105 were found to be significantly associated with AVG_DAILYCH4 by a t-test, while 1305 were significantly associated with AVG_DAILYCH4 by a regression analysis at p<0.05. Both the results of the t-test and regression analysis revealed that metabolites that were associated with enteric methane emissions in beef cattle were largely breed-specific whereas 4.29% to 6.39% CH4 associated metabolites were common across the three breed populations and 11.07% to 19.08% were common between two breed populations. Pathway analyses of the CH4 associated metabolites identified top enriched molecular processes for each breed population, including arginine and proline metabolism, arginine biosynthesis, butanoate metabolism, and glutathione metabolism for Angus; beta-alanine metabolism, pyruvate metabolism, glycolysis / gluconeogenesis, and citrate cycle (TCA cycle) for Charolais; phenylalanine, tyrosine and tryptophan biosynthesis, phenylalanine metabolism, arginine biosynthesis, and arginine and proline metabolism for Kinsella Composite. The detected CH4 associated metabolites and enriched molecular processes will help understand biological mechanisms of enteric methane emissions in beef cattle. The detected CH4 associated plasma metabolites will also provide valuable resources to further characterize the metabolites and verify their utility as biomarkers for selection of cattle with reduced methane emissions.

Analysis of liver miRNA in Hu sheep with different residual feed intake

Feed efficiency (FE), an important economic trait in sheep production, is indirectly assessed by residual feed intake (RFI). However, RFI in sheep is varied, and the molecular processes that regulate RFI are unclear. It is thus vital to investigate the molecular mechanism of RFI to developing a feed-efficient sheep. The miRNA-sequencing (RNA-Seq) was utilized to investigate miRNAs in liver tissue of 6 out of 137 sheep with extreme RFI phenotypic values. In these animals, as a typical metric of FE, RFI was used to distinguish differentially expressed miRNAs (DE_miRNAs) between animals with high (n = 3) and low (n = 3) phenotypic values. A total of 247 miRNAs were discovered in sheep, with four differentially expressed miRNAs (DE_miRNAs) detected. Among these DE_miRNAs, three were found to be upregulated and one was downregulated in animals with low residual feed intake (Low_RFI) compared to those with high residual feed intake (High_RFI). The target genes of DE_miRNAs were primarily associated with metabolic processes and biosynthetic process regulation. Furthermore, they were also considerably enriched in the FE related to glycolysis, protein synthesis and degradation, and amino acid biosynthesis pathways. Six genes were identified by co-expression analysis of DE_miRNAs target with DE_mRNAs. These results provide a theoretical basis for us to understand the sheep liver miRNAs in RFI molecular regulation.

Characterization of feed efficiency-related key signatures molecular in different cattle breeds

Feed efficiency is a major constraint in the beef industry and has a significant negative correlation with residual feed intake (RFI). RFI is widely used as a measure of feed efficiency in beef cattle and is independent of economic traits such as body weight and average daily gain. However, key traits with commonality or specificity among beef cattle breeds at the same level of RFI have not been reported. Accordingly, the present study hypothesized that signatures associated with feed efficiency would have commonality or specificity in the liver of cattle breeds at the same RFI level. By comparing and integrating liver transcriptome data, we investigated the critical signatures closely associated with RFI in beef cattle using weighted co-expression network analysis, consensus module analysis, functional enrichment analysis and protein network interaction analysis. The results showed that the consensus modules in Angus and Charolais cattle were negatively correlated, and four (turquoise, red, tan, yellow) were significantly positively correlated in Angus liver, while (turquoise, red) were significantly negatively correlated in Charolais liver. These consensus modules were found to be primarily involved in biological processes such as substance metabolism, energy metabolism and gene transcription, which may be one of the possible explanations for the difference in feed efficiency between the two beef breeds. This research also identified five key candidate genes, PLA2G12B, LCAT, MTTP, LCAT, ABCA1 and FADS1, which are closely associated with hepatic lipid metabolism. The present study has identified some modules, genes and pathways that may be the major contributors to the variation in feed efficiency among different cattle breeds, providing a new perspective on the molecular mechanisms of feed efficiency in beef cattle and a research basis for investigating molecular markers associated with feed efficiency in beef cattle.

Genome wide transcriptome analysis provides bases on hepatic lipid metabolism disorder affected by increased dietary grain ratio in fattening lambs

BACKGROUND

The liver is a principal metabolic organ and has a major role in regulating lipid metabolism. With the development of rapidly fattening livestock in the modern breeding industry, the incidence of hepatic steatosis and accumulation in animals was significantly increased. However, the molecular mechanisms responsible for hepatic lipid metabolic disturbances in a high concentrate diet remain unclear. The objective of this study was to evaluate the effects of increasing concentrate level in a fattening lamb diet on biochemical indices, hepatic triglycerides (TG) concentration, and hepatic transcriptomic profiles. In the present study, 42 weaned lambs (about 3 ± 0.3 months old) were randomly assigned to the GN60 group (60% concentrate of dry matter, GN60, n = 21) or GN70 group (70% concentrate of dry matter, n = 21) for a 3-months feeding trial.

RESULTS

No difference was observed in the growth performance or plasma biochemical parameters between the GN60 group and the GN70 group. The hepatic TG concentration was higher in the GN70 group than GN60 group (P < 0.05). Hepatic transcriptomic analysis showed that there were 290 differentially expressed genes identified between GN60 and GN70 groups, with 125 genes up-regulated and 165 genes down-regulated in the GN70 group. The enriched Gene Ontology (GO) items and KEGG pathways and protein-protein interaction (PPI) network of differentially expressed genes (DEGs) revealed that the majority of enriched pathways were related to lipid metabolism. Further analysis revealed that the fatty acid synthesis was up-regulated, while fatty acid transport, oxidation, and TG degradation were down-regulated in the GN70 group when compared with the GN60 group.

CONCLUSIONS

These results indicated that GN70 induced excess lipid deposition in the liver of lambs during the fattening period, with high synthesis rates and low degradation rates of TG. The identified mechanisms may help understand hepatic metabolism in lambs with a high concentrate diet and provide insight into decreasing the risk of liver metabolism disorder in animals.

Effect of breed and diet on the M. longissimus thoracis et lumborum transcriptome of steers divergent for residual feed intake

Improving cattle feed efficiency through selection of residual feed intake (RFI) is a widely accepted approach to sustainable beef production. A greater understanding of the molecular control of RFI in various breeds offered contrasting diets is necessary for the accurate identification of feed efficient animals and will underpin accelerated genetic improvement of the trait. The aim of this study was to determine genes and biological processes contributing to RFI across varying breed type and dietary sources in skeletal muscle tissue. Residual feed intake was calculated in Charolais and Holstein-Friesian steers across multiple dietary phases (phase-1: high concentrate (growing-phase); phase-2: zero-grazed grass (growing-phase); phase-3: high concentrate (finishing-phase). Steers divergent for RFI within each breed and dietary phase were selected for muscle biopsy collection, and muscle samples subsequently subjected to RNAseq analysis. No gene was consistently differentially expressed across the breed and diet types examined. However, pathway analysis revealed commonality across breeds and diets for biological processes including fatty acid metabolism, immune function, energy production and muscle growth. Overall, the lack of commonality of individual genes towards variation in RFI both within the current study and compared to the published literature, suggests other genomic features warrant further evaluation in relation to RFI.

Residual Feed Intake in Beef Cattle Is Associated With Differences in Hepatic mRNA Expression of Fatty Acid, Amino Acid, and Mitochondrial Energy Metabolism Genes

We evaluated the mRNA expression of genes involved in hepatic fatty acid, amino acid, and mitochondrial energy metabolism in crossbred beef steers with divergent low and high residual feed intake (RFI). Low-RFI beef steers (n = 8; RFI = - 1.93 kg/d) and high-RFI beef steers (n = 8; RFI = + 2.01kg/d) were selected from a group of 56 growing crossbred beef steers (average BW = 261 ± 18.5 kg) fed a high-forage total mixed ration after a 49-d performance testing period. At the end of the 49-d performance testing period, liver biopsies were collected from the low-RFI and high-RFI beef steers for RNA extraction and cDNA synthesis. The mRNA expression of 84 genes each related to fatty acid metabolism, amino acid metabolism, and mitochondrial energy metabolism were analyzed using pathway-focused PCR-based arrays. The mRNA expression of 8 genes (CRAT, SLC27A5, SLC27A2, ACSBG2, ACADL, ACADSB, ACAA1, and ACAA2) involved fatty acid transport and β-oxidation were upregulated (FC ≥ 2.0, FDR ≤ 0.05) in low-RFI, compared to high-RFI steers. Among those involved in amino acid metabolism, hepatic mRNA expression of a gene encoding for aminoadipate aminotransferase, an enzyme related to lysine degradation, was downregulated (FC = -5.45, FDR = 0.01) in low-RFI steers, whereas those of methionine adenosyltransferase I and aspartate aminotransferase 2, which both link amino acid and lipid metabolism, were upregulated (FC ≥ 2.0, FDR ≤ 0.05). Two mitochondrial energy metabolism genes (UQCRC1 and ATP5G1) involved in ATP synthesis via oxidative phosphorylation were upregulated (FC ≥ 2.0, FDR ≤ 0.05) in low-RFI beef steers, compared to high-RFI beef steers. The results of this study demonstrated that low-RFI beef steers exhibit upregulation of molecular mechanisms related to fatty acid transport, fatty acid β-oxidation, and mitochondrial ATP synthesis, which suggest that low-RFI beef steers have enhanced metabolic capacity to maximize capture of energy and nutrients from feeds consumed.

Proteomic analysis of adipose tissue revealing differentially abundant proteins in highly efficient mid-lactating dairy cows

The improvement of nutrient utilization efficiency in dairy cows represents an important task in view of the current rising demand for animal products and sustainable resource usage. In this perspective, the identification of appropriate markers to identify the most efficient animals for dairy production becomes a crucial factor. Residual feed intake (RFI), which represents the difference between predicted and actual intake, is used to define the efficiency of cows. In this study, subcutaneous adipose tissue (AT) was collected from five high efficient (HEF) and five low efficient (LEF) mid-lactation Holstein dairy cows, that represented subgroups of the 20% lowest RFI values (HEF) and highest 20% RFI values (LEF), out of a cohort of 155 cows that were examined for feed efficiency at the individual dairy barn at Volcani Institute, Israel. Adipose samples were examined for proteomic analysis by nano-LC/MS–MS and gene expression by RT-PCR. A total of 101 differential proteins (P ≤ 0.05 and fold change ± 1.5) and two protein networks related to feed efficiency were found between HEF and LEF cows. Among the enriched top canonical pathways, FAT10 signaling, EIF2 signaling, Sirtuin signaling, Acute phase response signaling, Protein ubiquitination and mTOR signaling pathways were related to feed efficiency in AT. Furthermore, abundance of transferrin (TF; FC = 78.35, P = 0.02) enriched pathways, including mTOR signaling, LXR/RXR and FXR/RXR activation was found in AT of HEF cows. Relative mRNA expression of RBM39, which is involved in energy metabolism, was decreased in AT of HEF versus LEF. The relationship found between the AT proteins and/or metabolic pathways and the feed efficiency demonstrates that AT may reflect metabolic adaptations to high efficiency, and suggests that these proteins together with their metabolic mechanisms are suitable candidates as biomarkers to identify efficient cows for dairy production.

Identification of Key Pathways Associated With Residual Feed Intake of Beef Cattle Based on Whole Blood Transcriptome Data Analyzed Using Gene Set Enrichment Analysis

We applied whole blood transcriptome analysis and gene set enrichment analysis to identify pathways associated with divergent selection for low or high RFI in beef cattle. A group of 56 crossbred beef steers (average BW = 261 ± 18.5 kg) were adapted to a high-forage total mixed ration in a confinement dry lot equipped with GrowSafe intake nodes for period of 49 d to determine their residual feed intake (RFI). After RFI determination, whole blood samples were collected from beef steers with the lowest RFI (most efficient; low-RFI; n = 8) and highest RFI (least efficient; high-RFI; n = 8). Prior to RNA extraction, whole blood samples collected were composited for each steer. Sequencing was performed on an Illumina NextSeq2000 equipped with a P3 flow. Gene set enrichment analysis (GSEA) was used to analyze differentially expressed gene sets and pathways between the two groups of steers. Results of GSEA revealed pathways associated with metabolism of proteins, cellular responses to external stimuli, stress, and heat stress were differentially inhibited (false discovery rate (FDR) < 0.05) in high-RFI compared to low-RFI beef cattle, while pathways associated with binding and uptake of ligands by scavenger receptors, scavenging of heme from plasma, and erythrocytes release/take up oxygen were differentially enriched (FDR < 0.05) in high-RFI, relative to low-RFI beef cattle. Taken together, our results revealed that beef steers divergently selected for low or high RFI revealed differential expressions of genes related to protein metabolism and stress responsiveness.

Transcriptome analysis reveals the potential roles of long non-coding RNAs in feed efficiency of chicken

Feed efficiency is an important economic trait and reduces the production costs per unit of animal product. Up to now, few studies have conducted transcriptome profiling of liver tissue in feed efficiency-divergent chickens (Ross vs native breeds). Also, molecular mechanisms contributing to differences in feed efficiency are not fully understood, especially in terms of long non-coding RNAs (lncRNAs). Hence, transcriptome profiles of liver tissue in commercial and native chicken breeds were analyzed. RNA-Seq data along with bioinformatics approaches were applied and a series of lncRNAs and target genes were identified. Furthermore, protein-protein interaction network construction, co-expression analysis, co-localization analysis of QTLs and functional enrichment analysis were used to functionally annotate the identified lncRNAs. In total, 2,290 lncRNAs were found (including 1,110 annotated, 593 known and 587 novel), of which 53 (including 39 known and 14 novel), were identified as differentially expressed genes between two breeds. The expression profile of lncRNAs was validated by RT-qPCR. The identified novel lncRNAs showed a number of characteristics similar to those of known lncRNAs. Target prediction analysis showed that these lncRNAs have the potential to act in cis or trans mode. Functional enrichment analysis of the predicted target genes revealed that they might affect the differences in feed efficiency of chicken by modulating genes associated with lipid metabolism, carbohydrate metabolism, growth, energy homeostasis and glucose metabolism. Some gene members of significant modules in the constructed co-expression networks were reported as important genes related to feed efficiency. Co-localization analysis of QTLs related to feed efficiency and the identified lncRNAs suggested several candidates to be involved in residual feed intake. The findings of this study provided valuable resources to further clarify the genetic basis of regulation of feed efficiency in chicken from the perspective of lncRNAs.

Chemical Group-Based Metabolome Analysis Identifies Candidate Plasma Biomarkers Associated With Residual Feed Intake in Beef Steers

We applied chemical group-based metabolomics to identify blood metabolic signatures associated with residual feed intake in beef cattle. A group of 56 crossbred growing beef steers (average BW = 261.3 ± 18.5 kg) were adapted to a high-forage total mixed ration in a confinement dry lot equipped with GrowSafe intake nodes for period of 49 d to determine their residual feed intake classification (RFI). After RFI determination, weekly blood samples were collected three times from beef steers with the lowest RFI [most efficient (HFE); n = 8] and highest RFI and least-efficient [least efficient (LFE); n = 8]. Plasma was prepared by centrifugation and composited for each steer. Metabolome analysis was conducted using a chemical isotope labeling (CIL)/liquid chromatography–mass spectrometry, which permitted the analysis of metabolites containing amine/phenol-, carboxylic acid-, and carbonyl-chemical groups, which are metabolites associated with metabolisms of amino acids, fatty acids, and carbohydrates, respectively. A total number of 495 amine/phenol-containing metabolites were detected and identified; pathway analysis of all these metabolites showed that arginine biosynthesis and histidine metabolism were enriched (P &lt; 0.10) in HFE, relative to LFE steers. Biomarker analyses of the amine/phenol-metabolites identified methionine, 5-aminopentanoic acid, 2-aminohexanedioic acid, and 4-chlorolysine as candidate biomarkers of RFI [false discovery rate ≤ 0.05; Area Under the Curve (AUC) &gt; 0.90]. A total of 118 and 330 metabolites containing carbonyl- and carboxylic acid-chemical groups, respectively were detected and identified; no metabolic pathways associated with these metabolites were altered and only one candidate biomarker (methionine sulfoxide) was identified. These results identified five candidate metabolite biomarkers of RFI in beef cattle which are mostly associated with amino acid metabolism. Further validation using a larger cohort of beef cattle of different genetic pedigree is required to confirm these findings.

Identifying the key genes and functional enrichment pathways associated with feed efficiency in cattle

Residual feed intake (RFI) is a measurement of feed efficiency, and is inversely correlated with feed efficiency. The differentially expressed genes (DEGs) associated with RFI vary substantially among studies, posing great challenges in finding the RFI-related marker genes. This study attempted to resolve this issue by integrating and comparing the multiple transcriptome sequencing data associated with RFI in the cattle liver, using differential, functional enrichment, protein-protein interaction (PPI) network, weighted co-expression network (WGCNA), and gene set enrichment analyses (GSEA) to identify the candidate genes and functional enrichment pathways that are closely associated with RFI. Four candidate genes namely SHC1, GPX4, ACADL, and IGF1 were identified and validated as the marker genes for RFI. Four functional enrichment pathways, namely the fatty acid metabolism, sugar metabolism, energy metabolism, and protein ubiquitination were also found to be closely related to RFI. This study identified several genes and signaling pathways with shared characteristics, which will provide new insights into the molecular mechanisms related to the regulation of feed efficiency, and provide basis for molecular markers related to feed efficiency in beef cattle.

Non-invasive metabolomics biomarkers of production efficiency and beef carcass quality traits

The inter-cattle growth variations stem from the interaction of many metabolic processes making animal selection difficult. We hypothesized that growth could be predicted using metabolomics. Urinary biomarkers of cattle feed efficiency were explored using mass spectrometry-based untargeted and targeted metabolomics. Feed intake and weight-gain was measured in steers (n = 75) on forage-based growing rations (stage-1, 84 days) followed by high-concentrate finishing rations (stage-2, 84 days). Urine from days 0, 21, 42, 63, and 83 in each stage were analyzed from steers with the greater (n = 14) and least (n = 14) average-daily-gain (ADG) and comparable dry-matter-intake (DMI; within 0.32 SD of the mean). Steers were slaughtered after stage-2. Adjusted fat-thickness and carcass-yield-grade increased in greater-ADG-cattle selected in stage-1, but carcass traits did not differ between ADG-selected in stage-2. Overall 85 untargeted metabolites segregated greater- and least-ADG animals, with overlap across diets (both stages) and breed type, despite sampling time effects. Total 18-bile acids (BAs) and 5-steroids were quantified and associated with performance and carcass quality across ADG-classification depending on the stage. Stepwise logistic regression of urinary BA and steroids had > 90% accuracy identifying efficient-ADG-steers. Urine metabolomics provides new insight into the physiological mechanisms and potential biomarkers for feed efficiency.

Genome-wide association study for meat tenderness in beef cattle identifies patterns of the genetic contribution in different post-mortem stages

The beef tenderization process during the post-mortem period is one of the most important sensorial attributes and it is well-established. The aim of this study was to identify the genetic contribution pattern to meat tenderness at 7-(LMD7), 14-(LMD14), and 21-(LMD21) days post-mortem. The heritabilities for LMD7 (0.194), LMD14 (0.142) and LMD21 (0.048) are well established in the population evaluated here. However, its genetic contribution in terms of genomic candidate regions is still poorly understood. Tenderness was measured in the Longissiums thoracis using Warner-Bratzler shear force in the three post-mortem periods. A total of 4323 crossbred beef cattle were phenotyped and genotyped using the Illumina BovineSNP50K. The percentage of the total genetic variance was estimated using the weighted single-step genomic best linear unbiased prediction method. The main candidate windows for LMD7 were associated with proteolysis of myofibrillar structures and the weakening endomysium and perimysium. Candidate windows for LMD14 and LMD21 were mapped in bovine QTLs for body composition, height and growth. Results presented herein highlight, the largest contribution of proteolysis related processes before 14-days post-mortem and body composition characteristics in later stages for meat tenderness.

Identification of candidate genes and enriched biological functions for feed efficiency traits by integrating plasma metabolites and imputed whole genome sequence variants in beef cattle

BACKGROUND

Feed efficiency is one of the key determinants of beef industry profitability and sustainability. However, the cellular and molecular background behind feed efficiency is largely unknown. This study combines imputed whole genome DNA variants and 31 plasma metabolites to dissect genes and biological functions/processes that are associated with residual feed intake (RFI) and its component traits including daily dry matter intake (DMI), average daily gain (ADG), and metabolic body weight (MWT) in beef cattle.

RESULTS

Regression analyses between feed efficiency traits and plasma metabolites in a population of 493 crossbred beef cattle identified 5 (L-valine, lysine, L-tyrosine, L-isoleucine, and L-leucine), 4 (lysine, L-lactic acid, L-tyrosine, and choline), 1 (citric acid), and 4 (L-glutamine, glycine, citric acid, and dimethyl sulfone) plasma metabolites associated with RFI, DMI, ADG, and MWT (P-value < 0.1), respectively. Combining the results of metabolome-genome wide association studies using 10,488,742 imputed SNPs, 40, 66, 15, and 40 unique candidate genes were identified as associated with RFI, DMI, ADG, and MWT (P-value < 1 × 10-5), respectively. These candidate genes were found to be involved in some key metabolic processes including metabolism of lipids, molecular transportation, cellular function and maintenance, cell morphology and biochemistry of small molecules.

CONCLUSIONS

This study identified metabolites, candidate genes and enriched biological functions/processes associated with RFI and its component traits through the integrative analyses of metabolites with phenotypic traits and DNA variants. Our findings could enhance the understanding of biochemical mechanisms of feed efficiency traits and could lead to improvement of genomic prediction accuracy via incorporating metabolite data.

Characterization and Duodenal Transcriptome Analysis of Chinese Beef Cattle With Divergent Feed Efficiency Using RNA-Seq

Residual feed intake (RFI) is an important measure of feed efficiency for agricultural animals. Factors associated with cattle RFI include physiology, dietary factors, and the environment. However, a precise genetic mechanism underlying cattle RFI variations in duodenal tissue is currently unavailable. The present study aimed to identify the key genes and functional pathways contributing to variance in cattle RFI phenotypes using RNA sequencing (RNA-seq). Six bulls with extremely high or low RFIs were selected for detecting differentially expressed genes (DEGs) by RNA-seq, followed by conducting GO, KEGG enrichment, protein-protein interaction (PPI), and co-expression network (WGCNA, n = 10) analysis. A total of 380 differentially expressed genes was obtained from high and low RFI groups, including genes related to energy metabolism (ALDOA, HADHB, INPPL1), mitochondrial function (NDUFS1, RFN4, CUL1), and feed intake behavior (CCK). Two key sub-networks and 26 key genes were detected using GO analysis of DEGs and PPI analysis, such as TPM1 and TPM2, which are involved in mitochondrial pathways and protein synthesis. Through WGCNA, a gene network was built, and genes were sorted into 27 modules, among which the blue (r = 0.72, p = 0.03) and salmon modules (r = -0.87, p = 0.002) were most closely related with RFI. DEGs and genes from the main sub-networks and closely related modules were largely involved in metabolism; oxidative phosphorylation; glucagon, ribosome, and N-glycan biosynthesis, and the MAPK and PI3K-Akt signaling pathways. Through WGCNA, five key genes, including FN1 and TPM2, associated with the biological regulation of oxidative processes and skeletal muscle development were identified. Taken together, our data suggest that the duodenum has specific biological functions in regulating feed intake. Our findings provide broad-scale perspectives for identifying potential pathways and key genes involved in the regulation of feed efficiency in beef cattle.

Distinctive gene expression patterns and imprinting signatures revealed in reciprocal crosses between cattle sub-species

Background

There are two genetically distinct subspecies of cattle, Bos taurus taurus and Bos taurus indicus, which arose from independent domestication events. The two types of cattle show substantial phenotypic differences, some of which emerge during fetal development and are reflected in birth outcomes, including birth weight. We explored gene expression profiles in the placenta and four fetal tissues at mid-gestation from one taurine (Bos taurus taurus; Angus) and one indicine (Bos taurus indicus; Brahman) breed and their reciprocal crosses.

Results

In total 120 samples were analysed from a pure taurine breed, an indicine breed and their reciprocal cross fetuses, which identified 6456 differentially expressed genes (DEGs) between the two pure breeds in at least one fetal tissue of which 110 genes were differentially expressed in all five tissues examined. DEGs shared across tissues were enriched for pathways related to immune and stress response functions. Only the liver had a substantial number of DEGs when reciprocal crossed were compared among which 310 DEGs were found to be in common with DEGs identified between purebred livers; these DEGs were significantly enriched for metabolic process GO terms. Analysis of DEGs across purebred and crossbred tissues suggested an additive expression pattern for most genes, where both paternal and maternal alleles contributed to variation in gene expression levels. However, expression of 5% of DEGs in each tissue was consistent with parent of origin effects, with both paternal and maternal dominance effects identified.

Conclusions

These data identify candidate genes potentially driving the tissue-specific differences between these taurine and indicine breeds and provide a biological insight into parental genome effects underlying phenotypic differences in bovine fetal development.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07667-2.

An examination of skeletal muscle and hepatic tissue transcriptomes from beef cattle divergent for residual feed intake

The selection of cattle with enhanced feed efficiency is of importance with regard to reducing feed costs in the beef industry. Global transcriptome profiling was undertaken on liver and skeletal muscle biopsies from Simmental heifers and bulls divergent for residual feed intake (RFI), a widely acknowledged feed efficiency phenotype, in order to identify genes that may be associated with this trait. We identified 5 genes (adj. p < 0.1) to be differentially expressed in skeletal muscle between high and low RFI heifers with all transcripts involved in oxidative phosphorylation and mitochondrial homeostasis. A total of 11 genes (adj. p < 0. 1) were differentially expressed in liver tissue between high and low RFI bulls with differentially expressed genes related to amino and nucleotide metabolism as well as endoplasmic reticulum protein processing. No genes were identified as differentially expressed in either heifer liver or bull muscle analyses. Results from this study show that the molecular control of RFI in young cattle is modified according to gender, which may be attributable to differences in physiological maturity between heifers and bulls of the same age. Despite this we have highlighted a number of genes that may hold potential as molecular biomarkers for RFI cattle.

A hierarchical dynamic model used for investigating feed efficiency and its relationship with hepatic gene expression in APOE*3-Leiden.CETP mice

BACKGROUND

Feed efficiency (FE) is an important trait for livestock and humans. While the livestock industry focuses on increasing FE, in the current obesogenic society it is more of interest to decrease FE. Hence, understanding mechanisms involved in the regulation of FE and particularly how it can be decreased would help tremendously in counteracting the obesity pandemic. However, it is difficult to accurately measure or calculate FE in humans. In this study, we aimed to address this challenge by developing a hierarchical dynamic model based on humanized mouse data.

METHODS

We analyzed existing experimental data derived from 105 APOE*3-Leiden.CETP (E3L.CETP) mice fed a high-fat high-cholesterol (HFHC) diet for 1 (N = 20), 2 (N = 19), 3 (N = 20), and 6 (N = 46) month. We developed an ordinary differential equation (ODE) based model to estimate the FE based on the longitudinal data of body weight and food intake. Since the liver plays an important role in maintaining metabolic homeostasis, we evaluated associations between FE and hepatic gene expression levels. Depending on the feeding duration, we observed different relationships between FE and hepatic gene expression levels.

RESULTS

After 1-month feeding of HFHC diet, we observed that FE was associated with vitamin A metabolism, arachidonic acid metabolism, and the PPAR signaling pathway. After 3- and 6-month feeding of HFHC diet, we observed that FE was associated most strongly with expression levels of Spink1 and H19, genes involved in cell proliferation and glucose metabolism, respectively.

CONCLUSIONS

In conclusion, our analysis suggests that various biological processes such as vitamin A metabolism, hepatic response to inflammation, and cell proliferation associate with FE at different stages of diet-induced obesity.

Plasma Carboxyl-Metabolome Is Associated with Average Daily Gain Divergence in Beef Steers

We applied an untargeted metabolomics technique to analyze the plasma carboxyl-metabolome of beef steers with divergent average daily gain (ADG). Forty-eight newly weaned Angus crossbred beef steers were fed the same total mixed ration ad libitum for 42 days. On day 42, the steers were divided into two groups of lowest (LF: n = 8) and highest ADG (HF: n = 8), and blood samples were obtained from the two groups for plasma preparation. Relative quantification of carboxylic-acid-containing metabolites in the plasma samples was determined using a metabolomics technique based on chemical isotope labeling liquid chromatography mass spectrometry. Metabolites that differed (fold change (FC) ≥ 1.2 or ≤ 0.83 and FDR ≤ 0.05) between LF and HF were identified using a volcano plot. Metabolite set enrichment analysis (MSEA) of the differential metabolites was done to determine the metabolic pathways or enzymes that were potentially altered. In total, 328 metabolites were identified. Volcano plot analysis revealed 43 differentially abundant metabolites; several short chain fatty acids and ketone bodies had greater abundance in HF steers. Conversely, several long chain fatty acids were greater in LF steers. Five enzymatic pathways, such as fatty acyl CoA elongation and fatty-acid CoA ligase were altered based on MSEA. This study demonstrated that beef steers with divergent ADG had altered plasma carboxyl-metabolome, which is possibly caused by altered abundances and/or activities of enzymes involved in fatty acid oxidation and biosynthesis in the liver.

Genomic Heritability and Genome-Wide Association Studies of Plasma Metabolites in Crossbred Beef Cattle

Metabolites, substrates or products of metabolic processes, are involved in many biological functions, such as energy metabolism, signaling, stimulatory and inhibitory effects on enzymes and immunological defense. Metabolomic phenotypes are influenced by combination of genetic and environmental effects allowing for metabolome-genome-wide association studies (mGWAS) as a powerful tool to investigate the relationship between these phenotypes and genetic variants. The objectives of this study were to estimate genomic heritability and perform mGWAS and in silico functional enrichment analyses for a set of plasma metabolites in Canadian crossbred beef cattle. Thirty-three plasma metabolites and 45,266 single nucleotide polymorphisms (SNPs) were available for 475 animals. Genomic heritability for all metabolites was estimated using genomic best linear unbiased prediction (GBLUP) including genomic breed composition as covariates in the model. A single-step GBLUP implemented in BLUPF90 programs was used to determine SNP P values and the proportion of genetic variance explained by SNP windows containing 10 consecutive SNPs. The top 10 SNP windows that explained the largest genetic variation for each metabolite were identified and mapped to detect corresponding candidate genes. Functional enrichment analyses were performed on metabolites and their candidate genes using the Ingenuity Pathway Analysis software. Eleven metabolites showed low to moderate heritability that ranged from 0.09 ± 0.15 to 0.36 ± 0.15, while heritability estimates for 22 metabolites were zero or negligible. This result indicates that while variations in 11 metabolites were due to genetic variants, the majority are largely influenced by environment. Three significant SNP associations were detected for betaine (rs109862186), L-alanine (rs81117935), and L-lactic acid (rs42009425) based on Bonferroni correction for multiple testing (family wise error rate <0.05). The SNP rs81117935 was found to be located within the Catenin Alpha 2 gene (CTNNA2) showing a possible association with the regulation of L-alanine concentration. Other candidate genes were identified based on additive genetic variance explained by SNP windows of 10 consecutive SNPs. The observed heritability estimates and the candidate genes and networks identified in this study will serve as baseline information for research into the utilization of plasma metabolites for genetic improvement of crossbred beef cattle.

Bovine hepatic miRNAome profiling and differential miRNA expression analyses between beef steers with divergent feed efficiency phenotypes

MicroRNAs (miRNAs) are small RNA molecules involved in regulation of multiple biological processes through modulating expression of their target genes. Here we employed RNAseq to profile liver tissue miRNAome of 60 steers from Angus, Charolais, and Kinsella Composite (KC) populations. Of these animals, 36 animals (n = 12 for each breed) were utilized to identify differentially expressed (DE) miRNAs between animals with high (n = 6) or low (n = 6) phenotypic values of residual feed intake (RFI), a common measurement of feed efficiency. At a threshold of fold-change > 1.5 and P-value < 0.05, we detected 12 (7 up- and 5 downregulated in low-RFI animals), 18 (12 up- and 6 downregulated), and 13 (8 up- and 5 downregulated) DE miRNAs for Angus, Charolais, and KC steers, respectively. Most of the DE miRNAs were breed specific, with bta-miR-449a and bta-miR-AB-2 being differentially expressed in all three breeds. The predicted target genes of the identified DE miRNA are mainly involved in cell cycle, cell death and survival, cell signaling, cellular growth and proliferation, protein trafficking, cell morphology, cell-to-cell signaling and interaction, cellular development, molecular transport, post-translational modification, as well as nutrient metabolism (lipids, carbohydrates, protein and amino acid). Our results provide insights into the bovine hepatic miRNAome and their potential roles in molecular regulation of RFI in beef cattle.

Exploring the Regulatory Potential of Long Non-Coding RNA in Feed Efficiency of Indicine Cattle

Long non-coding RNA (lncRNA) can regulate several aspects of gene expression, being associated with complex phenotypes in humans and livestock species. In taurine beef cattle, recent evidence points to the involvement of lncRNA in feed efficiency (FE), a proxy for increased productivity and sustainability. Here, we hypothesized specific regulatory roles of lncRNA in FE of indicine cattle. Using RNA-Seq data from the liver, muscle, hypothalamus, pituitary gland and adrenal gland from Nellore bulls with divergent FE, we submitted new transcripts to a series of filters to confidently predict lncRNA. Then, we identified lncRNA that were differentially expressed (DE) and/or key regulators of FE. Finally, we explored lncRNA genomic location and interactions with miRNA and mRNA to infer potential function. We were able to identify 126 relevant lncRNA for FE in Bos indicus, some with high homology to previously identified lncRNA in Bos taurus and some possible specific regulators of FE in indicine cattle. Moreover, lncRNA identified here were linked to previously described mechanisms related to FE in hypothalamus-pituitary-adrenal axis and are expected to help elucidate this complex phenotype. This study contributes to expanding the catalogue of lncRNA, particularly in indicine cattle, and identifies candidates for further studies in animal selection and management.

A Transcriptome Analysis Reveals that Hepatic Glycolysis and Lipid Synthesis Are Negatively Associated with Feed Efficiency in DLY Pigs

Feed efficiency (FE) is an important trait in the porcine industry. Therefore, understanding the molecular mechanisms of FE is vital for the improvement of this trait. In this study, 6 extreme high-FE and 6 low-FE pigs were selected from 225 Duroc × (Landrace × Yorkshire) (DLY) pigs for transcriptomic analysis. RNA-seq analysis was performed to determine differentially expressed genes (DEGs) in the liver tissues of the 12 individuals, and 507 DEGs were identified between high-FE pigs (HE- group) and low-FE pigs (LE- group). A gene ontology (GO) enrichment and pathway enrichment analysis were performed and revealed that glycolytic metabolism and lipid synthesis-related pathways were significantly enriched within DEGs; all of these DEGs were downregulated in the HE- group. Moreover, Weighted gene co-expression analysis (WGCNA) revealed that oxidative phosphorylation, thermogenesis, and energy metabolism-related pathways were negatively related to HE- group, which might result in lower energy consumption in higher efficiency pigs. These results implied that the higher FE in the HE- group may be attributed to a lower glycolytic, energy consumption and lipid synthesizing potential in the liver. Furthermore, our findings suggested that the inhibition of lipid synthesis and glucose metabolic activity in the liver may be strategies for improving the FE of DLY pigs.

Landscape of multi-tissue global gene expression reveals the regulatory signatures of feed efficiency in beef cattle

MOTIVATION

Feed efficiency is an important trait for sustainable beef production that is regulated by the complex biological process, but the mode of action behinds it has not been clearly defined. Here, we aimed to elucidate the regulatory mechanisms of this trait through studying the landscape of the genome-wide gene expression of rumen, liver, muscle and backfat tissues, the key ones involved in the energy metabolism.

RESULTS

The transcriptome of 189 samples across four tissues from 48 beef steers with varied feed efficiency were generated using Illumina HiSeq4000. The analysis of global gene expression profiles of four tissues, functional analysis of tissue-shared and -unique genes, co-expressed network construction of tissue-shared genes, weighted correlations analysis between gene modules and feed efficiency-related traits in each tissue were performed. Among four tissues, the transcriptome of muscle tissue was distinctive from others, while those of rumen and backfat tissues were similar. The associations between co-expressed genes and feed efficiency related traits at single or all tissues level exhibited that the gene expression in the rumen, liver, muscle and backfat were the most correlated with feed conversion ratio, dry matter intake, average daily gain and residual feed intake, respectively. The 19 overlapped genes identified from the strongest module-trait relationships in four tissues are potential generic gene markers for feed efficiency.

AVAILABILITY AND IMPLEMENTATION

The distribution of gene expression data can be accessed at https://www.cattleomics.com/transcriptome.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Potential Biomarkers for Feed Efficiency-Related Traits in Nelore Cattle Identified by Co-expression Network and Integrative Genomics Analyses

Feed efficiency helps to reduce environmental impacts from livestock production, improving beef cattle profitability. We identified potential biomarkers (hub genes) for feed efficiency, by applying co-expression analysis in Longissimus thoracis RNA-Seq data from 180 Nelore steers. Six co-expression modules were associated with six feed efficiency-related traits (p-value ≤ 0.05). Within these modules, 391 hub genes were enriched for pathways as protein synthesis, muscle growth, and immune response. Trait-associated transcription factors (TFs) ELF1, ELK3, ETS1, FLI1, and TCF4, were identified with binding sites in at least one hub gene. Gene expression of CCDC80, FBLN5, SERPINF1, and OGN was associated with multiple feed efficiency-related traits (FDR ≤ 0.05) and were previously related to glucose homeostasis, oxidative stress, fat mass, and osteoblastogenesis, respectively. Potential regulatory elements were identified, integrating the hub genes with previous studies from our research group, such as the putative cis-regulatory elements (eQTLs) inferred as affecting the PCDH18 and SPARCL1 hub genes related to immune system and adipogenesis, respectively. Therefore, our analyses contribute to a better understanding of the biological mechanisms underlying feed efficiency in bovine and the hub genes disclosed can be used as biomarkers for feed efficiency-related traits in Nelore cattle.

Genetic architecture of quantitative traits in beef cattle revealed by genome wide association studies of imputed whole genome sequence variants: I: feed efficiency and component traits

BACKGROUND

Genome wide association studies (GWAS) on residual feed intake (RFI) and its component traits including daily dry matter intake (DMI), average daily gain (ADG), and metabolic body weight (MWT) were conducted in a population of 7573 animals from multiple beef cattle breeds based on 7,853,211 imputed whole genome sequence variants. The GWAS results were used to elucidate genetic architectures of the feed efficiency related traits in beef cattle.

RESULTS

The DNA variant allele substitution effects approximated a bell-shaped distribution for all the traits while the distribution of additive genetic variances explained by single DNA variants followed a scaled inverse chi-squared distribution to a greater extent. With a threshold of P-value < 1.00E-05, 16, 72, 88, and 116 lead DNA variants on multiple chromosomes were significantly associated with RFI, DMI, ADG, and MWT, respectively. In addition, lead DNA variants with potentially large pleiotropic effects on DMI, ADG, and MWT were found on chromosomes 6, 14 and 20. On average, missense, 3'UTR, 5'UTR, and other regulatory region variants exhibited larger allele substitution effects in comparison to other functional classes. Intergenic and intron variants captured smaller proportions of additive genetic variance per DNA variant. Instead 3'UTR and synonymous variants explained a greater amount of genetic variance per DNA variant for all the traits examined while missense, 5'UTR and other regulatory region variants accounted for relatively more additive genetic variance per sequence variant for RFI and ADG, respectively. In total, 25 to 27 enriched cellular and molecular functions were identified with lipid metabolism and carbohydrate metabolism being the most significant for the feed efficiency traits.

CONCLUSIONS

RFI is controlled by many DNA variants with relatively small effects whereas DMI, ADG, and MWT are influenced by a few DNA variants with large effects and many DNA variants with small effects. Nucleotide polymorphisms in regulatory region and synonymous functional classes play a more important role per sequence variant in determining variation of the feed efficiency traits. The genetic architecture as revealed by the GWAS of the imputed 7,853,211 DNA variants will improve our understanding on the genetic control of feed efficiency traits in beef cattle.

Genes associated with body weight gain and feed intake identified by meta-analysis of the mesenteric fat from crossbred beef steers

Mesenteric fat is a visceral fat depot that increases with cattle maturity and can be influenced by diet. There may be a relationship between the accumulation of mesenteric fat and feed efficiency in beef cattle. The purpose of this study was to identify genes that may be differentially expressed in steers with high and low BW gain and feed intake. RNA-Seq was used to evaluate the transcript abundance of genes in the mesenteric fat from a total of 78 steers collected over 5 different cohorts. A meta-analysis was used to identify genes involved with gain, feed intake or the interaction of both phenotypes. The interaction analysis identified 11 genes as differentially expressed. For the main effect of gain, a total of 87 differentially expressed genes (DEG) were identified (P_ADJ<0.05), and 24 were identified in the analysis for feed intake. Genes identified for gain were involved in functions and pathways including lipid metabolism, stress response/protein folding, cell proliferation/growth, axon guidance and inflammation. The genes for feed intake did not cluster into pathways, but some of the DEG for intake had functions related to inflammation, immunity, and/or signal transduction (JCHAIN, RIPK1, LY86, SPP1, LYZ, CD5, CD53, SRPX, and NF2). At P_ADJ<0.1, only 4 genes (OLFML3, LOC100300716, MRPL15, and PUS10) were identified as differentially expressed in two or more cohorts, highlighting the importance of evaluating the transcriptome of more than one group of animals and incorporating a meta-analysis. This meta-analysis has produced many mesenteric fat DEG that may be contributing to gain and feed intake in cattle.

Transcriptome Analysis Identifies Candidate Genes and Pathways Associated With Feed Efficiency in Hu Sheep

In the genetic improvement of livestock and poultry, residual feed intake (RFI) is an important economic trait. However, in sheep, the genetic regulatory mechanisms of RFI are unclear. In the present study, we measured the feed efficiency (FE)-related phenotypes of 137 male Hu lambs, and selected six lambs with very high (n = 3) and very low (n = 3) RFI values and analyzed their liver transcriptomes. A total of 101 differentially expressed genes were identified, of which 40 were upregulated and 61 were downregulated in the low-RFI group compared with that in the high-RFI group. The downregulated genes were mainly concentrated in immune function pathways, while the upregulated genes were mainly involved in energy metabolism pathways. Two differentially expressed genes, ADRA2A (encoding adrenoceptor alpha 2A) and RYR2 (ryanodine receptor 2), were selected as candidate genes for FE and subjected to single nucleotide polymorphism scanning and association analysis. Two synonymous mutations, ADRA2A g.1429 C > A and RYR2 g.1117 A > C, were detected, which were both significantly associated with the feed conversion rate. These findings provide a deeper understanding of the molecular mechanisms regulating FE, and reveal key genes and genetic variants that could be used to genetically improve FE in sheep.

Liver transcriptome profiling of beef steers with divergent growth rate, feed intake, or metabolic body weight phenotypes1

Average daily gain (ADG) and daily dry matter intake (DMI) are key determinants of beef industry profitability. These traits together with metabolic body weight (MWT) are combined as component traits to calculate residual feed intake (RFI), a common measure of feed efficiency in beef cattle. Recently, there have been significant efforts towards molecular genetic characterization of RFI through transcriptomic studies in different breeds and tissues. However, molecular mechanisms of RFI component traits still remain predominately unexplored. Therefore, in the current study, we investigated the hepatic transcriptomic profiles and their associations with ADG, DMI, and MWT in Angus, Charolais, and Kinsella Composite (KC) populations through global RNAseq analyses. In each population and for each trait, 12 steers with extreme phenotypes (n = 6 low and n = 6 high) were analyzed for differential gene expression. These animals were from 20 beef steers of each Angus, Charolais, and KC breed population that were initially selected for a transcriptome study of RFI. At a false discovery rate <0.05 and fold change >1.5, we identified 123, 102, and 78 differentially expressed (DE) genes between high- and low-ADG animals of Angus, Charolais, and KC populations, respectively. For DMI, 108, 180, and 156 DE genes were identified between high- and low-DMI from Angus, Charolais, and KC populations, respectively, while for MWT, 80, 82, and 84 genes were differentially expressed between high- and low-MWT animals in Angus, Charolais, and KC populations, respectively. The identified DE genes were largely breed specific (81.7% for ADG, 82.7% for DMI, and 83% for MWT), but were largely involved in the same biological functions across the breeds. Among the most enriched biological functions included metabolism of major nutrients (lipids, carbohydrates, amino acids, vitamins, and minerals), small molecule biochemistry, cellular movement, cell morphology, and cell-to-cell signaling and interaction. Notably, we identified multiple DE genes that are involved in cholesterol biosynthesis, and immune response pathways for the 3 studied traits. Thus, our findings present potential molecular genetic mechanisms and candidate genes that influence feed intake, growth, and MWT of beef cattle.

Integrating genome-wide co-association and gene expression to identify putative regulators and predictors of feed efficiency in pigs

BACKGROUND

Feed efficiency (FE) has a major impact on the economic sustainability of pig production. We used a systems-based approach that integrates single nucleotide polymorphism (SNP) co-association and gene-expression data to identify candidate genes, biological pathways, and potential predictors of FE in a Duroc pig population.

RESULTS

We applied an association weight matrix (AWM) approach to analyse the results from genome-wide association studies (GWAS) for nine FE associated and production traits using 31K SNPs by defining residual feed intake (RFI) as the target phenotype. The resulting co-association network was formed by 829 SNPs. Additive effects of this SNP panel explained 61% of the phenotypic variance of RFI, and the resulting phenotype prediction accuracy estimated by cross-validation was 0.65 (vs. 0.20 using pedigree-based best linear unbiased prediction and 0.12 using the 31K SNPs). Sixty-eight transcription factor (TF) genes were identified in the co-association network; based on the lossless approach, the putative main regulators were COPS5, GTF2H5, RUNX1, HDAC4, ESR1, USP16, SMARCA2 and GTF2F2. Furthermore, gene expression data of the gluteus medius muscle was explored through differential expression and multivariate analyses. A list of candidate genes showing functional and/or structural associations with FE was elaborated based on results from both AWM and gene expression analyses, and included the aforementioned TF genes and other ones that have key roles in metabolism, e.g. ESRRG, RXRG, PPARGC1A, TCF7L2, LHX4, MAML2, NFATC3, NFKBIZ, TCEA1, CDCA7L, LZTFL1 or CBFB. The most enriched biological pathways in this list were associated with behaviour, immunity, nervous system, and neurotransmitters, including melatonin, glutamate receptor, and gustation pathways. Finally, an expression GWAS allowed identifying 269 SNPs associated with the candidate genes' expression (eSNPs). Addition of these eSNPs to the AWM panel of 829 SNPs did not improve the accuracy of genomic predictions.

CONCLUSIONS

Candidate genes that have a direct or indirect effect on FE-related traits belong to various biological processes that are mainly related to immunity, behaviour, energy metabolism, and the nervous system. The pituitary gland, hypothalamus and thyroid axis, and estrogen signalling play fundamental roles in the regulation of FE in pigs. The 829 selected SNPs explained 61% of the phenotypic variance of RFI, which constitutes a promising perspective for applying genetic selection on FE relying on molecular-based prediction.

Expression of genes related to lipid transport in meat-type ducks divergent for low or high residual feed intake

Objective

This study examined the effects of divergence in residual feed intake (RFI) on expression profiles of key genes related to lipid transport in the liver and duodenal epithelium and their associations with feed efficiency traits in meat-type ducks.

Methods

A total of 1,000 male ducks with similar body weight (1,042.1±87.2 g) were used in this study, and their individual RFI was calculated from 21 to 42 d of age. Finally, the 10 highest RFI (HRFI) and 10 lowest RFI (LRFI) ducks were chosen for examining the expression of key genes related to lipid transport in the liver and duodenal epithelium using quantitative polymerase chain reaction.

Results

In the liver, expression levels of albumin (ALB), CD36 molecule (CD36), fatty acid hydroxylase domain containing 2 (FAXDC2), and choline kinase alpha (CHKA) were significantly higher in LRFI ducks than in HRFI ducks (p<0.01); negative correlations (p<0.05) between expression levels of ALB, CD36, FAXDC2, and CHKA and RFI were detected in the liver. Additionally, ALB expression was strongly positively correlated (p<0.05) with CD36, FAXDC2, CHKA, and apolipoprotein H (APOH) expression in the liver. In duodenal epithelium, we found that mRNA levels of ALB, CD36, FAXDC2, and APOH were significantly higher in LRFI ducks than in HRFI ducks (p<0.01); RFI was strongly negatively correlated (p<0.05) with ALB, FAXDC2, and APOH expression, while ALB expression was strongly positively correlated with APOH expression (p<0.01) in duodenal epithelium. Furthermore, expression levels of both ALB and FAXDC2 genes were significantly associated with feed conversion ratio and RFI in both liver and duodenal epithelium (p<0.05).

Conclusion

Our findings therefore suggest that ALB and FAXDC2 genes might be used as potential gene markers designed to improve feed efficiency in future meat-type duck breeding programs.

The effect of breed and diet type on the global transcriptome of hepatic tissue in beef cattle divergent for feed efficiency

Background

Feed efficiency is an important economic and environmental trait in beef production, which can be measured in terms of residual feed intake (RFI). Cattle selected for low-RFI (feed efficient) have similar production levels but decreased feed intake, while also emitting less methane. RFI is difficult and expensive to measure and is not widely adopted in beef production systems. However, development of DNA-based biomarkers for RFI may facilitate its adoption in genomic-assisted breeding programmes. Cattle have been shown to re-rank in terms of RFI across diets and age, while also RFI varies by breed. Therefore, we used RNA-Seq technology to investigate the hepatic transcriptome of RFI-divergent Charolais (CH) and Holstein-Friesian (HF) steers across three dietary phases to identify genes and biological pathways associated with RFI regardless of diet or breed.

Results

Residual feed intake was measured during a high-concentrate phase, a zero-grazed grass phase and a final high-concentrate phase. In total, 322 and 33 differentially expressed genes (DEGs) were identified across all diets for CH and HF steers, respectively. Three genes, GADD45G, HP and MID1IP1, were differentially expressed in CH when both the high-concentrate zero-grazed grass diet were offered. Two canonical pathways were enriched across all diets for CH steers. These canonical pathways were related to immune function.

Conclusions

The absence of common differentially expressed genes across all dietary phases and breeds in this study supports previous reports of the re-ranking of animals in terms of RFI when offered differing diets over their lifetime. However, we have identified biological processes such as the immune response and lipid metabolism as potentially associated with RFI divergence emphasising the previously reported roles of these biological processes with respect to RFI.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5906-8) contains supplementary material, which is available to authorized users.

A comparative analysis of the transcriptome profiles of liver and muscle tissue in pigs divergent for feed efficiency

BACKGROUND

The improvement of feed efficiency is a key economic goal within the pig production industry. The objective of this study was to examine transcriptomic differences in both the liver and muscle of pigs divergent for feed efficiency, thus improving our understanding of the molecular mechanisms influencing feed efficiency and enabling the identification of candidate biomarkers. Residual feed intake (RFI) was calculated for two populations of pigs from two different farms of origin/genotype. The 6 most efficient (LRFI) and 6 least efficient (HRFI) animals from each population were selected for further analysis of Longissimus Dorsi muscle (n = 22) and liver (n = 23). Transcriptomic data were generated from liver and muscle collected post-slaughter.

RESULTS

The transcriptomic data segregated based on the RFI value of the pig rather than genotype/farm of origin. A total of 6463 genes were identified as being differentially expressed (DE) in muscle, while 964 genes were identified as being DE in liver. Genes that were commonly DE between muscle and liver (n = 526) were used for the multi-tissue analysis. These 526 genes were associated with protein targeting to membrane, extracellular matrix organisation and immune function. In the muscle-only analysis, genes associated with RNA processing, protein synthesis and energy metabolism were down regulated in the LRFI animals while in the liver-only analysis, genes associated with cell signalling and lipid homeostasis were up regulated in the LRFI animals.

CONCLUSIONS

Differences in the transcriptome segregated on pig RFI value rather than the genotype/farm of origin. Multi-tissue analysis identified that genes associated with GO terms protein targeting to membrane, extracellular matrix organisation and a range of terms relating to immune function were over represented in the differentially expressed genes of both liver and muscle.

Comparative Analysis of the Liver Transcriptome among Cattle Breeds Using RNA-seq

Global gene expression in liver transcriptome varies among cattle breeds. The present investigation was aimed to identify the differentially expressed genes (DEGs), metabolic gene networks and metabolic pathways in bovine liver transcriptome of young bulls. In this study, we comparatively analyzed the bovine liver transcriptome of dairy (Polish Holstein Friesian (HF); n = 6), beef (Hereford; n = 6), and dual purpose (Polish-Red; n = 6) cattle breeds. This study identified 895, 338, and 571 significant (p < 0.01) differentially expressed (DE) gene-transcripts represented as 745, 265, and 498 hepatic DE genes through the Polish-Red versus Hereford, Polish-HF versus Hereford, and Polish-HF versus Polish-Red breeds comparisons, respectively. By combining all breeds comparisons, 75 hepatic DE genes (p < 0.01) were identified as commonly shared among all the three breed comparisons; 70, 160, and 38 hepatic DE genes were commonly shared between the following comparisons: (i) Polish-Red versus Hereford and Polish-HF versus Hereford; (ii) Polish-Red versus Hereford and Polish-HF versus Polish-Red; and (iii) Polish-HF versus Hereford and Polish-HF versus Polish-Red, respectively. A total of 440, 82, and 225 hepatic DE genes were uniquely observed for the Polish-Red versus Hereford, Polish-HF versus Hereford, and Polish-Red versus Polish-HF comparisons, respectively. Gene ontology (GO) analysis identified top-ranked enriched GO terms (p < 0.01) including 17, 16, and 31 functional groups and 151, 61, and 140 gene functions that were DE in all three breed liver transcriptome comparisons. Gene network analysis identified several potential metabolic pathways involved in glutamine family amino-acid, triglyceride synthesis, gluconeogenesis, p38MAPK cascade regulation, cholesterol biosynthesis (Polish-Red versus Hereford); IGF-receptor signaling, catecholamine transport, lipoprotein lipase, tyrosine kinase binding receptor (Polish-HF versus Hereford), and PGF-receptor binding, (Polish-HF versus Polish-Red). Validation results showed that the relative expression values were consistent to those obtained by RNA-seq, and significantly correlated between the quantitative reverse transcription PCR (RT-qPCR) and RNA-seq (Pearson’s r  > 0.90). Our results provide new insights on bovine liver gene expressions among dairy versus dual versus beef breeds by identifying the large numbers of DEGs markers submitted to NCBI gene expression omnibus (GEO) accession number GSE114233, which can serve as useful genetic tools to develop the gene assays for trait-associated studies as well as, to effectively implement in genomics selection (GS) cattle breeding programs in Poland.

Systems Biology Reveals NR2F6 and TGFB1 as Key Regulators of Feed Efficiency in Beef Cattle

Systems biology approaches are used as strategy to uncover tissue-specific perturbations and regulatory genes related to complex phenotypes. We applied this approach to study feed efficiency (FE) in beef cattle, an important trait both economically and environmentally. Poly-A selected RNA of five tissues (adrenal gland, hypothalamus, liver, skeletal muscle and pituitary) of eighteen young bulls, selected for high and low FE, were sequenced (Illumina HiSeq 2500, 100 bp, pared-end). From the 17,354 expressed genes considering all tissues, 1,335 were prioritized by five selection categories (differentially expressed, harboring SNPs associated with FE, tissue-specific, secreted in plasma and key regulators) and used for network construction. NR2F6 and TGFB1 were identified and validated by motif discovery as key regulators of hepatic inflammatory response and muscle tissue development, respectively, two biological processes demonstrated to be associated with FE. Moreover, we indicated potential biomarkers of FE, which are related to hormonal control of metabolism and sexual maturity. By using robust methodologies and validation strategies, we confirmed the main biological processes related to FE in Bos indicus and indicated candidate genes as regulators or biomarkers of superior animals.