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

    探索DNA條形碼對哈搖蚊類群部分屬物種鑒定的有效性

    2020-10-14 06:10:46余海軍黃鑫張鈺瑩周晶晶田召志陳佳林王茜
    天津農學院學報 2020年3期
    關鍵詞:搖蚊種間類群

    余海軍,黃鑫,張鈺瑩,周晶晶,田召志,陳佳林,王茜

    探索DNA條形碼對哈搖蚊類群部分屬物種鑒定的有效性

    余海軍,黃鑫,張鈺瑩,周晶晶,田召志,陳佳林,王茜通信作者

    (天津農學院 水產(chǎn)學院 天津市水生生態(tài)及養(yǎng)殖重點實驗室,天津 300384)

    采集中國部分地區(qū)哈搖蚊類群,并對其進行了DNA條形碼的鑒別研究。研究結果顯示:種間平均遺傳距離大約是種內平均遺傳距離的2.8倍,種內與種間存在明顯的“DNA barcoding gap”?;贏utomatic Barcode Gap Discovery(ABGD)法分析遺傳距離與其分布頻數(shù)的關系;鄰接樹(Neighbor joining tree,NJ tree)分析結果指出,87.5%的物種能區(qū)分開,并與形態(tài)學鑒定一致。本研究為探索DNA條形碼在哈搖蚊類群中鑒別的有效性提供了基礎數(shù)據(jù),同時為填補中國搖蚊昆蟲的DNA條形碼數(shù)據(jù)提供豐富的材料。

    哈搖蚊類群;DNA條形碼;物種鑒定;遺傳距離

    搖蚊因大部分幼蟲體色為紅色而常被漁民們稱為“紅蟲”,其隸屬于雙翅目(Diptera)長角亞目(Nematocera),其同時也是水生態(tài)系統(tǒng)中豐富度極高的無脊椎動物類群[1],能生存于各種類型的水體[2]。由于不同的搖蚊類群對不同水域環(huán)境有一定的適應性和敏感性,因此可以不同搖蚊物種的分布為指標,監(jiān)測水體水質變化[3]。

    哈搖蚊群歸屬于搖蚊亞科,其類群包括哈搖蚊屬()、枝角搖蚊屬()、擬枝角搖蚊屬()、小搖蚊屬()、隱搖蚊屬()、擬隱搖蚊屬()、彎鋏搖蚊屬()和羅搖蚊屬()8個近緣屬[4]。該類群的形態(tài)鑒別特征以生殖節(jié)上的下附器退化為特征,其中、是上下附器都比較退化,其余6屬的上附器為細長狀、足靴狀、指狀和葉狀等形狀而未退化。盡管屬間的分類特征明顯,但是各屬之間的物種存在極大的相似性,從而導致形態(tài)學鑒定難度加大。

    2003年加拿大動物學家Hebert提出了DNA條形碼,即通過使用一段標準的簡短基因序列作為條形碼對物種進行鑒別[5]。近年來DNA條形碼已成功用于搖蚊昆蟲的物種鑒定,特別是在物種的區(qū)分上,Lin等人基于DNA條形碼對長跗搖蚊屬()121個形態(tài)種進行物種分析,結果表明,DNA條形碼有效區(qū)分了早期形態(tài)學研究分類中94.6%的長跗搖蚊屬種類[6]。Song等通過DNA條形碼對多足搖蚊屬()的93個形態(tài)種進行物種鑒別,結果顯示94.4%的多足搖蚊屬物種能通過DNA條形碼區(qū)分[7]。此外,在新種的發(fā)現(xiàn)上,DNA條形碼也起到重要的驗證作用,Song等人基于DNA條形碼對一個中國新紀錄屬進行了物種鑒別,結果顯示該屬單獨聚為一支,與形態(tài)學研究結果相符[8]。Qi等人通過DNA條形碼對一個在魚塘表面滑行的海洋物種進行了物種鑒別,結果顯示其與二叉搖蚊屬()聚為一支,驗證了其屬于搖蚊的形態(tài)學鑒定,并將該物種命名為Qi& Lin sp. n,證明了搖蚊不只是生存在淡水環(huán)境中,也可以在海洋環(huán)境中生存[9]。Lin等人通過DNA條形碼對sp. n.進行分子分析,結果顯示其能與其他多足搖蚊屬物種區(qū)分開,與形態(tài)學鑒定一致[10]。

    本研究通過對中國部分地區(qū)的哈搖蚊類群物種中8屬20種(總計40樣本),選用COI基因作為分子標記,進行了DNA條形碼研究,研究結果為國內搖蚊昆蟲的DNA條形碼數(shù)據(jù)庫提供基礎數(shù)據(jù),并對國內外的同類研究做一個補充。

    1 材料與方法

    1.1 樣品來源

    本研究所選對象為哈搖蚊類群物種中的8屬20種,均來自中國部分地區(qū)野外采集,采集方式為掃網(wǎng)和燈誘,標本均保存于75%~85%乙醇中。所有采集的樣本均在解剖鏡下解剖封片,并將頭部和胸部保存于85%乙醇中,用于后續(xù)分子研究,所有證據(jù)標本均保存于天津農學院水產(chǎn)遺傳育種實驗室。研究中的對比序列為30條,均從BOLD(Barcode Of Life Data)數(shù)據(jù)庫上下載。

    1.2 DNA提取、PCR擴增、序列對比和分析

    1.2.1 DNA提取

    本試驗采用的是試劑盒(QIAGEN DNeasy Blood and Tissue kit)提取法,將整個胸部和頭部經(jīng)蛋白酶處理后,按照試劑盒說明書提取樣本DNA。

    1.2.2 PCR擴增

    聚合酶鏈式反應(PCR)擴增使用COI基因通用引物LCO1490和HCO2198[11]。反應體系為25 μL,包含12.5 μL 2× EsMasterMix(CoWin Biotech Co,Beijing,China),上下游引物各0.625 μL,dd H2O 8.75 μL。PCR通過MasterCycler Gradient進行擴增,擴增程序如下,預變性98 ℃ 10s,94 ℃ 5 min,35個循環(huán):94 ℃ 1 min,51 ℃ 1 min,72 ℃ 1 min。最后延伸72 ℃ 5 min,10 ℃保存。PCR擴增產(chǎn)物經(jīng)1.5%瓊脂糖凝膠電泳檢測后送北京六合華大基因科技股份有限公司進行測序。

    1.3 數(shù)據(jù)對比和分析方法

    1.3.1 序列對比與序列特征分析

    將測序公司返回的測序序列導入SeqMan中進行序列編輯,將編輯好的序列導出為FASTA文件保存,然后導入MEGA7.0進行比對并轉換成氨基酸檢查序列檢查是否存在終止密碼子[12],并對序列特征及各堿基含量進行統(tǒng)計。

    1.3.2 基于遺傳距離分析

    基于DNA條形碼遺傳距離的分析方法有ABGD法。利用ABGD[13]軟件在線分析,將編輯好的FASTA文件在線提交到ABGD網(wǎng)站(http:// wwwabi.snv.jussieu.fr/public/abgd/),選擇Kimura- 2-parameter(K2P)模型,種內差異先驗值(Prior intraspecific divergence)為0.001到0.100,相對gap寬度值(relative gap width)為1.0,其余參數(shù)設置為默認值[6]。

    1.3.3 基于系統(tǒng)發(fā)育樹分析

    在建樹過程中,將比對好的FASTA文件導入MEGA 7.0中并選用K2P分子進化模型,節(jié)點處置信度為1 000,其他設置均為默認值。

    單倍型結構分析通過將比對好的NEXUS文件導入PopART軟件選用TCS模型分析[14-15]。

    2 結果與分析

    2.1 序列特征結果

    本試驗一共獲得40條COI序列。所有序列的長度均為658 bp,各序列特征所示,其中各堿基含量的平均值分別為A=27.5%,T=38.0%,C=18.1%,G=16.5%,序列存在堿基的偏倚性,C+G的含量顯著低于A+T的含量,特別是第三核酸位點A+T的含量高達88.4%。通過比對序列得出,C(保守位點)占63.5%,V(變異位點)占36.2%。大部分的變異位點都集中在第三核酸位點,其比例為84.4%。結果表明,在第三核酸位點上的堿基進化速率快。

    2.2 遺傳距離結果

    本試驗40個樣本,形態(tài)鑒定為種,種內遺傳距離為0~0.154 4,平均值為0.048 4;種間遺傳距離為0.068 8~0.173 3,平均值為0.137 8。

    ABGD能使用直方圖來展現(xiàn)種內、種間遺傳距離及其分布頻數(shù)的關系,比較種內和種間的遺傳距離,當種間最小遺傳距離大于種內最大遺傳距離時,則出現(xiàn)“空白區(qū)”,如圖1所示。

    圖1 搖蚊亞科物種遺傳距離頻率直方圖

    2.3 系統(tǒng)發(fā)育樹結果

    基于40條DNA條形碼序列,構建NJ樹,如圖2所示。

    圖2 哈搖蚊類群NJ樹

    注:不同顏色或形狀代表不同物種,下同

    3 討論

    3.1 基于遺傳距離閾值分析

    DNA條形碼研究中,種內與種間遺傳距離的差異是一項非常重要的指標,當種間遺傳距離大于種內遺傳距離時,就會形成一個明顯的空白區(qū),即“DNA barcoding gap”。本研究中種內平均遺傳距離為0.048 4,種間平均遺傳距離為0.137 8,后者大約是前者的2.8倍,形成了“DNA barcoding gap”。

    遺傳距離閾值是指物種在能形成“DNA barcoding gap”的前提下,使用物種種間遺傳距離與參比的遺傳距離閾值進行對比,當種間遺傳距離大于參比的閾值時,則可以將物種區(qū)分為不同物種。閾值法在DNA條形碼的物種鑒定中起著不可替代的作用,然而隨著DNA條形碼研究不斷深入,有學者發(fā)現(xiàn)不同物種DNA條形碼的閾值各不相同,沒有一個統(tǒng)一的標準閾值[16]。BOLD數(shù)據(jù)庫最初將遺傳距離1%作為生物物種的劃分閾值,由于其閾值寬泛,在一定程度上造成了物種的“鑒別過度”,導致同種異名的出現(xiàn),為了不使物種混亂,近年來BOLD系統(tǒng)默認鑒別閾值為3%[17]。

    將本研究的序列結果導入Speciesdentifier軟件中,設定閾值為1%~10%,物種條形碼聚類結果如圖3所示,形態(tài)鑒定的數(shù)量為20種,當閾值在7%~8%時,mOTU數(shù)量最接近形態(tài)鑒定的數(shù)量,我們推測哈搖蚊類群的遺傳距離閾值為8%~9%。關于搖蚊亞科中其他類群的閾值,Song等區(qū)分多足搖蚊屬的遺傳距離閾值為7%~9%[7],Lin等人區(qū)分長跗搖蚊屬的遺傳距離閾值為4%~5%[6]。綜上所述,本研究分析的哈搖蚊類群遺傳距離閾值,與其同亞科的類群遺傳距離閾值相比,沒有較大差異,在一定程度上證明了DNA條形碼能有效地對哈搖蚊群部分屬物種進行鑒定。

    圖3 基于閾值的物種聚類分析

    3.2 基于系統(tǒng)發(fā)育樹分析

    系統(tǒng)發(fā)育樹能直觀地反映物種之間的親緣關系,本研究中采用的是NJ法,其基于遺傳距離和最小進化原理,通過自舉檢驗法進行聚類,并通過統(tǒng)計給定樹分支的可信度,得出最優(yōu)樹[18]。在NJ樹中,87.5%的物種能區(qū)分開,與形態(tài)學鑒定結果一致。本研究顯示在哈搖蚊屬中(圖4),雙葉哈搖蚊()與參考序列分成了3支,但因不能查詢參考序列的證據(jù)標本,所以無法對該數(shù)據(jù)做出明確的定論;此外,在哈搖蚊屬中,TN01短葉哈搖蚊()與參考序列GBMIN55160-17長距哈搖蚊()聚為一支,且置信度為100%,據(jù)閆春財博士論文中描述,與在外觀形態(tài)上很相似,即使通過生殖節(jié)的比對也不易區(qū)分,但腋瓣上是否具有緣毛可區(qū)分二者[19],TN01經(jīng)過形態(tài)鑒定此物種應為(圖5),而非?;谝陨贤茢?,對研究中哈搖蚊屬所有個體的DNA條形碼數(shù)據(jù)在PopART軟件中使用TCS模型進行單倍型結構分析,結果如圖6所示:除與聚為一體外,其他個體都能單獨聚為一體,該結果表明TN01與GBMIN55160-17出現(xiàn)同物異名。由于數(shù)據(jù)庫中GBMIN55160-17參考序列不能提供證據(jù)標本與TN01比較,本研究無法做出準確的定論,在今后的研究中將進一步驗證。系統(tǒng)發(fā)育樹的結果在一定程度上表明,DNA條形碼能有效鑒別搖蚊亞科物種。研究者在進行搖蚊亞科昆蟲遺傳距離界定時,需要比對正確而無疑的DNA條形碼數(shù)據(jù),在今后的工作中,項目組將與條形碼的同行做更深入地交流,為中國搖蚊亞科DNA條形碼數(shù)據(jù)庫的完善提供參考。

    圖4 哈搖蚊屬物種NJ樹

    圖5 TN01樣本的證據(jù)標本(翅部分)

    圖6 哈搖蚊屬基于TCS haplotype network的分析

    4 結論

    本研究通過采集中國部分地區(qū)的哈搖蚊類群昆蟲8屬20種(共計40樣本),進行了DNA條形碼數(shù)據(jù)分析,運用ABDG法分析該數(shù)據(jù)集種內、種間遺傳距離及其分布頻數(shù)的關系,基于NJ樹結果得出數(shù)據(jù)集分離出20個物種,與傳統(tǒng)形態(tài)學鑒定一致,證明DNA條形碼能準確鑒別哈搖蚊類群昆蟲,同時為填補中國搖蚊昆蟲的DNA條形碼數(shù)據(jù)提供豐富的材料。

    [1] Milo?evi?D,Simi?V,Stojkovi?M,et al. Chironomid faunal composition represented by taxonomic distinctness index reveals environmental change in a lotic system over three decades[J]. Hydrobiologia,2012,683(1):69-82.

    [2] Nicacio G,Juen L. Chironomids as indicators in freshwater ecosystems:An assessment of the literature[J]. Insect Conservation and Diversity,2015,8(5):393-403.

    [3] S?ther O A. Zoogeographical patterns in Chironomidae(Diptera)[J]. SIL Proceedings,2000,27(1):290-302.

    [4] 王新華,鄭樂怡,紀炳純. 中國搖蚊亞科記述Ⅲ. 哈搖蚊屬(雙翅目:搖蚊科)[J]. 動物分類學報,1993,18(4):459-465.

    [5] Hebert P D N,Cywinska A,Ball S L,et al. Biological identifications through DNA barcodes[J]. Proc R Soc Lond B Biol Sci,2003,270(1512):313-321.

    [6] Lin X L,Stur E,Ekrem T. Exploring genetic divergence in a species-rich insect genus using 2790 DNA barcodes[J]. Plos One,2015,10(9):1-24.

    [7] Song C,Lin X L,Wang Q,et al. DNA barcodes successfully delimit morphospecies in a superdiverse insect genus[J]. Zoological Scripta,2018,47(3):311-324.

    [8] Song C,Lin W B,Zhang R L,et al. Two gynandromorphic species ofKieffer,1912(Diptera:Chironomidae),with DNA barcodes from Oriental China[J]. The Pan-pacific Entomologist,2017,93(2):1-10.

    [9] Qi X,Lin X L,Ekrem T,et al. A new surface gliding species of Chironomidae:An independent invasion of marine environments and its evolutionary implications[J]. Zoological Scripta,2019,48(1):81-92.

    [10] Lin X L,Yu H J,Zhang R L,et al.()sp.n.(Diptera:Chironomidae)from Gaoligong Mountains,Yunnan,China[J]. Zootaxa,2019,4571(2):255-262.

    [11] Folmer O,Black M,Hoeh W,et al. DNA primers for amplification of mitochondrial cytochrome coxidase subunit 1 from diverse metazoan invertebrates[J]. Molecular Marine Biology and Biotechnology,1994,3(5):294-299.

    [12] Kumar S,Stecher G,Tamura K. MEGA7:Molecular evolutionary genetics analysis version 7.0 for bigger datasets[J]. Molecular Biology and Evolution,2016,33(7):1870-1874.

    [13] Puillandre N,Lambert A,Brouillet S,et al. ABGD,automatic barcode gap discovery for primary species delimitation[J]. Molecular Ecology,2012,21(8):1864-1877.

    [14] Clement M,Posada D,Crandall K A. TCS:A computer program to estimate gene genealogies[J]. Molecular Ecology,2000,9(10):1657-1659.

    [15] Leigh J W,Bryant D. POPART:Full-feature software for haplotype network construction[J]. Methods in Ecology and Evolution,2015,6(9):1110-1116.

    [16] Collins R A,Cruickshank R H. The seven deadly sins of DNA barcoding[J]. Molecular Ecology Resources,2013,13(6):969-975.

    [17] Ratnasingham S,Hebert P D N. BOLD:The barcode of life data system[J]. Molecular Ecology Notes,2007,7(3):355-364.

    [18] Duarte S,Batista D,Barlocher F,et al. Some new DNA barcodes of aquatic hyphomycete species[J]. Mycoscience,2015,56(1):102-108.

    [19] 閆春財. 中印區(qū)哈搖蚊屬復合體系統(tǒng)學研究(雙翅目:搖蚊科)[D]. 天津:南開大學,2005.

    Exploring the utility of DNA barcoding in species identification ofcomplex

    YU Hai-jun, HUANG Xin, ZHANG Yu-ying, ZHOU Jing-jing, TIAN Zhao-zhi, CHEN Jia-lin, WANG QianCorresponding Author

    (Key Laboratory of Aquatic-Ecology and Aquaculture of Tianjin, College of Fisheries, Tianjin Agricultural University, Tianjin 300384, China)

    This study collected Harischia complex species in regions of China, using DNA barcoding for species identification. The results showed: the mean interspecific genetic distance was 2.8 times larger than the mean intraspecific genetic distance, with existing distinct “DNA barcoding gap”. Based on Automatic Barcode Gap Discovery genetic distance was analyzed, and it was related with distributed frequency. The results of Neighbor joining tree analysis pointed that 87.5%of the species were distinguished, and it kept conformed to morphological identification. Our results produced basal data for exploring the utility of DNA barcoding in species identification ofcomplex, and provided abundant material for filling up the Chinese DNA barcodes database of Chironomidae.

    complex; DNA barcoding; species identification; genetic distance

    1008-5394(2020)03-0057-05

    10.19640/j.cnki.jtau.2020.03.013

    S931.1

    A

    2019-06-26

    國家自然科學基金(31672264,31672324)

    余海軍(1994-),男,碩士在讀,主要從事分子系統(tǒng)學研究。E-mail:yhj18322040463@163.com。

    王茜(1971-),女,副教授,博士,主要從事水生動物及昆蟲分子系統(tǒng)學的研究。E-mail:wqgt1999@163.com。

    責任編輯:張愛婷

    猜你喜歡
    搖蚊種間類群
    干旱條件對鬼針草和醉魚草種間相互作用及生長的影響
    植物研究(2023年5期)2023-09-09 08:01:22
    三峽庫區(qū)支流花溪河浮游植物種間關聯(lián)及影響因子分析
    云南省雙江縣4種多足搖蚊翅膀形態(tài)特征及COI基因序列分析
    薏苡種質資源ISSR分子標記篩選及親緣關系分析
    沉睡搖蚊的生存智慧
    黑石頂自然保護區(qū)土壤動物功能類群對季節(jié)動態(tài)的響應
    常用基因序列在搖蚊科昆蟲系統(tǒng)發(fā)育研究中的應用進展
    江蘇省宜興市茶園秋季雜草種間生態(tài)關系及群落分類
    雜草學報(2015年2期)2016-01-04 14:57:55
    2011年春夏季黃海和東海微型浮游動物類群組成及其攝食的研究
    營養(yǎng)鹽輸入方式對3種微藻生長及種間競爭的影響
    精品福利永久在线观看| 精品国产一区二区三区四区第35| 在线观看舔阴道视频| 人妻久久中文字幕网| 三级毛片av免费| 亚洲aⅴ乱码一区二区在线播放 | 国产精品久久电影中文字幕 | 男女高潮啪啪啪动态图| 亚洲一区高清亚洲精品| 国产精华一区二区三区| 中文字幕人妻丝袜制服| 日本精品一区二区三区蜜桃| 久久精品国产a三级三级三级| 建设人人有责人人尽责人人享有的| av电影中文网址| 亚洲一区二区三区不卡视频| 午夜日韩欧美国产| 欧美最黄视频在线播放免费 | 国产精品国产av在线观看| 色在线成人网| 精品国产亚洲在线| 男男h啪啪无遮挡| xxxhd国产人妻xxx| 成人18禁在线播放| 中文字幕色久视频| 午夜老司机福利片| 超色免费av| 中国美女看黄片| 精品免费久久久久久久清纯 | 91麻豆av在线| 亚洲欧美激情综合另类| 日韩有码中文字幕| av超薄肉色丝袜交足视频| 黑丝袜美女国产一区| 深夜精品福利| 99精品欧美一区二区三区四区| 日韩成人在线观看一区二区三区| 婷婷成人精品国产| 国产精品国产高清国产av | 欧美黑人精品巨大| 免费看a级黄色片| 亚洲精品成人av观看孕妇| 一二三四在线观看免费中文在| 欧美乱色亚洲激情| 咕卡用的链子| 国产精品电影一区二区三区 | 国产亚洲精品第一综合不卡| 日本黄色视频三级网站网址 | 天堂动漫精品| 久久天堂一区二区三区四区| 新久久久久国产一级毛片| 中文字幕制服av| 女同久久另类99精品国产91| 久久精品亚洲精品国产色婷小说| 一级毛片高清免费大全| 精品一区二区三区av网在线观看| 精品电影一区二区在线| 亚洲国产精品合色在线| 岛国在线观看网站| 韩国av一区二区三区四区| 中文字幕av电影在线播放| www.熟女人妻精品国产| 少妇的丰满在线观看| 国产欧美日韩综合在线一区二区| 如日韩欧美国产精品一区二区三区| 亚洲av熟女| 精品国产乱码久久久久久男人| 老司机影院毛片| 国产欧美亚洲国产| 一本一本久久a久久精品综合妖精| 国产不卡一卡二| 精品久久久久久,| 嫩草影视91久久| 欧美成人午夜精品| 99国产精品一区二区三区| 亚洲中文日韩欧美视频| 黄片播放在线免费| 欧美国产精品一级二级三级| 大片电影免费在线观看免费| 亚洲avbb在线观看| 成熟少妇高潮喷水视频| 很黄的视频免费| 日本一区二区免费在线视频| 久99久视频精品免费| 欧美乱妇无乱码| 91国产中文字幕| 成人手机av| 久久天躁狠狠躁夜夜2o2o| 久久久久久久久久久久大奶| 大型av网站在线播放| 欧美激情高清一区二区三区| 18禁美女被吸乳视频| 欧美激情极品国产一区二区三区| 深夜精品福利| 国产欧美日韩精品亚洲av| 免费在线观看完整版高清| 欧美午夜高清在线| 在线永久观看黄色视频| 久久精品aⅴ一区二区三区四区| 一a级毛片在线观看| 亚洲欧美色中文字幕在线| 天天躁狠狠躁夜夜躁狠狠躁| 国产一区在线观看成人免费| 亚洲国产看品久久| 日韩 欧美 亚洲 中文字幕| 国产精品久久久人人做人人爽| 国产精品免费一区二区三区在线 | 人人妻人人澡人人爽人人夜夜| 一个人免费在线观看的高清视频| 90打野战视频偷拍视频| 久久久久久久午夜电影 | a级片在线免费高清观看视频| 中文字幕人妻丝袜一区二区| 国产精品二区激情视频| 在线天堂中文资源库| 国产有黄有色有爽视频| 51午夜福利影视在线观看| 中文字幕人妻丝袜制服| 亚洲性夜色夜夜综合| 大香蕉久久网| 黄色视频,在线免费观看| 五月开心婷婷网| 国产三级黄色录像| tube8黄色片| 女同久久另类99精品国产91| 日日爽夜夜爽网站| 国产激情欧美一区二区| 日日爽夜夜爽网站| 两性午夜刺激爽爽歪歪视频在线观看 | 欧洲精品卡2卡3卡4卡5卡区| 国产精品免费大片| 成人18禁在线播放| av天堂久久9| 免费在线观看日本一区| 两个人看的免费小视频| 亚洲中文日韩欧美视频| 国产成+人综合+亚洲专区| 一本大道久久a久久精品| 热99国产精品久久久久久7| 亚洲一区中文字幕在线| 亚洲午夜精品一区,二区,三区| 欧美精品啪啪一区二区三区| 超碰成人久久| 咕卡用的链子| 777久久人妻少妇嫩草av网站| 中文亚洲av片在线观看爽 | 久久99一区二区三区| 中国美女看黄片| 成人手机av| 水蜜桃什么品种好| av不卡在线播放| 美女国产高潮福利片在线看| 午夜激情av网站| 黄网站色视频无遮挡免费观看| 欧美+亚洲+日韩+国产| 国产在线一区二区三区精| 久久精品aⅴ一区二区三区四区| 亚洲熟妇熟女久久| 欧美日韩乱码在线| 久久精品91无色码中文字幕| 18禁裸乳无遮挡免费网站照片 | 日韩 欧美 亚洲 中文字幕| 日本欧美视频一区| 在线永久观看黄色视频| 1024香蕉在线观看| 成人免费观看视频高清| 国产av精品麻豆| 国产精品 国内视频| 国产欧美日韩一区二区精品| 大码成人一级视频| 黑人欧美特级aaaaaa片| av国产精品久久久久影院| 亚洲一区中文字幕在线| 男人舔女人的私密视频| 国产男靠女视频免费网站| 日本欧美视频一区| 国产不卡av网站在线观看| 国产一区二区三区综合在线观看| 91成年电影在线观看| 欧美精品人与动牲交sv欧美| 嫁个100分男人电影在线观看| 国产欧美日韩一区二区精品| 日本撒尿小便嘘嘘汇集6| 美女高潮到喷水免费观看| 极品教师在线免费播放| 99国产精品一区二区蜜桃av | 高潮久久久久久久久久久不卡| 搡老熟女国产l中国老女人| 亚洲国产精品合色在线| 老汉色∧v一级毛片| 成人国产一区最新在线观看| 精品一品国产午夜福利视频| 一本一本久久a久久精品综合妖精| 黑丝袜美女国产一区| av网站免费在线观看视频| 亚洲熟女精品中文字幕| 最新在线观看一区二区三区| 多毛熟女@视频| 在线国产一区二区在线| 又黄又爽又免费观看的视频| 女同久久另类99精品国产91| av中文乱码字幕在线| 国产一卡二卡三卡精品| 在线天堂中文资源库| 亚洲av日韩在线播放| 久久性视频一级片| 精品一区二区三区四区五区乱码| 成熟少妇高潮喷水视频| 亚洲 欧美一区二区三区| 又黄又爽又免费观看的视频| 美女视频免费永久观看网站| 国产野战对白在线观看| 黄色女人牲交| 久久人妻福利社区极品人妻图片| 视频在线观看一区二区三区| 精品熟女少妇八av免费久了| 一级片'在线观看视频| 国产高清videossex| 韩国av一区二区三区四区| 欧美成人午夜精品| 久久国产乱子伦精品免费另类| 久久精品国产亚洲av香蕉五月 | 精品午夜福利视频在线观看一区| 伦理电影免费视频| 亚洲av片天天在线观看| 成人特级黄色片久久久久久久| 欧美精品高潮呻吟av久久| 高清视频免费观看一区二区| 久久精品国产亚洲av高清一级| 亚洲免费av在线视频| 国产在线精品亚洲第一网站| 久久精品亚洲熟妇少妇任你| 黄色 视频免费看| 亚洲国产看品久久| 露出奶头的视频| 十八禁人妻一区二区| 国产单亲对白刺激| 国产成人av激情在线播放| 757午夜福利合集在线观看| 日韩欧美国产一区二区入口| 人妻一区二区av| 亚洲欧美一区二区三区黑人| 精品无人区乱码1区二区| 国产精品 欧美亚洲| 老司机深夜福利视频在线观看| 老司机靠b影院| 久久亚洲精品不卡| 大陆偷拍与自拍| 日韩有码中文字幕| 搡老乐熟女国产| 91在线观看av| 精品电影一区二区在线| 91精品三级在线观看| 丝袜人妻中文字幕| www.精华液| 国产精品久久久久成人av| 国产视频一区二区在线看| 日本欧美视频一区| 国产一区二区三区在线臀色熟女 | 久久中文字幕人妻熟女| 嫁个100分男人电影在线观看| 女人被狂操c到高潮| 韩国av一区二区三区四区| 久久人人97超碰香蕉20202| 日韩欧美三级三区| cao死你这个sao货| 可以免费在线观看a视频的电影网站| 成人精品一区二区免费| 日本黄色视频三级网站网址 | 欧美日韩一级在线毛片| 国产片内射在线| 国产亚洲精品一区二区www | 久久久久国内视频| 美女扒开内裤让男人捅视频| 欧美性长视频在线观看| 极品少妇高潮喷水抽搐| 啦啦啦 在线观看视频| 国产精品成人在线| 国产精品永久免费网站| 99国产精品一区二区三区| 日本五十路高清| 十八禁网站免费在线| 精品福利永久在线观看| 成人三级做爰电影| 99国产精品一区二区三区| 国产高清视频在线播放一区| 1024香蕉在线观看| 一级,二级,三级黄色视频| 久久久国产一区二区| 亚洲专区字幕在线| 国产精品综合久久久久久久免费 | 精品国产一区二区三区久久久樱花| 免费观看a级毛片全部| 色尼玛亚洲综合影院| 岛国在线观看网站| 搡老岳熟女国产| 俄罗斯特黄特色一大片| 免费在线观看黄色视频的| 久久青草综合色| 黄色 视频免费看| 十八禁网站免费在线| 如日韩欧美国产精品一区二区三区| svipshipincom国产片| 日韩免费高清中文字幕av| 99国产极品粉嫩在线观看| 90打野战视频偷拍视频| 一级作爱视频免费观看| 国产99久久九九免费精品| 欧美亚洲日本最大视频资源| 亚洲精品美女久久久久99蜜臀| 麻豆国产av国片精品| 精品国产乱子伦一区二区三区| 国产一区二区三区视频了| www.999成人在线观看| 欧美成狂野欧美在线观看| 亚洲欧美日韩高清在线视频| 每晚都被弄得嗷嗷叫到高潮| 久久久水蜜桃国产精品网| 欧美精品一区二区免费开放| 久久热在线av| 精品国产亚洲在线| 日本精品一区二区三区蜜桃| 欧美精品一区二区免费开放| 欧美大码av| 黑人巨大精品欧美一区二区mp4| 两个人免费观看高清视频| 99国产精品一区二区蜜桃av | 国产高清视频在线播放一区| 人妻一区二区av| 欧美精品一区二区免费开放| 我的亚洲天堂| 一进一出抽搐gif免费好疼 | 日韩欧美三级三区| tocl精华| 国产精品电影一区二区三区 | 国产黄色免费在线视频| 伦理电影免费视频| 久久香蕉精品热| 另类亚洲欧美激情| av超薄肉色丝袜交足视频| 美女扒开内裤让男人捅视频| a级毛片黄视频| 成年动漫av网址| 51午夜福利影视在线观看| 最新在线观看一区二区三区| 色婷婷久久久亚洲欧美| 午夜视频精品福利| 精品国产一区二区三区四区第35| 国产欧美亚洲国产| 免费女性裸体啪啪无遮挡网站| 成人国语在线视频| 免费日韩欧美在线观看| 天天影视国产精品| 啪啪无遮挡十八禁网站| 亚洲国产中文字幕在线视频| 757午夜福利合集在线观看| 欧美av亚洲av综合av国产av| 男女免费视频国产| 成人18禁在线播放| 黄网站色视频无遮挡免费观看| 日韩三级视频一区二区三区| 日韩一卡2卡3卡4卡2021年| 两性午夜刺激爽爽歪歪视频在线观看 | 99国产精品免费福利视频| 一区在线观看完整版| 啦啦啦视频在线资源免费观看| 欧美+亚洲+日韩+国产| 亚洲av电影在线进入| 无人区码免费观看不卡| 搡老熟女国产l中国老女人| 亚洲国产毛片av蜜桃av| 国产精品 国内视频| 人人妻人人添人人爽欧美一区卜| 亚洲全国av大片| 欧美人与性动交α欧美精品济南到| 欧美中文综合在线视频| 一二三四在线观看免费中文在| 欧美午夜高清在线| 波多野结衣一区麻豆| 王馨瑶露胸无遮挡在线观看| 一级a爱片免费观看的视频| av超薄肉色丝袜交足视频| 久久久久久久久免费视频了| 一进一出抽搐gif免费好疼 | 亚洲片人在线观看| 国产精品成人在线| 国产成人免费无遮挡视频| 亚洲成人免费av在线播放| 国产精品自产拍在线观看55亚洲 | 亚洲av成人不卡在线观看播放网| 精品国产乱码久久久久久男人| 国产精品免费大片| 热re99久久精品国产66热6| 国产免费现黄频在线看| 国产av又大| svipshipincom国产片| 变态另类成人亚洲欧美熟女 | 亚洲黑人精品在线| svipshipincom国产片| 动漫黄色视频在线观看| 欧美激情极品国产一区二区三区| 99久久精品国产亚洲精品| 亚洲一区中文字幕在线| 精品人妻熟女毛片av久久网站| 夜夜夜夜夜久久久久| 9热在线视频观看99| 99热只有精品国产| 桃红色精品国产亚洲av| 精品国内亚洲2022精品成人 | 欧美激情 高清一区二区三区| 视频在线观看一区二区三区| 怎么达到女性高潮| 欧美+亚洲+日韩+国产| 18禁美女被吸乳视频| 国内久久婷婷六月综合欲色啪| 亚洲精品中文字幕在线视频| 午夜福利在线免费观看网站| 欧美成人免费av一区二区三区 | av线在线观看网站| 亚洲精品中文字幕一二三四区| 国产精品 欧美亚洲| 午夜精品国产一区二区电影| 精品免费久久久久久久清纯 | 极品少妇高潮喷水抽搐| 久久精品国产亚洲av高清一级| av一本久久久久| 美女视频免费永久观看网站| 午夜福利免费观看在线| 大码成人一级视频| 日韩中文字幕欧美一区二区| 岛国毛片在线播放| 午夜成年电影在线免费观看| www.自偷自拍.com| 国产精品久久视频播放| 在线国产一区二区在线| 少妇被粗大的猛进出69影院| 69精品国产乱码久久久| 女性被躁到高潮视频| bbb黄色大片| 丁香六月欧美| 国产xxxxx性猛交| 在线观看午夜福利视频| 91精品国产国语对白视频| 好男人电影高清在线观看| 国产精品久久久久久人妻精品电影| 国产aⅴ精品一区二区三区波| 大型黄色视频在线免费观看| 国产精品 国内视频| 看免费av毛片| 波多野结衣av一区二区av| 这个男人来自地球电影免费观看| 看免费av毛片| 亚洲第一av免费看| 男人的好看免费观看在线视频 | 久久天堂一区二区三区四区| 三级毛片av免费| 国产亚洲av高清不卡| 大陆偷拍与自拍| 欧美精品av麻豆av| 99国产极品粉嫩在线观看| 人成视频在线观看免费观看| 亚洲精品久久午夜乱码| 精品亚洲成国产av| 99riav亚洲国产免费| 搡老熟女国产l中国老女人| 精品一品国产午夜福利视频| 亚洲色图 男人天堂 中文字幕| 日韩免费高清中文字幕av| 一进一出好大好爽视频| 久久这里只有精品19| 侵犯人妻中文字幕一二三四区| 国产精品自产拍在线观看55亚洲 | 新久久久久国产一级毛片| 大香蕉久久成人网| 午夜91福利影院| 91麻豆精品激情在线观看国产 | bbb黄色大片| 欧美精品亚洲一区二区| 精品久久久久久久毛片微露脸| 多毛熟女@视频| 新久久久久国产一级毛片| 又大又爽又粗| 咕卡用的链子| 中文字幕人妻丝袜一区二区| 久久久国产成人精品二区 | 成年版毛片免费区| 午夜福利免费观看在线| 搡老乐熟女国产| 一区二区三区精品91| 黑人操中国人逼视频| 免费在线观看完整版高清| 色综合欧美亚洲国产小说| 日韩制服丝袜自拍偷拍| 久久九九热精品免费| www日本在线高清视频| 久久香蕉国产精品| 成人18禁高潮啪啪吃奶动态图| 成年动漫av网址| 女人精品久久久久毛片| 妹子高潮喷水视频| 亚洲黑人精品在线| 亚洲av成人一区二区三| 亚洲美女黄片视频| 亚洲黑人精品在线| 成人永久免费在线观看视频| 丁香欧美五月| 一级,二级,三级黄色视频| 一本综合久久免费| 99国产精品免费福利视频| 免费久久久久久久精品成人欧美视频| 久久人人97超碰香蕉20202| 欧美色视频一区免费| 老司机在亚洲福利影院| 99国产综合亚洲精品| 精品国产乱子伦一区二区三区| 欧美黑人精品巨大| 亚洲精品乱久久久久久| 成人国产一区最新在线观看| 久久精品亚洲av国产电影网| 精品福利观看| 男女之事视频高清在线观看| 看片在线看免费视频| a级片在线免费高清观看视频| 美女扒开内裤让男人捅视频| 精品久久久久久,| 国产成人av教育| 999久久久精品免费观看国产| 一级,二级,三级黄色视频| 99国产精品一区二区三区| 久久国产精品大桥未久av| 窝窝影院91人妻| 91麻豆精品激情在线观看国产 | 日本撒尿小便嘘嘘汇集6| 女人被躁到高潮嗷嗷叫费观| 中国美女看黄片| 色婷婷久久久亚洲欧美| 黄色视频,在线免费观看| 国产1区2区3区精品| 一本一本久久a久久精品综合妖精| 欧美国产精品一级二级三级| 精品久久久久久,| 精品国产一区二区久久| 飞空精品影院首页| 在线观看免费日韩欧美大片| 变态另类成人亚洲欧美熟女 | 日日摸夜夜添夜夜添小说| 99re6热这里在线精品视频| 制服人妻中文乱码| 国产主播在线观看一区二区| 国产成人欧美在线观看 | 欧美乱妇无乱码| 变态另类成人亚洲欧美熟女 | 欧美日韩国产mv在线观看视频| 国产精品九九99| 久久青草综合色| 欧美黑人欧美精品刺激| 高清在线国产一区| 精品一品国产午夜福利视频| 777米奇影视久久| 日韩成人在线观看一区二区三区| 日本撒尿小便嘘嘘汇集6| av中文乱码字幕在线| 老司机午夜十八禁免费视频| 日韩欧美国产一区二区入口| 女警被强在线播放| 国产亚洲一区二区精品| 久久久久国内视频| 一二三四在线观看免费中文在| 交换朋友夫妻互换小说| av不卡在线播放| 日韩欧美免费精品| 中文字幕高清在线视频| 女性生殖器流出的白浆| 日日摸夜夜添夜夜添小说| 99国产精品免费福利视频| 国产精品一区二区免费欧美| 欧美午夜高清在线| 伦理电影免费视频| 1024香蕉在线观看| 9色porny在线观看| 中文字幕人妻丝袜一区二区| 91老司机精品| 亚洲久久久国产精品| 人妻丰满熟妇av一区二区三区 | 狠狠狠狠99中文字幕| 欧美 亚洲 国产 日韩一| 免费女性裸体啪啪无遮挡网站| 亚洲专区字幕在线| 免费观看人在逋| 欧美乱色亚洲激情| 性色av乱码一区二区三区2| 热99国产精品久久久久久7| 狠狠狠狠99中文字幕| 亚洲一区中文字幕在线| 满18在线观看网站| 亚洲午夜精品一区,二区,三区| 夜夜爽天天搞| 97人妻天天添夜夜摸| 久久天躁狠狠躁夜夜2o2o| 麻豆成人av在线观看| 亚洲五月婷婷丁香| 老汉色∧v一级毛片| 中文字幕人妻熟女乱码| 精品电影一区二区在线| 免费不卡黄色视频| 91成人精品电影| 国产亚洲欧美精品永久| 午夜福利乱码中文字幕| 国产成人啪精品午夜网站| 亚洲视频免费观看视频| 国产成人av激情在线播放| 天堂动漫精品| 99国产精品99久久久久| 国产有黄有色有爽视频| 久久久久国产精品人妻aⅴ院 |