999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

埃及品系尼羅羅非魚不同選育世代間遺傳潛力分析

2014-07-18 09:56:23王騰肖煒李大宇鄒芝英祝璟琳樂貽榮楊弘
江蘇農(nóng)業(yè)科學(xué) 2014年2期

王騰+肖煒+李大宇+鄒芝英+祝璟琳+樂貽榮+楊弘

摘要:選用25個(gè)微衛(wèi)星基因座對(duì)連續(xù)3個(gè)選育世代的埃及品系尼羅羅非魚進(jìn)行檢測(cè),通過分析世代間的遺傳差異來探討埃及品系尼羅羅非魚的保種現(xiàn)狀和遺傳潛力。結(jié)果表明,3個(gè)世代群體在25個(gè)微衛(wèi)星基因座中共檢測(cè)出105個(gè)等位基因,不同基因座上的等位基因數(shù)為2~7,3個(gè)世代群體的平均有效等位基因數(shù)(Ae)分別為2.59、2.52、230個(gè),平均期望雜合度(He)分別為0.547、0.545、0.518,平均觀測(cè)雜合度(Ho)分別為0.479、0.457、0.451,平均多態(tài)信息含量(PIC)分別為0.487、0.480、0.458,3個(gè)世代群體遺傳多樣性水平較高,世代間各遺傳參數(shù)隨世代數(shù)增長(zhǎng)呈略微下降趨勢(shì),但差異不顯著(P>0.05);3個(gè)世代間遺傳分化指數(shù)(Fst)在0.017~0.048之間,均小于0.05,遺傳分化程度相對(duì)較低;通過分子方差對(duì)遺傳變異來源分析表明,經(jīng)過選育群體世代間的遺傳分化較小,世代間變異僅占總變異的334%,分化主要來源于世代內(nèi)個(gè)體間的遺傳差異。目前,采用保種選育方法對(duì)埃及品系尼羅羅非魚世代間的遺傳信息較為穩(wěn)定,具有較高的遺傳潛力。

關(guān)鍵詞:埃及品系;尼羅羅非魚;選育世代;微衛(wèi)星;遺傳多樣性;遺傳潛力

中圖分類號(hào): S917.4文獻(xiàn)標(biāo)志碼: A文章編號(hào):1002-1302(2014)02-0184-04

收稿日期:2013-07-29

基金項(xiàng)目:現(xiàn)代農(nóng)業(yè)產(chǎn)業(yè)技術(shù)體系專項(xiàng)(編號(hào):CARS-49);中央級(jí)公益性科研院所基本科研業(yè)務(wù)費(fèi)專項(xiàng)(編號(hào):2013JBFM01);現(xiàn)代農(nóng)業(yè)人才支撐計(jì)劃(編號(hào):2130106)。

作者簡(jiǎn)介:王騰(1990—),女,山東臨沂人,碩士,研究方向?yàn)轸~類遺傳育種。E-mail:wangteng0416@163.com。

通信作者:楊弘(1966—),男,江蘇無錫人,研究員,主要從事魚類育種工作。E-mail:yangh@ffrc.cn。尼羅羅非魚(Orechromis niloticus)原產(chǎn)于非洲,具有生長(zhǎng)速度快、食性雜及抗病力強(qiáng)等優(yōu)點(diǎn),是聯(lián)合國糧農(nóng)組織(FAO)向世界各國推薦的主要養(yǎng)殖品種之一。自1978年引進(jìn)我國以來,羅非魚養(yǎng)殖在我國得到迅速發(fā)展,現(xiàn)已成為我國淡水養(yǎng)殖的重要品種。優(yōu)良的苗種是保證羅非魚成魚品質(zhì)的關(guān)鍵,因此,良種的選擇和培育是保障優(yōu)良苗種供給、實(shí)現(xiàn)我國羅非魚產(chǎn)業(yè)穩(wěn)定發(fā)展的決定性因素。為此,1999年,中國水產(chǎn)科學(xué)研究院淡水漁業(yè)研究中心從埃及農(nóng)業(yè)和農(nóng)墾部水產(chǎn)研究中心實(shí)驗(yàn)室引進(jìn)尼羅羅非魚基礎(chǔ)群體5 000尾,在南泉養(yǎng)殖基地采用閉鎖群體繼代選育法,以生長(zhǎng)速度和尾鰭條紋為主要選擇目標(biāo),2年1個(gè)世代,選擇總強(qiáng)度為30%左右。經(jīng)過多年選育,該群體生長(zhǎng)、經(jīng)濟(jì)性狀進(jìn)一步得到提高,現(xiàn)已成為1個(gè)新的優(yōu)良品系[1]。目前,群體選育法已經(jīng)較多地用于羅非魚的選育工作,往往重視經(jīng)濟(jì)性狀的選擇,在分子水平上對(duì)遺傳變化的跟蹤研究卻較少,忽略了群體遺傳型變化可能會(huì)導(dǎo)致遺傳多樣性的降低,從而降低遺傳效應(yīng),造成遺傳瓶頸效應(yīng)。

