何軍敏，黃錫霞，田可川，趙冰茹，柏妍，田月珍，徐新明，付雪峰

(1.新疆農(nóng)業(yè)大學(xué)動物科學(xué)學(xué)院，新疆烏魯木齊830052;2.新疆畜牧科學(xué)院畜牧研究所，分子遺傳重點實驗室，新疆烏魯木齊830011)

KAP16基因?qū)?xì)毛羊重要經(jīng)濟(jì)性狀的遺傳效應(yīng)分析

何軍敏1，2，黃錫霞1*，田可川2*，趙冰茹1，柏妍2，田月珍1，徐新明2，付雪峰2

(1.新疆農(nóng)業(yè)大學(xué)動物科學(xué)學(xué)院，新疆烏魯木齊830052;2.新疆畜牧科學(xué)院畜牧研究所，分子遺傳重點實驗室，新疆烏魯木齊830011)

為了分析KAP16基因的遺傳多態(tài)性，獲取相關(guān)毛性狀的DNA標(biāo)記，為細(xì)毛羊在分子育種上提供依據(jù)，采用重測序技術(shù)和PCR-SSCP技術(shù)，驗證KAP16候選基因的2個SNPs在中國美利奴(新疆型)群體中的遺傳多態(tài)性。同時，運(yùn)用SAS 8.1軟件進(jìn)行最小二乘方差分析方法，研究其與毛性狀的關(guān)聯(lián)性。結(jié)果表明，C41006938T突變點得到3種基因型為AA、AB和BB型，該群體中C41006938T突變點基因型頻率為0.22、0.47和0.31;A和B等位基因頻率為0.46和0.54，其中B等位基因為優(yōu)勢等位基因。T41007938C基因位點AA和AB基因型頻率為0.77和0.23;A和B等位基因頻率分別為0.89和0.11，A基因為優(yōu)勢等位基因。χ2檢驗表明，C41006938T和T41007938C位點均處于哈代溫伯格平衡狀態(tài)(P＞0.05)。C41006938T和T41007938C突變位點的多態(tài)信息含量分別為0.37和0.18，所以中國美利奴(新疆型)群體在該基因的遺傳變異處在中等水平。C41006938T不同基因型對體格大小和鑒定時體重有顯著影響(P＜0.05);T41007938C不同基因型對纖維直徑纖維直徑變異系數(shù)有極顯著影響(P＜0.01)，對剪毛量和彎曲有顯著影響(P＜0.05);C41006938T和T41007938C組合基因型對剪毛量有極顯著影響(P＜0.01)，對纖維直徑纖維直徑變異系數(shù)有顯著影響(P＜0.05)。

KAP16基因;SNPs;細(xì)毛羊;關(guān)聯(lián)性分析

我國是世界紡織業(yè)發(fā)達(dá)國家之一，尤其對羊毛加工的需求量很大，而細(xì)羊毛是主要的精紡原料，盡管其他天然纖維以及隨著工業(yè)發(fā)展涌現(xiàn)出的許多新型合成的纖維分享了細(xì)羊毛市場，但細(xì)羊毛因其絕緣性、吸濕性、彈性、抗靜電、抗皺性等方面的特色優(yōu)勢，始終作為理想精紡加工的首選原料，占紡織原料的17%［1］。

