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

    基于線粒體COⅠ的南海北部長棘銀鱸遺傳多樣性分析

    2020-06-30 08:08:48裴麗梅宮亞運黃柏炎王業(yè)磷
    海洋漁業(yè) 2020年3期
    關(guān)鍵詞:碣石瓊海陵水

    裴麗梅,宮亞運,黃柏炎,章 群,徐 示,王業(yè)磷

    (暨南大學(xué)生態(tài)學(xué)系,廣州 510632)

    長棘銀鱸(Gerres filamentosus)隸屬于鱸形目(Perciformes)銀鱸科(Gerreidae)銀鱸屬,廣泛分布于世界熱帶和亞熱帶海域和部分河口區(qū)[1-2],一般生活在淺海沿岸50 m以深的沙質(zhì)底部,以生活在沙子和泥漿底部的小型生物為食[3-4]。長棘銀鱸為小型優(yōu)質(zhì)魚類 ,我國見于南海和東海南部水域,歷史上北部灣曾盛產(chǎn),但近年來由于過度捕撈,野生資源已經(jīng)嚴(yán)重衰竭,資源密度由1961—1962年的 78.4 kg·km-2下降至 1992—1993年以后的 0.1 kg·km-2[5]。隨著長棘銀鱸等魚類資源的衰竭,北部灣海洋生態(tài)系統(tǒng)發(fā)生了明顯變化,如游泳動物和底棲生物食性魚類種類數(shù)日益減少,而浮游生物食性魚類種類數(shù)及其漁獲比例明顯增加[6]。雖然自1999年以來實行的休漁制度對北部灣部分魚類資源恢復(fù)有明顯效果,但是對衰竭型種類效果甚微[5],長棘銀鱸種質(zhì)資源保護、恢復(fù)迫在眉睫。

    了解種群的遺傳變異是評估資源現(xiàn)狀、制定野生群體保護與恢復(fù)策略和實施養(yǎng)殖群體種質(zhì)遺傳改良的重要前提[7]。目前國內(nèi)外對于長棘銀鱸的研究大多集中在形態(tài)比較、物種分類、生長發(fā)育和線粒體全基因測序等方面[7-13],尚未有長棘銀鱸種群遺傳相關(guān)的研究報道。線粒體DNA(mtDNA)具有進(jìn)化速度快、母系遺傳、幾乎不發(fā)生重組等特點[14-15],已被用作脊椎動物的系統(tǒng)發(fā)育和群體遺傳分析的有效工具[16-17]。其中細(xì)胞色素氧化酶Ⅰ(cytochrome oxidase subunitⅠ,COⅠ)序列進(jìn)化速率適中,是開展動物DNA條形碼研究首選分子標(biāo)記[18],用于種群遺傳學(xué)研究時,也可在一定程度上避免因線粒體控制區(qū)易變異飽和、有時不能有效檢測到種群結(jié)構(gòu)以及學(xué)術(shù)上對線粒體控制區(qū)序列插入和缺失狀況處理的共識不足[19-20]。近年來國內(nèi)外有不少使用COⅠ基因研究魚類遺傳變異的報道,如日本鯖(Scomber japonicus)[21]、 銀 鯧 (Pampus argenteus)[22]和烏賊(Sepia officinalis)[23]等。本研究采用線粒體COⅠ基因序列對南海北部7個長棘銀鱸群體進(jìn)行分析,旨在揭示中國長棘銀鱸主要分布區(qū)的遺傳多樣性和種群結(jié)構(gòu),為種質(zhì)資源保護和開發(fā)利用提供科學(xué)依據(jù)。

    1 材料與方法

    1.1 材料來源

    2014年7月至2019年10月間采自廣東碣石(JS)34尾、廣東陽江(YJ)34尾、海南瓊海(QH)14尾、海南陵水(LS)24尾、海南東方(DF)20尾、海南新盈(XY)26尾、廣西東興(DX)24尾共7個地點的176尾長棘銀鱸(表1),用95%的酒精處理并保存。所有樣品均為從近海作業(yè)的小型漁船上購買的野生魚類。

    1.2 基因組DNA的提取與PCR擴增

    DNA提取和PCR擴增條件參照樂小亮等[24]方法,PCR擴增引物參見WARD等[25],經(jīng)瓊脂糖凝膠電泳檢測合格的PCR擴增產(chǎn)物送華大基因公司測序。

    1.3 數(shù)據(jù)處理

    將測序后得到的序列在 Bioedit 7.0.9.0[26]軟件上進(jìn)行人工校對和多重比較。運用MEGA 7.0[27]軟件計算核苷酸組成、堿基變異情況、轉(zhuǎn)換和顛換值。通過 DnaSP 6.12.03[28]軟件計算單倍型數(shù)、單倍型多樣性(haplotype diversity,Hd)、核苷酸多樣性(nucleotide diversity,π)。使用Network 5.0[29]構(gòu) 建 單 倍 型 網(wǎng) 絡(luò) 圖。通 過Arlequin 3.5.2.2[30]軟件計算遺傳分化系數(shù)(Fst)和群體間基因流值(Nm);進(jìn)行AMOVA分子方差分析遺傳變異的分布情況,并以1 000次置換分析統(tǒng)計學(xué)顯著性;計算Tajima’sD和Fu’sFs值并進(jìn)行中性檢驗,得到SSD值、Rg值、τ值參數(shù)。根據(jù)公式T=(τ/2μk)×代時[31]計算種群擴張時間。其中,τ是種群擴張時間參數(shù);μ是序列的變異速率 [采用 1% ~3%/MY[32](million year,MY)];k是序列長度;T是種群擴張時間。

    2 結(jié)果與分析

    2.1 COⅠ基因片段序列特征及遺傳多樣性

    在長棘銀鱸COⅠ基因5′端652 bp序列中,T、C、A和 G含量分別為 29.0%、28.7%、23.3%和19.1%。A+T的含量(52.3%)略高于C+G的含量(47.8%),與大多數(shù)硬骨魚類相似[33]。共發(fā)現(xiàn)單倍型25個,變異位點155個,其中簡約信息位點87個。轉(zhuǎn)換和顛換比為2.41,表明序列未達(dá)到飽和,適合進(jìn)行種群遺傳發(fā)育分析[34]。7個群體整體呈現(xiàn)較高的單倍型多樣性(Hd=0.767±0.018)與較低的核苷酸多樣性(π=0.0550±0.003 7)。其中遺傳多樣性指數(shù)最高的是碣石群體(Hd=0.606±0.092,π=0.059 2±0.006 6),其次為陽江群體(Hd=0.323±0.102,π=0.027 0±0.010 7);瓊海群體遺傳多樣性指數(shù)最低(Hd=0.560±0.125,π=0.001 0±0.000 3),陵水群體次低(Hd=0.685±0.077,π=0.0015±0.000 3)。

    2.2 長棘銀鱸遺傳分化

    176條COⅠ序列中共檢測到25個單倍型,單倍型頻率最高的為 Hap1(60條),其次是Hap15(50條)。共享單倍型4個,其中 Hap1、Hap2和Hap3為所有群體共享;Hap15為碣石、陽江、新盈3個群體共享;其它單倍型為各群體特有。

    用 Arlequin 3.5.2.2軟件以 Kimura雙參數(shù)模型構(gòu)建的單倍型網(wǎng)絡(luò)中呈現(xiàn)非典型的星狀結(jié)構(gòu),沒有較為明顯的地理聚群和譜系結(jié)構(gòu)(圖1)。7個群體間的遺傳分化情況(Fst=-0.025 5~0.851 4,Nm=-68.993~56.971)詳見表 2。其中,碣石和陽江群體分別與瓊海、陵水、東方、新盈、東興5個群體間出現(xiàn)極顯著性遺傳分化(Fst>0.25,P<0.001)[35],且碣石和陽江群體分別與瓊海、陵水、東方、新盈、東興5個群體之間基因交流匱乏(0<Nm<1)[36-37]。碣石和陽江群體間雖然為隨機交配(Nm>4),但也出現(xiàn)了低等程度分化(0.05<Fst=0.111 1<0.15,P<0.05),其余組別之間無明顯分化。

    將7個群體通過4種不同方式進(jìn)行AMOVA分析:①7個地理群體;②大陸群組:碣石、陽江、東興;海島群組:瓊海、陵水、東方、新盈;③碣石和陽江為一組;瓊海、陵水、東方、新盈、東興為另一組;④碣石群體;陽江群體;其他群體(瓊海、陵水、東方、新盈)。AMOVA分析(表3)表明4種分組方式的主要變異來源均發(fā)生在群體內(nèi),變異比例分別為 74.13%、68.87%、61.69% 和64.76%,同時4種分組方法下4種群組間存在顯著性分化(Fst>0.25,P<0.001),原因可能與長棘銀鱸只進(jìn)行短距離洄游有關(guān)[5]。

    表1 南海長棘銀鱸群體的遺傳多樣性Tab.1 Genetic diversity of G.filamentosus in coastal waters of the South China Sea

    圖1 長棘銀鱸單倍型網(wǎng)絡(luò)圖Fig.1 Parsimony network of G.filamentosus hap lotypes

    2.3 長棘銀鱸種群歷史動態(tài)

    中性檢驗結(jié)果見表4。除陵水 Tajima’sD(-1.518 6,0.001<P<0.05)和 Fu’sFs(-3.293 6,P<0.05)、東 方 Tajima’sD(-2.590 6,P<0.001)、新 盈 Tajima’sD(-2.650 4,P<0.001)為顯著性負(fù)值外,其余群體Tajima’sD和Fu’sFs均為不顯著性負(fù)值或正值。結(jié)合陵水、東方、新盈群體SSD值及Rg值數(shù)值均較小且不具有統(tǒng)計顯著性,表明上述觀測值符合擴張模型。因Fu’sFs檢驗比Tajima’sD對近期擴張更敏感[38],再結(jié)合上述3個群體的核苷酸錯配圖(圖2-a~c)為明顯單峰以及單倍型網(wǎng)絡(luò)圖為非典型星型結(jié)構(gòu),可認(rèn)為陵水群體經(jīng)歷了近期的快速擴張,而東方和新盈群體經(jīng)歷過較為早期的擴張[39-40]。7個群體整體的 Tajima’sD=-0.765 9(P>0.1)、Fu’sFs=6.849 1(P>0.5),整體上未呈現(xiàn)擴張現(xiàn)象,結(jié)合核苷酸錯配圖(圖2-d)呈現(xiàn)明顯多峰,說明可能是種群經(jīng)歷多次擴張疊加的結(jié)果[40]。陵水、東方、新盈3個群體擴張參數(shù) τ分別為 1.068 4、0.949 2和1.175 8,以線粒體COⅠ基因進(jìn)化速率為(1%~3%)/百萬年[32],性成熟期 3年,估算出陵水、東方、新盈3個群體擴張大約分別發(fā)生在0.245 8~0.081 9、0.218 4~0.072 8、0.270 5~0.090 2百萬年之前,均處于晚更新世時期。

    表2 長棘銀鱸7群體的F st及N m值Tab.2 F st and N m among 7 G.filamentosus populations

    表3 長棘銀鱸種群結(jié)構(gòu)分子方差分析Tab.3 AMOVA analysis of G.filamentosus

    表4 長棘銀鱸中性檢驗和不對稱分析Tab.4 Neutrality tests and m ismatch analysis of G.filamentosus

    圖2 長棘銀鱸核苷酸不對稱分布Fig.2 Nucleotidem ismatch distribution of G.filamentosus

    3 討論

    3.1 長棘銀鱸的遺傳多樣性與種群歷史動態(tài)

    種群遺傳多樣性的高低很大程度上與物種的適應(yīng)能力、分化能力以及對環(huán)境影響的耐受力相關(guān)[41],而單倍型多樣性(Hd)和核酸多樣性(π)是衡量某一物種遺傳多樣性的2個重要指標(biāo)[42]。長棘銀鱸COⅠ序列的遺傳多樣性(Hd=0.767±0.018,π=0.055 0±0.003 7)呈現(xiàn)的高Hd值低π值的特點,與同海域基于COⅠ序列的銀鯧 (Pampus argenteus)[22]和 湖 棲 蝦 虎 魚(Gobiopterus lacustris)[43]等相似。南海北部[碣石(0.059 2±0.006 6)和陽江(0.027 0±0.010 7)]的遺傳多樣性明顯高于南海南部[瓊海(0.001 0±0.000 3)和陵水(0.001 5±0.000 3)],推測可能與碣石和陽江所處的海域有關(guān)。在冰期大陸架露出時,南海變成邊緣海,陽江和碣石群體瀕臨該邊緣海,種群相對穩(wěn)定,受瓶頸效果影響?。涣硪环矫?,碣石與陽江所處海域會同時受到南海暖流、黑潮及沿岸流的影響,使得碣石與陽江群體具有更強的分化能力和環(huán)境耐受力。而瓊海、陵水、東方、新盈的潮汐汊道港灣特征使內(nèi)灣水域處于比較穩(wěn)定的半封閉海洋環(huán)境[43];東興則位于半封閉的北部灣,海洋環(huán)境也相對穩(wěn)定,遺傳多樣性偏低。此外,冰期時海南島及北部灣所在的大陸架露出,間冰期時海平面上升[44-45],海南島和北部灣的長棘銀鱸群體發(fā)生擴張,也致使瓊海、陵水、東方、新盈和東興群體遺傳多樣性偏低。

    Fst值顯示陽江和碣石群體分別與瓊海、陵水、東方、新盈、東興5個群體間發(fā)生極顯著性遺傳分化(Fst>0.25,P<0.001),且 AMOVA分析結(jié)果顯示長棘銀鱸4種劃分方式的群組間都有顯著性分化。長棘銀鱸在河口生活,繁殖期遷移至沿海地區(qū)達(dá)到成熟的最后階段,然后進(jìn)行產(chǎn)卵,不進(jìn)行遠(yuǎn)距離洄游[3]。瓊州海峽是在全新世中期全球性海侵淹沒峽區(qū)原來的常態(tài)低地而形成,主要形成時段為(10 570±560)~(7 125±96)aBP[46-47]。瓊州海峽形成前,黑潮從巴士海峽進(jìn)入南海,同沿岸流一樣在沿廣東海岸流向西南過程中受阻,導(dǎo)致碣石和陽江群體與海南各群體及東興群體間發(fā)生遺傳分化;瓊州海峽形成后,黑潮從巴士海峽進(jìn)入南海,同沿岸流一樣在沿廣東海岸流向西南過程中受瓊州海峽和南海大陸架特殊地況的影響,也在一定程度上阻礙了碣石和陽江群體的遠(yuǎn)距離傳播,導(dǎo)致碣石和陽江群體與海南各群體及東興群體間有明顯的遺傳分化。

    大部分中國近海魚類在近期歷史上因冰期海平面變化而出現(xiàn)種群快速擴張[45]。長棘銀鱸Tajima’sD和Fu’sFs中性檢驗和核苷酸不對稱分析圖呈現(xiàn)單峰,顯示陵水、東方和新盈群體都經(jīng)歷過種群擴張,推算擴張時間均處于晚更新世時期。海南島為南海大陸架島嶼,在晚更新世構(gòu)造運動頻繁,火山活動活躍,使不少海灣底部抬升為階地;全新世氣候變暖,海平面上升,海水侵入河口、海灣、洼地等,使長棘銀鱸群體隨著海侵發(fā)生擴張[44]。

    3.2 長棘銀鱸種質(zhì)資源保護

    長棘銀鱸生活在淺海沿岸水域,是包括珊瑚礁、紅樹林等生境生態(tài)系統(tǒng)的重要組分[2-3]。近年來捕撈技術(shù)的發(fā)展與革新使得長棘銀鱸資源快速衰退[5-6],長棘銀鱸野生資源的保護迫在眉睫。碣石和陽江群體分別與瓊海、陵水、東方、新盈群體間有極顯著性分化,可以將碣石和陽江以及其他群體作為2個獨立的管理保護單位。其中碣石群體遺傳多樣性最高(Hd=0.606±0.092,π=0.059 2±0.006 6),應(yīng)作為長棘銀鱸種質(zhì)資源保護和開發(fā)利用的核心群體和優(yōu)先對象;瓊海(Hd=0.560±0.125,π=0.001 0±0.000 3)和陵水(Hd=0.685±0.077,π=0.001 5±0.000 3)群體遺傳多樣性較低,更容易受到環(huán)境變化的影響,建議及時采取措施,避免捕撈過度使其遺傳多樣性難以恢復(fù),重蹈大黃魚、帶魚等種質(zhì)資源衰退的覆轍。

    隨著分子生物學(xué)技術(shù)的進(jìn)步,特別是高通量測序技術(shù)的進(jìn)步和測序成本的下降,種群遺傳學(xué)研究開始進(jìn)入基因組時代。線粒體DNA有效群體大小不僅為雙親遺傳核基因組的1/4,可提供核基因組不能提供的母系遺傳信息,開展種群遺傳研究的成本遠(yuǎn)低于種群全基因組分析,仍是種群遺傳的首選分子標(biāo)記[48]?;诰€粒體COⅠ的長棘銀鱸種群遺傳研究發(fā)現(xiàn),南海北部瓊州海峽兩側(cè)的群體間存在極顯著的遺傳分化,可作為2個不同的管理保護單位;廣東沿海群體的遺傳多樣性明顯高于海南和廣西群體,其中碣石群體遺傳多樣性最高,可優(yōu)先保護。研究結(jié)果為長棘銀鱸種群的保護和開發(fā)利用措施的制訂與實施提供了科學(xué)依據(jù)。但受采樣范圍的限制,研究結(jié)果不能排除未采樣地區(qū)存在種群分化的可能,今后需要進(jìn)行物種全分布區(qū)的采樣調(diào)查。另外,相對于其他線粒體基因,COⅠ基因的進(jìn)化速度較慢,較適合發(fā)現(xiàn)歷史相對久遠(yuǎn)的種群擴張,對近期的種群擴張不太敏感,在后續(xù)有關(guān)長棘銀鱸種群遺傳的研究中,可考慮使用高多態(tài)性的微衛(wèi)星指紋圖譜分析技術(shù)以及基于二代測序的簡化基因組乃至全基因組分析技術(shù),以揭示雙親遺傳的核基因遺傳信息,反映更為近期的全球氣候變化與人類活動的影響,從而制訂更為精準(zhǔn)的保護與開發(fā)利用策略。

    猜你喜歡
    碣石瓊海陵水
    碣石山下歌聲起
    陵水17-2半潛式生產(chǎn)儲油平臺落物風(fēng)險分析
    文昌至瓊海高速預(yù)計11月通車!全長65.7km
    石油瀝青(2019年2期)2019-02-13 17:24:47
    碣石在秦漢時期歷史中的地位
    山西青年(2019年4期)2019-01-15 12:21:00
    在陵水錯把大風(fēng)吹過的地方當(dāng)成故鄉(xiāng)
    陵水陵水真奇美
    新教育(2018年35期)2018-02-27 00:03:03
    大禹治水曾到碣石山
    鄉(xiāng)音(2016年4期)2016-05-28 07:26:22
    海南瓊海辣椒收獲上市
    海南陵水:用豐厚文化家底為旅游添翼
    珠江口盆地瓊海凸起斷裂調(diào)整聚集成藏模式
    久久久精品欧美日韩精品| 久久久久久九九精品二区国产| 日本熟妇午夜| 免费av观看视频| 日韩欧美国产在线观看| 国产成人av教育| 久久久久国内视频| 露出奶头的视频| 免费电影在线观看免费观看| av中文乱码字幕在线| 一区二区三区激情视频| 亚洲欧美日韩东京热| 亚洲三级黄色毛片| 精品午夜福利视频在线观看一区| 精品熟女少妇八av免费久了| 久久亚洲真实| 制服丝袜大香蕉在线| 九色成人免费人妻av| 国内揄拍国产精品人妻在线| 99国产极品粉嫩在线观看| 精品国产亚洲在线| 又黄又爽又免费观看的视频| 国产精品乱码一区二三区的特点| 欧美潮喷喷水| 脱女人内裤的视频| av在线老鸭窝| 久久国产精品人妻蜜桃| 欧洲精品卡2卡3卡4卡5卡区| 国产v大片淫在线免费观看| 久久99热6这里只有精品| 免费在线观看成人毛片| 精品熟女少妇八av免费久了| 1024手机看黄色片| 一本久久中文字幕| 夜夜爽天天搞| 夜夜爽天天搞| 99久久无色码亚洲精品果冻| 制服丝袜大香蕉在线| 听说在线观看完整版免费高清| 亚洲熟妇熟女久久| 黄色女人牲交| 精华霜和精华液先用哪个| 亚洲专区中文字幕在线| 日日摸夜夜添夜夜添小说| 久久久久九九精品影院| 午夜福利免费观看在线| avwww免费| 在线a可以看的网站| 老女人水多毛片| 国产视频内射| 毛片一级片免费看久久久久 | 在线观看免费视频日本深夜| 国产在线男女| 亚洲人成网站在线播| 禁无遮挡网站| 深夜精品福利| www日本黄色视频网| 国产精品久久久久久亚洲av鲁大| 婷婷精品国产亚洲av| 少妇人妻精品综合一区二区 | 天堂av国产一区二区熟女人妻| 麻豆av噜噜一区二区三区| 天美传媒精品一区二区| 久久久久久久久久黄片| 精品久久久久久,| 两个人的视频大全免费| 亚洲午夜理论影院| 国产成人福利小说| 亚洲美女视频黄频| 久久欧美精品欧美久久欧美| 麻豆成人av在线观看| 亚洲在线自拍视频| 欧美中文日本在线观看视频| 男人的好看免费观看在线视频| 床上黄色一级片| 一a级毛片在线观看| 国产精品1区2区在线观看.| 午夜免费男女啪啪视频观看 | 午夜福利成人在线免费观看| 日韩 亚洲 欧美在线| 免费人成视频x8x8入口观看| 午夜免费成人在线视频| 国产精品亚洲一级av第二区| 中文字幕久久专区| 久久午夜亚洲精品久久| 日韩欧美一区二区三区在线观看| 国产精品98久久久久久宅男小说| 成人一区二区视频在线观看| 免费黄网站久久成人精品 | 国产精品电影一区二区三区| 91麻豆av在线| 天天躁日日操中文字幕| 露出奶头的视频| 久久久久久久久大av| 国产av在哪里看| 亚洲中文字幕一区二区三区有码在线看| 91在线观看av| 国产激情偷乱视频一区二区| 成人鲁丝片一二三区免费| 久久久国产成人免费| 久久午夜福利片| 我要搜黄色片| 长腿黑丝高跟| 香蕉av资源在线| 国产一区二区三区视频了| 三级男女做爰猛烈吃奶摸视频| 老熟妇乱子伦视频在线观看| 中文字幕av在线有码专区| 直男gayav资源| 精品日产1卡2卡| 日韩 亚洲 欧美在线| 午夜免费激情av| 日韩av在线大香蕉| 观看美女的网站| 国产三级黄色录像| 在线免费观看的www视频| 中文字幕人成人乱码亚洲影| 久久久久亚洲av毛片大全| 国语自产精品视频在线第100页| 搡老熟女国产l中国老女人| 日韩大尺度精品在线看网址| 亚洲精品乱码久久久v下载方式| 欧美黄色淫秽网站| 18禁黄网站禁片午夜丰满| 欧美成人一区二区免费高清观看| 国产一区二区在线av高清观看| 久久久久久久久大av| 国产伦精品一区二区三区视频9| 午夜影院日韩av| 国产高清视频在线播放一区| 夜夜看夜夜爽夜夜摸| 在线十欧美十亚洲十日本专区| 9191精品国产免费久久| 赤兔流量卡办理| 日日夜夜操网爽| 亚洲在线自拍视频| 久久久久久久亚洲中文字幕 | 国产一区二区三区视频了| 国产白丝娇喘喷水9色精品| 757午夜福利合集在线观看| 性欧美人与动物交配| 最新在线观看一区二区三区| 国产亚洲av嫩草精品影院| 欧美在线一区亚洲| 一个人免费在线观看电影| 成人午夜高清在线视频| 性色av乱码一区二区三区2| 舔av片在线| 国产精品女同一区二区软件 | 亚洲自拍偷在线| 永久网站在线| 中文亚洲av片在线观看爽| 波野结衣二区三区在线| 免费无遮挡裸体视频| av专区在线播放| 人妻丰满熟妇av一区二区三区| 老熟妇乱子伦视频在线观看| 韩国av一区二区三区四区| 91九色精品人成在线观看| 国产色婷婷99| 日日夜夜操网爽| 国产伦精品一区二区三区四那| 精品国内亚洲2022精品成人| 日本一本二区三区精品| 精品久久久久久久久久久久久| 97超级碰碰碰精品色视频在线观看| 好男人在线观看高清免费视频| 性色avwww在线观看| 欧美成人一区二区免费高清观看| 97热精品久久久久久| 嫩草影院入口| 亚洲欧美清纯卡通| 国产成人影院久久av| 99国产精品一区二区蜜桃av| 日韩欧美免费精品| 偷拍熟女少妇极品色| 欧美性猛交黑人性爽| 国产黄a三级三级三级人| 日本精品一区二区三区蜜桃| 国产精品,欧美在线| 国产伦精品一区二区三区视频9| 亚洲片人在线观看| 草草在线视频免费看| 3wmmmm亚洲av在线观看| 午夜福利欧美成人| 亚洲美女视频黄频| 国产精品久久久久久精品电影| 午夜a级毛片| 精品久久久久久久久亚洲 | 人人妻人人看人人澡| 91九色精品人成在线观看| 又黄又爽又刺激的免费视频.| 国内久久婷婷六月综合欲色啪| 真实男女啪啪啪动态图| 啦啦啦观看免费观看视频高清| 99热这里只有精品一区| 国产v大片淫在线免费观看| 欧美丝袜亚洲另类 | 91av网一区二区| av中文乱码字幕在线| 国产乱人视频| 亚洲精品亚洲一区二区| 国产乱人视频| 男人舔奶头视频| ponron亚洲| 人人妻,人人澡人人爽秒播| 成人国产一区最新在线观看| 午夜福利在线在线| av在线观看视频网站免费| 亚洲欧美日韩无卡精品| 嫩草影视91久久| 欧美又色又爽又黄视频| 又粗又爽又猛毛片免费看| 日韩人妻高清精品专区| 亚洲专区中文字幕在线| 午夜免费成人在线视频| 欧美在线黄色| 美女黄网站色视频| 国产欧美日韩精品亚洲av| 97人妻精品一区二区三区麻豆| 亚洲精品成人久久久久久| 亚洲国产精品999在线| 日韩成人在线观看一区二区三区| 免费在线观看影片大全网站| 在线观看av片永久免费下载| 亚洲av成人不卡在线观看播放网| 热99在线观看视频| 国产高清激情床上av| 99在线人妻在线中文字幕| 在线观看美女被高潮喷水网站 | 露出奶头的视频| 精品免费久久久久久久清纯| 国产精品久久视频播放| 青草久久国产| 一区二区三区免费毛片| 亚洲中文字幕日韩| 国产探花极品一区二区| 精品不卡国产一区二区三区| 国产午夜福利久久久久久| 嫩草影视91久久| 亚洲精品日韩av片在线观看| 九九在线视频观看精品| 国产精品综合久久久久久久免费| 色在线成人网| 久久久久亚洲av毛片大全| 午夜a级毛片| 国产免费av片在线观看野外av| 99在线人妻在线中文字幕| 免费大片18禁| 18美女黄网站色大片免费观看| 日韩亚洲欧美综合| 日韩 亚洲 欧美在线| 久久久精品大字幕| 中文字幕av在线有码专区| 伦理电影大哥的女人| 精品国产亚洲在线| 人人妻人人澡欧美一区二区| 久久99热这里只有精品18| 色哟哟·www| 两个人视频免费观看高清| 精品欧美国产一区二区三| 热99re8久久精品国产| 国产单亲对白刺激| 最近最新中文字幕大全电影3| 国产伦精品一区二区三区四那| 日本黄色视频三级网站网址| 亚洲片人在线观看| 一二三四社区在线视频社区8| 嫩草影院精品99| 少妇裸体淫交视频免费看高清| 日韩av在线大香蕉| 国产美女午夜福利| 蜜桃久久精品国产亚洲av| 精品国内亚洲2022精品成人| 美女被艹到高潮喷水动态| 欧美绝顶高潮抽搐喷水| 精品不卡国产一区二区三区| 成人国产一区最新在线观看| 日韩中字成人| 变态另类丝袜制服| 男插女下体视频免费在线播放| 狂野欧美白嫩少妇大欣赏| 欧美3d第一页| 最近最新中文字幕大全电影3| 国产精品野战在线观看| 亚洲色图av天堂| 国内少妇人妻偷人精品xxx网站| 亚洲自拍偷在线| 国产私拍福利视频在线观看| 久久久久久久久久黄片| 中文字幕熟女人妻在线| 国产欧美日韩一区二区三| 色哟哟·www| 国产高清有码在线观看视频| 男人的好看免费观看在线视频| 午夜激情福利司机影院| 精品久久久久久久久久久久久| 赤兔流量卡办理| 一个人看的www免费观看视频| 久久热精品热| 男女做爰动态图高潮gif福利片| 国产av不卡久久| 欧美性感艳星| 午夜免费激情av| 很黄的视频免费| 又爽又黄a免费视频| 亚洲午夜理论影院| 99久久九九国产精品国产免费| 夜夜看夜夜爽夜夜摸| 少妇人妻一区二区三区视频| 欧美在线黄色| 亚洲av.av天堂| 美女xxoo啪啪120秒动态图 | 99riav亚洲国产免费| 中国美女看黄片| 国模一区二区三区四区视频| 亚洲国产色片| 亚洲第一电影网av| 99久久精品一区二区三区| 亚洲经典国产精华液单 | 久久中文看片网| 国产一区二区在线av高清观看| 欧美激情在线99| 在线免费观看不下载黄p国产 | 久久精品夜夜夜夜夜久久蜜豆| 亚洲国产欧美人成| 一区二区三区免费毛片| 一进一出抽搐动态| 九色国产91popny在线| 亚洲国产精品合色在线| 男人舔奶头视频| 精品午夜福利在线看| 国产伦人伦偷精品视频| 亚洲人成电影免费在线| 国产aⅴ精品一区二区三区波| 欧美乱色亚洲激情| 欧美日韩综合久久久久久 | 一夜夜www| 可以在线观看的亚洲视频| 天堂影院成人在线观看| 一个人看的www免费观看视频| 成人特级黄色片久久久久久久| 国产主播在线观看一区二区| 亚洲狠狠婷婷综合久久图片| 亚洲欧美激情综合另类| 亚洲一区二区三区不卡视频| 自拍偷自拍亚洲精品老妇| 国产高清激情床上av| 免费大片18禁| 国产精品美女特级片免费视频播放器| 成人鲁丝片一二三区免费| 窝窝影院91人妻| 日本成人三级电影网站| 欧美在线黄色| a级毛片a级免费在线| 久久午夜亚洲精品久久| 国产成人欧美在线观看| 麻豆一二三区av精品| 亚洲一区二区三区色噜噜| 婷婷六月久久综合丁香| 九色成人免费人妻av| 18禁裸乳无遮挡免费网站照片| 精品午夜福利视频在线观看一区| 麻豆一二三区av精品| 91狼人影院| 嫩草影视91久久| 亚洲av美国av| 99久久无色码亚洲精品果冻| x7x7x7水蜜桃| 午夜福利成人在线免费观看| 免费av观看视频| 国内少妇人妻偷人精品xxx网站| 欧美午夜高清在线| 夜夜爽天天搞| 老司机福利观看| 国产精品人妻久久久久久| 国产精品一区二区性色av| 精品国产三级普通话版| 色综合站精品国产| 国产精品嫩草影院av在线观看 | 国产乱人视频| 男人舔女人下体高潮全视频| 99国产综合亚洲精品| 中亚洲国语对白在线视频| 欧美xxxx黑人xx丫x性爽| 婷婷丁香在线五月| 国产伦精品一区二区三区四那| 精品一区二区三区视频在线观看免费| 国产高清激情床上av| 搡老妇女老女人老熟妇| 精品欧美国产一区二区三| 淫妇啪啪啪对白视频| 在线天堂最新版资源| 精品一区二区免费观看| 天美传媒精品一区二区| 首页视频小说图片口味搜索| 亚洲天堂国产精品一区在线| 国产精品1区2区在线观看.| 美女大奶头视频| 99国产综合亚洲精品| 啪啪无遮挡十八禁网站| 精品99又大又爽又粗少妇毛片 | 在线十欧美十亚洲十日本专区| 精品久久久久久久久亚洲 | 午夜免费成人在线视频| 欧美一区二区国产精品久久精品| 给我免费播放毛片高清在线观看| 亚洲av电影不卡..在线观看| 欧美日本视频| 中文字幕av成人在线电影| 亚洲,欧美精品.| 少妇裸体淫交视频免费看高清| 中文字幕精品亚洲无线码一区| 人人妻人人看人人澡| 自拍偷自拍亚洲精品老妇| 欧美黑人巨大hd| 精品久久久久久久久久免费视频| 91字幕亚洲| 欧美在线一区亚洲| 国产色爽女视频免费观看| 国产伦在线观看视频一区| 神马国产精品三级电影在线观看| 亚州av有码| 高清毛片免费观看视频网站| 亚洲男人的天堂狠狠| 老司机深夜福利视频在线观看| 日韩人妻高清精品专区| 十八禁网站免费在线| 真人一进一出gif抽搐免费| 精品久久久久久久久久久久久| 一个人看的www免费观看视频| 亚洲成人中文字幕在线播放| 成人午夜高清在线视频| 免费看日本二区| 757午夜福利合集在线观看| 国产av不卡久久| 中文字幕免费在线视频6| 亚洲一区二区三区色噜噜| 欧美在线一区亚洲| 色尼玛亚洲综合影院| 在线观看av片永久免费下载| 一个人观看的视频www高清免费观看| h日本视频在线播放| 精品久久久久久,| 一本一本综合久久| 成人高潮视频无遮挡免费网站| 床上黄色一级片| 欧洲精品卡2卡3卡4卡5卡区| 91av网一区二区| 国产成人影院久久av| 欧美高清成人免费视频www| 久久久久久久久大av| 欧美3d第一页| 国产免费一级a男人的天堂| 免费人成视频x8x8入口观看| 小蜜桃在线观看免费完整版高清| 美女黄网站色视频| 蜜桃亚洲精品一区二区三区| 999久久久精品免费观看国产| 嫩草影院新地址| 欧美最黄视频在线播放免费| 精品国产三级普通话版| 99热6这里只有精品| 国产精品影院久久| 男女之事视频高清在线观看| 看免费av毛片| 草草在线视频免费看| 麻豆国产av国片精品| 久久99热6这里只有精品| 国产欧美日韩精品亚洲av| 国产精品爽爽va在线观看网站| 国产视频内射| 中文在线观看免费www的网站| 看片在线看免费视频| 最好的美女福利视频网| 亚洲欧美日韩卡通动漫| 两人在一起打扑克的视频| 亚洲在线自拍视频| 女同久久另类99精品国产91| 色哟哟·www| 日本黄色片子视频| 欧美日韩综合久久久久久 | 国产伦人伦偷精品视频| 校园春色视频在线观看| 乱码一卡2卡4卡精品| 精品久久久久久久久久免费视频| 黄色一级大片看看| 精品久久久久久,| 他把我摸到了高潮在线观看| 亚洲性夜色夜夜综合| 麻豆av噜噜一区二区三区| 欧美日本视频| 亚洲专区中文字幕在线| 国产欧美日韩一区二区三| 欧美黑人巨大hd| 国产精品久久久久久久久免 | 成人一区二区视频在线观看| av福利片在线观看| 国产精品久久视频播放| 成年女人看的毛片在线观看| 久久天躁狠狠躁夜夜2o2o| 每晚都被弄得嗷嗷叫到高潮| 最近在线观看免费完整版| 日本与韩国留学比较| 最新在线观看一区二区三区| 色综合亚洲欧美另类图片| 精品久久久久久久久久久久久| 国产av不卡久久| 麻豆一二三区av精品| 狠狠狠狠99中文字幕| 亚洲色图av天堂| 两个人视频免费观看高清| 午夜两性在线视频| 免费在线观看亚洲国产| 看十八女毛片水多多多| 麻豆国产av国片精品| 亚洲国产高清在线一区二区三| 最后的刺客免费高清国语| 日韩中文字幕欧美一区二区| 一本一本综合久久| 国产精品一区二区免费欧美| 日本精品一区二区三区蜜桃| 国产伦一二天堂av在线观看| 九九热线精品视视频播放| 91九色精品人成在线观看| 精品久久久久久久久av| 亚洲va日本ⅴa欧美va伊人久久| а√天堂www在线а√下载| 日本熟妇午夜| 亚洲精品久久国产高清桃花| 亚洲精品在线美女| 久久99热6这里只有精品| 深夜a级毛片| 亚洲第一电影网av| 亚洲熟妇中文字幕五十中出| 国产激情偷乱视频一区二区| 嫩草影院新地址| 少妇熟女aⅴ在线视频| 国产精品精品国产色婷婷| 亚洲经典国产精华液单 | 国产精品一区二区三区四区免费观看 | 午夜精品久久久久久毛片777| 色在线成人网| 一本精品99久久精品77| 国产精品自产拍在线观看55亚洲| 欧美日韩综合久久久久久 | 老司机福利观看| 一区福利在线观看| 国产激情偷乱视频一区二区| 亚洲性夜色夜夜综合| 国产视频一区二区在线看| 亚洲成av人片免费观看| 全区人妻精品视频| 成人性生交大片免费视频hd| 国产精品久久久久久亚洲av鲁大| 日本在线视频免费播放| 深夜精品福利| 中文字幕久久专区| 亚洲在线观看片| 18禁黄网站禁片免费观看直播| 国产精品自产拍在线观看55亚洲| 亚洲欧美日韩高清专用| 可以在线观看毛片的网站| 国产精品三级大全| 在线观看舔阴道视频| 日本精品一区二区三区蜜桃| 伊人久久精品亚洲午夜| 嫩草影院入口| 99国产精品一区二区三区| 久9热在线精品视频| 嫩草影视91久久| 级片在线观看| 九色成人免费人妻av| 午夜免费男女啪啪视频观看 | 日韩成人在线观看一区二区三区| 成人亚洲精品av一区二区| 天堂动漫精品| 亚洲午夜理论影院| 18禁裸乳无遮挡免费网站照片| 99久久久亚洲精品蜜臀av| 一卡2卡三卡四卡精品乱码亚洲| 免费av毛片视频| 日韩大尺度精品在线看网址| 两个人视频免费观看高清| 极品教师在线视频| 国产白丝娇喘喷水9色精品| 搡老岳熟女国产| 亚洲精品影视一区二区三区av| 国产毛片a区久久久久| 亚洲18禁久久av| 99热这里只有是精品在线观看 | 最新在线观看一区二区三区| 国产午夜精品久久久久久一区二区三区 | 欧美在线一区亚洲| 精品一区二区三区人妻视频| 在线看三级毛片| 88av欧美| 91在线观看av| 国产蜜桃级精品一区二区三区| 亚洲av电影不卡..在线观看| 欧美一区二区国产精品久久精品| 在线十欧美十亚洲十日本专区| 色精品久久人妻99蜜桃| 免费看a级黄色片| 日本三级黄在线观看| 欧美黑人欧美精品刺激| 日本熟妇午夜| 嫁个100分男人电影在线观看| 老熟妇仑乱视频hdxx| 又爽又黄a免费视频| 床上黄色一级片| 赤兔流量卡办理| 禁无遮挡网站| 日本五十路高清| 亚洲avbb在线观看|