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

晉大52\\57及其后代群體POD\\SOD同工酶遺傳多樣性及聚類分析

2011-01-01 00:00:00胡培培,李貴全
天津農(nóng)業(yè)科學(xué) 2011年1期

摘要:對(duì)晉大52(母本)、晉大57(父本)及其6個(gè)有代表性的后代SN-1、SN-2、SN-3、SN-4、SN-5、SN-6花莢期超氧化物歧化酶(S0D)和過(guò)氧化物酶(POD)兩種同工酶,進(jìn)行電泳圖譜檢測(cè)和分析,把SOD酶譜分為3個(gè)區(qū),A區(qū)有3條酶帶(Rf=0.27,0.41,0.44),B區(qū)有3條酶帶(Rf=0.55,0.59,0.62),C區(qū)有3條酶帶(Rf=0.68,0.74,0.79);POD同工酶電泳顯示有13條酶帶,其中3、5、11、13號(hào)酶帶出現(xiàn)頻率為100%,1、2、10號(hào)酶帶頻率均在85%以上,8種材料共出現(xiàn)78條酶帶,平均每種材料9.8條。后代品種與父本的平均遺傳距離為0.168 6,與母本的平均遺傳距離為0.387 6。

關(guān)鍵詞:大豆;同工酶;遺傳多樣性;聚類分析

中圖分類號(hào):S565.1文獻(xiàn)標(biāo)識(shí)碼:ADOI編碼:10.3969/j.issn.1006-6500.2011.01.002

Genetic Diversity and Cluster Analysis of Soybean POD, SOD Isozymes

HU Pei-pei,LI Gui-quan

(College of Agricultural,Shanxi Agricultural University,Taigu, Shanxi 030801,China)

Abstract:A population constructed from a cross between Jinda52 and Jinda57 were used as materials in this paper. By the analysis of POD electrophoresis spectrum, we concluded that there were 13 enzyme spectrum, and 78 enzyme spectrum in all 8 materials. SOD enzyme spectrum included 3 regions, the numbers of electrophoresis profile and the moving distance were basically the same in cross population, except for the color deepness in the profile, and so was SDS-PAGE electrophoresis, this is completely different from POD enzyme. The genetic distance between offspring and male is 0.168 6, between offspring and female is 0.387 6. These results were very significance to specific classifying and relationship researching.

Key words: soybean;isozymes;genetic diversity;cluster analysis

同工酶是由同一基因位點(diǎn)不同等位基因編碼的具有同一底物的蛋白質(zhì)分子, 在植物器官中普遍存在[1],具有組織、發(fā)育和物種的特異性。同時(shí)它也是較好的遺傳信息標(biāo)志,其電泳圖譜譜帶的多少及譜帶的多態(tài)性變化,反映了生物在蛋白質(zhì)水平上的多樣性,是生物發(fā)生遺傳分化的證據(jù),可作為遺傳多樣性的生化標(biāo)記。從20世紀(jì)60年代開始,同工酶被廣泛地應(yīng)用于大豆分類、遺傳[2-3]、變異[4]和親緣關(guān)系[5-7]研究中,獲得了一定的進(jìn)展。楊琪等[8]認(rèn)為大豆同工酶與雜種優(yōu)勢(shì)及產(chǎn)量性狀變異相關(guān),用酯酶與過(guò)氧化物酶同工酶分析后得出結(jié)論:具有互補(bǔ)帶的組合雜種優(yōu)勢(shì)高;類型相同的雜交組合F1代無(wú)互補(bǔ)酶帶;含秣食豆種質(zhì)組合F1代互補(bǔ)帶比率高;F2單株同工酶酶譜類型豐富的組合,產(chǎn)量性狀變異潛力大。缺失型酶譜組合的F1代,主要產(chǎn)量性狀均值最低。虞京葳等[9]利用PAGE法研究了大豆的分類,用酯酶同工酶分析大豆幼苗,結(jié)果酶帶明確,具有較為明顯的種的專一性,對(duì)栽培大豆(G.max)與野生大豆(G.soja)的酯酶分析同形態(tài)分類學(xué)研究相互驗(yàn)證結(jié)果較為一致。裴顏龍等[10]、鐘珍萍等[11]分別報(bào)道了對(duì)大豆亞屬間、種間的酶譜譜型分析,證實(shí)種群間有明顯的遺傳分化,指出染色體組型同酶帶譜型有關(guān)系,相近的材料譜型及酶活性相似,證明中國(guó)為野生大豆遺傳變異中心。