目前，角蛋白和角蛋白相關(guān)蛋白是模式生物中主要的結(jié)構(gòu)蛋白，它確定了絨毛的基本性質(zhì)，如細(xì)度、長度、強(qiáng)度、顏色和光澤等。在毛發(fā)的形成過程中，一系列的角蛋白和角蛋白關(guān)聯(lián)蛋白還決定了毛囊細(xì)胞的分化過程［2］。一個顯著的特征是角蛋白基因按家系聚類在染色體結(jié)構(gòu)中并緊密連接在一起，這種緊密的遺傳連鎖可能是共同作用調(diào)控它們表達(dá)的結(jié)果。羊毛生產(chǎn)性狀在品種間、品種內(nèi)和不同品系間的差異明顯受到遺傳的影響［3］。本文通過分析KAP16基因的遺傳多態(tài)性和遺傳效應(yīng)，旨在獲取相關(guān)毛性狀的DNA標(biāo)記，為細(xì)毛羊的育種在分子育種學(xué)上提供理論依據(jù)。為細(xì)毛羊的育種在分子育種學(xué)上提供理論依據(jù)。馮靜、艾買提和瑪依拉等人分別研究了綿羊KAP基因和KRT基因多態(tài)性與毛性狀的相關(guān)性［4－8］，趙志東、于輝和趙苗等對KAP基因遺傳變異與基因性狀及其對羊絨性狀的影響都有深入研究［9－11］。獲取與綿羊相關(guān)毛性狀的DNA分子標(biāo)記，并為細(xì)毛羊的育種在分子育種學(xué)上提供理論依據(jù)。羊毛纖維具有高度組織結(jié)構(gòu)，其中毛纖維的90%是由角蛋白中間絲(KRT-RF)和角蛋白關(guān)聯(lián)蛋白(KAPs)構(gòu)成，這兩種蛋白類型通過高度組織分化的方式形成規(guī)則的晶狀結(jié)構(gòu)，其中KIF蛋白形成纖維絲鑲嵌在KAP蛋白基質(zhì)中，這種規(guī)則的晶狀結(jié)構(gòu)需要高度協(xié)調(diào)的遺傳調(diào)控，然而，這種調(diào)控機(jī)理還不完全清楚。

1 材料與方法

1.1 試驗材料

試驗動物為新疆鞏乃斯種羊場提供的中國美利奴羊(新疆型)，采取416只周歲母羊的毛樣與血樣。其中鑒定時體重、毛彎曲評分、油汗、體格大小、剪毛后體重、剪毛量等數(shù)據(jù)由新疆鞏乃斯種羊場提供。毛樣在羊毛分析實驗室檢測羊毛的纖維直徑、纖維直徑變異系數(shù)、毛長、細(xì)度、彎曲和彎曲數(shù)，通過整理共得到832條有效數(shù)據(jù)。靜脈采血3～5 mL，用EDTA抗凝并冷凍保存?zhèn)溆?，用于提取DNA，并篩選有羊毛性狀鑒定記錄的個體樣本進(jìn)行試驗。

1.2 試驗方法

1.2.1 重測序技術(shù)及引物的合成取所測纖維細(xì)度數(shù)值中的兩端極值各15個樣本，建立2個DNA池(粗池和細(xì)池)，由諾禾致源生物公司運(yùn)用Illunima HiSeqTM2000/MiseqTM軟件進(jìn)行質(zhì)控分析。引物來自NCBI提供的綿羊基因序列并結(jié)合重測序結(jié)果，采用Primer Premier 5.0軟件設(shè)計，由上海生工生物工程有限公司合成，引物序列如表1。

1.2.2 PCR擴(kuò)增程序的設(shè)定見表2。

1.2.3 PCR產(chǎn)物測序經(jīng)過普通PCR擴(kuò)增和2%瓊脂糖凝膠電泳檢測將合格的PCR產(chǎn)物送到上海生工生物工程有限公司進(jìn)行測序。

1.3 數(shù)據(jù)統(tǒng)計與分析

對數(shù)據(jù)進(jìn)行χ2獨立性檢驗，多態(tài)信息含量(PIC)，遺傳純合度(Ho)，雜合度(He)和有效等位基因數(shù)(Ne)的計算

1.4 關(guān)聯(lián)分析模型

通過整理共得到有效數(shù)據(jù)832條，并對毛性狀進(jìn)行遺傳效應(yīng)分析。運(yùn)用SAS 8.1統(tǒng)計軟件檢測不同基因型與毛性狀的關(guān)聯(lián)性，并估計最小二乘均值。最小二乘方差分析數(shù)學(xué)模型:Yijk=u+αi+βj+eijk，式中，Yijk:個體表型記錄;μ:群體性狀均值;αi:標(biāo)記基因型效應(yīng)或標(biāo)記組合基因型效應(yīng);βj:群別效應(yīng)，eijk:隨機(jī)誤差。

表1 PCR擴(kuò)增的引物序列Table 1 PCR amplification primer sequence

表2 PCR擴(kuò)增條件Table 2 PCR amplification conditions

表3 KAP16的重測序結(jié)果Table 3 KAP16 sequencing result

