绵羊脂联素基因（ADIPOQ）多态性及其与生长及胴体性状关联性分析

2019-06-11安清明周辉通吴震洋罗玉柱JonHickford

中国农业科学 2019年10期

安清明，周辉通，吴震洋，罗玉柱，Jon G.Hickford

安清明1,2，周辉通2,3，吴震洋1，罗玉柱2，Jon G.Hickford2,3

（1铜仁学院，中国贵州铜仁 554300；2甘肃农业大学草食动物生物技术重点实验室，中国兰州 730070；3新西兰林肯大学基因标记实验室，新西兰林肯 7646）

【目的】检测绵羊ADIPOQ基因多态性及连锁现状，评估该基因突变对绵羊生长及胴体性状的影响，以期丰富绵羊相关重要经济性状的分子遗传研究基础。【方法】以8个不同品种的商品绵羊为研究对象，利用PCR-SSCP方法检测ADIPOQ基因Exon-1区和Exon-2区变异，利用GLMs模型进行评估该基因突变对绵羊生长及胴体性状的影响。【结果】绵羊ADIPOQ基因Exon-1区和Exon-2区共检测到13个突变位点，其中Exon-2区发现的c.46T/C突变导致编码氨基酸p.Tyr16His转变。Exon-1区等位基因A和B为优势等位基因，Exon-2区等位基因A和D为优势等位基因，且两个区域的等位基因均存在种群差异，大多数品种在两个区域中的变异为中度多态（0.25＜＜0.5），只有美利奴羊在Exon-2区为高度多态（＞0.5），特克塞尔、派伦代和丘陵陶塞特羊在Exon-2区为低度多态（＜0.25）；两段区域间存在的突变位点为高度连锁，且趋向于共同遗传（D’=0.952，r=0.365）。关联分析结果表明，ADIPOQ基因Exon-1区变异对绵羊生长性状存在性别差异，携带等位基因A的公羔具有较低的断尾重、断奶重和断奶前生长速度（＜0.05），而携带等位基因A的母羔却与生长性状无显著关联；携带等位基因B的公羔与生长性状无显著关联，但携带等位基因B的母羔却具有较高的断尾重（＜0.05）；同时发现基因型为BB的公羔个体具有更高的断尾重和断奶重（＜0.05）；胴体性状关联分析结果表明，携带等位基因A的群体具有较低的热胴体重、腰部瘦肉量、后腿瘦肉量和总瘦肉量（＜0.05），携带等位基因B的群体则具有较高的后腿瘦肉量、后腿瘦肉比例和较低的肩部瘦肉比例（＜0.05）；基因型为BB的个体均有较高的热胴体重、腰部瘦肉量、后腿瘦肉量和总瘦肉量（＜0.05）。【结论】绵羊基因的两段区域具有丰富的多态性，Exon-2区域中的c.46T/C为非同义突变。Exon-1区变异影响绵羊的生长性状和胴体性状，淘汰携带等位基因A的个体和选留存在等位基因B的个体、或留存BB基因型的个体和淘汰AA的个体，均可有效改善绵羊后代群体的部分生长性状和胴体性状。

绵羊；ADIPOQ基因；性别差异；生长性状；胴体性状

0 引言