本研究對(duì)晉大52(母本)、晉大57(父本)及其6個(gè)有代表性的后代SN-1、SN-2、SN-3、SN-4、SN-5、SN-6花莢期過(guò)氧化物酶(POD) 、超氧化物歧化酶(SOD)兩種同工酶[12]進(jìn)行了電泳圖譜檢測(cè)和分析,通過(guò)對(duì)酶譜進(jìn)行相似性和聚類分析,研究了晉大52 、晉大57雜交親本及其后代群體之間酶系統(tǒng)的多樣性,旨在對(duì)今后進(jìn)行大豆品種的分類及親緣關(guān)系的研究提供借鑒。

1材料和方法

1.1試驗(yàn)材料

晉大52(母本)、晉大57(父本)及其6個(gè)有代表性的后代SN-1、SN-2、SN-3、SN-4、SN-5、SN-6。

1.2過(guò)氧化物酶(POD)

花莢期取幼葉各0.5 g,樣品提取液(Tris 1.21%,用HCl調(diào)pH=8.0)1.5 mL;采用不連續(xù)垂直板聚丙烯酰胺凝膠電泳法[13-15]。染色采用醋酸聯(lián)苯胺法。測(cè)量酶帶遷移率(Rf=酶帶遷移距離/前沿指示劑距離)。

1.3超氧化物歧化酶(SOD)

花莢期取幼葉各0.5 g,提取液(Tris-Gly pH8.7)1.5 mL;采用不連續(xù)垂直板聚丙烯酰胺凝膠電泳法,分離膠濃度為10% pH8.9 ,濃縮膠濃度為2.5% pH6.7。染色:(1)在2.45 mmol·L-1氮藍(lán)四唑(NBT中),遮光,浸搖20 min;(2)在3.6 mmol·L-1 pH7.8磷酸緩沖液(含28 mmol·L-1 四甲基己二胺,28 μmol·L-1核黃素,10 μmol·L-1 KCN)中浸搖15 min;(3)浸在含0.1 mmol·L-1EDTA的50 mmol·L-1 pH7.8磷酸緩沖液中,在光照培養(yǎng)箱中照射20~25 min,清水漂洗,照相 [16-18]。

1.4 大豆種子SDS-PAGE

每份試樣取單粒種子脫殼,置于無(wú)菌的雙層濾紙間,用小鐵錘碾碎成粉末,將粉碎樣品置于1.5 mL的離心管中,用三氯甲烷∶甲醇∶丙酮(2∶1∶1)進(jìn)行脫脂;提取樣品時(shí)加入Samplebuffer(pH6.8,1 mol·L-1Tris-HCl 5.0 mL+20%SDS 2 mL+甘油8.0 mL+α巰基乙醇3.6 mL+溴酚藍(lán)20.0 mg,稀釋4倍);電泳采用聚丙烯酰胺垂直板非連續(xù)系統(tǒng);用0.25%考馬斯亮藍(lán)R-250染色。

2結(jié)果與分析

2.1晉大52、晉大57及其后代種子SDS-PAGE分析

采用10%的SDS-PAGE法分離大豆HMW-GS,由圖1可以看出,品系間譜帶數(shù)目的變異范圍在20條帶左右,品系間譜帶基本一致,差異較小,只是在帶的強(qiáng)弱上有區(qū)別,只有CK1 (父本晉大57)出現(xiàn)了特異帶Rf=0.27。

2.2超氧化物歧化酶同工酶分析

