MyoG和MEF2a基因多态性聚合效应对鸭屠宰性状的影响

赵忠海，李　辉，易恒洁，杨胜林，彭邦星，卜小雁

（贵州大学动物科学学院/高原山地动物遗传育种与繁殖教育部重点实验室，贵阳 550025）

MyoG和MEF2a基因多态性聚合效应对鸭屠宰性状的影响

赵忠海，李辉，易恒洁，杨胜林，彭邦星，卜小雁

（贵州大学动物科学学院/高原山地动物遗传育种与繁殖教育部重点实验室，贵阳 550025）

【目的】探讨MyoG和MEF2a基因聚合效应对鸭屠宰性状的影响，为进一步确定与鸭屠宰性状相关的分子遗传标记提供研究基础，为鸭屠宰性状的多基因聚合育种提供依据。【方法】试验以240只三穗鸭为研究素材，扩增MyoG和MEF2a基因并进行PCR产物直接测序以检测两基因所有外显子的单核苷酸突变（SNPs）位点。运用SPSS 18.0软件中的GLM统计模型对MyoG和MEF2a基因的SNPs所对应的不同基因型与三穗鸭屠宰性状进行关联分析，根据单基因关联分析结果，将对屠宰性状存在显著影响的MyoG和MEF2a基因的多态位点利用软件PHASE 2.0构建聚合基因型，再进行聚合基因型与屠宰性状的关联分析。【结果】在试验群体中一共发现8个SNPs，其中在MyoG基因中有6个SNPs被找到，MEF2a基因中找到2个SNPs位点，在所有突变位点中，其中MyoG基因的g.2977G＞C位点发生的G/C突变使密码子由GAG变为GAC，所编码的氨基酸由谷氨酸变成天冬氨酸；而MEF2a基因中的两个多态位点，g.47915G＞A位点发生的G/A突变使密码子由GAA变为AAA，编码的氨基酸由谷氨酸变成赖氨酸，g.47918G＞A位点的G/A突变引起的密码子由GAT变成AAT，所编码的氨基酸由天冬氨酸变成天冬酰胺。剩下的5个突变位点均属于同义突变，并未引起编码氨基酸的改变。此外，进行χ2适合性检验，除了MyoG基因的g.1131C＞T位点和MEF2a基因的g.47915G＞A、g.47918G＞A位点未处于Hardy-Weinberg平衡状态（P＜0.05）外，其他的突变位点均处于平衡状态。单基因关联分析结果表明，MyoG基因g.1131C＞T和g.2204G＞A突变分别对胸肌率、体重和全净膛重有着显著影响，其所对应的纯合子基因型CC、GG型为优势基因型。MEF2a基因g.47915G＞A/g.47918G＞A位点影响全净膛率，GA基因型个体属于优势基因型个体。通过挑选出与屠宰性状（胸肌率、体重、全净膛重和全净膛率）有关联的MyoG基因g.1131C＞T/g.2204G＞A位点与MEF2a基因g.47915G＞A/g.47918G＞A进行聚合效应分析，结果显示，聚合后的8种聚合基因型个体的全净膛率，在各基因型间无显著差异（P＞0.05），TTGAGA基因型的平均值最高，其次为CCGGGA基因型；其他3个指标各基因型间差异达到了显著，其中体重和全净膛重存在正相关，都是CCGAGA基因型的平均值最高，CTGGGA基因型的平均值次之；CCGGGG基因型的胸肌率平均值最高，其次是CCGGGA基因型。结果显示，单个基因的平均值最高的基因型分别为CC、GG和GA，在两基因聚合后在4个指标中CCGGGA基因型都不是最优的组合，说明两个基因间存在互作效应。【结论】两个基因间存在互作效应，所以用单个基因分子标记进行选育可能会顾此失彼，不能收到良好的效果，但是本研究的聚合优势基因型个体偏少，有待于进一步扩大样本进行验证分析，进行更多基因的聚合效应分析。

MyoG基因；MEF2a基因；屠宰性状；聚合效应

0　引言

【研究意义】畜禽屠宰性能由于与经济效益密切相关，在畜禽育种工作中一直备受重视，但屠宰性状属于数量性状，受多基因控制，易受环境影响，因此采用传统的选育方法进行改良进展缓慢。分子标记由于其结果不受环境影响，不存在等位基因显隐性表型关系，并且可进行早期选择，被认为是缩短育种周期，加快遗传进展的有效手段。但是基因之间往往存在互作效应，因此单基因分子标记往往会出现顾此失彼的现象，而多基因聚合的分子育种技术将有望弥补这方面的缺陷。因此，本研究选择在畜禽研究中被证明与畜禽屠宰性状相关的两个候选基因-肌细胞生成素（MyoG）基因和肌细胞增强因子2a（MEF2a）基因，进行两基因聚合效应对屠宰性状的影响研究，为实现畜禽的多基因聚合育种提供依据。【前人研究进展】由于MyoG基因和MEF2a基因在肌细胞分化和肌肉形成中发挥重要作用，而肌肉的发育与畜禽屠宰性状息息相关，所以这两个基因在畜禽研究中一直被作为屠宰性状的候选基因。MyoG基因是生肌家族因子中最重要的一员，在肌细胞的分化过程中起着重要的作用，能促进肌细胞的增殖和分化［1］。而MEF2基因家族（MEF2a、MEF2b、MEF2c和MEF2d）在骨骼发育、肌肉形成、肝脏纤维化和神经系统发育等多种生理过程中发挥作用，同时它也参与一些疾病的发生，如阿尔茨海默氏症和帕金森氏病等［2］。在骨骼肌发育的过程中，敲除microRNA可导致试验动物在胚胎期致死或骨骼肌发育不全、肌纤维形态异常、肌原细胞凋亡和成肌细胞死亡增加［3-4］。LIU等［5］通过构建小干扰RNA质粒载体获得敲除牛MyoG基因骨骼肌细胞的试验结果显示该基因的表达显著低于其他试验阳性基因。NEVILLE等［6］研究表明，MyoG通过控制成肌细胞的融合和肌纤维的形成而对肌肉的分化起关键作用，是唯一能在所有骨骼肌细胞中表达的基因，是骨骼肌分化所必须的基因。许多骨骼肌和心肌基因的调控区都存在MEF2的结合位点，能够与大多数肌肉特异基因的启动子或增强子直接结合，在所有的肌肉细胞类型中可作为肌源性基因表达的主要调节物［7-9］。LIU［10］等对鸭的MEF2a基因mRNA表达进行探究，发现在平滑肌中的表达量高于心肌和骨骼肌。JUSZCZUKKUBIAK等［11］对波兰荷斯坦奶牛MEF2a基因启动子区多态性及该基因mRNA表达水平进行研究，结果显示，启动子区域突变的不同基因型个体的背最长肌MEF2a基因mRNA水平显著高于其他基因型个体，因此认为MEF2a的核苷酸序列突变可能作为牛的生长性状指标分子标记。【本研究切入点】MyoG和MEF2a基因在人、小鼠及其他家畜方面已有大量研究，家禽方面近几年的研究也在逐步增多，但水禽方面的研究相对较少，迄今为止还未见到关于鸭的MyoG和MEF2a基因聚合效应分析的报道。虽然植物的多基因聚合育种已取得了初步成功，但动物的聚合育种还未见到成功的报道。本研究在前人研究的基础上进一步探讨MyoG和MEF2a基因之间的关系，探讨聚合基因型与鸭屠宰性状的关联性，为畜禽多基因聚合育种的研究提供基础。【拟解决的关键问题】以贵州省优良地方品种——三穗鸭为研究对象，首次将MyoG和MEF2a基因的聚合基因型与三穗鸭屠宰性状进行关联分析，进行多基因聚合效应的探讨，为家禽屠宰性状分子标记选育提供依据。