圖1 DNA的結(jié)果Fig.1 Results of DNA

圖2 P1的PCR產(chǎn)物Fig.2 PCR products of P1

圖3 P2的PCR產(chǎn)物Fig.3 PCR products of P2

2 結(jié)果與分析

2.1 KAP16的重測序結(jié)果

由表3可知，KAP16基因的C41006938T和T41007938C突變位點分別位于11號染色體的第5和4外顯子上。

2.2 DNA完整性檢測結(jié)果

圖4 P1的SSCP產(chǎn)物Fig.4 SSCP products of P1

圖5 P2的SSCP產(chǎn)物Fig.5 SSCP products of P2

提取基因組DNA模板用1%瓊脂糖凝膠電泳檢測(圖1)，實驗提取的基因組DNA模板質(zhì)量良好、純度高、不需要純化，可直接用作模板DNA。

2.3 KAP16基因的PCR擴(kuò)增結(jié)果

如圖2～3所示，PCR產(chǎn)物經(jīng)2.0%的瓊脂糖凝膠電泳檢測，可知KAP16基因P1、P2引物的PCR產(chǎn)物長度為216、196 bp與生物公司提供片段符合。結(jié)果表明，KAP16擴(kuò)增特異性良好，片段長度一致，引物的片段長度符合預(yù)期大小，可直接用于PCRSSCP的分析。

2.4 KAP16基因的PCR-SSCP多態(tài)性檢測結(jié)果

如圖4～5所示，P1位點存在3種不同的帶型，將其命名為AA、BB和AB型;P2位點存在，2種基因型，命名為AA、AB型。

2.5 KAP16基因位點多序列比

圖6 P1突變點不同基因型測序峰圖Fig.6 P1 mutation point of different genotype sequencing peak figure

圖7 P2突變點不同基因型測序峰圖Fig.7 P2 mutation point of different genotype sequencing peak figure

2.5.1 P1位點多序列比對比由圖6可知，在中國美利奴羊(新疆型)群體中發(fā)現(xiàn)此序列在第5外顯子編碼區(qū)第31 bp發(fā)生突變。通過對氨基酸序列的分析發(fā)現(xiàn)，其C31T突變點為非同義突變，由Arg (精氨酸)突變?yōu)門rp(色氨酸)。

2.5.2 P2位點多序列比對結(jié)果及峰圖由圖7可知，在中國美利奴羊(新疆型)群體中發(fā)現(xiàn)該序列在第四外顯子處存在編碼區(qū)第38 bp發(fā)生了突變。突變位點為T38C，通過對氨基酸序列的分析發(fā)現(xiàn)，其突變點由Val(纈氨酸)突變?yōu)锳la(丙氨酸)。

2.6 KAP16的遺傳多態(tài)性分析

由表4可知，該群體中C41006938T突變點的AB基因型為優(yōu)勢基因型，B基因是優(yōu)勢等位基因; T41007938C基因的AA基因型是優(yōu)勢基因型，A基因是優(yōu)勢等位基因?？ǚ綑z驗表明，C41006938T和T41007938C位點均處Hardy-Weinberg平衡狀態(tài)(P＞0.05)。中國美利奴(新疆型)群體中C41006938T和T41007938C突變點的PIC為0.37，0.18均屬于中度多態(tài)，因此中國美利奴(新疆型)群體的遺傳變異居于中等水平。

表4 基因頻率和基因型頻率的分布以及Hardy-Weinberg平衡的卡方檢驗Table 4 Distribution of gene frequency and genotype frequency and Hardy-Weinberg equilibrium card square test

表5 C41006938T和T41007938C基因型和群別對毛性狀的影響Table 5 The effect of C41006938T and T41007938C genotype and group on wool traits

表6 C41006938T和T41007938C不同基因型與毛性狀間的關(guān)系Table 6 C41006938T and T41007938C relationship between different genotypes and wool traits

2.7 KAP16基因位點與羊毛性狀的關(guān)聯(lián)性分析

