• <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種微藻生長及種間競爭的影響
    免费在线观看成人毛片| 毛片女人毛片| 久久精品国产亚洲av香蕉五月| 国产精品久久久久久亚洲av鲁大| 亚洲色图av天堂| 人人妻人人看人人澡| 精品久久国产蜜桃| 日韩中文字幕欧美一区二区| 国产精品不卡视频一区二区| 色av中文字幕| 久久草成人影院| 亚洲精华国产精华精| 村上凉子中文字幕在线| 我的老师免费观看完整版| 永久网站在线| 久久久久精品国产欧美久久久| 国产精华一区二区三区| 日韩欧美精品免费久久| 日日干狠狠操夜夜爽| 国产蜜桃级精品一区二区三区| 国产精品福利在线免费观看| 欧美一区二区精品小视频在线| 国内精品一区二区在线观看| 国产精品一区二区免费欧美| 国产精品永久免费网站| 国产69精品久久久久777片| 最新在线观看一区二区三区| 一本久久中文字幕| 97超视频在线观看视频| av国产免费在线观看| 久久热精品热| 91午夜精品亚洲一区二区三区 | 亚洲乱码一区二区免费版| 亚洲欧美日韩东京热| 亚洲在线观看片| x7x7x7水蜜桃| 午夜老司机福利剧场| 国产伦人伦偷精品视频| 97超视频在线观看视频| 丰满的人妻完整版| 精品人妻偷拍中文字幕| 国产精品一区二区性色av| 亚洲男人的天堂狠狠| 男人和女人高潮做爰伦理| 久久久久久久午夜电影| av国产免费在线观看| 精品一区二区三区人妻视频| 日韩,欧美,国产一区二区三区 | 成人二区视频| 国产精品伦人一区二区| 欧美黑人巨大hd| 国产一区二区在线av高清观看| 精品人妻视频免费看| 国产乱人视频| 在线观看av片永久免费下载| 国产精品一区二区性色av| 少妇高潮的动态图| 国产亚洲精品av在线| 91狼人影院| 看片在线看免费视频| 天堂网av新在线| 少妇人妻精品综合一区二区 | 亚洲国产精品成人综合色| 亚洲美女视频黄频| 国产蜜桃级精品一区二区三区| 麻豆成人av在线观看| 1000部很黄的大片| 国产一区二区三区视频了| 国产伦精品一区二区三区四那| 免费看日本二区| 很黄的视频免费| 很黄的视频免费| 久久九九热精品免费| 亚洲欧美日韩高清专用| 国产精品一区www在线观看 | 黄色配什么色好看| 欧美绝顶高潮抽搐喷水| 国产午夜精品论理片| 窝窝影院91人妻| 啦啦啦韩国在线观看视频| 欧美高清成人免费视频www| 91在线观看av| 日韩大尺度精品在线看网址| 日本一二三区视频观看| 亚洲av一区综合| 嫩草影院精品99| 国内毛片毛片毛片毛片毛片| 夜夜爽天天搞| 噜噜噜噜噜久久久久久91| 国产精品国产三级国产av玫瑰| 午夜福利欧美成人| 日本色播在线视频| 成人三级黄色视频| 国产伦精品一区二区三区视频9| 亚洲人成网站在线播放欧美日韩| 国产伦一二天堂av在线观看| 黄色欧美视频在线观看| 91麻豆av在线| x7x7x7水蜜桃| videossex国产| 18禁在线播放成人免费| 久久久久久伊人网av| 亚洲av免费在线观看| 久久久久久久久久成人| 国产高清三级在线| 俺也久久电影网| 永久网站在线| 99精品在免费线老司机午夜| 亚洲成人久久爱视频| 日本欧美国产在线视频| 国产精华一区二区三区| 亚洲av熟女| 男人舔女人下体高潮全视频| 国产精品乱码一区二三区的特点| 97超视频在线观看视频| 午夜福利欧美成人| 3wmmmm亚洲av在线观看| 午夜免费男女啪啪视频观看 | 97热精品久久久久久| 亚洲国产高清在线一区二区三| 国产精品99久久久久久久久| 老司机福利观看| 欧美xxxx性猛交bbbb| 午夜老司机福利剧场| 观看美女的网站| 在线观看一区二区三区| 亚洲在线自拍视频| 欧美成人免费av一区二区三区| 高清毛片免费观看视频网站| 成年女人永久免费观看视频| 搡女人真爽免费视频火全软件 | 亚洲av熟女| 色吧在线观看| av中文乱码字幕在线| 日日夜夜操网爽| 成人亚洲精品av一区二区| 春色校园在线视频观看| 黄色视频,在线免费观看| 午夜a级毛片| 欧美日韩精品成人综合77777| 亚洲成人中文字幕在线播放| 又紧又爽又黄一区二区| 免费一级毛片在线播放高清视频| 久久草成人影院| 亚洲四区av| 69人妻影院| 欧美精品国产亚洲| 色综合亚洲欧美另类图片| 成人国产麻豆网| 欧美高清成人免费视频www| 级片在线观看| 成人亚洲精品av一区二区| 老司机福利观看| 美女高潮的动态| 亚洲va在线va天堂va国产| 观看美女的网站| 欧美激情在线99| 99九九线精品视频在线观看视频| 欧美性猛交黑人性爽| 亚洲欧美日韩无卡精品| 又爽又黄a免费视频| 91久久精品国产一区二区成人| 搡老妇女老女人老熟妇| 欧美性猛交╳xxx乱大交人| 欧美日韩中文字幕国产精品一区二区三区| 欧美日韩精品成人综合77777| 国产午夜精品论理片| 国产伦精品一区二区三区四那| 亚洲av不卡在线观看| 在线观看午夜福利视频| av视频在线观看入口| 国产精品亚洲美女久久久| 午夜福利成人在线免费观看| 亚洲欧美激情综合另类| 有码 亚洲区| 日韩欧美三级三区| 国产欧美日韩精品亚洲av| 精品久久久久久,| 国产黄片美女视频| 亚洲欧美激情综合另类| 国产精品人妻久久久影院| 99视频精品全部免费 在线| 国产精品一区二区免费欧美| 在线免费十八禁| 亚洲成人中文字幕在线播放| 嫩草影视91久久| 变态另类成人亚洲欧美熟女| 人人妻人人看人人澡| 人妻丰满熟妇av一区二区三区| 成熟少妇高潮喷水视频| 亚洲成a人片在线一区二区| 亚洲精品456在线播放app | 久久久国产成人免费| 69人妻影院| 国产精品爽爽va在线观看网站| 变态另类丝袜制服| 97人妻精品一区二区三区麻豆| 中文字幕av成人在线电影| 中文字幕av成人在线电影| 日本-黄色视频高清免费观看| 麻豆久久精品国产亚洲av| 女同久久另类99精品国产91| 久久久久久久久中文| 成人美女网站在线观看视频| 最近最新中文字幕大全电影3| 国产综合懂色| 少妇的逼水好多| 中国美白少妇内射xxxbb| 免费在线观看日本一区| 精品久久久久久久人妻蜜臀av| 中文字幕高清在线视频| 国产 一区 欧美 日韩| 内地一区二区视频在线| 狂野欧美白嫩少妇大欣赏| 国产av不卡久久| 久久久成人免费电影| 国产精品久久久久久精品电影| 深夜精品福利| 国产国拍精品亚洲av在线观看| 一区福利在线观看| 一区福利在线观看| 日本 av在线| 老熟妇乱子伦视频在线观看| 欧美激情国产日韩精品一区| 欧美人与善性xxx| 久久国产精品人妻蜜桃| 又爽又黄a免费视频| 精品人妻熟女av久视频| 午夜亚洲福利在线播放| 两个人视频免费观看高清| av视频在线观看入口| 69av精品久久久久久| 天美传媒精品一区二区| 国产伦精品一区二区三区视频9| 99热网站在线观看| 久久精品国产清高在天天线| avwww免费| 精品久久久久久久久亚洲 | 99热这里只有是精品50| 国产 一区 欧美 日韩| 美女被艹到高潮喷水动态| 内地一区二区视频在线| 欧美中文日本在线观看视频| 色视频www国产| 无人区码免费观看不卡| 国产成人av教育| 国产精品美女特级片免费视频播放器| 精品欧美国产一区二区三| 亚洲人与动物交配视频| 午夜精品久久久久久毛片777| 99在线视频只有这里精品首页| 最后的刺客免费高清国语| 亚洲精华国产精华精| 亚洲黑人精品在线| 69人妻影院| 99久久无色码亚洲精品果冻| 女人十人毛片免费观看3o分钟| 一个人免费在线观看电影| 91在线精品国自产拍蜜月| 亚洲av成人av| 亚洲无线在线观看| 久久人人爽人人爽人人片va| a级毛片免费高清观看在线播放| 哪里可以看免费的av片| 男人舔女人下体高潮全视频| 少妇人妻一区二区三区视频| 日韩欧美一区二区三区在线观看| 有码 亚洲区| 99国产精品一区二区蜜桃av| 欧美一区二区国产精品久久精品| 亚洲精品成人久久久久久| 国产蜜桃级精品一区二区三区| 国产私拍福利视频在线观看| 午夜福利在线观看吧| x7x7x7水蜜桃| 精品人妻熟女av久视频| 久久久久九九精品影院| 深夜a级毛片| 成人三级黄色视频| 久久国产乱子免费精品| 美女大奶头视频| 在线观看免费视频日本深夜| .国产精品久久| 一级黄色大片毛片| 两性午夜刺激爽爽歪歪视频在线观看| 日本熟妇午夜| 动漫黄色视频在线观看| 在线观看av片永久免费下载| 成人国产麻豆网| 欧美日本亚洲视频在线播放| 国产精品日韩av在线免费观看| 亚洲黑人精品在线| 久久久色成人| 亚洲人成伊人成综合网2020| 久久国产精品人妻蜜桃| 黄色一级大片看看| 精品人妻偷拍中文字幕| 久久99热这里只有精品18| 一进一出好大好爽视频| 亚州av有码| 国产免费一级a男人的天堂| 国产精品久久电影中文字幕| 免费在线观看日本一区| 三级国产精品欧美在线观看| 久久精品综合一区二区三区| a在线观看视频网站| 91麻豆av在线| 日韩欧美在线乱码| 国内精品美女久久久久久| 色综合亚洲欧美另类图片| 国产精品自产拍在线观看55亚洲| 久久午夜亚洲精品久久| 亚洲图色成人| 真人做人爱边吃奶动态| 国产成人一区二区在线| 久久久久免费精品人妻一区二区| 欧美激情久久久久久爽电影| 麻豆成人av在线观看| 老熟妇仑乱视频hdxx| 美女免费视频网站| 久久久久久久久大av| 欧美性猛交╳xxx乱大交人| 日本a在线网址| 久久午夜福利片| 在线观看午夜福利视频| 色视频www国产| 啦啦啦啦在线视频资源| 亚洲欧美精品综合久久99| 久久精品国产鲁丝片午夜精品 | 91麻豆av在线| 日韩欧美在线乱码| 国产极品精品免费视频能看的| 神马国产精品三级电影在线观看| 少妇丰满av| 日韩中字成人| 国产男人的电影天堂91| 日本免费a在线| 欧美色视频一区免费| 蜜桃久久精品国产亚洲av| 可以在线观看毛片的网站| 色综合色国产| 中文字幕久久专区| 天堂av国产一区二区熟女人妻| 成人永久免费在线观看视频| 不卡一级毛片| 一区福利在线观看| 免费看美女性在线毛片视频| 国产精品三级大全| 国产精品一区二区三区四区久久| 久久久精品欧美日韩精品| 亚洲国产精品成人综合色| 一区二区三区免费毛片| 欧美色欧美亚洲另类二区| 成人一区二区视频在线观看| 中文亚洲av片在线观看爽| 久久久久国产精品人妻aⅴ院| 午夜福利高清视频| 国产真实伦视频高清在线观看 | 欧美不卡视频在线免费观看| 亚洲中文字幕一区二区三区有码在线看| 黄片wwwwww| 长腿黑丝高跟| 毛片一级片免费看久久久久 | 国产高清视频在线播放一区| 欧美+亚洲+日韩+国产| 哪里可以看免费的av片| 1024手机看黄色片| 欧美色欧美亚洲另类二区| 性欧美人与动物交配| 久久精品国产亚洲av天美| 大又大粗又爽又黄少妇毛片口| 亚洲,欧美,日韩| 国产伦在线观看视频一区| 欧美日韩综合久久久久久 | 久久精品国产亚洲av天美| 看黄色毛片网站| 亚洲精品一卡2卡三卡4卡5卡| 国产黄片美女视频| 国内精品美女久久久久久| 男女下面进入的视频免费午夜| 国产国拍精品亚洲av在线观看| 色哟哟·www| 我要看日韩黄色一级片| 日本 欧美在线| 亚洲欧美精品综合久久99| 欧美xxxx性猛交bbbb| 日日夜夜操网爽| 国产精品国产高清国产av| 国产精品一区二区性色av| 国产三级在线视频| 女生性感内裤真人,穿戴方法视频| 国产美女午夜福利| 中文字幕av在线有码专区| 国产乱人视频| 国国产精品蜜臀av免费| 性色avwww在线观看| 亚洲美女搞黄在线观看 | 亚洲成人久久爱视频| 成人精品一区二区免费| 成年女人看的毛片在线观看| 九九久久精品国产亚洲av麻豆| 十八禁网站免费在线| 亚洲av成人av| 男人的好看免费观看在线视频| 精品99又大又爽又粗少妇毛片 | 国产白丝娇喘喷水9色精品| 日韩,欧美,国产一区二区三区 | 国产一区二区三区视频了| 国产av在哪里看| 日韩欧美三级三区| 99九九线精品视频在线观看视频| 国产精华一区二区三区| 色尼玛亚洲综合影院| 少妇高潮的动态图| 男人舔奶头视频| 少妇的逼好多水| 麻豆成人午夜福利视频| 色尼玛亚洲综合影院| 国产精品99久久久久久久久| 亚洲真实伦在线观看| 亚洲熟妇熟女久久| 两个人视频免费观看高清| 婷婷丁香在线五月| 久久久精品大字幕| 国产精品国产高清国产av| 波野结衣二区三区在线| 一级a爱片免费观看的视频| 欧美激情在线99| 免费黄网站久久成人精品| 18+在线观看网站| 亚洲 国产 在线| 亚洲熟妇中文字幕五十中出| 亚洲精品影视一区二区三区av| 如何舔出高潮| 国产探花极品一区二区| 男女啪啪激烈高潮av片| 丝袜美腿在线中文| 窝窝影院91人妻| 天堂动漫精品| 亚洲最大成人手机在线| 变态另类成人亚洲欧美熟女| 国内毛片毛片毛片毛片毛片| 久久久久久久久久久丰满 | 日本五十路高清| 久久亚洲真实| 亚洲中文日韩欧美视频| 美女黄网站色视频| 色视频www国产| 国产精品亚洲一级av第二区| 国产国拍精品亚洲av在线观看| 夜夜夜夜夜久久久久| 亚洲av五月六月丁香网| 久久中文看片网| 很黄的视频免费| 国产精品1区2区在线观看.| 亚洲av成人av| 国产单亲对白刺激| 国产精品不卡视频一区二区| 亚洲一区高清亚洲精品| 中文字幕熟女人妻在线| 成人av在线播放网站| 亚洲中文日韩欧美视频| 久久欧美精品欧美久久欧美| 免费高清视频大片| 啪啪无遮挡十八禁网站| 国产高清有码在线观看视频| 欧美3d第一页| 一本久久中文字幕| 亚洲不卡免费看| 天天一区二区日本电影三级| 欧美三级亚洲精品| 国产精品免费一区二区三区在线| 日韩欧美精品免费久久| 国产在线精品亚洲第一网站| 国产精品精品国产色婷婷| 男插女下体视频免费在线播放| 别揉我奶头 嗯啊视频| 亚洲欧美日韩东京热| 男女做爰动态图高潮gif福利片| 国产成人福利小说| 午夜爱爱视频在线播放| 人妻丰满熟妇av一区二区三区| 在线观看av片永久免费下载| 国产高清不卡午夜福利| 免费av毛片视频| 99久久无色码亚洲精品果冻| 一本精品99久久精品77| 永久网站在线| 村上凉子中文字幕在线| 一个人观看的视频www高清免费观看| 亚洲国产精品久久男人天堂| 精品一区二区三区人妻视频| 亚洲性夜色夜夜综合| 97超级碰碰碰精品色视频在线观看| 一卡2卡三卡四卡精品乱码亚洲| 久久久久免费精品人妻一区二区| 日本黄色视频三级网站网址| 国产男人的电影天堂91| 欧美色欧美亚洲另类二区| 欧美不卡视频在线免费观看| 男女啪啪激烈高潮av片| 女生性感内裤真人,穿戴方法视频| 热99re8久久精品国产| 高清在线国产一区| 精品一区二区三区视频在线观看免费| 国国产精品蜜臀av免费| 亚洲一区高清亚洲精品| 免费看光身美女| 少妇丰满av| 亚洲国产精品成人综合色| 国产大屁股一区二区在线视频| 女同久久另类99精品国产91| 最好的美女福利视频网| 91av网一区二区| 国产精品一区二区三区四区久久| 嫩草影院精品99| www.www免费av| 亚洲欧美日韩高清在线视频| 动漫黄色视频在线观看| 久久人人爽人人爽人人片va| 国产成人aa在线观看| av在线老鸭窝| 欧美一区二区精品小视频在线| 女生性感内裤真人,穿戴方法视频| 99精品久久久久人妻精品| 欧美日本视频| 精品久久久久久久久亚洲 | 岛国在线免费视频观看| 国产精品一区二区免费欧美| 成人国产一区最新在线观看| 成人午夜高清在线视频| av在线老鸭窝| 成人毛片a级毛片在线播放| 露出奶头的视频| 如何舔出高潮| 伦精品一区二区三区| 亚洲专区中文字幕在线| 欧美激情在线99| 久久国产精品人妻蜜桃| 夜夜看夜夜爽夜夜摸| 国产成人福利小说| 狂野欧美白嫩少妇大欣赏| 日本黄色视频三级网站网址| 久久久久久久精品吃奶| 国产精品99久久久久久久久| 在线观看美女被高潮喷水网站| 欧美日韩国产亚洲二区| 99九九线精品视频在线观看视频| 久久久久久久久中文| 露出奶头的视频| 久久精品国产鲁丝片午夜精品 | 在线观看午夜福利视频| 亚洲va日本ⅴa欧美va伊人久久| 极品教师在线免费播放| 免费av观看视频| 淫秽高清视频在线观看| 成人精品一区二区免费| 久久久久久久午夜电影| 日日夜夜操网爽| 琪琪午夜伦伦电影理论片6080| 精品久久久久久久人妻蜜臀av| 国产精品日韩av在线免费观看| 淫妇啪啪啪对白视频| 精品午夜福利在线看| 久久精品国产鲁丝片午夜精品 | 国产真实乱freesex| 午夜福利高清视频| 别揉我奶头~嗯~啊~动态视频| 国产三级中文精品| 99热这里只有是精品在线观看| 精华霜和精华液先用哪个| 99久国产av精品| 午夜a级毛片| 国产av一区在线观看免费| 国产精品,欧美在线| 啦啦啦韩国在线观看视频| 国国产精品蜜臀av免费| 联通29元200g的流量卡| 久久久久久久久久黄片| 日本一本二区三区精品| 高清日韩中文字幕在线| 久久人妻av系列| 99热只有精品国产| 淫妇啪啪啪对白视频| 国产亚洲91精品色在线| 欧美最黄视频在线播放免费| 国产黄a三级三级三级人| 日韩一区二区视频免费看| 国产一区二区在线av高清观看| 99久久精品一区二区三区| 国产一区二区亚洲精品在线观看| 国产av麻豆久久久久久久| 色综合站精品国产| 免费av观看视频| 一级毛片久久久久久久久女| 欧美精品啪啪一区二区三区| 99国产精品一区二区蜜桃av| 国产精品1区2区在线观看.| 免费观看的影片在线观看| 成年女人看的毛片在线观看| 久久国内精品自在自线图片| 国产一级毛片七仙女欲春2| 在现免费观看毛片| 亚洲一区二区三区色噜噜| 99在线视频只有这里精品首页| 免费在线观看影片大全网站| 在线免费十八禁| 97碰自拍视频| 两人在一起打扑克的视频| 色在线成人网| 精品日产1卡2卡| 黄色女人牲交| 一本精品99久久精品77| 国产精品人妻久久久影院| 免费在线观看日本一区|