牛PRKAA2基因SNP检测及与其生长和肉质性状的关联分析

田万强,梅楚刚,付常振,辛亚平,张　松,王国庆,昝林森,3

(1 杨凌职业技术学院 动物工程分院，陕西 杨凌 712100；2 西北农林科技大学 动物科技学院，陕西 杨凌 712100；3 国家肉牛改良中心，陕西 杨凌 712100)



牛PRKAA2基因SNP检测及与其生长和肉质性状的关联分析

田万强1,2,梅楚刚2,付常振2,辛亚平2,张松2,王国庆2,昝林森2,3

(1 杨凌职业技术学院 动物工程分院，陕西 杨凌 712100；2 西北农林科技大学 动物科技学院，陕西 杨凌 712100；3 国家肉牛改良中心，陕西 杨凌 712100)

[摘要]【目的】 研究秦川牛、南阳牛、鲁西牛和安格斯牛共630头个体的单磷酸腺苷激活的蛋白激酶α2亚基(AMPKα2)基因PRKAA2的多态性，并分析其多态性与200头秦川牛生长及肉质性状的关联性。【方法】 利用PCR-SSCP、DNA测序、飞行质谱(MALDI-TOF)等技术，分析检测4个品种牛PRKAA2基因6个外显子的单核酸多态性(SNPs)；采用生物信息学方法结合统计分析软件，分析秦川牛各多态位点的遗传变异特性与体尺及肉质性状的关联性。【结果】 PRKAA2基因第2、4、6、7、8、9外显子中，只检测到第4外显子27 G>C和60 T>C 2个突变位点；关联分析表明，27 G>C位点与秦川牛背膘厚和眼肌面积显著关联(P<0.05)；60 T>C位点与秦川牛体斜长极显著相关(P<0.01)，与体高、尻长、腰角宽、坐骨端宽、眼肌面积显著相关(P<0.05)。【结论】 PRKAA2基因对秦川牛部分体尺及肉质性状有极显著或显著影响，可作为秦川肉牛新品种早期选择的候选基因和分子辅助标记。

[关键词]秦川牛；PRKAA2基因；单核酸多态性(SNPs)；体尺性状；肉质性状

在动物生产中，家畜自开始驯化就受到集中选择，特别是关于生长、发育、繁殖、行为和对疾病的抵抗力等性状[1]。近年来，随着分子育种技术的成熟，标记辅助选择(Marker assisted selection, MAS)方法已成为育种工作者对家畜重要经济性状进行标记选择进而提高育种效率的有效途径[2]。

哺乳动物单磷酸腺苷活化蛋白激酶(AMPK)是一个异源三聚体蛋白质，其组成亚基有α、β、γ 3种[3-6]。已知的AMPK亚型包括α1、α2、β1、β2、γ1、γ2和γ3共7种。从理论上来讲，α、β、γ的不同异构型可形成多达12种组合[5-8]。AMPKα2亚单位的编码基因为PRKAA2，该基因主要在心脏、骨骼肌和肝脏中表达[9]，且大部分位于细胞核内[10]。定位于核内的PRKAA2可以通过部分磷酸化转录因子来达到调节基因表达的作用[11-12]。

AMPKα2是参与AMPK激活的重要组成部分，是脂联素抑制肝脏葡萄糖输出的关键性靶点[13]，PRKAA2基因的缺陷或表达水平降低将直接影响AMPK的作用，从而导致高糖血症及糖尿病的发生。Viollet等[14]研究发现，PRKAA2基因敲除小鼠(PRKAA2-/-)表现为明显的高糖血症和胰岛素抵抗。Beall等[15]认为活化的AMPKα2亚基可能通过上调线粒体解偶联蛋白2(UCP2)的浓度来维持胰岛β细胞对葡萄糖的敏感性。Foretz等[16]和Viana等[17]研究发现，在正常小鼠的肝脏中短期表达活化形式的AMPKα2可以引起轻微的低糖血症，而在糖尿病小鼠模型中同样的处理可使原有的高糖血症有所缓解。