近年來,分子標(biāo)記技術(shù)已經(jīng)廣泛應(yīng)用于水產(chǎn)動(dòng)物的親緣關(guān)系鑒定、物種多樣性檢測(cè)、種質(zhì)資源的保護(hù)以及基因連鎖圖譜的構(gòu)建和基因定位等領(lǐng)域[2-5],在開展品種選育時(shí),對(duì)群體遺傳結(jié)構(gòu)影響的研究也有一些成果。Romana-Eguia等利用微衛(wèi)星標(biāo)記對(duì)尼羅羅非魚不同選育世代的研究、Patricia等利用RAPD技術(shù)對(duì)不同世代凡納濱對(duì)蝦群體遺傳多樣性的分析,結(jié)果均顯示選育世代遺傳多樣性有逐代降低的趨勢(shì)[6-7];頡曉勇等對(duì)吉富品系尼羅羅非魚選育群體遺傳變異分析發(fā)現(xiàn),隨著選育的進(jìn)行,群體遺傳參數(shù)逐漸降低,群體呈現(xiàn)出純化趨勢(shì)[8]。本研究擬采用微衛(wèi)星標(biāo)記對(duì)埃及品系尼羅羅非魚3個(gè)世代保種群遺傳變異進(jìn)行檢測(cè),初步評(píng)價(jià)埃及品系尼羅羅非魚在現(xiàn)有選育方式下的保種現(xiàn)狀及遺傳潛力,了解尼羅羅非魚在保種選育過程中遺傳型的變化,為制定更加合理有效的保種方案提供理論依據(jù),為其他品種的引種和保種工作提供參考。

1材料與方法

1.1試驗(yàn)材料

1.2方法

1.2.1基因組DNA的制備采用常規(guī)酚-氯仿法制備基因組DNA,并通過瓊脂糖凝膠電泳和紫外分光光度計(jì)檢測(cè)其濃度和純度,稀釋至100 ng/μL,作為PCR擴(kuò)增的模板。

1.2.2微衛(wèi)星引物的設(shè)計(jì)本研究中使用的引物均從已發(fā)表的尼羅羅非魚遺傳連鎖圖譜中選取[9],按照互不連鎖的原則,從不同的連鎖群中選擇具有多態(tài)性的25個(gè)微衛(wèi)星基因座,相關(guān)信息見表1。引物由上海申能博彩公司合成。

1.2.3PCR擴(kuò)增及電泳檢測(cè)反應(yīng)體系:10×Buffer(Mg2+ Free)2.5 μL,MgCl2(25 mmol/L)1.5 μL,微衛(wèi)星上下游引物(10 μmol/L)各1.5 μL,dNTP(2.5 mmol/L)2 μL,Taq-DNA聚合酶1 U,DNA模板(100 ng/μL)1 μL,補(bǔ)無菌去離子水至25 μL。反應(yīng)程序:94 ℃變性5min;94 ℃變性30 s,各引物 55~62 ℃ 復(fù)性30 s,72 ℃延伸35 s,30個(gè)循環(huán);最后72 ℃延伸10 min。