【研究意义】随着人们生活水平的逐步提高，绿色健康的饮食已成为日常生活的主题之一，而羊肉因其独特的营养风味越来越受到人们的偏爱。育种学家们也在根据市场需求逐渐把如何提高羊肉的产量和质量、满足人类追求绿色健康的饮食需求作为品种改良的主要参考因素，其中，从人类营养学角度来分析绵羊的胴体脂肪组成已逐渐成为一个研究热点。本研究通过检测影响绵羊生长、胴体性状的候选基因，评估该基因突变对不同性别绵羊生长性状、屠宰公羔的胴体性状的影响，以期丰富绵羊生长及胴体性状的分子遗传研究基础。【前人研究进展】脂联素，也称脂肪连接蛋白（adiponetin，ADIPOQ或Acrp30），是一种脂肪组织分泌的具有生物活性的脂肪因子，现阶段研究发现其主要与胰岛素耐受性调节、机体肥胖、心血管疾病和II类糖尿病等多种生理疾病相关[1-2]。脂联素由ADIPOQ基因编码，该基因最早在人类3号染色体上发现，总长约17kb，包含3个外显子和2个内含子（其中外显子1不参与蛋白编码），共编码247个氨基酸。编码的ADIPOQ蛋白分子量约为30kD，主要由信号肽、可变区、N-胶原蛋白三螺旋区和C-球状区组成[3-4]。研究发现，ADIPOQ蛋白主要通过与脂联素受体（ADIPOR-1和ADIPOR-2）蛋白结合，调节AMP激酶与PPAR配体的生物活性，从而间接调节脂肪酸氧化和糖类摄取[5]，也可通过调节PPAR-α信号因子，从而刺激肌肉组织及肾脏组织的脂肪酸氧化而减少甘油三酯的分泌，进而改变生物体对胰岛素的敏感性，从而参与到胰岛素的分泌调节中[6-7]。2015年，孟宪然等[8]应用RNA-Seq技术对山羊的转录组进行高通量测序分析，结果表明ADIPOQ基因是研究家畜肉品质的重要候选基因。其它一些研究也表明，ADIPOQ基因对家畜的生产性状有一定影响，如ADIPOQ基因启动子区核苷酸c.-67G/A和c.-892C/T变异对猪胴体性状及肉质有显著影响[9]，内含子g.1735A/G变异对猪肩部脂肪量有显著影响[10]；在牛ADIPOQ基因启动子区中，核苷酸c.-176A/G变异对胴体性状有一定影响，c.-199C/T、c.-34G/A变异对牛眼肌面积和背膘厚度也具有一定影响[11-12]。【本研究切入点】近年来，关于羊ADIPOQ基因多态性变异也有报道[13]，但相关变异对绵羊生长性状的性别差异和对公羔胴体性状的分子遗传机制研究却未见报道。【拟解决的关键问题】本研究以ADIPOQ基因作为不同性别罗姆尼羊生长性状、公羔胴体性状的候选基因，检测不同绵羊群体ADIPOQ基因相应区域的遗传变异，结合不同性别罗姆尼羊的生长性状及屠宰公羔的胴体性状，分析罗姆尼羊ADIPOQ基因变异对生长性状的性别差异和对公羔胴体性状的影响。

1 材料与方法

研究于2015年4—8月在新西兰林肯大学动物基因标记实验室完成。

1.1 试验样品表型

用于检测基因多态性的316份样本主要来自于新西兰南岛26个不同农场的8个绵羊品种（美利奴羊，萨福克羊，特克赛尔羊，考利代羊，罗姆尼羊，派伦代羊，丘陵陶赛特羊和杜泊羊），饲养过程严格按照不同品种分地域饲养，以保证品种间无基因交流。用于表型性状分析的试验样本均采自新西兰南岛单个牧场的17只优秀种公羊的子一代，共1 185只罗姆尼羊。羔羊在出生当天佩戴耳标，记录耳号、性别、出生日期、初生重和出生等级（单羔、双羔或三羔）。3周龄左右对羔羊进行断尾并测定断尾重，同时用FTA卡收集羔羊血样。3月龄断奶，记录并测定断奶重，再结合初生重计算所有羔羊断奶前生长速度。样本在后期的胴体生长性状测定中，由于不可抗拒因素（如羔羊死亡，耳标脱落等因素），部分羔羊缺乏相关表型性状数据，因此，用于结果分析的羔羊数量与最初血液采集总数量有一定偏差。

屠宰公羔的热胴体重及V-GR在屠宰厂直接测定，其它胴体性状在屠宰公羔时采用澳大利亚肉品和畜牧业有限公司开发的视频图像分析系统（VIASCAN○®Sastek）技术直接估算相关胴体肌肉性状（肩部瘦肉量、腰部瘦肉量、后腿瘦肉量、总瘦肉量（肩部、腰部、后腿瘦肉量总和）、肩部、腰部及后腿瘦肉比例）。

1.2 试验方法

1.2.1 基因组DNA提取羔羊用FTA卡（Whatman, Middlesex, UK）收集血液，父本公羊从颈静脉采血后滴于FTA卡，常温保存。采用“两步法”提取基因组DNA[14]。

1.2.2 引物设计及PCR扩增根据GenBank公布的绵羊ADIPOQ基因序列（Genbank Number：NC_019458.1），应用Primer 5.0在线软件设计特异性引物P1和P2（表1），分别用于扩增ADIPOQ基因Exon-1和Exon-2区序列，引物由Integrated DNA Technologies（Coralville, IA, USA）合成。

PCR反应体系：总体积20 μL，包括FTA卡1.2mm血样圆片1个，10×PCR Buffer缓冲液2 μL，5×Q Buffer缓冲液2 μL，3 μmol·L-1MgCl2溶液1.2 μL，150 μmol·L-1dNTPs 1.2 μL，0.25 μmol·L-1上下游混合引物1 μL，5U·μL-1Taq聚合酶0.1 μL和ddH2O加至20 μL。

