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

    長(zhǎng)江上游中華金沙鰍和短身金沙鰍線粒體遺傳多樣性研究

    2020-06-06 05:46:08龍安雨田輝伍汪登強(qiáng)陳大慶周湖海段辛斌
    淡水漁業(yè) 2020年3期
    關(guān)鍵詞:控制區(qū)金沙魚(yú)類(lèi)

    龍安雨,田輝伍,汪登強(qiáng),陳大慶,周湖海,3,段辛斌

    (1.上海海洋大學(xué),水產(chǎn)科學(xué)國(guó)家級(jí)實(shí)驗(yàn)教學(xué)示范中心,上海 201306; 2.中國(guó)水產(chǎn)科學(xué)研究院長(zhǎng)江水產(chǎn)研究所,農(nóng)業(yè)農(nóng)村部長(zhǎng)江中上游漁業(yè)資源環(huán)境科學(xué)觀測(cè)實(shí)驗(yàn)站,武漢 430223;3.華中農(nóng)業(yè)大學(xué)水產(chǎn)學(xué)院,武漢 430070)

    長(zhǎng)江上游指宜昌以上江段,該江段地形地貌復(fù)雜、河流水系眾多,魚(yú)類(lèi)區(qū)系與物種組成具有多樣性、特有性特征,有魚(yú)類(lèi)200多種,其中特有魚(yú)類(lèi)100多種[1]。近年來(lái),隨著經(jīng)濟(jì)和社會(huì)的發(fā)展,該江段大型水利工程相繼建設(shè),導(dǎo)致水文等環(huán)境條件發(fā)生改變,魚(yú)類(lèi)適宜生境不斷破碎、萎縮,種群數(shù)量急劇下降[1],魚(yú)類(lèi)資源和遺傳多樣性保護(hù)已成為長(zhǎng)江上游水電開(kāi)發(fā)過(guò)程中受關(guān)注的熱點(diǎn)之一。

    中華金沙鰍(Jinshaiasinensis)和短身金沙鰍(J.abbreviata)是長(zhǎng)江上游的特有魚(yú)類(lèi),是金沙鰍屬僅有的兩個(gè)種[2],其中中華金沙鰍也是“長(zhǎng)江上游珍稀特有魚(yú)類(lèi)國(guó)家級(jí)自然保護(hù)區(qū)”的保護(hù)對(duì)象之一[3]。兩種金沙鰍均為底棲、雜食性、喜激流、產(chǎn)漂流性卵魚(yú)類(lèi),其生活類(lèi)型是長(zhǎng)江上游最主要的生態(tài)類(lèi)群,具有代表性。兩種金沙鰍不僅有相似的形態(tài)特征和生活習(xí)性,而且分布區(qū)域基本重疊,主要分布在長(zhǎng)江上游干流、金沙江、雅礱江、嘉陵江和岷江[2,4]等水域。但是在資源量上,中華金沙鰍和短身金沙鰍有明顯差別,根據(jù)魚(yú)類(lèi)早期資源調(diào)查,中華金沙鰍規(guī)模大于短身金沙鰍[4],表明這兩種魚(yú)對(duì)環(huán)境的適應(yīng)力存在差異。因此研究群體遺傳多樣性有助于闡明近緣同域分布種類(lèi)的適應(yīng)性進(jìn)化。本研究采用線粒體DNA (mtDNA)控制區(qū)序列對(duì)長(zhǎng)江上游兩種金沙鰍遺傳多樣性進(jìn)行分析,旨在了解兩種金沙鰍的遺傳多樣性現(xiàn)狀,為長(zhǎng)江上游金沙鰍屬魚(yú)類(lèi)的資源保護(hù)提供科學(xué)依據(jù)。

    1 材料與方法

    1.1 材料

    2016-2018年在長(zhǎng)江上游8個(gè)樣點(diǎn)采集中華金沙鰍136尾,3個(gè)樣點(diǎn)采集短身金沙鰍66尾。樣本采集地見(jiàn)圖1,樣本量見(jiàn)表1。所有樣本經(jīng)形態(tài)鑒定[2]后剪取少量鰭條并保存于無(wú)水乙醇中備用。

    圖1 中華金沙鰍和短身金沙鰍采樣點(diǎn)及中華金沙鰍線粒體DNA控制區(qū)譜系分布Fig.1 Sampling sites and mtDNA clad distribution of J.sinensis and J.abbreviate

    1.2 基因組DNA 提取、PCR 擴(kuò)增及測(cè)序

    采用高鹽法提取基因組DNA[5]。線粒體控制區(qū)PCR擴(kuò)增采用申紹祎等[6]研究所用引物。PCR的擴(kuò)增體系為50 μL,包含2×PCR Mix(北京擎科新業(yè)生物技術(shù)有限公司)25 μL,上下游引物(10 pmol/L)各2 μL,DNA模板2 μL,滅菌去離子水加至50 μL。PCR擴(kuò)增條件:94 ℃預(yù)變性3 min;94 ℃變性30 s,55 ℃退火30 s,72 ℃延伸60 s,35個(gè)循環(huán);72 ℃再延伸8 min。PCR結(jié)束后,取2.5 μL產(chǎn)物,用1.0%瓊脂糖凝膠電泳檢測(cè)PCR,將擴(kuò)增效果良好的PCR產(chǎn)物送武漢天一輝遠(yuǎn)生物科技有限公司進(jìn)行雙向測(cè)序,測(cè)序引物與PCR引物相同。

    1.3 數(shù)據(jù)處理和分析

    使用Lasergene v7.1軟件包中的SeqMan對(duì)雙向測(cè)序結(jié)果進(jìn)行序列拼接,Clustal X[7]軟件進(jìn)行對(duì)位排列。用 MEGA v6.0[8]分析軟件統(tǒng)計(jì)序列的平均堿基組成、計(jì)算轉(zhuǎn)換/顛換比率。在DnaSP v5.0[9]軟件中獲得單倍型數(shù)(H)、變異位點(diǎn)數(shù)(S)、簡(jiǎn)約信息位點(diǎn)數(shù),計(jì)算單倍型多樣性指數(shù)(Hd)和核苷酸多樣性指數(shù)(Pi)等參數(shù),并計(jì)算群體間的K2P遺傳距離,。

    應(yīng)用Arlequin v3.11[10]軟件計(jì)算群體間的遺傳分化指數(shù)(FST)值,進(jìn)行分子方差分析(AMOVA),設(shè)1 000次重復(fù)隨機(jī)抽樣作顯著性檢驗(yàn)。用公式Nm=(0.5/FST)-0.5計(jì)算群體間基因流[11]。在MEGA v6.0[7]分析軟件中構(gòu)建單倍型鄰接(NJ)發(fā)育樹(shù),bootstrap驗(yàn)證各分支支持率,設(shè)置為1 000。同時(shí)采用MrBayes v 3.1.2[12]軟件對(duì)單倍型之間的系統(tǒng)發(fā)生樹(shù)進(jìn)行了貝葉斯推理(BI),最佳堿基替換模型由在線軟件IQtree (http://iqtree.cibiv.univie.ac.at/)獲得。使用Network5.0程序[13]構(gòu)建單倍型網(wǎng)絡(luò)圖,對(duì)單倍型進(jìn)化網(wǎng)絡(luò)關(guān)系進(jìn)行分析。

    在DnaSP v5.0軟件中計(jì)算中性檢驗(yàn)值Tajima′s D和Fu′s Fs。用Beast 2.4.3.0[14]進(jìn)行Bayesian Skyline Plots (BSP)分析種群歷史,MCMC設(shè)為600 000 000,所得結(jié)果由TRACER ver.1.4 讀取,各項(xiàng)參數(shù)的ESS均大于200。

    2 結(jié)果與分析

    2.1 序列變異

    序列比對(duì)后共獲得中華金沙鰍136條控制區(qū)序列,有效序列長(zhǎng)度為915 bp,A、T、C、G的平均含量為35.63%、30.19%、20.02%和14.15%,具有明顯的反G含量?jī)A向,Ts/Tv =23.63,轉(zhuǎn)換數(shù)明顯高于顛換數(shù)。共獲得66尾短身金沙鰍控制區(qū)序列,比對(duì)后得到有效長(zhǎng)度為908 bp,A、T、C、G的平均含量為35.82%、31.39%、19.49%和13.30%,Ts/Tv =12.87,轉(zhuǎn)換數(shù)也明顯高于顛換數(shù)。

    中華金沙鰍序列中共發(fā)現(xiàn)132個(gè)變異位點(diǎn),其中39個(gè)為單一突變位點(diǎn),93個(gè)為簡(jiǎn)約信息位點(diǎn)。136條序列共定義了131個(gè)單倍型,其中頻率在2以上的單倍型只有5個(gè),其余單倍型的頻率均為1??倶颖竞塑账岫鄻有?Pi)和單倍型多樣性指數(shù)(Hd)分別為0.018 2和0.999(表1),各群體之間的遺傳多樣性差別較小。短身金沙鰍序列僅檢測(cè)到2個(gè)變異位點(diǎn),定義了3個(gè)單倍型,Pi和Hd分別為0.000 5和0.411(表1)。

    表1 中華金沙鰍和短身金沙鰍采樣點(diǎn)、樣本量、遺傳多樣性及中性檢驗(yàn)值Tab.1 Sampling sites,sample size,genetic diversity parameters and neutral test of J.sinensis and J.abbreviate populations

    注:*表示顯著性檢驗(yàn)P<0.05,**表示顯著性檢驗(yàn)P<0.01。

    2.2 遺傳結(jié)構(gòu)

    基于K2P模型計(jì)算,中華金沙鰍種內(nèi)單倍型之間的遺傳距離為0.001 1~0.039 6,短身金沙鰍種內(nèi)單倍型之間的遺傳距離為0.001 1~0.002 2,這兩種魚(yú)種間的遺傳距離為0.035 7~0.051 5。

    中華金沙鰍群體間分子變異方差分析(AMOVA)結(jié)果顯示,來(lái)自群體間遺傳變異組成占-0.75%,來(lái)自群體內(nèi)的變異組成占100.75%,群體總的遺傳分化指數(shù)(FST)為-0.007 53,說(shuō)明中華金沙鰍群體遺傳變異主要來(lái)自群體內(nèi),群體間未出現(xiàn)顯著遺傳分化(表2)。

    短身金沙鰍群體間AMOVA分析結(jié)果也顯示群體間未發(fā)生顯著遺傳分化,其變異主要來(lái)源于群體內(nèi)(98.72%),且FST< 0.05 (FST=0.012 78)(表2)。

    分別計(jì)算兩兩群體間的FST,結(jié)果顯示,中華金沙鰍除NN及QJ與部分群體間的FST值大于0.05外,其余群體間FST均小于0.05(表3)。短身金沙鰍顯示,除BN與HJ群體間FST大于0.05外,其余群體間FST均小于0.05(表4)。

    表2 中華金沙鰍和短身金沙鰍群體間分子變異分析Tab.2 Analysis of molecular variance (AMOVA)among populations of J.sinensis and J.abbreviata

    中華金沙鰍和短身金沙鰍群體間基因流分別見(jiàn)表3和表4。這兩種魚(yú)群體間的Nm值(或絕對(duì)值)均大于4(除JJ與NN、BN與HJ外)表明群體間基因交流比較頻繁。

    利用Median Joining方法構(gòu)建中華金沙鰍單倍型網(wǎng)絡(luò)結(jié)構(gòu)圖(圖2),可見(jiàn)單倍型關(guān)系復(fù)雜,無(wú)明顯的中心單倍型。根據(jù)連接單倍型之間的突變步長(zhǎng),可將單倍型劃分為3個(gè)譜系:Clad I、Clad II、

    表3 中華金沙鰍群體間的遺傳分化(對(duì)角線下方)和基因流Nm值(對(duì)角線上方)Tab.3 Pairwise FST (below diagonal)and Nm values (above diagonal)among populations of J.sinensis

    表4 短身金沙鰍種群間的遺傳分化(對(duì)角線下方)和基因流Nm值(對(duì)角線上方)Tab.4 Pairwise FST (below diagonal)and Nm values (above diagonal)among populations of J.abbreviata

    Clad III譜系間的步長(zhǎng)為10和11,譜系內(nèi)相鄰單倍型之間突變步長(zhǎng)都低于10(除了Clad II內(nèi)的一個(gè)單倍型)。以犁頭鰍(Lepturichthysfimbriata,GenBank登錄號(hào):DQ105283.1)作為外類(lèi)群,構(gòu)建兩種金沙鰍單倍型NJ樹(shù)和BI樹(shù)(圖3)。兩棵樹(shù)顯示相似的拓?fù)浣Y(jié)構(gòu),中華金沙鰍和短身金沙鰍的單倍型分為兩個(gè)獨(dú)立分支,形成兩個(gè)單系群。中華金沙鰍的單倍型又分為3個(gè)具有較高支持率的支系。這3個(gè)支系包含的單倍型與中華金沙鰍的Network網(wǎng)絡(luò)結(jié)構(gòu)一致。

    2.3 種群歷史

    通過(guò)中性檢驗(yàn)Tajima′s D和Fu′s Fs等方法分析種群歷史[15],中華金沙鰍結(jié)果顯示,總?cè)后w以及3個(gè)譜系群體的Tajima′s D值均不顯著,而Fu′s Fs檢驗(yàn)中均為極顯著負(fù)值(表1)。Bayesian skyline plot (BSP)進(jìn)一步分析表明(圖4),中華金沙鰍總?cè)后w和3個(gè)譜系群體的擴(kuò)張時(shí)間分別約在距今0.12~0.17、0.1~0.12和0.12~0.17百萬(wàn)年(Ma)。短身金沙鰍群體由于序列變異位點(diǎn)和單倍型數(shù)量少,不作種群歷史分析。

    圖2 基于線粒體DNA控制區(qū)序列構(gòu)建的中華金沙鰍單倍型網(wǎng)絡(luò)結(jié)構(gòu)圖(短線代表增加一個(gè)突變步驟)Fig.2 Haplotype network of J.sinensis based on mtDNA control region sequences

    圖3 基于線粒體DNA控制區(qū)序列的中華金沙鰍和短身金沙鰍單倍型NJ樹(shù)(A)和BI樹(shù)(B)(以犁頭鰍作外類(lèi)群)Fig.3 Haplotype NJ phylogenetic tree (A)and BI phylogenetic tree (B)based on mtDNA control region sequence of J.sinensis and J.abbreviata (L.fimbriata was used as outgroup)

    圖4 中華金沙鰍的BSP分析Fig.4 BSP analysis of J.sinensis population

    3 討論

    3.1 群體遺傳多樣性

    本研究顯示,中華金沙鰍具有很高水平的線粒體遺傳多樣性,與Duan等[16]采用微衛(wèi)星標(biāo)記研究的結(jié)果一致,但該研究?jī)H對(duì)雅礱江里莊31尾樣本進(jìn)行了分析,本研究的采樣點(diǎn)和樣本量更為豐富,包含Duan等[16]研究的江段。長(zhǎng)江上游已開(kāi)展的一些鰍科魚(yú)類(lèi)同樣顯示出高水平的遺傳多樣性,如中華沙鰍(Hd=0.986,Pi=0.003 65)[17]、小眼薄鰍(Hd=0.958,Pi=0.004 20)[18]、紅唇薄鰍(Hd=0.907,Pi=0.003 15)[6]、長(zhǎng)薄鰍(Hd=0.916,Pi=0.004 50)[19]等,與長(zhǎng)江鯉科、鲿科等魚(yú)類(lèi)遺傳多樣性高低不同,高遺傳多樣性可能是長(zhǎng)江上游鰍科魚(yú)類(lèi)線粒體DNA共有特征,更快積累變異。但是短身金沙鰍遺傳多樣性水平卻明顯較低,低水平的遺傳多樣性通常是由于近親繁殖或小種群導(dǎo)致[20]。魚(yú)類(lèi)早期資源調(diào)查發(fā)現(xiàn),金沙江干流、長(zhǎng)江上游干流、雅礱江和岷江等均存在中華金沙鰍和短身金沙鰍產(chǎn)卵場(chǎng),表明它們是同域分布的,但所有產(chǎn)卵場(chǎng)的中華金沙鰍產(chǎn)卵規(guī)模均大于短身金沙鰍[4],說(shuō)明短身金沙鰍的種群規(guī)模比中華金沙鰍小。種群小容易近交,發(fā)生遺傳漂變,導(dǎo)致遺傳多樣性降低。

    3.2 種間及種群分化

    長(zhǎng)江大多數(shù)魚(yú)類(lèi)的種類(lèi)鑒別可以通過(guò)DNA條形碼(mtDNACOI基因)成功解決,但中華金沙鰍和短身金沙鰍例外[21]。這兩種魚(yú)之間Cytb基因的遺傳距離也很低(Cytb:0.008±0.001),系統(tǒng)發(fā)育關(guān)系上出現(xiàn)種間單倍型嵌套,無(wú)法獲得確定的種間區(qū)分[22]。本研究采用線粒體控制區(qū)分析的結(jié)果顯示,兩種金沙鰍的種間單倍型遺傳距離達(dá)到0.035以上,雖然與中華金沙鰍種內(nèi)單倍型之間遺傳距離存在部分重疊,但NJ系統(tǒng)發(fā)育樹(shù)可以將它們分開(kāi)(圖3)。這說(shuō)明線粒體控制區(qū)比COI和Cytb基因?qū)鹕出q種類(lèi)鑒別具有更高的區(qū)分度,但遺傳距離上仍然沒(méi)有明顯種間間隔,也說(shuō)明了這兩種魚(yú)物種分化時(shí)間比較晚。

    AMOVA分析顯示,中華金沙鰍和短身金沙鰍群體均未發(fā)生地理遺傳分化,可能與其產(chǎn)卵特征有關(guān)。這兩種魚(yú)均產(chǎn)漂流性卵,受精卵隨水漂流較長(zhǎng)的距離而發(fā)育,成體向各棲息地?cái)U(kuò)散,促進(jìn)了不同群體間的基因交流。長(zhǎng)江上游一些產(chǎn)漂流性卵魚(yú)類(lèi)的遺傳結(jié)構(gòu)也表現(xiàn)出類(lèi)似的同質(zhì)性,如長(zhǎng)鰭吻鮈[23]和異鰾鰍鮀[24]等。

    盡管未檢測(cè)到地理遺傳分化,但中華金沙鰍單倍型系統(tǒng)發(fā)育樹(shù)(圖2)顯示形成了三個(gè)具有較高支持率的譜系,表明中華金沙鰍發(fā)生群體內(nèi)遺傳分化。中華金沙鰍這一遺傳結(jié)構(gòu)特征符合Avise等[25]提出的線粒體系統(tǒng)地理學(xué)的第四類(lèi)模式,即具有弱的同域單倍型譜系分化,但沒(méi)有顯著地理分化,是由于缺少長(zhǎng)期亞群體隔離和高水平基因流,或群體擴(kuò)張和單倍型廣泛分布引起。長(zhǎng)江上游一些魚(yú)類(lèi)就具有相似的遺傳結(jié)構(gòu),如小眼薄鰍[18]和異鰾鰍鮀[24]等。短身金沙鰍單倍型沒(méi)有出現(xiàn)譜系分化,可能與其單倍型數(shù)少,單倍型間序列差異小有關(guān)。

    3.3 種群歷史

    BSP分析結(jié)果顯示中華金沙鰍群體在10萬(wàn)~17萬(wàn)年前發(fā)生擴(kuò)張,這一時(shí)期剛好是在第四紀(jì)冰期的廬山亞冰期(20萬(wàn)~23萬(wàn)年前)與大理亞冰期(1萬(wàn)~11萬(wàn)年前)之間的間冰期,間冰期相對(duì)于冰期溫暖[26],有利于物種生存,也有利于物種從冰期的“避難所”向外擴(kuò)散,發(fā)生種群擴(kuò)張。研究表明,長(zhǎng)江上游的一些魚(yú)類(lèi)像圓筒吻鮈[27]和銅魚(yú)[28]等,也受到冰期影響,其種群在歷史上也發(fā)生過(guò)瓶頸和擴(kuò)張現(xiàn)象。中華金沙鰍在歷史上經(jīng)歷了快速擴(kuò)張,群體內(nèi)可能積累了大量變異,與其較高的遺傳多樣性結(jié)果相吻合。

    3.4 保護(hù)建議

    本研究發(fā)現(xiàn)兩種金沙鰍均沒(méi)有發(fā)生地理遺傳分化,為單一的進(jìn)化顯著單元,建議就地保護(hù)。目前長(zhǎng)江上游已經(jīng)建成向家壩和溪洛渡等大型水電工程,正在建設(shè)或規(guī)劃建設(shè)更多梯級(jí)水電站,這將使中華金沙鰍和短身金沙鰍等的生境壓縮和碎片化,阻礙群體基因交流。因此,今后還需加強(qiáng)對(duì)長(zhǎng)江上游兩種金沙鰍等魚(yú)類(lèi)產(chǎn)卵場(chǎng)和遺傳多樣性的監(jiān)測(cè),特別是對(duì)群體小、遺傳多樣性較低的短身金沙鰍的監(jiān)測(cè)和保護(hù)。

    猜你喜歡
    控制區(qū)金沙魚(yú)類(lèi)
    基于MFCC和ResNet的魚(yú)類(lèi)行為識(shí)別
    基于OMI的船舶排放控制區(qū)SO2減排效益分析
    當(dāng)遠(yuǎn)方不再遙遠(yuǎn)
    如果……就……
    管好高速建筑控制區(qū)
    奇妙的古代動(dòng)物 泥盆紀(jì)的魚(yú)類(lèi)
    阿什河流域非點(diǎn)源污染優(yōu)先控制區(qū)識(shí)別
    山東日照劃定大氣污染物排放控制區(qū)
    魚(yú)類(lèi)是怎樣保護(hù)自己的
    歡樂(lè)小野人
    国产一区二区激情短视频 | 亚洲av男天堂| 操美女的视频在线观看| 成人国产麻豆网| 性少妇av在线| 99久久精品国产亚洲精品| 日韩中文字幕视频在线看片| 可以免费在线观看a视频的电影网站 | 久久久久久久久免费视频了| 久久久久人妻精品一区果冻| 国产视频首页在线观看| 青青草视频在线视频观看| 国精品久久久久久国模美| 丁香六月欧美| 亚洲精品久久午夜乱码| 国产av国产精品国产| 99久久99久久久精品蜜桃| 婷婷色综合www| 曰老女人黄片| 亚洲熟女毛片儿| 日韩中文字幕欧美一区二区 | 人人妻,人人澡人人爽秒播 | 精品少妇一区二区三区视频日本电影 | 秋霞在线观看毛片| 在线 av 中文字幕| 水蜜桃什么品种好| 日韩一区二区三区影片| 久久精品久久久久久久性| 亚洲国产欧美日韩在线播放| 久久久久久久大尺度免费视频| 女人被躁到高潮嗷嗷叫费观| 亚洲图色成人| 一级爰片在线观看| 热99久久久久精品小说推荐| 亚洲,欧美,日韩| 性高湖久久久久久久久免费观看| 精品国产露脸久久av麻豆| 亚洲精品aⅴ在线观看| 丝袜人妻中文字幕| 欧美日韩福利视频一区二区| 亚洲国产av影院在线观看| 中文精品一卡2卡3卡4更新| www.精华液| 亚洲情色 制服丝袜| 久久久久久久精品精品| 精品少妇内射三级| 在线 av 中文字幕| 日韩中文字幕欧美一区二区 | 日本欧美视频一区| 国产精品久久久人人做人人爽| a级片在线免费高清观看视频| 中文字幕人妻丝袜一区二区 | 亚洲国产欧美网| 一级毛片我不卡| 亚洲欧美中文字幕日韩二区| 黄网站色视频无遮挡免费观看| 久久人人97超碰香蕉20202| 国产亚洲av片在线观看秒播厂| 美女午夜性视频免费| 宅男免费午夜| 久久女婷五月综合色啪小说| 欧美精品一区二区大全| 香蕉国产在线看| 亚洲伊人久久精品综合| 国产成人啪精品午夜网站| 啦啦啦中文免费视频观看日本| 无限看片的www在线观看| 亚洲五月色婷婷综合| 999久久久国产精品视频| 欧美av亚洲av综合av国产av | 欧美中文综合在线视频| 国产一区二区 视频在线| 最近中文字幕高清免费大全6| 午夜日韩欧美国产| 在线免费观看不下载黄p国产| 免费黄色在线免费观看| 91精品国产国语对白视频| 国产精品三级大全| 18禁裸乳无遮挡动漫免费视频| 激情五月婷婷亚洲| 丝袜美足系列| 9热在线视频观看99| 国产精品一区二区精品视频观看| 国产亚洲一区二区精品| 韩国精品一区二区三区| 秋霞伦理黄片| 精品酒店卫生间| 交换朋友夫妻互换小说| 超色免费av| 国产不卡av网站在线观看| 精品亚洲成a人片在线观看| 欧美精品亚洲一区二区| 亚洲精品国产区一区二| videos熟女内射| 亚洲国产精品999| 在线天堂最新版资源| 两性夫妻黄色片| 一级毛片 在线播放| 欧美另类一区| 人人澡人人妻人| 国产免费视频播放在线视频| 亚洲精品国产色婷婷电影| 91老司机精品| 亚洲精品中文字幕在线视频| 久热这里只有精品99| 波野结衣二区三区在线| 一边亲一边摸免费视频| 蜜桃国产av成人99| 国产av码专区亚洲av| 国产亚洲欧美精品永久| 水蜜桃什么品种好| 久久鲁丝午夜福利片| 亚洲国产av新网站| 国产老妇伦熟女老妇高清| 午夜免费观看性视频| 国产黄色免费在线视频| 成人漫画全彩无遮挡| 国产精品久久久久久精品古装| 蜜桃在线观看..| 成人免费观看视频高清| 免费观看av网站的网址| 国产黄色免费在线视频| 婷婷色综合大香蕉| 伦理电影免费视频| 午夜福利视频精品| 黄色视频在线播放观看不卡| 极品少妇高潮喷水抽搐| 久久免费观看电影| 国产成人a∨麻豆精品| 久久久久久久大尺度免费视频| 高清视频免费观看一区二区| 不卡视频在线观看欧美| 男人操女人黄网站| 久久热在线av| 丝袜在线中文字幕| 男女床上黄色一级片免费看| 国产国语露脸激情在线看| 黑人猛操日本美女一级片| 91精品伊人久久大香线蕉| 久久午夜综合久久蜜桃| 国产成人av激情在线播放| 中文字幕av电影在线播放| 看免费av毛片| 综合色丁香网| 欧美日韩国产mv在线观看视频| 亚洲精品国产av成人精品| 在线免费观看不下载黄p国产| 女人被躁到高潮嗷嗷叫费观| 自线自在国产av| 亚洲av成人不卡在线观看播放网 | 最近中文字幕2019免费版| 日韩av不卡免费在线播放| 免费观看a级毛片全部| 久久久久久久久久久久大奶| 免费高清在线观看视频在线观看| 哪个播放器可以免费观看大片| 亚洲一级一片aⅴ在线观看| 十分钟在线观看高清视频www| 日本wwww免费看| 国产av国产精品国产| 女的被弄到高潮叫床怎么办| 欧美在线一区亚洲| 又粗又硬又长又爽又黄的视频| 欧美人与性动交α欧美软件| 午夜激情av网站| 18禁动态无遮挡网站| 一边亲一边摸免费视频| 最黄视频免费看| 国产极品天堂在线| 麻豆av在线久日| 免费观看av网站的网址| 女人高潮潮喷娇喘18禁视频| 欧美 日韩 精品 国产| 性色av一级| 午夜免费观看性视频| 丝袜美腿诱惑在线| 精品久久久久久电影网| 啦啦啦视频在线资源免费观看| 亚洲国产欧美网| 午夜福利视频在线观看免费| 欧美日本中文国产一区发布| 电影成人av| 日本vs欧美在线观看视频| 天美传媒精品一区二区| 欧美成人午夜精品| 欧美黑人精品巨大| 国产精品.久久久| 波多野结衣一区麻豆| 熟女少妇亚洲综合色aaa.| 午夜福利一区二区在线看| 久久精品亚洲av国产电影网| 熟女少妇亚洲综合色aaa.| 一本大道久久a久久精品| 女的被弄到高潮叫床怎么办| 国产成人午夜福利电影在线观看| 老鸭窝网址在线观看| 久久久久精品久久久久真实原创| 成人漫画全彩无遮挡| 国产男女内射视频| 免费日韩欧美在线观看| 精品人妻一区二区三区麻豆| 99精国产麻豆久久婷婷| 精品一区二区免费观看| 日韩欧美精品免费久久| 在线观看一区二区三区激情| 精品酒店卫生间| 色吧在线观看| 黄色怎么调成土黄色| 亚洲五月色婷婷综合| 久热这里只有精品99| 亚洲国产欧美一区二区综合| 国产精品欧美亚洲77777| 国产精品久久久久成人av| 高清在线视频一区二区三区| 一级毛片黄色毛片免费观看视频| 免费女性裸体啪啪无遮挡网站| 国产老妇伦熟女老妇高清| svipshipincom国产片| 国产人伦9x9x在线观看| 国产一区二区 视频在线| 久久久久精品国产欧美久久久 | 狠狠精品人妻久久久久久综合| 熟妇人妻不卡中文字幕| 精品亚洲成a人片在线观看| 日韩 欧美 亚洲 中文字幕| 久久久久久人妻| 国产精品蜜桃在线观看| 国产日韩欧美在线精品| 在线精品无人区一区二区三| 久久精品人人爽人人爽视色| 午夜福利网站1000一区二区三区| 免费在线观看黄色视频的| 久久这里只有精品19| 在线观看国产h片| netflix在线观看网站| 亚洲欧美一区二区三区国产| 国产精品国产三级国产专区5o| a 毛片基地| 亚洲国产成人一精品久久久| 欧美黑人精品巨大| 久久久国产精品麻豆| 国产高清不卡午夜福利| 肉色欧美久久久久久久蜜桃| 纵有疾风起免费观看全集完整版| 国产精品国产三级专区第一集| 成人国产麻豆网| 天堂8中文在线网| 久久久久网色| 日韩av在线免费看完整版不卡| 交换朋友夫妻互换小说| 久久久精品国产亚洲av高清涩受| 亚洲精品av麻豆狂野| 欧美成人精品欧美一级黄| 国产亚洲精品第一综合不卡| 免费观看av网站的网址| 在线观看三级黄色| 午夜久久久在线观看| 怎么达到女性高潮| 啦啦啦观看免费观看视频高清 | 国产精品1区2区在线观看.| 国产一卡二卡三卡精品| 亚洲欧美一区二区三区黑人| 国产激情欧美一区二区| 亚洲国产看品久久| 久久九九热精品免费| 久久久久久久久中文| 日韩成人在线观看一区二区三区| 一区二区三区精品91| 女人爽到高潮嗷嗷叫在线视频| 午夜福利欧美成人| 嫁个100分男人电影在线观看| 成人欧美大片| 自线自在国产av| 69av精品久久久久久| 久久人妻av系列| 少妇粗大呻吟视频| 亚洲aⅴ乱码一区二区在线播放 | 两性夫妻黄色片| 麻豆一二三区av精品| 人人妻,人人澡人人爽秒播| 国产成人精品久久二区二区91| 久久久久精品国产欧美久久久| 亚洲国产欧美日韩在线播放| 国产高清有码在线观看视频 | 久久中文字幕一级| 国产精品 欧美亚洲| 老汉色∧v一级毛片| 99国产综合亚洲精品| 精品久久蜜臀av无| 丝袜在线中文字幕| 成年人黄色毛片网站| 国产亚洲精品一区二区www| 在线观看免费视频网站a站| 欧美乱妇无乱码| 成人精品一区二区免费| 在线观看免费视频网站a站| 国产熟女午夜一区二区三区| 超碰成人久久| 亚洲性夜色夜夜综合| www.精华液| 国产一级毛片七仙女欲春2 | 好男人在线观看高清免费视频 | 日本 欧美在线| 黄色a级毛片大全视频| 纯流量卡能插随身wifi吗| 国产又爽黄色视频| 国产国语露脸激情在线看| 看黄色毛片网站| 丁香欧美五月| 人成视频在线观看免费观看| 99国产精品一区二区三区| 国产不卡一卡二| 一区在线观看完整版| 丝袜美足系列| 国产精品,欧美在线| 久久亚洲真实| 国产激情欧美一区二区| 午夜福利欧美成人| 午夜老司机福利片| 成人国产综合亚洲| 香蕉丝袜av| 久久热在线av| 波多野结衣av一区二区av| 国产麻豆69| 女同久久另类99精品国产91| 91av网站免费观看| 黄色毛片三级朝国网站| 午夜免费观看网址| 精品午夜福利视频在线观看一区| 国产一区在线观看成人免费| 欧美色视频一区免费| 欧美绝顶高潮抽搐喷水| 日韩欧美三级三区| 性欧美人与动物交配| 日本黄色视频三级网站网址| 香蕉丝袜av| 九色国产91popny在线| 国产精品久久久久久亚洲av鲁大| 满18在线观看网站| 窝窝影院91人妻| 日韩成人在线观看一区二区三区| 久久香蕉激情| 亚洲欧美日韩另类电影网站| 亚洲人成电影观看| 又黄又粗又硬又大视频| 久久精品影院6| 国产单亲对白刺激| 97碰自拍视频| 免费看十八禁软件| 亚洲成人免费电影在线观看| 亚洲九九香蕉| 亚洲精品av麻豆狂野| 国产一区在线观看成人免费| 99香蕉大伊视频| 欧美亚洲日本最大视频资源| 狠狠狠狠99中文字幕| 看免费av毛片| 一区二区日韩欧美中文字幕| 精品熟女少妇八av免费久了| 一进一出抽搐gif免费好疼| 人成视频在线观看免费观看| 黑人操中国人逼视频| 国产精品久久久人人做人人爽| 美女 人体艺术 gogo| 天天躁夜夜躁狠狠躁躁| 国产精品久久视频播放| 亚洲专区字幕在线| 久久久久久国产a免费观看| 欧美国产日韩亚洲一区| 亚洲,欧美精品.| 十分钟在线观看高清视频www| 一区二区三区激情视频| 国产高清有码在线观看视频 | 亚洲熟妇熟女久久| 亚洲第一电影网av| 国产成人系列免费观看| www.精华液| 色在线成人网| 一区二区三区精品91| 亚洲va日本ⅴa欧美va伊人久久| АⅤ资源中文在线天堂| 伦理电影免费视频| 国产私拍福利视频在线观看| 亚洲av电影不卡..在线观看| 男女午夜视频在线观看| 男人操女人黄网站| 日韩欧美国产一区二区入口| tocl精华| 一级黄色大片毛片| 免费女性裸体啪啪无遮挡网站| 欧美 亚洲 国产 日韩一| 女性生殖器流出的白浆| 国产片内射在线| 欧美在线黄色| 免费不卡黄色视频| www.自偷自拍.com| 久久草成人影院| 乱人伦中国视频| 亚洲成av人片免费观看| 一级a爱视频在线免费观看| 天天添夜夜摸| 亚洲色图综合在线观看| 在线天堂中文资源库| svipshipincom国产片| 村上凉子中文字幕在线| 国产成人一区二区三区免费视频网站| 久久久久久免费高清国产稀缺| 亚洲黑人精品在线| 又黄又粗又硬又大视频| 叶爱在线成人免费视频播放| 手机成人av网站| 日本免费a在线| 久久国产亚洲av麻豆专区| 欧美黄色片欧美黄色片| 91麻豆av在线| av天堂在线播放| 香蕉国产在线看| 国产视频一区二区在线看| 国产精品av久久久久免费| 亚洲色图综合在线观看| 18禁国产床啪视频网站| 丁香六月欧美| 国产私拍福利视频在线观看| 婷婷丁香在线五月| 国产男靠女视频免费网站| 好男人在线观看高清免费视频 | 国产av在哪里看| 两性夫妻黄色片| 亚洲精品一卡2卡三卡4卡5卡| 性色av乱码一区二区三区2| 手机成人av网站| 亚洲欧美一区二区三区黑人| 国产亚洲精品一区二区www| av有码第一页| 日日爽夜夜爽网站| 欧美激情高清一区二区三区| 国产亚洲精品一区二区www| 亚洲中文av在线| 久久久久久亚洲精品国产蜜桃av| 搡老岳熟女国产| 老司机午夜福利在线观看视频| 亚洲全国av大片| 大香蕉久久成人网| 久久精品91蜜桃| 女警被强在线播放| 18禁观看日本| 亚洲人成电影免费在线| 亚洲欧美日韩另类电影网站| 香蕉丝袜av| 欧美精品亚洲一区二区| 狠狠狠狠99中文字幕| 老司机深夜福利视频在线观看| 19禁男女啪啪无遮挡网站| 韩国av一区二区三区四区| 99国产综合亚洲精品| 一级毛片高清免费大全| 亚洲无线在线观看| 免费在线观看日本一区| 91av网站免费观看| 亚洲国产日韩欧美精品在线观看 | av视频免费观看在线观看| 欧美成狂野欧美在线观看| 一边摸一边抽搐一进一出视频| 久久久久九九精品影院| 亚洲欧美精品综合一区二区三区| 桃色一区二区三区在线观看| 国产精品亚洲美女久久久| 高清黄色对白视频在线免费看| www日本在线高清视频| 日韩欧美免费精品| 亚洲一区二区三区不卡视频| 99国产综合亚洲精品| 看黄色毛片网站| 高清黄色对白视频在线免费看| 欧美乱色亚洲激情| 精品久久久久久久毛片微露脸| 国产精品电影一区二区三区| 国产精品九九99| 99国产精品一区二区三区| 久久中文字幕一级| 一本综合久久免费| 久久人妻福利社区极品人妻图片| 免费观看精品视频网站| 亚洲久久久国产精品| 天堂影院成人在线观看| 国产精品久久久久久精品电影 | 亚洲国产欧美网| 日日干狠狠操夜夜爽| 日韩国内少妇激情av| 午夜福利在线观看吧| 亚洲精品一卡2卡三卡4卡5卡| 99在线视频只有这里精品首页| 国产av在哪里看| 女人被狂操c到高潮| 天天躁夜夜躁狠狠躁躁| 欧美成人免费av一区二区三区| 国产又爽黄色视频| 99久久精品国产亚洲精品| 婷婷六月久久综合丁香| 午夜久久久在线观看| 亚洲久久久国产精品| 成人欧美大片| av视频在线观看入口| 熟妇人妻久久中文字幕3abv| 欧美日韩乱码在线| 精品国产超薄肉色丝袜足j| 18美女黄网站色大片免费观看| 男人舔女人的私密视频| 久久精品91无色码中文字幕| 国语自产精品视频在线第100页| 午夜免费观看网址| 亚洲激情在线av| 中国美女看黄片| 亚洲精品av麻豆狂野| 国产精品久久久久久精品电影 | 亚洲视频免费观看视频| 日本撒尿小便嘘嘘汇集6| 国产亚洲精品久久久久久毛片| 老鸭窝网址在线观看| 男男h啪啪无遮挡| 欧美黑人精品巨大| 久久久精品国产亚洲av高清涩受| 黄色成人免费大全| 露出奶头的视频| 亚洲第一电影网av| 国产主播在线观看一区二区| 性少妇av在线| 国产成人欧美| av在线天堂中文字幕| 搡老熟女国产l中国老女人| 欧美黄色片欧美黄色片| 精品久久蜜臀av无| av欧美777| 麻豆一二三区av精品| 大陆偷拍与自拍| 亚洲三区欧美一区| 在线观看日韩欧美| 麻豆成人av在线观看| 日韩有码中文字幕| 大型黄色视频在线免费观看| www.精华液| 国产精品免费一区二区三区在线| 亚洲av片天天在线观看| 亚洲av第一区精品v没综合| 免费在线观看视频国产中文字幕亚洲| 夜夜爽天天搞| 国产一级毛片七仙女欲春2 | 视频在线观看一区二区三区| 精品电影一区二区在线| 亚洲成a人片在线一区二区| 国产成人系列免费观看| 亚洲成国产人片在线观看| 欧美+亚洲+日韩+国产| 在线av久久热| 日韩欧美一区视频在线观看| 最新美女视频免费是黄的| 91av网站免费观看| 97人妻天天添夜夜摸| 青草久久国产| 嫩草影视91久久| 日韩欧美国产在线观看| 啦啦啦 在线观看视频| 韩国av一区二区三区四区| 极品教师在线免费播放| 国产一级毛片七仙女欲春2 | 伊人久久大香线蕉亚洲五| 操出白浆在线播放| 91成人精品电影| 日韩欧美在线二视频| 欧美黑人欧美精品刺激| 一级毛片女人18水好多| 一级片免费观看大全| 黄网站色视频无遮挡免费观看| 精品欧美国产一区二区三| 国产成人系列免费观看| svipshipincom国产片| 99国产精品免费福利视频| 久久 成人 亚洲| 久久精品成人免费网站| 国产99白浆流出| 美女 人体艺术 gogo| 亚洲国产精品999在线| www.熟女人妻精品国产| 亚洲精品美女久久av网站| 一区二区三区国产精品乱码| 一夜夜www| 天天躁夜夜躁狠狠躁躁| 大码成人一级视频| 亚洲第一青青草原| 欧美 亚洲 国产 日韩一| 这个男人来自地球电影免费观看| 老熟妇仑乱视频hdxx| 99香蕉大伊视频| 国产男靠女视频免费网站| 精品日产1卡2卡| 午夜福利,免费看| 亚洲午夜精品一区,二区,三区| 黄片大片在线免费观看| 午夜精品国产一区二区电影| 日本 欧美在线| 色在线成人网| svipshipincom国产片| АⅤ资源中文在线天堂| 91在线观看av| 久9热在线精品视频| 国产男靠女视频免费网站| 日韩欧美在线二视频| 欧美激情高清一区二区三区| 久久国产精品影院| 国产极品粉嫩免费观看在线| 村上凉子中文字幕在线| 一级毛片精品| 51午夜福利影视在线观看| 日韩中文字幕欧美一区二区| 亚洲精品国产一区二区精华液|