1　材料与方法

试验于2014年7月采集样品，2014年10月至2015年6月在贵州大学动物科学学院高原山地动物遗传育种与繁殖教育部重点实验室进行相关操作。

1.1试验材料

17周龄三穗鸭240只（公母各半），相同饲养管理水平，翅静脉采血，肝素钠抗凝，置-70℃冰箱保存。屠宰测定按照中华人民共和国农业部制订的《家禽生产性能名词术语和度量统计方法》（2004）进行［12］。

1.2基因组DNA提取与引物设计

采用血液/细胞/组织基因组DNA提取试剂盒（北京天根）提取基因组DNA。根据NCBI上发布的家鸭MyoG基因编码区全序列（GenBank: NW_004676592.1）、家鸭MEF2a基因编码区全序列（GenBank: NW_004676438.1），应用软件Primer Premier 5.0及Oligo 6.0软件设计引物。引物由上海英潍捷基生物技术有限公司合成，引物详细信息见表1。

1.3PCR扩增

PCR扩增采用30 μL反应体系，包括2×Taq PCR Master Mix（北京天根）12 μL，上、下游引物各3.0 μL（浓度均为10 μmol·L-1），模板DNA（50 ng·μL-1）4 μL，ddH2O 8 μL。反应条件如下：94℃预变性4 min，94℃变性30—35 s，退火30—35 s（退火温度如表1所示），72℃延伸30—40 s，35个循环。扩增完成后，72℃再延伸10 min后，4℃保存。

1.4数据处理分析

采用DNAstar进行序列对比拼接，PHASE 2.0软件进行两基因聚合基因型的判定，Excel 2010进行数据整理，SPSS18.0软件中一般线性模型（GLM）进行关联分析，结果以均值±标准差体现。统计模型如下：

单基因效应关联分析数学模型：

Y=μ+Gi+S+Gi×S+e（i=1，2）

聚合基因型效应关联分析数学模型：Y=μ+G+e其中：Y为性状测定值，μ代表群体均值，Gi为MyoG或MEF2a基因型效应，S为样品性别，Gi×S为基因型与性别互作效应，G为MyoG和MEF2a基因的聚合基因型效应，e为随机残差。

2　结果

2.1PCR 产物检测

扩增目的片段，1.2%琼脂糖电泳检测PCR扩增目的产物片段，结果见图1、2。扩增产物大小与设计引物预期扩增片段长度一致。

2.2PCR扩增产物测序及序列分析

根据目的片段扩增结果，将所有PCR产物送往上海英潍捷基生物技术有限公司进行直接测序。测序结果经序列拼接、对比，并与GenBank上公布的基因参考序列矫正，发现MyoG基因共有6个SNPs（图3、4），分别位于外显子1的g.1131C＞T位点，外显子2的g.2186G＞A位点和g.2204G＞A位点，外显子3的g.2920G＞A位点、g.2962C＞T位点和g.2977G＞C位点，其中，g.2977G＞C位点发生的G/C突变使密码子由GAG变为GAC，所编码的氨基酸由谷氨酸变成天冬氨酸。MEF2a基因共有2个SNPs（图5），分别位于外显子11的g.47915G＞A位点和g.47918G＞A位点，g.47915G＞A位点发生的G/A突变使密码子由GAA变为AAA，所编码的氨基酸由谷氨酸变成赖氨酸，g.47918G＞A位点的G/A突变引起的密码子由GAT变成AAT，编码的氨基酸由天冬氨酸变成天冬酰胺。

表1　MyoG、MEF2a基因引物信息Table 1　The primer sequences and their information for PCR amplification of the MyoG and MEF2a genes in Sansui duck

图1　MyoG基因PCR扩增Fig. 1　The amplification result of MyoG from PCR

图2　MEF2a基因PCR扩增Fig. 2　The amplification result of MEF2a from PCR

图3　MyoG基因1131、2920和2977 bp处的碱基变异Fig. 3　The sequencing map of the MyoG gene at positions 1131， 2920 and 2977

图4　MyoG基因2186、2204 和2962 bp处的碱基变异Fig.4　The sequencing map of the MyoG gene at position 2186， 2204 and 2962

2.3MyoG基因和MEF2a基因遗传学分析