PCR扩增程序：94˚C预变性2 min，94˚C变性30 s，58˚C（P1引物）、62˚C（P2引物）退火30 s，72˚C延伸30 s，37个循环，延伸7 min，4˚C保存。PCR产物用1%琼脂糖凝胶电泳检测。

1.2.3 SSCP检测取琼脂糖检测过的PCR扩增产物2 μL于灭菌PCR离心管中，加入8 μL变性上样缓冲液（98%去离子甲酰胺、10 mmol·L-1EDTA、0.025%二甲苯氰、0.025%溴酚蓝），经金属浴105˚C热变性5 min后立即置于冰水混合物中冷却5 min，上样于14%非变性聚丙烯酰胺凝胶中（Acr:Bis=37.5﹕1），0.5×TBE缓冲液中电泳。电泳结束后进行银染显色，判断相应SSCP带型。P1和P2引物最佳SSCP检测条件见表2。

表 1 绵羊ADIPOQ基因的扩增引物信息

表2 SSCP电泳条件

1.3 数据处理

应用MEGA5软件比对等位基因序列；Popgen32.0软件计算等位基因频率、有效等位基因数（Ne）、纯合度（Ho）、杂合度（He）、并进行c2检验；PIC软件计算多态信息含量（）；SHEsis软件分析两个区域之间SNPs的连锁平衡性；

本研究利用Minitab（Version 16，Minitab Inc.，Pennsylavania）一般线性模型评估特定等位基因的存在/缺失和基因型对绵羊生长速度、胴体生产性能的影响。对单性状存在（1）/缺失（0）分析模型中，等位基因/基因型、父本、初生重/出生等级为固定因素，遵从下列最小二乘法方差模型：

Yijknm= μ+Gi+Mk+Fm+Xn+eijknm

其中，Y为表型，μ为群体平均值，Gi为父本，Mk为初生重/出生等级（取决于哪个因素对结果更具影响力），Fm为等位基因/基因型，Xn为因素互作效应，eijknm为随机误差。

在单个等位基因存在/缺失分析模型中，仅对群体频率大于5%的等位基因进行关联性分析。同时需对结果＜0.2的等位基因进行互作校正。分析基因型对生产性状的影响时，仅对群体频率大于10%的基因型进行关联性分析。所有分析值均采用“平均值±标准误”表示，＜0.05为显著水平，＜0.2为有影响趋势，＞0.2为无影响。

2 结果

2.1 PCR扩增及产物SSCP检测

P1和P2引物扩增绵羊ADIPOQ基因Exon-1和Exon-2区域，扩增产物经琼脂糖电泳检测与目的片段预测长度一致，并无特异性条带。两对引物PCR扩增产物SSCP检测结果见图1和图2。Exon-1区检测到A，B，C和D等4种等位基因。Exon-2区检测到A，B，C和D等4种等位基因。

图1 绵羊ADIPOQ基因Exon-1区SSCP 检测

图2 绵羊ADIPOQ基因Exon-2区SSCP 检测

表3 绵羊ADIPOQ基因Exon-1区和Exon-2区突变位点

#非同义突变 Non-synonymous

2.2 绵羊ADIPOQ基因等位基因序列比对

ADIPOQ基因起始密码子位于Exon-2区，因此Exon-1区序列不参与氨基酸编码。绵羊ADIPOQ基因不同区域序列比对结果见表3，Exon-1区共检测到7个SNPs突变，Exon-2区检测到6个SNPs突变，其中Exon-2区的c.46 T/C为非同义突变，导致第16位的中性络氨酸（Tyr）和碱性组氨酸（His）变异。

2.3 绵羊ADIPOQ基因多态性分析

2.3.1 绵羊ADIPOQ基因群体遗传学分析 Exon-1区，等位基因A和B在所有品种中均为优势等位基因，整体频率分别为55.70%和42.07%，除了在考利代羊中，A和B在的频率分别为39.02%和57.32%。在Exon-1区的4个等位基因中，C仅在罗姆尼羊中检测到，频率为4.93%，而D在罗姆尼、美利奴、派伦代和考利代羊中检测到，其余品种中没有检测到该等位基因，但其频率极低，频率分别为1.47%、0.70%、3.66%和3.56%，结果见表4。

在Exon-2区，A为优势等位基因且在所有品种中均有检测，基因频率为72.78%，但在特克赛尔和丘陵陶塞特羊中只检测到等位基因A。等位基因B是基因频率最低的，只有1.74%，且在罗姆尼、丘陵陶塞特、特克赛尔和派伦代羊中没有检测到。C在除派伦代、特克塞尔、萨福克和丘陵陶塞特羊的其它品种中均能检测到，基因频率为8.70%。D是仅次于A的等位基因，基因频率为16.77%，但在特克塞尔和丘陵陶塞特羊中没有检测到该等位基因，结果见表4。

