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

蘆筍種質(zhì)資源的SRAP遺傳多樣性分析

2016-04-09 08:13:15李書(shū)華于繼慶李保華包艷存山東省濰坊市農(nóng)業(yè)科學(xué)院山東濰坊261071
農(nóng)學(xué)學(xué)報(bào) 2016年3期

李 霞,李書(shū)華,楊 林,于繼慶,李保華,鄭 鑫,包艷存(山東省濰坊市農(nóng)業(yè)科學(xué)院,山東濰坊261071)

?

蘆筍種質(zhì)資源的SRAP遺傳多樣性分析

李霞,李書(shū)華,楊林,于繼慶,李保華,鄭鑫,包艷存
(山東省濰坊市農(nóng)業(yè)科學(xué)院,山東濰坊261071)

摘要:為了更深入地揭示蘆筍品種(系)的遺傳多樣性,本研究采用SRAP技術(shù)對(duì)國(guó)內(nèi)外12個(gè)蘆筍栽培種和雜交種的遺傳多態(tài)性進(jìn)行分析。結(jié)果表明,從170對(duì)引物中篩選出29對(duì)多態(tài)性強(qiáng)、條帶清晰的引物,每對(duì)引物產(chǎn)生12~36條譜帶,共產(chǎn)生551條譜帶,其中多態(tài)性條帶343條,占總帶數(shù)的61.21%。平均每個(gè)引物擴(kuò)增的DNA帶數(shù)是20.41條,多態(tài)性帶數(shù)是12.7條。采用NTSYSpc2.1軟件對(duì)SRAP擴(kuò)增結(jié)果進(jìn)行聚類分析,12份蘆筍材料的相似系數(shù)在0.39~0.97之間,平均相似系數(shù)為0.69。在Coefficient=0.69處劃等值線,可將12份蘆筍材料劃分為4類。第一類:‘魯蘆筍一號(hào)’、‘88-5’、‘格蘭德’、‘07-1’和‘06-6’。第二類:‘阿波羅’、‘澤西巨人’和‘阿特拉斯’。第三類:‘西班牙達(dá)瑪斯’。第四類:‘荷蘭全雄’、‘UC800’和‘西德全雄’。

關(guān)鍵詞:蘆筍;種質(zhì)資源;SRAP;遺傳多樣性

0 引言

相關(guān)序列擴(kuò)增多態(tài)性(sequence-related amplified polymorphism,SRAP)是一種新型的基于PCR的標(biāo)記系統(tǒng),由美國(guó)加州大學(xué)蔬菜作物系Li與Quiros博士于2001年提出[1]。該技術(shù)與其他分子標(biāo)記技術(shù)相比有明顯的優(yōu)越性,如實(shí)驗(yàn)操作簡(jiǎn)單、多態(tài)性高、在基因組中分布均勻、實(shí)驗(yàn)結(jié)果穩(wěn)定、引物具有通用性等。近年來(lái)SRAP技術(shù)已被用于黃花蒿[2]、麻櫟[3]、甜瓜[4]、蔥[5]、桂花[6]等作物的遺傳多樣性分析和親緣關(guān)系鑒定等方面。蘆筍(Asparagus officinalis L.),又稱石刁柏,是百合科天門冬屬植物,該屬約有300種,除美洲外,全世界都有分布。核型分析表明該屬植物有二倍體、四倍體和六倍體,共同的染色體基數(shù)x=10。目前,國(guó)內(nèi)外對(duì)蘆筍的研究多集中在同工酶[7-8]、繁殖栽培技術(shù)[9-10]、細(xì)胞學(xué)等方面的研究[11-12],而關(guān)于分子水平的研究較少,李媛媛[13]對(duì)蘆筍性別連鎖分子標(biāo)記及蘆筍性別決定基因克隆的研究概況作了綜述。周勁松等[14]2007年,研究表明STS標(biāo)記Asp1-T7在雄性DNA池與雌性DNA池間具有多態(tài)性,可以用來(lái)對(duì)蘆筍雄性兩性株S1群體進(jìn)行性別輔助選擇。周勁松等[15]2012年研究表明顯性標(biāo)記Asp-T7和Asp1-T7sp特異性好,擴(kuò)增效果穩(wěn)定,均能夠鑒別蘆筍雌、雄性別;一共測(cè)驗(yàn)了14株蘆筍兩性株自交后代S1代的雄株,從中篩選出2株超雄株。高武軍[16]利用RAPD技術(shù)對(duì)石刁柏雌雄株基因組差異進(jìn)行了分析,發(fā)現(xiàn)了2條雌性相關(guān)的RAPD片斷,長(zhǎng)度分別為928、1178 bp。蘇俊祥[17]利用RACE技術(shù)從石刁柏中分離了一個(gè)SUPERMAN類單鋅指蛋白的基因AoSUP,并對(duì)AoSUP的基因序列、蛋白結(jié)構(gòu)、雌雄株不同營(yíng)養(yǎng)組織表達(dá)差異進(jìn)行了初步分析。李霞等[18]以蘆筍(Asparagus officinalis L.)為材料,對(duì)影響SRAP-PCR擴(kuò)增結(jié)果的因素dNTPs、Taq DNA酶、Mg2+、引物、模板DNA進(jìn)行了正交優(yōu)化設(shè)計(jì),確立了適合蘆筍SRAP-PCR分析的最佳反應(yīng)體系。尋找蘆筍種質(zhì)間遺傳多樣性及親緣關(guān)系分析中適用的有效分子標(biāo)記技術(shù),是目前蘆筍品種鑒定、遺傳關(guān)系分析及種質(zhì)資源評(píng)價(jià)工作中急需解決的問(wèn)題。本研究就采用SRAP技術(shù)對(duì)12份蘆筍種質(zhì)資源進(jìn)行多態(tài)性分析,在DNA水平上探索蘆筍種質(zhì)資源的遺傳距離,以便為蘆筍種質(zhì)資源的有效利用在分子水平上提供參考。