测序结果经序列比对发现，在MyoG基因外显子扩增片段中存在6个SNPs位点，将6个SNPs位点所对应纯合子和杂合子的基因型分别定义为CC、CT、TT；GG、GA；GG、GA；GG、GA、AA；CC、 CT和GG、GC、CC型。对几种基因型分析，可知相对应的基因型中TT、GG、GG、GA、CC、GG为优势基因型，等位基因T、G、G、G、C、G则为优势等位基因。在MEF2a基因外显子11中有2个SNPs位点，这两个多态位点属于完全连锁平衡位点，将它们所对应纯合子和杂合子的基因型分别定义为GG、GA型。对基因型进行分析，可知相对应的基因型中GA为优势基因型，等位基因G则为优势等位基因。

图5　MEF2a基因47915和47918 bp处的碱基变异Fig. 5　The sequencing map of MEF2a gene at position 47915 and 47918

统计各基因型的个体数，分别计算MyoG基因和MEF2a基因中的SNPs位点的基因型频率和基因频率，计算结果如表2所示。χ2适合性检验结果表明，在MyoG基因中除g.1131C＞T外，该基因座其他多态位点均处于Hardy-Weinberg平衡状态（P＞0.05）。对MEF2a基因中的2个SNPs进行χ2适合性检验，结果表明，该基因这两个多态位点均未处于Hardy-Weinberg平衡状态（P＜0.05），这可能是由于在育种过程中的人为选择或遗传漂变造成的。

2.4MyoG基因与MEF2a基因多态性与屠宰性状的相关性分析

MyoG基因多态性与三穗鸭屠宰性状的显著性关系分析如表3所示。可以看出，公鸭的胸肌率在MyoG基因g.1131C＞T位点的不同基因型间差异显著，CC型显著高于CT和TT型，TT型显著高于CT型。MyoG基因g.2204G＞A的不同基因型对体重和全净膛重有显著影响，GG型均显著高于GA型。其余的性状在各个SNPs的不同基因型间和性别间差异不显著。MEF2a基因的2个突变位点g.47915G＞A/g.47918G＞A的不同基因型对全净膛率有显著影响（P＜0.05），且GA型显著高于GG型，其余的性状在两个SNPs的不同基因型、性别间无差异显著性，分析统计结果如表4所示。

表2　MyoG、MEF2a基因多态位点基因频率及基因型频率Table 2　Gene frequency and genotype frequency of SNPs in MyoG and MEF2a gene

2.5聚合基因型与屠宰性状间相关性分析

根据单个SNP位点分析结果，只选取与屠宰性状有显著影响的3个SNPs位点，采用PHASE 2.0软件，进行两基因聚合基因型的判定，选取的性状也是在单基因关联分析时存在显著差异的指标。由于聚合后试验群体个体数较少缘故，故不考虑性别因素。聚合后共得到8种聚合基因型（表5）。从体重指标看，试验群体中CTGGGG型体重除与TTGAGA型个体间差异不显著，与其他基因型个体间均呈现显著差异，且CTGGGA、TTGGGG与CTGGGG基因型间呈差异极显著（P＜0.01），TTGAGA与TTGGGG、CTGGGA间存在差异显著（P＜0.05）；从全净膛重指标看，几种聚合基因型中，除了CTGGGG基因型和TTGAGA基因型间差异不显著，与其他几种基因型间均差异显著；CCGAGA基因型的体重和全净膛重平均值都最高，但个体数量少，有必要扩大样本做进一步验证，其次为CTGGGA基因型，两个指标一致，这是由于两个指标呈正相关引起的。此外，基因型CCGGGA的胸肌率与TTGAGA基因型间表现出显著差异（P＜0.05），与CTGGGG、CTGGGA、TTGGGG、TTGGGA间均呈差异极显著关系（P＜0.01），CCGGGG聚合基因型的胸肌率平均值最高，其次是CCGGGA基因型；在全净膛率这个指标，TTGAGA基因型的平均值最高，其次为CCGGGA基因型，但各基因型间差异都没达到显著水平。

3　讨论

3.1MyoG基因多态性对肌肉的影响

在不影响肉质风味的前提下，动物能快速生长、肌肉沉积率高一直以来是畜禽育种的目标之一，研究证实，MyoG基因能在所有骨骼肌中表达，其表达程度直接影响畜禽的产肉性能［13］。WEINTRAUB等［14-15］研究表明，敲除MyoG基因的小鼠出生后会发生骨骼肌的发育严重缺陷。缺失MyoG基因的小鼠要比正常小鼠个体小。唐莹等［16］利用PCR-SSCP技术检测京海黄鸡MyoG基因外显子多态性，结果表明，MyoG基因多态与京海黄鸡生长性状相关。ZHANG等［17］再一次证实了MyoG基因多态性对京海黄鸡的生长性状有影响。XUE等［19］报道，肌细胞生成素基因单核苷酸多态性对猪的初生重和背标厚度有显著影响。王健等［20］对太湖鹅MyoG基因研究报道，在外显子1中第108位点处发现碱基C/T突变，推断该基因多态性对太湖鹅早期增重有显著影响。本研究中，共发现三穗鸭MyoG基因6个SNPs，其中g.2977G＞C位点发生的G/C突变使密码子由GAG变为GAC；所定义的几种基因型中，MyoG基因g.1131C＞T的不同基因型胸肌率有显著影响，CC型显著高于CT和TT型，TT型显著高于CT型，公鸭的胸肌率差异显著。MyoG基因g.2204G＞A的不同基因型对体重和全净膛重有显著影响，GG型均显著高于GA型。其余的性状在各个SNPs的不同基因型间和性别间差异不显著。

