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

    基于細(xì)胞色素B(Cyt b)基因?qū)βN嘴鲌不同群體遺傳多樣性的分析

    2021-09-03 02:48:36張桂寧方弟安薛向平張敏瑩馮曉婷楊雪軍
    中國農(nóng)學(xué)通報(bào) 2021年23期
    關(guān)鍵詞:淀山湖長(zhǎng)蕩湖太湖

    張桂寧,方弟安,,薛向平,張敏瑩,馮曉婷,楊雪軍

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

    0 引言

    翹嘴鲌(Culter alburnus Basilewsky)隸屬于鯉科(Cyprinidae)、鲌亞科(Cultrinae)、鲌屬(Culter),俗稱白魚、翹嘴白等,是鲌亞科中最大的一種,是太湖“三白”之一,與松江鱸(Trachidermus fasciatus)、黃河鯉(Cyprinus carpio)、松花江鮭(Siniperca chuatsi)被譽(yù)為中國四大名魚[1]。翹嘴鲌?jiān)谥袊植紡V泛,幾乎各個(gè)江流湖泊中都有存在。翹嘴鲌多2~3齡性成熟,生殖季節(jié)在5—7月,5月中旬為盛產(chǎn)期,具有溯河產(chǎn)卵特性,卵黏性,受精卵黏附于水生植物或砂石上發(fā)育[1]。翹嘴鲌生長(zhǎng)迅速,耐低氧,抗病力強(qiáng),適合高密度養(yǎng)殖,是長(zhǎng)江三角洲地區(qū)的主要名特淡水養(yǎng)殖品種[2]。鑒于翹嘴鲌的重要經(jīng)濟(jì)價(jià)值和生態(tài)地位,受到眾多學(xué)者的廣泛關(guān)注與研究,目前主要集中在生物學(xué)特征研究[3]、資源現(xiàn)狀[4]、育種[5]胚胎發(fā)育[6]、繁殖[7]、病害[8]和遺傳多樣性[9-11]等方面。

    研究遺傳多樣性對(duì)生物多樣性保護(hù)和生物資源可持續(xù)利用有重要意義,由于線粒體DNA(Mitochondrial DNA,mtDNA)具有母系遺傳、進(jìn)化速率快等特性,因此,常用作研究親緣物種遺傳多樣性的理想分子標(biāo)記[12]。在mtDNA的13個(gè)蛋白編碼基因中,細(xì)胞色素b基因(Cytochrome b,Cyt b)的結(jié)構(gòu)和功能較為清楚,堿基含量的富集反映出Cyt b基因在密碼子使用上的偏倚性[13-14]。Cyt b基因進(jìn)化速率適中,包含著從種內(nèi)到種間的較強(qiáng)進(jìn)化信號(hào),故被廣泛地用來進(jìn)行物種多樣性和系統(tǒng)進(jìn)化等方面的研究[15]。王利華等[16]探討基于CO I和Cyt b DNA條形碼鑒定6種鲌屬魚類(達(dá)氏鲌、海南鲌、尖頭鲌、蒙古鲌、擬尖頭鲌和翹嘴鲌)的可行性,通過計(jì)算其堿基含量、群體間的平均遺傳距離、種內(nèi)遺傳距離和構(gòu)建系統(tǒng)進(jìn)化樹,研究顯示,線粒體CO I和Cyt b基因可以作為DNA條形碼鑒定鲌屬魚類。王偉[17]應(yīng)用擴(kuò)增片段長(zhǎng)度多態(tài)性(AFLP)、內(nèi)部簡(jiǎn)單序列重復(fù)(ISSR)和DNA序列(COII和D-LOOP)分析技術(shù),對(duì)中國不同地理的7個(gè)群體翹嘴鲌遺傳結(jié)構(gòu)和多樣性進(jìn)行了分析,并對(duì)各種標(biāo)記的結(jié)果進(jìn)行了橫向比較,研究顯示AFLP標(biāo)記要略優(yōu)于ISSR標(biāo)記技術(shù),兩種序列均能提供一定的遺傳信息,但D-LOOP序列進(jìn)化速度更快,另外中國野生翹嘴鲌遺傳多樣性比養(yǎng)殖群體豐富,遺傳結(jié)構(gòu)合理。黃小彧[11]針對(duì)長(zhǎng)江水系不同江段的14個(gè)翹嘴鲌群體,應(yīng)用線粒體控制區(qū)研究其遺傳多樣性,發(fā)現(xiàn)不同江段的翹嘴鲌群體遺傳多樣豐富度不同,表現(xiàn)出高單倍型多樣性和低核苷酸多樣性共存的現(xiàn)象,遺傳多樣性較低[11]。

    長(zhǎng)江是中國第一大河,對(duì)其漁業(yè)資源遺傳多樣性的研究極其重要,太湖、淀山湖和長(zhǎng)蕩湖均是位于長(zhǎng)江下游的湖泊,是翹嘴鲌的良好棲息場(chǎng)所,近幾年,這些湖泊都有不同程度的富營養(yǎng)化,翹嘴鲌資源有所衰退[18]。翹嘴鲌是“太湖三白”之一,在淀山湖和長(zhǎng)蕩湖中亦為漁業(yè)資源的優(yōu)勢(shì)物種。近年來,基于線粒體DNA對(duì)長(zhǎng)江下游水域翹嘴鲌遺傳多樣性的研究較少。因此,本研究以長(zhǎng)江、長(zhǎng)蕩湖、太湖和淀山湖4個(gè)水域的翹嘴鲌為研究對(duì)象,通過mtDNA Cyt b基因核苷酸序列的比較分析,探討長(zhǎng)江下游4個(gè)水域翹嘴鲌群體的遺傳多樣性、遺傳結(jié)構(gòu)和群體間的進(jìn)化關(guān)系,以深入了解其遺傳背景和種質(zhì)資源現(xiàn)狀,為翹嘴鲌的野生資源保護(hù)提供基礎(chǔ)數(shù)據(jù)。

    1 材料和方法

    1.1 實(shí)驗(yàn)材料

    2019年5—8月,使用定制絲網(wǎng)和流刺網(wǎng)分別于淀山湖(DSH)、太湖(TH)、長(zhǎng)蕩湖(CD)和長(zhǎng)江(CJ)水域的固定采樣位點(diǎn)隨機(jī)捕獲收集翹嘴鲌尾鰭樣本,浸泡于無水乙醇中保存,其中淀山湖47尾,太湖40尾,長(zhǎng)蕩湖40尾,長(zhǎng)江84尾,共計(jì)211尾。

    1.2 DNA提取

    將鰭條樣品剪取約30 mg的組織樣品于離心管中,使用試劑盒法提取DNA(海洋動(dòng)物組織基因組DNA提取試劑盒,天根生化有限公司),并用濃度1%的瓊脂糖凝膠電泳檢測(cè)DNA的濃度和純度,符合后續(xù)擴(kuò)增要求的DNA樣品置于-20℃保存?zhèn)溆谩?/p>

    1.3 PCR和DNA測(cè)序

    PCR擴(kuò)增Cyt b基因所需引物序列為Cytb-F:5-GACTTGAAAAACCACCGTTG-3 Cytb-R:5-CTCCG ATCTCCGGATTACAAGA-3。

    PCR反應(yīng)體系:總體積25 μL,包括dNTP Mixture 12 μL,上、下游引物各1.25 μL(50 μmol/L),DNA 2 μL,ddH2O 8.5 μL。PCR擴(kuò)增程序:94℃預(yù)變形5 min;94℃變性10 s,55℃退火10 s,72℃延伸30 s,35個(gè)循環(huán);72℃延伸5 min;4℃結(jié)束。2%瓊脂糖凝膠電泳檢測(cè)PCR產(chǎn)物后,將條帶清晰擴(kuò)增有效的PCR產(chǎn)物送亦欣生物科技(上海)有限公司進(jìn)行雙向測(cè)序,測(cè)序引物與擴(kuò)增引物一致。

    1.4 數(shù)據(jù)分析

    采用MEGA 7.0.26軟件中的Align by clustalW軟件進(jìn)行Cyt b基因原始序列的比對(duì)、剪切。用DNASP v5[19]軟件計(jì)算群體單倍型,DNASP v5軟件進(jìn)行Tajima's D和Fu’s Fs的中性測(cè)試,以衡量群體是否發(fā)生了擴(kuò)張。采用Arlequin 3.5軟件[20]進(jìn)行遺傳多樣性、群體的分子方差分析(Analysis of molecular variance,AMOVA)、遺傳距離和基因流(Nm)值的計(jì)算。在MEGA 7.0.26軟件中以Kimura雙參數(shù)法(Kimura-2-parameter)為替代模型,采用NJ法(Neighbour-joining,NJ)構(gòu)建翹嘴鲌不同群體單倍型的分子聚類樹,從分支關(guān)系分析各群體之間的親緣遠(yuǎn)近關(guān)系。

    2 結(jié)果與分析

    2.1 序列特征

    PCR產(chǎn)物基因測(cè)序之后經(jīng)MEGA7.0.26軟件中的Align by Clustal W進(jìn)行原始序列的比對(duì)剪切,得到1088 bp的序列,與NCBI數(shù)據(jù)庫的翹嘴鲌Cyt b序列比對(duì),數(shù)據(jù)同源性達(dá)到99%以上。MEGA7.0.26軟件分析211個(gè)個(gè)體序列的堿基組成平均含量為T(27.1%),C(23.1%),A(27.9%),G(21.9%),A+T含量為55%,G+C為45%。表1中211個(gè)個(gè)體共有46個(gè)單倍型,單倍型多樣性在(0.849±0.004)~(0.946±0.015),4個(gè)群體中太湖群體的最高,其次是長(zhǎng)蕩湖群體,再次是長(zhǎng)江群體,淀山湖水域的翹嘴鲌單倍型最低。核苷酸多樣性在0.0083~0.174,整體而言核苷酸多樣性較高,其中淀山湖群體>長(zhǎng)蕩湖群體>長(zhǎng)江群體>太湖群體。MEGA7.0.26軟件中分析得到單態(tài)突變位點(diǎn)46個(gè),簡(jiǎn)約信息點(diǎn)1011個(gè),變異位點(diǎn)1057,變異率為97.15%,保守位點(diǎn)31個(gè)。DNASP軟件分析得到258個(gè)多態(tài)性位點(diǎn),約占核苷酸總數(shù)的23.71%。轉(zhuǎn)換與顛換比值(si/sv)為0.7。

    表1 不同群體的單倍型及遺傳多樣性參數(shù)

    2.2 單倍型分布

    在211個(gè)翹嘴鲌個(gè)體共有46個(gè)單倍型,其中Hap1、Hap2、Hap3、Hap4和Hap8在長(zhǎng)江和長(zhǎng)蕩湖群體中共享,其余41個(gè)單倍型分別在4個(gè)地理群體中獨(dú)有,占單倍型的89.13%,其中Hap5、Hap6和Hap10-Hap19在長(zhǎng)江群體中特有,Hap20-Hap34、Hap45及Hap46在淀山湖群體中獨(dú)有,Hap7、Hap9和Hap40在長(zhǎng)蕩湖中特有,Hap35-Hap39和Hap41-Hap44在太湖群體中獨(dú)有。

    通過Network軟件繪制單倍型網(wǎng)絡(luò)圖,如圖1所示,每個(gè)圓代表一個(gè)單倍型,其大小與總頻率成正比,圓圈大小表示該單倍型出現(xiàn)的次數(shù)。每個(gè)單倍型按地理分布用不同的顏色表示,4個(gè)地理群體的單倍型構(gòu)成一個(gè)網(wǎng)狀結(jié)構(gòu)圖,只有長(zhǎng)江群體的長(zhǎng)蕩湖群體有共享單倍型,長(zhǎng)江群體和淀山湖群體的單倍型散布在各個(gè)位置,其余群體的單倍型分布較均勻。圖中圓的面積代表出現(xiàn)次數(shù),由圖可知,Hap2出現(xiàn)次數(shù)最多為54次,包含長(zhǎng)江群體和長(zhǎng)蕩湖群體,其次是Hap41,出現(xiàn)23次,屬于太湖群體。Hap14出現(xiàn)22次,屬于長(zhǎng)江群體。

    圖1 4個(gè)地理群體翹嘴鲌的單倍型分布網(wǎng)絡(luò)圖

    2.3 遺傳分化分析

    由MEGA7.0.26軟件中的Data-distance計(jì)算出的4個(gè)地理群體的遺傳距離如表2所示,在表中可看出種內(nèi)遺傳距離的范圍是0.01~0.04,其中太湖群體的種內(nèi)遺傳距離是0.01,長(zhǎng)蕩湖群體的種內(nèi)遺傳距離最大,達(dá)到0.04。4個(gè)群體的種間遺傳距離0.052~0.152,太湖群體和淀山湖群體的種間遺傳距離最遠(yuǎn)是0.152。

    表2 不同翹嘴鲌群體遺傳距離

    4個(gè)地理的翹嘴鲌群體基因流和遺傳分化指數(shù)如表3和4所示,從表中顯示出兩兩群體間的遺傳分化指數(shù)范圍在0.197~0.646?;蛄髦档姆秶鸀?.152~1.010,其中長(zhǎng)江和淀山湖群體的基因流值最小為0.152,長(zhǎng)江與長(zhǎng)蕩湖群體的基因流值最大為1.019。

    表3 不同翹嘴鲌群體基因流(Nm)及遺傳分化指數(shù)(Fst)

    表4 4個(gè)翹嘴鲌地理群體間遺傳分化參數(shù)

    2.4 聚類分析

    應(yīng)用MEGA7.0.26軟件對(duì)211個(gè)翹嘴鲌個(gè)體的Cyt b的序列進(jìn)行聚類分析,得出4個(gè)翹嘴鲌群體46個(gè)單倍型的系統(tǒng)進(jìn)化樹如圖2所示。不同地理來源的單倍型總共分成7個(gè)大分支,Hap25、Hap28、Hap31和Hap33獨(dú)自各成一支,屬于太湖群體。Hap30和Hap34組成一支,亦屬于太湖群體,Hap45獨(dú)自為一支,屬于淀山湖群體,剩余一支包含很多單倍型,其中不同地理群體的單倍型散亂分布。

    圖2 4個(gè)地理群體翹嘴鲌的單倍型NJ系統(tǒng)發(fā)育樹

    2.5 4個(gè)翹嘴鲌群體錯(cuò)配分布圖

    用DNASPv5軟件分析得到如圖3所示的錯(cuò)配分布圖,由圖3可知,當(dāng)把4個(gè)水域的翹嘴鲌群體作為一個(gè)整體分析時(shí),其歧點(diǎn)分布圖大體呈現(xiàn)為單峰。

    圖3 錯(cuò)配分布曲線

    3 討論

    3.1 序列組成特點(diǎn)

    本研究中211個(gè)個(gè)體序列的堿基組成平均含量為T(27.1%),C(23.1%),A(27.9%)和 G(21.9%),A+T 的含量占55%,G+C占45%,表現(xiàn)出A+T的含量高于G+C的含量。與張久盤[21]文中A+T含量為56.8%,G+C含量為43.2%;馮曉宇[22]文中的A+T含量為56.6%,G+C含量為43.4%,研究結(jié)果相一致,均表現(xiàn)出堿基偏倚性,符合大多數(shù)硬骨魚的特征[22]。

    3.2 翹嘴鲌地理群體遺傳多樣性

    遺傳多樣性是生物進(jìn)化和適應(yīng)環(huán)境的基礎(chǔ),黃小彧[11]等對(duì)線粒體DNA序列遺傳變異的研究表明,若單倍型多樣度(Hd)高于0.5且核苷酸多樣度(Pi)低于0.005,是受瓶頸效應(yīng)后群體快速擴(kuò)張導(dǎo)致,表明該群體的遺傳多樣性水平低;如果單倍型多樣度(Hd)高于0.5且核苷酸多樣度(Pi)高于0.005,則該群體遺傳多樣性水平較高[23]。本研究中整個(gè)翹嘴鲌群體的遺傳多樣性水平較高,遺傳多樣性豐富,與徐丹丹[24](Hd:0.948,Pi:0.017),徐宇[25](Hd:0.969,Pi:0.2008)和張久盤[21](Hd:0.998,Pi:0.039)的研究結(jié)果一致。本研究中長(zhǎng)江群體表現(xiàn)出高的單倍型多樣性和高的核苷酸多樣性,與黃小彧等[11]的結(jié)果(Hd:0.866,Pi:0.00330)存在差異,可能是黃小彧等研究中使用的是線粒體控制區(qū),控制區(qū)是線粒體DNA的高突變區(qū);還可能是近幾年長(zhǎng)江的水域魚類生境優(yōu)化和資源恢復(fù)措施導(dǎo)致長(zhǎng)江翹嘴鲌的遺傳多樣性較高。

    遺傳多樣性的高低通常由遺傳距離的大小來反映,種內(nèi)個(gè)體間遺傳距離小,彼此間親緣關(guān)系近,說明遺傳多樣性低,反之,遺傳距離大則遺傳多樣性高[17]。本研究中4個(gè)地理群體的種內(nèi)遺傳距離的范圍是0.01~0.04,其中長(zhǎng)江群體和長(zhǎng)蕩湖群體的種內(nèi)遺傳距離相同且最大,遺傳多樣性較高,太湖群體的種內(nèi)遺傳距離最小,遺傳多樣性也較低,與王偉等[17]的研究結(jié)果一致。通常情況下,地理距離越遠(yuǎn),群體間遺傳距離也就越大,而徐丹丹等[24]在基于微衛(wèi)星標(biāo)記和線粒體Cyt b基因序列針對(duì)鲇(Silurus asotus)遺傳多樣性的研究一文中,發(fā)現(xiàn)不同鲇群體遺傳距離與其地理距離之間沒有相關(guān)性;沙航等[26]基于COI序列分析長(zhǎng)江中上游6個(gè)鰱(Hypophthalmichthys molitrix)地理群體遺傳多樣性也斧正了遺傳距離的大小與群體地理距離沒有明顯相關(guān)性。本研究中淀山湖距離長(zhǎng)蕩湖比距離太湖遠(yuǎn),但淀山湖群體與長(zhǎng)蕩湖群體的遺傳距離小于淀山湖群體與太湖群體的遺傳距離,這3個(gè)湖泊的翹嘴鲌群體之間遺傳距離與其分布的地理距離亦沒有相關(guān)性,這一結(jié)果與上述研究結(jié)果類似。

    群體遺傳學(xué)中通常以Fst值的大小作為衡量群體間遺傳分化程度的常見指標(biāo),F(xiàn)st<0.05表示群體間是低度分化;0.05<Fst<0.15 表示群體間是中度分化;0.15<Fst<0.25表示群體間是高度分化;0.25<Fst<1表示群體間有極大的遺傳分化[27]。本研究中長(zhǎng)江翹嘴鲌群體和長(zhǎng)蕩湖群體的Fst相對(duì)較小,為0.197,屬于高度分化,長(zhǎng)江群體、淀山湖群體和太湖群體間的Fst值較大,屬于極大分化,這有可能是研究水域存在地理隔離,翹嘴鲌群體間不存在基因交流。群體間平均凈遺傳距離(Da)是衡量群體間分化程度的重要指標(biāo),值越大,群體間的遺傳分化程度越大[24,28]。長(zhǎng)蕩湖和太湖群體的平均凈遺傳距離相對(duì)其它群體之間較大(0.48),遺傳分化程度高,與其遺傳分化指數(shù)較高的研究結(jié)果相一致。本研究中4個(gè)群體之間的基因流值都小于4,研究結(jié)果進(jìn)一步說明了4個(gè)翹嘴鲌群體之間基因交流較少,存在一定程度的遺傳分化現(xiàn)象,這與陳會(huì)娟等[29]的結(jié)果一致。

    3.3 翹嘴鲌群體歷史動(dòng)態(tài)

    一般要判斷某一群體有沒有發(fā)生過群體擴(kuò)張,有兩種方法:中性檢測(cè)和歧點(diǎn)分布分析[11,30-31]。本研究用這兩種方法來分析了翹嘴鲌群體的種群歷史動(dòng)態(tài)情況。當(dāng)檢測(cè)數(shù)值為負(fù)值,且具有顯著意義時(shí),則認(rèn)為該群體曾發(fā)生過群體擴(kuò)張現(xiàn)象,如果群體發(fā)生過擴(kuò)張,則歧點(diǎn)分布圖呈現(xiàn)出單峰;如果群體未經(jīng)歷過擴(kuò)張,呈現(xiàn)穩(wěn)定狀態(tài),則其歧點(diǎn)分布圖呈現(xiàn)為多峰[11]。對(duì)211個(gè)翹嘴鲌個(gè)體的Cyt b中性檢驗(yàn)得出Tajima's D值,除了長(zhǎng)蕩湖群體,其余群體都為負(fù)值,且長(zhǎng)江和太湖群體的Tajima's D表現(xiàn)出極顯著性差異P<0.01,表明群體在歷史上曾發(fā)生過群體擴(kuò)張。淀山湖群體的Tajima's D為負(fù)值,但不具有顯著性差異。除此之外,歧點(diǎn)分布圖顯示,4個(gè)翹嘴鲌地理群體的錯(cuò)配分布曲線中,將所有翹嘴鲌群體作為一個(gè)整體分析,錯(cuò)配曲線呈現(xiàn)單峰,說明該翹嘴鲌群體發(fā)生過群體擴(kuò)張,這與張力文等[23]的結(jié)果相一致。

    4 結(jié)論

    本研究對(duì)4個(gè)水域的翹嘴鲌進(jìn)行遺傳多樣性和遺傳結(jié)構(gòu)的分析,采樣點(diǎn)覆蓋長(zhǎng)江下游典型水域,樣本較多。結(jié)果顯示4個(gè)地理水域的翹嘴鲌群體,表現(xiàn)出高的單倍型和高的核苷酸多樣性,遺傳多樣性較為豐富,群體間都存在遺傳分化,分化程度不同。研究表明,線粒體DNA Cyt b可以作為評(píng)估翹嘴鲌的遺傳多樣性和種質(zhì)資源現(xiàn)狀的有效分子標(biāo)記,研究結(jié)果為翹嘴鲌種質(zhì)資源評(píng)估和保護(hù)提供了科學(xué)的數(shù)據(jù)支撐。后續(xù)研究將繼續(xù)增加野生翹嘴鲌群體的樣品數(shù)量和研究水域,采用多種分子標(biāo)記方法,全面而系統(tǒng)的評(píng)估長(zhǎng)江下游各個(gè)水域翹嘴鲌的種質(zhì)資源狀況和分析其遺傳結(jié)構(gòu)。

    猜你喜歡
    淀山湖長(zhǎng)蕩湖太湖
    豐水期長(zhǎng)蕩湖水質(zhì)評(píng)價(jià)與微生物群落結(jié)構(gòu)特征研究
    『長(zhǎng)蕩湖』八鮮宴
    美食(2017年9期)2017-09-15 08:52:47
    太湖思變2017
    玩具世界(2017年4期)2017-07-21 13:27:24
    太湖攬春
    寶藏(2017年2期)2017-03-20 13:16:42
    太湖
    中亞信息(2016年3期)2016-12-01 06:08:24
    蘆葦在淀山湖是站不穩(wěn)的〔外一首〕
    中國詩歌(2016年6期)2016-11-25 23:54:04
    太湖一角
    長(zhǎng)蕩湖河蟹產(chǎn)業(yè)現(xiàn)狀及發(fā)展對(duì)策
    昆山市淀山湖鎮(zhèn)河道水體生態(tài)修復(fù)分析研究
    淀山湖夏秋季浮游植物優(yōu)勢(shì)種生態(tài)特征分析
    亚洲国产成人一精品久久久| 少妇的丰满在线观看| 国产精品1区2区在线观看. | 激情在线观看视频在线高清 | 日韩免费av在线播放| av国产精品久久久久影院| 午夜成年电影在线免费观看| 欧美日韩一级在线毛片| 亚洲av日韩在线播放| 亚洲人成77777在线视频| 亚洲男人天堂网一区| 天天影视国产精品| 动漫黄色视频在线观看| 久久国产精品人妻蜜桃| 国产极品粉嫩免费观看在线| 亚洲av片天天在线观看| 国产精品久久久av美女十八| 日日爽夜夜爽网站| 午夜久久久在线观看| 精品免费久久久久久久清纯 | 在线观看66精品国产| 最近最新中文字幕大全电影3 | 久久免费观看电影| 高清视频免费观看一区二区| 亚洲全国av大片| 脱女人内裤的视频| 国产精品久久电影中文字幕 | 两个人看的免费小视频| 最新的欧美精品一区二区| 精品乱码久久久久久99久播| 大型黄色视频在线免费观看| 在线观看舔阴道视频| 国产亚洲欧美精品永久| 亚洲av第一区精品v没综合| 久热这里只有精品99| 国产精品亚洲一级av第二区| 日日夜夜操网爽| 最新在线观看一区二区三区| 啪啪无遮挡十八禁网站| 国产精品自产拍在线观看55亚洲 | 热99国产精品久久久久久7| 亚洲国产av影院在线观看| 国产成人精品久久二区二区91| 国产一区二区激情短视频| 久久精品国产a三级三级三级| 91九色精品人成在线观看| 香蕉国产在线看| 叶爱在线成人免费视频播放| avwww免费| 妹子高潮喷水视频| 亚洲va日本ⅴa欧美va伊人久久| 啦啦啦中文免费视频观看日本| 久久精品国产99精品国产亚洲性色 | 男人舔女人的私密视频| 精品福利观看| 一区二区三区国产精品乱码| 十分钟在线观看高清视频www| 欧美精品人与动牲交sv欧美| 天堂动漫精品| 日本五十路高清| 午夜激情av网站| 99精品在免费线老司机午夜| 国产单亲对白刺激| 桃红色精品国产亚洲av| 精品久久久久久久毛片微露脸| 亚洲国产看品久久| 美女视频免费永久观看网站| 成人永久免费在线观看视频 | 天天躁狠狠躁夜夜躁狠狠躁| 动漫黄色视频在线观看| 手机成人av网站| 亚洲精品国产色婷婷电影| 热99re8久久精品国产| 免费女性裸体啪啪无遮挡网站| 亚洲国产看品久久| 婷婷成人精品国产| 人妻久久中文字幕网| 成人黄色视频免费在线看| 精品人妻熟女毛片av久久网站| 精品人妻1区二区| 国产高清视频在线播放一区| 精品国产国语对白av| 两人在一起打扑克的视频| 日本精品一区二区三区蜜桃| 精品亚洲乱码少妇综合久久| 国产精品免费视频内射| 欧美 亚洲 国产 日韩一| 欧美中文综合在线视频| 国产麻豆69| 日韩熟女老妇一区二区性免费视频| 麻豆成人av在线观看| 亚洲av成人一区二区三| 欧美精品av麻豆av| 高潮久久久久久久久久久不卡| 美女福利国产在线| 午夜激情久久久久久久| 如日韩欧美国产精品一区二区三区| 亚洲九九香蕉| 亚洲,欧美精品.| 日韩成人在线观看一区二区三区| 男女午夜视频在线观看| 老司机深夜福利视频在线观看| 最近最新中文字幕大全电影3 | 自拍欧美九色日韩亚洲蝌蚪91| 久久精品91无色码中文字幕| 日韩视频一区二区在线观看| 俄罗斯特黄特色一大片| 欧美成人午夜精品| 免费观看a级毛片全部| 久久精品亚洲av国产电影网| 一级,二级,三级黄色视频| 亚洲精华国产精华精| 中文亚洲av片在线观看爽 | 欧美日韩福利视频一区二区| 国产激情久久老熟女| 国产精品免费大片| 啦啦啦中文免费视频观看日本| 一边摸一边抽搐一进一出视频| 男人舔女人的私密视频| 午夜福利,免费看| 久久久国产精品麻豆| 18在线观看网站| 国产麻豆69| 国产99久久九九免费精品| 无限看片的www在线观看| 欧美精品啪啪一区二区三区| 亚洲自偷自拍图片 自拍| 久久天躁狠狠躁夜夜2o2o| 波多野结衣av一区二区av| 国产精品影院久久| 变态另类成人亚洲欧美熟女 | 三级毛片av免费| 亚洲一区二区三区欧美精品| 女人爽到高潮嗷嗷叫在线视频| 十八禁高潮呻吟视频| 亚洲精品自拍成人| 男女之事视频高清在线观看| 国产日韩欧美在线精品| 在线观看人妻少妇| 男女午夜视频在线观看| 女人被躁到高潮嗷嗷叫费观| 麻豆乱淫一区二区| 国产精品久久久久久精品电影小说| 精品卡一卡二卡四卡免费| 国精品久久久久久国模美| 日日摸夜夜添夜夜添小说| 正在播放国产对白刺激| 精品国产乱码久久久久久小说| 高清毛片免费观看视频网站 | 露出奶头的视频| 久久精品熟女亚洲av麻豆精品| 国产不卡一卡二| 午夜福利欧美成人| 精品亚洲成国产av| 大型黄色视频在线免费观看| 婷婷成人精品国产| 一边摸一边做爽爽视频免费| 欧美日韩精品网址| www日本在线高清视频| 又黄又粗又硬又大视频| 亚洲国产欧美网| 亚洲 欧美一区二区三区| 最新美女视频免费是黄的| 熟女少妇亚洲综合色aaa.| 亚洲精品乱久久久久久| 性少妇av在线| 免费在线观看黄色视频的| 咕卡用的链子| 国产区一区二久久| 啪啪无遮挡十八禁网站| 狠狠婷婷综合久久久久久88av| 精品国产乱码久久久久久男人| 91av网站免费观看| 91成年电影在线观看| 巨乳人妻的诱惑在线观看| 欧美在线黄色| 午夜福利欧美成人| 两个人免费观看高清视频| 精品一品国产午夜福利视频| 久久狼人影院| 激情视频va一区二区三区| 水蜜桃什么品种好| 国产精品免费视频内射| 建设人人有责人人尽责人人享有的| 黑人巨大精品欧美一区二区mp4| 97人妻天天添夜夜摸| 成人永久免费在线观看视频 | 亚洲精品一卡2卡三卡4卡5卡| 中文字幕高清在线视频| 一本综合久久免费| 精品久久久精品久久久| 成年人午夜在线观看视频| 欧美性长视频在线观看| 我要看黄色一级片免费的| 亚洲va日本ⅴa欧美va伊人久久| 国产视频一区二区在线看| 亚洲熟女毛片儿| 国产亚洲av高清不卡| 麻豆乱淫一区二区| 精品国产亚洲在线| 成人国产一区最新在线观看| 伊人久久大香线蕉亚洲五| 亚洲欧美色中文字幕在线| 国产真人三级小视频在线观看| 日韩欧美一区二区三区在线观看 | 久久精品国产99精品国产亚洲性色 | 在线观看舔阴道视频| 美女国产高潮福利片在线看| 亚洲成人免费电影在线观看| 极品少妇高潮喷水抽搐| 美女午夜性视频免费| 在线永久观看黄色视频| 欧美黄色片欧美黄色片| 19禁男女啪啪无遮挡网站| 叶爱在线成人免费视频播放| 色精品久久人妻99蜜桃| 视频区图区小说| 国产麻豆69| 亚洲免费av在线视频| 少妇的丰满在线观看| 久久久国产一区二区| 色尼玛亚洲综合影院| 怎么达到女性高潮| 亚洲 欧美一区二区三区| 午夜福利影视在线免费观看| 女同久久另类99精品国产91| 国产亚洲一区二区精品| 午夜视频精品福利| 精品亚洲成a人片在线观看| 日韩精品免费视频一区二区三区| 日韩视频一区二区在线观看| av福利片在线| 天堂8中文在线网| 欧美人与性动交α欧美软件| 国产日韩欧美视频二区| 午夜视频精品福利| 久久久欧美国产精品| 黑丝袜美女国产一区| a在线观看视频网站| 操美女的视频在线观看| 国产成人av激情在线播放| 欧美日韩av久久| 亚洲美女黄片视频| 99精品欧美一区二区三区四区| 日本a在线网址| 女警被强在线播放| 精品第一国产精品| 色老头精品视频在线观看| 日韩制服丝袜自拍偷拍| 国产精品影院久久| 亚洲第一av免费看| 久久久欧美国产精品| 捣出白浆h1v1| 欧美精品人与动牲交sv欧美| 汤姆久久久久久久影院中文字幕| 亚洲一卡2卡3卡4卡5卡精品中文| 精品欧美一区二区三区在线| 最近最新中文字幕大全电影3 | 69精品国产乱码久久久| 极品人妻少妇av视频| 亚洲精品国产一区二区精华液| 夫妻午夜视频| 久久天堂一区二区三区四区| 国产午夜精品久久久久久| av有码第一页| 高清黄色对白视频在线免费看| 亚洲伊人久久精品综合| 久久久久精品国产欧美久久久| 精品高清国产在线一区| 欧美国产精品va在线观看不卡| 一二三四在线观看免费中文在| 老熟妇仑乱视频hdxx| 国产日韩欧美在线精品| 91麻豆精品激情在线观看国产 | 视频区图区小说| 亚洲欧洲日产国产| 国产99久久九九免费精品| 成人永久免费在线观看视频 | 精品久久久久久久毛片微露脸| 国产精品久久久久久人妻精品电影 | 大香蕉久久成人网| 欧美日韩成人在线一区二区| 国产麻豆69| 成人18禁在线播放| 国产1区2区3区精品| 一边摸一边抽搐一进一出视频| 热re99久久精品国产66热6| 亚洲少妇的诱惑av| 久久久欧美国产精品| 人人妻人人澡人人爽人人夜夜| 久久免费观看电影| 国产精品秋霞免费鲁丝片| 一级毛片电影观看| 脱女人内裤的视频| 五月天丁香电影| 国产一区二区三区视频了| 巨乳人妻的诱惑在线观看| 十八禁网站网址无遮挡| 亚洲精品自拍成人| 国产片内射在线| 久久天躁狠狠躁夜夜2o2o| 老司机深夜福利视频在线观看| 亚洲色图av天堂| 自线自在国产av| av网站免费在线观看视频| 国产精品一区二区免费欧美| 男人舔女人的私密视频| 国产一区二区三区在线臀色熟女 | 法律面前人人平等表现在哪些方面| 免费人妻精品一区二区三区视频| 免费在线观看影片大全网站| 国产成人精品久久二区二区免费| 一本综合久久免费| 欧美 亚洲 国产 日韩一| 啦啦啦中文免费视频观看日本| 精品一区二区三区视频在线观看免费 | 最近最新中文字幕大全电影3 | 一本久久精品| 亚洲va日本ⅴa欧美va伊人久久| 大香蕉久久成人网| 亚洲人成伊人成综合网2020| www.精华液| 成年版毛片免费区| 亚洲精品美女久久av网站| 精品熟女少妇八av免费久了| 亚洲免费av在线视频| 我要看黄色一级片免费的| 精品国产亚洲在线| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲第一青青草原| 午夜91福利影院| 亚洲自偷自拍图片 自拍| 精品国产一区二区三区四区第35| 一边摸一边抽搐一进一小说 | 99热国产这里只有精品6| 两人在一起打扑克的视频| 欧美乱码精品一区二区三区| 国产成人欧美| 不卡av一区二区三区| 一区二区av电影网| 首页视频小说图片口味搜索| 高清在线国产一区| 亚洲va日本ⅴa欧美va伊人久久| 国产野战对白在线观看| 久久中文字幕一级| 久久人人爽av亚洲精品天堂| 亚洲欧美精品综合一区二区三区| 国产免费av片在线观看野外av| 欧美亚洲 丝袜 人妻 在线| 国产97色在线日韩免费| 亚洲精品av麻豆狂野| 亚洲一区二区三区欧美精品| 精品国产一区二区三区久久久樱花| 国产成人av教育| 另类亚洲欧美激情| 视频区欧美日本亚洲| 久久久久国内视频| 欧美日韩亚洲国产一区二区在线观看 | 亚洲国产欧美网| 日韩一区二区三区影片| 亚洲男人天堂网一区| 欧美激情极品国产一区二区三区| 国产无遮挡羞羞视频在线观看| 涩涩av久久男人的天堂| 国产福利在线免费观看视频| 日日夜夜操网爽| 亚洲精品在线观看二区| 亚洲中文av在线| 99热网站在线观看| 精品少妇久久久久久888优播| 在线观看舔阴道视频| 亚洲伊人久久精品综合| 亚洲avbb在线观看| 精品福利观看| 精品一区二区三卡| 99国产综合亚洲精品| 成人国产av品久久久| 欧美久久黑人一区二区| 变态另类成人亚洲欧美熟女 | 亚洲国产欧美日韩在线播放| 日韩精品免费视频一区二区三区| 久久久久久久久久久久大奶| 成人影院久久| 黄色 视频免费看| 亚洲精品粉嫩美女一区| 青草久久国产| 中文字幕av电影在线播放| 亚洲专区中文字幕在线| 久久精品人人爽人人爽视色| 欧美日韩福利视频一区二区| 精品国产国语对白av| 18禁观看日本| 亚洲av日韩精品久久久久久密| 69精品国产乱码久久久| 一级毛片精品| 亚洲 国产 在线| 少妇裸体淫交视频免费看高清 | 18禁观看日本| 国产精品免费一区二区三区在线 | 精品国产超薄肉色丝袜足j| 国产欧美日韩一区二区三| 满18在线观看网站| 老司机午夜福利在线观看视频 | 久久亚洲精品不卡| 久久亚洲真实| 天堂8中文在线网| 久久久久久久国产电影| 欧美日韩亚洲综合一区二区三区_| 菩萨蛮人人尽说江南好唐韦庄| 高清欧美精品videossex| 国产不卡一卡二| 欧美精品人与动牲交sv欧美| 午夜激情av网站| 少妇粗大呻吟视频| 国产午夜精品久久久久久| 午夜精品久久久久久毛片777| 黄色a级毛片大全视频| 成在线人永久免费视频| 性少妇av在线| 激情在线观看视频在线高清 | 国产一区有黄有色的免费视频| 中国美女看黄片| 国产1区2区3区精品| 精品一区二区三卡| 后天国语完整版免费观看| 夜夜骑夜夜射夜夜干| 国产av一区二区精品久久| 国产欧美日韩综合在线一区二区| 欧美人与性动交α欧美软件| 日本黄色日本黄色录像| 大片电影免费在线观看免费| 亚洲全国av大片| 久久精品成人免费网站| 国产成人免费无遮挡视频| 美女午夜性视频免费| 久久午夜亚洲精品久久| 免费在线观看视频国产中文字幕亚洲| a级片在线免费高清观看视频| 热99re8久久精品国产| av免费在线观看网站| 美女高潮到喷水免费观看| 欧美另类亚洲清纯唯美| 久久久久久人人人人人| 大型av网站在线播放| 免费在线观看日本一区| 巨乳人妻的诱惑在线观看| 成人黄色视频免费在线看| 国产精品熟女久久久久浪| 一二三四社区在线视频社区8| 日韩中文字幕视频在线看片| 国产无遮挡羞羞视频在线观看| 少妇裸体淫交视频免费看高清 | www.精华液| 自线自在国产av| 久久免费观看电影| 国产成人精品久久二区二区免费| 久热这里只有精品99| 精品一区二区三区四区五区乱码| 精品少妇一区二区三区视频日本电影| 19禁男女啪啪无遮挡网站| 国产高清国产精品国产三级| 伊人久久大香线蕉亚洲五| 男女边摸边吃奶| 在线观看舔阴道视频| 欧美黑人欧美精品刺激| 婷婷丁香在线五月| 久久久久久久大尺度免费视频| 亚洲国产欧美一区二区综合| 亚洲av片天天在线观看| 欧美老熟妇乱子伦牲交| 欧美黄色淫秽网站| 黄色 视频免费看| 电影成人av| 欧美乱码精品一区二区三区| 人人妻人人爽人人添夜夜欢视频| 亚洲第一欧美日韩一区二区三区 | 国产一卡二卡三卡精品| 成人国产av品久久久| 欧美日韩中文字幕国产精品一区二区三区 | av欧美777| 国产麻豆69| 妹子高潮喷水视频| 亚洲av欧美aⅴ国产| 精品久久久精品久久久| 精品人妻熟女毛片av久久网站| 法律面前人人平等表现在哪些方面| 黑人巨大精品欧美一区二区蜜桃| 夜夜爽天天搞| 亚洲午夜理论影院| 精品少妇黑人巨大在线播放| 午夜激情久久久久久久| 欧美+亚洲+日韩+国产| 大香蕉久久成人网| 激情在线观看视频在线高清 | 人人妻人人澡人人爽人人夜夜| 国产一区二区 视频在线| 午夜福利欧美成人| 极品人妻少妇av视频| 每晚都被弄得嗷嗷叫到高潮| 亚洲av国产av综合av卡| 亚洲自偷自拍图片 自拍| 美女高潮到喷水免费观看| 正在播放国产对白刺激| 日韩视频在线欧美| 夜夜夜夜夜久久久久| 三上悠亚av全集在线观看| 高清毛片免费观看视频网站 | 18禁国产床啪视频网站| 成人三级做爰电影| 女人高潮潮喷娇喘18禁视频| 久久久国产成人免费| 亚洲欧美日韩另类电影网站| 狠狠婷婷综合久久久久久88av| 国产免费视频播放在线视频| 久久av网站| av一本久久久久| 99re在线观看精品视频| 欧美黄色淫秽网站| 乱人伦中国视频| 久久香蕉激情| 热re99久久精品国产66热6| a级毛片在线看网站| 免费日韩欧美在线观看| 叶爱在线成人免费视频播放| 脱女人内裤的视频| 亚洲男人天堂网一区| 久久久久精品国产欧美久久久| 99re在线观看精品视频| 最新在线观看一区二区三区| 国产成人欧美| 欧美精品一区二区免费开放| 亚洲欧美日韩另类电影网站| 电影成人av| 一边摸一边做爽爽视频免费| 桃红色精品国产亚洲av| 另类精品久久| 欧美人与性动交α欧美精品济南到| 天堂俺去俺来也www色官网| 久久亚洲真实| 男男h啪啪无遮挡| 另类亚洲欧美激情| 菩萨蛮人人尽说江南好唐韦庄| 久热爱精品视频在线9| 老司机影院毛片| 制服人妻中文乱码| 波多野结衣av一区二区av| 如日韩欧美国产精品一区二区三区| 免费在线观看视频国产中文字幕亚洲| 亚洲成av片中文字幕在线观看| av又黄又爽大尺度在线免费看| 久久婷婷成人综合色麻豆| 亚洲av日韩在线播放| 99精品在免费线老司机午夜| 亚洲精品久久成人aⅴ小说| 国产国语露脸激情在线看| 少妇 在线观看| 日韩欧美免费精品| 大片免费播放器 马上看| 老汉色av国产亚洲站长工具| 久久ye,这里只有精品| 久久狼人影院| 狂野欧美激情性xxxx| 国产在线一区二区三区精| 中文字幕人妻熟女乱码| 动漫黄色视频在线观看| 男女边摸边吃奶| 老熟妇乱子伦视频在线观看| 国产亚洲精品久久久久5区| 大片免费播放器 马上看| 国产午夜精品久久久久久| 天堂8中文在线网| av视频免费观看在线观看| 久久久国产一区二区| 亚洲国产毛片av蜜桃av| av片东京热男人的天堂| 日韩欧美国产一区二区入口| 精品熟女少妇八av免费久了| 亚洲精品中文字幕在线视频| 嫁个100分男人电影在线观看| 侵犯人妻中文字幕一二三四区| 欧美在线一区亚洲| av不卡在线播放| 黄频高清免费视频| 成人永久免费在线观看视频 | 国产一区二区 视频在线| 另类亚洲欧美激情| √禁漫天堂资源中文www| av在线播放免费不卡| 大码成人一级视频| 丁香欧美五月| 婷婷丁香在线五月| 美女午夜性视频免费| 大码成人一级视频| 在线观看一区二区三区激情| 国产成人欧美在线观看 | 纯流量卡能插随身wifi吗| 美女主播在线视频| 女人被躁到高潮嗷嗷叫费观| 久久久国产成人免费| 成人av一区二区三区在线看| 午夜精品久久久久久毛片777| 12—13女人毛片做爰片一| 别揉我奶头~嗯~啊~动态视频| 成人av一区二区三区在线看| 中文字幕人妻熟女乱码| 久久久久久久国产电影| 视频区图区小说| 人人妻,人人澡人人爽秒播| 日本vs欧美在线观看视频| 国产亚洲精品第一综合不卡| 在线av久久热| 日韩视频在线欧美| 多毛熟女@视频| 亚洲成人国产一区在线观看| 9色porny在线观看|