本試驗(yàn)分離的SOD同工酶掃描圖見圖2,從圖2中可以看出,被測(cè)的8個(gè)材料中,花莢期SOD同工酶共出現(xiàn)9條酶帶,根據(jù)遷移率(Rf)可將酶譜分為3個(gè)區(qū):A區(qū)有3條酶帶(Rf=0.27,0.41,0.44),B區(qū)有3條酶帶(Rf=0.55,0.59,0.62),C區(qū)有3條酶帶(Rf=0.68,0.74,0.79),這與前人的研究結(jié)果是一致的[19-22]。所不同的是C區(qū)的3條酶帶顏色深、酶帶寬、酶活性強(qiáng);B區(qū)的酶帶顏色較深,活性較強(qiáng),酶帶比C區(qū)的相對(duì)較窄;A區(qū)第二條帶(Rf=0.41)顏色深,酶活性強(qiáng),第一、三條帶顏色淡、活性較弱。晉大52、晉大57及其后代群體超氧化物歧化酶同工酶譜帶基本相同(沒(méi)有條帶的缺失及遷移率的變異),有很大的同源性,品系間的差異僅體現(xiàn)在某些酶帶的強(qiáng)弱上。

2.3過(guò)氧化物酶同工酶分析

2.3.1過(guò)氧化物同工酶酶譜特征 分析POD同工酶電泳掃描圖,計(jì)算各酶帶的Rf值,繪制POD同工酶模式圖(圖3),共出現(xiàn)13條酶帶。本研究根據(jù)遷移率將酶譜分成3個(gè)區(qū):Ⅰ區(qū)(慢區(qū))Rf值在(0.25~0.35)之間,共5條帶;Ⅱ區(qū)(中區(qū))Rf值在(0.50~0.54)之間,共3條帶;Ⅲ區(qū)(快區(qū))Rf值在(0.59~0.69)之間,共5條帶。它們?cè)诓煌牧现谐霈F(xiàn)的頻率也各不相同(見表1),其中3、5、11、13號(hào)酶帶出現(xiàn)頻率為100%,所有試材中均有這兩條酶帶,這為晉大52、晉大57及其雜交后代過(guò)氧化物酶同工酶的特征譜帶。第1、2、10號(hào)酶帶出現(xiàn)頻率均在85%以上,可認(rèn)為是該群體的基本穩(wěn)定譜帶。從基本酶譜可知,Ⅰ區(qū)、Ⅲ區(qū)的酶帶數(shù)多,顯色深,酶帶寬;而Ⅱ區(qū)酶帶數(shù)少,顯色淺,酶帶窄,這說(shuō)明慢、快區(qū)的酶活性強(qiáng),較穩(wěn)定,在研究中可能有較高的利用價(jià)值。

從圖3、表1中可以看出:各供試材料在花莢期葉片POD同工酶酶譜特征明顯不同,8種材料出現(xiàn)的POD同工酶酶帶總數(shù)為78條,每種材料平均為9.8條,品系間酶帶數(shù)的差異也較大,最少的有8條,最多的達(dá)到12條,而且在酶帶的組合形式上也表現(xiàn)各異。說(shuō)明花莢期POD同工酶酶帶比較豐富,并且不同材料間POD同工酶酶譜酶帶數(shù)、酶帶位點(diǎn)(Rf值)、和酶帶強(qiáng)弱(顏色深淺)均存在著不同程度的差異,同時(shí)也說(shuō)明不同材料的POD同工酶蛋白質(zhì)表達(dá)水平不同,體現(xiàn)了基因水平上的差異。

2.3.2POD酶譜的相似性根據(jù)POD同工酶電泳模式圖,將有帶量化為1,無(wú)帶量化為0(表1),計(jì)算成對(duì)品系的相似系數(shù)Sij和遺傳距離Dij,計(jì)算公式為:

Sij=2Nij/(Ni+Nj)Dij=1 - Sij