2.7.1 KAP16基因位點不同基因型的最小二乘方差分析由表5可知，C41006938T不同基因型對體格大小和鑒定時體重有顯著影響(P＜0.05)。群別對鑒定時體重、纖維細(xì)度、體格大小和纖維直徑有極顯著影響(P＜0.01)。T41007938C不同基因型對中國美利奴(新疆型)的纖維直徑變異系數(shù)有極顯著影響(P＜0.01);對剪毛量和彎曲有顯著影響(P＜0.05);不同群別對中國美利奴(新疆型)鑒定是體重、纖維細(xì)度、體格大小和纖維直徑有極顯著影響(P＜0.01)。

2.7.2 KAP16基因位點不同基因型與毛性狀的關(guān)聯(lián)性分析由表6可見，C41006938T位點BB基因型對AA和AB基因型鑒定時體重和體格大小差異顯著(P＜0.05);T41007938C位點AB基因型對AA基因型纖維直徑變異系數(shù)差異極顯著(P＜0.01)，對剪毛量和彎曲差異顯著(P＜0.05)。

2.8 KAP16基因位點不同基因型組合效應(yīng)與羊毛性狀的關(guān)聯(lián)性

2.8.1 KAP16不同基因型組合效應(yīng)最小二乘方差分析由表7可知，C41006938T和T41007938C對剪毛量有極顯著影響(P＜0.01);對纖維直徑變異系數(shù)有顯著影響(P＜0.05);不同群別對鑒定時體重、纖維細(xì)度、體格大小、剪毛后體重、纖維直徑有極顯著影響(P＜0.01)。

2.8.2 KAP16不同基因型組合效應(yīng)由表8可知，KAP16組合基因AABB型對剪毛量差異極顯著(P＜0.01);AABB基因型、ABAA基因型和ABBB基因型對纖維直徑變異系數(shù)差異顯著(P＜0.05)。

3 討論

3.1 KAP基因的多態(tài)性分析

Mclaren等［12］通過連鎖分析已經(jīng)將KAP1.1、KAP1.2、KAP1.3、KAP3.2定位到綿羊的21號染色體上。研究人員針對KAP基因家族的多態(tài)性做了許多研究，Parsons等［13］以KRTAP6和KRTAP8為候選基因，證實這兩個基因位點和羊毛產(chǎn)量及纖維直徑存在明顯的相關(guān)性，推測在KAP6和KAP8同一區(qū)域可能存在控制羊毛直徑的主基因。Henry等［14］利用超細(xì)的美利奴羊和羅姆尼羊的回交群體，測定了包括周歲體重、羊毛產(chǎn)量、平均纖維直徑在內(nèi)的多個性狀，同時用216個微衛(wèi)星標(biāo)記進(jìn)行了基因組掃描，找到了一個與羊毛纖維直徑相關(guān)的位點。Itenge-Mweza［15］在澳洲美利奴羊中報道KAP1.3 (B2C)基因存在9個獨特的SSCP帶型，序列分析顯示這9個帶型代表8個新的KAP1.3等位基因，在KAP1.3基因c.240G＞A突變位點的3種基因型分別為GG、AA和AG，在KAP1.3基因c.312G＞A突變位點的GG、AA和AG等3種基因型。KRT和KAP多態(tài)性位點可能導(dǎo)致基因部分改變的從而導(dǎo)致羊毛特征的變化，這些識別標(biāo)記可用于在羊毛產(chǎn)業(yè)的基因分子標(biāo)記輔助選擇。劉桂芬［16］選擇21號染色體上高硫蛋白家族中KAP1.1、KAP1.3中的部分序列和1號染色體高甘氨酸-酪氨酸蛋白KAP6.1的外顯子進(jìn)行研究，結(jié)果發(fā)現(xiàn)KAP1.1、KAP1.3所在區(qū)域中位點W08667與羊毛細(xì)度顯著相關(guān)。張亞妮［17］應(yīng)用PCR-SSCP技術(shù)分析了內(nèi)蒙古白絨山羊和遼寧絨山羊KAP基因與經(jīng)濟(jì)性狀的關(guān)系，KAP1.3基因S4位點AB基因型與遼寧絨山羊產(chǎn)絨量相關(guān)，KAP1.3基因S3位點AA基因型與內(nèi)蒙古白絨山羊體重相關(guān)，KAP1.3基因S2位點BB基因型與遼寧絨山羊體重相關(guān)。馮靜等［3］應(yīng)用PCR-RFLP方法分析綿羊KAP1.3基因遺傳多態(tài)性，結(jié)果表明KAP1.3基因存在2個SNP位點。王春昕［18］應(yīng)用PCR-SSCP技術(shù)對羊毛纖維組成蛋白中的KAP1.3基因部分序列進(jìn)行多態(tài)性分析，結(jié)果表明KAP1.3基因與產(chǎn)毛量和拉伸長度相關(guān)。