牛PRKAA2基因(GenBank登录号：AC_000160)位于20号染色体，基因全长71 908 bp，包含8个内含子和9个外显子。目前，对该基因多态性的研究多集中在人上。Horikoshi等[18]研究发现，PRKAA2多态性与日本人群胰岛素抵抗和Ⅱ 型糖尿病有关。刘靖星[19]以甘肃兰州人群为对象，研究了PRKAA2基因多态性与Ⅱ型糖尿病及脂联素、抵抗素(Resistin)的相关性，发现病例组和对照组中均存在AMPK-rs2796516位点G/A基因突变，Ⅱ型糖尿病患者AMPK-rs2796516 G/A基因多态性与胰岛素抵抗和血脂代谢明显相关(P<0.05)。宋娇等[20]在肉鸡上的研究发现，PRKAA2基因表达水平与肌内脂肪(IMF)沉积显著相关(P<0.05)，1日龄AA鸡PRKAA2基因表达水平显著高于56日龄(P<0.05)，IMF含量1日龄显著高于56日龄(P<0.05)，PRKAA2基因表达与IMF含量呈正相关。高丽南[21]研究发现，西杂牛、利杂牛、夏杂牛PRKAA2基因第4外显子55 513 bp处均发生了C>T突变，在第6外显子 59 233 bp处发生了A>G突变，但关联分析结果发现，这3种肉牛PRKAA2基因第4外显子和第6外显子突变位点不同基因型与肉牛个体脂肪含量、体质量、屠宰率等性状均无显著差异(P>0.05)。

本试验以秦川牛、南阳牛、鲁西牛和安格斯牛等4个牛品种作为研究对象，探究其PRKAA2基因单核苷酸多态性(SNP)，并分析其多态性与秦川牛体尺和肉质性状是否相关，以期寻求该基因与秦川牛体尺和肉质性状密切相关的遗传标记位点，为肉牛新品种(系)的培育提供理论依据和技术方法。

1材料与方法

1.1主要试剂和仪器

蛋白酶K购自德国默克公司；Tris饱和酚、EDTA、Tris、硼酸、甲醛、丙烯酰胺(Acrylamide)、N,N′-亚甲基双丙烯酰氨(Bis)、二甲基亚砜(DMSO)、N,N,N,N′-四甲基乙二胺(TEMED)、含染料Reaction MIX、DNA Marker、琼脂糖等，均购自天根生化科技(北京)有限公司；十二烷基磺酸钠(SDS)和去离子甲酰氨，购自北京鼎国生物技术有限公司。

PCR仪、凝胶成像系统(Bio-Rad公司，美国)，冷冻高速离心机(CENTURION，英国)，电泳槽、稳流稳压电泳仪、脱色摇床(六一仪器，北京)，Aquila Vet CE0344型兽用B超仪(Pie公司，荷兰)。

1.2牛血样的采集及其体尺和肉质性状的测量

从4个牛品种群体中随机抽取个体，每头牛颈静脉采血10 mL，ACD抗凝(V(血液)∶V(ACD)=6∶1)，轻微颠倒混匀，-80 ℃保存备用。共采取血液样品630份，其中秦川牛364份，血样采自陕西秦川牛业有限公司杨凌秦川肉牛良种繁育中心；南阳牛84份，血样采自河南省南阳市南阳黄牛繁殖场；鲁西牛102份，血样采自山东鲁西牛保种场；安格斯牛80份，血样采自陕西秦宝牧业发展有限公司眉县繁育场。从秦川牛中选取饲养条件相同的24～30月龄健康育肥阉牛200头，测量其8项体尺指标(体斜长、体高、腰高、尻长、腰角宽、胸深、胸围和坐骨端宽)，并采用B超仪测定3项肉质指标(背膘厚、眼肌面积和肌内脂肪)。

1.3血液DNA的提取及检测

采用酚-氯仿法[22]从血样中提取牛基因组DNA，用质量分数1.5%琼脂糖凝胶电泳检测DNA质量。DNA无降解的用紫外分光光度计测定其质量浓度，稀释至50 ng/μL，-20 ℃保存备用。

1.4牛PRKAA2基因外显子的多态性检测