率eat 2.08肉）1±瘦Leanmpercentage（%24.4率uscle tage1.51肌）5±胸Breastmpercen（%12.3率tage1.35肌uscle）腿Legmpercen（%12.06±率膛）2.96净（%全iscerated Evpercentage.41±75率膛2.56净tage（%）半Semi-evisceratedpercen.76±84率tage）2.25Traits宰屠Dressingpercen（%2±93.1状性重uscle肌14胸Breastmweight（g）122±htertraits 重ht肌uscle）9腿Legmweig（g8±11Ggeneandslaug重膛）64净（g全iscerated Evweight986±重ht8 yo膛）净（g08±6Semi-eviscerated半weig11析重）±85分体ghter（g联屠Slauweight1218关的状重）性体（g宰Weight1308±76屠与er点别）位性Gend态♀（24多型因因type8）（6基基G GenoCC3　MyoheassociationresultsbetweenthegenotypeofMT Table3　T点Ps表位SNg.1131C＞2.924±22.51.27b7±11.32.1211.16±5.20.14±735.944±81.92.880±94.220115±202±111911±1103733±111±12213021383±136）♂（442.750±23.21.40b2±11.71.8911.48±4.57.94±735.104±82.92.6693.82±187±1117114±02±110021624±111±11512721356±121T（68）2.230±22.60.897±10.91.4011.62±3.57.40±763.57.71±853.814±93.68113±109±11834±9101261±111±11812681357±154）♀（20（24）CT2.411±19.32a7.89±1.31.5911.42±4.03.32±781.433±88.42.659±96.2±2185193±12377±9101316±112±10413241375±113♂（4）2.41.05±22a1.49.46±101.258±11.53.28.72±763.38.16±863.578±94.013108±90±12941±8100270±111±10812771360±138T（24）3.072±22.61.853±11.31.8011.29±4.59.80±725.00.08±822.359±91.926114±233±1187999±326±911±8812631373±82）♀（56）（100TT2.247±22.51.15c4±10.61.7911.93±3.15.56±742.96.45±833.448±92.315109±223±12626±810749±911±11112721377±110）♂（442.68.59±22c1.59.02±111.787±11.54.04.58±734.20.68±822.826±92.122112±237±11611±810336±911±9712671374±93T（100）3.761±24.23.520±11.50.240±12.70.94.86±730.44.45±831.240±93.523109±92±1291957±83181±±11512121295±106♀（8）（12）GAAg.2186G＞1.992±23.31.038±12.31.9110.94±4.71.85±743.691±.8833.079±93.319146±179±122179±1114220±113±10814711575±192♂（4）2.71.91±232.54.79±111.031±12.10.87.19±740.37.57±830.806±93.427121±84±124331±1105060±111±17112981388±178T（12）2.788±22.61.675±11.21.5911.43±4.08.05±744.30.28±832.554±92.421114±185±11208±810134±911±8912591363±0972）♀（11）（228GG2.335±22.21.283±10.61.7511.61±4.15.25±744.578±83.13.407±93.718109±219±112724±1104447±111±14812931379±1486）♂（112.556±22.41.503±10.91.6611.52±4.08.15±744.403±83.23.0693.11±191±1119117±06±110162041±111±12312761371±124T（228）1.275±24.10.649±12.21.9111.86±3.35.00±763.81.95±842.856±93.624120±14±11514973±6187±110±19512001285±247♀（8）（16）GAAg.2204G＞0.874±22.50.882±10.50.0112.02±0.23.73±710.04.28±811.298±94.1±1791114±1089861±98976±±12911311200±120♂（8）1.285±23.31.190±11.41.1011.94±3.13.87±733.052±83.11.832±93.924106±99±107a11917±2731±110±14111661243±166AT（16）2.869±22.61.789±11.11.5811.49±4.00.89±734.21.17±832.503±92.421113±186±11907±710834±811±8312601364±862）♀（11）（224GG2.386±22.21.330±10.71.7811.56±4.16.45±744.634±83.33.452±93.718111±210±122241±1104066±111±14513111399±1442）♂（112.62.48±221.58.95±101.673±11.54.05.17±744.38.26±833.068±93.019112±198±1104b24±1101750±111±12012851381±119BT（224）

率eat 2.37肉）2±瘦Leanmpercentage（%24.2率uscle tage2.30肌）5±胸Breastmpercen（%11.7率tage0.32肌uscle）腿Legmpercen（%12.46±率膛）1.81净（%全iscerated Evpercentage.33±75率膛2.36净tage（%）半Semi-evisceratedpercen.45±85率tage）1.80Traits宰屠Dressingpercen（%0±94.3状性重uscle肌20胸Breastmweight（g）115±重ht肌uscle）5腿Legmweig（g2±12重膛）61净（g全iscerated Evweight982±重膛ht）3净（g半Semi-evisceratedweig13±611重体ghter）±81屠Slauweight（g1229重）体Weight（g1303±67别er性）edtable3Gend♀（16型因typetinu基Geno（56）onAA3　CA表Ps续点位SNg.2920G＞2.221±21.71.292±10.21.3411.49±3.05.94±743.121±84.44.307±93.319104±207±110514±1101542±111±12012611351±113）♂（402.46.43±221.70.65±101.217±11.72.69.05±752.87.71±843.713±93.619107±178±11404±9100133±111±10812521337±102T（56）2.865±22.51.786±10.81.5011.69±4.42.83±724.67.84±812.291±91.519108±166±1187993±0017±111±8612481365±105）♀（56）（100GA2.396±22.01.145±10.52.0711.50±5.41.95±735.995±82.62.549±94.517110±219±116445±1108868±111±18813361411±175）♂（442.624±22.31.512±10.71.7411.61±4.80.32±735.190±82.22.8292.86±189±1018117±27±110164439±111±14412871385±139T（100）2.928±22.51.568±11.51.7810.99±3.68.01±753.57.26±842.548±93.022120±223±11226±810853±811±9912741369±96）♀（48（84）GG2.279±23.11.338±11.31.8211.81±3.51.91±733.833±82.53.186±93.121117±221±121127±1102948±111±12812941392±161）♂（362.62.84±221.43.49±111.804±11.33.57.54±743.69.52±832.752±93.121119±227±11327±9100551±111±11012831379±125T（84）2.849±22.61.741±11.21.6211.47±3.88.81±734.16.13±832.552±92.520112±185±11302±810229±911±9012561358±962）♀（1132）（2CCTg.2962C＞2.308±22.21.309±10.61.7211.59±4.08.27±744.500±83.23.345±93.719110±209±112829±1104553±111±14912991385±1490）♂（122.568±22.41.544±10.91.6611.53±3.95.05±744.307±83.13.0393.16±191±1119117±09±110162241±111±12512781372±126T（232）0 0 0 0 0 0000000♀（0）（8）CT2.175±24.12.272±12.00.0912.13±4.95.13±774.59.55±851.863±92.430127±67±12547±510462±711±12212581363±159♂（8）2.17.15±242.27.02±120.093±12.14.95.13±774.53.55±851.863±92.430127±67±12547±510462±711±12212581363±159T（8）3.600±23.02.340±11.21.5211.79±4.07.57±752.78.76±841.681±91.431111±226±11177±809095±710±7511821293±81）♀（168）（2CCCg.2977G＞1.527±21.70.753±11.61.8710.14±4.34.62±733.975±82.10.961±92.217117±262±100904±1103021±111±18212621370±213）♂（122.777±22.41.728±11.31.7611.09±3.96.74±743.33.64±831.385±91.724113±230±1186988±106±911±12612161326±142T（28）2.582±22.01.263±11.31.6210.69±5.54.36±735.89.56±822.812±93.220117±220±11228±9100157±111±6813061401±67）♀（482）（7GC3.230±23.21.275±11.42.7611.75±3.71.35±744.166±83.23.184±94.319119±281±124242±1106568±111±14213191399±154）♂（242.782±22.41.237±11.32.0511.04±4.91.69±735.269±82.72.9093.59±198±1124114±0733±1102161±111±9513101400±99T（72）2.808±23.32.043±11.21.2912.14±1.94.17±742.46.50±832.352±92.219111±110±1273993±619±811±9012351341±111）♀（56）（140GG2.139±22.01.243±10.31.3111.75±4.33.34±744.824±83.33.621±93.819106±171±123229±1104754±111±15312981384±148）♂（842.461±22.6数1.64体otsignificant（P＞个9±本asn10.6样为字1.29数中11.91±号括，值3.53均平.27±的总74母公体3.99群验0±试ithoutanyidentificationtablethedifferencew83.4示），w0.053.23。“T”表ples93.17±.05）＞0（Pifference（P＜berofsam著19显8±不10异差eanthenum示表15识120±标何任parenthesesm无12±15），1015＜0.0bersbetween（Pum著26显40±1平水nwithdifferentsuperscriptlettersindicatesignificantdeanthen110.05示ecolumracketsm±134表母1273字同mbersinb不nu有1366±134，esitesandtraits，thesam具列同比对T（140）同paredtothesam间型因alearetotalaverageand基不的状aleandfem性一eansm同点同Differentgenotypescom位一）.“T”m0.05

