Abstract
OBJECTIVE: The objective of this study was to identify candidate genes that play important roles in skeletal muscle development in ducks. METHODS: In this study, we investigated the transcriptional sequencing of embryonic pectoral muscles from two specialized lines: Liancheng white ducks (female) and Cherry valley ducks (male) hybrid Line A (LCA) and Line C (LCC) ducks. In addition, prediction of target genes for the differentially expressed mRNAs was conducted and the enriched gene ontology (GO) terms and Kyoto encyclopedia of genes and genomes signaling pathways were further analyzed. Finally, a protein-to-protein interaction network was analyzed by using the target genes to gain insights into their potential functional association. RESULTS: A total of 1,428 differentially expressed genes (DEGs) with 762 being up-regulated genes and 666 being down-regulated genes in pectoral muscle of LCA and LCC ducks identified by RNA-seq (p<0.05). Meanwhile, 23 GO terms in the down-regulated genes and 75 GO terms in up-regulated genes were significantly enriched (p<0.05). Furthermore, the top 5 most enriched pathways were ECM-receptor interaction, fatty acid degradation, pyruvate degradation, PPAR signaling pathway, and glycolysis/gluconeogenesis. Finally, the candidate genes including integrin b3 (Itgb3), pyruvate kinase M1/2 (Pkm), insulinlike growth factor 1 (Igf1), glucose-6-phosphate isomerase (Gpi), GABA type A receptorassociated protein-like 1 (Gabarapl1), and thyroid hormone receptor beta (Thrb) showed the most expression difference, and then were selected to verification by quantitative realtime polymerase chain reaction (qRT-PCR). The result of qRT-PCR was consistent with that of transcriptome sequencing. CONCLUSION: This study provided information of molecular mechanisms underlying the developmental differences in skeletal muscles between specialized duck lines.