Abstract
Crisp grass carp (CGC), the most representative improved varieties of grass carp (GC), features higher muscle hardness after feeding faba bean (Vicia faba L.) for 90-120 days. DNA methylation, a most widely studied epigenetic modification, plays an essential role in muscle development. Previous studies have identified numerous differentially expressed genes (DEGs) between CGC and GC. However, it remains unknown if the expression levels of these DEGs are influenced by DNA methylation. In the present study, we performed a comprehensive analysis of DNA methylation profiles between CGC and GC, and identified important candidate genes related to muscle development coupled with the transcriptome sequencing data. A total of 9,318 differentially methylated genes (DMGs) corresponding to 155,760 differentially methylated regions (DMRs) were identified between the two groups under the CG context in promoter regions. Combined with the transcriptome sequencing data, 14 key genes related to muscle development were identified, eight of which (gsk3b, wnt8a, wnt11, axin2, stat1, stat2, jak2, hsp90) were involved in muscle fiber hyperplasia, six of which (tgf-β1, col1a1, col1a2, col1a3, col4a1, col18a1) were associated with collagen synthesis in crisp grass carp. The difference of methylation levels in the key genes might lead to the expression difference, further resulting in the increase of muscle hardness in crisp grass carp. Overall, this study can help further understand how faba bean modulates muscle development by the epigenetic modifications, providing novel insights into the texture quality improvement in other aquaculture fish species by nutritional programming.