张 娟，虎红红，邓占钊，母 童，顾亚玲，辛国省

(1宁夏大学农学院，宁夏 银川750021；2彭阳县畜牧技术推广服务中心，宁夏 彭阳756599；3宁夏大学生命科学学院/宁夏饲料工程技术研究中心，宁夏 银川750021)

静原鸡主要分布于宁夏彭阳县和甘肃静宁县，2006年被农业部认定为地方优良品种，已被列入国家级畜禽遗传资源保护名录[1]，根据羽色的不同，分为白羽、麻羽和黑羽3个群体。由于长期处于气候高寒湿润、地形起伏、牧草及昆虫众多的生活环境和优越的地理位置，静原鸡具有耐粗饲、抗逆性强、氨基酸含量丰富、肉质细嫩、营养和滋补价值高等特点，尤其多不饱和脂肪酸(二十碳五烯酸、二十二碳六烯酸)含量较高[2]，符合现今人们对高品质鸡肉的要求，属鸡肉中优选的绿色食品[3]。

腺苷酸激酶(Adenylate kinase,AK)是一种普遍存在于动植物体内的单体酶，对生物的生长和维持至关重要[4]，主要调节腺嘌呤核苷酸代谢，催化反应：ATP+AMP↔2ADP，在维持细胞能量平衡中起着不可或缺的作用[5]。在哺乳动物中，有几种AK 亚型具有组织特异性分布和独特的亚细胞定位[6]，其中AK1蛋白被认为是骨骼肌细胞质中的主要AK亚型[7]。猪心脏中胞质酶AK1蛋白的三维结构已经被Schulz等[8]和Sachsenheimer 等[9]准确测定。有研究表明，AK1基因缺失的小鼠尽管肌肉形成正常，但由于ADP的异常积累会导致骨骼肌松弛延迟[10]。此外，还有研究发现AK1蛋白磷酸化可能是质膜线粒体和KATP 通道之间传递信号所必需的[11]。目前关于AK1基因的研究集中在与溶血性贫血相关的罕见遗传疾病中，国内外将AK1作为肌苷酸影响基因及有关AK1基因结构及功能预测的研究鲜有报道。

本研究通过对宁夏优良地方品种静原鸡AK1基因CDS区序列进行扩增、克隆、构建其真核表达载体，运用生物信息学分析方法预测其理化性质、亲水性/疏水性、跨膜结构域、磷酸化位点、亚细胞定位、磷酸化位点、B细胞抗原表位、基因共表达、Cp G岛、蛋白的空间结构等，为后续深入研究AK1基因蛋白水平、代谢水平和细胞功能研究等的调控机理和表达模式提供科学依据，并为我国地方品种鸡的分子育种技术平台奠定基础。

1 材料与方法

1.1 试验材料与试剂

随机屠宰宁夏彭阳县朝那鸡繁育中心180日龄静原鸡公、母各15只(饲养管理条件相同)，采集腿肌、胸肌于液氮罐中，带回实验室于−80℃保存备用。试剂：Trizo，Hifair®Ⅱ1st Strand c DNA Synthesis Super Mix for qPCR(gDNA digester plus)试剂盒(上海翊圣生物公司)，DNA loading buffer(百泰克)，核酸染料(索莱宝)，p MD18-T、DH5α、质粒小提中量试剂盒、普通琼脂糖凝胶DNA 回收试剂盒(天根)，T 4连接酶(NEB)，EcoRⅠ、BamHⅠ(宝生物有限公司)。

1.2 试验方法

1.2.1总RNA 的提取及反转录 通过传统Trizo(Invitrogen,USA)法提取静原鸡腿肌、胸肌的总RNA，并对其RNA 的纯度、完整性、浓度进行检测。采用Hifair®Ⅱ1st Strand c DNA Synthesis SuperMix for qPCR(gDNA digester plus)试剂盒合成cDNA 后，于−20℃保存备用。

1.2.2AK1基因的引物设计 根据GenBank 上已公布的原鸡AK1基因序列，利用Primer 5.0软件设计特异性引物(表1)，并送至生工生物工程(上海)股份有限公司合成。