绵羊ADIPOQ基因区域的遗传特性及Hardy- Weinberg平衡检验结果见表5。在Exon-1区，杜泊羊和丘陵陶塞特羊有较高的纯合度；在Exon-2区，仅有派伦代羊有较高的纯和度。比较不同群体间的多态性信息含量（），大部分品种在两段区域中均为中度多态（0.25＜＜0.5），但美利奴羊在Exon-2区呈高度多态（＞0.5），特克塞尔、派伦代和丘陵陶塞特在Exon-2区呈低度多态（＜0.25）。在所检测的8个品种中，两段多态区域的检测等位基因数与有效等位基因数均不相同，均为观测等位基因数大于有效等位基因数。Hardy-Weinberg平衡检验结果表明，在所检测的两段区域中，大多数品种均处于平衡状态（＞0.05），只有罗姆尼羊在Exon-1区＜0.05、美利奴羊和杜泊羊在Exon-2区＜0.01，偏离了Hardy-Weinberg平衡。

表4 ADIPOQ基因不同区域在各品种中的等位基因频率

表5 不同绵羊ADIPOQ基因遗传特性及Hardy-Weinberg平衡检验

粗体表示有显著相关性（＜0.05） Statistically significant in bold (＜0.05)

2.3.2 绵羊ADIPOQ基因不同区域间连锁平衡分析常用于衡量连锁不平衡的两个参数为D’和r，其中区域内重组事件发生连锁不平衡的概率通过D’值反映，而连锁分析通过r值反映。SLATKIN等[15]研究分析认为D’值大于0.8则表明位点间强连锁且不平衡，ARDLIE等[16]分析认为r值大于0.33表明不同位点是紧密连锁且趋向于一个整体遗传。本研究中Exon-1区域中A、B的与Exon-2区域中的A、D最为普遍，基因频率最高，因此，将其做连锁不平衡分析发现，绵羊ADIPOQ基因不同变异区域之间D’为0.952，大于0.8，且r值为0.365，大于0.33，因此这两段区域处于强连锁不平衡状态，表明这两段区域的变异为高度连锁，趋向于整体遗传。

2.3.3 绵羊ADIPOQ基因变异与生长及胴体性状相关性分析在1 185只羔羊中，在Exon-1区域中检测到4个等位基因中等位基因C和D在新西兰罗姆尼羊中的频率小于5%，不能用于关联分析，因此，只对等位基因A和B进行关联分析。

对具有完整生长数据特征的550只公羔和537只母羔进行等位基因存在与缺失分析，结果表明，公羔中携带等位基因A的群体具有较低的断尾重（=0.019）、断奶重（=0.017）和断奶前生长速度（=0.025），等位基因B与公羔生长性状无显著相关。在母羔中等位基因A与生长性状无显著相关，但发现携带等位基因B的群体具有较高的断尾重（=0.019）；模型多因素矫正以上结果同样发现具有显著相关（＜0.05），结果见表6。

基因型与不同性别绵羊生长性状相关性分析表明，不同基因型与公羔的断尾重和断奶重有显著相关（＜0.05），携带BB基因型的群体较其它基因型群体有更高的断尾重和断奶重；但在母羔中未发现基因型与生长性状有显著相关，结果见表7。

关联分析时将所有公羔和母羔置于同一混合群体模型中，将性别作为影响因素进行GLMs模型分析矫正，其它参数相同，分析结果表明混合群体模型中等位基因和基因型均对羔羊生长性状无显著影响，结果未列出。

表6 ADIPOQ基因Exon-1区等位基因与罗姆尼公/母羔羊生长性状关联性分析

数值为平均值±标准误；粗体表示有显著相关性（＜0.05），斜体表示有影响趋势（0.05＜＜0.2）；n 表示样本数量

Estimates are given as mean±SE; statistically significant in bold (＜0.05) and trends in italics (0.05＜＜0.2); n: the number of observed sheep

表7 ADIPOQ基因Exon-1区基因型与罗姆尼公/母羔羊生长性状关联性分析

粗体表示相关性显著（＜0.05），斜体表示有影响趋势（0.05＜＜0.2）Statistically significant in bold (＜0.05) and trends in italics (0.05＜＜0.2)

对屠宰的456只公羔进行胴体性状相关性分析，结果表明等位基因A与热胴体重、腰部瘦肉量、后腿瘦肉量和总瘦肉量均有显著相关（＜0.05），携带等位基因A的群体具有较低的热胴体重、腰部瘦肉量、后腿瘦肉量和总瘦肉量；等位基因B与肩部瘦肉比例、后腿瘦肉量和后腿瘦肉比例有显著相关（＜0.05）；携带等位基因B的群体具有较低的肩部瘦肉比例和较高的后腿瘦肉量、后腿瘦肉比例；多因素矫正以上结果同样具有显著相关（＜0.05）。其它胴体性状与这两个等位基因均没有显著相关，结果见表8。

