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

    基于mt DNA COI基因的新疆部分地區(qū)葉蟬的分子鑒定

    2021-02-16 09:40:02張海燕陳光輝古麗扎爾阿不都克力木張秀英王玉濤
    關(guān)鍵詞:葉蟬條形碼遺傳

    張海燕,陳光輝,古麗扎爾·阿不都克力木,張秀英,王玉濤

    (1喀什大學(xué)生命與地理科學(xué)學(xué)院,新疆喀什 844000;2新疆帕米爾高原生物資源與生態(tài)重點(diǎn)實(shí)驗(yàn)室,新疆喀什 8440000)

    0 引言

    葉蟬屬半翅目,角蟬總科(Membracoidea),葉蟬科(Cicadellidea)昆蟲(chóng)[1]。全世界已知 43亞科、2345 屬[2],超過(guò)20000種[3],中國(guó)24亞科[2],約2000種[3]。葉蟬危害植物葉片、根莖部、韌皮組織造成直接損害,部分種類傳播植物病毒[4],蔡平等[5]統(tǒng)計(jì),植物病毒病媒介昆蟲(chóng)世界范圍內(nèi)已知的約397種,其中葉蟬133種,將近1/3,傳播86種病原物[6]。害蟲(chóng)種系的準(zhǔn)確鑒定影響防治方法的有效性[7],但部分葉蟬體型和體色易隨發(fā)育階段或生境改變[8],如廣頭葉蟬Macropsinae若蟲(chóng)的體色、警覺(jué)性隨蟲(chóng)齡變化[9],也說(shuō)明葉蟬在低齡期和若蟲(chóng)期較易防治,但若蟲(chóng)形態(tài)差異較小,鑒定較為困難。DNA條形碼技術(shù)補(bǔ)充了形態(tài)鑒定法的短板,能夠?qū)Σ煌l(fā)育歷期、不同形態(tài)、殘?bào)w昆蟲(chóng)有效鑒定,岳巧云等[10]基于COI基因成功鑒定玉米象Sitophilus zeamais幼蟲(chóng)。崔中翌等[11]基于CO1基因?qū)?43頭果蠅Drosophila melanogaster幼蟲(chóng)或卵成功鑒定。Caterino等[12]利用COI和18S rDNA基因鑒定了伴閻甲亞科Histeridae的幼蟲(chóng)。1993年Fang[13]首次用16SrRNA基因序列法研究了角頂葉蟬Dloctephaliane的系統(tǒng)發(fā)育。Johnson等[14]用16SrRNA、ND1和tRNA基因分析了紅額葉蟬屬Errhomus幾種葉蟬的種間親緣關(guān)系。Palomera、Bluemel[15-16]利用COI基因分別研究了玉米黃翅葉蟬Dalbulusmaidis和黃條脊冠葉蟬Aphrodes leafhoppers的遺傳多樣性。BOLD數(shù)據(jù)庫(kù)中記錄了1973種葉蟬,擁有條形碼的有1414種,4461條葉蟬的條形碼被記錄[17]。戴仁懷等[18]在2008年,基于28S rDNA D2、16S rDNA首次在國(guó)內(nèi)分析研究角頂葉蟬亞科進(jìn)化關(guān)系。倪俊強(qiáng)等[19]基于COⅡ基因分析了閩桂地區(qū)尖凹大葉蟬Bothrogonia acuminata5個(gè)地理種群的進(jìn)化地位。俞鵬飛等[20]采用引物步移法研究白邊大葉蟬Kolla paulula(Walker)線粒體基因組序列特征,發(fā)現(xiàn)白邊大葉蟬屬角蟬總科葉蟬科。喬利等[21]通過(guò)Cytb、COI基因,明確了信陽(yáng)地區(qū)4種葉蟬種類、遺傳關(guān)系。中國(guó)對(duì)于葉蟬的DNA條形碼研究大多集中于對(duì)成蟲(chóng)的研究,對(duì)葉蟬若蟲(chóng)的研究較少。昆蟲(chóng)幼蟲(chóng)的形態(tài)分類是個(gè)難題[22];COI基因進(jìn)化速率快、序列保守,分子標(biāo)記效果較好[23-24]?;诖?,本研究首先通過(guò)形態(tài)分類,對(duì)未知葉蟬分類識(shí)別;提取葉蟬基因組,基于mt DNACOI基因,對(duì)提取的基因組用通用引物PCR擴(kuò)增、測(cè)序;對(duì)目的基因序列用生物信息學(xué)軟件進(jìn)行相似性分析、計(jì)算種間距離、分析遺傳進(jìn)化地位,獲得不同種類葉蟬快速鑒定的DNA條形碼。南疆與多國(guó)接壤,較容易遭生物入侵,作為農(nóng)業(yè)大區(qū),植被種類多,易于葉蟬生存,但區(qū)域內(nèi)葉蟬分類識(shí)別基礎(chǔ)薄弱,研究葉蟬DNA條形碼的報(bào)道更為少見(jiàn)。本研究選用COI基因?qū)?種葉蟬若蟲(chóng)和2種成蟲(chóng)進(jìn)行分類鑒定,可豐富邊境地區(qū)葉蟬基因數(shù)據(jù)庫(kù)、為維護(hù)生物安全和害蟲(chóng)有效防治提供依據(jù)。

    1 材料與方法

    1.1 標(biāo)本來(lái)源

    本研究葉蟬樣本分布于新疆部分農(nóng)田,用捕蟲(chóng)網(wǎng)掃捕、黃板誘捕獲得(表1)。首先形態(tài)分類鑒定,初步分為6種葉蟬標(biāo)本,標(biāo)記為A、B、C、D、E、F,在離心管中加入99%的無(wú)水乙醇,將標(biāo)本放入離心管,冷藏于4℃冰箱中,實(shí)驗(yàn)于2020年5—10月在新疆帕米爾高原生物資源與生態(tài)重點(diǎn)實(shí)驗(yàn)室開(kāi)展。

    表1 葉蟬標(biāo)本信息

    1.2 儀器設(shè)備

    體視顯微鏡(Motic Cam2506 SMZ168)、離心機(jī)(BECKMAN COULTERTM AllegraTM X-22R Centrifuge)、FM100 雪花制冰機(jī)(YKKY)、振蕩器(IKA MS3 digital)、恒溫水浴鍋、PCR儀(Veriti TM 96-Well PCR 4375786新加坡制造)、伯樂(lè)Bio-RadPowerpacHC電泳設(shè)備、凝膠成像儀(Bio-Rad Molecular Imager Gel Doc XR)等。

    1.3 提取葉蟬DNA基因組

    由于葉蟬若蟲(chóng)體型較小,為保證基因組提取效果,參照鄒志文等Chelex100改進(jìn)法提取葉蟬DNA[25]。

    (1)取1只葉蟬,用雙蒸水在離心管中將樣本清洗數(shù)次,后將樣本置于濾紙上數(shù)分鐘直至干燥。

    (2)將干燥的葉蟬樣本放入離心管中(1.5 mL),對(duì)玻璃棒高壓滅菌后緊貼管壁充分研磨樣本。

    (3)吸取0.5 mL雙蒸水(已滅菌),加入裝有葉蟬樣本的離心管中,充分渦旋溶液至均勻,放入-70℃冰箱中冷卻4 min。

    (4)將離心機(jī)參數(shù)設(shè)置為12000 r/min,將溶液離心5 min,緩慢抽取上清棄去,只留下沉淀,重復(fù)此步驟3次。

    (5)將5%的Chelex-100溶液搖勻,抽取120 μL,加入到沉淀中,調(diào)整恒溫水浴鍋的溫度為56℃,將溶液搖勻后加熱2 h。

    (6)渦旋5~10 s,將離心管封口(防止溶液噴濺)放入100℃的恒溫水浴鍋,加熱10 min。

    (7)渦旋溶液 5~10 s,將離心機(jī)參數(shù)設(shè)置為10000 r/min,將溶液離心5 min,提取離心所得上清液,在-20℃條件下保存,作為DNA模板。

    1.4 基因組PCR擴(kuò)增

    引物由北京奧科生物公司合成,上游引物L(fēng)CO1490:5’-GGTCAACAAATCATAAAGATATTGG-3’,下游引物 HCO2198:5’-TAAACTTCAGGGTGAC CAAAAAATCA-3’[26-28]。使用通用引物對(duì)COI序列PCR擴(kuò)增,擴(kuò)增體系總體積設(shè)置為21 μL,其中2×Taq酶7.2 μL,ddH2O 12 μL,上下游引物各為0.4 μL,DNA模板1 μL。反應(yīng)條件及步驟為:95℃預(yù)變性5 min,95℃變性45 s,52℃退火45 s,72℃延伸1 min,循環(huán)40次,循環(huán)結(jié)束后72℃延伸7 min,4℃條件下保存。PCR產(chǎn)物送蘇州金唯智生物科技有限公司測(cè)序。

    1.5 處理數(shù)據(jù)、構(gòu)建系統(tǒng)進(jìn)化樹(shù)

    將擴(kuò)增的目的基因送專業(yè)測(cè)序公司測(cè)序,經(jīng)測(cè)序儀生成abi文件,首先用Chromas 2.6.6軟件打開(kāi)測(cè)序獲得的基因序列的峰圖,查找誤讀的堿基手動(dòng)校正。用DNAMAN4.0軟件進(jìn)行序列正反鏈拼接匹配校正。用NCBI數(shù)據(jù)庫(kù)的BLAST程序比較序列同源性,確定所測(cè)序列屬昆蟲(chóng)COI基因序列。下載Score較高、Query cover大于90%、Per ident高于80%的序列,用ClustalX 1.83軟件比對(duì)下載序列與測(cè)序所得序列。用MEGA7.0.14統(tǒng)計(jì)序列堿基組成、變異位點(diǎn)、核苷酸組成。選用Kimura 2-Parameter模型計(jì)算種間遺傳差異,用DAMBE7.2.136基于P距離分析序列堿基替換飽和性、檢測(cè)系統(tǒng)發(fā)育信號(hào)。用NJ鄰接法(Neighbour-Jioning)構(gòu)建系統(tǒng)進(jìn)化樹(shù),循環(huán)1000次,估計(jì)進(jìn)化樹(shù)中節(jié)點(diǎn)的置信度(bootstrap confidence level,BCL)。

    2 結(jié)果與分析

    2.1 葉蟬標(biāo)本初步形態(tài)鑒定

    參照表2中葉蟬成蟲(chóng)及若蟲(chóng)的形態(tài)特征及6個(gè)葉蟬樣本的外部形態(tài),初步劃分樣本所屬類群,經(jīng)分析所研究樣本分屬6個(gè)亞科,6種葉蟬背面、側(cè)面、腹面形態(tài)特征見(jiàn)圖1。

    圖1 6種葉蟬形態(tài)特征圖

    表2 6種葉蟬形態(tài)特征

    續(xù)表2

    2.2 序列分析

    本研究獲得葉蟬核酸序列共6條,其中4種葉蟬若蟲(chóng)(樣本A、B、C、D)的4條序列屬首次獲得,2種葉蟬成蟲(chóng)(樣本E、F)的2條COI序列在新疆首次報(bào)道。經(jīng)NCBI數(shù)據(jù)庫(kù)中的BLAST程序,發(fā)現(xiàn)與該序列同源性最高的是昆蟲(chóng)的COI基因序列,確定所測(cè)序列為COI基因序列。下載GeneBank中與該COI序列同源性較高的葉蟬的序列,將測(cè)序獲得序列與下載序列用ClustalX1.83軟件比對(duì)分析,非保守區(qū)域去除,測(cè)序拼接后僅以表6中的序列長(zhǎng)度進(jìn)行分析。COI基因核苷酸組成見(jiàn)表3,本研究6個(gè)樣本的COI序列的A+T含量分別為65%、72%、67%、69%、67%、68%,COI基因的A+T含量高于G+C,A+T堿基偏移性明顯,與昆蟲(chóng)mt DNA基因堿基構(gòu)成特點(diǎn)一致,成功對(duì)葉蟬樣本的COI基因序列在分子水平上監(jiān)測(cè)分析。

    表3 mt DNA核苷酸堿基組成表

    2.3 分子鑒定結(jié)果分析

    將獲得的葉蟬樣本的COI序列在GenBank中進(jìn)行比對(duì),可獲得相似度高低排列表,結(jié)果顯示均和葉蟬相關(guān)序列的相似性最高。COI序列覆蓋度大于90%、相似度大于95%,被認(rèn)定為達(dá)到鑒定標(biāo)準(zhǔn)[39]。COI序列相似度大于98%被認(rèn)為達(dá)到鑒定標(biāo)準(zhǔn)或水平[40]。本研究獲得的6條COI基因序列中,3個(gè)葉蟬樣本與NCBI中相關(guān)序列相似度高達(dá)98%,大于或等于鑒定標(biāo)準(zhǔn)98%,可以鑒定為該類昆蟲(chóng),鑒定依據(jù)和結(jié)果如表4所示。樣本A與KR564712.1條沙葉蟬Psammotettix confinis COI基因的相似度為99.02%,確定為條沙葉蟬。樣本B與MF716879.1大青葉蟬Cicadella viridis COI基因的相似度為99.09%%,確定為大青葉蟬。樣本E與JQ755806.1多態(tài)廣頭葉蟬Macropsis notate COI基因的相似度為100%,確定為多態(tài)廣頭葉蟬。樣本C與HQ929177.1紅帶鏟頭葉蟬Hecalus arcuatus的相似度只有85.5%,未達(dá)到鑒定標(biāo)準(zhǔn),經(jīng)形態(tài)特征鑒定初步鑒定為鏟頭葉蟬亞科Hecalinae鏟頭葉蟬屬Hecalus Stal葉蟬,由于該樣本為若蟲(chóng),尚未發(fā)育完全,形態(tài)特征不明顯,無(wú)法將其鑒定到種,故該葉蟬鑒定為鏟頭葉蟬屬葉蟬。樣本D與KR560629.1斑葉蟬族Erythroneura vitifex的相似度只有87.7%,未達(dá)到鑒定標(biāo)準(zhǔn),經(jīng)形態(tài)特征鑒定初步鑒定為阿小葉蟬屬Arboridia Zachvathin葡萄阿小葉蟬Arboridia kakogawana,(斑葉蟬屬Erythroneura和阿小葉蟬屬Arboridia均屬斑葉蟬族)。樣本F與JF737032.1六點(diǎn)葉蟬Macrosteles sexnotatus的相似度最高,達(dá)到91.92%,但未達(dá)到鑒定標(biāo)準(zhǔn),經(jīng)形態(tài)特征鑒定初步鑒定為殃葉蟬亞科Euscelinae二叉葉蟬屬M(fèi)acrostelesFieber六點(diǎn)葉蟬Macrosteles sexnotatus,形態(tài)特征與分子鑒定結(jié)果一致,該葉蟬為六點(diǎn)葉蟬。本研究中葉蟬樣本的形態(tài)特征鑒定結(jié)果與分子鑒定結(jié)果基本一致,樣本C、D、F的COI基因序列相似度與NCBI中相關(guān)葉蟬的相似度存在一定的差距,與樣本的蟲(chóng)態(tài)、新疆對(duì)葉蟬的條形碼研究較少數(shù)據(jù)庫(kù)中數(shù)據(jù)不夠全面有關(guān)。

    表4 葉蟬樣本COI基因鑒定依據(jù)及結(jié)果

    2.4 堿基替換飽和性分析

    啟動(dòng)DAMBE7.2.136軟件,P距離位于橫軸,轉(zhuǎn)換(s)和顛換(v)位于縱軸,繪制函數(shù)曲線圖,以此檢測(cè)堿基替換是否飽和[41],詳見(jiàn)圖2。由圖2可知COI序列組的轉(zhuǎn)換值隨著P距離增加而升,與P距離存在較強(qiáng)的線性關(guān)系,可知本研究中序列堿基替換未飽和,可進(jìn)行系統(tǒng)進(jìn)化地位分析[41]。

    圖2 葉蟬COI基因序列堿基替換飽和性

    2.5 葉蟬系統(tǒng)發(fā)育分析

    從NCBI數(shù)據(jù)庫(kù)篩選下載覆蓋度大于90%、同源性80%以上的葉蟬科昆蟲(chóng)的COI基因序列,與實(shí)驗(yàn)所得序列,用MEGA7.0.14構(gòu)建COI葉蟬部分屬種NJ分子進(jìn)化樹(shù)(圖3)。本研究分析的6種葉蟬在系統(tǒng)發(fā)育樹(shù)上的位置清晰,分類明確,各種在各分類階元間形成各自獨(dú)立的單系。本研究樣本的CO1序列和葉蟬類昆蟲(chóng)相關(guān)序列分為2個(gè)大支,屬間以亞科為分支聚類為6個(gè)小支,和已知葉蟬相關(guān)序列以較高置信值聚在一起,明顯將6種葉蟬準(zhǔn)確區(qū)分開(kāi)來(lái)。阿小葉蟬屬位于進(jìn)化樹(shù)的根部,是葉蟬科最為原始的類群[42]。

    圖3 與葉蟬相關(guān)的部分種類的COI基因系統(tǒng)進(jìn)化NJ樹(shù)(數(shù)字表示置信度,用紅色菱形標(biāo)注樣本)

    2.6 遺傳距離分析

    用MEGA 7.0.14軟件,基于K2P距離模型,計(jì)算6種葉蟬不同個(gè)體的COI基因序列間的遺傳距離,如表5所示。6種葉蟬種間遺傳距離為0.183~0.265,種間遺傳距離平均數(shù)為0.232;鑒定到種的樣本A、B、E種內(nèi)遺傳距離分別為0.0037、0.0019、0.0037,遠(yuǎn)低于0.01。種內(nèi)平均遺傳距離小于0.01,種間平均遺傳距離大于0.2,符合“10倍閾值定律”[40]且與形態(tài)分類結(jié)論一致[43],對(duì)6種葉蟬成功分類鑒定,因此,這些mt DNA序列可以作為鑒定葉蟬科昆蟲(chóng)的DNA條形碼[44]。

    表5 6種葉蟬COI基因種間遺傳距離

    3 討論

    3.1 葉蟬mt DNA堿基組成特征

    本研究中6個(gè)樣本的COI基因序列A+T含量明顯高于G+C,表現(xiàn)出明顯的A/T堿基偏嗜,與昆蟲(chóng)線粒體基因序列堿基組成特點(diǎn)一致[45]。高婭蓉等研究小葉蟬亞科5個(gè)族的COI基因片段序列時(shí)得出小葉蟬亞科5個(gè)族A+T含量為69.2%,含量較高[46],與本研究結(jié)果一致。確定了mt DNACOI基因序列可作為葉蟬的DNA條形碼,獲得COI基因序列6條,其中4條葉蟬若蟲(chóng)的COI基因序列為葉蟬科昆蟲(chóng)首次報(bào)道,2條葉蟬成蟲(chóng)序列為新疆首次報(bào)道。

    3.2 葉蟬mt DNA分子鑒定結(jié)果

    本研究中樣本C、D、F的COI基因序列與NCBI數(shù)據(jù)庫(kù)中相關(guān)葉蟬的COI基因序列相似度之間有一定差距,未達(dá)到鑒定水平。經(jīng)查證葡萄阿小葉蟬COI基因序列在各類網(wǎng)絡(luò)數(shù)據(jù)庫(kù)中均未記錄,故本研究中獲得的樣本D葡萄小葉蟬的COI基因序列與斑葉蟬族葉蟬相似度較低。另外新疆溫差較大、干燥少雨、光照強(qiáng),地理環(huán)境特殊,區(qū)域內(nèi)相關(guān)昆蟲(chóng)研究較少,獲得的數(shù)據(jù)未在條形碼數(shù)據(jù)庫(kù)中登記,造成新疆葉蟬科昆蟲(chóng)在GenBank數(shù)據(jù)庫(kù)中序列不全,相似度低;另一方面在實(shí)驗(yàn)過(guò)程中因機(jī)器或人工的中間環(huán)節(jié)較多,影響實(shí)驗(yàn)結(jié)果,而給鑒定結(jié)果造成一定誤差,數(shù)據(jù)被誤傳到數(shù)據(jù)庫(kù),造成同一物種相似度存在差異。

    3.3 葉蟬mt DNA基因的遺傳距離

    研究發(fā)現(xiàn)絕大多數(shù)物種間的COI基因序列之間存在相對(duì)較小的種內(nèi)遺傳距離和相對(duì)較大的種間遺傳距離,98%以上的種間遺傳距離大于2%,而種內(nèi)遺傳距離大多低于1%,很少超過(guò)2%[47-48]。本研究6種葉蟬種間遺傳距離平均數(shù)為0.232,樣本A、B、E種內(nèi)遺傳距離平均數(shù)為0.0031,遠(yuǎn)低于1%,達(dá)到鑒定標(biāo)準(zhǔn)。樣本E與F的遺傳距離最大為0.265,樣本A與C的遺傳距離最小為0.183,來(lái)自吐魯番的樣本D與來(lái)自喀什、克州的其他樣本的地理距離最遠(yuǎn),但遺傳距離并不是最大的,田小娟指出地理間距和遺傳差異之間無(wú)相關(guān)性[49],地理距離不會(huì)必然導(dǎo)致種群間的遺傳分化[49]。李樂(lè)對(duì)12種假眼小綠葉蟬的研究發(fā)現(xiàn),全國(guó)不同地理位置、海拔和氣候的假眼小綠葉蟬種群間遺傳變異不明顯[50],與本研究結(jié)果一致。本研究中地理環(huán)境并不是引起葉蟬遺傳結(jié)構(gòu)差異的主要原因,但本研究采集的樣本地域覆蓋度還較低,不能代表該區(qū)域的葉蟬遺傳結(jié)構(gòu)與地理環(huán)境的關(guān)系,今后將擴(kuò)大邊境區(qū)域及周邊采樣范圍,以期更全面評(píng)估該區(qū)域葉蟬遺傳多樣性、分布特點(diǎn)等。

    4 結(jié)論

    本研究基于mt DNACOI基因,首先通過(guò)形態(tài)分類法確定葉蟬的大致類群,再利用DNA條形碼技術(shù),獲得葉蟬COI基因序列,比對(duì)序列相似性,分析堿基含量、檢測(cè)系統(tǒng)發(fā)育信號(hào)、計(jì)算遺傳距離、構(gòu)建系統(tǒng)進(jìn)化關(guān)系,確定了6條mt DNACOI基因序列可作為葉蟬快速鑒定的DNA條形碼,其中4種葉蟬若蟲(chóng)的基因序列為葉蟬科昆蟲(chóng)首次研究報(bào)道。本研究獲得如下結(jié)論:(1)經(jīng)形態(tài)特征鑒定,6個(gè)樣本分屬6個(gè)亞科、6個(gè)屬、6種葉蟬;(2)葉蟬的mt DNA基因序列構(gòu)成特點(diǎn)與昆蟲(chóng)的mt DNA序列堿基組成特點(diǎn)一致;(3)通過(guò)mt DNA序列對(duì)比、遺傳距離分析、系統(tǒng)發(fā)育構(gòu)建,確定了本研究6種葉蟬的種類,與形態(tài)分類結(jié)果一致,獲得基因序列可作為本區(qū)域6種葉蟬快速鑒定的DNA條形碼。

    猜你喜歡
    葉蟬條形碼遺傳
    非遺傳承
    葉蟬預(yù)防早 來(lái)年葡萄生長(zhǎng)好
    創(chuàng)意條形碼
    還有什么會(huì)遺傳?
    還有什么會(huì)遺傳
    還有什么會(huì)遺傳?
    從條形碼到二維碼
    從條形碼到二維碼
    基于COI基因序列片段的柔突葉蟬屬部分近緣種分子鑒定與系統(tǒng)發(fā)育
    條形碼大變身
    国产美女午夜福利| 亚洲欧美日韩卡通动漫| 国产精品免费一区二区三区在线| 人体艺术视频欧美日本| 久久精品国产亚洲av天美| 男女视频在线观看网站免费| 日韩精品有码人妻一区| 你懂的网址亚洲精品在线观看 | 亚洲美女视频黄频| 赤兔流量卡办理| 极品教师在线视频| 亚洲高清免费不卡视频| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 观看美女的网站| 日本五十路高清| 天天一区二区日本电影三级| 亚洲人成网站在线播| 午夜福利在线观看吧| 国产精品,欧美在线| 国产精品久久久久久久电影| 日韩视频在线欧美| 国产麻豆成人av免费视频| 91久久精品国产一区二区三区| 黄片无遮挡物在线观看| 国产成人freesex在线| 99久久人妻综合| 久久九九热精品免费| 一级二级三级毛片免费看| 亚洲乱码一区二区免费版| 老女人水多毛片| 久久草成人影院| 国内精品宾馆在线| 久久久精品欧美日韩精品| 色噜噜av男人的天堂激情| 亚洲电影在线观看av| 中文字幕av在线有码专区| 九九爱精品视频在线观看| 国产精品久久久久久精品电影小说 | 国产高清视频在线观看网站| 久久精品夜夜夜夜夜久久蜜豆| 日本五十路高清| 中文资源天堂在线| 最近最新中文字幕大全电影3| 日韩高清综合在线| 国产一区二区亚洲精品在线观看| eeuss影院久久| 美女脱内裤让男人舔精品视频 | 午夜视频国产福利| 男人狂女人下面高潮的视频| 中国美女看黄片| 床上黄色一级片| 少妇人妻一区二区三区视频| .国产精品久久| 久久久精品大字幕| 婷婷色av中文字幕| 国产成人福利小说| 一区二区三区高清视频在线| 少妇高潮的动态图| 色尼玛亚洲综合影院| 欧美性猛交黑人性爽| 精品久久久久久久久亚洲| 亚洲最大成人中文| 在线播放国产精品三级| 久久午夜亚洲精品久久| 国产精品女同一区二区软件| 婷婷色综合大香蕉| 亚洲第一区二区三区不卡| 久久久久久久久久成人| 又黄又爽又刺激的免费视频.| 久久草成人影院| 99久国产av精品国产电影| а√天堂www在线а√下载| 亚洲一级一片aⅴ在线观看| 国产午夜福利久久久久久| 尾随美女入室| 精品国内亚洲2022精品成人| 成人鲁丝片一二三区免费| 欧美成人一区二区免费高清观看| 色哟哟·www| 久久婷婷人人爽人人干人人爱| 国产精品综合久久久久久久免费| 日韩国内少妇激情av| 国产午夜精品一二区理论片| av在线老鸭窝| 国产91av在线免费观看| 亚洲成人久久爱视频| 一个人观看的视频www高清免费观看| 亚洲美女视频黄频| 欧美三级亚洲精品| 熟妇人妻久久中文字幕3abv| 精品久久久噜噜| 国产v大片淫在线免费观看| 一本一本综合久久| 女人十人毛片免费观看3o分钟| 麻豆国产97在线/欧美| 亚洲精品成人久久久久久| 免费看a级黄色片| 欧美日韩综合久久久久久| 亚洲在久久综合| 亚洲精华国产精华液的使用体验 | 黄色日韩在线| 日韩欧美在线乱码| 亚洲欧美日韩高清在线视频| 男人狂女人下面高潮的视频| 最近的中文字幕免费完整| 午夜福利视频1000在线观看| 免费搜索国产男女视频| 我要搜黄色片| 免费看美女性在线毛片视频| 国产国拍精品亚洲av在线观看| 国产成年人精品一区二区| 亚洲内射少妇av| 国产黄a三级三级三级人| 亚洲综合色惰| 免费不卡的大黄色大毛片视频在线观看 | 亚洲av免费在线观看| 免费电影在线观看免费观看| 亚洲av成人av| 热99在线观看视频| 国产精品福利在线免费观看| 12—13女人毛片做爰片一| 国产精品一区二区在线观看99 | 大型黄色视频在线免费观看| 在线观看免费视频日本深夜| 日本熟妇午夜| 亚洲无线观看免费| 国产午夜精品一二区理论片| 啦啦啦观看免费观看视频高清| 高清在线视频一区二区三区 | 久久久久久大精品| 久久婷婷人人爽人人干人人爱| 免费观看a级毛片全部| 亚洲精品自拍成人| 欧美日韩一区二区视频在线观看视频在线 | 精品人妻视频免费看| 18+在线观看网站| 日韩在线高清观看一区二区三区| a级一级毛片免费在线观看| 日本一本二区三区精品| 乱码一卡2卡4卡精品| 国产乱人偷精品视频| 真实男女啪啪啪动态图| 精品日产1卡2卡| 天堂影院成人在线观看| 国产高清三级在线| 少妇的逼水好多| 欧美最黄视频在线播放免费| 乱码一卡2卡4卡精品| 国产在线精品亚洲第一网站| 日韩欧美国产在线观看| 国产成人精品一,二区 | www日本黄色视频网| avwww免费| 日日撸夜夜添| 日韩精品有码人妻一区| 菩萨蛮人人尽说江南好唐韦庄 | 欧美日韩在线观看h| 亚洲国产精品合色在线| 国产又黄又爽又无遮挡在线| 久久九九热精品免费| 欧美性感艳星| 国产成人精品久久久久久| 国产亚洲av片在线观看秒播厂 | 亚洲天堂国产精品一区在线| 一区二区三区四区激情视频 | 亚洲国产欧美在线一区| 亚洲av中文av极速乱| 国产av麻豆久久久久久久| 韩国av在线不卡| 国产精品久久久久久精品电影小说 | av国产免费在线观看| 免费看光身美女| 亚洲图色成人| 少妇人妻精品综合一区二区 | 有码 亚洲区| 一区福利在线观看| 国内精品美女久久久久久| 赤兔流量卡办理| 老师上课跳d突然被开到最大视频| 91精品国产九色| 免费观看的影片在线观看| 最近视频中文字幕2019在线8| 一区福利在线观看| 国产免费男女视频| 男女那种视频在线观看| 99热这里只有是精品在线观看| 免费搜索国产男女视频| 哪里可以看免费的av片| av免费观看日本| 亚洲无线观看免费| 久久人妻av系列| 婷婷色av中文字幕| 日本一二三区视频观看| 亚洲无线在线观看| 丰满的人妻完整版| 搞女人的毛片| 成熟少妇高潮喷水视频| 听说在线观看完整版免费高清| 99热只有精品国产| 最好的美女福利视频网| 亚洲中文字幕一区二区三区有码在线看| 精品久久久久久久久久久久久| 国产淫片久久久久久久久| 青春草亚洲视频在线观看| 丰满的人妻完整版| 三级国产精品欧美在线观看| 成人高潮视频无遮挡免费网站| 国产综合懂色| 青春草国产在线视频 | 国产毛片a区久久久久| 欧美潮喷喷水| 天堂影院成人在线观看| 国产69精品久久久久777片| 精品一区二区三区视频在线| 精品一区二区免费观看| 国产亚洲5aaaaa淫片| av在线蜜桃| 免费av不卡在线播放| 国产伦理片在线播放av一区 | 国产中年淑女户外野战色| 久久久色成人| 一本久久中文字幕| 日韩人妻高清精品专区| 在线a可以看的网站| 亚洲丝袜综合中文字幕| 欧美极品一区二区三区四区| 99热全是精品| 麻豆国产av国片精品| 免费人成视频x8x8入口观看| 精品午夜福利在线看| 此物有八面人人有两片| 成年女人看的毛片在线观看| 精品少妇黑人巨大在线播放 | 久久国产乱子免费精品| 日本免费一区二区三区高清不卡| 精品不卡国产一区二区三区| 国产成人一区二区在线| 国产精品麻豆人妻色哟哟久久 | av在线天堂中文字幕| 免费不卡的大黄色大毛片视频在线观看 | 美女高潮的动态| 在线观看午夜福利视频| 久久久久久国产a免费观看| 国产精品一区二区三区四区久久| 少妇人妻一区二区三区视频| 自拍偷自拍亚洲精品老妇| 欧美日韩综合久久久久久| 国产黄片美女视频| 两性午夜刺激爽爽歪歪视频在线观看| 精品久久久久久久人妻蜜臀av| 大又大粗又爽又黄少妇毛片口| 免费黄网站久久成人精品| 三级毛片av免费| 国产成人aa在线观看| 可以在线观看的亚洲视频| 久久韩国三级中文字幕| 免费大片18禁| 乱人视频在线观看| 天堂√8在线中文| 91久久精品国产一区二区三区| 久久精品影院6| 九九热线精品视视频播放| 国产精品av视频在线免费观看| 国产一区二区在线观看日韩| 亚洲精品久久久久久婷婷小说 | 久久久午夜欧美精品| 欧美一区二区精品小视频在线| 插逼视频在线观看| 熟女电影av网| 亚洲一区二区三区色噜噜| 高清午夜精品一区二区三区 | 22中文网久久字幕| 12—13女人毛片做爰片一| 亚洲欧美精品自产自拍| 欧美三级亚洲精品| 国产亚洲精品久久久久久毛片| 一级黄色大片毛片| 永久网站在线| 91精品一卡2卡3卡4卡| 啦啦啦韩国在线观看视频| 欧美变态另类bdsm刘玥| 午夜视频国产福利| 亚洲国产精品国产精品| 六月丁香七月| 欧美激情在线99| 成人特级av手机在线观看| 国产av不卡久久| 国产成人福利小说| 亚洲av免费在线观看| 在线免费观看的www视频| 丰满人妻一区二区三区视频av| 男女下面进入的视频免费午夜| 久久精品国产自在天天线| 欧美日韩在线观看h| 国产av在哪里看| 国产极品天堂在线| 少妇被粗大猛烈的视频| 亚洲欧美成人综合另类久久久 | 一个人免费在线观看电影| 成年女人永久免费观看视频| 久久久久久久久久黄片| 亚洲18禁久久av| 哪里可以看免费的av片| 欧美性猛交黑人性爽| 我要看日韩黄色一级片| 在线免费观看的www视频| 国产精品人妻久久久影院| 久久久久久久久久黄片| 99久久久亚洲精品蜜臀av| 国产乱人视频| av在线蜜桃| 亚洲激情五月婷婷啪啪| 国内少妇人妻偷人精品xxx网站| 亚洲精品自拍成人| 亚洲国产精品sss在线观看| 欧美日本亚洲视频在线播放| 最好的美女福利视频网| 91aial.com中文字幕在线观看| 亚洲成人精品中文字幕电影| 九九在线视频观看精品| 精品不卡国产一区二区三区| 亚洲欧美精品自产自拍| 亚洲电影在线观看av| 人妻夜夜爽99麻豆av| 国产爱豆传媒在线观看| 天天躁夜夜躁狠狠久久av| 亚洲av中文av极速乱| 波多野结衣高清无吗| 亚洲精品日韩av片在线观看| 26uuu在线亚洲综合色| 日本五十路高清| 成人国产麻豆网| 国内精品美女久久久久久| 在线观看免费视频日本深夜| 一级黄色大片毛片| 看片在线看免费视频| 亚洲无线观看免费| 国产探花极品一区二区| 国产黄片美女视频| 三级男女做爰猛烈吃奶摸视频| 亚洲国产高清在线一区二区三| 亚洲一区二区三区色噜噜| 国产精品一区二区性色av| 午夜福利在线在线| 97在线视频观看| 久久精品夜色国产| 少妇的逼水好多| 国产高清不卡午夜福利| 亚洲第一区二区三区不卡| 12—13女人毛片做爰片一| 91aial.com中文字幕在线观看| 久久热精品热| 黄片无遮挡物在线观看| 久久亚洲精品不卡| 久久久久久大精品| 深夜精品福利| 国产一区二区在线av高清观看| 村上凉子中文字幕在线| 亚洲国产精品sss在线观看| 精品国产三级普通话版| 热99在线观看视频| 热99在线观看视频| 国内揄拍国产精品人妻在线| 久久久久久久久久黄片| or卡值多少钱| 99九九线精品视频在线观看视频| 深爱激情五月婷婷| 性插视频无遮挡在线免费观看| 一个人观看的视频www高清免费观看| 久久韩国三级中文字幕| 日韩欧美精品免费久久| 久久午夜亚洲精品久久| 欧美色欧美亚洲另类二区| 一卡2卡三卡四卡精品乱码亚洲| 尤物成人国产欧美一区二区三区| 亚洲三级黄色毛片| 免费av观看视频| 久久久久久久久久久免费av| 夜夜看夜夜爽夜夜摸| 少妇人妻一区二区三区视频| 菩萨蛮人人尽说江南好唐韦庄 | 国国产精品蜜臀av免费| 欧美性感艳星| 午夜福利在线观看吧| 可以在线观看的亚洲视频| 国产精品伦人一区二区| 国产大屁股一区二区在线视频| 在线播放无遮挡| 国语自产精品视频在线第100页| 亚洲av中文av极速乱| eeuss影院久久| 又粗又硬又长又爽又黄的视频 | 亚洲真实伦在线观看| 国产精品免费一区二区三区在线| 欧美一区二区国产精品久久精品| 免费看美女性在线毛片视频| 中国美女看黄片| 美女 人体艺术 gogo| 久久精品国产亚洲av天美| 欧美日本视频| 成人亚洲欧美一区二区av| 欧美精品一区二区大全| 亚洲成人中文字幕在线播放| 国产精品一及| 听说在线观看完整版免费高清| 欧美zozozo另类| 午夜精品在线福利| 国产视频内射| 国产综合懂色| 日本熟妇午夜| 亚洲欧美精品自产自拍| 18禁黄网站禁片免费观看直播| 黑人高潮一二区| 一夜夜www| 春色校园在线视频观看| 精品久久久久久成人av| 色综合色国产| 国产真实伦视频高清在线观看| 菩萨蛮人人尽说江南好唐韦庄 | 老女人水多毛片| 成人综合一区亚洲| 国语自产精品视频在线第100页| 国产三级中文精品| 精品一区二区免费观看| 97在线视频观看| 丝袜美腿在线中文| 久久九九热精品免费| 国产一级毛片在线| 亚洲av中文字字幕乱码综合| 欧美日韩国产亚洲二区| 亚洲成人av在线免费| 男女下面进入的视频免费午夜| 99久久中文字幕三级久久日本| 亚洲七黄色美女视频| 国产黄色小视频在线观看| 给我免费播放毛片高清在线观看| 99久国产av精品| 看片在线看免费视频| 久久久久久久久久久免费av| 亚洲精品自拍成人| 国产三级中文精品| 最近最新中文字幕大全电影3| 色5月婷婷丁香| 黄色日韩在线| 村上凉子中文字幕在线| 国产成年人精品一区二区| 久久午夜福利片| 亚洲美女搞黄在线观看| 人妻夜夜爽99麻豆av| 综合色丁香网| 亚洲人成网站高清观看| 精品人妻视频免费看| 久久九九热精品免费| 亚洲欧美日韩东京热| 亚洲精品色激情综合| 女人被狂操c到高潮| 精品少妇黑人巨大在线播放 | 在线免费观看不下载黄p国产| 国产在线精品亚洲第一网站| 一个人看的www免费观看视频| 亚洲乱码一区二区免费版| 亚洲精品乱码久久久久久按摩| 91aial.com中文字幕在线观看| 亚洲av免费在线观看| 亚洲自偷自拍三级| 一区二区三区免费毛片| 国产91av在线免费观看| 免费看光身美女| 欧美日韩在线观看h| 一级二级三级毛片免费看| 能在线免费观看的黄片| 国产欧美日韩精品一区二区| 熟女电影av网| 久久久成人免费电影| 亚洲精品影视一区二区三区av| 嘟嘟电影网在线观看| 国产免费男女视频| 黄色一级大片看看| 亚洲乱码一区二区免费版| 观看美女的网站| 精品久久久噜噜| 欧美bdsm另类| 在线免费观看不下载黄p国产| 在线a可以看的网站| 18禁在线播放成人免费| 亚洲在线自拍视频| 国产不卡一卡二| 国产麻豆成人av免费视频| 亚洲av成人精品一区久久| 国内精品美女久久久久久| 又粗又硬又长又爽又黄的视频 | 国产成人精品一,二区 | 91午夜精品亚洲一区二区三区| 国产成人精品婷婷| 亚洲欧美日韩高清专用| 国产麻豆成人av免费视频| 美女被艹到高潮喷水动态| 国产免费一级a男人的天堂| 亚洲国产精品sss在线观看| 欧美一区二区亚洲| 一级av片app| 免费av不卡在线播放| 色哟哟哟哟哟哟| АⅤ资源中文在线天堂| 少妇人妻精品综合一区二区 | 国产精品日韩av在线免费观看| 中文字幕精品亚洲无线码一区| 午夜久久久久精精品| 寂寞人妻少妇视频99o| 国产av不卡久久| 日本撒尿小便嘘嘘汇集6| 午夜免费男女啪啪视频观看| 人人妻人人澡欧美一区二区| 日本免费a在线| 女的被弄到高潮叫床怎么办| 欧美激情久久久久久爽电影| 久久亚洲精品不卡| 六月丁香七月| 日本在线视频免费播放| 老熟妇乱子伦视频在线观看| 91aial.com中文字幕在线观看| 偷拍熟女少妇极品色| 啦啦啦啦在线视频资源| 寂寞人妻少妇视频99o| 午夜福利视频1000在线观看| 亚洲中文字幕日韩| 91久久精品国产一区二区成人| 中文字幕制服av| 国产白丝娇喘喷水9色精品| 草草在线视频免费看| 免费人成视频x8x8入口观看| 别揉我奶头 嗯啊视频| 伊人久久精品亚洲午夜| 日韩av不卡免费在线播放| av在线天堂中文字幕| 国模一区二区三区四区视频| 九九爱精品视频在线观看| 亚洲无线观看免费| 久久国产乱子免费精品| 日本爱情动作片www.在线观看| 亚洲精品国产av成人精品| av在线老鸭窝| 国内精品一区二区在线观看| 久久精品国产亚洲av香蕉五月| 午夜福利在线观看吧| 国产大屁股一区二区在线视频| 悠悠久久av| 日本一二三区视频观看| 亚洲三级黄色毛片| 人妻夜夜爽99麻豆av| 国产成人福利小说| 日韩 亚洲 欧美在线| 一卡2卡三卡四卡精品乱码亚洲| 伦理电影大哥的女人| 最新中文字幕久久久久| 禁无遮挡网站| 中文字幕久久专区| 日韩av不卡免费在线播放| 在线观看美女被高潮喷水网站| 国产精品嫩草影院av在线观看| 97超碰精品成人国产| 精品人妻视频免费看| 成人特级av手机在线观看| 国产精品一二三区在线看| 成人欧美大片| 99国产极品粉嫩在线观看| 日本爱情动作片www.在线观看| 精品人妻熟女av久视频| 在线播放无遮挡| 国产不卡一卡二| 男人的好看免费观看在线视频| 亚洲丝袜综合中文字幕| 精品一区二区三区人妻视频| 国产一区二区三区在线臀色熟女| 欧美3d第一页| av专区在线播放| 国产不卡一卡二| 国产成人影院久久av| 99热精品在线国产| 久久99精品国语久久久| 国产成人a∨麻豆精品| 午夜精品一区二区三区免费看| 自拍偷自拍亚洲精品老妇| 男女视频在线观看网站免费| 亚洲在线观看片| 最好的美女福利视频网| 欧美成人一区二区免费高清观看| 午夜精品一区二区三区免费看| 不卡一级毛片| 国产黄片视频在线免费观看| 丝袜喷水一区| 亚洲成人中文字幕在线播放| 亚洲无线在线观看| 久久久国产成人精品二区| 人人妻人人澡人人爽人人夜夜 | 欧美成人精品欧美一级黄| 日韩人妻高清精品专区| 在线免费十八禁| 日韩av不卡免费在线播放| 51国产日韩欧美| 蜜臀久久99精品久久宅男| 美女高潮的动态| 欧美日韩在线观看h| 国产单亲对白刺激| 91av网一区二区| 嘟嘟电影网在线观看| eeuss影院久久| 国产久久久一区二区三区| 男女做爰动态图高潮gif福利片| 九九爱精品视频在线观看| 日韩制服骚丝袜av| 国产又黄又爽又无遮挡在线| 亚洲av免费在线观看| 青春草视频在线免费观看| 色哟哟·www|