表1 静原鸡AK1基因CDS区引物序列Table1 Primer sequencefor CDSregion in AK1 gene of Jingyuan chicken

1.2.3AK1基因CDS区扩增、克隆及测序PCR扩增总体系为20μL：cDNA 样品1μL，上、下游引物各0.5μL，TaqPCR Master Mix 11μL，ddH2O 7μL。AK1基因CDS区PCR 扩增程序为：94℃预变性5 min；94℃变性45 s，退火40 s，72℃延伸45 s，共31个循环；72℃延伸10 min，4℃保存。PCR产物通过10 g/L 的琼脂糖凝胶电泳检测后，使用普通琼脂糖凝胶DNA 回收试剂盒对目的片段进行回收。

将回收纯化好的目的片段与pMD18-T 载体4 ℃下连接并过夜。将连接好的带有pMD18-T 载体的目的片段转化到DH5α 感受态细胞中，冰浴30 min，42 ℃热激90 s，冰浴3 min，加入无抗液体LB培养基，混匀后置于37℃摇床150 r/min 振荡培养45 min，涂布于含氨苄抗性的LB平板培养基上，37℃条件下培养12～16 h，随机挑取5个单菌落放入含氨苄抗性的LB液体培养基220 r/min 震荡培养8 h，菌液PCR 扩增。收集菌液送至生工生物工程(上海)股份有限公司进行测序。

1.2.4A K 1载体的构建 引物设计合成见“1.2.2”，扩增回收同“1.2.3”，选择p EGFPN1真核表达载体进行质粒提取。利用EcoRⅠ、BamHⅠ对AK1基因与pEGFP-N1质粒进行双酶切，双酶切体系为50μL，具体加入量分别为：11μL dd H2O，5μL 10×Buff er，2μLEcoRⅠ，2μLBamHⅠ，30μL 重组质粒(目的片段)。37℃酶切6 h 后使用普通琼脂糖凝胶DNA 回收试剂盒回收pEGFP-N1载体片段和AK1目的片段，AK1基因与pEGFP-N1载体的酶切回收片段经过T4连接酶连接后转化。AK1-pEGFP-N1载体进行双酶切鉴定后送至生工生物工程(上海)股份有限公司测序。

