Comparative Transcriptome Analysis of Key Genes and Pathways Activated in Response to Fat Deposition in Two Sheep Breeds With Distinct Tail Phenotype

Fat tail in sheep presents a valuable energy reserve that has historically facilitated adaptation to harsh environments. However, in modern intensive and semi-intensive sheep industry systems, breeds with leaner tails are more desirable. In the present study, RNA sequencing (RNA-Seq) was applied to determine the transcriptome profiles of tail fat tissues in two Chinese sheep breeds, fat-rumped Altay sheep and thin-tailed Xinjiang fine wool (XFW) sheep, with extreme fat tail phenotype difference. Then the differentially expressed genes (DEGs) and their sequence variations were further analyzed. In total, 21,527 genes were detected, among which 3,965 displayed significant expression variations in tail fat tissues of the two sheep breeds (P < 0.05), including 707 upregulated and 3,258 downregulated genes. Gene Ontology (GO) analysis disclosed that 198 DEGs were related to fat metabolism. In Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis, the majority of DEGs were significantly enriched in “adipocytokine signaling,” “PPAR signaling,” and “metabolic pathways” (P < 0.05); moreover, some genes were involved in multiple pathways. Among the 198 DEGs, 22 genes were markedly up- or downregulated in tail fat tissue of Altay sheep, indicating that these genes might be closely related to the fat tail trait of this breed. A total of 41,724 and 42,193 SNPs were detected in the transcriptomic data of tail fat tissues obtained from Altay and XFW sheep, respectively. The distribution of seven SNPs in the coding regions of the 22 candidate genes was further investigated in populations of three sheep breeds with distinct tail phenotypes. In particular, the g.18167532T/C (Oar_v3.1) mutation of the ATP-binding cassette transporter A1 (ABCA1) gene and g.57036072G/T (Oar_v3.1) mutation of the solute carrier family 27 member 2 (SLC27A2) gene showed significantly different distributions and were closely associated with tail phenotype (P < 0.05). The present study provides transcriptomic evidence explaining the differences in fat- and thin-tailed sheep breeds and reveals numerous DEGs and SNPs associated with tail phenotype. Our data provide a valuable theoretical basis for selection of lean-tailed sheep breeds.