我国是世界纺织业发达国家之一，尤其对羊毛加工的需求量很大，而细羊毛是主要的精纺原料，尽管其他天然纤维以及随着工业发展涌现出的许多新型合成的纤维分享了细羊毛市场，但细羊毛因其绝缘性、吸湿性、弹性、抗静电、抗皱性等方面的特色优势，始终作为理想精纺加工的首选原料，占纺织原料的17%。

目前，角蛋白和角蛋白相关蛋白是模式生物中主要的结构蛋白，它确定了绒毛的基本性质，如细度、长度、强度、颜色和光泽等。在毛发的形成过程中，一系列的角蛋白和角蛋白关联蛋白还决定了毛囊细胞的分化过程。一个显著的特征是角蛋白基因按家系聚类在染色体结构中并紧密连接在一起，这种紧密的遗传连锁可能是共同作用调控它们表达的结果。羊毛生产性状在品种间、品种内和不同品系间的差异明显受到遗传的影响。本文通过分析KAP16基因的遗传多态性和遗传效应，旨在获取相关毛性状的DNA标记，为细毛羊的育种在分子育种学上提供理论依据。为细毛羊的育种在分子育种学上提供理论依据。冯静、艾买提和玛依拉等人分别研究了绵羊KAP基因和KRT基因多态性与毛性状的相关性，赵志东、于辉和赵苗等对}P基因遗传变异与基因性状及其对羊绒性状的影响都有深人研究。获取与绵羊相关毛性状的DNA分子标记，并为细毛羊的育种在分子育种学上提供理论依据。羊毛纤维具有高度组织结构，其中毛纤维的90%是由角蛋白中间丝(KRT-RF)和角蛋白关联蛋白(KAPs)构成，这两种蛋白类型通过高度组织分化的方式形成规则的晶状结构，其中KIF蛋白形成纤维丝镶嵌在KAP蛋白基质中，这种规则的晶状结构需要高度协调的遗传调控，然而，这种调控机理还不完全清楚。

1 材料与方法

1.1 试验材料

试验动物为新疆巩乃斯种羊场提供的中国美利奴羊(新疆型)，采取416只周岁母羊的毛样与血样。其中鉴定时体重、毛弯曲评分、油汗、体格大小、剪毛后体重、剪毛量等数据由新疆巩乃斯种羊场提供。毛样在羊毛分析实验室检测羊毛的纤维直径、纤维直径变异系数、毛长、细度、弯曲和弯曲数，通过整理共得到832条有效数据。静脉采血3一5 mL,用EDTA抗凝并冷冻保存备用，用于提取DNA，并筛选有羊毛性状鉴定记录的个体样本进行试验。

1.2 试验方法

1.2.1重测序技术及引物的合成取所测纤维细度数值中的两端极值各15个样本，建立2个DNA池(粗池和细池)，由诺禾致源生物公司运用Illuni-ma HiSeqTM2000/MiseqTM软件进行质控分析。引物来自NCBI提供的绵羊基因序列并结合重测序结果，采用Primer Premier 5.0软件设计，由上海生工生物工程有限公司合成，引物序列如表to

1.2.2   PCR扩增程序的设定见表2}

1.2.3   PCR产物测序经过普通PC R打一增和2 %琼脂糖凝胶电泳检测将合格的PCR产物送到上海生工生物工程有限公司进行测序。

1.3数据统计与分析

对数据进行zX独立性检验，多态信息含量，遗传纯合度(Ho)，杂合度(He)和有效等位基因数(Ne)的计算

1.4关联分析模型

通过整理共得到有效数据832条，并对毛性状进行遗传效应分析。运用SAS 8.1统计软件检测不同基因型与毛性状的关联性，并估计最小二乘均值。最小二乘方差分析数学模型:Y_ijk = μ + α_i + β_j + e_ijk ,式中，:Y_ijk:个体表型记录; μ:群体性状均值; α_i:标记基因型效应或标记组合基因型效应; β_j:群别效应，e_ijk:随机误差。

2结果与分析

2.1   KAPl6的重测序结果

由表3可知，KAPl6基因的C41006938T和T41007938C突变位点分别位于11号染色体的第5和4外显子上。

2.2  DNA完整性检测结果

提取基因组DNA模板用1%琼脂糖凝胶电泳检测(图1)，实验提取的基因组DNA模板质量良好、纯度高、不需要纯化，可直接用作模板DNA2.3  KAP16基因的PC R扩增结果

