DNA methylation is a pivotal epigenetic regulatory mechanism in the development of skeletal muscles. Nonetheless, the regulators responsible for DNA methylation in the development of embryonic duck skeletal muscles remain unknown. In the present study, whole genome bisulfite sequencing (WGBS) and transcriptome sequencing were conducted on the skeletal muscles of embryonic day 21 (E21) and day 28 (E28) ducks. The DNA methylation pattern was found to fall mainly within the cytosine-guanine (CG) context, with high methylation levels in the intron, exon, and promoter regions. Overall, 7902 differentially methylated regions (DMRs) were identified, which corresponded to 3174 differentially methylated genes (DMGs). By using integrative analysis of both WGBS with transcriptomics, we identified 1072 genes that are DMGs that are negatively associated with differentially expressed genes (DEGs). The gene ontology (GO) analysis revealed significant enrichment in phosphorylation, kinase activity, phosphotransferase activity, alcohol-based receptors, and binding to cytoskeletal proteins. The Kyoto Encyclopedia of Genes and Genomes (KEGGs) analysis showed significant enrichment in MAPK signaling, Wnt signaling, apelin signaling, insulin signaling, and FoxO signaling. The screening of enriched genes showed that hyper-methylation inhibited the expression of Idh3a, Got1, Bcl2, Mylk2, Klf2, Erbin, and Klhl38, and hypo-methylation stimulated the expression of Col22a1, Dnmt3b, Fn1, E2f1, Rprm, and Wfikkn1. Further predictions showed that the CpG islands in the promoters of Klhl38, Klf2, Erbin, Mylk2, and Got1 may play a crucial role in regulating the development of skeletal muscles. This study provides new insights into the epigenetic regulation of the development of duck skeletal muscles.
Keywords: DNA methylation; RNA-seq; WGBS; embryonic ducks; skeletal muscle development.