摘要

为了研究五指山猪背最长肌与半腱肌的遗传差异,鉴定差异甲基化区域(DMRs)和差异甲基化基因(DMGs),采用全基因组重亚硫酸盐测序技术(WGBS)对五指山猪背最长肌和半腱肌进行DNA甲基化水平检测和DMRs分析,对DMGs进行基因本体(GO)功能、京都基因与基因组百科全书(KEGG)通路富集分析,并与已有的转录组数据进行整合分析,探讨基因组甲基化对基因表达差异的影响。结果显示:五指山猪背最长肌和半腱肌全基因组范围内胞嘧啶甲基化主要发生在CG序列上;在CG序列环境下,启动子和第一外显子区域的甲基化水平低于其他基因区域,该环境下五指山猪背最长肌和半腱肌间存在348 613个DMRs和17 484个DMGs, DMGs主要富集到61个GO条目和289个KEGG通路;在五指山猪背最长肌和半腱肌mRNA差异表达前10位的基因中,锌指结构家族成员4(ZIC4)、锌指结构家族成员1(ZIC1)、成对同源异型结构域蛋白转录因子1(PITX1)、碱性-螺旋-环-螺旋转录因子1(SIM1)、小清蛋白(PVALB)、SEC14类脂质结合因子5(SEC14L5)、丝氨酸蛋白酶抑制剂B家族成员11(SERPINB11)、肌球蛋白结合蛋白H (MYBPH)这8个基因存在DMRs, 5个基因(ZIC4、PITX1、SIM1、SEC14L5、MYBPH)存在启动子区甲基化差异。本研究为阐明五指山猪骨骼肌生长发育机制提供表观遗传学依据。