1 材料與方法

1.1供試材料

所有的試驗(yàn)材料均來(lái)自于山東省濰坊市農(nóng)業(yè)科學(xué)院試驗(yàn)基地,選取當(dāng)年萌發(fā)的幼嫩筍或葉片,洗凈、吸干水分,于-80℃冰箱中保存?zhèn)溆谩9?2份(表1)。

表1 供試12份蘆筍材料

1.2試驗(yàn)方法

1.2.1基因組DNA提取采用改良的CTAB法[19]提取蘆筍DNA,產(chǎn)物用1.0%的瓊脂糖凝膠電泳和紫外分光光度計(jì)檢測(cè),選出符合試驗(yàn)要求的DNA樣品,-20℃保存。

1.2.2 SRAP-PCR擴(kuò)增擴(kuò)增程序:94℃預(yù)變性5 min;94℃變性1 min,35℃復(fù)性1 min,72℃延伸1.5 min,共5個(gè)循環(huán);30 s后,94℃變性1 min,50℃復(fù)性1 min,72℃延伸1.5 min,共35個(gè)循環(huán);72℃延伸7 min,4℃保存。1.2.3電泳檢測(cè)擴(kuò)增反應(yīng)在美國(guó)ApolloATC-401擴(kuò)增儀上完成。擴(kuò)增產(chǎn)物采用6%變性聚丙烯酰胺凝膠電泳,電泳緩沖液為1×TBE,恒功率70 W,溴酚藍(lán)距離凝膠2~3 cm處結(jié)束電泳。電泳后銀染。蒸餾水沖洗后照相或者保存。

2 結(jié)果與分析

2.1基因組DNA的質(zhì)量檢測(cè)

將提取的12份基因組DNA用TE溶解后,1%瓊脂糖凝膠電泳和溴化乙錠染色檢測(cè),凝膠成像系統(tǒng)觀察,照相。結(jié)果顯示,蘆筍基因組DNA帶型表現(xiàn)整齊、清晰(圖1)。

圖1 12份蘆筍基因組DNA電泳結(jié)果

2.2 SRAP-PCR反應(yīng)體系的正交優(yōu)化

蘆筍SRAP-PCR反應(yīng)體系參照李霞[18]的正交優(yōu)化反應(yīng)體系。20 μL的反應(yīng)體積中包括0.2 mmol/LdNTPs 0.4 μL;1 U Taq DNA酶1.0 μL;1.5 mmol/L Mg2+1.2 μL;上下引物各30 ng,每個(gè)引物0.5 μL;60 ng模板DNA1.0μL;10×Buffer2.0μL;無(wú)菌雙蒸水13.4μL。

2.3引物組合篩選