表7 C41006938T和T41007938C不同組合基因型對毛性狀的影響Table 7 Effect of C41006938T and T41007938C different combination genotype on wool traits

表8 KAP16不同組合基因型與毛性狀間的關(guān)系Table 8 KAP16 and different combination of relationship between genotype and wool traits

3.2 國內(nèi)外對綿羊KAP基因家族的相關(guān)研究

目前，該領(lǐng)域的研究熱點之一是關(guān)于羊毛細(xì)度決定位點或相關(guān)基因的研究Mclaren等［12］通過連鎖分析已經(jīng)將角蛋白輔助蛋白KAP6.1，KAP7和KAP8定位到綿羊的第1號染色體上。Ponz等［19］通過片段定位方法尋找控制羊毛性狀的QTL，利用法國的合成品系INR401對40個微衛(wèi)星標(biāo)記進(jìn)行了分析，結(jié)果發(fā)現(xiàn)了在3號、7號、25號染色體上有控制毛長(SL)的QTL，在6號染色體上有控制毛纖維直徑(WFD)的QTL，在7號染色體上有控制毛纖維直徑纖維直徑變異系數(shù)的QTL。Cockett N E［20］認(rèn)為KAP6.n和羊毛直徑有關(guān)。Chen等［21］研究表明KAP8.1和KAP1.3基因表達(dá)模式相似，在羊皮膚中KAP1.3高度表達(dá)，KAP8.1表達(dá)較弱。此外，KAP8.1和KAP1.3的位點多態(tài)性與纖維直徑差異顯著。KAP8.1中BC基因型對纖維直徑有極顯著的影響，KAP1.3中AA基因型對纖維直徑有極顯著的影響。此外，KAP8.1與KAP1.3基因的最優(yōu)組合基因型為BC-AA組合基因型。

3.3 KAP16基因與中國美利奴羊(新疆型)毛性狀的關(guān)連分析

羊毛性狀是細(xì)毛羊重要的經(jīng)濟(jì)性狀之一，分子標(biāo)記成為研究重要表型的方法之一，并被廣泛用于動物的遺傳研究中。KAP16基因通過全基因組測序攻擊獲得8個SNPs突變位點，本文研究了其中的兩個突變點對中國美利奴(新疆型)的遺傳效應(yīng)分析，將繼續(xù)研究其他位點對中國美利奴(新疆型)的遺傳效應(yīng)。

C41006938T不同基因型對體格大小和鑒定時體重有顯著影響。群別對鑒定時體重、纖維細(xì)度、體格大小和纖維直徑有極顯著影響。C41006938T基因位點對鑒定時體重和體格大小差異顯著; T41007938C不同基因型對中國美利奴(新疆型)的纖維直徑變異系數(shù)有極顯著影響;對剪毛量和彎曲有顯著影響;不同群別對中國美利奴(新疆型)鑒定是體重、纖維細(xì)度、體格大小和纖維直徑有極顯著影響。T41007938C基因型對中國美利奴羊(新疆型)纖維直徑變異系數(shù)差異極顯著;對剪毛量和彎曲差異顯著。C41006938T和T41007938C對剪毛量有極顯著影響;對纖維直徑變異系數(shù)有顯著影響;不同群別對鑒定時體重、纖維細(xì)度、體格大小、剪毛后體重、纖維直徑有極顯著影響。KAP16組合基因型對中國美利奴(新疆型)剪毛量差異極顯著;對纖維直徑變異系數(shù)差異顯著。

在所研究的綿羊群體中，KAP16基因的2個突變位點均呈現(xiàn)出中度多態(tài)性，說明在這些位點處中國美利奴(新疆型)群體遺傳變異處于中等水平，變異程度相似，沒有太大的差異，同時也說明在上述各群體的培育改良中可以加強(qiáng)人工選擇的強(qiáng)度。