其中Nij為兩個(gè)品系共有的條帶數(shù),Ni和Nj分別為第i個(gè)和第j個(gè)品系各自的條帶數(shù)[23]。 Sij越接近1.0,所比較的兩個(gè)品系的相似程度越高。由此可以得到相似系數(shù)矩陣(表2)。由表計(jì)算,SN-1、SN-2、SN-3、SN-4、SN-5、SN-6與父本的平均遺傳距離為0.168 6,與母本的平均遺傳距離為0.387 6。

2.3.3POD聚類分析在POD酶譜相似系數(shù)的基礎(chǔ)上采用SAS軟件系統(tǒng)聚類平均法(UPGMA)[24-25]進(jìn)行聚類分析,并得到樹形圖(圖4)。

由樹形圖(圖4)可以看出,在類平均距離為3.82的水平上,所有品系被聚為兩大類:第一大類包括晉大57、SN-4、 SN-1、SN-2、SN-5、SN-3,說(shuō)明SN-4、 SN-1、SN-2、SN-5、SN-3與父本晉大57親緣關(guān)系較近。該類在類平均距離為2.50處分又可以為兩個(gè)亞類,第一亞類(晉大57、SN-4、SN-1)與第二亞類(SN-2、SN-5、SN-3)。父本晉大57被聚在第一亞類中,可見SN-4、SN-1與父本的親緣關(guān)系更為接近;第二大類包括晉大52和SN-6,說(shuō)明SN-6與母本的親緣關(guān)系最為接近。聚類結(jié)果在一定程度上反映了子代品系間及子代與父母本親緣關(guān)系的遠(yuǎn)近。

3討論

3.1同工酶酶譜變異及其意義

同工酶與基因的表達(dá)有著直接的聯(lián)系,它是生物遺傳多態(tài)性在分子水平上的表現(xiàn),同時(shí)在一定的發(fā)育階段和特定的組織器官中同工酶譜具有相對(duì)穩(wěn)定性,以及有一定的遺傳傳遞規(guī)律,因此被廣泛用于植物遺傳育種,物種起源與分類,雜種優(yōu)勢(shì)預(yù)測(cè)和抗病性篩選,植物生理生化等研究領(lǐng)域[26-31]。

本試驗(yàn)根據(jù)遷移率(Rf)將SOD酶譜分為3個(gè)區(qū):A區(qū)有3條酶帶(Rf=0.27,0.41,0.44),B區(qū)有3條酶帶(Rf=0.55,0.59,0.62),C區(qū)有3條酶帶(Rf=0.68,0.74,0.79),這與前人的研究結(jié)果是一致的。而且親本及其后代群體SOD同工酶譜帶基本相同,有很大的同源性,差異僅體現(xiàn)在某些酶帶的強(qiáng)弱上。POD同工酶電泳顯示有13條酶帶,可分為3個(gè)區(qū),Ⅰ區(qū)(慢區(qū))、Ⅲ區(qū)(快區(qū))酶帶數(shù)多,顯色深,酶帶寬;而Ⅱ區(qū)(中區(qū))酶帶數(shù)少,顯色淺,酶帶窄;其中3、5、11、13號(hào)酶帶出現(xiàn)頻率為100%,1、2、10號(hào)酶帶出現(xiàn)頻率均在85%以上,8種材料共出現(xiàn)78條酶帶,平均每種材料9.8條。說(shuō)明花莢期POD同工酶酶帶比較豐富,組合形式上也表現(xiàn)各異。不同品系的POD同工酶酶譜之間既具有某些相同點(diǎn),又各自具有自身的特征,反映出不同品系之間的相互聯(lián)系與區(qū)別。

試驗(yàn)證明,POD同工酶指標(biāo)在大豆品系中更具有多態(tài)性,體現(xiàn)了蛋白質(zhì)表達(dá)水平和基因水平上的差異。因此,對(duì)雜交后代不同品系之間POD同工酶的測(cè)定,再輔以其他區(qū)域的差異及酶活力的差異,可為大豆品種親緣關(guān)系的鑒定提供依據(jù),并為篩選優(yōu)良的雜種后代提供一定的科學(xué)依據(jù)。