率肉meat）3.46瘦Leanpercentage（%）1.8110uscle 6±胸肌率Breastmpercentage（%.54±tageuscle率）2.77Legm肌腿percen（%.12±11膛率Eviscerated）6.11percentage净5±70.4（%全膛）6.22率净（%Semi-eviscerate dpercentage 1±79.4）2.13半率屠宰TraitsDressingpercentage（%.67±91状性重uscle 21.6ht20肌）胸weig（g102±ghtertraits Breastm重uscle ht）29slau肌腿Legmweig（g107±eand2agen重膛）87净weight（g967±EF全Eviscerated heassociationresultsbetweenthegenotypeofM重膛ht7析净weig（g）90±9分半10联Semi-eviscerated关的重7状体ghter（g）性屠Slauweight60±812宰屠与点重ht位（g）±80态Weig1374多因基别type体）2a性4　MEFGenoTable4　T♀（20型pe表因基GenotyGG1.793±21.31.28.21±100.98.12±111.9972.87±1.885±82.50.74.78±9223102±11110±155996±7128±1116966±112±1811364）♂（202.600±21.51.49.37±101.96.12±114.47b71.66±4.6480.98±1.61.22±9220102±20109±119982±3309±1112763±112±1321369T（40）2.67.01±231.73.41±111.2911.59±3.065±74.73.257±84.02.56.68±9220116±15117±±801013939±811255±912±10113550）♀（10GA2.387±22.41.30.78±101.84.69±114.3574.55±4.8783.33±3.635±93.9182±1122121±1036±1244358±1114705±113±14613900）♂（102.52.74±221.55.10±111.5711.64±3.72a5±74.64.120±83.73.17.31±9319114±18119±±10410241849±1112480±112±1261372T（200）nificant（P＞asnotsig数体个本样为字数中号括均iden值，tificationtablethedifferencew平的总母公验0.05），withoutany体群试示表T”。“.05）＞0（P著显不异差ples表记fsam标ercaselettersindicatesignificantdifference（P＜无bero，.05）um＜0（P异ithdifferentsuperscriptloweanthen著显thesesm差示nwaren表母berinp不etraits，thesam字ecolumum写小文英同有具列同，alearetotalaverageandn状性paredtothesam同相比escomaleandfem对间eansm型enotyp不Differen因tg基同）.“T”m0.05A7918G＞/g.4type＞AGeno5G型g.4791因A、Ps/基G＞SNg.2204T、31C＞g.11析分otypeandslaughtertraits联关状性宰regatedgen屠与型因基nalysisofagg合5　聚表Table5　A状性TraitsGAGA2）cdTT（11225±67CGAGG）Bab TT（116±22A1369GGGG8）a TT（24A1439±4GAGG4）a CT（2±48A1442GGGG）d CT（82C±81168GAa GA1ACC（4）60±314GAGG4）Bab CC（41379±35AGGCbc GG）3BCC（4±21280ht（g）Weig重体cBC609±92a±19A1022Bab1033±39ABab±42A105473Cd831±a3A75±310Bab1034±31ABbc28A960±t（g）eighEvisceratedw重膛净全Bb0.629±.411c0.20BC4±10.4bcBC0.41.59±10c0.44CD3±10.1d.76D±09.18Aa0.3612.74±0.32Aa5±13.0Aa0.48.65±13e（%）ercentagusclepBreastm率肌胸Aa2.3375.95±0.75Aab4±74.6Ab1.5271.78±1.65Aab7±72.9A2.85.15±710.89Aab73.63±1.22Aab3±75.0Aab1.02.00±75Evisceratedrate（%）率膛净全osignificantdifference.eann样0.01），and数samelettersm体个本为字数中号1level（P＜括，著显不异差示表母字nificantdifferenceat0.0同相有eansig具1），＜0.0（Pitallettersm著显异差平水0.010.05），differentcap在示ers（P＜表母字.05lev写大同不5），＜0.0（P著nificantdifferenceat0ples异eansig显差平berofsam水0.05母ifferentsmalllettersm在示eanthenum表字，dthesesm同erow写小同Inthesam不标ersinparen行Numb