1.4.1引物的设计与合成通过查阅NCBI上提供的牛PRKAA2基因序列(GenBank登录号：AC_000160)，利用Primer 5.0软件分别针对外显子2、4、6、7、8、9设计P1、P2、P3、P4、P5、P6共6对特异性引物(见表1)，交由上海(生工)生物工程股份有限公司合成。

表 1　牛PRKAA2基因外显子的引物设计Table 1　Primers of PRKAA2 gene in cattle

1.4.2PCR扩增反应由15 μL体系组成：其中dNTPs(10 mmol/L)0.5 μL、TaqDNA 聚合酶(2.5 U/μL)0.15 μL、10×Buffer(含15 mmol/L Mg2+)的MIX 6.85 μL、5 μmol/L混合引物(上游引物和下游引物均为10 pmol/μL) 各0.3 μL、模板DNA (50 ng/μL) 1.0 μL、ddH2O 5.9 μL。PCR反应程序：95 ℃预变性5 min；94 ℃变性 30 s，退火30 s(退火温度见表1)，72 ℃ 延伸35 s，35个循环；72 ℃延伸10 min，4 ℃保存。PCR扩增产物用质量分数 1.5%琼脂糖凝胶电泳检测。

1.4.3PCR扩增产物的回收及纯化使用DNA回收试剂盒进行PCR扩增产物回收纯化，按照试剂盒具体说明进行操作，纯化产物送交上海(生工)生物工程股份有限公司测序。

1.4.4PCR-SSCP分析取PCR产物5 μL与10 μL变性缓冲液( 按照去离子甲酰胺9.9 mL，5 mol/L EDTA 100 μL，二甲苯青0.012 5 g，溴酚蓝0.012 5 g配制，调pH 为8.0)混合于PCR反应管，98 ℃变性10 min，迅速将管插入冰水混合物中冰浴5～10 min，然后加样于12%非变性聚丙烯酰胺凝胶中，先进行15 min 250 V高压电泳，再115 V电泳16 h。电泳结束后，将凝胶取出，用超纯水漂洗 2～3 次，然后加0.1%硝酸银溶液浸没凝胶，在脱色摇床上染色10～15 min；弃染色液，超纯水漂洗2～3次，加少量20 g/L NaOH与甲醛混合显色液(500 mL 20 g/L NaOH溶液和5 mL甲醛)预显色30 s，弃显色液，然后再加新鲜显色液至凝胶浸没，摇动直至显现清晰条带为止。

1.4.5飞行质谱(MALDI-TOF)选取质量和质量浓度(20～50 ng/μL)均符合 MALDI-TOF飞行质谱要求的DNA，与针对待测位点设计的引物(每个待测位点处设计3条特异性引物，其中2条用于 PCR扩增，1条用于延伸)一起，送上海(生工)生物工程股份有限公司进行MALDI-TOF分析。

1.5数据统计分析

运用遗传多样性分析软件(Genpop 32)统计各突变位点的等位基因频率和基因型频率，同时计算多态信息含量(PIC)、有效等位基因数(Ne)、遗传杂合度(He)和遗传纯合度(Ho)[23-24]。其中，PIC>0.5说明该突变位点处于高度多态；0.25

2结果与分析

2.1PRKAA2基因外显子2、6、7、8多态位点的检测

2.1.1PRKAA2基因外显子2、6、7、8 PCR扩增结果对4品种牛PRKAA2基因的外显子 2、6、7、8进行PCR扩增，所得产物用1.5%琼脂糖凝胶电泳检测，扩增片段长度与预期扩增产物片段长度一致，条带清晰特异性好(以秦川牛为例的结果见图1)，可直接进行SSCP 检测。

图 1　秦川牛PRKAA2基因外显子 2、6、7、8的PCR扩增结果 A.外显子2;B.外显子6;C.外显子7;D.外显子8；M.DNA Marker DL 2000；1～5.PCR扩增产物Fig.1　Amplified products of PRKAA2 exon 2,exon 6,exon 7 and exon 8 in Qinchuan cattle A.Exon 2;B.Exon 6;C.Exon 7;D.Exon 8；M.DNA Marker DL 2000；1-5.Products of PCR

