Background: Fatty acid composition contributes greatly to the quality and nutritional value of meat. However, the molecular regulatory mechanisms underlying fatty acid accumulation in poultry have not yet been cleared. The aims of this study were to characterize the dynamics of fatty acid accumulation in duck breast muscle and investigate its correlations with gene expression. Results: Here, we analyzed the fatty acid profile and transcriptome of breast muscle derived from Pekin ducks and mallards at the ages of 2 weeks, 4 weeks, 6 weeks and 8 weeks. Twenty fatty acids were detected in duck breast muscle, with palmitic acid (C16:0, 16.6%~ 21.1%), stearic acid (C18:0, 9.8%~ 17.7%), oleic acid (C18:1n-9, 15.7%~ 33.8%), linoleic acid (C18:2n-6, 10.8%~ 18.9%) and arachidonic acid (C20:4n-6, 11.7%~ 28.9%) as the major fatty acids. Our results showed that fatty acid composition was similar between the two breeds before 6 weeks, but the compositions diverged greatly after this point, mainly due to the stronger capacity for C16:0 and C18:1n-9 deposition in Pekin ducks. By comparing the multistage transcriptomes of Pekin ducks and mallards, we identified 2025 differentially expressed genes (DEGs). Cluster analysis of these DEGs revealed that the genes involved in oxidative phosphorylation, fatty acid degradation and the PPAR signaling pathway were upregulated in mallard at 8 weeks. Moreover, correlation analysis of the DEGs and fatty acid composition traits suggested that the DEGs involved in lipogenesis, lipolysis and fatty acid β-oxidation may interact to influence the deposition of fatty acids in duck breast muscle. Conclusions: We reported the temporal progression of fatty acid accumulation and the dynamics of the transcriptome in breast muscle of Pekin ducks and mallards. Our results provide insights into the transcriptome regulation of fatty acid accumulation in duck breast muscle, and will facilitate improvements of fatty acid composition in duck breeding.

Background: As a major economic trait in poultry, egg production efficiency attracts widespread interest in breeding and production. However, limited information is available about the underlying genetic architecture of egg production traits in ducks. In this paper, we analyzed six egg production-related traits in 352 F ₂ ducks derived from reciprocal crosses between mallard and Pekin ducks. Results: Feed conversation ratio (FCR) was positively correlated with feed intake but negatively correlated with egg-related traits, including egg weight and egg production, both phenotypically and genetically. Estimates of pedigree-based heritability were higher than 0.2 for all traits investigated, except hip-width. Based on whole-genome sequencing data, we conducted genome-wide association studies to identify genomic regions associated with these traits. In total, 11 genomic regions were associated with FCR. No genomic regions were identified as significantly associated with hip-width, total feed intake, average daily feed intake, and total egg production. Analysis of selective sweeps between mallard and Pekin ducks confirmed three of these genomic regions on chromosomes 13, 3 and 6. Within these three regions, variants in candidate genes that were in linkage disequilibrium with the GWAS leader single nucleotide polymorphisms (SNPs) (Chr13:2,196,728, P = 7.05 × 10 ^–14 ; Chr3:76,991,524, P = 1.06 × 10 ^–12 ; Chr6:20,356,803, P = 1.14 × 10 ^–10 ) were detected. Thus, we identified 31 potential candidate genes associated with FCR, among which the strongest candidates are those that are highly expressed in tissues involved in reproduction and nervous system functions of ducks: CNTN4 , CRBR , GPR63 , KLHL32 , FHL5 , TRNT1 , MANEA , NDUFAF4 , and SCD . Conclusions: For the first time, we report the identification of genomic regions that are associated with FCR in ducks and our results illustrate the genomic changes that occurred during their domestication and are involved in egg production efficiency.

Chickens play a crucial role as the primary global source of eggs and poultry, and the quality of rooster semen significantly impacts poultry reproductive efficiency. Therefore, it is imperative to comprehend the regulatory mechanisms underlying sperm development. In this study, we established transcriptome profiles of lncRNAs, miRNAs, and mRNAs in 3 testis tissues and 3 epididymis tissues from “Jing Hong No.1” roosters at 24, 35, and 64 weeks of age. Using the data, we conducted whole transcriptome analysis and constructed a ceRNA network. We detected 10 differentially expressed mRNAs (DEmRNAs), 33 differentially expressed lncRNAs (DElncRNAs), and 10 differentially expressed miRNAs (DEmiRNAs) in the testis, as well as 149 DEmRNAs, 12 DElncRNAs, and 10 DEmiRNAs in the epididymis. These genes were found to be involved in cell differentiation and development, as well as various signaling pathways such as GnRH, MAPK, TGF-β, mTOR, VEGF, and calcium ion pathways. Subsequently, we constructed two competing endogenous RNA (ceRNA) networks comprising DEmRNAs, DElncRNAs, and DEmiRNAs. Furthermore, we identified four crucial lncRNA-mRNA-miRNA interactions that govern specific biological processes in the chicken reproductive system: MSTRG.2423.1-gga-miR-1563-PPP3CA and MSTRG.10064.2-gga-miR-32-5p-GPR12 regulating sperm motility in the testis; MSTRG.152556.1-gga-miR-9-3p-GREM1/THYN1 governing immunomodulation in the epididymis; and MSTRG.124708.1-gga-miR-375-NDUFB9/YBX1 controlling epididymal sperm maturation and motility. Whole transcriptome sequencing of chicken testis and epididymis screened several key genes and ceRNA regulatory networks, which may be involved in the regulation of epididymal immunity, spermatogenesis and sperm viability through the pathways of MAPK, TGF-β, mTOR, and calcium ion. These findings contribute to our comprehensive understanding of the intricate molecular processes underlying rooster spermatogenesis, maturation and motility.

Background: Carcass traits are crucial for broiler ducks, but carcass traits can only be measured postmortem. Genomic selection (GS) is an effective approach in animal breeding to improve selection and reduce costs. However, the performance of genomic prediction in duck carcass traits remains largely unknown. Results: In this study, we estimated the genetic parameters, performed GS using different models and marker densities, and compared the estimation performance between GS and conventional BLUP on 35 carcass traits in an F₂ population of ducks. Most of the cut weight traits and intestine length traits were estimated to be high and moderate heritabilities, respectively, while the heritabilities of percentage slaughter traits were dynamic. The reliability of genome prediction using GBLUP increased by an average of 0.06 compared to the conventional BLUP method. The Permutation studies revealed that 50K markers had achieved ideal prediction reliability, while 3K markers still achieved 90.7% predictive capability would further reduce the cost for duck carcass traits. The genomic relationship matrix normalized by our true variance method instead of the widely used $$\sum {2p}_{i}(1-{p}_{i})$$∑ 2 p i ( 1 - p i ) could achieve an increase in prediction reliability in most traits. We detected most of the bayesian models had a better performance, especially for BayesN. Compared to GBLUP, BayesN can further improve the predictive reliability with an average of 0.06 for duck carcass traits. Conclusion: This study demonstrates genomic selection for duck carcass traits is promising. The genomic prediction can be further improved by modifying the genomic relationship matrix using our proposed true variance method and several Bayesian models. Permutation study provides a theoretical basis for the fact that low-density arrays can be used to reduce genotype costs in duck genome selection.