• <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)整聚集成藏模式
    欧美最黄视频在线播放免费| 久久久精品欧美日韩精品| 色播亚洲综合网| 99久久国产精品久久久| 成人一区二区视频在线观看| 色哟哟哟哟哟哟| 男女做爰动态图高潮gif福利片| 久久久久久久午夜电影| 好男人在线观看高清免费视频| 校园春色视频在线观看| 香蕉久久夜色| 久久欧美精品欧美久久欧美| 老汉色av国产亚洲站长工具| 欧美乱妇无乱码| 国产精品 欧美亚洲| 欧美日韩国产亚洲二区| 久久人人精品亚洲av| 国产精品一区二区三区四区久久| 欧美成狂野欧美在线观看| 熟女少妇亚洲综合色aaa.| 色综合婷婷激情| 亚洲国产高清在线一区二区三| 成人三级黄色视频| 欧美一区二区国产精品久久精品| 精品熟女少妇八av免费久了| 精品国内亚洲2022精品成人| 麻豆久久精品国产亚洲av| 草草在线视频免费看| 亚洲九九香蕉| 很黄的视频免费| 亚洲欧美日韩高清在线视频| 1000部很黄的大片| АⅤ资源中文在线天堂| 欧美日韩一级在线毛片| 日韩欧美国产在线观看| 又紧又爽又黄一区二区| 怎么达到女性高潮| 在线观看66精品国产| 精华霜和精华液先用哪个| 国产一区二区三区在线臀色熟女| 99热6这里只有精品| 两个人看的免费小视频| 国产99白浆流出| 国产主播在线观看一区二区| 老鸭窝网址在线观看| 欧美黄色淫秽网站| 搞女人的毛片| 男人舔奶头视频| 天堂动漫精品| 亚洲国产色片| 亚洲av成人精品一区久久| 欧美不卡视频在线免费观看| 精品久久久久久,| 亚洲无线在线观看| 在线观看免费午夜福利视频| 夜夜躁狠狠躁天天躁| 在线观看66精品国产| 午夜福利欧美成人| 男女之事视频高清在线观看| 又爽又黄无遮挡网站| 午夜日韩欧美国产| 欧美日本亚洲视频在线播放| 日本免费a在线| 老司机午夜十八禁免费视频| 国产亚洲精品av在线| ponron亚洲| 波多野结衣高清作品| 俺也久久电影网| 国产av一区在线观看免费| 在线十欧美十亚洲十日本专区| 最好的美女福利视频网| 嫩草影院精品99| 一进一出抽搐动态| 中文字幕av在线有码专区| 欧美日韩黄片免| 久久天堂一区二区三区四区| 亚洲精华国产精华精| 国产精品久久电影中文字幕| 757午夜福利合集在线观看| 欧美日本亚洲视频在线播放| 91在线精品国自产拍蜜月 | 久久这里只有精品19| 成人欧美大片| 亚洲专区中文字幕在线| 琪琪午夜伦伦电影理论片6080| 国产真实乱freesex| 人妻丰满熟妇av一区二区三区| 久久香蕉精品热| 热99在线观看视频| 国模一区二区三区四区视频 | 亚洲精品美女久久久久99蜜臀| 俺也久久电影网| 国产精品一区二区免费欧美| 国产极品精品免费视频能看的| 在线免费观看的www视频| 国产精品久久视频播放| 级片在线观看| 亚洲成人久久性| 老汉色av国产亚洲站长工具| x7x7x7水蜜桃| 欧美日韩瑟瑟在线播放| 精品国产美女av久久久久小说| 免费搜索国产男女视频| 九色国产91popny在线| 成人特级黄色片久久久久久久| 久久天堂一区二区三区四区| 在线观看免费午夜福利视频| 国产亚洲精品av在线| 亚洲真实伦在线观看| 噜噜噜噜噜久久久久久91| 免费在线观看成人毛片| 久久久久国产精品人妻aⅴ院| 巨乳人妻的诱惑在线观看| 精品免费久久久久久久清纯| 丰满人妻熟妇乱又伦精品不卡| 亚洲色图av天堂| avwww免费| 免费人成视频x8x8入口观看| 国产 一区 欧美 日韩| 国产亚洲av高清不卡| 宅男免费午夜| 香蕉国产在线看| 窝窝影院91人妻| 一进一出抽搐gif免费好疼| 19禁男女啪啪无遮挡网站| 操出白浆在线播放| 欧美日本视频| 搡老岳熟女国产| 亚洲欧美日韩卡通动漫| 亚洲在线观看片| 午夜福利高清视频| 成在线人永久免费视频| 欧美极品一区二区三区四区| 国产亚洲精品久久久久久毛片| 在线免费观看的www视频| 久久九九热精品免费| 国内精品久久久久精免费| 欧美一区二区精品小视频在线| 国产一区二区在线av高清观看| 99热6这里只有精品| 成人永久免费在线观看视频| 巨乳人妻的诱惑在线观看| 最近最新免费中文字幕在线| 精品久久久久久久久久久久久| 日本成人三级电影网站| 久久香蕉国产精品| e午夜精品久久久久久久| 亚洲自偷自拍图片 自拍| 亚洲国产看品久久| 青草久久国产| 色综合婷婷激情| 一级a爱片免费观看的视频| 嫩草影院入口| 中文字幕最新亚洲高清| 男人舔女人的私密视频| 日本免费a在线| 国产精品久久电影中文字幕| 每晚都被弄得嗷嗷叫到高潮| 精品国产乱码久久久久久男人| www国产在线视频色| 黄色日韩在线| 国产欧美日韩精品亚洲av| 亚洲成人久久爱视频| 国产熟女xx| 国产又色又爽无遮挡免费看| 国产亚洲av嫩草精品影院| 中文字幕久久专区| 亚洲欧美日韩无卡精品| 桃红色精品国产亚洲av| 丁香六月欧美| 日本与韩国留学比较| 日韩欧美免费精品| 国产又黄又爽又无遮挡在线| www.自偷自拍.com| 99国产精品一区二区蜜桃av| 欧美一级a爱片免费观看看| 日韩精品青青久久久久久| 欧美性猛交╳xxx乱大交人| 色综合欧美亚洲国产小说| 又紧又爽又黄一区二区| a级毛片a级免费在线| 婷婷丁香在线五月| 成人特级黄色片久久久久久久| 19禁男女啪啪无遮挡网站| 桃色一区二区三区在线观看| 成人国产综合亚洲| 国产成人福利小说| 精品久久久久久成人av| 大型黄色视频在线免费观看| 亚洲九九香蕉| xxx96com| 国产精品久久久久久亚洲av鲁大| 免费看十八禁软件| 中文字幕高清在线视频| 色吧在线观看| 欧美中文日本在线观看视频| 不卡一级毛片| 国产成人精品无人区| 国产综合懂色| 国产男靠女视频免费网站| 麻豆av在线久日| 国产欧美日韩精品一区二区| 国产精品av视频在线免费观看| 成年人黄色毛片网站| 国产亚洲精品一区二区www| 国产精品av视频在线免费观看| 99在线视频只有这里精品首页| 国产精品亚洲美女久久久| 在线观看一区二区三区| 午夜精品久久久久久毛片777| 国产亚洲精品一区二区www| 丰满的人妻完整版| 国产69精品久久久久777片 | 国产欧美日韩一区二区精品| 精品福利观看| 国产午夜精品久久久久久| 久久久久国产一级毛片高清牌| 亚洲人成网站在线播放欧美日韩| 精品久久久久久久毛片微露脸| 熟女少妇亚洲综合色aaa.| 99久久久亚洲精品蜜臀av| 免费观看人在逋| 亚洲人成网站高清观看| av欧美777| 最近在线观看免费完整版| 欧美乱色亚洲激情| 国产一区二区在线av高清观看| 人妻久久中文字幕网| 两个人的视频大全免费| 嫩草影院入口| 久久久久久久久中文| av欧美777| 国产美女午夜福利| 国产视频一区二区在线看| 一个人免费在线观看的高清视频| 成人18禁在线播放| 国产三级在线视频| 俺也久久电影网| 亚洲一区二区三区不卡视频| 免费电影在线观看免费观看| 国产视频内射| 日本免费一区二区三区高清不卡| 99国产精品一区二区蜜桃av| 丁香六月欧美| 他把我摸到了高潮在线观看| 非洲黑人性xxxx精品又粗又长| av福利片在线观看| а√天堂www在线а√下载| 91老司机精品| 国产精品av视频在线免费观看| 99久久国产精品久久久| 国产一级毛片七仙女欲春2| 制服丝袜大香蕉在线| 高潮久久久久久久久久久不卡| 日韩成人在线观看一区二区三区| 欧美成狂野欧美在线观看| 757午夜福利合集在线观看| 亚洲精品粉嫩美女一区| 亚洲人成网站高清观看| 国产一区在线观看成人免费| 国产欧美日韩一区二区精品| 激情在线观看视频在线高清| 国产三级黄色录像| 国产高清三级在线| 国产高清激情床上av| 黄色片一级片一级黄色片| 欧美日韩瑟瑟在线播放| 免费看美女性在线毛片视频| 国产精品一区二区免费欧美| 99re在线观看精品视频| 日本撒尿小便嘘嘘汇集6| 午夜两性在线视频| 一个人观看的视频www高清免费观看 | 国产精品一区二区精品视频观看| 欧美日韩中文字幕国产精品一区二区三区| 99久久成人亚洲精品观看| 午夜福利在线在线| 亚洲成人久久爱视频| 女同久久另类99精品国产91| 极品教师在线免费播放| 小蜜桃在线观看免费完整版高清| 久久天堂一区二区三区四区| 狂野欧美激情性xxxx| 人妻夜夜爽99麻豆av| 岛国在线免费视频观看| 亚洲国产日韩欧美精品在线观看 | 亚洲五月婷婷丁香| 国产精品一区二区精品视频观看| 久久中文字幕人妻熟女| 午夜日韩欧美国产| 亚洲中文av在线| 在线播放国产精品三级| 美女免费视频网站| 久久这里只有精品19| 日韩中文字幕欧美一区二区| 国内毛片毛片毛片毛片毛片| 亚洲人成网站高清观看| 国产蜜桃级精品一区二区三区| 亚洲成人久久爱视频| 又紧又爽又黄一区二区| 欧美激情在线99| 丁香六月欧美| 国产精品一区二区三区四区久久| 欧美日韩亚洲国产一区二区在线观看| 国产精品影院久久| 麻豆成人午夜福利视频| 亚洲人成网站在线播放欧美日韩| 免费大片18禁| 91老司机精品| 亚洲一区二区三区色噜噜| 国产午夜精品久久久久久| 黄色片一级片一级黄色片| 啦啦啦观看免费观看视频高清| av福利片在线观看| 99精品欧美一区二区三区四区| 欧美色视频一区免费| 国产精品一区二区免费欧美| 午夜福利在线观看免费完整高清在 | 一a级毛片在线观看| 欧美成人性av电影在线观看| 校园春色视频在线观看| 午夜两性在线视频| 成人特级av手机在线观看| x7x7x7水蜜桃| www日本在线高清视频| 成人精品一区二区免费| 久久人人精品亚洲av| cao死你这个sao货| 精品一区二区三区av网在线观看| 午夜福利在线观看免费完整高清在 | 成年版毛片免费区| 听说在线观看完整版免费高清| 午夜精品一区二区三区免费看| 天天躁狠狠躁夜夜躁狠狠躁| 国产成人aa在线观看| 亚洲人成网站在线播放欧美日韩| 日日干狠狠操夜夜爽| 看黄色毛片网站| 黄片小视频在线播放| 久久精品国产99精品国产亚洲性色| 成人av一区二区三区在线看| 麻豆av在线久日| 无限看片的www在线观看| 欧美日韩中文字幕国产精品一区二区三区| 精品国内亚洲2022精品成人| 18禁国产床啪视频网站| 欧美国产日韩亚洲一区| 亚洲片人在线观看| 欧美成人免费av一区二区三区| 欧美+亚洲+日韩+国产| 久久午夜亚洲精品久久| 精品一区二区三区视频在线观看免费| 十八禁人妻一区二区| 免费无遮挡裸体视频| а√天堂www在线а√下载| 欧美高清成人免费视频www| 国产三级黄色录像| 亚洲av成人一区二区三| 免费观看人在逋| 国产午夜精品论理片| 欧美性猛交╳xxx乱大交人| 黑人操中国人逼视频| 最近视频中文字幕2019在线8| 色综合婷婷激情| 噜噜噜噜噜久久久久久91| 国产成人精品久久二区二区免费| 亚洲av日韩精品久久久久久密| 99精品在免费线老司机午夜| 床上黄色一级片| 国产欧美日韩精品一区二区| 国产三级黄色录像| 亚洲成av人片在线播放无| 久久久久久大精品| 中文资源天堂在线| 午夜亚洲福利在线播放| 老熟妇乱子伦视频在线观看| 在线观看美女被高潮喷水网站 | 少妇的丰满在线观看| 2021天堂中文幕一二区在线观| 99久久精品国产亚洲精品| 波多野结衣高清无吗| 国产麻豆成人av免费视频| 18禁观看日本| 成年免费大片在线观看| 国产一区在线观看成人免费| 校园春色视频在线观看| 1024香蕉在线观看| 免费在线观看日本一区| 99在线人妻在线中文字幕| 一级毛片精品| 国产精品av久久久久免费| 亚洲专区字幕在线| 成人午夜高清在线视频| 久久国产精品影院| 国产精华一区二区三区| 99久久久亚洲精品蜜臀av| 亚洲人成电影免费在线| 人人妻人人澡欧美一区二区| 久久久国产成人精品二区| 日本五十路高清| 亚洲成a人片在线一区二区| 精品人妻1区二区| 久久久国产成人免费| 老司机在亚洲福利影院| 美女被艹到高潮喷水动态| 身体一侧抽搐| 午夜福利在线观看吧| 国产精品亚洲一级av第二区| 免费大片18禁| 黑人操中国人逼视频| 少妇裸体淫交视频免费看高清| xxx96com| 亚洲精品国产精品久久久不卡| 国内久久婷婷六月综合欲色啪| 久久99热这里只有精品18| 亚洲人成伊人成综合网2020| 国产99白浆流出| 成人性生交大片免费视频hd| 黑人巨大精品欧美一区二区mp4| 美女cb高潮喷水在线观看 | 三级国产精品欧美在线观看 | 男人的好看免费观看在线视频| 欧美日韩综合久久久久久 | 亚洲国产色片| 国语自产精品视频在线第100页| 亚洲五月天丁香| 中文资源天堂在线| 成人午夜高清在线视频| 国产精品99久久久久久久久| 欧美成人性av电影在线观看| 国产午夜精品论理片| 97超视频在线观看视频| 波多野结衣巨乳人妻| 美女高潮喷水抽搐中文字幕| 免费观看精品视频网站| 岛国视频午夜一区免费看| www.精华液| 在线免费观看的www视频| 国产成人啪精品午夜网站| 丰满人妻一区二区三区视频av | 成熟少妇高潮喷水视频| 大型黄色视频在线免费观看| 首页视频小说图片口味搜索| 97碰自拍视频| 国产精品亚洲av一区麻豆| 国内少妇人妻偷人精品xxx网站 | 国产精品乱码一区二三区的特点| 99久久99久久久精品蜜桃| 国产高清激情床上av| 最新在线观看一区二区三区| 男人舔女人下体高潮全视频| 91麻豆av在线| 国产蜜桃级精品一区二区三区| 精品日产1卡2卡| 久久99热这里只有精品18| 欧美极品一区二区三区四区| 在线视频色国产色| 亚洲天堂国产精品一区在线| 91麻豆精品激情在线观看国产| 久久久精品大字幕| 岛国视频午夜一区免费看| 看黄色毛片网站| 热99在线观看视频| 操出白浆在线播放| 久久精品人妻少妇| 麻豆成人午夜福利视频| 麻豆久久精品国产亚洲av| 国产黄色小视频在线观看| 午夜日韩欧美国产| 亚洲aⅴ乱码一区二区在线播放| 久久精品国产综合久久久| 欧美在线一区亚洲| 岛国在线观看网站| 精品国产乱码久久久久久男人| 亚洲专区字幕在线| a级毛片a级免费在线| 亚洲国产精品合色在线| 亚洲国产精品成人综合色| 亚洲成人免费电影在线观看| 中亚洲国语对白在线视频| 99久久久亚洲精品蜜臀av| 久久九九热精品免费| 一a级毛片在线观看| 国产成人一区二区三区免费视频网站| 国产又色又爽无遮挡免费看| 成熟少妇高潮喷水视频| 免费高清视频大片| 日本一本二区三区精品| 亚洲av日韩精品久久久久久密| 欧美一级a爱片免费观看看| 欧美性猛交╳xxx乱大交人| 在线观看舔阴道视频| 久久久久久久精品吃奶| 婷婷精品国产亚洲av| 90打野战视频偷拍视频| 日本一本二区三区精品| 国产综合懂色| 99国产精品一区二区三区| 天天躁日日操中文字幕| 国产精品一区二区精品视频观看| 国产真实乱freesex| 两个人的视频大全免费| 国产单亲对白刺激| 国产成人系列免费观看| 中文字幕最新亚洲高清| 国产精品,欧美在线| 免费在线观看亚洲国产| ponron亚洲| 久99久视频精品免费| 亚洲电影在线观看av| 美女cb高潮喷水在线观看 | 天天躁日日操中文字幕| 婷婷亚洲欧美| 偷拍熟女少妇极品色| 婷婷精品国产亚洲av| 免费看十八禁软件| 97碰自拍视频| 国产日本99.免费观看| 色综合亚洲欧美另类图片| 国产 一区 欧美 日韩| 成在线人永久免费视频| 91九色精品人成在线观看| 我要搜黄色片| 免费在线观看日本一区| 综合色av麻豆| 国产伦一二天堂av在线观看| 欧美成人性av电影在线观看| 麻豆成人av在线观看| 丝袜人妻中文字幕| 最好的美女福利视频网| xxx96com| 日韩欧美三级三区| 特级一级黄色大片| 黄片大片在线免费观看| 国产一区二区三区视频了| 老司机午夜福利在线观看视频| 午夜福利高清视频| 中出人妻视频一区二区| 亚洲国产欧洲综合997久久,| 国产精品电影一区二区三区| 91久久精品国产一区二区成人 | av欧美777| 中文字幕最新亚洲高清| 一二三四社区在线视频社区8| 国产在线精品亚洲第一网站| 一夜夜www| 欧美日韩黄片免| 嫁个100分男人电影在线观看| 99国产极品粉嫩在线观看| 国内久久婷婷六月综合欲色啪| 97超视频在线观看视频| 不卡av一区二区三区| 黄色视频,在线免费观看| 亚洲中文字幕日韩| 久久精品国产99精品国产亚洲性色| 男人舔女人下体高潮全视频| 成人精品一区二区免费| 国产精品九九99| 少妇熟女aⅴ在线视频| 成年人黄色毛片网站| 国产乱人视频| 国产人伦9x9x在线观看| 国产亚洲欧美在线一区二区| 黑人操中国人逼视频| 国产精华一区二区三区| 欧美激情久久久久久爽电影| 别揉我奶头~嗯~啊~动态视频| 美女免费视频网站| 99久久综合精品五月天人人| 国产一区二区在线av高清观看| 欧美乱码精品一区二区三区| 嫁个100分男人电影在线观看| 日韩欧美 国产精品| 欧美中文综合在线视频| 国产成年人精品一区二区| 欧美日韩瑟瑟在线播放| 俺也久久电影网| 亚洲国产精品999在线| 少妇的丰满在线观看| 国产成人福利小说| 99久久99久久久精品蜜桃| 毛片女人毛片| 又粗又爽又猛毛片免费看| 成年免费大片在线观看| 日韩欧美国产一区二区入口| 九色国产91popny在线| 亚洲欧美一区二区三区黑人| 亚洲人成网站在线播放欧美日韩| 亚洲精品456在线播放app | 成年人黄色毛片网站| www日本黄色视频网| 老司机福利观看| 视频区欧美日本亚洲| 成人高潮视频无遮挡免费网站| 18禁黄网站禁片午夜丰满| 国产精品久久久久久久电影 | 国产黄色小视频在线观看| 人妻夜夜爽99麻豆av| 国产麻豆成人av免费视频| 国产精品免费一区二区三区在线| 日本三级黄在线观看| 国产午夜精品久久久久久| 韩国av一区二区三区四区| 99在线视频只有这里精品首页| 国产高潮美女av| 精品久久久久久成人av| www.999成人在线观看| 激情在线观看视频在线高清| 色播亚洲综合网| 免费av不卡在线播放| 少妇的丰满在线观看| 亚洲片人在线观看| 久久久久久久久免费视频了| 女人高潮潮喷娇喘18禁视频| 欧美最黄视频在线播放免费|