2.1.2PRKAA2基因外显子2、6、7、8 的PCR-SSCP分析4品种牛PRKAA2基因外显子 2、6、7、8的 PCR扩增产物经SSCP检测，结果(以秦川牛为例的结果见图2)在各牛群中均未发现多态。

图 2　秦川牛PRKAA2基因外显子 2、6、7、8 的PCR-SSCP检测结果 A.外显子2;B.外显子6;C.外显子7;D.外显子8Fig.2　PCR-SSCP detection of PRKAA2 exon 2,exon 6,exon 7 and exon 8 in Qinchuan cattle A.Exon 2;B.Exon 6;C.Exon 7;D.Exon 8

2.2PRKAA2基因外显子4、9的多态位点检测

2.2.1PRKAA2基因外显子4、9 的PCR扩增及突变位点检测秦川牛、南阳牛、鲁西牛和安格斯牛4个群体中每个群体随机挑选10头个体的DNA样品，分别对其PRKAA2基因的外显子4、9进行PCR扩增，所得产物用1.5%琼脂糖凝胶电泳检测，PCR扩增片段长度与预期扩增产物片段长度一致，条带清晰、特异性好(以秦川牛为例的结果见图3)。扩增产物经胶回收、纯化后，送上海(生工)生物工程股份有限公司测序。测序结果用DNA STAR比对，在PRKAA2基因外显子 4上发现了2个潜在的多态位点，分别是第27位颠换G>C和第60位转换T>C。将这2个多态位点分别命名为27 G>C、60 T>C。而在第9外显子区域未发现突变位点。经分析，外显子 4上的27 G>C、60 T>C突变均未引起氨基酸的变化，属同义突变。27 G>C位点测序得到GG、GC 2种基因型(以秦川牛为例的结果见图4)，60 T>C位点测序得到TT、TC 2种基因型(以秦川牛为例的结果见图5)。

图 3　秦川牛PRKAA2 外显子4、9 PCR扩增结果 A.外显子4;B.外显子9；M.DNA Marker DL 2000；1～6.PCR扩增产物Fig.3　Amplified products of PRKAA2 exon 4 and exon 9 in Qinchuan cattle A.Exon 4;B.Exon 9；M.DNA Marker DL 2000；1-6.Products of PCR

图 4　秦川牛PRKAA2基因外显子4第27位点GG、GC基因型测序图(部分结果)Fig.4　Sequencing results of genotype GG and GC in PRKAA2-exon 4 at the 27th locus in Qinchuan cattle(Partial result)

图 5　秦川牛PRKAA2基因外显子4第60位点TT、TC基因型测序图(部分结果)Fig.5　Sequencing results of genotype TT and TC in PRKAA2-exon 4 at the 60th site in Qinchuan cattle(Partial result)

2.2.2PRKAA2基因外显子4多态位点飞行质谱分型检测及验证因PRKAA2基因第4外显子2个潜在突变位点距离太近，依靠PCR-SSCP等常规电泳试验技术难以分型，故采用飞行质谱技术对这2个位点进行进一步的基因分型检测。

(1)PRKAA2基因外显子4多态位点飞行质谱检测引物设计。对第4外显子27 G>C和60 T>C突变位点分别设计特异性引物(见表2)。

(2)质谱分型验证及检测结果。飞行质谱结果除检测到27 G>C和60 T>C突变位点通过测序确定的2种基因型外，还检测到2个位点测序没有发现的第3种基因型，即CC基因型，因此，飞行质谱法得到27 G>C和60 T>C 2个位点全部基因型(图6)。

表 2　牛PRKAA2基因外显子4多态位点飞行质谱引物设计Table 2　MALDI-TOF primer sequences of PRKAA2 gene exon 4 in cattle

图 6　秦川牛PRKAA2基因第4外显子多态位点飞行质谱图 A2-1为27 G>C位点；A2-2为60 T>C位点Fig.6　MALDI-TOF maps of SNPs of Qinchuan cattle PRKAA2 gene exon 4 A2-1 is 27 G>C locus and A2-2 is 60 T>C locus

2.3PRKAA2基因外显子4 SNP位点遗传多态性分析

