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

    基于線粒體D-loop 區(qū)的廣東沿岸和三沙海域綠海龜遺傳多樣性分析

    2022-08-18 04:23:40陳華靈葉明彬謝子強肖寶華廖寶林王付民王少鋒姜怡薇
    關(guān)鍵詞:惠東海龜種群

    陳華靈,葉明彬,謝子強,肖寶華,廖寶林,王付民,王少鋒,姜怡薇

    (1.廣東惠東海龜國家級自然保護區(qū)管理局,廣東惠州 516359;2.廣東海洋大學(xué)深圳研究院,廣東深圳 518108;3.深圳市碧海藍天海洋科技有限公司,廣東深圳 518108)

    綠海龜Chelonia mydas是一種古老的海洋生物,分布在熱帶和亞熱帶海域,是最常見的海龜種類[1]。研究指出,在自然環(huán)境中綠海龜具有強烈的產(chǎn)卵洄游特性,成年雌性海龜會準確地找回其出生地并進行繁殖產(chǎn)卵活動[2-3]。隨著時間的推移,種群之間會產(chǎn)生強烈的遺傳差異,尤其是在母系遺傳的線粒體基因(mtDNA)上[4]。隨著分子標記技術(shù)的發(fā)展,分子標記技術(shù)廣泛應(yīng)用于海洋生物的種群遺傳結(jié)構(gòu)和遺傳多樣性分析[5-6]。國外學(xué)者也通過相關(guān)分子標記技術(shù)開展了綠海龜種群結(jié)構(gòu)、聯(lián)通性和區(qū)域系統(tǒng)地理學(xué)等研究。20世紀90 年代,BOWEN,et al[2]通過研究綠海龜線粒體基因型提出,綠海龜種群有2 個主要的地理種群,一個為大西洋-地中海地理種群,另一個為印度-太平洋地理種群。此后,許多學(xué)者對大西洋、地中海、印度洋、西太平洋、中太平洋和東太平洋等相對較小距離上的綠海龜?shù)牡乩碜V系和種群結(jié)構(gòu)進行研究,并表明在不同海區(qū)的綠海龜種群中,有限的基因流動與偶然的長距離擴散事件是一種常見的現(xiàn)象[7-10]。

    綠海龜是在我國分布的5 種海龜中種群數(shù)量最多的物種,也是唯一一種在我國華南沿岸筑巢產(chǎn)卵的海龜,從我國山東到南海海域均有分布,現(xiàn)為我國一級保護動物[11-13]。由于受到漁業(yè)捕獲、非法貿(mào)易、產(chǎn)卵地生境破壞和海洋污染等威脅,我國綠海龜種群在過去幾十年中嚴重衰退[11]。因此了解我國綠海龜種質(zhì)資源和遺傳多樣性現(xiàn)狀是開展綠海龜種群保護和開展人工繁育的必要理論基礎(chǔ)。近年來,國內(nèi)包括香港、臺灣地區(qū)的學(xué)者紛紛對我國綠海龜種群結(jié)構(gòu)及其遺傳多樣性開展調(diào)查研究。楊文嘉[14]利用D-loop 在內(nèi)的多種分子標記對中國南海海域的175 只綠海龜遺傳多樣性進行了遺傳多樣性研究,初步探究了我國南海綠海龜種群遺傳結(jié)構(gòu)。魏文芝[11]采取問卷調(diào)查、實地調(diào)查和線粒體DNA 分子標記技術(shù)對中國南海綠海龜資源進行研究,結(jié)果顯示我國南海海域的綠海龜種群具有較高水平的遺傳多樣性。香港城市大學(xué)NG,et al[13]采集來自香港、廣東和臺灣海域的110 只綠海龜樣品,獲取了線粒體DNA 控制區(qū)的基因序列和單倍型信息,對其繁殖種群的遺傳多樣性進行了探究,并與整個東太平洋和東南亞海域的綠海龜種群聯(lián)系進行了分析和討論。

    廣東惠東海龜國家級自然保護區(qū)是目前我國大陸近岸唯一的綠海龜產(chǎn)卵地,自20 世紀80 年代成立保護區(qū)以來,開展了大量的綠海龜種質(zhì)資源保護工作,并在國內(nèi)首次人工繁育出綠海龜子代[15-16]。近10 年來,國內(nèi)雖然已有基于分子標記的南海綠海龜生物多樣性研究報道[13-14],但僅集中于較小范圍內(nèi)的地理種群研究。本研究利用D-loop 序列部分片段作為分子標記,對廣東惠東海龜國家級自然保護區(qū)、粵西海域和海南三沙海域綠海龜群體的遺傳多樣性進行比對,為探究我國綠海龜種群遺傳多樣性和開展綠海龜種質(zhì)資源保護奠定理論基礎(chǔ)。

    1 材料與方法

    1.1 樣本采集與基因組總DNA 的提取

    本研究所采集的綠海龜?shù)慕M織樣本均來自廣東惠東海龜國家級自然保護區(qū)救護或洄游產(chǎn)卵孵化的野生綠海龜,其中包括從海南三沙(七連嶼)海域救護綠海龜樣本25 個(編號S)、廣東粵西(硇洲島、東海島至徐聞一帶)海域救護的綠海龜樣本26 個(編號I)、廣東惠東海龜國家級自然保護區(qū)(巽寮灣、平海和港口)海域救護或洄游綠海龜產(chǎn)卵孵化的綠海龜樣本81 個(編號X),共計132 個樣本。采集樣品的綠海龜均為活體,剪取綠海龜前肢量上皮組織并保存在無水乙醇中,用安多福溶液對傷口消毒處理。記錄樣品編號及對應(yīng)的電子芯片號,用北京擎科生物科技有限公司生產(chǎn)的動物基因組DNA 提取試劑盒(TSP201-200)提取海龜DNA,并置于-20 ℃保存?zhèn)溆谩?/p>

    1.2 PCR 擴增與測序分析

    D-loop 基因的PCR 引物參照NORMAN,et al[17]的研究。上下游引物序列分別為TCR-5:TTGTACATC TACTTATTTACCAC(5′-3′)和TCR-6:GTAAGTAAAACTACCGTATGCCAGGTTA(5′-3′)。引物由北京擎科生物科技有限公司合成。PCR 反應(yīng)體系的反應(yīng)總體積50 μL:1.1×T3 Super PCR Mix(擎科生物)47 μL,上下游引物各1 μL,DNA 模板1 μL。PCR 反應(yīng)條件:98 ℃預(yù)變性2 min,98 ℃變性10 s,55 ℃退火10 s,72 ℃延伸20 s,共30 個循環(huán),最后72 ℃延伸5 min。將擴增好的PCR 產(chǎn)物進行瓊脂糖凝膠電泳(2 μL 樣品+6 μL溴酚藍),300 V 電壓下12 min,獲取鑒定膠圖,通過膠圖確定擴增條件是否單一,是否彌散,有無非特異性條帶。選取擴增質(zhì)量較高的PCR 產(chǎn)物送至北京擎科生物科技有限公司進行測序,獲得測序峰圖和堿基序列。

    1.3 序列分析

    測序獲得的堿基序列使用MEGA 6.0 軟件統(tǒng)計序列的堿基組成與含量、堿基信息位點數(shù)量,計算遺傳距離,并構(gòu)建基于K2-P 的鄰位連接(neighbor-joining,NJ)系統(tǒng)發(fā)育樹。采用DnaSP 5.10 軟件統(tǒng)計變異位點位置,單倍型分布等,采用Arlequin 3.1 軟件計算Fst值、Nm值、Tajima′s D 和Fu′s Fs中性檢驗以及進行AMOVA 分析。應(yīng)用DNAsp 軟件與PopART 1.7 軟件單倍型中介連接網(wǎng)絡(luò)(median-joining)分析。

    2 結(jié)果與分析

    2.1 序列多態(tài)性分析

    本研究擴增獲得了132 條D-loop 基因序列,排序去除兩端冗余序列后得到堿基片段在336~388 bp之間,絕大部分的堿基片段在384 bp 以上。在這132 條D-loop 基因序列中T、C、A 和G 堿基的平均含量分別為35.6%、14.4%、32.6%、17.4%,其中A+T 含量(68.2%)明顯高于C+G 含量(31.8%),堿基組成具有明顯的偏向性(表1)。

    表1 3 個地理群體D-loop 基因片段的堿基組成比例Tab.1 The nucleotide percentage composition of the D-loop fragments in 3 populations

    2.2 群體遺傳多樣性及單倍型分析

    在132 條D-loop 基因序列中,共檢測出32 個變異位點,總單倍型多樣性為0.779,總核苷酸多樣性為0.007 55;各群體的單倍型多樣性(Hd)在0.157~0.775 之間,核苷酸多樣性(π)在0.000 63~0.010 17 之間(表2)。根據(jù)GRANT,et al[18]的研究報道,種群單倍型多樣性及核苷酸多樣性臨界值分別為0.5 和0.005,二者的值越大,群體的多樣性程度越高。本研究的三沙地理群體(S)內(nèi)的單倍型多樣性和核苷酸多樣性均未超過臨界值,表明其群體內(nèi)的遺傳多樣性程度較低;粵西地理群體(I)的單倍型多樣性大于0.5,而核苷酸多樣性小于0.005,表明其群體內(nèi)遺傳多樣性程度一般;廣東惠東國家級海龜保護區(qū)群體(X)的單倍型多樣性和核苷酸多樣性均超過臨界值,表明其群體內(nèi)的遺傳多樣性程度較高。以132 個樣本為1 個總?cè)后w來算,其種群單倍型多樣性和核苷酸多樣性均超臨界值,總?cè)后w的遺傳多樣性程度較高。

    表2 3 個綠海龜群體D-loop 基因的遺傳多樣性指數(shù)Tab.2 Genetic diversity index based on D-loop sequences of 3 populations in C.mydas

    在132 條D-loop 基因序列中共檢測出19 個單倍型,其中共享單倍型有7 個(Hap_1、Hap_3、Hap_4、Hap_5、Hap_6、Hap_7 和Hap_9)。從單倍型網(wǎng)絡(luò)圖(圖1)可以看出,Hap_3 和Hap_5 作為中心向其他單倍型發(fā)散,其單倍型出現(xiàn)的頻率也是最高的2 個,分別為31.06%和33.33%。各單倍型之間的突變步數(shù)1~14 步不等。三沙地理群體(S)包含了Hap_4、Hap_5 和Hap_10 這3 種單倍型,粵西地理群體(I)包含了Hap_1、Hap_3、Hap_4、Hap_5、Hap_6、Hap_7、Hap_8 和Hap_9 這8 種單倍型,廣東惠東國家級海龜保護區(qū)群體(X)包含了除Hap_4、Hap_8 和Hap_10 之外的所有16 種單倍型。通過在NCBI(https://blast.ncbi.nlm.nih.gov)上進行的序列比對,并查閱DUTTON,et al[19]、JENSEN,et al[20]和楊文嘉[14]有關(guān)印度-太平洋地區(qū)綠海龜單倍型分類研究的結(jié)果,本研究獲得的單倍型對應(yīng)的相同的序列及對應(yīng)的單倍型如表3 所示。其中Hap_10、Hap_14 以及Hap_19 這3 個單倍型在NCBI 上比對不到完全相同的序列。Hap_18 的單倍型序列與楊文嘉的研究報道中CmC8 一致,其余的15 種單倍型序列信息與DUTTON,et al[19]和JENSEN,et al[20]的研究報道中的相關(guān)15 種單倍型信息一致。

    表3 19 個綠海龜單倍型與其他學(xué)者研究對比情況Tab.3 The 19 haplotypes of C.mydas were compared with other studies

    圖1 基于D-loop 基因序列的3 個綠海龜群體單倍型網(wǎng)絡(luò)中介圖Fig.1 Median-joining network based on D-loop sequences of 3 populations in C.mydas

    2.3 遺傳分化和基因流

    從群體間的遺傳分化結(jié)果(表4)可以看出,S 和I 群體間的Fst值最大,S 和X 群體間的Fst值其次、且大于0.15,I 和X 間的Fst值最小、且小于0.05。群體遺傳分化系數(shù)(Fst)是一個衡量群體間遺傳分化程度的參數(shù):Fst<0.05 表示分化較小,0.05<Fst<0.15 表示中等分化,0.15<Fst<0.25 表示高度分化,0.25<Fst<1.00 表示極高度分化[21]。基因流Nm值顯示S、I、X 這3 個群體之間存在基因交流現(xiàn)象。從3 個群體的種群AMOVA 分子變異結(jié)果分析,群體的遺傳變異主要來自群體內(nèi),占比達到88.3%,群體內(nèi)的占比僅為11.7%(表5)。本研究3 個群體間存在一定的分化,三沙群體和其他2 個群體存在高度的分化,粵西群體和海龜保護區(qū)群體存在較小的分化。

    表4 基于D-loop 基因的群體間Fst(對角線下)和Nm 值(對角線上)Tab.4 The Fst and Nm value based on D-loop sequences

    表5 3 個群體D-loop 序列的種群AMOVA 分析Tab.5 The AMOVA analysis of D-loop sequences in 3 populations

    2.4 基于Tajima′s D 和Fu′s Fs 的綠海龜群體中性檢測分析

    從3 個綠海龜群體的Tajima′sD和Fu′sFs中性檢測分析結(jié)果來看(表6),3 個群體的Tajima′sD和Fu′sFs中性檢測D值均為負值,三沙地理群體(S)的Tajima′sD和Fu′sFs中性檢測P值均小于0.05,粵西地理群體(I)的Tajima′sD的P值大于0.05 而Fu′sFs的P值大于0.05,廣東惠東國家級海龜保護區(qū)群體(X)的Tajima′sD和Fu′sFs中性檢測P值均大于0.05。Tajima′sD和Fu′sFs中性檢測結(jié)果說明三沙群體可能經(jīng)歷過快速種群擴張。

    表6 3 個綠海龜群體的Tajima′s D 和Fu′s Fs 中性檢測值Tab.6 Neutrality tests of 3 populations in C.mydas

    2.5 基于D-loop 序列片段對綠海龜群體和單倍型的系統(tǒng)進化分析

    以棱皮龜Dermochelys coriacea作為外源種(GenBank 序列號:AF121964.1),采用鄰位連接法(neighborjoining,NJ)構(gòu)建3 個綠海龜?shù)乩砣后w共的D-loop 系統(tǒng)發(fā)育樹(圖2)。結(jié)果顯示3 個群體的綠海龜存在較為復(fù)雜的交叉聚類現(xiàn)象。132 個樣本的系統(tǒng)發(fā)育樹主要分為3 個大的分支,其中第1 大分支包含了3 個群體共125 個樣本,第2 大分支包含了廣東惠東國家級海龜保護區(qū)群體(X)的6 個樣本,第3 個分支僅含有廣東惠東國家級海龜保護區(qū)群體(X)的1 個樣本。同時根據(jù)19 個單倍型的D-loop 序列和棱皮龜?shù)男蛄胁捎绵徫贿B接法構(gòu)建了單倍型的系統(tǒng)發(fā)育樹(圖3),單倍型的系統(tǒng)發(fā)育樹的聚類情況也與132 個樣本的系統(tǒng)發(fā)育樹的結(jié)果相互對應(yīng),一共分為3 大分支,其中Hap_11 單獨為一支(分支Ⅲ),Hap_2、Hap_15、Hap_16 和Hap_17 聚為分支II,其余的14 個單倍型聚為分支I。

    圖2 利用132 個樣品的D-loop 序列數(shù)據(jù)構(gòu)建的分子系統(tǒng)進化NJ樹Fig.2 The neighbor-joining tree constructed with D-loop gene sequence of 132 samples

    圖3 利用19 個單倍型D-loop 序列數(shù)據(jù)構(gòu)建的分子系統(tǒng)進化NJ樹Fig.3 The neighbor-joining tree constructed with D-loop gene sequence of 19 haplotypes

    3 討論

    物種群體遺傳多樣性是評估地理區(qū)域內(nèi)物種種質(zhì)資源豐富度和現(xiàn)狀的重要指標,也是開展物種保護和人工繁育的基礎(chǔ)[21-22]。本研究通過檢測群體單倍型數(shù)量、單倍型多樣性和核苷酸多樣性3 個指標來對3個綠海龜?shù)乩砣后w的遺傳多樣性進行評估分析。其中3 個綠海龜?shù)乩砣后w132 個樣本共定義了19 個單倍型,其中三沙、粵西和惠東海龜保護區(qū)的群體單倍型數(shù)量分別為3、7 和16 個,總單倍型多樣性為0.779,總核苷酸多樣性為0.007 55。NG,et al[13]在香港、廣東惠東海龜保護區(qū)和臺灣海域共采集110 只綠海龜并檢測出27 個單倍型,其中廣東惠州海龜國家級自然保護區(qū)檢測出CmP19 和CmP49 2 種單倍型;魏文芝[11]對采集于三沙的51 個野生綠海龜樣品的D-loop 序列進行分析,共定義了14 個單倍型,總單倍型多樣性為0.818,總核苷酸多樣性為0.007 72。楊文嘉[14]對采集于海南島附近的綠海龜樣本進行D-loop 序列擴增,獲得91 個D-loop 序列共定義了8 個單倍型,總單倍型多樣性為0.45,總核苷酸多樣性為0.003 5。JENSEN曾研究了全球4 878 只海龜?shù)倪z傳多樣性,共統(tǒng)計出144 種單倍型,群體平均單倍型多樣性為0.420,平均核苷酸多樣性為0.009 3[20]。綜合國內(nèi)的相關(guān)報道,地處西太平洋的華南沿岸和南海海域的綠海龜群體,單倍型多樣性水平高于全球平均值,具有較高的生物多樣性水平。

    目前,主流的綠海龜系統(tǒng)演化觀點仍遵循著20 世紀末BOWEN,et al[2]提出的2 大地理種群學(xué)說,并由后續(xù)全球各地區(qū)學(xué)者對較小地理區(qū)域的種群結(jié)構(gòu)和譜系進行深入研究。但由于缺乏對D-loop 單倍型命名的標準協(xié)議,全球已報道的綠海龜D-loop 單倍型就有700 多個,其中出現(xiàn)了大量的重復(fù)現(xiàn)象[20]。因此為了方便我們研究對比討論,本研究主要以DUTTON,et al[19]和JENSEN,et al[20]的研究報道中以CmP 命名的印度-太平洋種群的單倍型數(shù)據(jù)作為參考。JENSEN,et al[20]通過收集全球127 個產(chǎn)卵地的4 878 個綠海龜樣本的mtDNA 數(shù)據(jù)集,對種群的聯(lián)通性和基因流動障礙進行了深入的研究,并通過D-loop 序列在2 大地理種群的基礎(chǔ)上細分了11 個綠海龜線粒體DNA 分支譜系。本研究的3 個群體的19 個單倍型,其構(gòu)建的系統(tǒng)進化樹的3 大分支與JENSEN,et al[20]研究成果中的Ⅲ、Ⅶ和Ⅷ這3 個分支譜系是相互契合的。而在JENSEN,et al[20]的研究報道中,本研究對應(yīng)相同的單倍型包含了西太平洋中部、西太平洋南部、日本、印度-太平洋、西南印度洋、西北印度洋等不同地理區(qū)域的綠海龜個體,也從側(cè)面說明華南沿岸和南海海域這個較小地理范圍內(nèi)的綠海龜與整個印度-太平洋大區(qū)的綠海龜具有種群連通性。

    不同群體間的遺傳分化系數(shù)Fst值和基因流Nm值反映了其分化狀況的程度,但國內(nèi)對我國南海海域和華南沿岸兩大區(qū)域內(nèi)的綠海龜遺傳分化方面的研究仍鮮有報道。本研究的3 個地理群體間存在一定的遺傳分化現(xiàn)象,三沙和粵西2 個野生地理群體與廣東惠東海龜國家級自然保護區(qū)群體之間的Fst值分別為0.171 7 和0.043 4,三沙地理群體與粵西地理群體之間的Fst值為0.194 8。南海海域的綠海龜群體與華南沿岸綠海龜群體的Fst值均大于0.15,屬于高度分化的水平。由于綠海龜是一種洄游性海洋生物,在覓食、繁殖和海流等因素作用下,其遷徙活動的范圍十分廣袤。葉明彬等[23]采用衛(wèi)星追蹤的手段對惠東幼年海龜?shù)匿в我?guī)律及覓食地選擇進行研究,發(fā)現(xiàn)在惠東產(chǎn)卵的海龜在東海和南海之間作季節(jié)性定向洄游,覓食活動足跡主要沿我國華南沿岸一直到粵西再到越南沿岸。本研究的遺傳分化結(jié)果也側(cè)面說明了粵西地區(qū)由于地處惠東海龜洄游路線上,其地理群體與惠東海龜保護區(qū)群體存在更近的親緣關(guān)系。

    我國華南沿岸和南海海域是位于整個西太平洋綠海龜?shù)乩砘顒拥慕裹c位置,是西北太平洋、西南太平洋、中太平洋等多個地理譜系綠海龜洄游的交界區(qū)域,從遺傳育種的角度看,在該地區(qū)開展綠海龜?shù)倪z傳雜交育種有利于提升本地區(qū)海龜品種的遺傳多樣性水平。近年來,廣東惠東海龜國家級自然保護區(qū)開展了大量綠海龜人工繁育的技術(shù)研究,并在國內(nèi)首次全人工繁育出綠海龜子代[15-16]。本研究的綠海龜總體樣本具有較高的遺傳多樣性水平,尤其是廣東惠東綠海龜保護區(qū)圈養(yǎng)的綠海龜群體具有豐富的遺傳多樣性,其單倍型數(shù)據(jù)可為該保護區(qū)綠海龜人工選育工作提供指導(dǎo)依據(jù)。

    4 小結(jié)

    本研究對3 個綠海龜?shù)乩砣汗?32 個綠海龜?shù)木€粒體D-loop 基因部分序列進行了擴增、測序及分析,結(jié)果表明:3 個綠海龜?shù)乩砣后w的D-loop 基因片段序列的A+T 含量為68.2%,堿基組成具有明顯的偏向性。132 個樣本共定義了19 個單倍型,單倍型多樣性在0.157~0.775 之間,核苷酸多樣性在0.000 63~0.010 17 之間,本研究的綠海龜總體樣本具有較高的遺傳多樣性。3 個群體間存在一定的分化,其中三沙群體和其他2 個群體存在高度的分化,粵西群體和海龜保護區(qū)圈養(yǎng)群體存在較小的分化,3 個群體之間存在基因交流現(xiàn)象。采用鄰位連接法(NJ)對132 個序列片段和19 個單倍型分別構(gòu)建了系統(tǒng)發(fā)育樹,3 個群體的綠海龜存在較為復(fù)雜的交叉聚類現(xiàn)象,19 個單倍型聚類為3 個大的分支。

    猜你喜歡
    惠東海龜種群
    山西省發(fā)現(xiàn)刺五加種群分布
    海龜
    可怕:一塊塑料便可“殺死”一只海龜
    中華蜂種群急劇萎縮的生態(tài)人類學(xué)探討
    紅土地(2018年7期)2018-09-26 03:07:38
    惠東漁家
    嶺南音樂(2017年2期)2017-05-17 11:42:06
    海龜
    海龜
    惠東港口
    海洋世界(2016年12期)2017-01-03 11:33:00
    惠東漁歌的淵源、藝術(shù)特色及發(fā)展新思路
    樂府新聲(2016年4期)2016-06-22 13:03:14
    崗更湖鯉魚的種群特征
    国产精品福利在线免费观看| 91精品一卡2卡3卡4卡| av在线老鸭窝| 久久人妻av系列| 国产视频内射| 亚洲欧美日韩高清专用| 国产美女午夜福利| 日韩 亚洲 欧美在线| 亚洲不卡免费看| 久久人妻av系列| 啦啦啦观看免费观看视频高清| 亚洲av成人精品一区久久| 国产午夜精品一二区理论片| 亚洲经典国产精华液单| 国内精品一区二区在线观看| 精品国产一区二区三区久久久樱花 | 国产乱来视频区| 最近最新中文字幕大全电影3| 亚洲在线观看片| 中文字幕免费在线视频6| 又黄又爽又刺激的免费视频.| 可以在线观看毛片的网站| 久久精品国产99精品国产亚洲性色| 内地一区二区视频在线| 老女人水多毛片| a级一级毛片免费在线观看| 99久久中文字幕三级久久日本| 精品久久久久久电影网 | 99热这里只有精品一区| 青春草视频在线免费观看| 亚洲国产日韩欧美精品在线观看| 国产精品乱码一区二三区的特点| 欧美三级亚洲精品| 欧美最新免费一区二区三区| 国产午夜福利久久久久久| 精品人妻偷拍中文字幕| 哪个播放器可以免费观看大片| 日本午夜av视频| 亚洲国产欧美在线一区| 深爱激情五月婷婷| 亚洲最大成人手机在线| 韩国高清视频一区二区三区| 亚洲不卡免费看| a级毛片免费高清观看在线播放| 日本午夜av视频| 国产一区亚洲一区在线观看| 我要看日韩黄色一级片| 天天躁日日操中文字幕| 老司机影院成人| 久久久久网色| 能在线免费观看的黄片| 国产91av在线免费观看| 亚洲精品一区蜜桃| 最近最新中文字幕大全电影3| 精品国内亚洲2022精品成人| 天堂av国产一区二区熟女人妻| 尾随美女入室| 免费看光身美女| av天堂中文字幕网| 大香蕉久久网| 在线天堂最新版资源| 22中文网久久字幕| 久久欧美精品欧美久久欧美| 亚洲乱码一区二区免费版| 哪个播放器可以免费观看大片| 亚洲av男天堂| 久久久久久久午夜电影| 美女xxoo啪啪120秒动态图| 最近中文字幕高清免费大全6| 岛国在线免费视频观看| 午夜精品在线福利| 啦啦啦韩国在线观看视频| 国产高潮美女av| 久久精品影院6| 简卡轻食公司| 内射极品少妇av片p| 免费观看a级毛片全部| 国产黄片视频在线免费观看| 国产大屁股一区二区在线视频| 欧美日韩一区二区视频在线观看视频在线 | 国产片特级美女逼逼视频| 中文字幕久久专区| 国产亚洲5aaaaa淫片| 99九九线精品视频在线观看视频| 欧美成人午夜免费资源| 两个人的视频大全免费| 免费看美女性在线毛片视频| 欧美高清性xxxxhd video| 一级黄色大片毛片| 久久久精品94久久精品| 晚上一个人看的免费电影| 大香蕉久久网| 国产成人精品婷婷| 成人性生交大片免费视频hd| 熟妇人妻久久中文字幕3abv| 精品久久久久久久人妻蜜臀av| 亚洲,欧美,日韩| 成人亚洲欧美一区二区av| 老司机影院毛片| 看十八女毛片水多多多| 亚洲国产高清在线一区二区三| 国产精品国产三级专区第一集| 日本wwww免费看| 激情 狠狠 欧美| 亚洲欧美清纯卡通| 久久99热6这里只有精品| 亚洲国产精品国产精品| a级一级毛片免费在线观看| 日韩在线高清观看一区二区三区| 看免费成人av毛片| 人妻夜夜爽99麻豆av| 成人国产麻豆网| 色综合站精品国产| 性插视频无遮挡在线免费观看| 免费黄色在线免费观看| 天堂影院成人在线观看| 国产伦理片在线播放av一区| 欧美日韩一区二区视频在线观看视频在线 | 久久99热6这里只有精品| 国产伦一二天堂av在线观看| 日韩制服骚丝袜av| 人妻少妇偷人精品九色| 女人久久www免费人成看片 | 91狼人影院| 欧美bdsm另类| 简卡轻食公司| 伦理电影大哥的女人| 九九久久精品国产亚洲av麻豆| 男人的好看免费观看在线视频| 青春草视频在线免费观看| 蜜臀久久99精品久久宅男| 最近中文字幕高清免费大全6| av线在线观看网站| 日本一本二区三区精品| 亚洲av中文av极速乱| 少妇裸体淫交视频免费看高清| 欧美极品一区二区三区四区| 我的老师免费观看完整版| 欧美性猛交黑人性爽| 蜜桃亚洲精品一区二区三区| 欧美日韩一区二区视频在线观看视频在线 | 亚洲无线观看免费| 热99在线观看视频| 非洲黑人性xxxx精品又粗又长| 久久鲁丝午夜福利片| 精品一区二区免费观看| 韩国av在线不卡| 国产毛片a区久久久久| 99热网站在线观看| 亚洲成人久久爱视频| 高清视频免费观看一区二区 | 精品久久久久久久末码| 日韩欧美在线乱码| 国产精品一及| 女人被狂操c到高潮| .国产精品久久| 久久精品影院6| 干丝袜人妻中文字幕| 国产精品一区二区在线观看99 | 18禁在线播放成人免费| 免费黄网站久久成人精品| 日本-黄色视频高清免费观看| 丰满乱子伦码专区| 又爽又黄a免费视频| 久久久久久九九精品二区国产| 少妇裸体淫交视频免费看高清| 色尼玛亚洲综合影院| 一区二区三区高清视频在线| 美女内射精品一级片tv| 啦啦啦观看免费观看视频高清| 亚洲乱码一区二区免费版| 午夜福利在线观看吧| 欧美潮喷喷水| 亚洲精品日韩av片在线观看| 久久久国产成人免费| 国产午夜精品论理片| 国产成人91sexporn| 国产人妻一区二区三区在| 日韩国内少妇激情av| 97超视频在线观看视频| 国产精品蜜桃在线观看| 日韩av在线免费看完整版不卡| 日日摸夜夜添夜夜爱| 国产免费又黄又爽又色| 国产 一区精品| 日韩高清综合在线| 两性午夜刺激爽爽歪歪视频在线观看| 又粗又硬又长又爽又黄的视频| 韩国高清视频一区二区三区| 欧美三级亚洲精品| 免费搜索国产男女视频| 色5月婷婷丁香| 久久久久久久久大av| 亚洲成人久久爱视频| 免费无遮挡裸体视频| 日本猛色少妇xxxxx猛交久久| 久久久国产成人精品二区| 欧美区成人在线视频| 激情 狠狠 欧美| 中文字幕人妻熟人妻熟丝袜美| 日韩一本色道免费dvd| av在线播放精品| 亚洲av成人精品一区久久| 亚洲精品,欧美精品| 日韩欧美精品免费久久| 国产精品乱码一区二三区的特点| 变态另类丝袜制服| 国产综合懂色| 五月伊人婷婷丁香| 狂野欧美白嫩少妇大欣赏| 亚洲经典国产精华液单| 亚洲av一区综合| 老女人水多毛片| av播播在线观看一区| 国产亚洲精品久久久com| 国产伦一二天堂av在线观看| 精品久久国产蜜桃| 99久久人妻综合| 午夜日本视频在线| 毛片女人毛片| 看十八女毛片水多多多| 少妇裸体淫交视频免费看高清| 男女那种视频在线观看| 久久精品久久久久久噜噜老黄 | 婷婷色综合大香蕉| 三级国产精品片| 高清午夜精品一区二区三区| 久久6这里有精品| 国产激情偷乱视频一区二区| 色吧在线观看| 欧美另类亚洲清纯唯美| 亚洲自拍偷在线| 一本一本综合久久| 国产毛片a区久久久久| 国产精品日韩av在线免费观看| 亚洲精品乱码久久久v下载方式| 中文字幕av成人在线电影| 一区二区三区免费毛片| 大香蕉久久网| 听说在线观看完整版免费高清| 精品人妻一区二区三区麻豆| www日本黄色视频网| 麻豆成人午夜福利视频| 国产精品1区2区在线观看.| 我要看日韩黄色一级片| 天堂av国产一区二区熟女人妻| 菩萨蛮人人尽说江南好唐韦庄 | 国产精品综合久久久久久久免费| 免费观看的影片在线观看| 不卡视频在线观看欧美| 嫩草影院精品99| 又黄又爽又刺激的免费视频.| 成人三级黄色视频| 国产视频内射| 国产三级中文精品| 日韩亚洲欧美综合| 欧美日韩综合久久久久久| 国产精品精品国产色婷婷| 亚洲人成网站高清观看| 久久99热这里只有精品18| 国产精品一二三区在线看| 免费观看在线日韩| 好男人视频免费观看在线| 国产伦一二天堂av在线观看| or卡值多少钱| 波多野结衣巨乳人妻| 日韩欧美三级三区| 亚洲天堂国产精品一区在线| 久久亚洲精品不卡| 国产精品蜜桃在线观看| 人人妻人人澡欧美一区二区| 亚洲欧美一区二区三区国产| 我的老师免费观看完整版| 亚洲精品影视一区二区三区av| 免费av观看视频| 中文精品一卡2卡3卡4更新| 色吧在线观看| 久久精品国产鲁丝片午夜精品| 男女啪啪激烈高潮av片| 天堂√8在线中文| 国产精品嫩草影院av在线观看| 成人高潮视频无遮挡免费网站| 成人特级av手机在线观看| 久久久久久久亚洲中文字幕| 免费av毛片视频| 免费看日本二区| 日韩人妻高清精品专区| 国产片特级美女逼逼视频| 老司机影院毛片| 久久精品久久精品一区二区三区| 欧美+日韩+精品| 97在线视频观看| 色5月婷婷丁香| 有码 亚洲区| 真实男女啪啪啪动态图| 国产成人免费观看mmmm| 99久久人妻综合| 91av网一区二区| 三级毛片av免费| 国产精品一区二区性色av| 一边摸一边抽搐一进一小说| av天堂中文字幕网| 一二三四中文在线观看免费高清| 伊人久久精品亚洲午夜| 午夜亚洲福利在线播放| 亚洲成av人片在线播放无| 精品久久久久久久久av| 亚洲天堂国产精品一区在线| 菩萨蛮人人尽说江南好唐韦庄 | 国产极品天堂在线| 中文字幕精品亚洲无线码一区| 国产一级毛片七仙女欲春2| 国产精品久久视频播放| .国产精品久久| 欧美zozozo另类| 麻豆成人av视频| 国产成人福利小说| 国产淫语在线视频| 精品久久久久久久久av| 国产伦精品一区二区三区四那| 亚洲av免费在线观看| 亚洲国产精品国产精品| 七月丁香在线播放| 国产又黄又爽又无遮挡在线| 日日撸夜夜添| 丝袜喷水一区| 亚洲av免费在线观看| 日韩欧美国产在线观看| 亚洲一级一片aⅴ在线观看| 国产乱人偷精品视频| 26uuu在线亚洲综合色| 又爽又黄无遮挡网站| 2022亚洲国产成人精品| 国产亚洲一区二区精品| 国产午夜精品久久久久久一区二区三区| 热99re8久久精品国产| 亚洲欧美精品自产自拍| 中文乱码字字幕精品一区二区三区 | 欧美一区二区精品小视频在线| 国产亚洲精品久久久com| 99视频精品全部免费 在线| 18禁在线无遮挡免费观看视频| 国产精品一区二区三区四区免费观看| 一区二区三区四区激情视频| 精品酒店卫生间| av在线播放精品| 免费看光身美女| 日本免费a在线| 久久亚洲国产成人精品v| 丰满乱子伦码专区| 尾随美女入室| 一夜夜www| 亚洲在久久综合| 免费观看人在逋| 国产伦一二天堂av在线观看| 国产黄片视频在线免费观看| www.av在线官网国产| av在线播放精品| 日本熟妇午夜| 国产免费视频播放在线视频 | 精品国产三级普通话版| 国产乱人偷精品视频| 菩萨蛮人人尽说江南好唐韦庄 | 97人妻精品一区二区三区麻豆| 欧美丝袜亚洲另类| 亚洲自拍偷在线| 日本黄色片子视频| 亚洲精品成人久久久久久| 国产精品久久久久久久久免| 久久精品国产自在天天线| 国产精品久久久久久久久免| 免费播放大片免费观看视频在线观看 | 欧美+日韩+精品| 在线观看美女被高潮喷水网站| 国产精品99久久久久久久久| 岛国在线免费视频观看| 波野结衣二区三区在线| 国产在视频线精品| 两个人的视频大全免费| 午夜激情福利司机影院| 日本与韩国留学比较| 麻豆av噜噜一区二区三区| 亚洲真实伦在线观看| 麻豆精品久久久久久蜜桃| 国产三级中文精品| 又爽又黄无遮挡网站| av国产免费在线观看| 亚洲人成网站在线播| 草草在线视频免费看| 在线播放国产精品三级| 亚洲精品色激情综合| 亚洲在线观看片| 一边亲一边摸免费视频| 97超碰精品成人国产| 国产精品精品国产色婷婷| 亚洲,欧美,日韩| 日本爱情动作片www.在线观看| 看免费成人av毛片| 成人无遮挡网站| 国产熟女欧美一区二区| 在现免费观看毛片| 国产乱人视频| 99热这里只有是精品在线观看| 国产午夜精品久久久久久一区二区三区| 国产精品伦人一区二区| 亚洲美女视频黄频| 男女那种视频在线观看| 国国产精品蜜臀av免费| 欧美激情国产日韩精品一区| 免费搜索国产男女视频| 久久精品91蜜桃| 麻豆成人午夜福利视频| 久久99精品国语久久久| 欧美极品一区二区三区四区| 99久久中文字幕三级久久日本| 一级爰片在线观看| 国国产精品蜜臀av免费| 搞女人的毛片| 美女内射精品一级片tv| 大话2 男鬼变身卡| 精品无人区乱码1区二区| 村上凉子中文字幕在线| 亚洲精品乱码久久久v下载方式| 亚洲av男天堂| 国产欧美日韩精品一区二区| 国产精品女同一区二区软件| 黄色配什么色好看| 男插女下体视频免费在线播放| 国产成人精品久久久久久| 简卡轻食公司| 日韩欧美精品v在线| 午夜精品国产一区二区电影 | 全区人妻精品视频| 简卡轻食公司| 亚洲欧美日韩无卡精品| 国产久久久一区二区三区| 99热网站在线观看| av专区在线播放| 欧美一区二区亚洲| 黄色日韩在线| 日韩亚洲欧美综合| 2021天堂中文幕一二区在线观| 亚洲欧美精品综合久久99| 日本wwww免费看| 久久久精品大字幕| av在线天堂中文字幕| 日韩,欧美,国产一区二区三区 | 国内少妇人妻偷人精品xxx网站| 三级国产精品片| 26uuu在线亚洲综合色| 国产大屁股一区二区在线视频| АⅤ资源中文在线天堂| 欧美人与善性xxx| 99久久无色码亚洲精品果冻| 国内精品美女久久久久久| 99热全是精品| 亚洲欧洲国产日韩| 嘟嘟电影网在线观看| 男女下面进入的视频免费午夜| 级片在线观看| 色综合站精品国产| 欧美激情久久久久久爽电影| 99久久精品国产国产毛片| 91午夜精品亚洲一区二区三区| 免费看日本二区| 99热网站在线观看| 日本一本二区三区精品| 中国国产av一级| 色综合色国产| 亚洲四区av| ponron亚洲| 又爽又黄无遮挡网站| 国产精品久久久久久精品电影| 联通29元200g的流量卡| 视频中文字幕在线观看| 在线天堂最新版资源| 国产麻豆成人av免费视频| 美女大奶头视频| 亚洲国产精品国产精品| 久久久久久久午夜电影| 大香蕉久久网| 中文乱码字字幕精品一区二区三区 | 国产不卡一卡二| 观看免费一级毛片| 国产极品精品免费视频能看的| 51国产日韩欧美| av国产免费在线观看| 国产色爽女视频免费观看| 亚洲精品国产成人久久av| 国产精品国产高清国产av| 亚洲四区av| 亚洲真实伦在线观看| 国产久久久一区二区三区| 九色成人免费人妻av| 国产伦在线观看视频一区| 午夜福利在线观看吧| 久久热精品热| 免费黄网站久久成人精品| 久久精品影院6| 国产精品99久久久久久久久| 日本一二三区视频观看| 国产精品福利在线免费观看| 老师上课跳d突然被开到最大视频| 亚洲国产欧美人成| 白带黄色成豆腐渣| 国产精品久久久久久精品电影小说 | 久久久色成人| 成人三级黄色视频| 亚洲在线观看片| 亚洲av一区综合| 你懂的网址亚洲精品在线观看 | 菩萨蛮人人尽说江南好唐韦庄 | 亚洲国产欧美在线一区| 欧美激情国产日韩精品一区| 禁无遮挡网站| 黄片无遮挡物在线观看| 亚洲天堂国产精品一区在线| 中国国产av一级| 精品不卡国产一区二区三区| 免费在线观看成人毛片| 精品久久久久久电影网 | 色播亚洲综合网| 麻豆久久精品国产亚洲av| 成人二区视频| 黑人高潮一二区| 亚洲,欧美,日韩| 国产亚洲精品久久久com| 久久99蜜桃精品久久| 乱码一卡2卡4卡精品| 欧美性猛交╳xxx乱大交人| 精品国产露脸久久av麻豆 | 99久久成人亚洲精品观看| 综合色av麻豆| 在线a可以看的网站| 51国产日韩欧美| 精品久久久久久电影网 | 国产成人a∨麻豆精品| 国产精品综合久久久久久久免费| 久久久久久久国产电影| 变态另类丝袜制服| 网址你懂的国产日韩在线| 欧美xxxx性猛交bbbb| 国产伦精品一区二区三区视频9| 少妇的逼水好多| av卡一久久| 国产乱人视频| 精品人妻一区二区三区麻豆| 18禁在线播放成人免费| 国产一级毛片七仙女欲春2| 99久国产av精品| 久久精品影院6| 岛国毛片在线播放| 最近2019中文字幕mv第一页| 亚洲伊人久久精品综合 | 国产黄a三级三级三级人| 久久久久国产网址| 国产黄色视频一区二区在线观看 | 国产精品野战在线观看| 男女视频在线观看网站免费| 中文字幕av成人在线电影| 直男gayav资源| 国产精品精品国产色婷婷| 精品熟女少妇av免费看| 特大巨黑吊av在线直播| 国产成人福利小说| 黄色日韩在线| 麻豆av噜噜一区二区三区| 国产高清视频在线观看网站| 人妻制服诱惑在线中文字幕| 日韩 亚洲 欧美在线| 久久久成人免费电影| 色5月婷婷丁香| 日韩亚洲欧美综合| 亚洲欧美成人综合另类久久久 | 亚洲精品影视一区二区三区av| 天天躁夜夜躁狠狠久久av| 国产69精品久久久久777片| 国产一区二区在线观看日韩| 色尼玛亚洲综合影院| 日韩三级伦理在线观看| 天堂av国产一区二区熟女人妻| 国产亚洲最大av| 伦精品一区二区三区| 国产成人a∨麻豆精品| 永久网站在线| 午夜激情欧美在线| 五月伊人婷婷丁香| 国产av一区在线观看免费| 色吧在线观看| 天堂av国产一区二区熟女人妻| 在线观看美女被高潮喷水网站| 99久久九九国产精品国产免费| av卡一久久| 国内精品美女久久久久久| 久久精品夜色国产| 婷婷色综合大香蕉| 最近最新中文字幕大全电影3| 精品久久久久久久久av| 国产午夜精品论理片| 亚洲五月天丁香| 日本黄色视频三级网站网址| 国产老妇伦熟女老妇高清| 国产av在哪里看| 久久精品国产亚洲av天美| 国产精品一区二区三区四区久久| 午夜视频国产福利| 欧美性感艳星| 小说图片视频综合网站| 国产人妻一区二区三区在| av播播在线观看一区| 国产av一区在线观看免费| АⅤ资源中文在线天堂| 亚洲无线观看免费| 天天躁日日操中文字幕| 日本免费在线观看一区|