3.2MEF2a基因多态性

研究表明，多种生理过程与MEF2基因家族相关，如骨骼发育、肌肉形成、肝脏纤维化和神经系统发育等［2］。MEF2a基因属于MEF2基因家族一员，它在其中也扮演着至关重要的角色，NAYA等［21］报道，小鼠敲除MEF2a基因，大多数在第一周就会表现出明显心室扩张，肌纤维碎片等现象而死亡，能存活的个体到了成年因缺乏心脏线粒体也容易发生猝死。LIEB等［22］认为，缺少MEF2基因对冠心病（CAD）不存在影响，然而在2012年，LIU等［23］对不同人群的冠心病（CAD）病例进行分析，结果虽未能证实冠心病的产生与MEF2a基因的多态性有直接关系，但在1 008个研究病例中，有5个病例在MEF2a基因外显子11区域出现了21个碱基缺失，因此推测MEF2a基因的变异可能是造成CAD原因之一。2009年，周艳等［24］报道，MEF2a基因的3个SNPs位点与鸡的部分屠宰性状有显著相关，可推断MEF2a基因可作为影响鸡屠宰性状的候选基因。本研究显示，在三穗鸭MEF2a基因外显子11中发现的2个SNPs位点，分别定义为GG、GA两种基因型，分析表明，GA为优势基因型，等位基因G则为优势等位基因，与屠宰性状进行关联分析表明，MEF2a基因g.47915G＞A/g.47918G＞A的不同基因型对全净膛率有显著影响，GA型显著高于GG型，其余的性状在两个SNPs的不同基因型、性别间差异不显著。

3.3MyoG、MEF2a基因与三穗鸭屠宰性状的关联

研究表明，MyoG与MEF2a基因具有协同作用。常国斌等［25］研究报道，对鸡肌内脂肪影响最明显的优势单个基因，经两基因或三基因聚合后，并未表现出相应的优势情况；而在扩大样本数进一步研究发现，最佳聚合基因型与相应的单个基因的有利基因型结果一致，且整体效应要高于单个基因型，表现出一定程度的累加效应［26］。其他研究也表明，单个SNP所对应基因型与性状间关联分析结果在不同的品种间通常是不同的，没有多个位点联合分析准确［27-29］。本研究以三穗鸭为研究对象，在单基因关联分析时，3个位点的平均值最高的基因型分别是CC、GG和GA型，但是在聚合效应分析时，并不是CCGGGA基因型的平均值最高，这说明两个基因间存在互作效应，并没有表现出累加效应，这与常国斌等［25-26］的研究结果不一致。在基因间存在交互作用时，如果利用单基因进行分子标记辅助选择往往会顾此失彼，不能收到良好的育种效果，因此本研究的结果表明对影响同一性状的不同基因进行聚合效应分析是十分必要的。本研究结果所显示的平均值最高的聚合基因型个体数量都偏少，只有4个，因此有必要进一步扩大样本进行验证，其次本研究只进行了两个基因的聚合效应分析，证明两个基因间存在互作，然而影响屠宰性状的基因较多，所以有必要进行更多基因的聚合效应分析，同时在其它鸭品种中是否存在与本研究一致的聚合效应也需要做进一步验证，以期得到更真实有效的分子标记，从而为鸭屠宰性状的多基因聚合育种提供依据。再者，通过对不同模式生物的深入研究，对于人类医学进行肌肉相关疾病的研究、治疗也能提供一定的帮助或具有参考意义，与此同时，在进行人类肌肉相关疾病的研究时，也应该考虑从多基因聚合的角度去研究，也许能达到事半功倍的效果。

4　结论

本试验首次将鸭的MyoG基因和MEF2a基因进行聚合效应分析，结果单个基因的有利基因型聚合后并不是最佳基因型，说明基因间存在互作效应，但是最佳基因型在本群体中个体数偏少，有必要扩大样本做进一步研究和进行更多个基因聚合效应分析，才能为家禽的多基因聚合育种提供有效的分子标记。

References

［1］ 宋兴超， 魏海军， 杨镒峰， 陈秀敏， 薛海龙， 岳志刚.不同物种肌细胞生成素基因序列结构与功能特性的生物信息学分析.畜牧与兽医，2012， 44（7）:64-68.

SONG X C， WEI H J， YANG Y F， CHEN X M， XUE H L， YUE Z G. Bioinformatics analysis of different species of muscle cell gene structure and function characteristics. Animal Husbandry & Veterinary Medicine， 2012， 44（7）:64-68.（in Chinese）

［2］ DIETRICH J B. The MEF2 family and the brain: from molecules to memory. Springer， 2013， 352（2）:179-190.

［3］ O'ROURKE J R， GEORGES S A， SEAY H R， TAPSCOTT S J，MCMANUS M T， GOLDHAMER D J， SWANSON M S， HARFE B D. Essential role for Dicer during skeletal muscle development. Development Biology， 2007， 311（2）: 359-368.

［4］ HORAK M， NOVAK J， BIENERTOVA-VASKU J. Muscle-specific microRNAs in skeletal muscle development. Development Biology，2015， 410（1）: 1-13.

［5］ LIU C C， ZHAO D D， TONG H L， YE F， YANG Y， LI S F， JIA M Y，YAN Y Q. Impact of bovine skeletal muscle satellite cell differentiation by small interfering RNA targeting myogenin gene.Journal of Northeast Agricultural University （English Edition）， 2013，20（2）:32-37.

［6］ NEVILLE C M， SCHMIDT M， SCHMIDT J. Response of myogenic determination factors to cessation and resumption of electrical activity in skeletal muscle:a possible role for myogenin in denervation supersensitivity. Cellular and Molecular Neurobiology， 1992，12（6）:511-527.