2.3.1等位基因频率、基因型频率及H-W检验PRKAA2基因第4外显子 2个SNPs 的基因型频率和等位基因频率统计结果见表3和表4。从表3和表4可以看出，2个SNPs 在4个牛群体中均以普通纯合基因型频率最高，其次为杂合基因型，突变型最低，说明这2个位点普通纯合型为优势基因型。

表 3　PRKAA2基因第4外显子 27 G>C SNPs位点基因型及等位基因频率Table 3　Genotypic and allelic frequencies for 27 G>C locus in exon 4 of PRKAA2 gene

注：括号中的数据为观测个体数。表4同。

Note：Data in brackets are the number of observations.The same for Table 4.

表 4　PRKAA2基因第4外显子60 T>C SNPs位点基因型及等位基因频率Table 4　Genotypic and allelic frequencies for 60 T>C locus in exon 4 of PRKAA2 gene

2.3.22个SNPs 位点在群体中的遗传特性分析PRKAA2基因外显子4中2个SNPs 在群体中的纯合度、杂合度、有效等位基因数、PIC分析统计结果见表5。由表5可以看出，试验牛群在2个位点的群体基因杂合度在0.304 7～0.397 7，有效等位基因数在1.438 2～1.660 2，PIC为0.258 3～0.318 6，均介于0.25～0.50，表明4个牛群体在这2个位点均处于中度多态。经χ2检验，秦川牛、南阳牛、鲁西牛、安格斯牛在27 G>C和60 T>C 2个位点的χ2值均小于5.991，说明4个牛群体在这2个位点均处于H-W平衡状态(P>0.05)。从总群体来看，4个牛群体在这2个位点均处于 H-W平衡状态。

表 5　PRKAA2基因外显子4中2个SNPs位点在各群体中的遗传特性Table 5　Genetic indices of 2 SNPs for PRKAA2 gene exon 4 in every population

2.4秦川牛PRKAA2基因外显子4多态性与体尺和肉质性状的关联分析

2.4.1与体尺性状的关联分析考虑场地、年龄和基因型的效应，利用 SPSS 软件的 GLM 模型分析秦川牛PRKAA2基因多态性与其体尺性状的关联性。结果(表6)显示，PRKAA2基因第4外显子27 G>C多态位点与秦川牛8项体尺性状指标未达到显著相关(P>0.05)。60 T>C多态位点与秦川牛体斜长极显著相关(P<0.01)，与体高、尻长、腰角宽、坐骨端宽显著相关(P<0.05)。具体表现为体斜长性状上，CC基因型个体极显著高于TC和TT基因型个体(P<0.01)，TC、TT基因型个体之间差异不显著(P>0.05)；在尻长、腰角宽、坐骨端宽3项体尺性状指标上，均是CC基因型个体显著高于TC和TT基因型个体(P<0.05)，TC与TT基因型个体之间差异不显著。

表 6　PRKAA2基因外显子4中2个SNPs位点与秦川牛体尺性状的关联性Table 6　Association of 2 SNPs locus for exon 4 of PRKAA2 gene with growth traits in Qinchuan cattle 　cm

注：同列数据后标不同小写字母表示差异显著(P<0.05)，标不同大写字母表示差异极显著(P<0.01)。下表同。

Note:Data with a different superscript capital letters in one column differ significantly at the level of 0.01,small letters at the level of 0.05.The following table is same.

2.4.2与肉质性状的关联分析利用 SPSS软件的 GLM 模型分析秦川牛PRKAA2基因多态性与肉质性状的关联性，结果(表7)显示，PRKAA2基因外显子4上27 G>C与背膘厚和眼肌面积2项指标显著相关，GG、GC基因型显著大于CC基因型(P<0.05)；60 T>C位点与秦川牛眼肌面积显著相关，CC基因型显著大于TC基因型(P<0.05)，CC与TT、TT与TC基因型之间差异不显著(P>0.05)。

表 7　PRKAA2基因外显子4中2个多态位点与秦川牛肉质性状的关联性Table 7　Association of 2 SNP locus for exon 4 of PRKAA2 gene with meat quality traits in Qinchuan cattle

3讨论与结论