1.3 數(shù)據(jù)分析

使用PopGene 1.32軟件計(jì)算每個(gè)群體的等位基因數(shù)(observed number of alleles,Na)、有效等位基因數(shù)(effective number of alleles,Ae)、觀測(cè)雜合度(observed heterozygosity,Ho)、期望雜合度(expected heterozygosity,He)、群體間遺傳距離(genetic distance,Ds)和遺傳相似性系數(shù)(genetic similarity index,I)。利用ARLEQUIN 3.1軟件,計(jì)算群體遺傳分化指數(shù)Fst,估計(jì)群體間的遺傳分化,分析群體內(nèi)和群體間的分子方差(AMOVA)。根據(jù)Botstein等計(jì)算多態(tài)信息含量(polymorphism information content,PIC)[10],PIC計(jì)算公式為:PIC=1-(∑ni=1P2i)-(∑n-1i=1∑nj=i+12P2iP2j)。式中,n為某一位點(diǎn)上等位基因數(shù),Pi、Pj分別為第i和第j個(gè)等位基因在群體中的頻率,j=i+1。endprint

2結(jié)果與分析

2.1PCR擴(kuò)增結(jié)果

采用25對(duì)微衛(wèi)星引物對(duì)3個(gè)世代共84尾樣本進(jìn)行分析比較。由表2可見,25個(gè)基因座在3個(gè)群體中共檢測(cè)到了105個(gè)等位基因,大多集中在100~300 bp之間,不同基因座位上的等位基因數(shù)在2~7個(gè)之間,平均每個(gè)基因座上的等位基因數(shù)為4.2個(gè)。采用12%非變性聚丙烯酰胺凝膠電泳,得到清晰且重復(fù)性較好的產(chǎn)物。圖1所示為引物GM277分別在3個(gè)群體中的擴(kuò)增圖譜。

2.23個(gè)世代間的遺傳多樣性

由表3可見,3個(gè)世代的平均有效等位基因數(shù)(Ae)分別為2.59、2.52、2.30,平均期望雜合度(He)分別為0.547、0545、0.518,平均觀測(cè)雜合度(Ho)分別為0.479、0457、0451,平均多態(tài)信息含量(PIC)分別為0.487、0.480、0.458,3個(gè)羅非魚保種選育群體的遺傳多樣性均屬于中度偏高水平。其中,F(xiàn)1代遺傳多樣性最高,F(xiàn)3代最低,呈現(xiàn)出略微下降的趨勢(shì)。經(jīng)t檢驗(yàn)P>0.05,3世代間無明顯差異。

2.33個(gè)世代間的遺傳差異

3個(gè)羅非魚群體間的Neis遺傳距離(Ds)和遺傳相似系數(shù)見表4。在F1、F2、F3 3個(gè)世代中,F(xiàn)1和F3的遺傳距離最大,為0.077,遺傳相似系數(shù)最小,為0.926,說明兩者親緣關(guān)系較遠(yuǎn);F2與F3的遺傳距離最小,為0.040,遺傳相似系數(shù)最大,為0.961,說明兩者親緣關(guān)系較近。

Wright認(rèn)為,遺傳分化指數(shù)Fst≤0.05表示群體間遺傳分化微弱;0.05

3小結(jié)與討論

3.1埃及尼羅羅非魚各世代保種群的遺傳分析