基因型与胴体性状相关分析表明，基因型与热胴体重、腰部瘦肉量、后腿瘦肉量和总瘦肉量均有显著相关（＜0.05），与其它胴体性状无显著相关（＞0.05）；其中携带基因型BB的个体均有较高的热胴体重、腰部瘦肉量、后腿瘦肉量和总瘦肉量，结果见表9。

3 讨论

本研究所检测的绵羊包括肉用、毛用和肉毛兼用的8种来自新西兰不同农场的商业品种绵羊，共有13个核苷酸变异位点被发现，该结果表明绵羊ADIPOQ基因具有丰富的多态性。已有研究发现在人类[17-18]、猪[19]、牛[11]及山羊[20]中也发现该基因具有一定的核苷酸变异，进一步验证了绵羊ADIPOQ基因的遗传多样性。在发现的13个突变位点中，Exon-2区的c.46 T/C突变为非同义突变，导致第16位的中性络氨酸（Tyr）与碱性组氨酸（His）的变化。该核苷酸突变恰位于ADIPOQ蛋白信号肽区域且临近于信号肽与易变区的分界处，推测该氨基酸的突变可能导致ADIPOQ蛋白空间结构的改变，从而影响信号肽的相关功能，诸如可能影响ADIPOQ蛋白的转录后修饰信号位点的结合、信号肽工作的水准及进程。已有研究表明，在人类编码区相近位置，rs2241766 G/T的变异导致氨基酸变异与多个人群ADIPOQ基因表达水平具有关联性[21-22]，因此本研究发现的氨基酸突变值得进一步研究其具体生理功能，如其是否对绵羊生长、胴体组成及脂肪含量存在影响。

本研究所检测的等位基因中只有等位基因A和B较为普遍，等位基因C和D的数量较少，这可能由于本研究所选的罗姆尼羊均为商品羊，经过了牧民长期的自然选择，含有等位基因C和D的群体因为某些生产性能不符合当地的需求已被筛选淘汰。分析等位基因A和B对绵羊生长性状的影响发现，其对绵羊的生长性状存在性别差异，等位基因A的存在/缺失对公羔的断尾重、断奶重和断奶前生长速度均有显著影响，而对母羔却无影响；而等位基因B的存在/缺失却只对母羔的断尾重有显著影响，对公羔却无影响。这表明ADIPOQ基因变异对绵羊生长性能的影响存在性别差异。相关研究表明，母牛胎盘中ADIPOQ蛋白的表达水平与牛犊的初生重具有显著影响[23]，在非洲裔美国人群体中也发现ADIPOQ基因变异对表达水平和肥胖存在性别差异，且该变异仅与非洲裔美国人群体中的妇女肥胖有显著关联[24]。目前亦有其它研究发现一些与肌肉生长调控相关基因中也存在一定的性别差异，REISZ-PORSZASZ等[25]发现MSTN基因在雄性小鼠中可能对其肌肉生长具有抑制作用，HAN等[26]发现MSTN基因变异在绵羊中存在性别比例差异。LEE等[27]发现WFIKKN2基因变异与小鼠的不同生长时期的肌肉比重存在性别关联性差异，而WANG等[28]发现WFIKKN2基因变异与绵羊生长性状存在性别关联性差异。这与本研究结果类似，同时在实际生产中一般都是屠宰优质公羔肉用，选留优质母羔用于繁殖，因此有必要进一步加大群体数量，研究该变异对生长性状性别差异的有效机制，进一步证明本研究结果，最终用于指导提高屠宰公羔的生产效益。

表8 ADIPOQ基因启动子变异体对罗姆尼羔羊胴体肌肉性状的影响

粗体表示相关性显著（＜0.05），斜体表示有影响趋势（0.05＜＜0.2）；后退、腰部、肩部和总瘦肉量=后腿、腰部、肩部、总瘦肉量与总胴体重的比值，后腿、腰部和肩部瘦肉比例=后腿、腰部和肩部瘦肉与相应部位总重量的比值。下同

Statistically significant in bold (＜0.05) and trends in italics (0.05＜＜0.2); Leg, loin, shoulder and total yield means the ratio of leg lean meat, loin lean meat, shoulder lean meat and total lean meat to total weight of hot carcass. Proportion of leg, loin and shoulder yield means the ratio of leg lean meat, loin lean meat and shoulder lean meat to total weight of corresponding parts. The same as below