所用引物均購(gòu)于上海生工生物工程技術(shù)服務(wù)有限公司,上引物17個(gè),下引物10個(gè),隨機(jī)排列形成170對(duì)引物(表2)。

本試驗(yàn)用優(yōu)化后的體系對(duì)所有170對(duì)引物組合進(jìn)行擴(kuò)增,所用材料為12份蘆筍材料。結(jié)果顯示,有29對(duì)組合均能擴(kuò)增出多態(tài)性強(qiáng)、清晰度高的條帶帶(圖2)。

2.4蘆筍SRAP標(biāo)記遺傳多樣性分析

本試驗(yàn)用篩選出29對(duì)引物對(duì)12份蘆筍材料進(jìn)行重新擴(kuò)增,結(jié)果顯示,擴(kuò)增片段大小在100~1500 bp之間,每對(duì)引物組合產(chǎn)生12~36條譜帶,共產(chǎn)生551條譜帶,其中多態(tài)性條帶343條,占總帶數(shù)的61.21%。平均每對(duì)引物擴(kuò)增的DNA帶數(shù)是20.41條,多態(tài)性帶數(shù)是12.7條(表3)。不同的引物組合擴(kuò)增出的總帶數(shù)和多態(tài)性帶數(shù)的差異較大,產(chǎn)生帶數(shù)對(duì)多的引物組合是K1,擴(kuò)增出36條帶;產(chǎn)生帶數(shù)最少的引物組合有2個(gè),分別是E7、K10,擴(kuò)增出12條帶。產(chǎn)生多態(tài)性最高的引物組合是K1,多態(tài)性百分率為86.11%;多態(tài)性較低的引物組合是C3,多態(tài)性百分率為42.86%。引物組合F6擴(kuò)增的SRAP帶型(圖3)。

表2 SRAP隨機(jī)引物的序列

圖2 引物組合篩選擴(kuò)增圖譜

2.5聚類分析

對(duì)12份蘆筍材料的SRAP擴(kuò)增結(jié)果進(jìn)行數(shù)據(jù)統(tǒng)計(jì),每個(gè)材料的擴(kuò)增帶按有或無(wú)記錄,“有”附值為1,“無(wú)”賦值為0,從而得到原始數(shù)據(jù)矩陣。采用NTSYSpc 2.1軟件計(jì)算基因型間的Nei氏遺傳相似系數(shù)(表4)。UPGMA法聚類構(gòu)建樹(shù)狀圖(圖4)。由表4可以看出,12份蘆筍材料的相似系數(shù)在0.39~0.97之間,平均相似系數(shù)為0.69,遺傳多樣性相對(duì)偏低。

由圖4可知,在Coefficient=0.69處劃等值線,可將12份蘆筍材料劃分為4類。第一類:‘魯蘆筍一號(hào)’、‘88-5’、‘格蘭德’、‘07-1’和‘06-6’。第二類:‘阿波羅’、‘澤西巨人’、‘阿特拉斯’。第三類:‘西班牙大瑪斯’。第四類:‘荷蘭全雄’、‘UC800’和‘西德全雄’。

表3 29個(gè)引物組合產(chǎn)生的多態(tài)性條帶數(shù)

圖3 引物組合F6的擴(kuò)增產(chǎn)物

從聚類圖上可以看出,各品種的地理來(lái)源同聚類結(jié)果有一定關(guān)系。如來(lái)源于中國(guó)的‘魯蘆筍一號(hào)’和‘88-5’改良系(濰坊市農(nóng)科院自培品種)、來(lái)自美國(guó)的‘格蘭德’和‘07-1’,相似系數(shù)均高達(dá)0.97;在第一類中‘07-1’和‘06-6’均來(lái)自于美國(guó),相似系數(shù)為0.94。美國(guó)培育的‘澤西巨人’和‘阿特拉斯’(相似系數(shù)為0.81)聚為一小類后,又和‘阿波羅’聚為一類(相似系數(shù)為0.77)。在相似系數(shù)為0.69水平上‘西班牙大瑪斯’單獨(dú)聚為一類。從圖中還可以看出不同的地理來(lái)源,但同為全雄品種也聚為一類,如‘荷蘭全雄’、‘UC800’和‘西德全雄’聚為一類,相似系數(shù)為0.81~0.90,平均相似系數(shù)為0.85。

3 結(jié)論與討論