4 結(jié)論

本試驗首次通過重測序技術(shù)，發(fā)現(xiàn)KAP16基因在染色體11上第4外顯子SNP突變點為T41007938C，第5外顯子SNP突變點為C41006938T。通過分析可知，C41006938T不同基因型對體格大小和鑒定時體重有顯著影響; T41007938C不同基因型對纖維直徑纖維直徑變異系數(shù)有極顯著影響，對剪毛量和彎曲有顯著影響; C41006938T和T41007938C組合基因型對剪毛量有極顯著影響，對纖維直徑纖維直徑變異系數(shù)有顯著影響。

［1］田月珍.中國美利奴羊(新疆型)母羊羊毛細(xì)度相關(guān)候選基因的篩選及其關(guān)聯(lián)性分析［D］.烏魯木齊:新疆農(nóng)業(yè)大學(xué)，2012.

［2］汪玲.KAP7.1、KAP8.2基因在遼寧新品系絨山羊興盛期毛囊中的表達(dá)及分析［D］.大連:遼寧師范大學(xué)，2010.

［3］馮靜，宋天增，楊華，等.綿羊KAP1.3基因遺傳多態(tài)性分析［J］.西南農(nóng)業(yè)學(xué)報，2010，23(6):2074－2078.

［4］馮靜，楊華，宋天增，等.綿羊KAP1.1基因克隆與單核苷酸多態(tài)性分析［J］.西南農(nóng)業(yè)學(xué)報，2010，23(5):1720－1724.

［5］馮靜，宋天增，楊華，等.綿羊KAP6.1基因的遺傳多態(tài)性分析［J］.西南農(nóng)業(yè)學(xué)報，2012，25(1):276－278.

［6］瑪依拉.6個KRT基因的遺傳多態(tài)性及其與中國美利奴羊(新疆型)毛性狀的關(guān)聯(lián)性分析［D］.烏魯木齊:新疆農(nóng)業(yè)大學(xué)，2013.

［7］艾買提.5個KAP基因在中國美利奴羊(新疆型)的多態(tài)性及其與毛性的關(guān)聯(lián)性分析［D］.烏魯木齊:新疆農(nóng)業(yè)大學(xué)，2013.

［8］Cheng C T，Kuo C Y，David K A.KAPtain in charge ofmultiplemissions:Emerging rolesof KAP1［J］.World Journal of Biological Chemistry，2014，8.26，5(3):308－320.

［9］趙志東.KAP9.2基因的群體遺傳變異及羊絨生長期和休止期基因的表達(dá)研究［D］.咸陽:西北農(nóng)林科技大學(xué)，2012.

［10］于輝.絨山羊KAP基因家族與Hoxc13基因遺傳變異及其對產(chǎn)絨性狀的影響［D］.咸陽:西北農(nóng)林科技大學(xué)，2009.

［11］趙苗.絨山羊6個候選基因遺傳變異及其與經(jīng)濟(jì)性狀關(guān)系研究［D］.咸陽:西北農(nóng)林科技大學(xué)，2008.

［12］Mclaren R J，Rogers G R，Davies K P，et al.Linkage mapping of wool keratin and keratin-associated protein genes in sheep［J］.Mammali-an Genome，1997(8):938－940.

［13］Parsons Y M，Cooper DW，Piper L R.Evidence of linkage between high-glycine-tyrosine keratin gene loci and wool fibre diameter in a merino half-sib family［J］.Anim Genet，1994，25:105－108.

［14］Henry H M.A genome screen for QTL forwool traits in a Merino X Romney backcross flock［J］.Mammalian Genome，1998(8):938－940.

［15］Itenge M，Weza TQ，F(xiàn)orrestR H，etal.Polymorphism of the KAP1.1，KAP1.3 and K33 genes in Merino sheep［J］.Mol Cell Probes，2007，21(5－6):338－342.

［16］劉桂芬.新疆優(yōu)質(zhì)細(xì)毛羊遺傳多樣性及羊毛細(xì)度候選基因的分析［D］.烏魯木齊:新疆農(nóng)業(yè)大學(xué)，2005.