如图2-3所示，PC R产物经2.0%的琼脂糖凝胶电泳检测，可知KAPl6基因P1,P2引物的PCR产物长度为216 ,196娜与生物公司提供片段符合。结果表明，KAPl6打一增特异性良好，片段长度一致，引物的片段长度符合预期大小，可直接用于PCR-SSCP的分析。

2.4  KAPl6基因的PCR-SSCP多态性检测结果

如图4 }5所示，P1位点存在3种不同的带型，将其命名为AA,BB和AB型;P2位点存在，2种基因型，命名为AA,AB型。

2.5  KAP16基因位点多序列比

2.5.1   P1位点多序列比对比

由图6可知，在中国美利奴羊(新疆型)群体中发现此序列在第5外显子编码区第31娜发生突变。通过对氨基酸序列的分析发现，其C31T突变点为非同义突变，由Arg(精氨酸)突变为Trp(色氨酸)。

2.5.2  P2位点多序列比对结果及峰图

由图7可知，在中国美利奴羊(新疆型)群体中发现该序列在第四外显子处存在编码区第38娜发生了突变。突变位点为T38C，通过对氨基酸序列的分析发现，其突变点由Val ( _氨酸)突变为Ala丙氨酸)。

2.6  KAPl6的遗传多态性分析

由表4可知，该群体中C41006938T突变点的AB基因型为优势基因型,B基因是优势等位基因;T41007938C基因的AA基因型是优势基因型,A基因是优势等位基因。卡方检验表明，C41006938T和T41007938C位点均处Hardy-Weinberg平衡状态(P> 0.OS )。中国美利奴(新疆型)群体中C41006938T和T41007938C突变点的PIC为0.37 ,0.18均属于中度多态，因此中国美利奴(新疆型)群体的遗传变异居于中等水平。

2.7  KAPl6基因位点与羊毛性状的关联性分析

2.7.1   KAPl6基因位点不同基因型的最小二乘方差分析

由表5可知，C41006938T不同基因型对体格大小和鉴定时体重有显著影响(P <0.05)。群别对鉴定时体重、纤维细度、体格大小和纤维直径有极显著影响(P < 0.O1 ) o T41007938C不同基因型对中国美利奴(新疆型)的纤维直径变异系数有极显著影响(P<0.01);对剪毛量和弯曲有显著影响(P<0.05);不同群别对中国美利奴(新疆型)鉴定是体重、纤维细度、体格大小和纤维直径有极显著影响(P<0.01)。

2.7.2KAP16基因位点不同基因型与毛性状的关联性分析由表6可见，C41006938T位点BB基因型对AA和AB基因型鉴定时体重和体格大小差异显著(P < 0.OS ) ; 741007938 C位点AB基因型对AA基因型纤维直径变异系数差异极显著(P <0.01),对剪毛量和弯曲差异显著(P <0.05) o

2.8   KAPl6基因位点不同基因型组合效应与羊毛性状的关联性

2.8.1   KAPl6不同基因型组合效应最小二乘方差分析由表7可知，C41006938T和T41007938C对剪毛量有极显著影响(P<0.01);对纤维直径变异系数有显著影响(P <0.05);不同群别对鉴定时体重、纤维细度、体格大小、剪毛后体重、纤维直径有极显著影响(P<0.01)0

2.8.2   KAPl6不同基因型组合效应

由表8可知，KAPl‘组合基因AABB型对剪毛量差异极显著(P< 0.O1 ) ; AABB基因型,ABAA基因型和ABBB基因型对纤维直径变异系数差异显著(P<0.05)0

3讨论

3.1  KAP基因的多态性分析