表9 ADIPOQ基因启动子基因型对罗姆尼羔羊胴体肌肉性状的影响

ADIPOQ基因变异对罗姆尼羊的胴体性状有一定影响。本试验结果表明，等位基因A与腰部瘦肉量、后腿瘦肉量和总瘦肉量有一定相关，存在等位基因A的群体较缺失的群体平均低（0.40±0.19）%，（0.25±0.17）%和（0.84±0.16）%，这意味着对一个胴体重为18kg的羔羊来说，缺失较存在的个体相当于增加了151.2g左右的瘦肉，尽管看起来增加的量并不是很多，但对于养羊产业及现代居民对羊肉需求的营养水平来讲，这意味着更好的市场、更高的利润，尤其是腰部瘦肉量的提高，更具有价值，因为腰部瘦肉对羔羊胴体来讲是价值最高的部位。而等位基因B与肩部瘦肉比例和后腿瘦肉比例有一定关联性。进一步基因型分析结果还表明，不同基因型与热胴体重、腰部瘦肉量、后腿瘦肉量和总瘦肉量也相关，且本研究中发现的标记位点均与改善绵羊的瘦肉产量相关，而现代居民对肉制品的需求正在转变，越来越倾向于追求健康的饮食，过高脂肪含量的肉制品为消费者所不选择，缘由为过于肥胖的肉制品会导致肥胖、糖尿病或心脏类疾病[29]。本研究中的发现结果在绵羊上少有报道，但在其它动物上已有相关报道。如SHIN等[12]发现ADIPOQ基因突变与牛的背最长肌、脂肪厚度和眼肌面积等性状相关；CIESLAK等[9]发现ADIPOQ基因突变与猪的眼肌面积和腰部肌肉纤维长度等性状相关；OWECKI等[31]发现ADIPOQ基因突变可引起猪体内脂肪连接蛋白水平的改变。这进一步证明了本研究的结果，因此本研究中的相应等位基因和基因型可作为选育提高绵羊胴体性状的分子标记。

4 结论

本研究在8个绵羊群体的ADIPOQ基因中的2个不同区域共发现13个突变位点，其中Exon-2区的c.46 T/C突变为非同义突变，且这2个区域变异为高度连锁。罗姆尼羊Exon-1区等位基因与基因型均可影响绵羊生长性状和胴体性状，淘汰携带等位基因A的个体和选留存在等位基因B的个体、或选留BB基因型的个体，均可有效改善绵羊后代群体的部分生长性状和胴体性状。

[1] YANG Y, ZHANG F, DING R, SKRIP L, WANG Y, LEI H, HU D.gene polymorphisms and cancer risk: a meta-analysis., 2013, 61(2): 565-671.

[2] YANG W S, CHUANG L M. Human genetics of adiponectin in the metabolic syndrome., 2006, 84(2): 112-121.

[3] SCHERER P E, WILLIAMS S, FOGLIANO M, BALDINI G, LODISH H F. A Novel Serum Protein Similar to C1q, Produced Exclusively in Adipocytes., 1995, 270(45): 26746-26749.

[4] HSUEH W C, PL S J, MITCHELL B D, POLLIN T I, KNOWLER W C, EHM M G, BELL C J, SAKUL H, WAGNER M J, BURNS D K. Genome-wide and fine-mapping linkage studies of type 2 diabetes and glucose traits in the Old Order Amish: evidence for a new diabetes locus on chromosome 14q11 and confirmation of a locus on chromosome 1q21-q24., 2003, 52(2): 550-557.

[5] YAMAUCHI T, KAMON J, ITO Y, TSUCHIDA A, YOKOMIZO T, KITA S, SUGIYAMA T, MIYAGISHI M, HARA K, TSUNODA M. Cloning of adiponectin receptors that mediate antidiabetic metabolic effects., 2003, 423(6941): 762-769.

[6] KHARROUBI I, RASSCHAERT JEIZIRIK D L, CNOP M. Expression of adiponectin receptors in pancreatic beta cells., 2003, 312(4): 1118-1122.

[7] YAMAUCHI T, HARA K, KUBOTA N, TERAUCHI Y, TOBE K, FROGUEL P, NAGAI R, KADOWAKI T. Dual roles of adiponectin/Acrp30as an anti-diabetic and anti-atherogenic adipokine., 2003, 3(3): 243-253.

[8] 孟宪然，杜琛，王静，付绍印，郑竹清，张文广，李金泉. 基于RNA-Seq识别山羊肉品质候选基因. 畜牧兽医学报, 2015, 46(08): 1300-1307.