［7］ GONG H J， XIE J， ZHANG N， YAO L， ZHANG Y. MEF2A binding to the Glut4 promoter occurs via an AMPKa2-dependent mechanism. Medicine&Science in Sports&Exercise， 2011， 43（8）:1441-1450.

［8］ VARGAS M A， TIRNAUER J S， GLIDDEN N， KAPILOFF M S，DODGE-KAFKA K L. Myocyte enhancer factor 2（MEF2） tethering to muscle selective A-kinase anchoring protein （mAKAP） is necessary for myogenic differentiation. Cellular Signalling， 2012， 24（8）: 1496-1503.

［9］ SNYDER C M， RICE A L， ESTRELLA N L. MEF2A regulates the Gtl2-Dio3 microRNA mega-cluster to moudulate WNT signaling in skeletal muscle regeneration. Development， 2013， 140（1）:31-42.

［10］ LIU H H， WANG J W， SI J M， JIA J， LI L， HAN C C， HUANG K L，HE H， XU F. Molecuar cloning and in silico analysis of the duck（Anas platyrhynchos） MEF2A gene Cdna and its expression profile in muscle tissues during fetal development. Genetics and Molecular Biology， 2012， 35（1）:182-190.

［11］ JUSZCZUKKUBIAK E， STARZYŃSKI R R， WICIŃSKA K，FLISIKOWSKI K. Promoter variant-dependent mRNA expression of the MEF2A in longissimus dorsi muscle in cattle. DNA and Cell Biology， 2012， 31（6）:1131-1135.

［12］ 中华人民共和国农业部. 家禽生产性能名词术语和度量统计方法（NY/T 823-2004）. 2004-8-25.

Ministry of Agriculture of the People's Republic of China. Performance ferms and measurement for poultry （NY/T 823-2004）. 2004- 8-25.（in Chinese）

［13］ LIU Y Y， LI F N， KONG X F， TAN B， LI Y H， DUAN Y H，BLACHIER F， CHIEN-AN A. HU C A A， YIN Y L. Signaling pathways related to protein synthesis and amino acid concentration in pig skeletal muscles depend on the dietary protein level， genotype and developmental stages. PLoS ONE， 2015， 10（9）:1-21.

［14］ WEINTRAUB H， DAVIS R， TAPSCONTT S， THAVER M，KRAUSE M， BENEZRA R， BLACKWELL T K， TURNER D， RUPP R， HOLLENBERQ S. The MyoD gene family: nodal point during specification of the muscle cell lineage. Science， 1991， 251（4995）: 761-766.

［15］ HSATY P， BRADLEY A， MORRIS J H， EDMONDSON D G， VENUTI J M， OLSON E N， KLEIN W H. Muscle deficiency and neonatal death in mice with a targeted mutation in the myogenin gene. Nature， 1993， 364（6437）:501-506.

［16］ 唐莹， 王金玉， 张跟喜， 施会强， 张涛. MyoG基因外显子1多态性与京海黄鸡生长性状的相关性分析.中国畜牧杂志， 2013， 49（23）: 5-8.

TANG Y， WANG J Y， ZHANG G X， SHI H Q， ZHANG T. Relationship between polymorphisms of Exon 1 of MyoG gene and growth traits in Jinhai Yellow Chicken. Chinese Journal of Animal Science， 2013， 49（23）， 5-8.（in Chinese）

［17］ ZHANG G X， TANG Y， ZHANG T， WANG J Y， WANG Y J. Expression profiles and association analysis with growth traits of the MyoG and Myf5 genes in the Jinghai yellow chicken. Molecular Biology Reports， 2014， 41（11）: 7331-7338.

［18］ KNAPP J R， DAVIE J K， MYER A， MEADOWS E， OLSON E N，KLEIN W H. Loss of myogenin in postnatal life leads to normal skeletal muscle but reduced body size. Development， 2006， 133（4）: 601-610.

［19］ XUE H L， ZHOU Z X. Effects of the MyoG gene on the partial growth traits in pigs. Acta Genetica Sinica， 2006， 33 （11）:992-997.

［20］ 王健， 董飚， 侯庆永， 殷洁鑫. 太湖鹅MyoG基因多态性与体重的相关性分析.浙江农业学报， 2015， 27（1）:28-31.

WANG J， DONG B， HOU Q Y， YIN J X. Association analysis between of MyoG gene and body weight in Taihu goose. Acta Agriculture Zhejiangensis， 2015， 27（1）， 28-31.（in Chinese）

［21］ NAYA F J， BLACK B L， WU H， BASSEL-DUBY R， RICHARDSON J A， HILL J A， OLSON E N. Mitochondrial deficiency and cardiac sudden death in mice lacking the MEF2A transcription factor. Nature Medicine， 2002， 8（11）:1303-1309.

［22］ LIEB W， MAYER B， KÖNIG I R， BORWITZKY I， GÖTZ A， KAIN S， HENGSTENBERG C， LINSEL-NITSCHKE P， FISCHER M，DÖRING A， WICHMANN H.-E， MEITINGER T， KREUTZ R，ZIEGLER A， SCHUNKERT H， ERDMANN J. Lack of association between the MEF2A gene and myocardial infarction. Circulation，2008， 117: 185-191.

［23］ LIU Y， NIU W Q， WU Z J， SU X X， CHEN Q J， LU L， JIN W. Variants in Exon 11 of MEF2A gene and coronary artery disease: evidence from a case-control study， systematic review， and Meta-analysis. PLoS ONE， 2012， 7（2）:1-10.

［24］ 周艳， 刘益平， 蒋小松， 杜华锐， 朱庆. 优质鸡MEF2A基因的SNPs检测及其与屠体性状的相关研究.畜牧兽医学报. 2009，40（8）:1164-1170.

ZHOU Y， LIU Y P， JIANG X S， DU H R， ZHU Q. Study onassociation of single nucleotide polymorphism of MEF2A gene with carcassr taits in chicken. Acta Veterinaria et Zootechnica Sinica， 2009，40（8）， 1164-1170.（in Chinese）

［25］ 常国斌， 周琼， 雷黎立， 张学余， 王克华， 陈蓉， 栾德琴， 陈国宏.鸡肌内脂肪性状的多基因聚合效应分析.中国家禽， 2009， 31（19）: 25-28.

