Milk transcriptome biomarker identification to enhance feed efficiency and reduce nutritional costs in dairy ewes

In recent years, rising prices for high-quality protein-based feeds have significantly increased nutrition costs. Consequently, investigating strategies to reduce these expenses and improve feed efficiency (FE) have become increasingly important for the dairy sheep industry. This research investigates the impact of nutritional protein restriction (NPR) during prepuberty and FE on the milk transcriptome of dairy Assaf ewes (sampled during the first lactation). To this end, we first compared transcriptomic differences between NPR and control ewes. Subsequently, we evaluated gene expression differences between ewes with divergent FE, using feed conversion ratio (FCR), residual feed intake (RFI), and consensus classifications of high- and low-FE animals for both indices. Lastly, we assess milk gene expression as a predictor of FE phenotype using random forest. No effect was found for the prepubertal NPR on milk performance or FE. Moreover, at the milk transcriptome level, only one gene, HBB, was differentially expressed between the NPR (n = 14) and the control group (n = 14). Further, the transcriptomic analysis between divergent FE sheep revealed 114 differentially expressed genes (DEGs) for RFI index (high-FE_RFI = 10 vs low-FE_RFI = 10), 244 for FCR (high-FE_FCR = 10 vs low-FE_FCR = 10), and 1 016 DEGs between divergent consensus ewes for both indices (high-FE_consensus = 8 vs low-FE_consensus = 8). These results underscore the critical role of selected FE indices for RNA-Seq analyses, revealing that consensus divergent animals for both indices maximise differences in transcriptomic responses. Genes overexpressed in high-FE_consensus ewes were associated with milk production and mammary gland development, while low-FE_consensus genes were linked to higher metabolic expenditure for tissue organisation and repair. The best prediction accuracy for FE phenotype using random forest was obtained for a set of 44 genes consistently differentially expressed across lactations, with Spearman correlations of 0.37 and 0.22 for FCR and RFI, respectively. These findings provide insights into potential sustainability strategies for dairy sheep, highlighting the utility of transcriptomic markers as FE proxies.