3.2本試驗(yàn)的局限性

本試驗(yàn)利用同工酶分析技術(shù),從蛋白質(zhì)水平、遺傳背景方面來(lái)尋找大豆不同品系間的親緣關(guān)系,可能會(huì)存在一些不足之處。首先,在進(jìn)行酶提取、活性測(cè)定或電泳染色時(shí)酶活性與酶本身結(jié)構(gòu)會(huì)發(fā)生變化,因此操作結(jié)果的重復(fù)性不強(qiáng)。其次,電泳只能檢測(cè)編碼酶蛋白的基因位點(diǎn),對(duì)非結(jié)構(gòu)基因則無(wú)能為力,所檢測(cè)的位點(diǎn)數(shù)不可能很多(一般少于30個(gè)),一批等位酶位點(diǎn)變異并不一定代表整個(gè)基因組的變異,且有些隱避的變異可能無(wú)法用電泳檢測(cè)到,在要求結(jié)果準(zhǔn)確的檢測(cè)中,同工酶檢測(cè)技術(shù)只能用作輔助手段。再者,與大豆產(chǎn)量相關(guān)的性狀多屬數(shù)量性狀,而同工酶標(biāo)記位點(diǎn)較少,因此難于準(zhǔn)確標(biāo)定。最后,在聚類方法上,本試驗(yàn)采用二態(tài)性狀對(duì)POD酶帶進(jìn)行標(biāo)記,是以酶帶的有無(wú)作為分類的依據(jù),沒(méi)有酶量和酶活性的信息參與,這除了試驗(yàn)本身的限制外,考慮到酶帶的有無(wú)比深淺(表示酶活性和酶量)更具有較大的遺傳差異,前者是質(zhì)的差別,后者是量的差別。當(dāng)然,如果把酶量的差異考慮進(jìn)去,以及結(jié)合多酶系統(tǒng)、形態(tài)特征和分子標(biāo)記的研究,將會(huì)使大豆品種的親緣關(guān)系和多態(tài)性研究更趨于合理和真實(shí),對(duì)大豆的雜交育種也具有較好的指導(dǎo)意義。

4結(jié)論

通過(guò)對(duì)晉大52(母本)、晉大57(父本)及其6個(gè)有代表性的后代SN-1、SN-2、SN-3、SN-4、SN-5、SN-6進(jìn)行花莢期超氧化物歧化酶(SOD) 、過(guò)氧化物酶(POD)兩種同工酶的電泳圖譜檢測(cè)和分析,把SOD酶譜分為3個(gè)區(qū):A區(qū)有3條酶帶(Rf=0.27、0.41、0.44),B區(qū)有3條酶帶(Rf=0.55、0.59、0.62),C區(qū)有3條酶帶(Rf=0.68、0.74、0.79),譜帶基本相同,有很大的同源性;POD同工酶電泳顯示有13條酶帶,其中3、5、11、13號(hào)酶帶出現(xiàn)頻率為100%,1、2、10號(hào)酶帶頻率均在85%以上,8種材料共出現(xiàn)78條酶帶,平均每種材料9.8條,說(shuō)明POD酶帶豐富,變異程度高。利用聚類分析,6個(gè)后代與父本的平均遺傳距離為0.168 6,與母本的平均遺傳距離為0.387 6。

參考文獻(xiàn):

[1] Koehn R K, Hilbish T J. The adaptive importance of genetic variation [J]. Am Sci, 1987, 75 :134 - 141.

[2] 李軍, 錢波, 鄭師章, 等.大豆種子庫(kù)中同工酶水平的遺傳多樣性的初步研究[J].生態(tài)應(yīng)用學(xué)報(bào), 1998, 9(2) :145-149.

[3] 劉昭軍, 李希臣.同工酶技術(shù)在大豆育種中的應(yīng)用[J].黑龍江農(nóng)業(yè)科學(xué), 2001(6):35-37.