SRAP分子標(biāo)記技術(shù)不僅多態(tài)性高、重復(fù)性好、操作簡(jiǎn)單、在基因組中分布均勻、引物具有通用性,而且其正向引物可以與反向引物兩兩搭配組合,因此用少量的引物可組配得到多個(gè)引物組合,提高了引物的使用效率,降低合成引物成本。因而利用SRAP構(gòu)建品種指紋圖譜是一種經(jīng)濟(jì)實(shí)用的方法。本研究采用SRAP技術(shù)對(duì)國(guó)內(nèi)外12個(gè)蘆筍栽培種和雜交種的遺傳多態(tài)性進(jìn)行分析。從170對(duì)引物中篩選出29對(duì)多態(tài)性強(qiáng)、條帶清晰的引物,每對(duì)引物產(chǎn)生12~36條譜帶,共產(chǎn)生551條譜帶,其中多態(tài)性條帶343條,占總帶數(shù)的61.21%。平均每個(gè)引物擴(kuò)增的DNA帶數(shù)是20.41條,多態(tài)性帶數(shù)是12.7條。采用NTSYSpc2.1軟件對(duì)SRAP擴(kuò)增結(jié)果進(jìn)行聚類分析,12份蘆筍材料的相似系數(shù)在0.39~0.97之間,平均相似系數(shù)為0.69。在Coefficient=0.69處劃等值線,可將12份蘆筍材料劃分為4類。第一類:‘魯蘆筍一號(hào)’、‘88-5’、‘格蘭德’、‘07-1’和‘06-6’。第二類:‘阿波羅’、‘澤西巨人’和‘阿特拉斯’。第三類:‘西班牙達(dá)瑪斯’。第四類:‘荷蘭全雄’、‘UC800’和‘西德全雄’。試驗(yàn)結(jié)果與李可峰[20]的比較吻合,說(shuō)明蘆筍總體上差異不大,特別是遺傳基礎(chǔ)在逐漸變窄,與本研究結(jié)果相同。

為此,中國(guó)蘆筍育種除了作好現(xiàn)有優(yōu)良品種的保藏和選育外,還應(yīng)做好以下4方面工作:(1)繼續(xù)選育自有優(yōu)良品種。目前專門從事蘆筍育種研究的有江西省農(nóng)業(yè)科學(xué)院、山東省濰坊市農(nóng)業(yè)科學(xué)院、北京市農(nóng)林科學(xué)院等,育成的品種有15個(gè),在此基礎(chǔ)上繼續(xù)進(jìn)行蘆筍自有優(yōu)良品種的選育完全有可能實(shí)現(xiàn)。(2)加大兩性株的應(yīng)用。蘆筍雌雄異株,在自然情況下,蘆筍雄株和雌株的株數(shù)大體各占一半。蘆筍性別是由一對(duì)獨(dú)立基因控制的,蘆筍雄性兩性株的遺傳規(guī)律不同于雌、雄株。它們同株可以自交,其自交后代75%為雄株(1/3為同質(zhì)雄株,2/3為異質(zhì)雄株),25%為雌株。在整個(gè)自交選擇中,可以篩選出超雄株,繼而快速選育出全雄品種,蘆筍產(chǎn)量大幅度提升。(3)優(yōu)良品種與野生近緣種雜交。選育出具有野生近緣種抗病性強(qiáng)、耐貧瘠、抗風(fēng)沙等的品種。(4)加大對(duì)外合作力度,廣泛搜集國(guó)外具有優(yōu)良性狀的品種,拓寬親本選擇范圍,徹底改變中國(guó)蘆筍育成品種遺傳基礎(chǔ)狹窄的局面。

表4 12份蘆筍品種的Nei氏相似系數(shù)矩陣

圖4 基于SRAP標(biāo)記的12份蘆筍材料聚類圖

參考文獻(xiàn)

[1]Li G,Quiros C F.Sequence-related amplified polymorphism (SRAP),a new marker system based on a simple PCR reaction:its application to mapping and gene tagging in Brassica[J].Theor Appl Genet,2001,103:455-461.

[2]鄧蜻.黃花蒿SRAP-PCR反應(yīng)體系的建立與優(yōu)化[J].中草藥,2007, 38(1):125-128.

[3]葉青雷,曾憲云.麻櫟SRAP-PCR體系優(yōu)化與遺傳多樣性分析[J].生物技術(shù),2009,19(3):24-27.