在真核生物中，AMPK是一个重要的蛋白激酶，是哺乳动物激酶级联反应的核心组成部分，主要协调能量代谢的需要，通过阻断三磷酸腺苷(ATP)消耗和促进ATP再生的生物合成途径来防止哺乳动物细胞代谢应激，故有“细胞能量调节器”之称[4-7,25]。

本研究在秦川牛、南阳牛、鲁西牛和安格斯牛4个品种肉牛PRKAA2基因第2、6、7、8、9外显子未检测到突变位点，只在第4外显子中均发现了27 G>C和60 T>C 2个突变，且均未引起相应氨基酸的变化，属于同义突变。同义SNP是进化过程中对环境压力的适应和自然选择而形成的“密码子使用偏好”在基因组中的遗迹，有研究表明它能通过影响翻译效率[26]、mRNA的稳定性[27]和拼接控制来调节机体生理活动[28]，或者由于密码子替换导致翻译速度发生改变，进而通过影响蛋白折叠的动力学控制使蛋白质无法完成正确的折叠[29]。因此，同义突变从某种方面也能影响基因的表达和调节，进而影响个体表型性状。

在4个牛群体中，27 G>C、60 T>C 2个突变位点均以普通纯合基因型频率最高，突变型频率最低，2个位点的等位基因频率表现同一趋势，均是 G>C、T>C。对该基因的研究，单从检测得到的SNP数量与Zhang 等[30]的研究结果比较，本研究只在第4外显子检测到2个SNPs，与Zhang等[30]在PRKAA2基因扫描的外显子区域不同。

关联分析显示，27 G>C位点与秦川牛背膘厚和眼肌面积显著关联(P<0.05)，60 T>C位点与秦川牛体斜长极显著相关(P<0.01)，与体高、尻长、腰角宽、坐骨端宽、眼肌面积显著相关。本研究的部分结果与Lin 等[31]在猪上得到的结果相一致，但与高丽南[21]在牛上报道的结果有较大差异，这可能与试验牛品种差异有一定关系。

由本研究结果可知，PRKAA2基因可作为秦川牛标记辅助育种的候选基因，为秦川肉牛应用生物技术育种提供了新的思路和手段。然而因牛屠宰成本太高，本研究分析的数据和肉质性状相对偏少，还有待进一步扩大样本数量对其他肉质性状进行分析，并在细胞水平更深层次研究该基因功能及其相关调控机理。

[参考文献]

[1]Andersson L,Georges M.Domestic animal genomics:deciphering the genetics of complex traits [J].Nat Rev Genet,2004,5(3):202-212.

[2]Wu X L,MacNeil M D,De S,et al.Evaluation of candidate gene effects for beef backfat via Bayesian model selection [J].Genetica,2005,125(1):103-113.

[3]Hardie D G,Carling D.The AMP-activated protein kinase——fuel gauge of the mammalian cell [J].Eur J Biochem,1997,246:259-273.

[4]Hardie D G,Salt1 P,Hawley S A,et al.AMP-aetivated proteinkinase:an utrasensitive system for monitoring cellar energy change [J].Biochem J,1999, 338(3):717-722.

[5]Hardie D G,Carling D,Carlson M.The AMP-activated/SNF1 protein kinase subfamily:metabolic sensors of the eukaryotic cell [J].Annu Rev Biochem,1998,67:821-855.

[6]Hardie D G,Sakamoto K.AMPK:a key sensor of fuel and energy status in skeletal muscle [J].Physiology (Betllesda),2006,21:48-60.

[7]Wadley G D,Lee-Young R S,Canny B J,et al.Effect of exercise intensity and hypoxia on skeletal muscle AMPK signaling and substrate metabolism in humans [J].Am J Physiol Endocrinol Metab,2006,290(4):E694-E702.

[8]Yu H,Fujii N,Hirshman M F,et al.Cloning and characterization of mouse 5′-AMP-activated protein kinase γ3 subunit [J].Am J Physiol Cell Physiol,2004,286:C283-C292

[9]Stapleton D,Mitehelhill K I,Gao G,et al.Mammalian AMP-activated protein kinase subfamily [J].J Biol Chem,1996,271:611-614.