通過對(duì)埃及品系尼羅羅非魚3個(gè)世代保種選育群體進(jìn)行遺傳多樣性檢測(cè),對(duì)遺傳多樣性參數(shù)進(jìn)行比較,結(jié)果顯示,3個(gè)羅非魚選育群體平均多態(tài)信息含量(PIC)分別為0.487、0480、0.458,均處于中度偏高的遺傳多樣性水平[10],與Hassanien等研究埃及尼羅羅非魚5個(gè)野生群體遺傳多樣性得出的結(jié)果(Ae=8.4、Ho=0.702、He=0.869)[12]相比,這3個(gè)羅非魚選育群體的遺傳多樣性水平較低(Ae=2.47、Ho=0.462、He=0.537),這可能與引種后尼羅羅非魚為人工養(yǎng)殖群體,并且經(jīng)歷了選育有關(guān)。Macaranas等對(duì)尼羅羅非魚非洲野生群體和亞洲養(yǎng)殖群體進(jìn)行了遺傳差異分析發(fā)現(xiàn),養(yǎng)殖群體遺傳多樣性高于野生群體,推測(cè)是由于尼羅羅非魚養(yǎng)殖群體與莫桑比克羅非魚發(fā)生了漸滲雜交的原因[13]。本試驗(yàn)?zāi)崃_羅非魚遺傳多樣性雖然低于其他學(xué)者的研究結(jié)果,但其遺傳多樣性仍處于較高水平,3個(gè)保種世代群體的遺傳多樣性有逐代降低的趨勢(shì),世代間無明顯差異(P>0.05),這反映了埃及品系尼羅羅非魚保種群體的遺傳多樣性得到了很好的保存,保種選育效果較好。

埃及品系尼羅羅非魚從F1至F3代遺傳差異很?。‵st<0.05),通過分子方差對(duì)遺傳變異來源分析,經(jīng)過人工選育,群體世代間遺傳分化并不大,世代間變異僅占總變異的334%,分化主要來源于世代內(nèi)個(gè)體間的遺傳差異。F3與F2之間的分化小于F2與F1之間的分化,推測(cè)相鄰世代間的分化隨著世代增加可能呈現(xiàn)出逐漸降低的趨勢(shì),群體遺傳結(jié)構(gòu)逐漸趨于穩(wěn)定。從 3個(gè)選育世代間的遺傳距離和遺傳相似性結(jié)果可以看出,隨著選育的進(jìn)行,F(xiàn)2、F3與 F1之間的遺傳距離逐漸增大,而遺傳相似性系數(shù)逐漸減小,這表明埃及品系尼羅羅非魚經(jīng)過3代群體選育,人工選擇壓力在一定程度上改變了選育群體的遺傳結(jié)構(gòu),部分基因座上等位基因頻率發(fā)生改變,可以獲得具有更穩(wěn)定遺傳水平的選育良種。

3.2保種選育方法的可行性探討

保種就是有效地控制群體的近交和遺傳漂變,盡可能多地保存群體原有的基因和雜合度,群體規(guī)模是決定保種效果的基本因素[14]。為了避免種群中某些基因的丟失,生產(chǎn)單位多以大量的群體進(jìn)行保種和引種。2003年,廣西水產(chǎn)研究所從南京農(nóng)業(yè)大學(xué)無錫漁業(yè)學(xué)院引進(jìn)埃及品系尼羅羅非魚種苗51 360尾,經(jīng)多代選育已獲得較為理想的苗種繁育親本[15],其遺傳多樣性也得到了良好的保存[16]。從維持種群遺傳多樣性的角度來講,保留群體的規(guī)模越大越好,而在實(shí)際生產(chǎn)中,如果保留群體規(guī)模過大,就不可避免造成人力、物力資源的大量浪費(fèi)。本研究以雌魚600尾、雄魚300尾進(jìn)行保種,群體規(guī)模較小,親魚群體有效大小可達(dá)500尾以上,近交增量在0.1%以內(nèi),理論上可基本避免近親交配[17-18]。從微衛(wèi)星檢測(cè)結(jié)果來看,埃及尼羅羅非魚保種群體遺傳多樣性得到了很好的保存,本研究的保種選育方法具有一定的可行性,由于受選育保種數(shù)目所限,選育世代遺傳多樣性發(fā)生了一定程度的變化和下降。另外,本研究?jī)H通過微衛(wèi)星標(biāo)記對(duì)3個(gè)世代進(jìn)行了遺傳方面的分析和探討,還應(yīng)結(jié)合生長(zhǎng)指標(biāo)、血液生化指標(biāo)等其他生物學(xué)指標(biāo)對(duì)保種效果進(jìn)行進(jìn)一步分析,從而制定出更為合理有效的保種方案。