[4]Ferrio L M,Pico B,Cordava P F.Molecular diversity of a germplasm collection of sqash(Cucurbita moschata)determined by SRAP andAFLP makers[J].Crop Science,2004,44:653-664.

[5]李慧芝,尹燕枰.SRAP在蔥栽培品種遺傳多樣性研究中的適用性分析[J].園藝學(xué)報(bào),2007,34(4):929-934.

[6]李梅,侯喜林,郝日明.基于SRAP分子標(biāo)記的桂花品種親緣關(guān)系研究[J].園藝學(xué)報(bào),2009,36(11):1667-1675.

[7]李霞,牟萌,李書(shū)華,等.雌雄蘆筍過(guò)氧化物酶及酯酶同工酶的比較研究[J].江西農(nóng)業(yè)學(xué)報(bào),2009,11:34-36.

[8]洪達(dá),孫富從,姬艷克,等.蘆筍雌雄性別間同工酶酶譜的比較分析[J].河南師范大學(xué)學(xué)報(bào):自然科學(xué)版,2008,36(1):117-120.

[9]孫福來(lái).蘆筍無(wú)公害高產(chǎn)種植技術(shù)[J].中國(guó)種業(yè),2009,7:60-61.

[10]邱桂蘭,喬維鳳.蘆筍高產(chǎn)栽培技術(shù)[J].云南農(nóng)業(yè)科技,2009,2:36-37.

[11]高武軍,張小莉,尹為治,等.蘆筍性染色體的細(xì)胞學(xué)研究[J].河南農(nóng)業(yè)科學(xué),2008,10:118-122,140.

[12]張?jiān)獓?guó),倫君,王振華,等.蘆筍雄性系兩性株調(diào)查及鏡檢技術(shù)初步研究[J].山東農(nóng)業(yè)科學(xué),2006,4:10-12.

[13]李媛媛.蘆筍性別分子標(biāo)記研究進(jìn)展[J].中國(guó)蔬菜,2010,18:1-4.

[14]周勁松,湯泳萍,羅紹春,等.DNA分子標(biāo)記早期快速鑒別蘆筍雌雄株[J].分子植物育種,2007,5(3):363-366.

[15]周勁松,湯永萍,羅紹春,等.蘆筍超雄株的DNA分子標(biāo)記輔助篩選[J].園藝學(xué)報(bào),2012,39(11):2182-2188.

[16]高武軍.石刁柏雌雄株間細(xì)胞學(xué)、內(nèi)源物質(zhì)及分子標(biāo)記研究[D].開(kāi)封:河南大學(xué),2008.

[17]蘇俊祥.石刁柏SUPERMAN基因的克隆及初步分析[D].新鄉(xiāng):河南師范大學(xué),2012.

[18]李霞,于繼慶,李書(shū)華,等.利用正交設(shè)計(jì)優(yōu)化蘆筍SRAP反應(yīng)體系[J].農(nóng)學(xué)學(xué)報(bào),2011,3:5-9.

[19]Stewardt C N,Via L E.A rapid CTAB DNA isolation technique useful for RAPD fingerprinting and other PCR applications[J]. Biotechniques,1993,14:748-750.

[20]李可峰.石刁柏(Asparagus officinalis L.)種質(zhì)資源遺傳多樣性分析[D].濟(jì)南:山東師范大學(xué),2004.

Genetic Diversity Analysis of Asparagus Based on SRAP Markers

Li Xia,Li Shuhua,Yang Lin,Yu Jiqing,Li Baohua,Zheng Xin,Bao Yancun
(Weifang Academy of Agricultural Sciences,Weifang 261071,Shandong,China)

Abstract:This paper analyzed the genetic polymorphism of 12 asparagus cultivars and hybrids from home and abroad using SRAP method,so as to deeply reveal the genetic diversity of asparagus varieties.The results showed that 29 primers with high polymorphism and clear bands were screened from 170 primers,each of them amplified 12 to 36 bands and a total of 551 DNA fragments were obtained,343 bands showed polymorphism in all fragments,accounted for 61.21%.The average number of DNA bands produced by each primer was 20.41 and polymorphic bands were 12.7.The result of genetic similarity analysis of 12 asparagus genotypes showed that the Nei’s coefficient ranged from 0.39 to 0.97,and the average Nei’s coefficient was 0.69.A DNA molecular dendrogram was established for 12 genotypes which were divided into 4 groups according to the coefficient of 0.69.The first group included‘Lu Asparagus No.1’,‘88-5’,‘Gland’,‘07-1’and‘06-6’. The second group included‘Apollo’,‘Jersey Giant’and‘Atlas’.The third group included‘Spanish Damas’. The forth group included‘Holland all male’and‘West Germany all male’.