[10]Turnley A M,Stapleton D,Mann R J,et al.Cellular distribution and development expression of AMP-activated protein kinase isoform in mouse central nervous system [J].J Neurochem,1999,72:1707-1716.

[11]Dyck J R,Gao G,Widmer J,et al.Regulation of 5′-AMP-activated protein kinase activity by the noneatalytie beta and gamma subunits [J].J Biol Chem, 1996, 271: 17798-17803.

[12]Dyck J R,Kudo N,Barr A J,et al.Phosphorylation control of cardiac acetyl-CoA carboxylase by cAMP-dependent protein kinase and 5′-AMP activated protein kinase [J].Eur J Biochem,1999,262(1):184-190.

[13]Andreelli F,Foretz M,Knauf C,et al.Liver adenosine monophosphate-activated kinase-alpha2 catalytic subunit is a key target for the control of hepatic glucose production by adiponectin and leptin but not insulin [J].Endocrinology,2006,147:2432-2441.

[14]Viollet B,Andreelli F,Jorgensen S B,et al.The AMP-activated protein kinase alpha2 catalytic subunit controls whole-body insulin sensitivity [J].Clin Invest,2003,111(1):91-98.

[15]Beall C,Piipari K,Al-Qassab H,et al.Loss of AMP-activated protein kinase alpha2 subunit in mouse beta-cells impairs glucose-stimulated insulin secretion and inhibits their sensitivity to hypoglycaemia [J].Biochem J,2010,429(2):323-333.

[16]Foretz M,Ancellin N,Andreelli F,et al.Short-term overexpression of a constitutively active form of AMP-activated protein kinase in the liver leads to mild hypoglycemia and fatty liver [J].Diabetes,2005,54(5):1331-1339.

[17]Viana A Y,Sakoda H,Anai M,et al.Role of hepatic AMPK activation in glucose metabolism and dexamethasone-induced regulation of AMPK expression [J].Diabetes Res Clin Pract,2006,73(2):135-142.

[18]Horikoshi M,Hara K,Ohashi J,et al.A polymorphism in the AMPKα2 subunit gene is associated with insulin resistance and type 2 diabetes in the Japanese population [J].Diabetes,2006,55(4):919-923.

[19]刘靖星.AMPKα2基因多态性与2型糖尿病及脂联素、抵抗素相关性研究 [D].兰州：兰州大学,2009.

Liu J X.Correlation between AMPKα2 subunit gene polymorphism and type2 diabetes mellitus,adiponectin and resistin [D].Lanzhou:Lanzhou University,2009.

[20]宋娇,付睿琦,李颖颖,等.肉鸡PRKAA基因表达与MyHC基因表达和肌内脂肪沉积关系的研究 [J].中国畜牧杂志,2012,48(23):10-13.

Song J,Fu R Q,Li Y Y,et al.The study of relationship betweenMyHCgene expression and intramuscular fat deposition in broilerPRKAAgene expression [J].Chinese Journal of Animal Science,2012,48(23):10-13.

[21]高丽南.AMPK对肉牛脂肪沉积的影响及其基因多态性研究 [D].呼和浩特:内蒙古农业大学,2012.

Gao L N.Research on AMPK and its gene polymorphism influence on beef fat deposition [D].Hohhot:Inner Mongolia Agricultural University,2012.

[22]J萨姆布鲁克,D W拉塞尔.分子克隆实验指南(上、下册) [M].黄培堂译.北京:科学出版社,2002.

Sambrook J,Russell D W.Molecular cloning:a laboratory manual (Volume 1,2) [M].Translated by Huang P T.BeiJing:Science Press,2002.

[23]Nei M,Li W H.Mathematical model for studying genetic variation in terms of restriction endonucleases [J].Proceedings of the National Academy of Sciences,1979,76(10):5269-5273.

[24]Nei M,Roychoudhury A K.Sampling variances of heterozygosity and genetic distance [J].Genetics,1974,76(2):379-390.

[25]Carling D.AMPK [J].Curr Biol,2004,14(6):R220.

[26]Lavner Y,Kotlar D.Codon bias as a factor in regulating expression via translation rate in the human genome [J].Gene,2005,345(1):127-138.