CHANG G B， ZHOU Q， LEI L L， ZHANG X Y， WANG K H， CHEN R， LUAN D Q， CHEN G H. Genetic analysis of polygene pyramiding in intramuscular fat traits in chicken. China Poultry， 2009， 32（19）: 25-28.（in Chinese）

［26］ 常国斌， 刘向萍， 陈蓉， 栾德琴， 王克华， 张颖， 马腾， 周伟， 戴爱琴， 陈国宏. 鸡肌内脂肪性状候选基因的聚合效应及初步验证.中国农业科学， 2011， 44（20）: 4284-4294.

CHANG G B， LIU X P， CHEN R， LUAN D Q， WANG K H， ZHANG Y， MA T， ZHOU W， DAI A Q， CHEN G H. Pyramiding effect and preliminary verification of candidate genes for intramuscular fat traits in chickens. Scientia Agricultura Sinica， 2011， 44（20）:4284-4294.（in Chinese）

［27］ ZHANG W H， COLLINS A， MORTON N E. Does haplotype diversity predict power for association mapping of disease susceptibility？ Human Genetics， 2004， 115: 157-164.

［28］ CLARK A G. The role of haplotypes in candidate gene studies. Genetic Epidemiology， 2004， 27（4）:321-333.

［29］ BADER J S. The relative power of SNPs and haplotype as genetic markers for association tests. Pharmacogenomics， 2001， 2（1）: 11-24.

（责任编辑林鉴非）

The Effect of MyoG and MEF2a Gene Pyramiding on Slaughter Traits of Ducks

ZHAO Zhong-hai， LI Hui， YI Heng-jie， YANG Sheng-lin， PENG Bang-xing， BU Xiao-yan
（College of Animal Science， Guizhou University/Key Laboratory of Animal Genetics， Breeding and Reproduction in the Plateau Mountainous Region， Ministry of Education， Guizhou University， Guiyang 550025）

【Objective】 The aim of the present study was to explore the polymerization effects of MyoG and MEF2a genes on duck slaughter traits in order to provide a research foundation for further determining the molecular genetic markers related to duck growth traits， also provide a basis of polygene pyramiding breeding of slaughter traits of ducks.【Method】 A total of 240 individuals of Sansui ducks were selected as experimental material in the study， MyoG gene and MEF2a gene were amplificated and had PRC direct sequencing to detect the single nucleotide mutation （SNPs） of all exons of two genes. Base mutation （SNPs） was detected by direct sequencing of the PCR products. GLM statistical model of SPSS 18.0 software was used to analyze the association with different genotypes corresponding to the SNPs MyoG gene and MEF2a gene with Sansui duck slaughter traits. Based on the single gene association analysis results， the polymorphic sites of MyoG and MEF2a genes withsignificant influence on slaughter traits were employed to build polymerization genotype by using software PHASE 2.0. 【Result】The result showed that eight SNPs were found in MyoG gene and MEF2a gene， and six SNPs were found in MyoG gene and two SNPs were found in MEF2a gene. In all mutations， the G/C mutation in the g.2977G＞C SNP of MyoG gene resulted in the change of codon from GAG to GAC， and the coding amino acid changed from Glu to Asp； While 2 polymorphic site in the MEF2a gene， the G/A mutation at the g.47915G＞A SNP and the G/A mutation at the g.47918G＞A SNP led to codon change from GAA to AAA and GAT to AAT， and the coding amino acid from Glu/Lys and Asp/Asn. The other five SNPs belonged to synonymous mutations， which did not cause the variation of encoding amino acids. Besides， the SNPs fit with Hardy-Weinberg equilibrium except that g.1131C＞T of MyoG and g.47915G＞A，g.47918G＞A of MEF2a gene which were tested by χ2. The results of correlation analysis between polymorphism sites and slaughter traits showed that the SNP of g.1131C＞T and the SNP of g.2204G＞A in MyoG gene had significant influence over the breast muscle percentage， the body weight and eviscerated weight， and the correspondings to homezygote genotype CC and GG were dominant genotypes. The SNP of g.47915G＞A and g.47918G＞A in MEF2a gene affected the eviscerated weight， and the GA genotype individuals belong to dominant genotype individuals. The g.1131C＞T and g.2204G＞A in MyoG gene and g.47915G＞A and g.47918G＞A in MEF2a gene， which relating to slaughter traits （body weight， eviscerated weight， breast muscle rate and eviscerated rate） were selected and the multiple gene polymerizations （interaction） were analyzed， the results showed that after polymerization， the eviscerated rate of eight kinds of aggregated genotype individuals were not significantly different among different genotypes， the mean value of TTGAGA genotype was the highest， followed by CCGGGA genotype. The differences of other three indexes among different genotypes reached a significant level， and weight and eviscerated weight were positively correlated， and the CCGAGA genotype was the highest， followed by CTGGGA genotype； The average rate of chest muscle of CCGGGG genotype was the highest， followed by CCGGGA genotype. The result indicated that the highest main value of genotype of single gene was CC， GG and GA. After two genes combined， CCGGGA genotype in the four indicators was not the optional combination， which showed that there exist interactive effect between MyoG gene and MEF2a gene.【Conclusion】The results revealed that one single molecular marker breeding maybe not good and cannot obtain good result from the interaction of two genes. However， regnant aggregated genotype individuals was not more than enough， more samples should be selected to investigate the aggregated effect of more genes in further study， and to obtain effective molecular markers for poultry breeding.

MyoG gene； MEF2a gene； slaughter traits； polymerization effect

2015-08-27；接受日期：2016-07-12

教育部科学技术研究重点项目（211168）、贵州省科技厅农业重大专项［黔科合重大专项字（2012）6004号］、《三穗鸭国家标准》制定与养殖技术规程编制横向（H120183）、贵州省科技合作计划联合基金项目［黔科合LH字（2015）7677号］

联系方式：赵忠海，E-mail：andyzhzhao@163.com。通信作者李辉，Fax：0851-88298003；E-mail：ellenlihui@sina.cn