MENG X R, DU C, WANG J, FU S Y, ZHENG Z Q, ZHANG W G, LI J Q. RNA-Seq Approach for Identifying Candidate Genes of Meat Quality in Goats., 2015, 46(08): 1300-1307. (in Chinese)

[9] CIESLAK J, FLISIKOWSKA T, SCHNIEKE A, KIND A, SZYDLOWSKI M, SWITONSKI M, FLISIKOWSKI K. Polymorphisms in the promoter region of the adiponectin () gene are presumably associated with transcription level and carcass traits in pigs., 2013, 44(3): 340-343.

[10] DAI L H, XIONG Y Z, DENG C Y, JIANG S W, ZUO B, ZHENG R, LI F E, LEI M. Association of the A-G polymorphism in porcine adiponectin gene with fat deposition and carcass traits., 2006, 19(6): 779-783.

[11] MORSCI N S, SCHNABEL R D, TAYLOR J F. Association analysis of adiponectin and somatostatin polymorphisms on BTA1 with growth and carcass traits in Angus cattle., 2006, 37(6): 554–562.

[12] SHIN S, CHUNG E. Novel SNPs in the bovineandgenes are associated with carcass traits in Hanwoo (Korean cattle)., 2013, 40(7): 4651-4660.

[13] 刘重旭, 王凭青, 张宝云, 储明星, 邓腊梅, 谭颖, 樊奇. 贵州白山羊和古蔺马羊脂联素基因多态性及其与繁殖性能的关联研究. 中国农业科学, 2011, 44(9): 1916-1922.

LIU C X, WANG P Q, ZHANG B Y, CHU M X, DENG L M, TAN Y, FAN Q. Polymorphism of adiponectin gene and its relationship with reproductive ability in Guizhou White and Gulin Ma Goats., 2011, 44(9): 1916-1922. (in Chinese)

[14] ZHOU H, HICKFORD J. G. H. , FANG Q. A two-step procedure for extracting genomic DNA from dried blood spots on filter paper for polymerase chain reaction amplification., 2006, 354(1): 159-161.

[15] SLATKIN M. Linkage disequilibrium--understanding the evolutionary past and mapping the medical future., 2008，9(6): 477-485.

[16] ARDLIE K G, KRUGLYAK L, SEIELSTAD M. Patterns of linkage disequilibrium in the human genome., 2002, 3(4): 299-309.

[17] CHU H, WANG M, ZHONG D, SHI D, MA L, TONG N, ZHANG Z.polymorphisms are associated with type 2 diabetes mellitus: A meta‐analysis study., 2013, 29(7): 532-545.

[18] YUAN Y, JIANG H, KUANG J, HOU X, FENG Y, SU Z. Genetic variations ingene are associated with chronic obstructive pulmonary disease., 2012, 7(11): e50848.

[19] HOUDE A, MURPHY B, MATHIEU O, BORDIGNON V, PALIN M F. Characterization of swine adiponectin and adiponectin receptor polymorphisms and their association with reproductive traits., 2010, 39(3): 249-257.

[20] FANG X, DU Y, ZHANG C, SHI X Y, CHEN D X, SUN J J, JIN Q J, LAN X Y, CHENG H. Polymorphism in a microsatellite of the acrp30 gene and its association with growth traits in goats., 2011, 49(7-8): 533-539.

[21] YANG H, YE E, SI G, CHENG L M, CAI L Q, YE C F, ZHANG C, LU X M. Adiponectin gene polymorphism rs2241766 T/G is associated with response to pioglitazone treatment in type 2 diabetic patients from Southern China., 2014, 9(11): e112480

[22] GU H F. Biomarkers of adiponectin: plasma protein variation and genomic DNA polymorphisms., 2009, 4(4): 123-133.

[23] 沈留红, 江涛, 巫晓峰, 姜思汛, 肖劲邦, 曹随忠, 余树民, 邓俊良, 左之才, 彭广能, 马晓平, 钟志军, 任志华, 王娅, 胡延春. 奶牛胎盘脂联素、瘦素、内脂素与犊牛初生重相关性研究. 畜牧兽医学报, 2017, 48(01): 185-192.

SHEN L H, JIANG T, WU X F, JIANG S X, XIAO J B, CAO S Z, YU S M, DENG J L, ZUO Z C, PENG G N, MA X P, ZHONG Z J, REN Z H, WANG Y, HU Y C. The Correlation between adiponectin, leptin, visfatin in placenta and calf birth weight., 2017, 48(01): 185-192. (in Chinese)

[24] RIESTRA P, GEBREAB S Y, XU R, KHAN R J, BIDULESCU A, CORREA A, TEKOLA-AYELE F, DAVIS S K. Gender-specific associations betweengene polymorphisms and adiponectin levels and obesity in the Jackson Heart Study cohort., 2015, 16(1): 65.

