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

    中國北方常見芫菁分子物種界定(鞘翅目: 芫菁科)

    2022-10-19 06:58:04程海云段家充
    昆蟲學(xué)報 2022年9期
    關(guān)鍵詞:物種分析

    程海云, 段家充, 張 超, 潘 昭

    (河北大學(xué)生命科學(xué)學(xué)院, 生命科學(xué)與綠色發(fā)展研究院, 河北省動物系統(tǒng)學(xué)與應(yīng)用重點(diǎn)實(shí)驗(yàn)室, 河北保定 071002)

    芫菁科(Meloidae)隸屬于鞘翅目(Coleoptera)擬步甲總科(Tenebrionoidea),世界已知3亞科16族133屬近3 000種,廣泛分布于除新西蘭、南極和波利尼西亞群島之外的世界各地;中國記錄2亞科8族27屬近200種(亞種),各省市自治區(qū)均有記載(Bolognaetal., 2010; Pan and Ren, 2020)。該科昆蟲是一類重要的資源昆蟲,因其復(fù)變態(tài)等獨(dú)特的生物學(xué)特性而被學(xué)者廣泛關(guān)注,尤其是斑蝥素的藥用價值更是近年芫菁應(yīng)用研究的熱點(diǎn)之一(楊玉霞和任國棟, 2005; Bolognaetal., 2010)。然而,芫菁科部分物種的外部形態(tài)特征十分近似,在形態(tài)分類研究中經(jīng)常出現(xiàn)難以區(qū)分、錯誤鑒定的情況,如北方常見物種圓點(diǎn)斑芫菁Mylabrisaulica和蘋斑芫菁Mylabriscalida(Pan and Ren, 2020)、橫紋溝芫菁Hycleussolonicus(Panetal., 2017)等。分子系統(tǒng)學(xué)的發(fā)展,特別是DNA條形碼的提出和興起,為解決以上問題提供了可行的方案。Liu等(2016)基于線粒體COI基因片段序列,對西北豆芫菁Epicautasibirica、疑豆芫菁Epicautadubia和中華豆芫菁Epicautachinensis分類研究的有效性進(jìn)行了探討,提出后兩者是西北豆芫菁的次異名,然而該研究并未使用分子物種界定方法,僅是基于分支單系性得出的結(jié)論。

    為更加全面地評估分子物種界定方法在芫菁科分類研究中的可行性,本研究以中國北方常見芫菁科昆蟲為研究材料,選取線粒體COI基因片段和核CAD基因片段序列,分別采用自動條形碼間隔探索(automatic barcode gap discovery, ABGD)、廣義混合Yule溯祖模型(generalized mixed Yule coalescent,GMYC)、貝葉斯泊松樹進(jìn)程(Bayesian Poisson tree processes, bPTP)和貝葉斯系統(tǒng)發(fā)育和系統(tǒng)地理分析(Bayesian phylogenetics and phylogeography, BPP)方法進(jìn)行分子物種界定。在對比形態(tài)鑒定結(jié)果的基礎(chǔ)上,探討以上序列和分析方法界定的準(zhǔn)確性和適用性,以期為芫菁科的整合分類研究提供數(shù)據(jù)支持。

    1 材料與方法

    1.1 樣本選取

    本研究以中國北方常見芫菁科昆蟲為研究對象,選取58個樣本,經(jīng)形態(tài)鑒定,共涉及2亞科4族6屬18種;選擇擬步甲科(Tenebrionidae)黃粉蟲Tenebriomolitor為外群。實(shí)驗(yàn)樣本大多來自河北大學(xué)博物館館藏,少部分樣本序列下載自GenBank,詳見表1。

    1.2 形態(tài)鑒定

    選取樣本均使用Nikon SMZ1500立式顯微鏡進(jìn)行形態(tài)學(xué)鑒定。形態(tài)鑒定參考中國芫菁分類相關(guān)成果,如:Pan等(2010, 2017),潘昭等(2010a, 2010b, 2011, 2013), Li等(2020)以及Pan和Bologna(2021)。

    1.3 DNA的提取、PCR擴(kuò)增及測序

    從酒精泡制標(biāo)本中取后足肌肉組織,使用昆蟲基因組DNA提取試劑盒(倍沃醫(yī)學(xué)科技)提取總DNA。通過PCR擴(kuò)增COI和CAD基因片段,COI基因擴(kuò)增引物為TW-N1284(5′-AAACTAACARCCTTC AAAGYTGT-3′)和COI-R2329(5′-ACHGTAAAYATA TGATGGGCTCA-3′),反應(yīng)程序:94℃預(yù)變性3 min;94℃變性1 min, 48℃延伸1.5 min, 40個循環(huán); 72℃復(fù)性1.5 min;72℃最終延伸3 min。CAD基因擴(kuò)增引物為CD439F(5′-TTCAGTGTACARTTYCAYC CHGARCAYAC-3′)和CD688R(5′-TGTATACCTAG AGGATCDACRTTYTCCATRTTRCA-3′),反應(yīng)程序: 94℃預(yù)變性3 min; 94℃變性30 s, 55℃延伸30 s, 50個循環(huán); 72℃復(fù)性1 min; 72℃最終延伸3 min(Wild and Maddison, 2008)。PCR擴(kuò)增采用25 μL體系:上下游引物(10 μmol/L)各1 μL,去離子水8.5 μL, PCR Mix 12.5 μL,DNA模板2 μL。 PCR產(chǎn)物經(jīng)1%瓊脂糖凝膠電泳檢測,送通用生物系統(tǒng)(安徽)有限公司進(jìn)行雙向測序。

    1.4 基因序列處理和系統(tǒng)發(fā)育分析

    所得序列使用DNASTAR SeqMan v7.1.0進(jìn)行檢查和拼接,在NCBI中進(jìn)行BLAST分析,確定所得序列是否為目的序列, 使用MEGA X (Kumaretal., 2018)中的Clustal W進(jìn)行序列的比對。所有新得到的序列上傳至GenBank(登錄號見表1)。

    使用MEGA X分別對COI和CAD序列數(shù)據(jù)集進(jìn)行序列堿基含量分析及種內(nèi)和種間遺傳距離分析(Kumaretal., 2018)。使用PhyloSuite v1.2.2 (Zhangetal., 2020)將COI和CAD序列數(shù)據(jù)集串聯(lián),構(gòu)建3個序列數(shù)據(jù)集,即COI,CAD和COI+CAD串聯(lián)序列數(shù)據(jù)集。使用PhyloSuite內(nèi)嵌的MrBayes 3.2(Ronquistetal., 2012),使用貝葉斯推斷法(Bayesian inference, BI),分別基于3個數(shù)據(jù)集重建系統(tǒng)發(fā)育樹,運(yùn)行2 000 000世代,采樣頻率為每100代一次,確保分裂頻率分支頻率的標(biāo)準(zhǔn)偏差小于0.01(Ronquist and Huelsenbeck, 2003),以獲得高支持率的合意樹,生成合意樹之前舍棄25%的老化樣本,最佳堿基替代模型使用ModelFinder (Kalyaanamoorthyetal., 2017)進(jìn)行評估。 最后,使用Evolview v3 (Subramanianetal., 2019)對系統(tǒng)樹進(jìn)行可視化編輯。

    1.5 物種界定分析

    基于3個數(shù)據(jù)集,分別使用自動條形碼間隔探索(ABGD),廣義混合Yule溯祖模型(GMYC),貝葉斯泊松樹進(jìn)程(bPTP)和貝葉斯系統(tǒng)發(fā)育和系統(tǒng)地理分析(BPP)方法進(jìn)行物種界定分析,具體參數(shù)設(shè)置如下:

    ABGD:使用ABGD在線分析(https:∥bioinfo.mnhn.fr/abi/public/abgd/abgdweb.html)(金倩等, 2017),設(shè)置參數(shù)X(relative gap width)為0.5,其他參數(shù)為默認(rèn)值,將得到的序列劃分情況與形態(tài)鑒定的結(jié)果進(jìn)行比較分析。

    GMYC(金倩等, 2017):基于BEAST v1.10.4 (Drummond and Rambaut, 2007)構(gòu)建超度量二歧分支樹(有根),選擇嚴(yán)格分子鐘。MCMC 鏈為10 000萬代,每1 000代參數(shù)取樣一次,Burn-In默認(rèn)為鏈長10%,使用Tracer v1.7.2(Rambautetal., 2018)評估系統(tǒng)樹的收斂情況。將含有樹文件的*.trees.txt 文件導(dǎo)入到Tree Annotator v1.10.4生成所需系統(tǒng)樹。Burn In (as trees)根據(jù) Tracer v1.7.2分析結(jié)果(Burn-In 鏈長/采樣頻率)設(shè)置為1 000,最后生成所需最大分支可信度樹(maximum clade credibility tree),導(dǎo)入GMYC進(jìn)行分析。

    bPTP(Zhangetal., 2013):本實(shí)驗(yàn)使用PTP的更新版本bPTP進(jìn)行界定分析,將GMYC分析中得到的Newick系統(tǒng)樹上傳到PTP在線網(wǎng)站(https:∥species.h-its.org/ptp/)進(jìn)行分析,設(shè)置為有根樹,輸入外群名稱,馬爾科夫鏈代數(shù)為500 000,其他參數(shù)為默認(rèn)值。

    BPP(Yang, 2015):基于串聯(lián)基因數(shù)據(jù)集構(gòu)建的貝葉斯系統(tǒng)樹,在BPP v4.2.9中利用A11模板輸入指導(dǎo)樹(基于COI序列構(gòu)建的BI樹),運(yùn)行分析。

    2 結(jié)果

    2.1 基因序列特征

    COI基因片段序列經(jīng)拼接、比對、剪切處理后的長度為561 bp。使用MEGA X軟件中的Kimura雙參數(shù)模型(Kimura 2-parameter model)進(jìn)行序列堿基含量分析及種內(nèi)、種間遺傳距離分析(表2和3),結(jié)果顯示:序列的T, C, A和G平均含量分別為33.5%, 22.3%, 26.8%和17.5%,平均A+T含量顯著高于G+C含量,變異多發(fā)生在堿基第3位;Hycleussolonicus,Mylabrisaulica,Epicautamegalocephala,Mylabrissibirica,Lyttacaraganae,Meloeproscarabaeus,Meloeviolaceus和Epicautasibirica的種內(nèi)平均遺傳距離分別為0.92%, 0.71%, 0.83%, 0.30%, 1.66%, 1.34%, 0.11%和0.24%,其余10個物種樣本量不足,無法進(jìn)行種內(nèi)遺傳距離分析;18種芫菁的種間遺傳距離在9.49%~22.10%之間,平均值為18.36%。

    表2 基于COI基因片段序列計算得出芫菁科種內(nèi)平均遺傳距離和種內(nèi)遺傳距離范圍Table 2 Average intraspecific genetic distance and range of intraspecific genetic distance of Meloidaecalculated based on COI fragment sequences

    CAD基因片段序列處理后長度為730 bp。使用MEGA X 軟件中的Kimura 雙參數(shù)模型進(jìn)行序列堿基含量分析及種內(nèi)種間遺傳距離分析(表4和5),結(jié)果顯示:T, C, A和G平均含量分別為33.5%, 22.3%, 26.8%和17.5%,平均A+T含量顯著高于G+C含量,沒有堿基變異偏向性,密碼子的1, 2和3位均有發(fā)生變異;Hycleussolonicus,Mylabrisaulica,Epicautamegalocephala,Epicautasibirica,Lyttacaraganae,Meloeproscarabaeus和Meloeviolaceus的種內(nèi)平均遺傳距離為1.71%, 0.65%, 0.91%, 0.73%, 0.46%, 0.27%和0.06%,Mylabrissibirica沒有變異,其余10個物種樣品量不足,無法進(jìn)行種內(nèi)遺傳距離分析; 18種芫菁的種間遺傳距離是處于0.43%~22.33%之間,平均值為13.85%。

    表4 基于CAD基因片段序列計算得出芫菁科種內(nèi)平均遺傳距離和種內(nèi)遺傳距離范圍Table 4 Average intraspecific genetic distance and range of intraspecific genetic distance ofMeloidae calculated based on CAD fragment sequences

    2.2 系統(tǒng)發(fā)育

    利用PhyloSuite內(nèi)嵌的ModelFinder(Kalyaanamoorthyetal., 2017)進(jìn)行最佳堿基替代模型分析, 實(shí)現(xiàn)Akaike’sInformation Criterion (AIC)選擇最適模型,3個數(shù)據(jù)集的最適模型均為GTR+F+I+G4。

    對58頭芫菁樣本進(jìn)行BI分析,基于COI,CAD和COI+CAD串聯(lián)序列數(shù)據(jù)集分別建樹,顯示出高度相似的拓?fù)浣Y(jié)構(gòu)(圖1),聯(lián)合數(shù)據(jù)集分支支持度高于單基因建樹,高度支持這18個種的單系性。除Hycleussolonicus,Meloevariegatus,Meloebrevicollis和Epicautaobscurocephala這幾個分支的后驗(yàn)概率(posterior probability, PP)<1,其余種劃分PP均等于1。

    2.3 物種界定

    本實(shí)驗(yàn)對58個芫菁樣本(不含外群)的分子物種界定結(jié)果如圖1所示,具體如下:

    ABGD劃分結(jié)果包括初始劃分(initial partition)和遞歸劃分(recursive partition)兩種情況。 基于COI基因片段和COI+CAD串聯(lián)基因片段序列的劃分結(jié)果顯示,初始劃分較為穩(wěn)定,遺傳距離在0.001~0.1之間,被劃分為18個分子操作分類單元(molecular operational taxonomic unit, MOTU),支持18個形態(tài)種,且與BI樹拓?fù)浣Y(jié)構(gòu)一致。遞歸劃分先驗(yàn)值浮動較大,不夠穩(wěn)定?;贑AD基因片段序列的ABGD劃分結(jié)果顯示,屬級分類單元的劃分和形態(tài)劃分吻合,而物種階元上則將Meloeproscarabaeus和Meloeviolaceus劃分為一個MOTU,物種劃分紊亂,如圖2(A, B, C)。

    基于3個數(shù)據(jù)集的GMYC分析所劃分的MOTUs結(jié)果皆與形態(tài)種相符。單閾值分析,基于COI基因片段序列劃分的似然實(shí)體數(shù)量為19,置信區(qū)間為19~22,基于CAD基因片段序列劃分的似然實(shí)體數(shù)量為19,置信區(qū)間為1~49,基于COI+CAD串聯(lián)序列劃分的似然實(shí)體數(shù)量為19,置信區(qū)間為15~21,如圖3-5所示。

    bPTP分析結(jié)果表明,基于COI基因序列劃分得到21個MOTUs,將Mylabrisaulica,Meloeproscarabaeus和Lyttacaraganae各自分為了2個MOTUs;基于CAD基因序列劃分為17個MOTUs,將Meloeproscarabaeus和Meloeviolaceus歸為一個MOTU;基于COI+CAD基因串聯(lián)序列得到18個MOTUs,和形態(tài)學(xué)劃分結(jié)果一致,如表6-8所示。

    表6 基于COI序列通過bPTP方法得到的芫菁科物種劃分Table 6 Species division of Meloidae based on COI sequence by bPTP method

    表7 基于CAD序列通過bPTP方法得到的芫菁科物種劃分Table 7 Species division of Meloidae based on CAD sequence by bPTP method

    表8 基于COI+CAD串聯(lián)序列通過bPTP方法得到的芫菁科物種劃分Table 8 Species division of Meloidae based on concatenated COI+CAD sequence by bPTP method

    BPP是基于給定指導(dǎo)樹進(jìn)行分析,得到18個MOTUs的后驗(yàn)概率(PP)是0.98925,遠(yuǎn)大于其他劃分方式的PP,即支持將58個樣本劃分為18個獨(dú)立物種,支持指導(dǎo)樹給定的18個物種,與形態(tài)學(xué)劃分一致。

    3 討論和結(jié)論

    Liu等 (2016)首次應(yīng)用分子系統(tǒng)學(xué)方法對芫菁科昆蟲的物種劃分問題進(jìn)行了探討并提出了2個新異名,然而此工作依據(jù)的是外部形態(tài)、遺傳距離和系統(tǒng)樹拓?fù)浣Y(jié)構(gòu)的綜合結(jié)果,并未使用其他分子物種界定方法進(jìn)行分析。

    本研究對中國北方常見的芫菁科6屬18個形態(tài)種進(jìn)行分子物種界定,總體而言所得結(jié)果基本與形態(tài)劃分一致,僅基于CAD序列的ABGD分析結(jié)果和基于單基因的bPTP分析結(jié)果有所差異(圖1)。該結(jié)果一方面驗(yàn)證了形態(tài)種劃分結(jié)果,例如目前使用觸角末節(jié)形狀對具有相似鞘翅斑紋的圓點(diǎn)斑芫菁和蘋斑芫菁進(jìn)行區(qū)分(Pan and Ren, 2020),該形態(tài)界定結(jié)果得到了分子物種界定結(jié)果的支持;另一方面也顯現(xiàn)出對芫菁科昆蟲相對可靠的可用于分子物種界定的DNA序列和分析方法,即多基因片段串聯(lián)序列優(yōu)于COI序列優(yōu)于CAD序列,ABGD方法最簡便、準(zhǔn)確性相對較高且穩(wěn)定,在一定程度上可減少人為界定的誤差。在進(jìn)行ABGD和bPTP分析時,發(fā)現(xiàn)基于CAD的劃分對物種劃分不清晰,因此推測是由于CAD基因進(jìn)化速率較慢,種間差異小,所以導(dǎo)致種的劃分紊亂。

    此外,參考分子物種界定結(jié)果,基于COI基因種內(nèi)和種間的遺傳距離分析結(jié)果表明,在芫菁科中,當(dāng)遺傳距離在1.66%以下可認(rèn)為是同一物種的種內(nèi)變異,而大于9.49%則可認(rèn)為是不同的物種,而在1.66%~9.49%之間則需要綜合更多證據(jù)謹(jǐn)慎處理(表2和3)?;贑AD基因種內(nèi)和種間的遺傳距離和物種界定綜合分析結(jié)果表明,在芫菁科中,當(dāng)遺傳距離在1.71%以下可認(rèn)為是同一物種的種內(nèi)變異,而大于1.71%則有成為不同物種的可能,在1.71%~2.77%之間則需要綜合更多證據(jù)謹(jǐn)慎處理(表4和5)。

    綜上,基于COI和CAD的獨(dú)立和串聯(lián)數(shù)據(jù)集和ABGD, GMYC, PTP和BPP 4種分析方法均適用于芫菁科昆蟲的分子物種界定,但有一定的誤差,需多基因、多方法綜合考慮。此外,形態(tài)和分子物種界定方法的結(jié)合能夠有效提高芫菁科昆蟲物種劃分的準(zhǔn)確性。

    猜你喜歡
    物種分析
    物種大偵探
    物種大偵探
    吃光入侵物種真的是解決之道嗎?
    英語世界(2023年10期)2023-11-17 09:18:18
    隱蔽失效適航要求符合性驗(yàn)證分析
    回首2018,這些新物種值得關(guān)注
    電力系統(tǒng)不平衡分析
    電子制作(2018年18期)2018-11-14 01:48:24
    電咖再造新物種
    汽車觀察(2018年10期)2018-11-06 07:05:26
    電力系統(tǒng)及其自動化發(fā)展趨勢分析
    瘋狂的外來入侵物種
    中西醫(yī)結(jié)合治療抑郁癥100例分析
    .国产精品久久| 国产黄片美女视频| 日韩国内少妇激情av| 久久久久国产精品人妻aⅴ院| 日日摸夜夜添夜夜添av毛片 | 国产v大片淫在线免费观看| 国产精品伦人一区二区| 又爽又黄a免费视频| 日日摸夜夜添夜夜添小说| 搞女人的毛片| 我的老师免费观看完整版| 日本免费一区二区三区高清不卡| eeuss影院久久| 色吧在线观看| 国产黄a三级三级三级人| 国产一区二区在线av高清观看| 男女边吃奶边做爰视频| 18禁在线播放成人免费| 嫩草影院入口| 一级黄片播放器| 女人被狂操c到高潮| 一本精品99久久精品77| 国产成人一区二区在线| 99视频精品全部免费 在线| 日韩精品青青久久久久久| 一级a爱片免费观看的视频| 久久久成人免费电影| 色哟哟·www| 一级毛片久久久久久久久女| 久久久国产成人免费| 91精品国产九色| 午夜精品一区二区三区免费看| 欧美+日韩+精品| 久久精品91蜜桃| 最好的美女福利视频网| 国产伦精品一区二区三区视频9| 亚洲精品色激情综合| 在线观看美女被高潮喷水网站| 最近最新免费中文字幕在线| 国产淫片久久久久久久久| 老司机午夜福利在线观看视频| 亚洲在线观看片| 亚洲最大成人中文| av在线观看视频网站免费| 免费av毛片视频| 露出奶头的视频| 男女视频在线观看网站免费| 久久久国产成人精品二区| 老女人水多毛片| 国产亚洲精品综合一区在线观看| 亚洲最大成人中文| 亚洲美女黄片视频| 久久久久久国产a免费观看| 国产精品久久久久久av不卡| 又爽又黄无遮挡网站| 国内精品久久久久精免费| 不卡视频在线观看欧美| 22中文网久久字幕| 色在线成人网| 嫩草影视91久久| 久久精品国产自在天天线| 国产一区二区激情短视频| 免费搜索国产男女视频| 两人在一起打扑克的视频| 九色国产91popny在线| 亚洲精品久久国产高清桃花| 久久久久国内视频| 欧美绝顶高潮抽搐喷水| 蜜桃亚洲精品一区二区三区| 中出人妻视频一区二区| 国产精品久久电影中文字幕| 久久草成人影院| 午夜精品在线福利| 日韩亚洲欧美综合| 亚洲欧美激情综合另类| 99久久精品国产国产毛片| 国产麻豆成人av免费视频| 少妇人妻精品综合一区二区 | 亚洲av美国av| 亚洲精华国产精华精| 国产日本99.免费观看| 中文在线观看免费www的网站| 国产成人一区二区在线| 久久久久久久久中文| 高清日韩中文字幕在线| 欧美人与善性xxx| 国产在线精品亚洲第一网站| 欧美区成人在线视频| 亚洲乱码一区二区免费版| 国产精品,欧美在线| 国产伦人伦偷精品视频| 国产美女午夜福利| av专区在线播放| 免费av不卡在线播放| 男女啪啪激烈高潮av片| 亚洲av免费在线观看| 少妇裸体淫交视频免费看高清| 男人的好看免费观看在线视频| 麻豆av噜噜一区二区三区| 联通29元200g的流量卡| 国产精品久久久久久久电影| 免费av毛片视频| 精品免费久久久久久久清纯| 美女 人体艺术 gogo| 国产大屁股一区二区在线视频| 日韩国内少妇激情av| 免费在线观看影片大全网站| 精品午夜福利在线看| 五月玫瑰六月丁香| 精品一区二区三区视频在线| 日本一本二区三区精品| 黄色一级大片看看| www日本黄色视频网| 噜噜噜噜噜久久久久久91| 亚洲精品影视一区二区三区av| 一级黄片播放器| 中文在线观看免费www的网站| 窝窝影院91人妻| 一区二区三区激情视频| 久久婷婷人人爽人人干人人爱| 黄色视频,在线免费观看| 亚洲四区av| 亚洲国产精品久久男人天堂| 一级黄色大片毛片| 午夜福利在线在线| 女生性感内裤真人,穿戴方法视频| 亚洲欧美日韩卡通动漫| 精品午夜福利视频在线观看一区| 国产激情偷乱视频一区二区| 日本黄大片高清| 婷婷精品国产亚洲av| 国产男靠女视频免费网站| 亚洲精品色激情综合| 中亚洲国语对白在线视频| 成年免费大片在线观看| 色5月婷婷丁香| 一进一出抽搐gif免费好疼| 亚洲综合色惰| 日日干狠狠操夜夜爽| 一进一出抽搐动态| 九色国产91popny在线| 看十八女毛片水多多多| 国产免费一级a男人的天堂| 男女啪啪激烈高潮av片| 我要看日韩黄色一级片| 夜夜夜夜夜久久久久| 联通29元200g的流量卡| 我的老师免费观看完整版| 美女高潮喷水抽搐中文字幕| 午夜老司机福利剧场| 欧美另类亚洲清纯唯美| 午夜精品一区二区三区免费看| av天堂中文字幕网| 精品人妻视频免费看| 免费搜索国产男女视频| 欧美成人免费av一区二区三区| 国产伦一二天堂av在线观看| 欧美bdsm另类| 国产黄色小视频在线观看| 美女高潮的动态| 午夜日韩欧美国产| 久久久久久九九精品二区国产| 国产成年人精品一区二区| 美女cb高潮喷水在线观看| 色尼玛亚洲综合影院| 人妻制服诱惑在线中文字幕| av在线观看视频网站免费| 国产一区二区三区在线臀色熟女| 可以在线观看的亚洲视频| 日本一本二区三区精品| 麻豆久久精品国产亚洲av| 高清在线国产一区| 91在线观看av| 九九久久精品国产亚洲av麻豆| 性色avwww在线观看| 99久久精品一区二区三区| 国产精品无大码| 亚洲va在线va天堂va国产| 看片在线看免费视频| 久久精品国产清高在天天线| 午夜福利视频1000在线观看| 色5月婷婷丁香| 一进一出好大好爽视频| 久久99热6这里只有精品| 国产高清不卡午夜福利| 制服丝袜大香蕉在线| 老师上课跳d突然被开到最大视频| 中文字幕av成人在线电影| 国产精品女同一区二区软件 | 亚洲av中文av极速乱 | 久久精品91蜜桃| 18禁在线播放成人免费| 一本精品99久久精品77| 国产真实伦视频高清在线观看 | 此物有八面人人有两片| 18禁裸乳无遮挡免费网站照片| 18+在线观看网站| 成人国产麻豆网| 久99久视频精品免费| 又黄又爽又免费观看的视频| 九色成人免费人妻av| 一级av片app| 国产一级毛片七仙女欲春2| 国内毛片毛片毛片毛片毛片| 99热精品在线国产| 99热网站在线观看| 国产单亲对白刺激| 国产乱人视频| 午夜精品在线福利| 日韩,欧美,国产一区二区三区 | 日韩国内少妇激情av| 国产精品人妻久久久影院| 成人鲁丝片一二三区免费| 国产亚洲精品久久久久久毛片| 九九热线精品视视频播放| 欧美高清性xxxxhd video| 国产乱人视频| 日本欧美国产在线视频| 看十八女毛片水多多多| 亚洲av一区综合| 观看免费一级毛片| 黄片wwwwww| 麻豆成人午夜福利视频| 在线观看66精品国产| 嫩草影院精品99| 观看免费一级毛片| 亚洲国产精品sss在线观看| 国内精品宾馆在线| 国产黄片美女视频| 亚洲av第一区精品v没综合| 国内精品久久久久精免费| 精华霜和精华液先用哪个| 午夜精品久久久久久毛片777| 欧美日韩瑟瑟在线播放| 丰满乱子伦码专区| 国产色婷婷99| 免费av毛片视频| 成年免费大片在线观看| 大型黄色视频在线免费观看| 日韩亚洲欧美综合| 亚洲中文字幕日韩| 亚洲国产精品成人综合色| 人妻丰满熟妇av一区二区三区| 一区二区三区高清视频在线| 国产免费一级a男人的天堂| 小蜜桃在线观看免费完整版高清| .国产精品久久| 精品一区二区免费观看| 国产精品久久视频播放| 亚洲第一电影网av| 一区二区三区高清视频在线| 精品久久国产蜜桃| 一本久久中文字幕| 欧美极品一区二区三区四区| 免费无遮挡裸体视频| 日韩强制内射视频| 十八禁国产超污无遮挡网站| 精华霜和精华液先用哪个| 国产淫片久久久久久久久| 亚洲经典国产精华液单| 少妇裸体淫交视频免费看高清| 联通29元200g的流量卡| 欧美黑人欧美精品刺激| 别揉我奶头 嗯啊视频| 日本五十路高清| 性色avwww在线观看| 国产精品亚洲美女久久久| 亚洲av成人av| 十八禁国产超污无遮挡网站| 亚洲综合色惰| 精品99又大又爽又粗少妇毛片 | 成人性生交大片免费视频hd| 亚洲专区国产一区二区| 少妇猛男粗大的猛烈进出视频 | 精品国内亚洲2022精品成人| 少妇丰满av| 99热6这里只有精品| 免费观看在线日韩| a在线观看视频网站| 免费看光身美女| 国国产精品蜜臀av免费| 国产精品电影一区二区三区| 久久精品国产清高在天天线| 精品日产1卡2卡| 97碰自拍视频| 乱系列少妇在线播放| 欧美区成人在线视频| 亚洲七黄色美女视频| 亚洲精华国产精华液的使用体验 | 亚洲美女视频黄频| 中文字幕免费在线视频6| av中文乱码字幕在线| 亚洲国产色片| 亚洲一区二区三区色噜噜| 欧美xxxx黑人xx丫x性爽| 特级一级黄色大片| 国产精品女同一区二区软件 | 日韩av在线大香蕉| 欧美激情久久久久久爽电影| 一个人免费在线观看电影| 真人做人爱边吃奶动态| 亚洲av成人av| 中出人妻视频一区二区| 亚洲精品一卡2卡三卡4卡5卡| 黄片wwwwww| 国产激情偷乱视频一区二区| 日本免费一区二区三区高清不卡| 俄罗斯特黄特色一大片| 99在线视频只有这里精品首页| 国产精品野战在线观看| 亚洲国产精品sss在线观看| 欧美色视频一区免费| 真人做人爱边吃奶动态| 国产精品无大码| 女人十人毛片免费观看3o分钟| 国产久久久一区二区三区| 欧美区成人在线视频| 亚洲国产欧洲综合997久久,| 亚洲欧美激情综合另类| netflix在线观看网站| 亚洲精品亚洲一区二区| 精品免费久久久久久久清纯| 亚洲第一电影网av| 亚洲一区二区三区色噜噜| 国产亚洲91精品色在线| 麻豆av噜噜一区二区三区| 日本成人三级电影网站| 丰满人妻一区二区三区视频av| 国产精品久久久久久久电影| 欧美日韩中文字幕国产精品一区二区三区| 国产成人av教育| 不卡视频在线观看欧美| 搡女人真爽免费视频火全软件 | 日韩大尺度精品在线看网址| 久久国产精品人妻蜜桃| 国产乱人伦免费视频| 久久精品国产清高在天天线| 亚洲四区av| 亚洲精华国产精华液的使用体验 | 麻豆精品久久久久久蜜桃| 看免费成人av毛片| 国产成年人精品一区二区| 很黄的视频免费| 国产乱人伦免费视频| 色哟哟·www| 老司机福利观看| 精品一区二区三区视频在线观看免费| 国产一区二区亚洲精品在线观看| 亚洲成a人片在线一区二区| 国产午夜福利久久久久久| 亚洲精品久久国产高清桃花| av.在线天堂| 一个人观看的视频www高清免费观看| 人妻制服诱惑在线中文字幕| 免费搜索国产男女视频| av天堂中文字幕网| 99久久精品国产国产毛片| 亚洲成人免费电影在线观看| 精品欧美国产一区二区三| 亚洲电影在线观看av| 午夜免费激情av| 国产精品亚洲一级av第二区| 一边摸一边抽搐一进一小说| 中文亚洲av片在线观看爽| 97热精品久久久久久| 亚洲经典国产精华液单| 嫁个100分男人电影在线观看| 在线天堂最新版资源| 中国美女看黄片| 亚洲精品一卡2卡三卡4卡5卡| 亚洲五月天丁香| 国产白丝娇喘喷水9色精品| 久久久久久久精品吃奶| 麻豆成人av在线观看| 在线观看午夜福利视频| 欧美三级亚洲精品| 成人av在线播放网站| 国产高潮美女av| 国产精品一区www在线观看 | 午夜福利成人在线免费观看| 国产免费av片在线观看野外av| 日韩欧美在线二视频| 国产免费av片在线观看野外av| 国内久久婷婷六月综合欲色啪| 九九久久精品国产亚洲av麻豆| 高清在线国产一区| 五月玫瑰六月丁香| 波多野结衣巨乳人妻| 噜噜噜噜噜久久久久久91| 中国美白少妇内射xxxbb| 99九九线精品视频在线观看视频| 搞女人的毛片| 欧美成人a在线观看| 婷婷精品国产亚洲av| 成人无遮挡网站| 免费观看在线日韩| 能在线免费观看的黄片| 免费看a级黄色片| 久久精品国产99精品国产亚洲性色| 亚洲成人中文字幕在线播放| 久久久久久久久久久丰满 | 免费黄网站久久成人精品| 久久精品国产亚洲av香蕉五月| 国产在视频线在精品| 久久久久久久久久黄片| 亚洲av中文字字幕乱码综合| 欧美性猛交╳xxx乱大交人| 91狼人影院| 日本a在线网址| 国产白丝娇喘喷水9色精品| 欧美xxxx黑人xx丫x性爽| 夜夜看夜夜爽夜夜摸| 麻豆国产av国片精品| 老熟妇仑乱视频hdxx| bbb黄色大片| 最好的美女福利视频网| 国产精品久久电影中文字幕| 亚洲av美国av| 中文字幕av在线有码专区| 舔av片在线| 亚洲va日本ⅴa欧美va伊人久久| 免费人成视频x8x8入口观看| av.在线天堂| 国产精品一区www在线观看 | 最近最新中文字幕大全电影3| 88av欧美| 看黄色毛片网站| 亚洲狠狠婷婷综合久久图片| 成人高潮视频无遮挡免费网站| 夜夜爽天天搞| 欧美三级亚洲精品| 中文字幕精品亚洲无线码一区| 国产av在哪里看| 免费看日本二区| 久久这里只有精品中国| 人妻丰满熟妇av一区二区三区| 久久久久久大精品| 日韩,欧美,国产一区二区三区 | 亚洲精品亚洲一区二区| 91麻豆av在线| 丰满人妻一区二区三区视频av| 国产精品电影一区二区三区| 我要搜黄色片| 精品不卡国产一区二区三区| 午夜精品久久久久久毛片777| av天堂在线播放| 精品日产1卡2卡| 久久精品国产清高在天天线| 大又大粗又爽又黄少妇毛片口| 久久久精品欧美日韩精品| 看免费成人av毛片| 中文字幕免费在线视频6| 国产麻豆成人av免费视频| 日韩中文字幕欧美一区二区| 男女之事视频高清在线观看| 日韩在线高清观看一区二区三区 | АⅤ资源中文在线天堂| 91麻豆av在线| 色在线成人网| 久久99热6这里只有精品| 午夜福利高清视频| 国产亚洲精品久久久com| 三级男女做爰猛烈吃奶摸视频| 亚洲欧美日韩卡通动漫| 国产免费男女视频| 最近视频中文字幕2019在线8| 日韩一区二区视频免费看| 国产精品一区二区三区四区久久| 99久久精品一区二区三区| 制服丝袜大香蕉在线| 亚洲精品一卡2卡三卡4卡5卡| 好男人在线观看高清免费视频| 男人舔奶头视频| 老司机午夜福利在线观看视频| 香蕉av资源在线| a级一级毛片免费在线观看| 超碰av人人做人人爽久久| 伦理电影大哥的女人| 春色校园在线视频观看| 久久久午夜欧美精品| 国产视频内射| 午夜福利18| 如何舔出高潮| 国产高清不卡午夜福利| 久久久久久久久大av| 国产淫片久久久久久久久| 国产在线男女| 麻豆国产97在线/欧美| 精品一区二区三区av网在线观看| 精品一区二区免费观看| 精品免费久久久久久久清纯| 少妇高潮的动态图| 嫩草影院入口| 伦精品一区二区三区| 变态另类丝袜制服| 人妻丰满熟妇av一区二区三区| 国产精品人妻久久久影院| 窝窝影院91人妻| 成人特级av手机在线观看| 99久久精品一区二区三区| 舔av片在线| 久久久久久久久中文| 俄罗斯特黄特色一大片| 中国美白少妇内射xxxbb| 极品教师在线免费播放| 国产成人一区二区在线| 在线观看美女被高潮喷水网站| 国产探花极品一区二区| 天美传媒精品一区二区| 国产成人一区二区在线| 少妇裸体淫交视频免费看高清| 特级一级黄色大片| 狠狠狠狠99中文字幕| 日本熟妇午夜| 亚洲七黄色美女视频| 国产精品国产高清国产av| 亚洲黑人精品在线| 免费av不卡在线播放| 午夜老司机福利剧场| 成人特级黄色片久久久久久久| 最后的刺客免费高清国语| 夜夜看夜夜爽夜夜摸| 成年人黄色毛片网站| 亚洲成人久久性| 人妻少妇偷人精品九色| 男人舔奶头视频| 成年免费大片在线观看| 精品人妻视频免费看| 老司机午夜福利在线观看视频| 亚洲精品成人久久久久久| 91麻豆av在线| 男女做爰动态图高潮gif福利片| 尤物成人国产欧美一区二区三区| 免费电影在线观看免费观看| 国产 一区 欧美 日韩| 国产高清视频在线观看网站| 99久久精品一区二区三区| 极品教师在线免费播放| 亚洲欧美日韩东京热| 亚洲欧美日韩卡通动漫| 亚洲av二区三区四区| 国模一区二区三区四区视频| av.在线天堂| 人妻少妇偷人精品九色| 99国产极品粉嫩在线观看| 性欧美人与动物交配| 国产不卡一卡二| 国内久久婷婷六月综合欲色啪| 亚洲综合色惰| 天堂网av新在线| 综合色av麻豆| 免费看光身美女| 色精品久久人妻99蜜桃| 久久精品91蜜桃| 国产美女午夜福利| 午夜福利18| 内地一区二区视频在线| 色综合亚洲欧美另类图片| 国产aⅴ精品一区二区三区波| 亚洲精品色激情综合| 国产激情偷乱视频一区二区| av在线天堂中文字幕| 国产精品99久久久久久久久| 一区二区三区四区激情视频 | 嫩草影院新地址| 免费无遮挡裸体视频| 久久久久九九精品影院| 成人av一区二区三区在线看| 午夜激情欧美在线| 91av网一区二区| 国产av一区在线观看免费| 欧美精品国产亚洲| netflix在线观看网站| 91在线精品国自产拍蜜月| 变态另类成人亚洲欧美熟女| 在现免费观看毛片| av天堂在线播放| 校园春色视频在线观看| 国产精品国产三级国产av玫瑰| 国产伦精品一区二区三区四那| 欧美性感艳星| 欧美日韩瑟瑟在线播放| 99riav亚洲国产免费| 两个人的视频大全免费| 99热6这里只有精品| 久久精品国产自在天天线| 十八禁国产超污无遮挡网站| 哪里可以看免费的av片| 色综合婷婷激情| 欧美激情在线99| 亚洲色图av天堂| 国内精品宾馆在线| 欧美日本亚洲视频在线播放| 亚洲国产精品久久男人天堂| 天堂√8在线中文| 国产伦在线观看视频一区| 国模一区二区三区四区视频| 国内少妇人妻偷人精品xxx网站| 麻豆av噜噜一区二区三区| 久久99热这里只有精品18| 窝窝影院91人妻| 日本a在线网址| 久久精品91蜜桃| 少妇人妻精品综合一区二区 | 国产精品美女特级片免费视频播放器| 村上凉子中文字幕在线| av黄色大香蕉| 九九久久精品国产亚洲av麻豆| 美女被艹到高潮喷水动态| 变态另类丝袜制服| 日韩国内少妇激情av| 久久国内精品自在自线图片| 亚洲专区中文字幕在线| 老司机午夜福利在线观看视频| 日日撸夜夜添|