參考文獻(xiàn):

[1]夏德全 .論羅非魚產(chǎn)業(yè)發(fā)展[J]. 海洋與漁業(yè),2003(8):7-9.

[2]Naish K A,Skibinski D O F. Tetranucleiotide microsatellites loci for Indian major carp[J]. Journal of Fish Biol,1998,53(4):886-889.

[3]王婷婷,宋學(xué)宏,許愛國,等. 應(yīng)用微衛(wèi)星標(biāo)記分析4個(gè)黃顙魚群體的遺傳多樣性[J]. 江蘇農(nóng)業(yè)科學(xué),2012,40(4):41-45.

[4]Banks M A,Blouin M S,Baldwin B A,et al. Isolation and inheritance of novel microsatellites in chinook salmon(Oncorhynchus tschawytscha)[J]. The Journal of Heredity,1999,90(2):281-288.

[5]Kocher T D,Lee W J,Sobolewska H,et al. A genetic linkage map of a cichlid fish,the tilapia(Oreochromis niloticus)[J]. Genetics,1998,148(3):1225-1232.

[6]Romana-Eguia M R,Ikeda M,Basiao Z U,et al. Genetic changes during mass selection for growth in Nile tilapia,Oreochromis niloticus(L.),assessed by microsatellites[J]. Aquaculture Research,2005,36(1):69-78.

[7]Patricia D F,Pedrom M J. Assessment of the genetic diversity in five generations of a commercial broodstock line of Litopenaeus vannamei shrimp[J]. African Journal of Biotechnology,2005,4(12):1362-1367.

[8]頡曉勇,李思發(fā),蔡完其. 吉富品系尼羅羅非魚選育過程中遺傳變異的微衛(wèi)星分析[J]. 水產(chǎn)學(xué)報(bào),2007,31(3):385-390.

[9]Lee B Y,Lee W J,Streelman J T,et al. A second-generation genetic linkage map of tilapia(Oreochromis spp.)[J]. Genetics,2005,170(1):237-244.

[10]Botstein D,White R L,Skolnick M,et al. Construction of a genetic linkage map in man using restriction fragment length polymorphisms[J]. American Journal of Human Genetics,1980,32(3):314-331.

[11]Wright S. Evolution and the genetics of populations[M]. Chicago:University of Chicago Press,1978:1-30.

[12]Hassanien H A,Gilbey J. Genetic diversity and differentiation of Nile tilapia(Oreochromis niloticus)revealed by DNA microsatellites[J]. Aquaculture Research,2005,36(14):1450-1457.

[13]Macaranas J M,Agustin L Q,Ablan M A,et al. Genetic improvement of farmed tilapias:biochemical characterization of strain differences in Nile tilapia[J]. Aquaculture International,1995,3(1):43-54.

[14]Wright S. Evolution in mendelian populations[J]. Genetics,1931,16(2):97-157.

[15]曹謹(jǐn)玲,陳劍杰,楊弘,等. 埃及品系尼羅羅非魚的選育及其效果分析[J]. 水生生物學(xué)報(bào),2010,34(4):866-871.

[16]肖煒,楊弘,李大宇,等. 同種羅非魚(Tilapia)不同地區(qū)選育群體的遺傳多樣性分析[J]. 海洋與湖沼,2010,41(4):530-537.

[17]李思發(fā). 魚類繁育群體遺傳性能的保護(hù)[J]. 水產(chǎn)學(xué)報(bào),1988,12(3):283-290.

[18]Franklin I R. Evolutionary change in small populations[C]//Soule M E,Wilcox B A. Conservation biology:An evolutionary-ecological perspective. Sunderland:Sinauer Associate,1980:135-150.endprint

[5]Kocher T D,Lee W J,Sobolewska H,et al. A genetic linkage map of a cichlid fish,the tilapia(Oreochromis niloticus)[J]. Genetics,1998,148(3):1225-1232.