[4] 李軍, 陶蕓, 鄭師章, 等.同工酶水平上野生大豆種群內(nèi)分化的研究[J].植物學(xué)報(bào), 1995, 27(9) :669-676.

[5] 鐘珍萍.大豆種間的同工酶分析[J].大豆科學(xué), 1992, 11(4) :56-59.

[6] 張維民, 苗以農(nóng), 姚明昌.不同大豆品種的幾種同工酶初步研究[J].中國(guó)油料, 1988(4):329-335.

[7] 盧翠華, 何志鴻, 宋英淑.黑龍江省主要大豆品種同工酶酶譜分析[J].大豆科學(xué), 1990, 9(2):145-148.

[8] 楊琪, 王金陵, 楊慶凱, 等.大豆同工酶與雜種優(yōu)勢(shì)及產(chǎn)量性狀變異的相關(guān)研究[J]. 大豆科學(xué), 1992, 11(1):1-7.

[9] 虞京葳, 夏文勝, 黃偉英.大豆脂酶同工酶的初步研究[J].大豆科學(xué), 1983, 2(1):104-108.

[10] 裴顏龍, 王嵐, 葛頌, 等.野生大豆遺傳多樣性研究4個(gè)天然居群等位酶水平的分析[J]. 大豆科學(xué), 1996, 15(4):302 -309.

[11] 鐘珍萍, 劉德金, 陳啟鋒, 等.大豆種間的同工酶分析比較[J].大豆科學(xué), 1992, 11(4):329- 335.

[12] Wendel J F, Weeden N F. Visualization and interpretation of plant isozymes [M]. London: Dioscorides Press, 1989:5-45.

[13] 胡能書, 萬(wàn)賢國(guó).同工酶技術(shù)及其應(yīng)用[M].湖南:湖南科學(xué)科技出版社,1985.

[14] 喬富廉.植物生理學(xué)實(shí)驗(yàn)分析測(cè)定技術(shù)[M].北京:中國(guó)農(nóng)業(yè)科學(xué)技術(shù)出版社, 2002.

[15] 王金勝.農(nóng)業(yè)生物化學(xué)技術(shù)[M].太原:山西科學(xué)技術(shù)出版社, 1997:74-76.

[16] 中國(guó)科學(xué)院上海植物生理研究所.現(xiàn)代植物生理學(xué)實(shí)驗(yàn)指南[M].北京:科學(xué)出版社, 1999:314-315.

[17] 秦永梅, 于卓. 幾種高粱-蘇丹草雜交種的光合特性及同工酶酶譜差異[J].華北農(nóng)學(xué)報(bào),2004, 19(3):18-21.

[18] 羅廣華.超氧化物歧化酶在細(xì)胞內(nèi)的分布[J].植物生理學(xué)報(bào), 1985, 11(2):25-27.

[19] 秦國(guó)庫(kù), 苗以農(nóng).大豆不同品種不同節(jié)位葉片超氧化物歧化酶的研究[M]//大豆育種應(yīng)用基礎(chǔ)和技術(shù)的研究進(jìn)展. 南京:江蘇科學(xué)技術(shù)出版社, 1990.

[20] Salin M L, Bridges S M. Isolation and characterization of an Iron-Containing superoxide dismutase from water lily, Nuphar Luteum[J]. Plant physiol, 1982.

[21] 羅廣華.大豆和花生種子超氧化物歧化酶的同工酶研究[J].植物生理學(xué)報(bào), 1984, 10(2):37-40.

[22] Chevreau E, Manganaris A G, Gallet M. Isozyme segregation in five apple progenies and potential use form apconstruction [J].Theor Appl Genet, 1999, 98: 329-336.

[23] 裴喜慶, 薛河儒.SAS及應(yīng)用[M].北京:中國(guó)農(nóng)業(yè)出版社, 1998.

[24] 曲慶云, 趙曉梅.統(tǒng)計(jì)分析方法——SAS實(shí)例精選[M].北京:清華大學(xué)出版社,2004.