[25] REISZPORSZASZ S, BHASIN S, ARTAZA J N, SHEN R , SINHA- HIKIM I, HOUGUE A, FIELDER T J, GONZALEA- CADAVID N F. Lower skeletal muscle mass in male transgenic mice with muscle- specific overexpression of myostatin., 2003, 285(4): E876-888.

[26] HAN J, ZHOU H, FORREST R H, SEDCOLE J R, FRAMPTON C M, HICKFORD J G H. Effect of myostatin (MSTN) g+6223G ＞ A on production and carcass traits in New Zealand Romney Sheep., 2010, 23(7): 863-866.

[27] LEE Y S, LEE S J. Regulation of GDF-11 and myostatin activity by GASP-1 and GASP-2., 2013, 110(39): E3713-3722.

[28] WANG J Q, ZHOU H, FANG Q, LIU X, LUO Y Z, HICKFORD J G H. Effect of variation in ovine WFIKKN2 on growth traits appears to be gender-dependent., 2015, 5: 12347; doi: 10. 1038/ srep12347.

[29] VOLK M G. An examination of the evidence supporting the association of dietary cholesterol and saturated fats with serum cholesterol and development of coronary heart disease., 2007, 12(3): 228-246.

[30] OWECKI M, MICZKE A, KACZMAREK M, HOPPE-GOLEBIEWSKA J, PUPEK-MUSIALIK D, SLOMSKI R, BRYLL W, CYMERYS M, NIKISCH E, SOWIŃSKI J. The Y111 H (T415C) polymorphism in exon 3 of the gene encoding adiponectin is uncommon in Polish obese patients., 2007, 39(11): 797-800.

（责任编辑林鉴非）

Polymorphisms ofGene and Their Association with Growth and Carcass Traits in Sheep

AN QingMing1, 2, ZHOU HuiTong2, 3, WU ZhenYang2, LUO YuZhu2, Jon G. Hickford2, 3

(1Tongren University, Tongren 554300, Guizhou, China;2Gansu Key Laboratory of Herbivorous Animal Biotechnology, Gansu Agricultural University, Lanzhou 730070,Gansu, China;3Gene-Marker Laboratory, Lincoln University, Lincoln 7647, New Zealand)

【Objective】In this study, polymorphisms and linkage relationship of ovine(adiponectin) gene were investigated and their effects on some growth and carcass traits were estimated, so as to enrich the molecular genetic data for sheep. 【Method】Mutations in Exon-1 and Exon-2 ofgene were detected by PCR-SSCP in 8 commodity sheep populations, and the relationship between mutations in growth and carcass traits in New Zealand (NZ) Romney lambs was investigated using General Linear Models (GLMs). 【Result】In total, thirteen SNPs were detected in Exon-1 and Exon-2 regions of ovinegene, and the nucleotide substitution c.46T/C in Exon-2 resulted in amino acid change (p.Tyr16His). AlleleAandBwere the dominant allele in Exon-1, AlleleAandDwere the dominant allele in Exon-2, and there were difference of allele frequencies between these two regions. The majority population was moderately polymorphic in all regions (＜0.5), except Texel, Perendale and Dorset Down was low polymorphic in Exon-2 region (＜0.25), and there was a high linkage relationship of these mutation and tend to shared genetic linkages (D’=0.952,r=0.365). The association analysis showed that the mutations in Exon-1 region of ovinegene had different effects on growth traits of male and female lambs. In male lambs, individuals with possessing allele Ahad lower tailing weight, weaning weight and pre-weaning growth rate than those no-possessing (＜0.05), but no associations were detected in female lambs (＞0.05). In female lambs, individuals with possessing allele Bhad higher tailing weight than those no-possessing (＜0.05), however, no associations were detected in male lambs (＞0.05). And individuals with possessingBBhad higher tailing weight and weaning weight in male lambs. The carcass traits association analysis results showed that individuals with possessing alleleAhad lower hot-carcass weight, loin yield, leg yield and total yield than those no-possessing (＜0.05), and individuals with possessing alleleBhad higher leg yield, proportion of led yield and lower proportion shoulder yield than those no-possessing (＜0.05), individuals with possessingBBhad higher hot carcass weight, loin yield, leg yield and total yield (＜0.05). 【Conclusion】Exon-1 and Exon-2 of ovinegene had abundant polymorphisms, and SNP c.46T/C was non-synonymous. The mutations of Exon-1 of ovinegene might affect some growth traits and carcass traits, and selecting sheep with alleleBand genotypeBB, or eliminating ones with alleleAand AAcould improve some growth traits and carcass traits of Romney sheep.