［17］張亞妮.兩個絨山羊品種微衛(wèi)星多態(tài)性與產(chǎn)絨性能的分子標(biāo)記研究［D］.咸陽:西北農(nóng)林科技大學(xué)，2005.

［18］王春昕，張金新，金海國，等.KAP1.3基因與綿羊產(chǎn)毛性能的相關(guān)分析［J］.吉林畜牧獸醫(yī)，2009，30(4):5－6.

［19］Ponz R，Moreno C，Daniel A，et al.Assessment of genetic variation explained bymarkers for wool traits in sheep via a segmentmapping approach［J］.Mammalian Genome，2001，12:569－572.

［20］Cockett N E.Analysis of the sheep genome［J］.Physiol Genomics，2001(7):69－81.

［21］Chen H Y，Zeng X C，HuiW Q，et al.Developmental expression patterns and association analysisof sheep KAP8.1 and KAP1.3 genes in Chinese Merino sheep［J］.Indian Journal of Animal Sciences，2011，81(4):391－96.

(責(zé)任編輯 李潔)

Analysis of Genetic Effects of KAP16 Gene w ith M ain Econom ic Traits in Fine W ool Sheep

HE Jun-min1，HUANG Xi-xia1*，TIAN Ke-chuan2*，ZHAO Bing-ru1，BAIYan2，TIAN Yue-zhen1，XU Xin-ming2，F(xiàn)U Xue-feng2
(1.Animal Science of Xinjiang Agricultural University，Xinjiang Urumqi 830052，China;2.Xinjiang Academy of Animal Husbandry and Animal Husbandry Institute Laboratory of Molecular Genetics，Xinjiang Urumqi830011，China)

The genetics of KAP16 genewas analyzed to find out the polymorphisms and obtain DNA markers that associated with wool traits，to providemolecular bases formolecular breeding on fine wool sheep.In the research，re-sequencing and PCR-SSCPmethodswere used to analyze the genetic variations of2 SNPsmutation points in KAP16 gene and genetic polymorphism in Chinese Merino(Xinjiang type)population.At the same time，the SAS 8.1 least-squares analysis technique was used to study correlation with wool traits.The results showed that，C41006938Tmutation pointswere observed(named as AA，AB and BB，respectively)，the frequency of genotypeswere0.22，0.47 and 0.31，respectively;correspondingly，frequenciesofallele A and Bwere0.46 and 0.54，among them，B allelewas dominantallele.The frequency of genotypes of AA and AB at T41007938Cmutation pointswere0.77 and 0.23，respectively;frequenciesof allele A and Bwere 0.89 and 0.11，A allelewas dominantallele.χ2 test raveled that the C41006938T and T41007938C pointwere in the Hardy-Weinberg disequilibrium(P＞0.05).The polymorphic information content of C41006938T and T41007938C gene points were 0.37 and 0.18;therefore，Chinese Merino(Xinjiang type)population was reached themoderate level in genetic variation of the gene.Different genotypes on the physical size and identification has obvious effect on body weight(P＜0.05)on C41006938T;different genotype has a significant effect on fiber diameter variation coefficient(P＜0.01)on T41007938C，There is a significant quantity of shearing and bending effect(P＜0.05);C41006938T and T41007938C combination genotype has a significant effect on shearing quantity(P＜0.01)and fiber diameter variation coefficient(P＜0.05).

KAP16 gene;SNPs;Fine wool sheep;Correlation analysis

S814.8

A

1001－4829(2017)2－0458－08

10.16213/j.cnki.scjas.2017.2.036

2016－04－07

國家自然科學(xué)基金(31360543，31460593);農(nóng)業(yè)部體系項目(CARS-40);新疆農(nóng)業(yè)大學(xué)產(chǎn)學(xué)研聯(lián)合培養(yǎng)研究生示范基地項目(xjaucxy-yjs-20152011)

何軍敏(1991－)，女，新疆烏蘇人，碩士研究生，主要從事動物遺傳育種的研究，E-mail:1415230131@qq.com，Tel: 15299126423，*為通訊作者:黃錫霞(1963－)，教授，研究方向為動物遺傳育種，E-mail:au-huangxixia@163.com，田可川(1963－)，研究員，研究方向為動物遺傳育種，E-mail:tiankechuan@ 163.com。