[6]Romana-Eguia M R,Ikeda M,Basiao Z U,et al. Genetic changes during mass selection for growth in Nile tilapia,Oreochromis niloticus(L.),assessed by microsatellites[J]. Aquaculture Research,2005,36(1):69-78.

[7]Patricia D F,Pedrom M J. Assessment of the genetic diversity in five generations of a commercial broodstock line of Litopenaeus vannamei shrimp[J]. African Journal of Biotechnology,2005,4(12):1362-1367.

[8]頡曉勇,李思發(fā),蔡完其. 吉富品系尼羅羅非魚選育過程中遺傳變異的微衛(wèi)星分析[J]. 水產(chǎn)學(xué)報(bào),2007,31(3):385-390.

[9]Lee B Y,Lee W J,Streelman J T,et al. A second-generation genetic linkage map of tilapia(Oreochromis spp.)[J]. Genetics,2005,170(1):237-244.

[10]Botstein D,White R L,Skolnick M,et al. Construction of a genetic linkage map in man using restriction fragment length polymorphisms[J]. American Journal of Human Genetics,1980,32(3):314-331.

[11]Wright S. Evolution and the genetics of populations[M]. Chicago:University of Chicago Press,1978:1-30.

[12]Hassanien H A,Gilbey J. Genetic diversity and differentiation of Nile tilapia(Oreochromis niloticus)revealed by DNA microsatellites[J]. Aquaculture Research,2005,36(14):1450-1457.

[13]Macaranas J M,Agustin L Q,Ablan M A,et al. Genetic improvement of farmed tilapias:biochemical characterization of strain differences in Nile tilapia[J]. Aquaculture International,1995,3(1):43-54.

[14]Wright S. Evolution in mendelian populations[J]. Genetics,1931,16(2):97-157.

[15]曹謹(jǐn)玲,陳劍杰,楊弘,等. 埃及品系尼羅羅非魚的選育及其效果分析[J]. 水生生物學(xué)報(bào),2010,34(4):866-871.

[16]肖煒,楊弘,李大宇,等. 同種羅非魚(Tilapia)不同地區(qū)選育群體的遺傳多樣性分析[J]. 海洋與湖沼,2010,41(4):530-537.

[17]李思發(fā). 魚類繁育群體遺傳性能的保護(hù)[J]. 水產(chǎn)學(xué)報(bào),1988,12(3):283-290.

[18]Franklin I R. Evolutionary change in small populations[C]//Soule M E,Wilcox B A. Conservation biology:An evolutionary-ecological perspective. Sunderland:Sinauer Associate,1980:135-150.endprint

[5]Kocher T D,Lee W J,Sobolewska H,et al. A genetic linkage map of a cichlid fish,the tilapia(Oreochromis niloticus)[J]. Genetics,1998,148(3):1225-1232.

[6]Romana-Eguia M R,Ikeda M,Basiao Z U,et al. Genetic changes during mass selection for growth in Nile tilapia,Oreochromis niloticus(L.),assessed by microsatellites[J]. Aquaculture Research,2005,36(1):69-78.

[7]Patricia D F,Pedrom M J. Assessment of the genetic diversity in five generations of a commercial broodstock line of Litopenaeus vannamei shrimp[J]. African Journal of Biotechnology,2005,4(12):1362-1367.

[8]頡曉勇,李思發(fā),蔡完其. 吉富品系尼羅羅非魚選育過程中遺傳變異的微衛(wèi)星分析[J]. 水產(chǎn)學(xué)報(bào),2007,31(3):385-390.

[9]Lee B Y,Lee W J,Streelman J T,et al. A second-generation genetic linkage map of tilapia(Oreochromis spp.)[J]. Genetics,2005,170(1):237-244.

[10]Botstein D,White R L,Skolnick M,et al. Construction of a genetic linkage map in man using restriction fragment length polymorphisms[J]. American Journal of Human Genetics,1980,32(3):314-331.

[11]Wright S. Evolution and the genetics of populations[M]. Chicago:University of Chicago Press,1978:1-30.