1.2.5AK1基因生物信息学分析 采用BioEdit 软件对AK1基因的核苷酸序列进行比对，分析碱基组成、氨基酸数目及组成；相对分子质量和等电点运用在线网站expasy 的protparam 程序(http://web.ex pasy.org/protparam/)进行分析；蛋白质亲水性运用在线网站expasy 的protscale程序(https://web.ex pasy.org/protscale/)进行分析；蛋白质三级结构运用在线网站swiss-model(http://swissmodel.expasy.org/)进行预测；蛋白质跨膜结构域运用在线网站

TMHMM-2.0(http://www.cbs.dtu.dk/services/TMHMM-2.0/)进行预测；磷酸化位点运用在线网站NetPhos(http://www.cbs.dtu.dk/services/NetPhos/)进行预测；B细胞抗原表位运用在线网站immuneepitope的bcell程序(http://tools.immuneepitope.org/main/bcell/)进行分析；蛋白质二级结构运用在线网站psipred(http://bioinf.cs.ucl.ac.uk/psipred/)进行预测；基因共表达和GO功能注释分析通过在线网站string(https://string-db.org/)进行分析；CpG 岛运用在线网站novopro(http://www.novopro.cn/tools/cpg_islands.html)进行预测；基因的亚细胞定位通过在线网站psort Ⅱ(https://psort.hgc.jp/form2.html)进行预测。

2 结果与分析

2.1 总RNA 提取

静原鸡胸肌、腿肌总RNA 通过10 g/L的琼脂糖凝胶电泳检测，结果如图1所示，可清晰看到3个条带，无明显降解。核酸蛋白仪检测RNA 的D260nm/D280nm在1.8～2.2之间，完整性较好(RIN≥6.5)，可用于后续试验。

图1 静原鸡不同组织总RNA 提取Fig.1 Total RNA extr action fr om differ ent tissues of Jingyuan chicken

2.2 AK1 基因CDS区扩增及克隆

PCR 产物经10 g/L的琼脂糖凝胶电泳检测，经凝胶成像仪检测发现扩增片段条带整齐明亮且无引物二聚体等杂带，表明引物特异性较好，扩增条件合适。CDS区扩增产物条带长度大小约为741 bp，与预期目的片段长度大小相符(图2)，可以进行回收纯化；回收后产物连接pMD18-T载体进行克隆，菌液PCR 鉴定片段大小与预期目的片段一致。

2.3 AK1 基因加酶切位点扩增

AK1基因CDS区加酶切位点扩增产物条带长度大小约为759 bp，与预期目的片段长度大小相符(图3)，可以进行回收纯化及菌液扩增。

图2 静原鸡AK1 基因PCR 扩增电泳图Fig.2 Electrophoresis of PCR amplification of AK1 gene in Jingyuan chicken

图3 静原鸡AK1 基因加酶切位点PCR 扩增电泳图Fig.3 Electrophoresis of PCR amplification of AK1 gene with restriction site in Jingyuan chicken

2.4 AK1-pEGFP-N1载体的双酶切鉴定

重组质粒AK 1-p EGFP-N 1经EcoRⅠ和BamHⅠ37℃双酶切6 h 后，经10 g/L 的琼脂糖凝胶电泳检测得到4 033和759 bp的2个片段(图4)，一个为质粒pEGFP-N1本身的片段，为4033 bp，一个为目的基因AK1的片段，约为759 bp，酶切结果说明AK1基因与pEGFP-N1载体正确重组。将菌液通过甘油保存后送至生工生物工程(上海)股份有限公司测序。

图4 静原鸡AK1-p EGFP-N1载体酶切电泳图Fig.4 Restrictive enzyme digestion and electrophoresis of AK1-p EGFP-N1 vector in Jingyuan chicken

2.5 AK1基因CDS区生物信息学分析

2.5.1AK1基因编码区碱基组成分析 测序后对静原鸡AK1基因序列进行校正、拼接，AK1基因CDS区全长585 bp，共编码194个氨基酸(图5)。核苷酸组成分析(图6)发现AK1基因CDS区鸟嘌呤(G)核苷酸数量占比最多(33.68%)，胸腺嘧啶(T)核苷酸占比最少(16.41%)，腺嘌呤(A)和胞嘧啶(C)数目相当，且G+C含量(60.00%)高于A+T含量(40.00%)。编码区DNA 单链和双链相对分子质量分别为177991和356652。

图5 AK1 基因编码区全长序列与对应的氨基酸序列Fig.5 Thefull length coding sequence of AK1 geneand corresponding amino acid sequence

图6 AK1基因CDS区核苷酸组成Fig.6 Nucleotidecomposition in CDSregion of AK1 gene

2.5.2AK1基因编码氨基酸序列及蛋白的理化分析 静原鸡AK1基因编码蛋白质的各种氨基酸中赖氨酸(Lys)数目占比最多(11.3%)，天冬酰胺(Asn)和半胱氨酸(Cys)数目占比最少(1.0%)(表2)。氨基

表2 AK1基因编码蛋白的氨基酸组成Table 2 The composition of amino acid coded by AK1 gene

酸残基中正电荷残基总数(Arg+Lys)为32个，负电荷残基总数(Asp+Glu)为29个。理化性质分析结果显示，静原鸡AK1蛋白的原子组成为C962H1566N262O291S7，相对分子质量为21683，理论等电点为8.68，不稳定指数为21.56，脂肪指数为85.88。

2.5.3AK1基因编码蛋白的亲水性/疏水性 使用在线分析网站对静原鸡AK1基因编码蛋白对水分子的亲和力进行研究的结果(图7)表明，第15位缬氨酸(Val)疏水性最强(最高分数1.933)，第105位谷氨酸(Glu)亲水性最强(最低分数−2.589)，总体亲水性均值小于0(−0.469)，说明该蛋白具有较强的亲水区域。

图7 静原鸡AK1蛋白的疏水性/亲水性分析Fig.7 Hydrophobic/hydrophilic analysis for AK1 protein in Jingyuan chicken

2.5.4AK1基因编码蛋白的跨膜结构域和CpG岛预测 采用在线软件TMHMM-2.0的预测结果表明静原鸡AK1基因的编码蛋白不存在跨膜结构域，不是跨膜蛋白。CpG 岛预测显示，AK1基因编码蛋白存在2个CpG 岛。

2.5.5AK1基因编码蛋白的磷酸化位点 静原鸡AK1蛋白可能存在7个丝氨酸(Ser)磷酸化位点，分别位于氨基酸序列上第20、39、50、52、88、137、189位，8个苏氨酸(Thr)磷酸化位点，分别位于第3、24、36、87、136、146、153、179位，3个酪氨酸(Tyr)磷酸化位点，分别位于第35、118、155位(图8)。

2.5.6AK1基因编码蛋白的B细胞抗原表位预测利用immuneepitope在线软件预测AK1基因编码蛋白的抗原表位(图9)，结果发现静原鸡AK1蛋白存在10个B细胞抗原表位，分别位于1、17～25、38、48～56、85～88、97～111、122～124、134～147、155～158、176～181位氨基酸残基。

图8 静原鸡AK1蛋白磷酸化位点Fig.8 Thephosphorylation sites of AK1 protein in Jingyuan chicken

图9 AK1蛋白B 抗原表位预测Fig.9 Prediction of B epitope of AK1 protein

2.5.7AK1基因编码蛋白的二级和三级结构预测二级结构预测结果(图10)显示，静原鸡AK1蛋白含有8.76%的伸展链，46.91%的α–螺旋，44.33%的无规则卷曲。三级结构预测结果仍以α–螺旋和无规则卷曲为主(相似性为86.39%)(图11)。

2.5.8AK1基因共表达分析 通过string数据库对AK1基因进行共表达分析，发现在原鸡中AK1和AMPD1、PKM2存在共表达，共表达系数分别为0.116和0.063。PKM2与AMPD1共表达系数为0.067，AMPD1与AMPD3共表达系数为0.082。在其他有机体中AK1与AMPD3、AMPD1、NTPCR、PNPT1、CD39L1、PKM2、APRT、ADSL、ADK和TPK1均存在共表达。

图10 静原鸡AK1蛋白质二级结构Fig.10 The secondary structure of AK 1 protein in Jingyuan chicken

图11 静原鸡AK1蛋白质预测三级结构Fig.11 The predicted tertiary structure of AK1 protein in Jingyuan chicken

2.5.9AK1基因编码蛋白亚细胞定位 通过在线工具psortⅡ预测AK1基因编码产物的亚细胞定位，其分布在细胞质的可能性为52.2%，分布在细胞骨架的可能性为17.4%，分布在细胞核中的可能性为13.0%，分布在线粒体、过氧化物酶体、内质网中的可能性分别为8.7%、4.3%、4.3%。由此推断，AK1基因编码产物主要在细胞质中发挥生物学作用。

2.5.10AK1基因的GO富集分析 通过string数据库对AK1基因进行GO富集分析，发现AK1基因参与4个分子功能GO条目，参与11个生物学过程GO条目，参与1个细胞组分GO条目(表3)。

表3 AK1基因的GO注释分析Table 3 GO annotation analysis of AK1 gene

3 讨论与结论

随着生命科学的发展，为揭示复杂的基因组信息结构及遗传语言的根本规律就出现了生物信息学，生物信息学为大数据“插上翅膀”，不仅给生命科学、农学、遗传学、细胞生物学和生物医学带来巨大的推动力且具有诱人的前景，它已经渗透到人类基因组研究的各个层面，如基因组、转录组学、蛋白质组学、代谢组学、疾病潜在靶点基因预测以及系统发育分析的研究[12-13]。

腺苷酸激酶(AK)普遍存在于生物体，催化腺嘌呤核苷酸的相互转化：Mg2++ATP+AMP=Mg2++2ADP，在不同哺乳动物组织中已经鉴定出6种1型AK 同工酶(AK1)[14]，它们的主要功能是维持细胞的能量稳态，处理与细胞能量利用相关的代谢信号，并为核酸合成提供核苷酸[15-16]。AK1是一种小型胞质酶，在骨骼肌、红细胞和脑等高周转率细胞中表达，其活性依赖于Mg2+或Mn2+离子，在没有AK1的情况下腺嘌呤核苷酸大量减少[17]。鸡AK1基因约为6000碱基对，由7个外显子组成，Suminami等[18]在鸡AK1和牛线粒体腺苷酸激酶(AK2)之间的点阵图发现，AK1基因可能是通过缺失1个或多个外显子从AK2基因进化而来的。有研究表明鸡AK1的一级结构与猪AK1的一级结构有85%的相似性，因此，提出猪AK1的二级和三级结构[19]也将适用于鸡AK1，而且发现AK1外显子4相对应的区域与AK2A构象的相似性最高，已知该区域有助于ATP结合和催化功能[20]。本研究发现静原鸡AK1基因CDS区共编码194个氨基酸，理论等电点为8.68(＞7)，不稳定指数为21.56(＜40)，脂肪指数为85.88，推断该蛋白是碱性稳定蛋白[21]，蛋白的稳定性会直接影响生物的生命活动和细胞周期，因此，AK 1蛋白具有稳定的细胞周期和生命活动。AK1蛋白总体亲水性均值为−0.469(＜0)，说明该蛋白是可溶性蛋白。信号肽和跨膜结构域预测发现，AK1蛋白不存在信号肽和跨膜结构，说明该蛋白不是一个膜蛋白，不能通过膜受体来发挥作用，对其进行纯化可能较为简单。磷酸化位点分析发现AK1蛋白存在8个苏氨酸(Thr)磷酸化位点、7个丝氨酸(Ser)磷酸化位点和3个酪氨酸(Tyr)磷酸化位点，这些位点可以调节一些细胞的代谢和信号转导，该蛋白可能需要通过翻译后的磷酸化和去磷酸化影响细胞生长、分化以及基因表达[22]，这与Carrasco等[11]发现AK1磷酸化可能是在质膜线粒体和KATP通道之间传递信号所必需的的研究结果一致。蛋白二三级结构预测发现AK 1蛋白以α–螺旋和无规则卷曲为主，且二级结构为混合模型[23]，因此可能通过影响蛋白质肽链中配体和受体结合的活性[24]，进而影响蛋白的功能。赵骁等[25]等研究发现如果β–转角和无规则卷曲在蛋白中占比居多，则会出现抗原表位较多的现象，因此，AK1蛋白抗原表位预测发现其具有10个B细胞抗原表位。基因共表达分析发现，在原鸡中AK1和AMPD1、PKM2存在共表达，而AMPD1和PKM2都是肌苷酸(Inosinc acid,IMP)从头合成途径的关键基因，因此AK1可能存在与AMPD1与PKM2相似的功能，可将AK1基因作为与IMP相关的功能基因进行研究。Cp G 岛常位于基因转录调控区附近，静原鸡AK1基因存在2个CpG 岛，这为后期AK1基因转录区CpG 岛甲基化状态的研究十分重要。亚细胞定位发现AK1蛋白位于细胞质中，与基因进行GO细胞组分分析和Choo等[26]研究结果一致，由此推断该蛋白主要在细胞质中发挥生物学作用。

综上，本研究通过生物信息分析对前期转录组筛选的静原鸡不同部位与IMP相关的关键差异基因AK1进行分析并通过分子生物学的方法构建其真核表达载体AK1-pEGFP-N1。

静原鸡AK1基因CDS区全长585 bp，AK1蛋白含有194个氨基酸，是具有较强的亲水区域的非跨膜蛋白，AK1蛋白可能存在7个丝氨酸(Ser)磷酸化位点，8个苏氨酸(Thr)磷酸化位点，3个酪氨酸(Tyr)磷酸化位点，并且具有10个B细胞抗原表位，2个CpG 岛，预测其在细胞质中发挥作用，与AMPD1、PKM2、ADSL存在共表达，AK1蛋白二级结构和三级结构均以α–螺旋和无规则卷曲为主(相似性为86.39%)。