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

Abstract

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.