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

    采用6個(gè)微衛(wèi)星位點(diǎn)研究6個(gè)綿羊品種遺傳多樣性

    2016-06-23 13:49:50陳李鵬黃勇富趙永聚俄廣鑫娜日蘇
    生物技術(shù)通報(bào) 2016年5期
    關(guān)鍵詞:微衛(wèi)星雜合綿羊

    陳李鵬 黃勇富 趙永聚 俄廣鑫 娜日蘇

    (西南大學(xué)動(dòng)物科技學(xué)院,重慶 400716)

    采用6個(gè)微衛(wèi)星位點(diǎn)研究6個(gè)綿羊品種遺傳多樣性

    陳李鵬 黃勇富 趙永聚 俄廣鑫 娜日蘇

    (西南大學(xué)動(dòng)物科技學(xué)院,重慶 400716)

    利用6個(gè)微衛(wèi)星座位對(duì)中國(guó)烏珠穆沁羊、小尾寒羊、灘羊、昭通綿羊、呼倫貝爾羊和甘肅高山細(xì)毛羊6個(gè)綿羊品種,共計(jì)280個(gè)個(gè)體進(jìn)行遺傳多樣性分析。計(jì)算了6個(gè)綿羊品種間的多態(tài)信息含量、雜合度和遺傳距離,并進(jìn)行了主成分分析和UPGMA聚類(lèi)。發(fā)現(xiàn)了101個(gè)等位基因,平均雜合度0.599-0.691,平均多態(tài)信息含量0.609-0.680;甘肅高山細(xì)毛羊與其他綿羊品種遺傳分歧最大,而呼倫貝爾羊與烏珠穆沁羊間的分歧最小。結(jié)果表明6個(gè)綿羊品種均具有較高的遺傳多樣性,品種間的遺傳關(guān)系與其形成歷史、分化及地理分布基本一致。

    微衛(wèi)星;遺傳多樣性;綿羊

    我國(guó)地方綿羊品種眾多,遺傳資源十分豐富,是世界少有的寶貴基因庫(kù)[1],但受世界范圍內(nèi)商業(yè)品種的沖擊,導(dǎo)致市場(chǎng)競(jìng)爭(zhēng)力低的地方品種逐漸走向衰落。因此,對(duì)中國(guó)綿羊品種遺傳多樣性進(jìn)行評(píng)估和保護(hù),已是當(dāng)務(wù)之急。

    微衛(wèi)星又稱(chēng)簡(jiǎn)單串連重復(fù)序列,一般由2-6 bp組成核心序列,重復(fù)10-20次,首尾相接而組成[2]。它廣泛存在于真核生物的基因組中,可通過(guò)PCR擴(kuò)增和電泳分型,具有數(shù)量多、多態(tài)性豐富、保守性好、共顯性遺傳以及檢測(cè)方便快速等優(yōu)點(diǎn)[3-5],在度量品種遺傳多樣性、估測(cè)品種間遺傳距離及構(gòu)建系統(tǒng)發(fā)生樹(shù)等研究中顯示出巨大的優(yōu)勢(shì)。目前,有很多學(xué)者應(yīng)用微衛(wèi)星DNA標(biāo)記技術(shù)對(duì)綿羊的遺傳多樣性進(jìn)行了研究。如Sun等[6]利用15個(gè)微衛(wèi)星標(biāo)記評(píng)估了5個(gè)綿羊品種的遺傳多樣性,結(jié)果表明群體遺傳結(jié)構(gòu)與群體間的遺傳距離不一致。Agaviezor等[7]利用15個(gè)微衛(wèi)星標(biāo)記分析了4個(gè)尼日利亞綿羊品種共384只綿羊的遺傳多樣性,結(jié)果表明綿羊品種內(nèi)遺傳變異高于品種間遺傳變異。Yilmaz等[8]利用18個(gè)微衛(wèi)星標(biāo)記分析了9個(gè)綿羊品種的遺傳多樣性,發(fā)現(xiàn)歐洲綿羊品種的遺傳變異高于土耳其綿羊品種;Leroy等[9]利用21個(gè)微衛(wèi)星標(biāo)記分析了49個(gè)綿羊品種的遺傳多樣性,評(píng)估了法國(guó)部分綿羊品種的起源分化模式。?urkovi?等[10]利用30個(gè)微衛(wèi)星分析了18個(gè)綿羊品種的遺傳多樣性,評(píng)估了綿羊品種形成過(guò)程中自然和人工選擇的重要性。

    本研究選取呼倫貝爾羊、烏珠穆沁羊、小尾寒羊、灘羊和昭通綿羊等5個(gè)中國(guó)地方綿羊品種,以及甘肅高山細(xì)毛羊1個(gè)培育品種,在常染色體基因組水平上研究6個(gè)微衛(wèi)星位點(diǎn)的遺傳多樣性和遺傳分化關(guān)系,旨在為揭示其品種遺傳結(jié)構(gòu)和遺傳背景,以及下一步開(kāi)展資源保護(hù)與開(kāi)發(fā)利用提供科學(xué)依據(jù)。

    1 材料與方法

    1.1 材料

    1.1.1 實(shí)驗(yàn)樣本 本實(shí)驗(yàn)樣本采自6個(gè)綿羊品種分布的主產(chǎn)區(qū),每個(gè)品種隨機(jī)采集無(wú)親緣關(guān)系個(gè)體40只或48只,共計(jì)280只。所采血樣為8-10 mL/只,用檸檬酸葡萄糖抗凝,-20℃凍存。各品種的樣本信息,見(jiàn)表1。

    表1 6個(gè)綿羊品種的樣本信息

    1.1.2 微衛(wèi)星引物 根據(jù)相關(guān)微衛(wèi)星選擇標(biāo)準(zhǔn),并綜合考慮微衛(wèi)星位點(diǎn)在染色體中位置、引物序列組成特點(diǎn)、微衛(wèi)星位點(diǎn)多態(tài)性和產(chǎn)物大小等多方面因素,最后確定了分布于5條不同常染色體上的6對(duì)微衛(wèi)星引物,分別為 OarJMP29、OarAE129、MCM-527、ILSTS005、MAF209和OarFCB304(表2)。

    表2 6個(gè)微衛(wèi)星位點(diǎn)的信息

    1.2 方法

    1.2.1 DNA的提取 用基因組DNA提取試劑盒提取全血中的基因組DNA,-20℃保存。

    1.2.2 PCR擴(kuò)增 采用20 μL反應(yīng)體系,其中各組分的終濃度分別為dNTPs 0.2 mmol/L、Mg2+1.5 mmol/L、混合上下游引物0.5 mmol/L、Taq酶5 U/μL、D NA模板1 μL(約60 ng)。PCR反應(yīng)程序:94℃預(yù)變性5 min;94℃變性30 s,50-60℃退火30 s,72℃延伸30 s,經(jīng)過(guò)35個(gè)循環(huán);于72℃延伸7 min,4℃保存。PCR產(chǎn)物用ABI 3130 xl全自動(dòng)基因分析儀進(jìn)行分型檢測(cè),片段分析的上樣操作參照曾盛誠(chéng)[11]的文獻(xiàn)報(bào)道。

    1.2.3 數(shù)據(jù)的統(tǒng)計(jì)處理 利用Microsatellite Toolkit軟件計(jì)算各微衛(wèi)星座位的等位基因數(shù)、多態(tài)信息含量(PIC)、觀察雜合度(HO)和期望雜合度(HE)等參數(shù)[12],Arlequin 3.5.1.3軟件計(jì)算6個(gè)品種間遺傳距離(DS)[13],GENEPOP 3.4軟件進(jìn)行哈德溫伯格平衡檢驗(yàn)[14],F(xiàn)stat 2.9.3.2軟件計(jì)算群體近交系數(shù)(FIS)[15],MVSP3.1軟件進(jìn)行品種間遺傳分化的主成分分析和UPGMA聚類(lèi)分析[16,17]。

    2 結(jié)果

    2.1 等位基因分布

    6個(gè)綿羊品種微衛(wèi)星座位的等位基因分布見(jiàn)表3。由表3可見(jiàn) ,在OarFCB304座位檢測(cè)到的等位基因數(shù)最多,達(dá)到24個(gè),OarAE129座位檢測(cè)到的等位基因數(shù)最少,為9個(gè),6個(gè)綿羊品種在6個(gè)微衛(wèi)星座位中共檢測(cè)到101個(gè)等位基因,平均每個(gè)座位的等位基因數(shù)為16.8;同一個(gè)微衛(wèi)星座位,在不同品種間的等位基因分布也不同,如在OarJMP29座位上,UQ、ZT、HBR、STH、TAN、GSH檢測(cè)到的等位基因數(shù)分別為10、9、11、8、9、7。

    表3 6個(gè)微衛(wèi)星座位在6個(gè)綿羊品種中的等位基因數(shù)

    2.2 品種內(nèi)的遺傳變異檢測(cè)

    2.2.1 觀測(cè)雜合度和期望雜合度 6個(gè)綿羊品種的觀測(cè)雜合度(HO)和期望雜合度(HE)的檢測(cè)結(jié)果見(jiàn)表4。就座位而言,MAF209 座位的平均HO最高,達(dá)到0.764;OarAE129 座位的平均HO最低。就品種而言,6個(gè)綿羊品種的群體雜合度均較高,除ZT的HO為0.599外,其余都在 0.6-0.7之間,屬于高度雜合群體。

    2.2.2 多態(tài)信息含量(PIC) 由表4可知,就座位而言,6個(gè)微衛(wèi)星座位平均PIC值在0.505-0.755 之間,表現(xiàn)為高度多態(tài),MCM527座位平均PIC值最大,達(dá)到 0.755。表明所選6個(gè)微衛(wèi)星標(biāo)記均具有較高的多態(tài)性,能夠提供大量的遺傳信息,可作為有效的分子遺傳標(biāo)記用于群體遺傳多樣性分析。就品種而言,6個(gè)綿羊品種PIC均高于0.5。GSH為培育品種,平均PIC值為0.672,其余5個(gè)地方綿羊品種平均PIC值在0.609-0.680間分布,表明這6個(gè)綿羊品種的遺傳多樣性較高,具有很高的育種潛力。

    2.2.3 Hardy-Weinberg平衡檢驗(yàn) 本研究中Hardy-Weinbery 平衡檢驗(yàn)結(jié)果(表5)顯示,GSH在2個(gè)基因座處于哈溫不平衡、其余5個(gè)綿羊品種均在1個(gè)基因座處于不平衡狀態(tài)。

    2.2.4 群體內(nèi)近交系數(shù) 6個(gè)綿羊品種的近交系數(shù)(FIS)測(cè)定值(表6)顯示,6個(gè)綿羊品種的FIS范圍是0.041-0.166,6個(gè)微衛(wèi)星位點(diǎn)的平均FIS值在-0.005-0.254 7之間,矯正后的P值為0.001 4,除ZT的P值與校正后的P值相同外,其余5個(gè)品種的P值都大于校正后的P值,顯示近交程度不明顯。ZT的FIS值最大(FIS=0.166),其近交程度較高,通過(guò)校正其近交系數(shù)顯著偏離0,這可能與尚未建立ZT的保護(hù)區(qū)和保種場(chǎng)和處于農(nóng)戶(hù)自繁自養(yǎng)狀態(tài)有關(guān)。

    2.3 品種間的遺傳關(guān)系

    2.3.1 遺傳距離 基于SSR數(shù)據(jù)計(jì)算了6個(gè)綿羊品種間的Nei氏標(biāo)準(zhǔn)遺傳距離(DS),結(jié)果見(jiàn)表7。由表7可見(jiàn),UQ、ZT、HBR、STH和TAN 5個(gè)地方品種之間遺傳距離在0.012-0.156之間,HBR與UQ遺傳距離最小,STH與UQ和HBR兩個(gè)品種間的遺傳分歧較小,GSH與其他5個(gè)地方品種之間的遺傳距離均較大。

    表 4 6個(gè)綿羊品種在6個(gè)微衛(wèi)星位點(diǎn)上的遺傳變異參數(shù)

    表5 哈溫平衡檢測(cè)各品種各基因座P值

    2.3.2 聚類(lèi)分析 主成分分析是揭示群體間遺傳分化程度和遺傳關(guān)系的一種聚類(lèi)方法,研究群體間在多個(gè)標(biāo)記座位上的差異性,具有多變量的特征[18]?;谖⑿l(wèi)星數(shù)據(jù),利用MVSP3.1軟件進(jìn)行了主成分分析,結(jié)果如表8所示,在5個(gè)主成分中前3個(gè)主成分累計(jì)貢獻(xiàn)率達(dá)到76.795%,足以代表原始因子的大部分信息,據(jù)此繪制主成分分析散點(diǎn)圖(圖1)。其中第1主成分、第2主成分和第3主成分分別占總變異的38.02%、23.19%和15.58%。

    在前3個(gè)主成分中,第1主成分將GSH與其他5個(gè)地方品種分開(kāi);第2主成分將ZT從5個(gè)地方品種中分出;第3主成分將UQ、HBR和STH分為一類(lèi),與TAN分開(kāi)。同樣利用MVSP3.1軟件進(jìn)行了UPGMA聚類(lèi)分析,結(jié)果如表9所示,據(jù)此得到聚類(lèi)圖(圖2),與主成分分析結(jié)果基本一致。

    表 6 六個(gè)綿羊品種的近交系數(shù)

    表7 6個(gè)綿羊品種間的遺傳距離

    表9 UPGMA聚類(lèi)結(jié)果

    表8 MVSP主成分分析結(jié)果

    3 討論

    3.1 遺傳多樣性和遺傳分化分析

    本研究所用的6個(gè)微衛(wèi)星座位在6個(gè)綿羊品種中均檢測(cè)到9個(gè)以上等位基因,各座位的多態(tài)信息含量都很高。根據(jù)Bolstein等[19]提出的基因變異程度衡量指標(biāo),本研究涉及的6個(gè)微衛(wèi)星座位均屬高度多態(tài)基因座。因此,初步認(rèn)為這6個(gè)微衛(wèi)星標(biāo)記是普遍適合于綿羊遺傳分化的多態(tài)標(biāo)記。

    多態(tài)信息含量(PIC)是群體內(nèi)遺傳變異的量度指標(biāo),可用來(lái)描述微衛(wèi)星位點(diǎn)的變異程度。Ma等[20]利用30個(gè)微衛(wèi)星標(biāo)記分析了12個(gè)地方綿羊品種的遺傳多樣性,平均PIC在0.519-0.666之間;儲(chǔ)明星等[21]利用2個(gè)微衛(wèi)星標(biāo)記對(duì)5個(gè)綿羊品種的研究結(jié)果得出平均PIC為0.60;本研究中6個(gè)微衛(wèi)星均為高度多態(tài),其中以MCM527的多態(tài)性最高,為0.755。

    雜合度又稱(chēng)基因多樣性,反映群體在多個(gè)位點(diǎn)上的遺傳變異。一般認(rèn)為它是度量群體遺傳變異的一個(gè)重要參數(shù)。本實(shí)驗(yàn)所選取的6個(gè)微衛(wèi)星平均HE在0.544-0.795之間,6個(gè)綿羊品種的平均HO為0.599-0.691,表明這6個(gè)綿羊品種具有豐富的遺傳變異和較高的遺傳多樣性。呂慎金等[22]利用30個(gè)微衛(wèi)星標(biāo)記得出8個(gè)綿羊群體的平均雜合度為0.539-0.714。儲(chǔ)明星等[23]用5個(gè)微衛(wèi)星標(biāo)記揭示的STH平均雜合度為0.456 9,而本研究中6個(gè)微衛(wèi)星標(biāo)記得到STH的平均雜合度為0.671,表明本研究中所選擇的微衛(wèi)星座位具有較高的多態(tài)性。

    圖 1 6個(gè)綿羊品種前3個(gè)主成分的三維散點(diǎn)圖

    圖 2 6個(gè)綿羊品種間的UPGMA聚類(lèi)圖

    微衛(wèi)星DNA本身是選擇中性的,不受選擇的壓力,因此在一個(gè)理想群體中各等位基因在群體中的分布頻率應(yīng)該是穩(wěn)定的[24]。通過(guò)哈代溫伯格平衡檢驗(yàn)6個(gè)座位,發(fā)現(xiàn)MAF209、OarJMP29和OarFCB304三個(gè)座位在6個(gè)品種都處于平衡狀態(tài),ILSTS005座位在GSH偏離平衡,MCM527座位在ZT和GSH偏離平衡,OarAE129在UQ、HBR、STH和TAN等4個(gè)品種中偏離哈代溫伯格平衡。

    3.2 遺傳距離和聚類(lèi)分析

    該研究通過(guò)Arlequin軟件計(jì)算得到了Nei 氏標(biāo)準(zhǔn)遺傳距離,應(yīng)用MVSP 3.1進(jìn)行了主成分分析和UPGMA聚類(lèi)。通過(guò)這兩種方法,發(fā)現(xiàn)6個(gè)綿羊品種間的聚類(lèi)結(jié)果與遺傳距離基本一致。一般認(rèn)為群體分化時(shí)間越短,遺傳距離越小,HBR、UQ、STH和TAN等4個(gè)本地品種之間的遺傳距離較小,ZT與這4個(gè)品種之間的遺傳距離相對(duì)較大。而培育品種GSH與5個(gè)地方品種遺傳距離均較大,獨(dú)自聚為一類(lèi),聚類(lèi)結(jié)果與品種的形成史、分化及地理分布基本一致。如HBR與UQ遺傳距離最小,首先聚為一類(lèi),它們均來(lái)自?xún)?nèi)蒙古,生活環(huán)境極為相似,均是在當(dāng)?shù)亻L(zhǎng)期的自然和人工選擇情況下而逐漸育成的肉脂兼用型綿羊品種[25]。GSH是從1950年開(kāi)始,采用復(fù)雜育成雜交方法進(jìn)行培育,于1981年育成的培育品種[26],與其他5個(gè)地方綿羊品種遺傳距離最大,親緣關(guān)系較遠(yuǎn)。

    HBR、UQ、STH和TAN等4個(gè)品種均隸屬于蒙古羊系統(tǒng),起源于共同的祖先群體而親緣關(guān)系相對(duì)較近,有相關(guān)研究也證實(shí)了這一觀點(diǎn)。如仲濤等[27]采用21個(gè)微衛(wèi)星座位,運(yùn)用DA遺傳距離將9個(gè)地方綿羊品種進(jìn)行了分類(lèi),其中屬于蒙 古羊系的UQ、HBR、TAN和湖羊聚 為一類(lèi)。郭小雅等[28]分析了我國(guó)6個(gè)綿(山)羊品種的遺傳分化,進(jìn)一步證明了UQ和STH、TAN、湖羊、同羊有較近的親緣關(guān)系,共同歸屬于“蒙古羊”系。任春環(huán)等[29]分析了5個(gè)綿羊品種的遺傳多態(tài)性,結(jié)果表明TAN與STH的親緣關(guān)系較近。馬月輝等[1]采用30個(gè)微衛(wèi)星標(biāo)記,對(duì)12個(gè)中國(guó)地方綿羊品種和一個(gè)外來(lái)品種進(jìn)行遺傳多樣性研究,發(fā)現(xiàn)蘇尼特羊和UQ首先聚為一類(lèi),TAN和STH聚為一類(lèi),然后這4個(gè)品種作為一個(gè)分支聚為一類(lèi)。本研究中TAN與STH遺傳距離較小為0.076,STH與UQ遺傳距離最小為0.040,本研究中僅用了6個(gè)微衛(wèi)星位點(diǎn)分析6個(gè)綿羊品種遺傳多樣性,但與先前研究結(jié)果基本一致,說(shuō)明本研究采用的6個(gè)高度多態(tài)微衛(wèi)星在一定程度上能夠反映綿羊群體的遺傳結(jié)構(gòu)。

    昭通綿羊的起源有可能是云南本土或是由西藏的野生綿羊馴化而來(lái)[27]。蘭蓉等[30]對(duì)云南綿羊進(jìn)行了線粒體DNA和血液蛋白多態(tài)性研究,結(jié)果發(fā)現(xiàn)昭通綿羊和印度綿羊、尼泊爾綿羊關(guān)系較近,而與西藏綿羊關(guān)系較遠(yuǎn),但其起源還有待于深入研究。仲濤等[31]對(duì)中國(guó)部分綿羊品種的遺傳多樣性進(jìn)行了研究,表明昭通綿羊與其他藏系綿羊間遺傳分化較明顯,支持了昭通綿羊起源于云南本土的說(shuō)法。婁淵根等[32]認(rèn)為昭通綿羊因位于云貴高原,地理隔離阻礙了其與其他品種的基因交流,使其在分類(lèi)上自成一類(lèi)。本研究中昭通綿羊與其他5個(gè)綿羊品種遺傳距離較大,主成分分析和聚類(lèi)分析也得到相同的結(jié)果,研究結(jié)果支持昭通綿羊起源于云南本土的說(shuō)法。

    4 結(jié)論

    本研究選用的6個(gè)微衛(wèi)星為高度多態(tài)性位點(diǎn),在一定程度上可以反映綿羊群體的遺傳結(jié)構(gòu)。通過(guò)微衛(wèi)星標(biāo)記分析發(fā)現(xiàn)6個(gè)綿羊品種具有豐富的遺傳變異和較高的遺傳多樣性。品種間遺傳距離分析、主成分分析及UPGMA聚類(lèi),結(jié)果均一致說(shuō)明培育品種GSH與其他5個(gè)地方品種的遺傳分歧最大,HBR與UQ遺傳分歧最?。籗TH、TAN、HBR和UQ共同歸屬于蒙古羊系,親緣關(guān)系相對(duì)較近,進(jìn)一步證實(shí)了STH與TAN親緣關(guān)系較近的說(shuō)法,各綿羊品種間的遺傳關(guān)系與其形成史、分化及地理分布基本一致,昭通綿羊與其他綿羊品種間的遺傳距離較大,支持了其起源于云南本地的觀點(diǎn)。

    致謝:感謝中國(guó)農(nóng)業(yè)科學(xué)院北京畜牧獸醫(yī)研究所馬月輝研究員、四川農(nóng)業(yè)大學(xué)動(dòng)物科技學(xué)院仲濤副研究員在樣品采集中給予的幫助。

    [1]馬月輝, 呂慎金, 侯冠或, 等. 中國(guó)部分綿羊品種遺傳多樣性研究[C]//中國(guó)畜牧獸醫(yī)學(xué)會(huì)養(yǎng)羊?qū)W分會(huì)全國(guó)養(yǎng)羊生產(chǎn)與學(xué)術(shù)研究會(huì)議論文集. 中國(guó)草食動(dòng)物雜志出版社, 2005.

    [2]岳遠(yuǎn)瑞, 袁曉航, 周路, 等. 新疆7個(gè)綿羊品種MHC區(qū)段微衛(wèi)星遺傳多樣性的分析[J]. 石河子大學(xué)學(xué)報(bào):自然科學(xué)版, 2014, 32(4):438-443.

    [3]Botstein D, White RL, Skolnick M, et al. Construction of a genetic linkage map in man using restriction fragment length polymorphisms[J]. Am J Hum Genet, 1980, 32:314-331.

    [4]Lumsden JM, Loard EA, Montgomery GW. Characterisation and linkage mapping of ten microsatellite markers derived from a sheep x hamster cell hybrid[J]. Animal Genetics, 1996, 27:203-206.

    [5]Chen HJ, Yue YS, Fan XZ, et al, Study comparatively on genetic distance and genus analysis measures in Shandong indigenous chicken breeds[J]. Journal of Farming and Veterinarian, 2004, 35(1):33-36.

    [6]Sun W, Chang H, Tsunoda K, et al. Analysis of geographic and pairwise genetic distances among sheep populations[J]. Biochem Genet, 2010, 48(5-6):376-84.

    [7]Agaviezor BO, Peters SO, Adefenwa MA, et al, Morphological and microsatellite DNA diversity of Nigerian indigenous sheep[J]. Anim Sci Biotechnol, 2012, (1):38.

    [8]Yilmaz O, Cemal I, Karaca O. Genetic diversity in nine native Turkish sheep breeds based on microsatellite analysis[J]. Anim Genet, 2014, 45(4):604-608.

    [9]Leroy G, Danchin-Burge C, Palhière I, et al. How do introgression events shape the partitioning of diversity among breeds:a case study in sheep[J]. Genet Sel Evol, 2015, 47(1):48.

    [10]?urkovi? M, Ramljak J, Ivankovi? S, et al. The genetic diversity and structure of 18 sheep breeds exposed to isolation and selection[J]. J Anim Breed Genet, 2016, 133(1):71-80.

    [11]曾盛誠(chéng). 利用20個(gè)微衛(wèi)星標(biāo)記分析國(guó)內(nèi)外11個(gè)雞群體的遺傳多樣性[D]. 北京:中國(guó)農(nóng)業(yè)科學(xué)院, 200 9.

    [12]Park SDE. Trypanotolerance in West African Cattle and thepopulation genetic effects of selection[D]. Dublin:University of Dublin, 2001.

    [13]Excoffier L, Lischer HEL. Arlequin 3.5:A new series of programs to perform population genetics analyses under Linux and Windows[J]. Molecular Ecology Resources, 2010, 10:564-567.

    [14]Raymond M, Rousset F. GENEPOP(Version 1. 2):population genetics software for exact tests and ecumenicism[J]. Journal of Heredity, 1995, 86(3):248-249.

    [15]Goudet J. STAT(Version 1. 2):A computer program to Calculate F-Statistic[J]. J Hered, 1995, 86:485-486.

    [16]Nei M, Li WH. Mathematical model for studying genetic variation in terms of restriction endonucleases[J]. Proc Natl Acad Sci USA, 1979, 76:52-69.

    [17]Nei M, Li WH. Analysis of gene diversity in subdivided populations[J]. Proc Nat l Acad Sci USA, 1973, 70:3321-3323.

    [18]常爽. 9個(gè)綿羊群體遺傳多樣性的微衛(wèi)星分析[D]. 鄭州:河南農(nóng)業(yè)大學(xué), 2010.

    [19]Botstein D, White RL, Skolnick MH, Davis R. Construction of a genetic linkage map in man using restriction fragment length polymorphisms[J]. Am J Hum Genet, 1980, 32:314-331.

    [20]Ma YH, R ao SQ, Lu SJ, et al. Phylogeography and origin of sheep breeds in Northern China[J]. Conservation Genetics, 2006, 7(1):117-127.

    [21]儲(chǔ)明星, 程金華, 過(guò)緯. 微衛(wèi)星標(biāo)記OarAE101和BM1329在五個(gè)綿羊品種中的初步研究[J]. 遺傳學(xué)報(bào), 2001, 28(6):510- 517.

    [22]呂慎金, 楊燕, 候冠玉, 等. 運(yùn)用微衛(wèi)星標(biāo)記對(duì)中國(guó)地方綿羊品種的遺傳多樣性分析[J]. 畜牧獸醫(yī)學(xué)報(bào), 2008, 7:858-865.

    [23]儲(chǔ)明星, 王吉振, 王愛(ài)國(guó), 等. 小尾寒羊五個(gè)微衛(wèi)星基因座遺傳多態(tài)性研究[J]. 遺傳學(xué)報(bào), 2002, 29(6):502-506.

    [24]李延璐. 7個(gè)綿羊群體15個(gè)微衛(wèi)星DNA座位遺傳多樣性分析[D]. 合肥:安徽農(nóng)業(yè)大學(xué), 2009.

    [25]趙有璋. 中國(guó)養(yǎng)羊?qū)W[M]. 北京:中國(guó)農(nóng)業(yè)出版社, 2011.

    [26]杜立新. 中國(guó)畜禽遺傳資源志. 羊志[M]. 北京:中國(guó)農(nóng)業(yè)出版社, 2011.

    [27]仲濤, 馬月輝, 關(guān)偉軍, 等. 10 個(gè)綿羊品種的微衛(wèi)星 DNA 多態(tài)性研究[J]. 畜牧獸醫(yī)學(xué), 2008, 39(5):555-561.

    [28]郭小雅, 楊章平, 常洪, 等. 13我國(guó)6個(gè)綿(山)羊群體遺傳分化的微衛(wèi)星分析[J]. 畜牧與獸醫(yī), 2007, 39(5):16-20.

    [29]任春環(huán), 張子軍, 吳建平, 等. 5個(gè)綿羊群體微衛(wèi)星標(biāo)記的遺傳多態(tài)性分析[J]. 安徽農(nóng)業(yè)大學(xué)學(xué)報(bào), 2010, 37(3):436-439 .

    [30]蘭蓉, 洪瓊花, 張俊, 等. 云南不同地區(qū)綿羊血液蛋白多態(tài)性研究[J]. 草食家畜, 2002, 3(1):43-47.

    [31]仲濤. 中國(guó)部分綿羊品種的遺傳多樣性研究[D]. 北京:中國(guó)農(nóng)業(yè)科學(xué)院, 2007.

    [32]婁淵根, 劉丑生, 劉建學(xué), 等. 9個(gè)綿羊品種微衛(wèi)星DNA遺傳多樣性分析[J]. 中國(guó)畜牧獸醫(yī), 2011, 38(10):106-112.

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

    The Genetic Diversity of Six Chinese Sheep Breeds by Six Microsatellite Markers

    CHEN Li-peng HUANG Yong-fu ZHAO Yong-ju E Guang-xin NA Ri-su
    (College of Animal Science and Technology,Southwest University,Chongqing 400716)

    Six microsatellite loci were used to analyze the genetic diversity of 280 sheep from 6 Chinese bre eds of Ujumqin sheep,Smalltailed Han sheep,Tan sheep,Zhaotong sheep,Hulun Buir sheep,and Gansu Alpine Merino. We calculated the polymorphism information content,heterozygosity and genetic distance,analyzed the principal components and clustered UPGMA from these 6 sheep breeds. Altogether,101 alleles were found,the mean heterozygosity was 0.599-0.691,and mean polymorphism information content was 0.609-0.680. The Gansu Alpine Merino sheep had the most genetic divergence compared to other sheep breeds,while the Hulun Buir and Ujumqin ones were the least in genetic divergence. The above results showed that these 6 breeds had relatively high genetic diversity,and their genetic relationships were consistent with their formation history,differentiation and geographical distribution.

    microsatellite;genetic diversity;sheep

    10.13560/j.cnki.biotech.bull.1985.2016.05.012

    2015-07-28

    中央高?;究蒲袠I(yè)務(wù)費(fèi)專(zhuān)項(xiàng)資金(SWU114023)

    陳李鵬,男,碩士,研究方向:動(dòng)物遺傳;E-mail:clpshanxi@126.com

    黃勇富,男,研究員,研究方向:動(dòng)物遺傳多樣性、遺傳資源評(píng)價(jià)與創(chuàng)新利用;E-mail:hyf65@163.com

    猜你喜歡
    微衛(wèi)星雜合綿羊
    梅花鹿基因組微衛(wèi)星分布特征研究
    甘藍(lán)型油菜隱性上位互作核不育系統(tǒng)不育系材料選育中常見(jiàn)的育性分離及基因型判斷
    種子(2021年3期)2021-04-12 01:42:22
    數(shù)綿羊
    數(shù)綿羊
    奔跑的綿羊
    林麝全基因組微衛(wèi)星分布規(guī)律研究
    從翻譯到文化雜合——“譯創(chuàng)”理論的虛涵數(shù)意
    巧計(jì)得綿羊
    基于轉(zhuǎn)錄組測(cè)序的波紋巴非蛤微衛(wèi)星標(biāo)記研究
    中國(guó)“一箭雙星”成功將“遙感衛(wèi)星二十一號(hào)”與“天拓二號(hào)視頻微衛(wèi)星”發(fā)射升空
    河北遙感(2014年3期)2014-07-10 13:16:48
    国产精品野战在线观看| 99热只有精品国产| 久久精品91蜜桃| 啦啦啦韩国在线观看视频| 蜜桃久久精品国产亚洲av| 少妇人妻精品综合一区二区 | av视频在线观看入口| 日韩欧美在线二视频| 别揉我奶头 嗯啊视频| 国产伦一二天堂av在线观看| 久久久久久久久久成人| 99精品久久久久人妻精品| 69av精品久久久久久| 69人妻影院| 深爱激情五月婷婷| 天堂动漫精品| 中文字幕人妻熟人妻熟丝袜美| 亚洲人成网站高清观看| 看十八女毛片水多多多| 欧美bdsm另类| 亚洲精华国产精华精| 免费看日本二区| 欧美高清成人免费视频www| 1000部很黄的大片| 欧美一区二区国产精品久久精品| 给我免费播放毛片高清在线观看| 亚洲图色成人| 久久精品91蜜桃| 国产精品,欧美在线| 久久午夜福利片| 亚洲男人的天堂狠狠| 欧美激情国产日韩精品一区| 少妇被粗大猛烈的视频| 中国美女看黄片| 欧美日韩亚洲国产一区二区在线观看| 老熟妇乱子伦视频在线观看| 国产精品免费一区二区三区在线| 国产淫片久久久久久久久| 亚洲在线观看片| 一区二区三区激情视频| 亚洲精品粉嫩美女一区| 国产欧美日韩一区二区精品| 午夜免费男女啪啪视频观看 | 99热网站在线观看| 在线观看午夜福利视频| 啦啦啦韩国在线观看视频| 国产精品一区二区性色av| 欧美黑人欧美精品刺激| bbb黄色大片| 国产激情偷乱视频一区二区| 22中文网久久字幕| 免费观看人在逋| 亚洲一区高清亚洲精品| 99久久九九国产精品国产免费| eeuss影院久久| 两个人视频免费观看高清| 99久国产av精品| 一个人看视频在线观看www免费| 99国产极品粉嫩在线观看| 日韩强制内射视频| 国产亚洲欧美98| 91狼人影院| 久久人人精品亚洲av| av在线天堂中文字幕| 午夜激情欧美在线| 亚洲最大成人中文| 免费电影在线观看免费观看| 精品人妻熟女av久视频| 小蜜桃在线观看免费完整版高清| 国产精品无大码| 男女之事视频高清在线观看| 精品久久久久久久久久久久久| 亚洲国产精品合色在线| 国产高清激情床上av| 国产精品久久久久久av不卡| 人妻丰满熟妇av一区二区三区| 成人av一区二区三区在线看| 国产精品女同一区二区软件 | 一进一出好大好爽视频| 高清毛片免费观看视频网站| a级一级毛片免费在线观看| 国产爱豆传媒在线观看| 制服丝袜大香蕉在线| 在线免费观看不下载黄p国产 | 亚洲av熟女| 亚洲av不卡在线观看| 黄色一级大片看看| 国产精品人妻久久久影院| 真人一进一出gif抽搐免费| 免费在线观看日本一区| 亚洲精品在线观看二区| 亚洲av第一区精品v没综合| 如何舔出高潮| 成年女人看的毛片在线观看| 狂野欧美激情性xxxx在线观看| 日本撒尿小便嘘嘘汇集6| 国产成人一区二区在线| 99热网站在线观看| 狂野欧美白嫩少妇大欣赏| 色综合婷婷激情| 国产三级在线视频| 亚洲精华国产精华精| 真实男女啪啪啪动态图| 国产高清视频在线播放一区| 美女xxoo啪啪120秒动态图| 女人被狂操c到高潮| 日韩国内少妇激情av| 久久亚洲精品不卡| av女优亚洲男人天堂| 波多野结衣巨乳人妻| 最好的美女福利视频网| 69人妻影院| 色尼玛亚洲综合影院| 在线天堂最新版资源| 亚洲熟妇中文字幕五十中出| 性欧美人与动物交配| 韩国av一区二区三区四区| 午夜福利在线在线| 国产精品久久久久久av不卡| 国产精品电影一区二区三区| 久久国产乱子免费精品| 春色校园在线视频观看| 可以在线观看毛片的网站| 婷婷精品国产亚洲av在线| 我要搜黄色片| 欧洲精品卡2卡3卡4卡5卡区| 小说图片视频综合网站| 88av欧美| 三级毛片av免费| 欧美绝顶高潮抽搐喷水| 一a级毛片在线观看| 中国美女看黄片| 人妻丰满熟妇av一区二区三区| 波多野结衣高清无吗| 日本爱情动作片www.在线观看 | 欧美精品啪啪一区二区三区| 男女做爰动态图高潮gif福利片| 久久久久九九精品影院| 国产色婷婷99| 嫩草影视91久久| 国产伦人伦偷精品视频| 久久精品国产亚洲av天美| 欧美三级亚洲精品| 国国产精品蜜臀av免费| 国产单亲对白刺激| 国产精品1区2区在线观看.| 99久久中文字幕三级久久日本| 麻豆成人午夜福利视频| 欧美日韩国产亚洲二区| 久久精品国产自在天天线| 日韩欧美精品免费久久| 大又大粗又爽又黄少妇毛片口| 国产av在哪里看| 成人二区视频| 日韩精品中文字幕看吧| 国产91精品成人一区二区三区| av女优亚洲男人天堂| 欧美日韩国产亚洲二区| 国产精品av视频在线免费观看| 久久国内精品自在自线图片| 99久久精品一区二区三区| 亚洲乱码一区二区免费版| 色播亚洲综合网| 又黄又爽又刺激的免费视频.| 精品人妻偷拍中文字幕| 精品免费久久久久久久清纯| 久久精品国产鲁丝片午夜精品 | 99热这里只有是精品在线观看| 国产精品国产高清国产av| 亚洲一级一片aⅴ在线观看| 久久午夜亚洲精品久久| 在线播放无遮挡| 男女做爰动态图高潮gif福利片| 99在线视频只有这里精品首页| 久久久久免费精品人妻一区二区| 在线观看一区二区三区| 成年人黄色毛片网站| netflix在线观看网站| 久久精品夜夜夜夜夜久久蜜豆| 在现免费观看毛片| 小蜜桃在线观看免费完整版高清| 18禁黄网站禁片免费观看直播| a在线观看视频网站| 中文字幕高清在线视频| 久久精品久久久久久噜噜老黄 | 日韩欧美精品v在线| 久久久久国内视频| 少妇的逼水好多| 精品久久久久久久久亚洲 | 国产男靠女视频免费网站| 中亚洲国语对白在线视频| 免费人成在线观看视频色| 91av网一区二区| 日本一二三区视频观看| 国产精品不卡视频一区二区| 久久精品国产亚洲av香蕉五月| 一本一本综合久久| 国产伦在线观看视频一区| 高清日韩中文字幕在线| 亚洲人成伊人成综合网2020| 成年版毛片免费区| 亚洲熟妇熟女久久| 免费看日本二区| 日本黄大片高清| 校园人妻丝袜中文字幕| 精品久久久久久久久久免费视频| 亚洲乱码一区二区免费版| 桃红色精品国产亚洲av| 91狼人影院| 97超视频在线观看视频| 精品无人区乱码1区二区| 一个人看的www免费观看视频| 麻豆国产97在线/欧美| 校园春色视频在线观看| 国产免费一级a男人的天堂| 日本免费a在线| 亚洲av成人精品一区久久| 亚洲自拍偷在线| 春色校园在线视频观看| 看黄色毛片网站| 3wmmmm亚洲av在线观看| 国产真实乱freesex| 欧美激情国产日韩精品一区| 99久久中文字幕三级久久日本| 国产精品一区www在线观看 | 午夜精品久久久久久毛片777| 看免费成人av毛片| 国产成人福利小说| 亚洲久久久久久中文字幕| 一区二区三区高清视频在线| 两个人视频免费观看高清| 欧美潮喷喷水| 一区福利在线观看| 亚洲人与动物交配视频| 看十八女毛片水多多多| a级毛片a级免费在线| 尤物成人国产欧美一区二区三区| 波多野结衣高清无吗| 久久午夜福利片| 99热精品在线国产| 男女下面进入的视频免费午夜| 乱系列少妇在线播放| 人妻制服诱惑在线中文字幕| 美女黄网站色视频| 桃色一区二区三区在线观看| 国产精品久久久久久久电影| 久久精品综合一区二区三区| 亚洲图色成人| 亚洲乱码一区二区免费版| 亚洲最大成人手机在线| 小蜜桃在线观看免费完整版高清| 草草在线视频免费看| 欧美色欧美亚洲另类二区| av.在线天堂| 91在线精品国自产拍蜜月| 国产伦精品一区二区三区四那| 久久精品夜夜夜夜夜久久蜜豆| 国产成人av教育| 国产视频一区二区在线看| 一级毛片久久久久久久久女| 久久久精品大字幕| 午夜a级毛片| 韩国av在线不卡| 亚洲中文字幕一区二区三区有码在线看| 日韩欧美免费精品| 国产色婷婷99| 婷婷六月久久综合丁香| 国产中年淑女户外野战色| 国产亚洲精品久久久久久毛片| 听说在线观看完整版免费高清| 亚洲欧美日韩高清在线视频| 婷婷六月久久综合丁香| 国产亚洲欧美98| 一区二区三区高清视频在线| 久久久久国内视频| 日本熟妇午夜| 亚洲精品国产成人久久av| 真实男女啪啪啪动态图| 精品人妻视频免费看| 99精品久久久久人妻精品| 免费av不卡在线播放| av在线天堂中文字幕| 成人一区二区视频在线观看| 丰满的人妻完整版| 成年女人毛片免费观看观看9| 日韩人妻高清精品专区| 嫁个100分男人电影在线观看| 欧美精品啪啪一区二区三区| 国国产精品蜜臀av免费| 亚洲人成网站高清观看| 日韩中字成人| 国产视频内射| 午夜福利高清视频| 一本久久中文字幕| 亚洲av二区三区四区| 最近视频中文字幕2019在线8| 十八禁国产超污无遮挡网站| 国产高清激情床上av| 麻豆av噜噜一区二区三区| ponron亚洲| 给我免费播放毛片高清在线观看| 一个人看的www免费观看视频| 俄罗斯特黄特色一大片| 亚洲av中文字字幕乱码综合| 国产高清视频在线观看网站| 精品久久久久久久久久久久久| 欧美日韩精品成人综合77777| 韩国av在线不卡| 免费在线观看影片大全网站| 午夜老司机福利剧场| 亚洲最大成人中文| 成人国产综合亚洲| 偷拍熟女少妇极品色| 国产在线男女| 国产男人的电影天堂91| 午夜久久久久精精品| 精品久久久久久久久亚洲 | 婷婷色综合大香蕉| 麻豆av噜噜一区二区三区| 又粗又爽又猛毛片免费看| 国产成人影院久久av| 看十八女毛片水多多多| 精品免费久久久久久久清纯| 久久人妻av系列| 热99在线观看视频| 国产精品久久久久久久久免| 欧美日本亚洲视频在线播放| 热99re8久久精品国产| 悠悠久久av| 日本欧美国产在线视频| 在线观看免费视频日本深夜| 亚洲精品影视一区二区三区av| 日本熟妇午夜| 日韩欧美免费精品| 91在线观看av| 国国产精品蜜臀av免费| 日韩欧美在线二视频| 亚洲精品成人久久久久久| 51国产日韩欧美| 12—13女人毛片做爰片一| 国产精华一区二区三区| 欧美日韩综合久久久久久 | 色综合婷婷激情| 亚洲精品粉嫩美女一区| 精品一区二区三区av网在线观看| 国产91精品成人一区二区三区| 国内久久婷婷六月综合欲色啪| 国产一区二区在线观看日韩| 一级黄片播放器| 老熟妇仑乱视频hdxx| 男人和女人高潮做爰伦理| 日本黄色视频三级网站网址| 一本久久中文字幕| 中文亚洲av片在线观看爽| 一边摸一边抽搐一进一小说| 变态另类成人亚洲欧美熟女| 性插视频无遮挡在线免费观看| 久久久久国产精品人妻aⅴ院| 嫩草影院新地址| 特级一级黄色大片| 免费看美女性在线毛片视频| 国产伦在线观看视频一区| 18禁在线播放成人免费| 精品久久久久久久人妻蜜臀av| 久久精品影院6| 欧美一级a爱片免费观看看| 免费电影在线观看免费观看| 乱码一卡2卡4卡精品| 久久精品国产清高在天天线| 国产精品久久电影中文字幕| 少妇丰满av| 色av中文字幕| 国产 一区精品| 久久99热6这里只有精品| 国产三级在线视频| 99热6这里只有精品| 国产免费一级a男人的天堂| 真人一进一出gif抽搐免费| 在线a可以看的网站| 综合色av麻豆| 在线播放无遮挡| 国产精品精品国产色婷婷| 国产亚洲精品av在线| 亚洲性夜色夜夜综合| 99精品在免费线老司机午夜| 色哟哟·www| 亚洲在线自拍视频| 亚洲一级一片aⅴ在线观看| 日本在线视频免费播放| 香蕉av资源在线| 中文字幕人妻熟人妻熟丝袜美| 我的女老师完整版在线观看| 真实男女啪啪啪动态图| 亚洲18禁久久av| 国内精品久久久久精免费| 床上黄色一级片| 国产在线男女| 色综合站精品国产| 日韩欧美在线乱码| 亚洲av一区综合| 成人高潮视频无遮挡免费网站| 69人妻影院| 白带黄色成豆腐渣| 亚洲精品久久国产高清桃花| 日本黄色视频三级网站网址| 在现免费观看毛片| 国产爱豆传媒在线观看| 非洲黑人性xxxx精品又粗又长| 欧美+亚洲+日韩+国产| 性插视频无遮挡在线免费观看| 麻豆成人av在线观看| 午夜老司机福利剧场| 亚洲第一区二区三区不卡| 一本久久中文字幕| 我的老师免费观看完整版| 国产麻豆成人av免费视频| 99国产极品粉嫩在线观看| 欧美3d第一页| 男人舔女人下体高潮全视频| 免费黄网站久久成人精品| 久久久成人免费电影| 九色国产91popny在线| 长腿黑丝高跟| 少妇的逼水好多| 免费看美女性在线毛片视频| 亚洲国产精品合色在线| 国模一区二区三区四区视频| 国产视频一区二区在线看| 男人和女人高潮做爰伦理| 在线观看av片永久免费下载| 亚洲人与动物交配视频| 免费在线观看影片大全网站| 在线免费观看不下载黄p国产 | 别揉我奶头 嗯啊视频| 国产精品免费一区二区三区在线| 免费观看精品视频网站| 狂野欧美激情性xxxx在线观看| 伦精品一区二区三区| 丝袜美腿在线中文| 久久香蕉精品热| 99精品久久久久人妻精品| 久久热精品热| 精品无人区乱码1区二区| 美女大奶头视频| 99久国产av精品| 久久99热这里只有精品18| 久久久久久久亚洲中文字幕| 色哟哟·www| 欧美日韩国产亚洲二区| 99久久无色码亚洲精品果冻| 久久精品91蜜桃| 熟女人妻精品中文字幕| 人人妻人人澡欧美一区二区| 极品教师在线视频| 高清日韩中文字幕在线| 日韩一区二区视频免费看| 欧美一区二区亚洲| 国产精品一区二区三区四区免费观看 | 99热网站在线观看| 99久久久亚洲精品蜜臀av| 男人的好看免费观看在线视频| 国产女主播在线喷水免费视频网站 | 欧美3d第一页| 国产精品久久久久久精品电影| 黄色欧美视频在线观看| 亚洲狠狠婷婷综合久久图片| 国产白丝娇喘喷水9色精品| 成年人黄色毛片网站| 亚洲av中文av极速乱 | 美女黄网站色视频| 久久久久久久久中文| 岛国在线免费视频观看| 国产极品精品免费视频能看的| 非洲黑人性xxxx精品又粗又长| 搡老岳熟女国产| 一本精品99久久精品77| 露出奶头的视频| 看免费成人av毛片| 嫁个100分男人电影在线观看| 中文字幕av成人在线电影| 欧美中文日本在线观看视频| 国产免费av片在线观看野外av| 精品国内亚洲2022精品成人| 免费看美女性在线毛片视频| 小说图片视频综合网站| 无人区码免费观看不卡| 一区福利在线观看| 午夜精品一区二区三区免费看| 国产精品无大码| 亚洲国产精品久久男人天堂| 国产男人的电影天堂91| 老司机深夜福利视频在线观看| 性色avwww在线观看| 一进一出抽搐gif免费好疼| 国产亚洲欧美98| 黄片wwwwww| 黄色欧美视频在线观看| 国产黄色小视频在线观看| 国产熟女欧美一区二区| 久久久久性生活片| 国产精品av视频在线免费观看| 联通29元200g的流量卡| 欧美成人性av电影在线观看| 国内精品美女久久久久久| 夜夜夜夜夜久久久久| 精品一区二区免费观看| 久久久久久久久中文| 啪啪无遮挡十八禁网站| av在线蜜桃| 久久草成人影院| 国产精品国产三级国产av玫瑰| 亚洲成av人片在线播放无| 最近最新免费中文字幕在线| 亚洲专区中文字幕在线| 蜜桃久久精品国产亚洲av| 床上黄色一级片| 又黄又爽又免费观看的视频| 1024手机看黄色片| 一级a爱片免费观看的视频| 日本一本二区三区精品| 一个人免费在线观看电影| 欧美三级亚洲精品| 国产单亲对白刺激| 亚洲18禁久久av| 两性午夜刺激爽爽歪歪视频在线观看| 一区二区三区免费毛片| 日韩精品青青久久久久久| 内地一区二区视频在线| 色吧在线观看| 亚洲无线观看免费| 国产亚洲av嫩草精品影院| 久久人人精品亚洲av| 国产 一区 欧美 日韩| 国产精品爽爽va在线观看网站| 国内精品一区二区在线观看| 一区二区三区四区激情视频 | 男女那种视频在线观看| av在线天堂中文字幕| 亚洲中文字幕日韩| 久久久国产成人精品二区| 美女黄网站色视频| 特级一级黄色大片| 国产欧美日韩一区二区精品| 少妇丰满av| 国产激情偷乱视频一区二区| 亚洲精品影视一区二区三区av| 又粗又爽又猛毛片免费看| 色尼玛亚洲综合影院| or卡值多少钱| 国产精品久久久久久久电影| 欧美中文日本在线观看视频| 欧美精品国产亚洲| 色综合色国产| 精品人妻熟女av久视频| 黄色女人牲交| 尾随美女入室| 日本一二三区视频观看| 日本色播在线视频| 成人欧美大片| 一区二区三区免费毛片| 欧美日本亚洲视频在线播放| 成人美女网站在线观看视频| 国产高清有码在线观看视频| 91久久精品电影网| 国产精品久久电影中文字幕| 露出奶头的视频| 亚洲中文字幕日韩| 国产精品久久久久久久久免| 久久久久久久久久黄片| 成人永久免费在线观看视频| 国模一区二区三区四区视频| 午夜福利在线在线| 亚洲人成网站高清观看| 成人午夜高清在线视频| 91午夜精品亚洲一区二区三区 | 自拍偷自拍亚洲精品老妇| 欧美高清性xxxxhd video| 露出奶头的视频| 麻豆成人av在线观看| 国产真实乱freesex| 2021天堂中文幕一二区在线观| 亚洲第一电影网av| 欧美性猛交黑人性爽| 一级黄色大片毛片| 桃色一区二区三区在线观看| 亚洲av免费在线观看| 成人高潮视频无遮挡免费网站| 18禁黄网站禁片午夜丰满| 久久精品国产清高在天天线| 日韩国内少妇激情av| 国产精品美女特级片免费视频播放器| 免费av毛片视频| 欧美+日韩+精品| 国产人妻一区二区三区在| 精品一区二区免费观看| 熟妇人妻久久中文字幕3abv| 精品人妻熟女av久视频| 国产蜜桃级精品一区二区三区| 久久久久免费精品人妻一区二区| 久久久色成人| 国产伦精品一区二区三区视频9| 色5月婷婷丁香| 日韩亚洲欧美综合| 九九久久精品国产亚洲av麻豆| 国产伦一二天堂av在线观看| 成人性生交大片免费视频hd| 久久久久久伊人网av| 免费看光身美女| 97热精品久久久久久| 亚洲国产精品sss在线观看| 亚洲av中文字字幕乱码综合| 亚洲欧美清纯卡通| 成人三级黄色视频| 成人精品一区二区免费| 一边摸一边抽搐一进一小说|