Comparative Transcriptome Analysis Reveals the Key Genes Involved in Lipid Deposition in Pekin Ducks (Anas platyrhynchos domesticus)

There are differences in lipid deposition in fatty-type (FT) and lean-type (LT) ducks. Fatty ducks have a higher rate of sebum and abdominal fat, lower meat yield and hepatic lipid contents than LT ducks. However, the underlying changes in gene expression profiles regarding the lipid deposition between FT and LT ducks have not yet been clarified. To identify the differentially expressed genes in the liver, sebum, and abdominal fat between both ducks, we identified the gene expression profiles in the liver, sebum, and abdominal fat derived from FT and LT ducks by comparing the multistage transcriptomes. Our results showed that there were 622, 1536, and 224 differentially expressed genes (DEGs) in the liver, sebum, and abdominal fat between the FT and LT ducks, respectively. KEGG enrichment showed that the DEGs related to lipid metabolism were enriched in the biosynthesis of unsaturated fatty acid, glycerolipid and fatty acid metabolism in the liver; and were enriched in the fatty acid metabolism, fatty acid biosynthesis, glycerolipid metabolism, linoleic acid metabolism, and the PPAR signaling pathway in the sebum. There was no pathway related to a lipid metabolism enriched in abdominal fat. A gene functional analysis showed that the DEGs involved in adipogenesis were found to be upregulated. In contrast, those involved in lipolysis were downregulated in the liver and serum of FT ducks. The DEGs showed that ATP-binding cassette sub-family G member 8 (ABCG8), fatty acid synthase (FASN), and phospholipid transfer protein (PLTP) were highly expressed in the liver of FT ducks, and acyl-CoA synthetase long-chain family member3 (ACSL3), ACSL5, ACSL6, 1-acyl-sn-glycerol-3-phosphate acyltransferase alpha (AGPAT1), AGPAT9, ELOVL fatty acid elongase 6 (ELVOL6), fatty acid desaturase 1 (FADS1), FADS2, monoacylglycerol O-acyltransferase 1 (MOGAT1), serine/threonine kinase 17a (STK17A), and serine/threonine kinase 39 (STK39) were highly expressed in the sebum of FT ducks. A weighted correlation network analysis (WGCNA) of the DEGs showed ABCG8, FADS2, ACSL5, and ELOVL6 positively correlated with hepatic fatty acid synthesis, and AGPAT1, STK17A, STK32A, FADS1, and ACSL3 positively correlated with lipid deposition in the sebum. In summary, ABCG8 might be the key gene for the reduced hepatic lipid deposition in FT Pekin ducks, and FADS2, ACSL5, ELOVL6, AGPAT1, STK17A, STK32A, FADS1, and ACSL3 were the key genes for lipid deposition in the sebum of FT Pekin ducks. Our results provide new insights into the transcriptome regulation in lipid deposition of Pekin ducks and will be helpful for duck breeding.