• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

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

    2017-04-13 05:34:58何軍敏黃錫霞田可川趙冰茹柏妍田月珍徐新明付雪峰
    西南農(nóng)業(yè)學(xué)報 2017年2期
    關(guān)鍵詞:剪毛多態(tài)性基因型

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

    (1.新疆農(nóng)業(yè)大學(xué)動物科學(xué)學(xué)院,新疆烏魯木齊830052;2.新疆畜牧科學(xué)院畜牧研究所,分子遺傳重點(diǎn)實(shí)驗(yàn)室,新疆烏魯木齊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é)院畜牧研究所,分子遺傳重點(diǎn)實(shí)驗(yàn)室,新疆烏魯木齊830011)

    為了分析KAP16基因的遺傳多態(tài)性,獲取相關(guān)毛性狀的DNA標(biāo)記,為細(xì)毛羊在分子育種上提供依據(jù),采用重測序技術(shù)和PCR-SSCP技術(shù),驗(yàn)證KAP16候選基因的2個SNPs在中國美利奴(新疆型)群體中的遺傳多態(tài)性。同時,運(yùn)用SAS 8.1軟件進(jìn)行最小二乘方差分析方法,研究其與毛性狀的關(guān)聯(lián)性。結(jié)果表明,C41006938T突變點(diǎn)得到3種基因型為AA、AB和BB型,該群體中C41006938T突變點(diǎn)基因型頻率為0.22、0.47和0.31;A和B等位基因頻率為0.46和0.54,其中B等位基因?yàn)閮?yōu)勢等位基因。T41007938C基因位點(diǎn)AA和AB基因型頻率為0.77和0.23;A和B等位基因頻率分別為0.89和0.11,A基因?yàn)閮?yōu)勢等位基因。χ2檢驗(yàn)表明,C41006938T和T41007938C位點(diǎn)均處于哈代溫伯格平衡狀態(tài)(P>0.05)。C41006938T和T41007938C突變位點(diǎn)的多態(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 試驗(yàn)材料

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

    1.2 試驗(yàn)方法

    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è)計(jì),由上海生工生物工程有限公司合成,引物序列如表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)計(jì)與分析

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

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

    通過整理共得到有效數(shù)據(jù)832條,并對毛性狀進(jìn)行遺傳效應(yīng)分析。運(yùn)用SAS 8.1統(tǒng)計(jì)軟件檢測不同基因型與毛性狀的關(guān)聯(lián)性,并估計(jì)最小二乘均值。最小二乘方差分析數(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突變位點(diǎn)分別位于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),實(shí)驗(yàn)提取的基因組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位點(diǎn)存在3種不同的帶型,將其命名為AA、BB和AB型;P2位點(diǎn)存在,2種基因型,命名為AA、AB型。

    2.5 KAP16基因位點(diǎn)多序列比

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

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

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

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

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

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

    表4 基因頻率和基因型頻率的分布以及Hardy-Weinberg平衡的卡方檢驗(yàn)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基因位點(diǎn)與羊毛性狀的關(guān)聯(lián)性分析

    2.7.1 KAP16基因位點(diǎn)不同基因型的最小二乘方差分析由表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基因位點(diǎn)不同基因型與毛性狀的關(guān)聯(lián)性分析由表6可見,C41006938T位點(diǎn)BB基因型對AA和AB基因型鑒定時體重和體格大小差異顯著(P<0.05);T41007938C位點(diǎn)AB基因型對AA基因型纖維直徑變異系數(shù)差異極顯著(P<0.01),對剪毛量和彎曲差異顯著(P<0.05)。

    2.8 KAP16基因位點(diǎn)不同基因型組合效應(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為候選基因,證實(shí)這兩個基因位點(diǎn)和羊毛產(chǎn)量及纖維直徑存在明顯的相關(guān)性,推測在KAP6和KAP8同一區(qū)域可能存在控制羊毛直徑的主基因。Henry等[14]利用超細(xì)的美利奴羊和羅姆尼羊的回交群體,測定了包括周歲體重、羊毛產(chǎn)量、平均纖維直徑在內(nèi)的多個性狀,同時用216個微衛(wèi)星標(biāo)記進(jìn)行了基因組掃描,找到了一個與羊毛纖維直徑相關(guān)的位點(diǎn)。Itenge-Mweza[15]在澳洲美利奴羊中報道KAP1.3 (B2C)基因存在9個獨(dú)特的SSCP帶型,序列分析顯示這9個帶型代表8個新的KAP1.3等位基因,在KAP1.3基因c.240G>A突變位點(diǎn)的3種基因型分別為GG、AA和AG,在KAP1.3基因c.312G>A突變位點(diǎn)的GG、AA和AG等3種基因型。KRT和KAP多態(tài)性位點(diǎn)可能導(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ū)域中位點(diǎn)W08667與羊毛細(xì)度顯著相關(guān)。張亞妮[17]應(yīng)用PCR-SSCP技術(shù)分析了內(nèi)蒙古白絨山羊和遼寧絨山羊KAP基因與經(jīng)濟(jì)性狀的關(guān)系,KAP1.3基因S4位點(diǎn)AB基因型與遼寧絨山羊產(chǎn)絨量相關(guān),KAP1.3基因S3位點(diǎn)AA基因型與內(nèi)蒙古白絨山羊體重相關(guān),KAP1.3基因S2位點(diǎn)BB基因型與遼寧絨山羊體重相關(guān)。馮靜等[3]應(yīng)用PCR-RFLP方法分析綿羊KAP1.3基因遺傳多態(tài)性,結(jié)果表明KAP1.3基因存在2個SNP位點(diǎn)。王春昕[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)域的研究熱點(diǎn)之一是關(guān)于羊毛細(xì)度決定位點(diǎn)或相關(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的位點(diǎn)多態(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突變位點(diǎn),本文研究了其中的兩個突變點(diǎn)對中國美利奴(新疆型)的遺傳效應(yīng)分析,將繼續(xù)研究其他位點(diǎn)對中國美利奴(新疆型)的遺傳效應(yīng)。

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

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

    4 結(jié)論

    本試驗(yàn)首次通過重測序技術(shù),發(fā)現(xiàn)KAP16基因在染色體11上第4外顯子SNP突變點(diǎn)為T41007938C,第5外顯子SNP突變點(diǎn)為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è)部體系項(xiàng)目(CARS-40);新疆農(nóng)業(yè)大學(xué)產(chǎn)學(xué)研聯(lián)合培養(yǎng)研究生示范基地項(xiàng)目(xjaucxy-yjs-20152011)

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

    猜你喜歡
    剪毛多態(tài)性基因型
    單核苷酸多態(tài)性與中醫(yī)證候相關(guān)性研究進(jìn)展
    機(jī)械剪毛中羊只的保定及高效剪毛方式
    草食家畜(2019年4期)2019-02-16 22:52:54
    羊駝剪毛“三注意”
    剪毛對小尾寒羊行為和生理的影響
    馬鈴薯cpDNA/mtDNA多態(tài)性的多重PCR檢測
    GlobalFiler~? PCR擴(kuò)增試劑盒驗(yàn)證及其STR遺傳多態(tài)性
    新型雙刃靜音多插頭鼻腔剪毛器的研制與臨床應(yīng)用
    西安地區(qū)育齡婦女MTHFRC677T基因型分布研究
    BAMBI基因敲除小鼠的繁育、基因型鑒定
    CYP3A4*1G基因多態(tài)性及功能的初步探討
    欧美精品国产亚洲| 不卡视频在线观看欧美| 蜜桃在线观看..| 80岁老熟妇乱子伦牲交| 国产毛片在线视频| 80岁老熟妇乱子伦牲交| 精品亚洲乱码少妇综合久久| 一级毛片我不卡| 91精品伊人久久大香线蕉| 久久午夜福利片| 夜夜骑夜夜射夜夜干| 午夜日韩欧美国产| 亚洲,一卡二卡三卡| 中文字幕av电影在线播放| 啦啦啦在线免费观看视频4| 丝袜美足系列| 色婷婷久久久亚洲欧美| 亚洲少妇的诱惑av| 精品国产乱码久久久久久男人| 免费播放大片免费观看视频在线观看| 免费观看性生交大片5| 国产精品国产三级国产专区5o| 欧美变态另类bdsm刘玥| 丝袜人妻中文字幕| 亚洲成人手机| 精品亚洲乱码少妇综合久久| av.在线天堂| 少妇人妻久久综合中文| 一本久久精品| 高清在线视频一区二区三区| 熟女电影av网| 日韩三级伦理在线观看| av国产久精品久网站免费入址| 蜜桃国产av成人99| 伦精品一区二区三区| av.在线天堂| 成人国产av品久久久| 亚洲国产欧美网| 欧美日韩亚洲高清精品| 美女高潮到喷水免费观看| 制服诱惑二区| 成人漫画全彩无遮挡| 最近的中文字幕免费完整| 国产野战对白在线观看| 精品国产乱码久久久久久男人| 久久99一区二区三区| 久久久精品国产亚洲av高清涩受| 亚洲欧洲日产国产| 好男人视频免费观看在线| 国产深夜福利视频在线观看| 国产白丝娇喘喷水9色精品| 街头女战士在线观看网站| av又黄又爽大尺度在线免费看| 91精品伊人久久大香线蕉| 超碰97精品在线观看| 人妻 亚洲 视频| 成年美女黄网站色视频大全免费| 最近最新中文字幕免费大全7| 精品少妇久久久久久888优播| 少妇熟女欧美另类| 1024香蕉在线观看| av网站免费在线观看视频| 一级毛片 在线播放| 免费观看无遮挡的男女| 免费黄网站久久成人精品| 欧美人与善性xxx| 欧美亚洲 丝袜 人妻 在线| 亚洲精品国产一区二区精华液| 人妻人人澡人人爽人人| 免费黄频网站在线观看国产| 国产av码专区亚洲av| 久久久久精品久久久久真实原创| 久久青草综合色| 国产日韩欧美视频二区| 久久亚洲国产成人精品v| 日韩人妻精品一区2区三区| 黑人猛操日本美女一级片| 亚洲四区av| 老女人水多毛片| 日韩熟女老妇一区二区性免费视频| 欧美精品一区二区免费开放| 精品第一国产精品| 国产精品久久久久久精品古装| 午夜免费男女啪啪视频观看| 欧美日韩成人在线一区二区| 一区二区三区四区激情视频| 国产爽快片一区二区三区| 高清视频免费观看一区二区| 国产熟女午夜一区二区三区| 大片免费播放器 马上看| 亚洲国产成人一精品久久久| 成人亚洲精品一区在线观看| 一本色道久久久久久精品综合| 亚洲婷婷狠狠爱综合网| 在线观看人妻少妇| 黄频高清免费视频| 午夜激情久久久久久久| 国产精品久久久久久精品电影小说| 国产野战对白在线观看| 国产女主播在线喷水免费视频网站| 国产成人免费观看mmmm| 亚洲四区av| av在线app专区| videosex国产| 视频区图区小说| 91在线精品国自产拍蜜月| 欧美成人精品欧美一级黄| 久久av网站| 亚洲av国产av综合av卡| 久久综合国产亚洲精品| 日韩一区二区三区影片| 国产在线免费精品| 国产精品亚洲av一区麻豆 | 一二三四在线观看免费中文在| 日韩欧美一区视频在线观看| 男女下面插进去视频免费观看| 夫妻午夜视频| 日韩中文字幕欧美一区二区 | 晚上一个人看的免费电影| 国产探花极品一区二区| 日本vs欧美在线观看视频| 高清欧美精品videossex| 国产成人精品一,二区| 国产成人精品无人区| 多毛熟女@视频| 少妇熟女欧美另类| 成人亚洲欧美一区二区av| 中国三级夫妇交换| 在线免费观看不下载黄p国产| 99re6热这里在线精品视频| 日韩中字成人| videossex国产| 国产一区有黄有色的免费视频| 电影成人av| 美女xxoo啪啪120秒动态图| 免费看av在线观看网站| 精品国产乱码久久久久久男人| www.精华液| 麻豆精品久久久久久蜜桃| 日本黄色日本黄色录像| 日韩大片免费观看网站| 婷婷色av中文字幕| av.在线天堂| 丰满乱子伦码专区| 久久人人97超碰香蕉20202| 熟女电影av网| 综合色丁香网| 精品一区二区三卡| 日本午夜av视频| 精品卡一卡二卡四卡免费| 国产一区有黄有色的免费视频| 亚洲精品视频女| 黄片播放在线免费| 精品亚洲成国产av| 亚洲av.av天堂| 亚洲一区中文字幕在线| www.av在线官网国产| 女人精品久久久久毛片| 在线免费观看不下载黄p国产| 看免费成人av毛片| 午夜激情久久久久久久| 久久人人97超碰香蕉20202| 亚洲第一区二区三区不卡| 多毛熟女@视频| 欧美成人精品欧美一级黄| 中文字幕人妻丝袜制服| 啦啦啦中文免费视频观看日本| 日本91视频免费播放| 男女下面插进去视频免费观看| 黑人欧美特级aaaaaa片| 亚洲国产毛片av蜜桃av| a级毛片在线看网站| 亚洲色图综合在线观看| 涩涩av久久男人的天堂| 黄片无遮挡物在线观看| 我要看黄色一级片免费的| 少妇的逼水好多| 美女主播在线视频| 欧美日韩精品网址| 人体艺术视频欧美日本| 久久人人爽av亚洲精品天堂| 亚洲三级黄色毛片| 成人手机av| 王馨瑶露胸无遮挡在线观看| 久久精品国产a三级三级三级| 伦精品一区二区三区| 精品一区二区免费观看| 亚洲欧美一区二区三区久久| 日本欧美视频一区| 国产欧美日韩综合在线一区二区| 国产精品99久久99久久久不卡 | 日本午夜av视频| 亚洲四区av| 国产精品二区激情视频| 2018国产大陆天天弄谢| 久久精品国产亚洲av天美| 国产成人精品无人区| 久久99热这里只频精品6学生| videos熟女内射| 麻豆av在线久日| 国产成人免费无遮挡视频| 巨乳人妻的诱惑在线观看| 大香蕉久久成人网| 亚洲精品国产av成人精品| 免费在线观看完整版高清| 国产一区二区在线观看av| 久久精品夜色国产| 国产熟女午夜一区二区三区| 一二三四在线观看免费中文在| 亚洲成色77777| www.精华液| 九九爱精品视频在线观看| 日韩一本色道免费dvd| 久久97久久精品| 日本欧美视频一区| 色婷婷av一区二区三区视频| 亚洲av综合色区一区| 国产熟女欧美一区二区| 国产精品女同一区二区软件| 精品人妻在线不人妻| 久久这里只有精品19| 一边摸一边做爽爽视频免费| 免费观看av网站的网址| 卡戴珊不雅视频在线播放| 久久狼人影院| 99香蕉大伊视频| 成年人午夜在线观看视频| 亚洲久久久国产精品| 国产精品亚洲av一区麻豆 | 十八禁网站网址无遮挡| √禁漫天堂资源中文www| 日韩制服丝袜自拍偷拍| 中文欧美无线码| 色婷婷久久久亚洲欧美| 国产精品av久久久久免费| 亚洲av电影在线观看一区二区三区| 国产在线一区二区三区精| 这个男人来自地球电影免费观看 | videos熟女内射| 香蕉国产在线看| 超碰成人久久| 精品国产一区二区三区久久久樱花| 激情五月婷婷亚洲| 高清在线视频一区二区三区| 久久精品国产自在天天线| 欧美日韩av久久| av有码第一页| 成人影院久久| 日韩精品免费视频一区二区三区| 久久久久视频综合| 久久精品国产亚洲av天美| 超色免费av| 丰满少妇做爰视频| 黄片小视频在线播放| 日韩不卡一区二区三区视频在线| 老汉色∧v一级毛片| 高清黄色对白视频在线免费看| 午夜福利视频精品| 国产精品国产av在线观看| 日韩不卡一区二区三区视频在线| 999精品在线视频| 男女下面插进去视频免费观看| 久久 成人 亚洲| 欧美变态另类bdsm刘玥| 国产激情久久老熟女| 日韩制服丝袜自拍偷拍| 欧美日韩精品成人综合77777| 在线 av 中文字幕| 美女高潮到喷水免费观看| 国产高清国产精品国产三级| 男女无遮挡免费网站观看| 成人国产麻豆网| 久久精品夜色国产| 精品人妻偷拍中文字幕| 99久久精品国产国产毛片| 国产一区二区激情短视频 | 老汉色av国产亚洲站长工具| 高清av免费在线| 天天躁夜夜躁狠狠躁躁| 日韩中字成人| 一边摸一边做爽爽视频免费| 1024视频免费在线观看| 免费日韩欧美在线观看| 叶爱在线成人免费视频播放| 国产成人a∨麻豆精品| a级片在线免费高清观看视频| 中国三级夫妇交换| 国产野战对白在线观看| 亚洲国产av新网站| 18禁观看日本| 午夜av观看不卡| 777久久人妻少妇嫩草av网站| 大陆偷拍与自拍| 久久97久久精品| 18+在线观看网站| 中文字幕色久视频| 日本91视频免费播放| 久久精品久久久久久噜噜老黄| 一级黄片播放器| 纵有疾风起免费观看全集完整版| 美国免费a级毛片| 国产精品.久久久| 国产一区二区 视频在线| 精品国产一区二区三区四区第35| 久久久久视频综合| 天天躁狠狠躁夜夜躁狠狠躁| 日日撸夜夜添| 精品亚洲乱码少妇综合久久| 成人国产麻豆网| 美女视频免费永久观看网站| 一级片免费观看大全| 婷婷色综合www| 成年人免费黄色播放视频| 欧美少妇被猛烈插入视频| 欧美变态另类bdsm刘玥| 黑人猛操日本美女一级片| 天天影视国产精品| 成人国产麻豆网| 国产精品蜜桃在线观看| 亚洲精品久久久久久婷婷小说| 亚洲精品日本国产第一区| 人人妻人人澡人人爽人人夜夜| 亚洲 欧美一区二区三区| 18禁国产床啪视频网站| 国产激情久久老熟女| 一个人免费看片子| 日韩中文字幕视频在线看片| 九九爱精品视频在线观看| 国产精品偷伦视频观看了| 成年人午夜在线观看视频| 黄频高清免费视频| 免费在线观看视频国产中文字幕亚洲 | 中文天堂在线官网| 一级毛片我不卡| 夫妻性生交免费视频一级片| 日本爱情动作片www.在线观看| 国产伦理片在线播放av一区| 成人漫画全彩无遮挡| kizo精华| 制服人妻中文乱码| 宅男免费午夜| 两个人看的免费小视频| 国产毛片在线视频| 久久人人爽av亚洲精品天堂| 不卡视频在线观看欧美| 日本爱情动作片www.在线观看| 欧美激情 高清一区二区三区| 伦理电影免费视频| 国产精品国产av在线观看| 欧美国产精品一级二级三级| 日韩精品免费视频一区二区三区| 夫妻午夜视频| 国产午夜精品一二区理论片| 免费久久久久久久精品成人欧美视频| 巨乳人妻的诱惑在线观看| 成年女人毛片免费观看观看9 | 午夜福利影视在线免费观看| 中文字幕人妻丝袜制服| 老司机影院成人| 欧美日韩一区二区视频在线观看视频在线| 最近手机中文字幕大全| 免费观看在线日韩| 国产精品久久久久成人av| 久久女婷五月综合色啪小说| 久久精品国产综合久久久| 国产成人精品福利久久| 午夜福利视频精品| 欧美最新免费一区二区三区| 一级毛片我不卡| 久久久久久久久久人人人人人人| 又黄又粗又硬又大视频| 欧美激情高清一区二区三区 | 一个人免费看片子| 精品久久蜜臀av无| 男女无遮挡免费网站观看| 亚洲精品国产一区二区精华液| 日韩视频在线欧美| 男女国产视频网站| 亚洲第一青青草原| av国产精品久久久久影院| 国产探花极品一区二区| 黄色毛片三级朝国网站| 亚洲成人av在线免费| 91精品伊人久久大香线蕉| 午夜福利在线免费观看网站| 久久午夜福利片| 丝袜美腿诱惑在线| 热99国产精品久久久久久7| 天堂中文最新版在线下载| 国语对白做爰xxxⅹ性视频网站| 夫妻性生交免费视频一级片| 亚洲综合精品二区| 欧美精品一区二区大全| 黄色视频在线播放观看不卡| 午夜日本视频在线| 日韩一区二区三区影片| 制服诱惑二区| 精品99又大又爽又粗少妇毛片| 亚洲国产av新网站| av在线老鸭窝| 熟女少妇亚洲综合色aaa.| 国产色婷婷99| 成人影院久久| 久久久久国产网址| 国产免费视频播放在线视频| 国产av精品麻豆| 午夜福利影视在线免费观看| 成人毛片a级毛片在线播放| 在线 av 中文字幕| 国产精品国产三级专区第一集| 青春草国产在线视频| 国产精品成人在线| 纯流量卡能插随身wifi吗| 国产欧美日韩综合在线一区二区| 国产乱人偷精品视频| 日日撸夜夜添| 国产xxxxx性猛交| 国产激情久久老熟女| kizo精华| av国产精品久久久久影院| 90打野战视频偷拍视频| 咕卡用的链子| 亚洲av欧美aⅴ国产| 黄色视频在线播放观看不卡| 国产精品无大码| 秋霞在线观看毛片| 男人操女人黄网站| 久久久久久久国产电影| 18禁裸乳无遮挡动漫免费视频| 青春草视频在线免费观看| 欧美精品国产亚洲| 亚洲,欧美精品.| 美女xxoo啪啪120秒动态图| 国产欧美亚洲国产| 性色avwww在线观看| 制服丝袜香蕉在线| 777久久人妻少妇嫩草av网站| av在线观看视频网站免费| 久久久精品区二区三区| 国产麻豆69| 女人高潮潮喷娇喘18禁视频| 色哟哟·www| 1024香蕉在线观看| 美女福利国产在线| 丝袜美腿诱惑在线| 一二三四在线观看免费中文在| 亚洲人成电影观看| 女的被弄到高潮叫床怎么办| 亚洲av男天堂| 成年动漫av网址| 久久久精品94久久精品| 久久久久国产网址| 一本久久精品| 日韩av不卡免费在线播放| 秋霞在线观看毛片| 人人妻人人澡人人看| 精品亚洲成a人片在线观看| 在线观看三级黄色| 亚洲国产精品999| 日韩成人av中文字幕在线观看| 国产成人精品在线电影| 90打野战视频偷拍视频| 水蜜桃什么品种好| 中文字幕人妻丝袜一区二区 | 免费播放大片免费观看视频在线观看| 99热全是精品| 国产精品 国内视频| 天天躁日日躁夜夜躁夜夜| 亚洲精品美女久久av网站| 日韩不卡一区二区三区视频在线| 午夜av观看不卡| 欧美成人午夜精品| 免费人妻精品一区二区三区视频| 一二三四在线观看免费中文在| 国产免费福利视频在线观看| 一级毛片我不卡| 欧美激情高清一区二区三区 | 综合色丁香网| 国产乱来视频区| 人体艺术视频欧美日本| 久久精品国产综合久久久| 你懂的网址亚洲精品在线观看| 国产在视频线精品| 国产在线视频一区二区| 久久午夜综合久久蜜桃| av视频免费观看在线观看| 精品午夜福利在线看| 欧美日韩亚洲高清精品| 免费人妻精品一区二区三区视频| 一级黄片播放器| 久久99一区二区三区| 日本欧美国产在线视频| 午夜福利在线观看免费完整高清在| 国产综合精华液| av视频免费观看在线观看| 日本猛色少妇xxxxx猛交久久| 国产野战对白在线观看| 一区二区av电影网| 国产亚洲一区二区精品| 成人毛片60女人毛片免费| 美女福利国产在线| 亚洲精品久久成人aⅴ小说| 国产成人免费观看mmmm| 熟妇人妻不卡中文字幕| 久久人妻熟女aⅴ| 少妇的丰满在线观看| 亚洲精品国产av蜜桃| 国产极品天堂在线| 国产老妇伦熟女老妇高清| 免费高清在线观看日韩| 天美传媒精品一区二区| 久热这里只有精品99| 中文字幕色久视频| 国产在线一区二区三区精| 成人国产麻豆网| 国产成人欧美| 国产高清不卡午夜福利| 亚洲婷婷狠狠爱综合网| 捣出白浆h1v1| xxxhd国产人妻xxx| 日本免费在线观看一区| 下体分泌物呈黄色| 成人国产麻豆网| 美女午夜性视频免费| 成人18禁高潮啪啪吃奶动态图| 最近中文字幕2019免费版| 最近2019中文字幕mv第一页| 免费在线观看黄色视频的| 国产色婷婷99| 中文字幕av电影在线播放| 寂寞人妻少妇视频99o| 丝袜在线中文字幕| 国产亚洲欧美精品永久| 亚洲久久久国产精品| 热99国产精品久久久久久7| 黄片无遮挡物在线观看| 巨乳人妻的诱惑在线观看| 国产亚洲最大av| 91精品三级在线观看| 国精品久久久久久国模美| 黄色配什么色好看| 啦啦啦中文免费视频观看日本| 一个人免费看片子| 免费人妻精品一区二区三区视频| 国产片特级美女逼逼视频| 黄色一级大片看看| 欧美日韩亚洲国产一区二区在线观看 | 中文字幕av电影在线播放| 国产成人精品福利久久| 午夜福利视频精品| 国产熟女欧美一区二区| 99久国产av精品国产电影| 国产免费福利视频在线观看| 欧美日韩av久久| 日韩av在线免费看完整版不卡| 三级国产精品片| 波多野结衣一区麻豆| 天天躁夜夜躁狠狠躁躁| 日韩大片免费观看网站| 久久国产精品男人的天堂亚洲| 精品久久蜜臀av无| 18+在线观看网站| av天堂久久9| 国产男女超爽视频在线观看| 亚洲欧美一区二区三区国产| 亚洲综合色惰| 久久av网站| 老司机影院成人| 免费播放大片免费观看视频在线观看| 少妇的逼水好多| 国产在视频线精品| 观看av在线不卡| 免费少妇av软件| 卡戴珊不雅视频在线播放| 国产精品 国内视频| 极品少妇高潮喷水抽搐| 亚洲欧洲精品一区二区精品久久久 | 久久久久久久大尺度免费视频| 日韩制服骚丝袜av| 亚洲欧美色中文字幕在线| 国产片内射在线| 涩涩av久久男人的天堂| 最近最新中文字幕免费大全7| 一个人免费看片子| 激情五月婷婷亚洲| 在线精品无人区一区二区三| 日日摸夜夜添夜夜爱| 日韩,欧美,国产一区二区三区| 久久97久久精品| 一级毛片黄色毛片免费观看视频| 夜夜骑夜夜射夜夜干| 天天操日日干夜夜撸| 韩国高清视频一区二区三区| 成人手机av| 免费黄网站久久成人精品| 久久久久视频综合| 80岁老熟妇乱子伦牲交| 精品少妇久久久久久888优播| 久久精品人人爽人人爽视色| 最近手机中文字幕大全| 男女啪啪激烈高潮av片| 免费看av在线观看网站| 亚洲综合色惰| 久久婷婷青草| 国产片内射在线| 美女高潮到喷水免费观看| 美国免费a级毛片| 1024视频免费在线观看| 如何舔出高潮| 五月伊人婷婷丁香| 精品国产一区二区三区四区第35| 精品少妇内射三级| 亚洲精品久久久久久婷婷小说| 久久国产精品男人的天堂亚洲| 亚洲国产欧美网|