Feed efficiency in dairy sheep: An insight from the milk transcriptome

Aroa Suárez-Vega; Pilar Frutos; Beatriz Gutiérrez-Gil; Cristina Esteban-Blanco; Pablo G Toral; Juan-José Arranz; Gonzalo Hervás

doi:10.3389/fvets.2023.1122953

Feed efficiency in dairy sheep: An insight from the milk transcriptome

Front Vet Sci. 2023 Apr 3:10:1122953. doi: 10.3389/fvets.2023.1122953. eCollection 2023.

Authors

Aroa Suárez-Vega¹, Pilar Frutos², Beatriz Gutiérrez-Gil¹, Cristina Esteban-Blanco¹, Pablo G Toral², Juan-José Arranz¹, Gonzalo Hervás²

Affiliations

¹ Departamento de Producción Animal, Facultad de Veterinaria, Universidad de León, León, Spain.
² Instituto de Ganadería de Montaña (CSIC-Universidad de León), Grulleros, León, Spain.

Abstract

Introduction: As higher feed efficiency in dairy ruminants means a higher capability to transform feed nutrients into milk and milk components, differences in feed efficiency are expected to be partly linked to changes in the physiology of the mammary glands. Therefore, this study aimed to determine the biological functions and key regulatory genes associated with feed efficiency in dairy sheep using the milk somatic cell transcriptome.

Material and methods: RNA-Seq data from high (H-FE, n = 8) and low (L-FE, n = 8) feed efficiency ewes were compared through differential expression analysis (DEA) and sparse Partial Least Square-Discriminant analysis (sPLS-DA).

Results: In the DEA, 79 genes were identified as differentially expressed between both conditions, while the sPLS-DA identified 261 predictive genes [variable importance in projection (VIP) > 2] that discriminated H-FE and L-FE sheep.

Discussion: The DEA between sheep with divergent feed efficiency allowed the identification of genes associated with the immune system and stress in L-FE animals. In addition, the sPLS-DA approach revealed the importance of genes involved in cell division (e.g., KIF4A and PRC1) and cellular lipid metabolic process (e.g., LPL, SCD, GPAM, and ACOX3) for the H-FE sheep in the lactating mammary gland transcriptome. A set of discriminant genes, commonly identified by the two statistical approaches, was also detected, including some involved in cell proliferation (e.g., SESN2, KIF20A, or TOP2A) or encoding heat-shock proteins (HSPB1). These results provide novel insights into the biological basis of feed efficiency in dairy sheep, highlighting the informative potential of the mammary gland transcriptome as a target tissue and revealing the usefulness of combining univariate and multivariate analysis approaches to elucidate the molecular mechanisms controlling complex traits.

Keywords: RNA-Seq; dairy sheep; feed efficiency; mammary gland; sPLS-DA.

Grants and funding

This research was funded by the Junta de Castilla y León (JCyL, Spain; project CSI276P18) and the Spanish Research State Agency (Agencia Estatal de Investigación, Spain; project PID2020-113441RB-I00 and MCIN/AEI/10.13039/501100011033). Cofunding by the European Regional Development Fund (ERDF/FEDER) is also acknowledged. Research activities for the University of Leon were funded by the SMARTER project through the H2020 research and innovation program of the European Commission (Grant Agreement No. 772787).