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

    基于COI和16S rRNA基因片段鑒定廈門海域的仔稚魚

    2021-11-04 07:41:56徐春燕莊之棟沈長(zhǎng)春蔡建堤謝少卿
    漁業(yè)研究 2021年5期

    徐春燕,莊之棟,馬 超,劉 勇,沈長(zhǎng)春,蔡建堤,謝少卿

    (福建省水產(chǎn)研究所,福建省海洋生物增養(yǎng)殖與高值化利用重點(diǎn)實(shí)驗(yàn)室,福建 廈門 361013)

    魚類早期個(gè)體(包括魚卵、仔稚魚)發(fā)育過(guò)程是魚類生命周期中非常短暫、形態(tài)學(xué)和生理學(xué)等特征變化明顯的過(guò)渡時(shí)期[1],魚卵和仔稚魚調(diào)查研究是魚類生態(tài)學(xué)以及魚類生活史研究的一部分,同時(shí)也是魚類資源調(diào)查、早期補(bǔ)充機(jī)制和種群數(shù)量變動(dòng)、漁業(yè)資源可持續(xù)利用等研究的重要基礎(chǔ)[1-2]。目前,魚卵和仔稚魚種類的鑒定多是根據(jù)形態(tài)特征進(jìn)行區(qū)別,但由于個(gè)體較小、用于鑒別的形態(tài)特征較少,不同種類很難進(jìn)行準(zhǔn)確區(qū)分。與已知的成魚種類相比,已鑒別的魚卵和仔稚魚種類少之又少,魚類早期生活史的相關(guān)研究相對(duì)滯后,相關(guān)形態(tài)學(xué)資料和圖集十分有限,嚴(yán)重限制了鑒定工作的開展;另外國(guó)內(nèi)缺少相關(guān)的分類學(xué)家,鑒定者的水平不一,很難保證鑒定結(jié)果的準(zhǔn)確性[3-5]。

    2003年,加拿大分類學(xué)家Hebert P D N[6]首次提出DNA條形碼的概念,此后,DNA條形碼技術(shù)被廣泛應(yīng)用于各種動(dòng)物、植物及海洋生物的種類鑒定[7-9]。DNA條形碼技術(shù)利用DNA序列對(duì)物種進(jìn)行鑒定,突破了傳統(tǒng)分類學(xué)只能依靠外部形態(tài)特征鑒定物種的限制,對(duì)產(chǎn)卵場(chǎng)相互交疊和產(chǎn)卵期相互交錯(cuò)、形態(tài)特征近似、形態(tài)學(xué)分類困難的魚類早期個(gè)體種類的分類以及形態(tài)學(xué)分類上存有疑問(wèn)的種類的驗(yàn)證提供了新的方法和手段[1]。線粒體細(xì)胞色素C氧化酶亞基I(Cytochrome C oxidase subunit I,COI)基因是魚類常用的條形碼基因,被廣泛應(yīng)用于物種鑒定和多樣性分析[7,10-12]。然而,目前有研究指出單一的COI基因并不適用于所有魚類物種的鑒別[13],結(jié)合其他條形碼基因相互參照和驗(yàn)證是確保條形碼技術(shù)鑒定準(zhǔn)確性的必要手段[7,14-15]。

    本研究選取應(yīng)用比較廣泛的COI和16S rRNA基因,對(duì)廈門海域采集的仔稚魚樣品進(jìn)行種類鑒定,通過(guò)分析遺傳距離及系統(tǒng)進(jìn)化關(guān)系,探究其在仔稚魚種類鑒定中的適用性,同時(shí)可以為仔稚魚的準(zhǔn)確鑒定和分類學(xué)研究提供依據(jù)和基礎(chǔ)資料。

    1 材料與方法

    1.1 樣品采集與形態(tài)學(xué)鑒定

    在廈門海域設(shè)置12個(gè)采樣站位(圖1),分別于2019年4月22日和2020年9月14日采集仔稚魚樣品。參照GB/T 12763.6—2007《海洋調(diào)查規(guī)范 第6部分:海洋生物調(diào)查》[16],利用大型浮游生物網(wǎng)(網(wǎng)口內(nèi)徑80 cm,網(wǎng)長(zhǎng)280 cm,孔徑0.505 mm)進(jìn)行水平拖網(wǎng)采樣,樣品用95%酒精現(xiàn)場(chǎng)固定,帶回實(shí)驗(yàn)室進(jìn)行形態(tài)學(xué)和分子生物學(xué)鑒定。

    采集到的仔稚魚樣品首先在奧林巴斯SZX7體視顯微鏡下根據(jù)形態(tài)特征進(jìn)行鑒定和形態(tài)學(xué)歸類,參照Okiyama M的《An Atlas of the Early Stage Fishes in Japan》[17]及張仁齋等的《中國(guó)近海魚卵與仔魚》[2]進(jìn)行鑒定,之后根據(jù)初步的鑒定結(jié)果,每一種選取代表性個(gè)體進(jìn)行基于COI和16S rRNA基因的DNA條形碼分析。

    1.2 DNA提取、PCR擴(kuò)增及測(cè)序

    選取的仔稚魚樣品用QIAGEN DNeasy Blood & Tissue Kit試劑盒進(jìn)行DNA提取,DNA樣品置于4℃保存?zhèn)溆?。COI基因序列擴(kuò)增采用通用引物[18],F(xiàn)ish F1:5′-TCAACCAACCACAAAGACATTGGCAC-3′,F(xiàn)ish R1:5′-TAGACTTCTGGGTGGCCAAAGAATCA-3′;16S rRNA基因序列擴(kuò)增采用通用引物[19],16S AR:5′-CGCCTGTTTATCAAAAACAT-3′,16S BR:5′-CCGGTCTGAACTCAGATCACGT-3′。

    PCR反應(yīng)體系:2.5 μL 10×buffer,2 μL dNTPs(2.5 mmol/L),正反引物(10 μmol/L)各1 μL ,0.2 μL Taq酶,DNA模板1 μL,超純水補(bǔ)足至25 μL。反應(yīng)條件:94℃預(yù)變性5 min;94℃變性60 s,50℃退火50 s,72℃延伸60 s,30個(gè)循環(huán);最后72℃延伸5 min。用1%的瓊脂糖凝膠電泳對(duì)擴(kuò)增產(chǎn)物進(jìn)行檢測(cè),選取擴(kuò)增效果好的PCR產(chǎn)物送至上海生工進(jìn)行雙向測(cè)序。

    1.3 序列分析

    測(cè)序獲得的序列首先利用Chromas 2.6.5軟件查看原始序列質(zhì)量,質(zhì)量好的序列利用Seqman軟件進(jìn)行雙向拼接。將拼接所得COI和16S rRNA基因序列在NCBI(National Center for Biotechnology Information)網(wǎng)站進(jìn)行比對(duì),COI基因同時(shí)可以在BOLD(Barcode of Life Database)進(jìn)行比對(duì)。比對(duì)后分別根據(jù)相似度進(jìn)行種類鑒定,與庫(kù)中物種序列相似度≥99%的即鑒定為同一個(gè)種,92%~99%之間的鑒定為同一屬,85%~92%之間的鑒定為同一科[4],另外,兩種基因鑒定結(jié)果一致時(shí),最終將此種鑒定為一個(gè)物種,結(jié)果不一致時(shí),則將其鑒定到科或?qū)佟hb定后相關(guān)序列及信息遞交至GenBank獲取序列號(hào)。兩種基因分別用MEGAX軟件中Clustal W程序進(jìn)行多重比對(duì),比對(duì)后去除兩端冗余序列,保留共有序列區(qū)域進(jìn)行后續(xù)分析,統(tǒng)計(jì)序列的長(zhǎng)度、堿基組成、保守位點(diǎn)、變異位點(diǎn)、簡(jiǎn)約信息位點(diǎn)、單突變位點(diǎn),種內(nèi)和種間遺傳距離基于Kimura-2-parameter雙參數(shù)模型(K2P)進(jìn)行計(jì)算,分子系統(tǒng)進(jìn)化樹利用鄰接法(Neighbor-joining)構(gòu)建,以Bootstrap 1 000次重復(fù)抽樣檢驗(yàn)其分支置信度[20]。

    2 結(jié)果與分析

    2.1 測(cè)序結(jié)果及序列分析

    本研究共選取80尾仔稚魚樣品分別進(jìn)行COI和16S rRNA基因的條形碼分析,由于在DNA提取、PCR擴(kuò)增及測(cè)序過(guò)程中,均有可能出現(xiàn)結(jié)果不佳的情況,最終獲得64條COI基因序列和74條16S rRNA基因序列。COI基因序列平均長(zhǎng)度655 bp,16S rRNA基因序列平均長(zhǎng)度575 bp,序列均遞交至GenBank獲取序列號(hào),序列號(hào)見表1。

    COI基因序列的平均堿基組成為T:30.4%、C:28.0%、A:23.4%、G:18.2%,A+T含量(53.8%)高于G+C含量(46.2%),符合多數(shù)魚類COI序列的堿基組成特征。各密碼子堿基組成表現(xiàn)出明顯的差異性,第2個(gè)密碼子平均GC含量最高,為56.5%;第3密碼子次之,平均為42.7%;第1密碼子平均GC含量最低,為39.5%。全部位點(diǎn)中包括保守位點(diǎn)379個(gè)、變異位點(diǎn)276個(gè)、簡(jiǎn)約信息位點(diǎn)265個(gè)、單突變位點(diǎn)11個(gè)。

    16S rRNA基因序列的平均堿基組成為T:22.9%、C:25.4%、A:28.6%、G:23.1%,A+T含量(51.5%)高于G+C含量(48.5%),與多數(shù)魚類線粒體DNA的堿基組成特征一致。各密碼子堿基組成表現(xiàn)出明顯的差異性,第1個(gè)密碼子平均GC含量最高,為54.8%;第3密碼子次之,平均為45.7%;第2密碼子平均GC含量最低,為44.9%。全部位點(diǎn)中包括保守位點(diǎn)301個(gè)、變異位點(diǎn)298個(gè)、簡(jiǎn)約信息位點(diǎn)278個(gè)、單突變位點(diǎn)20個(gè)。

    2.2 序列比對(duì)及仔稚魚種類鑒定

    所得COI基因同時(shí)在NCBI和BOLD數(shù)據(jù)庫(kù)進(jìn)行比對(duì),16S rRNA基因在NCBI比對(duì),比對(duì)后依據(jù)相似度進(jìn)行鑒定,鑒定結(jié)果見表1。COI基因?qū)⒆兄婶~樣品鑒定為26個(gè)種類,其中19個(gè)種類鑒定到種、6個(gè)鑒定到屬、1個(gè)種類僅鑒定到科;16S rRNA基因?qū)⒆兄婶~樣品鑒定為29個(gè)種類,其中23個(gè)種類鑒定到種、6個(gè)鑒定到屬。

    表1 基于COI和16S rRNA基因鑒定的廈門海域仔稚魚

    兩種基因鑒定結(jié)果基本一致,僅有三處存在分歧。COI基因鑒定為鯡形目sp.(Clupeiformes sp.)的種類,在NCBI和BOLD中與潔白鯡(Escualosathoracata)和印度小公魚(Stolephorusindicus)相似性均≥99%,不能確定具體為哪一個(gè)種類,僅將其鑒定到鯡形目sp.;而NCBI數(shù)據(jù)庫(kù)中沒(méi)有印度小公魚16S rRNA基因,所以此種16S rRNA基因的比對(duì)結(jié)果僅與潔白鯡相似性≥99%,將其鑒定為潔白鯡;兩種基因比對(duì)鑒定結(jié)果不同,本研究?jī)H將其鑒定為鯡形目sp.。COI基因鑒定為小沙丁魚屬sp.(Sardinellasp.)的種類, 在NCBI和BOLD中與黑尾小沙丁魚(S.melanura)、 白腹小沙丁魚(S.albella)、 纟遂鱗小沙丁魚(S.fimbriata)和多耙小沙丁魚(S.jussieui)相似性均大于99%, 無(wú)法確定為哪一個(gè)種類, 僅鑒定到小沙丁魚屬sp.; 而其16S rRNA基因僅與纟遂鱗小沙丁魚相似性大于99%, 其余小沙丁魚屬種類相似性均低于99%, 將其鑒定為纟遂鱗小沙丁魚;本研究將其鑒定為小沙丁魚屬sp.。 COI基因鑒定為小公魚屬sp.(Stolephorussp.)的種類, 在NCBI和BOLD中與短背小公魚(S.bataviensis)和江口小公魚(S.commersonnii)相似性均大于99%, 無(wú)法確定為哪一個(gè)種類, 僅鑒定到小公魚屬sp.; NCBI數(shù)據(jù)庫(kù)中沒(méi)有江口小公魚16S rRNA基因, 所以此種16S rRNA基因的比對(duì)結(jié)果僅與短背小公魚相似性≥99%, 將其鑒定為短背小公魚; 本研究將其鑒定為小公魚屬sp.。

    續(xù)表1

    2.3 遺傳距離

    利用MEGAX軟件,基于K2P模型計(jì)算各分類階元的遺傳距離,結(jié)果見表2。

    各種類COI基因的種內(nèi)遺傳距離在0~0.005 4之間,平均種內(nèi)遺傳距離為0.001 5(SD=0.001 75)。多個(gè)種類的種內(nèi)遺傳距離為0,褐菖鲉的種內(nèi)遺傳距離最大。同屬不同物種間遺傳距離在0.085 1~0.238 5之間,棘鯛屬的黃鰭鯛和黑棘鯛的種間遺傳距離最小,小公魚屬的青帶小公魚和小公魚屬sp.的種間遺傳距離最大,平均種間遺傳距離為0.197 6(SD=0.048 71),約為種內(nèi)平均遺傳距離的132倍,符合Hebert P D N等[21]提出的“10×規(guī)則”,即序列的種間遺傳距離大于平均種內(nèi)遺傳距離的10倍以上時(shí)才能有效鑒別物種。同科不同屬的僅魚箴科的異鱗魚箴屬和下魚箴屬,屬間遺傳距離為0.216 4。目?jī)?nèi)科間遺傳距離最小值出現(xiàn)在鱸形目彈涂魚科和蝦虎魚科,為0.153 1;最大值出現(xiàn)在鱸形目羊魚科和魚喜科,為0.269 7;平均科間遺傳距離為0.238 7(SD=0.018 24)。

    16S rRNA基因種內(nèi)的遺傳距離在0~0.002 0之間,平均為0.000 3(SD=0.000 57),各種類的種內(nèi)遺傳距離明顯小于其COI序列的遺傳距離。多個(gè)種類的種內(nèi)遺傳距離為0,籃子魚屬sp.的種內(nèi)遺傳距離最大。同屬不同物種間遺傳距離在0.026 6~0.164 7之間,棘鯛屬的棘鯛屬sp.和黑棘鯛的種間遺傳距離最小,肩鰓鳚屬的斑點(diǎn)肩鰓鳚和肩鰓鳚屬sp.的種間遺傳距離最大,平均種間遺傳距離為0.089 2(SD=0.045 42),同樣超過(guò)種內(nèi)平均遺傳距離的10倍以上。同科不同屬的魚箴科的異鱗魚箴屬和下魚箴屬,屬間遺傳距離為0.182 3。目?jī)?nèi)科間遺傳距離最小值同樣出現(xiàn)在鱸形目彈涂魚科和蝦虎魚科,為0.059 1;最大值出現(xiàn)在鱸形目鳚科和銀鱸科,為0.227 5;平均科間遺傳距離為0.162 3(SD=0.033 63)。

    由表2可以看出,COI和16S rRNA基因遺傳距離均隨著分類階元的提高而增大,但增幅隨著分類階元的提高越來(lái)越小。

    表2 不同分類階元遺傳距離

    2.4 系統(tǒng)進(jìn)化樹

    基于64條COI基因序列和74條16S rRNA基因序列分別構(gòu)建系統(tǒng)進(jìn)化樹(圖2、圖3)。由圖可以看出,基于COI和16S rRNA基因序列的系統(tǒng)進(jìn)化樹結(jié)果基本一致,在種的水平上,聚類效果清晰,所有物種都分別聚為獨(dú)立的一支從而被有效區(qū)分,同一個(gè)種類的不同個(gè)體也都能聚在同一個(gè)分支。同屬的個(gè)體在基于16S rRNA基因構(gòu)建的系統(tǒng)進(jìn)化樹上聚類效果也很明顯,如棘鯛屬、銀鱸屬、屬、肩鰓鳚屬、小公魚屬的不同種類分別聚為一支,而基于COI基因構(gòu)建的進(jìn)化樹上,肩鰓鳚屬的斑點(diǎn)肩鰓鳚和肩鰓鳚屬sp.并未聚在一支。在科和目的水平上,兩條進(jìn)化樹的聚類效果都不清晰,聚類關(guān)系相互交叉、比較混亂。

    3 討論

    3.1 COI和16S rRNA基因在仔稚魚種類鑒定中的適用性及在系統(tǒng)進(jìn)化分析中的局限性

    Hebert P D N等[22]的研究表明,同一物種內(nèi)的遺傳距離多小于0.020,0.020常作為DNA條形碼物種劃分的閾值,另外指出只有種間遺傳距離大于種內(nèi)遺傳距離的10倍以上時(shí)才能有效鑒別物種[21]。

    本研究各個(gè)種類COI基因的種內(nèi)遺傳距離(0~0.005 4)均小于0.020,同屬種間的遺傳距離(0.085 1~0.238 5)均大于0.020,符合Hebert P D N等提出的種內(nèi)遺傳距離的標(biāo)準(zhǔn),且種間和種內(nèi)遺傳距離分布范圍沒(méi)有交叉,存在明顯的條形碼間隙;平均種間遺傳距離(0.197 6)為種內(nèi)遺傳距離(0.001 5)的約132倍;以上結(jié)果均表明,COI基因可以作為條形碼對(duì)仔稚魚進(jìn)行種類鑒定。16S rRNA基因種內(nèi)遺傳距離(0~0.002 0)遠(yuǎn)小于0.020,同屬種間的遺傳距離(0.026 6~0.164 7)遠(yuǎn)遠(yuǎn)大于種內(nèi)遺傳距離分布范圍,平均種間遺傳距離(0.089 2)大于平均種內(nèi)遺傳距離(0.000 3)的10倍以上,16S rRNA基因同樣也適合用來(lái)進(jìn)行仔稚魚的種類鑒定。另外,基于COI基因和16S rRNA基因構(gòu)建的系統(tǒng)進(jìn)化樹也直觀地反映出兩種基因都能將所有物種區(qū)別開來(lái),印證了兩種基因在仔稚魚種類鑒定中的適用性。

    COI和16S rRNA基因在同科屬間及同目科間的遺傳距離的分布范圍存在重疊,對(duì)基于COI和16S rRNA基因構(gòu)建的系統(tǒng)進(jìn)化樹的分析也表明,在科和目的水平上,聚類關(guān)系不是很清晰。這說(shuō)明,COI和16S rRNA基因不適合用來(lái)分析種以上高分類階元的系統(tǒng)進(jìn)化關(guān)系,這與其他諸多研究者的結(jié)論[12,23-24]相一致。

    3.2 COI和16S rRNA基因?qū)B門海域仔稚魚的鑒定結(jié)果分析

    本研究共選取80尾仔稚魚樣品進(jìn)行DNA提取、COI和16S rRNA基因的擴(kuò)增、測(cè)序、序列比對(duì)與分析。通過(guò)序列比對(duì),COI基因?qū)⒆兄婶~樣品鑒定為26個(gè)種類(4個(gè)種類的樣品PCR擴(kuò)增失敗),16S rRNA基因?qū)⒆兄婶~樣品鑒定為29個(gè)種類(1個(gè)種類的樣品PCR擴(kuò)增失敗),總的來(lái)說(shuō),16S rRNA基因通用引物擴(kuò)增效率高于COI基因。除去擴(kuò)增失敗的種類,兩種基因的鑒定結(jié)果基本一致,存在分歧的3處,多是由于COI基因比對(duì)結(jié)果與多個(gè)近緣種相似而無(wú)法確定具體哪個(gè)種類,而NCBI數(shù)據(jù)庫(kù)中16S rRNA基因序列數(shù)據(jù)不足,16S rRNA基因比對(duì)結(jié)果僅與其中一個(gè)種類相似性大于99%。目前,DNA條形碼數(shù)據(jù)庫(kù)并不全面,尤其是一些非經(jīng)濟(jì)種和少見種的COI和16S rRNA基因序列還十分欠缺。另外,數(shù)據(jù)庫(kù)中存在一些錯(cuò)誤鑒定后提交的序列,如果未能加以甄別則會(huì)干擾后續(xù)研究者序列的比對(duì)和物種鑒定[25]。因此,DNA條形碼基因序列庫(kù)的完整和準(zhǔn)確是物種鑒定的前提條件[26],除不斷補(bǔ)充DNA條形碼數(shù)據(jù)庫(kù)外,還必須嚴(yán)格開展序列的審查工作,保證數(shù)據(jù)庫(kù)中序列的質(zhì)量,從而更好地發(fā)揮其在物種鑒定中的作用。

    本研究中鑒定為鯡形目sp.、小沙丁魚屬sp.、小公魚屬sp.、鲾屬sp.、銀鱸屬sp1.、銀鱸屬sp2.、棘鯛屬sp.、籃子魚屬sp.的種類,經(jīng)COI或16S rRNA基因比對(duì)后,均與數(shù)據(jù)庫(kù)中多個(gè)種類高度匹配,相似性均大于99%,無(wú)法鑒定到種的水平。說(shuō)明有些魚類不僅形態(tài)上非常相似、無(wú)法區(qū)分,序列差異也較小,利用COI或16S rRNA基因均無(wú)法準(zhǔn)確鑒定,這也是利用DNA條形碼技術(shù)進(jìn)行種類鑒定的難點(diǎn)[27],針對(duì)這種情況,仍需要結(jié)合更多分子標(biāo)記(如線粒體Cytb基因、核基因分子標(biāo)記RAG2等)作進(jìn)一步研究。

    綜上所述,COI和16S rRNA基因均可以作為條形碼對(duì)形態(tài)學(xué)上難以區(qū)分的仔稚魚進(jìn)行種類鑒定,兩種基因鑒定結(jié)果基本一致,且多數(shù)能鑒定到種的水平,兩種基因結(jié)合使用可以相互印證或補(bǔ)充,從而提高鑒定結(jié)果的準(zhǔn)確性。由于本研究采集到的樣品數(shù)量有限,有的種類僅1條序列,相關(guān)數(shù)據(jù)不夠全面,下一步應(yīng)采集更多季節(jié)或月份的仔稚魚樣品,對(duì)DNA條形碼技術(shù)在仔稚魚種類鑒定中的應(yīng)用進(jìn)行更充分的研究和探討。

    91麻豆精品激情在线观看国产 | 考比视频在线观看| 丰满饥渴人妻一区二区三| 日韩制服丝袜自拍偷拍| 午夜福利一区二区在线看| 亚洲精品粉嫩美女一区| 欧美黑人欧美精品刺激| 国产在线观看jvid| 国产精品麻豆人妻色哟哟久久| 天天躁夜夜躁狠狠躁躁| 男女午夜视频在线观看| 久久精品国产a三级三级三级| 亚洲精品中文字幕在线视频| 韩国精品一区二区三区| 欧美精品高潮呻吟av久久| 午夜免费观看性视频| 国产精品久久久人人做人人爽| 老汉色av国产亚洲站长工具| 国产精品自产拍在线观看55亚洲 | 亚洲国产av新网站| 亚洲精品乱久久久久久| 国产亚洲欧美精品永久| 搡老乐熟女国产| 91精品伊人久久大香线蕉| 高清欧美精品videossex| 亚洲精品一卡2卡三卡4卡5卡 | 捣出白浆h1v1| 91av网站免费观看| 最新的欧美精品一区二区| 久久香蕉激情| 黑人操中国人逼视频| 嫩草影视91久久| 日韩三级视频一区二区三区| 国产麻豆69| 啦啦啦视频在线资源免费观看| 老汉色∧v一级毛片| 亚洲精品一二三| 1024视频免费在线观看| 精品一区二区三区av网在线观看 | 黄色 视频免费看| 亚洲男人天堂网一区| 亚洲一区二区三区欧美精品| 精品亚洲成a人片在线观看| 男人添女人高潮全过程视频| 亚洲精品自拍成人| 91大片在线观看| 一二三四在线观看免费中文在| 日韩电影二区| 一级,二级,三级黄色视频| 亚洲国产中文字幕在线视频| 秋霞在线观看毛片| 美女高潮喷水抽搐中文字幕| 如日韩欧美国产精品一区二区三区| 亚洲av电影在线观看一区二区三区| 久久久久国产精品人妻一区二区| 99精品欧美一区二区三区四区| 国产精品偷伦视频观看了| 成人影院久久| 欧美日韩亚洲高清精品| 天天添夜夜摸| av网站在线播放免费| 久久久久网色| 欧美午夜高清在线| 麻豆国产av国片精品| 久久久久久久精品精品| 久久人人97超碰香蕉20202| 欧美激情极品国产一区二区三区| 精品一区二区三区四区五区乱码| 黑人操中国人逼视频| 国产深夜福利视频在线观看| 午夜福利在线免费观看网站| 爱豆传媒免费全集在线观看| 老司机亚洲免费影院| 人人妻,人人澡人人爽秒播| 欧美黑人精品巨大| 精品免费久久久久久久清纯 | 久9热在线精品视频| 久久久久久久久免费视频了| 国产成人欧美| 国产免费现黄频在线看| 日本撒尿小便嘘嘘汇集6| 王馨瑶露胸无遮挡在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 一级片免费观看大全| 国产麻豆69| 欧美国产精品一级二级三级| 黄片播放在线免费| 美女扒开内裤让男人捅视频| 免费久久久久久久精品成人欧美视频| 99久久精品国产亚洲精品| 欧美黄色淫秽网站| 亚洲色图综合在线观看| 国产99久久九九免费精品| 热99国产精品久久久久久7| 乱人伦中国视频| av有码第一页| 精品熟女少妇八av免费久了| 淫妇啪啪啪对白视频 | 超碰97精品在线观看| 国产欧美日韩一区二区三区在线| 在线观看免费视频网站a站| 91成年电影在线观看| cao死你这个sao货| 中文字幕色久视频| 丁香六月欧美| 亚洲欧美日韩另类电影网站| 亚洲全国av大片| kizo精华| 国产亚洲精品久久久久5区| 成年人黄色毛片网站| 亚洲国产欧美在线一区| 热re99久久国产66热| 久久人人爽人人片av| 日韩精品免费视频一区二区三区| 亚洲少妇的诱惑av| 国产国语露脸激情在线看| 视频区图区小说| 精品第一国产精品| 国产成人免费无遮挡视频| 欧美日韩黄片免| 久久久国产欧美日韩av| 在线观看人妻少妇| 黄片小视频在线播放| www.999成人在线观看| 中文字幕另类日韩欧美亚洲嫩草| 亚洲色图 男人天堂 中文字幕| 亚洲国产欧美在线一区| 啦啦啦在线免费观看视频4| 欧美精品亚洲一区二区| 国产深夜福利视频在线观看| 99国产精品免费福利视频| 亚洲激情五月婷婷啪啪| 日韩大片免费观看网站| 一级毛片女人18水好多| 欧美日韩黄片免| 一本一本久久a久久精品综合妖精| 国产福利在线免费观看视频| 日韩制服骚丝袜av| 一级黄色大片毛片| 国产精品一区二区在线观看99| 久久久久国产一级毛片高清牌| 97人妻天天添夜夜摸| 人妻一区二区av| 久久精品久久久久久噜噜老黄| 人妻人人澡人人爽人人| 欧美老熟妇乱子伦牲交| 日韩大片免费观看网站| 国产一级毛片在线| 亚洲性夜色夜夜综合| 午夜久久久在线观看| 妹子高潮喷水视频| 满18在线观看网站| 久久99热这里只频精品6学生| 精品亚洲乱码少妇综合久久| 亚洲精品日韩在线中文字幕| 久久精品国产综合久久久| 国产xxxxx性猛交| 十八禁人妻一区二区| 在线观看免费高清a一片| 国产在线观看jvid| 十八禁人妻一区二区| 久久 成人 亚洲| 国产精品久久久av美女十八| 成人黄色视频免费在线看| 欧美日韩国产mv在线观看视频| 首页视频小说图片口味搜索| 亚洲精品一卡2卡三卡4卡5卡 | 亚洲色图 男人天堂 中文字幕| 精品卡一卡二卡四卡免费| 巨乳人妻的诱惑在线观看| 欧美亚洲 丝袜 人妻 在线| 国产免费视频播放在线视频| 精品视频人人做人人爽| 菩萨蛮人人尽说江南好唐韦庄| 免费日韩欧美在线观看| 91麻豆精品激情在线观看国产 | 一区福利在线观看| 国产成人一区二区三区免费视频网站| 老熟妇乱子伦视频在线观看 | 久久性视频一级片| 亚洲全国av大片| 日韩免费高清中文字幕av| 亚洲欧洲精品一区二区精品久久久| av国产精品久久久久影院| 色老头精品视频在线观看| 欧美+亚洲+日韩+国产| 国产日韩欧美在线精品| 咕卡用的链子| 青草久久国产| 中文字幕人妻熟女乱码| 2018国产大陆天天弄谢| 久久 成人 亚洲| 欧美日韩亚洲高清精品| www.999成人在线观看| 在线观看免费日韩欧美大片| 国产人伦9x9x在线观看| 777久久人妻少妇嫩草av网站| 日韩精品免费视频一区二区三区| 国产福利在线免费观看视频| 大码成人一级视频| 欧美 亚洲 国产 日韩一| 搡老岳熟女国产| 亚洲欧美精品自产自拍| 国产男人的电影天堂91| 美女国产高潮福利片在线看| 久久国产精品男人的天堂亚洲| av网站免费在线观看视频| 丁香六月欧美| 亚洲国产精品999| 成人国语在线视频| 美女视频免费永久观看网站| 电影成人av| 真人做人爱边吃奶动态| 人妻 亚洲 视频| 日本一区二区免费在线视频| 韩国高清视频一区二区三区| 美女脱内裤让男人舔精品视频| av天堂久久9| 涩涩av久久男人的天堂| 久久精品久久久久久噜噜老黄| a 毛片基地| 久久久国产精品麻豆| 欧美国产精品一级二级三级| 亚洲第一av免费看| 一区二区av电影网| 侵犯人妻中文字幕一二三四区| 91成年电影在线观看| 成人影院久久| 一区福利在线观看| 中文字幕色久视频| av天堂在线播放| 18在线观看网站| 亚洲欧美清纯卡通| 久久精品人人爽人人爽视色| 国产片内射在线| 国产深夜福利视频在线观看| 免费观看a级毛片全部| 亚洲午夜精品一区,二区,三区| 欧美激情极品国产一区二区三区| 中国国产av一级| 丁香六月天网| 精品一区二区三区四区五区乱码| 黄色 视频免费看| 精品福利永久在线观看| 美国免费a级毛片| 十八禁网站免费在线| 精品久久蜜臀av无| 午夜精品国产一区二区电影| 亚洲五月婷婷丁香| 亚洲精品一卡2卡三卡4卡5卡 | 亚洲综合色网址| 午夜精品国产一区二区电影| 亚洲五月婷婷丁香| 亚洲欧美日韩另类电影网站| 51午夜福利影视在线观看| 久久 成人 亚洲| 国产91精品成人一区二区三区 | 亚洲一区中文字幕在线| 狠狠婷婷综合久久久久久88av| 黄色怎么调成土黄色| 精品一品国产午夜福利视频| 777米奇影视久久| 日韩制服丝袜自拍偷拍| 一区二区三区精品91| av不卡在线播放| 国产精品九九99| 高清在线国产一区| 亚洲精品久久成人aⅴ小说| 日韩一卡2卡3卡4卡2021年| 一区二区三区激情视频| 一级,二级,三级黄色视频| 一本大道久久a久久精品| 成年人免费黄色播放视频| 欧美国产精品一级二级三级| 久久精品熟女亚洲av麻豆精品| 热99re8久久精品国产| 日韩中文字幕欧美一区二区| 视频在线观看一区二区三区| 国产精品一区二区精品视频观看| 国产亚洲精品久久久久5区| 真人做人爱边吃奶动态| 日本精品一区二区三区蜜桃| 久久影院123| 国产真人三级小视频在线观看| 热99re8久久精品国产| 欧美日韩黄片免| 亚洲精品自拍成人| 水蜜桃什么品种好| 亚洲一码二码三码区别大吗| 免费日韩欧美在线观看| 51午夜福利影视在线观看| 亚洲av片天天在线观看| 99re6热这里在线精品视频| 成年女人毛片免费观看观看9 | 亚洲欧美清纯卡通| 美女脱内裤让男人舔精品视频| 窝窝影院91人妻| 精品乱码久久久久久99久播| 国产精品一区二区在线不卡| 国产精品久久久人人做人人爽| 精品欧美一区二区三区在线| 国产在线观看jvid| 国产区一区二久久| 涩涩av久久男人的天堂| 日本一区二区免费在线视频| 999久久久国产精品视频| 亚洲五月婷婷丁香| 亚洲熟女毛片儿| a级毛片在线看网站| 免费日韩欧美在线观看| 亚洲国产精品一区二区三区在线| 又黄又粗又硬又大视频| 中文字幕人妻丝袜制服| 一区二区日韩欧美中文字幕| 午夜成年电影在线免费观看| 两个人免费观看高清视频| 狂野欧美激情性xxxx| 日本五十路高清| 在线观看免费午夜福利视频| 在线亚洲精品国产二区图片欧美| 美女主播在线视频| 老熟妇仑乱视频hdxx| 久久九九热精品免费| 久久久久国产精品人妻一区二区| 精品免费久久久久久久清纯 | 久久国产亚洲av麻豆专区| 国产欧美日韩精品亚洲av| 国产高清国产精品国产三级| 黄片小视频在线播放| 精品一区在线观看国产| 国产免费av片在线观看野外av| 亚洲avbb在线观看| 国产淫语在线视频| 美女主播在线视频| 秋霞在线观看毛片| 99香蕉大伊视频| 国产精品1区2区在线观看. | 国产又爽黄色视频| 人妻久久中文字幕网| 国产精品麻豆人妻色哟哟久久| 99久久国产精品久久久| 人人妻,人人澡人人爽秒播| 视频区图区小说| 三上悠亚av全集在线观看| 欧美成狂野欧美在线观看| 亚洲熟女毛片儿| 日韩大码丰满熟妇| 亚洲三区欧美一区| 两性夫妻黄色片| 精品少妇一区二区三区视频日本电影| 色94色欧美一区二区| 久久影院123| 久久久国产成人免费| e午夜精品久久久久久久| av在线播放精品| 在线亚洲精品国产二区图片欧美| www.av在线官网国产| 极品人妻少妇av视频| 国产成人精品久久二区二区免费| 一本大道久久a久久精品| 99九九在线精品视频| 一区二区三区激情视频| 美女高潮喷水抽搐中文字幕| 欧美亚洲 丝袜 人妻 在线| 91麻豆av在线| 亚洲欧美清纯卡通| 深夜精品福利| 777米奇影视久久| 午夜两性在线视频| 少妇裸体淫交视频免费看高清 | 免费在线观看完整版高清| 亚洲综合色网址| 久久人人97超碰香蕉20202| 又大又爽又粗| 一本大道久久a久久精品| 欧美黄色片欧美黄色片| 久久久久国产一级毛片高清牌| 日韩中文字幕欧美一区二区| 精品免费久久久久久久清纯 | 国产免费福利视频在线观看| xxxhd国产人妻xxx| avwww免费| bbb黄色大片| 亚洲熟女毛片儿| 亚洲七黄色美女视频| 亚洲精品av麻豆狂野| 人妻 亚洲 视频| e午夜精品久久久久久久| 婷婷色av中文字幕| 人成视频在线观看免费观看| 99久久99久久久精品蜜桃| 免费在线观看完整版高清| 亚洲第一欧美日韩一区二区三区 | 欧美日韩精品网址| 叶爱在线成人免费视频播放| 免费在线观看日本一区| 久久亚洲精品不卡| √禁漫天堂资源中文www| 精品人妻在线不人妻| 欧美精品av麻豆av| 十八禁网站免费在线| 日日爽夜夜爽网站| 叶爱在线成人免费视频播放| 亚洲国产精品成人久久小说| 国产成人影院久久av| 丰满饥渴人妻一区二区三| 后天国语完整版免费观看| 国产一级毛片在线| 老司机午夜十八禁免费视频| 少妇精品久久久久久久| 欧美另类亚洲清纯唯美| 热99久久久久精品小说推荐| 一级毛片女人18水好多| 天天添夜夜摸| 曰老女人黄片| 一区二区三区精品91| 美女视频免费永久观看网站| 国产一区二区在线观看av| 精品亚洲成a人片在线观看| 久久久久国产一级毛片高清牌| avwww免费| 成人影院久久| 国产一区二区激情短视频 | 男女边摸边吃奶| 日韩中文字幕欧美一区二区| 久久精品久久久久久噜噜老黄| av天堂在线播放| 精品欧美一区二区三区在线| 欧美午夜高清在线| 啦啦啦在线免费观看视频4| 人人妻人人澡人人爽人人夜夜| 丰满少妇做爰视频| 国产精品久久久久久人妻精品电影 | 国产97色在线日韩免费| 久久精品熟女亚洲av麻豆精品| 亚洲欧洲精品一区二区精品久久久| 深夜精品福利| 精品国产国语对白av| 亚洲精品中文字幕一二三四区 | 久久中文看片网| 脱女人内裤的视频| 久久人妻福利社区极品人妻图片| 国产精品一区二区精品视频观看| 亚洲天堂av无毛| 韩国精品一区二区三区| 成在线人永久免费视频| 不卡一级毛片| 国产精品av久久久久免费| 亚洲 欧美一区二区三区| 国产一区有黄有色的免费视频| 99久久精品国产亚洲精品| 中文字幕制服av| 亚洲色图 男人天堂 中文字幕| 国产亚洲av高清不卡| 亚洲av欧美aⅴ国产| av电影中文网址| videosex国产| 国精品久久久久久国模美| 蜜桃在线观看..| 免费久久久久久久精品成人欧美视频| av在线老鸭窝| av片东京热男人的天堂| 亚洲专区中文字幕在线| 成年人黄色毛片网站| 一边摸一边做爽爽视频免费| 母亲3免费完整高清在线观看| av不卡在线播放| 丝袜美腿诱惑在线| 99国产精品99久久久久| 在线天堂中文资源库| 人人妻人人添人人爽欧美一区卜| 80岁老熟妇乱子伦牲交| 亚洲精品国产一区二区精华液| 国产又爽黄色视频| 亚洲第一欧美日韩一区二区三区 | netflix在线观看网站| 国产免费现黄频在线看| 啦啦啦视频在线资源免费观看| 99国产精品免费福利视频| 久久午夜综合久久蜜桃| 人妻 亚洲 视频| 中文精品一卡2卡3卡4更新| av天堂在线播放| 天天躁狠狠躁夜夜躁狠狠躁| 别揉我奶头~嗯~啊~动态视频 | 久久人妻熟女aⅴ| 91九色精品人成在线观看| 精品国产一区二区三区久久久樱花| 夫妻午夜视频| 国产成人系列免费观看| 女人精品久久久久毛片| av在线app专区| 99热网站在线观看| 精品亚洲乱码少妇综合久久| 视频区图区小说| 久久久久久久国产电影| 69精品国产乱码久久久| 欧美日韩亚洲高清精品| 国产欧美日韩一区二区精品| 国产精品久久久久久精品古装| 超色免费av| 黄色视频在线播放观看不卡| cao死你这个sao货| 日韩人妻精品一区2区三区| 色老头精品视频在线观看| 香蕉丝袜av| 美国免费a级毛片| 国产精品av久久久久免费| 精品久久久精品久久久| 精品国产乱码久久久久久男人| 亚洲专区字幕在线| 成年美女黄网站色视频大全免费| 又黄又粗又硬又大视频| 天堂中文最新版在线下载| 欧美性长视频在线观看| 亚洲av日韩在线播放| www日本在线高清视频| 精品福利永久在线观看| av免费在线观看网站| 亚洲精品久久久久久婷婷小说| 国产精品一二三区在线看| 中文欧美无线码| 亚洲精品中文字幕在线视频| 中文字幕人妻熟女乱码| 91字幕亚洲| e午夜精品久久久久久久| 波多野结衣av一区二区av| 可以免费在线观看a视频的电影网站| av在线老鸭窝| 成人三级做爰电影| 精品一品国产午夜福利视频| av一本久久久久| 一本一本久久a久久精品综合妖精| 人人妻人人添人人爽欧美一区卜| 亚洲国产日韩一区二区| 精品高清国产在线一区| 天堂中文最新版在线下载| 97人妻天天添夜夜摸| 国产精品偷伦视频观看了| 纯流量卡能插随身wifi吗| 亚洲国产av影院在线观看| 亚洲欧美色中文字幕在线| 国产精品免费大片| 少妇被粗大的猛进出69影院| 精品第一国产精品| 亚洲精品美女久久av网站| 精品少妇内射三级| 99国产综合亚洲精品| 色婷婷久久久亚洲欧美| 中亚洲国语对白在线视频| 国产亚洲一区二区精品| 一区在线观看完整版| 狠狠精品人妻久久久久久综合| 老司机午夜十八禁免费视频| 亚洲国产毛片av蜜桃av| av网站在线播放免费| 久久天堂一区二区三区四区| 黄色怎么调成土黄色| 成人亚洲精品一区在线观看| 亚洲专区字幕在线| 亚洲黑人精品在线| 天堂俺去俺来也www色官网| 欧美日韩一级在线毛片| 国产亚洲精品第一综合不卡| 人妻久久中文字幕网| 亚洲视频免费观看视频| xxxhd国产人妻xxx| 一个人免费在线观看的高清视频 | 亚洲国产欧美日韩在线播放| 一区二区av电影网| 国产在线观看jvid| 丰满人妻熟妇乱又伦精品不卡| 亚洲第一青青草原| 亚洲精品av麻豆狂野| 久久人妻福利社区极品人妻图片| 下体分泌物呈黄色| 91国产中文字幕| 黑人巨大精品欧美一区二区蜜桃| 成人国语在线视频| 久久久精品免费免费高清| 91老司机精品| 久久国产精品大桥未久av| 亚洲一区中文字幕在线| 大香蕉久久成人网| 成人国产一区最新在线观看| av又黄又爽大尺度在线免费看| 国产又色又爽无遮挡免| 人人妻人人添人人爽欧美一区卜| 人人妻人人澡人人爽人人夜夜| 老司机亚洲免费影院| 精品福利永久在线观看| 超碰97精品在线观看| 日本撒尿小便嘘嘘汇集6| 日韩精品免费视频一区二区三区| 久久国产精品人妻蜜桃| 亚洲av日韩精品久久久久久密| 欧美性长视频在线观看| 十八禁人妻一区二区| 91字幕亚洲| av网站在线播放免费| 欧美黑人欧美精品刺激| 丰满少妇做爰视频| av超薄肉色丝袜交足视频| 女人被躁到高潮嗷嗷叫费观| 精品久久蜜臀av无| 99精品欧美一区二区三区四区| 热re99久久国产66热| 欧美xxⅹ黑人| 激情视频va一区二区三区| 精品一区二区三区四区五区乱码| 亚洲国产看品久久| 亚洲欧美日韩高清在线视频 | 桃红色精品国产亚洲av| 嫩草影视91久久| 国产精品免费大片|