[12]Hassanien H A,Gilbey J. Genetic diversity and differentiation of Nile tilapia(Oreochromis niloticus)revealed by DNA microsatellites[J]. Aquaculture Research,2005,36(14):1450-1457.

[13]Macaranas J M,Agustin L Q,Ablan M A,et al. Genetic improvement of farmed tilapias:biochemical characterization of strain differences in Nile tilapia[J]. Aquaculture International,1995,3(1):43-54.

[14]Wright S. Evolution in mendelian populations[J]. Genetics,1931,16(2):97-157.

[15]曹謹(jǐn)玲,陳劍杰,楊弘,等. 埃及品系尼羅羅非魚的選育及其效果分析[J]. 水生生物學(xué)報(bào),2010,34(4):866-871.

[16]肖煒,楊弘,李大宇,等. 同種羅非魚(Tilapia)不同地區(qū)選育群體的遺傳多樣性分析[J]. 海洋與湖沼,2010,41(4):530-537.

[17]李思發(fā). 魚類繁育群體遺傳性能的保護(hù)[J]. 水產(chǎn)學(xué)報(bào),1988,12(3):283-290.

[18]Franklin I R. Evolutionary change in small populations[C]//Soule M E,Wilcox B A. Conservation biology:An evolutionary-ecological perspective. Sunderland:Sinauer Associate,1980:135-150.endprint

主站蜘蛛池模板: …亚洲 欧洲 另类 春色| 三上悠亚精品二区在线观看| 无码电影在线观看| 免费毛片全部不收费的| 日日噜噜夜夜狠狠视频| 91精品久久久无码中文字幕vr| 欧美a√在线| 激情视频综合网| 自偷自拍三级全三级视频 | 亚洲日韩欧美在线观看| 亚洲三级色| 亚洲一区免费看| 五月综合色婷婷| 日韩精品毛片| 国产欧美另类| 亚洲大学生视频在线播放 | 日韩免费无码人妻系列| 久久精品人妻中文系列| 国产精品福利导航| 澳门av无码| 视频一区亚洲| 精品一区国产精品| 91在线激情在线观看| 一本二本三本不卡无码| 国产综合色在线视频播放线视| 网久久综合| 免费国产一级 片内射老| 国产va免费精品观看| 日韩福利视频导航| 国产成人精品在线| 最新精品久久精品| 亚洲国产精品美女| 一级毛片在线免费看| 国产流白浆视频| 无码aaa视频| 亚洲国产综合自在线另类| 亚洲日本一本dvd高清| 青青草原国产一区二区| 亚洲欧美成人综合| 国产综合精品日本亚洲777| 亚洲天堂视频在线免费观看| 黄片一区二区三区| 欧美一区国产| 欧洲亚洲欧美国产日本高清| 欧美有码在线| 亚洲第一区在线| 蜜芽国产尤物av尤物在线看| 国产精品刺激对白在线| 亚洲丝袜第一页| 高清久久精品亚洲日韩Av| 欧美午夜在线视频| 动漫精品中文字幕无码| 国产男人天堂| 亚洲AV成人一区国产精品| 67194亚洲无码| 久久久久无码精品| 超碰色了色| 久久精品国产国语对白| 精品视频福利| 91成人试看福利体验区| 欧洲一区二区三区无码| 国产毛片片精品天天看视频| 一级毛片不卡片免费观看| 最新国产成人剧情在线播放| 成年A级毛片| 免费网站成人亚洲| 综合久久久久久久综合网| 国产凹凸一区在线观看视频| 日韩AV无码免费一二三区| 97成人在线观看| 亚洲va视频| 真实国产乱子伦高清| 97视频在线观看免费视频| 伊人久久福利中文字幕| 99久久无色码中文字幕| 国产亚洲精品97AA片在线播放| 亚洲第一黄色网址| 亚洲av无码成人专区| 亚洲av成人无码网站在线观看| 欧美在线视频a| 蜜臀av性久久久久蜜臀aⅴ麻豆 | 在线日韩日本国产亚洲|