Mclaren等通过连锁分析已经将KAPl.1、KAPl.2 、 KAPl.3 , KAP3.2定位到绵羊的21号染色体上。研究人员针对KAP基因家族的多态性做了许多研究，Parsons等以KRTAP6和KRTAP8为候选基因，证实这两个基因位点和羊毛产量及纤维直径存在明显的相关性，推测在KAP6和KAP8同一区域可能存在控制羊毛直径的主基因。Henry等利用超细的美利奴羊和罗姆尼羊的回交群体，测定了包括周岁体重、羊毛产量、平均纤维直径在内的多个性状，同时用216个微卫星标记进行了基因组扫描，找到了一个与羊毛纤维直径相关的位点。Itenge-MWPZar在澳洲美利奴羊中报道KAPl.3(B2C)基因存在9个独特的SSCP带型，序列分析显示这9个带型代表8个新的KAPl.3等位基因，在KAPl.3基因c.2406 > A突变位点的3种基因型分别为GG ,AA和AG，在KAPl.3基因c.3126 > A突变位点的GG,AA和AG等3种基因型。KRT和KAP多态性位点可能导致基因部分改变的从而导致羊毛特征的变化，这些识别标记可用于在羊毛产业的基因分子标记辅助选择。刘桂芬选择21号染色体上高硫蛋白家族中KAPl.1、KAPl.3中的部分序列和1号染色体高甘氨酸一酪氨酸蛋白KAP6.1的外显子进行研究，结果发现KAPl.1、KAPl.3所在区域中位点W08667与羊毛细度显著相关。张亚妮应用PCR-SSCP技术分析了内蒙古白绒山羊和辽宁绒山羊KAP基因与经济性状的关系，KAPl.3基因S4位点AB基因型与辽宁绒山羊产绒量相关，KAPl.3基因S3位点AA基因型与内蒙古白绒山羊体重相关，KAPl.3基因S2位点BB基因型与辽宁绒山羊体重相关。冯静等应用PCR-RFLP方法分析绵羊KAPl.3基因遗传多态性，结果表明KAPl.3基因存在2个SNP位点。王春听应用PCR-SSCP技术对羊毛纤维组成蛋白中的KAPl.3基因部分序列进行多态性分析，结果表明KAPl.3基因与产毛量和拉伸长度相关。

3.2国内外对绵羊KAP基因家族的相关研究

目前，该领域的研究热点之一是关于羊毛细度决定位点或相关基因的研究Mclaren等通过连锁分析已经将角蛋白辅助蛋白KAP6.1 , KAP7和KAP8定位到绵羊的第1号染色体上。Ponz等通过片段定位方法寻找控制羊毛性状的QTL，利用法国的合成品系INR401对40个微卫星标记进行了分析，结果发现了在3号、7号、25号染色体上有控制毛长(SL)的QTL，在6号染色体上有控制毛纤维直径(WFD)的QTL，在7号染色体上有控制毛纤维直径纤维直径变异系数的QTL}  Cockett N E认为KAP6.n和羊毛直径有关。Chen等研究表明KAPB.1和KAPl.3基因表达模式相似，在羊皮肤中KAPl.3高度表达，KAPB.1表达较弱。此外，KAPB.1和KAPl.3的位点多态性与纤维直径差异显著。KAPB.1中BC基因型对纤维直径有极显著的影响，KAPl.3中AA基因型对纤维直径有极显著的影响。此外，KAPB.1与KAPl.3基因的最优组合基因型为BC-AA组合基因型。

3.3  KAPl6基因与中国美利奴羊(新疆型)毛性状的关连分析

羊毛性状是细毛羊重要的经济性状之一，分子标记成为研究重要表型的方法之一，并被广泛用于动物的遗传研究中。KAPl6基因通过全基因组测序攻击获得8个SNPs突变位点，本文研究了其中的两个突变点对中国美利奴(新疆型)的遗传效应分析，将继续研究其他位点对中国美利奴(新疆型)的遗传效应。

C41006938T不同基因型对体格大小和鉴定时体重有显著影响。群别对鉴定时体重、纤维细度、体格大小和纤维直径有极显著影响。C41006938T基因位点对鉴定时体重和体格大小差异显著;T41007938C不同基因型对中国美利奴(新疆型)的纤维直径变异系数有极显著影响;对剪毛量和弯曲有显著影响;不同群别对中国美利奴(新疆型)鉴定是体重、纤维细度、体格大小和纤维直径有极显著影响。T41007938C基因型对中国美利奴羊(新疆型)纤维直径变异系数差异极显著;对剪毛量和弯曲差异显著。C41006938T和T41007938C对剪毛量有极显著影响;对纤维直径变异系数有显著影响;不同群别对鉴定时体重、纤维细度、体格大小、剪毛后体重、纤维直径有极显著影响。KAP16组合基因型对中国美利奴(新疆型)剪毛量差异极显著;对纤维直径变异系数差异显著。

在所研究的绵羊群体中，KAPl6基因的2个突变位点均呈现出中度多态性，说明在这些位点处中国美利奴(新疆型)群体遗传变异处于中等水平，变异程度相似，没有太大的差异，同时也说明在上述各群体的培育改良中可以加强人工选择的强度。

4 结论

本试验首次通过重测序技术，发现KaPl6基因在染色体11上第4外显子SNP突变点为T41007938C，第5外显子SNP突变点为C41006938T。通过分析可知，C41006938T不同基因型对体格大小和鉴定时体重有显著影响;T41007938C不同基因型对纤维直径纤维直径变异系数有极显著影响，对剪毛量和弯曲有显著影响;C41006938T和T41007938C组合基因型对剪毛量有极显著影响，对纤维直径纤维直径变异系数有显著影响。