Key words:Asparagus officinalis L.;Germplasm;SRAP;Genetic Diversity

收稿日期:2015-10-15,修回日期:2015-12-16。

通訊作者:楊林,男,1980年出生,山東煙臺(tái)人,農(nóng)藝師,推廣碩士,主要從事蘆筍栽培研究。

作者簡(jiǎn)介:第一李霞,女,1976年出生,山東萊蕪人,副研究員,碩士,主要從事蘆筍遺傳育種的研究。

通信地址:261071山東省濰坊市勝利東街1921號(hào)山東省濰坊市農(nóng)業(yè)科學(xué)院,Tel:0536-2118589,E-mail:55654687@qq.com。 261071山東省濰坊市勝利東街1921號(hào)山東省濰坊市農(nóng)業(yè)科學(xué)院,Tel:0536-2118589,E-mail:80408239@qq.com。

基金項(xiàng)目:山東省科技發(fā)展計(jì)劃專項(xiàng)計(jì)劃項(xiàng)目“蘆筍超雄株技術(shù)研究與示范”(2012GZX21024)。

中圖分類號(hào):S644.6

文獻(xiàn)標(biāo)志碼:A論文編號(hào):cjas15100011

主站蜘蛛池模板: 亚洲国产清纯| 为你提供最新久久精品久久综合| 毛片a级毛片免费观看免下载| 日本欧美一二三区色视频| 亚洲Va中文字幕久久一区| 高潮爽到爆的喷水女主播视频| 香蕉99国内自产自拍视频| 国产亚洲视频免费播放| 男人的天堂久久精品激情| 亚洲国产精品不卡在线| 欧美色亚洲| 国产在线98福利播放视频免费| 国产欧美视频在线观看| 成人免费黄色小视频| 国产欧美视频在线| 欧美成人午夜在线全部免费| 亚洲国产理论片在线播放| av在线5g无码天天| 日本在线亚洲| 特级精品毛片免费观看| 国产福利不卡视频| 欧美中文一区| 久久一本精品久久久ー99| 欧美精品高清| 久草网视频在线| 国产亚洲高清视频| 国产一区自拍视频| 老司机午夜精品网站在线观看 | 精品久久综合1区2区3区激情| 幺女国产一级毛片| 亚洲成人动漫在线观看| 99久久精品视香蕉蕉| 韩日无码在线不卡| 孕妇高潮太爽了在线观看免费| 99久久精品免费视频| 久久亚洲国产最新网站| 好久久免费视频高清| 国产自在线播放| 伊人久久大香线蕉影院| 久久综合色天堂av| 成人午夜网址| 欧美全免费aaaaaa特黄在线| 精品国产成人三级在线观看| 99九九成人免费视频精品| 美女无遮挡拍拍拍免费视频| 国产伦片中文免费观看| 欧美日韩久久综合| 国产91在线|日本| 久久这里只精品国产99热8| 久久人与动人物A级毛片| 欧美国产中文| 性欧美在线| 免费人成黄页在线观看国产| 国产97视频在线| 九九热免费在线视频| 成人永久免费A∨一级在线播放| 国产精品国产主播在线观看| 国产国产人成免费视频77777| 波多野吉衣一区二区三区av| 国产成人一区免费观看| 久久黄色免费电影| 亚洲系列中文字幕一区二区| 精品久久久无码专区中文字幕| 999国产精品永久免费视频精品久久 | 无码丝袜人妻| 亚洲午夜综合网| 一级毛片基地| 欧美日韩一区二区在线免费观看 | 亚洲第一区欧美国产综合| 婷婷色中文网| 日韩欧美国产区| 亚洲成a∧人片在线观看无码| 无码精品国产dvd在线观看9久| 国产黄色免费看| 97在线公开视频| 亚洲乱码在线视频| 中文字幕在线看| 色网站在线免费观看| 国产一二三区视频| 中文字幕在线日韩91| 日韩一区二区在线电影| 亚洲成人在线免费观看|