[27]Chamary J V,Hurst L D.Evidence for selection on synonymous mutations affecting stability of mRNA secondary structure in mammals [J].Genome Biology,2005,6(9):R75.1-R75.12.

[28]Eskesen S T,Eskesen F N,Ruvinsky A.Natural selection affects frequencies of AG and GT dinucleotides at the 5′and 3′ends of exons [J].Genetics,2004,167(1):543-550.

[29]Purvis I J,Bettany A J E,Santiago T C,et al.The efficiency of folding of some proteins is increased by controlled rates of translationinvivo:a hypothesis [J].Journal of Molecular Biology,1987,193(2):413-417.

[30]Zhang Q,Zhao S,Chen H,et al.SNP discovery and haplotype analysis in the bovinePRKAA2 gene [J].Mol Biol Rep,2011,38:1551-1556.

[31]Lin L,Flisikowski K,Schwarzenbacher H,et al.Characterization of the porcine AMPK alpha 2 catalytic subunitgene (PRKAA2):genomic structure,polymorphism detection and association study [J].Anim Genet,2010,41:203-207.

DOI：网络出版时间:2016-06-0816:2010.13207/j.cnki.jnwafu.2016.07.001

[收稿日期]2014-12-05

[基金项目]国家863计划项目(2013AA102505)；国家自然科学基金项目(31272411)；陕西省自然科学基金项目(2007C127)；陕西省科学技术发展计划项目(2011K01-09)

[作者简介]田万强(1969-),男,陕西扶风人,副教授，博士，主要从事家畜遗传育种与繁殖研究。E-mail:twqiang2003@163.com [通信作者]昝林森(1963-)，男，陕西扶风人，教授，博士，博士生导师，主要从事动物遗传育种及生殖生理调控研究。 E-mail:zanlinsen@163.com

[中图分类号]S823.9+2

[文献标志码]A

[文章编号]1671-9387(2016)07-0001-09

SNPs of bovinePRKAA2 gene and its association with growth and meat quality traits

TIAN Wanqiang1,2,MEI Chugang2,FU Changzhen2,XIN Yaping2,ZHANG Song2,WANG Guoqing2,ZAN Linsen2,3

(1AnimalEngineeringDepartment,YanglingVocational&TechnicalCollege,Yangling,Shaanxi712100,China;2CollegeofAnimalScienceandTechnology,NorthwestA&FUniversity,Yangling,Shaanxi712100,China;3NationalBeefCattleImprovementCenter,Yangling,Shaanxi712100,China)

Abstract:【Objective】 Polymorphism of PRKAA2 gene on 630 individuals from four cattle breeds (Qinchuan,Nanyang,Luxi and Angus) was detected to analyze its genetic structures and varieties.The association of different genotypes with growth and meat quality traits of 200 Qinchuan cattle was also analyzed. 【Method】 SNPs located in exon 2,exon 4,exon6,exon 7,exon 8 and exon 9 were detected in PRKAA2 gene of four cattle breeds by PCR-SSCP,random DNA sequencing method and MALDI-TOF technology.Genetic variation of SNPs and the association of different genotypes with growth and meat quality traits of Qinchuan cattle were also analyzed by bioinformatics method and SPSS software.【Result】 Two SNPs (27 G>C and 60 T>C) located in only exon 4 were detected in 6 exons of PRKAA2 gene.The 27 G>C SNP was associated with back fat thickness (BF) and loin muscle area (LMA) (P<0.05) significantly.The 60 T>C SNP was associated with body length (BL) (P<0.01) extremely significantly,and with withers height (WH),rump length (RL),hip width (HW),pin bone width (PBW) and LMA (P<0.05) significantly.【Conclusion】 PRKAA2 gene had potential effects on growth and meat quality traits of Qinchuan cattle and could be used for candidate gene and marker-assisted selection (MAS) in early selection of new varieties of Qinchuan cattle.

Key words：Qinchuan cattle;PRKAA2 gene,SNPs;growth traits;meat quality traits

网络出版地址:http://www.cnki.net/kcms/detail/61.1390.S.20160608.1620.002.html