[25] 袁志發(fā), 周靜芋.多元統(tǒng)計(jì)分析[M].北京:科學(xué)出版社, 2000.

[26]Jobin-Decor M P , Graham G C , Henry R J , et al. RAPD and isozyme analysis of genetic relationships between Carica papaya and wild relatives[J].Genetic Resources and Crop Evaluation , 1997, 44:471-477.

[27] 趙義, Becker H C.同工酶分子標(biāo)記研究中國(guó)和歐洲栽培的遺傳差異[J].作物學(xué)報(bào), 1998, 24 (2):213-220.

[28] 范傳珠, 王述民, 馬緣生, 等.多花菜豆基因庫(kù)種子繁殖更新后遺傳完整性分析[J]. 中國(guó)農(nóng)學(xué)通報(bào), 1999(4):5-8.

[29] 張宗文.紅花品種資源的同工酶遺傳多樣性及分類研究[J].植物遺傳資源科學(xué), 2000, 1 (4):6-13.

[30] 王中仁.植物等位酶分析[M].北京:科學(xué)出版社, 1996.

[31] Spinosa A, Pignone D. Assessment of genetic variation in a working collection of Vigna vexillata by isozyme and RAPD analysis [J]. Genetic Resources and Crop Evaluation, 1998, 45:374-354.

主站蜘蛛池模板: 色成人综合| 亚洲香蕉在线| 夜夜操狠狠操| 国产日韩精品欧美一区灰| 欧美激情视频一区二区三区免费| 综1合AV在线播放| 亚洲国产成人超福利久久精品| 久久久久人妻精品一区三寸蜜桃| 亚洲国产欧美目韩成人综合| 亚洲一区二区三区中文字幕5566| 狠狠躁天天躁夜夜躁婷婷| 伊人久久久大香线蕉综合直播| 国产精品一区在线观看你懂的| 欧美日韩在线国产| 久久综合结合久久狠狠狠97色| 四虎影院国产| 九月婷婷亚洲综合在线| 欧美午夜视频| 青青青国产精品国产精品美女| 欧美精品啪啪一区二区三区| 中文字幕免费在线视频| 亚洲精品欧美重口| 亚洲精品中文字幕无乱码| 成人午夜在线播放| 国产女人18毛片水真多1| 青草91视频免费观看| 五月天天天色| 国产噜噜在线视频观看| 一级全黄毛片| 2022国产无码在线| 色欲综合久久中文字幕网| 97国产精品视频自在拍| 亚洲日韩精品无码专区| 久久午夜影院| 久久五月视频| 精品乱码久久久久久久| 2022精品国偷自产免费观看| 伊人大杳蕉中文无码| 青青草91视频| 欧美国产综合色视频| 亚洲男人在线| 国产视频 第一页| 日韩免费毛片视频| 精品国产中文一级毛片在线看| 亚洲中文字幕23页在线| 无码中文字幕加勒比高清| 欧美亚洲国产日韩电影在线| 漂亮人妻被中出中文字幕久久| 国产99精品久久| 麻豆精品久久久久久久99蜜桃| 精品成人一区二区三区电影 | 欧美日本中文| 亚洲无卡视频| 1769国产精品免费视频| 一区二区三区精品视频在线观看| 久久精品人人做人人综合试看| 精品国产www| 精品久久久久久中文字幕女| 成人午夜精品一级毛片| 亚洲视频免费在线看| 九九久久精品免费观看| 色网站免费在线观看| 99热在线只有精品| 欧美成人精品一级在线观看| 粉嫩国产白浆在线观看| 九一九色国产| 日本道综合一本久久久88| 国产成人综合亚洲网址| 中文字幕不卡免费高清视频| 亚洲成年网站在线观看| 国产成人精品一区二区三区| 操美女免费网站| 四虎成人精品| 亚洲无码电影| 在线观看国产精美视频| 久久午夜影院| 国产一二视频| 欧美在线黄| 久久精品国产亚洲麻豆| 特级欧美视频aaaaaa| 欧美一级高清片久久99| 呦视频在线一区二区三区|