• <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
    崗更湖鯉魚的種群特征
    久久综合国产亚洲精品| 国产成人精品福利久久| 国产亚洲5aaaaa淫片| 七月丁香在线播放| 大香蕉97超碰在线| 久久久久久久国产电影| 久久鲁丝午夜福利片| 汤姆久久久久久久影院中文字幕| 777米奇影视久久| 97热精品久久久久久| 久久精品国产亚洲网站| 亚洲欧美成人精品一区二区| 只有这里有精品99| 18禁在线播放成人免费| 一级黄片播放器| 久久精品熟女亚洲av麻豆精品| 你懂的网址亚洲精品在线观看| 最新中文字幕久久久久| 亚洲aⅴ乱码一区二区在线播放| 亚洲av综合色区一区| 人人妻人人看人人澡| av网站免费在线观看视频| 亚洲精品乱久久久久久| 国内精品宾馆在线| 精品一品国产午夜福利视频| 日日啪夜夜撸| 欧美少妇被猛烈插入视频| 亚洲国产最新在线播放| 久久久久久久亚洲中文字幕| 国产高潮美女av| 国产一区亚洲一区在线观看| a 毛片基地| 成人综合一区亚洲| 亚洲av欧美aⅴ国产| 中文精品一卡2卡3卡4更新| 97热精品久久久久久| 亚洲av在线观看美女高潮| 少妇高潮的动态图| 久久久成人免费电影| 91精品国产国语对白视频| 国产无遮挡羞羞视频在线观看| 亚洲精品一二三| 亚洲精品亚洲一区二区| 国产精品欧美亚洲77777| 成年人午夜在线观看视频| 亚洲国产毛片av蜜桃av| 国产女主播在线喷水免费视频网站| 国产精品人妻久久久影院| 婷婷色综合www| 国产大屁股一区二区在线视频| 网址你懂的国产日韩在线| 3wmmmm亚洲av在线观看| 插阴视频在线观看视频| 欧美一级a爱片免费观看看| 亚洲国产高清在线一区二区三| 日日撸夜夜添| 日本欧美国产在线视频| 成人综合一区亚洲| 只有这里有精品99| 少妇猛男粗大的猛烈进出视频| 狠狠精品人妻久久久久久综合| 不卡视频在线观看欧美| 久久精品久久久久久噜噜老黄| 国产高清有码在线观看视频| 极品教师在线视频| 亚洲欧美日韩东京热| 一区在线观看完整版| 91在线精品国自产拍蜜月| 一级av片app| 伦精品一区二区三区| 亚洲熟女精品中文字幕| 国模一区二区三区四区视频| 日韩视频在线欧美| 欧美亚洲 丝袜 人妻 在线| 国产精品一区二区在线不卡| 国产欧美日韩一区二区三区在线 | 国产欧美日韩一区二区三区在线 | 欧美日本视频| 黄片wwwwww| 亚洲欧美日韩卡通动漫| 丝袜喷水一区| 欧美精品一区二区大全| 偷拍熟女少妇极品色| 日本黄色片子视频| 久久久久久久久久久免费av| 日韩中文字幕视频在线看片 | 国产精品偷伦视频观看了| 国产精品女同一区二区软件| 亚洲不卡免费看| .国产精品久久| 我的女老师完整版在线观看| 国产 精品1| 国产成人午夜福利电影在线观看| 深夜a级毛片| 免费观看性生交大片5| 亚洲国产成人一精品久久久| 身体一侧抽搐| 国产视频首页在线观看| 国产久久久一区二区三区| 黄色日韩在线| 九色成人免费人妻av| 51国产日韩欧美| 中文资源天堂在线| 日本av免费视频播放| 午夜福利在线在线| 久久久久久久精品精品| 中文字幕免费在线视频6| h日本视频在线播放| 日本vs欧美在线观看视频 | 一个人免费看片子| 欧美精品一区二区大全| 免费大片18禁| 欧美高清性xxxxhd video| 成人免费观看视频高清| 欧美精品国产亚洲| 日本黄色片子视频| 夫妻性生交免费视频一级片| 欧美另类一区| 建设人人有责人人尽责人人享有的 | 永久网站在线| 亚洲国产精品一区三区| 精品人妻一区二区三区麻豆| 亚洲激情五月婷婷啪啪| 卡戴珊不雅视频在线播放| 国产久久久一区二区三区| 亚洲成色77777| 免费观看a级毛片全部| 亚洲欧美日韩无卡精品| 99精国产麻豆久久婷婷| 亚洲av成人精品一二三区| 青青草视频在线视频观看| 97热精品久久久久久| 日本黄色片子视频| 成人综合一区亚洲| h日本视频在线播放| 中文欧美无线码| 一边亲一边摸免费视频| 欧美xxxx黑人xx丫x性爽| 欧美3d第一页| 成年女人在线观看亚洲视频| 在线天堂最新版资源| 国产精品av视频在线免费观看| 国产精品爽爽va在线观看网站| 亚洲精品日本国产第一区| 国产午夜精品一二区理论片| 少妇人妻 视频| 国产精品久久久久成人av| av专区在线播放| 久久国产乱子免费精品| 久久久久网色| 岛国毛片在线播放| 久久久精品94久久精品| 久久久久久久久久久丰满| 午夜福利在线观看免费完整高清在| 国产精品三级大全| 国产久久久一区二区三区| 国产精品一区二区性色av| 久久ye,这里只有精品| 欧美日韩亚洲高清精品| 高清黄色对白视频在线免费看 | 一本久久精品| 欧美+日韩+精品| 久久久成人免费电影| 亚洲欧美中文字幕日韩二区| 国产伦在线观看视频一区| 亚洲怡红院男人天堂| 热99国产精品久久久久久7| 欧美老熟妇乱子伦牲交| 国产精品国产三级国产av玫瑰| 乱系列少妇在线播放| 亚洲最大成人中文| 国产亚洲91精品色在线| 美女视频免费永久观看网站| av网站免费在线观看视频| 嫩草影院新地址| 精品亚洲成国产av| 久久99蜜桃精品久久| 精品久久久久久久末码| av专区在线播放| 日韩中字成人| 午夜激情福利司机影院| 久久国产精品大桥未久av | 亚洲欧美日韩另类电影网站 | 久久影院123| 国产爽快片一区二区三区| 日本黄大片高清| 深爱激情五月婷婷| 黄色欧美视频在线观看| 五月伊人婷婷丁香| a 毛片基地| 黄色一级大片看看| 狂野欧美激情性xxxx在线观看| 成年av动漫网址| 中文天堂在线官网| 性色avwww在线观看| 另类亚洲欧美激情| 我的女老师完整版在线观看| 特大巨黑吊av在线直播| 日本-黄色视频高清免费观看| 我要看日韩黄色一级片| 精品一区二区三卡| 五月伊人婷婷丁香| 免费在线观看成人毛片| 成年美女黄网站色视频大全免费 | 777米奇影视久久| 亚洲精品中文字幕在线视频 | 国产精品国产三级专区第一集| 午夜福利视频精品| 午夜视频国产福利| 日日撸夜夜添| 国产爽快片一区二区三区| 男的添女的下面高潮视频| 一二三四中文在线观看免费高清| 成人综合一区亚洲| 99热网站在线观看| 久久久久久久国产电影| 国产极品天堂在线| 日韩亚洲欧美综合| 免费av不卡在线播放| 久久久精品94久久精品| 亚洲av成人精品一区久久| 国产午夜精品久久久久久一区二区三区| 一区二区三区免费毛片| 国产在线视频一区二区| 91久久精品国产一区二区三区| 久久久久久人妻| 中文字幕av成人在线电影| 天天躁日日操中文字幕| 伊人久久精品亚洲午夜| 我要看日韩黄色一级片| 精品一品国产午夜福利视频| 亚洲怡红院男人天堂| 国产乱来视频区| 精品熟女少妇av免费看| 综合色丁香网| 97超碰精品成人国产| 一级毛片 在线播放| 超碰97精品在线观看| 99精国产麻豆久久婷婷| 日韩成人av中文字幕在线观看| 亚洲第一av免费看| 亚洲综合色惰| 熟妇人妻不卡中文字幕| freevideosex欧美| 国产片特级美女逼逼视频| 女人久久www免费人成看片| 精品一区在线观看国产| 国产免费又黄又爽又色| 国产成人91sexporn| 韩国高清视频一区二区三区| 精品99又大又爽又粗少妇毛片| 交换朋友夫妻互换小说| 99久国产av精品国产电影| 日韩,欧美,国产一区二区三区| 色哟哟·www| 亚洲av在线观看美女高潮| 国产精品偷伦视频观看了| 国产精品嫩草影院av在线观看| 精品酒店卫生间| 久久久精品94久久精品| 精品视频人人做人人爽| www.色视频.com| 亚洲精品中文字幕在线视频 | 少妇被粗大猛烈的视频| 亚洲成人一二三区av| 永久免费av网站大全| 99久久综合免费| 久久av网站| 国产成人精品婷婷| 老师上课跳d突然被开到最大视频| 精品人妻视频免费看| 国产精品久久久久久精品电影小说 | 啦啦啦在线观看免费高清www| 亚洲av不卡在线观看| 国内少妇人妻偷人精品xxx网站| 欧美xxxx黑人xx丫x性爽| 久久精品国产亚洲av涩爱| 亚洲内射少妇av| 极品少妇高潮喷水抽搐| 国产伦精品一区二区三区视频9| 少妇熟女欧美另类| av网站免费在线观看视频| 18禁裸乳无遮挡免费网站照片| 亚洲在久久综合| 91久久精品国产一区二区三区| 国产成人a区在线观看| 日韩不卡一区二区三区视频在线| 成年免费大片在线观看| 国产黄色免费在线视频| 乱系列少妇在线播放| 欧美日韩一区二区视频在线观看视频在线| 王馨瑶露胸无遮挡在线观看| 久久久久久久久久成人| 久久人人爽av亚洲精品天堂 | 久久6这里有精品| 亚洲精品一二三| 九色成人免费人妻av| 精品亚洲成国产av| 色吧在线观看| 新久久久久国产一级毛片| 久久精品国产自在天天线| 免费在线观看成人毛片| 国产精品国产三级专区第一集| 成年人午夜在线观看视频| 精品国产露脸久久av麻豆| 国产免费视频播放在线视频| 99国产精品免费福利视频| 精品一区在线观看国产| 国产精品国产av在线观看| 久久久国产一区二区| 一个人免费看片子| 免费不卡的大黄色大毛片视频在线观看| 亚洲精品国产av蜜桃| 国产黄频视频在线观看| 国产精品国产三级国产专区5o| 22中文网久久字幕| 高清欧美精品videossex| 国产男女超爽视频在线观看| 国产av国产精品国产| 麻豆成人av视频| 亚洲精品国产色婷婷电影| 亚洲三级黄色毛片| 97超视频在线观看视频| 极品少妇高潮喷水抽搐| 日韩大片免费观看网站| 少妇精品久久久久久久| 日日撸夜夜添| 亚洲国产毛片av蜜桃av| 亚洲av福利一区| 国产深夜福利视频在线观看| 男人舔奶头视频| 一区二区三区免费毛片| 国精品久久久久久国模美| 插阴视频在线观看视频| 精品少妇黑人巨大在线播放| 日韩大片免费观看网站| 夜夜骑夜夜射夜夜干| 久久精品国产亚洲av涩爱| 日本猛色少妇xxxxx猛交久久| 嘟嘟电影网在线观看| 熟女电影av网| 五月玫瑰六月丁香| 亚洲国产精品成人久久小说| 欧美一级a爱片免费观看看| 国产精品福利在线免费观看| 国模一区二区三区四区视频| 亚洲三级黄色毛片| 一区二区三区精品91| 日本av手机在线免费观看| 欧美成人精品欧美一级黄| 国产精品久久久久久久电影| 免费观看无遮挡的男女| 亚洲人成网站高清观看| 免费人妻精品一区二区三区视频| 91精品国产国语对白视频| 观看免费一级毛片| 亚洲av中文av极速乱| 国产精品久久久久久av不卡| 联通29元200g的流量卡| 国产一区二区三区av在线| 亚洲四区av| 一个人看视频在线观看www免费| 亚州av有码| 国产精品久久久久久精品古装| 国产有黄有色有爽视频| 亚洲国产精品一区三区| 老师上课跳d突然被开到最大视频| 亚洲精品久久久久久婷婷小说| av国产久精品久网站免费入址| 亚洲欧美成人精品一区二区| 97超视频在线观看视频| 欧美 日韩 精品 国产| 日韩av不卡免费在线播放| 黄色视频在线播放观看不卡| 久久热精品热| 国产视频内射| av免费观看日本| 女的被弄到高潮叫床怎么办| 最黄视频免费看| 国产精品蜜桃在线观看| 精品久久久久久电影网| 一级爰片在线观看| 麻豆国产97在线/欧美| 成人一区二区视频在线观看| 国产精品三级大全| 久久久欧美国产精品| 在线观看国产h片| 国产毛片在线视频| 免费观看无遮挡的男女| 欧美日韩综合久久久久久| 国产免费一区二区三区四区乱码| 国产av精品麻豆| 日本一二三区视频观看| 亚洲精品aⅴ在线观看| 国产精品99久久99久久久不卡 | 欧美日本视频| av在线播放精品| 伦理电影大哥的女人| 卡戴珊不雅视频在线播放| 人人妻人人澡人人爽人人夜夜| 欧美xxxx黑人xx丫x性爽| 伦理电影大哥的女人| 在线观看国产h片| 久热久热在线精品观看| 女的被弄到高潮叫床怎么办| 一级毛片黄色毛片免费观看视频| 国产精品一区二区在线不卡| 黄色欧美视频在线观看| 在线观看一区二区三区| 午夜激情福利司机影院| 国产有黄有色有爽视频| 校园人妻丝袜中文字幕| 久久国产精品男人的天堂亚洲 | 一本色道久久久久久精品综合| 校园人妻丝袜中文字幕| 久久这里有精品视频免费| av免费在线看不卡| 亚洲av不卡在线观看| 国产极品天堂在线| 热99国产精品久久久久久7| 伦理电影大哥的女人| 成人18禁高潮啪啪吃奶动态图 | 九九久久精品国产亚洲av麻豆| 熟女人妻精品中文字幕| 大陆偷拍与自拍| 高清不卡的av网站| 高清午夜精品一区二区三区| 亚洲精品日本国产第一区| 91午夜精品亚洲一区二区三区| 韩国av在线不卡| 少妇人妻一区二区三区视频| 国产久久久一区二区三区| 18禁裸乳无遮挡动漫免费视频| 免费人成在线观看视频色| 成人影院久久| 黄色配什么色好看| 亚洲va在线va天堂va国产| 永久免费av网站大全| 久久精品国产亚洲网站| av免费观看日本| 久久精品国产鲁丝片午夜精品| 亚洲精品色激情综合| h视频一区二区三区| 欧美性感艳星| 三级经典国产精品| 五月天丁香电影| 亚洲av国产av综合av卡| 极品教师在线视频| 成人国产av品久久久| 三级国产精品片| 精华霜和精华液先用哪个| 亚洲av在线观看美女高潮| 五月玫瑰六月丁香| 中文资源天堂在线| 大片免费播放器 马上看| 制服丝袜香蕉在线| 久久99热6这里只有精品| 最新中文字幕久久久久| 久久国产精品男人的天堂亚洲 | h日本视频在线播放| 我要看黄色一级片免费的| 亚洲欧美日韩东京热| 久久人人爽人人爽人人片va| 老女人水多毛片| 免费av不卡在线播放| 久久久久久伊人网av| 最近最新中文字幕大全电影3| 又大又黄又爽视频免费| 中文字幕亚洲精品专区| 看非洲黑人一级黄片| 狠狠精品人妻久久久久久综合| 亚洲av成人精品一二三区| 亚洲成人一二三区av| 国产 一区精品| 91精品国产九色| 精品国产乱码久久久久久小说| 免费看日本二区| 深夜a级毛片| 国产精品久久久久久精品古装| 国产黄频视频在线观看| 色视频www国产| 久久99蜜桃精品久久| 性色av一级| 精品久久久久久久久亚洲| 欧美成人精品欧美一级黄| av播播在线观看一区| 在线免费观看不下载黄p国产| freevideosex欧美| 秋霞在线观看毛片| 中文字幕av成人在线电影| 日韩亚洲欧美综合| 久久国产精品大桥未久av | 噜噜噜噜噜久久久久久91| 日韩一区二区视频免费看| 又粗又硬又长又爽又黄的视频| 久久久久久久亚洲中文字幕| 我的女老师完整版在线观看| 97在线人人人人妻| 亚洲av男天堂| 涩涩av久久男人的天堂| 色婷婷av一区二区三区视频| av黄色大香蕉| 最近最新中文字幕免费大全7| 国产亚洲av片在线观看秒播厂| 亚洲av.av天堂| 精品亚洲成国产av| 18禁裸乳无遮挡动漫免费视频| 免费黄频网站在线观看国产| 成人国产av品久久久| 亚洲国产毛片av蜜桃av| 国产精品国产三级国产av玫瑰| 乱码一卡2卡4卡精品| 联通29元200g的流量卡| 国产亚洲91精品色在线| 一级毛片aaaaaa免费看小| 久久ye,这里只有精品| 一级爰片在线观看| 亚洲内射少妇av| 亚洲国产精品999| 国产高清有码在线观看视频| 免费大片黄手机在线观看| 91精品国产国语对白视频| 在线精品无人区一区二区三 | 中文字幕久久专区| 亚洲国产欧美人成| 日韩成人伦理影院| 青青草视频在线视频观看| 黄色一级大片看看| 我要看日韩黄色一级片| 久久久欧美国产精品| 一区二区三区四区激情视频| 99精国产麻豆久久婷婷| 男人和女人高潮做爰伦理| 高清黄色对白视频在线免费看 | 国产黄色视频一区二区在线观看| 高清黄色对白视频在线免费看 | 性色av一级| 哪个播放器可以免费观看大片| 免费看日本二区| 男女下面进入的视频免费午夜| 99久久人妻综合| 亚洲综合色惰| 亚洲第一av免费看| 国产精品免费大片| 18+在线观看网站| 深爱激情五月婷婷| 久久鲁丝午夜福利片| 在线观看免费高清a一片| .国产精品久久| 国内精品宾馆在线| 国内少妇人妻偷人精品xxx网站| 免费在线观看成人毛片| 国产一区二区三区av在线| 国产爽快片一区二区三区| 久久精品久久久久久噜噜老黄| 日韩视频在线欧美| 日日摸夜夜添夜夜爱| 中文字幕av成人在线电影| 男人狂女人下面高潮的视频| 亚洲欧美成人综合另类久久久| 亚洲精品国产av蜜桃| 熟女av电影| av一本久久久久| 久久鲁丝午夜福利片| 日本爱情动作片www.在线观看| 国产成人精品福利久久| 天天躁日日操中文字幕| 欧美三级亚洲精品| 精品人妻熟女av久视频| 欧美精品一区二区大全| 超碰97精品在线观看| 国产午夜精品久久久久久一区二区三区| 一级黄片播放器| 男女下面进入的视频免费午夜| 精品人妻一区二区三区麻豆| 国产成人精品婷婷| 国产又色又爽无遮挡免| 嫩草影院新地址| av福利片在线观看| 精品一区二区三区视频在线| 亚洲真实伦在线观看| 观看免费一级毛片| 插逼视频在线观看| 22中文网久久字幕| 在线亚洲精品国产二区图片欧美 | 国产免费又黄又爽又色| 成人影院久久| 一个人看的www免费观看视频| 在线免费十八禁| 蜜臀久久99精品久久宅男| 国产91av在线免费观看| 色婷婷av一区二区三区视频| 深爱激情五月婷婷| 高清午夜精品一区二区三区| 久久久久性生活片| 九九久久精品国产亚洲av麻豆| 亚洲三级黄色毛片| 午夜免费鲁丝| 日本免费在线观看一区| 日本黄色片子视频| 日日撸夜夜添| 日日摸夜夜添夜夜添av毛片| 午夜福利视频精品| 22中文网久久字幕| 成年av动漫网址| 久久人人爽人人爽人人片va| 一个人看视频在线观看www免费| 最近手机中文字幕大全| 欧美亚洲 丝袜 人妻 在线| 欧美极品一区二区三区四区| 五月玫瑰六月丁香| 十八禁网站网址无遮挡 | 99热这里只有是精品在线观看| 菩萨蛮人人尽说江南好唐韦庄| 精品视频人人做人人爽|